{ "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": [ "# SEIR model" ], "cell_type": "markdown", "metadata": {} }, { "cell_type": "code", "execution_count": 36, "metadata": {}, "outputs": [], "source": [ "import epistochmodels" ] }, { "cell_type": "code", "execution_count": 37, "metadata": {}, "outputs": [], "source": [ "pop = 1000\n", "beta = 0.7\n", "gamma = 0.2\n", "e = 0.3\n", "model = epistochmodels.SEIR(pop, beta, gamma, e)" ] }, { "cell_type": "code", "execution_count": 38, "metadata": {}, "outputs": [], "source": [ "inits = 990,10,0,0\n", "model.initialize(*inits)" ] }, { "cell_type": "code", "execution_count": 39, "metadata": {}, "outputs": [ { "output_type": "stream", "name": "stdout", "text": [ "CPU times: user 221 ms, sys: 75.3 ms, total: 296 ms\nWall time: 299 ms\n" ] } ], "source": [ "%%time\n", "res = model.run(0,1000)" ] }, { "cell_type": "code", "execution_count": 40, "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": 41, "metadata": {}, "outputs": [ { "output_type": "display_data", "data": { "text/plain": "
", "image/svg+xml": "\n\n\n\n \n \n \n \n 2020-12-23T17:55:41.526700\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", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABAAklEQVR4nO3dd3hUVfrA8e+Zkkx6JxVI6C10QZogNhAV10WxLIuKsr9VUVdXRde1rXUtKGsDEYVdBRQrCiJSLIh0pHcCSUjvPZmZ8/tjhh5a2p0k7+d58mTm3HPvfQcm75w599xzlNYaIYQQzYPJ6ACEEEI0HEn6QgjRjEjSF0KIZkSSvhBCNCOS9IUQohmxGB3AmYSHh+v4+HijwxBCiEZl/fr12VrriOq2eXTSj4+PZ926dUaHIYQQjYpS6uDptkn3jhBCNCOS9IUQohmRpC+EEM3IWZO+UmqmUipTKbX1uLJQpdQSpdQe9+8Qd7lSSk1VSu1VSm1WSvU+bp/x7vp7lFLj6+flCCGEOJNzael/CIw4qWwysFRr3R5Y6n4OMBJo7/6ZCLwDrg8J4EmgP9APePLIB4UQQoiGc9akr7X+Ccg9qXg0MMv9eBZw7XHls7XLb0CwUioauAJYorXO1VrnAUs49YNECCFEPatpn36k1jrN/TgdiHQ/jgWSj6uX4i47XfkplFITlVLrlFLrsrKyahieEEKI6tR6nL7WWiul6mx+Zq31dGA6QN++fWt03PSCcj5efdphqiiliA32IdjXerQsPMCbqEAbJqWOlnlZTPhYzXhbTJhMqrpDCSFEo1LTpJ+hlIrWWqe5u28y3eWpQMvj6sW5y1KBYSeVr6jhuc8eXGE5/1m+97Tba7KEgLfFhM1qxmpWgKJzdACXd43iss6RRAXZah6sEEI0IHUui6gopeKBb7TW3dzPXwZytNYvKqUmA6Fa64eVUqOAe4ArcV20naq17ue+kLseODKaZwPQR2t98rWCE/Tt21fXxx25doeT1PwyisrtgOtDILOonMyiiqN1tIZKu4OyKiflVQ7KqxyUVTmwOzUOh2ZNUi4HsksASIwNIiLAm4FtwxjTJ45gX686j1kIIc6VUmq91rpvddvO2tJXSs3B1UoPV0ql4BqF8yLwiVJqAnAQuMFdfSGuhL8XKAVuA9Ba5yql/gWsddd75mwJvz5ZzCZah/mdVBp0XsfQWrM3s5jvt2fw0+4sknNLeXZnJlOW7GZw+3AubBNG15ggLogPQSnpGhJCeIZzaukbpb5a+vVl++FCpi7dw/JdmVTYnQAoBaG+XgT5WukRF8zFnVowtEMEQT7WsxxNCCFq5kwtfUn69UBrTVZRBXPXJpNZVI7WkJxXxrbUAnJKKrGYFDf3b8WYPnF0jQnCLBeJhfA4ztJS7Lm5OHJysOfm4iwtRZdX4Kwod/0uL0NXVoLDCdqJdjjB6URrJzg1OB1op+uxdjrcZdXUPbK/e/uRur4XXED4xDtrFHutunfE+VNK0SLQxr2XtD+h3OHUbErO4/Uf9jB71UFmrzpIXIgPo7pHc23PWDpHBxoUsRDNiz07m6r0DKrSDlO+eTP27BzsuTk4cvOOJnldXn72A5lMKLMZTCbXY/dvTCZXt67ZDCaFUiYwm11lp9RVKJP52GN3XWdpSb28dmnpGyQlr5S3lu9j9YEc9me5/nOjg2yM6RPH/w1ti5+3fB4LURtaa6pSD1OxaydlW7aAhqrkZEo3bsSelnasotWKJTwcS2go5tBQ1++wMCyhIZhDw7CEhWIOCcHk54fytmGyeaNsNkze3iirZ3bTSveOh9uSUsCn65PZlV7E6gO5xAb78IdesVyQEMrQDtWugyBEs1V56BDlO3ZSuX8fzvIK0BpnaSmOggIchQU4Cwpdj/PycOTnH9vRYkF5eeE/aBC+fftgbdkSS1gY3h07YrI1rWHX0r3j4RLjgkiMc40eWrojg6cXbOfN5XthOUwYnMAjIzrhZZEJUUXTprWmct8+ynfuciXvwiLs2dlUJiVRmZSEs7gY7XTiLCw8tpPFlcJMvr6YAwNdP8FBWKKiMAcFYevUEe/27bElJmLy9jbolXkWael7qIKyKsZOW8XO9CJig3144bpELpJWv2iC7Hl5ZP/nTYpXrKDq8OETtpl8ffGKj8crIQFzkKthZG3ZEt++ffFu1xaTj48RIXs86d5ppKocTr7dnMZL3+0kraCcPq1D+POA1ozuWe20RUI0Gva8PPLmzKF823ZKVq1Cl5biP3w4/kOH4tu7F+bgYEyBgdI6ryHp3mmkrGYT1/aK5YquUUxdtof3fz7AfXM3UVhuZ9yFrY0OT4hzprWmbMMGst+dhj0ri6rkZJwlJShfX/wvuojQP4/Dt3fvsx9I1Jq09BuRKoeT2z9cy6p9Obz7pz5c2iXy7DsJYQCtNfaMDCoPHKB4xQoKFnyDI9d1E76tWzdsXbsScuNYbJ07Gxxp0yTdO01IVlEFY979lUO5pUy5oSfX9pKuHuE5yjZtInvadCp276YqNdVVaLHg27Mn/sOH49OrJ769ehkbZDMgSb+JyS6uYOy0VezLKsHPy8yKhy4mIkD6PkXDq0pPp2zzZsp37KD4h6VU7NmDKSAAW+fO+A8dine7tvj06YvZ/+S5rkR9kqTfBJVW2rn+3VVsO1xIRIA33947mBYBTWussfBMJavXkDtzJuW7d59wk5NX27YEXH4ZITfdhLVFCwMjFJL0myinU/P5xlT+/unv+FjNLJg0iHYtAowOSzRRpRs2kPnaa5StW48lMhLfPn3wbt8O7w4d8Bs4UIZPehAZvdNEmUyKMX3iWH8wjzlrDjF22m989teBxIfLV2lRe7qqiry58yhevpzKQ4eoSkkBIHzSPYRNmNDk7mJtLqSl30Qs3ZHBhFmuf6s3buzJNT1iZB5/USOVKSlkv/U2BV99BU4nymp1XYTt3p2ga0djCQszOkRxFtK900x8tzWd//vfegA6RQXw9i29aRPhb3BUojHQWlO0aBH5X3xJycqV4HTiFR9P2B0TCBgxUi7ENjKS9JuR3JJK3lq+lw9WHsCp4amru3DroASjwxIeLuvNt8h+802sMTEEjr6GkBtuwBodbXRYooakT78ZCfXz4p9XdeHyLpFM/nwLzy/cyeD2EbRrIS1+cSpdWUn+V1+R/eabBI4aRczL/3bN8y6aLPnfbaL6twnjk78MwMfLzOVTfiS7uOLsO4lmQ9vtlKxeQ9ItfyL9n09gDg6mxd8flITfDMj/cBMWEeDNM6O74tTwn6V7jA5HeIjin35iz0VDOTR+PBU7dxL15BO0W7ZUunOaCeneaeJG94zlt/25fLT6EH8Z2paYYBlL3VyVbdlC0k03g92OJSqKqH8945qiOEGu+TQn0tJvBv5yURtMSvHX/62n0u40OhzRwByFhRye/ChJN4zF5OND4JUjabPga0Kuv14SfjMkSb8ZiA/344HLO/B7SgHPL9xhdDiiAWmnk0O3T6Dgyy8J+sMfaLdsKbGvvYY5QO7cbq4k6TcTf7moDZd3ieTjNYcoqbAbHY5oAJXJyRwafyvlW7cSMHIEMc8/J8leSNJvLpRyTdlQaXfy/fZ0PPn+DFF7RUuXsn/0tZSuXUvU008T+9prRockPIQk/Wakb3woAH+b9zu3fbgWp1MSf1OUN2cOKfdMwis2hoSvvyJk7A0yJYc4SpJ+MxLq58Xndw0kNtiHFbuymLnygNEhiTrkKChg72WXk/70M/gNGEDr//4XW4cORoclPIwk/Wamd6sQfn74YlqF+vLstzvYl1VsdEiiDjiKS0j+611UJScTeNVVxLz6CubgYKPDEh5Ikn4zZDIp3rrZtQj1ze/9RnmVw+CIRE05S0vJevtt9g4bRtmGDcS89CKxr7yMJSTE6NCEh5Kk30wlxgXx5NVdyCisYPgrK6hyyPj9xijrP2+SPfU/AEQ8+ABBo0cbHJHwdJL0m7FbB8YTG+zD4YJy3l2xz+hwxHmqysgkf948fPv3p+O6tYTfeafRIYlGoFZJXyn1N6XUNqXUVqXUHKWUTSmVoJRarZTaq5Sap5Tyctf1dj/f694eXyevQNSYUoqfH76YthF+zP7toHTzNDKHJz+Cs7SUiPvuNToU0YjUOOkrpWKBe4G+WutugBm4EXgJmKK1bgfkARPcu0wA8tzlU9z1hMFMJsXjV3Uhq6iCaT/uNzoccY4cRUWUrl1H8E034tu7t9HhiEaktt07FsBHKWUBfIE0YDgw3719FnCt+/Fo93Pc2y9RMnjYIwzrEEGLAG+m/LCbrzalGh2OOAfZ774LTifBfxxjdCiikalx0tdapwKvAIdwJfsCYD2Qr7U+cp9/ChDrfhwLJLv3tbvrn7LYplJqolJqnVJqXVZWVk3DE+dBKcXUm3oBcP+8TTIpmwfLmjqVHYndyX1/Jv7DL8anW1ejQxKNTG26d0Jwtd4TgBjADxhR24C01tO11n211n0jIiJqezhxji5sE8bUm3qhNdwy4zcZzeNBHEVF5H/5JUl/+hPZb7+DsloJu/MOWjz4oNGhiUaoNvPpXwoc0FpnASilPgcGAcFKKYu7NR8HHOkvSAVaAinu7qAgIKcW5xd17KrEaHamFfL2in28tGgnj1/VxeiQmr3C7xaT9thjOEtLUT4++F00hMhHH5UpkUWN1aZP/xBwoVLK1903fwmwHVgOHOloHA985X78tfs57u3LtMz65VFMJsXDIzpxSacWzPjlAAWlVUaH1GxprSn++RdS778flCL29Sl0XLeWVtOnS8IXtVKbPv3VuC7IbgC2uI81HXgEeEAptRdXn/377l3eB8Lc5Q8Ak2sRt6hH4wfGA3Dje79hl24eQ+TO/IDkO+/E5O9Pq9mzCBwxAmU2Gx2WaAJqtVyi1vpJ4MmTivcD/aqpWw5cX5vziYZxUYcI+rYOYd3BPNYk5TKwbbjRITUbFQcOcPDmW3Dk5QGQ8OWXeMXFnmUvIc6d3JErqjVlbE8Abp25VkbzNJDKpCSSxlyPIy8P/4svpv0vP0vCF3VOkr6oVstQX67rHUulw8kHMgVzvStZvYZ9I0biLCkh+rnnaPnO21jC5RuWqHuS9MVp/Wt0NwJtFl5bspuCMrmoW1+01hwa7xrjkPDlFwT/8TqDIxJNmSR9cVp+3hZev7EnFXYnc9ccMjqcJitnxgwAfPr0wdapk8HRiKZOkr44o+GdIundKpgXFu3EIcsr1rnMV14h61XX+rUt333H4GhEcyBJX5zVsI4tAPjXN9sNjqRpKVq6lJwZ7+M3cABtFi3EHBBgdEiiGZCkL87q9sEJtAjw5sNfk6Sbp444y8vJnDIFgOjnn5cbrkSDkaQvzsrf28LKycPpEh3I5M+3kJxbanRIjVrF/v3sHjiIyr37aDl9GtaoKKNDEs2IJH1xTqxmE5OGtwNg5d5sg6NpvJwVFRwYcz26tJQWkx/Bb8gQo0MSzYwkfXHOrugaha+XmR93y5TXNVW4cJEr4T/8MGG33oosKSEamiR9cc5MJsVlXSJZtDWdvZnFRofT6DgrKsh++20sUVGE/nmc0eGIZkqSvjgvk4a3B+DNZXsMjqRxsefmcvihh6lKTqbFgw+gLLWa9kqIGpOkL85Luxb+9IgL4stNh9maWmB0OI1C+a7dHPjjGIqWLSPigQcIuvpqo0MSzZgkfXHenrjatUTfa0t2GxyJ50t/9jkOjB6Nrqqi9exZhE+80+iQRDMn3zHFeevTOoT4MF9p6Z+Bs7SUQ7dPoGzTJkwBAcS+9iq+vXsbHZYQ0tIXNTNuQDyZRRUyZv80sqdNp2zTJgJGjKDd8uX49TtliQkhDCFJX9TIkPauaX+XbM8wOBLPkzt7NjnTpuHTtw9xr0/B7O9ndEhCHCVJX9RI+xb+BNgs/EdG8ZxAOxxkTf0PALGvvmpwNEKcSpK+qBGlFJd2jiSvtIrU/DKjw/EYef/7H87iYqJffAFrZKTR4QhxCkn6osaOTMsw4+f9BkfiGdKfeYaMF15E+fgQePnlRocjRLUk6YsaaxPhT0SANx+sTKLC7jA6HEPlzfuEvI/nYG3VijZffoHJ19fokISoliR9USsTh7QB4M1lew2OxDhlW7aS/uSTmENDiZ83F6/WrY0OSYjTkqQvauX2wQn4epn5z7K9zXJlLV1VRdL11wMQ/8k8LCEhBkckxJlJ0he1YjYp/nShq2W7KTnf2GAMkPHCiwCEjv8zXnFxBkcjxNnJHbmi1u4a1pbpP+1n/vpk+rRuPi1de14eefPmYevWjRaPPGJ0OKIOVVVVkZKSQnl5udGhnJHNZiMuLg6r1XrO+0jSF7UW7OtFoM3C9sOFRofSYJwVFRx+8EFwOIh68gmUSb40NyUpKSkEBAQQHx/vsWseaK3JyckhJSWFhPNYblPeqaJOjOoew+8pBUz/aZ/RoTSInBkzKPl1FeF3341PYqLR4Yg6Vl5eTlhYmMcmfHDdKxMWFnbe30Yk6Ys68bdLXfPsf7gyydhAGoAjP5/s/7yJd/v2REy6x+hwRD3x5IR/RE1ilKQv6kSLQBv3Dm/H4YJyfm/iF3Qzp7wOQOhttxkbiGjynnvuObp27Ur37t3p2bMnq1evrvUxpU9f1Jmx/VoxddlevtiYSo+WwUaHUy9K168nf948rLGxBF/3B6PDEU3YqlWr+Oabb9iwYQPe3t5kZ2dTWVlZ6+NK0hd1JjbYh46RAaw+kIvWulF8PT4fWmtSH3oIgJbvTTc4GtHUpaWlER4ejre3NwDh4eF1ctxaJX2lVDAwA+gGaOB2YBcwD4gHkoAbtNZ5ypUB3gCuBEqBW7XWG2pzfuF5burXkqcWbGdfVgntWvgbHU6dKl66FPvhNMLuvBPvNm2MDkc0kKcXbKvzkWldYgJ50r0C3elcfvnlPPPMM3To0IFLL72UsWPHMnTo0Fqfu7Z9+m8A32mtOwE9gB3AZGCp1ro9sNT9HGAk0N79MxF4p5bnFh5ocPsIwPWH0pRou52UeyYBEPaXiQZHI5oDf39/1q9fz/Tp04mIiGDs2LF8+OGHtT5ujVv6Sqkg4CLgVgCtdSVQqZQaDQxzV5sFrAAeAUYDs7XWGvhNKRWslIrWWqfVOHrhcdq559lvSmP2HUVFJN/pSvQh48Zh9m9a32DEmZ2tRV6fzGYzw4YNY9iwYSQmJjJr1ixuvfXWWh2zNi39BCAL+EAptVEpNUMp5QdEHpfI04Ejk4rHAsnH7Z/iLjuBUmqiUmqdUmpdVlZWLcITRpk4pA05JZW8tbxpTMKW9sQTlG3aRPi9k4h87FGjwxHNxK5du9iz59giRZs2baJ1HUzmV5ukbwF6A+9orXsBJRzrygHA3ao/r1m4tNbTtdZ9tdZ9IyIiahGeMMqRuXheXrwLu8NpcDS1U7RiBUWLvsMrPp6Iu+5qchenhecqLi5m/PjxdOnShe7du7N9+3aeeuqpWh+3NhdyU4AUrfWRgaPzcSX9jCPdNkqpaCDTvT0VaHnc/nHuMtHEhPh58ejITrywaCc704voFhtkdEg1Ys/NJeX//gpA1DNPGxyNaG769OnDr7/+WufHrXFLX2udDiQrpTq6iy4BtgNfA+PdZeOBr9yPvwb+rFwuBAqkP7/purhTCwB+259jcCQ1V/DFl4BrBk2/fv2MDUaIOlLbcfqTgI+UUl7AfuA2XB8knyilJgAHgRvcdRfiGq65F9eQTbmdsQlrFepaOeq3/TncMaRxDm8sWrIE7y6daTF58tkrC9FI1Crpa603AX2r2XRJNXU1cHdtzicaD5vVzMhuUaw5kGt0KDVmz8zEp1cv6ccXTYrMvSPqTf+EUHJKKnlh4Q6jQzlvOe/PpOrwYXx69DA6FCHqlCR9UW/+2Me1ktT0n/c3qoXTHQUFZL35JpYWLQi+cazR4QhRpyTpi3oTYLPyxFVd0LpxTbmcP38+uqyMyEcnY/LyMjocIeqUJH1Rr24fnMDgduG8sGgnezOLjQ7nrOzZ2WRPfw+AgBEjDI5GNGdms5mePXse/XnxxRfr5Lgyy6aod49d2Zkrp/7Mp+uTeXRkZ6PDOaPUB/+Os6CAmJdflgu4wlA+Pj5s2rSpzo8rLX1R77rEBNKndQjTftxPfmnt5wOvL1UZmZSuXo05Ipygq68yOhwh6oW09EWDGNg2jPUH83ho/mbe+3N1o3yNl//JJwDEPP+8wZEIj7JoMqRvqdtjRiXCyDN315SVldGzZ8+jzx999FHGjq39wAJJ+qJBPHh5R77adJgl2zPYdriArjGeNTWDo6iIggULwGrFb/Bgo8MRot66dyTpiwbz3B+6Me79Ncz8JYlXb/Cc8e+VBw+S/H9/perQIVpOnyZ9+eJEZ2mRNzbSpy8azJD2Efh6mflsQ4rRoZwge9p0Kg8coMXDD+N/0UVGhyNEvZKkLxrU6J4xABzILjE4EpfyHTso+Pxzgq+/nrDbZToo4TmO9Okf+ZlcR3NASfeOaFDX9Y5jzppk3l6+l5evN76LJ+uNqQBy563wOA5H/dzFLi190aAuiA+lc3Qg321NN3z4ZslvqylesQJrbCw+XY1bEk+IhiRJXzS4l8d0p6jCzifrks9euR5lvvYaAC2nvWtoHEI0JEn6osF1jQnEYlJ8vsG4hdPs2dmUb95M+F134d2unWFxCNHQJOmLBqeUonWYLzvTi6gyaA3doh9+AMCnT29Dzi+EUSTpC0Pc1K8VACv3Zhty/rLfN4PFgt/AgYacXwijSNIXhri6h2vo5r+/22XI+YuWLCFg+HC5EUs0O5L0hSEiA210jAxge1ohKXmlDXruwsXf4ywuxle6doQH8/f3r5fjStIXhvnXtd0AeHj+5gY7p7O0lNT77gMg6I9/bLDzCuEpJOkLw/RLCAXg1305OJ26Qc5ZuGgRAFFPPYm5nlpSQngyuSNXGOovQ9sw7cf97EwvoktMYL2fL2/eJ1iiowmugylqRfPw0pqX2Jm7s06P2Sm0E4/0e6ROj3mupKUvDHVVouuC7o60wno/V+H331O+eTMBFw+TC7ii2ZKWvjBU15hAooNsvPjdTv7YJ67ezqOdTg5PfhSA8EmT6u08oukxqkVeX6SlLwxlMikGtA0jq6iiXhdOL1q8GF1aSsR992IJCam38wjh6STpC8NNvKgNAC99V7f9pkc4KyrIeP4FAEJuvrleziFEYyFJXxiuU1QgQztEsGR7Bj/vyarz46dMmoQ9K4vIf/wDc5BnLdMoxOkUF9fPN19J+sIjTBvXhwBvC89+s6NOj+ssL6fkp58xBwcTcou08oWQpC88gs1qpk2EH7syiiittNfZcfPmzAUg8tHJKJO83YWQvwLhMe6+2DXF8eJt6XV2zNzZs1G+vgRedVWdHVOIxkySvvAYl3aOJMjHysq9OXVyvKqMTOxpaYTdeivKbK6TYwrR2NU66SulzEqpjUqpb9zPE5RSq5VSe5VS85RSXu5yb/fzve7t8bU9t2haTCbFBfGhrE3KRevaTcvgrKhg/6hRAPgPvaguwhOiSaiLlv59wPFX314Cpmit2wF5wAR3+QQgz10+xV1PiBMM6xjBwZxSttfyDt3c2bNxFhfjP2wYtm7d6ig6IRq/WiV9pVQcMAqY4X6ugOHAfHeVWcC17sej3c9xb79Eyb3w4iSjEqNdv6f+gr2Gq2rlf/4FWa++hnfnzsRNfUO6dkSjZDab6dmzJ926dePqq68mPz+/To5b25b+68DDwJG/zjAgX2t9ZPhFChDrfhwLJAO4txe4659AKTVRKbVOKbUuK6vux2wLzxbi50WbcD8Avt+ecd77a63JfPVVAKKfeRrl5VWn8QnRUHx8fNi0aRNbt24lNDSUt956q06OW+Okr5S6CsjUWq+vk0jctNbTtdZ9tdZ9IyIi6vLQopH45t7BAMz85cB571vw+Rc4cnKIeuopfBIT6zo0IQwxYMAAUlNT6+RYtZlwbRBwjVLqSsAGBAJvAMFKKYu7NR8HHIk0FWgJpCilLEAQUDfDNEST4utloXWYLxsO5eF0akymc+8FLPhmAQBB1/2hvsITzUz6889TsaNupwjx7tyJqMceO6e6DoeDpUuXMmHChLNXPgc1bulrrR/VWsdpreOBG4FlWutbgOXAGHe18cBX7sdfu5/j3r5M13aIhmiyburXCqd2LbByrpylpZRt3ETQtddikm4d0ciVlZXRs2dPoqKiyMjI4LLLLquT49bH1MqPAHOVUs8CG4H33eXvA/9VSu0FcnF9UAhRrYs7tuDFRTtZcyCHwe3Dz2mfjJdfRpeXE+geqilEXTjXFnldO9KnX1payhVXXMFbb73FvffeW+vj1snNWVrrFVrrq9yP92ut+2mt22mtr9daV7jLy93P27m376+Lc4umqWNUAO1a+PPuj/vPacrlypRU8ufMxeTri9/gQQ0QoRANw9fXl6lTp/Lqq69it9d+ihK5I1d4rFeu70Glw8lnG1LOWrfk15UAxL4+RVbFEk1Or1696N69O3PmzKn1sWTlLOGxerYMZnC7cGb9msQt/VsRF+J72rqV+/YB4DdIWvmiaTh5auUFCxbUyXGlpS882iMjOlFa6eCl73adsV7RD0uxde0qN2IJcRaS9IVHS4wLomWoDxsO5p22Tum6dVSlpuLdvn0DRiZE4yRJX3i8cRe2JjW/jJziimq3533yCQAh4/7UkGEJ0ShJ0hcer1cr10Lmc9YcOmVb8U8/Ufj1AgJHjcKna9eGDk00YY3hNqKaxChJX3i8vq1DiA324dUlu08od1ZUkHL3PQCE/99fjAhNNFE2m42cnByPTvxaa3JycrDZbOe1n4zeER5PKcXVPWJ498d9/LIn++jNWmUbNqCrqoh6+mnpzxd1Ki4ujpSUFDx90kebzUZcXNx57SNJXzQKtw6M590f9/Gn91ez+9mReFlMFP/0M5hMBI660ujwRBNjtVpJSEgwOox6Id07olGICrJxU79WALyzYh8Ve/aQ+8EH2Dp3xuzvb3B0QjQekvRFo/Hctd3w8zIzY+FGkv56NwCht99mcFRCNC7SvSMaDZNJ8eQ1Xan62904s5KJefUVgmRyNSHOi7T0RaMyYMbz9M7aw4+t+uB9+QijwxGi0ZGkLxqNkl9/pXjFCoq79OS1HmN47tsdRockRKMjSV80Godud60clPji0/j4+/Lf3w6SkldqcFRCNC6S9EWjUJmcDIBP3z7YOnRg6k29ABj80nJKK2s/x7gQzYUkfdEoZL/5FgAt7r8fgKEdIrhziGsc9dvL9xkVlhCNjiR94fG03U7BwoUA+PTpc7T8sSs7A7Bwa5ohcQnRGEnSFx4v9aGHoKqKsDvvPGFVLKUUT17dhf1ZJXy5MdXACIVoPCTpC4+WO/u/FC36DmtMDBEP/O2U7Rd1iADg/nmbWLwtvaHDE6LRkaQvPFb+F1+S8fzzAMS98061a9+2jfBn1u39AHjo098bND4hGiNJ+sIjOQoKSH/ySQDiP5uPrWOH09Yd2iGCUd2jKSy3sy+r+LT1hBCS9IWHKvz+e3RlJa0+mHlOi6NMHNIGgPnrU+o7NCEaNUn6wiMVfPY5lshIfC+88Jzqd48LAmDJ9gzsDmd9hiZEoyZJX3icqoxMyjZtwrdfv2r78aujlOKei9uxN7OYMe+uqucIhWi8JOkLj5P2+OMABF8/5rz2u//S9phNik3J+WxJKaiP0IRo9CTpC49SuGQJJT//jCUyEr9+/c5rX4vZxCd/cXUH3Tt3Y32EJ0SjJ0lfeJS0Rx8DIP7jj85ad97OeczYMuOEsj6tQ4kJsnEgu4SH58sQTiFOJklfeITKlFQOjvszzuJi/IYMwRobe8b6M7fO5NnVz/LGhjdInJVI4qxE8svzAVj292HEBNn4ZF0KezNlCKcQx5OkLwynHQ5S7rqL0rVrscbGEvf6lLPuM2X9qXWGzBtC4qxEKpzFzJno6uYZMXUxSdklACQXJbMufd3RD4mF+xfW7QsRohFQWmujYzitvn376nXr1hkdhqhnhx9/nIL5nxH0x+uIee65M9YtqChg8NzBAITZwgjzCaNdcDsWHjgxgScEJXCg4AAA9uR7GX3Rfr498M0px7so7iKuaXsNV8RfUUevRgjjKaXWa637VrutpklfKdUSmA1EAhqYrrV+QykVCswD4oEk4AatdZ5yjb17A7gSKAVu1VpvONM5JOk3fdrhYGfXbgB0WLcOs79ftfUcTgcZpRlc8dmx5BwfGM+CPywAYMaWGbyx4Y1axzM4djDvXPpOrY8jhJHOlPRr071jBx7UWncBLgTuVkp1ASYDS7XW7YGl7ucAI4H27p+JgPxlCYoWLwYg8p+PnzbhA/T8b88TEv4DfR7g06s/Pfp8QrcJPHLBIwyMGXjG8y0Zs4S1t6xlUq9J1W7/JfUXEmcl8s4meXuKpqnOuneUUl8Bb7p/hmmt05RS0cAKrXVHpdQ09+M57vq7jtQ73TGlpd+0OcvL2XvpZTiys+n4+yZM3t5UOapYeXglk5a5kvKvN/3KwDknJvKXhrzElW2urPaYWmvK7GX4Wn2Plr25Zi5vb5iNOX0SW54aebReSVUJXmYvHvrxIZYlL2NQ7CBWpq48ut9rw17jstaX1fXLFqLenamlb6mjE8QDvYDVQORxiTwdV/cPQCyQfNxuKe6yE5K+Umoirm8CtGrVqi7CEx4qb85cHNnZRD72GCZvbwB6/6/3CXVOTvjLb1hOuE/4aY+plDoh4QPc0+9GNm5vy9LyTKYu3cO9l7RHKYW/lz8Abww/1i1UUlXCxzs+ZurGqTyw4gFevuhlRiSMqNXrFMKT1Hr0jlLKH/gMuF9rXXj8Nu36GnFeXyW01tO11n211n0jIiJqG57wUM6SEnJnzUJZrYTccjMA23O2n7b+f0f+lyVjlpwx4Z/Ja2N7un4v2U1mUflp6/lZ/biz+51M7D4RgHm75lFmL6O0ShZgF01DrVr6SikrroT/kdb6c3dxhlIq+rjunUx3eSrQ8rjd49xlohlKe+pp7OnphN/1V5TZDMDYb8YC8OSAJxnTYQylVaXHWu0VxWC21vh8QT5Wvr5nENe8uZJRU39h0X1DCPf3Pm39Sb0m8dvh31iXsY5+H514Z/CXo7+kbXDbGscihJFq3NJ3j8Z5H9ihtX7tuE1fA+Pdj8cDXx1X/mflciFQcKb+fNF0VaWmUrhgAZjNhN9zDwBpxcfeCmM6uObcOZrwC1LhhVh4NvKUY52P7nHB/OPKzmQVVdD32R+Y8fN+HM7TfxG9IOqCasuv/epa7lh8R61iEcIotWnpDwLGAVuUUpvcZY8BLwKfKKUmAAeBG9zbFuIarrkX15DN22pxbtFI2bOz2XvJpQCU/GkUhVVFFFYWcuXnrguzp4yq+eL/4Pc57icangqCJ/PhwI/QejCYz+8tfOdFbdifXcKcNYd49tsdpBWU88+rulRbd2L3iYTYQugQ0oE+kX0od5QzaM4gAFanr2Zl6ko6h3Um1BZ6XjEIYSS5OUs0GO10cmD0aCr27GVJT8V7I0xw0tTJm/+8GVWaC76h8OVd8PvHpz9g7/FwzdQaxWJ3OGn3j0UAfHRHfwa1O/drBTtzd3L9gutPKJszag4JQQmsTV9LanEqt3S+pUZxCVEX6uXmrIYgSb/pcFZUcOCa0VQePMjSHoppV5pPqfN8VjZXdxgDG2afeoDBf4OAaFj08LEykwWeyKlxTL/uy+bm91bTJTqQBZMGYzad29z9AP0+6keZvey027uHd2fmiJl4m09/3UCI+iJJXxgu//MvSHvsMQ6Fw9/vMDMovAcVKWuIdDhc200m3s3IOnXHQffDZU8fe/7D0/DLa6fWO9kVz0OPm8AWBKZTP2COeOjT3/l0fQrjLmzNv67tdp6vClanreaO78/ev39Z68t4deir57wojBC1IUlfGKbMXsbURf/kmoe+BWDsZDNaKTaXBaHStxyreNdv8HY1SyM+Vc1iKJWlsGYa/PDUuQUxfgFE9wRb4CmbyqscXPLqj6Tml9E/IZTZE/rhbTn9h8TpFFcWYzFZsJqsrEpbxV9/+Gu19SzKwg/X/0CYT9h5n0OIcyVJXxhj8yckbvwXExY7uGKDZvZwE9/0N7HpwCGOptWHD4CXP1i8XM+3fwWxfcAvAixn6RrZ+S3MvRm6/RG2fnb2eO5c5jr2SdIKyhjzzipS88sY1C6Mj+44t3V5z8SpnczfPZ8Inwh8rD7M3DKTVWnHlnG8tNWlrM1Yy4obVmAx1ck9kkIcJUlfNKzSXPh3AokJreiQonn2vw4Oh8D9fzGzJem4m7J9QuGRA/UTg6MKijNhykkjc/6RAVbbKdWdTk37xxfhcGo2PXEZwb5edR5SUWXRKXcYA9zY8UZ6teh12qklhDhf9TXhmhCnStsM/05gTEwU0TmayZ+6+uz7J+acmPAB7vih/uIwWyEo1tU99I/0Y+XPRULK+lOqm0yK/03oD8Dot1aSX1pZ5yEFeAXw09ifTimfu2suj/z8CN/u/7bOzynEyaSlL+qO0wnPhPBRoD8vhoUy9R07UfkQ//GH+PTuDwdXQdZO6P3nM15crRflBfDicXM5PZEHplPbPJM/28zcta4Pp1Hdo/nbpe1p1yKgXkJyaidf7/uaj3d8zI7cHQD8PPZngm3B5JXnkVSYxPqM9dyRKDeCifMj3Tui/j0VxBYvL2YGB/KDny999jh5ZL4TvyFDaPXedKOjO+apoGOPn8w/5T4BrTUvfreTD1cmUWF3ArDkbxfRPrJ+Ev8RIz4bQWrx6Wclef3i1/EyeTF311zu6nEX9y67l67hXXnj4jdkRJA4hSR9Ub9Ksln7ny7cHh2Jxa559BMniQc11tatiJ87F0tIiNERHlOWDy+1dj2+dSHEDzpt1Y9WH+QfX2ylX0Io8yZeWK/JVWvN7O2zeWXdK+e134zLZ9A/un89RSUaK0n6on64u0yW+/pwb6RrRtR/zHHQI8n1nmqz8Fu827QxMsLqpayHGcNdjydtgLDTT542Zclu3li6B7NJMaBNGB/edgEWc/1dCtNaszx5OV3CuhDlF4VTO/k+6XtSi1NJK0nDYrKwMXMjQ+OG8s7v73BDhxv454B/1ls8onGSpC/qxbb5f+avhevIc8+SOeJgCLd/nIWtSxdazXwfc3CwsQGeycKHXWP9j/fATvjxRbjsGddNXbiS8Evf7eLdH/cB0CLAm4/v7F9v/fzn45Zvb2Fz9mZZ4lGcQpK+qDO55blMWT+FL/d+ebTMt1zzt91x9Pj2IADtfvwRa2QLgyI8R1rD08Gn3973dhj12tE+f601j32xlTlrDgEwbVwfruga1QCBnl5+eT5D5g0B4Lr213F9h+vpFn7+dxWLpkeSvqgzT/36FJ/tOXYjVGy2Zsp7rmGZyteX2FdeJmD4cKPCOz+FafDmBdDpSihIgYMrT1/3of3gF8b329KZ+F/XkM/HR3XmtkEJ5zVnT137KeUn7ll6DxpNqC2U/438Hy0DW559R9GkSdIXtaM1bPqIn8LiuHv5samPf96fTPaqflQdSibysUcJGTeuaYwkObwJpg89sSysHUxyJfu9mUVc+pprvL2P1cyN/Vry+KguhiX/SkclO3J3MPH7ididdv5+wd+5rv11MtlbMyZJX9TM4Y2uOWueDubRiDC+8fcDYEZaBv3KKji8KpjCQ75EPv44oX9qYlMJO52QstY1PcT0YcfKrb4w+RAFe37ly92VvLzeSXGFnRBfKxv+eZmhH3qbMjcxbtE4AK5pew1PDXgKay1WGxONlyR9cX7c0yikWMyMbBl7wqZBpWW8fTiLXZ/GHC1r99OPWFt4eB9+bZRkw8vVj/DRfe9g0JYrOVxQzsxb+zK8U+1W96qtHTk7+GDbByw6sIj2Ie25qdNNjIwfeXQReNE8SNIXZ5a1G966AK6eCm0vhtcTAUhMaHVCtck5udxSWMz+dRdRsXcvtsREWr033bNH6dQVrV3ffN67+IzV8jrdTMiNxo+k+d/2//HS2peOPo/0jWTxHxdjbug7oYUhJOmL6hWkwoGf4Mv/O1qkcc1tf09kBJttrj7h1X2ewtz5Khy/byN31myKFi/Gu307Er7+umn04Z8PpxO0E9Cw/0f46I+nVEm1xpN64VP0G37tKXf8NqRfU39l2uZpbMjccLRsRPwI+kb25boO11FQUYC32ZsAL+OHn4q6JUlfHFOaC15+8O4QyN51tFgDN8dEstX7xIt/3177Da2CWpMz8wMy//1vMJnwbtuWVrM+xBIqa8OStBI+vBLCO1JwwX0ELbrrhM1L+n9AVOIldIsNNOwDstxezgtrXuDzPZ9Xu/2OxDu4teutBHkHVbtdND6S9IVLRRG8EHdq+fDHmervzXtb3jta1C2sGzOumIGf1Y/UBx6kcOFCsFhot+R7rNHRDRh0I+N0UP7hddgOrTih+DdTL7JaDGJvzDVEtIhiTJ84bNaG7Wpxaiff7v+W9RnrTxh2C+Br8WX+NfNpGSDDPZsCSfrCZc8Px7oj2l8ON80Dpx1tttJ9dnfg2CyPANrh4MC1f6Bizx4A2n6/GK9Wrao7sqhG0f61BMy+9JTyJGckV1U+R+8OrYkK9MZmNXNlYjT94kMxNdCwzypnFVaTlazSLL4/+D3/XvtvhsUN45Vhr2A1yYifxk6SvnCNQHmjB1QWwyNJ4BNClaOKf/zyDxYlLQLgzsQ7ubf3vQCU79xJ2j+foHzLFpTVSseNG1AWWeGpRipLXf/uix6GbV8cLR7u8yllTjNpBeWAa4qHCxJCsZgU3eOC+dOFrWq0dGNNzNo2i1fWvUKIdwgjE0ZyaetLaRnQkkjfyOZ33aYJkKTf3BVlwKsdXI/7/xVGvkhBRQHXfXUdmWWZR6v9etOveCVnkj1tGoVfLwAg6NpriXnxBSOibpoqS+H547rHbppLduxwFm9L54ftGWw4lE9BWdXRzVGBNoJ9rVTanVzeNYq/X96hXiZ801rzS+ovfLXvK5YdWkaV0xVDx5COjGozihHxI4jyi5IPgEZCkn5zVpSO/vh6UrK28feoGLxjerExa9PRzeO6jGOITyItft2D1/rtFP/4IwC2rl0Jv+duAi4+8xBFUQOVpTD3Jti/wvW8w0gY9ojrRjityS2zM/OXA2xJLcDbYiK/tIo1SblHd48L8cHf20KonxeXd4nEz9vC1T1i6uwaQXZZNs+seoaNmRvJr8g/Wh5gDaCoqoiXh77MiPgRdXIuUT8k6Td15YXgtMP/roPoHjDyZVc3QkUhH6x8mtdCq5nPXmsuTvLnsaCx5M/8AF3latl5t29PzMv/xtapUwO/iGaoKB1e7XhqeUxvuOSf0PbYHEYlFXY+Wn2QvZnF5JdWsfpA7gnfCFqH+dKndQjmWrTEbVYzCeF+eFlMdIsNonN0AA5dwcbMjSQXJTNv1zz25u8FwMvkRb/oftiddiocFQyIHkC38G4MiBkgC717AEn6TdHWz2H+bacUFyvFZm9vVvramB0UeLQ83i+WCYl3cmXIYAqXLCHr2eddNxwBlhYtaPHww/j174clIqLBXoIA7BWw6i3Xj7c/FB4Gh3t9XndX3OmUVtopq3Tw4a9JrNqXQ1pBOTX9e3ZoTUZhxQllXhYTbcL9CLRZaRnqS1SgFyH+ik0FC9lauBhvi5W00kOYlBmHtgNgM/swqtUtXJYwlD5R3TCpuuuKspqVdC+dI0n6Tc26D+Cb+08ocgJvhAQxMzgIpTUBpdB7n6anOZIhSSGQmoEjJ+dofVNgIMHXXUfYHRMwBwaivLwa9jWI6mkNyath5hWu5y26wO3fgdkb9iyGTlfV2/rCpZV28kqr2J1RxL7MYrYfLqSk0k5GYQV7MoooqXRUs5cdsICpHIv/Dqwhq7H4JrleitMC2oK9pC2Okg7Yi7qCcqDtNbsfwGpWhPh6EeRjPWFyO18vMwnh/kbeB1cverQMZtyFrWu0ryT9pqA0F/IPwsc3QnG6qxtn/AL0lvns004m7f8E86E0xi910u3gqf+n1thYvNq2wbd3b/wGDsSWmCitJk/mnv/oFL5hcNHD0OdWsNoaLBytNfmlVWjA7nSyO70YRzW5w6kd7M7fwrbsnRwu24tGk1y2kTJn/gn1gi2xdPK/hGJ7DrG2RFr59MasTj9UVGtNSaWDvJJKdxzHzn04v5zcksq6eqke45LOLXhmdM3WR5Ck31hVlnJo3xJiN8/HvONrAPZZLfzmYyNuxDu031pE6vPP4l9+4m5ebdvi178//hcPw9a1KyabDZOvb8PHL2pHa0j6GZY9C0VpkH/oxO19boMrXwaTxbUeQOFh17QPQXGuDwcAi/HTK2ut+S3tN77a9xXf7v+WYXHDKKgsYGPmxqN1Yv1jua/3fQyMGSh3BtcBSfqeyF4JG2dD6kYoTIGITvweHk9R9g5W7ltISat+rMncSKrFTGwO3LzSgdOkKAe6HtSEF7kO4wR+6KUY1Gs08UEJBI64Aq/WNftKKBoBRxUc+g1mXXVu9aO6Q69x0Pc28KBplp3ayY6cHZTaSzlQcIA3N75JXkUeZmWmTXAbgr2DCbAGEOoTSuuA1iQEJRz9MIjxj6GFbxOe1bUOSNJvaEf+Tff+AN6BENwSqsrAaacqYys7175Nu7Tt+FSWUqwUP9t8WOflzfde/vTfpWmVpfGp0gQ5nMSkKcLyj3XDlEYEYFVmKnHy9RVB+A4cwMPDnqzTC2aiEcg7CMv+5foAKEiGwQ9A64Gu917+QagohIztsHW+q35QK4hKhA6XQ5thEBgHZs8ZZeNwOtiSvYWfUn5iT94eCisLKaoqIrs0m7yKvFPqtwpoRZvgNtiddnLKcjArM0HeQQR6BRLoHUigVyBB3kF0DO1Iu+B2hNnCUErh1M5m8bfiUUlfKTUCeAMwAzO01qcdnlDjpJ+1C+beAoHREBBz4m/fMHA6XF+VzVaI6QXh7V1/LNpZ/UUyrV1DIt0tJe10ostLXV+9HVVQUcQhRxkF+7cTvX4BWymj0OJNmbZz0GqlyGkiOheKlYlcs5k26ZqWWRofpxmH00FUHlicJ57SEh6CslhQNj98evQg4JYb8e2WKFPjivNTlA47v4FvHzx1W+/xMOJFsPq4uoUcdqgqhcwdrm+f9goIiXcdI2uX628mKBbi+kFJJoS2AVswRHY59dhOh+vvqQ6+XRRUFJBUmERJZQkO7WBbzjZ25+1mf/5+vC3ehPuE43A6KKwspLCykIKKAooqi3DoYxeew2xhVDgqKHeU0y+qH71a9CIuII5wn3C6hnVtcjONekzSV0qZgd3AZUAKsBa4SWu9vbr6NU762XtcraDCNEqL0sgoyyTTBJkWMxVKEVtux8sJ2RYz3/j40THFi3Z5dhzaTrHyIi7TTGC+xmTXR6fR1WjsyozdZMa7UGOpqt2/W2mQN+ltQ7CaLFQF+9G762VYzV749u2LrUtnTLaGu0gnmomyfMjdD7sXw48ntbXC2kPeAVfj5nQsPuCocP9NHCcgBmzHhgfjqIT8ZNAOCG17YkMqsht0uw58Ql3XIkLbgF/Yub8GrV0fRmhAua5ZaOexYa5Hq2kKKwvZlruDnXm72Z2/j2CfUIoqi9iWvY19BftOfAleAYTZwgixhRBqCyXUFnr08ZFyX4sv3hZvbGYb3mZvbBbXb7Myo5TCpEwoPGNYqScl/QHAU1rrK9zPHwXQWld7n39Nk/6OtYvJePBh7NqO86Q3qE8FR/vDTyczCJLDFZXVNFIUYDdBeghUWBUBJivRPi0IxoSy+uNs1ZqYuM60CmyFUgpfsy8o8GrdGrO/e/UipTAHycUqYSCnE3Z8Bbu/h6oSV5nFx7UWsMkEZi8Ibu26d8AWBCEJrlZ9cQZk7XBdk7J4Qc4+17KS9uNGEygzBLcCNOQlHSuvKncNOz1ZcCvXMpRnorVr/qLSnBPP5RXg+nCpKj37a7bYXN2ttkBKTWYylSbNpNlscpKrNLlKk6c0Oe7H+QqcNcjfJg0mXLnCdNKPwtXFofSRMnXCNtOR7cAQ/wQeuv6r8w+AMyf9hu7UiwWSj3ueAvQ/voJSaiIwEaBVDWd09PcPZX9cCN5mGz5mb2xmb7yPfCr7+FDcIgCnxYSPxYcQnxBs7Tqgu7bHavHGx+pDqLkMv6Lko/OPAJiVmUi/SKL9olFKUep+kwV4BTSLPkLRxJhM0PUPrp/zERjt+jmi7XDod+e575+xDYrd8z05Kl33JOTsPbd9vfxd3bO2INe3B+10HauyxPWN4UwtbIcdKgpc04uXF+KrHcQD8cCA0+2iNYU4yNVV5Go7ZTgpx0mF1q7f7scajRNwuAeSOgGnu8zpLnPoM2xzP9fHlTuBKL+oc/t3OU8N3dIfA4zQWt/hfj4O6K+1vqe6+o32Qq4QQhjoTC39hm6ipgLHr9IQ5y4TQgjRABo66a8F2iulEpRSXsCNwNcNHIMQQjRbDdqnr7W2K6XuARbjul4xU2u9rSFjEEKI5qzB787QWi8EFjb0eYUQQjR8944QQggDSdIXQohmRJK+EEI0I5L0hRCiGfHoWTaVUlnAwVocIhzIrqNwGkJjixck5obQ2OIFibkhnCne1lrratc+9eikX1tKqXWnuyvNEzW2eEFibgiNLV6QmBtCTeOV7h0hhGhGJOkLIUQz0tST/nSjAzhPjS1ekJgbQmOLFyTmhlCjeJt0n74QQogTNfWWvhBCiONI0hdCiGakSSZ9pdQIpdQupdRepdRko+OpjlJqplIqUym19biyUKXUEqXUHvfvECNjPJ5SqqVSarlSartSaptS6j53uSfHbFNKrVFK/e6O+Wl3eYJSarX7/THPPc23x1BKmZVSG5VS37ife3q8SUqpLUqpTUqpde4yj31fACilgpVS85VSO5VSO5RSAzw5ZqVUR/e/75GfQqXU/TWJucklfffi628BI4EuwE1KqS7GRlWtD4ERJ5VNBpZqrdsDS93PPYUdeFBr3QW4ELjb/e/qyTFXAMO11j2AnsAIpdSFwEvAFK11OyAPmGBciNW6D9hx3HNPjxfgYq11z+PGjXvy+wLgDeA7rXUnoAeuf2+PjVlrvcv979sT6AOUAl9Qk5i11k3qB9eSl4uPe/4o8KjRcZ0m1nhg63HPdwHR7sfRwC6jYzxD7F8BlzWWmAFfYAOuNZmzAUt17xejf3CtJrcUGA58g2uNbI+N1x1TEhB+UpnHvi+AIOAA7oEsjSHmk+K8HFhZ05ibXEuf6hdfjzUolvMVqbVOcz9OByKNDOZ0lFLxQC9gNR4es7urZBOQCSwB9gH5Wmu7u4qnvT9eBx7GtYY2QBieHS+ABr5XSq1XSk10l3ny+yIByAI+cHejzVBK+eHZMR/vRmCO+/F5x9wUk36ToF0f3R43nlYp5Q98BtyvtS48fpsnxqy1dmjXV+I4oB/QydiITk8pdRWQqbVeb3Qs52mw1ro3ri7Vu5VSFx2/0QPfFxagN/CO1roXUMJJ3SIeGDMA7us51wCfnrztXGNuikm/MS++nqGUigZw/840OJ4TKKWsuBL+R1rrz93FHh3zEVrrfGA5ru6RYKXUkVXjPOn9MQi4RimVBMzF1cXzBp4bLwBa61T370xc/cz98Oz3RQqQorVe7X4+H9eHgCfHfMRIYIPWOsP9/LxjbopJvzEvvv41MN79eDyufnOPoJRSwPvADq31a8dt8uSYI5RSwe7HPriuQezAlfzHuKt5TMxa60e11nFa63hc79tlWutb8NB4AZRSfkqpgCOPcfU3b8WD3xda63QgWSnV0V10CbAdD475ODdxrGsHahKz0Rcl6ulCx5XAblz9t/8wOp7TxDgHSAOqcLU8JuDqv10K7AF+AEKNjvO4eAfj+uq4Gdjk/rnSw2PuDmx0x7wVeMJd3gZYA+zF9TXZ2+hYq4l9GPCNp8frju1398+2I39vnvy+cMfXE1jnfm98CYQ0gpj9gBwg6Liy845ZpmEQQohmpCl27wghhDgNSfpCCNGMSNIXQohmRJK+EEI0I5L0hRCiGZGkL4QQzYgkfSGEaEb+H4jzhfNoB4psAAAAAElFTkSuQmCC\n" }, "metadata": { "needs_background": "light" } } ], "source": [ "plot(res[0],res[1])\n", "legend(['S','E','I','R']);" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [ { "output_type": "execute_result", "data": { "text/plain": [ "inf" ] }, "metadata": {}, "execution_count": 17 } ], "source": [ "res[0][-1]" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ] }