Skip to content

Instantly share code, notes, and snippets.

@cjtu
Created February 25, 2019 23:33
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save cjtu/d4a2cf4687564f871528d60ee0015e29 to your computer and use it in GitHub Desktop.
Save cjtu/d4a2cf4687564f871528d60ee0015e29 to your computer and use it in GitHub Desktop.
PHY 530 homework 2
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "import numpy as np\nimport pandas as pd\nimport matplotlib.pyplot as plt\n\n\nn1 = 1 # vacuum (~air)\ndn2 = {\"glass\": 1.5, \"diamond\": 2.4, \"graphite\": 2.7}\n\ndef Fresnell_TE(n1, n2, theta):\n \"\"\"\n Return intensity of transverse electric component of the Fresnell\n equations (perpendicular reflected / S-polarized).\n \"\"\"\n n = n2 / n1\n num = np.cos(theta) - np.sqrt(n**2 - np.sin(theta)**2)\n denom = np.cos(theta) + np.sqrt(n**2 - np.sin(theta)**2)\n I = (num/denom)**2\n return I\n \n \ndef Fresnell_TM(n1, n2, theta):\n \"\"\"\n Return intensity of transverse magnetic component of the Fresnell\n equations (parallel reflected / P-polarized).\n \"\"\"\n n = n2 / n1\n num = -(n**2) * np.cos(theta) + np.sqrt(n**2 - np.sin(theta)**2)\n denom = (n**2) * np.cos(theta) + np.sqrt(n**2 - np.sin(theta)**2)\n I = (num/denom)**2\n return I\n\ndef main():\n dTE = {}\n dTM = {}\n theta = np.radians(np.linspace(0, 90, (90*20)+1))\n for n2 in (\"glass\", \"diamond\", \"graphite\"):\n dTE[n2] = Fresnell_TE(n1, dn2[n2], theta)\n dTM[n2] = Fresnell_TM(n1, dn2[n2], theta)\n\n # Save as DataFrame\n df_te = pd.DataFrame(dTE)\n df_tm = pd.DataFrame(dTM)\n df_te['theta'] = np.degrees(theta)\n df_tm['theta'] = np.degrees(theta)\n\n # Init plots\n fig, axs = plt.subplots(nrows=2, ncols=1, figsize=(8,8))\n plt.subplots_adjust(hspace=0.4)\n\n df_te.plot(x='theta', ax=axs[0])\n axs[0].set_ylabel(\"intensity\")\n axs[0].set_title(\"Transverse Electric\")\n\n df_tm.plot(x='theta', ax=axs[1])\n axs[1].set_ylabel(\"intensity\")\n axs[1].set_title(\"Transverse Magnetic\")",
"execution_count": 74,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "main()",
"execution_count": 75,
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAfUAAAHwCAYAAAC/hfaiAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3XlYVdX6wPHvYp5BwBFRcUScFaWcyynNrGww01KvQ3PdrO6tbpY/m61bWdmg3qxMK60sh8yynM0JZ3BCBZlEAZnHw1m/PzYgIiIohwP4fp7nPOecvdfe+z0n4z1r7TUorTVCCCGEqP1srB2AEEIIIaqGJHUhhBCijpCkLoQQQtQRktSFEEKIOkKSuhBCCFFHSFIXQggh6ghJ6kIIq1JKTVFKbaima01QSq2pjmsJYQ2S1IWoAkqpjBIPs1Iqu8T7cdaOz9qUUjGlvpMMpdQHVXj+1kqpK066obX+Sms9vKquK0RNY2ftAISoC7TWbkWvlVKRwBSt9brLlVdK2WmtTdUR29WyQIzDtdYbqvB8lVIbvnMhrpXU1IWoBkqp15RS3yulvlVKpQPjlVI3KqW2K6VSlFLxSqkPlVL2heXtlFJaKfWQUipCKXVeKfVhifO1VUptUkqlKqUSlVJLCrcvUEq9Veraq5VSTxa+bqqUWq6UOqeUOqWUeuwKMd6glNqjlEpTSiUopd4pUb5Pifj3KaX6V9F3FaSUWqeUSlZKHVFK3VVin4tS6n2l1OnCz75JKeUIbCrcX9QK0LOwWX9T4feaDLxUuqlfKdWpxLXOKKX+VRWfQQhrkaQuRPW5E1gCeALfAybgKcAX6APcAjxU6pgRQA+gG0aSHVy4/XVgNVAPaArMLdy+BLhPKaUAlFI+wM3A90opW2AVsAvwA4YAzymlBpUT40fAO1prD6A18EPhef2BFcArgDfwPPBT4fWumlLKHfgD+BpoAIwD5iml2hUWeR/oDIQUXvdFwAz0B6PFpPCxq7B8b+AwUB94u9S1PIF1wEqgMdAW2HAt8QthbZLUhag+W7TWK7XWZq11ttZ6l9Z6h9bapLU+CcwDBpQ65k2tdarWOhIj4XQt3J4PtAAaa61ztNZbC7dvAOyBGwvf3wts1lonADcAHlrrN7TWeVrrCOB/wH2Xi7HwOm2UUj5a63St9Y7Ccg8CK7TWawvL/gbsx/hhcjmrCmv1RY9JZZQZBRzTWn9d+L2EAj8Ddxf+KJkIPKm1jtdaF2itt2it88u55mmt9aeFZbPLuFa01nqO1jpXa52mtd5ZzrmEqPEkqQtRfaJLvlFKBRY2jZ9RSqUBszBq7SWdKfE6Cyi6d/8MRvLerZQ6qJSaAKC1NmPUsMcWlrsfWFz4ujnQrGRiBf4FNLpcjMAkIAg4qpTaqZQaUeJcY0ud6wagSTmff6TW2qvEY2EZZZoDfUqddwxGTboh4ACcKOcapZX+PCX5AxGVOJcQNZ50lBOi+pTunf05sB0Yo7XOUEo9C4ys0Im0jgemABTey/5DKbVJa30K+BajVvwe0B1YXnhYNHBca92+ojFqrY9iNOfbAPcAPyql6hWea6HW+pGKxFsJ0cCfZfVQL6yp5wGtgLDy4q7A9qJr3Xk1QQpRU0lNXQjrcQdSgUylVHsuvZ9+WUqpe5VSfoVvUzCSVwFA4f3kVIzm/F+11mmF5f4G8pRSzyilnJRStoUdxXqUc50HlFK+hS0AqYXXMQOLgDuVUkMKz+OklLpJKVVeTb0iVgAdlFL3K6XsCx+9lFLttNYFwJfAB0qpRoXX7VPYufAsoJVSLSt5rWZKqceVUg5KKQ+lVK9rjF8Iq5KkLoT1PANMANIxau3fV+LYEGCXUioT+Al4TGt9usT+b4HBGJ3eACgczjUC6AVEAomF1/Uo5zojgMOFveHfxWhVyCu8x38nMAM4B5wu/Dzl/U1Zoy4ep76sdAGtdSowDBgPxGPcfngTcCws8jRGx7dQIBl4A1Ba6/TCcjsKm+2Dy4mj5LWGAHdh/Cg4xqV9GoSoVZTWV5yvQQghhBC1gNTUhRBCiDpCkroQQghRR0hSF0IIIeoISepCCCFEHSFJXQghhKgjat3kM76+vrpFixbWDkMIIYSoFqGhoYla6/oVKVvrknqLFi3YvXu3tcMQQgghqoVSKqqiZaX5XQghhKgjJKkLIYQQdYTFkrpS6gul1Fml1KHL7FdKqQ+VUhFKqQNKqe6WikUIIYS4HljynvqXwMfA15fZPxxoU/gIAT4tfK60/Px8YmJiyMnJuZrDr2tOTk40bdoUe3t7a4cihBDiGlksqWutNymlWpRT5Hbga21MPr9dKeWllGpcuKRkpcTExODu7k6LFi1QSl1lxNcfrTVJSUnExMQQEBBg7XCEEEJcI2veU/fDWM+4SEzhtkrLycnBx8dHEnolKaXw8fGRFg4hhKgiX/36LtP/N4wsK/1dtWZSLysDl7lknFJqmlJqt1Jq97lz58o+mST0qyLfmxBCVI19x1byccKXHDPHYco7b5UYrJnUYwD/Eu+bAnFlFdRaz9NaB2utg+vXr9D4+xpj4sSJ/PDDD9YOQwghhAWFx27nka0v4msyc3ezt/HwaGyVOKyZ1FcADxb2gr8BSL2a++lCCCGENR1LDOehPx7C3WzCPe5+7howxGqxWKyjnFLqW2Ag4KuUigFeAewBtNafAb8CI4AIIAuYZKlYqsurr77K4sWL8ff3x9fXlx49ely0f9asWaxcuZLs7Gx69+7N559/jlKKDz/8kM8++ww7OzuCgoL47rvv2LhxI0899RRgNJFv2rQJd3d3a3wsIYQQl3HifARTfx2PQ0EebWKG4Nv1DtydrDeayJK938deYb8GHqvq6/7fyjDC49Kq9JxBTTx45bYO5ZbZvXs3P/74I3v37sVkMtG9e/dLkvrjjz/Oyy+/DMADDzzAqlWruO2223jrrbc4deoUjo6OpKSkAPDuu+8yd+5c+vTpQ0ZGBk5OTlX6mYQQQlybU6mnmLz6fmzyc/in7scT2cP4s3cLq8YkM8pVkS1btnD77bfj7OyMu7s7t9122yVl1q9fT0hICJ06deKvv/4iLCwMgM6dOzNu3Di++eYb7OyM31l9+vRh+vTpfPjhh6SkpBRvF0IIYX2n004zZfV4dF4G85zbMyv+Lm4ObECr+m5WjavOZYor1agtxWh4uLycnBweffRRdu/ejb+/PzNnziweSrZ69Wo2bdrEihUrePXVVwkLC+P555/n1ltv5ddff+WGG25g3bp1BAYGVsdHEUIIUY6Y9Bgmr3mQvJwUvsj34kDQG5w7EMHkvtaf70Nq6lWkb9++rFy5kpycHDIyMli9evVF+4sSuK+vLxkZGcU94s1mM9HR0dx0003Mnj2blJQUMjIyOHHiBJ06deLf//43wcHBHDlypNo/kxBCiIvFZ8Qz5beJZGUlMj+tgNb3LWX+9gQCG7nTu5WPtcOrezV1a+nZsyejRo2iS5cuNG/enODgYDw9PYv3e3l5MXXqVDp16kSLFi3o2bMnAAUFBYwfP57U1FS01jz99NN4eXkxY8YM1q9fj62tLUFBQQwfPtxaH00IIQSQkJnAP36bSFpmAvMT0wl8YBVbzjpx5Ew679zduUbM+6Gu1Gxc0wQHB+vS66kfPnyY9u3bWymiCzIyMnBzcyMrK4v+/fszb948unev+evU1JTvTwghaqozmWeYsnYyienRzItLoPM9S6D1ICYu3Mmh2DS2Pn8Tjna2Frm2UipUax1ckbJSU69C06ZNIzw8nJycHCZMmFArEroQQojyxWXEMXntZM5nxPNZXDydh70HrQcRcTadDUfPMX1IW4sl9MqSpF6FlixZYu0QhBBCVKHo9GimrJ1CetY55sfG0unGp6H7AwD8b0skjnY2jAtpZuUoL5CkLoQQQpQhKi2KyWsnk5ObyoKY0wS1vxsGvgBAcmYeP+2JYXT3pvi4OVo50gskqQshhBClnEw9yZS1UzDlZ/G/6NO08+sNt30IhZ3hFv0dRa7JzOS+LawbaCmS1IUQQogSIs5HMOX3KWA28UVsHK09W8GYRWDnAEBWnokvt51icPsGtG5Qs6bvlnHqQgghRKGjyUf5x9p/YKPNfBEXT2sHTxj/IzhdGKL83c5ozmfl88jA1laMtGxSU7eQmTNn4ubmRlpaGv3792fw4MFWjScyMpKRI0dy6NAhq8YhhBA11YFzB3hk3SM42zrwv/hzNDcrmPgzlFhGNc9kZsHmk/QK8KZH83pWjLZsktQtbNasWdYOQQghxBVsj9/Ok389ia+jN/POJtE0MwUmrgKfVheVW7E/jrjUHF4f3clKkZZPmt+r0Ouvv067du0YPHgwR48eBWDixInFU8LOmjWLnj170rFjR6ZNm1Y8X/zAgQN5+umn6d+/P+3bt2fXrl2MHj2aNm3a8NJLLxWf/7333qNjx4507NiRDz74ADBq4O3bt2fq1Kl06NCBoUOHkp2dDUBoaChdunThxhtvZO7cudX5VQghRK3xZ9SfPLruUZq6NuGrlDyaJkXC2CXQpOtF5cxmzWcbT9C+sQcD29a3TrBXUPdq6muehzMHq/acjTrB8LfKLRIaGsp33313VUuvAjg4OLBp0ybmzJnD7bffTmhoKN7e3rRq1Yqnn36ayMhIFi5cyI4dO9BaExISwoABA6hXrx7Hjx/n22+/Zf78+dx77738+OOPjB8/nkmTJvHRRx8xYMAAnnvuuar9ToQQog74JeIXXt72Mh19OvBJcjaesXvhnq8goP8lZf84nEDE2Qw+HNutRkwJWxapqVeRzZs3c+edd+Li4oKHhwejRo26pMzlll4Fist36tSJDh060LhxYxwdHWnZsiXR0dFs2bKFO++8E1dXV9zc3Bg9ejSbN28GICAggK5djV+UPXr0IDIyktTUVFJSUhgwYABg/IgQQghxweLDi3lp60v0atiT+dmOeJ74C0a+D0GX/v3WWvPJhhM083ZhRMdGVoi2YupeTf0KNWpLKu+XW3lLrwI4OhqTF9jY2BS/LnpvMpnKXdq1ZHlbW1uys7PRWtfYX5JCCGFNWms+2/8Zn+z/hEHNBjE7LR+HQ9/AoJehx8Qyj9l+Mpn90Sm8dkdH7Gxrbn245kZWy/Tv35/ly5eTnZ1Neno6K1euvGj/5ZZercz5f/75Z7KyssjMzGT58uX069fvsuW9vLzw9PRky5YtACxevLiSn0gIIeoeszYze9dsPtn/CaNajeLdbHsc9n4D/Z6Ffs9c9riP/jpOfXdH7u7RtBqjrby6V1O3ku7duzNmzBi6du1K8+bNL0m4l1t6tTLnnzhxIr169QJgypQpdOvWjcjIyMses3DhQv7xj3/g4uLCsGHDKv2ZhBCiLskryOPFLS+yNnIt49uP57nULGx2fgw3PAY3v3TZ43aeSmbbiSReurU9TvY1Y+GWy5GlV4V8f0KIOi8tL41/rv8nu87sYnqP6UxMPIfa8Dr0mGTcRy/nduW4Bds5eiaDzf+6CWeH6k/qlVl6VZrfhRBC1GkJmQlM/G0iexP28ma/N5mUlmkk9M73wa3vlZvQd0UmszUiiYcHtLRKQq8saX4XQghRZ51IOcHD6x4mLTeNuYPn0js2HH7/DwTdAbfPBZvy67Zz1h3H182BcSHNqyniayM1dSGEEHXSnoQ9PLDmAUxmE1/e8iW944/B6meg7XAYPR9sy6/XhkYlsyUikWn9a0ctHSSpCyGEqIPWRa1j6u9T8XHyYdHwRbQ/tR1WPgWth8A9XxavuFaeD9Ydx8fVgfE31I5aOkhSF0IIUYdorfny0JdM3zCdQJ9Avh7+NU2PrIVV/4Q2w+C+xWDvdMXzhEadZ/Nxo5bu4lB77lTXnkiFEEKIcuSb83l9++v8ePxHhjYfyut9X8dpzyL49Vloewvc+zXYOV75RMB7fxzFu5bV0sHCNXWl1C1KqaNKqQil1PNl7G+mlFqvlNqrlDqglBphyXhqow0bNjBy5Mgy902ZMoXw8HAA3njjjeoMSwghapS0vDQeWfcIPx7/kamdpvLOgHdwCv3aSOjtRlQqoW+NSGRrRBKP3dQaV8faVfe1WFJXStkCc4HhQBAwVikVVKrYS8BSrXU34D7gE0vFUxOYTKYqPd+CBQsICjK+UknqQojrVXR6NON/HU9oQiiv9nmVJ7s/ic2OebDmOWh3q7FASwUTutaa2b8doYmnE+NCmlk48qpnyZp6LyBCa31Sa50HfAfcXqqMBjwKX3sCcRaMx+JeffVVAgMDGTJkCGPHjuXdd99l4MCBvPjiiwwYMIA5c+awcuVKQkJC6NatG4MHDyYhIQGAmTNn8sADD3DzzTfTpk0b5s+fX3zejIwM7r77bgIDAxk3btxFS7bu3r2b559/nuzsbLp27cq4ceMA+Oabb+jVqxddu3bloYceoqCgoPq/ECGEsLB9Z/cxbvU4krKTmDdkHne0vgP+/gR++zcEjqxwp7gia8POsD8mlX8OblvjZ48riyXbFfyA6BLvY4CQUmVmAr8rpZ4AXIHBZZ1IKTUNmAbQrFn5v5ze3vk2R5KPXF3ElxHoHci/e/273DK7d+/mxx9/LHPp1ZSUFDZu3AjA+fPn2b59O0opFixYwOzZs/nvf/8LwIEDB9i+fTuZmZl069aNW2+9FYC9e/cSFhZGkyZN6NOnD1u3bqVv377F137rrbf4+OOP2bdvH2DMEPf999+zdetW7O3tefTRR1m8eDEPPvhglX4vQghhTb+e/JUZW2fQyLURcwfNpYVHc9g4G9a/Du1Hwd1fgK19hc9XYNa8+/sxWtV3ZXR3PwtGbjmWTOplTdFTek7ascCXWuv/KqVuBBYppTpqrc0XHaT1PGAeGNPEWiTaa7RlyxZuv/12nJ2dAYrXSQcYM2ZM8euYmBjGjBlDfHw8eXl5BAQEFO8rOt7Z2ZmbbrqJnTt34uXlRa9evWja1FhEoGvXrkRGRl6U1Ev7888/CQ0NLZ5fPjs7mwYNGlTp5xVCCGspMBfw8b6PWXBwAd0bdGfOTXPwcvSEP2bAto+gy1gY9fEVx6GX9tOeGCLOZvDJuO41eiW28lgyqccA/iXeN+XS5vXJwC0AWuu/lVJOgC9w9moveqUataWUN4e+q6tr8esnnniC6dOnM2rUKDZs2MDMmTOL95VeKrXofemlVa90b15rzYQJE3jzzTcr8xGEEKLGS89L5/nNz7MpZhN3tbmLF0NexEHZGkPWQr+EnlNh+OwrzhRXWq6pgA/WHaeTnyfDa/B66VdiyZ8iu4A2SqkApZQDRke4FaXKnAYGASil2gNOwDkLxmQxffv2ZeXKleTk5JCRkcHq1avLLJeamoqfn9Gs89VXX12075dffiEnJ4ekpCQ2bNhQqZXc7O3tyc/PB2DQoEH88MMPnD1r/DZKTk4mKirqaj6WEELUGKdST3H/6vvZFruNl0Je4pUbX8EBBT9NMxJ63+kw4p1KJ3SAJTtOE5uSzb9uaXdJBas2sVhNXWttUko9DqwFbIEvtNZhSqlZwG6t9QrgGWC+UuppjKb5ibq2LRtXqGfPnowaNYouXbrQvHlzgoOD8fT0vKTczJkzueeee/Dz8+OGG27g1KlTxft69erFrbfeyunTp5kxYwZNmjTh2LFjFbr+tGnT6Ny5M927d2fx4sW89tprDB06FLPZjL29PXPnzqV589o13lIIIYpsitnEvzf9G3sbe+YPnU9wo2DIz4FlE+HYGhj0CvSbflXnTsvJ5+O/IrixpQ99W/tWbeDVTJZerUIZGRm4ubmRlZVF//79mTdvHt27d6/QsTNnzsTNzY1nn33WwlFeqqZ8f0IIUZrWmi8OfcGcPXNo592OOTfNoYlbE8hNh+/GwamNMOJd6DX1qq/x1pojfLbxBCsf70unppdWxqytMkuv1q5R9TXctGnTCA8PJycnhwkTJlQ4oQshhLhUVn4WM7fNZE3kGm5pcQuz+szC2c4ZMs7C4nvgzEG44zPoOvaqrxGdnMUXW04xurtfjUzolSVJvQotWbLkqo8t2WFOCCGud5GpkTy94WlOpJzgqe5PMbnjZONed/JJWDQa0s/A2G+h7bBrus7stUexsYHnhrWrositS5K6EEKIGuWPqD+YsXUG9jb2fDbkM3o36W3siNtr1NDNBTBhJfhXvDNxWUKjzrNyfxxP3tyaxp7OVRC59dWZpK61rtU9Fq2ltvWpEELUXfnmfOaEzuGr8K/o7NuZ/w78L41cC4eXRfwJ3z8ALj7wwE/g2+aarqW15rXV4dR3d+ShAa2qIPqaoU4kdScnJ5KSkvDx8ZHEXglaa5KSknByuvIyhEIIYUnnss7x7MZn2XN2D2MDx/Jc8HPYF80Gt/97+OVRqN8exi0Dj8bXfL3VB+PZezqFt+/qVOsWbSlPnfgkTZs2JSYmhnPnauUQd6tycnIqnq1OCCGsYdeZXTy38TmyTFm81e8tbm1pTJGN1rDtQ/jjZWjRz1gL3enaO7Pl5Bfw1pojBDZy5+4e/lc+oBapE0nd3t7+oulWhRBC1HxmbeaLQ1/w8d6P8Xf3Z8HQBbSu19rYWZBvLJsa+iV0uBPu/LzCK61dybxNJ4k5n803k0Owtalbrbt1IqkLIYSoXRKzE3lh8wtsj9/OsBbDmHnjTNwc3Iyd2SmwbAKc3AD9noGbXrqqWeLKEp2cxdz1EYzo1Ii+bWr3RDNlkaQuhBCiWm2L3cYLW14wxqHfOJPRbUZf6A91PhKWjIGkCLh9LnQbX6XXfm11ODZK8dKtQVV63ppCkroQQohqkW/O5+O9H/PFoS9o7dWa/w3934XmdoDoXfDtfWDOhweWQ0D/Kr3+xmPnWBuWwHPD2tHEq24MYStNkroQQgiLi0mP4d+b/s2BxAPc0/Ye/tXzXzjZlRh5c+gn+PkRcG9s9HC/xiFrpeWaCpi5IowAX1em9Ku7fbAkqQshhLCo3yJ/Y9a2WQC8O+BdhrUoMQuc2QybZsOGN8H/BrhvCbj6VHkM/9tyilOJmXw5qSeOdrZVfv6aQpK6EEIIi0jLS+PNHW+y6uQqOvt25u3+b9PUvcQQ2tx0WP4wHFkFXcbCyA/AvurnzYhLyeajPyMYGtSQge0aVPn5axJJ6kIIIarcrjO7eHHLi5zLOsejXR9laqep2NmUSDnJJ+Hb+yHxGAx7E254BCwweZjWmldWhKHRzBhZNzvHlSRJXQghRJXJK8jjwz0f8nX41zTzaMai4YvoVL/TxYVO/AXLJhlJ/IGfoOVAi8WzNuwMf4Qn8MLwQPy9XSx2nZpCkroQQogqcTT5KC9seYHj548zpt0YpveYjot9iUSqNWz/BH5/CeoHGvfPvS3XaS01O5+XfwkjqLEHk/vW3c5xJUlSF0IIcU0KzAUsCl/Eh3s/xMPBg7mD5tK/aanhaHlZsOppOPAdtB8Fd3wKjm4WjWv2b0dIzMhlwYRg7GyrZvKamk6SuhBCiKt2MvUkL299mf3n9nOz/8280vsVvJ28Ly6UdAKWPggJYXDTf6Dfs1U2Q9zl7IpMZvGO00zuG0Dnpl4WvVZNIkldCCFEpRWYC/g6/Gs+3vsxTnZOvNH3DUa2HHnpSpnhK+DnR8HWHsb/AK0HWzy2XFMBL/x0ED8vZ6YPaWvx69UkktSFEEJUyomUE8zYOoODiQe52f9mZtw4A1/nUvOoF+TDupnw98fgFwz3fAle1bMi2qcbThBxNoOFE3vWqWVVK+L6+rRCCCGumsls4suwL/lk3ye42rsyu/9sbmlxy6W187Q4o3d79HboNQ2Gvg52DtUSY1hcKh//FcHtXZtwU2DdHpNeFknqQgghruho8lFe3vYy4UnhDG0+lBdDXsTHuYyZ305uhB8nGx3j7vofdLq72mLMM5l5Zul+6rk68H+jOlTbdWsSSepCCCEuK9uUzaf7P+XrsK/xdPS8dJrXIgX5sOEt2Pxf8G0LE1dD/XbVGuuHfx7nyJl0/jchGC+X6mkZqGkqlNSVUt5a62RLByOEEKLm2BK7hde2v0ZsRiyj24zm6e5P4+VURk/y81Hw4xSI2QndHoDhb4ODa7XGuj86hU83nuDuHk0Z1L5htV67JqloTX2HUmofsBBYo7XWFoxJCCGEFSVmJ/L2zrf5LfI3AjwDWDhsIcGNgssuHLYcVjwFaLj7C+h4V7XGCpCTX8Azy/bTwN2Rl2+r+1PBlqeiAwXbAvOAB4AIpdQbSqkrjhNQSt2ilDqqlIpQSj1/mTL3KqXClVJhSqklFQ9dCCFEVTJrM0uPLmXU8lH8efpPHuv6GD/c9kPZCT0vE1Y8AcsmQv228PBmqyR0gPf+OEbE2QzevqszHk72VomhpqhQTb2wZv4H8IdS6ibgG+BRpdR+4Hmt9d+lj1FK2QJzgSFADLBLKbVCax1eokwb4AWgj9b6vFLq+uuqKIQQNcDR5KO8tv019p3bR69GvZhxwwxaeLYou/CZg/DDPyDxOPR7Bga+YIxDt4JtJxKZv/kk94c0o3/b+laJoSap6D11H2A8Rk09AXgCWAF0BZYBZU2q2wuI0FqfLDzHd8DtQHiJMlOBuVrr8wBa67NX9zGEEEJcjdTcVD7e+zFLjy3F08GT1/q8xqhWoy4dpgZQYIJtc2D9m+DiAw/+Ai0HVH/Qhc5n5jH9+/0E+Ljy0q3trRZHTVLRe+p/A4uAO7TWMSW271ZKfXaZY/yA6BLvY4CQUmXaAiiltgK2wEyt9W8VjEkIIcRVMmszy48vZ86eOaTmpTKm3Rge6/oYno6eZR+QdMJY+zxmJwTdASPfBxfvsstWA601z/90gKTMXBZM6IOLgwzmgoon9Ze01ktLblBK3aO1Xqa1fvsyx5S1MG7pDnZ2QBtgINAU2KyU6qi1Til1rWnANIBmzZpVMGQhhBBlOXjuIG/seINDSYfo3qA7L4a8SDvvyww/0xp2f2GsrGZrb4w973iXRdY+r4xvd0azNiyB/4xoT0e/y/wQuQ5VNKk/Dywtte0FjKb3y4kBSs4J2BSIK6PMdq11PnBKKXUUI8nvKllIaz0Po6MewcHB0vNeCCGuQlJ2EnP2zGF5xHLqO9fnrX5vMSJgRNlN7WDMDPfL43DiT2h1M9w+Fzwyqo/YAAAgAElEQVSaVG/QZYg4m86sVWH0a+N73SypWlHlJnWl1HBgBOCnlPqwxC4PwHSFc+8C2iilAoBY4D7g/lJlfgbGAl8qpXwxmuNPVjx8IYQQV5JbkMviw4uZf2A+OaYcJnWYxENdHsLV/jJjybWGg8vg1+fAlAsj3oWeU6xeOwdjsZYnvt2Hi4Md/72nCzY21o+pJrlSTT0O2A2MAkJLbE8Hni7vQK21SSn1OLAW4375F1rrMKXULGC31npF4b6hSqlwoAB4TmuddHUfRQghRElaa36L/I0PQj8gLjOOAU0HMD14Oi09W17+oNQYWDUdjq+Fpj3hjs/At3X1BX0Fr64K53B8Gv+bEEwDDydrh1PjqIrMI6OUstNaX6lmXi2Cg4P17t27rR2GEELUaPvO7uOd3e9w4NwB2tVrx7M9n+WGxjdc/gCzGfZ8Cb+/DGYTDJoBIQ+DjW21xXwlv+yL5anv9vFQ/5a8MOL66e2ulArVWl9m9p+LXan5fanW+l5gr1KqZPZXGMPXO19DnEIIIapYTHoMH+z5gLWRa6nvXJ9ZvWcxqtUobMtLzkknYMWTELUFAvrDbR+Cd826V308IZ0XfjpIrxbePDuseueUr02u1Pz+VOHzSEsHIoQQ4uql5KSw4OAClhxZgp2NHY90eYSJHSbiYu9y+YPMBfD3XFj/Otg6GMm8+4M14t55SZm5Jh5ZvAcXB1s+ur8b9rYVnQz1+lNuUtdaxxe+TASytdbmwulhA4E1lg5OCCFE+TLzM/k6/Gu+CvuKbFM2o1qN4vGuj9PQ9QqLmsSGwqqnIX4/tBsBt/63RvRsL01rzQs/HeTkuQy+mRxCQ7mPXq6KDmnbBPRTStUD/sToPDcGGGepwIQQQlxeXkEeS48uZf7B+STnJDOo2SCe6PYErbxalX9gdgr8OcsYe+7WEO5eCB3urHG18yKLtkexYn8czw5tS+/WvtYOp8araFJXWusspdRk4COt9Wyl1F5LBiaEEOJSJrOJlSdW8un+T4nPjCekcQhPdXuKTvU7lX+g1nBgKfz+H8hKMjrB3fQiOHlUT+BX4e8TScxaGc6gwAY8OrDm9MCvySqc1JVSN2LUzCdX8lghhBDXyKzNrItax9x9czmZepKOPh35v97/x41NbrzyweeOwerpELkZ/IJh/I/QuIvlg74G0clZPLo4lOY+Lrx/X1cZj15BFU3MT2HMILe8cKx5S2C95cISQggBRjL/Pep3Pt//OREpEbT0bMn7A99nULNBl58JrkhuBmz+L2z7CBxcjPnau08Em5rd0Swz18TUr3dTYNYsmNDzul9OtTIquvTqJoz76kXvTwJPWiooIYS43hWYC/gj6g8+2/8ZJ1JP0NKzJbP7z2Zo86HlD08DY8z5ge9h3UzIOANdxsKQV8Gt5i9NajZrnlm6n2MJ6Xw5qRcBvpeZ9U6UqaJLr7YFngValDxGa32zZcISQojrU4G5gLWRa/n8wOecTD1JK89WvNP/HYY0H3LlZA4QsxvW/Mvo3e7XA8Z8A/49LR94Ffnwr+P8FnaGl25tL+ujX4WKNr8vAz4DFmBM5yqEEKIKmcwm1pxaw/yD8zmVeorWXq15Z8A7DG0+FBtVgebytHijZn7gO3BrZEzv2nlMjW9qL+mXfbF8sO44o7v7yUItV6miSd2ktf7UopEIIcR1KNuUzfLjy/kq7CviMuNo7dWadwe8y5DmQyqWzPOyYPtc2Pw+mPOh73ToNx0c3S0ffBXacTKJ55YdoFeAN2+O7nTl/gKiTBVN6iuVUo8Cy4Hcoo1a62SLRCWEEHVcam4q3x75liWHl3A+9zxd63flhZAX6N+0f8WSubkA9i2G9W9AejwEjoShr9W46V0rIuJsOlO/3o2/tzPzHuiBo13NmW++tqloUp9Q+PxciW0aKGepHyGEEKWdyTzDovBFLDu2jGxTNv2b9mdyx8l0b9i9YifQGo6thXWvwLkjxhC1u7+A5r0tG7iFnEvPZeLCXTjY2fDlpF54uThYO6RaraK932vfTz8hhKhBjiYfZVH4IlafWo3WmuEBw5nUcRJt67Wt+ElidsMfL0PUVvBuBfd+De1H1djZ4K4kK8/E5K92kZSRx/cP3YC/dznz1IsKqWjvdxdgOtBMaz1NKdUGaKe1XmXR6IQQohYzazObYzazKHwRO87swMnWiXvb3suDHR7Ez82v4idKPA5/vQrhv4BrfWOe9u4TwLb2jt/OM5l5bPEeDsWmMu+BYDo39bJ2SHVCRZvfFwKhQFH7TgxGj3hJ6kIIUUpWfhY/R/zMkiNLiEqLoqFLQ/7Z/Z/c3fZuPB09K36i5JOwcbYx5tzOGQY8D70fr3Wd4EorMGueWbaf9UfP8eboTgwOusLiM6LCKprUW2mtxyilxgJorbOVdE0UQoiLxGXE8e2Rb/nx2I+k56fTuX5n3un/DoOaD8LephK16pTTRjLft8Sojd/wKPT5Z62YPOZKtNa8/MshVu6P4/nhgYzt1czaIdUpFU3qeUopZ4zOcSilWlGiF7wQQlyvzNrMjvgdLD26lL+i/0KhGNJ8COODxtOlfiXnV0+Nhc3vwp5Fxn3ynlOM4WnujSwTvBW8s/Yoi3ec5uEBrXh4wBVWlBOVVtGkPhP4DfBXSi0G+gCTLBWUEELUdCk5Kfxy4heWHVtGVFoU9RzrMbHDRMYGjqWRayWTcGosbPsQdi8EbYbuD0C/Z8GzEvfda4F5m07wyYYTjO3VjH/f0s7a4dRJFe39/rtSKhS4AVDAU1rrRItGJoQQNYzWmgOJB1h6dCm/nfqNPHMeXet35eF+DzOk+RAcbR0rd8KkE7D1A9j3rZHMu94P/Z+Des0t8wGs6Mutp3jj1yPc2rkxr93RUSaXsZCK9n7/U2s9CFhdxjYhhKjTMvMz+fXUryw7uozDyYdxsXPhzjZ3ck/be2jnfRU1zoQw2PwehP0ENvbQYwL0frJOJnMwEvrMleEMDWrI+/d2xVaWUbWYcpO6UsoJcAF8lVL1MGrpAB5AEwvHJoQQVqO1Zu/ZvSyPWM7ayLVkm7JpW68tM26Ywa0tb8XV/ipWD4vZbSyFevRXcHCDGx83Hu51t/f3V9simbkynCFBDfn4/u442NWeuehroyvV1B8C/omRwEO5kNTTgLkWjEsIIawiITOBlSdX8nPEz0SlReFi58KIgBHc0foOutTvUvlmY7MZTvxp3DM/tQmc68HAF6HXVHDxtsyHqCG+/juSV1aEMSSoIXMloVeLcpO61noOMEcp9YTW+qNqikkIIapVfkE+G2I2sPz4crbGbcWszfRo2IOpnaYypPkQXOyvYqaz/GxjfPnfn0DiUXBvYqxpHvwPcHSr+g9RwxQ1uQ9uLwm9OlW0o9xHSqneXLqe+tcWiksIISyqqNPb6pOr+e3Ub5zPPU8D5wZM7jiZ21vfTnOPq7y/nXEOdi0wHlmJ0KgzjJ4PHe6s1TPAVZTWmrnrI3j392NSQ7eCinaUWwS0AvZxYT11DZSb1JVStwBzAFtggdb6rcuUuxtjhrqeWuvdFQtdCCEqLzI1ktWnVrP65Gqi06NxsHFggP8A7mx9J72b9MbW5ipXCDt7BP7+GA4shYJcaDvcmP2teZ9aOzd7ZWmteXPNEeZtOsmd3fyYfXdn7G0loVenio5TDwaCtNa6oidWStli3HcfgjGt7C6l1AqtdXipcu7Ak8COip5bCCEqIzE7kTWn1rD65GrCksJQKHo17sXUTlMZ3Hww7g5XOe1qQb7R6W3nfIjcbEzl2m28MQOcb+uq/RA1XIFZ85/lB/luVzQP3ticmbd1wEZ6uVe7iib1Q0AjIL4S5+4FRGitTwIopb4DbgfCS5V7FZgNPFuJcwshRLlSc1NZH72eNafWsD1+O2Ztpr13e54NfpZbWtxCQ9dr6HGefgZCv4LQhcZa5p7NYNArxiIrrj5V9yFqiTyTmaeX7mP1gXieuLk104e0lXHoVlLRpO4LhCuldlJielit9ahyjvEDoku8jwFCShZQSnUD/LXWq5RSktSFENckNTeVv07/xe9Rv7M9fjsmswk/Nz8md5zMyJYjaenV8upPrjVEbYNd8+HwSjCboPVgGPk+tBkKV9tsX8ulZufzyDehbDuRxIsjApnWX6Z+tabKTBNbWWX9TCtuvldK2QDvAxOveCKlpgHTAJo1k8n/hRAXnM85X5zId8bvxKSNRP5A+wcY2mIoHXw6XFutMSsZDi6D0C/hbDg4eUHIw0Yvdp/rO4HFpmQzaeFOTiVm8t69XRjdvam1Q7ruVbT3+8arOHcM4F/ifVMgrsR7d6AjsKHwf7hGwAql1KjSneW01vOAeQDBwcEVvq8vhKibzmWdY0PMBn6P/J1dZ3ZRoAto6taUBzs8yNAWQwnyDrq2RG42w6mNsHcRHF5ldHxr3BVGfQwd7wKHqxjiVseExaUyaeEusvMK+GpSL3q39rV2SIIrzyi3RWvdVymVTolaNkYtXGutPco5fBfQRikVAMQC9wH3F+3UWqdiNOsXXWsD8Kz0fhdClKa15kTKCdZHr2dD9AYOJB4AoJl7MyZ1nMTQ5kMJ9A689vu4qTGwdzHs+8ZY/tTJC3pMNBZYadTp2j9IHbHh6FkeW7wHT2d7fnikN+0a1e713euSK00+07fwudL/xbTWJqXU48BajCFtX2itw5RSs4DdWusVVxOwEOL6YDKb2Ht2L+uj17P+9HpiMmIA6OjTkce7Ps5A/4G0rVcFHbLyc+DYGmO50xN/ARoCBhgd3wJHgr3TtX+YOkJrzcKtkbz+62HaNXRn4aSeNPSQ76cmUZUYpVYjBAcH6927pTIvRF2UmpvK3/F/szF6I5tiNpGWl4a9jT0hjUO4yf8mBjQdcG291ouYzRC11ZjxLfwXyE0DDz/oOg66jYN6La79GnVMrqmA/yw/xA+hMQwNash7Y7ri5ljRblniWiilQrXWwRUpK/9FhBBWY9ZmDicdZkvsFrbEbuFA4gHM2oynoycD/Qcy0H8gfZr0ubppWsuSEG4k8oPLIC3WWFSl/SjofI9RO79Oe7Bfydm0HB76JpS9p1N4clAb/jmojYxBr6EkqQshqlVyTjLb4raxNXYr2+K2kZyTDEAHnw5M6TSFfn796OjbETubKvrzlBoLh340ZnpLOAjK1hiKNmQWtBshnd6uYF90Cg8t2k16jolPx3VneKfG1g5JlEOSuhDCovIK8th/bj874newNXYrYUlhaDT1HOvR2683fZr0oXeT3vg4V+GkLamxRrN6+M8QXThZpV8wDH/HmIPdrX7VXauO0lrzzfYoXl11mAYejvz4SG/aNy6vb7SoCSSpCyGqVIG5gMPJh9kev52d8TvZe3YvOQU52CgbOvt25tGuj9LXry9BPkHYqCqcF7ysRN6wE9z8EnQYfd2PKa+M9Jx8XvjpIKsOxDOgbX3eH9MVb1cHa4clKkCSuhDimhQNN9txZgc74new+8xu0vPTAWjt1Zq72t5FSKMQghsFX/0c65eTGgPhK8pO5EF3Xnfzr1eF8Lg0Hluyh6ikTJ4b1o5HBrSS++e1iCR1IUSlFJgLOJ5ynNCEUPYk7CE0IZSknCQA/Nz8GNpiKL0a9aJX4174OlfxhCRaw5kDcHQNHFltvAZjDPnNMyDoDknkV0lrzfe7onllRRiezvZ8O/UGQlpef/PY13aS1IUQ5coryONQ4iH2nDUS+P6z+4tr4g1dGhLSOISQxiH0atSLpu4WmCbUlAdRW+DIr0YyT4sBFPiHGJ3dAkdK0/o1Ss7M44WfDrA2LIG+rX354L6u+Lo5WjsscRUkqQshLpKam8rBxIPFtfBDiYfIM+cBEOAZwLCAYXRv0J0eDXvQxK2JZYLITDImgjn6K0SsM8aR2zlDq5vhphegzTDp7FZFNhw9y3M/HCA1K58XRwQypW9LaW6vxSSpC3EdM5lNHD9/nIOJB9l/bj8Hzh0gMi0SAFtlS6B3IGMCx9CjQQ+6NeyGt5O3ZQIxF0DsHiOBR6yD2FBAg2sD6HCHMfSs5UCwd7bM9a9D2XkFvLXmMF/9HUXbhm58NakXQU2kd3ttJ0ldiOtIQmYCBxMPcuDcAfaf28/h5MNkm7IB8HbyprNvZ25rdRudfDvRpX6Xqpv0pSzpCXDiTyOJn/gLss8DCpoGw8DnjbHkTbqDTRX2kBcA7Dl9nueW7efEuUwm9w3guWHtcLKXiXfqAknqQtRRidmJhCeFczjpMOFJ4YQlhZGQlQCAvY097b3bc1ebu+jk24nO9Tvj5+Z37fOolycvC6K3w8mNRhIv6uTm2gDaDofWg4zmdRcLtQYIsvJMvLv2GAu3naKxhxPfTA6hbxtZXa0ukaQuRC2ntSYhK8FI3slGEj+cdJiz2WeLy7TwaEH3ht3p7NuZzvU7E+gdiIOthccdm/KMZvRTm4xHzE4oyAMbO2jaCwa9bNTGG3aS2ng12BqRyPM/HSA6OZsHb2zOv24JlLnb6yD5LypELZJvzicqNYrjKcc5dv4Yh5ONBF401aqNsiHAI4CQxiG092lPe+/2BHoH4ubgZvngzAVw5qCxDvmpTRD1N+RnAgoad4GQh4351ZvdAI7VEI8AIDUrnzfXHOa7XdEE+Lqy9KEb6RUgrSF1lSR1IWogrTXnss9x/LyRvIueT6aeJN+cD4CdsqOlV0v6N+1Pe+/2BPkE0bZeW8veBy8pPwfi9sDpv40EHr0TclONffUDjdXOAvpD8z7SpG4FZrPmh9AY3vrtCKnZ+Tw8oBX/HNxG7p3XcZLUhbCylJwUTqae5GTqSSJSIoqTeEpuSnGZBs4NaOPdht5NetOmXhva1mtLgGeA5ZvQS8pOMWZtK0ricXuM5nQwknjH0dC8t5HI3RtVX1ziEodiU3np50Psi06hZ4t6zLq9o8zbfp2QpC5ENTBrM2cyz3Ay9SSnUk8ZSTzFeH0+93xxOWc7Z9p4tWFQs0HFybuNVxu8nLyqOWAzJJ+AmN0Qs8tI5glhgDbuiTfpZjSnN7vRaE6XmniNkJqVz7u/H+WbHVH4uDrw3r1duLObhTtAihpFkroQVSgjL4PT6ac5nX6aqNQoTqWd4mTKSSLTIouHjgF4OnrS0rMlNze7mQDPAAI8A2jp2ZImbk2qdpGTispMgtjdRhKP3W10cMspbEp3cAf/nhB0u5HE/XrIcqU1TJ7JzKLtUXz013HSsvOZcGMLnh7SFk9ne2uHJqqZJHUhKik1N5Xo9GhOp50mKj2K6LRoTqefJjo9urjDWpEmrk0I8AygR8MetPRqSYBHAC29WlpuEpeKyMuChENG4i5K4ucjjX3KBhp0MJYn9Qs2xoz7tgUbuQ9bE2mtWXPoDG//doSopCz6tvblxRHtZRKZ65gkdSFKySvIIz4zntiMWOIzjOfYjFgjkaefJrWoM1ihRq6NaObejJv8b6K5R3OauTfD38Ofpm5Nq6/T2uVkpxg90s8cgPj9xiPxGGizsd+9iZG4g/9hJPEmXcHB1boxiwoJjTrP66vD2XM6hXYN3flyUk8GtK0vTe3XOUnq4rqTY8ohPjOeuIw44jLjiMuIIzYjlriMOOIz4i8a3w3GdKkNXRrSzKMZw5oPo5lHM/zd/Wnu0Rw/Nz+c7Jys9ElKyTgHZ/ZfSN7xB+D8qQv73ZsYQ8uC7oDGnY374h4WmrtdWExYXCofrDvOH+EJ1Hd35O27OnF3D39sZb52gSR1Ucdk5WeRkJVgPDITOJt19pL3RcuEFrFTdjR0bYifmx+9/XrTxLUJTdyMh5+bHw1cGmBnU4P+V8nLhLNH4GwYnD1sdGA7Gw6Z5y6UqRdgJO7uDxiJvFEXWQClljt6Jp0P1h1jzaEzeDjZMX1IWyb3DcBVJpARJci/BlErZJuyScpOIiknicTsRJKyk4oT9tmss8UJu2hJ0JI8HT1p4NKAhi4NCfIJopFrI/zc/IqTdn3n+tjWxHvGBfmQdKJE8g43Xp+PArRRxt7FGE7Wdhg0CIJGnY21xZ2rube8sJiIsxnM+fM4qw7E4epgx5OD2jC5b4B0ghNlkqQurCbblM35nPOXJOvE7ESScpIu2p6Zn3nJ8TbKBl8nXxq4NKCFZwtCGocYydu1IQ1djEd9l/o429Xwlb2yUyApwrjXnXgMEo8bj+STUDjRDMoWfFpD467QdZyRwBsGgVcLmWK1jjoYk8pnG0+w5lA8Tva2PDqwFVP7tcTLpRrnJhC1jiR1USWyTdmk5KRwPvc8Kbkpxa9Tc1M5n1O4rfBxPsfYnlOQU+a5PBw88HH2wdfZl/be7fF19sXH2QcfJ5/i7UWva1SzeHnMBZAaYyTrpOMlkvcxyEi4UM7GDrxbgW8bCBxh1MIbBBk90O1ryL17YTFaa7ZGJPHZxhNsiUjE3dGOhwa0YkrfAHzcHK0dnqgFaslfRGFpuQW5pOelk5aXRnpeOhl5GRe9L37kp1/8vrBMbkHuZc/t4eBBPad6eDl60dClIe3qtaOeUz08HT2p51ivOFH7Ovvi7eRdvbOkVaX8HEiJguRTRge15FNGbfv8KaPJvKjWDeDkBfXbQeshRgL3bWs86jUHW2lWvd6YCsysDUvgs40nOBibSn13R14YHsj9Ic1wd5J/D6LiJKnXIgXmAnILcosfOaac4tfZpmyy8rPINGWSlZ9V/D7LlEVmfiZZpqzi99n52Zdszy+ZcMpgp+zwcPTAzd4Ndwd33B3caeDSAA8HD9wd3PFy9DIeTl7Uc6xX/NrDwaP21KavRGujM1pKNKSeLpW8T0FaLMX3usGYtMW7BTTsAIEjwTvgQvJ28QEZenTdO5+Zx3e7ovlmexSxKdkE+Lry5uhO3NnNT+ZoF1fFon9tlVK3AHMAW2CB1vqtUvunA1MAE3AO+IfWOsqSMZWUmJ1IdHo0Zm3GXDhut+i1RqO1vuS1GTNaazSF+7SmQBdgMpvIN+djMpsuPHQZ2y5TLseUQ15BHjkFFxJ1rimXnIIL201mU6U/o5OtEy72LrjYuRQ/uzm40cClQfF7Z3tnIznbu+PmYCTtomRd9HCydar7418L8o3EnBpTmLijIeV04XO0sb10i4RrfaOneYu+RtKuF2A8e7eUxC0uKywula+2RfLLvjhyTWZubOnDjJHtGRLUSIamiWtisaSulLIF5gJDgBhgl1JqhdY6vESxvUCw1jpLKfUIMBsYY6mYStsUs4lXtr1i8evYKBvslB12Nhc/7G3ssbexx87GDkdbRxxtHXGxc6GeYz0c7RyLtznaOuJk54SDrQNOtheeS5ZxtnPGxd4FVzvXC8nazrlm9uq2hvwcSI+H9DOQHgdp8YXv4yE11kjc6fEXJmUp4toAvPyNHuWBI8DT33h4+UO9FuDobpWPI2qfnPwC1oad4ZvtUeyKPI+zvS139WjKhBtb0K6R/DsSVcOSNfVeQITW+iSAUuo74HagOKlrrdeXKL8dGG/BeC7Rp0kfPh/yOTbKBoUqflZKXbLNRtlctF0phQ0XthUl56JkXZy8lZ0kVksy5RlN4plnIT3BSNjpZyAt7kIST4uD7ORLj7VzAvfG4OFnrCxWlKw9m4JnM+NZOqeJaxQel8bS3dEs3xtLanY+/t7O/GdEe+4N9sfTRe6Xi6plyaTuB0SXeB8DhJRTfjKwpqwdSqlpwDSAZs2aVVV8xtAn14ZVdj5RBbSGvAwjUWcUJuvi14XvS77OSS3jJArcGhgJ26sZ+PcyZlPzaGwsCVr02slLmseFRaTl5LNiXxxLd0dzICYVB1sbhnVsxJhgf3q38sFGmtiFhVgyqZf1r1aXsQ2l1HggGBhQ1n6t9TxgHkBwcHCZ5xA1kNkMOSmQfR6yko3n7OTC16WfzxuPzEQosZrZRZy8jGTtWt/ofOZ2k/G66OHeyEjkbg3Bto50zhO1Rp7JzMZj5/hlXyzrDieQk28msJE7r9wWxB1d/ajnWktHdYhaxZJ/+WIA/xLvmwJxpQsppQYD/wEGaK0vPy5KVD9zAeSmQU6aUSMuel38nFrifYnXRUk8J+XSe9RFlI2RpF28wdnbSMYNOxidy1zrFybvBuDqa7x28QU7+aMoahazWbMzMplf9sXx68F4UrPz8XZ14J4e/tzdoymdm3rW/Q6mokaxZFLfBbRRSgUAscB9wP0lCyilugGfA7dorc9eegpRIVqDKcd45OcYc4PnZUB+1oXXeVkV316UnPMyrnxtOydw9AAnjwvPnv7gXO9Cwr7ouZ7xcPKSmdBErWQ2a/ZGp/DboXhWHYgnPjUHFwdbhgY15PZufvRt7Yu9rfzbFtZhsaSutTYppR4H1mIMaftCax2mlJoF7NZarwDeAdyAZYW/Zk9rrUdZKqZL5GUZyUubjcSozYU1S11i2+W2l9pmLjCGRJnzocBU+JwPZpPxKN6XX+p9qbL52WDKvZCkL3pd+D6/1L5yJn4pk7IBe1djiU0Hl8JnNyPRejQBR8+Lk7STZ6nE7Xlhm9SexXUgv8DMjpPJrA07w9qwM5xNz8XeVtGvTX2eHx7IkKCGuDjILR9hfRb9V6i1/hX4tdS2l0u8HmzJ61/RoR/+n737Do+jvBY//j3aXfVeXOXeuzG2MRiMqaETuKGFnhDzCyHJJSEJqZcLKYRLcgMJNwmQQAg9QIJNCSSAAQM2lsHGuIFcJRdZlqyulbac3x+zkleyysrWalXO53nmmZl33pk9u17r7Lwz876w9OsxDcEhTi9icR7nbmt3+JQQOhtOc5ql3QngTjpU7k4AT/h6opOgw5O1p2k5NLkT7QYxYzpR1+hnxWcHeHVDCf/eVEJlvY8kj4vFk/I4a/oQTpk8iHTr7c30MgP7p+WIBXDe/wLinL1KnJPsmpaby+VQeVt1EafPbpfbScwuD8S5wpY9rba5w8pDdY0xMaWqbC2tZfmW/bz1aSmrtpXTGAiSnujm9KmDOWvaEBZNzLOe3kyvNng4SnYAACAASURBVLCTet5EZzLGDEj1jQHe33aANzeXsvzT/RSVO09eTBiUyrUnjGLxpEHMH5Nt18hNnzGwk7oxZkBp8AdYu6uC97aW8f62Mj7adRBfQEnyuFg4PocbF41j8aQ88rOSYx2qMUfEkroxpt/yB4J8sqeK97Ye4P2tZazeUY7XF0QEpg/L4EsLx7BwfC7zx2Rbs7rpFyypG2P6jWqvj492VbBm50HW7DzIR7sOUtsYAGDS4DQunzeSE8blcNyYHOui1fRLltSNMX2SqlJ8sJ41Ow9SsLOcgh0H2VJSjSrECUweks7Fc/I5bmw2C8bmkJuaEOuQjYk6S+rGmD6hpMrLx8WVrC+u4OPdlawvrqSsthGA1AQ3x4zM5KzpQ5g7KptZIzJIs8fNzABkSd0Y06uoKiVVDWzcWxlK4pV8vLuS0mqnkyVXnDBhUCqnTh7EzBGZHDsyi0lD0mwccmOwpG6MiaHaBj+fllSzeV81W/ZVs2lvFVtKqqmo8wFOVxDj81I5aUIuM4dnMCM/k6lD00mKt5vajGmLJXVjTNTVNwbYdqCGraW1FO6vYcu+Kjbvq2ZXeR0aGncxJd7FxCFpnD19KFOGpjF5SDrThqWTkmB/poyJlP1vMcZ0C1WltLqBwlIneW/dX8PW0hq2ldayu+LQcLpxAqNzU5g2LJ3/mJPP5CFOAs/PSrJxxo05SpbUjTERCwSVPRX1FJXXsbO8jl3ldewqc+Y7DtRS3eBvrpsc72JcXirzRmdxed4IxualMm5QCqNzUuyZcGOixJK6MaZZMKiU1Tayt7I+lLzr2Vley67yenaVOWfcvoA21/e4hPysZEZmJ3PMyEzG5aU606AUhqQn2ljixvQwS+rGDCDVXh97K73srqhnb4WXPRX17Akl8L2VXvZWeGkMBFvsk5HkYWR2MtOGZ3D2jKGMynaS+MicZIZmJNld58b0IpbUjenjVJUqr5/Sai/7qxrYX91AaXUD+6u97K9uCJU528Kbx8F5PGxwWgLDMpOYmZ/JWdMTGZaRxNCMRIZlJjEiK9l6XjOmD7Gkbkwv5AsEOVjXSHltI+U1jZSHlstqnPmBGid5NyXrBn/wsGMkeuIYlJZIXloCEwenceL4XIZlJjE0M4nhmYkMzUhiUFoCbhuBzJh+w5K6MVGkqtT7AlTW+5ypzte8XFHnc5J1TSNltY3NSbyspoEqr7/dY2Yme8hNTWBQWgJzRmYxKC2BQWmJDEpPIC9sOS3Bbde0jRlgLKkb04FgUKnzBajx+qlpcKbaBj/VXj9VoeRc5T2UqMOnpu3hN5a15nEJWcnxZKfEk5Maz7Rh6eSkxJOdkkB2ajzZYduyU+LJTPLYmbUxpl2W1E2/oao0+IPUNwao8wWobwxNvgB1jX68vgB1jYHm5FzjDSXoUKJ2ygPUeH3UNjj1ahv9zZ2jtCdOID3JQ0bYNCwzqXk5PbHltuYp2UN6op1NG2O6jyV1ExWqij/oJNkGX4DGQJAGX9BZ9wdo9B9abvAFw7YHQuXO5A0l57rGQCgp+6n3hSfrQHOyrvcFOk3A4VxxQmqC+9CU6CYjyUN+ZhIpCS5SEzykJrpJDS2nJLhIS3STEu8mJcHdnJhT493WaYoxplewpB4jqkpQnc48AkEloM48GFoOtijj0LIq/oAzD4TXbV6mRZk/qPgCQfzBIL5AaLlpHlR8/iC+oOIPBPEFnDr+YBCfX/EFD9VtKvcHlMZAEH9o/0a/s1940m4MJetgFxJsexLccSTHu0jyuEiKd6Zkj5uM5HiGhpUleVwkx7tIDM2b63tcJMe7SYqPI8njJineRWqCm7RENwnuODtLNsb0KwM6qb/48R7uemVz89mdqqKAKigamhPaHr4eVi+0TFvbwvYhbD0YSui9iccluOPicLsEjyuued0TWne7wpbjhOR4d4u6iW4X8e44EtxxJHhcJLjjiHfFkeCJI8EdWncfWm4qb94ntN+hfZz9LekaY0zkBnRSz01NYP7obBAQBBEQCM1D6wK0tw2ak87h5c620O4ttsWJEBcnuERwxYErLg5XnFPuinOm5uWmumHb3WHbDx0nbJ9Wx4oPJWV3nBDvduZuV1yo3DmeJU9jjOn7BnRSXzA2hwVjc2IdhjHGGNMt7NkYY4wxpp+IalIXkbNEZIuIFIrIbW1sTxCRp0PbV4nI6GjGY4wxxvRnUUvqIuIC7gfOBqYCV4jI1FbVvgwcVNXxwP8Cv4xWPMYYY0x/F80z9flAoapuU9VG4CngwlZ1LgT+Elp+FjhN7I4tY4wx5ohEM6kPB4rC1otDZW3WUVU/UAnYnWvGGGPMEYhmUm/rjLv109mR1EFElohIgYgUlJaWdktwxhhjTH8TzUfaioERYev5wJ526hSLiBvIAMpbH0hVHwAeABCRUhHZ2Y1x5gIHuvF4pn32WfcM+5x7hn3OPcM+ZxgVacVoJvXVwAQRGQPsBi4HvtiqzlLgWuB94AvAG6od996tqnndGaSIFKjq3O48pmmbfdY9wz7nnmGfc8+wz7lropbUVdUvIjcDrwIu4M+qukFE7gAKVHUp8CfgryJSiHOGfnm04jHGGGP6u6j2KKeqLwMvtyr7SdiyF7gkmjEYY4wxA4X1KBe6Vm96hH3WPcM+555hn3PPsM+5C6STS9jGGGOM6SPsTN0YY4zpJwZ0Uu+sb3pzZERkhIi8KSKbRGSDiHwzVJ4tIv8Skc9C86xYx9ofiIhLRD4SkRdD62NCYyl8FhpbIT7WMfZ1IpIpIs+KyObQ9/p4+z53PxG5JfQ34xMReVJEEu373DUDNqlH2De9OTJ+4NuqOgVYAHwt9NneBryuqhOA10Pr5uh9E9gUtv5L4H9Dn/NBnDEWzNG5F/inqk4GZuF83vZ97kYiMhz4BjBXVafjPDV1OfZ97pIBm9SJrG96cwRUda+qfhharsb5Azicln39/wX4fGwi7D9EJB84F3gotC7AqThjKYB9zkdNRNKBRTiP4KKqjapagX2fo8ENJIU6I0sG9mLf5y4ZyEk9kr7pzVEKDad7DLAKGKyqe8FJ/MCg2EXWb/wG+C4QDK3nABWhsRTAvtfdYSxQCjwcuszxkIikYN/nbqWqu4F7gF04ybwSWIN9n7tkICf1iPqdN0dORFKB54D/VNWqWMfT34jIecB+VV0TXtxGVfteHx03MAf4vaoeA9RiTe3dLnRPwoXAGGAYkIJzebQ1+z53YCAn9Uj6pjdHSEQ8OAn9cVV9PlRcIiJDQ9uHAvtjFV8/sRC4QER24Fw+OhXnzD0z1HwJ9r3uDsVAsaquCq0/i5Pk7fvcvU4Htqtqqar6gOeBE7Dvc5cM5KTe3Dd96G7Ky3H6ojdHKXRd90/AJlX9ddimpr7+Cc1f6OnY+hNV/b6q5qvqaJzv7xuqeiXwJs5YCmCf81FT1X1AkYhMChWdBmzEvs/dbRewQESSQ39Dmj5n+z53wYDufEZEzsE5s2nqm/5nMQ6pXxCRE4F3gPUcutb7A5zr6s8AI3H+A1+iqoeNyme6TkQWA7eq6nkiMhbnzD0b+Ai4SlUbYhlfXycis3FuRowHtgHX45wU2fe5G4nIfwOX4TxB8xFwA841dPs+R2hAJ3VjjDGmPxnIze/GGGNMv2JJ3RhjjOknLKkbY4wx/YQldWOMMaafsKRujDHG9BOW1I0ZgEKjjt0UWl7cNMJbF/a/TkSGRSc6Y8yRsqRuzMCUCdx0FPtfh9OVpzGmF7Hn1I0ZgESkaVTCLYAPpz/zA8B0nEE0rlJVFZFjgV8DqaHt1+F0T/sIsBuoB44HvgOcDyQB7wE3qv1xMabHWVI3ZgAKjZ73oqpOD/VG9wIwDadf7XdxkvQq4C3gQlUtFZHLgM+p6pdEZDlOD3YFoeNlN/WmJiJ/BZ5R1WU9+66MMe7OqxhjBoAPVLUYQETWAqOBCpwz9385XXHjwhkSsy2niMh3ccbAzgY2AJbUjelhltSNMQDhfWkHcP42CLBBVY/vaEcRSQT+D5irqkUicjuQGK1AjTHtsxvljBmYqoG0TupsAfJE5HhwhtMVkWlt7N+UwA+ISCqHRtQyxvQwO1M3ZgBS1TIReVdEPsG52a2kjTqNIvIF4D4RycD5e/EbnKb1R4A/iEjTjXIP4ozKtwNnWGNjTAzYjXLGGGNMP2HN78YYY0w/YUndGGOM6ScsqRtjjDH9hCV1Y4wxpp+wpG6MMcb0E5bUjTHGmH7CkroxxhjTT1hSN8YYY/oJS+rGGGNMP2FJ3RhjjOknLKkbY4wx/YQldWOMMaafsKRujDHG9BOW1I0xxph+wpK6McYY009YUjfGGGP6CUvqxhhjTD9hSd0YY4zpJyypG2OMMf2EJXVjjDGmn7CkbowxxvQTltSNMcaYfsKSujHGGNNPWFI3xhhj+glL6sYYY0w/YUndGGOM6ScsqRtjjDH9hCV1Y8yAJiKviciVsY7DmO5gSd2YLhCRmrApKCL1YesDPjGISLGIeEUkq1X5ehFREcmPVWyhOH4qIo+El6nqmar6eIxCMqZbWVI3pgtUNbVpAnYB54eVHZYYRMTd81F2TRRi3AlcHnb8OUB8N7+GMaYNltSN6UahM8GnReRJEakGrhKR40VkpYhUiMheEblPRDyh+u7QGeyNIlIoIgdF5L6w400UkbdFpFJEDojIE6Hyh0Tkrlav/ZKIfCO0nC8ifxeRUhHZLiJf6yTGBSLyoYhUiUiJiPxPWP2FYfGvFZFFnXwMfwWuCVu/Bni0VawXhI5VLSK7ROTHrbZfHyo/ICI/CLUALA6L/0kReSy0/yehHw5N+7b53kXkPOC7wJWhlpU1ofIVInJd2P43isjmsGPP6uT9GtN7qKpNNtl0BBOwAzi9VdlPgUbgfJwfzUnAPOA4wA2MBT4Fbg7VdwMKvABkAKOB8qbjAn8Dvhc6ViKwMFR+auj1JbSeA9QDgwEXsBb4Ac4Z8vhQ3dM6iHE1cEVoexpwXGh5BFAGfC5U9yzgAJDTzmdSDCwGCoEJofe3O/S+FMgPi3966JizQsc8L7RtBlANnAAkAP8L+IHFYfHXh2JyAf8DrAhti+S9P9Iq5hXAdaHlK4Ai4FhAgInAiFh/12yyKdLJztSN6X4rVHWZqgZVtV5VV6vqKlX1q+o24AHg5Fb7/EJVK1V1B7AcmB0q9+EkxKGq6lXVd0PlywEPcHxo/VLgHVUtARYA6ar6c1VtVNVC4E+ENYm3jjH0OhNEJEdVq1V1VajeNcBSVX01VPefwDqc5N6Rx0L7ngV8DOwL36iqb6jqJ6FjrgOeCvtMLgH+oarvqWoD8KM2jv9WKKYATstA0+cVyXvvyA3AXaq6Rh2fqmpRhPsaE3OW1I3pfi2SgIhMDjWN7xORKuAOILfVPuFJrw5IDS1/Gyd5F4RuNrsWQFWDwNM4Z5YAXwSarumPAkaGmssrRKQCp9l5SHsxAtcDU4EtIvKBiJwTdqwrWh1rATCsk8/gUeBK4FpaNb2HPpPjRWR5qIm8EieZNn0mw8LjU9Va4GCrQ7T+vFK68N47MgLYGmFdY3qdXn8TjzF9kLZa/yOwErhMVWtE5FbgvIgOpLoXJ+ERupb9LxF5W1W3A08CL4rIr4E5wN9DuxUBn6nqlEhjVNUtwOUiEodzpvxc6A72IuBhVf1qJPGGHW+biOwBzqDl9fUmTwH3AGepqldEfsehHzJ7cZIzofedAmQdfog2dfbeW//btLX/uAhfy5hex87UjYm+NKASqBWRKcCNke4oIpeKyPDQagVOUgoAqOrq0HEfAF5W1apQvfeBRhH5togkiohLRGaIyLEdvM7VIpIbagGoDL1OEKdp+yIROSN0nEQROUVEOjtTB7gO51p2fRvb0oDyUEJfQMvm8b8Bnw/dvBeP07IRqc7eewkwWkSknf0fAr4rIseIY4KIjOjC6xsTU5bUjYm+b+M0Q1fjnLU/3YV9jwNWi0gt8DzwNVXdFbb9SeB04ImmAlX1A+cA83FuEjsQet30Dl7nHGBT6G74e3BaFRpD1/gvAn4MlOI8xvdtIvjboaqFqrqmnc1fBX4Rer0fAM+E7fcxcAtOct+Dc6NeGdAQwWt29t6fxrmBrlxEPmhj/yeBX4bqVeF85pG2EhgTc013zhpjTK8kIuk4rRSj7KY1YzpmZ+rGmF4n9Bx7soikAr8CPrSEbkznLKkbY3qji3Ca3otxHum7osPaxhjAmt+NMcaYfsPO1I0xxph+wpK6McYY00/0uc5ncnNzdfTo0bEOwxhjjOkRa9asOaCqeZHU7XNJffTo0RQUFMQ6DGOMMaZHiMjOSOta87sxxhjTT1hSN8YYY/qJqCV1EfmziOwXkU/a2S4icp+IFIrIxyIyJ1qxGGOMMQNBNK+pPwL8jjaGXQw5G5gQmo4Dfh+ad5nP56O4uBiv13skuw9oiYmJ5Ofn4/F4Yh2KMcaYoxS1pK6qb4vI6A6qXAg8qk7vNytFJFNEhoaGmuyS4uJi0tLSGD16NO0PvmRaU1XKysooLi5mzJgxsQ7HGGP6vIA/gMvtitnrx/Ka+nCcsYubFIfKuszr9ZKTk2MJvYtEhJycHGvhMMaYbvKDRy/kcw9Nx+/3x+T1Y5nU28rAbfZZKyJLRKRARApKS0vbPpgl9CNin5sxxnSf/b79COB2x+aJ8Vgm9WJgRNh6Ps4ADodR1QdUda6qzs3Li+j5+17juuuu49lnn411GMYYY6JMVSmTOvI0KWYxxDKpLwWuCd0FvwCoPJLr6cYYY0xvsKe8mhKPMsSTG7MYotY+ICJPAouBXBEpBv4L8ACo6h+Al4FzgEKgDrg+WrH0lDvvvJPHH3+cESNGkJuby7HHHtti+x133MGyZcuor6/nhBNO4I9//CMiwn333ccf/vAH3G43U6dO5amnnuKtt97im9/8JuA0kb/99tukpaXF4m0ZY4yJwNrNK6mLi2NkRuxuPI7m3e8djn8cuuv9a939uv+9bAMb91R16zGnDkvnv86f1mGdgoICnnvuOT766CP8fj9z5sw5LKnffPPN/OQnPwHg6quv5sUXX+T888/nrrvuYvv27SQkJFBRUQHAPffcw/3338/ChQupqakhMTGxW9+TMcaY7vXZzlUATBsxO2YxWI9y3WTFihVceOGFJCUlkZaWxvnnn39YnTfffJPjjjuOGTNm8MYbb7BhwwYAZs6cyZVXXsljjz3WfHPFwoUL+da3vsV9991HRUVFzG66MMYYE5mSik0AjMtfELMY+l2m6OyMOlqchof2eb1ebrrpJgoKChgxYgS3335786NkL730Em+//TZLly7lzjvvZMOGDdx2222ce+65vPzyyyxYsIB///vfTJ48uSfeijHGmC5SVap8e3CpMix3SszisDP1bnLiiSeybNkyvF4vNTU1vPTSSy22NyXw3Nxcampqmu+IDwaDFBUVccopp3D33XdTUVFBTU0NW7duZcaMGXzve99j7ty5bN68ucffkzHGmMjsrqjH66pgUNCFxxW7Hjr73Zl6rMybN48LLriAWbNmMWrUKObOnUtGRkbz9szMTL7yla8wY8YMRo8ezbx58wAIBAJcddVVVFZWoqrccsstZGZm8uMf/5g333wTl8vF1KlTOfvss2P11owxxnRiXVEllZ4G8t05MY1DOms27m3mzp2rrcdT37RpE1OmxK65o0lNTQ2pqanU1dWxaNEiHnjgAebM6f3j1PSWz88YY/qqe174gOfLrufszKn8+OK/deuxRWSNqs6NpK6dqXejJUuWsHHjRrxeL9dee22fSOjGGGOO3r6iAqrT4hiVNT6mcVhS70ZPPPFErEMwxhjTw/yBIHWV6yANRg6eFdNY7EY5Y4wx5ih8tr+GJPcuAEYMiW0LrSV1Y4wx5iisK6ogLv4AcQojMkbHNBZL6sYYY8xRWFdcQW18DcPjEoh3xcc0FrumbowxxhyFzTv3EJcSZExi7EcRtaQeJbfffjupqalUVVWxaNEiTj/99JjGs2PHDs477zw++eSTmMZhjDH9SV2jn2DpZooy3SyIcdM7WFKPujvuuCPWIRhjjImST3ZXke/5jG1xcYzOmxHrcOyaenf62c9+xqRJkzj99NPZsmULANddd11zl7B33HEH8+bNY/r06SxZsqS5v/jFixdzyy23sGjRIqZMmcLq1au5+OKLmTBhAj/60Y+aj//rX/+a6dOnM336dH7zm98Azhn4lClT+MpXvsK0adM488wzqa+vB2DNmjXMmjWL448/nvvvv78nPwpjjBkQ1hVVkJHg3Pk+ZmhE/cNEVf87U3/lNti3vnuPOWQGnH1Xh1XWrFnDU089dURDrwLEx8fz9ttvc++993LhhReyZs0asrOzGTduHLfccgs7duzg4YcfZtWqVagqxx13HCeffDJZWVl89tlnPPnkkzz44INceumlPPfcc1x11VVcf/31/Pa3v+Xkk0/mO9/5Tvd+JsYYY/io6CDDkkoBGJ05NsbR2Jl6t3nnnXe46KKLSE5OJj09nQsuuOCwOu0NvQo0158xYwbTpk1j6NChJCQkMHbsWIqKilixYgUXXXQRKSkppKamcvHFF/POO+8AMGbMGGbPdsbvPfbYY9mxYweVlZVUVFRw8sknA86PCGOMMd1HVSnYcZAGTxVpuMhJjG2/79Afz9Q7OaOOJhFpd1tHQ68CJCQkABAXF9e83LTu9/s7HNo1vL7L5aK+vh5V7TAeY4wxR6f4YD011ZXszfExJmFwr/iba2fq3WTRokX8/e9/p76+nurqapYtW9Zie3tDr3bl+P/4xz+oq6ujtraWv//975x00knt1s/MzCQjI4MVK1YA8Pjjj3fxHRljjOlIwc5yxstutnvcjE4dEetwgP54ph4jc+bM4bLLLmP27NmMGjXqsITb3tCrXTn+ddddx/z58wG44YYbOOaYY9ixY0e7+zz88MN86UtfIjk5mc997nNdfk/GGGPaV7DjIFMSinjZ7WZ0bu8Y6dKGXjX2+RljzBE46zdvczH389usjfx60T2cMSY6J09dGXrVmt+NMcaYLqqs97GlpJpEz24AxmSOi3FEDkvqxhhjTBd9tOsgqkqFlhIHjEwfGeuQAEvqxhhjTJet2XmQIXGV7BQfwz3pMR/IpYkldWOMMaaLCnYc5MycUrbGexiXPjrW4TSLalIXkbNEZIuIFIrIbW1sHykib4rIRyLysYicE814jDHGmKPlCwRZW1TB8Wl72enxMGHQrFiH1CxqSV1EXMD9wNnAVOAKEZnaqtqPgGdU9RjgcuD/ohWPMcYY0x027a2i3hcgzVWIX4RxudNiHVKzaJ6pzwcKVXWbqjYCTwEXtqqjQHpoOQPYE8V4+qTly5dz3nnntbnthhtuYOPGjQD8/Oc/78mwjDFmwPpgezkAFY07ABifOT6G0bQUzaQ+HCgKWy8OlYW7HbhKRIqBl4GvRzGemPP7/d16vIceeoipU53GD0vqxhjTM1ZuK2diTjzbG8twIYzJGBPrkJpFM6m31Qlu655urgAeUdV84BzgryJyWEwiskRECkSkoLS0NAqhdo8777yTyZMnc8YZZ3DFFVdwzz33sHjxYn7wgx9w8sknc++997Js2TKOO+44jjnmGE4//XRKSkoAuP3227n66qs59dRTmTBhAg8++GDzcWtqavjCF77A5MmTufLKK1sM2VpQUMBtt91GfX09s2fP5sorrwTgscceY/78+cyePZsbb7yRQCDQ8x+IMcb0M4Gg8sH2Ms4bVsVWj4sRCdm95s53iG43scVAeGe4+RzevP5l4CwAVX1fRBKBXGB/eCVVfQB4AJwe5Tp60V9+8Es2l28+ushbmZw9me/N/16HdQoKCnjuuefaHHq1oqKCt956C4CDBw+ycuVKRISHHnqIu+++m1/96lcAfPzxx6xcuZLa2lqOOeYYzj33XAA++ugjNmzYwLBhw1i4cCHvvvsuJ554YvNr33XXXfzud79j7dq1gNND3NNPP827776Lx+Phpptu4vHHH+eaa67p1s/FGGMGmk17q6jy+jkxrYTX6j1MzJoQ65BaiGZSXw1MEJExwG6cG+G+2KrOLuA04BERmQIkAr33VLwDK1as4MILLyQpKQmgeZx0gMsuu6x5ubi4mMsuu4y9e/fS2NjImDGHmm2a9k9KSuKUU07hgw8+IDMzk/nz55Ofnw/A7Nmz2bFjR4uk3trrr7/OmjVrmvuXr6+vZ9CgQd36fo0xZiBaua0MgFG6jSK3m7PzZsY4opailtRV1S8iNwOvAi7gz6q6QUTuAApUdSnwbeBBEbkFp2n+Oj3Kzug7O6OOlo7CTklJaV7++te/zre+9S0uuOACli9fzu233968rfWwfU3rrYdW7ezavKpy7bXX8otf/KIrb8EYY0wnVm4rY0xuCvsrPiYowvjsibEOqYWoPqeuqi+r6kRVHaeqPwuV/SSU0FHVjaq6UFVnqepsVX0tmvFE04knnsiyZcvwer3U1NTw0ksvtVmvsrKS4cOd+wX/8pe/tNj2wgsv4PV6KSsrY/ny5V0ayc3j8eDz+QA47bTTePbZZ9m/37mKUV5ezs6dO4/kbRljjAkJBJVV28tZMCaLwoqtQO+68x2sR7luM2/ePC644AJmzZrFxRdfzNy5c8nIyDis3u23384ll1zCSSedRG5ubott8+fP59xzz2XBggX8+Mc/ZtiwYRG//pIlS5g5cyZXXnklU6dO5ac//SlnnnkmM2fO5IwzzmDv3r1H/R6NMWYg27S3imqvn1OH1FGID7fE9Zo+35vY0KvdqKamhtTUVOrq6li0aBEPPPAAc+bMiWjf22+/ndTUVG699dYoR3m43vL5GWNMb/bQO9v46UubWHtxDT9a80N2D5rA3//jlai/bleGXo3mjXIDzpIlS9i4cSNer5drr7024oRujDGm92u6np5Z+QGfxXuYmdP7ToYsqXejJ554BefDUgAAIABJREFU4oj3Db9hzhhjTO/SdD39vJlDqdqzhj1uN5fktO75PPbsmroxxhjTifW7K6n2+lkwJpstZU733JOzJ8c4qsP1m6Te1+4N6C3sczPGmM6986nThcqivDo+pQGwpB41iYmJlJWVWYLqIlWlrKyMxMTEWIdijDG92juFB5g+PJ2sqk1sjo8nJz6d3KTcznfsYf3imnp+fj7FxcX05n7he6vExMTm3uqMMcYcrqbBz4c7D3LDSWNhzxtsiY9nci+8ng79JKl7PJ4W3a0aY4wx3WXl1jL8QWXRhFx8731EYXw8V/fSpN4vmt+NMcaYaFlReIBETxzHjspkW+l6fNI7r6eDJXVjjDGmQ29/VspxY3JIqNrJFq0HYFL2pBhH1TZL6sYYY0w7dlfUs620lpMm5EJxAZvjPSTGxTMqbVSsQ2uTJXVjjDGmHSs+Cz3KNjEPilezJTGJidmTcMW5YhxZ2yypG2OMMe14+7MDDE5PYMKgVHT3ajbHJzCpl15PB0vqxhhjTJv8gSDvfFrKogl5iN/L3tJNVEmQSVm983o6WFI3xhhj2rRm50GqvH5OmzII9q5jg8dJmdNyp8U4svZZUjfGGGPa8Mbm/XhcwokTnOvpnyTE4xY3E7Mmxjq0dllSN8YYY9rw+ub9HDcmh9QENxSvZkNKBpOyJxHvio91aO2ypG6MMca0squsjsL9NZw6eRAAweICNnhcTM+dHuPIOmZJ3RhjjGnljc0lAE5Sr9rDzvoSaggwLaf3Xk8HS+rGGGPMYV7fvJ+xeSmMzk1xrqfHO03uvfkmObCkbowxxrRQ2+Bn1bZyTp3kNL2zayUbkpJJciUyNmNsbIPrhCV1Y4wxJsyKwgM0BoKcOiWU1He+xyepWUzJmYo7rncPbmpJ3RhjjAnzr40lpCe6mTc6G7xV+PZ9zOY4f69veocIk7qIZB/JwUXkLBHZIiKFInJbO3UuFZGNIrJBRJ44ktcxxhhjuoMvEOTfm0o4fcpgPK44KP6ArR4XDRpgek7vvvMdINJ2hFUishZ4GHhFVbWzHUTEBdwPnAEUA6tFZKmqbgyrMwH4PrBQVQ+KyKAuvwNjjDGmm3ywvZyKOh+fmz7EKdj5HusTkwB6/eNsEHnz+0TgAeBqoFBEfi4inXWpMx8oVNVtqtoIPAVc2KrOV4D7VfUggKrujzx0Y4wxpnv985N9JHlcLJqQ5xTsfJ91mYPJTsxmRNqI2AYXgYiSujr+papXADcA1wIfiMhbInJ8O7sNB4rC1otDZeEmAhNF5F0RWSkiZ3UxfmOMMaZbBIPKqxv2sXhSHknxLvB5YXcB6+LdzMqbhYjEOsRORdT8LiI5wFU4Z+olwNeBpcBs4G/AmLZ2a6OsdbO9G5gALAbygXdEZLqqVrR6/SXAEoCRI0dGErIxxhjTJR8VVbC/uoGzmpre93xIufrZEajlokGzYxtchCJtfn8fSAc+r6rnqurzqupX1QLgD+3sUwyEt1XkA3vaqPOCqvpUdTuwBSfJt6CqD6jqXFWdm5eXF2HIxhhjTORe3bAPj0s4ZXLTo2zvsi4xAYDZef0rqf9IVe9U1eKmAhG5BEBVf9nOPquBCSIyRkTigctxzu7D/QM4JXS8XJzm+G1diN8YY4w5aqrKPz/Zx8LxuaQnepzCne+zNmsY7jg3U3OmxjbACEWa1Nt6HO37He2gqn7gZuBVYBPwjKpuEJE7ROSCULVXgTIR2Qi8CXxHVcsijMkYY4zpFhv2VLGrvI6zpoWa3v2NsGsla5NTmJozlUR3YmwDjFCH19RF5GzgHGC4iNwXtikd8Hd2cFV9GXi5VdlPwpYV+FZoMsYYY2Ji2bo9uOOEzzUl9d1r8Plq2RCo4bI+0vQOnd8otwcoAC4A1oSVVwO3RCsoY4wxpqcEg8rSdXtYNDGPrJTQWOnblrMpIZEG9TO7j9wkB50kdVVdB6wTkcdDzenGGGNMv1Kw8yB7K71876zJhwq3LWftoDFAbZ+5SQ46b35/RlUvBT4SkfDH0QSn9XxmVKMzxhhjomzput0keuI4Y+pgp8BbBcWr+WjKfIa7MslL7jtPXXXW/P7N0Py8aAdijDHG9DRfIMjL6/dx2pTBpCSEUuLO9whqgAJfBScPOz22AXZRh3e/q+re0OIBoEhVdwIJwCwOf+bcGGOM6VPeLTxAeW0jF8wadqhw23IKk1Ko8Ncyb8i82AV3BCJ9pO1tIFFEhgOvA9cDj0QrKGOMMaYnLF23h7REN4snhTWxb3+L1UMmAfTbpC6qWgdcDPxWVS8C+saT+MYYY0wb6hr9vPrJPs6ePoQEt8sprC6B/RspSElleOpwhqUO6/ggvUzEST00cMuVwEuhskiHbTXGGGN6nVfW76O2McAXjg3r0XzbcoLA6sbSPneWDpEn9W/i9CD391CvcGNxeoAzxhhj+qRn1xQzKieZeaOzDhV+9iqfpQ+m0tf3rqdDhGfbqvo2znX1pvVtwDeiFZQxxhgTTUXldby/rYxvnzHx0JCqAT8U/puCMceAdxtzB8+NbZBHINKhVycCtwKjw/dR1VOjE5YxxhgTPc+uKUYELj42/1Bh8WrwVrI6MYHh7r53PR0ivy7+N5whVh8CAtELxxhjjImuYFB57sNiFo7LZXhm0qENn71KQFysrivm1JGnxS7AoxBpUver6u+jGokxxhjTA1ZuL6P4YD23njmp5YZPX2PDyDlUNZZwwrATYhPcUYr0RrllInKTiAwVkeymKaqRGWOMMVHwzOoi0hLch0ZkA6gshv0beDd3OIKwYOiC2AV4FCI9U782NP9OWJkCY7s3HGOMMSZ6ymoaeHn9Pq6YP4KkeNehDZ+9BsD7WsvUnKlkJWa1c4TeLdK738dEOxBjjDEm2v62ppjGQJArF4xqueHT16jOHMnHFYV8afqXYhNcN4io+V1EkkXkRyLyQGh9gojYIC/GGGP6jGBQeWLVLuaPzmbi4LRDGxprYdubrBo5m4AGWDh8YeyCPEqRXlN/GGgEmu4cKAZ+GpWIjDHGmCh4p/AAu8rruHLByJYbCv8Nfi/vpaaS4klhZl7fHVU80qQ+TlXvBnwAqlqPM6a6McYY0yc8tnInOSnxnDV9SMsNm5ahSdm8V7WN+UPm44nzxCbAbhBpUm8UkSScm+MQkXFAQ9SiMsYYY7rR3sp6Xt9UwqXzRhwavAXA3wCfvsqOCYvZXbuHhcP6btM7RH73++3AP4ERIvI4sBBn+FVjjDGm13v0/Z0AfHF+q6b37W9DQxVvZQ+Bajgp/6QYRNd9Ir37/TURWQMswGl2/6aqHohqZMYYY0w3qG3w8/jKnZw1fQgjspNbbtz4AsSn8WbdbiZlTeqTXcOGi/Tu99dVtUxVX1LVF1X1gIi8Hu3gjDHGmKP13IfFVHn9fPnEVl2rBPyw5WUOTjiVtQfWsXjE4pjE1506TOoikhjqOS5XRLLCepMbDfTtnzPGGGP6vWBQ+fOK7cwekcmxo1p1KLPrPagr451BYwhqkFNGnBKbILtRZ2fqNwJrgMmhedP0AnB/ZwcXkbNEZIuIFIrIbR3U+4KIqIj0vXHujDHG9Fqvb97PjrI6bjipjT7U1j8LnhSWByoYlDSIKTlTej7AbtbhNXVVvRe4V0S+rqq/7cqBRcSFk/jPwHmufbWILFXVja3qpeGMzb6qS5EbY4wxnXjonW0Mz0zirGmtHmPzN8DGF2icdDbv7l3JOWPPIU4ifSCs94roHajqb0XkBBH5oohc0zR1stt8oFBVt6lqI/AUcGEb9e4E7ga8XYrcGGOM6cCanQdZtb2c6xeOxu1qle4KXwdvBatHzqbOX9cvmt4hwrvfReSvwDhgLYfGU1fg0Q52Gw4Uha0XA8e1Ou4xwAhVfVFEbo00aGOMMaYz979ZSFayhytaP8YGsP4ZSM7hX75Skt3JzB8yv+cDjIJIn1OfC0xVVe3Csdvqca55fxGJA/4XuK7TA4ksAZYAjBzZxj+OMcYYE+aT3ZW8sXk/t545kZSEVqmuoRq2vIJv9pW8XvQmJ484mUR3YmwC7WaRXkD4BBjSaa2WioERYev5wJ6w9TRgOrBcRHbgPAO/tK2b5VT1AVWdq6pz8/LyuhiGMcaYgeb/lheSluDm6uNHH75x04vg97I6fzoVDRV8bvTnejy+aIn0TD0X2CgiHxDWPayqXtDBPquBCSIyBtgNXA58MWzfytBxARCR5cCtqloQcfTGGGNMK5+VVPPKJ/v42uLxZCS10Y/7+mcgcySv1e4k2Z3c57uGDdeVbmK7RFX9InIz8CrgAv6sqhtE5A6gQFWXdvWYxhhjTGf+b/lWEt0uvnRiG4+xVRTB1jfxnfRtXi96jcUjFvebpneIvJvYt47k4Kr6MvByq7KftFN38ZG8hjHGGNPks5JqXli7my+fOIbslPjDK6x9AsBpei9+pl81vUMnSV1EVqjqiSJSTdhNbjg3wamqpkc1OmOMMaYLfvXapyTHu/nq4vGHbwwG4aPHYOzJvFq2lhRPCguH95+md+jkRjlVPTE0T1PV9LApzRK6McaY3mRdUQX/3LCPG05q5yx9+3Ko3IV39hW8tuM1Tht5GgmuhB6PM5r6fvc5xhhjDPA/r24hOyWeG04a23aFD/8KSVksT06hxlfDBeM6ute7b7Kkbowxps97r/AAKwoPcNPicaS2fi4doK4cNr8IMy9n6Y5XGJIyhHlD5vV8oFFmSd0YY0yfFgwqv/znZoZlJHLVglFtV/rorxBo5MC0C3hvz3ucN/a8ftHXe2v97x0ZY4wZUP6xdjfriiv51pmTSPS4Dq8QDMAHD8GoE3m5upCABjh/7Pk9H2gPsKRujDGmz6pt8HPXK5uZlZ/BxccMb7vSllegchc6fwlLty5les50xma2c929j7Okbowxps/6/fKt7K9u4CfnTyMurq0hR4AP/gjp+XySO4otB7fw+fGf79kge5AldWOMMX1SUXkdD7yzjc/PHsaxo7LarlSyEba/DfO+zN8KnyfJncS5Y8/t2UB7kCV1Y4wxfdLPXtqES4TvnT25/Uof/BHciVRNv5hXtr/CuWPPJTU+teeC7GGW1I0xxvQ5r28q4Z8b9nHzqeMZmpHUdqXqElj7JMy8lGX73sUb8HLJxEt6NtAeZkndGGNMn1Lb4OcnL2xg4uBUvtJeRzMAq34PgUb0+G/w7KfPMj1nOlNzpvZcoDFgSd0YY0yf8ut/fcruinp+ftEM4t3tpDFvJaz+E0y9kIJABYUVhVwyqX+fpYMldWOMMX3I+uJKHn53O188biRzR2e3X3H1n6ChCk68hUc3PEpWQhbnjDmn5wKNEUvqxhhj+oRGf5DvPvcxOakJfO+sDm6O89XDyt/DuFPZkZLJ8uLlXDb5sn41bnp7IhpP3RhjjIm1e1//lE17q3jwmrlkJHnar7jmEajdDyfewmObHiM+Lp7LJl3WY3HGkp2pG2OM6fXW7DzI75dv5ZJj8zlj6uD2KzbWwju/gtEnUTF0Ji8UvsB5484jNym354KNITtTN8YY06vVNfq59W/rGJqRxE/O7+Tu9VV/hNpSuOxxnt7yNN6Al6unXN0zgfYCdqZujDGmV/vZS5vYfqCWey6ZRVpiB83u3kp4916YcCZ1Q2fw2KbHOGn4SYzPGt9zwcaYJXVjjDG91rJ1e3h81S5uXDSW48fldFz5/fvBWwGn/JAnNz9JRUMF/2/W/+uZQHsJS+rGGGN6pR0Havn+8+uZMzKTWz83qePKVXvgvd/C1Aupy5vIXzb8hYXDFzIzb2bPBNtLWFI3xhjT63h9Ab72xIe4XcJvvzgHj6uTdPX6HRD0w+n/zdNbnuZgw0G+OuurPRNsL2JJ3RhjTK+iqty+dAMb9lTxq0tmMTyznb7dm+z+ENY9CQtuoi5tMI9seIQThp3ArLxZPRNwL2JJ3RhjTK/y6Ps7eWp1ETctHsdpUzp4fA1AFf75fUjJg5O+zV82/IVybzk3zb6pZ4LtZaKa1EXkLBHZIiKFInJbG9u/JSIbReRjEXldREZFMx5jjDG927uFB7jjxY2cPmUQt57ZyXV0gI+fhqKVcOqPOKCNPLzhYc4YdcaAPEuHKD6nLiIu4H7gDKAYWC0iS1V1Y1i1j4C5qlonIl8F7gZ6rNufFwpf4Kcrf4qIIAhxEocgiLS/3FQvTpzfQ03LbnHjjmt78oin7W2hfTxxHjwuDwmuhMMndxtlbWxzx1mXA8aYvm3HgVpuevxDxuWl8JvLjyEuTjreoa4cXv0B5M+DY67h/lV34gv6+M85/9kzAfdC0cwE84FCVd0GICJPARcCzUldVd8Mq78SuCqK8RxmbMZYrph8BUENomjzXLXVMkFU9VCd8GWUYDCIX/34gj78QX/z1OBvoDZYi18PlfmCvsPqNZUresTvxSUuEt2JJLuTSfYkk+xOJsmdRJInySkLK0/2ONtalyW7k0nxpJAan0pafBoJroRu/LSNMaZ95bWNfOmR1YjAQ9fMIzUhgvT02o+dZ9PPv5fCym08/9nzXDH5Ckamj4x+wL1UNJP6cKAobL0YOK6D+l8GXoliPIeZkTeDGXkzevIl26Wq+IN+vAEvDYEGGgONeANeZ+73tlxvp7zeX0+9v546fx31Pmde6a1kr39vc3mdrw5f0BdRTPFx8c0JPtXjzMOXU+NTSfOkHbacmZBJRkIGSe4kRDr5pW2MGfDqGv1c/8hqdlfU89gNxzEyJ7nznba/A2sfg4X/iQ6ayt3/upEUdwo3zrwx+gH3YtFM6m39NW/zVFRErgLmAie3s30JsARg5Mj++QtMRPC4nGb4NNKi+lq+gM9J/K1+ANT6aqnx1VDTWEONr4bqxmqqG6upaayh2ucsl9aVNi/X++s7fJ34uHgyEzJJT0gnMyGzOdmHL7e1bpcSjBk4fIEgNz3+IeuLK/jDVccyr6PhVJs01MDSr0PmKDj5e7y641Xe3/s+t82/jazErOgH3YtF869nMTAibD0f2NO6koicDvwQOFlVG9o6kKo+ADwAMHfu3CNvozYAeFweMlxOAj0a/qCfWl8tVY1VzT8EqhqqqGiooKKhgsqGSiobK6nwOuvbK7c3l/vV3+5xMxIyyEnMITsxm+zEbHKSclrOE3Oc7UnZJLuTrTXAmD4qGFS+99zHLN9Syi8unsGZ04ZEtuNrP4SDO+D6l6khyN2r72ZqzlQun3R5VOPtC6KZ1FcDE0RkDLAbuBz4YngFETkG+CNwlqruj2IsJgrcce7ms+uuUFVqfbVOwm+ooNJb2fxDoKKhgnJvOeXecsrqy/j04KeU7S2jurG6zWMluhIPS/55yXkMShrkzJMHkZeUR05SjrUAGNOLNCX05z/cza1nTuSK+RG2wm75pzO06sJvwqgT+N0Hd3Gg/gC/PfW3uOJcUY25L4jaXzlV9YvIzcCrgAv4s6puEJE7gAJVXQr8D5AK/C10trVLVS+IVkyHaaiB+oMgAhLnTDQth80J295U3laZnTFGRERIjU8lNT6V4anDI9rHF/A5id5b1pzww5N/ubec/XX72Vi2kTJvGUENttg/TuLIScxpkfDbSv5ZiVnNTzYYY6IjGFS+//x6/rammG+eNoGbT50Q2Y61B5xm98HT4ZQfsq50HU9ufpJLJ13KtNxp0Q26jxDVvtWaPXfuXC0oKOieg334qPMF6S7iApcH4jzgckOcO2zZE9oWKm+qFxe+T2hyJ4I7IWye1Go9fN5WWQJ4kiE+2ZkPsB8bgWDASfL1+ymtK2V/3X5K60sprSulpK6E0rpSSutLKfeWH7avO87N4OTBDEkZwtCUoQxNGcqQlCEt1lPjU2PwrozpHwJB5Yd/X89Tq4v4xqnjueWMiZFdQgsG4PEvwI4VsGQ59TljuXTZpTQEGnj+guf79f9LEVmjqnMjqTuw2yNHngAX/A40CKgz16DTQ5FqG2XBVnU1rDzg9Dsc8LWcB30QCM2D/kPLzdv94Ks/VC/QCP4G8HsPzQNt3moQIWmZ4ONTnKlpuWlbfGpYvZRDZQnpkJgOCWnOckKaUx7Xe89mXXGu5jNxOhjUyRfwcaD+QIvkX1JXwr7afeyr3ceHJR9SUldCQAMt9kv1pDYn+fD5kJQhDE8dzuDkwdYMaEwbGvwBvvX0Ol5av5ebT+lCQgd4627Y+gacfy8MnsZvVv2CHVU7+NOZf+rXCb2rBnZSzx3vTL1dMBhK9t5D89aJP3zZ5wVfHTTWHpqHLzfN68rBF9rWWOcst2q2bpu0TPKtk35iemg5tJ6U1WrKdFoTYszj8jA0dShDU4e2WycQDFBaX9qc6PfW7m0x/+TAJxxsONhiH7e4GZo6lOGpwxmeOpz8tHzyU/OblzMTMu3mPjPgVHt93PjXNby3tYwfnDOZJYvGRb7zZ/+Gt34Js74Ic65l5d6VPLH5Ca6cciXzh86PXtB90MBufjctqTo/DBprnQTfUAMN1aGp0pl7q0Lrobm3Mmw5bJvf2/FreZIPT/RJWZCY2caPgCxIzoGU3F7xY6C1en89JbUl7Kndw56aPeyu2U1xdTG7a3azu2b3Yc38ye5khqcNb5Ho81PznXlavnX6Y/qd0uoGrn/kAzbtrebu/5jJfxybH/nOZVvhodMgfTh8+V8cCNTxhaVfIC0+jWfOf4YkdyeDvfQD1vxujowIeBKdqaN260j4Gw/9GKivcG5IrD8I3qblsLL6CjhQGFoud1oj2hOf5iT3lFxIzoWUnNA8r2VZSp6z7Ek8uvcRgSR3EqMzRjM6Y3Sb22t9tYcl+uLqYoqqi3h/z/t4A4d+AAnCkJQhjEwfyai0Uc48fRSj0keRn5qPx+WJ+vsxpjt9sruSrzxawMG6Rh685lhOndzJAC3h6srhiUsBgUsfJeBO4Ltv3kytr5YHznxgQCT0rrKkbqLDHQ/uHCfBdoWqc4+BNyzp15VDXRnUHYDaMqgtdZYri2HvWueO2PZ6yYtPDZ3l50HqYEgdBGlDnHnqEEgb7JSnDHJijoIUTwoTsyYyMWtiG29XKfOWsbtmN0XVReyq2sXOqp3sqtrFKzteafEoX5zEMSxlGKPSDyX7kWnOfFjqMHtkz/Q6L328l2//bS1ZyfE8+/9OYPrwLjz+6m+Ap6+Cil1wzVLIGcf9H97H6n2ruXPhnW3+fzKW1E1vIxK6SS8Z0odFto+qcxmgrsxJ8HUHWs5rDzg/BA7ugKJVTnlbkrJDCT+U6JsSfvP6ECem+JRufLtCblIuuUm5h40qpapUNFQ4Sb76ULLfWbWTtaVrqfXVNtd1i5v8tHzGZoxlTMYYxmaObV5O8XRfvMZEwh8I8pt/f8bv3ixkzshM/nj1XPLSunBZKRiEF26Gne/CxQ/BqON5bcdrPLj+QS4afxGfH//56AXfx1lSN32fSOiafCbkRHDzTcDnJPnqfVBT4kzVJVCzD2r2O+VlW531ti4FJGY41/fSh4WmNpYT0o/6UUIRISsxi6zELGYPmt1iW9MZflOS31m1kx1VO9hWuY23i99u0WPfoORBjM04lOTHZoxlbOZYchJz7IY90+32VXr5xlMf8cH2ci6dm8+dn59OgrsLT4OowkvfgvXPwKk/hpmXsL50PT9Y8QNm583mhwt+GL3g+wG7Uc6Y9qg6zf81+50EX10C1Xugqmna7cxrSg7fNz61ZaJPG3poOXMkZI5wngyIAl/QR1F1EdsrtrO9ajvbKraxrXIb2yu3U+eva66XFp92KMlnjGVc5jgmZE5gSMoQS/bmiLyxuYRvP7OOBn+Qn35+OhfP6cINceD8n3v1h7DyfjjxFjjtv9hbu48rXrqCRHciT5z7BNmJEfQN38905UY5S+rGHC1/o5P0wxN96+XqvYc/LpiYGUrwoSljhJPsm5aTsrq14yBVpaSupDnBb6/czrbKbWyr2EaZt6y5XqonlfGZ4xmfNZ7xmeOZmDWR8ZnjB/xAGaZ91V4fv3hlM0+s2sXkIWncf+Uc/n97dx5dZ1kncPz7e++aNHvSljZpuiVNKmXpFKVsCsWFQafAAQ9w4AgjICpzxJVRRhyc0TN6VMBBcRlUGNyP4hErR6YWEBBBCgjS9qZZ2qbpln1P7vrMH897b+5N09I0N0l78/sc7nmX+9zLm/e8vb9nf1bOn+TYcWNg813w/P1w9ofhkq/QNdrNjX+4ka6RLh659BFWlkxiGFwO0aCu1IkmHoOhdtu5r7cV+vZC7960/VY7f0A6f6EN8sVLxkr3xUvsylSlyyC/LGtBvy/cR3NvM029Tezs2UlTbxONPY30R/pTacqD5dSW1lJTUpPa1pTUkO87hmUyVc56trGDz/767+zvG+Hm85fzqXfXEfRNcvKlRBw2fQJeeRjeegtc+jX6Iv3c9MRNtA608r13fY+1C9ZOzx9wEtCgrtTJxhjby7+v1Qb43r1jwb53rz0/2pf5mUARlC6F0uVQttxuS5fZ/aIqOz3xlC7J0DHSQVNPE429jalA39zbnDEMr7KgktqSWmpKa6grq6OutI7qwmqdVS/HdQ2G+eofQvxyaxsr5s/ja1edwbqlx1GbEwvDb26Fbb+BCz4NGz5PX6Sfj/zxI4S6Q3zr4m9x7uJzs/8HnEQ0qCuVi0b7xkr3PbuhZxd077L7vXsyO/U5XluqHx/sk/uB459WM2ES7BvYR2NvI409Ntg39Taxu293qoNenjeP2pJaVpWtor60nrqyOlaVrtJSfQ6IxRP89K+tfP2JBoYjcW6+YAUff2ft5EvnAIMd8Ivr7KiUd/0nnPcxOkc6uXXzrezq28U33vENLqq+KPt/xElGg7pSc00ibtvuxwf75P5ob2b6efNtgC+vgfIV7rYGylYc95C9SDxCc28zDT0NNHQ3pLbJKnxBWFK4JFWaryuro76snoX5C7Vj3kniucZOvvz4DnYc6Oe8mnK+uPFUahYcZ4fPg2/Az66xQ07BYIzHAAAS50lEQVSv+A6cegUHBg9wy+ZbaB9u576L7pvzJfQkDepKqUwjPW6Q3+0GfDfYdzXbHv3pChfboYHlNWnbGtuWP8kJeowxHBw6SENPA6HuEDt7dtLQ3UDrQGsqTZG/iPqyelaVrqK+zJbqVxav1NnzTiCvtvbwtScaeL65i8qSPD7/3tVcsmYKoyRe+4VtQw8WwbU/g8Vreb3jdW5/6nbC8TAPXPzAYcM45zIN6kqpYxcZgu4W6GpyX2n7I2nz1ovHdtjLCPYroWwlFFfZZYSP0VB0iMaeRkLdoVSJvrGnMdVW73W8rCxeSV1ZHavLVtvSfVkdRf6ibP/16ij+3tbH/U828n/bD1E+z89tF9Vw3frqyY07TxcehMc/A6/91K6SedUPoWgRm1o28e9//ncW5C/g/g33U1N6Eiy0NYM0qCulsmO4Oy3gN48F++4WiAyOpfMEbICvqIWKOqhYZffLa465/T6eiNM60Jqqut/RvYOG7gY6R8ZmAKwsqKSutI768nrqS+upL6vXcfVZZozh2cZOvvunZp5v7qIw6OWWC1bwwfOXUxCYQufL1hfgt7fZ5+gd/wpv/wxh4tz78r38ZMdPOGvhWdxz4T06dHICGtSVUtPLGDvpTirYN9pFeTp32qr99DH5RVVusF+Vtl1lp909hmDcOdJJQ/dYkA91h9jTvweD/e1KVt8nX3VldSwvXo7P0er7yRiOxNj02gEeen432w/0s7AowE3nL+fat1VTGJzCvRzthy1fhJcehOJquPwBWH4BLb0t3PHMHTT0NHD96uv55LpPapPLEWhQV0rNnljYttd37nRfjWPbtAVq8BdOHOzLVrxp2/1wdJjG3kZCXSFCPSEauhvY2bOTcDxsv9rxU1NaY4N8aV2qzb7Af/y9/nNVw8EBfvriHh59dR8DozFqFxRwywUruGzt4uOvZgeb8Xvj17D5C7YT5/qPwEX/RtQX4Mfbf8wDf3uAPG8eXzr/S7y96u3Z+4NykAZ1pdSJxxg7s15GoHf3+/eNpROPHX43PthX1NpZ9o4glojR2t+aUaIPdYfoCfek0lQXVqd63ScD/oL8BXOu+r69f5TfvX6Ax17bz2t7e/F7HC497RSuW7+Us5aWTv1+tL4AT9wJ+16GU06D994LS97KG51vcPfzd9PQ08CFSy7krvV3sSB/QXb+qBymQV0pdXIJD9iq/GSw72iw+93NmePv581Pa6+vHWu3L1k64WQ7xhjah9tTve9D3aHDet+XBctSpflkx7ylRUtzbvKc9oFRtuxo53ev7ecvLV0YA6cuLuLyMyu5cl0VZfOysPRw6wvwzNehabNd72DDXXDGNewfPsT9r97P71t+T0VeBXeefScXV1885zJTx0uD+jH64/ZDfPvpJgS7Ipbd2vG07n+pY5HM/aTMz2Uek/zcuO9xxG49juARwUnfOmSc8ziCyOHnbdr0z9vvc8Z9Z/o5r0fweQSv47j7Dh5H8KWOM9/zOoLX46TO+zyi/wjVzIrH7MQ640v2nQ12mF6Sx2+r7ctrDg/4+YcvADIYGWRnz85U7/tQd4jGnkaiiSgAQU+Q2tJaW6ovrae+vJ7aktqTavIcYwzb9vfzZKidLTsO8VqbnZFwWXk+G8+sZOMZi6lZkIXmiEQCmrfAc/fBnucgvxzOuQ3O/jCHooM8vP1hfh76OY44XLf6Om4+7WYK/dOzmFGu0qB+jJ5uaOcHz+3CGDAYu03fBxh3nLxfdj+ZxqQdp39PWvq09xLGbuPGEE8YEgnj7tv3Ms9lpp1tHjcz4RsX8CfKCPg8Dn6Pg9/r4PM4BLx2P3nOP+444KabKE0g7Tg9TWDc93g9zmzfIjVThrrcDnqNh3fUS4wtPUt+xViAT1bnl9faKXbTOmZFE1F29e1KVd0nO+elT56ztGjpYZ3yKvIqZvovn5AxhpbOIV5o6eKFlm5eaOmiYyCMCJy5pISL6xewoX4hqxcVZidzPtgBrz4CL//IznJYuAjO/Risu4Hm4YM8tO0hNrVsImESbFy5kdvOvI1T5p0y9f/vHKRBPYelB/uxDACHnct4380wROMJYglDLJ4gGjfEEglicTN2zt3G4oao+97kPpOZLhpPEInZVziWIBLPPBdJncveM+hxhKDXIejzEEhufR6CPoeg124D7jbo87jvp53zetzz49ONHQfS0gR9HryO1mCcUOJR6NkzccAfHhseh+O1s+qlAn7aMLz8chBJTZ4T6h7rkBfqDrFvcKwPQEVeRcZ4+vrSeqqLqnFkejOYg+EY2/b18XpbH39r6+Wvu7rpGLAdBRcWBVi/opzzaiq4qG4B8wsD2fmfjvZB6PfwxqPQ8pTNPC27AM76IMM1F7O57WkebXyUV9pfIegJckXtFXzgLR+gqnCSS7CqDBrU1UklkTBE4jbAR2ITZwbGZwQi48/Hk+njjEYTjEbdbSxOOJo8P/ZeOJZME2fU/Y7jlcxI5Pm95Pkd8n1egn4P+T4PeX735fOQ727Tj4M+D/nu5/J8XvL8h6fL83lwHM00ZMVIjw3wGQG/0Y67T2+79xdC2TI7T35qwRx3v3gJ/fFhO54+bahdc29z5tz3pbWsLF7JsuJlLC9azrLiZVQVVk16qJ0xhv19ozQeGqDx0CChgwO83tZLU8cgyZ/vypI8zlpWyvoV5axfUc6y8vzsZDSNsfeneQs0bYFdf7L3qbga1lzB8JoreTZ8kC2tW3im7RmGokMsLVrK5TWXc2XtlTrmPEtOmKAuIpcA3wQ8wIPGmK+Mez8A/C+wDugCrjbG7D7ad2pQV9MhmbFID/zJDEEy8CczAeFYgvC4dCORBCPROCORGCPROMMRm3Y4EnfPx1PnjycDEfA6mcHe7zlq5iE4QSZi/Dbf7yXP5yHot00Xc7q2IRF32+7dyXWS8+b37LalfneoHODOrLckc6Gc4ioiBQtpkTihcBehXjsd7u7+3RmT53jFS1VhFcuLl6eC/dKipSyatxgTK2Jf7yhtPcPs7R5hT9cwTR2DNB0aYCgST33H/MIAp1cWc3pVCadXFXNaVTEVBVkqicejcGgb7NsKbVth93N2tUCAspVEa9/N9uq1vBQf4KVDL7H14FYiiQilgVIuXHIhG1duZN3CdXP7WZoGJ0RQFxEPsBN4F9AGvARca4zZnpbmo8DpxpgPi8g1wBXGmKuP9r0a1NXJLp4wY4E+FexjhwX/VKYgLWMwdt5NH02kMhLJ7xuOxpnsP2uPIzbAp2UGMjIMx7DN+FxaLUPy/ZO2v0MiYefHT583P7W/O3MqXQBxoOAUKK4kXriYzmA5ISfADhOnMTJEa6yXg9FO+hLtGMaCtUl4MNFSEu5rnjOfRfMWs6K0ivr5lZy+aAlrFlVQmo1e6tFRO4ywqwk6Qna0QUcIDm2H2AgG6CuYz97Fp7GzYhk7An4ahvbT0NPASGwEgJqSGtYvWs+G6g2sXbAWrzO1pX7VkZ0oQf0c4G5jzHvc488BGGP+Ky3NE26av4iIFzgIzDdHuSgN6kodnTGGcCwxlhmITpwxSM88JM+PRNMzDQlG0zIdo9FExv5k+TwyLjPgJc+X7Ojoyej46POI2/nRk9YRUtI6UXpS6QJeByc50sMd+eGIHS0iqfO458U9T0afk4ThsA6qCWP7q0Tjthkno1knrSbHjPYRGD5AYOgAeaMHKQgfoiTaQUW8g4V0sVi6yJPIYfcjjEOjv5imQDH7g/m0B/y0+6DTiXGQMD2J8GGfKfAEqfAXUx4opiJQQrGvgAJPkELHT4H7KnS8+OIxvNFhvJFhPJFhvJEhGOklMnSIyFAH4cgAURGGHIcex6E3r4ie/BK6A/nsd2BfdIDB2HDq/zvPN4+60jpWl69m3cJ1rFu4jrLg4aMK1PSYTFCfzqxVJbA37bgNOPtIaYwxMRHpA8qBTpRSx0VEUp0Ap6tFM5GwGYdkjcFEtQoZGYpxNRDp++Fogv6R6BH7SSS3Jxq/10l1ysz3l1Gct5CiEh9FQR9FeT6K8rx2P+hlgW+YxU4fFc4ApaaPQLiLwHAna4Y6WDPUaZcfHe6CvkE7Zj8yyLAI+7xe2r0eOjweOj0eujwOnZ4uOjweQl4P/Y7DgOMQO9bq7gAQyAPyMk57xUtJMEBJoJDFBYtZV1BFZUElVYVV1JbWUllQOe0d/1R2TGdQn+gpG18CP5Y0iMiHgA8BVFdXT/3KlFJT4jiSasefCcYYonFzxE6TtrSdLGHjlrJtqdskS+FmrPSdMGasdJ+c28EtzSdL/MkSvs8zNgIiOToi4HWmt/NiIkF+ZJDa8AC14QGIDtl2/3gUElHb6zweAxPHOD7CGAaJMZCIMZSIEPX4ifmC9iUOcWOr+f2OH5/HR8ATwO/xk+/NpzRYSoGvQNvBc8R0BvU2YEnacRWw/whp2tzq92JgXAMVGGO+D3wfbPX7tFytUuqEJSL4vbb6nSz1CTuhOY5dazz45kvNChB0XyfGiHk1m6azPuUloFZElouIH7gGeGxcmseAG9z9q4Anj9aerpRSSqkjm7aSuttG/i/AE9ghbT80xmwTkf8AthpjHgN+ADwiIk3YEvo103U9SimlVK6b1jEIxpjHgcfHnftC2v4o8P7pvAallFJqrtDujEoppVSO0KCulFJK5QgN6koppVSOOOkWdBGRDmBPFr+yAp3sZqbovZ4Zep9nht7nmaH3GZYaY+YfS8KTLqhnm4hsPdbp99TU6L2eGXqfZ4be55mh93lytPpdKaWUyhEa1JVSSqkcoUHdnX5WzQi91zND7/PM0Ps8M/Q+T8Kcb1NXSimlcoWW1JVSSqkcMaeDuohcIiINItIkIp+d7evJFSKyRESeEpEdIrJNRG53z5eJyGYRaXS307Xc95wiIh4ReVVENrnHy0XkRfc+/8JdUElNgYiUiMivRCTkPtfn6POcfSLyCfc34w0R+ZmIBPV5npw5G9RFxAN8G/hH4C3AtSLyltm9qpwRAz5ljFkNrAduc+/tZ4EtxphaYIt7rKbudmBH2vFXgXvd+9wD3DQrV5Vbvgn8wRhTD5yBvd/6PGeRiFQCHwPOMsaswS4Edg36PE/KnA3qwNuAJmNMizEmAvwcuGyWryknGGMOGGNecfcHsD+Aldj7+7Cb7GHg8tm5wtwhIlXAe4EH3WMBNgC/cpPofZ4iESkC3o5dVRJjTMQY04s+z9PBC+SJiBfIBw6gz/OkzOWgXgnsTTtuc8+pLBKRZcBa4EVgoTHmANjADyyYvSvLGfcBdwAJ97gc6DXGxNxjfa6nbgXQAfzIbeZ4UETmoc9zVhlj9gFfB1qxwbwPeBl9nidlLgd1meCcDgXIIhEpAH4NfNwY0z/b15NrROR9QLsx5uX00xMk1ed6arzAPwDfMcasBYbQqvasc/skXAYsBxYD87DNo+Pp83wUczmotwFL0o6rgP2zdC05R0R82ID+E2PMo+7pQyKyyH1/EdA+W9eXI84DNorIbmzz0QZsyb3Erb4Efa6zoQ1oM8a86B7/Chvk9XnOrncCu4wxHcaYKPAocC76PE/KXA7qLwG1bs9KP7ZDxmOzfE05wW3X/QGwwxhzT9pbjwE3uPs3AL+d6WvLJcaYzxljqowxy7DP75PGmOuAp4Cr3GR6n6fIGHMQ2Cside6pi4Ht6POcba3AehHJd39DkvdZn+dJmNOTz4jIpdiSjQf4oTHmy7N8STlBRM4HngX+zlhb753YdvVfAtXYf8DvN8Z0z8pF5hgRuRD4tDHmfSKyAltyLwNeBa43xoRn8/pOdiJyJrYzoh9oAf4ZWyjS5zmLROSLwNXYETSvAjdj29D1eT5GczqoK6WUUrlkLle/K6WUUjlFg7pSSimVIzSoK6WUUjlCg7pSSimVIzSoK6WUUjlCg7pSc5C76thH3f0Lkyu8TeLzN4rI4um5OqXU8dKgrtTcVAJ8dAqfvxE7ladS6gSi49SVmoNEJLkqYQMQxc5n3gmswS6icb0xxojIOuAeoMB9/0bs9LQPAfuAEeAc4DPAPwF5wPPArUZ/XJSacRrUlZqD3NXzNhlj1riz0f0WOBU7r/afsUH6ReBPwGXGmA4RuRp4jzHmgyLyNHYGu63u95UlZ1MTkUeAXxpjfjezf5VSyvvmSZRSc8BfjTFtACLyN2AZ0IstuW+2U3HjwS6JOZGLROQO7BrYZcA2QIO6UjNMg7pSCiB9Lu049rdBgG3GmHOO9kERCQIPAGcZY/aKyN1AcLouVCl1ZNpRTqm5aQAofJM0DcB8ETkH7HK6InLqBJ9PBvBOESlgbEUtpdQM05K6UnOQMaZLRP4sIm9gO7sdmiBNRESuAv5bRIqxvxf3YavWHwK+KyLJjnL/g12Vbzd2WWOl1CzQjnJKKaVUjtDqd6WUUipHaFBXSimlcoQGdaWUUipHaFBXSimlcoQGdaWUUipHaFBXSimlcoQGdaWUUipHaFBXSimlcsT/A7O2SuFiV1QoAAAAAElFTkSuQmCC\n",
"text/plain": "<Figure size 576x576 with 2 Axes>"
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "import scipy.constants as sc\n\ndef reduced_mass(m1, m2):\n return (m1*m2)/(m1 + m2)\n\ndef freq(J, m, re):\n num = sc.hbar * J # kilograms meters^2 / second\n denom = 2*np.pi*m*re**2 # kilograms meters^2\n return num / denom # Hz\n\nre = 1.13e-10 # meters\nm_C = 12.0107 # g/mol\nm_O = 15.999 # g/mol\nm_gmol = reduced_mass(m_C, m_O)\nm = m_gmol / (sc.Avogadro * 1000) # kg\n\nf = freq(np.array((1,2,3)), m, re)/1e9 # GHz\n\nprint(\"Frequency J=1 to J=0: {:.4} GHz\".format(f[0]))\nprint(\"Frequency J=2 to J=1: {:.4} GHz\".format(f[1]))\nprint(\"Frequency J=3 to J=2: {:.4} GHz\".format(f[2]))\nprint(\"Differences: {:.4}, {:.4} \".format(f[2]-f[1], f[1]-f[0]))",
"execution_count": 109,
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": "Frequency J=1 to J=0: 115.4 GHz\nFrequency J=2 to J=1: 230.8 GHz\nFrequency J=3 to J=2: 346.1 GHz\nDifferences: 115.4, 115.4 \n"
}
]
}
],
"metadata": {
"kernelspec": {
"name": "python3",
"display_name": "Python 3",
"language": "python"
},
"language_info": {
"name": "python",
"version": "3.7.2",
"mimetype": "text/x-python",
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"pygments_lexer": "ipython3",
"nbconvert_exporter": "python",
"file_extension": ".py"
},
"gist": {
"id": "",
"data": {
"description": "PHY 530 homework 2",
"public": false
}
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment