{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Dual Phosphorylation Cycle" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "from ecell4.prelude import *" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/markdown": [ "Takahashi K, Tanase-Nicola S, ten Wolde PR, Spatio-temporal correlations can drastically change the response of a MAPK pathway. *Proceedings of the National Academy of Sciences of the United States of America*, **6**(107), 2473-8, 2010. [10.1073/pnas.0906885107](https://doi.org/10.1073/pnas.0906885107). PubMed PMID: [20133748](https://www.ncbi.nlm.nih.gov/pubmed/20133748)." ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "citation(20133748)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "@species_attributes\n", "def attrgen(radius, D):\n", " K | Kp | Kpp | KK | PP | K_KK | Kp_KK | Kpp_PP | Kp_PP | {\"radius\": radius, \"D\": D}\n", "\n", "@reaction_rules\n", "def rulegen(kon1, koff1, kcat1, kon2, koff2, kcat2):\n", " (K + KK == K_KK | (kon1, koff1)\n", " > Kp + KK | kcat1\n", " == Kp_KK | (kon2, koff2)\n", " > Kpp + KK | kcat2)\n", "\n", " (Kpp + PP == Kpp_PP | (kon1, koff1)\n", " > Kp + PP | kcat1\n", " == Kp_PP | (kon2, koff2)\n", " > K + PP | kcat2)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "radius, D = 0.0025, 1.0\n", "ka1, kd1, kcat1 = 0.04483455086786913, 1.35, 1.5\n", "ka2, kd2, kcat2 = 0.09299017957780264, 1.73, 15.0\n", "\n", "m = NetworkModel()\n", "m.add_species_attributes(attrgen(radius, D))\n", "m.add_reaction_rules(rulegen(ka1, kd2, kcat1, ka2, kd2, kcat2))" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "K | {'D': 1.0, 'radius': 0.0025}\n", "Kp | {'D': 1.0, 'radius': 0.0025}\n", "Kpp | {'D': 1.0, 'radius': 0.0025}\n", "KK | {'D': 1.0, 'radius': 0.0025}\n", "PP | {'D': 1.0, 'radius': 0.0025}\n", "K_KK | {'D': 1.0, 'radius': 0.0025}\n", "Kp_KK | {'D': 1.0, 'radius': 0.0025}\n", "Kpp_PP | {'D': 1.0, 'radius': 0.0025}\n", "Kp_PP | {'D': 1.0, 'radius': 0.0025}\n", "\n", "K + KK > K_KK | 0.04483455086786913\n", "K_KK > K + KK | 1.73\n", "K_KK > Kp + KK | 1.5\n", "Kp + KK > Kp_KK | 0.09299017957780265\n", "Kp_KK > Kp + KK | 1.73\n", "Kp_KK > Kpp + KK | 15.0\n", "Kpp + PP > Kpp_PP | 0.04483455086786913\n", "Kpp_PP > Kpp + PP | 1.73\n", "Kpp_PP > Kp + PP | 1.5\n", "Kp + PP > Kp_PP | 0.09299017957780265\n", "Kp_PP > Kp + PP | 1.73\n", "Kp_PP > K + PP | 15.0\n" ] } ], "source": [ "show(m)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "session = Session(model=m, y0={\"K\": 120, \"KK\": 30, \"PP\": 30})" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "ret = session.run(60.0)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEGCAYAAACKB4k+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dd5wU5f3A8c9393qv9HIHIoRebSCigAoYUaOJ0RhBEozBRGPUqPlFMTGJpGhMsUUULBHRKFgQJAQUYqFL7x2uwR3Xy5bn98fsLQfcHXtlb29vv++85jUzz8zOfCce+915npnnEWMMSimlFIAt0AEopZRqPTQpKKWU8tKkoJRSykuTglJKKS9NCkoppbzCAh1AU6SlpZmMjIxAh6GUUkFl3bp1x40x6bVtC+qkkJGRwdq1awMdhlJKBRUROVjXNq0+Ukop5aVJQSmllJcmBaWUUl5B3aaglPJdVVUVe/fupaysLNChqBYSExNDz549iYiI8PkzmhSUChF79+4lKSmJ3r17Y7NpJUFb53a7ycnJYc+ePfTt29fnz/ntL0NEXhaRXBHZUqPsjyKyQ0Q2ich7IpJUY9vDIrJHRHaKyFX+ikupUFVWVkb79u01IYQIm81G+/btKSsr4/PPP8fXzk/9+dcxB7j6jLKlQH9jzEBgF/AwgIj0BW4G+nk+86yI2P0Ym1IhSRNCaLHZbIgIX375JQcP1vkU6umf8VcwxpjPgPwzyj4xxjg9q18CXTzLk4F5xphKY8x+YA9wgb9i25ldzO8WbaesynnunZVSKsiJCMXFxT7tG8ifDXcAH3uWOwOHa2w74ik7i4hMF5G1IrI2Ly+vUSc+UlDGi5/tY8vRokZ9Ximl2qqAJAUR+SXgBN6oLqplt1orwIwxLxpjhhtjhqen1/qW9jkN7GI1ZXx9+GSjPq+Uapy4uDjv8qJFi+jVqxeHDh065+cOHDhA//79vev//Oc/GTp0KAUFBbXuP2XKFN555x0A8vPzGTJkCK+88kqDjxOKWjwpiMjtwDXAreZUy8cRoGuN3boAx/wVQ3p8JJ2Tovn6iCYFpQJh2bJl/OQnP2Hx4sV069bttG1jxozhwIEDdX72tdde429/+xuffPIJycnJ9Z6nsLCQq666iunTpzN16tRGHyeUtOgjqSJyNfAL4DJjTM2Hpd8H/iUiTwGdgF7Aan/GMqhroiYFFbIe/2Ar2441b/Vp304JPPbNfufcb+XKlfzwhz9k0aJF9OzZs0HnmD9/Pk8++STLli0jLS2t3n1LSkqYMGECt9xyC3fddVejjxNq/PlI6pvAF0BvETkiItOAvwPxwFIR2SgizwMYY7YC84FtwGJghjHG5a/YAAZ1SeJwfjknSir9eRqlVA2VlZVMnjyZBQsW0KdPnwZ99uDBg9x999188skndOjQ4Zz733fffYwaNYqf/exnTTpOqPHbnYIx5ru1FM+uZ//fAr/1Vzxnqm5X2HS0kMt7t2up0yrVKvjyi94fwsPDueSSS5g9ezbPPPOMt/yVV17xru/Zs4eJEycSERFBZmYm7733HgDp6emkpKQwf/78s77oa3PFFVewcOFC7r//ftq1O/VvvKHHCTUh+9DygC6JiGhjs1ItyWazMX/+fNasWcPvfvc7b/nUqVPZuHEjGzduZPjw4SxatIiNGzd6EwJYXTZ8/PHHPP/887zxxhu1Hf40N998M3fddRcTJ0487XHMhh4n1IRsNxdxkWH0ahenSUGpFhYTE8OHH37IpZdeSvv27Zk2bZrPn01PT2fx4sWMGTOGtLQ0rrqq/s4P7r33XrKysrj++utZtGhRo48TSkL2TgGsKqRNRwp9fv1bKdU8UlJSWLx4MU888QQLFy5s0GczMzN5//33ueOOO/jqq6/Ouf+sWbPo2rUrt912G263u9HHCRUSzF+Iw4cPN00Zee21Lw/yqwVbWPng5XRNiWnGyJRqfdatW8ewYcMCHYZqYevWreOzzz5j3LhxDBgwAAARWWeMGV7b/iF9pzCoSyKAPpqqlFIeIZ0U+nRIIMJuY9ORwkCHopRqhBkzZjB48ODTpldeeSXQYQW1kG1oBogIs9G3UwIbtbFZqaD0j3/8I9AhtDkhfacAVhXSlqOFuNzB27ailFLNRZNC1yTKqlzsyNYeU5VSKuSTwqjz0hCBpdtyAh2KUkoFXMgnhXYJUYzonsKizVmBDkUppQIu5JMCwMQBHdiVU8KeXN9GJlJKNU5jx1OoT0ZGBsePHwesZ/IzMzPZsGFDg45R3zgLdY3N4IsxY8ZQ/S7VgQMH6NWrF0uWLGlQbAB2u53BgwfTv39/brrpJsrKyuotbwpNCsCEAR0BWLQ5O8CRKBUa6htPoS4zZ85kzpw5dW7ftGkTN954I2+99RZDhgxpdGx1jbNQ39gMc+bMYebMmXUe88iRI1x11VX8+c9/blSXGtHR0WzcuJEtW7YQERHB888/X295U4T0I6nV2idEMbx7Mos2Z/HTsb0CHY5S/vfxQ5C9uXmP2WEATHjynLvVNp7ClClTiIqKYuvWreTk5PDUU09xzTXX+Hzq7du3c/vtt/Paa69xwQXW8O4zZ85k7969HD16lMOHD/Pggw/ywx/+sN7j1DXOQn1jM5xLdnY23//+93niiSe49tprASuJvPfee1RWVrJ//35uueUWHnvsMZ+Od+mll7Jp0yafyxtKk4LHxAEd+fWH29ibV0LP9Lhzf0Ap1WDV4ymsWLHirPEUDhw4wKeffsrevXu5/PLL2bNnD1FRUT4dd/Lkybz++uuMGjXqtPJNmzbx5ZdfUlpaypAhQ5g0aRKdOnWq9RjV4yxs2LDhrHEW7rvvPn7wgx80qqvt6oRw0003nVa+evVqtmzZQkxMDCNGjGDSpEkMH15rzxNeTqeTjz/+mKuvvtqn8sbQpOAxYUAHfv3hNj7enMXdV+jdgmrjfPhF7w91jacA8O1vfxubzUavXr3o0aMHO3bswG63c9tttwHWL+6IiAj+8pe/AFYVVGpqKgDjxo3jpZde4qqrrsJut3uPOXnyZKKjo4mOjubyyy9n9erVXHfddbXGVt84C7WNzXDixAnGjh0LWG0NVVVVLFiwALCqoKr7GRo3bhyvvfYaU6ZMISbmVB9r48eP98Z/ww03sGrVqjqTQnl5OYMHDwasO4LqnmXrKm8KTQoeHROjGdotiY82Z2tSUMpPqsdTGDduHL/73e945JFHvNtE5LR9RYQBAwawceNGwKoOysjIYMqUKWcd9+9//zs/+tGP+PGPf8wLL7xQ7zHrUj3OwqhRo2jXrh233nqrd9vNN9/MqFGjmDhxIsuXLyc+Pp7U1FRvbHPmzOHAgQO1tis8+OCDvP7669x0000sXLiQsLCwBsdW3Xbga3lTaENzDRMHdGR7VhH78koCHYpSbVb1eApvvPEGs2efGozx7bffxu12s3fvXvbt20fv3r19PqbNZuPNN99k586dPProo97yhQsXUlFRwYkTJ1ixYgUjRoyo9zjV4yw88sgjZz0ldO+99zJ27Fiuv/56qqqqfI4N4OmnnyYhIYFp06Z5u+pfunQp+fn5lJeXs2DBAkaOHNmgY/qLJoUarh3UiTCb8PqXTXtETilVv9rGU+jduzeXXXYZEyZM4Pnnn/e5PaFaZGQkCxcu5P333/f2iXTBBRcwadIkLrroIn71q1/V2Z5QU33jLNQ1NsO5iAhz584lKyuLBx98EIBRo0Zx2223MXjwYL71rW+dsz2hpYT0eAq1uXfeBv6zPZcvHr6C+KjwZj22UoHUmsdTmDJlCtdccw033nhjsx1z5syZxMXFcf/99zfbMZvLnDlzWLt2LX//+9/9fi4dT6GJpo3qQUmlk/lrjwQ6FKWUanHa0HyGAV0SGZGRzJzP9zPlkgzstrobf5RSzaO+l9Iaq7ZG382bN3ufZqoWGRnZ4sNxTpky5awG85pPM9VU8ymrlqBJoRbTRmXyo9fXs3RbDlf373DuDyilgkLNp5lam5pPMwWSVh/VYnzfDnRJjublVfsDHYpSSrUovyUFEXlZRHJFZEuNshQRWSoiuz3zZE+5iMhfRWSPiGwSkaH+issXdpsw5ZIMVh/I52sdlU0pFUL8eacwBzjzneuHgGXGmF7AMs86wASgl2eaDjznx7h88p0RXUmMDuevy3YHOhSllGoxfksKxpjPgPwziicDcz3Lc4HrapS/aixfAkki0tFfsfkiPiqc6aN7sGxHro7hrJQKGS3dptDeGJMF4Jm385R3Bg7X2O+Ip+wsIjJdRNaKyNq8vDy/Bnv7JRmkxEbw9NJdfj2PUqHCH+Mp1KUp4yD4evzMzEwGDx7M0KFD+eKLL+otDxatpaG5tuc+a32rzhjzojFmuDFmeHp6ul+DiosM487RPfh0Vx7rDp5506OUaqzGjKdQlxUrVtTaH1K1+sZBaKo//vGPbNy4kSeffJI777zznOXBoKUfSc0RkY7GmCxP9VCup/wI0LXGfl2AYy0cW61uu7g7/1y5j6eX7ub1H1wY6HCUahazVs9iR/6OZj1mn5Q+/OKCX5xzv4aMp9CUcQeg9nEQVqxYwaOPPkpqaio7d+5k9OjRPPvss9hsNuLi4rjzzjtZvnw5ycnJzJs3D19+fI4ePZo9e/b4XN6atfSdwvvA7Z7l24GFNcq/73kK6SKgsLqaKdBiIsL40WU9WbXnOF/sPRHocJQKatXjKSxYsKDO8RQ++ugjfvSjH1FRUQFY4w688cYbbNy4kbfffpuGdG1z3333MWrUqLO6wl69ejV//vOf2bx5M3v37uXdd98FoLS0lKFDh7J+/Xouu+wyHn/8cZ/O88EHH3i7kPClvDXz252CiLwJjAHSROQI8BjwJDBfRKYBh4DqUScWAROBPUAZ0Lz3eE30vYu68/Kq/Tzx0Tbev3uUvuWsgp4vv+j9oaHjKUDd4w5ceOGFVFZWUlJSQn5+vndcgVmzZnmHvKxtHASwOsrr0aMHAN/97ndZtWoVN954Izabje985zsAfO973+OGG26o93oeeOABnnjiCdLT00/r8bWu8mBwzqQgIjcBi40xxSLyf8BQ4AljzPr6PmeM+W4dm856j9tYvfLN8CHegIgKt/OLCX24Z95G3ll3mO+MaFodqFKhqqHjKdRXXt01xYoVK5gzZ06tXWXUNg5Cfcc8U31jHIDVdlBbJ351lQcDX6qPfuVJCKOAq7AeJQ34ewQt7dpBnRjWPZk/LtlFcYUj0OEoFbQaOp5CU8cdqG0chNWrV7N//37cbjdvvfWWdxhPt9vtfWLpX//611nDe4YCX5KCyzOfBDxnjFkIRPgvpNZJRHj0mr4cL6nkH8v3BjocpYJaQ8ZTaI5xB84cB+Hiiy/moYceon///mRmZnL99dcDEBsby9atWxk2bBj//e9/TxuwJ1ScczwFEfkQOAqMA4YB5cBqY8wg/4dXP3+Mp3AuP5//NR98fYyl942me2psi55bqaYIxvEU/DHuwIoVK/jTn/7Ehx9+eNa2uLg4Skra1siL/hhP4dvAEuBqY8xJIAV4oLkCDjYPXt2biDAb/7dgC8E8QJFSStXmnA3NxpgyEckFRgG7AadnHpLaJ0Txi6t786uFW3l3/VG+NaxLoENSKujVNZ5CbeMONNWYMWMYM2ZMrdtqu0uYMWMG//vf/04ru+eee5r9RbjWwpenjx4DhgO9gVeAcOB1oHWMMh0At17YnQUbj/Gbj7ZxWe900uIiAx2SUj5xu93YbK2lI4PgUD3eczBqyDjS1Xz567geuBYoBTDGHAPiG3ymNsRmE568YQBllS5+/cG2QIejlE9iYmLIzs5u1BeFCj5ut5vs7GwcjoY9LenLy2tVxhgjIgZARLR1FejVPp4Zl5/H0//ZxbWDOjGub/tAh6RUvXr27Mm2bds4duzYOZ+/V22Dw+HwdjgYHh7u02d8SQrzReQFrO6sfwjcAfyzsUG2JXeN6cnirdn84t+bWNx1NOnxWo2kWq+IiAjOO+883n77bQoKCnz+klDBzeFwkJycTJcuvrV/nvORVAARGQ9cidWb6RJjzNImRdlMAvFI6pl25RTzzb+t4uKeqbwyZYT+AlOtXklJCbt37/b2LaTatqioKHr16nVat+X1PZLqU99HniTQKhJBa3N++3h+OekbPLpwK69+cZDbL8kIdEhK1SsuLo4hQ4YEOgzVStXZ0CwixSJSVGNeVHO9JYNs7W67qDuX907nt4u2szO7ONDhKKVUo9WZFIwx8caYhBrzhJrrLRlkayci/OHGQSREhXPX6+u0bySlVNA65yOpItKttqklggsm6fGR/P2WIRzML+P+t7/Wt52VUkHJl/cUPqoxLQP2AR/7M6hgdVGPVB6e0IclW3N47lPtNE8pFXx86ebitGGDRGQoEFyDjragaaMy2Xj4JH9aspP+nRIZfb5/x5FWSqnm1OD33T2D64zwQyxtgogw61sD6dUunhlvrGdXjjY8K6WChy9tCvfVmO4XkX8BeS0QW9CKjQzj5akjiIqwM/WVNeQVVwY6JKWU8okvdwrxNaZIrLaFyf4Mqi3onBTN7NuHc6K0kh+8upbyKte5P6SUUgHm0xvNrVVreKP5XD7Zms2dr69jbJ92PPe9YYTbtYdKpVRgNWmQHRFZKiJJNdaTRWRJcwbYll3ZrwO/ntyf/2zP5YG3v8btDt4krJRq+3zp5iLdM+IaAMaYAhFp58eY2pzbLupOUbmDPy7ZSXxUOL+e3E/7SFJKtUq+JAWXiHQzxhwCEJHugP7cbaAfj+lJUbmDFz7bR0yEnYcm9NHEoJRqdXxJCr8EVonIp5710cB0/4XUNokID03oQ2mVkxc+24fTbfi/Sd/QxKCUalV8eXltseeFtYuwus7+mTHmeFNOKiI/A36AdcexGZgKdATmASnAeuA2Y0xVU87T2ogIv5ncnzCbjdmr9uNyGx77Zl9NDEqpVsOXhmYBrgaGGmM+AGJE5ILGnlBEOgM/BYYbY/oDduBmYBbwtDGmF1AATGvsOVozEeGxb/blB6MymfP5AR5+dzNOlw6PqJRqHXx5PvJZ4GLgu571YqCpI1mHAdEiEgbEAFnAFcA7nu1zgeuaeI5WS0T45aRvcPfl5zFvzWF+/MZ6Khz6HoNSKvB8SQoXGmNmABVgPX0ERDT2hMaYo8CfgENYyaAQWAecNMY4PbsdATrX9nkRmS4ia0VkbV5e8L5YLSLcf1VvHvtmX5Zuz+H7s1dTWK5dbiulAsuXpOAQETueJ45EJB1odH2HiCRjvRGdCXQCYoEJtexa6xNOxpgXjTHDjTHD09ODv7O5qSMzeebmIWw4XMCNz33OoRNlgQ5JKRXCfEkKfwXeA9qJyG+BVcDvmnDOccB+Y0yeMcYBvAtcAiR5qpMAugDHmnCOoHLtoE7MveMCcosrue7Z/7F6f36gQ1JKhahzJgVjzBvAg8Dvsap7rjPGvN2Ecx4CLhKRGE8j9lhgG7AcuNGzz+3AwiacI+hc0jONBTNGkhQdzq0vfclbaw4FOiSlVAiqs+8jEUmp74PGmEb/nBWRx4HvAE5gA9bjqZ059UjqBuB7xph6uxcNhr6PGqqwzMHdb65n5e7jfPeCrjz2zX5EhdsDHZZSqg2pr++j+pLCfqx6/doeojfGmB7NF2LjtMWkAOB0ufnz0l08t2IvAzon8tz3htIlOSbQYSml2ohGJYVg0FaTQrVPtmbz8/lfIwKzvjWQCQM6BjokpVQb0KReUj0HuFZE/uSZrmne8FRdruzXgQ9/OorMtFjuemM9D7+7ibIq57k/qJRSjeTLG81PAvdgNQZvA+4Rkd/7OzBl6Z4ay9s/uoS7xvRk3prDXPPXVaw/VBDosJRSbdQ5q49EZBMw2Bjj9qzbgQ3GmIEtEF+92nr10Zk+33OcB97ZRFZhOdNH9+Rn43sRGaaN0Eqphmly9RGQVGM5sekhqca45Lw0Ft97KTcN68rzn+5l0l9XseaAvtOglGo+viSF3wMbRGSOiMzF6pKiKS+vqSaIjwpn1o0DeWXqCMqrXNz0/Bc8/O5mCsu0iwylVNP59PSRiHQERmA9nvqVMSbb34H5ItSqj85UWunkL//ZxexV+0mJjeDBq/pw47Au2GzaFbdSqm6NfU9haH0HNcasb4bYmiTUk0K1LUcLeez9raw7WMCgLonMvLYfQ7olBzospVQr1dik4Aa2AtVdkdb8+WmMMVc0a5SNoEnhFGMMCzYe5feLdpBbXMk3B3XigSt70y1VX3pTSp2uvqRQ38hrPwe+BZRjdT/xnjGmxA/xqWYgIlw/pAvj+3bghU/38s+V+1i8JYvbLsrgx5f3JC0uMtAhKqWCgC+PpGZiDbAzGTgI/M4Ys7EFYjsnvVOoW3ZhBU8t3ck7644QFW5n6sgMpl/ak8SY8ECHppQKsCY9kmqM2Y/VY+knwAXA+c0bnvKHDolR/OHGQXzys8sY+432/GP5XkbN+i9/WrKT/NI2NfS1UqoZ1dem0ANr7OTJwGGsKqQPjTEVLRde/fROwXfbs4r423938/GWbKLC7Nx6YTfuGJVJp6ToQIemlGphTWlo3oR1l1DEGSOhGWOeauY4G0yTQsPtyS3m2eV7Wfj1MQT45qBO/ODSTPp10ncSlQoVjU0KM6ljSEwAY8zjzRJdE2hSaLwjBWW88r8DzFt9iNIqFxdkpDBlZAZX9m1PmN3XF92VUsFIu85WdSosd/DWmkO8+sVBjhSU0zExiu+M6Mp3RnSlY6JWLSnVFmlSUOfkchv+uyOX1748yMrdeQhwRZ923DS8K5f3bkdEmN49KNVWNPY9BRVC7DZhfN/2jO/bnsP5Zcxbc4j5a4/wn+25pMRGcN3gzlw/pDP9OydgDa2tlGqL6mtTuMcY84yIjDTG/K+F4/KJ3in4l9Pl5rPdebyz7ghLt+XgcBl6pMcyeVBnrhnUkZ7pcYEOUSnVCI1taN5ojBksIuuNMfX2gxQomhRazsmyKj7eks2CDUf5ar/VXXefDvFMGtCRCQM6cF67+ABHqJTyVWOTwpvAxUA6sLfmJqy+j3SQnRCVVVjOx5uz+WhzFusOWqPA9UiP5cq+HRjftx2DuyZj155alWq1Gt3QLCIdgCXAtWduM8YcbLYIG0mTQuBlFZbzn205LNmaw5f7TuB0G5Jjwrm8dzsu653Opb3SSYmNCHSYSqkamvz0kYhEcKp7i53GmFYxoosmhdalsNzByt15/Hd7Lst35lJQ5kAEBnZOZFSvNEael8bQbslEhesQokoFUpOSgohcBrwKHMCqOuoK3G6M+ayZ42wwTQqtl8tt2Hy0kM925fHprjw2Hj6Jy22ICrcxrHsyF2amclGPVAZ1TdRxppVqYU1NCuuAW4wxOz3r5wNvGmOGNSGgJOAloD/WW9N3ADuBt4AMrAT0bWNMQX3H0aQQPIorHKzen8+qPcf5cl8+O7KLMAYiwmwM6pLI8IwUhnVLZki3JFK1m2+l/KqpSWHTmY3KtZU1MKC5wEpjzEueqqkY4BEg3xjzpIg8BCQbY35R33E0KQSvk2VVrN6fz5oD+aw5UMCWo4U43dbfYkZqDIO7JjGoaxIDuyTRr1OCVjkp1YyamhRexvo1/5qn6FYgzBgztZHBJABfAz1MjZOLyE5gjDEmyzMm9ApjTO/6jqVJoe0or3Kx6chJNhw+yfqDBWw8fJLc4krAerGuV7s4+ndOpH+nBPp2SqRPx3gSonRsCKUao6lJIRKYAYzCalP4DHjWGFPZyGAGAy8C24BBwDrgHuCoMSapxn4FxpizBhoWkenAdIBu3boNO3gw4A9BKT/JLqzg6yMn2XykkC3HCtlytJDjJafGguiSHE2fDvH07hBP7w4J9GoXR4/0WG2jUOocWlXfRyIyHPgSGGmM+UpEnsHqmvsnviSFmvROIbQYY8gtrmRbVhHbjhWxI7uYndlF7M0rxeWperLbhO4pMfRIj6Nnu1h6psfRIy2WzLRYUmIjtIsOpWh9fR8dAY4YY77yrL8DPATkiEjHGtVHuQGITbViIkL7hCjaJ0Rxee923vJKp4v9x0vZlVPC7pxi9uSWsDevhE935eJwnfrRkxAVRkZaLN1TY+meEkO31Bi6pVhT+4QofeFOKQKQFIwx2SJyWER6e55oGotVlbQNuB140jNf2NKxqeAUGWanT4cE+nRIOK3c6XJzuKCcA8dL2Xe8lP3HSzh4ooyvD5/ko03HcNe4SQ6zCZ2SoumcFE3nZM88KZqOSVF0TIyiY2I0sZHaf6Rq++r9KxcRO/CkMeaBZj7vT4A3PE8e7QOmYo0XPV9EpgGHgJua+ZwqxITZbWR6qo4uP2Obw+Um62QFh/LLOJRfxtGTZRzOL+dIQRmrdh8np7iCM2tW4yPDaJ8YRYeEKNolRNIuPor2nnl6fCTp8ZGkxUUQFxmm1VQqaNWbFIwxLhEZJiJimrHxwRizEaitPmtsc51DqfqE221W9VFqTK3bq5xucooqOHaynKzCCrIKK8gpqiC7sILsogr27yslt7jitOqpapFhNtLirASRGhdJamwEKZ4pOTaClJgIkmPDSY6JICkmgsTocK26Uq2GL/fDG4CFIvI2UFpdaIx5129RKRVgEWE2uqbE0DWl9qQB4HYbTpY7yCuuJLe4grziSo6XVHK8pIrjxZUcL60ip6iCbceKyC+tosrlrvNYCVFh3gSRFBNOQnQ4CVHhJEaHEx8V5lkPIz4qjPgoqywuMoz4yHBiI+06hKpqNr4khRTgBHBFjTIDaFJQIc1mE+8dQO8O9XcdboyhtMpFfkkVBWWnppNlDs9URWG5g5PlDgrLHRw9WU5RuZOicke9yaRaVLiNuEgrUcRGhhEbEUZspJ2YyDBiwu3ERoYRHWEnJtxuzSPCiImwExVuJybCKosOtxMVbiMq3CqPDrcTGWbThBNizpkUGvuSmlLqFBHxfmnXVWVVlwqHi6IKB0XlTkoqnRRXOCiucFJS4aS40pqXVlnbSr2TixOlVRzML6O8ykVppZOyKpf3rfGGCLOJJ1HYiAyzE1k9D7MRGWYjIuzUekSNsgi7jXDPPCLMRrhdvGXhds92u1Ue7lkO8y4LYTbP3G4jzCaE1dxCTV8AABpwSURBVCiz26z97DYhzCbahtOMzpkUPH0dPQe0N8b0F5GBwLXGmCf8Hp1SyvvLvTnGMapyuimvclHmcFrzKhflDhcVDhcVDjdlVU4qHW4qnC7Kq6yyCqe1vcrp9q5XOd1UOt1UOFwUVzg57qyiyumiyuWmylljcrlrbXdpbnabeBNEzbm1bMNmgzCblUTsItg8+9hsgl2sz9vk1Geql625tV3E+qy1jLVNBJsN7zabWMvVn7OJ9Tlr/1Pr1cs197eWa3wGamwHvGXW+jc6JjCoa9I5/p9pOF+qj/4JPAC8AGCM2SQi/wI0KSgVZCI8v+ITabkuQowx3uRQ5XTjdFkJpcrlxukyOFxuHC43Tnf1svHs46LSWUWly0Gly4HD6aLKXUWV04nD7cDhduNwVVHlduFyO3G6nThcTpzGicvtwuEpdxs3TuPC5XbhNi6cxonbs+7AUGlcuI0bl3FhjBu3243bU2aMwY3bKvfMDZ7JuDHe/7kwxoCYGvsYwHjKa6x75oK1P1RPnLF+aructR0uTL+SV7r+rNn/e/mSFGKMMavPuD1zNnskSqmAMsZQ6aqkwllBhavCu1zpqjx9clpzh9vhLatyVVmTu8q77HA7Ts3dVThd1V/mp8odbof1ZV7L3G3O3ZbSJA1oKrGJDZvYsIvdu1yzTJCzyqsnwXN3gM3zq9+GzXNym4SB4P28tXRq2QrTjvX1K9gQax+xMa5bt2b+P8TiS1I4LiI98aQnEbkRyPJLNEopnzhcDkocJZQ4SihzlFHiKKHUUUqZo8yaO8soc5SdNi93lltzRznlznIqXBXW3FnhTQRNEWYLI8IWQaQ9knB7OBG2CMLt4YTbrOUIewThtnBiwmMIt4UTZgsj3BZ+2nJt8zBbGHaxe5fDbeHYbXbCJMw7D7NZyzaxWds9X97Vn6u5bBc7Npu1XnO/6rKaX/52sYdce4UvSWEGVgd2fUTkKLAfq6dUpVQjGWMocZRQWFlIYVUhhZWFFFUWWfOqIoqrir3z6qnEUeKdV7p864+y+ks4JiyG6LBo75QWnUZ0WDRRYVHesqiwKCLtkUSHRRNpjyTSHukti7JHERlmlUXYIrzL4bZwq8we4fl1q4KdL08f7QPGiUgsYDPGFPs/LKWCizGGwspCTlScIL8inxPlJzhRcYKCigLyK/I5WXnSmlecpKCygKLKIpym7lrYSHsk8RHx3ikxKpEu8V2Ii4gjPjye2PBY4iLirHl4HDHhMaeWw2KsROD5Ra5UQ/jy9FEq8BhW19lGRFYBvzbGnPB3cEoFmjGGoqoiskuzyS3LtabyXPLK8qypPI/j5cc5UXECp/vsL3lBSIpMIiUqhaSoJDITMxkSNYSkyCSSIpNIjEwkMSLRmkcmkhCRQEJkApF2HX1OBYYv1UfzsMZQ+JZn/VasYTPH+SsopVqKw+0gpzSHYyXHOFZ6jKySLLJKszhWeoyc0hxyynIod5af9bmUqBTSotNIj07nvKTzSItOIy06jdToVFKjUkmJSiElOoXEiETsNh3fQQUPn95oNsb8psb6EyJynb8CUqq5lTvLOVR0iEPFhzhcfJhDRYc4UnKEI8VHyC7NxmVc3n0FIT06nQ6xHeiV3IvRXUbTPqY97WPb0z6mPe1i2pEenU64XatlVNvkS1JYLiI3A/M96zcCH/kvJKUazhhDblku+wr3sa9wH/sL93Og8AD7i/aTW3b60BwpUSl0ie/CoPRBTMycSJf4LnSK60Tn2M50iO2gX/gqpNWZFESkGOsxVAHuA173bLIBJVjtDEq1uMLKQnYV7GJXwS72nNzD7oLd7Du5j2LHqWcg4sPjyUzM5MIOF9I9oTvdE7vTLb4b3eK7ERcRF8DolWrd6kwKxphmeKleqcYzxpBTlsO2E9vYkb+D7fnb2ZG/g+zSbO8+iZGJ9ErqxcQeE+mZ1JOeiT3pkdSD1KjUkHu+XKnm4NNQUp7+jjJq7q9dZ6vmVlRVxJa8LXx9/Gu2Ht/KluNbOFFhPeRmExsZCRkMaTeEPil96J3cm/OTzyctOk2//JVqRr48kvoyMBDYClS/d65dZ6smMcZwsOggG3I3sDFvIxtzN7KvcB9gNfZmJmYysvNI+qX2o29qX85PPp+Y8Ib1LqqUajhf7hQuMsb09Xskqk0zxrC/cD9fZX/Fupx1rMtZx/Hy4wAkRCQwuN1gJvWYxMD0gfRP7a/1/koFiC9J4QsR6WuM2eb3aFSbklOawxdZX/DFsS9Ynb3amwQ6xHbgwo4XMrTdUIa1H0ZmYqZ2kaBUK+FLUpiLlRiygUqsp5GMMWagXyNTQcfhcrAhdwMrj65k5ZGV7C3cC0BqVCoXdryQCzteyIgOI+gS10XbAZRqpXxJCi8DtwGbOdWmoBRgPR668uhKVhxewaqjqyh1lBJmC2NY+2Fcd951XNzpYs5PPl+TgFJBwpekcMgY877fI1FB43j5cf576L98cvAT1mavxWVcpEWncXXG1YzuMpqLOl6kjcJKBSlfksIOz0hrH2BVHwH6SGqoKagoYOnBpSw5sIQ12WswGDISMpjSbwpju42lX1o/bRdQqg3wJSlEYyWDK2uU6SOpIaDCWcGKwyv4YN8HfH70c5zGSUZCBncOupMru1/JeUnnabWQUm2ML+MpTPXHiUXEDqwFjhpjrhGRTKweWVOA9cBtxpgqf5xb1c0Yw5bjW3hvz3t8vP9jShwltItpx239bmNS5iRtH1CqjfPl5bVX8I4qfYox5o4mnvseYDuQ4FmfBTxtjJknIs8D04DnmngO5aPiqmI+2PsB7+x+h90Fu4myRzG++3iuPe9aRrQfod0/KxUifKk++rDGchRwPXCsKScVkS7AJOC3wH1i/fS8ArjFs8tcYCaaFPxuZ/5O3tzxJov2L6LcWU6/1H786qJfMSFzAvER2v2VUqHGl+qjf9dcF5E3gf808bx/AR4Eqr91UoGTxnjHJzwCdK7tgyIyHZgO0K1btyaGEZpcbhefHvmU17e/zprsNUTZo5jYYyLf7v1t+qX2C3R4SqkA8qlDvDP0Ahr9bSwi1wC5xph1IjKmuriWXc+qsgIwxrwIvAgwfPjwWvdRtat0VfL+3veZu3UuB4sO0jG2I/cNu48bet1AYmRioMNTSrUCvrQpVI+rUC0b+EUTzjkSuFZEJmJVRyVg3TkkiUiY526hC02solKnlDpKmbdjHq9ue5X8inz6pfbjT5f9ibHdxhJma8zvAqVUW+VL9VGzViwbYx4GHgbw3Cncb4y5VUTexhrVbR5wO7CwOc8bikqqSvjXjn/x6rZXKaws5JJOlzCt/zRGdBihTxAppWpV38hr9VYRGWMONXMsvwDmicgTwAZgdjMfP2RUOCuYt2Mes7fM5mTlSS7rchl3DryTAekDAh2aUqqVq+9O4SNODcdZzQDpQDugyc8oGmNWACs8y/uAC5p6zFDmdDtZsGcBz218jtzyXEZ2GslPhvyEfmnaeKyU8k19w3Ge9rNSRDKwfs2PA37n16hUgxhjWHV0FU+te4o9J/cwOH0ws0bPYniH4YEOTSkVZHxpaO4F/BK4EPgz8FNjjMPfgSnfHCg8wJNrnuR/R/9H1/iuPD3macZ2G6ttBkqpRqmvTaE/VjLoB/wBmGaMcbVUYKp+ZY4yXtz0InO3zSXKHsX9w+/nlj63EG4PD3RoSqkgVt+dwtfAYay2hQuAC2r++jTG/NS/oam6fHbkM5748gmySrOY3HMy9w67l7TotECHpZRqA+pLCk3t20g1s+Plx3ly9ZMsObCEnok9eXXCqwxpNyTQYSml2pD6GprntmQgqn6L9y/mia+eoMxRxt2D7+aO/ndoVZFSqtnp66ytXEFFAb/96rcsObCEAWkDeGLkE/RI6hHosJRSbZQmhVbsy6wveWTlIxRUFvDTIT9lav+p2i2FUsqv9BumFXK4HTy78Vlmb55NRmIGz457lj4pfQIdllIqBJxzUF0ROV9ElonIFs/6QBH5P/+HFpqyS7OZungqL21+iRt63cC8SfM0ISilWowvI63/E6sDOweAMWYTcLM/gwpVq7NW850Pv8Pugt38cfQfmXnJTGLCYwIdllIqhPhSfRRjjFl9xhuyzrp2Vg1njGHu1rk8vf5puid055WrXtHGZKVUQPiSFI6LSE88YyqIyI1All+jCiFVrioe/+Jx3t/7PuO7j+c3I39DbHhsoMNSSoUoX5LCDKyRzvqIyFFgP/A9v0YVIvIr8vnZ8p+xPnc9dw26i7sG3aV9FimlAsqXQXb2AeNEJBawGWOK/R9W23eo6BB3Lr2TvPI8/jD6D0zInBDokJRSyqdeUiOBbwEZQFj1L1ljzK/9GlkbtvX4Vn687McYY3j5qpcZmD4w0CEppRTgW/XRQqAQWAdU+jectu/zY59z7/J7SY5M5oXxL5CRmBHokJRSysuXpNDFGHO13yMJAcsPLefnn/6czMRMnhv3HO1i2gU6JKWUOo0v7yl8LiI6uG8TLT24lPtW3Efv5N68fNXLmhCUUq1SfYPsbAHcnn2misg+rOojAYwxRivCfbR4/2IeWvkQ/dP689y454iPiA90SEopVav6qo86A4NbKpC2atnBZTy08iEGpQ/i2XHP6jsISqlWrb6ksN8Yc7DFImmDPj/2OQ989gD90vrx3LjntMsKpVSrV19SaCci99W10RjzlB/iaTM25G7g3uX3kpmYybNjn9WEoJQKCvUlBTsQh9WGoBpg78m9zPjPDNrHtOeF8S+QGJkY6JCUUson9SWFLH+8oCYiXYFXgQ5YDdkvGmOeEZEU4C2sl+QOAN82xhQ09/n97Xj5cX78nx8TGRbJC+NfIC06LdAhKaWUz+p7JNVfdwhO4OfGmG8AFwEzRKQv8BCwzBjTC1jmWQ8qZY4y7l52NwWVBfx97N/pFNcp0CEppVSD1JcUxvrjhMaYLGPMes9yMbAd60mnycBcz25zgev8cX5/cRs3D618iO352/nD6D/QL7VfoENSSqkGqzMpGGPy/X1yEckAhgBfAe2NMVmec2cBtb7dJSLTRWStiKzNy8vzd4g+e+7r51h+eDkPjniQMV3HBDocpZRqFF/eaPYLEYkD/g3ca4wp8vVzxpgXjTHDjTHD09PT/RdgAyw/tJznv36eyT0nc0ufWwIdjlJKNVpAkoKIhGMlhDeMMe96inNEpKNne0cgNxCxNdSBwgM8suoR+qb25VcX/0rHQ1BKBbUWTwpifWvOBraf8a7D+8DtnuXbsXpnbdXKHGXcu/xewm3h/GXMX4i0RwY6JKWUahJfekltbiOB24DNIrLRU/YI8CQwX0SmAYeAmwIQW4P8Yc0f2Fe4jxfGv0DHuI6BDkcppZqsxZOCMWYVdT/u6pcnnvxh6cGl/Hv3v5nWfxoXd7o40OEopVSzCFhDczDLLs1m5ucz6Z/anxlDZgQ6HKWUajaaFBrI5Xbx8MqHcbgdzBo9i3BbeKBDUkqpZqNJoYHm7ZzH2py1PHzBw3RL6BbocJRSqllpUmiAoyVHeWb9M4zqPIrrzguqF66VUsonmhR8ZIzh8c8fRxAevehRfR9BKdUmBeKR1KC0cO9Cvsj6gl9e+MvW8fipywGOMnBWgqPcmjsrTs1dleCsAteZk8MzVYHbaS27nbVMLmsyrtPLjPvU3LhqLPsyGc/kWcecKj9r2Zy+7J1jLUP928/c17tY87OcY/+a+9TUXPvUsX+du/hyHF8157H8qFmvuY0ZdS+Mm9nsh9Wk4IMT5Sf445o/MrTdUL7d+9tNO5gxUFEIZSegLB/KC6DiJJSftMorC6GiCCqLoaoEKkugqhiqyqCq1EoAjjJwO5rn4gBsYWAL98ztnnn1smddPMtiB5vNmovtVJnYrMkebs0Rz7YaywiIePYVz7ovy5z6bPXTzN47tRrlZ20/c9+a+3Nq+5n7nLl/bZprnzr3r3MnH/bxUdDc7QZLnC2s20V+OawmBR/8bcPfKHOU8dglj2GTemrcHOVQeAROHoLCw1CUBcXHoDgbSnKgJA9K8+r/Qg+LhqgEiIyHiDhrntAFImIhIgbCa05REBYF4dHW3DtFgD0SwiLBHuGZh1tf/PYIa7l63WYPoi8HpZS/aVI4h60ntvLu7nf5ft/v0yOxh1VYkgu52yB3B+TtgBN7IH8fFB0949MCce0gvgPEtYf2AyAuHWLTISYVolMgOtkzJUFkgvWFrpRSAaJJoR7GGGZ99XuSw2O5s7QS3rwFjm2wfv1Xi06G1F6QcSmk9IDk7pDUDRK7WsnAru8xKKWChyaF2hTnwLYFfLx9HhvkODPzThC/6w+eL/+R0HEwdOgP6d+w7gS0+kUp1UZoUqjp5CF4/6ew/1PKMTzVrRvfCE/kukm/t5JBdHKgI1RKKb/SpFCtKAvmfhPKCuDS+/lXtJ2cXW8wa9wz2NsPC3R0SinVIjQpgPVU0KvXQukJ+P5Citudz8v/vppRnUcxTBOCUiqE6BvNlcXw2vVw8jDcOh+6DOO1ba9RVFXET4b8JNDRKaVUi9I7hS3vQs5muGU+dL+EkxUneXXbq4zvPp6+qX0DHZ1SSrUovVPY/5n1DkGvKwF4ecvLlDnKmDFYx0lQSoWe0E4KxlhJIXM0iHC8/Dhv7niTST0m0TOpZ6CjU0qpFhfaSSFvB5TmQuZlAMzePBuH28Fdg+4KcGBKKRUYoZ0U9n1qzXtcRk5pDvN3zufantfq4DlKqZAV2klh/6eQnAFJ3Zi9ZTZu42b6wOmBjkoppQImdJOCywkHVkHmZWSXZvPOrneYfN5kusR3CXRkSikVMKGbFLK+hsoi6HEZL21+CWMMPxz4w0BHpZRSARW6SWG/1Z6wLyWDf+/+N9f1uo7OcZ0DHJRSSgVWq0sKInK1iOwUkT0i8pDfTrT/U3Z26MPUlfcRHx7PnQPv9NuplFIqWLSqpCAiduAfwASgL/BdEWn+14odFWzIXsfUmCoi7BHMmTCHDrEdmv00SikVbFpVUgAuAPYYY/YZY6qAecDk5j7JF5tfY3p6IqmRSbx69aunRlRTSqkQ19qSQmfgcI31I54yLxGZLiJrRWRtXl5eo07SISadYbY45lw1h45xHRsfrVJKtTGtrUO82oYwM6etGPMi8CLA8OHDTS37n1Nmn+t4vs91jfmoUkq1aa3tTuEI0LXGehfgWB37KqWUamatLSmsAXqJSKaIRAA3A+8HOCallAoZrar6yBjjFJG7gSWAHXjZGLM1wGEppVTIaFVJAcAYswhYFOg4lFIqFLW26iOllFIBpElBKaWUlyYFpZRSXpoUlFJKeYkxjXr/q1UQkTzgYCM/ngYcb8ZwAkmvpXVqK9fSVq4D9FqqdTfGpNe2IaiTQlOIyFpjzPBAx9Ec9Fpap7ZyLW3lOkCvxRdafaSUUspLk4JSSimvUE4KLwY6gGak19I6tZVraSvXAXot5xSybQpKKaXOFsp3Ckoppc6gSUEppZRXSCYFEblaRHaKyB4ReSjQ8TSEiLwsIrkisqVGWYqILBWR3Z55ciBj9IWIdBWR5SKyXUS2isg9nvJgvJYoEVktIl97ruVxT3mmiHzluZa3PN3BBwURsYvIBhH50LMelNciIgdEZLOIbBSRtZ6yYPwbSxKRd0Rkh+ffzMX+uo6QSwoiYgf+AUwA+gLfFZG+gY2qQeYAV59R9hCwzBjTC1jmWW/tnMDPjTHfAC4CZnj+OwTjtVQCVxhjBgGDgatF5CJgFvC051oKgGkBjLGh7gG211gP5mu53BgzuMYz/cH4N/YMsNgY0wcYhPXfxj/XYYwJqQm4GFhSY/1h4OFAx9XAa8gAttRY3wl09Cx3BHYGOsZGXNNCYHywXwsQA6wHLsR62zTMU37a311rnrBGPFwGXAF8iDVMbrBeywEg7YyyoPobAxKA/XgeDPL3dYTcnQLQGThcY/2IpyyYtTfGZAF45u0CHE+DiEgGMAT4iiC9Fk91y0YgF1gK7AVOGmOcnl2C6e/sL8CDgNuznkrwXosBPhGRdSIy3VMWbH9jPYA84BVPld5LIhKLn64jFJOC1FKmz+UGiIjEAf8G7jXGFAU6nsYyxriMMYOxfmVfAHyjtt1aNqqGE5FrgFxjzLqaxbXs2uqvxWOkMWYoVnXxDBEZHeiAGiEMGAo8Z4wZApTixyqvUEwKR4CuNda7AMcCFEtzyRGRjgCeeW6A4/GJiIRjJYQ3jDHveoqD8lqqGWNOAiuw2kmSRKR6dMNg+TsbCVwrIgeAeVhVSH8hOK8FY8wxzzwXeA8rYQfb39gR4Igx5ivP+jtYScIv1xGKSWEN0MvzNEUEcDPwfoBjaqr3gds9y7dj1c+3aiIiwGxguzHmqRqbgvFa0kUkybMcDYzDaghcDtzo2S0orsUY87AxposxJgPr38Z/jTG3EoTXIiKxIhJfvQxcCWwhyP7GjDHZwGER6e0pGgtsw1/XEehGlAA13EwEdmHV+/4y0PE0MPY3gSzAgfULYhpWne8yYLdnnhLoOH24jlFYVRCbgI2eaWKQXstAYIPnWrYAj3rKewCrgT3A20BkoGNt4HWNAT4M1mvxxPy1Z9pa/W89SP/GBgNrPX9jC4Bkf12HdnOhlFLKKxSrj5RSStVBk4JSSikvTQpKKaW8NCkopZTy0qSglFLKK+zcuyilRKT68T+ADoALq+sBgDJjzCUBCUypZqaPpCrVQCIyEygxxvwp0LEo1dy0+kipJhKREs98jIh8KiLzRWSXiDwpIrd6xlrYLCI9Pfuli8i/RWSNZxoZ2CtQ6hRNCko1r0FYYxEMAG4DzjfGXAC8BPzEs88zWGMTjAC+5dmmVKugbQpKNa81xtOdsYjsBT7xlG8GLvcsjwP6Wt0/AZAgIvHGmOIWjVSpWmhSUKp5VdZYdtdYd3Pq35sNuNgYU96SgSnlC60+UqrlfQLcXb0iIoMDGItSp9GkoFTL+ykwXEQ2icg24EeBDkipavpIqlJKKS+9U1BKKeWlSUEppZSXJgWllFJemhSUUkp5aVJQSinlpUlBKaWUlyYFpZRSXv8P2jOuVqthj5UAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret.plot(y=[\"K+K_KK\", \"Kp+Kp_KK+Kp_PP\", \"Kpp+Kpp_PP\"], legend=True)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "ret = session.run(60.0, ndiv=100, solver='gillespie')" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEGCAYAAACKB4k+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3de3wV1bnw8d8TLmJECiggKhSwGIuigBGkIqLSitGK9kCrWCqpFbS21Xq899TLkXrk7VHraW3VasHjFbBqvCDqy2tqtRUaNOWmURBEkJtiVUC5mOf9Y2ZvJmHPzuzL7Jmd/Xw/n3ySPXv2zDPJZNas9axZS1QVY4wxBqAs6gCMMcbEhxUKxhhjkqxQMMYYk2SFgjHGmCQrFIwxxiS1jTqAXOy///7ap0+fqMMwxpiisnDhwg9VtVuq94q6UOjTpw91dXVRh2GMMUVFRN7ze8+aj4wxxiRZoWCMMSbJCgVjjDFJRZ1TMMYEt2PHDlasWMG2bduiDsUUSHl5OYcccgjt27cP/BkrFIwpEStWrKBz585UVFRQVmaNBK1dY2MjGzZsYPny5QwYMCDw50I7M0TkTyKyUUSWeJb9WkTeEpFFIvKEiHT2vHeNiCwXkQYROSWsuIwpVdu2baNHjx5WIJSIsrIyevTowbZt2/jb3/5G0MFPwzw7ZgBjmi17EThCVY8E3gauARCRAcDZwOHuZ34vIm1CjM2YkmQFQmkpKytDRHjttdd47z3fXqhNhNZ8pKovi0ifZste8Lx8DRjn/jwWeFRVtwMrRWQ5MBT4exix3fj0UpZ98Gny9dhBBzFhWO8wdmWMMZETET777LNA60Z52/BD4Dn354OA9z3vrXGX7UFEJotInYjUbdq0Kecglq37lJr6tTlvxxhjWoNICgUR+QWwC3gosSjFaikbwFT1HlWtVNXKbt1SPqXdouu/fTgzpwxn5pThDOjZKattGGMy17Fjx+TPc+bMoX///qxevbrFz61atYojjjgi+fqPf/wjQ4YM4eOPP065/qRJk3jssccA2Lx5M4MHD2b69OkZb6cUFbxQEJHzgNOBc3V35mMN0Muz2sHAB4WOzRhTGPPmzeOnP/0pc+fOpXfvpk23o0aNYtWqVb6ffeCBB/jtb3/LCy+8QJcuXdLu55NPPuGUU05h8uTJVFdXZ72dUlLQLqkiMga4CjhBVb2dpZ8CHhaR24ADgf7AgkLGZkwpaZ5Xy4cBB3bi+m8f3uJ6f/3rX7nggguYM2cOhxxySEb7mDVrFrfccgvz5s1j//33T7vuli1bOPXUU5kwYQIXXXRR1tspNaEVCiLyCDAK2F9E1gDX4/Q22gt4UUQAXlPVC1V1qYjMApbhNCtdrKpfhhWbMSYa27dvZ+zYsdTW1nLYYYdl9Nn33nuPn/zkJ7zxxhsccMABLa5/2WWX8aMf/Yif//znOW2n1ITZ++icFIvvS7P+r4BfhRWPMWa3IHf0YWjXrh3f+MY3uO+++7jjjjuSy6dPn558vXz5cqqqqmjfvj19+/bliSeeAKBbt2507dqVWbNm7XGhT+Wkk06ipqaGyy+/nO7duyeXZ7qdUmOdlo0xBVNWVsasWbP4xz/+wc0335xcXl1dTX19PfX19VRWVjJnzhzq6+uTBQI4QzY899xz3HXXXTz00EOpNt/E2WefzUUXXURVVVWT7piZbqfU2DAXxpiCKi8v55lnnuH444+nR48enH/++YE/261bN+bOncuoUaPYf//9OeWU9IMfXHrppaxbt46zzjqLOXPmZL2dUmI1BWNMwXXt2pW5c+cydepUampqMvps3759eeqpp/jhD3/I/PnzW1x/2rRp9OrVi4kTJ9LY2Jj1dkqF1RSMMQWzZcuW5M+9evVi5cqVe6xTW1u7x7I+ffqwZElyGDWOOuoo1q71f+h0xowZTV5Pnz49+XMm2ylFVlMwxhiTZDUFY0zRuvjii3n11VebLLvkkkv2eFDNBGeFgjGmaN15551Rh9DqWPORMcaYJKsp4IyU+r27d4/SbUNpG2NKVckXCmMHNR2he9k6ZzwYKxSMMaWo5AuFCcN6NykAvDUGY4wpNZZTSCHRnPS9u//Ow/NbHuvdGBNMtvMppNOnTx8+/PBDABYuXEjfvn154403MtpGunkW/OZmCGLUqFHU1dUl99G/f3+ef/75jGIDaNOmDYMGDeKII45g/PjxbNu2Le3yXFih0MzYQQclJ96xWdmMCUe6+RT83HDDDXs8lOa1aNEixo0bx8yZMxk8eHDWsfnNs5BuboYZM2Zwww03+G5zzZo1nHLKKdx6661ZDamx9957U19fz5IlS2jfvj133XVX2uW5KPnmo+a8zUnWlGRareeuhvWL87vNAwbCqbe0uFqq+RQmTZpEhw4dWLp0KRs2bOC2227j9NNPD7zrN998k/POO48HHniAoUOHAk4hsmLFCtauXcv777/PlVdeyQUXXJB2O37zLKSbm6El69ev5wc/+AFTp07ljDPOAJxC5IknnmD79u2sXLmSCRMmcP311wfa3vHHH8+iRYsCL8+UFQrGmIJJN5/CqlWr+Mtf/sKKFSs48cQTWb58OR06dAi03bFjx/Lggw8yYsSIJssXLVrEa6+9xtatWxk8eDCnnXYaBx54YMptpJtnwW9uhiASBcL48eObLF+wYAFLliyhvLycY445htNOO43Kysq029q1axfPPfccY8aMCbQ8G1YoGFOKAtzRh8FvPgWA7373u5SVldG/f3/69evHW2+9RZs2bZg4cSLg3HG3b9+e3/zmN4DTBLXffvsBMHr0aO69915OOeUU2rRpk9zm2LFj2Xvvvdl777058cQTWbBgAWeeeWbK2NLNs5BqboaPPvqIk08+GXByDTt27ODJJ58EnCaogQMHJmN74IEHmDRpEuXl5cltfvOb30zG/53vfIdXXnnFt1D4/PPPGTRoEODUCBIjy/otz4UVCsaYgknMpzB69Ghuvvlmrr322uR77myMTV4PHDiQ+vp6wGkO6tOnD5MmTdpju7/73e+48MIL+fGPf8zdd9+ddpt+EvMsjBgxgu7du3Puuecm3zv77LMZMWIEVVVVvPTSS+y7777st99+ydhmzJjBqlWrUuYVrrzySh588EHGjx9PTU0Nbdu2zTi2RO4g6PJcWKLZGFNQifkUHnroIe67b/dkjLNnz6axsZEVK1bw7rvvUlFREXibZWVlPPLIIzQ0NHDdddcll9fU1PDFF1/w0UcfUVtbyzHHHJN2O4l5Fq699to9egldeumlnHzyyZx11lns2LEjcGwAt99+O506deL8889HVQF48cUX2bx5M59//jlPPvkkxx13XEbbDIvVFIwxBZeYT2HkyJHJhG5FRQUnnHACGzZs4K677gqcT0jYa6+9qKmp4YQTTqBHjx4ADB06lNNOO43Vq1fzy1/+0jef4JWYZ6GqqorHH3+8yXvTpk2jurqaiRMn8sgjj1BWFuy+WkS4//77Of3007nyyis5/PDDGTFiBBMnTmT58uVMmDChxXxCoUii1CpGlZWVmugDHIZE76OZU4aHtg9jCmXhwoUcffTRUYeR0qRJkzj99NMZN25c3rZ5ww030LFjRy6//PK8bTNfZsyYQV1dHb/73e9C39fChQt5+eWXGT16dDLPISILVTVlKWTNR8YYY5Ks+cgYE7l0D6VlK1XSd/HixcneTAl77bVXwafjnDRp0h4Jc29vJi9vL6tCsELBGFMyvL2Z4sbbmylK1nxkjDEmKbRCQUT+JCIbRWSJZ1lXEXlRRN5xv3dxl4uI/I+ILBeRRSIyJKy4jDHG+AuzpjADaP7M9dXAPFXtD8xzXwOcCvR3vyYDfwgxLmOMMT5CKxRU9WVgc7PFY4H73Z/vB870LP9fdbwGdBaRnmHFZowxJrVC5xR6qOo6APd7d3f5QcD7nvXWuMv2ICKTRaROROo2bdoUarDGmPwKYz4FP7nMgxB0+3379mXQoEEMGTKEv//972mXF4u4JJpTDfqR8qk6Vb1HVStVtbJbt24hh2UT7hgThmzmU/BTW1ubcjykhHTzIOTq17/+NfX19dxyyy1MmTKlxeXFoNBdUjeISE9VXec2D210l68BennWOxj4oMCx7cE7f7PN3Wxak2kLpvHW5rfyus3Duh7GVUOvanG9TOZTyGXeAUg9D0JtbS3XXXcd++23Hw0NDYwcOZLf//73lJWV0bFjR6ZMmcJLL71Ely5dePTRRwly8zly5EiWL18eeHmcFbqm8BRwnvvzeUCNZ/kP3F5IxwKfJJqZojRhWG9mThnOzCnDk7OxGWOyl5hP4cknn/SdT+HZZ5/lwgsv5IsvvgCceQceeugh6uvrmT17NpkMbXPZZZcxYsSIPYbCXrBgAbfeeiuLFy9mxYoVyTGOtm7dypAhQ3j99dc54YQTuPHGGwPt5+mnn04OIRFkeZyFVlMQkUeAUcD+IrIGuB64BZglIucDq4HErBNzgCpgObANyG8dzxjTRJA7+jBkOp8C+M87MGzYMLZv386WLVvYvHlzcl6BadOmJae8TDUPAjgD5fXr1w+Ac845h1deeYVx48ZRVlbG9773PQC+//3v853vfCft8VxxxRVMnTqVbt26NRnx1W95MWixUBCR8cBcVf1MRP4DGAJMVdXX031OVc/xeWuP57jVGZXv4gDxGmOKWKbzKaRbnhiaora2lhkzZqQcKiPVPAjpttlcujkOwMkdpBrEz295MQjSfPRLt0AYAZyC05XUniMwxmQl0/kUcp13INU8CAsWLGDlypU0NjYyc+bM5DSejY2NyR5LDz/88B7Te5aCIM1HX7rfTwP+oKo1InJDeCEZY1q7TOZTyMe8A955EKZMmcLw4cO5+uqrWbx4MSNHjuSss84CYJ999mHp0qUcffTRfOUrX2HmzJn5O+giEaRQWCsidwOjgWkishfx6cpqjCkiW7ZsSf7cq1cvVq5cCcATTzzBcccdx+23377HZ7p375523oFRo0YxatSoPZY3b05KPKNQW1tLeXm57wX/pptu4qabbmrpUHxHdg1jxNdCCnJx/y7wPDBGVf8FdAWuCDUqY4wxkWixpqCq20RkIzACeAfY5X43xpi88Lu7TjXvQK78ahbQtCaTcPHFF/Pqq682WXbJJZfk/UG4uAjS++h6oBKoAKYD7YAHgXjMMm2MCayxsTHwvMLGceedd0YdQtYaGxsz/kyQs+Ms4AxgK4CqfgDsm/GejDGRKi8vZ/369VldKEzxaWxsZP369ezcuTOjzwVJNO9QVRURBRCRfbIJ0BgTrUMOOYRly5bxwQcftNj/3rQOO3fuTA442K5du0CfCVIozHJ7H3UWkQuAHwJ/zDbI1uDh+aupqV+bfD120EF7jInUfB2vVOsbE7b27dvzta99jdmzZ/Pxxx8HvkiY4rZz5066dOnCwQcfHGj9IInm/xaRbwKf4uQVrlPVF3MLs7jV1K9l2bpPGdCzk+9Aed51vGxgPROljh07Mn78eN55553k2EKmdevQoQP9+/dvMmx5OoHGPnILgZIuCGD3MNqJi/3MKcP53t3+Y6Un1vFKt74xhdCxY0cGDx4cdRgmpnwLBRH5DGdOA6Hp3AaCM1xRSQ0b6h1Ge0DPTk1eZypRuCS2a7UGY0xc+BYKqmo9jDwmDOudl4u3zdFgjImzIM8ppLxiqapNQ9aMN7mcKp8ATQsXa0oyxsRNkJzCs56fOwB9gQbg8FAiKmLe5HKuTUzGGBOFIL2PmkwbJCJDgOKadLSAUiWXjTGmWGQ885qqvi4ix4QRTCmypLMxJk6C5BQu87wsw5l5bVNoEZUQSzobY+ImSE3B2wtpF06O4c/hhFNaLOlsjImbIDmFGwsRiDHGmOi1OEqqiLwoIp09r7uIyPPhhmWMMSYKQZqPurkzrgGgqh+LSPcQYypZlnQ2xkQtSKHwpYj0TjysJiJfpemwFyYPLOlsjImDIIXCL4BXROQv7uuRwOTwQipNzZPOVmswxkQhSKJ5rvvA2rE4g+H9XFU/zGWnIvJz4Ec4NY7FQDXQE3gU6Aq8DkxU1R257KdYWa3BGBOVIIlmAcYAQ1T1aaBcRIZmu0MROQj4GVCpqkcAbYCzgWnA7araH/gYOD/bfRS7CcN6M3PKcGZOGZ5y/CRjjAlLkDmafw8MB85xX38G5DqTdVtgbxFpC5QD64CTgMfc9+8HzsxxHwXjnWfBGGOKWZCcwjBVHSIib0Cy91H7bHeoqmtF5L+B1cDnwAvAQuBfqrrLXW0NkHI0ORGZjJvT6N07+iaVfM6zYIwxUQtSKOwUkTa4PY5EpBvQmO0ORaQLMBZntNV/AbOBU1OsmrKHk6reA9wDUFlZGXkvqHzNs2CMya/Zb89mzrtzkq+r+lUx/tDxEUZUHII0H/0P8ATQXUR+BbwC3JzDPkcDK1V1k6ruBB4HvgF0dpuTAA4GPshhH8aYEjfn3Tk0bG4AoGFzQ5MCwvhrsVBQ1YeAK4H/wmn7P1NVZ+ewz9XAsSJS7iaxTwaWAS8B49x1zgNqctiHMcZQ0bWC6WOmU9G1IupQika6OZq7el5uBB7xvqeqm7PZoarOF5HHcLqd7gLewGkOehZ4VESmusvuy2b7xhhjspcup7AQp11fUrynQL9sd6qq1wPXN1v8LpB1V1djjDG58y0UVLVvIQMxxphCsAR0eoFmXhORM3CGtwCoVdVnwgvJGGOCy/Qin0hAV3StSCairVDYLcgTzbcAl+Akg5cBl4jIf4UdmDHGBJFNLyNLQPsLUlOoAgapaiOAiNyPkwi+JszAjDGFVczNKomLfPXc6qhDKXpBnlMA6Oz5+SthBGKMiZb16zcQrFD4L+ANEZnh1hIWktvDa8aYmPI2qzRsbqB6bjXVc6uZ/XYujyYVViLuRAHX0nLTVJChsx8RkVrgGJzuqVep6vqwAzO7eedW8LJ5FkxYqvpVJX8upmSsN+6KrhXJ137LYXdhkVgv1+Ms5mY4SP/w2pBmi9a43w8UkQNV9fXwwjIJfgPs2TwLJkzjDx2fvJAVUzu9N+4gy8Mo/Iq9d1O6mkIdsBTY5L72PsSmOENdm5D5Dbhns7OZKBT6Ljjs/YVV+BVz4jtdTuHfgU9whreeDnxbVU90v6xAiNjYQQclJ+BZtu5TaurXRhyRKQWFTkZb8rvw0j3RfDtwu4j0xZlgZ56IvAfcrKr1hQrQpNZ8TmdjsuG9E080ebSk0HfB3v3lu/0/F62t1pQQJNG8UkRqgL2BicChgBUKRebh+aub1CasuclA0/bv5glYL2/Pnage+Ioq+e13Mfb+7uo21FG3oa7JMm+sfgVoJhf2QuUq0iWa++HMnTwWeB94FPiVqn6R9yhM6Grq17Js3acM6NnJktSmicSduJ90PXcKKarkd7qLceJ35y04/Ho9NZfNhb0QtbR0NYXlwCKceQ0+BXoDP3amQABVvS20qEwoBvTsxMwpw625yWTEr+dOPsW9G2dLF+NMez1BfHt1pSsU/pPdU2J2LEAsxpgSVezdOFuTdInmGwoYh4lIrrkGy1WYfIlDN85U7f9+OYJ85lea15T8YitEPifQ0Nmm9co112C5ivgqZJNMqota2E1A+b44+7X/B30yOhuJY6jbUAdAZY9K33W9+zus62E57TcdKxRMzrkGy1XEUyGbZJr3ugl7f2Ekv4PkTvKZX/HGXNmjMqNC9KqhV+UlhlTS9T66RFXvEJHjVPXV0CIweVHIp5u9TUaJWoKJh1TPHYTZJNP8bj3RiynsJqBCJL/DFtdjSFdTqAbuAH4LNB8HycSId3ykoE04iUIkm4u6t8loQM9OvuMzmdwFaWv28jZDhD3wW1y6qpr8SlcovCkiq4BuIrLIs1wAVdUjQ43MBJbp083ei3i2F/VEk5EJV6qHodLxa4YI48GvuN7pmtyk6310jogcADwPnFG4kEzY/AbZM/HU0sNlQcRl1FO/mk8heviYYNImmt15E44SkfY4w1sANKjqztAjM8a0On41nzB7+JjMtNj7SEROAP4XWIXTdNRLRM5T1ZdDjs1EoBiH47ZnJYpLSzUfa5aKVpAuqbcB31LVBgARORR4BDg6252KSGfgXuAInKemfwg0ADOBPjgF0HdV9eNs91HKsr2wZ5OwjoPW8KxE0GYVY8IWZI7mdokCAUBV3wba5bjfO4C5qnoYcBTwJnA1ME9V+wPz3NcmQ7nMszBhWG9mThnOzCnDi66baSLxXWxxJ3jnDfCy5hNTaEFqCnUich/wgPv6XGBhtjsUkU7ASGASgKruAHaIyFhglLva/UAtEN4TGq1U855IuXQ9DZs1+zSVj4RyUIVI5AYZMsLET5BC4SLgYuBnODmFl4Hf57DPfjhTfE4XkaNwCphLgB6qug5AVdeJSPdUHxaRycBkgN69S/cCEkQ+up6GqTU0+xSjQiRygwwZYeIpyCQ723HyCvkaKrstzsNwP1XV+SJyBxk0FanqPcA9AJWVldrC6iUtiq6nmeYzoh4iI0htpbXVaAqRyLVkcfGKYuyjNcAaVZ3vvn4Mp1DYICI93VpCT2BjBLGZHHhrIvNXbmb+ys3Ji2lcL6Te2opfzOnW8a4XRNznDTCm4IWCqq4XkfdFpMJNYJ8MLHO/zgNucb/XFDo2kxtvzaT5+EiJ9+MoUVtJF3OqdVKt1xKbN8DEXdpCQUTaALeo6hV53u9PgYfch+LexRlnqQyYJSLnA6sB+08pYpkOveHH7846jEH5gsTcvEnOb710TU5xmDfAGD9pu6Sq6pfA0ZKYgzNPVLVeVStV9UhVPVNVP1bVj1T1ZFXt737fnM99muLk7arZsLkhWUAkmnQg3kl0yLxrsDFRCtJ89AZQIyKzga2Jhar6eGhRmVbriufv5uUPXgBgm+6iZ9tvAE0H1vPeZa9q/yn779Mn5Z2136B8+XoqO2h3Xu/+vMu8SfTEOs7x7JVVPMYUQpBCoSvwEXCSZ5kCVii0cs0vdvlIFr/8wQts09WUS2/KOqyjXfk/99jf6x8/R9tO9XTaux3bdDUfbg2+z3w9lR20O6/fcu9nvOts276LDzOOxpjCCdIl1Ro+S1Dzi10+k8Xl0pv51X9ucufv3V/XHkv5su1GBnT7OsvW9eYrXw4NvO185TKCducNsp53nWHTbbJDE29BBsQ7FPgDzsNlR4jIkcAZqjo19OhMZFIlVMN8Orpdl/mUf9XJF7TZvI4BXb/O9DHTQ3l+oRAzx/mNZfSFvE8H7ZV8bUNEm7gJMvbRH4FrgJ0AqroIODvMoEz8eMdUCiOx600oh/3UayGS1H5jGXXQXsmaT1W/qmRBYE/6mrgIUpctV9UFzTog7QopHhNTYT0d7TfHb0vrr2r/qXtxDTb7m/fOfVX7T+l58FBm/ijzobVS1QD8HkBLdTzemo899WviKEih8KGIHIKTXEZExgHrQo3KlIRMx+Dxvr+t7G22lb1N9dy3ku81v8B+3OZlPmmzgOq5nZrMXfyFvA9tsou5+SQxuTyA1tqGzzCtQ5BC4WKcsYYOE5G1wEqckVKNyUmmd8re9b917zQ+abMA8L8wf9JmgVMAcHiTuYuHTf+3nOL21gByeQDNBgQ0cRSk99G7wGgR2QcoU9XPwg/LmPS6fDmSLl+OZPqY4WkvzB20V9omKb8npgs1RlHUAwIa01yLiWYR2U9E/gf4K1ArIneIyH7hh2ZM+PyemPZbbkxrF6T56FGcORQSde5zcabNHB1WUCZ+Mr1zbp7YBaie22mPLplBpHtiON1+v5D3afyiZ4tPN/uNReRdHqTrqHffuXQxtVyDiVKgJ5pV9SbP66kicmZYAZl4aF4IeBO1fm343s941wfnSd5l6z6lcXtPvtI2+MNoQZ4Y9vImgg8u/xo7dx0FXzYdynub7qJ8r+APkQVNiHv3nU0X00ThN3+lM+zXsL5dW0WuwQq54hLkP+MlETkbmOW+Hgc8G15IJg6a97LxJmr92vC9n/Gun7wo7AAExg4M/lxANl1hU3UF9V6Yyvdqm9H4Q5kkxLOdUtNbwA3r2zV54WwNuQZLqBcX30JBRD7D6YYqwGXAg+5bZcAW4PrQozORyuYCl+ozUcwA15w3huq5/k8wR/WEcRx+R2GyhHrx8C0UVHXfQgZiiot3UvZinD0s1cU/27mLwy5I8jXqqzFBBGpYdcc76uNd34bOLl3ei2Uxzh7md/HP5gnjbAuSoPI16qsxQQUZEO9PwJHAUqDRXWxDZ7dCQXvPeC+exTh7WD6Hl8hmW5kMLJivUV+NCSpITeFYVR0QeiQmcrn2njEtCzpPgzFRCVIo/F1EBqjqstCjMQWXqnaQaXI5DsM/xyGGIHJNKKfKL1iXT5NPQQqF+3EKhvXAdpzeSKqqR4YamQkslyEZcq0dhN2mXiwxFIJffiGOXT4znbOiFAu2uB5zkELhT8BEYDG7cwomRrwX9mwSv9n2rU/sJ+okcxxiKIR0+YW4dfn0FlRBmsniWLCFLa7HHKRQWK2qT4UeicmJ31ANCX4zgcW9ucUUr0RBlen63ln+ID530GGIW2EOwQqFt0TkYeBpnOYjwLqkFpvmTygntObmFlOcrBtutIIUCnvjFAbf8iyzLqkx5n2wzLssl2YiUzxSDSAI8bnjbqkW0Jq74eYr13Lj00sBuP7bh+c9xiDzKYTSEV1E2gB1wFpVPV1E+uKMyNoVeB2YqKo7wth3sQqSUPa767caQWnwa7v3DgiY6jNBC4vmFyk/fhe8bGoBcWpKyjU5HCTX4t2Hd3BE799w/srNDOvbNcejSS3Iw2vTcafi9FLVH+a470uAN4HEmTMNuF1VHxWRu4DzgT/kuI9WJUhCuVSSriY1vy6vfhfzTJtnvBe1dPwueJnWAuLWlJSP5HBLuRbvPryDI3r/honlYQjSfPSM5+cOwFnAB7nsVEQOBk4DfgVcJiICnARMcFe5H7gBKxT2kOkY/6b1yvbJaK9skrqZJpBzEVVTUroaQSYJ8ebbCfK38u7Dq1CDJgZpPvqz97WIPAL83xz3+xvgSiAx6N5+wL9UdZf7eg2QshgUkcnAZIDevaNvH41KqfTNN6nl68nouN2Jx0WQGkGQ313zmlUxPMUefKaR3bich08AAA/VSURBVPoDWZ81InI6sFFVF4rIqMTiFKvu0WQFoKr3APcAVFZWplynFFgzUWnL111j8zvxQrffZ1LbKbSWuosGrcUUsmaVD0FyCol5FRLWA1flsM/jgDNEpAqnOaoTTs2hs4i0dWsLB5NjE5UxJjOFrjVkU9sJs9BK1zMozoVXvgVpPsrrvAqqeg1wDYBbU7hcVc8Vkdk4s7o9CpwH1ORzv8aY9Ardfp9pbSfsQsuvZ1CpDWKYbua1tL9tVV2d51iuAh4VkanAG8B9ed6+MSYLmfatD0u6pi6vXGoQUSZ44yJdTeFZdk/HmaBAN6A70CbXnatqLVDr/vwuEHxGd2NMQWQ6jlEh+MVgyfLcpZuOc6D3tYj0wbmbHw3cHGpUxphYiVuyNJcutnEdnTQueYsgieb+wC+AYcCtwM9UdWfYgRlH0NnQjDH+eYcgTwnnejFO1ZwVdJtxylukyykcgVMYHA78H+B8Vf2yUIEZh82GZkxwfsnyIE8J5+t5D6+g24xT3iJdTeGfwPs4uYWhwFDnwWOHqv4s3NBMgg1kZ0zuwkwix+minqt0hUKuYxsZY4wpMukSzfcXMhBjjMm3uCRvi0lZ1AEYY+LJe0EtRmMHHVRUYw7FRTZjHxljWrk49YbJVmtq5y8kKxRiyLqhmqjZBbV0tdh8JCKHisg8EVnivj5SRP4j/NBKV6IbKtiw2MaYwgpSU/gjcAVwN4CqLhKRh4GpYQZW6qwbqjEmCkESzeWquqDZsl0p1zTGGFPUghQKH4rIIbhzKojIOGBdqFEZY4yJRJDmo4txZjo7TETWAiuB74calTHGmEgEmWTnXWC0iOwDlKnqZ+GHZYwxJgpBRkndC/g3oA/QNjH+kar+Z6iRGWOMKbggzUc1wCfAQmB7uOEYY4yJUpBC4WBVHRN6JMYYYyIXpFD4m4gMVNXFoUdTwuwpZmNMHKSbZGcJ0OiuUy0i7+I0HwmgqnpkYUIsDTaZjjEmDtLVFA4CBhUqEGNPMRtjopeuUFipqu8VLBJjjDGRS1codBeRy/zeVNXbQojHGGNMhNIVCm2Ajjg5BGOMMSUgXaGwLowH1ESkF/C/wAE4iex7VPUOEekKzMR5SG4V8F1V/Tjf+8+7uumw+LHdrweOg8rqQB8tmR5HOfyOjDGFlW5AvLBqCLuAf1fVrwPHAheLyADgamCeqvYH5rmv42/xY7De7a27fnHTi18LSmbehBx+R8aYwkpXUzg5jB2q6jrcUVZV9TMReROnp9NYYJS72v1ALXBVGDHk3QEDofpZmH5acpG3FgBQ1a+K8YeO3+OjJdPjKMXvyBgTP741BVXdHPbORaQPMBiYD/RwC4xEwdHd5zOTRaROROo2bdoUdohZ89YCGjY3NCkgjDEmriKbo1lEOgJ/Bi5V1U8TA+21RFXvwRnKm8rKSg0vQn9NagGywfk+txpkA1W6D4n6QKIWUD23OuVnW3UeIc5KPcdR6sdv0oqkUBCRdjgFwkOq+ri7eIOI9FTVdSLSE9gYRWxBeJ8+9qqT7dTJdubMrfa94Bf0yeVS/+dvfvwJ773ifP/qiN25jlL6vSRyPAcMjPb4S/38jKmCFwriVAnuA95s9qzDU8B5wC3u95pCx5aJZC4g0UY+Zjqzp5/AHNmafN/vgl+wPEJc/vmj4j1+r6+O2H0BKtUcRxxyPKV+fsZUFDWF44CJwGIRqXeXXYtTGMwSkfOB1cCeWdmYG09HxmtHiDJx7L37SvzDRf3Pny/Z3Fkmjj9OMZndsj0/W9vvPUbHE2SO5rxS1VdUVVT1SFUd5H7NUdWPVPVkVe3vfg890d0qebt/HjDQOblaizh2bY1jTKWgtf3eY3Q8kSWai02uCeKGzQ1Up8k15FU+7o792uMTorqTiWPNJ5OY8nVHGKM7y5QKEV+czoV8HG9MjscKhYBySRB71y2ah9T82uPB2n9zka929Ey3k6pZMZ1cL3J+8WUaRyHlcsytKD9ihUIGsk0Qjz90fMoH12LPr8YRhzuzYpavO8JMtuO9aAVpVszHRS5VfJnGUUi5HnNM7vRzZYVCUJ+th62bdv/B/e4i1i9ueR0TnlzuRPPRBJDp37+Q50umzYrei1yqOLP9fbUUR5RNY63kwp6Lgieai9bWTbDD6W7qmwgaOG73Rag1JL+KUS6J9lyTfZn+/YvlfPGLM6zkaIySrqXIagqZaL8PTEpzF1FZvfuOxruO351PkDuiMO6a4p6k9Apy5+/XDTcbLd0dp+P398/X+kGE8bdNF6ff7yvXfEGh7tjTnV/e4/GKWw0wz6xQKARvW+V7rzhfix9r+mStdznsPpHCSGCFlQT067GUywUiSBt0GO3U3m0UU+LQ71wLO6nb/HceRr4gjKZBv3PHL3a//1OvfJw7ERYqVigUSuLOx3tyep+sbX7Cw+4TIYy7pjCSgH4Xn1wvEEHu/PP9kFoYd/GFkupcCzup6/19hSWX8zPdzVWqc8fveNL9n6b6bDbnTsQ3JFYopNHk2QR2UEH73Dfqd7LlciLls5tfrhfXbD+fTbNH4m6qEF0b49AcEKSpw7u8EBfqoPL1t8olSZ2Pm6ts/k+bN0O1dP4030eqc+85d7qZU2/JLP4ArFBIY84//0TD1rVU0J6KHTupatdl95uFvCC1JM7d/ILKtJnMe4y5HnNLf8tc79yCFtotFTxBmjqy+V2EXeDl82/VkkI/L5DJuZNYPxFTkPPC79xL/BwCKxTS2bqJih07md6uF7Qjf/+EYQh7jJ9CyOROLl93wUH+lrk2BwQptIMWPJk0dQRRiKaKQtdYCpWkzvTcgcybayNoxrRCoSWJHkdeQU/yMGsTcX4yNB8KVRPL5oIVtJdN82NIV2inazLwbi/fv4tizp1ELdvCLpPzoqXPhsAKhbCEXZtoDU1GfuJYE0sI2ssml2PwWzduvwuTuXydFyGeC1YoNEtMzWYLczqWw74H5JZcLkSVOV93GdncfeRyJ9tSG3acEqTNBY0tl2OI8/FDeHeruZyfmX42qpxgEZwXVig060Y5Z+dGGra2o2LfA6igPVW6TzRx5ftk9rvLyObuI5c7We/7hexHb4LJJHGaz7vVXM7PTD8b55poDIhqJNMc50VlZaXW1dXltpHEHat7x109o9JZPKluj/cKJpcnoItJazueYmd/j5IhIgtVtTLVe1ZTiKMgzzK0Bq3teIqd/T0MNiCeMcYYDysUjDHGJFmhYIwxJskKBWOMMUlWKKSyY+vubnnGGFNCrPdRc/t0c74r1ofZGFNyrFBobt8DnK8x06OOxBhjCi52hYKIjAHuANoA96pq/gcM95jNFubIVpjr9M9u2NxARdeKMHdpjDGxFaucgoi0Ae4ETgUGAOeIyIAw9zlHttLAjuTriq4VVPWrCnOXxhgTW3GrKQwFlqvquwAi8igwFliWz51Mm/Vt3tq2DoAG/YIK6cB0ay4yxph41RSAg4D3Pa/XuMuSRGSyiNSJSN2mTZty3mGFdKDqwONz3o4xxrQGcaspSIplTUbsU9V7gHvAGRAvm51c9d2ns/mYMca0enGrKawBenleHwx8EFEsxhhTcuJWKPwD6C8ifUWkPXA28FTEMRljTMmIVfORqu4SkZ8Az+N0Sf2Tqi6NOCxjjCkZsSoUAFR1DjAn6jiMMaYUxa35yBhjTISsUDDGGJNkhYIxxpgkKxSMMcYkiWpWz3/FgohsAt7L8uP7Ax/mMZwo2bHEU2s5ltZyHGDHkvBVVe2W6o2iLhRyISJ1qloZdRz5YMcST63lWFrLcYAdSxDWfGSMMSbJCgVjjDFJpVwo3BN1AHlkxxJPreVYWstxgB1Li0o2p2CMMWZPpVxTMMYY04wVCsYYY5JKslAQkTEi0iAiy0Xk6qjjyYSI/ElENorIEs+yriLyooi8437vEmWMQYhILxF5SUTeFJGlInKJu7wYj6WDiCwQkX+6x3Kju7yviMx3j2WmOxx8URCRNiLyhog8474uymMRkVUislhE6kWkzl1WjOdYZxF5TETecv9nhod1HCVXKIhIG+BO4FRgAHCOiAyINqqMzADGNFt2NTBPVfsD89zXcbcL+HdV/TpwLHCx+3coxmPZDpykqkcBg4AxInIsMA243T2Wj4HzI4wxU5cAb3peF/OxnKiqgzx9+ovxHLsDmKuqhwFH4fxtwjkOVS2pL2A48Lzn9TXANVHHleEx9AGWeF43AD3dn3sCDVHHmMUx1QDfLPZjAcqB14FhOE+btnWXNznv4vyFM+PhPOAk4BmcaXKL9VhWAfs3W1ZU5xjQCViJ2zEo7OMouZoCcBDwvuf1GndZMeuhqusA3O/dI44nIyLSBxgMzKdIj8VtbqkHNgIvAiuAf6nqLneVYjrPfgNcCTS6r/ejeI9FgRdEZKGITHaXFds51g/YBEx3m/TuFZF9COk4SrFQkBTLrF9uRESkI/Bn4FJV/TTqeLKlql+q6iCcu+yhwNdTrVbYqDInIqcDG1V1oXdxilVjfyyu41R1CE5z8cUiMjLqgLLQFhgC/EFVBwNbCbHJqxQLhTVAL8/rg4EPIoolXzaISE8A9/vGiOMJRETa4RQID6nq4+7iojyWBFX9F1CLkyfpLCKJ2Q2L5Tw7DjhDRFYBj+I0If2G4jwWVPUD9/tG4AmcArvYzrE1wBpVne++fgynkAjlOEqxUPgH0N/tTdEeOBt4KuKYcvUUcJ7783k47fOxJiIC3Ae8qaq3ed4qxmPpJiKd3Z/3BkbjJAJfAsa5qxXFsajqNap6sKr2wfnf+H+qei5FeCwiso+I7Jv4GfgWsIQiO8dUdT3wvohUuItOBpYR1nFEnUSJKHFTBbyN0+77i6jjyTD2R4B1wE6cO4jzcdp85wHvuN+7Rh1ngOMYgdMEsQiod7+qivRYjgTecI9lCXCdu7wfsABYDswG9oo61gyPaxTwTLEeixvzP92vpYn/9SI9xwYBde459iTQJazjsGEujDHGJJVi85ExxhgfVigYY4xJskLBGGNMkhUKxhhjkqxQMMYYk9S25VWMMSKS6P4HcADwJc7QAwDbVPUbkQRmTJ5Zl1RjMiQiNwBbVPW/o47FmHyz5iNjciQiW9zvo0TkLyIyS0TeFpFbRORcd66FxSJyiLteNxH5s4j8w/06LtojMGY3KxSMya+jcOYiGAhMBA5V1aHAvcBP3XXuwJmb4Bjg39z3jIkFyykYk1//UHc4YxFZAbzgLl8MnOj+PBoY4Az/BEAnEdlXVT8raKTGpGCFgjH5td3zc6PndSO7/9/KgOGq+nkhAzMmCGs+MqbwXgB+knghIoMijMWYJqxQMKbwfgZUisgiEVkGXBh1QMYkWJdUY4wxSVZTMMYYk2SFgjHGmCQrFIwxxiRZoWCMMSbJCgVjjDFJVigYY4xJskLBGGNM0v8HCNgby+GW2F0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret.plot(y=[\"K+K_KK\", \"Kp+Kp_KK+Kp_PP\", \"Kpp+Kpp_PP\"], legend=True, step=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 }