{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Tyson1991" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "from ecell4.prelude import *" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This model is described in the article:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/markdown": [ "Tyson JJ, Modeling the cell division cycle: cdc2 and cyclin interactions. *Proceedings of the National Academy of Sciences of the United States of America*, **16**(88), 7328-32, 1991. [10.1073/pnas.88.16.7328](https://doi.org/10.1073/pnas.88.16.7328). PubMed PMID: [1831270](https://www.ncbi.nlm.nih.gov/pubmed/1831270)." ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ref = citation(1831270)\n", "ref" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/markdown": [ "Abstract: The proteins cdc2 and cyclin form a heterodimer (maturation promoting factor) that controls the major events of the cell cycle. A mathematical model for the interactions of cdc2 and cyclin is constructed. Simulation and analysis of the model show that the control system can operate in three modes: as a steady state with high maturation promoting factor activity, as a spontaneous oscillator, or as an excitable switch. We associate the steady state with metaphase arrest in unfertilized eggs, the spontaneous oscillations with rapid division cycles in early embryos, and the excitable switch with growth-controlled division cycles typical of nonembryonic cells." ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from IPython.display import display, Markdown\n", "display(Markdown(\"Abstract: {}\".format(ref.abstract)))" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "with reaction_rules():\n", " YT = Y + YP + M + pM\n", " CT = C2 + CP + M + pM\n", "\n", " ~Y > Y | 0.015 / CT\n", " Y > ~Y | 0.0 * Y\n", " CP + Y > pM | 200.0 * CP * Y / CT\n", " pM > M | pM * (0.018 + 180 * ((M / CT) ** 2))\n", " M > pM | 0.0 * M\n", " M > C2 + YP | 1.0 * M\n", " YP > ~YP | 0.6 * YP\n", " C2 > CP | 1000000.0 * C2\n", " CP > C2 | 1000.0 * CP\n", "\n", "m = get_model()" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "C2 + CP + M + pM > Y + C2 + CP + M + pM | (0.015 / (C2 + CP + M + pM))\n", "Y > | (0.0 * Y)\n", "CP + Y + C2 + M > pM + C2 + M | ((200.0 * CP * Y) / (C2 + CP + M + pM))\n", "pM + C2 + CP > M + C2 + CP | (pM * (0.018 + (180 * pow((M / (C2 + CP + M + pM)), 2))))\n", "M > pM | (0.0 * M)\n", "M > C2 + YP | (1.0 * M)\n", "YP > | (0.6 * YP)\n", "C2 > CP | (1000000.0 * C2)\n", "CP > C2 | (1000.0 * CP)\n" ] } ], "source": [ "show(m)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEJCAYAAACZjSCSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdd3hUVfrA8e+Z9EoqNYRApAcIvRcFpCggKs2+KFh/6tpdV+zd1VXXhopYQSwUKQooRVBaqAGkJySQkIT0nsmc3x8ngQApk2Rqcj7PkyfOzJ1735gw7z3tPUJKiaZpmtZ4GewdgKZpmmZfOhFomqY1cjoRaJqmNXI6EWiapjVyOhFomqY1cjoRaJqmNXJWSwRCiHlCiBQhRGwVrwshxLtCiKNCiL1CiF7WikXTNE2rmjVbBPOBsdW8Pg5oX/Y1G/jQirFomqZpVXC11omllBuFEBHVHDIJ+FKqFW1bhBABQogWUsqk6s4bEhIiIyKqO62maZp2sZiYmDQpZWhlr1ktEZihFZBQ4XFi2XOXJAIhxGxUq4Hw8HB27NhhkwA1TdMaCiFEfFWv2XOwWFTyXKX1LqSUc6WUfaSUfUJDK01omqZpWh3ZMxEkAq0rPA4DTtspFk3TtEbLnolgGXBL2eyhAUBWTeMDmqZpmuVZbYxACLEAGAGECCESgWcANwAp5UfASmA8cBTIB/5hrVg0TdO0qllz1tCMGl6XwL3Wur6maZpmHr2yWNM0rZHTiUDTNK2Rs+c6Aq2u8s5Cwhbw8IeIISAqm4mrOQSTCY79Bi5uEBgB/mHgov/ZOTQp1e8sOwk6TwCvAHtHZHX6L9KZ7PwK/nwX0g6ffy6sLwx/Ai4bqROCI1rzNPz1v/OPXdyh/10w4klw97ZfXNqlpITDv8CG1+D0LvXcykeh62ToNwtaNdxyaLpryFnkJMOKh8HVA0Y+A//4Ba7+r3r+m+vgiwmQc8beUWoVbflIJYG+d8CtP8PE96DrtSqZfzgITmy0d4RaRb88CQumQ366+l3N+h2iZ8DBn+GTK2Ddy2AqtXeUViGcbfP6Pn36yEZZYuLXp2DLB3DfDgiOPP+8sRh2fgFr5qiuoinzoc1Au4WplTmwDBbdAp2ugqlfgsHl/GsnNsKy+yHjBAz5J1wxBwz6nsyuinLhzQ7QYQxcO1d15Z17LQdWPgZ7voWIoXDtJ+Dfwn6x1pEQIkZK2aey1/RfnzPIS4Md86DblAuTAICru2q23vEbuPvAF1fD1o/tE6emZMTBT7MgrA9c9+mFSQCg7TC45y/ofRtsehu+u0l9EGn2c3AZlORB/zsvTAIAHn4w+UO45kM4FQMfD4XEGPvEaSU6ETiDv96HkgIY+nDVxzTrArPXQfsxsOoxWPVEg23GOryja8FYCNd8BG5elR/j5qW69sa9DodXwbwxkJlQ+bGa9e3+FoLaQev+VR8TfQPMWgdu3jD/Kji43HbxWZlOBI6uIAO2fQJdJkFox+qP9WwC076GAffA1g/h+9tUAtFs6+QW8G1+aevtYkKoO9Abf1BJ4LMr4cx+28SonZcRB3F/qA/6miZcNO2kWt/NuqiW3JaGsY2KTgSObuvHUJwDwx4173iDAca+AmNeVoNcX16jkolmOye3QPgA82dxXTYSZq4CJMwbB3GbrBqedpE9CwEB3aebd7xvKNy6XI3//PKEGp9zsrHWi+lE4Ohif4R2l0PzqNq9b+C9MOVzOL1Tfbhk68KuNpGZAFkJ0GZQ7d7XrCvcvgb8msNXk+HAUuvEp13IZFLdQu2GQ0Drmo8v5+6tJgH0vQM2vwNL7obSEuvFaWU6ETiy0hJIPw4te9bt/V0nw00/Qlai6nZIPVzze7T6SdiqvocPqP17A1rDzF/U73vRrbD9M8vGpl3q5J+QGQ/RN9b+vQYXGP8mXP4U7FkAC2ZAcZ7lY7QBnQgcWUYcmIwQ0qHu52g7DP6xAoxFMO9KSGyEU29tKf5PcPeDpl3r9n7vILh5iZrGuOIhWP+q03c7OLS936nfV6er6/Z+IWD4YzDhHbUa+ctJah2Ck9GJwJGlHVHfQ9rX7zwtesDtv4JngFp4dmRt/WPTKndyC7TuW78yEu7eatC/xw2w/hW1kFDPALOO5H3q91XfVd69b1NdRUl7Yd5Y1Qp3IjoROLLyUhLBl9X/XEHt4PbV6lwLpsGe7+p/Tu1CBRmQcgDCazk+UBkXN7jmAxj8AOz4rGwGWGH9z6tdKCMeAtpY5lydJ8DNP0FOkuqKTTlomfPagE4EjuzsEfBparmiV75N4bYVaiBz8WzY/K7udrCkhO2ArNv4QGWEgNHPw5hX1IKnr6+DwizLnFtTK4YL0iHQQokAVBHIf6xUXbrzxqgWohPQicCRpR2tf7fQxTz91bz1rpNVQbRf/6VmTmj1d/IvMLhCq96WPe/Ae+DaT9VA9OfjVVVMrf7KF/AFhFv2vM27qRlgPqFqzMAJFp7pRODI0g5bPhGAKlx33Tzof7eqX/TjTDWYrNXPyS3QIto6VUW7T4Ebv1cTCD4bDamHLH+NxiYzXn23VNdQRYFtYOZqNS140c2w/VPLX8OCdCJwVHlnVbM12AqJAM4vPBv9POxfDF9dqxee1YexSNWhsWbBv8jLVbdDabHqg47/y3rXagwyT6rv1kgEAD7Bqups+yvVgP/a5xy2K1YnAkd11kIzhqojhBqMLO92mDdW17upq6Q9UFoErS00PlCVFj3UoL9PiOp22L/YutdryDJPgquX+n9pLe4+MO2bsgKDb8HiO1XFYAejE4GjstTUUXN0n6JmO2Sfhk9HqQ81rXYyyroZ6rPmw1yBEaoPumVPNZtID/rXTUacGh+w9oZOLq6qwODl/1brFr6+FgoyrXvNWtKJwFGlHVa7WVmr2XqxtsNg5q9qsHPeODi82jbXbShyk9V3v+a2uZ53ENyyFLpcowb9Vz4KpUbbXLuhyDxp2RlD1REChj8Kkz9WY0nzxpy/eXAAOhE4qrNHISjy0lr21tSsC9yxFkLK1hps+8R213Z2Ocng5qNq19uKmydc/zkMuh+2fwILZ6gpkZp5Mk9afsZQTXpML2t9J6nW9ynH2NdAJwJHlXZEfSDbmn8LuG0lXDYaVj6itu/Tq1prlpMEfs1sv2+0wQBXvgBXvw1Hf1OtOSdb1WoXhVlQmGn7RACq9X37apXIP79K7WZnZzoROKLSErWNoS36myvj4QszFpyfXrrwBn2nWZOcM+Bnx+0L+8yEGxepfu9PRsKpnfaLxRlYe8ZQTcr3NWgepaaXbvqvXcd5dCJwROXF5qw1ddQcBhcY96qqrnhkjbrT1DOKqpaTBL7N7BvDZaPUnaaLu1p4tn+JfeNxZOX98/ZoEZTzbaqml3a9FtY+A0vvs9t6Hp0IHFF5jSFbzBiqSb9Z6k4zMx4+uQISttk7IseUa+cWQblmXWDW79CiO3x/K2x4Q88oqoy9WwTl3Lzgus9g+OOw+2s1JTgvzeZh6ETgiMqnjlqi2JwlXDZKDSK7+6i9WvcstHdEjqUoB4pzbTdjqCa+oXDLMrXj1roX4YeZUJxv76gcS+ZJcPdVs6/szWCAy/8F18+D07vgk8ttvmWpTgSOKM3CxeYsIbSjutMMH6AWxaz+tx5ELpdj46mj5nDzhMkfwajn1KKzz52vNLJVZcbbZg1BbURdV7ZyvAQ+Ha22mrURnQgc0dkj9hsoro53ENz0E/SbDX++B99cr8tSgGMmAlAfckMehBkL4exxmDtCl6UoZ4+po+Zo1RtmrYOmneG7m9TGRDYoCqkTgSPKiIOgCHtHUTkXNxj/Bkx4F078ocYNnKjuulWUJwJfB0sE5TqOVV17Hv7wxdV6C0wpyxKBnccHquLfQpWLL9+YaNHNUJht1UvqROBojMWQmwJNarGRtj30vlX9sRbnqemKjbnmja1XFddF006qay/yCrUF5rL7G2/F2YIMKMp2zBZBOTdPtTHRmFfg0Cr4dOT5sUMr0InA0eScBiT4t7J3JDUL7w+zN6hSu9/fBmvmNM4yB/ZYVVwXXgGqm2joI7DzC6fcUtEizs0YcuBEAKprb+A9qpRIfjrMvdxqpV+smgiEEGOFEIeEEEeFEE9U8nq4EGKdEGKXEGKvEGK8NeNxClmn1PcmTpAI4Hwzts9M2PwOfD0ZclPtHZVt2WtVcV0YXGDk06oiZtoR+HgYHF9v76hsqzwR2KrOUH21HQp3blBTgz39rXIJqyUCIYQL8D4wDugCzBBCdLnosH8Di6SUPYHpwAfWisdpZJclAv8w+8ZRG67uqsTBpA/UOoOPh8LJrfaOynbsvaq4LjpfDbPXqV20vpoMG99oPDvVZTrAYrLaahKmikJaahvUi1izRdAPOCqlPC6lLAYWApMuOkYC5SmuCXDaivE4h6yy1bvO0iKoqOeNqjyyqyfMHw9/vd84FjM5wqriughpr8ocdL0Wfn8Rvp2quiAausyTauDc04GmZ5vDii1OayaCVkDFmgSJZc9V9CxwkxAiEVgJ/F9lJxJCzBZC7BBC7EhNbeDdDlmn1B+ou4+9I6mbFt1h9npoP0bth/zdTQ1/iqmjrCquCw9fuO5TuOo/cGKD6ipq6KvHM+LVjCFn6MqzEWsmgsr+L198ezgDmC+lDAPGA18JIS6JSUo5V0rZR0rZJzQ01AqhOpDsU6oZ6My8AmD6NzDmZTj8i/pwSXSMcrsW52iriutCCOh7B8z8BYQBPh+nxnsaaldRxgnnGR+wEWsmgkSg4hzIMC7t+rkdWAQgpfwL8ASsuG+cE8g65RwzhmoiBAy8V/VrSmDelQ3zw8VRF5PVRavecOdG6DhezQBbMM0udW+sqtQI6Sccp3yLg7BmItgOtBdCtBVCuKMGgy8uvH0SGAkghOiMSgQNvO+nBtmJzjk+UJWwPnBXhQ+Xb65vWLOKGlIiANWam/qlqjp7fAN8OAiOrbN3VJaTdRJMJToRXMRqiUBKaQTuA34FDqJmB+0XQjwvhJhYdtjDwCwhxB5gAXCblI1hdLEKxfmqP70htAgq8gpUHy5XvQVxm+DDgQ1nK0xHX1VcF0KoqrOzflfjVV9NVkncATddr7Wzx9X34Ej7xuFgakwEQogpQgi/sv/+txDiJyFEL3NOLqVcKaXsIKWMlFK+VPbcHCnlsrL/PiClHCyl7CGljJZSNpBPhzoqnzrq7GMElREC+t6uBpJ9msK3U9Q+uyUF9o6sfpxhVXFdNY9Sv6/et6puvU9HQuohe0dVP2ePqu+6RXABc1oET0spc4QQQ4AxwBfAh9YNq5EqX+XZ0FoEFZXXyx9wL2ybCx8PV6V3nZWzrCquK3dvmPAOTP9W3ah8PAy2fuy8Yz3px9TUUZ8GPumklsxJBOW1hq8CPpRSLgXcrRdSI5btZKuK68rNE8a+DDcvVrNuPh0F619zzvIUzrSquD46XQV3/wURQ2HVY2oFuTPuWHf2KAS1a/i/r1oyJxGcEkJ8DEwFVgohPMx8n+Nx9D7O8vISDblFUFHkFXDPn9B1Mqx/WXU9nDlg76hqxxlXFdeVXzO48Xu4+r+QsF0NJO/62rkWDZ49qruFKmHOB/pU1IDvWCllJhAEPGrVqKxh7yJV+qD8w9YRZSeq/nNXD3tHYjtegWpB05T5alX1x8NUuQNnaR0466riuhIC+vwD7t4MzaJg6b1qJpgztA6MRSpOnQguUWMikFLmAynAkLKnjID16qFaS5MwlQQ+Hwvpx+0dTeWyTjX8bqGqdJ0M925TNXB+f1Ft13d6t72jqpkzryquj6C2qtjguNch/k/4YKDa58CRxw7STwBSzxiqhDmzhp4BHgeeLHvKDfjamkFZRZtBcOsyKMqFeeMcczOV7AaymKyufEJUy2Dql+oD9pMr1LRFR91vtyGsKq4PgwH63wl3/wmteqp9DuaPd9yZRedmDOlEcDFzuoYmAxOBPAAp5WnAOadItOql9gQFtYz+yFr7xnOxrAZQXsISukyCe7dCz5vUtMUPBjjmuoOGtpisroLawi3LYNL7kPo3fDgYfn/J8aYGpx9T34N0IriYOYmguGyRlwQQQjhpNbQyTTurmip+LeGb6+C35x2jP7owC4pzGneLoCKvQJj4rup+cPVU6w6+u8mxNlLRieA8IVTivne76ubb+LrjJfCzR9W0US8nqzpqA+YkgkVls4YChBCzgLXAJ9YNy8qC2sKs36DXLfDHf+DLiaoioT2Vf8A11jGCqkQMgbs2wcg5qgX3v76w4Q0oKbR3ZOfXPwS2tW8cjsQ3FK77RO2q5eKuEviCG8r65+3s7HHdGqiCOYPFbwI/AD8CHYE5Usr3rB2Y1bl5wcT3YPJcSNpbNtj1qf0Gu7KccEMaW3F1h6EPw33boP1oWPcivN8PDiyz79TFA0ugRQ8IcPD9pe2h3Qi4azOMek7tgPZ+P1jzjNU3Ya+WnjpaJbPWA0gp10gpH5VSPiKlXGPtoGyqxzS45y+1/+6Kh1XrwIqbRFcpW7cIahQQrgaSb1mm9mtYdLPadzdxh+1jyTwJp2JUN4hWOVd3GPIg/F8MRF0Hm/8L7/VWs4tKS2wbS1GOKgeiB4orVWUiEELkCCGyK3zPrvjYlkFaXUBruOkn1UIobx2sfRaK82wXQ9YpVQu+IRUvs5Z2w+HOP1Tpg/TjaiHaolsh9bDtYjiwVH3vco3trums/FvA5I/gjt/Vqt4VD8H7/WH/Etu16NJ1sbnqVJkIpJR+Ukr/Ct/9Kz62ZZA2IYQaM/i/GOg+FTa9rfqjdy8AU2nN76+v7FNqPrqLq/Wv1RC4uELv2+D+XTD8cTi6Fj7oD4vvsk1/9P6ybqEgPT5gtrDeaqLGjIXg4gbf36oWEP690voJQRebq5Y56wjCK/uyRXB24RsK13wAM1erGQZL7lJ/rEfWWvePNStRzxiqCw9fuPxf8MAeGHAP7F+suh9+uhNS/rbONTNPwqkdujVQF0JAx3Fq7cE1H6oum4UzYO4INeZjrZuus+VTR9tZ5/xOzpzbzxUV/tsTaAscArpaJSJHEd4fZq2D/T+pKabfXAetB8DwRyFypGWLVpWWqKZr636WO2dj4xMCY16CgffBn+9BzOew9ztVLG3A3dBmsOV+Z+XdQl11IqgzgwtE3wDdpqjf08Y31JhP8GUw6P+g+3RVnNBSzh5TEzHcvKo8pLi4mGPHjpGf76ALGM3k7e1NZGQk7u7m1wYVtd0HpmwvgjullHfWMj6L6NOnj9yxw8aDg8Zi2PkFbPqvGtRt1Vv9sXaaYJmunN+eV9NYp36pFlNp9ZeXBls+hB2fqc1+mneDfneqwV0P3/qd+5ORaperOzdaJlZNtQQOLFUDykl7wCtIrUvoM7P+3W9Zp1SF29COcMuSKg87ePAgAQEBNGvWDIPBOetqmkwmkpKSSEhIIDw8nJYtW557TQgRI6XsU9n7av3TSil3An3rHqoTcnVXOzbdv0sNUOafhe9vg3d6qLGE3JS6n/v4BvjjLfVHr5OA5fiEwMin4Z8HVLXMUiMsuw/e7KAKpcVtrls3REac7hayBoMLRF0LszfArT9D26Hw1/vwbk+YfzXs/BIKMmt/3qxTMP8qVQrkiqerPTQ/P9+pkwCAwWCgRYsWuLq6snjxYrKzzZvXU2OLQAjxUMXrAL2AYCnlmDpHWw92aRFczFQKh39Rd5xxf4BwUfPbe0yH9leqqY3myEtTy/E9/dVOUOa+T6s9KSFhqyqbvH+x+mDwDlH91R3HqW4/n+Cq319qhJ3zYd3Lql/7vu0QGGGr6Bun7NOw8yvVdZR+DFw8oO0wuGyU+vdW074C5Ukg/6yaFdi6+vvXmJgYevfubeEfwj5iYmLYvHkzEyZMoG1b1aKqrkVgTr9GxbpCRtSYwY/1jtSZGVxU33Onq9SA5J4Fqsz14V/Uaso2g1VCCOurSlpc3BVRlAPHfld3PAUZcNMPOglYmxAQPkB9jX0VjqyGv5er2T+7vlLHBLdX9aj8W4F/S7XrWO4ZVUri2O+qjk7EUBjzsk4CtuDfEkY8DsMfg1M7IfYHOPwr/PK4+vIOUV1+zaPUjDtXT/V19iic3KLWebh6mJUEGiIpJeZ2/dd6jMDeHKJFUBlTKcRvVn+oR9ZAWnkFRqE+NLwCwOCq7kyT90JpsaqnM+YViJ5hz8gbN2ORWpCWsBUStkHyPrXwyFSh/pSbj5p/Pvxxlfz17lb2lX5cJebTu9TvK+Wg+vdUzuAKLaJV0o++UW2PagZHaREkJyfz4IMPsn37djw8PIiIiODZZ5/lgQceIDs7GxcXF5566immTZtW5TliYmLYtGkTEyZMoF07NVOqXi0CIcQaYErZpjQIIQKBhfbqGnJYBhfVbG07TM1eyUxQH/hn9kPKAVX+WpaqhNF3FnQar7oj9LoB+3L1gIjB6qucyQR5qarl5tes4e5H7KyC2l04DbTUCCV5qtppSQH4NnXaFraUksmTJ3PrrbeycOFCAHbv3k1WVhZffvkl7du35/Tp0/Tu3ZsxY8YQEGCZAnrmfAqFlieBskAzhBBNLXL1hiygtfrqdJW9I9Fqy2BQCcCvEe085sxcXMGlCXg2sdgpn/t5PwdOW7aAQpeW/jwzofpZ9+vWrcPNzY277rrr3HPR0dEXHNOyZUuaNm1KamqqxRKBWZvXV1xAJoRoQ1lJak3TNM1yYmNja+ye2rZtG8XFxURGWq5chjktgqeATUKIDWWPhwGzLRaBpmmag6npzt1ekpKSuPnmm/niiy8sOs21xkQgpfylbBHZAEAA/5RSplksAk3TNA2Arl278sMPP1T6WnZ2NldddRUvvvgiAwYMsOh1zak1JICxQC8p5c+AtxBC10LQNE2zsCuuuIKioiI++eT83l/bt29nw4YNTJ48mVtuuYUpU6ZY/LrmtC0+AAYC5XMcc4D3LR6JpmlaIyeEYPHixaxZs4bIyEi6du3Ks88+y8aNG9m4cSPz588nOjqa6Ohodu/ebbHrmjNG0F9K2UsIsQvOzRoyv5qRpmmaZraWLVuyaNGiS55/+unqS2TUhzktghIhhAvnN68PBey0n6OmaZpmaeYkgneBxUBTIcRLwCbgZatGpWmaptmMObOGvhFCxAAjUbOGrpFSHrR6ZJqmaZpNVJkIhBBBFR6mAAsqvialTLdmYJqmaZptVNciiEGNC1RWYUsCNe75JoQYC7wDuACfSilfreSYqcCzZefcI6W8oeawNU3TNEupMhFIKeu1LVDZAPP7wGggEdguhFgmpTxQ4Zj2wJPAYF3DSNM0zT7MWqMshJgohHiz7OtqM8/dDzgqpTwupSwGFgIXb8E1C3hfSpkBIKWsx1ZfmqZpzi85OZnp06cTGRlJly5dGD9+PIcPH8bLy4vo6Gi6dOnCXXfdhclkucmb5qwsfhV4ADhQ9vWAEOIVM87dCkio8Dix7LmKOgAdhBCbhRBbyrqSKoththBihxBiR2pqqhmX1jRNcz7lZahHjBjBsWPHOHDgAC+//DJnzpwhMjKS3bt3s3fvXg4cOMCSJVXvv1xb5iwoGw9ESylNAEKIL4BdqC6d6lQ1tnDx9dsDI4Aw4A8hRFTFstcAUsq5wFxQG9OYEbOmaVrdrXpCbXpjSc27wbhLhkkvUFUZ6ri4uHOPXV1dGTRoEEePHrVYaOaWr6tY9Nrcot+JQOsKj8OA05Ucs1RKWSKlPAEcQiUGTdO0RsecMtT5+fn89ttvdOvWzWLXNadF8AqwSwixDnWXP4yaWwMA24H2Qoi2wClgOnDxjKAlqBpG84UQIaiuouNmxq5pmmYdNdy528OxY8eIjo5GCMGkSZMYN26cxc5tzoKyBUKI9UBfVCJ4XEqZbMb7jEKI+4BfUdNH50kp9wshngd2SCmXlb12pRDiAFAKPCqlPFv3H0fTNM15VVeGunyMwBqq7BoSQvQq/wJaoLpxEoCWZc/VSEq5UkrZQUoZKaV8qey5OWVJAKk8JKXsIqXsJqVcWP8fSdM0zTlVVYY6Pj7eqtetrkWwA9gPlE/TqTj4K4ErrBWUpmlaY1RehvrBBx/k1VdfxdPTk4iICP773/9a9brVJYKHgeuAAtQagMVSylyrRqNpmtbIVVWGOjY21mrXrLJrSEr5tpRyCHAfavbPb0KIRUKIaKtFo2maptlcjdNHy6Z1LgVWo1YLd7B2UJqmaZrtVFd9tB1qyuck1CDxQuAlKWWhjWLTNE3TbKC6MYKjwF5UayAbCAfuUXvZg5TyLatHp2mapllddYngec6XhPC1QSyapmmaHVRXhvpZG8ahaZqm2Ym5tYY0TdM0KxNCcPPNN597bDQaCQ0N5eqrza3+Xzc6EWiapjkIHx8fYmNjKSgoAGDNmjW0anVx9X7Lq27W0ANSyneEEIOllJutHommaZqDeG3ba/yd/rdFz9kpqBOP93u8xuPGjRvHihUruP7661mwYAEzZszgjz/+sGgsF6uuRfCPsu/vWTUCTdM07Zzp06ezcOFCCgsL2bt3L/3797f6NaubNXRQCBEHhAoh9lZ4XqDqxXW3amSapml2Ys6du7V0796duLg4FixYwPjx421yzepmDc0QQjRHlYqeaJNoNE3TNCZOnMgjjzzC+vXrOXvW+pX5q92PoGzfgR5CCHfOl5Y4JKUssXpkmqZpjdTMmTNp0qQJ3bp1Y/369Va/Xo0b0wghhgNfAnGobqHWQohbpZQbrRybpmlaoxQWFsYDDzxgs+uZs1XlW8CVUspDAEKIDsACoPqNNTVN07Rayc29tNL/iBEjGDFihFWva846ArfyJAAgpTwMuFkvJE3TNM2WzGkR7BBCfAZ8Vfb4RiDGeiFpmqZptmROIrgbuBe4HzVGsBH4wJpBaZqmabZTYyKQUhahxgl02WlN07QGSNca0jRNa+R0ItA0TWvkqk0EQggXIcQbtgpG0zStsZJSMmTIEFatWnXuuUWLFjF27FirX7umlcWlQojeQgghpZTVHatpmqbVnRCCjz76iClTpnD55ZdTWlrKU089xS+//GL1a5sza2gXsFQI8T2QV/6klPInq0WlaZpmR8kvv0zRQcuWofbo3Inm//pXtcdERUUxYcIEXnvtNfLy8rjllluIjIy0aByVMRKbPLgAACAASURBVCcRBAFngSsqPCcBnQg0TdMs7JlnnqFXr164u7uzY8cOm1zTnOmj/6jpGE3TtIakpjt3a/Lx8WHatGn4+vri4eFhk2vWOGtICNFBCPGbECK27HF3IcS/rR+apmla42QwGDAYbDep05wrfQI8CZQASCn3AtOtGZSmaZpmO+YkAm8p5baLnjNaIxhN0zTN9swZLE4TQkSiBogRQlwPJFk1Kk3TtEbs2Weften1zGkR3At8DHQSQpwCHgTuMufkQoixQohDQoijQognqjnueiGEFEL0MStqTdM0zWLMmTV0HBglhPABDFLKHHNOLIRwAd4HRgOJwHYhxDIp5YGLjvNDVTbdWtvgNU3TtPozZ9ZQsBDiXeAPYL0Q4h0hRLAZ5+4HHJVSHpdSFgMLgUmVHPcC8DpQWIu4NU3TNAsxp2toIZAKXAdcX/bf35nxvlZAQoXHiWXPnSOE6Am0llIur+5EQojZQogdQogdqampZlxa0zRNM5c5iSBISvmClPJE2deLQIAZ7xOVPHeuXpEQwgC8DTxc04mklHOllH2klH1CQ0PNuLSmaZpmLnMSwTohxHQhhKHsayqwwoz3JQKtKzwOA05XeOwHRKG6m+KAAcAyPWCsaZpmW1UmAiFEjhAiG7gT+BYoLvtaCPzTjHNvB9oLIdoKIdxRi9CWlb8opcySUoZIKSOklBHAFmCilNI2xTU0TdMcSHVlqF1cXIiOjiYqKoopU6aQn59v0WtXmQiklH5SSv+y7wYppWvZl0FK6V/TiaWURuA+4FfgILBISrlfCPG8EGKi5X4ETdM051dehvqhhx6isLCQvLw8nnrqKd5//328vLzYvXs3sbGxuLu789FHH1n02uYsKEMI0R2IqHi8OWWopZQrgZUXPTenimNHmBOLpmmatf2x6DBpCbkWPWdIa1+GTu1Q7THmlKEeOnQoe/futWhsNSYCIcQ8oDuwHzCVPa3LUGuapllBdWWojUYjq1atsviuZea0CAZIKbtY9KqapmkOrKY7d2uqrAx1QUEB0dHRKrahQ7n99tstek1zEsFfQoguF68I1jRN06zj4jLU5WME1mJOIvgClQySgSLU+gAppexutag0TdM0mzEnEcwDbgb2cX6MQNM0TWsgzEkEJ6WUy2o+TNM0TbOEi8tQ5+ZadgbTxcxJBH8LIb4FfkZ1DQHmTR/VNE3THJ85icALlQCurPCcnj6qaZrWQJizH8E/bBGIpmmaZh/mLCj7nApVQ8tJKWdaJSJN0zTNpszpGqq4V4AnMJkLq4hqmqZpTsycrqEfKz4WQiwA1lotIidXkJtDyoljpMYdpyAnm+LCAkoKi/APDaVZu/Y0j2yPT0CgvcPUtAZDSokQlW1/opnLrKJzF2kPhFs6EGdlLCkhcf9eju3cTtzuGDLPJJ17zeDigruXN67u7uRmpINUPWydBg/nytn/h5unp73CbvRMpaWkxB0n8WAsKSeOkZVyhuzUMxTkZOPm4Ymblxe+gUF0GjyczkNG4OVXY8Fdzcoykk+TfOQQycePknLiGHmZ6RTk5lKUl4unrx9BLVsR2CKMjgOHENGjl73DtagRI0Zw/Phx4uPjzyW9a665hrVr11pkaqk5YwQ5XDhGkAw8Xu8rOzFpMnHq7wPEbljLka2bKS4owNXdg/BuPeg2cgzN2l5G07bt8PT1O/dLKy4sICXuOMd3bmf7sh9JS4hn0sNPEdC8hZ1/msajMC+X4zHbOLLtT+L37aGksAAAv+BQApo1p033Xnj5+2MsLqaksIC0hHjWzZ/Lxq/n0XHQMEbPug9Xd3c7/xSNh8lUSkLsPo7t3MqJXTvITFY3Wa7uHoRGtKVp28vw9PXD08eH/KxM0k+f4uj2v4hdt5pOg4cz4pY7GlTrOyAggM2bNzNkyBAyMzNJSkqq+U1mMqdryM9iV3NyhXm5xK5bw+7VK8g6k4ybpxcdBw6hff9BtO7aHTd3jyrf6+7pRVinroR16kp41+6seOd1vv7Xg1z7xLO07NDZhj9F41JqNHJidwyx61ZzYlcMplIjvkHBdBk6grAu3Qjr1BXfoOAq358Sd5x9v//K7l9X4Obpxajb77Zh9I1TStxxDmz8nb83byAvMwNXN3daR3Wn1/hJhHWOIrhVawwuLpW+11hSwrYl37NtySLidscw7v8epl3PvrWOYd38uaTEH6/vj3KBpm3acflts6s9Ji4ujrFjx9K/f3927dpFhw4d+PLLLwGYPn06CxcuZMiQIfz0009ce+217N+/3yKxVZkIhBDVdv9IKU9aJAInkJ2awvaff2L/+rWUFBXSqlNXBk25kfZ9B9apeyeiRy9uevW/LHj6UbYt/YFrHn3aClE3brnpZ9m9eiWx61aTl5mBT0AgPcdNoEP/wbS4rAPCYM4urdA0oh0jZ96Ni5s7McsXEx7VnQ79B1s5+sbHWFLCkS2b2PXrcpKOHMLg4kq7Xn3oPPRy2kb3xs3DvH9nrm5uDJpyAx0HDmXZWy/z++cf07ZHb7N/347g0KFDfPbZZwwePJiZM2fywQcfADBy5EhmzZpFaWkpCxcuZO7cubzwwgsWuWZ1LYIVqC6hiqMwEggFmgKVp+QGJPNMMtuWLGL/ht8AQafBw+g1biLN2l1W73M3adqcyN79+fvPDZhKS6u8w9Fq58yJY8QsX8yhv/7AZDLRrldful0xhrbRvXFxrcuQmDJ0xi2c+ns/qz96l2ZtI2nStLkFo268CnNz2b16Bbt++Zn8rEwCW7RkxC2z6DLs8nqNywSHtWbAtdNY+d6bxMfuIaJ7z1q9v6Y7d2tq3bo1gwerm42bbrqJd999FwAXFxeGDBnCd999R0FBARERERa7ZpX/MqSU3So+FkJEoMYGRgEvWywCB5Sbkc6WHxew97dfMbi40H3UOPpOvA7/kFCLXie8WzR7f/uFpKOHadVRdw/Vx6lDB9m6+DtO7NqBu5cX0WOupufYCQQ0s8wHtourG1c/8BhfPf4Ay995nRkvvIHBoJN3XeVlZrD955/Yu/YXSgoLaBvdm17jJ9GmW7TF7t7b9x+M5/y57F27qtaJwJ4ungFV8fH06dOZPHnyJbWI6sucweL2wFNAf+A/wP1SyhKLRuEgigsL2L7sR3YsX4zJaKTH6HH0v2ZqtX3I9REe1R2E4OS+3ToR1FHSkUNsWvglJ2P34OXnz5DptxA95io8vH0sfq0mTZszZMat/PbZB6SfSiSkdRuLX6Ohy8/KZNuyH9mzeiWlxhI6DRpG34nXEdqmrcWv5ermRtfhI9m1atm57kFncPLkSf766y8GDhzIggULGDJkCD///DOgNqV58sknmTFjhkWvWd0YQRQqAXQFXgdul1KWWvTqDkJKyaE/N7Lh63nkpp+lw8ChDJl+M4HNW1r1ul5+/jRrexnx+3Yx8HrL/mIbutSTcWz+7iuO7diKl38Tht98Oz1GjbP6lNyQcPXhn5uRrhNBLRTl57Nj+U/ELF+CsbiYzkNHMOC66Vb/N9Z95Fhili8mdt0a+k+eatVrWUrnzp354osvuPPOO2nfvj133333uUQghOCRRx6x+DWraxHsARJQYwX9gH4VmyhSyvstHo0dnE08ydrPPiDxQCxN20Yy4Z9P2HQWT5tuPdixfDHFBfm4e3nb7LrOKjf9LJsXfU3s+rV4eHkzZPot9Bw3AXdPL5tc3zdQtQ7zMtJtcj1nV2o0smfNKrb8tJCC7Cw6DBzK4Kk3EtQyzCbXD2rZitZdurHv91/pN+l6pxg0NhgMfPTRRxc8t379+kqPtVR56uoSQYOuJVRqLGHbkh/Yuvg73Dy9GD3rPqKuGG3zft823XuybekPJByIJbJ3P5te25mUFBWyfdmPbP/5J2RpKb2vuob+k6fi5Wvb2c0+gap7ITf9rE2v64yO79rO+i8/I+N0Iq27dmfYDbfR/DLb7wXcfdRYVrz7BvF7dxER3dvm13cG1Q0Wf2HLQGzpzPGj/PLB26QlxNNp8HAuv3UW3k0C7BJLyw6dcXX3IH7fLp0IKiGl5MjWzaz/6jNy0lLpMHAoQ2fcarFB4Npyc/fAw8dHrRTXKpV++hTrv5jLid0xBLZoyTWPzaFdr752KwNxWb9BePn5s3/j7w6fCCIiIoiNjbX5des+n84JlRqNbF28iC0/LcQnIJBrHptj9w9fV3d3WnXqQvxe621M7azOnkrg93kfcjJ2L6Ft2jL+3ocJ6xJl77DwDQzWXUOVKCkqZOviRWxf9hOu7u4Mv/l2eo69GhdXN7vG5ermRkjrNma14kwm0wWbxjsjk6n2Owo3mkRw9lQCq/73FmeOH6Hz0Mu54rY78fT1tXdYgOoe2vj1PHLS0/ALCrF3OHZXUlzEtsWL2Lb0R9w8PRg58266jx7rMNM1fQKDyM3QXUMVHYvZxu+ff0R2agpdhl7OsJtmOtQsHXdvH7JTkqs9xtvbm+TkZJo3b+60ycBkMpGcnExJSQlCCLNbYY0mEcTtjiEr9QwT/vkEHQYMsXc4F2jTLRqAk/v20HX4SDtHY18nY/ewZu7/yDyT5JAfKAC+gUEkHrR9890R5Waks27+XA5v2URwWDjTnnnVIVptF/Pw9qaoIL/aYyIjI/n77785ffq0U1czLSkp4ejRo0gp8fc3b1GeOesIOgAfAs2klFFCiO7ARCnli/UL17Z6jZtI5yEj7DYWUJ3Q8Ag8ff04fehgo00Ehbm5bPj6M2LXrSGwRUumPP0S4VE97B1WpXwDg8hNT2/U5Y+llOz7fTUbv56HsaSYwdNupu/Ea+3eDVQVDx8fivLzqj3G3d2dbt26ERMTw6ZNm2wUmXUYDAauvPJKgoPNWwNlTovgE+BR4GMAKeXess3snSoRCIPBIZMAqNh8AgIpyM22dyh2cXT7FtZ++j752Vn0m3Q9A66fUW0BP3vzCQzGVGqkICcbb/8m9g7H5rJSzrB67nuc3LebsC5RjJ71fwS1bGXvsKrl4e1DUX4+0mSqdgqpEII+ffrQs2dPjEajDSO0LDc3t1p1b5mTCLyllNsuuvNx3v9DDsrD24eivOrvWBqagpxsfv/8Y/7evIHQiHZMfuJZmrWNtHdYNfItm0Kal5HeqBKBlJK9a1ex4at5IASj7riH7iPHOsXcfA8vb5CS4sJCPLxrXq/j4uKCSyOq/2VOIkgTQkRStieBEOJ6wHKFsDVANV3zszLtHYbNHN+5ndUfv0tBTjYDr7+B/pOn1qsonC35lC0qy81It0ppBEeUnZbK6o/fJX7vLtp078mVs/8P/9Cm9g7LbB4+amJIUX6eWYmgsTHnX969wFygkxDiFHACuMmqUTVCHt4+ZJw+Ze8wrK64sID1X37Kvt9+JSQ8gmuffI6mEe3sHVat+AYGATSKmUNSSg5uWs/v8z6itNTIyNvvocfocU43NlL+4a/GCSxbPLIhMGdjmuPAKCGED2CQUuZYP6zGx8PHl8IaBrOc3enDf7Pyf2+SlXKGvpOuZ9CUG3F1c8zBxer4lCWCvPSGvZagMDeXtZ++z6G//qBlxy6Mu+efTrujnntZEcKaBowbK3NmDXkA1wERgGv5nYCU8nmrRtbIePr4UJyf1yBnophMpWxdvIi/fliAX3AI0555hbDOjjfF0Fyubm54+vk36NXFCfv3svL9t8jPzGDI9FvoO+k6h1nHURee5YmgkY3DmcucrqGlQBYQAxTV5uRCiLHAO6hNbD6VUr560esPAXegBp9TgZlSyvjaXKOh8PD2wVRairGoqEFtap+dlsLK997k1N8H6DR4OKPuuMcqJaJtzTcwiLzMhpcISo1G/vphAVuXLCKweQtmPP+GXeoDWVp5i6BYtwgqZU4iCJNSjq3tiYUQLsD7wGggEdguhFgmpTxQ4bBdQB8pZb4Q4m5Uuetptb1WQ1D+4ViYn9tgEsGRbX/y60fvIE0mxt33MF2GXm7vkCzGNzCowbUIslLOsOK9N0g6/DdRl4/m8ttm26yqq7WdHyOoflFZY2VOIvhTCNFNSrmvlufuBxwtG2NACLEQmAScSwRSynUVjt9CIx6E9vA533R19jITxuJi1n/1GXtWr6BZu/Zc/cBjTtu3XBWfgCDSEhvOtt3nk7bkqgceo9OgYfYOyaI89BhBtarbmCYWMJUd8w8hxHFU15AApJSyew3nboXaz6BcImqXs6rcDqyqIpbZwGyA8PDwGi7rnDwaSB9mRvJplr/9Gilxx+h91TUMveFWh11tWh++QUHkZaTXuEDJ0RlLStj49Tx2/fIzzSPbc9UDj9utsqs1ubq74+LmphNBFaprEbQCoutx7spGPGWlBwpxE9AHGF7Z61LKuagprPTp06fSczi7cy0CJ/5DPbx1M79++A4Gg8EhKrtak09gENJkIj87y+FqIZkrK+UMP7/9KmeOH6HX+EkMu/G2Bpm0y6nVxXm8HfM2Ukoe6vOQvUOqFWtOJKkuEZyo58BtItC6wuMw4PTFBwkhRqG2xBwupazVYHRD4uFdtuAlzzI7DtlSqdHIxm8+Z+fKpTS/rAMTHnzCqRYb1cX5tQTpTpkIju/czqr//QeTycTEh/9F+36D7B2S1ZWv3t+YuIXi0mKnSgS56WdZ+b//MPSGW2lxWUeLn7+6RNC0bFZPpaSUb9Vw7u1AeyFEW+AUMB24oeIBQoieqBpGY6WUKeaF3DB5nmsRONdgVk56Gsv/+zqnDx2g59gJDL95ZoO+qyx3wZaVTlAWo5zJVMpf33/Llp++IzSiHRP++YTV9w12FOUVSFMLUskrzsNoMuJqcPzV7AkH9rH8v69RUlhIXkaGVa5R3f8FF8CXyrt4aiSlNAoh7gN+LTvXPCnlfiHE88AOKeUy4I2ya3xf1uQ5KaWcWJfrOTtnXPBS/gdaXFjA+PsfpfPgSnv2GiQfJ1xdXJCTzcr33iRuz06iLh/NFTPvcujifpbm7u1DYV4uWUVZACTlJdHar3UN77IfKSU7li/mj2/nE9C8JVPnvExwmHXGSKtLBEn1XTQmpVwJrLzouTkV/ntUfc7fkLi6ueHq5n7uD9XdxR0vV8ecuielZNeqZaz/6jMCmrVgytMvEdK6jb3DsimfAFXJ1lp3aJZ25sQxlv3nZfIyzjJ69n10H1nrGeFOz9Pbh8y0M+ceJ+QkOGwiKCks5JeP3uHwX3/Qvt8gxtz9oFVrJFWXCBrW8lYnUF4zffaa2XQO6syzg561d0iXKCkuYs3c/3Hwj3VE9unPuHsfahALxGrLxdUNL/8mTtEiOLhpPas/ehdPf3+mPfeaVfqYnYG7tw9F+efH4BJzEu0YTdUyk5NY+p+XOJtwkqE33EbfiddZvdpAdYmgce6QYkce3j4U5uZyOOMwUjre5KjstFSWvvkiKXHHGTT1RgZMnubUUyfr6+JFZY5WHsRUWsofC75gx88/EdY5ign/fMJh9+SwBQ9vb0oKCs89TshJIDutAGOJiaAWjnEzE7dnJ8vfeQ2B4NonnyWiRy+bXLfKRCClbFjLJp2Ah48P2TnpGIONxGXHUWosZf3Xh2jSzJs+4yLsGtupvw+w7K2XMRYXcc2j/yayd3VLQhoH38Cgc5vYpyXmsOStXQy/oSPt+zSzc2SqYNzyd14jfu8uosdcxYhbZjlNmW9r8fDxwVRcgjBBoHcQCTkJ/PrpfnLTC7nlpUG4uNnvpkZKyc6VS9nw1TyCW4cz6ZF/23Q9R+P+y3AwHj6+pKWrrR4KjAWsWbCPY1vScXU30H1EGO5e9vl17Vu3mrWffECTpk2ZOucVgsMcs1/V1nwCg0mJPwFA0tEsivKNrPlsP6ZSScf+9luUlX46kSWvP09WSgpX3nk/3a640m6xOJLyLkwvkztRIVGkJmeSEqd2BTy0LZkug+0ze8pYXMzaT99n/4bfaN9vEGPv/afNS3voROBAPLx9KExQfZhdk4dw7EQ6bboFE7/vLIe3nyFqmG23AzSZStn49TxiViylTfeeXP3g43iWbfChqdXF+ZmZmEpLyUjKw83ThaZt/Fk7/wCmUhOdB9n+gyVudwzL33kdg6srU+a8RFinrjaPwVGVJ4IWLiGE+4VTsF2N7/iHerFr9Uk6D2yBMNi2ay8vM4Ol/3mJpMN/M/D6Gxh43XS7dLc23g5eB+Tp44OxoICwzI4MPnEtbm2LGH93d4Jb+XBw8yVr8ayqKD+PJa89T8yKpfQcN4Frn3hWJ4GL+AYGIaWJ/KxM0pPzCWzuw1X3dqd15yB+//Jvko9n2SwWKSU7V/3MT68+h39oU256+W2dBC5SnghCXQJp7deaNilRBLfxZsCkdmSeyefE3jSbxpMaf4JvnnqI1LgTTPjnEwyacoPdxtx0InAg7t4+yEIjI+Knku2VRsaQWAwGQZchLUmJzyE1wTZ7AmWlJLPg6UeJ37eb0bPu44rb7sTQiPZvNde5LSvTz5KRnEdQC2/c3F0YOzsKVw8XDv5pmx1dS41GfvvsQ9bN/5h2vfsx/fnXG/zK7roon34ZbAigaXEYIflh+HWWRPYMxT/Ek52/xttsksaxmG0smPMY0mRi+nOv0WHAEJtctyo6ETgQD28fhEnindeEzLAE4gqOA9ChX3NcXA0c2GT9VsHpwwf55qmHyc04y3X/ep7uoxrffHNznRbqDnL9vjXkZxUT2Fzdcbp7uhIZHcrRmBSMJaVWjaEoP4/Frz3HnjUr6TvxOiY9/K8GUzra0soXbQbgizzqh8REUZtUDC4GokeFc+ZENknHrNuKKx8UXvrGiwS1bMWNL71Fs3aXWfWa5tCJwIGUN12FLMInxJX4bFXqydPHjcheoRzedoaSYut9sPz950YWPf8vPLy8ueHF/xAe1cNq12oIjpWqeehrtv8CcMEUxA79m1FcYCQ+1nrrDLJTU1jw9KMk7N/LlXfdz7Ab/9Gop/PWRLqrVq2/9CI1tohk/xOcRpUS7zSoBZ4+buxabb3S4qbSUn6b9xHrvviEyD79mfbMq/gGBVvterWh/2ocSKm7+i5lESEt/DmVe4ri0mIAugxpSXGBkWM7LV+SSUrJtqU/sOKd12ke2Z4ZL75JUMswi1+noUkwJlNqkEQY1aDwXwXrz70W1ikIb393Dm1Jtsq1k48e5punHiI3XbXcul2uZwbVJNdQAIBXjhcZSfmktThOQo6qlO/m7kKXoS2J35dGXpbla18WF+Sz5PXn2bN6BX0mXMvEh550qA2odCJwIFmirOCcLCI8vDkmaTq3+rFl+wCaNPWyePeQqbSUtZ+8zx/fzqfT4OFc/9SLePs3seg1GqqTOQkYvV0IMzXFZCjllYPPcyj9EAAGg6B9v2bEx56lMK/Eotc9umMr3z33JK7uHsx44U3dcjNThklNFTUlCYQAInMuWF3csX9zpIQj289UcYa6yU0/y3fPPknc3l2MuuNeht800+Fabo4VTSOXLlX/pJuXkXYhEQCcyFbz1IUQdBrYgqSjWWSlFljkesWFBSx5/Xn2/vYL/SdPY/x9D+Pq7m6RczcG8TnxGPy9KMjKILC5D1JIdpzZce71jv2aYyqVHI2xXCtu968rWPbmS4S0DueGF9/UazpqIa04nRIXE0VpxTSPbEKL0FBO5pzvCgpq4UNouB+HtlquFZeWEM+3/36EjKRTTH5sDj1Gj7PYuS1JJwIHklKqCpj5+EMbf1XErXycANQdCwKL/KHmZWbw3bNPELd3F6Nn38eQ6Tc73F2KIyswFpCcl4xnYBOK8jNp2qoJTb2aEpsWe+6YkNa+BLbw4fC2+v++pMnExm/n89u8D2nbqw9T57zilPsg2FNKfgrFriaKcwtoFuFPuF84mUWZ5BSfn43XsX9z0hJyOXu6/vuCJBzYx8JnHsNkKmXac6/Rtmefep/TWvS/fAdyukTdOXr5SPzc/Qj2DCYuK+7c635BnoR1DOTQlqR6TXNLP53It/9+hPTTiVzz6NONshJlfZX3LTcJbI7JmENAM0+iQqLYl3Z+a28hBB37NyPpaBbZaXVvxZUaS1j1/ltsX/oDPUaPZ9IjTzlU/7KzSM1PRbq6Ik1FhLT2O1d5tPx3CdC+bzOEQXB4a/26hw799Qc/vvQ0Pk0CueGFN2nm4HtW6ETgQM7kZALg5qlmBkU0iSAuO+6CYzoNaE52WiFJR+s2zS3pyCEWzHmMkqJCps15hXa9+tYr5sbqZLbqUgjwbgmU4tPERLfQbsRnx5+rdw9q6q8QcKCOCwKLC/JZ/NrzHNy0niHTb2Hk7XdjMOg1HXWRUpCCweAFsojQ8MoTgbe/O607B3F4WzLSVLebrZ2rlrH8nddpflkHpr/whlOs6dCJwIFknS0GXHB1VYOLEf4RF3QNAbTr2RRXDxcOban9YqXju7az6IV/4eHtzYwX3qD5ZR0sEXajVP578XVrAYCbeyHdQroBsD9t/7nj/II8aRMVzIHNSZQaTbW6huq+e5KTsXsYc/eD9J881aGqmzqb1PxU3IQPUERAM2/C/NTMuIqJAKDjgGbkZhRx6khmrc4vpeSPb+ezbv5cLuszgOufehEvXz9LhW9VOhE4iAJjAWR7gPBAGNSU0Qj/CNIL0y+4w3TzcOGynmWLlWqxpmD/ht9Y8voLBLUMY8bzbzSa7Qmt5WTOSYI8gxDGshlWMocuwV0QiAu6hwC6DmtFQXYxJ/aYX8Ig80wyC+c8RnpSIpMfm0PUCL2HU32lFaThJv0wuJRgMAh83HwI9gy+JBG07RGKm6dLrcbiSo1Gfv3wHbYt/YHuo8Yy4aEnnGrihU4EDiIxJ5HAgmYIgwcmo5rHXNmAMUDHgS0oLizl+J5Us8694+ef+OWDt2ndpRvTntGDjJYQlxVHG/825Oeo/ZnzstLxc/ejbZO2lySC8K7B+AV5ErvRvI1QUuKOs3DOoxTm5TL16ZcdepDRmaTkpWIw+oIsPvdca7/Wl/z7cnN3IbJXU47GpJg19bekwy1d9AAAGaJJREFUuIhl/3mJ/RvWMvD6Gxh1x71O132nE4GDSMhJIKCgKS5unhQVqH2LI5pEAFwyTtCqfQB+QZ4c3Fx995CUko3fzmfD1/Po0H8wk594Fncv621315iczDlJG/82ZJ81gHAh56y62+8W0o19afsuGMw3GARdh7Xk1KFMMpKr35M68WAsi557EuHiwvTnXqdF+8a5m5il5ZXk4ZrniTB5Umo8vzlNVEgU+1L3kV+Sf8HxPa5ojbGolNgN1SfvwrxcfnxpDsd37WDk7feownFO2H2nE4GDSMxJJKCgKe5e3hTlqalrYb5huAiXC2YOAQiDIGp4KxL/ziD5ROWDxiZTKWvmvsf2sqbqVQ8+hqubm7V/jEYhrySPtII0wn3akJ1WiId3E3IrJIL0wnSS8i5M0p0HtcTgIojdeKrK8x7fuZ0fX5qDd0AgM55/Xa8RsKCU/BRC8sIQwgNpMmIsVq2C4a2HU2wqZkvSlguODwnzJbxrMHvXJVbZBZuXmcGi554k6cghrn7gMaKvHG/1n8NadCJwEAkZifgVBeHt50dRnrprdHNxI8wvjKOZRy85Pmp4Kzx93Ni+PO6S10qNJax45w32/b6a/pOnOWVT1ZGVzxhqXhKONEl8AoLJSVeJICo0CoC9aXsveI+3vzuRPUM5tCW50npRBzetZ+mbLxLcOpzpz72Gf4jjzzRxJqn5qYTkhYHBA1DF+gB6N+uNn5sf6xPWX/KeXmPCKcgpqbSKbHZqCgufeYyM5NNMfnwOHQcOtWr81qYTgYNIOZOFwIBPgP+5P1KAAS0G8OfpP8ktvnCBi7unK9GjW3Ny/1nOnMg+93xJYSFLXn+Bw1s2MfymmWqhmBM2VR1ZfI7qU/ZOCwGgSdPQc11DHQI64G5wJzY19pL3RQ1vRVG+8ZIyIbtXr2Tl//5Dy46dmfL0y7rEhxWkFKgWgVeA6hot/zfmZnBjcKvBbEjcgEleOKurZfsAmrX1Z/fak5hKz7929lQCC555jIKcbKb8+0Wb7StsTToROAApJdlnVL+lb5D/uRYBwNXtrqaotIg18WsueV+3EWF4+LiyfaUqQ1GUn8ePr8whfu9urrzzfvpMuNY2P0AjE58VDxLObC8muJUvwWHNyT2bhpQSNxc3OgV3umTAGKDFZQGEdQpk+/ITFOSqromtS77nt88+oF3PPlz75HPnauZrlpWan0pwXisCm6uJEhVvtka0HkF6YfolvzMhBL2uVN1/x3apiRlnThzju2cex2Q0Mu2ZV2nZobPtfggr0onAAexM2YkpQ+0a6h8SgLGkGGOJmq3QI7QH4X7hLD++/JL3uXu6Ej0qnPh9Z4nfn8ii5/9F0pFDXPXAo3qfWis6mXOSLkV9yTidT/crwvALDqXUaKQgR7XMuod052D6QYwm4wXvE0IwdGoHigtL2bLkGH8s+IJNC76g46BhTHz4KdzcPezx4zQKqWnp+JQ0oWmY6nKreLM1pNUQXIQLGxI2XPK+tj1CCGjmTcyqeBIOqIF8Vw8Ppj/3GqFt2tosfmvTicABLPh7AaFFrfAJdMfbXy1AKS67YxFCcHW7q9mWvI2k3Ev7KruPCMPNo4BlbzzN2cQEJj36b6fvr3R08dnxdEsehqePGx36NsMvWNWUL+8eigqJosBYwOGMw5e8N6ilD92Gt2TP6m/YtuR7uo0cw/j/exgXV719uDXlJKmkHBIeBEBR/vlZQk08mtCrWS/WJay75H3CIOg/sR0pcfv54aU5+AQEMv251wlsYdv9w61NJwI7O5N3hrXxa2lT2oGgFr54+KjNTQov6h4C/r+9M4+So6r3+OdXvU7PTM+efUL2FRKWQGIACWsAWUQUiMpyQFARl/c8T+DxDgf1uIuKwkNBzANZ1ABKiKyGiKyBhC0bZM8wySQzmX3rrer3/qjqmelZQkhmMmb6fs6pc6uqb1f9bn/71q/qd2/dy9+3/b3H99ub95JqW0KivYFZZ17PhGPMkBEDTe2eJgqrypl58ij8QR/5xW5bQdoRzB05l5AvxIPrH+zxXce2aar+O3b8HaLD55uG/EOAqtJW5cb4h48vAzJDQwALxixgc8PmjGGpOw+wjVTr31CNctrV/0O0tGzAbT7UGEcwyCzZuISjKk/Fqo0wekphxyxl8bbOxuHyaDnHDDuGJ7c8mdE/vbaygj/d+h0cO8b4Odfy/kofuzZ/vNfiDR+PpkQTYypmIeI2/gLkexeGdBfS0pxSFk1bxLKty9jSsKXju25vrp/y/ssvMGnuBcTjc1mzou/upIb+4YktTyA1OVgFKaIlbkN8D0dQvgCAFyszw0PrX1rB0tt/wLDxEyga80VefrSSZHxgpx8dDIwjGEQSdoK3VmxjXsX5TD5+OMeedUSnI2jN/KOeN+E8tjZuZX3degD2bN3Mn2+7CcdxuPS2H3PeDWcRLQ3z7L1raWtK9DiXoX/YVrOd6dXzyJ8GeUXuCKCRaAGWz9fRhRTg6iOvJhKIcNc7dwHu3A9//cn32LjyFU65/Bou+Na1TDi6jJeXbOrX8e8NmexqrOKFhzYwvn4WM44eSyCcg4jVEXpNMzY6lgkFE3h2+7MdN1tvP7uMp++8nfIZR3LJrT/grGuOo7GmnZeXbDpkk9wfKowjGCTUUR5bupxjN55DdLLF6VdNRyzpCA11v2NZOG4hASvAXzf9le3vvuXOUJVutBo7jmCOn7OvO4p4W4qnf7umo1eKoX9wHGXd6x+y4o7thOwIM08d2fGZWBZ5xSUdoSGAonARl8+4nOd3PM+7O1ax5Pu3dAweN+e8ixBLOOtLMxk9tZDl929g23v7Pw6RYf9ormvnwZ+8xNSquUz8ZCEnXzIFESEYycloI0hz6dRLebv6bR7e8DCvPfoIL/zht0ycM4+LbnTfyB89tYhjzxrL+pd38c+HPsjoUnq4Y1qoDiGOo9TtamHTqmrWvVZJvDFEXdGHXHfD5/H5XJ/c4Qi6PREUhAo4f+L5rFq+jJy1b1A2ZiwX3XxbR3wa3Lchz7hqBv9YvJ6//PBNzvnyUQw7InroCjjEsJMOe7Y3sXNjPetX7qSlOkFTThOlCxs5btZpGXnziks7QkNprphxBU+s+hNP/fAH5LTCBd++hUlz5nZ87g/4OPers3jil2/z7L1rOf2K6UyaM8y893GAqCrNtTEq1tWy5Z0aKj+oI0g+kXPrOfuCTr1CkbyOt/e7smjaIl758GVeue8+dldGmHnK6Zx53dczGvLnXTgREWH1MztoqY+z8NqZBMOH/2X08C/BvyGOo7Q1xmmsbqehuo36PW3U7GimpqKZZNxGxaGiYAM7pr3HtRcsIhjqHPohHMkDej4ROLbNaRXjib5XSlVpO/OvvzTDCaSZdNwwoqVhnv7tGh7/2Vuc9LlJTD9pVIejMWRi2w5tjQla6mI018doqolRv7uV+t1t1FW1YicdFKUubxfrp7/EVy76PAvGLuhxnPySUvZs3ZSxb++6D1j4UilxO8awK87PcAJpgmE/5319NsvufI/n7lvHB2/s5pRFU8kvNhPP9IY6Sqw1SWtjnKa9MZprYzTubaduVwt7P2wh3ub2DmqN1LNxxGqCM9u48/xfZBwjFIkQb+/5RBBvbWX+K7nsroyw9Ujly1+6tkdvLrGEeZ+eSH5JmBcf2chjP13NiRdPonxG8WHtwI0j+JgkEzatDXHaGuO0Nibci0hDnFZvaa6L0VIfo+tLiuqzacmvpaJ4I3tyt9NQtpPPHv1pbpn2cwpCmW+R+kMhLJ8vwxHU7drJM3f9gqrNHzD55E/y0siX+fZrN7K4YDFTi3sOSjbsiCiX/PfxPHffOl58ZCOrn9nBrFPLmXHyKEI5Q09yVSUZt0nGbBKxFIl2m0R7inh7ikR7ilhbknhbinhrklhLklhrkvaWJG2NiV5Hl7RzY7REatk5bBsVeRuoju5gVvlMfjbvVsqjvY//k19SypY3X0dVcWybN55YwqtLHqakfCyvn9DCI1V3s/X1Om48/kYCvswxn3Lyglz8X8ey5p87eX3pVh7+7kpmnjiKKXOHUzY2/7C+wPSGqmKnHE+vtGapDs3SWsXbXO1iza5e7c0J2poSOHZmfF4DNu35Dewp3MGHwzdTFd3CEWNHcPHUizl73NlYknkTlF9axo417/DW009yzMJPIZbFpjde5YXFv6O9qZHpl1/MA3W/4qaXbuK2+bdRmtPzhmvmyaPJKw7zz4fe58nfvMuICVGOO2cc5dOL8fkPv5suGchGDxE5G7gD8AG/V9Ufd/s8BDwAHAfUApeq6vZ9HXPOnDm6atWqfWU5YFSVpvo2du+upba6mfqaFlrrE7Q1JEk0OaSaQeM9RXasFLFQCy3BBhoDe2kO1dMSrKcxXENjTg3kpphaMoVZpbM4quwo5o+aT24gt087/vdLn2f0tBlMmXcSNRXbefuZJ/H5/ZxxzfVMO/EUqlqq+OJTX6QuXse548/lihlX9OoQ1FF2rKvlnX9UsPODBiyfMHxclFGTCxkxoYBoWQ7RkjD+4MB0X0xXeDvpkEr2TFNJGzvhrSfszjTRmSYSKeKxJMlEikQ85V7wEzapuE0q7mAnFGc/mkNUHOxggkQgRjzQSpu/mSZfHS3+BtqCTR2atYTqyY/kMq5gHJMKJzF/1HzmjpxLfnDfE4y89dQTrLj/XuZf8gXWLH+O5toapp+0gDOvuwEJ+Pn1W79m8brFzC6bzfWzr2fOiDkEfT3Hq2+qbee1x7ew9d0anJRSODxC+fRiSsvzKCvPJ1qWQzDsGxDnoKo4jmJ/lFZJVxs76ZBM2B1pKuGQjKdIxJMk4l30ittddFXsuMJ+hNftQIJUIEE80Ep7oIUWXyONvlrago0dmjWF9mIHEowrHMe04mlMK57G/FHzmVw0uc/jNtfu5bl7fsP2d1YzauoMcvLz2bJqJWVjx3HWl7/BiElTeHD9g9y+6nYCvgBXzrySq2Ze1WudtVMOG16tYvXT22mpjxMI+Rg9tYjy6UUUjcilcHiEvMIQYg2+MxeR1ara65jmA+YIRMQHbATOBCqBN4FFqrq+S57rgVmq+hURuQy4SFUv3ddxD9YR7Gndw7u71rK1soKGPW201zhofZBQcx6RtiKCduYjeVugiZZgAy2helqDjbQGG4kHW3AiCaxcxZ+n5OaFKQgXUBAsoCSnhJJwCSU5JYzMG8mYvDFEg9GPVXEX/+dXqdvpTZYhwhGzj2XB1V8jr7ikI8/u1ir+uOF+lm75GzE7xrSi6UwsnMSEggmMyhtNXiCPSCBC2BfGJz5aq1JUr2mnfluMpl3JjIrojwi+sDselxUELAVLUVEURdVB1Q15aXqxcRfHTbEFtQW8RWwLcQ78zsjBIeVLkLISJC03Tfm8dV+cpJUg6Yt3LlachD/WkSZ87ST9cayQEghbBMN+coO55AXzyA/kkx/MpzBUSEGogOJwMWWRMspyyhiRO6LHU9pHoapsXPkqy375IwBGTZvJMed9hiNmH5eh+/KK5/n+yu/Smmwh4o8wb+R8JhZOYFhkOGU5w8gN5BL0BQlaQZyYUL22jeo17TRVJrATnfVU/BDIE3whkABIQN1uH6IobijLcRRVVyvHplOzDL3w9LKQdMqBX7B61Sq97mmUTK/7YiSsOEl/jIQvvbRjhRQrDIGQRW44l4g/QjQYJT/oaRYupDhUTFG4iOG5wxmZO5KScAm+j/kuhuM4rPvXCl584F7sVJITPrOIWQvPzwgFVTTt4O737uQfFc8T8oWYUjiV6SUzmFgwiYJwAdFglNxALgEJgG3RsDlO7aYYtRtjxBo6u5iKBf4cwZcDVhjEB/i8OoanmXbRzNPL8fQiXddsYdrCYs48o2eIcX8YLEfwCeA2VV3obd8MoKo/6pLnWS/PayLiB3YDZboPow7UEfzuxtuR6jGk/FHUl/kqvz9RRzBWQyhWTTC+l0C8jkCsjkCiCUttUAG6LDpw3l2B9qBFPGARSjoEkw7WviQSRayYO6uZ2OzPrZZthYiHR5MMlpAMlpAKFGL7cnB8OThWCBUfKj73H9xhFYh3FREUcVKI5w0sTSHqbouT9NZTiJN0P3NSWE4K0SRie6mTwnKSiCbdz5yku51e1OlyUfJS7bbdoUeX9Yy8A38XlpbGFqguCpHfmiJ3n/3MFbGSYCW8meg+Wi9FSAZLieWMIRksxvZHSfmjOL4wjhXEsYIgPhTp1MyLTQpOh06itqeL3amXp5WV1s3p1C2tizjdNfNSJ+HmsRPusbwz0pF21aZbqt3zHVq90uspS1CBgN13JRMrBVYckRRIqttRej9Hyl9AIjSMZGgYiVApti+C48vF9uWg4ncXq1uYVh0Ex03V6aaRmzZNsfna9285oLLvyxEMZMB4NNB1DrhKoLsr68ijqikRaQRKgIzuFyJyHXAdwNixYw/ImFBhDondVYTjOwiSIiwOYaeFoNOMRcr9D+Z4C3ne0in5oeg1nK4KESAX6P4QsT8PFbamSDlxHGwctTvuEF2UrpUtSBtB2oGdCIKodzcv4l2E06nl2efdMXqfi3Sup4++z8Lh95Zwpzn0WN3nvo/LQOnWVQoRCCCM6Xa2j9JLFRTF1gQpjaM4OOqgGc6hU7MACQLsQaj2bLDAEcTpqkPvqftkIohY4CPjzr/v38gCCQGhPjN1v7AeKIeqfrmp+5taknkB3J/6paqkNIGjKWxSOGqDd1/fvX5ZCCEaCdGYVsHVK12H7Mw6lqFbH3Xs6DkX98Mv0ZOBdAS9/azd9d6fPKjqPcA94D4RHIgxV918/YF8zWAwGIY8A9m8XQl07WIxBtjVVx4vNFQA1A2gTQaDwWDoxkA6gjeBySIyXkSCwGXA0m55lgJXeuufBV7YV/uAwWAwGPqfAQsNeTH/G4BncbuP/kFV14nI94BVqroUuA/4o4hsxn0SuGyg7DEYDAZD7wzo20Wq+hTwVLd9t3ZZjwGfG0gbDAaDwbBvDr9X4AwGg8HQrxhHYDAYDFmOcQQGg8GQ5RhHYDAYDFnOgA46NxCISA2w4wC/Xkq3t5azAFPm7MCUOTs4mDIfoaq9Trh82DmCg0FEVvU11sZQxZQ5OzBlzg4GqswmNGQwGAxZjnEEBoPBkOVkmyO4Z7ANGARMmbMDU+bsYEDKnFVtBAaDwWDoSbY9ERgMBoOhG8YRGAwGQ5aTNY5ARM4WkQ9EZLOI3DTY9gwEIlIuIitEZIOIrBORb3r7i0XkeRHZ5KVFg21rfyIiPhF5W0SWedvjRWSlV94/e8OgDxlEpFBEHhWR9z2tP5EFGv+H959eKyKPiEh4qOksIn8QkWoRWdtlX6+6isuvvevZeyJy7MGcOyscgYj4gLuAc4AZwCIRmTG4Vg0IKeDbqjodmAd8zSvnTcByVZ0MLPe2hxLfBDZ02f4J8EuvvPXANYNi1cBxB/CMqk4DZuOWfchqLCKjgW8Ac1T1SNxh7S9j6On8f8DZ3fb1pes5wGRvuQ64+2BOnBWOADgB2KyqW1U1AfwJuHCQbep3VLVKVd/y1ptxLxCjcct6v5ftfuDTg2Nh/yMiY4BPAb/3tgU4DXjUyzLUyhsFPok7lweqmlDVBoawxh5+IMebyTACVDHEdFbVf9Fzhsa+dL0QeEBdXgcKRWTkgZ47WxzBaODDLtuV3r4hi4iMA44BVgLDVbUKXGcBDBs8y/qdXwHfgY4Z30uABlVNedtDTesJQA2w2AuH/V5EchnCGqvqTuDnQAWuA2gEVjO0dU7Tl679ek3LFkcgvewbsv1mRSQPeAz4lqo2DbY9A4WInAdUq+rqrrt7yTqUtPYDxwJ3q+oxQCtDKAzUG15c/EJgPDAKyMUNjXRnKOn8UfTr/zxbHEElUN5lewywa5BsGVBEJIDrBB5S1ce93XvSj41eWj1Y9vUzJwIXiMh23HDfabhPCIVeCAGGntaVQKWqrvS2H8V1DENVY4AzgG2qWqOqSeBxYD5DW+c0fenar9e0bHEEbwKTvV4GQdyGpqWDbFO/48XH7wM2qOovuny0FLjSW78SeOJQ2zYQqOrNqjpGVcfhavqCqn4BWAF81ss2ZMoLoKq7gQ9FZKq363RgPUNUY48KYJ6IRLz/eLrMQ1bnLvSl61LgCq/30DygMR1COiBUNSsW4FxgI7AFuGWw7RmgMp6E+3j4HvCOt5yLGzdfDmzy0uLBtnUAyr4AWOatTwDeADYDS4DQYNvXz2U9Gljl6fw3oGioawx8F3gfWAv8EQgNNZ2BR3DbQJK4d/zX9KUrbmjoLu96tga3R9UBn9sMMWEwGAxZTraEhgwGg8HQB8YRGAwGQ5ZjHIHBYDBkOcYRGAwGQ5ZjHIHBYDBkOf6PzmIwZCciku66BzACsHGHdwBoU9X5g2KYwdDPmO6jBsN+ICK3AS2q+vPBtsVg6G9MaMhgOABEpMVLF4jIiyLyFxHZKCI/FpEviMgbIrJGRCZ6+cpE5DERedNbThzcEhgMnRhHYDAcPLNx50Q4CrgcmKKqJ+AOjf11L88duGPnHw9c7H1mMPxbYNoIDIaD5031xnkRkS3Ac97+NcCp3voZwAx3qBwAoiKSr+68EQbDoGIcgcFw8MS7rDtdth0665gFfEJV2w+lYQbD/mBCQwbDoeE54Ib0hogcPYi2GAwZGEdgMBwavgHM8SYaXw98ZbANMhjSmO6jBoPBkOWYJwKDwWDIcowjMBgMhizHOAKDwWDIcowjMBgMhizHOAKDwWDIcowjMBgMhizHOAKDwWDIcv4fB3VKtJS+DqwAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "run_simulation(100.0, model=m, y0={'CP': 0.75, 'pM': 0.25}).plot(legend=True)" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.1" } }, "nbformat": 4, "nbformat_minor": 1 }