{ "metadata": { "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.8.6-final" }, "orig_nbformat": 2, "kernelspec": { "name": "python3", "display_name": "Python 3", "language": "python" } }, "nbformat": 4, "nbformat_minor": 2, "cells": [ { "source": [ "# SIR model" ], "cell_type": "markdown", "metadata": {} }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [], "source": [ "from epistochmodels import SIR, SIR_Dem" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [ { "output_type": "stream", "name": "stdout", "text": [ "Populating the interactive namespace from numpy and matplotlib\n" ] } ], "source": [ "%pylab inline" ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [], "source": [ "pop = 1000\n", "beta = 0.7\n", "gamma = 0.2\n", "model = SIR(pop,beta,gamma)" ] }, { "cell_type": "code", "execution_count": 20, "metadata": {}, "outputs": [], "source": [ "model.initialize(990,10,0)" ] }, { "cell_type": "code", "execution_count": 21, "metadata": {}, "outputs": [ { "output_type": "stream", "name": "stdout", "text": [ "CPU times: user 2.66 ms, sys: 0 ns, total: 2.66 ms\nWall time: 2.67 ms\n" ] } ], "source": [ "%%time\n", "res = model.run(0,1000)" ] }, { "cell_type": "code", "execution_count": 22, "metadata": {}, "outputs": [ { "output_type": "display_data", "data": { "text/plain": "
", "image/svg+xml": "\n\n\n\n \n \n \n \n 2020-12-23T18:14:51.271293\n image/svg+xml\n \n \n Matplotlib v3.3.3, https://matplotlib.org/\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n\n", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABBgElEQVR4nO3dd3wUdf7H8dd3N71XAiRAQi8BgVBEihRRQKTYsGM/y1kPxXaKop7lbKh3P7EiemIBAUFEehGQ3kMJPRBICOkhbff7+2MWCZ1sNpnd5PN8PPaxM7Mz33lnknx2dnbmO0prjRBCiNrBYnYAIYQQ1UeKvhBC1CJS9IUQohaRoi+EELWIFH0hhKhFvMwOcD5RUVE6Pj7e6eULCgoIDAx0XSAXcMdMILkqyh1zuWMmkFwV5Ypca9asOaq1jj7ri1prt30kJSXpyliwYEGllq8K7phJa8lVUe6Yyx0zaS25KsoVuYDV+hx1VQ7vCCFELSJFXwghahEp+kIIUYtcsOgrpb5QSqUrpTaXmxahlJqjlNrpeA53TFdKqXFKqRSl1EalVMdyy4x0zL9TKTWyan4cIYQQ53Mxe/pfAQNOm/YMME9r3QyY5xgHGAg0czzuB/4LxpsE8BLQFegCvHTijUIIIUT1uWDR11ovBo6dNnkoMMExPAEYVm76144vkFcAYUqpesBVwByt9TGtdRYwhzPfSIQQQlQxpS+il02lVDwwQ2ud6BjP1lqHOYYVkKW1DlNKzQDe0Fovdbw2DxgN9Ab8tNavOqb/Eziutf73WdZ1P8anBGJiYpImTZrk9A+Xn59PUFCQ08tXBXfMBJKrotwxlztmAslVUa7I1adPnzVa605ne63SF2dprbVSymX9M2utxwPjATp16qR79+5d4TYy84v5fOkeSnMOMPiS9jSLCSLAxz2uQ1u4cCHO/ExVTXJVjDvmcsdMILkqqqpzOVsJjyil6mmt0xyHb9Id0w8CDcrNF+eYdhBjb7/89IVOrvuC9mYWMn7xbsrsmk83/YHVohjeIZb7ezWmWZ0gjA8nQghR+zh7yuZ04MQZOCOBaeWm3+E4i+dSIEdrnQbMBq5USoU7vsC90jGtSiQ1Cid57ABe7+HP/92WxMDEuvy0JpUr31tM//cWM33DoapatRBCuLUL7ukrpb7D2EuPUkqlYpyF8wbwg1LqHmAfcKNj9l+BQUAKUAjcBaC1PqaUGguscsz3itb69C+HXcrbaqF+kIXeiXUZkFiX5wYd59dNaXwwbyePfreO71ft5/qkOPq2jCHU37sqowghhNu4YNHXWt98jpf6nWVeDTx8jna+AL6oUDoXqh/mz709GzPysnj+9es2vvlzH3+kZAJQJ9iXQW3r8eSVzQnxkzcAIUTN5R7fblYjb6uFF69pzTMDW7JoRwY7juSx/kA2Xy3by+IdGcx98nIsFjnmL4S4MK01ZboMm92GTdsos5dh07azjp9rvvLTS2wlrMxbyfYN22kU0ogBCa4/s73WFf0TfLws9G8dQ//WMQCMnbGVz5fuYcLyvdzVPcHkdEKIC8kqyiK/NB+7tmPTNux249mmbX9NO1FM/5pH2/+atrFgI/m78ykuK6bIVkSxrZjismKOlx0ntySXnOIcckty/xo+Xnb8jCJt07aq+eGOwaCEQVL0q9ILV7fizz2ZvPzLVgpLbNzdPQF/H6vZsYSokbTWHC87Tk5xDunH08k8nsmFrhkqsZewN3cvu7J3sSt7FynZKZUPsuTMST4WH0J9QwnxCSHUN5S6AXVpHt6cAK8AvCxeeFm8sCorVosVL+WF1WLFqqynTD99vCLzJa9N5uo+V+NtrZpDzVL0HZRSTLy7K499v563Z2/n/xbu4rcnehEb5m92NCE8Qqm9lEP5hzh6/Ohfj8zjmWw6uokf5v5ATkkOucXGnnNucS5luqzC61AoYoNiaRrWlCvjr6R+YH0synJKET0xblGWM6aVH1+3Zh3du3bHz8sPX6vvXw+rxdydvTSvtCor+CBF/xThgT58fXcX/j17Ox8tSOH2z/5kzpOXY5Vj/KKWObEnnl2cTYmthFJ7KSX2EkptpZTajUeZvYyisiLWpa9jbfpaUrJSKLGXnNKOVVkJsgRR/3h9Qn1DqRdYjxCfEOPhazzXCahDpH8kXur85ciqrMQGx+Lv5ZodsXSfdOJD413SlieRon8Wo65qQYMIf0ZP3sQLUzfz+vBEuaBL1EglthI2ZGxgSeoSduXsIqMwg6ziLLKKsii2FV9UGz4WHzrU6cCtrW6lSVgTogOiifKPIso/ijDfMBYvWuyWV77WVlL0z2FE54ZsPZTLhOX7aFoniHt6yJe7wrMVlBYwa88skjOTOVRwiEP5h0jNS6XEXoK3xZv40HhiAmJoFt6MCL8Iwv3CCfMNw9fqi7fF23hYvU8OO8YbBjckwDvA7B9PXCQp+ufxwuDWJKfl8eZv2xiQWFeO7wuPUFRWxK7sXWw9tpW1R9aSedy4HmVL5hZyS3Lxs/qREJpAQmgCveJ60aFOBzrX7UywT7DJyUV1kKJ/Ht5WC2OHJXLV+4u54p1F3Nszgb9d3oQgX9lswn3kl+STfCyZ7zO/57nvniOvJO+v17yUF43DGuPv5U+riFZc3/x6rmh0BV4W+RuureQ3fwEt6gbzw9+68eH8nXw4P4Upaw8y6/GecuWuMF1afhrLDi3j9T9fNw7RKG96N+xNq4hWRPpHkhSTRFxQnOlnowj3IkX/InRJiGDiPV35YfUBnp2yiXFzd/LC4NZmxxK1kNaavbl7+WDtB8zbPw+AAK8A/tntn3jt82Jw78EmJxTuTop+BdzYqQFztx7hu5X7uaNbPA0j5csrUfW01mw8upGfdvzEHwf/ION4Bt4Wbx665CEub3A5zcKa4W31ZmHqQrOjCg8gRb+C7u/VmCU7j3LjJ8uZ9vfuxIT4mR1J1GBHjx/l1pm3cqjgEApFz7iePBD3AD1je1IvqJ7Z8YQHkqJfQZ3iI/jizs7cM2EVd325iul/746X1dnbEghxdruydzF2xVjWHFkDwJAmQ3i689OE+oaanEx4Oin6TujWJJIxQ9rw9E8bmb7hENd2jDM7kqgBcktymbdvHt8mf8v2rO34Wn156JKH6B7bnbZRbeUCQeESUvSdNLxDLJNW7mfM9C1c3a4evl5yhoRwTlZRFv9a+S9m7ZkFQFxQHPe2vZfrml1HXLDsUAjXkqLvJG+rhUf6NeOuL1fx7ORNvDuivdmRhAfRWnO44DBr09fy9qq3ySzKpEdsD3rH9ebK+CsJ9ws3O6KooaToV0Lv5tEMbV+fKesOckvXhnSKjzA7knBzR48fZWrKVKbvms6enD0ABPsE80LXFxjRcoTJ6URtIEW/EpRSjB7QkmnrD/HId+v46JYOJDWSwi/OtDNrJwsOLGDi1olkF2fTLqodz3R5hg51OtA8vLlcISuqjfylVVL9MH/evr4dT0/eyKgfN/Lb4z3l+L74y8zdM3lz5ZtkFWcB0LVeV0Z1GkXLiJYmJxO1lRR9F7ihUwOOl9p4cdoW5iWnM6itnD9d2x0uOMzYFWNZnLqYxqGNGd5sOIMSBtEiooXZ0UQtJ0XfRUZ0bsC4eSm8NjOZzvERRAf7mh1JmGRPzh5un3U7eSV5PNz+YW5pdQshPiFmxxICALmqyEV8vay8NjyRg9nHGfHJcuz289/vU9RMv+35jRt+uYGC0gL+0+8/PHDJA1LwhVuRPX0XuqpNXd68ri2jJ29iTvIRrmpT1+xIohqk5acxd/9clh5cyrJDy2gX3Y6xl42lcVhjs6MJcQYp+i52Xcc43p69g5/XHpSiX4MV2gr5eP3H/L73d3bn7AagcWhjbmt1Gw9c8oB0lyDclhR9F/OyWhjcrh5fLdvL5oM5JMbKP39NUmor5cstX/LZoc8oSi2iS70uDG86nN4NetfKm2wLzyNFvwrc2KkBXy3byz0TVjH14e7UC5XbLNYERwqO8OKyF1l2aBmt/VvzyhWvyNk4wuPIF7lVoHX9EKb/vTvpecWMnbHV7DiiknJLcvn3qn8zfNpwlh1axhNJT/BgnQel4AuPJHv6VaRdXBjXtKvP71sPk5KeR9M6ctNpT/Pr7l9ZdmgZM3fPpEyXkRSTxFOdnqJNVBsWLlxodjwhnCJFvwo9dVULZm1O4+MFu3hPOmTzKL/s+oXnlj4HQP9G/RnWdBi94nqZnEqIypOiX4UaRAQwILEef6QcRWst/aF7iOeWPMcvu3+hVUQr3u/zPvWD6psdSQiXqdQxfaXUE0qpLUqpzUqp75RSfkqpBKXUn0qpFKXU90opH8e8vo7xFMfr8S75CdzcpY0jSM8rZl9modlRxHkU24p5d827XPPzNfyy+xfqBtblqwFfScEXNY7TRV8pFQs8CnTSWicCVuAm4E3gPa11UyALuMexyD1AlmP6e475arwuju6W/9yTaXIScS7FtmKunXYtX27+EquyMrrzaGYOn0mAt9z4XtQ8lT17xwvwV0p5AQFAGtAX+Mnx+gRgmGN4qGMcx+v9VC043pEQFYiPl4Ux07dSUmY3O444TXJmMtdOu5b9efv5W7u/MXXYVG5rfRs+Vh+zowlRJZTWzvcRo5R6DHgNOA78DjwGrHDszaOUagDM0lonKqU2AwO01qmO13YBXbXWR09r837gfoCYmJikSZMmOZ0vPz+foKAgp5d3lS83F7MotYyH2/vSKqjILTKdzl221emqKlepLuXDIx+yp9i4kckVIVcwNHyo6bkqwx0zgeSqKFfk6tOnzxqtdaezvqi1duoBhAPzgWjAG5gK3AaklJunAbDZMbwZiCv32i4g6nzrSEpK0pWxYMGCSi3vKnlFpbrR6Bn65vHL9fz5882Oc1busq1O5+pcdrtdr0xbqYdNHaYTv0rULyx9Qe/P3a/tdrupuVzBHTNpLbkqyhW5gNX6HHW1MmfvXAHs0VpnACilpgDdgTCllJfWugyIAw465j/oeBNIdRwOCgVqxYHuIF8vbunakP/9uZ9Lw3zpY3agWuro8aM8vfhpVh1eRahvKP+94r/0iO1hdiwhqlVljunvBy5VSgU4js33A7YCC4DrHfOMBKY5hqc7xnG8Pt/xjlQrPHFFcwC2HZPj+mZYdGARg6YMYt2RddzX9j5+GfaLFHxRKzm9p6+1/lMp9ROwFigD1gHjgZnAJKXUq45pnzsW+RyYqJRKAY5hnOlTa0QH+9I8Jog9OXLqZnWbu28uzy19jhCfEMb3Hy9dHotarVIXZ2mtXwJeOm3ybqDLWeYtAm6ozPo8Xf/WMXy8YBcHs48TGyadsFWHVYdX8eySZ2kU0oj3er9Hg5AGZkcSwlTS4Vo1uq5jHABvzNpmcpLaodReyjOLnyHcL5wP+n4gBV8IpOhXq8bRQXSpa2X25sPkF5eZHafGm79/PunH03my05PEBsWaHUcItyBFv5r1iPWixGbnlw2HzI5So01NmcqoRaOI9IukZ2xPs+MI4Tak6FezJmFWABZsSzc5Sc2VXpjOy8teRqEYf+V4Ar0DzY4khNuQol/NAr0V13WMY+GODIpKbWbHqZF+2P4DNm1j5vCZNA9vbnYcIdyKFH0TDO8QS0mZnVmb08yOUuPsztnNxK0T6RnXU764FeIspOib4LImkUQE+rBkx9ELzywu2tx9c7l15q0UlhXyj07/MDuOEG5JbqJiAotFcVmTSJbIzVVc5unFTzNrzyyahjXlpW4v0ThULsAS4mxkT98kvZpFk5FXzMbUHLOjeLw5++Ywa88sWka05IfBP9C+TnuzIwnhtqTom+SK1jF4WxUfLUgxO4pHs2s7b696m1DfUN7u9TbeVm+zIwnh1qTomyQi0Ichl8QyZ+sR0vOKzI7jsX7c/iNpBWk8mfQk8aHxZscRwu1J0TfRrZc2BOCDuTtNTuKZtNb8b9v/aBXRiuFNh5sdRwiPIEXfRB0bhhPk68X/Vu6nFvUy7RKltlIemvcQu3N2M7TpUPkyXIiLJEXfZLdd2gitITktz+woHuXJRU+y9OBSbml5Cze1qFW9dAtRKVL0TXbnZfEALN6ZYW4QD7IrexdLUpfQLrodz3Z9FqvFanYkITyGFH2T1Q31o0VMMD+sPmB2FI+gtebfq/+Nl8WLdy5/x+w4QngcKfpuoFN8OLszCjhwTO6qdSFTU6b+dVinbmBds+MI4XGk6LuBm7sYZ/Hc8cVKk5O4t2NFx3hlxSvEBcXxSIdHzI4jhEeSou8GEmNDGda+PnuOFpBTWGp2HLdks9sYPm04ZfYynu36rFyEJYSTpOi7iWEdjDs7bTuca3IS92LTNn7b8xv9fuzHsaJjdKvXjcvqX2Z2LCE8lnS45iZa1g0BYNvhPLo2jjQ5jXvYkLGB1w+9Tvr+dGICYhjdeTQjWozAyyJ/tkI4S/573ERMiC8hfl7M3nKYkY7TOGuro8eP8tKyl1iSuoQwaxiv93idfg37EeAdYHY0ITyeFH03oZTi6nb1+G7lAfYcLSAhqnbe4i8tP427Z99NemE6t7a6lXb57RjYZKDZsYSoMeSYvht58PKmALw9e5vJSczz4boPSc1P5csBXzK6y2j8Lf5mRxKiRpGi70YaRgbQs1kUy3Zl1sq+eLKKspi3fx4DEwbSLrqd2XGEqJGk6LuZgYn1yC4sZX8tvFBraspUCssKua/tfWZHEaLGkqLvZtrFhQKw/kC2uUFMMGXnFNpFt6NZeDOzowhRY0nRdzMt6wbj723l+1W1qy+eRQcWsTd3L0ObDDU7ihA1mhR9N+NltRAT4suyXZkUlpSZHadabD+2neeWPkedgDoMbyY3QxGiKknRd0MP9m4CwKq9WSYnqXpaaz5c9yG5Jbm80fMNvC3SvYIQVUmKvhsa0KYeSsGafTW/6H+55UsWpS7ivrb30bluZ7PjCFHjVeriLKVUGPAZkAho4G5gO/A9EA/sBW7UWmcp4352HwCDgELgTq312sqsv6YKDfAmqWE4c7ce4cn+zc2OU2U+3/Q57699n/6N+kuvmcKtlJaWkpqaSlFRUbWvOzQ0lOTk5Iua18/Pj7i4OLy9L/4TcmWvyP0A+E1rfb1SygcIAJ4D5mmt31BKPQM8A4wGBgLNHI+uwH8dz+Is+rWK4c3ftjFt/UGGto81O47L7crexftr3wfg+a7Pyz1uhVtJTU0lODiY+Pj4av/bzMvLIzg4+ILzaa3JzMwkNTWVhISEi27f6cM7SqlQoBfwuSNAidY6GxgKTHDMNgEY5hgeCnytDSuAMKVUPWfXX9Nd29Eo9BOX7zM5iettP7adkb+NBGDKkClE+ksHc8K9FBUVERkZ6dY7I0opIiMjK/xpRDl75adSqj0wHtgKXAKsAR4DDmqtwxzzKCBLax2mlJoBvKG1Xup4bR4wWmu9+rR27wfuB4iJiUmaNGmSU/kA8vPzCQoKcnr5qlCRTF9uLmb1kTI+6htQ5X981bWtCmwFPJP6DF54cW+de2nj38YtclWUO+Zyx0zgmblCQ0Np2rRpNScy2Gw2rNaLv+9zSkoKOTk5p0zr06fPGq11p7MuoLV26gF0AsqAro7xD4CxQPZp82U5nmcAPcpNnwd0Ot86kpKSdGUsWLCgUstXhYpk+mbFXt1o9Ay9P7Og6gI5VMe2Kior0j2/66kTv0rUU3dOvahl3PF3qLV75nLHTFp7Zq6tW7dWX5DT5ObmVmj+s2UFVutz1NXKnL2TCqRqrf90jP8EdASOnDhs43hOd7x+EGhQbvk4xzRxDu1iwwDYkJptag5XKLYVc/uvt5NVnMXozqMZ2lQuwhLiQl577TXatGlDu3btaN++PX/++eeFF7oAp4u+1vowcEAp1cIxqR/GoZ7pwEjHtJHANMfwdOAOZbgUyNFapzm7/tqged0gfKwWNqXmXHhmN/fqildJPpbMyNYjua31bWbHEcLtLV++nBkzZrB27Vo2btzI3LlzadCgwYUXvIDKnr3zCPCt48yd3cBdGG8kPyil7gH2ATc65v0V43TNFIxTNu+q5LprPF8vKy3rBbPOw/vh+WTDJ0xNmUrD4IaM6jzK7DhCeIS0tDSioqLw9fUFICoqyiXtVqroa63XYxzbP12/s8yrgYcrs77aqEF4ADM3pVFQXEagr+fd86bUVsqErROID4ln8pDJZscRosJe/mULWw+59t7VreuH8NI15z+J4corr+SVV16hefPmXHHFFYwYMYLLL7+80uuWK3LdXFKjcAC2pnnmDdOf/+N58kryeKrzU/hYfcyOI4THCAoKYs2aNYwfP57o6GhGjBjBV199Vel2PW/XsZYZ3iGWsTO3smJXJp3jI8yOUyFHCo4wa88sEiMT6Va/m9lxhHDKhfbIq5LVaqV379707t2btm3bMmHCBO68885KtSl7+m4uPNCHlnVDWL470+woFWLXdl7/83UAnr/0eelITYgK2r59Ozt37vxrfP369TRq1KjS7cqevgfo1jiSL/7YQ1GpDT/vi79owyxaa55Y8ATzD8zn5pY30ybSvD0lITxVfn4+jzzyCNnZ2Xh5edG0aVPGjx9f6Xal6HuAjo3C+OIP+G3zYYZ1cP9+eJYeXMr8A/MZlDCIZ7o849aXsgvhrpKSkli2bJnL25XDOx7gqjZ1CfX35v25O8yOckF2bee5pc8R7hvOPy/9JxYlf2JCuBP5j/QA3lYLN3VpwN7MQpbvcu9j+3P3zSW7OJv72t1HkI/79bciRG0nRd9D3N+zMQALt6dfYE7zpOWn8Y9F/yA+JJ4RLUaYHUcIcRZS9D1EZJAvLesGk3w4z+wo5/T26rcBeLrz03JOvhBuSoq+B7kkLozFOzLYc7TA7Chn+GH7D8zZN4c7Wt9Bz7ieZscRQpyDFH0Pctulxjm6K9zsnP1vtn7D2BVj6VCnA48nPW52HCHEeUjR9yCJsSGE+nszL/mI2VH+Umov5btt3+Fr9eWjfh/JRVhCuEhV3XhGir4HUUrRt2Ud5ianU1hSZnYcADZlbGJ/3n6e6vQUIT4hZscRQlyAFH0PMyCxLgDb3eQL3TdWvgFAn4Z9TE4ihLgYckWuh2ldz9ibXr47kw4Nw03N8sP2H0g+lkyITwjR/tGmZhGiysx6Bg5vcm2bddvCwDdc2+ZFkj19D9MgIoCezaJ467ftpOUcNy1HcmYyY1eMJS4ojhnDZ0hXC0J4CNnT90D39WzMkp1HWbE7k+Ed4kzJMG7dOPy9/Jk4aCLhfuZ+4hCiSpm0R15VZE/fA3VrEomXRTFnqzln8UzcOpGlB5cyss1Iovxdcws3IUT1kKLvgbytFuLC/fl102FKbfZqXff69PW8teotgr2Dub317dW6biFE5UnR91ADEusBsONI9Z7F8/aqt/FSXvw89Gc5RVOIKpSfn18l7UrR91CD2xlF/8Cxwmpb59dbvmbj0Y08nvQ4MYEx1bZeIYTrSNH3UA0jAwBYvPNotaxPa82nmz7FoiwMaTKkWtYphHA9OXvHQ4X4edM4KpCd1XB4J78knw/WfkB2cTZjuo2Rs3WE8GBS9D1Yl4QIfl53sErXsTdnL9dNv44SewnNwpsxMGFgla5PCFG15PCOB0uICqS4zM7/LdpVJe1rrXly0ZOU2Et4t/e7TL5mMgHeAVWyLiFE9ZCi78Hu6BYPwPjFu7Hbtcvbn39gPjuzdvJoh0fp36i/XHUrRA0gRd+D+ftYeeHqVhwrKHH5TdO11ry96m18rb7clXiXS9sWQlyY1Wqlffv2JCYmcs0115Cdne2SdqXoe7g7L4snLtyf37Ycdmm7U3ZO4WD+QTrU6YCXRb76EaK6+fv7s379ejZv3kxERAQff/yxS9qVou/hvKwWBrWtR0p6vks7YJuzbw4Ab/SsWf2OCOGJunXrxsGDrjlpQ3bhaoDbL23E+MW7ef3XbXx4c4dKt1dUVsTKwyu5o/UdRPpHuiChEJ7rzZVvsu3YNpe22TKiJaO7jL6oeW02G/PmzeOee+5xybplT78GaBARwOXNo/llwyGWuuBirVl7ZlFqL6VHbA8XpBNCOOP48eO0b9+eunXrcuTIEfr37++Sdiu9p6+UsgKrgYNa68FKqQRgEhAJrAFu11qXKKV8ga+BJCATGKG13lvZ9QvDi9e0pt87i3hs0jqWP9sPHy/n3s8zCjP4cN2HtAhvwaX1LnVxSiE8z8XukbvaiWP6hYWFXHXVVXz88cc8+uijlW7XFXv6jwHJ5cbfBN7TWjcFsoATn0nuAbIc099zzCdcpEl0EGOHJZJZUMKyXc7v7X+++XMyjmfwwqUvyCmaQriBgIAAxo0bxzvvvENZWeXvjV2poq+UigOuBj5zjCugL/CTY5YJwDDH8FDHOI7X+ympKi51Q1Icvl4WFu9wrugXlRXx4/Yf6dugL+3rtHdtOCGE0zp06EC7du347rvvKt2W0tr5i3qUUj8B/wKCgVHAncAKx948SqkGwCytdaJSajMwQGud6nhtF9BVa330tDbvB+4HiImJSZo0aZLT+fLz8wkKCnJ6+apQ1Zn+vbqIzON2/tWzYlfO5ufnM6N4Bn/k/8HDdR6mpX/LKkpYMe74OwT3zOWOmcAzc4WGhtK0adNqTmSw2WxYrdaLnj8lJYWcnJxTpvXp02eN1rrTWRfQWjv1AAYD/3EM9wZmAFFASrl5GgCbHcObgbhyr+0Cos63jqSkJF0ZCxYsqNTyVaGqM322ZLduNHqG3ns0v0LLzZk/Ryd+laiv+fmaKkrmHHf8HWrtnrncMZPWnplr69at1RfkNLm5uRWa/2xZgdX6HHW1Mod3ugNDlFJ7Mb647Qt8AIQppU58QRwHnDi59KDjTQDH66EYX+gKF7qytdHP/ejJGyu03LSsaQDc1uo2l2cSQrgPp4u+1vpZrXWc1joeuAmYr7W+FVgAXO+YbSQwzTE83TGO4/X5jnck4UINIgLw8bKwYvcxbLsWw5hQ4zG2DhQeOznjD3cY0z++lO0r/8PCvIV0iunEdc2uMy27EO7EE8qTMxmr4uKs0cAkpdSrwDrgc8f0z4GJSqkU4BjGG4WoAo/2bcq/f9+BdeI1JyfaiuGthDNnzkjmu5VvQ3AQY7uPxWq5+GOJQtRUfn5+ZGZmEhkZ6bZnsWmtyczMxM/Pr0LLuaToa60XAgsdw7uBLmeZpwi4wRXrE+fx61P8feV4/n6Rfwe5FsXcAH/6FxQS9+N9sH8ZtL0B+r4A4fFVGlUIdxUXF0dqaioZGRnVvu6ioqKLLuR+fn7ExcVVqH3phqGmKC0CL19YOf6UyRmP7CI6MsoYKciEtxsbw62ugRaDmJe7nZw9U7gxNw/SlxmvbfrReIw59YwAIWoLb29vEhLO8sm4GixcuJAOHSrfncq5SNH3dFpDxnb4T9czXhpY+hZd/zjCmCGOoh8YeUYhnzvvYQDaFJec2fb+FdBQrsoVoiaRvnc82aH18HLYmQX/wWUwJoe4FklM33DonF/2rDmyhsWpixkQP4A1l0813hDG5MCzBwEFP90N2fur+IcQQlQnKfqe7IsBp453fQDunAkxbQDo2SyKYwUl/Lb57H3tP7/0eQCGNh166gu+QYCG3IPwfltj2tbpxtk+77aBjzobj/x0V/40QohqIEXfU5UUQpmj//wbJhh76APfhPiTPWNe3bYeABOW7z1j8fySfNIK0ujXsN/Ze9Mc+cvJ4Y86ww+3G8O5qXB0h/H45XHj8JIQwmNI0fdUrxsFHasPtBl21lkig3y5t0cCK3YfIzkt96/pWUVZ3DjjRuzazs0tbz57+wm94ImtxvBRx60Yg2JOnWf7TOPw0phQKf5CeAgp+p6oOO/k8ANLzzvr/b2Ms3X+NevkTSDGrRvHgbwDDGkyhC51zzi79qTQWLhnDkS1gMHvw6gdJ4/7/+O0e/Jm76voTyGEMIEUfU/07xbGc6d7ILrFeWeNDvalTrAvi3dkcKygBK01P+34iWCfYF7r8dqFLzxp0AX+vhI6nXZz9OAYo/jfPdsY/2IA/HgnZOyQvX4h3JgUfU8zJhRKC4zhK1664OxKKZ4ZaPSY+cA3a5ixewYA/Ru55i48NLwUvPwhLw22/AwfdzYO+Xx/G9jtrlmHEMJlpOh7irJimHzvyfGb/gd+oRe16LUd4+jZLIqVe44xYfP/AHi+6/Ouy/a3xWdOS/4FvrgSivNdtx4hRKVJ0fcUr9YxrpIFaNQdWgyq0OKvDE0ESxHbj22jcWhjfKw+rssW3dw41PNiFjy+6eT01FXwr1iY+7Ls9QvhJqToe4KUuaeO3zkTKtgJVEJUIHUTfgVLGS2973RdtvIsFghraLwBjNp5cvrSd40+fYQQppOi7+5spfCNo7vjUTuNgupEr3/Hio5R6rceCluxYWekazOeTVAdeCkbYh037/nqauP7iPRt511MCFG1pOi7u486nxwOquN0M++ufpcSewnDG/2NbYfzSM8tckG4C1AK7psHHUeenPafrvBeWzi4turXL4Q4gxR9d1daaDw/VrE7YZWXV5LH/APzaRXRils7Guflvz9v5wWWcqEh4+DJbcZ3EQA5++HTPsaef+au6sshhJCi7/a0hjbXQngjp5uYvGMyeSV5vHTZS7SoG0z/1jFMW3eQ4jKbC4NeQEg9uOtXeOYANL3i5PQPO8KnfeXcfiGqiRR9d7b/TyhIh8zK7ZVP3jmZtlFtaRNpdMR2S9eGFJTYztkRW5XyC4HbJhvH+7vcb0w7uAZW/Lf6swhRC0nRd1d2u3GeO0CPJ51uZsvRLezN3ctV8Vf9Na1Xs2gCfKz8+/ftZOQVVzapc5SCQW8bX057+cHsZ2H7b+ZkEaIWkaLvrr659uRwm+FON/Pdtu8AGJRw8rx+q0VxV/d4Dhw7zuu/JjvdtksE1YGRxlXCTP87lB43N48QNZwUfXeUfQB2LzCG753n1CmaADnFOUzbNY2r4q8iOiD6lNee7N+Cjg3DWL4r85w3Wak2DTrDHdOgIAPWfGVuFiFqOCn67uj9ROO5bjuI6+RUEwWlBTw07yEAbmh+5v3orRbF9UkNOJxbxJS1B52O6jKNe0Od1rD6S7OTCFGjSdF3N/kZJ4fvm+90M9dOu5aNGRsZGD+QrvXOvH8unLzJyqszt3K8zA3Onml6BRzdDis/NTuJEDWWFH13k7rKeB75C1i9nWoiOTOZQwWHiA+J563L3zrnfKEB3nx9dxeyCkt5cG4hNrvJhb/rA8bzr6OMM5eEEC4nRd/dTHLcyarB2ffOL8b3278H4MVuL15w3l7No2kUGQDAh/Or8YKtswmNhavfNYa/uBLGhJKweyJk7YWCo6ZGE6KmkKLvTsaU6yrZy9epJmx2G5N3TiYpJonOdTtfeAFg5qM9AfjPgl0cyjb57JnO98DN3xu3ZrR40Wj/T/DBJfB2E3ijoXHox+wvnoXwYFL03UVZufPlH3L+0Ma8/fMAGNx48EUvE+TrxZhufpTY7HyyyA26RWgxwLg144uZrOn4FsQ53ryKcoxDP3uXmJtPCA8mRd9dbJlqPHe6B+q0dKoJrTUTtk4AYED8gAotGx9qZcgl9Zm06gA5haVOrb8q5IW0gHvnGr2L9n3BmPjtDUbvo0KICpOi7w6O7YafHV0SdLrb6WZeXv4yGzM20jC4IUE+QRVefuRl8RSX2Rln9rH9c+n1lPEoK4KxUXJjFiGcIEXfHYzrYDwHRkPdRKeaWJG2gsk7J9MivAUzhs9wqo2kRuFc0iCMr5fvZXeGm97msPdzJ4d/HCnH94WoICn6Ztu14OTwUylON/PuauOslwkDJ6CcvIIXYNxN7Sm1afq+s4hx1dn98sWyWOC5NGM4ebpR+MtKzM0khAeRom+mjB0wcZgxXIkvb2fvnU3ysWRaRbQi0DuwUpEaRQZy+6WN8Pe28u6cHXy3cr/55++fzicAbp9qDG+dBtMeMjWOEJ7E6aKvlGqglFqglNqqlNqilHrMMT1CKTVHKbXT8RzumK6UUuOUUilKqY1KqY6u+iE81sflTqmMbuFUE4WlhYxaNAqAD/p84IpUjB2WyLoX+5MQFcizUzbR48357lf4m/QxvtwF44bx5T8xCSHOqTJ7+mXAP7TWrYFLgYeVUq2BZ4B5WutmwDzHOMBAoJnjcT9QuztQT119cvifmU53qjZ68WgArmt2HfWC6rkiGQB+3lZ+fKAbTaIDScspotU/fzO/Y7azuWO68TzpVsgz4f4AQngYp4u+1jpNa73WMZwHJAOxwFBggmO2CcAwx/BQ4GttWAGEKaVcV6U8zWf9jOfHN4HVy6kmthzdwsLUhfh7+fPPS//pwnCGqCBf5v2jN72aR1NiszN9wyGXr6PSGl8Og/4NpQXw3+5QUmh2IiHcmnLF3ptSKh5YDCQC+7XWYY7pCsjSWocppWYAb2itlzpemweM1lqvPq2t+zE+CRATE5M0adIkp3Pl5+cTFFTxUxerUn5+PnX1YTqt+QcAC3tPc7qtd9LeYW/JXl6JfYVwr/BK5zrXtios1Ty39DjZxZq3e/kTHVB9XwVd1O9Qazqtfoyggn2UWQNJq3cFu5o6f+qry3JVM3fMBJKrolyRq0+fPmu01mfvoldrXakHEASsAa51jGef9nqW43kG0KPc9HlAp/O1nZSUpCtjwYIFlVq+KiyYP1/rl0KMR9omp9vJKMzQiV8l6r4/9HVNrgtsq+nrD+pGo2foByaudsn6LlaFfofbZ5/ctr88UWWZtHbTvy03zKS15KooV+QCVutz1NVK7bIppbyBycC3WuspjslHThy2cTynO6YfBBqUWzzOMa1WqXt4XrkR587JB5i7by4AL3R9obKRLso1l9SnQ8Mw5mw9wsQV+6plnRXW/ErjxusAqz+Holxz8wjhhipz9o4CPgeStdbvlntpOjDSMTwSmFZu+h2Os3guBXK01mnOrt8j5R6i5fYPjeG7ZzvdTFFZERO3TiTCL4IesT1cFO7CPrk9iTK75tUZW7G729k8J/iFwPBPjOE3GkDyL+bmEcLNVGZPvztwO9BXKbXe8RgEvAH0V0rtBK5wjAP8CuwGUoBPgdp3cvW7rYzn6FbQ8FKnmtifu5/O33Zmf95+xnQbg7eTfe47o06wH68Pb0txmZ3/rdxfbeutsHYjoLmj76Hvb5OzeoQopzJn7yzVWiutdTutdXvH41etdabWup/WupnW+gqt9THH/Fpr/bDWuonWuq0+7QvcGq/8F+YPr3C6mTVH1gBwR+s76NOwT2VTVdjVbetRL9SPF6Zu5oWpmyizuWH/N0rBzZNggGN/4/MrpYM2IRzkitzq8nIYAMfCO1SqmZ92/ESEXwRPJD3hglAVFxrgzc8PdQfgmxX76fnWAkrK3LTwX/qgMZy9D2Y/D7YyczMJ4Qak6FeHKX/7a3Bby0edbmZr5lY2Ht1Is/BmeFmcO7ffFeqG+pH8ygAiA31Iyyli2Md/sP5Atml5zuuJrdCkL6z8BKY+IB20iVpPin5VKyuBjY5rDUZ8S4lvhNNNvbriVQCe6vSUK5JVir+PlTX/7M+Lg1uzNS2XYR//wedL95gd60yhsXDrZOM+BZt+hPG94dA644YsQtRCUvSrWvL0k8Mtr3a6mTJ7Gbuyd9GhTgdaRDjXT09VuLtHAjMe6UFYgDdjZ2xlxe5MsyOdyWKBgW/BkI8gY7tR+N9oCOu+kSt4Ra0jRb+qHdliPI/a6XT/OgDvrH6HwrJChjUd5ppcLpQYG8rMR3sSFuDNnV+udM+++K1e0PF2+PtKaHeTMW3aw/BWY5g4HPYtNzefENVEin5V0to4pBCbBEF1nG5mXfo6vk3+FoBrmlzjqnQuFRvmz+QHLwOg7zuLWJZy1ORE5xDWEK79BP55FBr3hrLjsGs+TBgMy/9j/M4ytstduUSNJUW/Kv32LOQcgOYDK9XMK8tfIcQ3hN+v+x1vS/Wdl19RTaKD+O+tSQC8O2eHyWkuwOoNd0yD5w7B/YvAXgaznzXOsvq4C7wSDllueuWxEJUgRb8qlB6HX5+CVZ8a490fc7qpxamLSclO4ZaWt7i06+Sq0qdlHYZcUp/V+7Lo9Opc5mw9Ynak8/MJhPrt4aEVxp5/qyEnX/ugHcwcJcVf1ChS9KvC7/+EleONvccR34KXj1PNZB7P5OF5DwNwQ/MbXJmwSr1+bVv6tIjmaH4x9329mn/NSjY70oXVaWXs+Y+YCC9mwbWfgk+w8cb9QTvjAq+Da8xOKUSlmXeyd02VtffkHv7w8dBqsNNN/XeDcZ+ZN3q+QXRAtAvCVY8gXy++vKsLxwpK6PzaXD5ZtJvDOUU8N6gVMSF+Zse7MIsF2t0IcZ0g74jxvczqL+DTvqCs8MBSsxMK4TTZ03el7P3wwSUnxy8Z4XRT49aO4/vt33NLy1u4urHzp3qaKSLQh9+f6EWv5tFMW3+Iq95fzM4jeWbHungRjaFRNxj8Ljy5FZLuAm2DT/sSc3gh2G1mJxSiwqTou9J/jLNXuOzRk/dvdcLhgsN8teUrQn1DebLTky4KZ44m0UF8fXcX/nNrR7ILSxnwwRJemraZnOMe1hdOSH245n14cDlEJNBq23vwSS9YMwFy3fCOYkKcgxR9V9n0E5TkQf0OcOVYp5ux2W08Mv8RtNZ8P/h7fK2+LgxpnkFt6zHt4e5EB/kyYfk+hn60lN+3eGDvlzGt4YE/2NJ6lPGF/S+PGr2nftILFr0FqXLcX7g3KfqVpTW81xYm32OM3zbl/POfx9HjR2k/sT3bjm2jX6N+xAbFuiike7ikQRjLn+3Lg72bsDezkPsnruHxSes4kltkdrSKsVjIqNMTHlkDD/wBV7wMFi9Y8Dp81hcWviF9/Ai3JV/kVtbG7yHH0bf8ffMhwPm+de6ebdzXtVVEK8Z2d/7TgjtTSjF6QEtGdotn1I8bmLr+EFPXH+KyJpE0igwgyNeLO7rFExfuj6rEFczVQinj7md1E6HH41B4DKY/Agv/BTtmw+0/g3+Y2SmFOIUU/co4uhN+dvSg+Vwa+AQ43dTM3TPZk7OHjnU6MmHgBBcFdF91Q/345t6uTFt/kMlrD7I7I5+d6flk5BXz6ZI9hPp781i/Znjl2tBau/8bABhv+Dd+Dcs+hLkvwfttoe8L0OH2Sv1tCOFKUvSdpbWxRwdw03eV+qeev38+zyx5hij/KMb1HeeigJ5haPtYhrY/eRhryc4MvlmxjxW7j/HKjK0AjFn+K/f0SODy5nXo1iQSq8WN3wAsVmOvPzDK6Ntn1tOw9mvjTmnRLeGSm8A32OyUohaTou+sFf+BzZOhz/PQcpDTzezI2sEzS54h1DeUz6/6nFDfUBeG9Dw9m0XTs1k0drtmzf4svpy9GhUUxadL9vDpkj1EBfnSMMKf3i3qcFf3eIL93LRbig63GbdsXPCacY7/kc3G9HmvQJ/n4JKb5dCPMIUUfWekzIXZzxnDlz1SqaZeWPoCFmXhp2t+om5gXReEqxksFkXn+AgKWvjQu3dHXi0o4feth1my8yiLtmewdn82787ZwYhODUiKD+eGpDj3OwQUGAWD34OrXjeO96dtgOUfwW/PGA/vAAiKgSZ9jL+jiMZmJxa1gBT9itqzGL67BQKj4dafwNvfqWa01oxaNIrkY8ncnXi3FPwLCA/0YUTnhozo3BCbXTNx+V4mrtjH96sP8P3qA7z7+w7aNwgjLMCbhpEBDG5bnwYRbvJlsLe/cTOX0FhoMdDozmHNV5CTCmVFxieB1V9AxzuMHlmbD4Bg+XsQVUOK/sU4sMroR6dRN6P7XS8/o2fGUOdPqXx8wePMPzAfgAcvedBVSWsFq0VxZ/cE7uyewJHcIr5atpfth/PYkpZDRl4xRaV23vptOxGBPgxrH8u1HWNpVS/kjDZMoZTRvUNcp5PTdsyGH+4wjv2v/RqC6xt7/Z3vgcRrzckpaiwp+heStgE+v+LUaV0frFTBf/GPF5l/YD51A+syZcgU/Lw8oD8aNxUT4sfoAS3/GrfbNQu2p7PtcB7frtjHF3/s4Ys/zryNY9M6QbSLC6VFTDB9WtYhISoQL4sy55NB86uMLp4PbzL+3jZ8B/uWGo/VX0DTftD5PvANqv5sosaRon8+BZnGlZan6z3a6SZX5a/i58yfAZgyZArBPnImhytZLIp+rWLo1yqG+3o2Zs2+LFbtPXbKPCVldpbszGDmxjSmlB3kX7O2AcaNYK5LiqNH0yi6JDh/vYVzwa1GF8/120PSSCjOh2XjYPMUmDsGNkyCfi9BsyuNu4AJ4ST56zmXspKTe/h9XoBeo4wOtirxDzdhywS+zvwagGlDp0nBr2I+Xha6NYmkW5PIM14bdVULbHbNnqMF/Lknk51H8lmwPZ1x83Yybt5O2sWF0r9VDG3jQgnxP3mGUJi/Nw0iAiiza0ptF767llUpLM4cSvINMs7y6f0szBoNGyfBpJvBNwSaXmH03trm2krdglPUTlL0z6akAKbcD8d2w+D3odNdxnQnC35WURb/WPQPVh1eRYJvAp8N+Yw6Ac7fPlG4htWiaFoniKZ1jMMmY2jDkdwiPpy/k2//3M/G1At0mvf7rItaT+t6ITSKDKBZTDBRQSfvrdAoMpBOjcLx8TJ6Q/G2nqVXFKVg0FtwxRhI/gV2zoa9f8CWKbD6S+PQUIOuENf5orIIIUW/vOJ8ow+dHb8Z482uPFnwK0hrTZm9jEWpi3hi4RMADEwYyADbACn4biwmxI9Xh7Xl+UGtySsuJeVIPiWOPXqtYVdGPkWlNnbv2UPjhITztrXtcB5FpXbsWrP9cB6ztxzGfp4ueeqH+tGkThDRwb4E+njRODoQyyl78l2hfleoZ6fl/u9ovf9/BO99AYCU+kPYb2vD115xaMup/9b+PlZC/LwJ9fcmxN/L8eyNr5cFRdV/Uiiza0rKLuJTkUVhUbjHGVc1mBT9E+x2+Kgz5Dm6yW17Iwz/xLmmtJ07Zt3BhowNf027rtl1jLlsDAsXLnRBWFHV/H2s+PtYqRN86pfsfVoab9gLFx6kd+9mFWqzqNRGQXEZAAeyjnMw6zh7juYDYLNDSkY++zML2J1RQEZe8V9vNmfXHmhPPTIZ5/MhlxycSVM1nYzDH7PN3pAMwlhlb8FOeyx7dD3CVD7e2Niu49Bm9LN4kZ+KALwsCqtFnXy2Wk4d/+vZUqmjW/n5hQStX+x8A1XkRK7Lm0fz7KBWLm9fij6ArQw+7W0U/DbD4drPnDqUY7PbSCtI49rp13K87DiJkYkMbjKYvg36esT9bUXV8vO24udtBSAyyJf2DcLOOW9RqY3Ckou9SctN5NtK2DPjXdrak+lWkI7l2EauLT7zDl8lPmEciu7JYf+mZPvEkOPXgMygpmhVdaVg9+7dNG58/gvPtNbYtfGpwGa3G8827RjXp053jJddxHcq53NUFxIV6X59Ip3IFRVUNd2q186ib7dD6kqjD/xju43T5AqPgn84XP9lhb4cO1xwmFWHV7EibQXTd03/a3qXul347MrP5KOqcEr5N4iL48PG+r3w6f2iMWorNW7onrrSOGwZEAGHN+Gz9mviMxYSX/JLuUWDjKuC63c0bhYTEGVcTRwYBQGRYPU1biHppIUqld69mzq9fFVZuHAhvXt3uvCM1ayqc9Weom8rg90LYM5LkL7l5PSwRpDQC1pdA62HnbPgl9hKmL9/Pn8e/pMSWwlZRVmkZKeQVpAGgJ/Vj56xPWkU0oi2UW0ZmDBQCr4wj9UbopoajxPaXg/9XzaGj2dDzgHI2A57lxoXiCX/ctamAON+AVYfo12rz6nDMW2gSV+od4nRqZxXzbjxT01Vo4u+spcZ3R9vngIrPoaicmdjXD7a6PI2rMEF28kpzuGan68hqzgLgJiAGLwsXtQNrEufBn1oF92Ofg37yUVWwnP4hxmPum2NNwOAolzITzc+9RYcNZ4Lj4GtxHiUFRufIGwlJ59LC2HfMtjyc7m2I4w+hYLqQFAMTbJKwHuDMd0/3Fivn2P9/uHgE3jhvHYbaLvxRiMqpdqLvlJqAPABYAU+01q/4fKVHFoPU+6nZ+YuWGx8cUZcF+j2EDTubfzBnWMvvNReysaMjezK3kVWURbbs7YzZ98cABqFNOJ/V/+PEJ+Qsy4rhEfzCzEeVPBQjN1uHCY9vBGO7jDeOPKPGM+pK6mfcwhSp557eZ9g4w0CfdqbimPYXmoUfIDQBhDdwvhEEdvRGFdW4/CTxcsxbD3PtJPDFluJcUpWLftEXq1FXyllBT4G+gOpwCql1HSt9VaXrigwCqKacSAgkUYd+xsfP2MSwWKhoLSA1amLOZB34JRFckty2ZCxgXXp6zhedvyU10J9Q+nboC+vdH/FpTGFqBEsljMPJZWzZMECel/WCY5nOR7ZxnNRtvFJ4sQbhMUKFu9yh5AcwxbHuLbDsV2Qvg32LAFbcaVi9wJYgtGX1omHtx94B578FPLXJ5Iw483JYjXeSE55nD6t3LjVu4LLVP0nmere0+8CpGitdwMopSYBQwGXFv0d9kKe9iugwJ5O4N7vYe/J144UHiG/NP+sy4X7hhMTEMONLW7ksvqXEeUfhY/VB38v53rSFEJg7En7BhuPsIauadNWBofWGoektM1x+MdmdIx44lCQ3TH+1+unTtudspPGDesbN7gvKzr5XFJovCEd2228QRVlG4exqklPiy+sj4FWQ+Cq11zevtLVeANnpdT1wACt9b2O8duBrlrrv5eb537gfoCYmJikSZMmVXg9GaUZTM+eTllZGV5ep76vKRSNfBuRFJCEj/Ip/wL+quq74s3PzycoyP06zpJcFeOOudwxE9SMXMpeitV2HKXtKG0r92y7yPGLW8ZiL0UXZhJIIflBTUhtMMSpn61Pnz5rtNZnPwVIa11tD+B6jOP4J8ZvBz461/xJSUm6MhYsWFCp5auCO2bSWnJVlDvmcsdMWkuuinJFLmC1Pkddre5L8w4C5U+XiXNME0IIUQ2qu+ivApoppRKUUj7ATcD0CywjhBDCRar1i1ytdZlS6u/AbIxTNr/QWm+5wGJCCCFcpNrP09da/wr8Wt3rFUIIUf2Hd4QQQphIir4QQtQiUvSFEKIWkaIvhBC1SLVekVtRSqkMYF8lmogCjroojqu4YyaQXBXljrncMRNIropyRa5GWuvos73g1kW/spRSq/W5LkU2iTtmAslVUe6Yyx0zgeSqqKrOJYd3hBCiFpGiL4QQtUhNL/rjzQ5wFu6YCSRXRbljLnfMBJKroqo0V40+pi+EEOJUNX1PXwghRDlS9IUQohapkUVfKTVAKbVdKZWilHrG7DwnKKX2KqU2KaXWK6VWm5jjC6VUulJqc7lpEUqpOUqpnY7ncDfJNUYpddCxzdYrpQZVc6YGSqkFSqmtSqktSqnHHNNN3V7nyWX29vJTSq1USm1w5HrZMT1BKfWn43/ye0fX6mZn+koptafctmpfXZlOy2dVSq1TSs1wjFfttjrX3VU89YHRZfMuoDHgA2wAWpudy5FtLxDlBjl6AR2BzeWmvQU84xh+BnjTTXKNAUaZuK3qAR0dw8HADqC12dvrPLnM3l4KCHIMewN/ApcCPwA3Oab/H/CgG2T6CrjerG1VLt+TwP+AGY7xKt1WNXFP/6+br2utS4ATN18XDlrrxcCx0yYPBSY4hicAw6ozE5wzl6m01mla67WO4TwgGYjF5O11nlym0oZ8x6i346GBvsBPjunVur3Ok8l0Sqk44GrgM8e4ooq3VU0s+rHAgXLjqbjBP4ODBn5XSq1x3ADencRordMcw4eBGDPDnObvSqmNjsM/1X7Y6QSlVDzQAWNP0W2212m5wOTt5ThcsR5IB+ZgfPLO1lqXOWap9v/J0zNprU9sq9cc2+o9pZRvdWZyeB94GrA7xiOp4m1VE4u+O+uhte4IDAQeVkr1MjvQ2Wjjc6Vb7AkB/wWaAO2BNOAdM0IopYKAycDjWuvc8q+Zub3Oksv07aW1tmmt22PcA7sL0LK6M5zu9ExKqUTgWYxsnYEIYHR1ZlJKDQbStdZrqnO9NbHou+3N17XWBx3P6cDPGP8Q7uKIUqoegOM53eQ8AGitjzj+Ye3Ap5iwzZRS3hiF9Vut9RTHZNO319lyucP2OkFrnQ0sALoBYUqpE3fqM+1/slymAY5DZFprXQx8SfVvq+7AEKXUXozD0H2BD6jibVUTi75b3nxdKRWolAo+MQxcCWw+/1LVajow0jE8EphmYpa/nCisDsOp5m3mOMb6OZCstX633Eumbq9z5XKD7RWtlApzDPsD/TG+b1gAXO+YrVq31zkybSv3pq0wjptX67bSWj+rtY7TWsdj1Kn5WutbqeptZfY311XxAAZhnM2wC3je7DyOTI0xziTaAGwxMxfwHcZH/1KMY4b3YBxLnAfsBOYCEW6SayKwCdiIUWjrVXOmHhiHbjYC6x2PQWZvr/PkMnt7tQPWOda/GXjRMb0xsBJIAX4EfN0g03zHttoMfIPjDB8zHkBvTp69U6XbSrphEEKIWqQmHt4RQghxDlL0hRCiFpGiL4QQtYgUfSGEqEWk6AshRC0iRV8IIWoRKfpCCFGL/D9zmh9NL47HWQAAAABJRU5ErkJggg==\n" }, "metadata": { "needs_background": "light" } } ], "source": [ "plot(res[0],res[1])\n", "legend(['S','I','R']);\n", "grid();" ] }, { "cell_type": "code", "execution_count": 23, "metadata": {}, "outputs": [], "source": [ "pop2 = 1000\n", "alpha = 0.05 # Mortality rate\n", "model2 = SIR_Dem(pop2, alpha,beta,gamma)" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [], "source": [ "model2.initialize(990,10,0)" ] }, { "cell_type": "code", "execution_count": 25, "metadata": {}, "outputs": [ { "output_type": "stream", "name": "stdout", "text": [ "CPU times: user 37.6 ms, sys: 7.04 ms, total: 44.6 ms\nWall time: 43 ms\n" ] } ], "source": [ "%%time\n", "res2 = model2.run(0,100)" ] }, { "cell_type": "code", "execution_count": 26, "metadata": {}, "outputs": [ { "output_type": "display_data", "data": { "text/plain": "
", "image/svg+xml": "\n\n\n\n \n \n \n \n 2020-12-23T18:14:58.365187\n image/svg+xml\n \n \n Matplotlib v3.3.3, https://matplotlib.org/\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n\n", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABH3ElEQVR4nO3dd3hURdvA4d+kEwIJNZQgofceOkLoXSwgioL15fPFjg2xYUGwYAdeUVCsINjoCEhEEOm9hx56DSSkZ74/zia7m2zalmySfe7ryrVnT52TTZ49Z87MM0prjRBCCM/g5e4CCCGEKDwS9IUQwoNI0BdCCA8iQV8IITyIBH0hhPAgPu4uQG4qVqyow8PD7d4+Pj6e0qVLO69ARZynnS/IOXsKOeeC2bJly0WtdSVby4p00A8PD2fz5s12bx8VFUVkZKTzClTEedr5gpyzp5BzLhil1PGclkn1jhBCeBAJ+kII4UEk6AshhAfJs05fKTULGASc11o3Nc0rD8wFwoFjwJ1a6ytKKQV8DAwAbgD3a623mra5D3jZtNu3tNaznXsqQgjhHCkpKcTExJCYmOi2MgQHB7Nv375c1wkICCAsLAxfX9987zc/D3K/Bj4DvrGYNw5YpbWerJQaZ3r/AtAfqGf6aQ9MB9qbviReAyIADWxRSi3QWl/Jd0mFEKKQxMTEUKZMGcLDwzGuZQvf9evXKVOmTI7LtdZcunSJmJgYatWqle/95lm9o7VeA1zOMnsIkHGlPhu41WL+N9rwLxCilKoK9AVWaK0vmwL9CqBfvksphBCFKDExkQoVKrgt4OeHUooKFSoU+G7E3iaboVrrM6bps0Coabo6cNJivRjTvJzmZ6OUGg2MBggNDSUqKsrOIkJcXJxD2xc3nna+IOfsKQr7nIODg4mLiyu049mSlpbG9evX81wvMTGxQL8bh9vpa621Uspp+Zm11jOAGQARERHannaqqWnp3Pn5evpWKcX/eVDbXmnL7BnknF1v3759uVatFIa8qncyBAQE0KpVq3zv197WO+dM1TaYXs+b5p8CalisF2aal9N8l9hw9DJbT1xl0kb3PYQRQghHTZw4kSZNmtC8eXNatmzJhg0bHN6nvVf6C4D7gMmm198t5j+mlJqD8SA3Vmt9Rim1HHhbKVXOtF4f4EX7i527TnUquGrXQghRKDZs2MCiRYvYunUr/v7+XLx4keTkZIf3m+eVvlLqR2A90EApFaOUeggj2PdWSh0CepneAywBjgDRwBfAGACt9WXgTWCT6ecN0zyXUErRvYGRdmJaVLSrDiOEEC5z7tw5KlasiL+/PwAVK1akWrVqDu83zyt9rfXdOSzqaWNdDTyaw35mAbMKVDoHNAsLYfWBC7y77ABjIusW1mGFECXM6wv3sPf0Nafus3G1srw2uEmu6/To0YP33nuP+vXr06tXL4YPH063bt0cPnaJ7ZE7tnf9zOlu760mfNxi1hy8QHq6jAkshCj6goKC2LJlCzNmzKBSpUoMHz6cr7/+2uH9Fuksm456NsKf9zcncfzSDQBGzdpI5TL+bHypl5tLJoQoLvK6Inclb29vIiMjiYyMpFmzZsyePZv777/foX2W2Ct9gMYVvLPNO389ifBxi91QGiGEyL9Dhw5x6NChzPfbt2+nZs2aDu+3RF/peynFsckDuZaYQqCvNyNnbmT9kUsApKdrvLyKbm87IYRni4uL47HHHuPq1av4+PhQt25dZsyY4fB+S3TQz1A2wEhG9MN/2vPYj9tYvPMM0/86zKPd5QGvEKJoatWqFf/884/T91uiq3eyUkrxYv+GAKzcd87NpRFCiMLnUUEfIKxcIADbTlx1b0GEEMINPC7oW0pOTXd3EYQQolB5ZND/aHhLABbuOO3eggghRCHzyKBfqYzRrfn1hXvcXBIhhChcHhn0O9etCMC1xFTCxy3muXk7MDJICCFEyeaRQT+reVtiqPXiEncXQwghMlWtWtUl+/XYoP/zfztmm3f+muTfF0KUbB4b9NvULM+xyQPZ/2Y/Imoaaf5PXrnh5lIJIYRreUSP3NwE+HpzZ9sabD5+hU9WRTP7wXbuLpIQoihZOg7O7nLuPqs0g/6T817PBTz2St/SoOZG3VmZAI//DhRClHAS5YBAP+PXsGjnGT4b4ebCCCGKFjddkbuKXOln8eu2GHcXQQghXEaCvknGmLq/bD3l5pIIIYTrSNA3+eqBdvj7eBF9Ps7dRRFCCM6cOeOS/UrQt5CUms6Z2ERW7JW0y0KIkkmCvoVGVcsCMOWPA24uiRBCuIYEfQtLn7wZgP1nr7u5JEII4RoS9LPo0ziUOpVKu7sYQgjhEhL0s6hbOYjjl26QkiYDrAghSh4J+lmU9vchNV2zLvqiu4sihBBOJ0E/i/a1ygMQcyXBzSURQniykJAQWrZsSdOmTRk8eDBXr151yn4l6GfRpFowAH8fuuDmkgghPFmpUqXYvn07u3fvpnz58kydOtUp+5Wgn0UpP28Alu+RtvpCiKKhY8eOnDrlnGwBknAtF1prlFLuLoYQwo3e2fgO+y/vd+o+G5ZvyAvtXsjXumlpaaxatYqHHnrIKceWK30bHugcDsDm41fcWxAhhMdKSEigZcuWVKlShXPnztG7d2+n7NehK32l1NPAw4AGdgEPAFWBOUAFYAswUmudrJTyB74B2gCXgOFa62OOHN9VutavxFfrjrH39DXahpd3d3GEEG6U3ytyZ8uo079x4wZ9+/Zl6tSpPPHEEw7v1+4rfaVUdeAJIEJr3RTwBu4C3gE+1FrXBa4AGfckDwFXTPM/NK1XJGW04Jm35aSbSyKE8HSBgYF88sknTJkyhdTUVIf352j1jg9QSinlAwQCZ4AewHzT8tnArabpIab3mJb3VEW0wjxjUJXdp665uSRCCAGtWrWiefPm/Pjjjw7vS2mt7d9YqSeBiUAC8AfwJPCv6WoepVQNYKnWuqlSajfQT2sdY1p2GGivtb6YZZ+jgdEAoaGhbebMmWN3+eLi4ggKCrJr2/uXxQMwpoU/7aoWj+fdjpxvcSXn7BkK+5yDg4OpW7duoR3PlrS0NLy9vfNcLzo6mtjYWKt53bt336K1jrC1vt3RTClVDuPqvRZwFZgH9LN3fxm01jOAGQARERE6MjLS7n1FRUVh9/bLFgMwbUcSz9/dy+4yFCaHzreYknP2DIV9zvv27aNMmTKFdjxbrl+/nq8yBAQE0KpVq3zv15HqnV7AUa31Ba11CvAL0BkIMVX3AIQBGY1LTwE1AEzLgzEe6BZJO17t4+4iCCGE0zkS9E8AHZRSgaa6+Z7AXmA1MNS0zn3A76bpBab3mJb/qR2pW3Kx4EDfzOkiXEwhhIsUh/97e8pod9DXWm/AeCC7FaO5phdGtcwLwFilVDRGs82Zpk1mAhVM88cC4+w9dmGpGOQPwLzNMli6EJ4kICCAS5cuFenAr7Xm0qVLBAQEFGg7h55Qaq1fA17LMvsI0M7GuonAMEeOV9jeHdqMB7/ezPM/7+TOtjXcXRwhRCEJCwsjJiaGCxfcl4MrMTExz4AeEBBAWFhYgfZbPJqluEm3+pXdXQQhhBv4+vpSq1Ytt5YhKiqqQA9o80vSMOTC20txf6dwSvt5F+nbPCGEyC8J+nm4FJ9MfHIaU1dHu7soQgjhMAn6ebitVTUA3v/joJtLIoQQjpOgn4ceDUPdXQQhhHAaCfoFkJwqg6ULIYo3Cfr5EOBr/JomLNzj5pIIIYRjJOjnw8qx3QBYf7jIZo0QQoh8kaCfD2HlAgE4ejHezSURQgjHSNAvoLR0aa8vhCi+JOgXUJ3xS9xdBCGEsJsE/Xxa/lRXdxdBCCEcJkE/nxpUce+ACkII4QwS9AugeVgwAGdjE91cEiGEsI8E/QJ4sLORda/DpFVuLokQQthHgn4BdKxTwd1FEEIIh0jQL4DQsuYBDS7HJ7uxJEIIYR8J+gVUPaQUIEFfCFE8SdAvoP/cbNTrrz8iKRmEEMWPBP0C6lS3IgB/7Dnr5pIIIUTBSdAvoDqVggD4+9BFN5dECCEKToJ+AXl7KXcXQQgh7CZB3w53tA4jtKy/u4shhBAFJkHfDsv3nOXctSRib6S4uyhCCFEgEvTtMCwiDICog+fdXBIhhCgYCfp26FDb6Jn75Jzt7i2IEEIUkAR9O/RtUsXdRRBCCLtI0LdTKV9vdxdBCCEKTIK+nbrWryiBXwhR7EjQt1PzsBASUtJITElzd1GEECLfJOjbqXxpP0ASrwkhiheHgr5SKkQpNV8ptV8ptU8p1VEpVV4ptUIpdcj0Ws60rlJKfaKUilZK7VRKtXbOKbhHRtC/GJfk5pIIIUT+OXql/zGwTGvdEGgB7APGAau01vWAVab3AP2Beqaf0cB0B4/tVjeVDwTg2KUbbi6JEELkn91BXykVDHQFZgJorZO11leBIcBs02qzgVtN00OAb7ThXyBEKVXV3uO7W1g5I6/+Ez9uc3NJnCM1PZWUtBS01sQlx9FsdjOazW7G1cSr/HXyL0YsHsGlBEknLURxp7TW9m2oVEtgBrAX4yp/C/AkcEprHWJaRwFXtNYhSqlFwGSt9VrTslXAC1rrzVn2OxrjToDQ0NA2c+bMsat8AHFxcQQFBdm9fV7uXxYPwNf9SrvsGAXhyPk+fvxxm/Pblm7LpvhNme9HVRhF26C2dh3DFVz9GRdFcs6ewZFz7t69+xatdYStZT4OlMkHaA08rrXeoJT6GHNVDgBaa62UKtC3itZ6BsaXCREREToyMtLuAkZFReHI9nm5+fAGDpy97tJjFIS95zt9+3Q4bnuZZcAH+ObSN3xz6RsAQgND+WPoH3gp97UHcPVnXBTJOXsGV52zI/+tMUCM1nqD6f18jC+BcxnVNqbXjAQ1p4AaFtuHmeYVW/Uql+HKjWTsvVsqCqZvn860HdMAmNV3Fp2rdQZgWP1hPND0AQCCfIPYOWonXcO6Wm177sY5WnzTgnPx5wq30EIIu9l9pa+1PquUOqmUaqC1PgD0xKjq2QvcB0w2vf5u2mQB8JhSag7QHojVWp9xqPRuFlauFClpmnPXkqgSHJD3BkVIXHIcCakJmQG/bkhd2lZpS9sqbbmWfI3SPqXx9vJmbJuxmdtM7TmV49eOE+gTyO+Hf+fjrR8D0Gt+L7aP3I63l3RWE6Koc/S+/HHge6XUTqAl8DZGsO+tlDoE9DK9B1gCHAGigS+AMQ4e2+38fIxf33++2ZzHmkXLnot76PhjR3rM65E5b/7g+ZnTZf3K5hjAa5atSaXASjzc7GE+7fFp5vyW37Z0WXlF4dp1YRep6anuLoZwEUfq9NFabwdsPSzoaWNdDTzqyPGKmkHNq/Lyb7vZdSrW3UXJt8TURO5afJfVPHuv0iNrRLJj1A5afNMCgGazm7Hrvl1OKacoPFprbqTewMfLhyG/DeFUnFHr2q5KOzae3QjA4tsWE1YmzK3Pb4RzOBT0PV1IoNFBq214OTeXJG9p6Wn8fvh3XvvnNQBaVGpB69DWtKvSzqFqGS/lxe+3/s6Q34YAEDk3ktV3rsZouCWKqnSdzs1zbuZa8rUc18kI+AADfx2YOT2s/jBGNx9NldKSbbY4kq9tB7WrVd7dRchVWnoaCakJPLfmucyADzC913TGthlLl+pdHD5G7eDaLLh1AQCXEi/xv53/c3ifwrXuWXxPjgH//ib357rtvIPz6D2/NwmpCcQmFZ+7XGGQK30H1alUmiW7zqK1LnJXt3HJcfSY14OE1ITMeaGBoSy/Y7nTH7rWCq7FiqEr6D2/N9O2T+OuBndRLqDo3wGVdP/b8T/qlatHz5uMGtdms5tZLW9aoSlTIqew5dwWxq8dzzNtnuH+pvfzTMQzmeukpadx8vpJPt/5OYevHmbf5X0AtPu+HQCTbp5E7eDaNK7QuJDOSjhCgr6DbipfmtiEFC7GJVOpTNEZLP168nU6/djJat7zbZ9nZOORLjtmldJVaFS+Efsu76Pvz30Z3Xw0H2/9mJaVWjKm5RhaVGpBoG+gy44vzNJ1Op/v/Jxp26dlzmtbJXunuh8H/QhAtaBqDK4z2Oa+vL28CQ8OZ9LNkwCITYqlyxzzHeKLf79otf6/I/6ltG/R6LAospPqHQddSzQGR9949LKbS2Jt/+X9mdPTe00n6s4olwb8DDP7zgQgITUhs0nn9gvbGb1iNO1/aE98SrzLy+Dp4pLjaPFNC6uAD7DprNHRrkrpKjSv1Jwdo3bYtf9g/2Am3TyJaqWrUa10tWzLO/zQgVm7Z7m0BdCK4yu4Z/E9jP97PB9u+dBlxymJ5ErfQb0bhzI96jBpRayD1oPLHwTgz2F/UimwUqEdt4xfGV5s9yKTNk6yuTxybiRzB8/lwWUP0qxiMz7p8UmRqxYrrhJTE1l6dCmv/vNq5rz+tfqjtWbZsWUA9KnZh/e7ve/w73xQ7UEMqj0o8/2HWz4k+mo0a2LWZL7/cMuHfNHnCzpU7eDQsTJorekyp4vVs4idF3cCRkuyVpVbOeU4JZ0EfQdVKWt0ynrix23c0iL7VU9hi0uO42js0cz3hRnwM4xoNIK2VdpStXRVgvyM3CFaa5p/05zEtMTMlj5RMVE0/6Y5IxuP5Pm2zxd6OUuatt9bV99YNp+d0GkCAd4BLutA93SbpwGjWimjCS/Af/74Dz7Kh22jzIkJU9NTSUlPoZRPKeJT4nOtCjoVd4p+P/ejtn9tHv8me36oUj6lSEhNYNTSUdQNqcv3A75Ho6V6KRcS9B1ULaSUu4uQ6esLX/P4j+Z/DHc2qatXrp7Ve6UUy+5YRr+f+2Vb99u93/Lt3m/55+5/KONXhs+2fUaryq3oXL0zN1JuUMqnVI5XpgnpCTbng1H3XNavbIm8k4hPiafDD+Yr6DahbayWL7ptkdX7wgqCXsqL7SO3E5cSl1nvn6pTsz1AttSwfMPM6sjZ/WaTptNoWbklA38ZyJl4o9P+kaQjVtvcWf9OXun4CmB+OB19NZr2P7TPtv/3u71P3/C+jp9cCWF3ls3CEBERoTdvtr+3a2ElaQoftxiAiJrlmP/fTnms7Rp7L+1l+KLhVvOKemoErTXf7fuOdze9m+e6E7tM5JY6t2S+//fMv3y+43M2nzP+PuYNnkdoYChl/cry+J+P8/epvwFoUqEJcwbZn6m1KPpq+Vd8cPaDHJe/0uEV7mxwZyGWKGeLjyxm3N/j8l4xn8a1G8c9je6xmpeUlsSamDWMjRqbw1bk2Wnw5LWTDPh1AAAbRmzId4ODL3d9SfTVaCbfPDnPdTPuWoJ8g1hz1xp8vXxJSE3gpbUv0S2sG0PqDrFa35H4pZTKMcumBH0neHrudn7dZvRiPDZ5YB5rO989S+5h5wWjbrNyqcp82fdLqgVVw9+76LQmyk1G1U9+3F7vdn459EuB9v9S+5e4q+Fdea9YDFxOvEy3ud2yzffx8mFArQFM7DLRDaXK3ew9s3l/8/s2l4WXDefYtWM5bvtKh1cYVn8YD/78IJvjN7Pl3i34efvZXDddp/Pd3u8IDw5n5fGV/BXzF5cTzQ0sfhz4IzXL1uSPY38wa/csktKSeKXDK0zbMY29l/Zm298HkR/Qu2Zvq3laa97e8DbB/sF8vvNzq2Vf9Pki87lC1v+9hYcXMn7t+BzP01Lryq35qt9XrPlrjQT9giqsoH/6agKdJv8JwM4JfSgb4OvyY2bIWoe6c9TOYludcTT2KAsOL6BPzT5sPLuR9ze/z+e9PufcjXNWDyctlfYqTXx69hZBEzpO4HLiZT7Z9gkArSq3Ik2nZX45Dq0/lPkH5xerO4HU9FRafWt+WFmcU15orUlITSDQN5D4lHgCfQJJTU/lxPUTBPsHc/7GeeqF1MPX2/hfsvd/+WjsUW757Za8V8yBj5cP20aan0fc9vttRF+Nzte2y+9YTtXSVYm+Gs3tC24v0HGndJuC3zE/CfoFVZg5uL/fcJyXft3Nb492pmWNkEI5JsB3e7/jnU3v0KFqB7roLtzX975CO3Zh2n5+O29veDuzY9DqO1dTsVTFbJ9xuk4nITUhsw47v/+k73V7j9439S6y1WEp6Sm0/tY8rPQbnd7gtnq3ubFEhcuR/+UjV48w5HfrqpOWlVqy/cL2zPfl/MvxVb+vqBZUjX4/96N8QHmrv5uG5Rvi5+2XedGQ4cnWT2Y2Tc4qyDeIpLQkUtJTMufN6juL1pVbs/3CdpYcWcJjrR4jxD+EVSdW8XTU05TxLcP1lOsADC8/nJcHv2zXOUvQLwSrD5znga82UadSaVY9UzjHXH5sOc/+9axx/DtXs3vDbhloIotrydd4/q/nWXd6Xb72d0udW3iy9ZNUDqzspBI6x8R/JzLngHFHMilsEoN6Dspji5LF0f/l5LRkktOSeWr1U2g0M/vOZMXxFYyNGpvj3fGamDU8uip7jshnI55leIPhpOt0An0DiU2KJdg/mLWn1pKSlkI66Ty1+qls2226ZxMBPnmnYLd86G3v3VxuQV9a7zhJ13pG08jDFwqv81FGwAeoWKpioR23OCnrV5b/9bbOBZSWnsalxEtUDqzMtvPbGLV0VOayBYcXsODwAl7t+CrD6g8r7OLmKOPqb+vIraxbk78vMGHm5+2Hn7cfX/b9MnNe75q9cw2qXcO68mbnN0lOS+bNf98E4PPen9OpmnVjjWD/YACrPFZTe07N/ML49ZZfqVuubr7Luu7udXyy9RMibtiM2Q6ToO8k3l6FV48efSWa2xaYb+0/7/15LmuLrLy9vDOv5FtVbpX5j7/9/HYeWfkI8SnxvLH+Dc7EneGJ1k+4s6gA/Bb9G4uPLKZ15db4ehXe8yIBt9a9FTDuANN0Wr6bvnYN68ofd/zBteRrBQr4YFyovNzhZaKiogpY2vyRNAxOlDGoysnLN1x2jHPx56wC/tq71ma78hD2aVm5JevuWkfTCk0B+GLXFzSb3YxNZzdx39L76DWvFyeuneDr3V/zx7E/AOPOoMdPPbiSeMXu46br9GzZKj/Y/AHNZjej2exmvLLOaI9eN6RgwUM4T4BPQIH7OlQNqkqD8g1cVCL7yZW+E5UP9OPstUTmbYlhbO/6LjnGR1s/ypxuXbl15q2lcA5vL29+HPQji44sykwklpHSAqzzyvOXebLrXGP84Cqlq7Dk9iW5XpHfSLnB6+tfZ8nRJYC5VykYHYksq+0svdzBvod6QliSoO9Efz7bjcavLueTVYdcEvT3XNrDoiNGT8t1d6+jrF9Zpx9DGAbVHkRsUiyTN5o73fgoH1J17knEzsaftWplk+GnQT/RqEIjElMTs/UatUx9bRnwn4t4jiF1h5CankqFUhXsPRUhrEjQd6JAP/Ov87dtp7i1VXWn7v+uReYORhLwXe+eRvdwT6N7OH/jPOUCyuHr5ZvZpjw5PZm1p9ZyLekat9W7jYsJF/nj2B85Jpq7c9GdzOo7y6p36piWYzIzYS65bUlmj1CAHaN2yNCEwiUk6DvZ7Afbcd+sjWw4etmpQf9IrDn3SHHulFMcWTbfzKjX9ff2zxyYBIzWUyMajaC0b2leXvcytYJrcTT2KPXK1ePQlUOAdTXR5ns34+/tz/AGw0lMTaRaUDV2jtrJjgs7aFGpRbHtYCeKPgn6TtaxdsZtuPP6P2itMzNTftHnC6ftVzjfkLpDsuVQSUxNtMqAOXfQ3Mxu+uUDzMNtKqVoWblloZRTeC4J+k6W0YIHnHel9t2+7zKnnZWbXBSeAJ8Adt23i01nN9G0YlNK+RSdzKzC80iloYv8uPGEU/ajtc7MQvlU66ecsk/hHm2rtJWAL9xOrvSLOMvsk56Ua8WdMlJlA/wyphO3T/sHgCnDWtCnSSg+Xl6U8rPO0XMpLok2b63ksxGtGNTc/YPpCJETCfou8Gyf+rz/x0ESktOyBYeCOBV3KnN60W2LrOp/PdWumFjG/bKT2Q+2Y130Rf48lIwz0w0duRBn9T4j4AM8M28HzMt9+8d+2EZkg8oE+cu/liia5C/TheZvjWFkh5p2b7/w8EIAZvSeQc2y9u+nJBn82VoAIt5amTnv+tebmHV/25w2yVPMlRs8PHsz+89ed7h8AE1fW86+N/rx18HztA0vT4Wg4jGugSh83284zu5Tsaw/fIljl4ye/M/0rs/jPevlsaX9JOi7QAdTC57D5+PyWDNnWmumbp9q7E8e3jI96jDvLNtvc9mf+89zMS6JirkE1zUHLzBq1kZCy/qzYXwvwLoaJ6ujkwaw+9S1zC+Zrx9oS1xSKo/9sM3m+mMi6wAwLeowAI1eXZa57PDbAwo1N5PIn4zP/+ikAW5pIvvw7M2s3Hcu2/wpKw4S2cB1WV4l6LtAm5rlAPj6n2NMuKWJXfuYtmNa5rQnt9nWWlPrxSVW8354uD0jvtxgNS/irZU2Ry27npjC0l1nef5nIw/6uWtJuQb721tXZ9LtzVBK0SwsmGOTB5Kalo6Pt9HmoVejUBq+YgT0Y5MHciM5NbNT3o3kVGKuJLBgx2mrfXZ550/eH9aCZmHBpKSm08Z0l7Ll5V5yF+Amk5buy5yu9eKSQh/xLvK91ZlX9rYM/mwtwxv4ObXqMoMEfRewDNJJqWn4+xSsXn/h4YX8b4eRDnjOwOIxqpOrfL7GekDsDrXL06luRaIn9sdLKZb/GcV/Vxr/POHjFhM9sT87T8XSvHow/x65zL0zN9jabTYfDm/Bba3CbC7LCPgAAb7eRE/sn3nlbtkLO9DPh0/ubsU97W9i+Ix/M+efiU3kHtOXVKUy5iCfEfzfHNKEkR3D81VOYe3HjSdodVMIDasYPdTT0zWz1h2lWfVg2tc2p644fTWBkTM38P3DHbhj+j+cuppgtZ9/Dl+kY+0KVv+7aemaKzeS+WHDCQa3qMaJyzdoF17eoed0aemagZ/8nRnwywb4MPmO5vRuHIqv6e8s46Jk7oFk3rH7SDmToO9ip68mUqti/rPznbx20moszSYV7btTKAkOnL3O5KVGlU7zsGAWPGbOV54RiEv5KAY0q8KSXWcBqPvS0hz3N2VYC8qX9mPFvnP8sOEEN9eryBejIjh6MZ5GVfOf1sLyS8CW9rUrMP+RjgQF+DDok7Wkpps76l24npRt/Vd+38MdbcKsvkBKgvu/2kjUgQv0ahTKl/c5Lzf81RvJ/LT5JJfjU/jfX0Z1WsaV+shZG1gXfQmA/W/2I8DXm7R0nTmcaYdJq6z2NWVYC56Zt4MRXxhfyhmDIE1dHc17yw9krvfBioOZ0/ve6Gd34K8z3nzXWrdyECvHZh/v+I0hTXj19z283invAVfsUbL+yoqQ94e14Nl5O/hlawzP9Ml/elXL/CtFaRCPwnYmNoG+H60BICTQ1yrgZ/XZ3a2pvWtJjssD/byZcEsT7mhjXMl3b1iZt28zj05UkICfXxHhRkur3a/3zawOyk3jV5cXehWDK12OTybqwAUAVu47x+mrCXh7KRbuOM1DXWpxMS6Z8qX98vWsIz1dczUxHTA3jc2qyzt/EnPF+uq94SvLGNe/IaV8bQfoh7rUYmDzqkarLJPDF+Jzrf4D6+c1GepWDuLdoc35et0xklLT+Hyk9Zdcerrmm/XHrObZCvgAozqGM6pjuMvy6Tsc9JVS3sBm4JTWepBSqhYwB6gAbAFGaq2TlVL+wDdAG+ASMFxrfczR4xdV3RsYI2kdPGdfixBPTbh1OT6Z37ef4vWFezPn/f5o51y38fJSHH57AFuOX+HOz9cT6OdN46plaRYWzKiO4QW603K2AF9vDr89gGsJKbSduJLUdM3WV3oTUsqXLSeuUMrXm0GfGg+L95+9lllNUdy1fnOF1fuMK22Atxab69PfurUpfRqHMuWPgzzZqx7VQqw7r+0/e41+H/0NwIrLW1m884zN42UN+Bky7hQBvhgVwad/HmJnTCx73+ibeWf1+6OdGTI199HI+jYJZfme7A9dM0Sfj7Nq3hs+bjFh5Uqx9oUePPPTDn7eGpO57Kle9RgWUSPX47mSM670nwT2ARl/re8AH2qt5yil/gc8BEw3vV7RWtdVSt1lWm+4E45fJJUv7QeQ6x9KVhlNNNuEtvHIgJ+WrrMFi89GtKJmhbyDtreXol2t8kXyatnbS1GutB/Rbw+wmt823LrfRb+P/rYKRnm5cD2J4FK+Fqk/ip4tL/eyeWWe4eXfdvPyb7sBmLv5JB8Nb0n3hpUJLuXLhetJmQEfsAr4U0e0JriUL53rVrB60P9c3wb8t1sd3vvjANNNLaky9G4cSu/GodnK0KJGCMcmD+TwhTh6TjEGSfjfva2pUymIOpWCUMp4Tnc2NjFb9VBuYq4k2LxreKqXa8bayC+Hgr5SKgwYCEwExirjKUgPYIRpldnABIygP8Q0DTAf+EwppXRRHpndAUopmlUPZvfp2LxXNslIy/tx949dVawirXGW2+ZejUI9onfrrgl9aDbBGIlr6PT1LHny5jy3uRKfTNuJ5mD6RI+6jC1ANaKr7Dtzjf4fmwN1hSB/xg9oyNtLbDe3zeqpudsB4wFn87AQm+vMe6Sj1RfmsckDSTc9N/EyVRe90K8hvl6KT/6MBiDSdOedmzqVgtj9el+uxCdTo3xgtuVVggM4NnkgUQfOUz+0DNVCSqG15kJcEj5eXvy27RTXE1P5cOVBG3s37HujX57lcDXlSMxVSs0HJgFlgGeB+4F/tdZ1TctrAEu11k2VUruBflrrGNOyw0B7rfXFLPscDYwGCA0NbTNnjv2tV+Li4ggKCrJ7e0f9d2U8CakwvVcgpXxyr7vUWvPECWM81k9rfmrX8dx9vo54bFU8cSnG9IBavtxa1xc/77zre4vzOVuavj2RDWfTCA1UdK7uwy+HUuh1kw/3Ns7epDMuLo7H1mb/3czsE+iW/gBXk9LZeCaNH/YnW81vE+rN462Mh5HXkjXLjqZwez1fridr0jUkp8GLa41qmXZVvNl4Ns3m/qf2DOSf4/F8H22c29f98lddl5aumX8ohfZVvAkr44VPIf5ubqRoxqwyN8kM8IYG5b15uk3+H8468rfdvXv3LVprm0/P7b7SV0oNAs5rrbcopSLt3U9WWusZwAyAiIgIHelAQ9WoqCgc2d5R/mv+ICE1hdiydejf7qZc152xcwacgDEtxhDZMtKu47n7fO118vIN4patBqBnw8pMK0Dv2uJ6zll162a0MDkTm8gvh4xvv5UnUll5wjxS18DmVZk6ojU/L/0TyF6HXaVha5pUK9zhM/efvcb9FlUwGSoG+TH/qV5WTSBvsbF9n8gkpkUd5rm+DXhv+QGGtglj/eFLvLHIeKbToXZ5BvbuSOmoKJKDy/F/3WpTt3KZfJevZ48Cn5LTHO2l0RrStMbHSxW4v42r/rYdqQzsDNyilDqG8eC2B/AxEKKUyvgyCQMyEsicAmoAmJYHYzzQLbE+uLMFADeSbV/BWPp0m3F174lJ1UZatKWf6UA6heJMKcWZ2MRc11m88wzh4xbzzF9GwP/PzbU4Nnkg7w41kvIN/GSty8tpKWude4aDb/Vn88u98xXkKgT588qgxgT4evPKoMY0qlqW+zuFZy6fMcp8sfresBYFCvjuppTCy0vh6+1VpDpY2h30tdYvaq3DtNbhwF3An1rre4DVwFDTavcBv5umF5jeY1r+Z0mtz8/QpqZR75iex2leuHEhc7pK6SouLVNRlNFRZesrvd1cEvfa/HKvzOn+TfP+O+jVyHgoeadFS5DC/JeyfKYAsOjxLhybPNDhB8teXopjkwdydNIAygbkPMC8sI8r2um/AMxRSr0FbANmmubPBL5VSkUDlzG+KEq0sgE++Pt4cTaPK7ge84x70P9r/n+FUawi5c/95tZNGS2ePFXFIH+r1kfLdp/lke+2sHJsN7TW9P5wTeaym8oHZmv9A0ZKgZVju7r8irjbe6szp0d2qMmbtzZ1+jGK0tVxSeKUoK+1jgKiTNNHgHY21kkEPKq3kVKKWhVLc+RifI7rHLlqTjMwuvnowihWkbErJpYHv94MGGkQhLV+TatYfQlkTNuq6130eJfM9v69PjB/OURP7J9nD+KCiktK5bjp7uyXMZ1ofVM5p+5fuFbRbeBbQtQLLcOBXFL2PrD8AQB61+yNn7dnXemOnGWuy7+1pfMGkfdETasHc3O9itnm28ri6PCxXlsOGD2lJeAXPxL0XaxhlTKcuppA7I0Um8svJ14G4P1u7xdmsYqEKmWN5mvuSm1b0nz7UHuOThrAnRHmxHGPfLeVhHw0JMiPqaujrTobfX5vG6fsVxQuCfou1qSa0VF575lrua7naT1wf90WkzloiQR851FK8e7QFlbVQj2nRDm837E/bbdKQFa7YunM/EKieJGEay7W2BT0X1+4h2VPdbVadv7GeQDubXRvoZfLna4lpvD0XCPJVd8m2bvFC+eoGhzAmdhETufRkCAnV+KTuf+rjeyIyd6rfMXYbjIwTDElQd/FKpcxqjD2n72O1tq4qtUarhzl1uVGsO9Ww3a2vZImLV1bpZYFsmUjFM6z/sWemdUxp68mEJ+Uyv/+OsL7w4x2/VnvsO6asZ5/jxjVjevG9aCzRZK0DLtf78uN5FQJ+MWYBP1C0EQd47wO5uGv/mHmkCqweiLs/pky9ZpwHc8YDvHoxXi6vx9lPW/SANsrC6cZ2Kwqi3edscpymZHxcVDzqrw4oBEXridxa5Ysk1kD/gOdw3ltsDG2gwz6XrzJp+dqy8az2N8Y65YTgCmtziUvL06nXmdQ7UFuK5qrJaWm8epve5i7+WS2ZeEVAqUuvxBMubMFi3fZTke8aOcZFmVJVdyzYWVW7T+f+X5omzAm3ta0wKO/iaJLgr6r/TvV5uxlQUYWv/Ak++pbi4MGL2cfbGLaPa0Z0KyqG0rjmQJ8valSNoCz1/L3d/b5yDacu56UeaX//jDpP1HSSNAvZOcrdaTyhfV8HWw84B257ivY9DO8eMLNJXOtv5/vzo8bT9CzUWV3F8XjRD0XyeBP1zL7wXaZg5RorRk1ayNeSvHXwQssfqILoWUD8PH2onpIKWlGW4JJ0HclizwoW4b8yas/rWfPyVocm7yMs7ON4foCtYakWJgQDK9dhRLyj7bl+JXM6Yzmg8/3a+iu4ni0AF9vVmQZmk8pxbcPtc9xGwn4JZdnNQ4vbMdNw6c1HESL5q3Yo2sB8Mf+owDcej3Oev3zeykp3jSlxrXMmCiEcD8J+q4yIRi+NrVO8Quyyn/yzAYjs3jnAVNh5G/mbaZ3govRsOpNY/v10wqxwM61/eRVAF4b3Ni9BRFCWJGg7wrJWRKstTeyZx6dNADvIPPVfJewrlCnO7xsbi3BZ23gb1NKhuUvurqkTpealm7VVV+qCYQoWiTou8KURubpvpOgemvACICBNb4BoFu1/gT5mYZC8/GHtg/b3teFA7bnF1F1X1qaOX13uxq5rCmEcAcJ+q6QZOq2PmYDdByTOXvtKfPIRotWWadkoP97UNc8iAZlTUmzFo11VSmdZvepWLafvMo/0ebhjl8b3JhJtzd3Y6mEELZI6x1nS083T1e2bq3y35X/NVZJDQKyVHt4ecG9P5vf37gM79aC42th00xo+5CLCuyYLcevcMf0f6zmNQgtwwOda7mpREKI3MiVvjNpDW/Yzi+ekGoeyHp45S9QyshFk6NAiwyGi8fCT/flvK6bnL+emC3gAyx76mY3lEYIkR9ype9Mr4eYp6tbJxKbu38uABGhEVQMCEBreG7eDuqGBjEmsq7t/b1wHN6paUzv/c355S0grTXp2rhHmbPpJON/3ZW5bOFjXWhavaw8uBWiiJOg7yzXLUYoqtoSHlphtXjKlimAMVjKzhPGoBa/bDsFGPlNMrJxWikVYnTYyvgyWTYe+r3t3HLnk9aaWi8usbnMFUPyCSFcQ/5TnWH/EphS35iu2gL+7y+jjt7kdNzpzOkKpSrQvYF1KoJV+86TI6VggKkJ579Tjfb7ls8NCsnrC3PuOCYBX4jiQ670HXXtNMy52/x+6FfZVnlw+YPZ5v3wn/aM+MIYI/bFX3YxtE0YvjkFzwYDYMmz5vdvlIMJ2Qe2cBWtNV//cwyA3x/tTI3ygSSmpOHr7UWlMv6FVg4hhOPkEs1Rlw5bvy9fO9sqp+KMapx1d5tzlneqU5FNL5mbaL69ZB+paem2H+4GV4ex+6GGRd79Jc85Vu4CyEjNqxS0qBFC+dJ+VAspJQFfiGJIgr6jZpvy4d82A8buy5Yw7cBlo3NVkwpNKOtX1mpZpTL+dKlbEYCv1h2j7ktLqTN+Cem2An/ZqvDQcqjd3Xi/cQb8+z/nnosNczae4LEftgGwYXxPlx9PCOFaEvQdEb3KPF2/L5Stlm2VoQuHAvBIi0ds7uK7h9tzU/lAq3m1x9t+YArA3XPM08tegLTU/Je3gPaduca4X4wWOn4+XrYfNgshihUJ+o7Y9p15ulRItsXaIrVyt7Ccx8H9T9fsVUJHL8bbWBPwDYCXLFoKLXg8z2IWRFq65tHvtzJ1dTT9P/47c/7Bt/o79ThCCPeQoO+IPb8Yry+dtbm4+TdGGgJfL99c26+P7FCTbx5sx2+PdsbHNOB01vFkrfgGwCumlAc7fjBa9KQmF7j4GbTWrNh7jm0nrlBn/BIW7zrDe8vNOX8y8uELIYo/ab3jDL6lcl28+LbFuS4H6Fq/EgDbX+tD09eWAxA+bnHOIxh5+1q/f6sSvp2+zV95TaLPx9Hrg79yXefI2zJ4uRAliVzp2yv+kvFataXNxXHJxgAp9zS6h6pB+R8TNsjfh4pB5lYxExfvs1quteZaYgrR56+T+MIZKFs9c1nnf0YaqSDSUiEl9zFRT11NyDHgf/tQO/4bWYfoif3x8pIetkKUJHKlb68T643Xrs/aXLz3ktGZqU5InQLvevPLvej/8d/sO3ONL9cepUfDynSqW5GUtHTqWaQuBjg2eS/89R6sfsuYETUZ/ppsTL983kjbbCE2IYUWr/+R7ZjvDW2Or7cXg1tUw9tLcXO9SgUutxCi6JOgb6/l443XwAo2F0/dPhWAeiH17Nr90idvzhyMZMSXG3Jcz1inMWHqI9b6P2UO+AD/fJr5pfTIt1tYtudsZhPRDPvf7Ed8UioVgqTNvRCeQIK+PS4dhqvHjembOtpcZev5rQC0rNzS7sPMGd2Bu2b8m691Y3RlNqfXJ8LroHnmn29Se0kD0i1q8dZa5LyfMLgxAb7eBPh6211GIUTxYnedvlKqhlJqtVJqr1Jqj1LqSdP88kqpFUqpQ6bXcqb5Sin1iVIqWim1UynV2lknUej2LTRP23jIejHhYrZ59uhQuwIbx/ekeoj5QfEzveuz5/W+7H+zHy1qhFitPzT5NWJ1IK+k3J85r5E6gSJ7rp5jkwdyv+S8F8LjOHKlnwo8o7XeqpQqA2xRSq0A7gdWaa0nK6XGAeOAF4D+QD3TT3tguum1+MloOTN2v83F3X8yes0+0eoJhw9VuWwA68b1sLns90c7W72Pioqi97pvOZ+UxCEdxhy/t1jsb1RD7e86jbU+7alX0Z/29bJ3IhNCeAa7g77W+gxwxjR9XSm1D6gODAEiTavNBqIwgv4Q4Btt9Fj6VykVopSqatpP8XLtNHj7QZkq2RadjTe32f9P8/8UZqkAI1XC5GX7mfNXnNX8hmvGYDWO19h9NnsQCyFKNmXZa9TunSgVDqwBmgIntNYhpvkKuKK1DlFKLQIma63XmpatAl7QWm/Osq/RwGiA0NDQNnPmzMFecXFxBAUF2b19TlptfZ7S8SdZe/OP2ZbNvDCT7Te2c3f5u+lUppPTj52brOebdPQf0q+dov+V77Ktm+odaLP8xY2rPuOiTM7ZMzhyzt27d9+itY6wuVBr7dAPEARsAW43vb+aZfkV0+sioIvF/FVARG77btOmjXbE6tWrHdo+R6+VNX6yiE+O102/bqr7zu+r09PTXXPsXOR4vunpWl8+qvXWb7U+uMJc/p9HF2bxXMJln3FRkZam9YWDxmdoUuLP2QY554IBNusc4qpDnbOUUr7Az8D3WmtTTgLOKaWqmpZXBTJGCDkF1LDYPMw0r3hJvJbjoqVHjTb09za6t2gNG6gUlAuHVvdCvV7mpG0758Ce39xZMpGb1GRj7ITPIozR09ZPc3eJRAngSOsdBcwE9mmtP7BYtADIGMX7PuB3i/mjTK14OgCxujjW52+eabx2fT7bor9PGQnKhtYfWpglKrgGFsnT5t2XfUwAUTTMiLR+v/xFmNIIrzT78ywJ4ciVfmdgJNBDKbXd9DMAmAz0VkodAnqZ3gMsAY4A0cAXwBgHju0+KycYr7WzZ83ceGYjDcs3JMCnGKQgfu2qefrT1rDD/mcnwkmunjSS52X8nN9jzLccje36abr+PQz2/m57H1mlpcLGL2Dl63D1BJzd7fxyi2LFkdY7a4Gc6jCyjbZhqmd61N7juV1aCvxtcUNTrZXV4oNXDnI95TqtKxeT7gdKwfgz8LYpL9DS56HFXe4tkydLT4OPmtpe1vR242dCsHne0nHQeIj5fcIV2P4DdDT9iyXfMMZb2PqNeZ21Fn+/o36H2pFOK74oPiThWn592Qui3ja/9ytttXjtqbUA3FH/jsIslWP8AmGMKcVDYqxU87jTG+Vtz3/K4sp8QizcahotLeuwnO+EG6lBlr9kJN17u6p1wM/qmyHyeXsoCfr5dWa7edqyasRk4eGF1Cxbk/rl6hdakZyickNo/19j+tPWELPFveXxRF9YdL4b9KER3DN+QmpYr9vybi5WaAvH18JpYxhLEq6Yl6//zHjoa6nhIAgw3SX4WTQB/LSY3JUKp5LcOwX18KpsqRcuJVwi+mo0t9S5xU2FclDvN2DDdGP6yx7w/FEIzOHKUzjXyU1wyvRFO2Ie1O+T5yanq/Wj4qVNsH4qhHeBhU/aXvGxLVCxrjGdmgw63RiAJz3NfGexYy7U7QmlK9rehyhx5Eo/P9JSjNe2D0NY9v4OOy7sAKDnTcV04HAfP2jzgPn9u7WMwCBca8ETMLOX+X0+Aj7A5QoRULUF7JpnHfCfOWC+om/7H3PAB+Mz9jU1MPDyhooNjOlfR8N7dYzqPeE86Wmw7Xu4fi7vdbNKTba+e3MyCfr5sfZD49Xbz+bijPb5bULbFFaJnG/wR+YhGAG+HuS2oniEc3th62zz+3EnC7Z9nSwXGGP3GWlBnj8GL8bAwPdz3/6uH6zfT77JeB4Qd972+qJgPmsLv4+BKfUL/r80tR28E071mIV5r2sHCfr5sXqi8XrzMzYX/3niTyoEVCDYP9jm8mLD2xf+s9qYPvEPHP/HveUpaa6fMzfHnG6RkntCLASULdi+IscZr7fNMLbPyKPk5QX+ZfLevmJdY7tXL5vnrf8M3q8H5/cZZXy3jnwJ2OPiIbhs8ZD82N+mJri2EzRaSY6HK0cBqBf9pUuKJ0E/L+9ZDIJio94z+ko0yenJ1Ctn32ApRU711tDB1IXiq/7mYSGF46bYeMj/bLR9+/LxN4J2i+GOlcnLGx5Zaz1vWgfj9cZF40tg98+OHaMkSLxm9I3I+NJe8RrMHZl9vZQEowc1QPkso+ZNa280pbXl3F6jP8Xbpi/vXhP4u4tr+s7Ig9zcXDkO8aYrndtm2FxlzgHjgxlYe2Bhlcr1+r5tpGe4fhp+vAseXuHuEhV/x9ebp1vcDbdOtzkWg1tUaWZ8gVw7DR80Ms+v1BAu7If5D0LTQmiKvPkr4y6ldiQor6LTmODUFusWVgDrPjJeJwTD4E/g+hmImmS9zshfISgUvrsdjq8z5r1dFV69YqTXyE3bh0lb75qWdBL0c/Nxc+PVJyDHK6q5B+YCMKTOEJvLiyWl4BnTLX7MRqN+8rFN7i5V8fZVP+O1y1jo9Zp7y5KTstVgzL/GlX6/ydDhv+YOYROC4ek9ULoynN0FV49ZfxHcuGzcfRz9Gxr0s96v1kY1kZcP/PaI8ZDTyweOroEnd0CZUONKetFT2ct07y9Qu7urzjhvc+6B/YtyX2ehjXEzxp8x+sEAPLDEuEv4aZTxPq+A/9I580N3F5Cgnx8v2s4Lt/Cw8aCleaXmRSvBmrM8vAq+7AkXDxrVPKVtjwcsbIjZYgTRfQth6XPm+UU14Geo3Mi46s9w6/+MQA3wYRPrdec/CP3fNc7x2N+29xdcA2JzeUhtq8rL0ne3AxDU5kO4UBUqNcjjBLLQ2nhgXqYq1O9rvSwlAXb+BLExxvOre+bBuT1QqpzxzCMlIXvAH/gBRDwIacnGl9+XNlrslQs3B/wMjYdAt3HWY1hb6jAGOj9lfAG6mAT9nFhmNPTO/ms6cPkA49cao1I90OSBbMtLhLAIeGCZcZX6Xm2XX4GUCGd2wuc32172wLLCLYsztLzbSDkyLYdB7pZmTzxoJbeAb8urV4yH0Wkp8O1tmV8mEVueNhK4j/gpe/DOSVoqzOwNp7ea5w3/zqg+mhSWff2MlCQAdXtB9Erz++eOWF/0+Pgb/x8t7zU6yY1ebczLTaRF0O/wKPR7O/f1XUQe5OYkxlSdcZftgUbe2/QeAEG+QfSq2cvmOiXCTR3M08teMG7jnTDwTomUmpRzwJ8QCzU72l5W1FVuCC+fN/fleGq3MVSor8XVbJXmxoPLF2OMZcoL/CxaEXV8DFDGF1+/d4x5T++BJrcbV87la8NLZ42AD0ZLsvsXwWNWYyzBD3fm7+8vMRZmdLMO+ABz77Ud8LOyDPjjz+R8l3vrVBjzT94BH4xq04ye1m4K+CBX+jk7sd64JWw4INuis/Fn2XDWyFmzfsT6bMtLFKWMZn1vlIctXxs/GWp1gxFzjVtj02041Vqb/9F8S8NLpwu7xK6VHA9TGhlVATXaGc3zKtaDfz6FFa+Y1+v+Mqx+y5jOYSzlYsXH3+jLMfgj87zxp40qkKxVGf5l4DUbnYv6mpo+1+wIHUxVRsO+yr6epYr14ME/iFn+MWGnFhvzXg8xrvhrR1oH213z4eeHoOtzsOY9i7KXgv7vZK97H/49XDwAq96AWz6DBY/BsNlGq6XFzxjpSfq8ZfNOvzgrWWfjLGmpxtP4HIxcajTVerJ1Dt3fSxovb2P8gDXvWs8/+hdMzDJOsOWVVUo8XDhQ8HrYomzR05AUC7Ny6T2bUUXR9Vmj7jc/V4HFkVLZA74r3NSe6HqjCWvYxgjQYFzxlw2DiPvhz7es17cM+C3uhiHTjM+jzX1Gk8kDS4x697AIYJC5/01riyaYbR924Qm5lwR9W9403cpVaZ5t0dHYo5yNP0uj8o14uFnJ/cPIpuOjcG439JpgdBGPXmn9zwVGPej1cxDaBKo0hT9ehp8fhkdyeMhX1N24bNzF1Olu1DHPvTfnB5YZxu43V1EoVXIDvjvc/IwRxDOalV6LyR7wM1RvY6SPztpRzS8QmhXxQY5cTIJ+VqkWoxJl9Ho0OXHtBLf8ZiRVm9hlYmGWyv1KhcDdFs83buoA7R+BpS9Atxegko1WGH+8DGd3wqx+xsO3Y+sg2qLNf5sHjLreiwfAt5Tx4C6rW6dDyxHGFVr0Cusc8hm0hstHYM8vEH4zbJ5l1DfXuhnO7IB1H8N9C6FGe4i/ACgIrp77+aYmGTmIbAmsAAPeNwaeuf1z2L/EqOqpWEI66BVlZasZdeKnt1mPLOZXBkJugv+uKzr9H4ooCfpZ/Wm6fbypIzS07nD1xa4vAKgeVL3k9MB1ROmKMHRmzstHRxn/mCfWGz9ZbfnK+MnNb/81tR6xeKBcuhLcMx9mdKNm+AjY/rgxKpSt/WeYPTj78rt+gDkjoGpL+L+/jHnRK+G7XDoi3f4FNL/TmG5qeo7R6p7cz0E4X7VWxoPgf6dDu9FQrqa7S1RsSNDPKinOeB35q9VsrTW/Rf9GKZ9SLL19qRsKVgxVawV9Jxlju2ao1AhqdYVNX4K2kcnzmYPZ225bBnwwrtZnGMNV1jqWJXFYQcwZYbye2Q4ftzTq36/Z6JPR+w3jLqFSQ+OORxQNwWHmh8Mi3yToW9LauDqs0cGobrCw7bwxYMWtdW8tmR2xXKXjGOMnqwFZHgqnJhvpf8HcOSg1Cd6qbEw3HWrcVWz91mhlkVWXp6H5cKPuvarpWUxKolGnrpRRPRQ1CZKuG/nj595rvb0pyRUAdXsbD/Ua3SJVBaLEkaBvKWPEoSy94hYcXsBLa18C4O6GdxdyoTxERsC3mudv3TsUjGDccoTRDrtUOdYvn0fHWmWgQf/s21t2JPMLhD5vmt+/dA7+nmJ8Aez5FTaYhiHsNcH4AhGihJKgn+HcHvN0H+tbxoyAf0e9O6gVnMPDPVF4vLwzk3ElBVSGBpEF34dvAPQwPldqtIeer2Yb91iIkkh65GaY3sl4vfNbq3FJp+8whhFsW6UtEzpNcEPBhMspJQFfeAwJ+lk1No9zm67TmbbdyMHzeqfX3VUiIYRwGgn6YE6udpN1bpRX1hnd6u+odwc1ytTIupUQQhQ7EvQBdpia/d1uHiglLjmOBYcXAPBqx1fdUSohhHA6z36Qe2wtfG3qgFWvj9GjD7iUcInInyIBeLfru3gp+W4UQpQMnhvNYk+ZAz7AbZ8DcDXxambAv7n6zfQL72djYyGEKJ48M+ivfB0+bGxM1+lh5AAPLM+6U+u4ea45H/q0XtOkI5YQokTxrOodrWHJs0YKAIDACpy+bSp/H11COum8vcEY2ODRlo8ysrGNke6FEKKYK9lB3zK3S1oqfHtrZmrcq+VuYs8d03jkZ+uh1+5rfB+PtHikEAsphBCFp+QG/WXjifx3KvyFkUvFIqXvnsHvc9fuT2CVkROmUflG7Lu8j/ub3M+jLR91U4GFEML1Cj3oK6X6AR8D3sCXWuschod3wLk98O9U8/voFaQBhyvWZmmrIXy5+xMA3uj0BnVD6tKsUjOnF0EIIYqiQg36SilvYCrQG4gBNimlFmit9zr1QD4B0GgwnxDCqdAAdl/ay4m4k0AqRP9MjTI1uLfRvdxWz8agHUIIUYIV9pV+OyBaa30EQCk1BxgCODXoH/RK43HvC5yO3wHHoUrpKrSr0o5WlVvRJ7wP9cvZGOVJCCE8gNJaF97BlBoK9NNaP2x6PxJor7V+zGKd0cBogNDQ0DZz5swp8HEupFxgwdUF+KT5MKDCACr5VnLOCRRxcXFxBAUFubsYhUrO2TPIORdM9+7dt2itI2wtK3IPcrXWM4AZABEREToyMtKu/QxjGFFRUdi7fXHkaecLcs6eQs7ZeQq7c9YpwDJzWZhpnhBCiEJQ2EF/E1BPKVVLKeUH3AUsKOQyCCGExyrU6h2tdapS6jFgOUaTzVla6z15bCaEEMJJCr1OX2u9BFhS2McVQgjhqQnXhBDCQ0nQF0IIDyJBXwghPIgEfSGE8CCF2iO3oJRSF4DjDuyiInDRScUpDjztfEHO2VPIORdMTa21zVQERTroO0optTmnrsglkaedL8g5ewo5Z+eR6h0hhPAgEvSFEMKDlPSgP8PdBShknna+IOfsKeScnaRE1+kLIYSwVtKv9IUQQliQoC+EEB6kRAZ9pVQ/pdQBpVS0Umqcu8vjCkqpGkqp1UqpvUqpPUqpJ03zyyulViilDpley7m7rM6klPJWSm1TSi0yva+llNpg+qznmlJ2lyhKqRCl1Hyl1H6l1D6lVMeS/DkrpZ42/U3vVkr9qJQKKImfs1JqllLqvFJqt8U8m5+rMnxiOv+dSqnW9h63xAV9i8HX+wONgbuVUo3dWyqXSAWe0Vo3BjoAj5rOcxywSmtdD1hlel+SPAnss3j/DvCh1roucAV4yC2lcq2PgWVa64ZAC4zzL5Gfs1KqOvAEEKG1boqRgv0uSubn/DXQL8u8nD7X/kA9089oYLq9By1xQR+Lwde11slAxuDrJYrW+ozWeqtp+jpGIKiOca6zTavNBm51SwFdQCkVBgwEvjS9V0APYL5plRJ1vgBKqWCgKzATQGudrLW+Sgn+nDFSvpdSSvkAgcAZSuDnrLVeA1zOMjunz3UI8I02/AuEKKWq2nPckhj0qwMnLd7HmOaVWEqpcKAVsAEI1VqfMS06C4S6q1wu8BHwPJBuel8BuKq1TjW9L4mfdS3gAvCVqVrrS6VUaUro56y1PgW8D5zACPaxwBZK/uecIafP1WlxrSQGfY+ilAoCfgae0lpfs1ymjfa4JaJNrlJqEHBea73F3WUpZD5Aa2C61roVEE+WqpwS9jmXw7iqrQVUA0qTvQrEI7jqcy2JQd9jBl9XSvliBPzvtda/mGafy7jtM72ed1f5nKwzcItS6hhGlV0PjLruEFM1AJTMzzoGiNFabzC9n4/xJVBSP+dewFGt9QWtdQrwC8ZnX9I/5ww5fa5Oi2slMeh7xODrpvrsmcA+rfUHFosWAPeZpu8Dfi/ssrmC1vpFrXWY1joc4zP9U2t9D7AaGGparcScbwat9VngpFKqgWlWT2AvJfRzxqjW6aCUCjT9jWecb4n+nC3k9LkuAEaZWvF0AGItqoEKRmtd4n6AAcBB4DDwkrvL46Jz7IJx67cT2G76GYBRz70KOASsBMq7u6wuOPdIYJFpujawEYgG5gH+7i6fC863JbDZ9Fn/BpQryZ8z8DqwH9gNfAv4l8TPGfgR47lFCsYd3UM5fa6AwmiVeBjYhdG6ya7jShoGIYTwICWxekcIIUQOJOgLIYQHkaAvhBAeRIK+EEJ4EAn6QgjhQSToCyGEB5GgL4QQHuT/AZDGrFB2DVHdAAAAAElFTkSuQmCC\n" }, "metadata": { "needs_background": "light" } } ], "source": [ "plot(res2[0],res2[1])\n", "legend(['S','I','R']); \n", "grid() " ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ] }