Skip to content

Instantly share code, notes, and snippets.

@Diimu
Created September 12, 2019 20:15
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 Diimu/363fc88d9dfd454e67fb91a2c6adb0f2 to your computer and use it in GitHub Desktop.
Save Diimu/363fc88d9dfd454e67fb91a2c6adb0f2 to your computer and use it in GitHub Desktop.
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgIAAAHnCAYAAADKJDVgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdd5xcdb3/8ddne7bvZjeb3c0mm7rpCUmIhJooJSAQpQiogKBGFK9erw1sKFcvth/lXgVEAQEVjBRFek1ACQQSQnonpPdkS9q2z++PmYQlpGyyM3umvJ+Pxzx25syZ7/l8OGHOZ77ne77H3B0RERFJTilBByAiIiLBUSEgIiKSxFQIiIiIJDEVAiIiIklMhYCIiEgSUyEgIiKSxAIvBMzsHjPbZGbz2iwrNrPnzWxp+G9RkDGKiIgkqsALAeCPwMQDll0HvOju/YEXw69FREQkwiwWJhQys2rgCXcfGn69GBjv7uvNrByY6u41AYYoIiKSkGKhR+Bgytx9PUD4b7eA4xEREUlIaUEH0BFmNhmYDJCVlTW6Z8+eR/X5Pc2wYVcr3bNTyIqz/xKtra2kpMRqHRdZyZQrJFe+yZQrJFe+yZQrBJfvkiVLtrh7aUfaiNXD30YzK29zamDTwVZy97uAuwBqamp88eLFR7WRZZvqOf3mV7jt0pFMGlnZ4aA709SpUxk/fnzQYXSKZMoVkivfZMoVkivfZMoVgsvXzN7raBuxWq49DlwZfn4l8I9obKQ0LwuATXV7o9G8iIhIzAu8EDCzB4HpQI2ZrTGzzwM/B84ws6XAGeHXEZeflUZmWgqbG1QIiIhIcgr81IC7X3aItz4W7W2bGaV5mWyq2xPtTYmIiMSkwHsEglZekMX6WhUCIiKSnJK+EOhe0IUN6hEQEZEklfSFQEW4RyAWJlYSERHpbElfCHQvyKKxuZVtOxuDDkVERKTTJX0hUF4QuoRQ4wRERCQZJX0h0L2gCwAbVAiIiEgSSvpCoGJ/j8DugCMRERHpfElfCHTNzSQtxXRqQEREklLSFwKpKUZZvuYSEBGR5JT0hQDsm1RIpwZERCT5qBAgdAmhBguKiEgyUiHA+9MMa1IhERFJNioEgPKCLuxtbmX7rqagQxEREelUKgRoO6mQxgmIiEhyUSFAaIwAwPodGicgIiLJRYUAUFEYml1wnXoEREQkyagQAEpzM8lIS2HtdhUCIiKSXFQIACkpRo/CLqzevivoUERERDqVCoGwyqIurFGPgIiIJBkVAmFVxdkqBEREJOmoEAjrUdSFbTsb2bm3OehQREREOo0KgbAeRdkA6hUQEZGkokIgrKoodAnhGg0YFBGRJKJCIGxfj8DqbSoEREQkeagQCCvJzSArPUWnBkREJKmoEAgzM3oUZWsuARERSSoqBNroobkEREQkyagQaEOFgIiIJBsVAm1UFWVTu7uJuj1NQYciIiLSKVQItLHvyoFVWzVOQEREkoMKgTZ6dQ0XArqEUEREkkRMFwJm9g0zm29m88zsQTPLiub2qktyAHh3y85obkZERCRmxGwhYGaVwNeAMe4+FEgFLo3mNnMz0yjNy2SlCgEREUkSMVsIhKUBXcwsDcgG1kV7g9Vds3lPYwRERCRJmLsHHcMhmdnXgZ8Bu4Hn3P0zB7w/GZgMUFpaOnrKlCkd3ubdc/cyZ0sLt03I7nBb0dTQ0EBubm7QYXSKZMoVkivfZMoVkivfZMoVgst3woQJM919TEfaSItUMJFmZkXAJKA3sAP4m5l91t3/tG8dd78LuAugpqbGx48f3+HtzvdlvPrsYo4fdzI5mTH7n4epU6cSiXzjQTLlCsmVbzLlCsmVbzLlCvGdbyyfGjgdeNfdN7t7E/AocGK0N1rdNTRgcOVWjRMQEZHEF8uFwCrgBDPLNjMDPgYsjPZGq0tCpwRWbtE4ARERSXwxWwi4+xvAw8AsYC6hWO+K9nZ7qUdARESSSOyeBAfc/Qbghs7cpi4hFBGRZBKzPQJB6t01Rz0CIiKSFFQIHESvrtms1FwCIiKSBFQIHER1SQ6b6/dSr7sQiohIglMhcBD9uoUmhVi+WacHREQksakQOIh9hcCyTQ0BRyIiIhJdKgQOoldxNumppkJAREQSngqBg0hLTaF3SQ7LNtUHHYqIiEhUqRA4hH7dctUjICIiCU+FwCH065bHqm272NPUEnQoIiIiURORQsDMeplZSfj5CWb2LTP7ZCTaDkq/brm0OryrGQZFRCSBdbgQMLMfAi8Br5vZT4FbgRLga2Z2a0fbD0r/8JUDS3V6QEREElgk7jVwGTAIyCZ0x8Du7r7LzNKA2RFoPxC9S3JIMV1CKCIiiS0ShcAed28EGs1subvvAnD3ZjNrjED7gchKT6WqOJvlKgRERCSBRaIQKDSzCwAD8sPPCb8uiED7genfLZeluoRQREQSWCQKgWnAeeHnr7R5vu913OrXLY9pSzbT1NJKeqousBARkcTT4ULA3a+KRCCxaGD3PJpanBWbd1LTPS/ocERERCKuw4WAmf3X4d5395s7uo2gDCrPB2Dh+joVAiIikpAi0d+dF36MAb4MVIYf1wCDI9B+YPqU5pCRmsLC9XVBhyIiIhIVkTg18BMAM3sOGOXu9eHXPwb+1tH2g5SemkK/brks3KABgyIikpgiOQKuJ9D2csFGoDqC7QdiYHmeegRERCRhRbIQeACYYWY/NrMbgDeA+yLYfiAGl+ezuX4vWxr2Bh2KiIhIxEWsEHD3nwFXAduBHcBV7n5TpNoPyr4Bg4vW6/SAiIgknkjMI7Cfu88CZkWyzaANDF8tsGhDHSf3Lwk4GhERkcjSLDlH0DU3k255mSzQOAEREUlAKgTaYWB5vk4NiIhIQlIh0A6DyvNYuqmexubWoEMRERGJqKgVAmb2gpk9bWbnRmsbnWVYZQFNLc6SjeoVEBGRxBLNHoErgB8AvaK4jU4xvLIQgDlragOOREREJLKiUgiYWRFQ4u4z3f230dhGZ6oq7kJBl3TmrNkRdCgiIiIRFbFCwMymmlm+mRUD7wD3mlnc3nCoLTNjeI8C9QiIiEjCiWSPQIG71wEXAPe6+2jg9I40aGaFZvawmS0ys4VmNi4ikR6D4T0KWLKxnj1NLUGFICIiEnGRLATSzKwc+BTwRITavA14xt0HAiOAhRFq96gNqyykudV13wEREUkokSwEbgSeBZa5+5tm1gdYeqyNmVk+cCpwN4C7N7p7YCfph/coADRgUEREEkvEphh297/R5rbD7r4CuLADTfYBNhMaazACmAl83d13dijQY1RekEVJboYKARERSSjm7pFpyOyXwE+B3cAzhLry/9Pd/3SM7Y0BXgdOcvc3zOw2oM7df9hmncnAZIDS0tLRU6ZM6WAWh3fzzD1s3d3Kz07Ojup22qOhoYHc3Nygw+gUyZQrJFe+yZQrJFe+yZQrBJfvhAkTZrr7mA414u4ReQCzw38/Sej2w8XAOx1orzuwss3rU4AnD7X+gAEDPNpufm6x977uCW/Y0xT1bR3Jyy+/HHQInSaZcnVPrnyTKVf35Mo3mXJ1Dy5f4C3v4PE7kmME0sN/zwEedPdtHWnM3TcAq82sJrzoY8CCjrTZUSN7FtLq8I7mExARkQQRyULgn2a2CBgDvGhmpcCeDrb5H8CfzWwOMBL4nw621yGjqooAmLlye5BhiIiIREwkBwteZ2a/IHQev8XMdgKTOtjmbEKFRUwoyE6nf7dcZq5SISAiIokhYoWAmaUDlwOnmhnANODOSLUfK8ZUF/HknPW0tjopKRZ0OCIiIh0SyVMDdwCjgdvDj1HhZQllVM8i6vY0s3xzQ9ChiIiIdFjEegSA4919RJvXL5nZOxFsPyaM7hUaJ/DWe9vpX5YXcDQiIiIdE8kegRYz67vvRXhmwYSbmL93SQ7FORnMfE/jBEREJP5Fskfg28DLZrYCMKAXcFUE248JZsaonkXMUiEgIiIJIJJXDbxoZv2BGkKFwCJ33xup9mPJmOoiXli4kW07GynOyQg6HBERkWPW4ULAzC44xFt9zQx3f7Sj24g1Y8LjBGa8u42JQ7sHHI2IiMixi0SPwHmHec+BhCsEhvcopEt6Kq+v2KpCQERE4lqHCwF3T7hxAEeSkZbCmOoipi/fGnQoIiIiHRLJqwaSygl9urJ4Yz1bGxJyGISIiCQJFQLHaFzfrgC8vqJD91YSEREJlAqBYzSssoCcjFSmr9gSdCgiIiLHLJLzCGBmJwLVbdt19/sjuY1YkZ6awvG9izVOQERE4lokbzr0ANAXmM37Mwo6kJCFAMC4Pl2Zungzm+r20C0/K+hwREREjlokewTGAIPd3SPYZkw7oU9onMD0FVuZNLIy4GhERESOXiTHCMwDkuqi+qGVBRR0SefVpRonICIi8SmSPQIlwAIzmwHsv6bO3c+P4DZiSmqKcUr/EqYt2Yy7Y2ZBhyQiInJUIlkI/DiCbcWN0waU8sSc9SxcX8/givygwxERETkqETs14O7TgEVAXvixMLwsoZ02oBSAaUs2BxyJiIjI0YtYIWBmnwJmABcDnwLeMLOLItV+rOqWn8Wg8nymLdkUdCgiIiJHLZKnBr4PHO/umwDMrBR4AXg4gtuISacNKOUPr66gYW8zuZkRnZpBREQkqiJ51UDKviIgbGuE249Zpw0opbnVeW2Zrh4QEZH4EskD9TNm9qyZfc7MPgc8CTwVwfZj1uheReRkpDJV4wRERCTORKwf292/bWYXAicBBtzl7o9Fqv1YlpGWwin9S3lp4SZaJzkpKbqMUERE4kNET2i7+yPAI5FsM16cOaSMZ+ZvYO7aWkZUFQYdjoiISLt0+NSAmf0r/LfezOraPOrNrK7jIcaHjw7sRmqK8dyCDUGHIiIi0m4dLgTc/eTw3zx3z2/zyHP3pJlhpzA7g7HVxTw3f2PQoYiIiLRbJOcReKA9yxLZmUPKWLqpgRWbG4IORUREpF0iedXAkLYvzCwNGB3B9mPeGYPLAHh+gXoFREQkPkRijMD1ZlYPDG87PgDYCPyjwxHGkR5F2QypyFchICIicSMSYwRucvc84FcHjA/o6u7XRyDGuHLm4O7MXLWdjXV7gg5FRETkiCJ5amCGmRXse2FmhWb2iY42amapZva2mT3R0bY6w8eHl+MOT8xZH3QoIiIiRxTJQuAGd6/d98LddwA3RKDdrwMLI9BOp+jXLZfB5fn88511QYciIiJyRBG918BBlnVowiIz6wF8HPhDR9rpbOePrGD26h2s3rYr6FBEREQOK5KFwFtmdrOZ9TWzPmZ2CzCzg23eCnwHaO14eJ3n48PKAfjnHPUKiIhIbDN3j0xDZjnAD4HTw4ueA37m7juPsb1zgXPc/StmNh74lrufe8A6k4HJAKWlpaOnTJlyrOFH3E9f383eFvjvk7pEpf2GhgZyc3Oj0nasSaZcIbnyTaZcIbnyTaZcIbh8J0yYMNPdx3SkjYgVAvsbNMt19w7PqGNmNwGXA81AFpAPPOrunz3Y+jU1Nb548eKObjZi/vjvd/nxPxfw/DdOpX9ZXsTbnzp1KuPHj494u7EomXKF5Mo3mXKF5Mo3mXKF4PI1sw4XApGcWfBEM1sALAi/HmFmtx9re+5+vbv3cPdq4FLgpUMVAbHonOHlpBg89vbaoEMRERE5pEiOEbgFOAvYCuDu7wCnRrD9uNItL4vxNd14ZNYaWloj2+siIiISKZEsBHD31QcsaolQu1MPHB8QDy4e3YONdXt5ZenmoEMRERE5qEgWAqvN7ETAzSzDzL5FHF3/Hw0fG1RGcU4Gf3vrwPpIREQkNkSyELgGuBaoBNYAI8Ovk1ZGWgqfGFnJ8ws2sm1nY9DhiIiIfEjECgF33+Lun3H3Mnfv5u6fdfetkWo/Xl08pgdNLc4/ZmvQoIiIxJ4OzfwHYGb/BxxyNJy7f62j24hng8rzGVZZwF/fXM3nTqzGzIIOSUREZL8OFwLAWxFoI6F9+iM9uf7Rubz13naOry4OOhwREZH9InFq4HR3vw8odPf7DnxEoP24N2lkBXlZadz32sqgQxEREfmASBQCo82sF3C1mRWZWXHbRwTaj3vZGWl8akwVz8zbwKa6PUGHIyIisl8kCoE7gWeAgYRuMtT2odMGYZef0IvmVucvM1YFHYqIiMh+HS4E3P1/3X0QcI+793H33m0efSIQY0KoLsnhtAGl/OWNVTS1xNXNFEVEJIFF8vLBL5vZyWZ2FYCZlZhZ70i1nwiuPLEXm+r38vS8DUGHIiIiAkT2pkM3AN8Frg8vygD+FKn2E8H4Ad3oU5rD76YtJ9J3fRQRETkWkZxZ8JPA+cBOAHdfB0T+/rtxLCXF+NKpfZi/ro5/L0v6uZZERCQGRLIQaPTQz1wHMLOcCLadMD5xXCXd8jK5c9ryoEMRERGJaCEwxcx+BxSa2ReBF4A/RLD9hJCZlsrVJ/fmX8u2MG9tbdDhiIhIkovkYMFfAw8DjwA1wI/c/X8j1X4i+fRHepKXmaZeARERCVwkewRw9+fd/dvu/i3gJTP7TCTbTxT5Wel8dlwvnpy7nqUb64MOR0REkliHCwEzyzez683sN2Z2poV8FVgBfKrjISamL57Sh+z0VG59cWnQoYiISBKLRI/AA4ROBcwFvgA8B1wMTHL3SRFoPyEV52Rw9cm9eXLOehaurws6HBERSVKRKAT6uPvn3P13wGXAGOBcd58dgbYT2hdO7kNeVhq3PL8k6FBERCRJRaIQaNr3xN1bgHfdXSe+26EgO50vnNyH5xZsZO4aXUEgIiKdLxKFwAgzqws/6oHh+56bmfq8j+Dqk6spzE7nF88s0myDIiLS6SJx06FUd88PP/LcPa3N8/xIBJnI8rLS+frH+vOvZVuYunhz0OGIiEiSiejlg3JsPntCL/qU5PDTJxfozoQiItKpVAjEgPTUFK47eyDLN+/koRmrgg5HRERiVGurs3LLTp6cs55fPrMoIm2mRaQV6bAzBpdxQp9ibnlhKeePrKSgS3rQIYmISICaW1pZtrmB+WvrmL+ujnnralm4ro76vc0ApKVYRLajQiBGmBk/+Phgzv/Nv/h/zy3mxklDgw5JREQ6yd7mFpZsaGDu2lrmratl/tpaFm2oZ29z6HRxVnoKg8rzmXRcBUMrChhSUcCA7rlk3dTxbasQiCFDKwu4Ylw1901fyYWjejCiqjDokEREJMJ2N7awcEMd89fWhg78a+tYsrGe5tbQlWP5WWkMqSjginG9GFJRwJCKfPqU5pIaoR6AA6kQiDHfPHMAT81dz/cem8s/rj2JtFQN4xARiVcNe5tZuL6OuWv2/dKvY9nmBlrCB/2i7HSGVhYwuaYPQysLGFpRQFVxF8yic9A/GBUCMSYvK50bzhvCtX+Zxf3T3+Pqk3sHHZKIiLRD7e4m5ocP9vPWhX7tv7tlJ/umiCnNy2RYZQFnDSljSGUBwyoLKC/I6tSD/sGoEIhB5wzrzviaUv7fc4s5Y3AZVcXZQYckIiJtbN/ZuP9gP39tHW8u38WmZ57b/35FQRZDKguYNKKSYT3yGVpRQLf8rAAjPrSYLQTMrAq4H+gOtAJ3ufttwUbVOcyMn35iKBNvfZXvPDyHP3/hI6RE6dyQiIgc3ub6vcxbWxt6rAud01+7Y/f+96uKu9AzL4UrT+kX7t7Pp2tuZoARH52YLQSAZuCb7j7LzPKAmWb2vLsvCDqwztCjKJsfnjuI7z4yl/unr+RzJ+kUgYhINLk7G+r2MG9tXfiXfujAv7Fu7/51+pTkMKpXEVeM68WwytDo/YLsdKZOncr48f0CjP7YxWwh4O7rgfXh5/VmthCoBJKiEAD41Jgqnpm3gZ8/s4hTB5TSpzQ36JBERBKCu7Nm++4P/Mqft7aWrTsbAUgx6Fuay4l9S/b/yh9ckU9eVuLN8RKzhUBbZlYNHAe8EWwkncvM+PmFwznzllf4xpR3+NuXxpGRpqsIRESORmur8962XR/q3q/dHbp5blqK0b8sj48O7BY66FcWMKg8j+yMuDhEdpjF+h3vzCwXmAb8zN0fPeC9ycBkgNLS0tFTpkwJIMLom7Ghmdtn72VidRqXDgydd2poaCA3Nzl6CJIpV0iufJMpV0iufIPKtdWdDTudlXWtvFfbwsq6VlbVt7I7NBkfaQY98lLolZ9CdX4KvQpS6JGbQkZqx8ZhBZXvhAkTZrr7mI60EdOFgJmlA08Az7r7zYdbt6amxhcvXtw5gQXgh3+fxwOvv8fvrxjDGYPLwuejxgcdVqdIplwhufJNplwhufLtjFybWlpZtqmBeWtrmb8udF5/wbo6dje1AJCZlsLgitCI/aGV+aHZ+MryotKzGtS+NbMOFwIx2+9hoQsr7wYWHqkISAbf//ggZq3azrf+9g5Pfu3koMMREelUe5tbWLoxPAXv2lrmratj0fq6/VPwZmekMqQin0uOr2JYuHu/b2mOJmVrh5gtBICTgMuBuWY2O7zse+7+VIAxBSYrPZXffnoU5/7fv7j2z7O4dlDs9uSIiHTEnqYWFq6vY96696fhXbKxnqaW0PdeXmYaQyrzuWJcL4aGR+73LsmJ2hS8iS5mCwF3/xegvdpGdUkO/+9TI/jSAzO5tyWVMz7qgc9IJSLSETvDU/DOW1vL3LV1zF9Xy9JN70/BW5idzrDKAj5/ch+GVuYzrLKAqqJsza0SQTFbCMjBnTWkO986cwC/fm4Jd0xbzlfi9LpVEUk+dXuawrfUrQ0f+GtZ0WYK3pLcDIZWFnDG4DKGhM/rVxZ27rz7yUiFQBy6dkI/Xp2znF89u5j+3fI4Y3BZ0CGJiHxAQ6Pzr6Vb2kzDW8vKrbv2v989P4uhlQWcNyJ0W91hPQrolpepg34AVAjEITPj88My2ZOWzn88OIs/f+EERvcqCjosEUlSWxr2vj8TX/iGO2u272bf1C89irowtKKAi0b32H9OvzQvfqbgTXQqBOJURqpx9+eO56I7XuPqP77Jw9eMo39ZXtBhiUgCc3c21u1tMylP6MC/oW7P/nWqu2YzsqqQE0ubOf+UUQypyKcoJyPAqOVIVAjEsZLcTB74/Ee44I7XuOKeGTzy5ROpKOwSdFgikgDcnbU7du8/2O+bjW9LQ2jefQtPwXtCn+L9s/ENrsgnPzwF79SpUzm5f0mQKUg7qRCIc1XF2dx31Vgu+d10PvuHN3hw8gmUxeitLkUkNrW2Oqu27Qpdo78udFvdeetq2bErNAVvaorRv1supw0oZVhlfngK3nxyMnUISQTaiwlgcEU+9151PFfeM4PL7npdxYCIHFJLq/Pulob9d9ibF56Nr35vaA7e9FSjpnseE4d03/9Lf2D3PLLSUwOOXKJFhUCCGFNdzH1Xj+XKe2Zw6V2v8+AXT6B7gYoBkWTWnil4B5XnM+m4iv231I3WFLwSu1QIJJAx1cXc//mxXHnPm1xy13Tuv3osvbrmBB2WiHSCvc0tLNnQ8P4gvnV1LFxfR+NBpuAdWlnAME3BK2EqBBLM6F7FPPD5sVz1xze58I7X+ONVYxlaWRB0WCISQXuaWliwvu4Dl+t9YArerDSGVhRwZXgK3qGVBVR31RS8cnAqBBLQcT2LePiaE7nynhlc8rvp3Hn5aE7pXxp0WCJyDBrCU/DOXfP+QL5lm9+fgrcoO52hlQV84ZQ++++y17M4WxPzSLupEEhQ/brl8uhXQsXAVfe+yY2ThvLpj/QMOiwROYza3U37p9/d90v/3TZT8JbmZTKssoCzhpQxJPxLv6IgSwd96RAVAgmsLD+LKdeM4z/+8jbfe2wu89fVcsN5QzQQSCQGbKrfw4J1dcxft2/u/TpWbXt/Ct6KgiyGVBbwiZGVDK3MZ2hFAd10NZBEgQqBBJeflc49nzueXz27mDunLWfpxgZ++5lRmt5TpJPsu0Z//ro6nl7SyB/fncH8dXVsrt+7f52q4i4MqyzgkuOrwqP38+maq/9HpXOoEEgCqSnGdWcPZFB5Ht99ZA5n3/Yqt1wyQuMGRCKssbmVJRvrWbC+Lvxrv5aF6+tpCF+jn2rQv2wPp/QvYUhF6IA/qDyfgi7pAUcuyUyFQBKZNLKSmu55fO3Bt7n87hl86dQ+fPPMGp0qEDkG9XuaWLi+nvnravd38S/d9P7I/eyMVAaV53PBqEqGVOQzuLyA9YtncebHTg04cpEPUiGQZAZ2z+cf157MT59cwO9eWcH0FVv51UUjqOmuGxaJHMqmuj3Mb/Mrf/66Ot5rc0vdktwMBlcUcFpNKYPL8xlSkU911xxSDrhcb+syDeqT2KNCIAl1yUjlZ58cxin9S/jeY/M49/9e5doJ/fjK+H7qHZCk1tLqvLd1Z+ga/XV1+3/p77vRDkCvrtkMLs/n4tE9GFIRutFOt7xMjdyXuKVCIIlNHFrO8dXF3PjEAm59YSlPz93ATRcOY1TPoqBDE4m62t1NLFpfx6IN9SxcX8fCDfUs2VC/f/rdtBSjf1ke42tKw137+Qxqc3c9kUShQiDJdc3N5LZLj+O84RX84O/zuOD217hgVCXfnThQNy6ShLDvV/7C9fUs2hCadnfh+nrW7ti9f52i7HQGledz2dieDCrPY1B5Pv3LcslM0412JPGpEBAATh9cxgl9u/Lbl5dx96vv8sy8DVw7oR+fP7m37jomceNIv/JTU4y+pTmM7lXEZ0/oxcDyPAaXq2tfkpsKAdkvNzON704cyKXHV/GzJxfyq2cXc//0lXz1o/25ZEyVxg9IzGhuaWXl1l0sDh/wQ7/0j/wrv1+3XBW2IgdQISAf0qtrDnddMYbpy7dy8/OL+eHf53Hn1OV87WP9+ORxPVQQSKdpbXXWbN/N4o31LAk/Fm+oZ8XmnTS2hO6qp1/5Ih2jQkAOaVzfrkzpM45Xlm7h5ucW891H5nLL80u56qRqLvtITw2akohxdzbW7Q0d8DfUs3hjPUs31rNkY8P+bn2AysIuDCjL5bSaUmrK8hhQlqdf+SIdpEJADsvMOG1AKaf2L2Hqks38/pUV3PT0Iv73xaVccnxPLh/Xi94lOUGHKXHC3dm6s5GlGxtYsrGel+fv5TcLX2PJxnrq9jTvX680L5OasjwuG9uTmu659E8tV2EAACAASURBVC/Lo3+3XPJUfIpEnAoBaRczY0JNNybUdGPe2lr+8OoK7p++knv+/S4n9CnmsrE9OWtId/0yEyA0Un/1tl0s39zAsk0Nbf7upHZ30/71ctJhSA/j/JEV+3/hDyjLoygnI8DoRZKLCgE5akMrC7j10uO4/pxBPDxzDQ+9uYqvPzSbgi7pnDeinHOHV3B8dTGpKTo/m+h2NTazYvNOlm9uYPmmBpZtbmD5pp28u+X9c/gQ+oXftzSH80aU07c0l37dcqkpy2P+zOlMmDAuwAxERIWAHLOy/CyundCPL5/Wl+krtvLQm6t5eOYa/vT6KrrlZXLOsHI+PrycUT2LVBTEsaaWVtZs383KLaED/Htbd/Lu1l0s39TwgVH6KRYaaNq3NIfxA0v3H/D7luRSkH3wLv0FGswnEjgVAtJhKSnGSf1KOKlfCTv3NvPiok08OWcdf5mxij++tpLC7HROG1DKhJpunDqglGJ1+8acppZW1m7fzbtbd7Jyy07e27qLd7fsZOXWnazZvpuWVt+/bm5mGtUl2YypLuKS0ir6dQsd8Ht1zdYEPCJxSIWARFROZhrnj6jg/BEV1O9p4uXFm5m6eBPTFm/mH7PXYQbDKgv4SO9ixvbuyvHVRRRmqzCINndn285GVm/fzZrtu1izfTert+1i9fbdrAof7JvbHOxzMlKpLslhaGUB5w2voFfXbHqX5FBdkkPXnAxdlieSQFQISNTkZaXvLwpaW525a2t5efEmXlu+lfumv8fvX30XgJqyPEb1KmRoZQHDKguo6Z6nX5ZHqbU1NBp/Q+0e1u7YxeptoQN+2wP/rsaWD3ymMDudqqJsBlfkc86wcqpLckIH+645lOTqYC+SLGK6EDCzicBtQCrwB3f/ecAhyTFKSTFGVBUyoqqQ/zwd9jS1MGdNLTPe3cob727jqbkbeHDGaiB0s5cBZXkMrgjNBLd7UzM9NzdQVZxNemryTWbU1NLKpvq9bKjdzfraPWzY96gL/V1fu4dN9XtoavEPfC4vM40exdlUd83h5H6lVBV3oUdRNj2KutCjqIsuxRMRIIYLATNLBX4LnAGsAd40s8fdfUGwkUkkZKWnMrZ3MWN7F/NVQl3Xa7bvZu7aWuatrWXu2lqmLdnMwzPXAHDbrGmkpxo9i7PpWZxNRWEXKgq7UBn+W16QRWleZlxcvtjY3Ert7ibq9jSxbWcjWxv2sqWhka0NjWxp2MvCd/dw++Lp+5e3vdxun6z0FMoLutA9P4uxvYvpXpBFeUEWZflZVBZ2oaoo+5AD9ERE2orZQgAYCyxz9xUAZvYQMAlQIZCAzIyq4myqirM5Z1j5/uV1e5p4+JlXyK+q2X+J2todu5m9egfbd334AJmdkUpxTgZdczIozsmgOCeTvKw0cjJTyc5IIzczjeyMVHIy0+iSkUpGagqpKUZaipGWmkJaiu1/3eqh6+Fb3WlpdVr2/Q0/dje2sLsp9NjT1MKuxhZ2N77/vGFvc+iAv7tp/4G/dncTe5paPxT3PoXZ6XSxVnpmw8Du+XTNzaBrTiZl+Zl0L8gKHfDzu5DfJU1d9yISEbFcCFQCq9u8XgN8JKBYJCD5Wen0KUxl/OgeH3pvV2Mz63bsYd2O3ayv3c2Whka27Qw9tu5sZHPDXhZvqKd+bzO7Gls+MPI9mjLTUuiSkUpeVhr5WekUdEmnb2kuBV3Sye+SFv6bTn5Weqhoyc2gJDeT4pwM0lNTmDp1KuPH69p6EekcsVwIHOznzge+yc1sMjA5/HKvmc2LelSxowTYEnQQnSSZcoXkyjeZcoXkyjeZcoXg8u3V0QZiuRBYA1S1ed0DWNd2BXe/C7gLwMzecvcxnRdesJIp32TKFZIr32TKFZIr32TKFeI731gegv0m0N/MeptZBnAp8HjAMYmIiCSUmO0RcPdmM/sq8Cyhywfvcff5AYclIiKSUGK2EABw96eAp9q5+l3RjCUGJVO+yZQrJFe+yZQrJFe+yZQrxHG+5t45I6lFREQk9sTyGAERERGJsoQoBMxsopktNrNlZnZd0PFEkplVmdnLZrbQzOab2dfDy4vN7HkzWxr+WxR0rJFkZqlm9raZPRF+3dvM3gjn+9fwANK4Z2aFZvawmS0K7+Nxibxvzewb4X/H88zsQTPLSqR9a2b3mNmmtpcyH2p/Wsj/hr+35pjZqOAiP3qHyPVX4X/Lc8zsMTMrbPPe9eFcF5vZWcFEfewOlm+b975lZm5mJeHXcbVv474QaDMV8dnAYOAyMxscbFQR1Qx8090HAScA14bzuw540d37Ay+GXyeSrwML27z+BXBLON/twOcDiSrybgOecfeBwAhCOSfkvjWzSuBrwBh3H0poEPClJNa+/SMw8YBlh9qfZwP9w4/JwB2dFGOk/JEP5/o8MNTdhwNLgOsBwt9ZlwJDwp+5PfzdHU/+yIfzxcyqCE2Fv6rN4rjat3FfCNBmKmJ3bwT2TUWcENx9vbvPCj+vJ3SgqCSU433h1e4DPhFMhJFnZj2AjwN/CL824KPAw+FVEiJfM8sHTgXuBnD3RnffQQLvW0IDlLuYWRqQDawngfatu78CbDtg8aH25yTgfg95HSg0s3LixMFydffn3L05/PJ1QvO/QCjXh9x9r7u/Cywj9N0dNw6xbwFuAb7DBye8i6t9mwiFwMGmIq4MKJaoMrNq4DjgDaDM3ddDqFgAugUXWcTdSuh/rH2T8ncFdrT5gkmUfdwH2AzcGz4N8gczyyFB9627rwV+TeiX03qgFphJYu7btg61PxP9u+tq4Onw84TM1czOB9a6+zsHvBVX+SZCIXDEqYgTgZnlAo8A/+nudUHHEy1mdi6wyd1ntl18kFUTYR+nAaOAO9z9OGAnCXIa4GDC58YnAb2BCiCHUBfqgRJh37ZHov67xsy+T+i05p/3LTrIanGdq5llA98HfnSwtw+yLGbzTYRC4IhTEcc7M0snVAT82d0fDS/euK+rKfx3U1DxRdhJwPlmtpLQaZ6PEuohKAx3J0Pi7OM1wBp3fyP8+mFChUGi7tvTgXfdfbO7NwGPAieSmPu2rUPtz4T87jKzK4Fzgc/4+9enJ2KufQkVte+Ev696ALPMrDtxlm8iFAIJPRVx+Pz43cBCd7+5zVuPA1eGn18J/KOzY4sGd7/e3Xu4ezWhffmSu38GeBm4KLxaQuTr7huA1WZWE170MUK32U7IfUvolMAJZpYd/ne9L9+E27cHONT+fBy4IjzC/ASgdt8phHhlZhOB7wLnu/uuNm89DlxqZplm1pvQILoZQcQYKe4+1927uXt1+PtqDTAq/P91fO1bd4/7B3AOoRGqy4HvBx1PhHM7mVCX0hxgdvhxDqHz5i8CS8N/i4OONQq5jweeCD/vQ+iLYxnwNyAz6PgilONI4K3w/v07UJTI+xb4CbAImAc8AGQm0r4FHiQ0/qGJ0IHh84fan4S6j38b/t6aS+hqisBz6GCuywidG9/3XXVnm/W/H851MXB20PFHIt8D3l8JlMTjvtXMgiIiIkksEU4NiIiIyDFSISAiIpLEVAiIiIgkMRUCIiIiSUyFgIiISBJTISAiIpLEVAiIiIgkMRUCIvIhZvZJM5t9wKPVzM42s4Z2fL7azJ4K33t+iZld3xlxi8jR04RCInJEZjYZ+AwwAahz99zDrJtCaOrvn7j742aWCfwFeNbd7+qUgEWk3VQIiMhhmdkA4CXgRHdfZWYNRygEzga+4O4XtllWDkxz9wHRj1hEjoZODYjIIYXvfPkX4FvuvqqdHxsEfOD+7B664Up++MZgIhJDVAiIyOH8NzDf3R86is+0AB/oMQjfbTCb0D3qRSSGqBAQkYMys/HAhcBXj/KjU4Fzwgf/fc4AZrl7a2SiE5FIUSEgIh9iZkXAvcAV7l5/mPVeNLPKtsvc/R3gbeDG8DplwM3A96IXsYgcq7SgAxCRmHQN0A2444M/7Llp35Pw1QH9gG1tVzCz64AxwGfN7GVCPQq9gNvN7BPuvjK6oYvI0dBVAyJyTMxsKHC1u/9X0LGIyLFTISAiIpLENEZAREQkiakQEBERSWIqBERERJKYCgEREZEkpkJAREQkiakQEBERSWIJM6FQYWGh9+vXL+gwombnzp3k5OQEHUbUKL/4lsj5JXJuoPzi3cyZM7e4e2lH2kiYQqCsrIy33nor6DCiZurUqYwfPz7oMKJG+cW3RM4vkXMD5RfvzOy9jrahUwMiIiJJTIWAiIhIElMhICIiksRUCIiIiCQxFQIiIiJJTIWAiIhIEkuYQmDTLqe5pTXoMEREROJKwhQCu5qdh95cHXQYIiIicSVhCoGsVLjl+SXU7WkKOhQREZG4kTCFQHFWCtt2NfLbl5YFHYqIiEjcSJhCICMVLhzVg3v/vZJVW3cFHY6IiEhcSJhCAODbZ9WQlmrc9PTCoEMRERGJCwlVCJTlZ3HNaX15et4G3lixNehwREREYl5CFQIAXzylD+UFWfz0yYW0tnrQ4YiIiMS0hCsEumSk8t2JA5m7tpZH314bdDgiIiIxLeEKAYDzR1QwoqqQXz27iF2NzUGHIyIiErMSshBISTF+dO4gNtbt5c5pK4IOR0REJGZFtRAws4lmttjMlpnZdQd5/1Qzm2VmzWZ20QHvXWlmS8OPK49226N7FXPu8HLuemU562t3dyQNERGRhBW1QsDMUoHfAmcDg4HLzGzwAautAj4H/OWAzxYDNwAfAcYCN5hZ0dHGcN3ZA2l1+PnTi44+ARERkSQQzR6BscAyd1/h7o3AQ8Cktiu4+0p3nwMceLegs4Dn3X2bu28HngcmHm0APYqyuebUPvxj9jpmvLvt2LIQERFJYGlRbLsSaHsXoDWEfuEf62crD1zJzCYDkwFKS0uZOnXqhxoakup0zTK++Zc3+MmJWaSYtTOE2NLQ0HDQ/BKF8otviZxfIucGyk+iWwgc7Ijb3gv72/VZd78LuAugpqbGx48ff9DGWrqt5yt/nsXarN5cPq66nSHElqlTp3Ko/BKB8otviZxfIucGyk+ie2pgDVDV5nUPYF0nfPZDzh7anRP7duXXzy1h287GY21GREQk4USzEHgT6G9mvc0sA7gUeLydn30WONPMisKDBM8MLzsmZsaPzx9Cw95mfv3c4mNtRkREJOFErRBw92bgq4QO4AuBKe4+38xuNLPzAczseDNbA1wM/M7M5oc/uw34b0LFxJvAjeFlx2xAWR5XjqvmwRmrmLe2tiNNiYiIJIxojhHA3Z8Cnjpg2Y/aPH+TULf/wT57D3BPJOP5zzP68/g7a/nRP+bx8DUnkpISnwMHRUREIiUhZxY8lPysdL4zcSCzVu3gMd2HQEREJLkKAYCLRvVgZFUhNz29iPo9TUGHIyIiEqikKwRSUoyfnD+ErTv38r8vLg06HBERkUAlXSEAMKKqkE+NruLef69kycb6oMMREREJTFIWAgDfmVhDblYaP/j7PNzbO8+RiIhIYknaQqBrbibXTRzIjHe38cgsDRwUEZHklLSFAMCnxlQxulcR//PUQrZrxkEREUlCSV0IpKQYP/3EUGp3N/HLZ3WrYhERST5JXQgADCrP5+qTqnlwxmpmvqdbFYuISHJJ+kIA4D9PH0B5QRbff2wezS2tQYcjIiLSaVQIADmZadxw3hAWbajnj6+tDDocERGRTqNCIOysIWV8dGA3bn5+Cet27A46HBERkU6hQiDMLDTjYKs7N/5zQdDhiIiIdAoVAm1UFWfzHx/tzzPzN/DSoo1BhyMiIhJ1KgQO8MVT+tC/Wy4/eGweDXubgw5HREQkqlQIHCAjLYWfXzic9XV7+PWzi4MOR0REJKpUCBzE6F5FXHFCL+6bvpJZq7YHHY6IiEjUqBA4hG9PHEj3/Cyue2QOjc2aW0BERBKTCoFDyM1M42efHMqSjQ3cOW150OGIiIhEhQqBw/jowDLOH1HBb15axrJN9UGHIyIiEnEqBI7gR+cNJjszlesemUtrqwcdjoiISESpEDiCktxMfvjxwbz13nb+/MZ7QYcjIiISUSoE2uGCUZWc0r+EXzyzWNMPi4hIQlEh0A5mxv98chgtrc4P/z4Pd50iEBGRxBDVQsDMJprZYjNbZmbXHeT9TDP7a/j9N8ysOrw83czuM7O5ZrbQzK6PZpztUVWczTfPHMCLizbx99lrgw5HREQkIqJWCJhZKvBb4GxgMHCZmQ0+YLXPA9vdvR9wC/CL8PKLgUx3HwaMBr60r0gI0lUn9WZ0ryJ+/PgCNtXtCTocERGRDjtsIWBmvcysJPz8BDP7lpl9sp1tjwWWufsKd28EHgImHbDOJOC+8POHgY+ZmQEO5JhZGtAFaATq2rndqElNMX510XD2NLXwvcfm6hSBiIjEvUMWAmb2Q+Al4HUz+ylwK1ACfM3Mbm1H25XA6jav14SXHXQdd28GaoGuhIqCncB6YBXwa3ff1p6Eoq1PaS7fPquGFxZu4rG3dYpARETiW9ph3rsMGARkEzoYd3f3XeFf6bPb0bYdZNmBP6EPtc5YoAWoAIqAV83sBXdf8YEPm00GJgOUlpYyderUdoTVcX3c6V+Ywg8efQfbtISirOiPuWxoaOi0/IKg/OJbIueXyLmB8pPDFwJ7wl36jWa23N13QeiXu5k1tqPtNUBVm9c9gHWHWGdNuMAoALYBnwaecfcmYJOZ/RsYA3ygEHD3u4C7AGpqanz8+PHtCCsyeg/bydm3vcI/N+Rx95VjCJ3RiJ6pU6fSmfl1NuUX3xI5v0TODZSfHH6MQKGZXWBmFwL54ef7Xhe0o+03gf5m1tvMMoBLgccPWOdx4Mrw84uAlzx04n0V8FELyQFOABYdRV5R17skh2+fNZCXFm3ikVk6RSAiIvHpcD0C04Dzws9fafN83+vDCvccfBV4FkgF7nH3+WZ2I/CWuz8O3A08YGbLCPUEXBr++G+Be4F5hE4f3Ovuc9qfVue46sRqnp23gZ/8cz4n9yuhe0FW0CGJiIgclUMWAu5+VUcbd/engKcOWPajNs/3ELpU8MDPNRxseaxJSTF+edFwJt72Ctc9Ood7P3d81E8RiIiIRNIhCwEz+6/DfdDdb458OPGnuiSH704cyE/+uYC/vrmaS8f2DDokERGRdjvcGIG88GMM8GVCl/pVAtcQmiBIwq4cV82Jfbty4xMLeG/rzqDDERERabdDFgLu/hN3/wmhuQNGufs33f2bhGb669FZAcaDlBTj1xePIC3F+MZfZ9Pc0hp0SCIiIu3SngvgexKa2W+fRqA6KtHEsYrCLvz3J4Yya9UO7py2POhwRERE2uVwVw3s8wAww8weIzTZzyd5f1pgaWPSyEpeWLiJW19YyqkDShneozDokERERA7riD0C7v4z4CpgO7ADuMrdb4p2YPHqp5OGUpqXyX/+dTa7G1uCDkdEROSw2jU3rrvPcvfbwo+3ox1UPCvITufXF49gxead/PzphUGHIyIiCaq1NTI3vov+JPlJ6KR+JVx9Um/um/4e05ZsDjocERFJMIs31HP2ba9GpC0VAlHynYk1DCjL5dt/e4ftO9tzawYREZEjm7e2lst+/zrbd0Xm2KJCIEqy0lO55ZKRbN/VyHcfmUPoFgoiIiLH7u1V2/n0718nKy2FKV8aF5E2j7oQMLMXzOxpMzs3IhEksCEVBXx34kCeW7CRP72xKuhwREQkjr2xYiuf/cMbFGZnMOWacVSX5ESk3WPpEbgC+AHQKyIRJLirT+rN+JpS/vuJBSzaUBd0OCIiEof+tXQLV947g7KCLKZ8aRw9irIj1vZRFQJmVgSUuPtMd/9txKJIYPtmHSzoks5X//K2LikUEZGj8tKijVx935tUd83hr5PHRfxOt0csBMxsqpnlm1kx8A5wr5nphkNHoSQ3k1s+NZLlmxu48Yn5QYcjIiJx4um56/nSAzOpKcvjwS+eQGleZsS30Z4egQJ3rwMuAO5199HA6RGPJMGd3L+EL53alwdnrObJOeuDDkdERGLc399ey1cffJvhPQr58xc/QlFORlS2055CIM3MyoFPAU9EJYok8c0zBzCyqpDrHp3D6m27gg5HRERi1F/fXMU3pszm+Ooi7r96LPlZ6VHbVnsKgRuBZ4Fl7v6mmfUBlkYtogSWnprC/112HDh8/aG3adJdCkVE5AD3T1/Jdx+Zyyn9S/njVWPJyWzPbYGOXXvuNfA3dx/u7l8Jv17h7hdGNaoEVlWczc8uGMasVTu45fklQYcjIiIx5M5py/nRP+ZzxuAyfn/FaLLSU6O+zfYMFvxleLBgupm9aGZbzOyzUY8sgZ0/ooJLxlRx+9TlvLx4U9DhiIhIwNydnz+9iJ8/vYjzRlRw+2dGkZkW/SIA2ndq4MzwYMFzgTXAAODbUY0qCfxk0hAGds/jG3+dzdodu4MOR0REAtLS6nzvsXncOW05n/lIT269ZCTpqZ038W97trRvhMI5wIPuvi2K8SSNrPRU7vjsaJpbnGv/PIvGZo0XEBFJNo3NrXz9obd5cMYqvjK+Lz/9xFBSU6xTY2hPIfBPM1sEjAFeNLNSYE90w0oOvUty+OVFw5m9egc36ZbFIiJJZXdjC5MfeIsn5qzn+rMH8p2JAzHr3CIA2jdY8DpgHDDG3ZuAncCkaAeWLM4ZVs5VJ1Vz779X8tRczS8gIpIManc3cfndbzBtyWZuumAYXzqtb2CxHPGaBDNLBy4HTg1XKtOAO6McV1K5/uxBzF69g+88PIdB5fn0jtCNJEREJPZsrt/LFffMYNmmen5z2Sg+Prw80Hjac2rgDmA0cHv4MSq87IjMbKKZLTazZWZ23UHezzSzv4bff8PMqtu8N9zMppvZfDOba2aRnVw5hmSkpfCbT48iLdX48p9msqdJ9yMQEUlEa7bv4uI7X2Pllp384crjAy8CoH2FwPHufqW7vxR+XAUcf6QPmVkq8FvgbGAwcJmZDT5gtc8D2929H3AL8IvwZ9OAPwHXuPsQYDzQ1M6c4lJlYRduvWQkizfW86N/zAs6HBERibBlm+q56I7pbNvZyJ++MJbTBpQGHRLQvkKgxcz2n7wIzyzYnp+sYwnNRrjC3RuBh/jw2IJJwH3h5w8DH7PQ+YczgTnu/g6Au29194T/mTy+phv/MaEfU95aw1/eWBV0OCIiEiFvr9rOxXdOp7nV+euXxjG6V3HQIe3XnnkLvw28bGYrAAN6AVe143OVwOo2r9cAHznUOu7ebGa1QFdCcxW4mT0LlAIPufsv27HNuPf10wcwZ20tNzw+j5rueYzuVRR0SCIi0gEvLdrItX9+m9K8TO6/eizVMTYOzNz9yCuZZQI1hAqBRe6+tx2fuRg4y92/EH59OTDW3f+jzTrzw+usCb9eTqgn4SrgWkKnIHYBLwI/cPcXD9jGZGAyQGlp6egpU6YcMZd4sLPJ+cn03TS2wA3jsijKSqGhoYHc3NygQ4sa5RffEjm/RM4NlF+0vbqmiXvnN1KVl8J/jc6iIDOylwdOmDBhpruP6Ugbh+wRMLMLDvFWXzPD3R89QttrgKo2r3sA6w6xzprwuIACYFt4+TR33xKO5SlCgxQ/UAi4+13AXQA1NTU+fvz4I4QUP/oNq+eTt/+bB1Zk8tDkcbz2r1dIpPwONHXqVOUXxxI5v0TODZRftLg7t09dzt3zFnNyvxLuvHw0uVG+edCxOlxU5x3mPQeOVAi8CfQ3s97AWuBS4NMHrPM4cCUwHbgIeMnd950S+I6ZZQONwGmEBhMmjZruefzqohFc+5dZ/OSf8zlDZwhEROJCS6tz4z/nc9/095g0soJfXTSCjLTOmzL4aB2yEAhfHXDMwuf8v0roFsapwD3uPt/MbgTecvfHgbuBB8xsGaGegEvDn91uZjcTKiYceMrdn+xIPPHo48PLmbu2L3dOW07GkAzGBx2QiIgc1p6mFv5rymyemruBL5zcm++dM4iUTp4y+GhFtZ/C3Z8Cnjpg2Y/aPN8DXHyIz/6J0CWESe3bZ9Uwf10tDyzYwnmrtjOqp7oGRERiUd2eJibf/xavr9jG988ZxBdP7RN0SO0Su30VAkBqivF/lx1HUVZosqFNdbrNg4hIrNlYt4dP3Tmdt1Zu59ZLRsZNEQAqBOJCYXYGXxuVRf2eZr54/1uaeVBEJIYs3VjPBbe/xuptu7jnc8fzieMqgw7pqLSrEDCzE83s02Z2xb5HtAOTD6rKS+GWS0YyZ20t3/rbO7Tnsk8REYmu15Zv4YI7XmNvcysPTR7HqTEyW+DRaM9Nhx4A+gKzeX9GQQfuj2JcchBnDenOd84ayC+eWUT/bnl8/fT+QYckIpK0Hpm5husenUN11xzuvep4ehRlBx3SMWnPYMExwGDXT9CYcM1pfVi6qZ5bXlhC3245nDu8IuiQRESSirtz24tLufWFpZzYtyt3fHY0BV3Sgw7rmLWnEJgHdAfWRzkWaQcz46YLhrFq6y6+OeUdqoqyGVFVGHRYIiJJobG5lesencOjs9Zy4age3HTBsJieI6A92hN9CbDAzJ41s8f3PaIdmBxaZloqd14+mpLcTL54/1usr90ddEgiIgmvdncTV94zg0dnreUbpw/g1xcPj/siANrXI/DjaAchR68kN5O7PzeGC29/jS/e/xZTvjSO7IzYnL5SRCTe/f/27js+qjJ7/PjnpDeSEBJCCSWQ0HsJIsgCVmzIiiJrR7FiXev6s+76ta1lZS1rQRELYEcXFQRRRKR3pERqCCTUBAKpc35/zMWNMQlDmcxMct6v17wy9859hnN4cidn7n3uc7fuOciotxewaXcBz13clT/3SPF1SCfMEUsZVf0eWAPUcx6/OOuMj7VrFMvYv3RnVXY+t01cSpnLhnEYY8yJtjxrH8Ne/okd+YWMH5VRq4oA8KAQEJGLgfm4ZwC8GJgnIsO9HZjxzOB2yTx8bgemr87hsS9W2WWFxhhzqRR8ewAAIABJREFUAn29cgcj/vMz4SFBfHLjyZzcOtHXIZ1wnhxLfgDoraq5ACKSBHwLfOTNwIznruqXyrZ9h3h99kaa1o/kugGtfR2SMcYEtMN3D3zmm7V0axbPa1f0pGG9CF+H5RWeFAJBh4sAx25sRkK/c/+Q9mTvK+T/pq6hSXykXVZojDHHqLCkjPs/WcGnS7YxtFsTnrqwCxGhwb4Oy2s8KQS+dm4L/IGzPIIKNxIyvhcUJDx7cVdy9xdy56RlNKwXQUZqgq/DMsaYgLJzfxHXT1jI4i37uOuMNtw8KA0R/7574PHyZLDg3cBrQBegK/Caqt7r7cDM0YsIDeb1K3qRkhDJ6HcWkpm739chGWNMwFidnc8FL81h9fZ8Xrm0B2MGp9f6IgA8PMSvqh+r6p2qeoeqfurtoMyxi48KY/zVGYQGB3HluAXk7re7FRpjzJFMX53D8Fd/osylfHTDyQzp3NjXIdWYKgsBEfnR+blfRPLLPfaLSH7NhWiOVrOEKMZd1Ys9BcVc/dYC8gtLfB2SMcb4JVXl1e9/5boJC0lvGMOUMf3o1DTO12HVqCoLAVXt7/ysp6qx5R71VDW25kI0x6JLSjyvXNaDtTv2M3q83brYGGMqKiwp464Pl/PkV2s4p3NjJl3fl4axtfPKgOp4Mo/ABE/WGf8zsG1Dnr24K/M37WHM+0soLXP5OiRjjPEL2/MOMeK1n/l4cRa3n5bO2JHda/WVAdXxZIxAx/ILIhIC9PROOOZEG9qtKY+d35Fvf8nh3o9X4LLZB40xddyCTXs4b+wcMnP285/Le3L7aW3qxKDAqlR5+aCI3A/8DYgsNyZAgGLcVxGYAHF535bsKSjh+W/XER8Vyv87p32d/qU3xtRNqsp787bwyJRVNEuI4oPRfUhPrufrsHyuykJAVZ8AnhCRJ1T1/hqMyXjBraemsfdgMW/+uJGE6DBuHpTm65CMMabGFJWW8fDnq5i4YCuD2ibxwiXdiYsM9XVYfsGTCYXmi0icquYBiEg8MFBVP/NuaOZEEhEeOrcD+w4W88w3a4mLDOWyk1r4OixjjPG6nPxCbnh3EUu27GPMoDTuOL0NwUF2VPQwTwqBh8vPHaCq+0TkYcAKgQATFCQ8c1FX8gtLefDzlUSEBjO8Z+26i5YxxpSXubeMu8f+SEFRKa9c2qNOzQ/gKU8GC1a2jd34PkCFBgfx8qU96Nc6kXs+WsaUZdm+DskYY044VeX9eVt4Yn4hUWHBfHpTPysCquBJIbBQRJ4TkdYi0kpEngcWefLmInKWiKwVkUwRua+S18NFZJLz+jwRaVnh9eYickBE7vLk3zOeOTwVca+WCdwxaSlfr9zu65CMMeaEOVTsnh/gb5+uoH2DYKbc3J+2jWxQYFU8KQRuwX2lwCRgMnAIuPlIjUQkGHgJGAJ0AEaKSIcKm10D7FXVNOB54KkKrz8PfOVBjOYoRYYFM+6q3nRNieOWD5Yw45ccX4dkjDHHbeOuAoa9PIdPlmRx66np3NkznLgoGxRYHU9uOlSgqvfhHiDYS1X/pqoFHrx3BpCpqhtUtRiYCAytsM1QYLzz/CPgVHGuaxORC4ANwCoPczFHKSY8hLdHZdC+cSw3vruYH9bt9HVIxhhzzL5euYPzx/7IjvxC3rqqN3ee3oYgu1T6iDyZWfBkEVkNrHaWu4rIyx68d1Nga7nlLGddpduoaimQBzQQkWjgXuBRD/4dcxxiI0J5Z1QGrZKiuW7CQub+utvXIRljzFEpLXPxf1N/4YZ3F9EqKZovb+nPwLYNfR1WwBDV6meaE5F5wHBgiqp2d9atVNVOR2h3EXCmql7rLF8OZKjqLeW2WeVsk+Us/4r7SML9wHxVnSwijwAHVPWflfwb1wHXASQlJfWcPHmyZ1kHoAMHDhATE+O1988vVp6cf4hdh5Tbe0TQoUHNTrXp7fx8zfILXLU5Nwj8/PYVunh5WRHr9roY3DyEke3CCC13aWCg53ckgwYNWqSqvY7nPTwa/a+qWyvMROfJHWyygGblllOAikPUD2+T5UxdHAfsAfoAw0XkaSAecIlIoar+u0Jcr+HMcti2bVsdOHCgJ+kEpFmzZuHt/Pr2LeKyN+bxryUFvH5FLwa0SfLqv1deTeTnS5Zf4KrNuUFg5/fzht384/0lFBQJL4zoxgXdKx50Duz8aoongwW3isjJgIpImDOC/xcP2i0A0kUkVUTCgEuAKRW2mQJc6TwfDsxUt1NUtaWqtgReAP6vYhFgTrykeuG8P7oPrZJiuPadhXy3JtfXIRljzB+4XMpL32Vy6RvziI0M4fMx/SotAoxnPCkEbsB9lUBT3N/gu+HBVQPOOf8xwDe4C4fJqrpKRB4TkfOdzd7EPSYgE7gT+MMlhqZmNYgJ54PRfWiTHMN1ExYyfbVdTWCM8R+5+wu5Ytx8nvlmLUM6NeLzm/vRxu4XcFyOeGpAVXcBlx7Lm6vqVGBqhXUPlXteCFx0hPd45Fj+bXPs4qPCeO/ak7hi3HxufHcRY0d2t4k4jDE+N3v9Tu6YtJQDRaU8+efOjOjdzG6gdgJUd/fBsUCVIwlV9VavRGT8QlxkKBOuyeCqcfMZ88ESnncp53dt4uuwjDF1UEmZi+enr+OV738lLSmG90efZEcBTqDqjggsrLEojF+KjQjlnWv6MOqtBdw2cQkFRaWMzGju67CMMXVI1t6D3DZxKYs27+WS3s14+LyORIbV7FVNtV11hcBpqnq5iNymqv+qsYiMX4kJD2H8qAxufG8R93+ygn0HS7hxYGtfh2WMqQO+WbWDuz9chkvhxZHd7aikl1Q3WLCniLQARolIfRFJKP+oqQCN70WGBfPa5b04v2sTnvp6DU989QtHmn/CGGOOVWFJGQ99vpLrJyyiRQP3BEFWBHhPdUcEXgW+BlrhvslQ+REZ6qw3dURYSBAvjOhGXGQo//l+A3kHS3h8WGe7p7cx5oRalZ3HbROXkpl7gGv6p3LPWW0JD7FTAd5UZSGgqi8CL4rIK6p6Yw3GZPxUUJDw2NCOxEeFMnZmJvmFJTw/opvtpMaY41bmUt6YvYF/TltL/agwJlyTwSnpNTepWV3myeWDN4pIfyBdVd8SkUSgnqpu9H54xt+ICH89oy3xUWH8/cvV5B1awCuX9SQ2wu7uZYw5Ntv2HeKvk5fy84Y9nNWxEU/8uTP1o8N8HVadccRCQEQeBnoBbYG3gDDgXaCfd0Mz/uya/qnER4Zy78fLufjVubx1dW8ax0X6OixjTICZsiybBz5dgculPDO8C8N7ptjcADXMk5kFhwHnAwUAqpoN2AWchgt7pvDW1b3J2nuIYS/9xJod+b4OyRgTIPILS7h94hJu/WAJ6Q1jmHrbKVzUyyYI8gVPCoFidQ8RVwDnFsHGAHBKehKTr++Lolz0ylzmZO7ydUjGGD83J3MXQ16YzRfLt3PHaW2YfH1fWjSwPy2+4kkhMFlE/gPEi8ho4FvgDe+GZQJJhyaxfHpTP5rER3LluPl8sjjL1yEZY/xQQVEpD362kkvfmEd4SBAf3tCX205LJyTYkz9Fxls8GSz4TxE5HcjHPU7gIVWd7vXITEBpEh/Jhzf25YYJi7hz8jKy9h7ilsFpdpjPGAO4bxl890fuz4Zr+6dy15ltiQi1K478wRELAQDnD/90ABEJFpFLVfU9r0ZmAk5sRChvX53BfR8v57np68jMPcDTw7vYzm5MHXawuJSnv17L2z9tokWDKCZf35feLW1OOn9S3U2HYvnf7Yen4C4EbgbuBpYCVgiYPwgLCeLZi7uSlhzDM9+sZfPuAl67ohfJsRG+Ds0YU8MWbtrDXR8uY9Pug1x1ckvuOastUWEeff80Nai6HpkA7AXmAtfiLgDCgKGqurQGYjMBSkS4aWAaaUkx3D5pKef/+0dev6IXXVLifR2aMaYGHCou49lpa3lzzkZS6kfyweiT6Nu6ga/DMlWoboRGK1W9SlX/A4zEPZfAuVYEGE+d0bERH994MiFBQVz06ly+WJbt65CMMV724/pdnPnCD7zx40b+ktGcr28bYEWAn6uuECg5/ERVy4CNqrrf+yGZ2qR941imjOlHl5Q4bvlgCc9OW0uZy25YZExts7egmLs+XMZlb84jOEiYeN1JPD6sM9HhdirA31XXQ11F5PAMMQJEOssCqKrGej06Uys0iAnn3Wv78OBnKxk7M5PlWXm8MKKbTSFqTC2gqnyxfDuPTllF3qESbh7UmlsGp9sg4QBS3U2HrBfNCRMeEsxTF3aha7N4Hp2ymnPH/sirl/Wkc0qcr0MzxhyjbfsO8eBnK5m5JpeuKXG8e20f2je274iBxmZxMDVGRLi0Twsm39AXVeXCV39i8oKtvg7LGHOUylzK23M2csZz3zP31908eG4HPrmpnxUBAcpO3pga161ZPF/eegq3frCEez5ezpKtexkUZ+MGjAkES7bs5cHPV7JyWz4D2iTx+AWdaJYQ5euwzHGwQsD4REJ0GONHZfDc9LW89N2v/BQbRJuuBbRMtPnGjfFHewuKefqbtUxcsIWkmHDGjuzOuV0a2+yhtYCdGjA+Exwk3H1mO16/ohe7Drk458XZfL50m6/DMsaU43IpkxdsZfCzs5i8cCuj+qUy469/4ryuTawIqCXsiIDxudM7JPPYyZFM3BzBbROXMnv9Lh49v6NddmSMj63OzufBz1eyaPNeeresz98v6ES7RjYOoLaxT1rjFxpEBvHB6JN4ccZ6xn6XyeItexk7sjsdm9hVBcbUtLyDJTz/7TrembuJ+lFh/POirlzYo6kdAailvHpqQETOEpG1IpIpIvdV8nq4iExyXp8nIi2d9aeLyCIRWeH8HOzNOI1/CAkO4s4z2vLetX0oKCpl2Es/8facjajaQEJjakJpmYsJczcx8J/f8c7cTfylT3Nm/nUgw3umWBFQi3ntiICIBAMvAacDWcACEZmiqqvLbXYNsFdV00TkEuApYASwCzhPVbNFpBPwDe6bH5k64OTWiUy99RTu/mg5j3yxmhlrcnl6eBcax0X6OjRjaq3Z63fy9y9Xsy7nAH1bNeCh8zrY5YB1hDePCGQAmaq6QVWLgYnA0ArbDAXGO88/Ak4VEVHVJap6eGL6VUCEiIR7MVbjZxrEhPPmlb14fFgnFm7ayxnP/8Ani7Ps6IAxJ9iOAhfXjl/A5W/Op7DExauX9eT90TYxUF3izTECTYHys8VkAX2q2kZVS0UkD2iA+4jAYRcCS1S1yIuxGj90eAKifq0T+euHy7hz8jKmrcrh8WGdaBBjdaExxyPvUAn/nrmecT8eIjKshPuGtOPqfi0JD7FJZesa8dY3LBG5CDhTVa91li8HMlT1lnLbrHK2yXKWf3W22e0sdwSmAGeo6q+V/BvXAdcBJCUl9Zw8ebJXcvEHBw4cICYmxtdheM2R8nOp8vXGEj5ZX0JUKFzVMZweyYEz1rWu918gq225lbiUGZtL+WJDMQdL4KRkZUSHKOLDa+fV5LWt/yoaNGjQIlXtdTzv4c1P0iygWbnlFKDifWgPb5MlIiFAHLAHQERSgE+BKyorAgBU9TXgNYC2bdvqwIEDT2T8fmXWrFnU9fwGD4Krt+dzx6SlvLhkP+d0SeCR8zqSVM//jw5Y/wWu2pKby6V8tnQbz05bx7Z9xQxok8S9Z7Vl57oltSK/qtSW/vMmb5aAC4B0EUkVkTDgEtzf7subAlzpPB8OzFRVFZF44L/A/ao6x4sxmgDjvq1xf+48vQ3TV+Vw2nPfM3nhVhs7YEwVVJVZa3M5Z+yP3Dl5GfWjQ3nv2j68MyrDLs81gBePCDjn/MfgHvEfDIxT1VUi8hiwUFWnAG8CE0QkE/eRgEuc5mOANOBBEXnQWXeGquZ6K14TOMJCgrj11HTO7tyI+z9ZwT0fLefzpdv4v2GdadHApig25rDlWft4Yuoa5m7YTfOEKF4c2Z1zOzcmKMguBTT/49WTrKo6FZhaYd1D5Z4XAhdV0u4fwD+8GZsJfGkN6zHpur68P38LT361hjNf+IHbT2vDqH6phIXUzvOdxnhidXY+L3y7jmmrc0iIDuOR8zrwlz4tbL8wlQqc0VbGVCIoSLjspBac1j6ZBz9fyZNfreHDhVt59PxO9E9P9HV4xtSotTv2868Z65i6Ygf1IkK447Q2jOrfknoRob4OzfgxKwRMrdAoLoLXr+jFjF9yePSL1Vz25jzO7tyIB87pQNN4m4jI1G6Zufv514xMvlyeTXRYCLeems41/VOJi7QCwByZFQKmVjm1fTL90hJ5/YcNvDQrk5lrchkzKI3RA1rZ9dGm1vl15wH+PTOTz5duIyI0mJsGtmb0Ka2IjwrzdWgmgFghYGqdiNBgbjk1nWE9mvL4f3/hn9PW8eGiLO49qx1DOjWyOdNNwFuRlcfLszL5etUOIkKCGT2gFdcPaE1CtBUA5uhZIWBqrZT6UbxyWU9mr9/JP778hZveW0yP5vE8cE57erZI8HV4xhwVVeXnDXt4eVYms9fvol5ECGMGpXHVyS1tpk1zXKwQMLXeKelJTL0tkY8WbeXZaeu48JW5nNWxEfcOaUdqol1uaPyby6XMXJPLS7MyWbJlH4kx4dw3pB2X9mlugwDNCWGFgKkTgoOEEb2bc17XJrwxeyP/+f5Xvv0lh7/0ac7Ng9JIjo3wdYjG/E5hSRmfLdnGuDkbWZdzgJT6kfz9gk5c1DOFiFAb72JOHCsETJ0S5YyoHpnRnBe+Xcf787YwacFWLu3TghsGtqJhPSsIjG/l7i/k3bmbeXfeFvYUFNO+cSzPj+jKeV2aEBJs8wCYE88KAVMnJdUL5/FhnbnhT60ZO3M94+du4v35m7mib0uuH9DKzrmaGrcqO49xP27ii2XZlLhcnNoumWv6p3JSqwQb4Gq8ygoBU6c1S4ji6eFduXFgGmNnrOeN2Rt49+fNXH5SC0b1T7VTBsariktdTFu9gwlzNzNv4x6iwoIZmdGMq/ql2vgVU2OsEDAGSE2M5rkR3bhpUBovzljP67M38NacTfy5R1OuG9CKVkm19zampuZt3XOQiQu2MGlBFrsOFJFSP5L7h7TjkozmNgmQqXFWCBhTTlrDGF4c2Z2/ntGG137YwIeLspi0cCtDOjXihj+1pktKvK9DNAGqzOW+C+B787bw3dpcBBjcLplLT2rOn9KT7EZAxmesEDCmEi0aRPP4sM7cflob3pqzkQk/b2bqih1ktEzgypNbcmbHZBu4ZTyyaVcBHy/O4uNFWWTnFdKwXji3DErjkozmNLHpr40fsELAmGok1QvnnrPacePA1kycv5Xxczdx8/uLaRwXwWUnteCS3s1sYKH5g/zCEqYu385Hi7JYuHkvQQL905N48NwOnNYhmVArIo0fsULAGA/Uiwhl9IBWjOqfysw1uYz/aRPPfLOWf81Yz3ldmjAyoxk9W9S30d11WGmZi59+3c3Hi7P4ZtUOCktctE6K5t6z2jGse1MaxdnAU+OfrBAw5igEBwmnd0jm9A7JZObuZ/xPm92HfRdn0Topmot7NePPPVJIqmdHCeqCMpcyf+Mevlyezdcrd7C7oJjYiBCG90zhwh4pdGsWb8Wh8XtWCBhzjNIa1uPvF3TiviHt+O/y7UxauJUnvlrDM9+sZXC7hlzcqxkD2iQRFmKHgWsTl0tZsnUvXyzbztQV28ndX0REaBCntk/mvC6NGdi2oc38ZwKKFQLGHKfo8BAu7t2Mi3s3IzN3P5MXZvHJ4iymrc4hPiqUIZ0a0VzLOMWlBNvI8IBUXOri5w27mb46h29/yWF7XiFhIUEMbJPEeV2bcGr7hkSF2cepCUz2m2vMCZTWsB5/O7s9d5/ZltnrdzJlaTafL83mYHEZb62ZwbldmnBOl8Z0bxZvl4v5ubyDJXy3Npf3lxYy5rvpHCgqJSI0iAHpSdx9ZltO75BsN/0xtYIVAsZ4QWhwEIPbJTO4XTKHissY+/F3ZJbE8+7Pmxk3ZyOJMeGc1r4hp3dIpl9aoh1K9gNlLmXltjxmr9/JD+t3sWjzXspcSmyYcG7XFOsrU2tZIWCMl0WGBZPROIR7BvYi71AJs9bmMm11Dl8u387EBVuJDA3mlPRETm3fkP7pSTS1a8trTPa+Q7/94Z+TuYt9B0sA6NgklusHtOK0Dsns+3Upgwd18XGkxniPFQLG1KC4yFCGdmvK0G5N/3DeedrqHABaJUbTLy2R/umJnNSqgU05e4KoKht2FbBg4x7mb9rDgk172LrnEADJseGc1j6ZU9IT6ZeWSGK5uSFmbbBTOKZ2s0LAGB8JCwliQJskBrRJ4rGhHVmfe4DZzjfTjxdnMeHnzQQJdG4aR88WCfRqWZ+eLerbjZA8dLC4lNXZ+Szduo+Fm/aycPMedh0oBiAhOozeLetz1cmpnJKeSHrDGLvMz9RZVggY4wdEhDbJ9WiTXI9r+qdSXOpi6dZ9/Lh+J/M27uH9+e6xBQAp9SPp2aI+3ZvF07FpHO0bxxITXrd35UPFZazN2c+KrH0sy8pjRVYe63P341L36yn1IxmQnkTv1AR6t0ygdVK0/eE3xuHVTw8ROQv4FxAMvKGqT1Z4PRx4B+gJ7AZGqOom57X7gWuAMuBWVf3Gm7Ea40/CQoLISE0gIzUBcF++tnp7Pos272XR5j3M/XU3ny/NBkAEWjaIpkPjWDo0iaVD41haJUXTND6y1t0Poai0jI27Cli7Yz/rcw6wNmc/63L2s2XPQdT5o58QHUaXlDjO7NSILk3j6JISR0M7imJMlbxWCIhIMPAScDqQBSwQkSmqurrcZtcAe1U1TUQuAZ4CRohIB+ASoCPQBPhWRNqoapm34jXGn4WFBNGtWTzdmsVzTf9UVJWc/CJWZeexOjufVdn5LN+2j/+u2P6/NsFBNG8QRavEaFolxZCaGEWT+Egax0XSJD7CL697L3MpuwuK2JFXyObdB9my5yBbDv/cc5DteYd++5YfHCSkJkbTsUksw7o3pV2jenROiadJXIR92zfmKHjzkyADyFTVDQAiMhEYCpQvBIYCjzjPPwL+Le49eCgwUVWLgI0ikum831wvxmtMwBARGsVF0CguglPbJ/+2Pu9QCetz9rNhVwEbdhawYecBNuwq4Lu1uZSU6e/eIz4qlCZxkTSKiyAhOoyE6DDqR4WREB1K/agw4qPCiAoLJjIsmKiwYKJCQ4gICyIsOKjaP7Qul1JU6qK41EVRaRlFpS72F5aSd6iE/MIS90/nsbugmJz8InL3F5KbX8TOA0WUuX4fZ2JMOM0TIslITaBZQhRpDWNokxxDamI04SF2KZ8xx8ubhUBTYGu55SygT1XbqGqpiOQBDZz1P1do29R7oRpTO8RFhtKrZQK9Wib8bn1pmYvteYVk7zvE9rxCtu07xPa8Q2zfV8iO/ELWbM9nd0ExRaWuI/4bIhAsgoi7IAkSUJeL4JlfU1Lm+kPBUd37xEeGkhwbQcPYCNom16NhbDjJsREkx0bQPCGK5glRRNfx8Q/GeJs397DKvjJU/ISoahtP2iIi1wHXASQlJTFr1qyjDDFwHDhwwPILYP6UXzwQL9Ax3lkA3MN4IikqUw4UK/uLlYISKCpTisp+/7PEBarOA/ejqLiM0FAhNCiEkCAIDYLQIPnteUSIEB0qRIdCVIgQFSpEhkCQCOACDjoPoAjYCTk7IaeG/28q40995w2Wn/FmIZAFNCu3nAJkV7FNloiEAHHAHg/boqqvAa8BtG3bVgcOHHiiYvc7s2bNwvILXJZf4KrNuYHlZ8CbQ4oXAOkikioiYbgH/02psM0U4Ern+XBgpqqqs/4SEQkXkVQgHZjvxViNMcaYOslrRwScc/5jgG9wH3ccp6qrROQxYKGqTgHeBCY4gwH34C4WcLabjHtgYSlws10xYIwxxpx4Xh2Fo6pTgakV1j1U7nkhcFEVbR8HHvdmfMYYY0xdV7tmGzHGGGPMUbFCwBhjjKnDrBAwxhhj6jArBIwxxpg6TFQ9mwXM34nIfmCtr+PwokRgl6+D8CLLL7DV5vxqc25g+QW6tqpa73jeoDbN3blWVXv5OghvEZGFll/gsvwCV23ODSy/QCciC4/3PezUgDHGGFOHWSFgjDHG1GG1qRB4zdcBeJnlF9gsv8BVm3MDyy/QHXd+tWawoDHGGGOOXm06ImCMMcaYoxQwhYCIBIvIEhH50llOFZF5IrJeRCY5dzisrN39IpIpImtF5MyajdozleT2nhPvShEZJyKhVbQrE5GlzqPinR39RiX5vS0iG8vF3q2Kdlc6/bteRK6sbBt/UEl+s8vlli0in1XRzu/7T0Q2icgKJ8aFzroEEZnu9Mt0EalfRVu/778q8ntGRNaIyHIR+VRE4j1t62+qyO8REdlW7nfv7CranuV8DmWKyH01G7lnqshvUrncNonIUk/b+hMRiReRj5zfxV9EpK/X9j1VDYgHcCfwPvClszwZuMR5/ipwYyVtOgDLgHAgFfgVCPZ1Lh7kdjYgzuODynJztjvg69iPMb+3geFHaJMAbHB+1nee1/d1Lp7kV+G1j4ErArX/gE1AYoV1TwP3Oc/vA54K1P6rIr8zgBDn+VOV5VdVW397VJHfI8BdR2gX7HxetgLCnM/RDr7O52j7AHgWeCgQ+w8YD1zrPA8D4r217wXEEQERSQHOAd5wlgUYDHzkbDIeuKCSpkOBiapapKobgUwgw/sRe65ibuC+a6M6gPlAiq/iO16V5eehM4HpqrpHVfcC04GzTnR8x6u6/ESkHu7f00qPCASwobj3Oah63wuI/quMqk5T1VJn8WcCeP87DhlApqpuUNViYCLufg8Yzt+Ji3F/mQooIhILDADeBFDVYlXdh5f2vYAoBIAXgHsAl7PcANhXbmfNAppW0q4psLXcclXb+VLF3H7jnBK4HPi6irYRIrJQRH4Wkcp+IfxBVfk97hx9O7vIAAAE8UlEQVR6fV5EwitpFwh9B9X0HzAMmKGq+VW0DYT+U2CaiCwSkeucdcmquh3A+dmwknaB0n+V5VfeKOCrY2zrD6qKcYyz/42r4vBybei/U4AcVV1/DG19rRWwE3jLOe34hohE46V9z+8LARE5F8hV1UXlV1eyaWWXP3i6nU9UkVt5LwM/qOrsKl5vru4Zs/4CvCAirb0R57GqJr/7gXZAb9yHr+6trHkl6/ym78Cj/htJ9d9G/Lr/HP1UtQcwBLhZRAZ42M7v+89RZX4i8gBQCrx3tG39SGUxvgK0BroB23EfPq8o4PuPI+9//tx/IUAP4BVV7Q4U4D4V4Imj7ju/LwSAfsD5IrIJ9+Gpwbi/hcWLyOEpklOA7EraZgHNyi1XtZ2v/CE3EXkXQEQeBpJwn3+ulKpmOz83ALOA7l6O92hVmp+qbnfOfBQBb1H56Rp/7zuovv8a4M7rv1U1DoD+Kx9jLvAp7pxyRKQxgPMzt5KmgdB/VeWHM8DqXOBS5xSdx239SWUxqmqOqpapqgt4ncDd/6rrvxDgz8Cko23rJ7KALFWd5yx/hLsw8M6+5+sBEUc5eGIg/xtw9iG/Hyx4UyXbd+T3gwU34IeDBSvJ7VrgJyCymu3rA+HO80RgPX44mKeK/Bo7PwV3UfdkJdsnABudPOs7zxN8nYcn+TnLNwDjA7n/gGigXrnnP+E+1/gMvx+w9HQg9l81+Z0FrAaSjratr3PyML/G5ba5A/c4qoptQ5zPy1T+N1iwo69z8rQPnDy/D/D+m437hkLgHuD5jLf2PZ8ne5T/MeX/mLTCPZAuE3dRcPhD9XzgsXJtHsA9+nUtMMTXOXiYW6kT81Ln8ZCzvhfwhvP8ZGCFs4OuAK7xdQ5Hkd9MJ+aVwLtATMX8nOVRTv9mAlf7OgdP83OWZ1X8YAm0/nP2sWXOYxXwgLO+ATADd/Ey4/CHTKD1XzX5ZeI+x3p4/3vVWd8EmFpdW396VJPfBOd3bjkwhf8V5r/l5yyfDaxzPosCJj/ntbeBGypsH2j91w1Y6PTTZ7j/qHtl37OZBY0xxpg6LBDGCBhjjDHGS6wQMMYYY+owKwSMMcaYOswKAWOMMaYOs0LAGGOMqcOsEDDGGGPqMCsEjDHGmDrMCgFjzB+IyLBy93Q//HCJyBAROeBB+5YiMtW5n/06Ebm/JuI2xhw9m1DIGHNEzt3ZLgUGAfmqGlPNtkHAAuBRVZ3i3F3yfeAbVX2tRgI2xnjMCgFjTLVEpA3uaaFPVtUtInLgCIXAEOBaVb2w3LrGuOd+b+P9iI0xR8NODRhjqiQiobi/zd+lqls8bNYe9xzuv1H3vdNjRSTsBIdojDlOVggYY6rzd2CVqk48ijZlwO+OGIiIAFG4b6hljPEjVggYYyolIgOBC4ExR9l0FnC288f/sNOBxarqOjHRGWNOFCsEjDF/ICL1gbeAK1R1fzXbzRCRpuXXqeoyYAnwmLNNMvAc8DfvRWyMOVYhvg7AGOOXbgAaAq/8/os9Txx+4lwdkAbsKb+BiNyH+/7ol4nId7iPKLQAXhaRC1R1k3dDN8YcDbtqwBhzTESkEzBKVe/0dSzGmGNnhYAxxhhTh9kYAWOMMaYOs0LAGGOMqcOsEDDGGGPqMCsEjDHGmDrMCgFjjDGmDrNCwBhjjKnDrBAwxhhj6jArBIwxxpg67P8DEN6C1W57gYkAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 576x576 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"Z0 = 50\n",
"Zl = np.arange(1, 150, 0.01)\n",
"gamma = (Zl-Z0)/(Zl+Z0)\n",
"\n",
"plt.figure(figsize=(8,8))\n",
"plt.subplot(211)\n",
"plt.plot(Zl, -20*np.log10(1-abs(gamma)**2))\n",
"plt.grid()\n",
"plt.ylim(0, 10)\n",
"plt.xlim(0, 150)\n",
"plt.xlabel('Zl, $\\Omega$')\n",
"plt.ylabel('Reflection loss, dB')\n",
"\n",
"Zl = np.arange(40, 60, 0.01)\n",
"gamma = (Zl-Z0)/(Zl+Z0)\n",
"plt.subplot(212)\n",
"plt.plot(Zl, -20*np.log10(1-abs(gamma)**2))\n",
"plt.grid()\n",
"#plt.ylim(0, 10)\n",
"plt.xlim(40, 60)\n",
"plt.xlabel('Zl, $\\Omega$')\n",
"plt.ylabel('Reflection loss, dB')\n",
"plt.savefig('reflection_loss.png')"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.6.8"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Display the source blob
Display the rendered blob
Raw
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment