{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Action Potentials in Neurons" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "from ecell4.prelude import *" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Hodgkin-Huxley Model" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/markdown": [ "HODGKIN AL, HUXLEY AF, A quantitative description of membrane current and its application to conduction and excitation in nerve. *The Journal of physiology*, **4**(117), 500-44, 1952. [10.1113/jphysiol.1952.sp004764](https://doi.org/10.1113/jphysiol.1952.sp004764). PubMed PMID: [12991237](https://www.ncbi.nlm.nih.gov/pubmed/12991237)." ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "citation(12991237)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "Q10 = 3.0\n", "GNa = 120.0 # mS/cm^2\n", "GK = 36.0 # mS/cm^2\n", "gL = 0.3 # mS/cm^2\n", "EL = -64.387 # mV\n", "ENa = 40.0 # mV\n", "EK = -87.0 # mV\n", "Cm = 1.0 # uF/cm^2\n", "\n", "T = 6.3 # degrees C\n", "Iext = 10.0 # nA\n", "\n", "with reaction_rules():\n", " Q = Q10 ** ((T - 6.3) / 10)\n", "\n", " alpha_m = -0.1 * (Vm + 50) / (exp(-(Vm + 50) / 10) - 1)\n", " beta_m = 4 * exp(-(Vm + 75) / 18)\n", " ~m > m | Q * (alpha_m * (1 - m) - beta_m * m)\n", "\n", " alpha_h = 0.07 * exp(-(Vm + 75) / 20)\n", " beta_h = 1.0 / (exp(-(Vm + 45) / 10) + 1)\n", " ~h > h | Q * (alpha_h * (1 - h) - beta_h * h)\n", "\n", " alpha_n = -0.01 * (Vm + 65) / (exp(-(Vm + 65) / 10) - 1)\n", " beta_n = 0.125 * exp(-(Vm + 75) / 80) \n", " ~n > n | Q * (alpha_n * (1 - n) - beta_n * n)\n", "\n", " gNa = (m ** 3) * h * GNa\n", " INa = gNa * (Vm - ENa)\n", " gK = (n ** 4) * GK\n", " IK = gK * (Vm - EK)\n", " IL = gL * (Vm - EL)\n", " ~Vm > Vm | (Iext - (IL + INa + IK)) / Cm\n", "\n", "hhm = get_model()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "Vm > m + Vm | (1.0 * ((((-0.1 * (Vm + 50)) / (exp((-(Vm + 50) / 10)) - 1)) * (1 - m)) - (4 * exp((-(Vm + 75) / 18)) * m)))\n", "Vm > h + Vm | (1.0 * ((0.07 * exp((-(Vm + 75) / 20)) * (1 - h)) - ((1.0 / (exp((-(Vm + 45) / 10)) + 1)) * h)))\n", "Vm > n + Vm | (1.0 * ((((-0.01 * (Vm + 65)) / (exp((-(Vm + 65) / 10)) - 1)) * (1 - n)) - (0.125 * exp((-(Vm + 75) / 80)) * n)))\n", "m + h + n > Vm + m + h + n | ((10.0 - ((0.3 * (Vm - -64.387)) + (pow(m, 3) * h * 120.0 * (Vm - 40.0)) + (pow(n, 4) * 36.0 * (Vm - -87.0)))) / 1.0)\n" ] } ], "source": [ "show(hhm)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEGCAYAAACO8lkDAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO2dd5xcZbn4v8+ULQkphCSUhE4ooSeh2S4iFrwqiqBIFQsiWFDvVdF7LZfrzwJXhSsiSC8KXEBBVFQQUVRIA0ISSgKkl02ySXazbWbOeX9/nHNmzs6cMzO7O+V9N+/388lnd2ZnM8++c877vE8XpRQWi8VisYRJNFsAi8ViseiHVQ4Wi8ViKcEqB4vFYrGUYJWDxWKxWEqwysFisVgsJaSaLUAtmDx5stpvv/2aLYbFYrEYxYIFCzYrpaZE/WxUKIf99tuP+fPnN1sMi8ViMQoRWRn3M+tWslgsFksJVjlYLBaLpQSrHCwWi8VSglUOFovFYinBKgeLxWKxlGCVg8VisVhKaJpyEJG9ReQJEXlRRJaIyOf95yeJyJ9EZJn/dddmyWixWCw7K820HHLAl5RShwEnApeJyEzgq8DjSqkZwOP+Y4tFO17f3MMLa7Y3WwyLZvx20Xo6ezLNFmPENE05KKXWK6UW+t93Ay8C04DTgdv9l90OvL85EloazVPLNtOfdZotRtW89eq/8N6fPNVsMapmzdZeZnz9d0YptOueWM5+X/0tpsyd6ezJcNkvFvLJO8wvytUi5iAi+wHHAs8Auyul1oOnQICpMb9zsYjMF5H5mzZtapSoljrx8oZuzrv5Ga58ZGmzRRm1LFy1jayj+PFjrzRblKq56g8vA9CTMePQkMm5ALy4vqvJkoycpisHEdkFeAC4XClV9YoqpW5USs1RSs2ZMiWyNYjFIDZ09QOwYktPkyUZOjnHbbYIVTGhPQ3Auu39TZZk6HT3Z5stQlX0+ZZvryHKrBxNVQ4iksZTDHcrpR70n94oInv6P98T6GiWfJbG0dXn3fzt6WSTJRk6W3vN2Liy/qnWFGUWpqsv12wRqqJvFCiFgGZmKwlwM/CiUuqHoR89DFzof38h8FCjZbM0nm29XgAvlWi6MTtkugw51eZcTyk4rhn++zCmWQ6jgWZ2ZX0jcD7wgog85z/3NeB7wH0i8nFgFXBWk+SzNJCcv2H1GnhzmRJEzzjeGueMVA5mWA6mXAvV0DTloJR6CpCYH7+tkbJYmk+wX/UOmLEJhOnPmuGmMdmtZIoPfzQpB/NseMuoxPW1Q8bAjcuUDSHrr61JlkM66Z0fTVnj0eRWssrBogWun8cepAKahCkbl4nKoSXpbVH9OTPW2MTrNw6rHCxa4CizLIdwUZYpp8Ug5pA1ZI0B0ilfORjiujMx2B+HVQ4WLQj2WlM2rvAeYMrGFcQaTNrA8paDIQo4XMht0jpHYZWDRQuCG8kUszx845tiOZjoVkolvJjDgCFr7Ia0gykKLQ6rHEYxvRlzMn+CmyrrmLFxhTcBUzausFvJlF5FgR7rN+TQENa7phwa4rDKYZQyb0UnM7/xB/62zIy+U8FNZcpGG95bTamKDSwHpWDAmM3WW2hTTuGOtRwsuvP86m0APP6iGd1HglRWU06IgzYBQzJpwq4wUzYu05SDssrBojvj27wma0HPIt0JNgHHVUYEpcNupb6M/vJCQQGDOUVlebeSIUH/8Bqbcl3EYZXDKCXhB/K6Dak4Ns0cdw0MSDsGpt8G1o4J1wQUZbEZYlHGYZXDKKXPoGA0FPnwDdgIwpuAjZPUj8BCM+GagGKL0gyZ47DKYZQSuA0Scd2rNCPsDx8wwIXgGngKH5RhZUhsJ7DQTLgmwMzrIg6rHEYppl2YpuWHuwYHdwEGDHF5FFJZzZIXzLku4rDKYZQSmLSmFJUN3mz1l9nEfPbBrjD91xgKcRJTNtrB9S9mrHEcVjmMUoJiMhM2WjAvkOcY6KJRBloOSpl1HYfjOqascRxWOYxSTAvkmZytZMoJ0Q2Jacpma1q20qDYmSGHhjischilmHZTDS4e0v+mMvGE6BhoORTqHEyR16xDTjmschilBBuBKaeXwScu/W+qYH0TYoYyA2/jGtOSBMy4LsIHBlNauQeXsYgZa1wOqxxGKYHbw5Rca9OCpcEJcUxLyphNQCnyysGEU61TlKRgQrNApRQJgbZU0pjrIg6rHEYpebeSAadw8JRZ0LvfhJsq2Kja0kkjLB3wFFpb2rccjFDA3tdAoZlgPTiuIiFCazphhAIuh1UOo5TA7WGO5aBoz7s89Jc52KfGtJhzQnSVNzwnlRAjDg0F68wcV5irvNY1ramEEQq4HFY5jFICC3wgZ4Y57ihoTwcuD/1vqvDGlcm5g7KXdMVVChGM2biCNc4fGgyQOXArtabMsSjjsMphlGJaSp15lsPgjcsEl4e3cQlt6aQRlkN+jdPmxEnc/BonjLjvymGVwyhlUAdOA1xLrqtIJYR0Uoy4qVSRP9yEU63r4vnDjbEcvK/tLSnAjEOOk1/jpBHKrBxWOYxSBrWjMOCU6CpFMiFelocRG1fRqdaANXYCt1LajDhJcA2PSZtjUbp5t5K1HEY9OwZyRpy8iwm7lUyQ33FB/CwPEzYBJ5StBGZYDspXwK0pMzJpTAxIK6VIJMTPYtNf3nJY5VCBo7/9R95y1RPNFmPIDK7U1P8i9TYufHPcDHkhvHGZsNn6Lg9DNq68AjbJdadCrjsDrolyWOVQAcdVbOoeaLYYQyZsOZgQLHX8QJ4pN1UhldUcf7ibL9Ayw3LIx3UMc90lBL/OQf9rohzaKgcReZeIvCwiy0Xkq82WJ2fABhvGMayvvKs8t1KLIb5at8itZNIaG2M5uGamsoofkDbhkFMOLZWDiCSB64DTgJnAR0RkZqPlGNQMzoCbKYxrWiqrq0gKxvhqjfWHG2Q5lNQ5GLDZui4kg1RWA5RZObRUDsDxwHKl1GtKqQxwD3B6o4UI3/C9hs1k9lol+O0oDNkICmmWBsgbqpAGQzauYI3TSSOGQOXXOG2e686mstaPacDq0OM1/nMNZcdAQSH0Z/S/MMM4rjLKH+64XpZHazpphJVmYvVuvu+PcZaDSYcc33VniHu0HLoqB4l4blB/AhG5WETmi8j8TZs21UWIcApob9Y8y8EslweF/HAjNgHz6hy8jQtjqncLysGcQ46rFIlEoZbEhNY1ceiqHNYAe4ceTwfWhV+glLpRKTVHKTVnypQpdREiGwpCm1ArEMazHMxyeQQ5+Ea4PIpjDgZYDoU6B8PiOmlzDjmuUiR9ywHMyBSMQ1flMA+YISL7i0gLcDbwcKOFyIULyQw4zYZxVOjEZcDGFaSyGhOQ9kU061RLvu+PGW4l72trOuENzzFE5kRIOZiczppqtgBRKKVyIvIZ4A9AErhFKbWk0XKYbDm4rmKXVvM2LmPqHIrdSkZsXEFX1iQ5V5FzXFJJXc+HhVTWpEE+/PwaD2r5kW6uUMNES+UAoJT6HfC7ZsqQDRULZA0zDx3XrC6nrhvO8tB/rUsrpPWXOWw5gCezzsohcCslEuY0snP9oH9bKgii639dxKHvlaEB4cK3jGNWYMlVXpfTlqQZlZr5mIMhvZUCl0c6mSCZECNkDs8aAP0VWuC6SxhmOST9rDvQf43LYZVDGcLBpKxhH7LjFgK8JmxcjqvyKYBZRw1q/6EjgXxBUZkJJ8RwLQno7woLLIdkAv/QYMIaF1JZQf81LodVDmXImexWUkHdgBk3lVJBZak/PEdzmQe5PIwZnuNtXG2GnGqDuI6I38rdgDXOV6EbssblsMqhDDk37FYy60NWqjB0xIRTrRPkhwe+Ws03grxyMGh4TqHzrRlrHMR1kn4rdxPco+FCQ9B/jcthlUMZMrlQZ1PDTgCO36vIFLdSweVhxhzpgj/crH5QiZDloPsaO4NiDmZYDq4KAug2ID2qCVsOWcMC0kE7CmO6nBp24iq2HEzwLRfn4OteN1Bw3WGMdeaWBP31XuNyWOVQhnCcwbSYQ1Cpac6plny2Eujvqx2cZmmIAi7Kwde9h5XrFhSwKdexikgXNpXYOgcR+Q1F/YzCKKXeVxeJNCJsLRjpVgo2Ls1PiBBkK0Fbyox2FEEyVdIgl4cyznLwvpqWdZcIF8Fpfh2Xo1wR3NUNk0JTTM5WclWhy2lXX7bZ4lREBT1p0qa5lbw0y3AHX13JT4Lz11h3y8EJr3HKjIC0q9TgVFbNr+NyxCoHpdSTjRRER7KOudlKXkDaHJeHY1xAupBm2ZpKsnlHpskSVaaQSROcavXeuAbHdcyxzpKJULqw5tdxOSq2zxCRGcB38SaytQXPK6UOqKNcWhAoh3RSjLMcTHMrlWR5aL4RDHJ5GFLVrfJjQs2wHMIxB1PqdVylSBt0HZejmoD0rcD1QA54K3AHcGc9hdKFoCvrmJYU2ZxZ2UpuuM7BhJvKLQxmB/0DeYMrpM2oJXH9OofCqVbvjStQwPmAtCFrnBAhlRASor8FXI5qlEO7UupxQJRSK5VS3wJOqa9YehC0zBjbkjTOrRRsBKacaoOeNG2GpAAOrpA2Z40HpwvrfU0Xp7L25xzth+c4vnUm+Qwr/a+LOKrpytovIglgmd9Gey0wtb5i6UHWP7q0GagcwpWaJpy48vKmzSgeUqFTrSlrHPT9aUmaMR9hkFsplUApL4OwJRU1KFIPvMQK73tT4n1xVGM5XA6MAT4HzAbOBy6sp1C6kHVc0knvZjKt8V4+W8kQt1K43Qfo37BscCaNKWvsue5MmXE8OJXVHIsyIZ52MKV1TRwVLQel1Dz/2x3ARfUVRy9yjksqkaAllTAzIC1CMilkHNfz6Sf0PXE5+cpSw1wefsGTCWscxKEgmJuh90YbVsDhorJxzRSqAq7f3BA8l+6oTGUVkR8rpS6PK4bbWYrg0kkhnUwY1T5DKZXP/sl3OXVc2hLJJksWT3iGNBigHNzBaZbgyRwMWNKRoM4BvM1W9zVWKnqNdWbQGo9iyyHISNppi+E8t1KCdFKMqpAeXL1b8OEHikJHghNXKpkgZcDwnEImzeAupzorh2BmBhhiORSlsoL+7sbgkAPmJIPEERtzUEot8L89Rin1ZPgfcExjxGsuOUeRSgotKbMC0vnZuwmMqjhOhAN5mp+4CoNozJmPEBRogRmWQ3H7DNA/UWGw607/NS5HNQHpqODzR2ssh5YElkOLYUVwxbN3Qf+Nyxl04tJ/eM7gCmlzTrUFBax/EL2wxuZ0OQ2aG4JXT6L7NVGOcjGHjwDnAPuLyMOhH40DttRbMB3Iusp3K5kVkM5bDoZsXEqpfPUumGI5MMh9APor4HAmTVta/zbjYevMmDV21SDLobNHb3nLUS7m8A9gPTAZ+J/Q893AonoKpQtetpJ5AWknfFMZEOBVoRgJmGGOF5/CwQyXRzjm0Ke5cnAiAtL6K7TQocEA66wc5RrvrQRWAic1Thy9KASkE0YFpJUvqtdHR39zPJyyCBhRWeqoQnC3zZC4jiqK62zt1btZYLh9hgmHHBjsVjJlCFQcFWMOInKGiCwTke0i0iUi3SLS1Qjhmk2QytqSEu0vyjB5yyGcSaPxqTYcIwEzLAelwpaOIZ1kQ8FSE4bnuOH+VYa4lVQ4IG3AGpejmvYZPwDeq5R6sd7C6EYhIJ0YNDJUdwrZSmacuArzmM1KswyfwkF/y6E4I0z3NR4UczCkzXjxdaG7vOWoJltp486oGKCQypo2rH1GdLaSvhdpYRPwHpvQnjloTwJmBKRLgv4GnGqdcEaYAWsMRe0zDLiOy1GN5TBfRO4Ffg0MBE8qpR6sm1SakHVddkmnSKcMC0iHspVMMMfDgUfwTlxbdugrLwx2H5jQSTZcMwBmnGpVSGaTAtLBoaHND0irUHzKJKpRDuOBXuAdoecUMPqVQzgg7ZjzIecrSxNixCxbVexWMiEgHXYf5Kt39V1jtyjo7/X90VdeGJyokDLAPQpFQf/QwUzn7gRxVNN4b6dqthcm5yhSCaHF78Gbc70Ate7k3TRixkSq4mwlEwLSg9okGOAPD64JCVk7Gc1PtYPHhJqhHIq7soK5yqGabKWDReRxEVnsPz5KRP5jJG8qIleJyEsiskhEfiUiE0M/u0JElovIyyLyzpG8z0jJOi7plGc5BI9NwLiAdCjwCEFAWl95oTBIHszIpAnPnwBT4iTe14Q/PMc7NOirgKEwlwTMSVSIo5qA9M+BK4AsgFJqEXD2CN/3T8ARSqmjgFf8/x8Rmen/34cD7wJ+KiJNU7lZx5sHm1cOhowKjQpI6+yrDbeiAIzYBFy3YOm0JM1zK7UZULgXHsUKZlTOD4pFGeDSLUc1ymGMUmpu0XO5kbypUuqPSqng/3gamO5/fzpwj1JqQCn1OrAcOH4k7zUSco5LKpkg7Z8ATGm+F+7Kmk6KN/VL4xNiSbDUgCwPV6l8nYMJnWTdGMtB5x5WJRalAbGo4nRhGN2Ww2YRORB/poOInInXVqNWfAz4vf/9NGB16Gdr/OdKEJGLRWS+iMzftGlTDcUpkHFUvvEemOdWMmXqV0mFdMgfritOka9e+zUONbEDMyyHSItSY3lhcLZSoa+Z3jLHUU220mXAjcChIrIWeB04t9IvichjwB4RP/q6Uuoh/zVfx7NC7g5+LeL1kTuEUupGXy7mzJlTl10k57r5YT9goHII93gxwK0U5Q/XNZAXbn8N+rf8UCWncBMsh6g11vsedIq6soK5lkM1ykEppU4VkbFAQinVLSL7V/FLp5b7uYhcCLwHeJsqHBHXAHuHXjYdWFeFjHUh5xS6sgLG9FcKZyuB/qfacFYKDG5kp6tyCLsPIKg41nmNva9Ra6wrTsQa677RqpC70YTWNeWoxq30AIBSqkcp1e0/d/9I3lRE3gV8BXifUqo39KOHgbNFpNVXQDOA4nhHw8g4br5COnhsAuFsJdD/xBVVoAV6n7jCWSmgf8VxSUDagGaB4bRQ0F8BQ1G2kiFDoOIoN8/hULysoQkickboR+OBthG+70+AVuBPvj/xaaXUJUqpJSJyH7AUz910mVKqaVdvznFJJxK0pIKYg74+8DDRjez03QSK/eEmpN+qkG8Z9K84Lq5zMKFZoFusgFP6u+6iYg46y1yOcm6lQ/DcPhOB94ae7wY+OZI3VUodVOZn3wG+M5L/vxY4rvdBh91K5sQcvK/JkA9fZ9O22B9ugq82XCENwfQ6ndfY+5ooqc3Qd41LYw4JtvSMKFGyrhTPJQmuY50VcDnKzXN4CHhIRE5SSv2zgTJpQaAIUknJ57Gb0nyvEJD2Hus+dCSqtxLofVNFuTxMsByKBxTpvMaOWwjuQpBYobG8EZX+oLcCLkc1MYfVfhVzh4hsFJEHRGR65V8zm5y/waaTYmCdQ1RAWt8LtKRltwGWQ3g2ApgQ9Pe+mmQ5qFCLEgj6Qekrb2mWoP7u0XJUoxxuxQsU74VXc/Ab/7lRTWAlBPMcwJyYQ3FAWvdAXumpVv8sD69ld+Gx9kH/kriOAZZDpHWmr7zFRXs7Q4X0VKXUrUqpnP/vNmBKneVqOlk3cCuZF3OICj7qfEIs7a2k/4krXCEN5riVTMoIK7XOdL+Ova/Fqaw6t64pRzXKYZOInCciSf/fecCWegvWbHK+ldCSlHwnVtOUgyntKJyiIjhTAtJSsnHpu8albiX90yzdoqB/myHXcXBZpJIJkgmzRgyHqUY5fAz4ELABr23Gmf5zo5p8QDphXhFcSbaS9ua499UkX63XYK3w2Nu49FZmUJourPOp1i2OOeiugIvcuaB/vK8c1VRI9yql3ld3STQjiC+kkkJLysyYw+BsJX0v0NJBNAZ0ko3YuHT23xdbk4mEl4Wn82bruKVBf8dV+SFculG8xuBZaDpfF+WIXWERea+IbAJeEJE1IvKGBsrVdALLocXgmEMhMKb3JlDSW8kAy6HErWSI5VAcJ9FZAStVlMqq+QwKpyjWB2ZbDuXU73eANyul9gQ+CHy3MSLpQS5vOSSMizkUbwS6m+NxdQ46u8KUKt1os47Kr71uFKdZghktP4pP4aDvxD23yJ0L+qc4l6OccsgppV4CUEo9A4xrjEh6EGQrpQ3srRTVPsNxFTlN5S9U73pfTQhIR6Wygr5xqeLaF9DfcnAiaklAf8shWZzirPEhpxzlYg5TReSLcY+VUj+sn1jNJ1znYF5AushyCJnjKQN8tamEkNB8QFFUDj54cZL2Fv06yRbXvoD+WWylnW/1jkUVz58AXwFrfMgpRznl8HMGWwvFj0c1QYV0KiEk/X/GuZVC8xzA22zHtjZNrFicopvKG1CU1HYTgOgcfNBXoRVbk+AN/NH5VFvaeE9vyyHaOtN7jctRrrfStxspiG4ELqSgdUY6KcZkK0W5lUDfE1e+YZlJp9qIHHzQ1xVWnN4M+gfRSzLCdA9Ix1hnPQP6Ngssh34+Bk0IAtLpRKAcEsa4ldwiH74pN1XJYBeNT1xRqaygbzuK4vRm0N9ycNzBLprCaFM9FVpxvQ7onwxSDqscYsiFurJCkI1ixoccla0E+p5qiyfBgf5jN11V6lsG/de42HLQ2R/uNd4rPC6MNtXzPiyu1wF/jTVVZpWwyiGGwK3UkipYDqYohzi3kq6nxCjloHsKYLFbSXfrLB9DS5oTcygN+uttOcTVkuh6TVSionIQkd1F5GYR+b3/eKaIfLz+ojWX4AMNOrJ6ysGMmEO85aDnRRro3LDLQ/+AdFwOvp5rXFxoCPpbDo5b7LozxHJIFFvAespbiWosh9uAP+C17AZ4Bbi8XgLpQvCBBifCdFKMqXMoLnjSPliqgsywwuWoe5txRxU33tM76B8ZLNU8ruO4gzvfmlIEVzr3Wk95K1GNcpislLoPcAGUUjnAzL92CATB59akd0GmkwljJsGVtmfWPVjqyVV8EtdVmUFUhbTm1llMXMdaDrUjqghutAeke0RkN0ABiMiJwPa6SqUBwcYUWA4tRgWkva+lRXB6bgSBvKmiecG6KjMoLdDS3TqL7Riq8Ro7rhoUI2nV3HIobj0P3hpncm5+TrpJVNOV9Yt4k+AOFJG/4w36ObOuUmlAxuiYw+CTeJuBlkNrSv9TrUlFcI4qVQ6BdaaKXGS6kHMV7TL4wAD6uu5UzBqDd10E35tCReWglFooIv8CHAII8LJSKlt3yZpMJud6bRz8D9qkmEO4uhv0v6mK5QW/QEtTZQb+PIcol4emaxx3qnWV14q+JaWfcnCVGnRNtCQTiMZtVeLWGLxEhVGnHHyOB/bzXz9LRFBK3VE3qTRgIOfmP1jwLIcdhlQ65hzP5ZFPZdV8PkJUsNTrg6+nvBDR90fzVNa4WQPgucJaUvplteccRTKUpCAitGmcxRYV1ynUZjhMIN0UuYZLReUgIncCBwLPUQhEK2BUK4dMzh10w7QYVOeQc9WgBnu696QJqtEHbVyaB/JyRcHSNs2bwkW2z8hbOy7j2pohVXkcd7DlAHrHoqLawBSquvWUuRzVWA5zgJnKxIjKCBjIOXk/MngBaVM+4Jzjki5yeYjou3FFn2r1TgH0Nq6CAg4mq+m6cbkR7TN0r5zPue6gawL0tiij2sAE1pnO8bM4qrElFwN71FsQ3Si2HEyqdCw+1XpdTvWVvxBzKKx3WzpJTuMZFDmndOPSOac9KiCdd3noqtAU0cpB0+s42j2qdyyqHNVYDpOBpSIyFxgInhztc6UzzmDloPOJpZicWzpjV2f5y95UOZddNJxBkYtwebRqXJsRXQSnv+VQssYaK+DgkBO+9/KWg6YKuBzVKIdv1VsIHRnIDg5I63zyLqa4eAjQOpCXH8kaESztzzrs0lpt3kTjKI7rgN7+8Ki+P22aB9EdRw3KCAO9DzmBlbvTWA5KqSfr9eYi8m/AVcAUpdRm8ZKtrwHeDfQCH1VKLazX+5fDZMsh66gIy0HnjcuTq3gQDeh7U0UHS/W9RspZDtrKrCKsM41jfwXLwZw1Lkc1jfdOFJF5IrJDRDIi4ohI10jfWET2Bt4OrAo9fRoww/93MXD9SN9nuMRZDibE5SMtB503rqhNQGN/uFIqZo31tS6jm8JpbjnEXccau8GAQem3hYC0nmtcjmqcuT8BPgIsA9qBT/jPjZQfAV/Gb8vhczpwh/J4GpgoInvW4L2GzIDj0hLKVmpN610BGybruIPaDkAwy1ZP2YsD6KD3iSuqaA/0dt1Ft5PWux1F1HWhcxZbtHvUXLdSVZE+pdRyIKmUcpRStwInj+RNReR9wFql1PNFP5oGrA49XuM/13AGsk6J5QBmKIcol0erzpaDE53PDnoGS4ONtjTmoO8aR81z0L1wz3HiLGA95Y2cmaF5P6hyVBPp6xWRFuA5EfkBsB4YW+mXROQxolNgvw58DXhH1K9FPBfpxxGRi/FcT+yzzz6VxBkyUTEH8D/kdr0rHbPO4Bx88OTf3qdn15OcGx14BD3dSkExZJRC6+zRT14oyJyOcnlounFFuRu94kg95Y1LyQY9r+NKVKMczsezMD4DfAHYG/hgpV9SSp0a9byIHAnsDzzvN/uaDiwUkePxLIW9Qy+fDqyL+f9vBG4EmDNnTs0DAZmcS6tBVcZhHLfUrdSWStCh6yYQE9wFPTeuqOAu+NaZrhtXUUsV0P+a9txK5iRW5CIODW2a99wqR1nlICJJ4DtKqfOAfuDbI31DpdQLwNTQe6wA5vjZSg8DnxGRe4ATgO1KqfUjfc/hMJBz82Y36L1ZFROVg6+zy8NR0ZsA6LlxRWWlgN5jN7OuG+kGA31bO3gB6cHPaX0dR7iVUskEyYRoe2goR1nloJRyRGSKiLQopTINkOd3eGmsy/FSWS9qwHtGksm5+XbdoP8pK0wu0q2k74krMuagcUC6YDlErbF+8gJkc2rQ9Qx6d5ItZIQVyeyP3dSxzXjWKXUrgWc96HrvlaMat9IK4O/+qb4neFIp9cNaCKCU2i/0vQIuq8X/O1L6ss6gFrtmWQ5uZMxB19NLdFaKvr7a+JiDvqfaXISrMZUQEqJn3x83aGIncYkK+rXADup1Sly6Gl8X5ahGOazz/yWAcfUVRw9yjksm5zI2VJlrkuWQcRRjWj8t/IcAACAASURBVEpPibpeoNEN1vQ91Ua5D8C3HDS9PqKSFETEG/ijsQJOp+ItSt2UQzYilRX0zrAqRzUV0t8GEJHx3kPVXXepmkyvvyGNaTHTcihuGgiFC1RPczxaXtDzVJuLCUi3pZI4riLrlPa2ajY5xy2JkYC+bWEC5VDiCtO4ODIo2iu+v1rTCS2v40pUUyE9R0ReABYBL4jI8yIyu/6iNY/eAe+DHGQ5aBwgLSYTMbwl2Gx1nGaXKesP10/eXJxvWeMDRK5oHnOAri6PYExvSRsYjWNR2QgLGIJEBf3krUQ1x5tbgEuVUvv58YHLgFvrKlWT6cl4E98GWQ4aX5TFZJzBabig92abdVzSqVKXh9dHR7/1zsX6ljVf40Tp7a6v5VDa4RT0tigdRw2aoxKgczJIOapRDt1Kqb8FD5RSTwGj2rWUtxxaTLUc4t00Om62WcelJcbloaMyjvMt6zyONc7VpavlkHcrlVzH+irgqMQK0HeNKxEbcxCRWf63c0XkBuCXeNXKHwb+Un/RmofxlkMZ5aDjTZUpqikJ0DWQF7g8wpMCYfBMZt3IOdFuJV0th8D9WVJLovkhJ04Bd/Xr2Z2gHOUC0v9T9Piboe/1b006AnoD5WBozGGgqEYDQm4lDTeurOOyS1vppahr+m2gHEoUsM6uu4j5E6DvgKK4gHR4CJRuZHKDOzkHmOpWilUOSqm3NlIQnejJu5VCXVlHjeWgn/yZiPkT4LfA1vCmyjjeGsa67jTcbIvnige0phJ09+eaIFF54gLSOt+HAxH3HejdrbccFVNZRWQicAGwX/j1SqnP1U+s5hJlOSQTQjop2lsOrqvIuSrWV6uj/F7MIcatpOFGm7cc4oKlGiq0OLdSWzrJ5h2NaH4wNAp1DgYdcmKUQ6um7tFKVFME9zvgaeAFwLy/cBj0ZkotBzDjBJCJDeTpe1NF1TmAvus9EOdW0rhwbyDnMHFMS8nzumaEZXKe5zrOraSjRTmQc0riUBBYwPqtcSWqUQ5tSqkv1l0SjQiUw5iWwcvTqvGkr4CBuFNtSt9TbSYXU6CV1tvlUexf1tly8NpNRKWyJrW8pgvZSjFtVXS0KOMOOZpawJWoJpX1ThH5pIjsKSKTgn91l6yJ9AzkSCel5INu1fQkGyY4obQXWz0an2pNS7OMs850bmTXn40/1eoob96tFOu600/m4madAW2pJFlH5duumEI1lkMGuApvSE/w1ynggHoJ1Wy6+3ODqqMDdA2QhilYPdFpljreVAO5eOWg46m2YsxBw1PiQEwmja6WQ3xAWt+MsIGcy9ix0fsGePde1L6iK9VI+kXgIKXU5noLowvb+rLsGumf1TPtL0ygHNqLmpIVUlk1vKmybomlA0GrY/3WOz6VVXe3kjmWQ6Bgi6/jdDAfQUOZYy2H0MHMJOVQjVtpCd5shZ2Gbb0ZJkSMAjUhX7kv6/no20viJXoWDzmuIuO4JZsA+A3LNJMX4pVDq8auu/6imegBrakkOVflp5jpQl/GkydSoWk6HyEuW0nn2oxyVKPGHLz50U8AA8GTozmVdXtflkljSy2H9pYkfRre+GGCm6rUraRnKmuwkUYpBy9bSS95wVvDZEJK22ekEojop4ChdLJhQPi6iCqSaxZ95a4LTQv3YuscNHbplqMa5fBr/99Ow9beDAdO2aXk+fZ0is6eviZIVD1BjUbxTdWS9DYu3S7QYBNoi3Ir+VkeurUZ7804jEknS2QKmgXqdkLMOS6Oq/JurzDhOSVjWxstWTz9+evCnLYqfVmn5FAGehfulaOaeQ63i0g7sI9S6uUGyNR0tvVmI91K7S1J+jL6pVaG6YvJVhIRLesG+mJiJOCdapXyGt0VpzQ2k75sLjJGAnpmWAXKKq5/FehX1d2XcUhIadAf9J2P0JvJlaS/g97NAstRzTyH9wLPAY/6j4/xR4aOSnKOS3d/joljSpXDmHQyH/DVlR0DnvIaG3OR6naBlnUraZr905uJPiGCnoV75RSwrsNzgklvURajjvMRHFfRn42OnencLLAc1TgZvwUcD2wDUEo9B+xfR5maSpdfdDUx1nLQ+wMOisbGt8c0stPsAs27laIKtDT11fZmnJKAf4COCrgnODBEZMoELg/tLIesE7nRgp5rHFzHY1vjXXe6HXIqUY1yyCmlthc9Z1Y1xxDY2uv1mYlqNTDGgIB0V1+WVEKiT4ka+sPjUm+h0OVUt9qSvnKWg4YKOLAmTXJ59GXiZ0TruMb5WF/kGuub4lyOapTDYhE5B0iKyAwR+V/gH3WWq2ls6vYSsibvUhqdG9Pipf1lNNtgw3T35xjXloo2x9P6meOBpTOuLSp1WE/LoSeTi1UOremkdqNYAwW8SznLQbM17h7wruModGxHEQwIGxNzKAO03jeiqEY5fBY4HC+N9ZdAF3B5PYVqJh2+cpg6vlQ5BKcCnV1LXf1Zxke4xMDvDqnZBbpjwBuCEjfPAfQ7cfUOxFsOrUn9qugLbqXooD/ol4O/oz9eObRqWOcQNSAsoMVQ5VBNtlIvXuuMr9dfnObT0dUPwNRx0ZYDQG82xwSiN+Bm09mTiYyXgJ4VxwXLoYzLQ7NTYld/lvERlg54AV7dkhaCjatszEG362Igy5QI6x30dCsVYn2l14WucZ1KlBsTWjYjSSn1vtqL03w2dQ/QkkxEp7L6J1ndbv4wm7oHmL7rmMiftaWTbOvVq3d/OeWga354V190qjN4qZdbc3qucZRbSWfL4YDJpbVG4Mms2ym8q8+zgKOui5aUngWolShnOZwErMZzJT0D6JNoXkc6ugeYMq410mcf5Lbr7Fbq6B5g1r67Rv5MxyyPrv4sLclEbMdQ0MsczzkuPRmnjOtOP7fStl5v44rrFwYaWg79uUhXI3gKWLeNdruvHKIsytZRqBz2AN4OfAQ4B/gt8Eul1JJGCNYsOrr7mRLhUoKCW0nXjKVMzqWzJxPpEgNoSekXLN3Wk42sKQE9T1xBqvP4MhuXdmvcm6EllYhJF9ZvjR1Xsa0vy6QIZQZ6zr0OrotIyyGp3yGnGmID0kopRyn1qFLqQuBEYDnwFxH5bMOkawJrt/YxbWJ75M/yMQdNLYctPX4wfVxb5M91nPrV2ZuJ7GMFevpqO/013rWMzDpaDruOSUdawzq67rb2ZnBcFXtIa015loNS+mTUb+vNkJDoxIpEQrS0dipRNltJRFpF5AzgLuAy4FrgwUYI1gyyjsuarX3su1u0z749HWQr6dlCo6MrUA5xloN+p9rOnnLKQb86h8IaRytgHdd4S0+Gie3l11gnmTfviE8nB09mpSCn0fCcjq4BJu/SSjIR7X1vSekXJ6lErHIQkdvx6hlmAd9WSh2nlLpSKbW2Fm8sIp8VkZdFZImI/CD0/BUistz/2Ttr8V7Vsm5bHzlXsd9uYyN/rrvlsH671xRwjwkxG5eGp5ctOwYqKgedNq6N3X42W0SqM+hpnW3s6i97TYBeCjioNYq3HAKLUh+ZN3T1s/v46DWGwNrR67qoRLmYw/lAD3Aw8LmQSSqAUkqNH+6bishbgdOBo5RSAyIy1X9+JnA2Xl3FXsBjInKwUqohq7pyize2Yp8Yy0F35bC8YwcA+0+OVm66zcB2XcW67f284/A9In9emEGhj8yB5RC3Eei2xgDrt/dz+F7Rt2siIaSTopUCLlgO0YeGfCwq60RmYDWDjV39TN812h0No8xyUEollFLj/H/jQ//GjUQx+Hwa+J5SasB/rw7/+dOBe5RSA0qp1/HiHMeP8L2qZuWWHoBYy6FN82ylZR07mDaxPXbaVGvSu0B18dVu7hkgk3Njb6pCloc+672xa4AxLcnYTakl6VXR6zIvOJNz2bxjINZyAD+IrtHGVdly0C+I3tE9wNSKloM+8lZDs6Z7HAy8WUSeEZEnReQ4//lpeOmzAWv850oQkYtFZL6IzN+0aVNNhFq5pZe2dCLWZx+UxuuarfTKxh3M2D06NxwKJ/Gso8fGtWar5waLSwBIJcQbnqPRTdXRXcF9oFn67Ua/qHPPcspBs1Ptxq4B2tKJWAWsW4bVQM6hsyfDHmWuC93WuBrqZpOJyGN46bDFfN1/313xsqCOA+4TkQOIrqWI3MmUUjcCNwLMmTOnJrvdii297DNpDImYoFIqmaAlqV8FLHgtEpZt7Oath0yJfU3ev5xzIidWNZq1gXKIsRyC4Tm6bAIAq7f2sdfE8qdw8JRD3MyHRrJ+u6cc9pgQ7/LQbTb665t72G+3sbEDnoKYgy6b7fptwRqXsxz0WuNqqJtyUEqdGvczEfk08KDy/BtzRcQFJuNZCnuHXjodWFcvGYtZ1dnDvjEupQBdB/7MX7mVnKs46cDdYl+jW4+XVZ1ejCfOcoAgNVSPm0opxasdOzhz9vTY1xROtQ5o0GLltU1eHGq/mDga6HeqfW3TDg6fNiH257q5G5f5sb6Dppax2jXMYqtEs46PvwZOARCRg4EWYDPwMHC2n0K7PzADmNsIgVxXsXJLL/tOir+JwAtK62g5/PWVTaSTwuyY6mjQL/tn6foupu/aHtmRNUCnm2pDVz87BnJlN4GCdaaHzC9t6GZMS5K9Y1qqgF7ptwM5h1WdvRwYk1QB+hVHvrKxG4AZ5a6LlH6V85VoVqj/FuAWEVkMZIALfStiiYjcBywFcsBljcpU2tjdz0DOZd8yFyV4lkOvJifZAMdVPLJoHf9y8NTInv0BLZrVDSxd1xWbRROgUzuKVzZ6J8Rym0A+w0qTjeulDV3M2H1crKsUgvRbPeR9taMHV8GBZU/hemWxLdvYzZ4T2ioecoIeV6bQFMtBKZVRSp2nlDpCKTVLKfXn0M++o5Q6UCl1iFLq942SKUhjLWd+g9d8T7dspT8t3cDGrgE+OCsydp8n76vV4JTY1Z9lxZYeDt8r3n0Aga+2+fICLFy5lYTAYWUUWjiu02wcV7F0XRcz9xxX9nU6WQ4LVm0FYNY+lS1gHdYY4LnV2ypex7q57qqh+VFJTQjSWPedVN5yGKPZqNCc4/Ljx5axz6QxvH3m7mVfq5Pl8PSrW1AKjt9/UtnX6VQ8NPf1TmbuNT62XTfola304vouuvpznLB/fBwK9CqOXLCikynjWsvWDOi0xhu7+lmxpZcTD6h0HZsXkLbKwWflll5SCSmbiQLewB+d3Eq3/n0FL23o5orTDiWVLP9x5gPSTvPlf/KVTYxpSZY9IYLvq9VgE+gZyLFw1daKG61OOfhPLd8MUDZJAfQ51Tqu4qnlmznxgN1iM5VAr7jOU8u8Na6ogDVZ46FglYPPyi29TN+1veIGOyatT7bS3Nc7+f6jL/H2mbvzriOiq4zD6LJxZXIuv1+8gZMPmVIxpVYXf/gTL3cwkHMrWmc6jYT8/eINHDFtfNm6DNDHdbdw1VY278jwjkprnNanIePvF29grwltHDGtQuxMk0POULDKwWfFlspprKBPttKSddv51J3z2XvSGK4+6+iyJ60AXbI8/vxSB509Gc6avXfF17amkgxo4A+/f8Eapoxr5bj9KrsPoPlrvLxjB8+v3sZ7jtqr4mtbUwkyGmy0985bzZiWJCeXqdUBfQ45Hd39PPlKB6cduWfF+89aDoailGLVlt6KwWgI6hyaeyMtXrudc296hrZ0ktsuOi52KlkxOvSVV0px699fZ4/xbbx5xuSKr9ehkd1rm3bwl5c3cd4J+8Z23Qxo0SRYevNTr9OSSnBWmZqMAB0C0lt7Mvzm+XV84NhpZbN+QB/r7O6nV5F1FOeesE/F1+pinQ0Fqxzw2kZ3D+TYxwDL4U9LN/KhG/7J2JYU9158UlXWToAOk9X+8eoWnnm9k0+ffGBFFx54LoRmbwI/eWI5LakE51S1CTR/jddu6+PBhWs449hp7BbT9jqMDq67G/76GhnH5cI37FfxtTpYZ9t7s9zxzxWccuhUDpgSn3YbEChgXfqaVYNVDsDKzurSWMFPZc06Df+QlVLc8OSrXHznfA6augsPXvqG2O6xcbQkm3tTDeQcvvHQYqZNbOfDx1V2KUHzfbUvrNnOr55dy0Vv3C+2EVwYHVx33/3diwB85pSDqnp9sy2H1Z293Pr31/nAMdM4ePfyabcA6aRnvTXTorzm8WVs78vyb+84pKrX6+IKGwp69LttMvk01qrcSt6S9Wcb1ztnIOfwH79azP8tWMO/HrUnV5959LDeu9kpgD/583Je3dTDbRcdR1u6Ovmbmco6kHP49/ufZ/IurVx6cnUbbbP7/jz5yiYeWbSez79tBtPLVEWHaWZXVqUUX31wEamE8KV3VrfRNrvn1qI127jjnyv48HH7MLNCEWdAuDtBtdd+s7HKAS9TSYSqbqbCTIdcQ5TDpu4BPn3XAuav3Mrn3zaDy0+dUVXwOYpmFmg9tWwzP3liOR+cNZ2TD5la9e81s+3ADx59mZc2dHPzhXOqjus0s0BrU/cAX7rvOQ7ZfRyfPvnAqn+vmenCP//ba/x9+Ra+84EjyvbYKqZZymHHQI7P/fJZpo5r5avvOrTq3xs01bB88pg2WOWApxz2mtBelUZvDw38KZ/ZPHIWr93OxXfMp7M3w3XnzOJfj9pzRP9fsxrvrdnay+fueZYZU3fhyvcfPqTfbVa20p1Pr+Tmp17nwpP25W2HlU+tDNOsQsNMzuXz9zxLd3+OX3zyxCGdTltTSRx/BkWlgHsteXTxBr77+5d495F7cM7xleM5YVrTjQ/wuq7iK/cvYlVnL/dcfBITxlTfWLElZDmYglUOeGms+1RouBcQWA71nunwyKJ1/Nv/Pc+kMS3cf8kbOKJMl8pqaUawdFtvhgtvmUvOcbn+vNllez9F0ZoqDCgarsU0VJ54qYNvPrSYUw6dyn++Z+aQfjeVEBLS2E1AKcUVD77AP17dwtVnHV2V3z5M+NDQKFfps6u2cvm9z3LM3hP54YeOGfJn61V1N9Y6+/6jL/HbF9ZzxWmHVqzsL6bQD6r5KcPVYgPSeAGxauINUP9RoUopfvinV/jML57liL0m8NBn3lQTxQDePIpEA4fn9GcdPn77fFZv7eOmC4/jwCqyOopp9GCXxWu3c9kvFnLYnuP5348cW1VGVRgRabib5kePLeOBhWu4/NQZZduJx9Foi3J1Zy+fuH0+U8a18vML5gzLB9/ocax3/HMFN/z1Nc4/cV8ufssBQ/59azkYSG8mx+YdGfau0nJoT3tLVo9ah6zjcsWDL3D/gjV8aM50rnz/EfkTR61oTSUbcoE6ruJzv3yWhau2ct05s4Z80goIpy3WO5C3blsfH799HhPb09zy0eNix61WojXVuPTb++at5trHl3HW7Ol8/m0zhvV/DI6T1HcGxfbeLB+9dS45V3HrR49nchWptlF4cz4as8Z/WrqRbz28hFMPm8o33ztzWBZsa5PcjSNhp1cOqzu9aWTVupXa826l2rbQ6BnIcdkvFvKXlzfxhVMP5nNvO6gubpRGVGoqpfjWw0v449KNfPO9M3n3kcOPlTRq4+ruz/Kx2+bRM+Bw/6dPqthyohwtDcqwevKVTVzxqxd484zJ/L8zjhx+okKD0iwzOZdL7lrAqs5e7vz4CWXnYlSiUVlsz6/exmd/uZAjp03g2mFYkgE61GYMlZ1eOQTTyKq1HOrhVtrak+Gjt87lhbXb+d4ZR3L2EINzQ6ERG9dP//Iqdz69kovfcgAXvXH/Ef1fjQjwZh2Xy37xLMs7dnDrRcdx6B7VpSfG0YhMmiXrtnPpXQs4ePdx/PTcWaSHuWlBY4ZABSmr/3xtCz/68NGceMDI0jlaG3DIWbWll4/fPo8p41q56cLjhhwvC6PbFMZqsMrBVw5VWw7p2iqHzp4M5970DK9t2sEN58+p2NhtpNR74/q/+au56g8v8/5j9hpSql8c9d64lFL8568X89dXNvH9Dx7Jm2eU7+tTDfWOOazd1sdFt85jfHuaWz96XMV2E5VohMvj2seX8+DCtVx+6gw+cOzQ4yLFtNR5eM7Wngwfvc1zf9120fFVFUCWQ7cZFNWw0wekV3f2sktril2rTEvLZyvVQDmEFcPPL6i/YoD6blx/ebmDrz74Am86aDI/OPPostPHqqXeU7+uf/JV7pm3msveeiAfPq42Fls9/eHd/Vk+dus8+jIOt150XNmh9tVS72Dpr55dw48ee4UzZk0bdlykmHr2KhrIOXzqzgWs6ezjxvPnDCuRohhrORjI6s5e9p40pmp/bWBajtRy2NY7WDG85eCRn1iroV7B0kVrtnHp3Qs5ZPdxXH/erIqtuKulkK1U+xPXH5Zs4AePvsz7jt6LL729uurcaqjX3Ous43Lp3Qt5ddMObrvo+BG7vwKCtir1uC4WrtrKV+5/gRMPmMT3zjiqZnE0L1up9teEUoov37+IuSs6ufYjxw47kaIYE9tn7PSWw6rOXvaZVH1lZls6gcjI6hyCFM9XG6wYoD4B6TVbe/nYbfOYNLaF2y4auZsjTL1uqmUbu/nivc9x9N4T+cGZR9XEygloqUMnWaUU33x4CX9btpnvfOAI3lRFR9tqqVcn2Y6ufi65cwG7T2jl+nNn1+zAAPVrFvijx5bx0HPr+Pd3HsL7jq7c7rxa6mU55OoYJ9qplYNSylcO1TewExHa00l6B4bn78w5Lp/1Uzyv+fAxDVUMAK01Lh7qzzpcctcCBrIut110PFNHkOUTRT2yPLr6s1x85wLaW5L87LxZNU+RrYfl8PO/vcYvnlnFpSfXzv0VUI/iyEzOs3K6+3PceP4cdh3bUrP/G+qTkn3/gjVc+/gyPjRnOpcOof1INRSu49rde30Zhw/+7J/c8c8VNfs/w+zUymFT9wADOXdIygFgXFtqWMEwpRTfeHiJlzf93sM5bQQpnsOlNV07y0Epxdd+9QKL13bx47OPGVFqYhyFYGltbirXVXzhnudY3dnLT8+dzZ4Tqrcaq6XWp9rHX9zId3//Eu85as+qu4AOhXqcav/rkSXMX7mV7595FIftWRv3V5haz/lYsHIrVzy4iDcetBvf+cDw04LjqHW6sFKKrzywiEVrtrFXHa5h2MmVw+qtXqbS9CEqhwntabb3ZYf8fsHp79MnH1hV3/p60JKs3an2jn+uzGegDKX/0FCotVvpx4+9wuMvdfCN986smT+5mFoOk399cw+X3/sch+81nqvPqk2Qv5haB6TvnbeKu55exafeckBNXTNhapl119Hdz6V3L2DPCe389JzZI0oLjqPW1/Etf1/Bw8+v40tvP5hT65TIslMHpIeaxhowoT3Ntr7MkH7nyVc28T2/ydiXq2xNXA9a07U51c59vZMrH1nKqYdN5XOn1CYDJYpatsD+w5INXPvn5Zw1ezrnn7jviP+/OGo1H6FnIMen7pxPKiH87LzZdasQr2Uq67OrtvKfv17Cm2dM5ss1SGWOI1AOI+25lXVcLrt7IV19OR689PghNdMbCrWcwvj0a1v4f797kXfM3L3qVvLDYadWDqcfPY0TD9iNKUMs4Z/Q3sLabX1Vv37F5h4++4uFHLz7OK46s7p5z/WiFpbDhu39XHr3QvaeNIYffviYupxmA2rVW2l5hx+Anj6BK99/RF0/g1q4lZRSfPmBRSzv2MEdHzuh6tkMwyHv8hjhddHR3c8ld3kB6GvPPrauHV5bfUWZcdwRtZj5zm9fZN6KrVxz9jF1cX8FJBLiNwsc2Rqv397HZ36xkH13G8P/fKg+lmTATu1WSiSEPSe0D7kkfkJ7mq4q3Uo7BnJ88o75JBLCzy+YM+x+PbVipNlKAzmHT9+9gL5MjhvPn834GmYmRVGL4qGu/iwX3+EFoK+v4wk8oBaWw01/e53fLlrPl991aE0zk6JorUEqaybnculdC9nel+WG82ofgC6mFifxBxeu4bZ/rODjb9qf04+ZVivRYqnJvXfXQvoyDjeeP7umWYFR7NTKYbhUG3NwXcUX7n2O1zb3cN05s6pu0VFPRlo89O3fLOXZVdu4+qyjmTHE1tDDYaTZSkEAelVnL9edM4u9hjBQZriM1HL4x/LNfPf3L3LaEXvwqWF0AB0qtaglufKRpcxfuZUfnHl01dPRRsJILcrFa7dzxYNe/cUVp9XP/RVmpP2gvv2bpTy3ehv/86GjOWhq/e89qxyGwYT2NDsGcmQrnA5/8sRy/rR0I19792G88aD6nv6qZSSnl3vmrsoH1BuVaTXS3krXPL6Mx1/q4D/fM5MTRtjPp1pGkma5blsfn/3lsxwwZReuOqsxLsiRnsLvm7c630urXgHoYkYS4N3Wm+HTdy9g0tgWfnLOrGE30xsqI7n37pu3On/vveuIxtx7VjkMg4l+0Kqca+mJlzr40WOv8P5j9uJjb9yvQZJVZrinl2dXbeUbD3mBxnqkU8aRTAjppAxL5j8u2cA1jy/jg7Omc8FJ9QtAF9OSSuC4asgFSgM5h0vvXshAzuWG82ezS4NckImEkErIsDau51Zv4z9+vZg3HTS5oYkWwx2eE1jzG7b389NzZw27ZfhwGG6G1QtrtvMfDy3mjQft1tB7zyqHYRDME94WoxxWbunh8/c8y6F7jOe7NWwZUAtaUgmyjsJ1VdW/482xXsjuE1r534/UN9AYRUty6Ceu5R07+OJ9z3PU9Al85wP1DUAXM9xT7ZWPeG6Dq886qib9fIbCcHpudXR7FdBTx7cOazDSSBhuQ8b//fNynnh5E994z0yO3WfXeogWy3A6Im/t8aycyWNb6h7kL6YpykFEjhGRp0XkORGZLyLH+8+LiFwrIstFZJGIzGqGfJXY0292tnZracZSX8Zr2iUi3HDe7IaNXayWfGpolTeV1856Idv6Mtxw3hwmjqlvoDGKoc4L9iqg59OaStQ1BTSO4RSVPbBgTb42oFFugzBDbYGdyXkpoNv6MnWpgK7EcNyNf3m5gx8//gpnHDuN8+qYyhzHUPuaOa7i94iGzQAADkNJREFU8/c+R0fXAD89bza7NdDKgeZZDj8Avq2UOgb4hv8Y4DRghv/vYuD65ohXnv0mjwU8CyFM0LP+5Y3dXHP2MexT5ejRRhJ0le2psv3HlY8sZe7rnXzvjKMaEmiMYiiuMNdVfPHe51i1pZfrzm1MALqYoQbRl6zbztd+5QVH/71JNTBD9Ydf+chS5q1oXAC6mKGu8erOXi6/9zkO2X1cXSqgq2GobqVrHl/GX1/ZxDffN5Nj9p5YR8miaZZyUEBwRU0A1vnfnw7coTyeBiaKSOOPURWYOq6VtnSClVt6Bz3/sydf46HnvKrFkw+Z2iTpyhOk0vYMVN5sfzl3FXf80ws0vv/Y+qf6xTGUm+qqP77MYy928B//etiIB8oMl6FYDlt7Mnz6roXsOqaF//1I44KjxQwl/fbOp1dy59Mr61oBXYmhZFj1ZnJ8+u4FOK7iZ0205oeigB9bupFrH1/GmbOnc04dh3+Vo1lJ95cDfxCRq/EU1Bv856cBq0OvW+M/t774PxCRi/GsC/bZp7GLJyLsO2ksK0LK4TfPr+P7j77Ee4/eq65ViyNll1bvxthRwXKYt6KTbzy0mH85eApfqWOlazVUOx/hV8+u4fq/vMpHjt+nae1JoPrajIGcw6fuWsCGrn7uufjEEQ+UGQnt6SS9mcrW5D+Wb+ZbDy/hlEOn1rUCuhLVVnV7luTzLF3XxU0Xzslb/c2gtcoBRUvWbedz9zzLkdMm8N91LtgsR92Ug4g8BuwR8aOvA28DvqCUekBEPgTcDJwKRK1CZORUKXUjcCPAnDlzqo+u1oiD9xjHM69tQSnF4y928KX7nue4/Xblqhq3f641u7R6wfRyymHF5h4uuXMBe+86hmubEIAuppre/QtXbeUrD3iumf86/fCmJgFUE5BWSvG1Bxcz9/VOrjn7GGY1ODhazPi2dMWNa8XmHi79xUIOmDyWa84+pqnXRbX9oK7648s8umQD//memZxyaP2HaZWjmoB0R1c/n7h9PhPa09x04ZyGx8vC1E05KKVOjfuZiNwBfN5/+H/ATf73a4C9Qy+dTsHlpBWnHjaV3zy/jkvuWsDjL3Ywc6/x/PyC5n6Y1TC2tXzMoaO7nwtumYsCbrpwTj4zq5lUavmxbGM3H7ttHnuMb+On59ancdpQCNwW5WZ+/OixZTywcA2XnzqjIdW5lRjXlmLTjoHYn2/s6ue8m59B8K6LelfnVqKaFth3Pr2S6//yKuecsI8W6eS7tKbZUUYBd/Vn+djt89jel+X/LjmJ3Wvc/n6oNOsuWgf8i//9KcAy//uHgQv8rKUTge1KqRKXkg68+8g9edNBk/nT0o2884g9uPsTJzQlk2eojGvzzgNd/aVpuNt7s3z0lnls3jHALR89jgManE4ZR2s6QX+M+2B1Zy/n3fwM6WSCOz9+PJManDUTxa7+dbC1J7o54w1Pvsq1jy/jrNnTazY2c6SMb4+3HLb1Zrjg5rls7clw+8eOZ9/dmueaCQiss7jr4oEFa/jPXy/m1MOm8u33NdeSDNh1TJqtvdHp7z0DOS66dR4vb+jmunNncfheExosXSnNijl8ErhGRFJAP37sAPgd8G5gOdALXNQc8SoTbEY5VzX9pDoUgqKfTd2DT4mbdwxw/s1zebVjBz+/cE5TsiPiGN+WZmPXjpLnX9/cw/k3P0NfxuG+S07SYtOCQpHk1t5S5XDzU6/nZzN874P61MDEzSjZ1D3ABbfM5fXNPdx20XEcNV2P6yJQwJ0RCvih59by7/c/zxsP2o2fnDNLm/tz4pg0fVmH/qwzyMOwYyDHJ2+fz7OrtnLdObN4qybJLE1RDkqpp4DZEc8r4LLGSzQ8RLzqXZOY0J6mJZUYpBzWbuvjwlvmsmZrLzd/dA5vntHY6XSVmDqulaeWbx703JJ127nwlrm4Cu7+xIk1m6dcCwLrpbOncEpUSnHVH17mp395lXcdvgc/+nBzffbFTBrTwrbeDFnHzW+ma7b2cv7Nc9mwvZ+bLpzDGzRpAQOe625ca6rkkHPr31/n279Zygn7T+LG8/Vy8waehe192bxcm3cMcNGt81i6vosffuiYpgwAi2Onbtm9MyIi7DG+jTV+y/GnX9vCpXcvJJtzuf2i4xvWf2goTB3fRnd/jr6MQ3tLkkcWrePL9y9iYnuaOz5+Ql0m0I2E9nSSsS1JNnb1A9Ddn+WrD7zAb19Yz0eO34f/fv8RWikGgOm7jsFVXjv2vSeN4R/LN/PZXz5L1nG56xMnMHvf5gbMo5gyrjWvHPqzDv/1yFJ+8cwq3nn47lxz9rFaKQYgH0NYu62P3ce3sWjNNi69eyGbdwzw8wtmNz1gXoxVDjshR0wbzz+Wb+a/H1nKrf9Ywb67jeHnF8xpeMuGatnfTz+ct6KTPyzZwN3PrGL2vrty3Tmz2GNCc4N2UYgIM3Yfx0sbupi3opOv3L+IlZ29fPW0Q/nUWw7QxpUUZvokr1hwybou7pu/muueWM4BU3bhZ+fN1k75Buw/eSwvbuhi6bouvvzA8yxe28Ul/3Ig//aOg5tWL1KOYB2Xruti/opOrv7DK0zepYV7Lj5JKzdugHieHLOZM2eOmj9/frPFMIZHF6/nkrsWAvDhOXvz9fccVve5DCNh844BTvru42QdRULg42/an39/56H5dEYd+e7vX+SGJ18DYK8Jbfzow8doaZUF9GcdZl/5J3oyXvbPGbOmceXpRzR9/kg5bvrba/z3b18EYLexLXz3jCN5x+FR2fN64LiKN3//z6zb7lmUpx62O1edeVTDW4+EEZEFSqk5kT+zymHn5LnV29ilNaXtqbCYRxdv4OnXtnDm7OkcMa35mRyV2N6X5ZrHljF1fCsXnLQvY1r03WQDnni5g0df2MDpx+ylVXwhjv6sw48ee4W2VJKPvXH/uo34rCXzVnTy62fXcsqhUznl0KlNtyKtcrBYLBZLCeWUg752ucVisViahlUOFovFYinBKgeLxWKxlGCVg8VisVhKsMrBYrFYLCVY5WCxWCyWEqxysFgsFksJVjlYLBaLpYRRUQQnIpuAlcP89cnA5oqvGl3Yv3nnwP7NOwcj+Zv3VUpFtmEeFcphJIjI/LgKwdGK/Zt3DuzfvHNQr7/ZupUsFovFUoJVDhaLxWIpwSoHuLHZAjQB+zfvHNi/eeegLn/zTh9zsFgsFksp1nKwWCwWSwlWOVgsFoulhJ1aOYjIu0TkZRFZLiJfbbY89UBE9haRJ0TkRRFZIiKf95+fJCJ/EpFl/lf9JsiPABFJisizIvKI/3h/EXnG/3vvFZHmzWasAyIyUUTuF5GX/M/6pJ3gM/6Cf00vFpFfikjbaPucReQWEekQkcWh5yI/V/G41t/PFonIrJG8906rHEQkCVwHnAbMBD4iIjObK1VdyAFfUkodBpwIXOb/nV8FHldKzQAe9x+PJj4PvBh6/H3gR/7fuxX4eFOkqh/XAI8qpQ4Fjsb720ftZywi04DPAXOUUkcASeBsRt/nfBvwrqLn4j7X04AZ/r+LgetH8sY7rXIAjgeWK6VeU0plgHuA05ssU81RSq1XSi30v+/G2zSm4f2tt/svux14f3MkrD0iMh34V+Am/7EApwD3+y8ZbX/veOAtwM0ASqmMUmobo/gz9kkB7SKSAsYA6xlln7NS6q9AZ9HTcZ/r6cAdyuNpYKKI7Dnc996ZlcM0YHXo8Rr/uVGLiOwHHAs8A+yulFoPngIBpjZPsprzY+DLgOs/3g3YppTK+Y9H22d9ALAJuNV3pd0kImMZxZ+xUmotcDWwCk8pbAcWMLo/54C4z7Wme9rOrBwk4rlRm9crIrsADwCXK6W6mi1PvRCR9wAdSqkF4acjXjqaPusUMAu4Xil1LNDDKHIhReH72U8H9gf2AsbiuVWKGU2fcyVqep3vzMphDbB36PF0YF2TZKkrIpLGUwx3K6Ue9J/eGJic/teOZslXY94IvE9EVuC5Ck/BsyQm+u4HGH2f9RpgjVLqGf/x/XjKYrR+xgCnAq8rpTYppbLAg8AbGN2fc0Dc51rTPW1nVg7zgBl+dkMLXjDr4SbLVHN8f/vNwItKqR+GfvQwcKH//YXAQ42WrR4opa5QSk1XSu2H95n+WSl1LvAEcKb/slHz9wIopTYAq0XkEP+ptwFLGaWfsc8q4EQRGeNf48HfPGo/5xBxn+vDwAV+1tKJwPbA/TQcduoKaRF5N96pMgncopT6TpNFqjki8ibgb8ALFHzwX8OLO9wH7IN3o52llCoOfBmNiJwM/JtS6j0icgCeJTEJeBY4Tyk10Ez5aomIHIMXgG8BXgMuwjv8jdrPWES+DXwYLyPvWeATeD72UfM5i8gvgZPx2nJvBL4J/JqIz9VXkj/By27qBS5SSs0f9nvvzMrBYrFYLNHszG4li8ViscRglYPFYrFYSrDKwWKxWCwlWOVgsVgslhKscrBYLBZLCanKL7FYLAEishteszOAPQAHr3UFQK9S6g1NEcxiqTE2ldViGSYi8i1gh1Lq6mbLYrHUGutWslhqhIjs8L+eLCJPish9IvKKiHxPRM4Vkbki8oKIHOi/boqIPCAi8/x/b2zuX2CxFLDKwWKpD0fjzZQ4EjgfOFgpdTxeFfNn/ddcgzd74Djgg/7PLBYtsDEHi6U+zAv62ojIq8Af/edfAN7qf38qMNPregDAeBEZ58/dsFiailUOFkt9CPfzcUOPXQr3XQI4SSnV10jBLJZqsG4li6V5/BH4TPDAb55nsWiBVQ4WS/P4HDDHHwa/FLik2QJZLAE2ldVisVgsJVjLwWKxWCwlWOVgsVgslhKscrBYLBZLCVY5WCwWi6UEqxwsFovFUoJVDhaLxWIpwSoHi8VisZTw/wE7dyOUDQL46AAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret = run_simulation(100, ndiv=1000, model=hhm, y0={'Vm': -75})\n", "ret.plot(y=\"Vm\", ylabel=\"Membrane Potential\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## FitzHugh–Nagumo Model\n", "\n", "FitzHugh R, Mathematical models of threshold phenomena in the nerve membrane. *Bull. Math. Biophysics*, **17**, 257—278, 1955." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "a = 0.7\n", "b = 0.8\n", "c = 12.5\n", "Iext = 0.5\n", "\n", "with reaction_rules():\n", " ~u > u | -v + u - (u ** 3) / 3 + Iext\n", " ~v > v | (u - b * v + a) / c\n", "\n", "fnm = get_model()" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "v > u + v | (((-v + u) - (pow(u, 3) / 3)) + 0.5)\n", "u > v + u | (((u - (0.8 * v)) + 0.7) / 12.5)\n" ] } ], "source": [ "show(fnm)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "scrolled": true }, "outputs": [], "source": [ "ret = run_simulation(200, ndiv=501, model=fnm)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEICAYAAABBBrPDAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOydd5Qb13m3n4u6vfddLjspdlKkKapLVreqVdxrnOhzi524JI6TOIkdl8QliRw3OW5x3GTLtmRJlqxiWZREUWKn2Hvbzu0F/X5/XMwuuAsMBsAMMITwOwcHu8Ds3PvuzLy9CCklBRRQQAEFFJAIjlxvoIACCiigAHujICgKKKCAAgrQRUFQFFBAAQUUoIuCoCiggAIKKEAXBUFRQAEFFFCALgqCooACCiigAF3kTFAIIWYJIf4ohNgnhNgjhPhonGOEEOI+IcRhIcQuIcSFudhrAQUUUMBrGa4crh0CPi6l3CaEKAe2CiGelFLujTnmJmBh9HUR8K3ouy7q6urknDlzLNhyAQUUUEB+YuvWrX1Syvp43+VMUEgpO4HO6M8jQoh9QCsQKyhuB/5XqqrAl4QQVUKI5ujfJsScOXPYsmWLVVsvoIACCsg7CCFOJPrOFjEKIcQcYA2wedpXrcCpmN9PRz8roIACCiggS8i5oBBClAEPAn8lpRye/nWcP4nbc0QIca8QYosQYktvb6/Z2yyggAIKeM0ip4JCCOFGCYmfSCl/HeeQ08CsmN/bgI5455JS3i+lXCelXFdfH9fNVkABBRRQQBrIZdaTAL4H7JNSfi3BYQ8D74pmP20AhpLFJwoooIACCjAXucx6uhR4J7BbCLEj+tmngXYAKeW3gceANwCHgXHgvTnYZwEFFFDAaxq5zHp6nvgxiNhjJPCh7OyogAIKKKCAeMh5MLuAAgoooAB7oyAoCgDgleP9vHT0bK63kTUc6Brhqb3dud5G1tAxOMHDOzt4rQwqGxgL8NCOM4Qjrw16h31BHtvdiT8UtuT8uYxRFGATjAdC3PPtTQD84D2v4+oLGnK8I2vhC4a585svMBYI819vWc3tq/O7NCcYjnDnN1+ka9jH4HiAd108J9dbshxvvn8TB7tHOd43zkevXZjr7ViOD/zfVl44fJb3XDKHf75tmennL1gUBfBkjGb9b4/vz3ut85Xj/YwFlOb1X08dIpLnWufR3jG6hn0A3Pf0YXxBa7ROu2BgLMDB7lEAvvPcEYYmgjnekbUIhSO8cmwAgD8d7GU8EDJ9jYKgKIBT/eMAfO6O5ezvGuHlY/053pG1ODMwAcCnbrqAo31jPLkvv11QJ6PX95M3LKZv1M/DO+KWIuUNjvQqIfHJGxYzHgjz85dP5nhH1uJo3xiBcIQv3rmCRz9yGSUe8x1FBUGhg3/47W4u+7dn2HlqMNdbsRTdw36qStzcfWEb5UUufrI5vx+sjiEfDgHvuWQObdXFfP/5Y7nekqU4cXYMgLeub2dJcwXfe/5YXluNmqC4dWULF8+r5UcvHs/rWMW+TtXQYvWsKkuEBBQERUKM+IL8dPNJTg9M8LcP7spr90T3sI/G8iKKPU7uurCNx1/ton8skOttWYaOwQkaK4oocjt56/p2Nh/rn2Sm+YhT/eOUe11Ul7h5x4Z2DnSPsOv0UK63ZRmO9o3hcTporS7mXRfPpmPIx8ZD+dvWpzvqVmyrLrZsjYKgSIBNR84SkUoL2981wsbDfbnekmXoHvbRWFkEwNsuaicQjvCrraeS/NX5i47BCVqq1EN154WtOAT8auvpHO/KOpwZnKCtpgQhBLeuasHrcvDAlvy9voNjQapL3TgdgmuWNFJT6uGXW/L3+g5PhHAIKLXImoCCoEiIziElpf/y9QuoK/Pyvy8ez+2GLET3sJ/Gci8AixrLed2can7+yqm8dU90DftoqlCCsbmymCsW1fOrrafz1j0x7AtRWayYSEWRmzesaObhnR15G9QeDYQo8yp6PS4Ht69u4cm93QzkqZU84gtS5nXhcOjWL2eEgqBIgJ4RH06HoKmiiLetn8UzB3omg775BCklvaN+6qOCAuDutW0c7R1jZ566J8b8IcqLprSvN6+bRWceuyfG/KFztM171rYx4gvxxJ6uHO7KOoz6QpQVuSd/v2ftLALhCA/tOJPDXVmHEV+IimJ38gMzQEFQJED3sJ+Gci8Oh+At69sBeHBb/pmv/lCEcERSFsM4b1rRjNfl4Nd5SC/AuD9Mscc5+Xu+uyfGA2FKvVPXd8O8Wtqqi/PW/TTqD1EeQ+/SlgqWt1bwQJ5e32FfkPKigqDICbqHfTREteyWqmIumV/Lr7edyTt3zES0nqDEPcU4K4rcXL+siYd3dhAIRXK1NUsgpWQscK6G7XE5uGN1K3/Y25WX7olRf+gcQeFwCO5ZO4sXDp/NSyt51DfletJwz9pZ7O0c5tUz+WclD/tCVBRZWztdEBQJ0DvipyHqxwa4c00bJ/vH2XJiIIe7Mh/jUT91rIYNKsg7OB7kjwd6crEty+APRYjI+PQGw5LHXs2/LvbK9XQuvXetVdXo+eiOGfWHzrGQAW5f3YLbKfjN9vyjd3iiYFHkDL0j5/rtb1zeRInHyYN5lh0zEa3iLJ6WMXH5gjrqy715534aj1pQ0xnnspYKFjSU8dD2/CpGi0TkDNcTQFt1Ca+bU81vd+Rf/yctuBuLqhIPVy1u4Hc7O/IuaUHFKAoWRU4wOi3gWep1cdPyZh7d1ZlX2SITAeVainU9AbicDu5Y3cIz+3vyyh2jtTeYXpgkhOCO1S28fLyf0wP5447RLMbpjBPgttWtHO4ZZW/n9AnE5y+klMqiiEPvHatb6Rnx513zy2FfkIqCRZF9hCMSfyhCifvcm+2uC1sZ8Yf4Qx51HR2ftCicM76788I2gmHJI7vyR8vWLIoS70x6teaAD+/MH3rH/FHBGIfem1c043KIvGrpMREME5HMcD0BXLOkgVKPM6/cbZGInKHUWoGCoIgDjXmWTnu4NsyrpaWyiN/kkTsmUYwCYElzBUuaK3hwW/48WBrjjFecNKumhLWzq/PK/aTRG0/Drin1cOWieh7e2ZE3nQdGfYnpLXI7uWF5E79/tStvvAL+UAQpZ1rIZiOngkII8X0hRI8Q4tUE318lhBgSQuyIvj6TjX1pWud05ulwCG5Z1cLGQ30MjedHR0qfRqt7pqAAuG1VCztODeZNdsxEgmur4fbVLRzoHpnsn3O+Y8yvxWTiM5LbVrfQOeTj5eP50QhyREcwgnI/jfhCPJsnSRpaVqLHZS0rz7VF8UPgxiTHbJRSro6+PpuFPelqnbesbCYUkXlTrDTpiknAOG9Z2QzAY7vzIxtoLKDPOG9e0YzTIXgoT9wxozquJ4DrljZSkkfuGE0RKEqg+Fwyv5a6Mk/eXF9/WNGb14JCSvkcYDtVJpFFAbCitZL2mhJ+lyd+ez3XEyh3zKq2Sh7NE0GhF5MBqC3zcsXCOh7ecSYv3DEavYk07BKPi+uXNvLY7i7LpqNlE4Gw0rC9CRiny+nglpUtPL2/h2Hf+e8V0CwKrzOPBYVBXCyE2CmE+L0QIuHoJiHEvUKILUKILb29mbViGNfROoUQ3LyymRePnM2LDquT6bEJNDCAm1c2s+v0UF50WJ28tgk0bFBB7Y4hH9tOnv81M76gYiSJNGyA29e0MjQRZOPB87/xpRFXzO2rWwiEIjzx6vnvFXituJ6SYRswW0q5Cvg68NtEB0op75dSrpNSrquvr89o0WRa5y0rmwlHJI/nwY02mR6rEwx7wwrlfsoHq2LS1eZOTO81SxrwOB154V4MRjVst47Geen8OsqLXDyeB/QaYZyrZ1XRWlWcF9dXs6Be04JCSjkspRyN/vwY4BZC1Fm9bjKtc2lzBXPrSvMibXQ8GMLjcuDU6TzZVl3CmvYqHtl5/gsKLdvF605865cXublkQS2P7+k674vRjDBOj8vBtUsaeWpfN6Hw+d2yZZJeHcEohOD6ZY08d6hvMh55vsIIvWbA1oJCCNEkhBDRn9ej9mt5tUwyrVMIwS0rm3np6Fn6Rv1Wb8dSTATCum4nDTevaGZv5zBHo9PDzlcYfbBuXNbEqf4J9nWOZGNblsEfNkbvDcsaGRwPnvdjcI1q2DcsayIQivCng+d3x+DXhOtJCPEzYBOwWAhxWgjxPiHE+4UQ748ecjfwqhBiJ3Af8BaZBRVvsnpXx499w7ImIhKe2X9+p9n5gxGKdLRrDTdF3U9PnefzpYPhCC6HSNq7/9qljTgE5707xqhgvGJRPV7X+e9uM8o4XzenhppSz2uG3kyR66ynt0opm6WUbillm5Tye1LKb0spvx39/r+llMuklKuklBuklC9mY1/Jcs9B9QZqriziqURV2sOd8NCH4Ee3wqEnrdimKQiEI4ZustaqYpY0V/DU3gSCcbwfHvkY/PAW2PMbk3dpHgIhY/TWlXlZN6eGPyRiJP4RePzT8IObYfv/mbxL8xA0qGGXeFxcsaieP+ztju9uC07A05+DH7wBNn8HbOqSCxi0oJwOwbVLGnhmf0/8DsnhIDz3FUXvxq9CxJ4uOX8hRpE7TARCCIGupi2E4NoljWw81DezynPoNPzPtbD7Qeg7DD9/O3TvtXjX6SEQihj2b163pIEtJ/pn9n6aGITvXQ/bfwwDx+FXfwanXjF/syYgGI7oBnZjccOyJvZ3jXC8b1q2V3BCKQCbvw3Dp5VCcOgpC3abOVLROG9Y1kTnkG/mPO2gD/7vbtj4FRjphN//Dez+lRXbzRiT6aIG6R3xhdg0vfdTOAQPvAue+RyMdMHTn4VXvmvFdjNGIUaRQ0wEw3hdDqLhkYS4dmkjE8EwLx6JSSuMhOGX7wXfELzvCXj/RnC6YdN/W7zr9BBIgXFeu7SRiOTc1uNSwm8/CAPH4B2/hg+8CKX1SguzIYxaUKD89sBM98SjH4eO7fDmH8OHXoaaefCnL5m9VVMQCEVwCHSTFTRcu6QBp0PMdLc98zk48Tzc+V348FZoXgXPfsGWVkUqgvHSBXWUeJwzr++L98GBx+Cmf4e/3Apzr4DnvqysDJshFcGYCQqCIg4CoQheV/IA74Z5NZR5XTwZ647Z8n04/TLc/BX1QJXWwfI7lTsm6LNw1+lB0WrsNljeUkljhffcOMXuX8KBR+G6z8Lcy6GoAla+GQ4/CWP269LpT8GCaqsuYUlzxbmC8fBTsOMncMUn4YKbweWFNe+E06/AwAmLdp0+gikIxqoSD+vn1PB07PU99TJs+gasfS+sfBM4HLDuz6D/KHTH7byTU6SSLlrkdnLV4nqeinW39eyDZ78IS2+H9feCELD+/8FYL5zIiuc7JbwmYhR2hd8g8/S6nFy5qJ6n90VvtIlB5cedd5VilhoW3QTBcTizxbI9p4tUXDEOh+CaJY386UCvukGDE/DUvyiBeNEHpg5ccitEQnDiBYt2nT6CYZnSQ3XFojq2nhhQaZThEDzx91A9VwkKDUtuU+9HnjZ5t5kjFcEIcOXieg52j9I15FMWwxN/D+XNcP3npg5afLN6P/i4ybvNHKm6Yq5a1EDPiJ+D3dFsvqf+BdzFcPPXlJAAmH81uIrg4BNWbDkjFOoocgijAU9QjKRnxM+hnlHls/YPwfX/OnWTAcy+GBBw3H6MMxVaAa5aVM9YIMz2kwPwyveUj/76zytNU0PzavVgndxkwY4zQyAUTo1xLqwnGJZsOnJWWU+9+5X15JoaakXtfCipU9q3zaBcbcmtYw1XLFTFqhsP9Srr6fTLcOUnwVs+dVBZPdQttmUcSnO1uQxe48sXqbKsjYd6FT0Hfw+XfER5AjR4SqHlQvW/sBkKMYocwqhFAcrPCbB533F46ZtwwS3QtOLcg4qrof4COLPV5J1mjlQsCoAN82txCNh0sAte+hbMuVy5nGLh8kDrWlsyzmBY4nYl99drWDunmmK3k+cO9ijfdcNSZTHFQgho3wCnNpu828yhkhWM03tBUzl1ZV6eO9Sn/PJV7bD6HTMPbHudcrfZLE6RSgwKoLmymIUNZaqeYuNXoaQWLnr/zAPb1kLnLgjZq21PwfWUQ/hDxrWwtuoS5taV4tj1cxXAvuxj8Q9sWg7de0zcpTnwp2hRVBS5WTWrCrn3N8qauOQv4x/YtEL5e22WVphKlhco9+LF82vx7XsCevYqbTNekkPzauW399urQC+VGAUo9+IVC+s4e/AlJfg2fFAJ/ulovRAm+mHQXnGZVK8vwOUL6+k8thd58HFY9z7wls08qHUthP22i8sUXE85hD8UTimL4NL5NWw4+xsiLRcqzSMeGpcrxjphr0ZzwXDqD9al8+u4cvC3hGsWwILr4h/UsBSCYzB4PPNNmohUXW0AVyys47rxRwmXNMDyu+If1LBEvfceyHCH5iItehfVc2foMcKuElj9tvgHNa1U7zZTflJR8jRcvqiOt/AEUjhVoD4ebEwvFFxPOUGqD9etVceYL85wfO5bEx/UuFy926yeIlVTHeDa+kEudBziSNtd58YmYtEYbfRrQ3pTcbUBXN0a5vWO7exrvDW+dg3QuFS924yRBEKp03v5LDe3Ojaxr/4mKKqMf1DDEkBAl8007BTcxhoumlXKXc6N7K+6Eiqa4x9UPQfcJba8vh5n8lT+TFEQFHGQSowCYE3/44zIYh4Nb0h80CTjtNeDFQzJlLWR5T0PE5ROHpKXJz6o/gL13mMzQZEGI2k/+RBOIfmtuDrxQVVzwF1qP3rTUARqT/werwjyQPiqxAd5y6BmLnTvzmyDJiMdekuOP021GOUXoSsSH+RwKuFos+c3HYsxHRQERRykxEyCE3gOPsLLRZey8bjOvIbyJiiusd+NFo6kFNwlEsH16i/ZWXIRT57UiT94y5QWZjMNLNXgPYDY/QCHi1fy6JmSxN1kHQ5ouMB+giINnz07f06vt50Hu+r1u8k22i/ulmpWGwA7f86ou5af9s2bnAgYFxq9NgrgB8LhgqDIFVSMwqCf8+Dj4B+mc/at7Dw9GL9vDKgAaNNyW5rqHmcKPt1Tm2G0m7Ozb+Zg96j+7PCGZfZjnKlqnD37oXc/fbNvpnPIx+mBicTHNiyxpastJXoHT8LJF+mbdwdjgQj7u3SC800roP8Y+O3TUThlDXtiAA49weD82wlKp0r7ToTG5SqAP2KfRoJpKQJpoCAo4iClh2vPb6Cskepl1+IPRdjXOZz42IZlKg/fRplAKVsUex8Cp5fKVbcAsP2U3oO1FM4esVVFesoP1r6HAUHtujsBeOW4ThvuhmUw3gej9ukonHKywv5HAai5SMXbdNuONy4DpK2UgZQF48E/QCRE9evehEPAK8f17mfNfWwfK8o2richxD1CiPLoz/8ghPi1EOJCy3eWQ/iDBl1PIT8cfhoW38SFc2oB2HoiyY0WHFd9kWwAKaVysxllJJGIYpwLrmXFvDYcAradHEx8fMNSkGHos08mUDAcwZ3Kg7X3IWjfwPx5C6kocukLChsGtFNmJPsfhYalNM5ZSmtVcRJ6owkaXfaJU6SsCOx/BMqaKJ17EUuaK9hi6PraxyuQTkwmHRhZ4R+llCNCiMuAG4AfAd+ydlu5heF//rGNEBiFxTfTXFlMS2WR/pxlmzGSUET5Wg377M9sheEzsPQ2Sr0uFjdVJDHV7Zf5lFJLi77DiiksvR2HQ7BqVtXMzqqxaIjS27Mv842ahJSynsb7VduVxW8A4MLZ1ew8paMIVLWDt9I29zOkKBiDE0rRu+AN4HCwbnY1O04NEo4kiEEUV0NFm/3otYnrSeuhfTPwLSnlQ0CCHMH8gGGL4sBjKmVursqWuHB2Ndv0LIr6aEqhTUz1lKs69z0EDjcsuhGAC9ur2HFykEiiB6tmPji90GOfByulArR9D6n3aCX28tZKDnaPzGwrr6GsXnXOtRG9KTHOg4+DjKhmh8CK1go6hnycTTTFUQilDNhIw06pgPTos6rW5wLlRl3eWsl4IMyx6W3lY9G4zFaCItWC2XRhZIUzQojvAG8CHhNCeA3+3XkLQxaFlHDg9zD/9eAuAmBVWxUdQz76p89r0OApiaYU2uNG04baGNY4D/4B5lwGxVUArGmvZsQf4miiB8vpgvrFtrIoUtLADj2lCq0q2wBY2VpJMCw5oBfgbVhqm+sLEEilCeL+R6G8BVrWAIpxAuw+o2NFaYzTJnG3lFwx+x9RFtEclea9ok3R+2oyevsOKLezDWCbGAVKQDwB3CilHARqgE/q/4kxCCG+L4ToEULEVUmEwn1CiMNCiF3ZiI2EwhHCEZk866lzJ4x0TJrpAEuaKwD0A9o20khSsiiGTqsHZME1kx8taVaN4vTptU8KZSgcISIN0usbVk3gYug1xjiXq0ypSAKrI8swnC4a8sORZ5QbJlq8pdGryzibliv3q00q8A0rAlIqxWfhtZNFlAvqy/C6HMkFRSQEfQdN2nFmCIRTrwtKB0lXkFKOAz3AZdGPQsAhk9b/IXCjzvc3AQujr3vJQmzEcO8UbbzpwusnP9IY596OJJlP/UchMJ7RPs2A0bGRABz5o3qfP8U4FzSU4XKI5IJxtAvG+hIfkyUEwynEZI49pxhCDL1t1cVUl7jZrRenaFwGoQl1jW0Awxr2qc0q0WLBtZMfVRS5mVNbkkQwRhtg2kQZMCwoevbCWI/yCEThcjpY0lzBqx1JFAE4/+jNEEaynv4J+Fvg76IfuQFThgRLKZ8DdNIMuB34X6nwElAlhEhQY28O/EGDE6OO/UnlkZfVT35UW+alscKbhHEuBaRKk80xUrIojjyt5hJoPY1QDfPm15clFxRgiwcrNXqfAU8ZzLpo8iMhBMtbK/UZZ5PGSOzhtw+GDVbeH30WhBNmX3rOx8tbK3n1jJ7iY69WHoYF49Fn1fu8q875eHlrBXvODCeOu9UuUHE3m1zfQCiCV2dks1kwssIbgduAMQApZQdQrvsX5qEVOBXz++noZzMghLhXCLFFCLGlt7c37QUNWRSBcaWBzb1yxldLmyvYm8wVA7ZgnJqGnfTBioTVgzX/9TM6py5pLmdfp47P3kb0TllQBupGjjytkhSm9Xa6oKmcI72jiTNj6hYrhmsDesMRSThiMEZx9E/Qtk5NKIzB0pYKzgxOMOJLUFjpKVHzOGzEOI0pAn+E2oWT8ScNy1sqGfGHEhdWOl2qAt8G1xeigtEOFgUQkKpvgQQQQpRau6VzEO+JjvuESinvl1Kuk1Kuq6+vj3eIIUxZFDoxilObIRyYoY2AilMc7hnFH0rgo66eA65iW2Q+aRp2UldMxw5VwRpjpmtY0lxB17CPgUQB/LJ6KG2wBSMx7FY8ewQGjseld359Gf5QhI7BBIzEXQR1C22hYRu+vhOD0LEt7v08v1613D7aq5cJtNw2tRSGBEUooNKA510146sFDYreI7061eY2irvZKZj9QDTrqUoI8RfAU8B3rd3WJE4Ds2J+bwM6rFwwEFYMXveff+xP4HBB+8UzvrqguYJQRCZ+sBxO22gkhmgFOPqMep931YyvJgP4XTpWVNNyewgKo66nY39S73EEhcZIDvfYn5EYpvf48yotdt5VM77SBIUu42xaruZS+HTugSwgEpGEIgZcbadfUfGYeVfN+GqeEXobl8FoN4ym77kwC7YRFFLKrwC/Ah4EFgOfkVJ+3eqNRfEw8K5o9tMGYEhK2Wnlgj4jMYpjz0HrurgDTubVKYPrxFn752IHQlpwN4kr5sQmlfYZOx4yivlRxpk097xnv5o5nUMYTgc+8SKUNUHNvBlfGWKcjctg6KQaZJVDGLagjj6rOt+2rpvx1ezaElwOkYTeaEA7x1ZySvQKx8zJjEBNqYeqEnfilG+YirvZoF4m5V5tacKQKJJSPiml/KSU8hNSyifNWlwI8TNgE7BYCHFaCPE+IcT7hRDaLMLHgKPAYZQV80Gz1k6EpDfbxCB0bJ8sspuO9toSAI6f1clqsklPII1WXaEYCauRpnGsJ4CmiiI8Tgcn9ehtXK6mg/UfyWS7GcPwfOETm9Sc8zg9/qtLPdSWepJbFJBzZcBwTObYc4reOLM23E4H7bUlSVxPUcaZY/eTNsQnaSLK8efVRMIEszbm1ZVyNJnrCXJ+fQH8WWrh4Ur0hRBiBBUPEJwbFxCAlFJWxP3DFCCl1Jn0oxYBPpTpOqkgadbTqc3KTE8gKCqK3NSWepJbFKAerJg8/WwjaMSH3bUbAiMw+5K4Xzsdglk1xRw3Sm/94nS3mzE0xqnb62nwpJpE2P7RhIfMry9L7ooBxUgS/N+yAUOup/F+VR+z8k0JD0lKb2WbYro5di8aojcUUPGYde9LeMi8+jKeO6jjViqtg7LGnAsKrVdbTl1PUspyKWVFzHtF7O+W7yxHSKpln9qsslpaE4w8RZnrx/t0NOymqKneuTPdbZoCQ6b6yU3qvT3xUKbZtaWc0LMo6haD0wNdu9LZpmnQGIluE8QTUXpnx7egQLnbjuhp2OXNavZIjq9vcNKi0EvMeFm961zf+fVlHO8bTzybQgjlfsox4wwaqQvq3AkhH7RflPCQefWl9Iz4E2d6gVJ+cmxBaVmLtii4E0K0x3tZvrMcwR/t45Mw6+nUy9C8UqUFJsCc2lJ9i6KkBqpmQ+eOTLaaMQz57E+8CJXtM9IIY9FeU8LJ/vHEQ31cHvVgdeSWXkMa58kXVVuHhqUJD5lfX0r/WIDB8QSZXkJAy+qcX9+prCcd19OpzSoxoyVx04N59aUEwhH9WRyNy1Srlhy28jB0fU+9pN5nJRaM8+oMxt16cxt3S6lgNkMYWeHRmNfTqJjB763cVC6hq2WHQ6qDatt63XPMri2lY8iXuHkcqH46OWacSQezSwknX9LVrkFZUOOBMH2jCRgnKJ9w546cTgczJhg3Ke3akVgLb6tWSoIu42xerbrI5nAWh98Q49ys+lklUXwATvbrWcnLVYO9HLbQN2Yhv6RS1MsbEx4yv17Rqy8olqsU+bOH09mqKUi5qWcGMJL1tCLmtRBYDzxv+c5yhIBeQKz7VZVWN0tfUMypUw/dKb0Hq2W1Sikc1ytMtxbBZA9W/1HV5kDHLQFKUECSTK+WNSoLKJeMJB8H1rkAACAASURBVNmDNXZW+euTCMa26mIATg/oXd81qgVIDt0xSa9vKKAUnyTXtzVK75lEtSNgi9kUSZMVpFSCUceaAKP0ah0HcheXsZWgmA4p5TbgdRbsxRbQ1cI0f+6sxP5NUBYFJMl8al6t3nPonkj6YJ3arN4TZDxp0OjVjVO0ROnt2J7SHs1EUo3zzBb1nuT6zjJiUUQ7sNKxLaU9momk17drt/LXJ1F8Gsu9OB1CXzA2LFEppzkUjEktqP6jMNarG58AKPG4qC5xc0bv+mpxNzs/vyYiYdaTBiHEx2J+dQAXArmvNLEIUxZFHNfDqc0qUKnjr4cpjfOMrsapMc4dcQu7soGkGueZbarfUd0i3fO0VRcjBJzQs6Dql6geOR07YPld6W45IyR9sM5sU8yueZXueSqKXZR7XfqCorINSmrtwUgSJmYk99eDapbXVFGkzzjdxaoPkp01bE3xSUIvKKsiYfU9qLhb0wo4k0vFx2DBrAkwskJ5zMuLilXcbuWmcgmt9Ubcf/7pl5X2FSe/PhY1JR48Tgddwzo964urla/UBowkoc++Y5uyfHT89aCEal2Zl+4hHX/8ZEDbzhbFVqi/ADz6XWqEELRUFesLCiFyHodKSu/pV6ByFlQk77PZWl2s74qBnLfySJqxeHoLeMrVNU6ClspiOgaTxJda16r7OUct5Q3FoEyCkRjFv8S8Pi+l/ImUMncROouRMEYx0q1y7JMEsgEcDkFDhZeuoSQPVvPq3DISvayYUEA99K1rDJ2rscJL90iS26JljUpPzFFmjG7diJRKMLYaG3nSUOGl1wi9Pfty1lI+afC+Y/uUiywJmiuL6Bo2QO/QqZwVkk5ZjAkUm84dypJ3JGesLVVJLApQlezBsZx1grZVjEII8aQQoirm92ohxBPWbit38IciCAEuxzTmqeXEG3ywmioMPlg5DGgHoi2oRTwLqWePyurQSZuMRVNFEV16FgWoh9Q/nLNZDboa9uBJGD9rmN7GiiK69SxGUIqADOfMHaOb1TYxoBofGryfG8q99Az7E6dAA7RFQ5ent6S4U3OgyzjDQdWoMYlbUUNjRREj/hDjAZ30V62W6szWVLdqCgzVBZkEIyvURyfbASClHAAarNtSbhEIqYlRM5in5iLSiuWSoLHSACNpi/bWyeGDpRufgBQ07CJ6RpLQm+MHa7KterwHS9uTQXobK7z0jvoTtxuPPVeuGUm8a6xZslqsLAkayovwhyIM+3QYZ8tqVZNx+pVUt2oKdH32PftUGxmDgrG+3AtA34hOynfNPFWRnqvra7S3lQkwskI4tsBOCDGbBK2+8wH+RBOjOneqYF2RsaJ0Yxr2GlXlrQXZsoxgOJK4GOvMNlVdXDXb0LmaKoroHwskbq8OyjfsrcgZvX49V1vHNpXF0rDM0LkaK4oIR2Ti+egAFS1Q0TaVLZdl6CYraIpPs0FBUaEYp667zV2sFKlcCQo9i0KLjaVgQQH06NHrcCgLNMeKgF0Exd8DzwshfiyE+DHwHFPT7vIO/lAErzuOj7Njh+GHCpRGMhEM65uunlL1YOWIcepaFB3blAWQJHCvoTHKSHr0rChHtPVJjhinNjYyrqvtzHZVeBanMV48aIykO5l7cdb6nNILCWIUHTuUElBSY+hcDeVFQJLrC8r9dGZbTgK8ulltnTtUxX31XEPn0iyK3mRWcts6VZGegziUrQSFlPJxVErsL4AHgLVSyjyOUcQZRj/WpxrFGfRvgmpXDHBWr1oZFCM5sy0nrQCURRHnFghEA3QG3TAANaXqwdLVsEHVKPTsycnsgoQWVCSiGEkK9NZHGWdSRjLrInXvDJ1OZaumQJeRaIFdg9AsiqTuxbbXqQBvzz7D5zYLullAHduhZZWhQDZMCQpD7lQZzkn2oq1aeAilft0IXCil/B1QIoRInvpzniLuDNrO1Py5ALWaoDDCOINjOeltn7BFcfce1SE3BQtKE4zJBcV6de4cmOsJLaiBYxAYVRaFQWjXdyBRvycNWjFbDqyKQDhBYsZ4vwpkp3B968qiPvtRo3G37LufEqbHhgLqnk7lfi7x4HQIfdcTTAXwT76UylZNga3SY4FvAhcDWkvwEeAblu0ox4gbo9ACf6kwkjJNwzaggUFOGEkwUTxGy4XX2mUbQK1RQdG2DhC5oTeRYEyD3uoSg/Q2rVCjb3MkKNzxXG0pZvABVBS5cDpEcsFYPVcVGuZCUCRyPfXsjWbwGafX4RBUl3joH9PpIAuq5Xjdoqkuy1mErVxPwEVSyg8BPpjMejLmyD0PEYgXo+jcoR6A4qr4fxQHGuPUbZQHUNWupqnliJHEtyheVf7cylkzv0uAmjKDjLOoUrV7OJ0DekMJXG3dr6qkgvolhs9VbpRxOt3RuEz241D+YCR+xpMmKFJwpQohqC5xMzCehHEKoazkHGjYgVAEl0PgmG5BaYpACvQCVJe4E3cIjkX7xXByc9bjMlPpsfaYcBcUQjiJZjoJIeqB3PUSthj+UHhmXnLnzpRvMsOuGCFg1uum2ilkEcFwAouie4+qojYYyAYo97pwO0VyVxtEA7yvZL3wLqGrretVqFsI7iLD51Iapzu5xgmqt1DXrqwHPAPhSPxWNN17oKLVcCBbQ1WJxxjjnH2JmmY4bOnU4hlI6Frs3gPuEsOBbA1VJe7kigAoev1DWR8Fa7f02PuA3wANQojPozrHfsGMxYUQNwohDgghDgshPhXn+/cIIXqFEDuirz83Y109+ILTYhT+EVWMlYJbAqDE48TrciQXFADtl6g1Bk+luNvMEFfDjkTUg5UivUrj9DBgiN6L1YOV5UK0QCJXW/erU91PU4BhemddpDrJnsluvUxCi6Jn71T30xSgBKMBeudcpt5PvJDyGplA10JuWGI4kK1BCUYjikC0aeaJ7LqfbOV6klL+BPgb4ItAJ3CHlPKXmS4ctVK+AdwELAXeKoSINy3mF1LK1dHX/2S6bjL4gmGKYl1PvQfUu84gm3gQQlBT6jH2YGlD3o9vTGmNTBFXAxs8rgK7aTDOmlKPMYtCYyQ5oHeGW3FiQLWdSFEwgpqf3W/INbFBNRs8lmV6w3EERTio7ukU72dQgtEQ42xaqepljmd3GkFcRUDKKQs5RSjXkwF6q9qVhXbyxZTXyASBUASnQ+Cc7mqzAAkFhRCiRnsBPcDPgJ8C3dHPMsV64LCU8qiUMgD8HBs0G/SHIucKCs2cbDDuv9ZQUeTWH6eooWGZKm7L9oMVljMtCq1NdDqMs8TD0IQBxlnZptwAWWac/lB4JuPU6G00VnEfixqjrpiiShVIPfZcymtkAn8wPFMR6DsEkWD6FpQReh1OJRyzbVHEs5BHu2GiPy16q4zSK4SyKk5syupgrkAi17EF0FtlK7Al+j79ZYYN3QrE+lpORz+bjruEELuEEL8SQiSMrgoh7hVCbBFCbOntTb8Lui8Ypij24erZp7JWquakfK6KYhfDEwbqIxwOmHNpbhjn9FTgrleV9ptCYFdDeZGLEb0WD7GYe4Uas5rFAKA/FEfD7oq6v9IQjFVGgrsa5lyuXE8BneFOJiMuvZri05i6RVFV6mZgLKjf70nD7Euh72BWGwT641lQmiKQhgVVVeLGH4roT6rUMPsSGO3Kah8z3YJZk5FwFSnlXCnlvOj79Nc8E9aOZy9NvwN/B8yRUq4EngJ+pLPf+6WU66SU6+rr69Pe1AzXU/ceaLggZf8mQHmRm2EjFgXAnCtg6KTKb88S4vqwu1+Fmvm6ozETobzInZqg8A9NZeBkAYreaa6n7t1QUgdliUdjJkJ5kYvRVOiNhLKaRhmXkXTvUf2YahemfL6qYg+BcIQJI4xzTtSdmkWrIiG9kKbryWCtDMC8q9T7kWdSXidd+O0gKGIhhLhNCPGV6OsWk9Y+DcRaCG1AR+wBUsqzUkqtEOG7wFqT1k4IXzBCkXuaRWGw/890VKSkYUcfrCxaFTPcbKBSCdPQrkExzuEJo4IxGqfIojsmoeupaXlKGV4ayovcTATDkz2VdNG+ARzuHNA7XTDuUdPZDLYqiUVZkZpzNuo3cE83r1JDr7JIb1xB0bNXDRtLMcMLoLLYDWAsTlEzDyrb4eizKa+TLhImZ1gAI5XZXwI+CuyNvj4qhPiiCWu/AiwUQswVQniAtwAPT1s7dqLKbYClfQGklPhCMRbFWJ+aGZ1GfAKgojgFi6L+Aiitz2qAdwbjDIyptufpCsZiN6OBEBG9jqoayptUoVJW6Z1mQUUi0HswLTcbQJlXMc4xI4zTU6qKDbOoCMTNAurZm5bbCZTiAxizopwuZUUdfjprfvu4jLP71bSsCUjx+goB869SgjFL7XjiJitYBCOrvAG4Tkr5fSnl91HtPG7OdGEpZQj4MPAESgA8IKXcI4T4rBDituhhHxFC7BFC7AQ+Arwn03X14A9FkJIpQaH1q0lXUBS5GZ4w6NMVQj1YR57JWn2BYpwxGmffIfVerz/6NBEqilxICaN6jRBjMfdKFacIJaleNwn+6e1Zhs+o9ilp0lseZZwpuds6d2Rt/sgM16JvSGV4peGvhynGaciiADXid/BE1vz2MwRjOJR2hhdMWVAjRumdd7Wat5KlOemBUJxkBYtgdJXYkuRKsxaXUj4mpVwkpZwvpfx89LPPSCkfjv78d1LKZVLKVVLKq6WUlo6S8gcVg54pKNLUwIpdRCSMBQwGbBder4a/d1nvt49E5OTsjUlogiLJjOxESJlxLrwOguNZ82P7g9NcMX3R1Oe06VWuCcNW44JrVZ+rLPmxZ2ic2v2coYZt+PouuFa9H34qrfVSxQzXU/8R1bojGxYFROMUAo78Ma31UoUtgtkx+CKwXQjxQyHEj1BZT6YU3NkNvugshckYRc9eKKpSbpI0UKExEqN++/nXAAIOPZnWeqlgsoFarIbdd0C1sqiZn9Y5NcZpKCUYVMDT6c0KvRCHcU4KxsVpna8iVcHYulalQWeJXn9wGiPRRnYamBkdD2Wp0lszV/nuDz+d1nqpYobrqe+gek9TEZi0oIzSW1KjYjNHsyQobJIeC4CU8mfABuDX0dfFUsqfW72xXEBLgytyxVgUDUvTCnTCFOMcMiooyupVvn0WGIlmPZ2jYfceUA93GoFOiBWMBh8sT4kK4h/6Q1rrpQIp5cwYRe8BKK5Wjd3SQFkqPntQ9QULroXDT2bFvTijhUffIZXqnUIPr1iUe9X1Nex6AkXv8Y1ZcS/OaPqYoaAoTdXVBrDgGtW3bWIgrTVTwQyLYqTbMremXsHdhdoLaEZlKZ0CWqKf5R18sa4nKaOCIr34BECpVz2kusOLpmPh9arz5tjZtNc1Av906wnUg5Wmdg2xrieDghEUvWcPw9kjaa9rBMGwRErOrczW6M1QERjxp0jv+NmpiWsWYkbBXd9BqFuQVqo3xArGFOidf03UvWh91fKMdNHeg6pi2luW1vlSjskALLpJzac4ZL27bYageO7f4b7VliQP6N0xW1B1C1+Jvr4a8/qK6TuxASYtCrcDhjtUnn9GgkLdaONGYxSgGAnScj+21st+UuMMhxSzrks9v15DSumTGrLkx9YE4zmmeu+BtAPZkEZMBpTGKRxZsaJmWFB9B9PWrmFK8Unp+s69HFxFcOD3aa9rFDNdi5nR63QIit1O4xYjKPdiaQMceDTtdY1ixkgEjd40FR896AmKjwNDwATwA+DWaED5ainl603fiQ0wJSicGQeyAYrdmkWRgqBoWaMKwA48lva6RqAxzskHa+C4au1Qn75FUeJJg97a+WoWucWMZFIwahbUeD+M92XESNISFCU1agbJwcfTXtcIwhFJKCKnNM7gBAycyMhi9LqceFwO41lAoNKC518D+x+xPE32nBiFlMrVlsH1BaX8jKXiEXA4YPGNyqIIGSjUywAzsrxMoDcR9Cqz/0NKeRkqhXUW8LQQ4gEhhPExUecZfCHN9eTIqMeThinGmeKNtuQWpXEGJ9JeOxl8kzGK6C0wmQGUgaBwK8Y5kYqgAFhyq8o/tzBtdMqCitLbmzm9HqcDp0OkTu/iN6g02YETaa+dDIHpFuPZI4DMyGIE1U4+JQ0b1P08fMbytNFzXDEjnRAYMYXelBQBgMU3q7VPWNu7zRcITyqj+IYVzRnSmwhGgtnHgIeAP6Aa+VkjsmwAzaLwuqIWRVlTWhWdGko8abieAJberjq4Wuh+mrQotBttknEuSPucRR51Oxlq8RCLpbcrv+5+68x1f+y1halAZwauJyEEJW5n6td32R3qfe9Daa+dDDNaUGcY2NVQ6nWl5noCWHSjyqbb90hGayfDORq2dj9nYCGDotdweqyGeVeq+Rf7rfUK+GI7K5zNLLU9GfSC2fOEEJ8WQmwG/gXYCVwgpXzAkp3YAOe6ntKvYNVQHLUoUtY451yusnEsZCT+GRbFQdXqoCj9MhlNw07JggI1y7iq3VrGOX2ect/BaAZQe0bnLfI4UxeM1XMUzVZe3+muxb6DgFBuvgxQ7HYaa5IXi5Ia1bJl3+8yWlsPwXCEcEROadgZ1gRpKEtHMLqLVSxq38OWNr2cCIQneYxZ9CaCnkVxGHgT8DiwCWgHPiiE+JgQ4mOW7CbH0DTDMjcq5zyD+ASk6bMHNT7zgpuV396itMIZrpgMA3+QgYYtBCy9Q/XJsSitcFIwumM0zgwygDSkxThBWRVntlg2rMofz6Konp3SFL94UIIxjdTeJbcqrbfHmi48E7FKHihXqrcyrWaPsShORxEAWH63anFu0egAKSUTsQ1M+w6qZo/VcyxZT+8p+Sxqsl0EKAPKp73yDpqJWTZ6FEI+NYAlA7idDjxOR+qMExTj9A9bVuU5lR4bTQXuPZixmQ7qwUqLcS69QwXTLQpqz8jy6jtgivZV7HambkGBcreBZVaFdg2KYxmJKfQ68KV1P9+u3E+7rHFIaFb7lIZ9UPnrM8wAKnY7U/cIACy6ATzlsDvjGW9x4Y+Np0JUEZirlEwLoBfM/mcp5b8kelmymxxjzK9uiJKz0dbEzZkJCohqJOkwkrlXqireXdbUNp5jUUwG/kxgJJ40LAqA1guhajbstIreaHqsy6FmVw+eyiiQraE4XQ27Zp5yP+22iHHGCopIBPoOmyYY09KwyxqUO2bXLywpNpwUFJMxN3MUH6/bMZn4kRLcxSqIv/dhS7wCMxUB6zKewHivp9cExgIhitwOHF27Ve53Gj37p6PE4zTe6ykWLg+sfJMK8FqQDTTlinHGBLLN0rDToFcIWP02lf00eDLjfUzHOYzk7CFAZhTI1lDsdqanYYOit3Onau1uMrRrUOJxqkaAoQnTFIG0LEaAVW9R2U8WdAw+RzD6htQQIRMygIrdzkklI2WsuFvVYlnQaeEcek2ogUqGgqCIwZg/pKoxu3apRmJOV8bnVBZFuozk7aqp2asPZryP6fDFFqBNdo3NXAMryYTeVW8FJOz4Wcb7mI7JB8vjVNommGhRpMtI7gGnB7b/JON9TMc59JoY6CxK16IAlRbsrVBWhcmYjFF4nMp6AvPoTfd+nnuVGh2w46cZ72M6znG1DZ5QbttcWBRCiI9G3y+1bHWbYcwfosTtVIIiw/iEhhJPmj5sUK6vxhWW3mglHmc08FeRceAPMmSc1bNVK+4dPzHdPXGOht13UFVH16bX/DAWaccoQGUDLX6DYpwmF2edw0hMSo0FxTjTtijcxSpWsfch8I9kvJdY+GItRhNqgjQUuR34QhFjowKmw+lSyt7Bx2HoTMZ7icVEbLq3idc3EfQsivdG379u2eo2w1ggzAJ3rzJdTYhPgKqlSMsVo2HN21Whkjbb2SSMxT5YvQdMK/0vdmdI7+p3KA3J5GKlSUHhdilGUj0HXN6Mz6tcMRkItTXvgIl+0yvxJ6bTW1ILpbUZnzft4K6Gte9RNUImx6I0xlniid7PDrdSPDJEsdtJOCIJhtOsKl/7HlUjtP3HGe8lFto9d64ikFnqsx70BMU+IcRxYLEQYlfMa7cQYpdlO8ohxvwh1ssoQ55tjiFVkomGDbDyzSrff/O3TdmPholAiGK3E4dDqBvNBLcTaK6nDCZ8Lb1NBfE3f8eU/WjQ9jTpejJB24QMgrsa5r9edXN9+bum7EfD+KQrxhGl1xxtU6M3LQ0bVC+kljWKXhNbeoxPV3xqF5iSAaSln/rSjVPUzFUtTLb9r6mT784JZvfsV8XBxdWmnX869LKe3opqL34YuDXmdUv0Pe8w5g+xJrxT/dNNerCKXBmY6qDcE6vfqtIKx/pM2RMoi6LE41R1C6PdptFbkm7WkwZ3Max7rwri9x8zZU+gGInLIfAQUsHsDFqzxKI4E9ciqNbj6+9VFlTHDlP2BFOCsUQrHjWR3oicKmBMGUIoevsOmDpP+5w6it590JDezI3pmBQUmTzD696rgvgm9vfSrLoityNa82UOvYmgG8yWUnZJKVcBnUzVT3RIKU1pUiOEuFEIcUAIcVgI8ak433uFEL+Ifr9ZCDHHjHUTIeQfZ9XEy2rymkkdGL1ux2Q7hbRx0fsh7IctPzBlT6ButBKvE7qjPa3SnAI2HRkF/zS87s8VA335flP2BEpQFHucqqV5JJRxMaUG5bOPGJsTnggXvgs8ZfDSN03ZE8BEIOqa8PeBbzDtueDTMck4Axnc08vuVK4wE63GSQ1b+FUPLTvRu+gm1Xnghf80zYqaTFZwOaJdkHMoKACEEFcCh4BvAN8EDgohrsh0YSGEM3rOm4ClwFuFENOf3vcBA1LKBcB/AP+W6bp6WDvxIkVyQqW1mQSP0zFZs5A26herdtyvfBeCPlP2pQL3rpjmh+YxzozprWhRzGTbj1W8yARMaBaUCc0eY6FV32dEc3EVrHmnym4b7jBlX+PBEB6XA2dfZnPfp6M4U1cMqOrwdX+m4jImVWpryknpULT5oUkatin0Ol1wyUfUnBmTKrU1wVjq64zOfTfHlZoIRtJjvwZcL6W8Ukp5BXADimlnivXAYSnlUSllAPg5cPu0Y25HzcQA+BVwjRAWNFsHkJI3hx6iz9uuei2ZBK/bBEEBcOlHlYto6w8zPxdKI1EWxR7V36mixZTzelwOAuE0s0RicclfqiLATeZo2ePBsGrS2LNPtTow0WcPKXYIjocN71fa5vP/acKuYjqL9pgsKLTGj5lajRs+qFqQ/+nfTdhVTBbQgNbs0SzXk6I3I9cTqKSF0np4/msm7CpGUAxHB36ZZEElghFB4ZZSHtB+kVIeBMyoE29FTczTcDr6WdxjpJQh1HyMuKkbQoh7hRBbhBBbent7U9+Nf4SAt5ZD89+j3B4mwetyZu56ApU2Oudy2PhVVVmcIcZjNeyGZea52qK9hdL2YWtoXglLboNN3zCl4FAL3tOzTwU60xz3Oh0avRkrA9VzFDPZ+gNTCg6nru8+xaDSHPc6HZpgzCiADyr2tv5e2PMbFYzNEBPRaX6OvmjGU828jM8JMfRmKhjdxUo4HnkGTm/JeF8zBWPuLYotQojvCSGuir6+C2w1Ye14nGm6GmrkGPWhlPdLKddJKdfV19envpuiCtb83VNcfI+5/Q49Lkf6lZ3T8fp/gLEe5YLKEGP+EMVaO/UMu+TGwjTGCXD1p1Uq5fOZG7BTjHOPado1TDXdM0UZuPJv1PufMvewTgSjMZmefab6r72TFpQJ9/TFH1btuJ/9Qsan8mnXt3e/qlA2qeeRd9L1ZML1Xf8XavrdE3+fcaxiNNpuyNt/QJ0zg3EIRmBEUHwA2AN8BPgosBd4vwlrn0YNRNLQBkx30E4eI4RwAZWAddNt1EKmns7rchAMS8KZBDs1tG9QsYqNX4XRNKymGEwEw7Q5+1TjQZPiE2Ay42xYoqqXX75fTeDLAOOBMNUuvzqPFfRmakEBVLbBuvepynQtySBNTATClLqIKgLmJCpArCJggqAorVUu1b0PZZwBNRYIq8Bux3ZoXJ753qLQXE8ZWxQA3nKl/Jx6SbUgzwCjPtVFQnRsgxbrZ8kZGVzkl1J+TUp5p5TyjdHJd2Z0uXoFWCiEmCuE8ABvAab/9x4G3h39+W7gGZmx8zu70LqVmsI4AW74gnI9PfmPGZ1mzB9mUTDqUWw270Yz1aIAuPafVNfRxz6ZkRY2EQizNBI101vXmrM3pmZwm3Z9r/ikihk98lcZVaePB8IscZxUgc6215mzN6aub9oFaNNx6UdURtDv/zajOoMRX5D53miq96z15uyNKdeTaV6BNe9U8YQnP5PRBMtRf5AGT0BlPLWuM2dvOshZr6dozOHDwBPAPuABKeUeIcRnhRC3RQ/7HlArhDgMfAyYkUJrd3jN1LBB+SIv/Qjs/BkcS7+52kQgxMKJXeAuheZV5uwNky0KUFr21Z9Wo2Ez0MLGgyGWBveo1h0mMk6P2YKxtBZu+Dyc2gzbfpT8+AQY9YdYGYn6/mddZM7eAI/TZMXHXQw3fFHFyjJIDx71h1grooqAidfXlDqKWDhdcNO/Kcv2mX9N+zSj/hDr3YcBCW15LCgApJSPSSkXSSnnSyk/H/3sM1LKh6M/+6SU90gpF0gp10spj+Zyv+nAY6apruHyT6iW3A99CCYGU/5zKSXjwTDzRrdC+0WmND/UYDojAVVH0rgCHv0EjHSndYoxf5gl41tVD6+iCtO2ZrpgBNUccc7l8OQ/pV10OOILsjqwVd0nVbOS/4FBWELvBTfDBbfAM5+DzvSaPoz4QqwPbVNV/U0rTNuaW7MYzbKgQI1KXfc+lahxYlNapxjxhbhUblddrtsvNm9vCaArKIQQTiHEly3fRR7DdFcMgKcE7vofVe358IdTdslMBMPM5zS1E8dVUzoTYaoPW4PTBXferxrJPfi+lMdLSinxTvQwa3yPYkomwrQsr1gIAbf/t0rl+OW706qdCU2MsHh8m+n0up0qhhcIm3h9hYBb71NFeA++F+SXuQAAIABJREFUDwJjKZ9iYsLHKl+0WNbErEVLBCPAdZ9Vvah+c29aHRfGfEHWB15WI2Y9JebuLQ6SVWaHgbWW1S68BuCd9HGafKPNWg/X/rOaQ5xibvbwRIi7nc8RwaFGVJoIyx6sxqVw81fVLIOn/imlP/UFI9whnkMg1SQ9E2GJBQUqXfaOb6t5FY99PGVl4KrAc7hlwHx6tRhFyORQYWktvPE7qiX6g3+ecrxije8lyiLDqlDTRGgxqKCZigCAtwzu/j6M9sAv3pHycKO54ztpDHeqkatZgBHX03bgISHEO4UQd2ovqzeWL9BuNFM1bA0XfxiW3wVPfxZe+R/DfzY62MPbnM/Q1Xo9lDeZuiXTYzKxWPN21d7jxa/Dc8YN3eGRId7teoKu2vWmDCuKhWWCEeCCN8AVfwPb/y+llEqfz8d7HY/QV7rI1MAuxLhSzWacoFwyb/iyqth+9K+NC0cpeXvo1wx6mpRFYSIsvb6ta+GOb8LJTfDbD6QkHO8Zf4AxZ+XUSF2LYcQ5XQOcBV4f85kEfm3JjvIMXrcFricNQigtLDAGj35caSUbPpg0xbdq42cpZYKDq/4Sc+qxp2B6cHc6bvoy+EdVIDAchCs/BQ59fcez8UtUiwFeWPoRzBWLsemxFigCoAL5/mF46RuqR9UNX0gaUwo9/3XmOzr546JPcbXZ6d5WWVAa1v8FjHTBxq+o63vrfUmLI0NbfsgqcYQnZ3+G60x0OwE4HQKHsMCi0LD8LjWW96l/Us/vnferinU97PktG+QO/tD0Ya7PgtsJDAgKKeV7kx1TQGJM+uwzmVmgB6cb7vkR/PrP4YlPq34yN30ZyuIUHUoJz3yOusO/5L/Dt3NZs3n55hosFxQOB9z+DeWH/tO/qbz5W/4TKqcX9Ufxwn1U7/gOPwldQ2ubedk/Giy1oEAJ/Ru/pFKEX/qGarly232Jhy5t+zGlz3+eR8IXEZp1renbcbuiMQqr6AVVWOoqgj/+q6oDeeO3ExdJ7nsE5+N/w8bwck61m+tm0+BxmdDYUw+X/ZWi94m/g+9eA3d8I3EK9/EXkL/9ILsi89gz621cb92uzoGRpoCLhBBPCyFejf6+UgjxD9ZvLT8wWUdhlUYCqsnaPf8Lr/9H2PcI3LdauSpOblZtxMf64MDj8KNbYeNXOTb7br4WuofKYnOqV2Oh0WuJq02D06WExU1fVinCX1+rLKpjG1W7j/F+OPw0/OQeePIf6Zl1E/8cejcVFtBrqWtCgxBw4xfgjm+p6YvfuAge+jAc+aO6tuP9qmDtgXfBwx9mpPUKPhF8P+VW0GuVzz4WQsCVn4Q3/0QNsfrWJfDgX6jZ06M9KtPv5Evwmw/AL95OsHYJHw5+hPJic9qyTIfb6bD2+QXV6+sdD6pOv9+9RsUt9j+qrCvfEJzeCo/8NfzoFmR5E38R+DjFxUXW7ikGRlxP3wU+CXwHQEq5SwjxUyD9JODXEKYsCgsZJyhN+4pPqOD0s19Sg442/fe5x5TUwc1f5bnANUQO7KWiyLy0WA2Wa9gahICL7lU+6ee+ogbDTI/TeCvhus/xUukbCR7aZQm9lltQsVj9Nph3larM3/6TmVPT3CVw1d+xq+kd+I7spLzIfEHhcjpwiCxcX4Alt6jUz41fVbTufmDaZorh4g9zZMlHGfrmK5RbcH1B3dNZoXf+6+FDL8ML/6Xu5X2/O/d74YR176N77Sfo+c+tVFmgCCSCkf9siZTy5WmJT+aNaspzWJIeq4f6xXDPD5Qlcfx51Zvf4VKft18M7iKGnz4EYAkjMbWlhRHUzFWm+o1fhBMvqLoDbR727EvBU8LwS2p8SoUV9DqzTG9Fi8r+uu6zcOJFNV9DSpUlNfcK8JYxvLsTwDLGqXUIzgpKa5U1dc0/Kiui9wBEglA9F2ZfAiU1DBxR6aWVFloUllpQsSiqULRe+beq6LJ3v6rgrmyFeVdDSQ1nz6jW+1Ul1tAbD0bupD4hxHyizfiEEHejBhkVYACmt/AwiuLqhKmvw74gxW7nJFM3E5bHZBKhqAIW3xT3qxGf0mssEYxmt/AwvHCpsqbiZPkMjgcBqCqxRuP0OLOkYcfCXQzzr1avaRgYU/TWlFrDOC2PUcSDywNzL1evadCub7VF1zfudgwc8yHgfuACIcQZ4Bjwdkt3lUfQmHFGg09MxtBEkIpi67RNyKKGbQAD4wE8Lsdkgzcz4XAI3E6RfUaig/4xlZNvKeO00fXtHw8AUF1qDeNUFoV9WswNROm16vrGg5Gsp6PAtUKIUsAhpRyxflv5g6kCJRs9WGNBqi0yWyfrRrJtUejg7GiAulIPVtWNmjLF0EScHQtQ5nVNWrNmIycWhQ4GxqKCwsJ72k7XVxMU2XQ9Gcl6qhVC3AdsBJ4VQvyXECLu8KACZsIVbXkQMqPNuEk4O+anrsxrybldTgdOh7CuriANnB3zU2sRvZAj14QO+scClmqbdqS3vMg12ZfJbLhdWYxRGIDmarPKtRgPRv6zPwd6gbtQrb57gV9Yual8QtaDnQZwdjRAXZmFjMTpsJ1F8VpjnNWvIXoHxq29vl67WVDj1grGeDCyUo2U8nNSymPR178CVVZvLF/gckQtChv5OPtGrdew7WRB9Y8FqLVSMNrNZz8WoNZCxpnVLCAD6B8LWOZ2AlVkaKfrOzBuLb3xYERQ/FEI8RYhhCP6ehPwqNUbyxc4HQIhIGSTG208EGI8ELaUcbqd9nmwpJT0jVrnagP7+eyz4nqyyfWFLNBrM8HYN+qnvty6+zkeEgoKIcSIEGIY+H/AT4FA9PVz4K+zs73zH0II3A6Huf3sM8DZURUIs5JxuhwO2wjGsUAYfyhiMeN02ibYKaXkrMUWhd2Cu93DfhosZJxumykC3cN+GitsIiiklOVSyorou0NK6Yq+HFLKjCa/CCFqhBBPCiEORd+rExwXFkLsiL4yGzKbQ7icwjaMs29UpU5a6ppwCdukE/aOKHottShcDkIZjC01EwPjQQKhCE2V1rV3sFOMIhiOcHbMT2OFxfTa5PkF6B720VCevfYdYKyOAiHESmBO7PFSyky6x34KeFpK+SUhxKeiv/9tnOMmpJTWTw63GHby6XYPqyE4Vj5Ybod96O0cUnOJWyxknG6HsB29TVYyThvdz70jfqTEWsFoI4tiPBBixBey9PmNh6SCQgjxfWAlsAfQ/luZthm/Hbgq+vOPgGeJLyjyAm6nIGiT4G7HoBIULVXFlq3hdjpsE7zvjNJrJSOxU0GWpgi8ViyKLo1eiy0KuwjGnmFlIWfb9WTEotggpVxq8rqNUspOACllpxCiIcFxRUKILajeUl+SUv7W5H1kBYpx2uNG6xicoMjtsLT83+W0j4atMZLmSusEo8sp8PntUTfSOWQ9vbaykIeyYCHbyKLQFAE7up42CSGWSin3pnJiIcRTEHdOzN+ncJp2KWWHEGIe8IwQYreU8kiC9e4F7gVob29PZauWQzFOe2icnUM+WiqLLatSBlV0Zx8LaoKqEjfFHmuqlMFeFlTXkA+HwNI6GTtZUJpgtNqCsgu9pweirtQq+wmKH6GERRfgR418l1LKlXp/JKVMODVFCNEthGiOWhPNQE+Cc3RE348KIZ4F1gBxBYWU8n5UTyrWrVtnj6sahZ189mcGJyx1OwF4nMI2LUs6h3yWatcQdS3a5PqeHpigqaIIl4XFWHai92T/OGVel6UWsp0silMD4wgBrdXW3tPTYeRu+j7wTuBG4Fbgluh7JngYeHf053cDD00/QAhRLYTwRn+uAy4FUrJq7AI7meqnB8ZptVhQuBz2yQI6cXaM2TXWjot02ej6njg7xuzaJKM0M4TbaZ+CypP948yqKbHUQvZE64Kk0RneFuJk/zhNFUWW9fFKBCOC4qSU8uFoVfYJ7ZXhul8CrhNCHAKui/6OEGKdEEKbPrME2CKE2An8ERWjOC8FhUqPzf1NNuIL0jcaYE6dxYzEJqZ6OCI51T/B7DprBYXHZoxzdq3VgtE+FuPJ/nHaa6y2GBWbDNvgGp/un2CWxYpPPBhxPe2PTrT7Hcr1BGSWHiulPAtcE+fzLcCfR39+EViR7hp2gtsmPvvjfeMAzLWYcdolXbRjcIJAOMJcizVsl8MejHPUH6JvNEC7xYLC43QQtIHFGIlITvWP8/oLEuXCmAPNjRcMS7KsyM/Aif4xLltQn/V1jQiKYpSAiJ3jnWl67GsKbptoYEf7RgGYW1dm6Tp2saBOnFWC0XJXjMsuisAYAHOsFow2Sc44MziBPxSxnF53tAN0IByhmNxJimFfkO5hPwsarH1+48HIPIr3ZmMj+Qy7xCiO9IziEFjumrALvYd61OiU+Q0WMxKbWFAavYsarWUkbqeDcEQSiUgcDutiA8mQTXoh9/3aDvcoRW+hHQWFEOIHRMegxkJK+WeW7CgP4XI6GAvkPs9+T8cw8+vLKHJbqxW5beKaONg9Qk2ph3oL23eAur52sKAOdI3idoqsBLMBgpEIXkfuNOxD3RrjLLd0HXeM6ymX0ASFLS0K4JGYn4uANwId1mwnP+F22KPX056OYTbMq7F8HbdNXE/7u0ZY1FhmaUYMRNMnbXB9D3aPML++zPI5BbGt873WTNQ1hAPdIzSUe6m0eICP5nrKtdW4v3OEIrfDnsFsKeWDsb8LIX4GPGXZjvIQdnDFnB310zXsY1lLpeVr2SFdNByRHOga4U3rZlm+ltsGTR+llOw+M8RlC+osX2tKw84tzbtOD7Gi1fr72S707jw9yIrWSpw5cPelo3osBOxV+mxz2CG4u6djGIBlLRk1/jUE5bPPLb0HukYYD4RZ0279jC2300FE5jZ9smPIR++IP0v0ahp27ugd9gU50jvK6lnZub6Q23HGwXCEV88MsaotNzPjjMQoRjg3RtFFHjfwswJ2SCecEhTZ0cByrX1tOzkAwJpZcTvYmwpXjGvCmSOf/fYs0msHDXvXqSGkhFVZEBTa9c1ldfbB7hH8oQgrs0BvPBhxPVkbKXoNQBUo5VbD3n1mkLbqYsv9uWCP4O72k4PUlXmYZXExFkzNRc+lxrn95CBel4MLmq1/XF2TWUC5o3fn6UEgO4LCDtd356khAFbbzaIQQui6l6SUJ83fTn7C5cxtSwspJZuP9nPFouwU6sS2PLA6kJwI208OsKa9Oivra8HdYCgC2e3+PIkdp5T/2upANpxbV5ArbD85yLz6UiqLs6H45D6YvfPUINUl7qwoPvGgZ1E8inI5xT5pEqgHGiCHlSfnGXI9+ORg9yhnxwJcPL82K+u5YloeaA9ZNjEwFuBo3xh3r2vLynqumHTRXMAXDLP7zBDv2jA7K+tN+exzQ28kItl2coCrF1tbka1h0tWWw2f4pWNnWTs7O4pPPOiNQl0hpVwZfV+BagT4AjAK/FW2NpgPcDlETs3WTUf6ALgkS4Ii18G/5w8rei+amx16PTnOs3/leD+BUIRLs5DxBLGMMzf0vtoxRP9YgMv/f3tnHh1HdTX43+1Fi2V53y1v8ooxBowwGLDZjG2WhA/CGsIQEmBIQiBDICGTmTk533fmzAQI+UJgyAeEkB34IIQkJnjBxgtLvIE3bC2WZVu2NsuLJGtXv/mjq+RG9Cap61V15/3O0VF3qdV136uqd9+97757p+tqr2VRuHQ/H6xv5kB9Mwun60/dYZPQThWR6SLyMvB3YCswWyn1c6cFyySCAXd99h/sq2fCsFwKhuqJv3bbNbGupI7BuUEtETFw2jXhVojsuuI6sgI+LtCwRwYiXDEuWRQbSsMTgUu0KQp3LYoNZXWAvvZGI6aiEJE51p6JNwjvm5ijlHpRKdWhTboMIehzL01xe2eID8vruahQ300WuSFLN6GQYl1JHQunj9AWb+72zt11JXXMnzyMAVl6dr9lubyYva6kjjPHDWKEwzvubQI+d11tG0qOMn5ILoUOZ32OR7w7aztwiPBaxXxgfqR/TCn1oLOiZQ5u+uw/Kq+nsbWTq2aP1nbOYMC98Mk91Q3UNbZxmSb/Nbi7c/fIiRZKa5u0bCy06V68d6G9Da0dbDtwnHsXFWo7Z1bAtpD1K8a2zi7eLzvKtXPHurY+AfEVhcnllCIiZ5y60xSv/LSaAVl+rWZr0Oeeoli5uwYRWDRDpwXl3gx7xe5qAK44Q6NidHEi8O6eGjpDisVnaJz4uJgU8IOyehrbOll6ZrSq0vqIqSiUUr/WKUgmE4zw6epMUxwKKVburuHSGSMdTwQYSTDgjutJKcXfdhzhginDtBaftwdON9Zklu+oYtaYfKaO1Jco7vREQL9iXL6jmrGDczhX48azgIsbDJfvrCI/J8BF0/QEZsTC+aBrg2uLYZsrjlHb2KZ9NhJwyaIormlkX90prp07Tut5gz53FrOrTraw5cBxrps7Vut5T08E9La3sbWD9aV1XD1nrNb05m6lLOnoCrHq0xquOmO09tKnPXFFUYjIzSKyW0RCIlIU53PLRKRYRMpE5DGdMqaS7qgYzeF1r2+tZGB2gCVn6jPTwb0H66/bj+ATuHqOXsV42hWjt73Ld1QBcM1ZehWFPRHQbUG9s6ua9s4Q1+pWjC5NfNbureVkS4f29kbDLYtiF3AjsD7WB0TEDzwLXA3MBm4Xkdl6xEstbuTGOdXWyfKdVVx71lht0TA2brS3syvE61srWTRjpLZoGJvuxV2NUTFKKV7dfIhzJw6hUKPbCdyLenptyyEKR+YxT0Piw0jsiYAb7R2Zn82lmjIqxCOZfRQzRORdEdllvZ8rIv+jPydVSu1RShUn+Nh8oEwpVa6UagdeAa7vz3ndwo0Z9ju7qmlu79K2OzmSgAs7d9eV1FHT0MZt5+tPbBx0YeDcdvAEpbVN3Ha+vmgnGzdSWpTVNrG54ji3Fk3QHv1jTwR0WlC1Da2sLa7jS/MKup8nN0lGgheAHwAdAEqpHcBtTgplMZ5weK5NpXUs7TgdFaPvRnt9ayWThg+gaJLz2UR74oZi/OOmQ4wYmM2VGqN/bNywoF7dfJC8LD/XaV6PgcgKd/qu739uOUTAJ9w4T//Ex43r+8a2w3SFFLe4MNGLRjKKYoBSalOPY52J/klEVovIrig/yVoF0aYNMe9MEblPRLaIyJa6urokT6EH+0bTNSMpq23kw/J6bj6vwJXYa90PVk1DK2uLa7npvAItSfF6onuG3dTWyd92VPGFs8eR50KJue6JgKbgjI6uEG9sq+SKWaMYma8/66LfJ/hEn8UYCile23KI+ZOHaXcrxiKZu+yoiEzFGqRF5CagKtE/KaUW91O2SiDSri4gTglWpdTzwPMARUVF7tfhjCDo1xsu+qv3K8gK+Lh9vjv1pXTvzP7jpoN0hRS3uuCGAf25nl7fcojm9i7X2qvbtbh8RxVHm9q5/QL36qXprLGyrrSO/UdP8dCV07WcLxmSURTfIjwAzxKRw8B+4CuOShVmMzBdRKYAhwm7u76s4bwpR2e2zRPN7byxrZIbzhnPcM2LujY6LajWji5+++EBrpg1iikupTjQmeupK6R46f0K5k0cwrkT9bsVQa9rUSnFixvLmToyj0tdTIqX5fdpmwj8csN+xgzK0R7NFo+Edrq1mLyYcHrxWUqpS5RSFf05qYjcICKVwAJguYissI6PE5G3rfN2Ag8AK4A9wGtKqd39Oa9bnK6Q5fyN9sdNh2jtCHH3JZMdP1csdC7u/vnjw9SfaueehVMcP1csdLraVn1azcFjzdy7UF8Ki57oDBf9x/5j7DrcwD0LC7XunehJwC9a2runqoGNZUe566LJZAXcX8S2SaYUajbwJWAyELB93kqpf+3rSZVSbwJvRjl+BLgm4v3bwNt9PY9X0GVRdHSF+M2HFVw8bTizxjhfGzsW3a42h9sbCile3LifM8cNYkGheztXde5UfmHDfiYMy2WJiykdfD7B79MzcL64YT/D8rK44Vx341iCmoqPvbRxP7lBP192yW0ci2RU1luEw1I7gVMRP4Yk0TXjfHtnFVUnW/naxe7NriHC9eTwYue6kjrKapu4Z+EUVxOmde9Udngg2XbwOFsPHOdrF0/Rlhk3FgGfOG4xltc18e7eGr5y4SStKWiiEfT7HPcI1Da28tYnR7i5qEBLyeLekMwaRYFSapnjkmQwp7NtOnejhUKKZ9eWMX3UQG2Vv2Khayf68+vLGTMox5UQ0UgCmiyK59eVk58T4GaNmWJjkeX3Ob4G9cKGcoI+H3dqqtwXj6BfHJ8I/Or9CjpCIe52eaIXjWQsig9E5CzHJclgsjTs7Fz5aTUlNU08cMU0V325oMeC2lJxjA/L67ln4RRXQmIj0ZFmvLi6kXd2V3PXgskMdCEkticBv7MWxeETLby+tZJbz5/gSkhsTwIORz0dP9XObz6o4Lq541wLyohHzDvO2okdsj5zt4iUA22E9zcopdRcPSKmP07n71dK8fM1ZUwZkef67Br0+Ox/9m4pIwZmcccF7s82RYSAwz77n68pJS/Lz9cv8cZs0+lw0efeKwPg/sumOnaO3uC06+ml9/dzqr2Lb18xzbFz9Id4U5PxwDm6BMlknJ5hr9lby+4jDTxx01zXfdfgfHbRrQeOs6H0KP/9mlnkZrnru7ZxcoZdWtPI8p1VfOPSqQzNy3LkHL0l6GC4aNXJFl7bXMlN501g/JBcR87RW5x0PZ1s7uDl9yu45qwxzBid78g5+ks8RbFfKXVAmyQZjJOlMpVSPL2mjIKhufyLy5EhNk6nGf/Zu6UMz8viKx7wXdsEHfTZP7O2jNygn3tcDIntSdDBcNH/WFdOSCm+6RFrApy1oH71wX4a2zp54HLvbLDrSTxFMUpEHo71R6XUUw7Ik5EEHAwX3VB6lO2HTvC/b5jjuq/exskNWR8fPM76kjoeu3qW9qy48Qj6fY5YFPvqmvjr9iPcu6iQYR6xJiDss3fifq5paOUPmw7ypXkFTBg2IOXf31fCijH117ehtYOXNu5nyezRzB7nXkh7IuI9aX5gINFzLhl6gVPhokopnlxZzPghudx0njeSh4GzPvun3y1l6ICgJyJhInHKNfHMmjKyA35XN9hFwynX03Pv7aMrpPjW5d7y1Qf9PpraEqa46zW/2lhBQ2snD3ooXUc04imKqv5sqjOcJuhQuOiK3TXsqDzJ4zfNdb0CVk8Cfkl5ezdXHGNtcR3fXzbLlWR48Qj4Ur/Yube6gT9/cpj7FhZqr7GRCCdcT5XHm/nDPw5y07wCJg73jjUBzliMx06188KGcpaeOZo54wen9LtTTTxfhbEkUoQTxdm7QoqfrCxm6sg8bvTI2kQk4SiR1LVXKcXj7+xlVH42X71ocsq+N1U4YVE8uaKEgdkBvuEhX72NEwPnv68uBYHvXOW92bUTFvL/W1tGc3snjyyZmdLvdYJ4iuJKbVJkOKcLn6Tuwfrzx4cprW3iu0tmeqKwSU9SnfLgveI6Nlcc58Erp3sm0imSVC92bj1wjNV7arj/0qkMGeCdtQmbgE9SunhfUtPIn7ZVcteCSYwd7I1Ip0iCgdRe38MnWvjNRwf40rwCpns00imSmCOMUuqYTkEyGdtnnyqLor0zxE9XlzBn/CCWuZjzJx7BFIaLhkKKx1cUM2n4ANdSaycikEKfvVKKH79TzIiB2dx98eSUfGeqyQr4UmohP7mimLysAN+8zFtrEzZBX2oXs3+2ugQUfOeqGSn7Tifx3lQ0Q0nljPOVzQepPN7Co0tnub4LOxYBX+rCRf+64wh7qhp4+KoZnons6klWCn3260rq2LT/GA9eOc1TkV2RBFI4cH588DgrP63hvkWFntkn0pNUPr9ltY28vrWSOxdM8sw+kUR486nLQAIpCq9rbu/k6XfLmD9lGIumj0iBZM6QKouioyvEU6tKOGPsIL7ggV3nsQikyGcfCikef6eYCcNyXan/nSypGjjD1tNeRgzM4mse2XUejVRajD9ZWUJu0O+pfSKJMIpCE1kp8tn/+oMDHG1q43tLZ7qaMTURqRpIXt18iAP1zTy6dIZnrSdIXRTQ8p1VfGpZT16qR9CTVF3fDaVH+aj8GA9cPs1zkWyRpMpi3H7oBH/fVc29iwpdKyzWF7x7J2YYAb/Q0c/wyZMtHfxi3T6umDWKosnDUiSZM6RiBtbS3sXT75Zy/uShrmfETUQqBs6OrhA/WVnMrDH5fPFs70WyRZKK8OdQSPHEimIKhua6WuY0GcIWY/8VxeMr9jIsL8tTu+yTwSgKTaRiIHlhfTknWzrSIpwuKwXhoi9/UEFtYxvfWzbL09YTWPUZ+jlwvrblEBX1zTyyZKYncnbFI+j30dHP8Oe/76pm5+GT/LfFMzy3D6gnqdhguLH0KO+X1fPA5dM8kQG4N7iiKETkZhHZLSIhESmK87kKEdkpIp+IyBadMqaa7ICPtn4oirrGNn65cT9fOHucp7f62/Q3LfPJ5g6ee6+MK2aN4nyPW0/Q/30jtvV03qShXHmGt60nsFxt/VCMnZb1NGP0QM/kKItHll/oCIVQqm9tVkrx+Iq9jB+Syx0Xett6ioZbFsUu4EZgfRKfvVwpdY5SKqZCSQeyA37aOvo+kDyzppT2rhAPp0k4XX+jYv5j/T4aWtNjMxLY+0b63t5ff1hBTUMb308D6wn6byG/vrWS8qOneHTpLM9bTxCe+CgV3ujaF97ZVc2OypN8Z/F0z1tP0XDF/lFK7QHS4oFIFdlBH22dXX3630PHmvnDpoPcev4ETxY1iUZWwEdja99y49Q2tPLS+/v5YppYT2BHefVt4DzZ0sFz7+3jspkjmT/F+9YThMOf+xrl1drRxb+vLmXexCEsTgPrCSLr3it6O853doV4YmUx00cN5MZ53snJ1hu8vkahgJUislVE7nNbmP6QHfDR1kfXxE9XleAT4SGPJw6LJDvg73N7n15TSmeXShvrCfq3eP/8+n31HYQXAAARJUlEQVScbOng0aXpYT1BuOZIX/fJ/PbDA1Q3tKbF2pONna+tL21+Y1sl5XWneGSp99eeYuGYRSEiq4Fo24Z/qJR6K8mvuVgpdURERgGrRGSvUiqqu8pSJPcBTJzoPR9gdsDPieb2Xv/f3uoG3vzkMPctKmT0oBwHJHOGnKCPto7eW1AH6k/xyqZD3Hr+BCanifUEfXfF1Da28tLGCr5w9jjOHOftxHCRBH19iwJqaO3g2ffKuHTGSC4sHO6AZM7QXVOml5Mf23o6Z8IQlswe7YRoWnBMUSilFqfgO45Yv2tF5E1gPjHWNZRSzwPPAxQVFTlb5b4P9NWi6E4Md2n6bM6BvlsUP11VQsAvnk+73JMsv/Spvc+sKaOjK8R308h6gvDAGbJ89r2ZJb+wvpwTzellPUFkTZneDS2/++gAVSdbeeqWc9LGeoqGZ11PIpInIvn2a2AJ4UXwtCQ72PuB0+uJ4eKRE/TR2kuLYk9VA29tP8JXL5qSVtYTQE6Wv9ftPXSsmT9uOsgtaWY9wemBszdWlB25d+3csZ5Pq92TvtSUaWzt4Nm1ZSycPoIFU9PHeoqGW+GxN4hIJbAAWC4iK6zj40Tkbetjo4GNIrId2AQsV0q944a8qSAn0DtXTDokhotHTrD3A+eTK4rT0noCyLUmAqFezDjttacHr0gv6wnCmQagd4ri2bVltHWmn/UEfasp88KG/Rxv7uB7S2c5JZY23Ip6ehN4M8rxI8A11uty4GzNojlGOOop+YfKTgz3b9ef6dnEcPHICfpo7UV7t1Qc4929tTy6dCaDBwQdlMwZcoLhUJi2zlBSadC7154WFjJmcHpZTxC+vgCtHSHykxDfLkp0S1EBhSMHOixd6gn2UjHWN7Xxyw3lXHPWGM4qSC/rKRqedT1lGr3x2UcmhrvVw4nh4pET8NMVUkk9WOGiROlrPUHYogBoSdKK8nJRomTItSYvLe3JtdcuSpRua082vVUUz67dR0tHFw9flV5rMbEwikIT4cXs5B4qOzHcd6+a6enEcPHItmacySjH90rq2FTh7bTaibBn2Mkoiq0HjrN6Tw3/dVFh2q092fRGMZZ6vChRMtjPYWsSm2YPn2jhdx8d4KbzCpg2Kv2sp2ik5yiUhmQH/HR0qYQ7Oz+bGM67abUTYbtiEq1TpEta7UQk216lFE+sCKfVvvti76bVTsQAy73W3J54U+WTK4sZkBXgGx4tSpQMeb2woH62ugSAhxan31pMLIyi0IQ9w04UNWEnhnt06UxPp9VORE4guYHzbzuruosSpav1BBEz7AQDSbqk1U6EvQ6TqL2fHDrBit013LuwkGEeLUqUDLZiPJVAMZbVNvH61kq+cmH6FCVKhvR9MtOM7IDtion9YNmJ4YomDeWKWemR2iAW2cHEpnpHV4inVhYzc7T302onIhmLIlzSNZwYzutptRORrOvpiRV7GZ6XxdcXpq/1BHQr9UQW1FOriskN+vnW5em59hQLoyg0YScCi+ez704Md3X6pDaIxekooNgDSaT1lK6pDWzsGXY8xfj3XdXsOhy2ntIxMVwkp11Psa+vnVb7m2mYVrsneUm0d2flSd7eWc3XF6ZXUaJkMIpCE90WRYyBxE4Md/nMkWmRVjsR2QkW/2zrad7EIWmRVjsRtqst1gy7syvET1alT1rtROQkcLXZazHjh+RyR5pbTwADbIuiLbaieHzFXoYOCHJvmltP0TCKQhOno4Ci32inE8Ol/+YciLAoYgyc6ZZWOxG5WfGjnuzEcN9Ng6JEyWBbFLHau2J3NdsrT/LQ4und90I6Y7vaYq1RfLivng2lR/nmZdPIz0m/fUCJMIpCE/FcT3ZiuHRKq52Ibp99FMVoW0+XzhjJBWmUGC4e2XEW7zMlMVwkdhhzNFdMV0jx5MoSpo7M48YMsJ4A/D4hJ+iL2l67KNGYQTncuWCSC9I5j1EUmogXZ28nhkuntNqJsNsbzdVml3RNt8Rw8Ti9RvH562snhvvespkZYT3BadditPv5T9sqKatt4pElMwn4M2eIycsKcKrt8xbF6j21fHzwRMZYT9HInKvocQbnhs3Rk80dnzl+sD6c2uC2+emXGC4e3eGxPSyK2sZWfrlxP9elYWK4eMQKj41MDHfR1BFuiOYIPp+QG/TT0sMV09YZtp7mFgxm2ZxoVQbSlwHZ/s9ZFF0hxZMripkyIo+bz0vPokTJYBSFJuwY8mOnPluT4qlVxeG02mmYGC4e+Tlh10RPxfjsmjLau0J8N01KnCZLToxw0RetxHCZZD3ZDMj6/MD5h38c5PCJFh5dmjnWk01eVuBz4bF/2X6Y4ppGHr5qRkZZTz3J3JZ5DFtR1EcoCjut9t0XT2FUmqXVTsTQAVn4fcLRptPttUu63lJUkDYlXZPFb82wI8u/1je18eKGcq6eM4a5BUNclM4ZcoL+zyjGprZOnllTxoLC4VwyLXOsJ5ueirG9M8RTq0qYPXYQ15411kXJnMcoCk0MyAqQE/Rx7FRb97EnVxSTnx3g/kWZtTkHwq6J4XlZ1DWebm93Wu00TQyXiDGDc6huaO1+byeGyzTrySY/J0BDy2mL8aWN+6k/1c6jGbQWE0le9mfXKF7dfJBDx1p4dFl6Z1FIBqMoNDI8L5tjp8IPlp1W+/7LpqZlWu1kGJmfTV1TWFEUVzfy5ieHueuiyWmbGC4RowdlU3MyrCgyMTFcT8YPyaXyeAsAx0+188L6cpbMHs28iUNdlswZBuUGOW65UpvbO3l6TRnzJw/jshkjXZbMeYyi0MiwvCyOnWqzihLtZVR+NndflHmbc2xG5md3WxRPrChmYFZ6FiVKlrGDc6myFEUmJobryYRhAzh8vAWlFM+t20dTeyePZOBajM3k4QM4dKyZjq4QL39QQV1jW0ZFssXDKAqNhBVFOys/rWFzxXG+feX0pIrcpCsjB4YVxUfl9azeU8N9iwoZmsaJ4RIxZnAONQ2t7Dp8kte3VnLngsxKDNeTgqG5NLZ1sutwAy+/X8GN5xYwY3S+22I5xpQRA+kMKbYfOsFza/dx5axRFGVAFoVkcKsU6hMisldEdojImyISdaVPRJaJSLGIlInIY7rlTDWThw9ge+VJHntjB7PG5HP7+RPcFslRbJ/9Q698zPghudyzsNBtkRxl7OAcOkOK+3+3lcG5wYyLZOtJwdCwErz/d1vx+yQjI7sisQMwvvH7bTR3dPGDa85wWSJ9uGVRrALmKKXmAiXAD3p+QET8wLPA1cBs4HYRma1VyhRznVVf4nhzB/96/ZyMDqcDuG5uuL01DW38z+vOyGjrCeiO9Kk83sIjaVrStTecZUVyHT7RwgNXTEvLkq69wV5rqmts46sXTc7YtadoiFLJFwt3RACRG4CblFJ39Di+APiRUmqp9f4HAEqp/5PoO4uKitSWLVucELdfhEKKZ9eWMWvsIK7KkFQOiXhp436GD8zi+nMyI5VDIv60rZKWji6+PH/iP4Xv+r3iWvZWN3LfwsKMj/yBcBDKB/vq+cZlU7vLo2YKIrJVKVUU9W8eUBR/BV5VSv2ux/GbgGVKqXus93cCFyilHojxPfcB9wFMnDjxvAMHDjgruMFgMGQQ8RSFY0niRWQ1EG0P/w+VUm9Zn/kh0An8PtpXRDkWU6sppZ4HnoewRdFrgQ0Gg8EQFccUhVJqcby/i8hdwHXAlSq6WVMJRK72FgBHUiehwWAwGJLBrainZcD3gS8qpZpjfGwzMF1EpohIFnAb8BddMhoMBoMhjFurMc8A+cAqEflERH4BICLjRORtAKVUJ/AAsALYA7ymlNrtkrwGg8HwT4srhWyVUtNiHD8CXBPx/m3gbV1yGQwGg+HzZFZ8l8FgMBhSjlEUBoPBYIiLURQGg8FgiIvrG+6cQETqgL7uuBsBHE2hOKnCyNU7jFy9w8jVOzJRrklKqag50zNSUfQHEdkSa3eimxi5eoeRq3cYuXrHP5tcxvVkMBgMhrgYRWEwGAyGuBhF8Xmed1uAGBi5eoeRq3cYuXrHP5VcZo3CYDAYDHExFoXBYDAY4mIUhYVXyq6KyAQRWSsie0Rkt4g8ZB3/kYgctnJjfSIi1yT6LgdkqxCRndb5t1jHhonIKhEptX4P1SzTzIg++UREGkTkO271l4i8JCK1IrIr4ljUPpIwT1v33A4RmadZrqgliUVksoi0RPTdLzTLFfPaicgPrP4qFpGlmuV6NUKmChH5xDqupb/ijA3O319KqX/6H8AP7AMKgSxgOzDbJVnGAvOs1/mES8XOBn4EPOJyP1UAI3ocexx4zHr9GPBjl69jNTDJrf4CFgHzgF2J+ohwXrO/E669ciHwD81yLQEC1usfR8g1OfJzLvRX1GtnPQfbgWxgivXM+nXJ1ePvPwH+l87+ijM2OH5/GYsizHygTClVrpRqB14BrndDEKVUlVJqm/W6kXDmXC/XEb0e+LX1+tfAv7goy5XAPqWUa+UNlVLrgWM9Dsfqo+uB36gwHwFDRGSsLrmUUitVOEszwEeEa75oJUZ/xeJ64BWlVJtSaj9QRvjZ1SqXiAhwC/BHJ84dR6ZYY4Pj95dRFGHGA4ci3lfigcFZRCYD5wL/sA49YJmQL+l28VgoYKWIbJVw6VmA0UqpKgjfyMAoF+SyuY3PPrxu95dNrD7y0n33NcKzT5spIvKxiKwTkYUuyBPt2nmlvxYCNUqp0ohjWvurx9jg+P1lFEWYXpVd1YGIDATeAL6jlGoAngOmAucAVYRNX91crJSaB1wNfEtEFrkgQ1QkXNzqi8B/Woe80F+J8MR9J58vSVwFTFRKnQs8DPxBRAZpFCnWtfNEfwG389kJidb+ijI2xPxolGN96i+jKMJ4quyqiAQJ3wi/V0r9CUApVaOU6lJKhYAXcMjkjocK1wtBKVULvGnJUGObs9bvWt1yWVwNbFNK1Vgyut5fEcTqI9fvOzldkvgOZTm2LddOvfV6K+G1gBm6ZIpz7bzQXwHgRuBV+5jO/oo2NqDh/jKKIoxnyq5a/s9fAnuUUk9FHI/0Ld4A7Or5vw7LlSci+fZrwguhuwj3013Wx+4C3tIpVwSfmeW53V89iNVHfwH+ixWdciFw0nYh6EBilCQWkZEi4rdeFwLTgXKNcsW6dn8BbhORbBGZYsm1SZdcFouBvUqpSvuArv6KNTag4/5yeqU+XX4IRwiUEJ4N/NBFOS4hbB7uAD6xfq4BfgvstI7/BRirWa5CwhEn24Hddh8Bw4F3gVLr9zAX+mwAUA8MjjjmSn8RVlZVQAfhGd3XY/URYdfAs9Y9txMo0ixXGWEftn2f/cL67Jesa7wd2AZ8QbNcMa8d8EOrv4qBq3XKZR1/Gbi/x2e19FecscHx+8vszDYYDAZDXIzryWAwGAxxMYrCYDAYDHExisJgMBgMcTGKwmAwGAxxMYrCYDAYDHEJuC2AwZDOiIgdmggwBugC6qz3zUqpi1wRzGBIISY81mBIESLyI6BJKfWk27IYDKnEuJ4MBocQkSbr92VWsrjXRKRERP6viNwhIpskXN9jqvW5kSLyhohstn4udrcFBkMYoygMBj2cDTwEnAXcCcxQSs0HXgS+bX3mZ8BPlVLnE97t+6IbghoMPTFrFAaDHjYrK8+OiOwDVlrHdwKXW68XA7PDKX0AGCQi+Spce8BgcA2jKAwGPbRFvA5FvA9x+jn0AQuUUi06BTMYEmFcTwaDd1gJPGC/EZFzXJTFYOjGKAqDwTs8CBRZld0+Be53WyCDAUx4rMFgMBgSYCwKg8FgMMTFKAqDwWAwxMUoCoPBYDDExSgKg8FgMMTFKAqDwWAwxMUoCoPBYDDExSgKg8FgMMTFKAqDwWAwxOX/A+Csy/fc+3fEAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "scrolled": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAFzCAYAAAD2RpSkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dd3hUZdoG8PtJQkKHBEInJBTp0iJgo4kI6IptV6youFjXvoryiYpddy2rWLCsvZcVFWkCVlpQeg0hQKiBhFACafN8f8wQJslMMpPMzHvmzP27rlycmXln5mYIuXPae0RVQUREFAhRpgMQEZF9sFSIiChgWCpERBQwLBUiIgoYlgoREQUMS4WIiAImxnSAUGratKkmJyebjkFEFFaWLVu2T1UTfRkbUaWSnJyMtLQ00zGIiMKKiGz1dSw3fxERUcCwVIiIKGBYKkREFDAsFSIiChiWChERBQxLhYiIAoalQkREAcNSISKigGGpEBFRwLBUiIgoYFgqREQUMBE19xdRJCtxKA4fK8bOvKNYu/Mgtuw7gsz9R7D3UAFyjxTiaFEJShyKEoei2KEoLnGcWHYoHKqoHxuDpg3i0K5JXfRo1Qj9kuPRrWVDNK0fh+goMf1XJAtgqRCFieISBzL3H8G89Xsxe80epG3NDXmGQwXFOFRQjC37jmDBhmy/nvvkRT1xSb82qBXNDSR2JqpqOkPIpKamKmcpJqtwOBRbc/Lx6dLteO2nzabjGJHYIA7f3HI6WjWuYzoKVUJElqlqqi9juaZCFAQHjxXh3d8y8dzcjQj2721dWjTAiG7NMbxbc3Rv1cgym6EO5Bdi5urdeH7uRuw5WOBxTPahApz21Lwy9704tjfG9G4diogUBFxTIfLD4YJifL9yJ56dtQH7DhcG/PX/ltoGVw1MRvdWDRFlkXIIhkPHinD120vw57YDlY77+5kpmHRutxClIm/8WVNhqRC5FJc48NPGbEz+Zg12HDga0Ne+eUgHXHt6ChIbxAX0de1kz8FjGPDEj14fv2VoB/zznC4hTETHsVS8YKlEtvzCYizNzMVD36xG5v78gLzmpaltcf2ZKejYrD5E7LtmEWqqirHTFmHxlhyPj2c8MdrWa3JWw1LxgqVif8UlDnyxLAsTv1pV49e6bVhHXHN6ChLqxQYgGVWXquLqt5fgl037Kjy2/tGRqF0r2kCqyMJS8YKlYg8HjxVh6vx0vP5TRo1e5+G/dMPF/dqgQe1aAUpGwbZpzyGc/fzPFe5fN2Uk6sSyXIKFpeIFSyU8FBSXYP2uQ5g6Px2z1+6p9uv0aN0QT154Mnq2aRTAdGQFB/IL0XvKnAr3Zz51roE09sdS8YKlYh0lDsWmvYfw3OyNNSqOv/RqhUfO785NVBEqv7AY3SbPKnPff685BUO7NDOUyJ5YKl6wVELvSEEx0rbm4u7PVmDfYc/nKlTlhsHtccOgDiwO8iotMweXvLawzH1bnhzNgycChKXiBUslePILi7F8+wE8P2cjlmb6P33IbcM64tL+SWjVqDZ/EFC1JU/8vsztD8YPwBmdmhpKYx8sFS9YKjXncCi25eRjYcZ+PPH9OhwqKPb5uS0a1sZzl/bCgJQmljnrm+xn78Fj6O92vkuHxHr48e4h5gLZAEvFC5aK/xwOxYY9h/D0zPV+TSD4r7/2wqgeLVAvjjMBkRnl11q4E7/6WCpesFSq5nAoZq/djRs/+MOn8YkN4vDYBT0wrEszzj5LlvP0zPV4dcGJyTo3PT6K36fVwFLxgqVSlqoiY98RXPfOUmz18Qzz167si3O6t+B+Dwob5c9tWTLpLDRrUNtgovDDUvEikkulsNiBDbsP4ZOl2/Dh4m0+PeeR87vjqoHtOB0Ghb0Sh6LDAzNKb//0zyFo16SewUThhaXiRaSUSn5hMVZm5eHbFTt9LpD+KQl44+pUNKrDs8vJvtz3s/z8z6FIalLXYJrwwVLxwq6l4nAoMvYdxtMzN2COjycSPnvJybikXxtuxqKIw2LxH0vFCzuVSkFxCT5Py8L//W91lWOTEurihbG90adtY5YIEcoWy9JJw3lJgiqwVLwI91LJPlSAuz5b7nG2VnfnntwSdw7vhA6JnI6dyBv3YuFRYZVjqXgRjqWyekce7v9qFVbtyPM65tyTW+LeczojKaEuS4TID+7FwvNYvGOpeBEupZK+9xDeW7gV7y3c6nUM94kQ1Vz5o8JYLJ75UypGT3cWkbcBnAdgr6r28PD4FQDuc908DOAmVV3heiwTwCEAJQCKff0LW1lG9mFMeH8Z0vce9vj4ncNPwvVnpvAsdaIAiY4S/HLvUJz5zHwAzjUXFkvNmP7p9A6AlwG85+XxLQAGq2quiIwCMA3AALfHh6pq5TsYwsDhgmJM/t9qfPXnjgqP1YuNxje3no6OzRoYSEZkf20T6uLRC3rgQddBL68u2IybhnQwnCp8GS0VVf1ZRJIrefx3t5uLALQJdqZQUlU8+t06vP3blgqP3TasI+4YfhJPPCQKgasGtsOLczdh3+ECPD1zPS4fkMRztqrJ+D4VV6l852nzV7lx9wDooqrXu25vAZALQAG8rqrTvDxvAoAJAJCUlNRv61bv+ylCpaC4BDd98Afmrd9b4bEf7x6MDon1DaQiIu649yysdtT7UioiMhTAKwDOUNX9rvtaqepOEWkGYA6Af6hqxYtXuzG5o15VsXrHQYyZ+isc5T7yRnVqYdn/DUcMD2kkMkpVkXI/d9yXFzY76n0hIicDeBPAqOOFAgCqutP1514R+RpAfwCVlooJ+w4X4KPF2/DcnI0VHjv35JZ4+bI+PIKLyCJEBB+MH4Ar31oMAPhx3R6c1bW54VThxdKlIiJJAL4CcJWqbnS7vx6AKFU95FoeAWCKoZheTfxyJT5Zur3C/Rf2aY3nL+1tIBERVcX9SpHj303j2oqfjG7+EpGPAQwB0BTAHgAPAagFAKr6moi8CeBiAMd3hBSraqqItAfwteu+GAAfqerjVb1fqDZ/rcw6gPNf/q3C/QPbJ+Djvw/kmglRGOD+lRPCap9KKAW7VFQVHyzeVnpo4nEJ9WLx631DUTfW0iuGRORm3+ECpD42FwDwyhV9MbpnS8OJzLHVPpVwct5Lv2LNzoNl7ltwzxAkN+V1G4jCTdP6cagfF4PDBcW4+cM/In5txVc83CgAjhaWIHni92UK5T+X9UHmU+eyUIjC2OpHzildLn/Ne/KMpVJDRwqK0XXyzDL3LZl0Fs7v1cpQIiIKpK9vPq10eVfeUYNJwgNLpYYue2NRmdsL7x/G618T2UifpPjS5VOfnGcwSXhgqdTA+t0HsTLrxJT0M+84Ey0b1TGYiIiCYcuTo0uXPZ1zRiewVGpg5Au/lC5Pu6ofurRoaDANEQWLiGBUjxYAgP/8uMlwGmtjqVTTfLd5uy7q2xojurcwmIaIgu3VK/uVLvd6ZLbBJNbGUqmG/MJiXPvO0tLbz/2NZ8cTRYLXr3IWS97RIkTSOX7+YKn4qaC4BN0mzyq97b6tlYjs7Ry3LRLuE0/SCSwVP6gqRrntR3lrXCqnXCGKMAvvH1a6nF9YbDCJNbFU/LBm50Fk7DtSepuzlxJFHvcjPN23WpATS8UP5730a+nysv8bbjAJEZm0/tGRpcsH8gsNJrEeloqPCopLSpfbJ9ZDk/pxBtMQkUm1a0WXLveeMsdgEuthqfhowYbs0uXpt55hMAkRWcHGx0aVLucdLTKYxFpYKj664f1lpcv14zi5M1Gki4058eOT562cwFLx0/CuzUxHICKL2PDYiX0rRwp4JBjAUvHbTUM6mo5ARBYRF3Ni30r3h3gkGMBS8Yn7mbN92jY2mISIrGaN2zVXHA6eZc9S8UF+4Ykjv6KieLIjEZ1Qz20fa/sHeJY9S8UHvDAPEVVm7l2DTUewDJaKD56ZucF0BCKysI7N6pcuj3n510pG2h9LxQez1+4xHYGILO6JC3sCAFa4XbgvErFU/DCcc30RkReXD0gqXf5w8VaDScxiqfjh/N6tTEcgIgvr2boRAGDS16sNJzGHpeKHXm0amY5ARBb27T9OTOGUcyQyJ5pkqfihecPapiMQUZjo+2hkTjTJUvFDXAw/LiKq3MqHR5iOYBR/SvqhsMRhOgIRWVzD2rVKl696a7HBJGawVPzA6a2JyBf/PKczAOCXTfsMJwk9looftmQfqXoQEUW8W4aemHh2c/Zhg0lCj6Xih9/SI++3DiKqmbP+/ZPpCCHFUvEDz6wnIl+5X2slkrBU/LB+9yHTEYgoTLhfa+XmD5dVMtJeWCpEREEyYVB7AMCMVbsNJwkdlooPBqQkmI5ARGHogdFdS5ePFZVUMtI+WCo+uLBPa9MRiCjMdXlwpukIIWG0VETkbRHZKyIeZ18Tp/+ISLqIrBSRvm6PjRORTa6vccHMOaJ7i2C+PBHZ2Jw7B5mOEFKm11TeAVDZIRKjAHRyfU0A8CoAiEgCgIcADADQH8BDIhIfrJAJ9WKD9dJEZHOdmjcoXV64eb/BJKFhtFRU9WcAOZUMGQPgPXVaBKCxiLQEcA6AOaqao6q5AOag8nIiIjLusjcWmY4QdKbXVKrSGsB2t9tZrvu83U9EZDnpj48yHSFkrF4q4uE+reT+ii8gMkFE0kQkLTs7u8aBHA6Pb0NE5FVM9IkftVPnpxtMEnxWL5UsAG3dbrcBsLOS+ytQ1WmqmqqqqYmJiTUONH2Fx7chIqrU8atCPjtrg+EkwWX1UpkO4GrXUWADAeSp6i4AswCMEJF41w76Ea77gu6uz5aH4m2IyGbcrwppZ6YPKf4YwEIAnUUkS0TGi8iNInKja8gMABkA0gG8AeBmAFDVHACPAljq+priui9oGtaOAQBw6xcR1dQ/Pv7TdISgMX3012Wq2lJVa6lqG1V9S1VfU9XXXI+rqt6iqh1Utaeqprk9921V7ej6+m+wsz5/ae9gvwUR2dzlA5IAAN/aeDO61Td/WcawLs1MRyCiMPfEhT1NRwg6loqPRDwdcEZEVD0XTP3NdISgYKkQEYXQbWd1AgAs337AcJLgYKlUw5GCYtMRiChM3XX2SaYjBBVLpRoe+36d6QhEZAPX/neJ6QgBx1Kpho+XbDMdgYjC2BWuo8Dmb6j5LB9Ww1Lxw42DO5iOQEQ28LiNjwJjqfjh7hH23hZKRKH30DceLycVtlgqfqgVzY+LiALjzE5NAQDvLtxqOElg8ackEZEB748fYDpCULBUqil972HTEYjIJr5YlmU6QsCwVKrpGhseCkhEZtzz+QrTEQKGpeKnOrWiAQBZuUcNJyGicLfxMftdEZKl4qcP/27P7aBEFHqxMSd+BG/PyTeYJHBYKn7qmxRvOgIR2dCZz8w3HSEgWCpERAbNumOQ6QgBxVKpgWNFJaYjEFGY69yiQemywwaXlmWp1MCkr+11JiwRmTXo2fDfBMZSqYEv/7DPseVEZM4dw53XWLHDUaUslWp45Yq+piMQkY3cMdw+8wqyVKphdM+WpiMQkU3NXL3LdIQaYakQEVnIjR/8YTpCjbBUaqiw2GE6AhHZwIrJI0xHCAiWSg09NJ1HgBFRzTWqW6t0ubgkfH9ZZanU0MdLtpuOQEQ2c8bT4XtoMUulml4c29t0BCKymduGdQQA7D54zHCS6mOpVNOY3q1NRyAim7lrRGfTEWqMpUJEZEHLtuaajlAtLJUA4BFgRBRoF7/6u+kI1cJSCYCHv11jOgIR2cS8uwebjlAjLJUA+GjxNtMRiMgm2ifWNx2hRlgqNfDy5X1MRyAiG3th7kbTEfzGUqmB805uZToCEdlQ3dhoAMALczcZTuI/lgoRkcX88eDZpiNUG0slQI4W8iqQRBQYtWtFly6H29UgWSoBcvsnf5qOQEQ2dFGYHVpstFREZKSIbBCRdBGZ6OHx50Vkuetro4gccHusxO2x6aFNXtHstXtMRyAiGzn3ZOd1m5ZvP1DFSGsxVioiEg1gKoBRALoBuExEurmPUdU7VbW3qvYG8BKAr9wePnr8MVU9P2TBy/n8xlNNvTUR2djLl4Xn0aUm11T6A0hX1QxVLQTwCYAxlYy/DMDHIUnmh1OSE0xHICIbEpHS5fzCYoNJ/GOyVFoDcJ83Pst1XwUi0g5ACoB5bnfXFpE0EVkkIhcELyYRkVmpj801HcFnJktFPNzn7TCHsQC+UFX3Q6ySVDUVwOUAXhCRDh7fRGSCq3zSsrOza5a4Cut2HQzq6xNRZLlxsPPHWn4YHV1qslSyALR1u90GwE4vY8ei3KYvVd3p+jMDwAIAHjdAquo0VU1V1dTExMSaZq7UqBd/CerrE1FkmTiqi+kIfjNZKksBdBKRFBGJhbM4KhzFJSKdAcQDWOh2X7yIxLmWmwI4HcDakKT2YFiXZqbemogixJGC8NivYqxUVLUYwK0AZgFYB+AzVV0jIlNExP1orssAfKKq7pvGugJIE5EVAOYDeEpVjZXK1Mv7mnprIooQpzweHvtVYky+uarOADCj3H2Ty91+2MPzfgfQM6jh/FAnNrrqQURE1XDDoPZ4/eeMsNmvwjPqA6zsChURUc2E234VlkqAPTtrg+kIRGQj4Xa+CkslwF5ZsNl0BCKyqUHPzDcdoUoslQCZdlU/0xGIyKbGnuI8+2Lf4ULDSarGUgmQEd1bmI5ARDb15EWWOS6pSiwVIiKLc9+vYvXrq7BUgiAj+7DpCERkUxPeX2Y6QqVYKkEw8gVO10JEgdU+sR4AYO46a1+7iaUSQKe2bwIAKCxxGE5CRHYz47YzTUfwCUslgN4cl2o6AhHZlPt1662MpRJA9eKMznpDRBHi/YWZpiN4xVIJkhKLH6FBROHrwW/WmI7gFUslSO7/aqXpCERkM3PvGmw6QpVYKkHyWVqW6QhEZDMdm9U3HaFKLJUA++H28DhCg4jC26Y9h0xH8IilEmBdWzY0HYGIIsDZz/9sOoJHLBUiojDywGhrX1+FpRJE037mNPhEFFjjz2hvOkKlWCpB9MSM9aYjEJHNREdZe3JJlkoQvHPtKaYjEFEEeGi69c5XYakEwZDOzUxHIKII8P6iraYjVMBSISIKM1Y+CZKlEmQzV+82HYGIbKaDaxp8K2KpBNmNH1j7gjpEFH7crwS5btdBg0kqYqkEyWMX9DAdgYgiwKgXrXVRwCpLRUTuFJE2oQhjJ1cObGc6AhHZ2JUDk0xH8MiXNZWGAGaJyC8icouINA92KCIiqtydw08yHcGjKktFVR9R1e4AbgHQCsBPIjI36MlsZFHGftMRiMhmmtSPMx3BI3/2qewFsBvAfgA8EcMPY6ctMh2BiGzs06XbTEco5cs+lZtEZAGAHwE0BfB3VT052MHs4J/ndDYdgYgiwH1frjIdoZQvayrtANyhqt1V9SFVXRvsUHZxy9COpiMQkY3dP8p6Mxb7sk9loqouD0UYIiLy3cgeLUxHqIDnqYTI3LV7TEcgIptpG1+3dLmoxGEwyQkslRC5/r000xGIyGai3KbBt8qMxSyVIHvlir6mIxBRBPhosTWOAGOpBNnoni1NRyAiG7tvpLV21rNUiIjCWGpyvOkIZRgtFREZKSIbRCRdRCZ6ePwaEckWkeWur+vdHhsnIptcX+NCm7x6Xvpxk+kIRGQznVs0KF0+XFBsMImTsVIRkWgAUwGMAtANwGUi0s3D0E9Vtbfr603XcxMAPARgAID+AB4SEWvVtQf/nrPRdAQispmGtWuVLt/3xUqDSZxMrqn0B5CuqhmqWgjgEwBjfHzuOQDmqGqOquYCmANgZJBy1tiM2840HYGIIsD3q3aZjmC0VFoD2O52O8t1X3kXi8hKEflCRNr6+VyIyAQRSRORtOzs7EDk9lu3Vg2NvC8RRYZLU9tWPShETJaKeLhPy93+FkCya66xuQDe9eO5zjtVp6lqqqqmJiYmVjtsoFjlBCUiso+uLRtUPShETJZKFgD3em0DYKf7AFXdr6oFrptvAOjn63Otavy7PAmSiAIrqcmJM+sLiksMJjFbKksBdBKRFBGJBTAWwHT3ASLifpLH+QDWuZZnARghIvGuHfQjXPdZVstGtQEAP280swmOiOwrKeFEqTzyrdk5f42ViqoWA7gVzjJYB+AzVV0jIlNE5HzXsNtEZI2IrABwG4BrXM/NAfAonMW0FMAU132WNfOOQaYjEJFNtXGbA8z0mfUxJt9cVWcAmFHuvsluy/cDuN/Lc98G8HZQAwZQozq1qh5ERFQNtWtFIzY6CoUW2GfLM+oN2J6TbzoCEdlMr7aNTEcAwFIx4sxn5puOQEQ209Ztv4rD4fFg2JBgqYTQ7Wd1Mh2BiGzK/doqb/26xVgOlkoI3TGcpUJEweF+BNjjM9ZVMjK4WCohJHLinE1Vc6unRGQ/7ueqmMRSMWTKd2aPJScie2kTX8d0BAAsFWP++1um6QhEZCN1Y42eIVKKpRJiSycNNx2BiGwoLqbsj/NdeUeN5GCphFhigzjTEYjIhmKjy/44v+OT5UZysFQMWr/7oOkIRGQTUVGCWtEnDgZavMXMzFUsFYNGvvCL6QhEZCOx0VH4W2oboxlYKgY8cWFP0xGIyIbiakUjLibaaAaWigGXD0gyHYGIbCg2Oiqir6dCAIotMKsoEdlDbEwUCotP/Ewxcf0mlophV721xHQEIrKJuJgoFLiVym2f/BnyDCwVQ7q1bAgAWJix33ASIrKL42sqtwztAAA4kF8U8gwsFUO+vuU00xGIyGZiXWsqtw41N3ktS8UQ00doEJH9xLnWVOrEmvv5wlKxgJfnbTIdgYhsYO/BAjQ0fOlylooF/Gv2RtMRiCjM5R0tQsa+I+hd7rLC+YXFIc3BUjHozwfPNh2BiGxi9Y48AMDJbRqXuf/VBZtDmoOlYlB8vVjTEYjIJpZvPwAAOLmNc03l+DxgL81LD2kOlopFfJ623XQEIgpjK7MOILlJXTSu6/xl9ZHzexjJwVKxiH9+sdJ0BCIKYyu256FX2xObvkxNLMlSMeyXe4eajkBEYW7PwWPYffBYmf0pMdFmfryzVAxrm1DXdAQiCnMrXPtTyh/5ZQJLxULmr99rOgIRhaGVWXmIjhJ0a+m5VI4VhW7mYpaKhVz7zlLTEYgoDK3IOoDOzRt4PZP+sxAeCMRSsYBZdwwyHYGIwpSqYmVWHnpVsunrqR/WhywPS8UCOrdoYDoCEYWprfvzkXe0CL3KnfToLr+Qm78i1k8GLqpDROFr/gbnvtjyZ9IDwMjuLUIdh6ViNePe5kW7iMg36XsP4ZmZGzCwfQK6eNjiMbZ/25BnYqlYxLy7B5uOQERh5FhRCW796E/UjY3Gi2P7ICpKKow5vWPTkOeKCfk7kkftE+ubjkBEYWTKd2uxfvchvHPtKWjesLbHMbUMnADJNRULCvWsokQUXr5buRMfLd6GGwa3x5DOzUzHKYOlYiHiWnt9emboDv8jovCybX8+7v9yFfokNcY9IzqbjlOB0VIRkZEiskFE0kVkoofH7xKRtSKyUkR+FJF2bo+ViMhy19f00CYPjjWPnGM6AhFZWGGxA7d+/AdEgJcu6+PX5i1VDWKyE4ztUxGRaABTAZwNIAvAUhGZrqpr3Yb9CSBVVfNF5CYAzwC41PXYUVXtHdLQQVY39sQ/h6pCpOKONyKKXE/PXI+VWXl47cp+aBPv37yBqie2hgSTyTWV/gDSVTVDVQsBfAJgjPsAVZ2vqvmum4sAmJnL2YDL31hsOgIRWcjctXvw1q9bMO7UdhjZw//zT/KOFgUhVUUmS6U1APcJabJc93kzHsAPbrdri0iaiCwSkQu8PUlEJrjGpWVnW//EwuFdnTvdFmbsN5yEiKxi54GjuOeLFejeqiHuH921Wq9xJETXqjdZKp5WxDxu9BORKwGkAnjW7e4kVU0FcDmAF0Skg6fnquo0VU1V1dTExMSaZg66N65ONR2BiCykuMSB2z7+E0XFDrx8eV/UruV50siqHDxq/1LJAuB+umcbADvLDxKR4QAmAThfVQuO36+qO11/ZgBYAKBPMMOGivt+lNwjhQaTEJEVvDB3E9K25uKJi3oipWm9ar9OYYkjgKm8M1kqSwF0EpEUEYkFMBZAmaO4RKQPgNfhLJS9bvfHi0ica7kpgNMBuO/gt4U+j84xHYGIDPppYzamLkjH31LbYEzvyvYOVK12rdD8uDdWKqpaDOBWALMArAPwmaquEZEpInK+a9izAOoD+LzcocNdAaSJyAoA8wE8Ve6osbD2/KW9TEcgIoNyjxRi8jercd07S9ExsT4ePr97jV+zflxoDvY1Ok2Lqs4AMKPcfZPdlod7ed7vAHoGN505F/Zpgzs/XWE6BhGFWFGJAx8s2ooX5m7CoWNFuHJgO9w5/KQypxtUV8tGdQKQsGqc+8vivliWhUv6RcyR1EQRa8GGvXjs+3VI33sYZ3RsigfP6xbQay1Fe5hwMhhYKhZ3z+crWCpENrY5+zAe/34d5q3fi+QmdfHG1akY3rVZ2J78zFKxqHVTRqLr5JmmYxBRkOQdLcJ/ftyEd3/PRJ1a0XhgdBeMOy0ZcTHVO2TYk/wQnZvijqViUXViT3xjHS4oDtlONiIKrhKH4pOl2/Dv2RuRm1+IS1Pb4u4RnZHYIC7g7zVn7Z6Av2ZV+JMqDPR4aBYynzrXdAwiqqHfN+/DlG+d10Hpn5KAyed1Q4/WjYL2fp+lba96UICxVCzsxbG9cfsny03HIKIa2rY/H4/PWItZa/agdeM6eOWKvhjVo0XQ95v8lh766Z5YKhY2pndrlgpRGDtcUIyp89Px1i9bEBMtuGfESbj+zPbVnmqluprWjw3Ze7FUwsTU+em4ZWhH0zGIyAcOh+KLP7Lw7KwNyD5UgIv6tMa9I7ugRSPPl/0NtvtGdgnZe7FULK5p/VjsO1yIZ2dtYKkQhYG0zBw88u1arNqRhz5JjTHtqn7okxQf8hxHCk4c+VXTKV78wVKxuMUPDEeHB2ZUPZCIjFFVrNl5EK//nIFvV+xEi4a18cKlvXF+r1aICtFJh+VN+npV6XJsTOhm5GKpWJz7WbBb9x9BuybVn6WUiAJHVbEyKw8zVu/CD6t2Y1tOPuJionDbsI64cUiHgEytUhP/W6Wj7HMAABgXSURBVF5h0veQYKmEkcHPLuChxUQGORyK5VkH8MOqXZixajd2HDiKmCjBaR2b4uYhHTCiewsk1AvdTnFv0jJzSpdfvaJvSN+bpRIG5t8zBEP/tcB0DKKI5HAolm3LxYxVuzBz9W7syjuGWtGCMzsl4o7hnXB2t+ZoXNd8kbi75LWFpcvDXFeTDRWWShhwvzBPfmGx8dVqIrsrcSiWZubgh1W78MPq3dh7qACxMVEYfFIi7h3ZGcO6NEejOrVMx/To0LET16J/4sKeAZ32xRf86RRmuk3m2fVEwVBc4sDiLTmYsWoXZq3ZjX2HCxEXE4WhnZthVM8WGNalGRrUtmaRuOv58OzS5cv6t61kZHCwVMLE++P746q3lpiOQWQrRSUOLNy8Hz+s3oVZa/Yg50gh6tSKxrCuzTC6R0sM6ZyIemE07577Jcjn3jXIyEzH4fNpRbgzOyWWLjscauwwRaJwV1jswG/p+zBj1S7MXrsHeUeLUC82Gmd1bY7RPVtg8EnNykzoGk7cL0HesVngrsXiD5ZKGEp9fC7+ePBs0zGIwsaxohL8uslZJHPW7cGhY8VoEBeDs7s1x6ieLXFmp6Yhnzol0OZv2Fu6HOojvtyxVMLIQ3/phke+XYsct1VcIvLsWFEJFmzIxg+rd+HHdXtxuKAYDWvH4JzuLXBuz5Y4rWOTkO/EDpYSh+La/y4tvX12t+bGsrBUwsi1p6fgkW/Xmo5BZFlHCorx08ZszFi1C/PW70V+YQni69bCeSe3xKieLXFq+yYhPbs8VG78YFnp8uieLRATbe7vyFIJU397fSE+u+FU0zGIjDpSUIy0rblYlLEfizL2Y1VWHoodiqb1Y3Fhn9YY3bMlBqQkGP0hG2wZ2YfLXIzrqYtPNpiGpRJ2bhjUHq//nIElW3KqHkxkM4cLipGWmYNFGTnOEtmRhxKHIiZK0KttY9wwuD3O6JiI/ikJZaY4squC4hIM+/dPpbc/mTAQDQ0f9sxSCTP3j+6K13/OAOCce8jEIYNEoXLoWJHbmkgOVrtKpFa0oFebxrhpcAcMbN8Efds1jriTgg8eK8LJbuek3HX2SRjYvonBRE6R9a9gM8Of+wk/3j3EdAyigDl4rAhpmTlY7LYm4lCgVrSgd9vGuHmIq0SS4sP2sN9AyNx3BEPcpm7q1Kw+/jHMGpfGYKmEocnndcOU79Zic/YR01GIauTgsSIs3ZKDxVucJbLarUT6tI3HrUM7YgBLpIyZq3eX2TEPAHPuGmwoTUUslTB03RkpmPKd8yiw4hKHrXdCkr3kHT1eIs7NWWt2OkskNjoKvZMa49ZhnTAwJQF9WCIe3fTBMvywenfp7fpxMZY7Z42lEuY6TvqBc4GRZeXlF2FJZg4WZ+zHoi37sWbnQairRPokNcY/hnXCgPYJ6JsUH/YnHwaTqiLl/rIX6zslOR6f3XCq5farslTC1HvX9cfVb3MuMLKWvPwiLN6yv3Rz1tpdrhKJiULfpMa4bVgnDGzfBH2SGrNEfLRl35EKl774ZMJAS+yU94SlEqYGnXRiLrAD+YWWu54D2Z+qYseBo0jLzMXSzBws25qLDXsOlZZIv6R43H6Ws0R6t2WJVMdl0xZhYcb+0tu3DO2Au87ubOnDpVkqNtB7yhxuAqOgK3Eo1u06iLTMHKRtzUVaZi52HzwGwLltv09SY4zq0RID2yegF0ukRrIPFeCUx+eWuW/F5BFoVNf6U++zVMLYL/cOxZnPzDcdg2zqSEExlm8/gLTMXKRtzcEfW3NxpLAEANCyUW2ckpKA1HbxSE2OR5cWDS3923M4GTttIRZlnDi5+ZzuzfH6VakGE/mHpRLG2ibULV1esiUH/VMSDKahcLf34DGkbXVuykrLzMXaXQdR4lCIAJ2bN8CFfVvjlOQEpCYnoHXjOqbj2k5aZk6ZywADwKcTBmKARfedeMNSsYm/vb6Qm8DIZw6HIj37sHMtxLU5a1tOPgCgdq0o9G7rPFs9NTkefZLiLXvpXDs4dKyozNUaj1twzxAku11KPFywVMLc5idGo8MDM6oeSBHtWFEJVu3Ic+5Qz8xF2tZc5B11Xsu8af1Y9GsXj6tPbYd+7eLRvVUjW87kazUOh+L5uRvx0rz0Mve3aFgbM24/Ewn1wvPgG5ZKmHPfjn3uf37B97edaTANWcXxiReXbHF+rczKQ2GJAwDQPrEeRnZvgX7J8TglOQHJTepa7lwHu5u9ZjcmvL+swv1f3nQa+rWLN5AocFgqNjBlTHdM/mYN1uw8aDoKGXIgvxBLM3OxOGM/lmTmlE53Eh0l6Nm6Ea45PRmp7eLRr108mtSPMx03IqkqPlqyDZO+Xl3hsScu7InL+re1RbkbLRURGQngRQDRAN5U1afKPR4H4D0A/QDsB3Cpqma6HrsfwHgAJQBuU9VZIYxuKVefmozJ36wBAOw4cJQ7USPA3kPHsHRLLpa4TjRcv/sQAOf5Ib3bNsYtQztiQIrzJMN6cfzd0aTsQwX4+3tpWL79QIXHXhzbG2N6tzaQKniMfbeJSDSAqQDOBpAFYKmITFdV90sbjgeQq6odRWQsgKcBXCoi3QCMBdAdQCsAc0XkJFUtCe3fwnpOf2oed9jb0I4DR7Fky34scU2+mOGaTLRubDT6tYvHuT1bYkD7Jji5TSOeH2IBqooPFm3Fg65f9sr76PoBOLVDE1usmZRn8leY/gDSVTUDAETkEwBjALiXyhgAD7uWvwDwsjj/FcYA+ERVCwBsEZF01+uVPR4vgqx+5Bz0eChiV9ZsJ3PfkdLpTpZsyUFW7lEAQIPaMeifnIBLU9uif0oCerRuhFqcUNQyPE2p4u7rm09Dn6Tw3mdSFZOl0hrAdrfbWQAGeBujqsUikgegiev+ReWea691SD/Vd9vEcfOHy/DKFf0MpqGauu6dpcjYdwRN6sWif0oCxp+Rgv4pCTzJ0IKKSxy4+LWFWOFh8xbgPHnxpcv6RswRdSZLxdP/DPVxjC/Pdb6AyAQAEwAgKSnJn3xh5+YhHfDKgs2YsWp31YPJ0p64qCea1o9Fh8T6ttxEYgev/7QZT/6w3uNjA1IS8NBfuqNrywYR9+9nslSyALR1u90GwE4vY7JEJAZAIwA5Pj4XAKCq0wBMA4DU1FSPxWMX947sglcWbAYArN6Rhx6tGxlORNVl1RloI90vm7Jx1VueZwdPbRePF8b2Rpv4uh4fjxQmS2UpgE4ikgJgB5w73i8vN2Y6gHFw7iu5BMA8VVURmQ7gIxF5Ds4d9Z0AcB54N+e99Ct32BMFwJ/bcnHhK797fCy+bi3Mu3sI4sP0RMVgMFYqrn0ktwKYBechxW+r6hoRmQIgTVWnA3gLwPuuHfE5cBYPXOM+g3OnfjGAW3jkl9OWJ0eXXsynqMTBnbhE1bBsay4uftVzkQDhM2OwCaJq6y1CZaSmpmpaWprpGEGXPPH70mWurRD55tdN+3DlW4u9Pv7dP86I2E3KIrJMVX2aKplnRdnQr/cNxRlPc0p8oso4HIpZa3bjpg//8Dpm0uiuuOrUdjz3xw8sFRty31E46sVf8MPtnA+MCAAKix146of1ePu3LV7HjD8jBTcMao9mDWuHMJl9sFRs6umLe+K+L1dh3S7OB0aRbfWOPJz30q+Vjpk0uiuuOT2Z+yADgKViU5eekoT7vlwFAHhx7ibcPryT4UREoZGVm48Lpv6OfYcLKh33y71Dy1zojgKDpWJjo3u2wIxVu/H83I0sFbKt3XnHcOErv2FX3rFKx/37r71wUd/WEXcyYqixVGzslSv6lR4J9sbPGfj7oPaGExHV3N6DxzDwyR/hqOLA1Yv6tsbjF/REnVjuZA8llorNDTopET9vzMbjM9axVCgs7co7ilOfnFfluOFdm+HJi05GYgNeL8YklorNvXdd/9K1lXd+24JrTk8xnIiocpuzD+Osf/9U5bhrTkvGxFFdeLivxbBUIkD/lAQs2ZKDh79dy1Ihy1mamYO/vlb1VSsmje6K8WekIIqzNFsaSyUCfHbDqaVrK68sSMfNQzoaTkSR7ItlWbjn8xVVjrPjVREjAUslQpzVpRl+XL8Xz8zcwFKhkCkuceDuz1fgm+UeJxEv493r+mPwSYkhSEXBxFKJEG9dc0rp2spfX/sdn994muFEZEc7DhzF6U9VvVMdAL655XT0ats4yIko1FgqEeThv3TDw9+uxdLMXNNRyCamzk/Hs7M2VDmuaf1YzLlzMKeIjwCcpTjCcAZjqq7tOfk48xnfJiq9bVhH3Hn2STzR0CY4SzF5NfeuQRj+3M8AgLyjRWhUh9eEoIpUFQ9PX4N3F271afz8e4YgpWm9IKeicMBSiTAdmzUoXe71yGyurRAAYO3Ogxj9n198GnvTkA6495zOXAshj1gqESjjidFo/4Dz6pA8ITLylDgUE95Lw4/r9/o0fvEDZ6E5p4EnH7FUIlBUlKBXm0ZYkZXHEyIjwO/p+3D5m96vaOhuypjuuPrU5OAGIlvjjvoIxp329lNQXIKx0xbhz20HqhzbomFtzL5rEBrW5n41qhx31JNPfp84DKe5zimYuXoXRvZoaTgR+cuftZDXruzLf2MKOpZKBGvVuA5EAFXgxg/+4NqKxTkcignvp2Huuqr3hfRo3RBf3HgaJ1ukkOPmL+JmMItav/sgRr7g2xFZb16diuHdmgc5EUUqbv4iv6x4aAR6PTIbAHDFm4vw4fUDDSeKPKqKp2duwGs/ba5ybJv4Oph1xyDUi+N/X7IeflcSGtWphWtPT8Z/f8vEb+n7sSorDz3bNDIdy9aOFpZg2L8XVHkJXAB48qKeuKx/UghSEdUcN39RKffNYOmPj0JMdJTBNPbi69ULAWDJpLPQrAHPCyHr8GfzF0uFyuD+lcBYsf0Axkz9rcpx553cEi9d1odnp5OlsVS8YKlUrcSh6OA62x5gsfhq7to9uP69qr+3Hjm/O8adlhz8QEQBxB31VG3RUYJf7h1aOhtt8sTvWSwezF+/F9e+s7TKcbzwFEUalgpV0DahLj4YPwBXvuU8qY7FAqzKysNfXv61ynGz7xyEk5o3qHIckV1x8xd59dUfWbjrsxPXEo+kYtl54GjpbAOVmXf3YLRPrB+CRETmcPMXBcRFfdsgOkpw+yfLAdh7jeVIQTF6PjwLjip+x/rq5tPQNyk+NKGIwhBLhSo1pndrtGhYG5dOWwTAWSwrJo9Ao7rhPQmhw6G48NXfsWJ75RMvvnJFX4zuyfmyiHzFzV/kk9wjhejz6JzS2+PPSMGD53UzmMh/k75ehQ8Xb6t0zD0jTsKtwzqFKBFReOAhxV6wVGrO/TwWwNonSb66YDOenrm+0jGje7bA1Mv78jwRokqwVLxgqQTG+HeWlrlqYJ1a0Vj36EiDiZy+Wb6jdP+PN0kJdfHj3YNRy6JFSGRFLBUvWCqBc6yoBF0enFnh/lDuyP9s6Xbc++XKKsetengEGvBCVETVZvlSEZEEAJ8CSAaQCeBvqppbbkxvAK8CaAigBMDjqvqp67F3AAwGkOcafo2qVv4rKlgqwTBv/R5c907Fz/TekZ1x85COAXufrfuPYMi/FsCXb9dF95+FFo04dxZRoIRDqTwDIEdVnxKRiQDiVfW+cmNOAqCquklEWgFYBqCrqh5wlcp3qvqFP+/LUgmenzZmY9zbS7w+3rFZfUy9vC86t/B+YmBhsQP/+3MHJv1vFYpKfP++nH/PEKQ0redXXiLyXTiUygYAQ1R1l4i0BLBAVTtX8ZwVAC5xlcw7YKlYkqpi8LMLsC0nPyiv/5derfDCpb0RHcUd60ShEg4nPzZX1V0A4CqWZpUNFpH+AGIBuF/B6HERmQzgRwATVbUgaGnJZyKCn+8dWnp754GjGPKvBSgsdvj1Os0axOH1q/qhD080JAorQVtTEZG5AFp4eGgSgHdVtbHb2FxV9fjT4/iaDIBxqrrI7b7dcBbNNACbVXWKl+dPADABAJKSkvpt3bq12n8nIqJIZIk1FVUd7u0xEdkjIi3dNn/t9TKuIYDvAfzf8UJxvfYu12KBiPwXwD2V5JgGZ/EgNTU1cg51IyIywNTB+tMBjHMtjwPwTfkBIhIL4GsA76nq5+Uea+n6UwBcAGB1UNMSEZFPTJXKUwDOFpFNAM523YaIpIrIm64xfwMwCMA1IrLc9dXb9diHIrIKwCoATQE8Ftr4RETkCU9+JCKiSvmzT4VzVRARUcCwVIiIKGBYKkREFDAsFSIiChiWChERBQxLhYiIAoalQkREAcNSISKigGGpEBFRwETUGfUikg3A32mKmwLYF4Q4gWDlbADz1RTz1YyV81k5G1AxXztVTfTliRFVKtUhImm+Tk8QalbOBjBfTTFfzVg5n5WzATXLx81fREQUMCwVIiIKGJZK1aaZDlAJK2cDmK+mmK9mrJzPytmAGuTjPhUiIgoYrqkQEVHAsFTKEZFnRWS9iKwUka9FpLGXcSNFZIOIpIvIxBBl+6uIrBERh4h4PTJDRDJFZJXrapkhuyqZH/lC/tm53jdBROaIyCbXn/FexpW4XW10eghyVfp5iEiciHzqenyxiCQHO5Mf2a4RkWy3z+v6UGVzvf/bIrJXRDxeUlyc/uPKv1JE+loo2xARyXP77CaHKpvr/duKyHwRWef6f3u7hzH+f36qyi+3LwAjAMS4lp8G8LSHMdEANgNoDyAWwAoA3UKQrSuAzgAWAEitZFwmgKYGPrsq85n67Fzv/QyAia7liZ7+bV2PHQ7hZ1bl5wHgZgCvuZbHAvjUQtmuAfByqL/X3N5/EIC+AFZ7eXw0gB8ACICBABZbKNsQAN8Z/OxaAujrWm4AYKOHf1+/Pz+uqZSjqrNVtdh1cxGANh6G9QeQrqoZqloI4BMAY0KQbZ2qbgj2+1SXj/mMfHYuYwC861p+F8AFIXrfyvjyebjn/gLAWSIiFslmlKr+DCCnkiFjALynTosANBaRlhbJZpSq7lLVP1zLhwCsA9C63DC/Pz+WSuWug7Oly2sNYLvb7SxU/McwSQHMFpFlIjLBdJhyTH52zVV1F+D8DwWgmZdxtUUkTUQWiUiwi8eXz6N0jOsXnjwATYKcy9dsAHCxa9PIFyLSNgS5/GH1/6unisgKEflBRLqbCuHapNoHwOJyD/n9+cUEMli4EJG5AFp4eGiSqn7jGjMJQDGADz29hIf7AnIYnS/ZfHC6qu4UkWYA5ojIetdvTVbIF7TPDqg8nx8vk+T6/NoDmCciq1R1c2ASVuDL5xHUz6wSvrzvtwA+VtUCEbkRzjWqYUFP5jtTn50v/oBz+pPDIjIawP8AdAp1CBGpD+BLAHeo6sHyD3t4SqWfX0SWiqoOr+xxERkH4DwAZ6lrw2I5WQDcfyNrA2BnKLL5+Bo7XX/uFZGv4dyMEZBSCUC+oH12QOX5RGSPiLRU1V2uVfi9Xl7j+OeXISIL4PwNLlil4svncXxMlojEAGiE0GxWqTKbqu53u/kGnPshrSSo32814f4DXFVniMgrItJUVUM2J5iI1IKzUD5U1a88DPH78+Pmr3JEZCSA+wCcr6r5XoYtBdBJRFJEJBbOnadBP0rIFyJST0QaHF+G88ADj0efGGLys5sOYJxreRyACmtWIhIvInGu5aYATgewNoiZfPk83HNfAmCel192Qp6t3Pb18+HcLm8l0wFc7TqKaSCAvOObQE0TkRbH942JSH84fx7vr/xZAX1/AfAWgHWq+pyXYf5/fqaOPLDqF4B0OLchLnd9HT/qphWAGW7jRsN5tMRmODf9hCLbhXD+5lAAYA+AWeWzwXmkzgrX15pQZfM1n6nPzvW+TQD8CGCT688E1/2pAN50LZ8GYJXr81sFYHwIclX4PABMgfMXGwCoDeBz1/fmEgDtQ/iZVZXtSdf32QoA8wF0CVU21/t/DGAXgCLX9954ADcCuNH1uACY6sq/CpUcNWkg261un90iAKeF+LM7A85NWSvdft6NrunnxzPqiYgoYLj5i4iIAoalQkREAcNSISKigGGpEBFRwLBUiIgoYFgqREQUMCwVIiIKGJYKkUEikux+vQ0RuUdEHjYYiahGWCpERBQwLBUiIgoYlgqRWcUo+/+wtqkgRIHAUiEyaw+AZiLSxDU78nmmAxHVREReT4XIKlS1SESmwHnFvS0A1huORFQjnKWYiIgChpu/iIgoYFgqREQUMCwVIiIKGJYKEREFDEuFiIgChqVCREQBw1IhIqKAYakQEVHA/D8S1WVf8qKYAgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret.plot(x='u', y='v', xlabel='u', ylabel='v', figsize=[6.0, 6.0])" ] } ], "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.7.1" } }, "nbformat": 4, "nbformat_minor": 1 }