{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "import numpy as np" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Nicholson & Bailey (1935) introduced this model to study oscillations in a system in which a parasitoid controls the host population.\n", "\n", "$$\n", "{N_{t+1}=\\lambda N_{t} \\exp \\left[-a P_{t}\\right]} \\\\\n", "{P_{t+1}=b N_{t}\\left(1-\\exp \\left[-a P_{t}\\right]\\right)}\n", "$$\n", "\n", "with the possible extension for logistic growth of prey:\n", "$$\n", "\\lambda_t = \\operatorname{Exp}\\left(r\\left(1-\\frac{N_{t}}{K}\\right)\\right)\n", "$$" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [], "source": [ "def step(N, P, λ, a, b):\n", " return (\n", " λ * N * np.exp(-a * P),\n", " b * N * (1 - np.exp(-a * P))\n", " )" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "def multistep(N0, P0, λ, a, b, tmax):\n", " N, P = np.empty((2, tmax))\n", " N[0], P[0] = N0, P0\n", " for t in range(1, tmax):\n", " N[t], P[t] = step(N[t-1], P[t-1], λ, a, b)\n", " return N, P" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAw0AAAIPCAYAAAAxXjJfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAWJQAAFiUBSVIk8AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3Xl8nWd54P3frdWWLC+SvK9kT7MQ20lIWAoN07KUsiVM006XaWn70pJhSqGdQoGGvl1oB1qGdKGl7wBT+pIAYWlaYKZpgAANJLGTmGxOQux4jy3LlrVY67nnj3Ok8xxZ25Geo3OO9Pt+Pvqc53nOs9yy5eRcuu/rukKMEUmSJEmaTE25ByBJkiSpshk0SJIkSZqSQYMkSZKkKRk0SJIkSZqSQYMkSZKkKRk0SJIkSZqSQYMkSZKkKRk0SJIkSZqSQYMkSZKkKRk0SJIkSZqSQYMkSZKkKRk0SJIkSZqSQYMkSZKkKRk0SJIkSZqSQYMkSZKkKRk0SJIkSZpSXbkHsBiFEPYBy4H9ZR6KJEmSFrZtwJkY4wvmchODhvJYvnTp0tZLL720tdwDkSRJ0sL1xBNPcPbs2Tnfx6ChPPZfeumlrbt27Sr3OCRJkrSA7dy5k927d++f633MaZAkSZI0JYMGSZIkSVMyaJAkSZI0JYMGSZIkSVMyaJAkSZI0JYMGSZIkSVMyaJAkSZI0Jfs0SJIkqWQymQydnZ10d3czMDBAjLHcQ6paIQQaGxtpaWmhtbWVmpr5+/2/QYMkSZJKIpPJcPDgQfr6+so9lAUhxkh/fz/9/f309vayefPmeQscDBokSZJUEp2dnfT19VFXV8e6detobm6e19+OLzSZTIbe3l6OHTtGX18fnZ2dtLe3z8uz/VuTJElSSXR3dwOwbt06WlpaDBjmqKamhpaWFtatWwfk/3zn5dnz9iRJkiQtKgMDAwA0NzeXeSQLy+if5+if73wwaJAkSVJJjCY9O8OQrhACwLwmlfs3KEmSJFWR0aBhPhk0SJIkSZqSQYMkSZKkKRk0SJIkLRAxRkYyNk9T+gwaJEmSFoBHD3ex8w/v5pUf+SYnuuevqo6mF0IghMDWrVvp7++f8Jxt27YRQmB4eHieRzczBg2SJEkLwEf+z146ewfZf7KPrz16tNzD0QQOHDjARz/60XIPY1YMGiRJkqrc8e5+7n26Y2z/+TMT/zZb5bNq1SpaW1v5kz/5Ezo6Oqa/oMIYNEiSJFW5rzx0pCCXobN3sIyj0USampp4//vfz5kzZ/jgBz9Y7uEUzaBBkiSpisUYuXP3oYJjHT0GDZXo7W9/O+effz5/+7d/y1NPPVXu4RTFoEGSJKmKPX70DE8e6y445kxDZaqvr+dDH/oQQ0ND/O7v/m65h1OUunIPQJIkSbN3567D5xw72VMd1ZO2/e6/lHsIM7b/Qz+Zyn1uuukmrr/+er70pS/xne98h5e+9KWp3LfUnGmQJEmqUkMjGb7y8ARBgzMNFe0jH/kIAO9617uIsTr6ahg0SJIkValv7T0xFiCsXd5IbU0AoLt/mIHhkXIOTVO4/vrruemmm7j//vv53Oc+V+7hzIjLkyRJkqpUMgH6Tds38YVdh+jILU061TvEuhW15RrajKS15KcafehDH+IrX/kK73nPe3jTm95U7uFMy5kGSZKkKnS6b5B/e+L42P6NOzbSvqxhbL+jSvIaFqvzzz+f3/iN32Dfvn3cdttt5R7OtAwaJEmSqtBde44yOJIB4MpNK7hwbQutzfmgwbyGyveBD3yAlStX8kd/9Ef09PSUezhTMmiQJEmqQnfuyi9NunHHJgDaljWOHevsdaah0rW2tvLe976XU6dOcfLkyXIPZ0oGDZIkSVXmhyd6ePjgaQDqawOvf+EGANqSMw02eKsK73jHO9i2bVu5hzEtE6ElSZKqTHKW4YZL1rAqFyy0uTypIk1VVrWxsZF9+/bN42hmx5kGSZKkKpLJRL70UL43w+jSJIDWZcmZBpcnKT0GDZIkSVXkvmdPcrSrH4BVTfW84uI1Y++1NSdzGpxpUHoMGiRJkqpIcmnSG67aSENd/uNcW0HJVYMGpcegQZIkqUr0DAzztUePje0nlyZBYU6DMw1Kk0GDJElSlfjaD45ydmgEgIvWLuPyjcsL3k+WXDWnQWkyaJAkSaoSX9xdmAAdQih4f/mSOuprs8d6B0fozwUY0lwZNEiSJFWBQ6f6uO/ZbAOwmgBv3L7xnHNCCHaFVkkYNEiSJFWBLyVmGV564WrWLl8y4XmtzS5RUvoMGiRJkipcjJEvFvRmOHeWYVT7MmcalD6DBkmSpAq3+8Ap9nX0AtDSWMerLls36bkFy5Msu6qUGDRIkiRVuDsTS5N+8sr1LKmvnfTcwgZvLk9SOgwaJEmSKlj/0Aj//MiRsf0bd26a4uzCBm/ONCgtBg2SJEkV7O4nnudM/zAAW1qbuHrrqinPb7N6kkrAoEGSJKmC3bnr0Nj2m3dsPKc3w3g2eFMpGDRIkiRVqOPd/dz7dMfY/o07pl6aBIWJ0J3ONFSEEELBV21tLe3t7dxwww384z/+Y7mHNyN15R6AJEmSJvaVh44wkokAXPuCVja3Nk17TbLkaoc5DRXl93//9wEYGhpi7969fPnLX+Yb3/gGu3bt4s///M/LPLqpGTRIkiRVqDt355cm3TSDWQYYV3LV6kkV5dZbby3Y/7d/+zd+/Md/nI9+9KO84x3vYNu2bWUZ10y4PEmSJKkCPXakiyePdQOwpL6G11wxeW+GpGWNdTTUZT/i9Q9l6BscLtkYNTevfOUrueSSS4gx8sADD5R7OFMyaJAkSapAd+7K92Z41WXraFlSP6PrQgiFFZRcolTRYswuP5suwb3cDBokSZIqzNBIhq88nA8aZpIAnVTQq8Fk6Ip19913s3fvXkIIXHPNNeUezpTMaZAkSaow39p7YuzD/rrlS3jJBe1FXd/aXCVlV29dUe4RzNytXXO/RS6nIZkIHWPkne98J1u3bp3z/UvJoEGSJKnCfPGhfAL0G7dvpLamuKUr7TZ4q0gf/OAHgexSpJUrV/Kyl72Mt771rfzcz/1cmUc2PYMGSZKkCnK6b5C7Hz8+tn/jjo1F36NgeZI5DRVjNH+hGhk0SJIkVZC79hxlcCQDwJWbVnDh2pai75FcntRZyWVXU1jyo/lhIrQkSVIFuXNXfmlSsQnQo5xpUNoMGiRJkirED0/08PDB0wDU1wZe/8INs7pPmzkNSplBgyRJUoX4YqID9A2XrGFV4sN/MdqWJaonVfLyJFUNgwZJkqQKkMlEvrQ735vhzbNcmgTY3E2pq7qgIYSwP4QQJ/k6Nsk1Lw4hfDWE0BlC6Ash7Akh/GYIoXaK57wuhPDNEEJXCKEnhPD9EMIvlu47kyRJi9l9z57kSFc/AKua6vmxi9fM+l7jm7tVc9WehSDGWPV/B9VaPakL+OgEx3vGHwghvAG4E+gH7gA6gZ8C/gJ4CfCWCa65BbgNOAl8BhgEbgI+FUK4Isb47nS+DUmSpKxkAvQbrtpIQ93sf7fb1FDHkvoa+ocyDA5n6BkYpmVJfRrD1CJVrUHD6RjjrdOdFEJYDnwCGAFeEWN8MHf8/cA9wE0hhJtjjLcnrtkGfJhscHF1jHF/7vgfAA8A7woh3BljvC/Nb0iSJC1ePQPDfO3R/IKJ2VZNSmprbuTw6bMAdPYOGjRoTqpueVKRbgJWA7ePBgwAMcZ+4H253V8fd80vA43AX44GDLlrTgF/nNt9W6kGLEmSFp+vP3qMs0MjAFy0dhmXb1w+53u2J5YodZjXoDmq1pmGxhDCzwFbgF5gD3BvjHFk3Hk35F6/PsE97gX6gBeHEBpjjAMzuOZr486RJEmas+TSpDfv2EQIYc73bE0kQ3dadlVzVK1BwzrgH8Yd2xdC+KUY47cSxy7OvT41/gYxxuEQwj7gMuA84IkZXHM0hNALbAohNMUY++byTUiSJB061cd9z54EoCbAm7ZvTOW+BWVXeyy7qrmpxqDhk8C3gceAbrIf+G8Bfg34Wgjh+hjjI7lzV+ReJ+tRPnp8ZeLYTK5pzp03ZdAQQtg1yVuXTHWdJElaPJJlVl964WrWLl+Syn1t8KY0VV3QEGP84LhDjwJvCyH0AO8CbgXeNMPbjc79FVMDazbXSJIknSPGyBcfygcNN+5IZ5YBxpVdNadhQSlH+daqCxqm8HGyQcOPJo6NzhasOPd0AJaPO290uz13zckprjkz3YBijDsnOp6bgdgx3fWSJGlh233gNPs6egFoaazjVZetS+3erc3l7wodQiDGSCaToaZmodffmT+jQUMauS8ztZD+9o7nXpsTx/bmXi8af3IIoQ54ATAMPDvDa9bn7n/IfAZJkjRXd+7OJ0C/9or1LKmftO9s0ZIzDeVKhG5szAYuvb29ZXn+QjX65zn65zsfFlLQcH3uNRkA3JN7ffUE5/8o0AT8e6Jy0nTXvGbcOZIkSbPSPzTCPz9yZGz/xp1z782QlMxpKFfJ1ZaWFgCOHTtGd3c3mUym6jsjl8vojE13dzfHjmV7eoz++c6HqlqeFEK4DDgaY+wcd3wr8Je53c8k3voC8KfAzSGE2xLN3ZYAf5g752/GPeaTwO8At4QQPplo7rYKeG/unI+n8x1JkqTF6u4nnudM/zAAW1qbuGbbqlTvn6ye1Fmm5Umtra309vbS19fHoUOHpr9AM9bU1ERra+u8Pa+qggbgLcDvhhC+AewjWz3pfOAngSXAV8l2cwYgxngmhPCrZIOHb4YQbifb6fn1ZEurfgG4I/mAGOO+EMJvAx8DHgwh3AEMkm0Utwn4iN2gJUnSXBX2ZtiY+vr0tnF9GmKM87oGHqCmpobNmzfT2dlJd3c3AwMDzjTMQQiBxsZGWlpaaG1tndc8kWoLGr5B9sP+drLLkZqB08B3yPZt+Ic47icxxvjlEMLLgd8DbiQbXDwD/BbwsfHn5665LYSwH3g38Atkl3E9Drwvxvjp0nxrkiRpsTje3c+9T3eM7b95e7pLkwCW1NfS3FBL7+AIQyORM/3DrFhan/pzplNTU0N7ezvt7e3z/mylp6qChlzjtm9Ne+K5130XeG2R19wF3FXssyRJkqbzTw8fYSST/b3ltdta2dLWVJLntC1rpLczW7vlZM9AWYIGLQwLKRFakiSpKnwhsTTpxp3p9WYYr7W5/BWUtDAYNEiSJM2jx4508eSxbgCW1Nfw2ivWl+xZ7cvKX0FJC4NBgyRJ0jy6c1e+A/SrLltHy5LSLRlKzjSUq8GbFgaDBkmSpHkyNJLhnx7JBw037kg/ATqpoOyqMw2aA4MGSZKkeXLvUyfGlgmtXd7ISy4obUWhtoKZBoMGzZ5BgyRJ0jy5c3c+AfqN2zdSW1PavgltywwalA6DBkmSpHlwum+Qux8/PrZ/U4mXJgG0NueXJ53sMadBs2fQIEmSNA/u2nOUwZEMAFduWsGFa1tK/szxXaGl2TJokCRJmgdfTCxNKnUC9Kj2RCK0JVc1FwYNkiRJJfbDEz08dOA0APW1gZ964YZ5ee6q5nw511N9g2RyXailYhk0SJIklVhyluHHLl5T0D+hlBrramlZUgfASCbSdXZoXp6rhcegQZIkqYQymciXdid6M+ycn6VJoyy7qjQYNEiSJJXQfc+e5EhXPwCrmur5sYvXzOvzkw3erKCk2TJokCRJKqFkb4Y3XLWRhrr5/fjV6kyDUmDQIEmSVCK9A8N8/dFjY/vzVTUpqd0Gb0qBQYMkSVKJfO3RY/QNjgBw4ZplXL5x+byPoWCmweVJmiWDBkmSpBK5c1eiN8POTYQQ5n0MbYmu0DZ402wZNEiSJJXAoVN93PfsSQBqArxp+8ayjKMtuTzJBm+aJYMGSZKkEkiWWX3phatZu3xJWcaRnGk42evyJM2OQYMkSVLKYox88aFEb4Yd5ZllAGcalA6DBkmSpJTtPnCafR29ACxrrOMnfmRd2caSbO5mToNmy6BBkiQpZcneDD95xXqWNtSWbSyrkkFD3yAjmVi2sah6GTRIkiSlqH9ohH9+5MjY/o075783Q1J9bQ0rltYDECOc6nO2QcUzaJAkSUrR3U88z5n+YQC2tDZxzbZVZR5RYV6DS5Q0GwYNkiRJKfpiomrSm3dsLEtvhvGSeQ0dNnjTLBg0SJIkpeRE9wDfeurE2P6bt5d3adIoG7xprgwaJEmSUvKVhw+PJRpfu62VLW1NZR5RlmVXNVcGDZIkSSn5wq581aQbd5avN8N4yeVJJ51p0CwYNEiSJKXgsSNdPHmsG4Al9TW89or1ZR5RXtuyRFdocxo0CwYNkiRJKUgmQL/qsnW0LKkv42gKtdrgTXNk0CBJkjRHQyMZvvJwsmpSZSRAjzKnQXNl0CBJkjRH9z51go7ch/G1yxt56QXtZR5RoWT1pJO9Lk9S8QwaJEmS5ujO3fkE6Ddu30htTfl7MyQVzDS4PEmzYNAgSZI0B2f6h7j78eNj+zdV2NIkgFVNDYz2mDvdN8TQSKa8A1LVMWiQJEmag93PnWIw9yH8knUtXLi2pcwjOldtTWBVU3624VSfsw0qjkGDJEnSHOx+7tTY9nXntZVxJFNLVlAyGVrFMmiQJEmag10H8kHDjq2ryjiSqbVZdlVzYNAgSZI0SyOZyMMHTo/t76zgoKE90eCtwwZvKpJBgyRJ0iztPdZN7+AIkC21umHFkjKPaHI2eNNcGDRIkiTNUnJp0s6tqwihskqtJtngTXNh0CBJkjRLySToHVsqd2kSFOY02KtBxTJokCRJmqVdzxXONFSytkROw0lzGlQkgwZJkqRZONE9wIHOPgAa6mq4bMOKMo9oaq3ONGgODBokSZJmYXcin+HKjStoqKvsj1Xty0yE1uxV9k+3JElShdpdRUuTAFqbLbmq2TNokCRJmoVkPkMlN3UbtXJpPTW54k7d/cMMDmfKOyBVFYMGSZKkIg0OZ9hzuGtsv9IrJwHU1ISC2QaXKKkYBg2SJElFeuxI19hv6re0NrG6pXGaKypDYdlVlyhp5gwaJEmSilRNpVaTbPCm2TJokCRJKlKyclI15DOMSpZddXmSimHQIEmSVIQYY+FMQxXkM4xqX2YFJc2OQYMkSVIRDp8+y/Nnsh+4mxtquXhdS5lHNHM2eNNsGTRIkiQVYfeB02PbV21ZSe1oHdMqkMxp6DSnQUUwaJAkSSrC7ipdmgRWT9LsGTRIkiQVodqauiW1JXIaXJ6kYhg0SJIkzVDf4DCPHz0ztr+9mmcaXJ6kIhg0SJIkzdCeQ12MZCIAF65Zxoql9WUeUXHa7AitWTJokCRJmqFqbeo2avnSOupyids9A8P0D42UeUSqFgYNkiRJM7S7ivMZAEIINnjTrFR90BBC+PkQQsx9/cok57wuhPDNEEJXCKEnhPD9EMIvTnPfXwwh3J87vyt3/etK811IkqRKF2Nk14HqnmmAccnQ5jVohqo6aAghbAZuA3qmOOcW4C7gcuAzwCeADcCnQggfnuSaDwOfAtbnzv8McAVwV+5+kiRpkXm2o5fTfUMArGyq57z25jKPaHYsu6rZqNqgIYQQgE8CJ4GPT3LONuDDQCdwdYzx7THGdwJXAj8E3hVCuH7cNS8G3pV7/8oY4ztjjG8Hdubu8+HcfSVJ0iJSsDRpyyqyH0WqT7LBmzMNmqmqDRqAdwA3AL8E9E5yzi8DjcBfxhj3jx6MMZ4C/ji3+7Zx14zu/1HuvNFr9gN/lbvfL81x7JIkqcrsXgBLk4CCnAZnGjRTVRk0hBAuBT4E/I8Y471TnHpD7vXrE7z3tXHnzOUaSZK0wO0aN9NQrdpt8KZZqCv3AIoVQqgD/gE4ALx3mtMvzr0+Nf6NGOPREEIvsCmE0BRj7AshNAMbgZ4Y49EJ7vd07vWiGY511yRvXTKT6yVJUmXoOjvEU89nUyhrawIv3LyizCOaPRu8aTaqLmgAPgBsB14aYzw7zbmj/6K7Jnm/C2jOndc3w/MBVs5sqJIkaSF4+ODpse1L17fQ1FCNH6GyLLmq2aiqn/gQwrVkZxc+EmO8L41b5l5jkdfN6PwY484JH5qdgdhR5DMlSVKZFDR1q+KlSTC+5Ko5DZqZqslpSCxLegp4/wwvG50ZmGwOcXnu9cwMz59uJkKSJC1A1d7ULamw5KozDZqZqgkagGVkcwkuBfoTDd0i8Pu5cz6RO/bR3P7e3Os5OQghhPVklyYdijH2AcQYe4HDwLLc++NdmHs9J0dCkiQtTCOZyEMLpHISWHJVs1NNy5MGgP9vkvd2kM1z+A7ZQGF06dI9wEuAVyeOjXpN4pyke4Cfz13zyRleI0mSFqi9x7rpHRwBYO3yRjauXFrmEc3NssY6GmprGBzJcHZohL7B4arO0dD8qJqfkFzS869M9F4I4VayQcOnY4x/n3jrk8DvALeEED452qshhLCKfOWl8Y3hPk42aPi9EMKXR3s15Bq6vZ1s8DI+mJAkSQtUsj9DNTd1GxVCoG1ZA0e7+oHsbENTa9V8JFSZVNPypKLFGPcBvw20Ag+GEP4qhPAXwB7gfCZIqI4x/jvw57n394QQ/iKE8FfAg7n7vDvZKE6SJC1syXyGal+aNKrVvAYVacGHlTHG20II+4F3A79ANlB6HHhfjPHTk1zzrhDCHuAW4NeADLAb+O8xxn+el4FLkqSKsOvAwkmCHpWsoNRpV2jNwIIIGmKMtwK3TvH+XcBdRd7z08CEQYUkSVocTnQP8NzJPgAa6mq4bMPyaa6oDskKSh0mQ2sGFvTyJEmSpLlI5jNcuXEFjXW1ZRxNetps8KYiGTRIkiRNYvcCXJoENnhT8QwaJEmSJlHQ1K3KO0En2eBNxTJokCRJmsDgcIZHDnWN7e/YurKMo0mXDd5ULIMGSZKkCTx2pIvB4QwAW1qbWNOypMwjSk+rOQ0qkkGDJEnSBHYfOD22vWPLwpllAGg3p0FFMmiQJEmawEJs6jZqfHO3GGMZR6NqYNAgSZI0ToyRB5/rHNtfSJWTAJoaallSn/0YODCcoXdwpMwjUqUzaJAkSRrnSFc/z5/JLttpbqjl4rUtZR5RukIItDW7REkzZ9AgSZI0zq7E0qSrtqykrnbhfWQqqKBkMrSmsfD+BUiSJM3RQu3PkFTQq8Gyq5qGQYMkSdI4C7UTdFJrYnlSZ6/LkzQ1gwZJkqSEvsFhHjtyZmx/x+aFGTS0J5YndTjToGkYNEiSJCXsOdTFSCZbgvTCNctY0VRf5hGVhg3eVAyDBkmSpISCpUkLNJ8BoM0GbyqCQYMkSVLCQm7qltTWbPUkzZxBgyRJUk6MsaDc6kJNgoZxJVfNadA0DBokSZJy9nX0cqpvCICVTfWc195c5hGVTmvBTIPLkzQ1gwZJkqScXeP6M9TUhDKOprTaCkquDhJjLONoVOkMGiRJknJ2Hzg9tr1jy8oyjqT0ljbU0txQC8DQSORM/3CZR6RKZtAgSZKUs3uR5DOMal1m2VXNjEGDJEkS0HV2iKeOdwNQWxN44aaFPdMAhUuULLuqqRg0SJIkAQ8fPM3osv5L17fQ3FhX3gHNA8uuaqYMGiRJkhi3NGkBN3VLsuyqZsqgQZIkicJO0Au5qVtSa0EFJZcnaXIGDZIkadEbyUQeKqictDiChvbETEOHMw2agkGDJEla9J56vpuegWzJ0TUtjWxatbTMI5ofyQZvVk/SVAwaJEnSopds6rZz6ypCWLhN3ZLaliWqJ7k8SVMwaJAkSYteMp9hsSxNgnHVk1yepCkYNEiSpEVvsTV1G1VQPcnlSZqCQYMkSVrUOnoG2H+yD4CG2hou37i8zCOaP+NzGjKZWMbRqJKlGjSEEN4UQqhN856SJEmllJxluGLTChrrFs9Hmca6WlpyTexGMpEz/UNlHpEqVdozDXcCz4UQ/iCEsCXle0uSJKVud6LU6mLpz5DUZtlVzUDaQcNfA03A+4AfhhDuCiG8LiyWEgSSJKnqFHaCXlnGkZSHZVc1E6kGDTHGW4ANwC8DDwI/CXyF7OzDB0IIG9N8niRJ0lwMDmd45NDia+qWVFB2tceyq5pY6onQMcb+GOOnYozXA1cCfwMsA24F9oUQvhRCeHXaz5UkSSrW40fPMDCcAWBz61LWLF9S5hHNv4Kyq840aBIlrZ4UY3w0MfvwS8DzwOuBfwkh7AshvDuE0FzKMUiSJE2moKnbIpxlgHFlV81p0CRKXnI1FxT8AvAOYCMQgEeANuDPgCdDCFeVehySJEnjFTR1W4RJ0ACtzXaF1vRKFjSEELaHED4OHAE+DlwE/D2wI8a4g+zsw+8C7cDHSjUOSZKkyRQmQS/OoKHdBm+agbo0bxZCaAJ+Bvh/gJ1kZxWeIBs0fDrGeGb03BhjD/BnIYTNwFvTHIckSdJ0jpw+y9GufgCaGmq5ZF1LmUdUHm3NJkJreqkGDWRnFVqAEbI9G/46xvjNaa45DCy+rCNJklRWyXyGqzavpK625Ku2K5IlVzUTaQcN3cBHgE/EGI/N8Jq/Bj6b8jgkSZKmlMxnWIxN3Ua1mwitGUg7aNgaY8wUc0FuydKZaU+UJElKkfkMWasSMw2n+gYZyURqa+zLq0Jpz8PdHUL4halOCCH8XAjhnpSfK0mSNGNnB0d47Ej+d5bbF2En6FH1tTWsWFoPQCbC6T5nG3SutIOGVwDbpjlnK/DylJ8rSZI0Y3sOnWY4EwG4YM0yVjY1THPFwtZmXoOmUY6Mn6XAcBmeK0mSBMCuAzZ1S0o2eOswr0ETSDunASBOdDCEEIAtwGuBgyV4riRJ0ozsfu702PaOrYt3adKoZAUlG7xpInOeaQghZEIIIyGEkdyhW0f3k19kZxeeBa4Cbp/rcyVJkmYjxmjlpHHaluV7Nbg8SRNJY6bhXvKzCz8KHAD2T3DeCHAS+DeynaElSZLm3f6TfWMfjFcsree89mVlHlH5tTe7PElTm3PQEGN8xeh2CCEDfDLG+Adzva8kSVIp7CootbqSGsuLjmvw5vIknSvtnIYXAKenPUuSJKlMkkGDS5OyksuTbPCmiaQaNMQYn0vzfpIkSWl76IBN3cZrK0iENmjQueYUNIQQPkA2n+GvYoyduf3+vfSQAAAgAElEQVSZiDHG/3cuz5YkSSrWmf4h9j7fDUBNgBdutnISjJ9pcHmSzjXXmYZbyQYNdwCduf2ZiIBBgyRJmlcPHzhNzJVvuXT9cpobS1F9vvq02txN05jrv5Qfy70eGLcvSZJUccxnmNiqpnpCgBjhVN8QwyMZ6mrL0QNYlWpOQUOM8VtT7UuSJFWS3eYzTKiutoaVS+s51TcEZAOH1S2N01ylxcQQUpIkLQojmcjDB/JFHp1pKFSQ12DZVY1TsoV8IYRmYCVQO9H7McYDEx2XJEkqhaePd9M9MAzA6pZGNq1aWuYRVZa25gaeyW1bdlXjpT7TEEL4+RDCo8AZsrkO+yb4enYO9//TEMK/hRAOhhDOhhA6QwgPhRB+P4TQNsk1Lw4hfDV3bl8IYU8I4TdDCBMGNLlrXhdC+GYIoSuE0BNC+H4I4RdnO25JklReBfkMW1YRgk3dktqWWXZVk0t1piGE8J+B/wmMAN8GDgLDaT4DeCewG/hX4DjQDFxHtnLTr4UQrosxHkyM6Q3AnUA/+SpPPwX8BfAS4C0TfB+3ALcBJ4HPAIPATcCnQghXxBjfnfL3JEmSSswk6Km1NVt2VZNLe3nSu4FTwEtjjE+kfO9Ry2OM/eMPhhD+CHgv8B7gN3LHlgOfIBvEvCLG+GDu+PuBe4CbQgg3xxhvT9xnG/BhssHF1THG/bnjfwA8ALwrhHBnjPG+En1/kiSpBB5K5DPs2Gp/hvEsu6qppL086QLgCyUMGJgoYMj5XO71wsSxm4DVwO2jAUPiHu/L7f76uPv8MtAI/OVowJC75hTwx7ndt81q8JIkqSxO9gywr6MXgIbaGi7bsKLMI6o87YnlSR3mNGictIOGTrLLgMrhp3KvexLHbsi9fn2C8+8F+oAXhxCSNcWmuuZr486RJElVYHdiluHyjctZUj9pWuOi1ZpYntRp9SSNk/bypH8GXhFCCDGO9lssjRDCu4FlwArgauClZAOGDyVOuzj3+tT462OMwyGEfcBlwHnAEzO45mgIoRfYFEJoijH2TTPGXZO8dclU10mSpHSZzzC9gkRoZxo0TtpBw3uA7wIfDyG8K8bYk/L9k94NrE3sfx34zzHGE4ljo3OPXZPcY/R4cmHjTK5pzp03ZdAgSZIqQ7Kpm0HDxNrMadAU0g4aPk/2g/SvAD8bQngaOD3BeTHG+Mq5PCjGuA4ghLAWeDHZGYaHQgivizHunuFtRmutFTMrMuNrYow7J7xBdgZiRxHPlCRJszQ0kuGRg4kkaDtBTyjZ3K3D6kkaJ+2g4RWJ7WbgqknOS23pUozxeeBLIYTdZJcU/S/g8tzbo7MFk2U7LR933uh2e+6ak1Ncc2a2Y5YkSfPn8SNnGBjOALBp1VLWLF9S5hFVppVL66kJkIlwpn+YweEMDXWpt/RSlUr1JyHGWDPDr9Szj2KMzwGPA5eFENpzh/fmXi8af34IoQ54Adk+Eslmc1Nds55sMHRounwGSZJUGcxnmJmamlBQdvVUn0uUlLfQwscNudeR3Os9uddXT3DujwJNwL/HGJNzcFNd85px50iSpAq3y3yGGUs2eHOJkpLmLWgIIawKITTP8R6XhBDWTXC8JtfcbQ3ZIGD0vw5fADqAm0MIVyfOXwL8YW73b8bd7pPAAHBLrtHb2PjJNo8D+Phcvg9JkjR/HkrMNJjPMDUbvGkyqeY0hBBeCbwK+JPRD+4hhDVkE6RfCgyHEP4qxvhbs3zEq4H/HkK4F/gh2ZyDtcDLyZZNPQb86ujJMcYzIYRfJRs8fDOEcDvZXhKvJ1ta9QvAHckHxBj3hRB+G/gY8GAI4Q5gkGyjuE3AR+wGLUlSdThy+ixHurItpJbW13LJupYyj6iyWXZVk0k7Efq/AJfHGH8ncezDwMuAp4EW4L+GEL4XY/zcRDeYxt3A3wEvAV5ItlRqL9kE6H8APhZj7ExeEGP8cgjh5cDvATcCS4BngN/KnX9OUnaM8bYQwn6yZV1/geyMzOPA+2KMn57FuCVJUhkkS61etXkldbULbWV2upJlV08606CEtIOGFwLfGt0JISwl+xv6f40xviqE0AL8AHgbUHTQEGN8FHj7LK77LvDaIq+5C7ir2GdJkqTKYRJ0cZJlV0+a06CEtMPtNcCRxP6LyP5m/1MAMcZusl2jLz7nSkmSpJTtPpDvz2DQMD1zGjSZtIOGAWBpYv9lZHsy3Js4dgZoTfm5kiRJBfqHRnjscL4V0/YtK8s4murQnshp6DCnQQlpBw37gBsS+zcCT8cYDyeObSZb0UiSJKlk9hzqYjiTTV08f3UzK5saprlCrYmSq529Lk9SXtpBw6eBK0II3w8hfBu4Avj/x52zg3wDNUmSpJIwn6F4BdWTXJ6khLQTof8GuA74aSCQTST+09E3QwjXApcCn035uZIkSQUMGorX3pxMhDZoUF6qQUOMcQj42RDC27K7sXvcKc8C24H9aT5XkiQpKcbIQwds6las5UvrqKsJDGciPQPD9A+NsKS+tvgbDZ2F+qXTn6eqUZJixTHGMxMEDMQYO2KMj8QYuya6TpIkKQ3PnewbW16zfEkd569eVuYRVYcQwtwrKP3Lu+GPN8D//r0UR6ZyK0nQEEJYHUJ4Wwjhf4QQ/n7c8Wtz/RskSZJKIrk0acfWVdTUhDKOprrMKWjoOQEPfAJiBr731zDYl/LoVC6pBw0hhLeSXX70V2Q7RP9S4u21wH3Az6b9XEmSpFG7EkuTdro0qSjtiQZvHcU2eDv4vfx2zEDHUymNSuWWatAQQvhx4O+Ap4A3kU2MHpPr6PwY8MY0nytJkpS02yToWZvTTMOB7xXuH38ihRGpEqRdPem/AUeBl8cYz4QQtk9wzh7g+pSfK0mSBEB3/xB7n8+mVtYEeOFmm7oVo6DsarEVlA5+v3D/+OMpjEiVIO3lSVcD/xxjPDPFOYeAdSk/V5IkCYCHD54mZnu6ccm65TQ3pv070oWtrXmWvRqGzsKRhwuPOdOwYKQdNDQAvdOcsxIYSfm5kiRJgP0Z5qptWbJXQxE5DYd3Q2ao8JhBw4KRdtCwH9g5zTkvwo7QkiSpRAwa5mbWMw0Hv3fusTOHoN9K+wtB2kHDV4CXhRDeMtGbIYRfAq4E7kz5uZIkSWQykYcPnB7bt6lb8QpyGooJGg58f+Ljx5+c44hUCdIOGv4MOAB8NoRwB7mE5xDCLbn9vwOeBm5L+bmSJEk8fbyH7oFhIFs6dHOrraGK1dY8i+VJmUzhTMOGRC0ck6EXhFQzg2KMp0IILwf+F5CcbfhY7vXbwM/GGKfLe5AkSSpa4dKklYRgU7ditS6bRcnVjr35ZUhN7XDJ6+DIQ9l98xoWhNTLCcQYDwCvCCFcSXamoQ3oAr4XY9yV9vMkSZJGmc8wdy2NdTTU1jA4kqFvcISzgyMsbaid+qJkf4bNL4K1l+X3nWlYEEpWgyzGuIdsTwZJkqR5sfuAQcNchRBobW7g2Jl+AE72DrCpoWnqi5L9GbZcB2suze8707AgpJ3TAEAIYWsI4eoQws4QwpZSPEOSJCmps3eQfR3ZFdD1tYHLNqwo84iqV9EN3g7cl9/ech2s2AL1zdn9vg7oOZHyCDXfUgsaQgjtIYQ/DyEcBZ4Fvg/cD+wLIRwJIfz3EEJrWs+TJElK2p1YmnT5xhUsqZ9mSY0m1dpcRF5D9/Nwan92u7YR1r8QampgzSX5c1yiVPVSCRpCCBcCDwL/FVhLtnnbceBEbnsd8FvAgyGE89J4piRJUtKu5NIkS63OSXuiwVvHdBWUklWTNu6Euty1LlFaUOYcNIQQaoB/BLYA3wL+A7Asxrg+xrgOaAF+ArgX2AZ8Zq7PlCRJGs8k6PS0FTPTkOzPsOVF+e01P5Lfdqah6qUx0/ATwNXA54BXxhjviTGO/XTFGAdijHcDNwBfAF4UQvjxFJ4rSZIEwNBIhj2HEk3dDBrmpLWYBm/JfIbN1+W3nWlYUNIIGm4EBoD/EmOMk52Ue+8WYAi4KYXnSpIkAfDE0TP0D2UA2LhyKWuXLynziKpbe/MMlycN9sGxRLHMzdfmtwtmGp6AyT8mqgqkETTsAL4bY5w2LT7GeBz4Tu4aSZKkVLg0KV0zToQ+vAsy2Q7crL4EmhI1b5athaW5v4vBbug6VIKRar6kETRsBh4r4vzHgK0pPFeSJAkwaEjbjEuuHhzX1C0phHNnG1S10ggalgOnpz0r7zTZ5GhJkqRU7DZoSFVbYnnSlDMNyU7QW6479/2CvAaToatZGkFDA9myqjOVyV0jSZI0Z0e7znKkK9u9eGl9LZes83eTc5WcaejoGWDCtNVMBg4+kN8fP9MAJkMvIGk1dzOzRZIklcXu5/ILHl64eQV1tan1rl20mhpqaazL/jkODGfoG5zg98MnnoCBrux28xponaAVl2VXF4y0/lXdGkIYmckX8IGUnilJkmQ+QwmEEAoavE2Y15AstbrlRdkchvFWJ7pCn9gLmWIWp6iSpBU0hCK/JEmSUlHQCdqgITUFydC9E5RdTTZ12zxBPgNkqym1rM9ujwxA574UR6j5NOegIcZYM4uv2jQGL0mSFrf+oREeP9I1tr99s0FDWpJlVyecaTg4TRL0qGRewwnzGqqVi/4kSVLV+sHhLoZGsqmV561uZlWztVbSkqygdM5Mw5kjcPpAdrtuKay7cvIbWXZ1QTBokCRJVasgn2GLswxpKlyeNG6mIVlqdeNOqJsiWEvmNZgMXbUMGiRJUtX67jMdY9vmM6SrbarlSQcT+QxbJii1muRMw4Jg0CBJkqrSwc4+vpMIGl5yQXsZR7PwJHMazmnwVtDU7fqpb7T64vz2yWdgeIKkalU8gwZJklSVPv/gQUZ7jr3swnY2tzaVd0ALTLLkakdP4oP+QA8c+0FuJ8Cma6a+UeMyWLk1u50ZzgYOqjoGDZIkqeqMZCKfe/DQ2P7N12wp42gWpklnGg4/CDHXb2HNpbB05fQ3c4lS1TNokCRJVedbTx3n2Jl+ILv2/sd/ZG2ZR7TwFCRCJ3MaCvozTJPPMCpZdtVk6Kpk0CBJkqrOZ+8/OLZ9485NNNT5kSZtyZKrnb2DxNG1YAeLyGcY5UxD1fNfmCRJqirHz/Rzz5PHx/Z/+prNZRzNwrW0oZamhmw/3sGRDN0Dw5AZgYMP5E+arnLSKGcaqp5BgyRJqiqf33WIkUz2t97XvqCV81cvK/OIFq5zlig9/xgMdmcPLFuXT3CeTvuFELIBCKf2w2BvugNVyRk0SJKkqpHJRG5/4MDY/s9c6yxDKbUWLFEaOLc/Qwgzu1FdI7RdkN8/8WRKI9R8MWiQJElV499/eJKDnWcBWL6kjtdcvr7MI1rY2hMVlDp6BovrzzBewRIl8xqqjUGDJEmqGp9NzDK8eccmltTXlnE0C985ZVeTQcNMKyeNMhm6qhk0SJKkqnCyZ4D/89ixsf2bXZpUcm2JBm8DHc/BmVxvjPomWHdFcTczGbqqGTRIkqSq8MXdhxkaySZAX7V5JZesW17mES18bYmZhmUnduXf2LgTauuLu5kzDVXNoEGSJFW8GGPB0iQToOdHsnrS6lMP598oNp8BoPUFUJubueg+Cn2dcxyd5pNBgyRJqngP7D/FsyeyZTqbG2p53ZUbyjyixSGZ07C1d0/+jZn2Z0iqqYXVF+f3raBUVQwaJElSxUuWWX39VRtobqwr42gWj/ZcTkMzZ9k8tC93NMCma2Z3w4IlSuY1VBODBkmSVNG6zg7x1R8cHdu/+ZotZRzN4jK6PGl7zTPUkskeXHs5LFkxuxtadrVqGTRIkqSK9pWHD9M/lP3Aeun65Vy5aZYfWFW00eVJV9fszR+czdKkUSZDVy2DBkmSVLFijHz2/oNj+z9z7WbCTLsQa84a62ppaazj6pAIGjZfN/sbji+7GuPs76V5ZdAgSZIq1p5DXTxx9AwAS+preMNVG8s8osVndXMt22ueyR+Yy0zDik3Q0JLdPnsKep6f2+A0bwwaJElSxUomQL/2ivWsWFpkbwDN2Y7GwzSHAQAGmtbDijmUuw3BJm9VyqBBkiRVpN6BYf7p4SNj+z9zrQnQ5bAjPDW23dm6PfvBfy7WXJLfPm7Z1Wph0CBJkirSXY8coXdwBIAL1izj6q2ryjyixemykfxswKGWK+d+Q8uuVqWqChpCCG0hhF8JIXwphPBMCOFsCKErhPCdEMJbQwgTfj8hhBeHEL4aQugMIfSFEPaEEH4zhFA7xbNeF0L4Zu7+PSGE74cQfrF0350kSUr67AP5BOibrzEBuixi5LyzPxjbfWbJ5XO/p2VXq1K1dUZ5C/A3wFHgG8ABYC3wZuDvgdeEEN4SYz4VP4TwBuBOoB+4A+gEfgr4C+AluXsWCCHcAtwGnAQ+AwwCNwGfCiFcEWN8d6m+QUmSBE8cPcMjB08D0FBbw5t3bCrziBaproO0DJ4AoCcu4ak4h3yGUcmZhhNPQiYDNVX1e+xFqdr+hp4CXg9sijH+pxjje2KMvwxcAhwEbiQbQAAQQlgOfAIYAV4RY3xrjPG3gauA+4CbQgg3Jx8QQtgGfJhscHF1jPHtMcZ3AlcCPwTeFUK4vrTfpiRJi9vt9+cToH/isrVj/QI0zw58f2zzocwFnOjLzP2ezauhqS27PdgDXQenPl8VoaqChhjjPTHGu2KMmXHHjwEfz+2+IvHWTcBq4PYY44OJ8/uB9+V2f33cY34ZaAT+Msa4P3HNKeCPc7tvm9t3IkmSJtM/NMKXHjo8tm8CdBkduG9sc1e8iJM9g3O/Zwg2eatCVRU0TGMo9zqcOHZD7vXrE5x/L9AHvDiE0DjDa7427hxJkpSyrz16lDP92f+db2lt4vrz2so8okXsYH6m4YHMxZzsHUjnvpZdrTrVltMwoRBCHfALud3kh/2Lc69PMU6McTiEsA+4DDgPeGIG1xwNIfQCm0IITTHGvmnGtWuSty6Z5LgkSYtesgP0T1+zmZoaE6DLor8Lnn8MgJEYeDhzAUt7U5hpAJOhq9BCmWn4EHA58NUY4/9OHF+Re+2a5LrR4ytncc2KSd6XJEmz9MMTPdy/rxOA2prAW3aaAF02Bx8AsrVlnohb6WUpnb2DZDJx6utmwuVJVafqZxpCCO8A3gU8Cfx8sZfnXov56Z/xNTHGnRPeIDsDsaOIZ0qStCjckSizesMla1izfEkZR7PIHfze2OaemuwiiUyE02eH5p6Yvjqx6KJjL4wMQ23Vfyxd0Kp6piGE8HbgfwCPAz8WY+wcd8p0swLLx51XzDVnihiqJEmaxuBwhjt3HRrb/5lrUyjvqdk7kA8anmnM92c42ZNCXsPSlbB8Y3Z7ZBA6n537PVVSVRs0hBB+E/hL4FGyAcOxCU7bm3u9aILr64AXkE2cfnaG16wHmoFD0+UzSJKk4vzr489zMrdmfv2KJbz8ojVlHtEiNjIEh/OpmclO0CdLktdgMnSlq8qgIYTw38g2Z3uYbMBwfJJT78m9vnqC934UaAL+PcaYDJmnuuY1486RJEkpuf2BfG+Gt1y9mVoToMvn2B4Yyv1+dMVmwoqNY2+lUnYVTIauMlUXNIQQ3k828XkX8MoYY8cUp38B6ABuDiFcnbjHEuAPc7t/M+6aTwIDwC25Rm+j16wC3pvb/TiSJCk1Bzv7+PbT2f+lhwD/8WoToMsq0dSNzS+ibVm+On1namVXk8nQzjRUuqrKOAkh/CLwB2Q7PH8beEcI5/wWYn+M8VMAMcYzIYRfJRs8fDOEcDvZTs+vJ1ta9QvAHcmLY4z7Qgi/DXwMeDCEcAcwSLZR3CbgIzHG+5AkSalJJkD/6IWr2bSqqYyjUTIJmi3X0XY6n/jc4UzDolRVQQPZHASAWuA3JznnW8CnRndijF8OIbwc+D3gRmAJ8AzwW8DHYoznVEGKMd4WQtgPvJts/4cassnW74sxfjqV70SSJAEwPJLh87vyQYMJ0GUWY+FMw5braBvKBw2daeU0tF9MtihlhM4fwlA/1Fstq1JVVdAQY7wVuHUW130XeG2R19wF3FXssyRJUnG+sfcEz5/JLnlpX9bIKy9dW+YRLXKn9kNPrr5M43JY8yO0Pv/82NupdYVuaILWF2QrJ8UMdDwF66+c/jqVRdXlNEiSpIXljkQC9E07N1Ff68eTsjqYmGXYdDXU1NLeXILlSWCTtyriv0pJklQ2x7r6uefJfBHEm69xaVLZHUjmM1wPQOuyEixPAsuuVhGDBkmSVDaff/AgmVx24fXntbGtvbm8A1LhTMPmFwHQ1pyvnpRKc7dRJkNXDYMGSZJUFplM5I4H8wnQN5sAXX5nT+V/4x9qs8uTgFVN9WOnnD47xPBIJp3nuTypahg0SJKksvjOMx0cOnUWgJVN9bzqsnVlHpE4+EB+e/2V0JCd+amrrRkLHGKEU31D6Tyv9XyoyQUkXQeg/0w691XqDBokSVJZJDtAv2n7RpbU15ZxNAIK+zNsvq7grdbmEuQ11DVA2wX5/RN707mvUmfQIEmS5l1HzwD/+ni+jOfPXLuljKPRmIL+DC8qeCvZFbp0eQ0mQ1cqgwZJkjTv7tx1iKGRbAb0ji0ruWhtS5lHJIYH4fCD+f1xMw3tiQpKJ1OtoJTIazjxZHr3VaoMGiRJ0ryKMXLHA8kEaGcZKsKxPTDcn91euRWWry94O7k8yZmGxcegQZIkzavv7+vk2Y5eAFoa63jdleunuULzoqA/w3XnvF1QdrVkvRqsoFSpDBokSdK8uv3+fAL066/aQFNDXRlHozEH7stvb37ROW+3lWp50qptULc0u93zPPSeTO/eSo1BgyRJmjen+wb56qPHxvZNgK4QMRY2dZtupiHN5Uk1tbD64vz+CWcbKpFBgyRJmjdffugwg8PZxmCXb1zO5RtXlHlEAqDzWeg9kd1uXAGrLz3nlJKUXB1lk7eKZ9AgSZLmRYyR25MJ0Nc4y1AxkrMMm6+FmnM/IhZUT+pJO2gwGbrSGTRIkqR58fDB0zx5rBuApfW1vOGqDWUekcYk8xm2nJvPAOOqJznTsOgYNEiSpHlx+/35WYbXXbmeliX1ZRyNCiSbum0+N58BYGVTAzUhu911doihkUx6zx8/0xBjevdWKgwaJElSyfUMDHPXniNj+zdfu7mMo1GBvk7o2JvdrqmDjTsnPK22JhTMNpxKc7Zh+YZsLgVAfxd0H03v3kqFQYMkSSq5f3r4CH2DIwBcuGYZO7asKvOINObg/fnt9S+EhqZJT00GDR1p5jWEYF5DhTNokCRJJXf7A/neDDdfu4UQQhlHowIF/RkmXpo0Kll2Nf0KSjZ5q2QGDZIkqaQeO9LFnkNdADTU1vDm7RvLPCIVmKY/Q1JrQYO3FHs1gMnQFc6gQZIklVQyAfrVl69jVWKJi8pseAAO787vTxM0tJdqeRK4PKnCGTRIkqSSOTs4wpcfPjy2bwJ0hTn6CIzkZgxWvQCWrZny9NaC5UlpzzQkg4YnIZNidSbNmUGDJEkqmX/5wVG6+4cB2NbWxPXntZV5RCpQ0J9h6lkGgLZSNnhrbofmXNAyfBZO70/3/poTgwZJklQyt9+fT4D+6WtMgK44B2aezwDQVsoGb2AydAUzaJAkSSXxzPFuHnzuFAB1NYGbdm4q84hUIMbCJOhpKicBtC3LL0862ZPy8iQYlwxtXkMlMWiQJEklkUyA/g+XrmV1S+MUZ2venfwh9HVkt5eshPaLpr0kuTwp9ZKr4ExDBTNokCRJqRsYHuHO3YfG9k2ArkAF/RleBDXTfywsWJ6Udk4DWHa1ghk0SJKk1P2fx57nVN8QABtXLuVlF64u84h0joPfy2/PIJ8BYPmSeupqsnkp3QPDDAyPpDum1RfntzueguESBCaaFYMGSZKUumQH6P949WZqa0yArjhFJkED1NSEgj4bqS9RWrIcVuRmpTLD0PnDdO+vWTNokCRJqXruZC/ffeYkADUB3nK1CdAVp7cDTj6d3a6phw3bZ3xp6Zco2eStEhk0SJKkVN3xQD4B+uUXrWbDyqVlHI0mlKyatGE71M/876igV4PJ0IuGQYMkSUrN0EiGz+9KJkBvKeNoNKkDyXyGFxV1aVvzfJZdNWioFAYNkiQpNfc8eZwT3dkPkqtbGrnhkjVlHpEmVGR/hqTWUuY0gDMNFcqgQZIkpSbZAfotOzdRX+tHjYoz1A9HHsrvby5upqE9sTypoxQ5De0XQcj93HQ+C0Nn03+Giua/ZEmSlIojp8/yradOjO3/9DX2ZqhIRx6CkdyH/bYLYFlx5XCTXaE7e0uwPKl+KbSel9uJcGJv+s9Q0QwaJElSKj7/4CEyMbv9kgva2NrWXN4BaWLJ/gxFLk2CwuVJJameBC5RqkAGDZIkac5GMpHPPZivmnTzNSZAV6yC/gzFLU2CwuVJJameBOOSoS27WgkMGiRJ0px9++kTHD6dXXu+qqmen7hsbZlHpAllMinMNCSqJ5VieRI401CBDBokSdKc3X5/fpbhxh2baKyrLeNoNKmTT8PZU9ntpa3QfmHRtyjo01Cy5UmWXa00Bg2SJGlOTnQPcPcTz4/t33ytCdAVq6A/w3UQQtG3aGmso742e13f4AhnB0fSGl1e63lQmwtOzhyC/q70n6GiGDRIkqQ5+cKuQwznMqCv2baKC9a0lHlEmlRBf4bi8xkAQgiFDd5KsUSptj5benXU8SfTf4aKYtAgSZJmLcbIHQ/kezOYAF3hDtyX395SfD7DqJI3eINxeQ0mQ5ebQYMkSZq1+549yf6TfQC0LKnjtVesL/OINKme49lmaQC1jbBh+6xvNT95DSZDVxKDBkmSNGvJBOg3XrWRpQ0mQFes5NKkDduhrnHyc6fRviy5PMmyq4uBQYMkSZqVU72DfP3RY2P7JkBXuIIk6NnlM4wqbPBm2dXFwKBBkiTNyhcfOszgSAaAKzet4LINK8o8Ik2pIEbGF1AAACAASURBVGi4fk63Si5PKllOw4otUJ/rKt7XAT0nSvMczYhBgyRJKlqMkdvvNwG6agydhaOP5PdnWTlpVFtipqGjVDkNNTWw5pL8vkuUysqgQZIkFe2fHjnC08d7AGhqqOX1V20o84g0pcO7ITOU3W6/CJpa53S7ZMnVzlJ1hQaXKFUQgwZJklSUZ453854v/mBs/807NrKssa6MI9K0kqVW5zjLANCarJ5UquVJYDJ0BTFokCRJM9Y3OMxv/ONu+nJdgLe1NfHfXn3JNFep7JKVk+aYzwDQnmzuVqrlSQCrk8uTnGkoJ4MGSZI0IzFG3vflR3nq+eyypMa6Gv76P+2kZUl9mUemKWUy44KG2Td1G1U40zBAjHHO95xQwUzDE1Cq52haBg2SJGlG7njgIF/cfXhs/w/ecBk/smF5GUekGenYC/1d2e2mdmg9b863bG6opbEu+zGyfygzNvOUupZ1sGRldnuwG7oOleY5mpZBgyRJmtZjR7r4wD89NrZ/445N/Mer7ctQFZL5DFuugxDmfMsQQkEFpZKVXQ3h3NkGlYVBgyRJmtKZ/iF+4x93Mzic7clw8doW/vCNlxNS+PCpeXAg3aVJo9oSXaE7StXgDcZVUDIZulwMGiRJ0qRijPzO5/fw3Mk+ILss5a9/bgdLG2rLPDLN2MFEU7fNaQYN8zDTAJZdrRAGDZIkaVL/87v7+fpjx8b2P3TjlZy/elkZR6SidB+DU/v/b3v3HV9VfT9+/PW52Ttkskcg7CVDhrIEcVXUVq1fR2ut1bbuau342v60u+5qtWpb19e9AcUNiIqogCwJG0IICdl73/v5/fE5SW6Sm5CEu/N+Ph73cfY5n3w43Hve57PMfGgkDJjitlMnOVVP8mgPStLtql+QoEEIIYQQLm0+XMpfV7W+2f3BnGGcO0UGcQsoh51KGQZOg9DwzvftoRSn6kmeHavBqaShaA84PNToWnRJRmIRQrQorqpnze5CwkIU4wfEMyIlhtAQebcgRF9UWt3A9c9vpslhuricPDiB/z1n3HGOEn7HzV2tOmtb0uDBNg3RSRDbH6ryoanOlJwkj/Tc9YRLARc0KKUuBBYAU4EpQBzwvNb68i6OmQvcAcwGIoF9wJPAw1prl+GqUuo7wG3ASUAI8C3wqNb6Gff9NUL4h6Kqev697gD/tyG7Tbd54aE2RqfHMq5/PGMHxDNuQBzjB8STGO2+N1VCCP/jcGhueWULR8vrAIiPDOWRS6cRESrtGAKOc0mDm4MG596TPFrSAKa0ocqqJlewU4IGHwi4oAHz8D8FqAKOAF0OQ6mUOg94HagDXgZKgHOBB4BTgItcHHM98DBQDDwHNAAXAk8rpSZprW9z1x8jhC8VVNbxxCcHeO7LbOoaHR22NzQ52JFbwY7cijbr+8dHMm5AHOMGmGBi/IA4hidLqYQQweLRtftYu7uwZfn+i6cyJCnahykSvdJQDfnbWpcHz3Tr6ZNjvRk0jIcDa8x8QRaMO9ez1xMdBGLQcAsmWNiHKXFY09mOSql44N+AHViotd5orf8dsBq4UCl1idb6JadjhgP3YoKLGVrrQ9b6PwBfA7cqpV7XWjt1eixEYCmoqOOxTw7w/JfZ1De1DRbG9o9jYGIUWXkV5FlvGdvLr6gjv6KONU4PFRGhNkanx7UGE/3jGT8gnoRoGSlWiECyfn8R93+4p2X52gUZLBmf7sMUiV7L3QSOJjOfOs5U83Gj5BinNg2erJ4E0u2qHwi4oEFr3RIkdKN/6AuBVODZ5oDBOkedUuoO4GPgZ8BLTsdcBUQAf28OGKxjSpVSfwH+C/wUkKBBBJz88joe+2Q/L3x1uKW/9WYTBsZz4+JMTh+Xjs1m/m+V1TSQlVdJVl4Fu/IryMqrZPexyg7HAtQ3OdieW8723PI26wcmRLZUbWoOJkakxBBik/7dhfA3BRV13PjiFqxmDJw8PIlfLh3j20SJ3mszPsMst58+yRuDuzWTAd58LuCChh46zZq+52LbOqAGmKuUitBa13fjmHfb7SNEQDhaVstjn+znpa9yaLC3feCfNCiBGxdnsmRcWodAPDE6nDkjk5kzMrllXZPdwaHianbmVbIrr4KsPBNM5Fe4LpU4Wl7H0fI6Vu8qaFkXGWZjTHocY/vHt6nmlBAlpRJC+EqT3cH1L37TMkhXSmw4D196klQ7DGQeGp+hWZvqSVUNaK09N+BfqlPwWrwPmuohNKLz/YXbBXvQ0HyH7Wm/QWvdpJQ6CEwAMoCsbhyTp5SqBgYrpaK11jVdXVwptamTTV22wxDCXXLLanl0zT5e3XikQ7AwZXACNy3JZNGYjsFCV0JDbIxKi2NUWhzLnLpeLK1uIMsqjdiVV0FWfgV7jlW5LJWoa3Sw9Ug5W4+0LZUYlBjFuAEmmJg2LLHHaRNC9N59H+7hq4MlACgF/7jkJNLjI32cKtFrDjvkfN267IGShujwUKLDQ6hpsNNgd1BV30RcpIde/kTEQuIwKMs2Va6K90H6BM9cS7gU7EFDgjUt72R78/rEHh4TY+3XZdAghK/klNTw6Nr9vLYph0a7brNt6pBEblqSycLRqW59IO8XE87ckSnMHZnSsq7J7uBgUTU7rdIIU8WpgmMVruu+5pbVkltWy0dZplTijAnpPHrZdKnKJISHrd51jH+t3d+yfMuS0ZwyKqWLI4TfK8iCeutRJjYd+o3wyGWSYsKpaagFTGmDx4IGMFWUyrLNfEGWBA1eFuxBw/E0P4noLvfq5TFa6+kuT2BKIKb14JpCdMvh4hoeWbOP1zcfaelbvdn0Yf24aXEm8zJTvPb2PjTERmZ6HJnpcZw3tXV9SXUDu/Iq2gQTe49VdSgNef/bY9y18lvuWjZBShyE8JAjpTXc8vLWluX5o1O5ftEoH6ZIuEWbqkmzTPGRByTHRnCk1AoaqhsYnhLjkesApjH0HqumuDSG9rpgDxqaSwsSOtke326/5vkU65jiLo6pcLFNCJ84VFTNI2v28cY3udjbBQszh/fjpsWjOWVUst88eCfFhDN3VApznd5kNtodHCisZld+BR9lFbBy61EAnv0imyH9ovnJ/AxfJVeIoNXQ5OC6F76hvLYRMN0pP3DxlJbOEEQAO+y5Qd2cJXtrgDeQxtA+FuxBw25gBjAaaNO+QCkVCowAmoAD7Y5JsY75ot0xAzBVk44crz2DEN5wsKiah1fvZfmWox2ChVkjkrhpSSZzMvwnWOhKWIiNMf3jGNM/jnMnm7YSzYHDn1dlMTAxinMmD/BlEoUIOn9ZlcXWnDIAQm2KRy47ieRYaVwaFA57thF0s2Sv9qAk3a76UrAHDauBy4AzgRfbbZsPRAPrnHpOaj7mFOuY9t2qnuW0jxA+s6+gikfW7GP5llzaxQrMyUjmpiWZzM5Idn1wALDZFPdeNJljFXUtDTNveWULafERzBzu3n7Gheir3tmWx9PrD7Us//qssUwfJv+/gkLFUSg/bOZDo2DAZI9dKsmbA7ylZIIKAW2H0kNm8LpwD1aHEm0Eez9qrwFFwCVKqRnNK5VSkcCfrMV/tTvmKaAeuN4a6K35mH7Ab63FxzyUXiG6tK+gkhtf/IbTH/iEN79pGzCcOiqFV66dw4vXzA7ogKFZRGgIT1wxnZGp5gehocnB1c9sZH9hlY9TJkTgO1BYxa9ebx0p+IwJ6fz4VM80lBU+4FzKMHgGhHiucXKK0wBvRZ6unhQaAclO7W0Kd3n2eqKNgCtpUEqdD5xvLfa3pnOUUk9b80Va69sAtNYVSqmfYIKHtUqplzAjPS/DdK36GvCy8/m11geVUr8EHgI2KqVeBhowA8UNBu6T0aCFt+3Or+Th1Xt5Z3seul3JwrzMFG5ekhmUbwgTo8N5+kcnc8Gj6ymqqqe8tpErn/qKN352CqlxUoVCiN6obbDz8+c3U1VvRgoemhTN3RdOCYhqjKKbcpzaMwxxf1erzrw6wBuYKkpFu818QRYMctnnjPCAgAsagKnAD9uty7A+ANnAbc0btNZvKaUWAP8LfA+IBPYBvwAe0rr9IxhorR9WSh2yzvMDTInMTuAOrfUzbv1rhOhCVl4FD6/ey6rt+R22LRyTyo2LM5k2tJ8PUuY9Q5KiefLKGXz/8Q3UNtrJKanl6me+5sVrZhMdHohfYUL41u+X72BXfiUA4aE2Hr1smgysGGwOO73b9GAjaOg4wJvHpY2DnW+ZeWkM7VUB94urtb4TuLOHx3wOnN3DY1YCK3tyjBDu8u3Rch76eC/vf3usw7bTxqZx4+JMpg5JdHFkcJo8OJF/XnoSP3l2Iw4NW4+Uc+OL3/D4FTNkDAcheuCVjTm8uulIy/Kd505g4qDOOhgUAam+CvJ3WAsKBs/06OVSnBrOe7xNA0hjaB8KuKBBiGC2I7ecf3y8lw93dgwWloxL56bFmUwa3Dd/4BePS+cP503kjrfMj+FHWQXcueJb/nCejOEgRHfsyq/g98t3tCxfcNIg/ufkIT5MkfCI3I2moTCYLkqjPPuCKcmbXa6CdLvqQxI0COEHsvIquO+D3S0jITtbOj6dGxdnyttA4PLZwzhSWstjn5iRa/9vQzZDkqK4Zv5IH6dMCP9WWdfIz5/bTF2jGUAxMy2WP18wUQLuYNRmfAbPtmeAjm0atNaeva/6jYCQCLDXQ2Ue1JRAdPC16fNHEjQI4UOHiqq5/8M9rNx2tEMD57Mm9ueG0zIZPzDe9cF91O1njOFoWS0rrDEc/rJqFwMSojh3ykAfp0wI/6S15tdvbOdAUTUAUWEhPHrZNGkTFKzatGeY4/HLRYaFEBsRSlV9E00OTUVtEwnRHmwjExIKqaMhf7tZLtwFw+Z67nqihXxjCOED+eV1/OPjvbyyMafNoGxKwdkTB3DD4lGM7S/Bgis2m+IeawyHL60xHG59ZSvp8ZGcPELeNgnR3rNfZPPOtryW5b9+dxKZ6XE+TJHwGIcdjmxsXfZwz0nNkmPDW3rjKq6u92zQAKaKUnPQULBTggYvCfZxGoTwKyXVDfz5nZ0suGcNL351uE3AsGRcOqtunMcjl02TgOE4zBgOMxiVFgtAg93BT57dyL4CGcNBCGdbcsr40zutjUUvnTWU808a5MMUCY869i00mJ6xiBsAiUO9ctk27Rq83hha2jV4i5Q0COEFVfVN/PfTg/z70wMtb2Oazc5I4pdnjGX6sODuOtXdEqLDeOrKmR3GcHjz5zKGgxAAZTUNXPf8Zhrt5uXExEHx/P47449zlAho2etb54fMMsXXXpDsNMCbNIYOXhI0COFBdY12ntuQzaNr93cY9Gby4AR+ecYYTh2VIo0Re2lIUjRPXTmTix//gtpGO0dKa/nxM1/zkozhIPo4h0Nz6ytbyS2rBSAuMpRHL51OZFiIj1MmPKa2DD5/sHXZi1V2kn1a0rATtPZagNSXSfUkITygye7gpa8Os+jetfzpnaw2AcOotFgeu3way687hXmZqRIwnKBJgxN45LKTaB6uYduRcm544Rua7A7fJkwIH3p83QE+3tXaG9u9F01haHK0D1MkPO7D35vehABi0mDyxV67tNcHeEsYAuGmeiq1pVDVsZty4X4SNAjhRg6HZuXWo5z+wDp+/cZ28srrWrYNSozi3oum8P7N8zlz4gAJFtzotLHp/PH8iS3LH+8q4M6V3+JiwHchgt6XB4q594PdLcs/mTeCMyb092GKhMcd+AQ2P9O6fM69EOW9Kq/JTgO8tS9V9wilpF2DD0j5vRBuoLVm7e5C7nl/NzvzKtpsS4mN4IbTRnHJyUOICJWqAZ5y2SwzhsO/1poxHJ7bcJjB/aL56QIZw0H0HYWV9dzw4jctnSxMH9aP288c6+NUCY9qqIGVN7YujzsXxp/n1SQ4V08q8kabBjBBw5GvzXxBFoxc5J3r9mESNAhxgr46WMI97+/i60OlbdbHR4Zy7YKR/OiU4VK/3kt+uXQMuaWtYzj87d1dDEyMYpmM4SD6ALtDc9NL31BQaR7akmLC+eelJxEWIpUKgtqaP0PpITMfmQBn3+v1JDhXT/JKSQO0awy9s/P9hNvIk4wQvbQjt5x7P9jN2t2FbdZHhYXwo1OGc+38kZ7vq1q04WoMh9te2Up6XASzMpJ9nDohPOvBj/awfn8xYGpvPPj9qQxIiPJxqoRHHdkEGx5tXT7jLxDn/apobbpc9UabBpDqST4gQYMQPXSgsIr7PtzTZrAkgLAQxaUnD+W600aRFhfpo9SJ5jEcvvfYevYVVLWM4fDGz+cyKk0GtBLBae3uAh5eva9l+YbTMpk/OtWHKRIe19QAy68DbXX6kLEQpl7mk6SkOLVp8ErvSdC2pKFwFzgcYJNSNU+S3BWim46W1fKr17Zx+gPr2gQMSsF3pw1i9a0Lueu8iRIw+IGE6DCe/tHMlvEaKuqauPKprymorDvOkUIEnqNltdzy8paW5VNHpXDT4kwfpkh4xWf3Q6H1hj0sGs79h8+6He0X7Vw9qR6HwwudUMSkQrRVgtxQBeU5nr9mHydBgxDHUVxVzx9W7mThPWt5eWNOm1Gcz5iQzvs3z+f+i6cyJEm6M/Qng/tF8+QPZxIdbhqfHymt5cdPb6Smoek4RwoROBqaHFz3wmZKaxoBSI+P4MFLphJik97ZgtqxnbDOqe3C4v8H/Yb7LDnhoTbiI03lFYeGstpGz19UKRnkzcskaBCiExV1jdz/wW7m372GJz8/SINTv/+njkph+XWn8PgVMxidLlVe/NWkwQk8cum0ljEctufKGA4iuPz9vV18c7gMgBCb4p+XTmtTVUQEIYcdVlwPDuvBfPDJcPJPfJsm2ne76sUelJpJY2iPk6BBiHbqGu08sW4/8+9ew0Or91HdYG/ZNnVIIi9cPYvnrp7FlCGJPkyl6K5FY9P40/mTWpY/3lXA/1shYziIwPfejjz++9nBluXbzxjDzOFJPkyR8IoN/4LcTWY+JByWPQw233fn3bbbVWkMHYykIbQQlka7g1c25vDQx3s5VtH2Lcno9FhuWzqG08eny6BsAejSWUM5UlrDo9YYDs9/acZw+NlCGcNBBKaDRdX88tVtLctLxqVxzfwMH6ZIeEXJAVj9p9bl+bdDmn+Mw+Hcg5Jvul2VoMHTJGgQfZ7doXl721Hu/3AP2cU1bbYNSYriF6ePZtmUQVJHOMDdtnQMuWW1LN9ixnD4+3u7GJgYyXlTB/k4ZUJ0n9aat7bkcueKnVTWm/Y5g/tFcd9FU+WFRrDTGlbcCE21Zjl9Ipx6s2/T5MS5elKxtwZ4S3UKmIp2g70JQuTR1lMkZ0WfpLXmm5wyVmw5yjvb8yisbPsFlxYXwQ2LM/n+jCGEh0otvmBgsynuvtCM4bDhgBnD4ZevbiM9PpLZMoaDCADHKur43ze381FWQcu68BAbj142TcaE6Qs2PwuHPjXzygbn/RNC/OffPcVpgDevdbsalQjxg6AiF+wNpiQmdbR3rt0HSdAg+pRd+RWs2HKUlduOklNS22F7QlQYP1s4kh/OGU5UuO/riAr3iggN4fHLZ3DhY+vZa43hcI2M4SD8nNaaNzbnctfKb6moa+39y5QwTGHyYGlfFfQqjsIHd7Quz70BBp7ku/S44JMB3sCUNlTkmvmCnRI0eJAEDSLoZRdXs3LrUVZsPcqeY1Uu90mJjeB/Th7C1fMySIjynzc3wv0SosN46kczueDR9RRW1lNR18QPn/yaN6+bK2NsCL+TX17Hb9/czupdBW3W/2DOMH515lhiIuRnPOhpDe/cCvUVZjkpAxb+xrdpcqFt70leDBrSxsH+j818QRZMON971+5j5NtGBKVjFXW8vS2PFVuPsjWnzOU+8ZGhnDVxAMumDmR2RrK0WehDBveL5qkrZ3Lx419Q02Ant6yWq57+mpevmSMPYcIvaK15bdMR/vD2TiqdSheGJEVx9/emMGekVKnrM759A3aval1e9jCERfkuPZ1o23uSl9o0QLvG0NLtqifJr6MIGmU1Dby7I58VW46y4WAxrnrUjAoLYcn4dJZNGcj80SlEhEoVpL5q4qAEHrlsGlc/sxG7Q7Mjt4LrX9jMv38wg9AQaccifCe/vI7fvLGNNbsL26z/4Zxh3C6lC31LdTGsur11ecZVMPxU36WnC8mxPug9CaTbVS+Sbx4R0Krrm/go6xgrthxl3d5CGu0dI4WwEMWC0amcO2UgS8alyw+uaLFoTBp/PG8iv31zOwBrdhfyu+Xf8pcLJkpPNMLrtNa8uukIf2xXujA0KZq7L5wsDfb7ovd/AzVFZj5+ECy5y7fp6UKbNg3eDBpSxwAK0FCyHxrrIEyqmnqCPD2JgFPfZOeT3YWs2HqUj7KOUdfYcXRfpWD2iGSWTR3IWRP7kxgd7uJMQpgxHHLLanhkjRnD4cWvDjMkKYqfLxzl45SJviSvvJZfv76dT/a0LV24cu5wbj9zDNHh8nPd5+z5ALa93Lr8nQcgMt536TmOJKff2dKaBuwO7Z1qv+Ex0G84lB4E7YCiPTBgsuev2wfJt5AICHaH5ov9xazYmst7O/Lb9CDibMqQRJZNGch3Jg8gPV7eNIjuuW3pGI6Uto7hcPd7uxmUGCVjOAiP01rz6kardKG+9XttWHI0d39vMrOkdKFvqquAt29pXZ50MYw+w3fp6YbQEBuJ0WGU1TSitQkcUpwaR3tU2ngTNICpoiRBg0dI0CD8ltaazYfLWLn1KG9vy+u0YdXo9FiWTRnIuVMGMiw5xsupFMFAKRnDQXjf0bJafv3GdtY5lS4oZUoXfnmGlC70aR/dCRVHzHx0Mpz5N58mp7uSY8Ipq2kETLer3gsaxsHud8y8NIb2GPlGEn5Fa82u/EpWbD3Kyq1HOVLacSwFMP2TL5sykGVTBzK2v/8W14rAEREawuNXzODCf7Udw+H1n80lM13GcBDuo7Xm5a9z+NM7WVQ5lS4MT47m7guncPKIJB+mTvjcoc9h439bl8+6G2IC4+VFcmwE+wurASiurge89N0pjaG9QoIG4Reyi6tZscWMpbC3oPOxFL4z2XSRetKQRGmoKtwuISqMp686mQse+ZyCljEcvuLxK2YwaXCCr5MngkBuWS2/fn0bn+4talmnFPxo7gh+ecYYGVSyr2ushRU3tC6PORsmfs936emhZF8N8Nam21UJGjxFggbhdXWNdg4VV3OwsJp9BVV8lHWMrUfKXe4rYykIbxuUGMWTTmM4HC2v44JHP+fWpWO4Zn6G3IOiV7TWvPR1Dn9uV7owIiWGey6czIzhUroggLV/Mz0AAUTEwzn3magyQPis29XkUWALBUcTlB82bUL8uNF4oJKgQXhEk91BblktB4pMcHCwqPWTW+a6ylEzGUtB+NrEQQk8dvl0fvbcJqob7DQ5NH9/bxdrdxfwwPenMjDR/wZWEv7rSGkNv359O5/ta1u68ONTRnDrUildEJaj38D6h1uXl/4R4gf6Lj29kBTT2oah2JsDvIWGQ3ImFFqlDIW7YchM712/j5CgQfSa1prCynoTGFifA4XVHCyq4nBJjcsxEzojYykIfzN/dCqrbprHTS9tYYs1qviXB0s488F1/OW7k/jO5MD6MRfep7Xmha8O85d3sqhusLesz0iJ4Z6LJjN9mJQuCIu9EZbfANq6T4bPg2k/9G2aeiEl1kdjNYBp19ASNGRJ0OAB8mTmIw1NDgoq6oiOCCU6LASbH1d5KK9ttIKCKg4WVnOgqLqlepHzD2F32RQMSYpmREoMI1JiGDcgnqXj02UsBeF3hiXH8OpP5/Dw6n38c/VeHBoq6pq4/oVvWLOrkDuXjScuMszXyRR+KKekhl+/sY3P9xW3rFMKfjIvg1+cPprIMCldEE4+fxCOmUEmCY2CZQ8FVLWkZkm+atMApl3Dt2+YeWnX4BESNPjInmOVnPyXj1uWo8NDiA4PJTbCTGMiQoiJCCUmPJTocGs+onkf53WhxIQ7rY8IISY8lMgwW48aCtc12skuruFgUVWHKkW9fVuQHh9hBQaxjEiJtqYxDE2KJjzU1qtzCuFtYSE2fnH6aOZnpnDzy1taevR6ffMRvj5UwgPfn8r0Yf18nErhLxwOU7rw11XtShdSY7jnwilyr4iOCnfDJ3e3Lp92ByRl+C49JyDZuXpStRerJ0G7HpSk21VPkKDBT9Q02KlpsFPkuuOgHrMpTMDRLvgwgYUJNGw2RU5JDQcKqzlaXovufm2iFnGRoWSkxpJhlRo0f4anxBArVYxEILI3QfFeyNtq3lbFDYCJ32PGcFNd6c7l3/LGN7kAHC6p4eLHv+CG00Zx/aJRhIZIMNyX5ZTU8KvXt7F+f2vpgs0qXbhFSheEKw47LL8e7NbLuUHTYfbPfJumE5Ds6+pJzaSkwSPkqc5HwkNtpMRGUF3fRG1jz6v4HI9DQ2V9kzXC6IlF++GhNkYkWwFBqpk2BwlJMeHS9akIXE315sclb2vr59i30NSusf4Hd8DYs4mf9kPuv2gRC8akcsdbO6isa8Lu0Dz40V7W7Snkwe+fxNDkaN/8LcJnHA7N819m89d3d1HjVLowMjWGey6awrShUrogOvH1f+DIV2beFgbL/gm2wA0ufdblKkC/4aZqV1MtVB2D6uKAGd8iUEjQ4COj0+PYeMcSAOwOTW2jner6Jqrrm6hpsFNV30RNQxPV9db6Bjs19U1UNTRRU2+nuqF1X3OcWdd8bEOTo0fpsSkY3K+1nUFGamupwcCEKL9ucyFEtzRUQ/4OyN8GeVuskoRd4Gg8/rGORti53HwShnLetCuY+ePvcfOqAr46aEaQ3ny4jLMf+pS7lk3gu9MGSTDdR+SU1HD7a9v44kC70oX5GdyyREoXRBdKs+Gju1qX590K6eM73z8AJEaHY1PmxWV5bSONdgdh3iqBtYVA6hjz/Q6mMXTMqd65dh8hQYMfCLEpYiNC3Vqdp9HukQglTwAAIABJREFUaAkoOgQfDU0tgUVzoDAkKUq6NhXBo7bMCg62Qp41Ld4LupvBdPwgGDDF/AAd+rz1TSCYPsDX/JmBa//Ky6NOZ9X0M/jFN2nUO2xU1Tdx66tbWbO7gD+fP4mEaGkkHawcDs1zX2bzt3alC6PSYrnnwsmcJKULoitaw8qboNGMnkzqOJj3C9+myQ1CbIp+0eEtVZNKqxtIi4/0XgLSxrcGDQVZMFyCBneSoCFIhYXYSIiykRAlDy0iyFUVWIGBVXqQvw1KD3X/+H4jTIAwYLKZ9p8Csalt9zm2EzY/C1tfhDrT/Sragdr7PufwPksT03i5cT6PV51Cjk7n7W15bM4u5b6LpzJnpBSPB5vs4mpuf20bX1qlTGBKF65dMJKbFmdK6YI4vi0vwIE1Zl7Z4Lx/QmhE18cEiOTY1qCh2OtBgzSG9iQJGoQQgUFrKD/iVIJgfSrzune8skHKaCtAmAL9J0P/SRCVePxj08fDWX+DJXfCrrdh09Nw6NOWzWE1BVzOa1we8Rqf2ifykv00PiyfzqX/2cBPF4zkliWjpcewAFdQWcf73x5j1bY8vjxYjMOp44jMtFjuuWgKU4d0414SovIYvP+b1uXZP4fBM3yXHjfzeberzaQxtNtJ0CCE8D9aQ8mBtsFB3laoLTn+sWAaFKaNaw0QBkyB9AkQHnNi6QqLhEkXmk/xflP6sOUFqC5o2WVeyA7mheygWMfxhn0eL32yiM/2FvHgJVMZmRp7YtcXXpVfXsd7O/JYtSOfrw+VdOhhzqbgpwtGcqOULoieWHUb1JWb+X7DYdFvfZocd0uO9aNuV7UOyPEu/JUEDUII/1BfCfvXwN73Ye+HpveL7giNgv4T2wYIqeMg1MODBSaPhNPvMn2q73kPNj0D+z4CzJNlsqrkJ6Gr+EnoKr4qHMO/H17MtLOu5KLZo6WRtB/LLavl3e15vLsjn03ZpS73UQpOHp7Eb88exxQpXRA9sXM5ZK1oXT73oRN/meFnfNqDUvxAiIiH+goTmFXmmXXCLSRoEEL4TtE+EyTseR+y1x+/J6OIeFOtyDlASB4FIT78KgsJg3Hnmk9ZDnzznPlUHGnZ5WTbbk5mNxXvPc3aL5Yw7fybSciY7rs0izYOF9fwrlWisDWnzOU+zYHCOZMHcMaE/qR7s562CA41JfDOba3L034AGQt8lx4PcR7grcTbYzUoZUobcr40ywU7JWhwIwkahBDe01QP2Z/Dng9MsFByoPN9IxNh0LS2AULicLD5cduAxCGw6Dew4HbYvxo2PY3e/R5KNwEQr2pYVLECnl1BZdIk4uZeBRMvhMh4Hye87zlYVM2q7Xm8uyOPHbkVLvexKZgzMpmzJppAITUuOBqqCh/54I7WqoxxA+D0P/o2PR6S1GaANy9XT4J2QUMWjFri/TQEKQkahBCeVZEHez8wnwNroaGLYc/7T4LMM2D0GWZk1EAd5MgWApmnQ+bpqMpjNG5+nor1/yW5vrX0Ia5kO7x9C/r9O1ATL4BpV5rGkFJ1yWP2FVTx7nZTopCV5zpQCLUp5o5K4eyJ/Tl9fHqb+tlC9Nq+j2HL863L59zfvU4YAlCKU/WkIm9XTwJpDO1BEjQIIdzLYYfcza3VjvK3db5vWAxkLITRSyFzaXAWI8elE7bgFyTPu5lvPn2bY2sfZ5FjAxHKlD6oxurWKk1p42HaD2HyxRCd5OOEBz6tNXuOVbWUKOw55jpgDQtRzMtM5SwrUEiM9nB7GNG31FfByptblyd8F8ae7bv0eJhzoO316kkg3a56kAQNQogTV1sG+z821Y72fQg1xZ3v22+EKUnIXGoG3gmSvsmPy2bjpAXLKJp+Bre9/CkpB97ikpA1jLG1lj5QsBPe+xV8+HsYv8wEEMNPldKHHtBaszOvgne357NqRx4HCqtd7hceamPB6FTOntSf08amy5g2wnNW/9EMCgkQ1Q/Outu36fGwtl2u+qJ6knNJwy5wOPy7WmsAkaBBCNFzWkPhLlOSsPcDOLwBtN31vrZQGDa3tdpR8qg+/RCcEhvBQ1ct5rkvM1n29tlMaNzNJSFr+E7IBqKV9QNrr4ftr5pPUoZpMDnpYkgY5NvE+ymtNdtzy1m1PZ93d+SRXVzjcr+IUBuLxqRx1qT+LB6XTmyE/AQKDzu8Ab58vHX5zL93HDwyyKQ4tWk4XFLDba9u5dr5GWSmx3knATEpEJMK1YXQVAtlh8z3qDhh8o0phOiexlo4+KnVJeoHUHa4831j0kxJwuilkLFIGvq2o5TiitnDmD0iiZteiuP2vNH8oekKloV8weVhaxiPUwPxkgPw0Z3mM2Q2TDgfxp8XnFW5ekBrzZacMt7dkc+q7XkcKa11uV9UWAinjUvj7IkDWDgmlRgJFIS3NNbBihto7oaZUaebqodBLj4yjLS4CAoq63FoeG3TEV7bdITFY9O4dsFIZg7v5/lup9PGwcFCM1+QJUGDmyjdfrQa4XFKqU3Tpk2btmnTJl8nRYiuleWYAGHP+3BwnXlr05mBJ7WWJgyYKsXB3VTfZOe+D/bwxLrWQGGCOsTvBn7NrKqPUPWVrg8cOgcmXADjlkH8AC+l1rdqG+zsOFrOu9vzeW9HHkfL61zuFxsRyuJxaZw1cQALRqcSFR6gDepFYPv4j/DpvWY+PBau+xISBvs2TV6y/Ug5f3x7J18d6jgg57ShiVy7YCSnj0vHZvNQ8PDur+DLx8z8aXfA/F965joBYvr06WzevHmz1vqE+vqWoMEHJGgQfklr0xahcJcZpGzPB1Dwbef7h8fByEUmSBh1OsSley+tQeizvUXc+uoWjlW01gEenxLKE9NzGJzztgnaXFYBU6b6V3MAEcD/DlprymoayS6pIbu4msPFNRwqruFwSTXZxTUUVHZePzouMpTTx6Vz9qQBnJqZIiM0C9/K2wb/XgQO0+EB59wHM6/2bZp8YFN2KU+s288HO491GFE9IzWGa+ZlcMG0QUSEuvn/66anYeVNZn7i9+DCJ917/gAjQUMAU0ptmjYketqmPy6GuP4Qm95xGpsOEXF9uu638ACHw4yQWXrQVHspaZ4egNJDZhTNriRntjZiHjrH86Mu9zGl1Q389s3tvLsjv2VdWIji1qVj+O6YCPodfp+wrOVw6FPQDhdnUDDsFFOFyU8DCIdDc6yyjuxiExhkF9eQXVJjBQjVVNY1dftcCVFhLB1vAoW5o5Ld/+AhRG/Ym+A/p0HeVrM8dC5c+U6fLn3dX1jFfz49wOubcmmwt/3uSo2L4KpTRnDprKHu65Ag5yv47+lmPm08/PwL95w3QEnQEMCUUpumDbBN23RNbNc7hkV3DCTi0iG2v9O0P0Ql9ekvI9GOvQnKc9oGA80BQulBaHJdpcOlkHDTe0/mGaZ9gtQL9TitNa9uPMKdK7+lpqFjyUJsRCgjoqs5y7aJhfbPGFu3DRsdAwitbDQNnkvIpAuwjV8GsWneSD4ADU0OcstqOWSVFmRbpQWHimvIKamhvslVwHN8oTbF4H5RzM5I5uxJA5gzMpmwEPnuE37mswdMGySAkAj42XpIGeXTJPmLgoo6nlp/iOc2ZHd4QRAbEcqls4Zy1Skj6J9wgiOu11XA34aYeVso/DavT7/kkqAhgHU7aOguW6hpeNo+oIhNswIOa11MWp/+TxNUGuugLLs1MGguMSg9aBooO7r/traN8FjTJeqgaaZEYcQCiHDTfSp65GBRNTe/vIWtOWVd7pdCOWeGfMU5ti+ZZcvCpjp+p9tRbA2ZxNfRC9iVtJDwuDT6xYSTFBNGUkwESTFh9IsOJykmnH4x4cRFhB63oWJ1fVNLMJDdrhrR0bJaHL38aYkKC2FYcjRDk6IZnhLD0KRohiVHMywphoGJkYRKkCD8WdE++Ndc0wMawJI74dRbfJkiv1RZ18hLX+Xw388Okl/R9kVWWIji/KmDuOZEe1x6YKJ5gQbw8w1tx2/oYyRoCGBKqU3TJo2btun1f0BlPlTlQ+UxqLI+lflm2pM3wt0VndxaahGbbvqMjoyHyATziXCaj4yHyESzLkR6HPG6+sq2wUBLcHAQKnJp6ZGjp6KSIGmEKTXoZ02TMsy6mFSpEudHGu0Onlh3gJVbj1JU1UBpTQP2Lp7GUykzAUTIl5ysdrkOILTiC8d43nHM5n37TEro2LNVqE2ZoMIKJEwwEUZNvd1qb1BD0Qn0v94vOoxhyTFWMBDN0Ob55GhSYyM837OKEJ7gcMDT58Dh9WZ5wBS4erX8fnahocnB8i25PLHuAHsLOg6+uGRcc49LvRjs8vmLTEceAKfcDDOugsShffI3ToKGANathtBaQ125UyBxzAou8p0CiwKzrq7c84kOi3EKJDoLMJrXu1gXGtm3/qNqbYK+hhpotD4N1da0BhqrW7c1VJvuTBuroaqwNUCoLuz99eMGOAUEw50ChBEmUBQBSWtNRV0TpdUNlNQ0UFJlps3LpdUNlFQ3UlrTAJX5zKj5lNMc65mpdrsMIJq0jfWOCVYAMYMy3NOPulIwID6SoVYJwdDkaIZbgcHQ5GjiI2UgNRFE7E3me3v7q/DJ3806Wyhcsxb6T/JlygKGw6FZs7uAxz7Zz9eHSjtsnz6sH9fOz2BJT3pc+vD38Pk/2q6LGwBDZsHQ2Wbaf3KfCOokaAhgbu89qbG2bWBRVdCuBMOaVhfS67fTJ8oW1kmAEQ/hMaBCTLsMFWK+bG0h1nxIu/lQUDan+fbbe7OvtV9TbccH+k4f8p3WN3ZynMuGqm6ibJAwxEWJwQjoN9zkqRBAk91BRcFhmna8ReTelcQXbHS5nx0bG22TWdk0i7cbph03gAgLUQzpZ1UdSnaqRpQczeB+0dJ7kQg+DTVQvBcK90DRbijcDUV7oHg/OBrb7jvvNlj8O9+kM8Btyi7h8U8O8MHOYx22ZaTGcO38DM4/qRs9Lh3eAE+eSZfPPWHRMGi6FUTMhiEzzbNJkJGgwcOUUoOBPwBnAslAHvAWcJfWumMY3LNz+6bLVXsT1BS1llZUFZhSivoKM60rN42Hmufrndb5Ktjoy0LCTQDgHBA0BwiJQ6V9iuid8lzIWgHfvgk5X7rcRdtCqR86j+Jh53A49TQK7VGUVNUTHhrCcKu0YEBCFCGe6mNdCF+qLe0YGBTutga07MZvYcpouPZTCDvBxrx93L4C0+PSG5s79riUFhfBVaeaHpe6LLnM3w67VkHOBsj5Gho6GfemhYL0CW1LI4KgSpMEDR6klBoJrAfSgOXALuBkYBGwGzhFa118AucPrHEaHA5oqHIKJiq6CDDar7fm7Q2+/iu8LyTcvMUIjzHTsKjW+fBoU+UrPNraZq2LTGwNEOIHmdIQITyl/AjsXA7fvgVHvnK9jy0MMhaacSAyFkJMCoRGeC+NQniC1uYFWtHutgFC4W6oLuj5+eIHmWAhfQLMutY8aAq3OFZRx1OfH+L5DdlU1nfscemyWUP5UXd6XHLY4di35mXJ4Q1m2txQuitBUKVJggYPUkq9DywFbtRaP+y0/n7gFuBxrfVPT+D8gRU0uENjXbtAoqx1ubHODFrlsLdO28w3mao+LfO92ddhljs7TtshNMp6iI9q+0Df5kHf1ToXD/9hMQH3pSL6uLIcK4B4E3JdV2FqERYD0UmmfUx0kmlc7zzfZmqtj0iQrqGF9znspqe5wnalBkV7jj8uTXvKZkp6U8eYACF1TOt8hHvaA4nOVdY18sKXh3ny84NtBsEEU13ygpNMj0uj0nrwb1Gea0ohDn8Jh7+AYzuOX7U4LAYGTzfVmYbOgsH+X6VJggYPUUplAPuBQ8BIrVvvHqVUHKaakgLStNbVvbxG3wsahBCBozTbBBA734JcN31PKZspSesssOgs6AiLcs/1RfBy2E3HE6WHOgYGRXtbuz/trpAISMlsFxiMgeSRUsrmB+qb7CzfcpQn1h1gn8sel9L56YIMZvSmx6X6SjiysbU04sjGnldpGjrbtDn0oypNEjR4iFLqauDfwBNa62tdbG8uhViitf64l9eQoEEIERhKD5kAImul6e63ttSUzHlLaJRTQGEFFyHhrR0mdOj8wKmjg/bLzh0tuOwoIdSUhrRZDrH2C227r7KZ10dAy0zLQ0Ivl3t1jDZVbXo1PdHjnab2BuvTaH0anKYNpqS3zT7t97WWHS6Obdmvk3172+lERHzHwCB1NCQOk6qhAcDh0KzeZXpc2pjtuselny4YyeKxad3vcanDRXpbpWmgKYVoLo1In+TT2gcSNHiIUuoe4DbgNq31fS62/xO4Dvi51vpfvbyGBA1CiMDkcJhqHbUlJoCoKTXzNSVtp7WlTutKu/G2TggPiUnrGBikjDGDn/rR22DRexsPlfD4ugN86KLHpcToMMLdOChkmi5iit7FFL2LyY5dZJJNCF0HrjVE8q3KZMQl95AyZo7b0tJd7goapNJ1R80V0zob/KB5feLxTqSU6iwqmJKVlcX06Sf0byeEEAEgHEgHnWa1IWpqbXPUZt4Out2yownpuU10j4KQMDMmUGikqUbUPG8LwdQszgPW+jaZwqNimxwUVdZTVttA8zvxPA9c50OneRv9iFL1xFBPNHVEq3psHYKIKqCIplevIjQi2gMp6lpWVhbA8BM9jwQNPedUJtxrttraWvvmzZu3uiNBgrHWdJdPUxE8JD/dS/LTvSQ/3SfI8tIO1PkyAUGWnz4XMPlZB3SrL/48n/0pw4EetvzvSIKGjppLEjprCh/fbr9OdVYM1FwCcaLFRMKQ/HQvyU/3kvx0L8lP95G8dC/JT/eS/PQ/0v9dR7ut6ehOtmda0z1eSIsQQgghhBA+J0FDR2us6VKlVJv8sbpcPQWoBTZ4O2FCCCGEEEL4ggQN7Wit9wMfYOp/Xddu811ADPBsb8doEEIIIYQQItBImwbXfg6sBx5SSi0GsoBZwCJMtaT/9WHahBBCCCGE8CopaXDBKm2YATyNCRZuBUYCDwFztNbFvkudEEIIIYQQ3iWDuwkhhBBCCCG6JCUNQgghhBBCiC5J0CCEEEIIIYTokgQNQgghhBBCiC5J0CCEEEIIIYTokgQNQgghhBBCiC5J0CCEEEIIIYTokgQNQgghhBBCiC5J0OAGSqnBSqknlVJHlVL1SqlDSqkHlVL9enieJOu4Q9Z5jlrnHeyptPsTpVSyUupqpdSbSql9SqlapVS5UuozpdSPlVLdvl+tPNSdfPI9+Xf4E3fmg7vu80CllLqyi7xs/ti7ea4+c38qpS5USj2slPpUKVVh/Y3PHeeYuUqpVUqpEqVUjVJqm1LqZqVUSC+uP14p9YpSqkApVaeU2q2UukspFdX7v8p3epKfSqlMpdSvlFKrlVI5SqkGpdQxpdRypdSiHl53+HHu/Zfc8xd6Vw/z0+154M573R/0MD+f7sZ36sfdvG5Q3p/+JtTXCQh0SqmRwHogDVgO7AJOBm4CzlRKndKdEaSVUsnWeUYDq4GXgLHAj4BzlFJztNYHPPNX+I2LgH8BecAa4DCQDnwX+A9wllLqIt39EQnLgQddrK9yQ1oDyQnng7vu8wC3Bbirk23zgNOAd3twvr5yf94BTMH8XUcw32udUkqdB7wO1AEvAyXAucADwCmY74luUUrNwnyfhgGvATmYf6ffA4uVUou11vU9/Ht8rSf5+Ufg+8BOYBUmL8cAy4BlSqmbtNYP9fD6W4G3XKzf0cPz+Ise3Z8Wt+SBO+91P9KT/HwLONTJtiuADHr2nQrBd3/6F621fE7gA7wPaOCGduvvt9Y/1s3zPG7tf3+79Tda69/z9d/qhbw8DfOFaWu3vj8mgNDA97p5rkPAIV//Tb7+uCsf3HWfB+sH+MLKh2Xe/HcJhA+wCMgEFLDQyqfnOtk3HigA6oEZTusjMUGrBi7p5nVDMA/Lbf5dMCXsr1nrf+3r/PFwfl4JnORi/QKgwcrnAd287nDrWk/7Og98mJ9uywN33uv+9OlJfnZxjkSgxsqblG4eE5T3p799pHrSCVBKZQBLMQ8Aj7Tb/P+AauAKpVTMcc4Tg4mqq63jnP3TOv8Z1vWCltZ6tdZ6pdba0W59PvCYtbjQ6wnr49x1nwcrpdREYDaQC7zj4+T4Ha31Gq31Xm39sh/HhUAq8JLWeqPTOeowbzABftbNSy8AxgHrtNYrnM7lAG63Fn+qlFLdPJ9f6El+aq2f1lp/42L9J8BaIByY6/5UBo4e3p/u5M573W+4KT+vAKKAN7TWRW5KmnADqZ50Yk6zph+4eNCtVEp9jnnYmg10VS9vDuY/yAda68p253EopT4ArsFE8MFeRakzjda0qQfHRCilLgeGYh5st2EeILpV7zyInGg+uOs+D1bXWtP/9vDekvuzo+Z77T0X29Zh3j7OVUpF6ONXK+r0XFrrA0qpPZjqoBnA/l6mN5D15jsVYKBS6logGSgGvtBab3NryvyfO/LAnfd6sPmJNX2iF8fK/elBEjScmDHWdE8n2/diHqZG0/XDVHfOg3WePkcpFQr8wFp09QXbmf7A/7Vbd1Ap9SPrTVtfcaL54K77POhYjWkvBxyYdjc9IfdnR53ea1rrJqXUQWAC5kE/q7fnsuzF3LOj6WNBg1JqGLAY82C6roeHn259nM+3Fvih1vqwWxLo/9yRB+6814OGUmoOMAnYo7Ve04tTyP3pQVI96cQkWNPyTrY3r0/00nmC1d+AicAqrfX73TzmKcyPYn8gBvMl9Dim3uO7SqkpHkinP3JHPsj92bmLMX/3u1rrnB4cJ/ena+681+S+dUEpFQE8D0QAd2qtS7t5aA2mYfV0oJ/1WYDptGIh8HEfqKLozjyQ+9O1a6zpv3t4nNyfXiBBg2c115U90bqS7jpPwFFK3Qjciumt54ruHqe1vstqI3FMa12jtd6htf4ppuFuFHCnRxLsZ7yUD332/qT1B+7xnhwk92evufNe63P3rdWN5/9heuZ5Gbi3u8dqrQu01r/XWm/WWpdZn3WYUsYvgVHA1Z5It7/wch70xfszAfMipgF4uifHyv3pHRI0nJjmNwEJnWyPb7efp88TVJRS1wH/wPSAskhrXeKG0zY3qJ7vhnMFsp7kg9yfLiilxmMakR7BdGfpDn39/nTnvSb3rRMrYHgO043nK8Dl7mj8q7VuorVqXp+8b3uZB3J/dnQ5EI0bG0DL/eleEjScmN3WtLO2BpnWtLM6te4+T9BQSt2M6TlqByZgcNeAVwXWtK8XU/YkH+T+dK23DaC70tfvz07vNatt0whMw93udAgh963FyrsXgUuAF4BLrYcpdym0pn31voWe54E77/Vg0dwAukclt90g96ebSNBwYpob6SxV7UYrVkrFYYqAa4ENxznPBmu/U6zjnM9jwxSvOV8vqCmlfoUZ3GYLJmAoOM4hPTHHmvalL2JXepIP7rrPg4ZSKhJTXc4B/NeNp+7r9+dqa3qmi23zMW8h13ezN5lOz2V1IzwayCbI81opFY4Zl+Ii4FngCg/00DXbmgZ1Xh5HT/PAnfd6wLMGYpyCaQC91s2nl/vTTSRoOAFa6/3AB5jGi9e123wXJqp9Vmtd3bxSKTVWKdVmhEStdRWmnmkMHesyX2+d/30d/CNCo5T6Habh8yZgcVdFlEqpMCs/R7ZbP0EpleRi/2GY0gswxfRBraf50Fl+9uY+7wMuwjS0W9VZA2i5P3vlNaAIuEQpNaN5pRWk/cla/JfzAUqpaCufh7Y71yeYXmfmK6WWOe1vA/5uLT7mg/75vcZq9PwmcB4muP1R+26TXRyTYOXngHbrZ1kBSPv9TwNusRaD+r7tTR50lp/04l4Pcs3tw7rsZlXuT99SQfx96RXWA8F6IA1YjvmRmoUZU2EPMFdrXey0vwbQWqt250m2zjMa8wbiK8zAROdhqizMtR7egpZS6oeYxk924GFc1+U8pLV+2tp/OHAQyNZaD3c6z53ArzFvyA8ClcBI4BzMaJurgAu01g2e+Dv8RU/zobP8tLb16D4PdkqpT4FTMSMNr+xkn+HI/YlS6nzgfGuxP3AG5o3fp9a6Iq31be32fw2oA14CSoBlmC4qXwMudn7QV0otxOTlJ1rrhe2uPQvzfRpmHXsY02vVDOBzzIuJgHqT25P8VEo9hRkVugh4FNeNatc6v9lVSl2J6d3rGa31lU7r12K6AF2LaccDMJnW8QZ+p7VuftgNGD3Mz7X0MA86y0+na3f7Xg8EPf3/bh0TDxzF/D8ddJyXhVfSh+5Pv6P9YFjqQP8AQzA3cR6m1X82pgFvkot9tcl2l+dJso7Lts6TBzwJDPb13+ilfLyzOX+6+Kx12n+4te5Qu/MswNTf3QWUYQYxKgQ+xIz3oHz9t3opP3uUD53lp9P2bt/nwfzBBPMayAFCuthP7k/drf/XHe43TJW3VUAppurbdszbwg75jelOsc13Q7vt44FXMQ/O9Zgg9y4gytd54+n8xDxAHe879c5257/SWv90u/U/Bt7GjAxfZeXlYUwvTPN8nS9eys8e50Fn+dmbez0QPr38//4za9uL3Th/n7o//e0jJQ1CCCGEEEKILkmbBiGEEEIIIUSXJGgQQgghhBBCdEmCBiGEEEIIIUSXJGgQQgghhBBCdEmCBiGEEEIIIUSXJGgQQgghhBBCdEmCBiGEEEIIIUSXJGgQQgghhBBCdEmCBiGEEEIIIUSXJGgQQgghhBBCdEmCBiGEEEIIIUSXJGgQQgjhU0qptUop7et0CCGE6JwEDUIIIdxCKaV7+LnS12kWQgjRPaG+ToAQQoigcZeLdTcDCcA/gLJ227ZY0x8A0R5MlxBCiBOktJYSYSGEEJ6hlDoEDANGaK0P+TY1QgghekuqJwkhhPApV20alFILrSpMdyqlZiil3lNKlSulSpVSryulhlj7ZSilXlJKFSqlapVSa5RSUzq5TrRS6jdKqS1KqWqlVJVS6gul1P++KYaBAAACcElEQVR44+8UQohAJkGDEEIIfzYT+NSa/zfwFfBd4GOl1FhreTDwLPAOsAD4UCkV63wSpVQi8BnwF8AOPAk8A6QCLyil/uT5P0UIIQKXtGkQQgjhz84GLtdaP9+8Qin1X+AqYD1wn9b6z07bfgf8Afgxph1FsweBk4Bfaa3vdto/EngL+K1S6jWt9RaEEEJ0ICUNQggh/NlnzgGD5RlrWg78rd22Z63p1OYVSqlk4HJgo3PAAKC1rgN+BSjgUnclWgghgo2UNAghhPBnG12sO2pNt2it7e225VrTwU7rZgIhgFZK3enifGHWdFxvEymEEMFOggYhhBD+rNzFuqbOtmmtm5RS0BoIACRb05nWpzOxXWwTQog+TaonCSGECHbNwcUDWmvVxWeRT1MphBB+TIIGIYQQwe4rwAHM83VChBAiUEnQIIQQIqhprQuA54EZSqnfKaU6VM1VSo1USo3wfuqEECIwSJsGIYQQfcH1QCamO9YrlFKfAceAgZgG0DOB/wEO+iyFQgjhxyRoEEIIEfS01hVKqQXANZiuVb8HRGICh73ALcCHvkuhEEL4N6W19nUahBBCCCGEEH5M2jQIIYQQQgghuiRBgxBCCCGEEKJLEjQIIYQQQgghuiRBgxBCCCGEEKJLEjQIIYQQQgghuiRBgxBCCCGEEKJLEjQIIYQQQgghuiRBgxBCCCGEEKJLEjQIIYQQQgghuiRBgxBCCCGEEKJLEjQIIYQQQgghuiRBgxBCCCGEEKJLEjQIIYQQQgghuiRBgxBCCCGEEKJLEjQIIYQQQgghuiRBgxBCCCGEEKJLEjQIIYQQQgghuvT/AXp2GBqO9i7FAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "image/png": { "height": 263, "width": 390 } }, "output_type": "display_data" } ], "source": [ "λ, a, b = 2, 0.05, 0.5\n", "N0, P0 = 40, 15\n", "\n", "N, P = multistep(N0, P0, λ, a, b, tmax=20)\n", "plt.plot(N, label='N')\n", "plt.plot(P, label='P')\n", "plt.xlabel('Time')\n", "plt.ylabel('Desnity')\n", "plt.xlim(0, None)\n", "plt.ylim(0, None)\n", "plt.legend();" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python [conda env:scipy]", "language": "python", "name": "conda-env-scipy-py" }, "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.7" } }, "nbformat": 4, "nbformat_minor": 2 }