{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Reversible\n", "\n", "This is for an integrated test of E-Cell4. Here, we test a simple reversible association/dissociation model in volume." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "from ecell4.prelude import *" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Parameters are given as follows. `D`, `radius`, `N_A`, `U`, and `ka_factor` mean a diffusion constant, a radius of molecules, an initial number of molecules of `A` and `B`, a ratio of dissociated form of `A` at the steady state, and a ratio between an intrinsic association rate and collision rate defined as `ka` and`kD` below, respectively. Dimensions of length and time are assumed to be micro-meter and second." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "D = 1\n", "radius = 0.005\n", "N_A = 60\n", "U = 0.5\n", "ka_factor = 0.1 # 0.1 is for reaction-limited" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "N = 20 # a number of samples" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Calculating optimal reaction rates. `ka` and `kd` are intrinsic, `kon` and `koff` are effective reaction rates." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "import numpy\n", "kD = 4 * numpy.pi * (radius * 2) * (D * 2)\n", "ka = kD * ka_factor\n", "kd = ka * N_A * U * U / (1 - U)\n", "kon = ka * kD / (ka + kD)\n", "koff = kd * kon / ka" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Start with no `C` molecules, and simulate 3 seconds." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "y0 = {'A': N_A, 'B': N_A}\n", "duration = 3\n", "opt_kwargs = {'legend': True}" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Make a model with effective rates. This model is for macroscopic simulation algorithms." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "with species_attributes():\n", " A | B | C | {'radius': radius, 'D': D}\n", "\n", "with reaction_rules():\n", " A + B == C | (kon, koff)\n", "\n", "m = get_model()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Save a result with `ode` as `obs`, and plot it:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEGCAYAAABiq/5QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dd3xb9dn//9clyfKIHY/YcfZeJGQ5JoSywvwFWgKU0RQKtDdfUkpboLQUaG9KB6VQWrihAwilkEIJpKFhQ5owy8zeMdk7jh0n3kuSr98fko0JHopt+djS9QQ9dM7ROUfvEyWXjz865/MRVcUYY0zscDkdwBhjTOeywm+MMTHGCr8xxsQYK/zGGBNjrPAbY0yM8TgdIByZmZk6ZMgQp2MYY0y3smLFikOqmnX08m5R+IcMGcLy5cudjmGMMd2KiOxqark19RhjTIyxwm+MMTHGCr8xxsSYbtHGb4wxTqitrWXbtm1UVlY6HaVFSUlJDB8+HK/XG9b6VviNMaYZ27ZtIy0tjdGjR+Nydc0Gkrq6Og4ePMjWrVsZO3ZsWNtE9EhEJE1EFohInohsEpGTRCRDRBaLyJbQc3okMxhjTFtVVlaSnZ3dZYs+gMvlIjs7m8rKSvLz88PbJsKZHgLeVNUxwERgE3A78JaqjgTeCs0bY0yX1JWLfj2Xy4WI8NJLLxEIBFpfP1JBRKQncBrwBICq1qpqMXAhMDe02lzgokhlWPH6k3z6rz9EavfGGNOlVFdXU1VV1ep6kfxRNgwoBJ4UkVUi8jcR6QFkq+oBgNBz76Y2FpHZIrJcRJYXFha2KYBsXMioDf9HTXXX/mLGGGNasnDhQkSEvLy8FtcTkbD2F8nC7wFygEdUdTJQwTE066jqHFXNVdXcrKwv3XEclrgTvk06Zax7659t2t4YY7qCefPmccopp/Dcc891yP4iWfj3AntV9dPQ/AKCPwgOikhfgNBzQaQCjDvlQvZLbxLWPhOptzDGmIgqLy/nww8/5Iknnuiwwh+xyzlVNV9E9ojIaFX9DDgL2Bh6XAPcG3p+KVIZXG43uwZ9nZN2Pcq+7RvoP2xcpN7KGBPlfvXKBjbuL+3QfY7t15O7Lmi5Lr344ovMmDGDUaNGkZGRwcqVK8nJyWnX+0b66+ofAv8UkbXAJOAeggX/HBHZApwTmo+Y4edeT0CF3Usei+TbGGNMRMybN49Zs2YBMGvWLObNm9fufUb0Bi5VXQ3kNvHSWZF838Z69x/K6h7TGLn/JXy19xPnje+stzbGRJHWzswjoaioiLfffpv169cjIgQCAUSE3//+92F/kduUrn+BakeYfDWZFLP+3X85ncQYY8K2YMECrr76anbt2sXOnTvZs2cPQ4cO5YMPPmjXfmOi8B8//VIKyMC18imnoxhjTNjmzZvHxRdf/IVll1xyCc8++2y79hsTffV44rxsG/h1Ttz9BHu2rmPgiPFORzLGmFa9++67X1p24403tnu/MXHGDzDyqzfhx8X+Rf/ndBRjjHFUzBT+zD6DWJN2FscXvEppcZHTcYwxxjExU/gB0s64kR5SzcbX/+p0FGOMcUxMFf6Rk05lU9xYBm55moDf73QcY4xxREwVfoCqnNn014Osfbtjbn02xpjuJuYK/4SzrySfLLzL7U5eY0xsirnC74nzsnPkVYyrXctny992Oo4xxrTI7XYzadIkJk6cSE5ODh999FG79xlzhR9g/MybKCaZqrd/73QUY4xpUWJiIqtXr2bNmjX87ne/44477mj3PmOy8PdISWPToCuZVPkx29d/2voGxhjTBZSWlpKe3v5hymPizt2mjL3oVioemsuRRffC8QudjmOM6ereuB3y13XsPvuMh/Na7qC4qqqKSZMmUV1dzYEDB3j77fY3UcfkGT9AakYWa/tdyqTSd9iztYM/TGOM6SD1TT15eXm8+eabXH311ahqu/YZs2f8ACMvvB3fI89z4NV7GHhz+/u4NsZEsVbOzDvDSSedxKFDhygsLKR37yaHKw9LzJ7xA2T2GciarJlMPrKIfds3OB3HGGNalJeXRyAQoFevXu3aT0yf8QMM//pd+B97hQMv/oL+t7zgdBxjjPmC+jZ+AFVl7ty5uN3udu0z5gt/Zr/BfNxvFifuf5odGz5l6LgTnY5kjDENAoFAh+8zppt66o297BeUSyIlr93ldBRjjIk4K/wEr/DZOOQ7TKr8mLxlS5yOY4wxEWWFP2TCpbdxiDQCi3+J1tU5HccYYyLGCn9IUnIq2467gXG161i9pH3jWRpjTFdmhb+RKV//ETtdg+j98a+pqa50Oo4xxkSEFf5GPHFeyqbfTX89yMr5v3M6jjHGREREC7+I7BSRdSKyWkSWh5ZliMhiEdkSem5/j0MdaPxpF7Iq6StM2DaHQ/t3OR3HGBPj8vPzmTVrFsOHD2fs2LGcf/75bN68uV377Iwz/jNUdZKq5obmbwfeUtWRwFuh+S4l6+v3E4ePHfNvczqKMSaGqSoXX3wx06dPZ9u2bWzcuJF77rmHgwcPtmu/TjT1XAjMDU3PBS5yIEOLBow4nhX9vskJxW+w6dNFTscxxsSod955h7i4OK6//vqGZZMmTeLUU09t134jfeeuAv8REQUeU9U5QLaqHgBQ1QMi0mRPQyIyG5gNMGjQoAjH/LIJV/yWA39cROKiH1M76XS88QmdnsEY03Xct/Q+8g7ndeg+x2SM4bapzbcsrF+/nilTpnToe0Lkz/hPVtUc4Dzg+yJyWrgbquocVc1V1dysrKzIJWxGj5Q0Ck67hyF1e1jxrN3Ra4yJHhE941fV/aHnAhFZCEwFDopI39DZfl+gIJIZ2mPimZezYtWz5Ox8gj1brmDgyIlORzLGOKSlM/NIGTduHAsWLOjw/UbsjF9EeohISv00cC6wHngZuCa02jXAS5HK0BEGX/kwNeKl9F8/oC4CnSUZY0xzzjzzTGpqanj88ccbli1btoz33nuvXfuNZFNPNvCBiKwBlgKvqeqbwL3AOSKyBTgnNN9lZfYZRN74WxlXu5ZlC+53Oo4xJoaICAsXLmTx4sUMHz6ccePG8ctf/pJ+/fq1a7+tNvWIyGXAm6paJiL/C+QAd6vqypa2U9XtwJfaRlS1CDirjXkdccLFN7Fmy2tM2PhH9mw5z5p8jDGdpl+/fsyfP79D9xnOGf+doaJ/CvD/EbwE85EOTdHFictF/6ufoEa8VD3///D7ap2OZIwxbRZO4a9v2P4q8IiqvgR4Ixepa8rsN5htU3/NKP9mlj1zp9NxjDGmzcIp/PtE5DHgcuB1EYkPc7uoM+X8a1mecha5Ox8nb/lbTscxxnSCum7QTfuxZgyngF8OLAJmqGoxkAHceuzRosPI/5lDoasXaa9+l5Ki9t02bYzp2pKSksjPz+/Sxb+uro78/Hx8Pl/Y27T65a6qVopIAXAKsAXwh55jUmp6JgcveJwhL32djU9cw8SfvI64YvIXIGOi3vDhw1m9ejX79+9HRJyO0yyfz8fu3btRVVxh1KNwruq5C8gFRgNPAnHAM8DJ7czabY3Kmc4nebcwbfP9fDLvN0y70u7sNSYaeb1eEhMTWbx4MampqWEVVaeUlZWRnZ1NYmJiq+uGcxQXAzOBCmi4GzelXQmjwImzfsaqpJOZsvkhNn78htNxjDERcvzxxzN9+nQSEhIQkS75cLlcjBw5kpkzZ4b1m0k4XTbUqqqGOlqrvws35onLxbDr/kH+w6fSZ9FsDvRbQt/Bo52OZYzpYCJCTk4OOTk5TkfpMOGc8c8PXdWTJiLXAUuAx1vZJiakpmei35xHnPqo/McsKstLnI5kjDGtarXwq+ofgAXACwTb+X+hqn+KdLDuYtCoSWw//WGG+neQ9+hVaBf+9t8YYyDM6/FVdbGq3qqqP1HVxZEO1d1MPPNylg7/ITnl7/HJEz9yOo4xxrSo2cIvImUiUtroubTxfGeG7A5O/NavWJpxASfte4pP51tnbsaYrqvZwq+qKaras9Fzz8bznRmyOxCXi5zv/Z3VidPI3fBbVv3nGacjGWNMk1pt6hGRQU09OiNcd+OJ8zLq+/PZFjeSMR/+iI2fvOl0JGOM+ZJw2vhfa/R4C9gO2IXrzUhKTqXXdQspdGcx6I1vs3nlu05HMsaYLwjnqp7xjR4jCQ6f+EHko3VfvbIHkHDta5S4etLn5SvYtu4TpyMZY0yDY77/ODQAywkRyBJVevcfilzzClUkkP7CZezYuMzpSMYYA4TXV88tjWZdBEfgKoxYoijSb8ho9nzrJeSZmaTNv5itFz/HiImnOB3LGBPjwjnjT2n0iCfY1n9hJENFk4EjxuO7+lWqSaD3wsvIW7bE6UjGmBgnqup0hlbl5ubq8uXLnY7RLvm7t+B/8gIy6g6z7azHGX+a/ew0xkSWiKxQ1dyjl4dzOediEUlrNJ8uIos6OmC06zNoJAnXLeKguy+j3/oOy19+1OlIxpgYFU5TT1Zo5C0AVPUI0DtykaJXZr/B9LrxbbbEjyN35W188o87rW8fY0ynC2uw9cY3bInIYKDrtw91UT3TejHilkWsSDmDadsfZtmfr6a2ptrpWMaYGBJO4f858IGIPC0iTwPvA3dENlZ0i09IYvLNL/Bxv2uYevgVtv7hLA4X7HM6ljEmRoRzA9ebBC/hfB6YD0xR1bDb+EXELSKrROTV0PxQEflURLaIyPMi4m1r+O7M5XZz0uyHWZ57P8NqP6Pmr6ezdc2HTscyxsSAcL7cFWAGkKOqrwBJIjL1GN7jJmBTo/n7gAdDdwEfAa49hn1FndyvzWbPRf/GTYCB/76QpQsesHZ/Y0xEhdPU81fgJOCbofky4C/h7FxEBgBfBf4WmhfgTIIDuwDMBS46hrxRaeTk0/Dc8AGbE8Yzdf2vWP7QN2w0L2NMxIRT+E9U1e8D1dBwVU+4zTP/B/wUqD+F7QUUq6o/NL8X6N/UhiIyW0SWi8jywsLov1E4o3d/xt66mI8HzWZK8WKK/jiNLavedzqWMSYKhVP4fSLiJnQlj4hk8Xkhb5aIfA0oUNUVjRc3sWqTVwip6hxVzVXV3KysrDBidn9uj4eT/ud+Np77NF6tYciLF/HJUz8j4Pe3vrExxoQpnML/MLAQ6C0ivyXYM+c9YWx3MjBTRHYCzxFs4vk/goO21/cRNADYf6yho93xJ19Awo2fsrbnaUzb+Re23nsyu/JWOh3LGBMlwrmq558Em2t+BxwALlLVf4Wx3R2qOkBVhwCzgLdV9UrgHeDS0GrXAC+1MXtUS83IIudH/2b5lN+T7d9L33nn8PHcn+GrrXE6mjGmm2tpzN2M+gdQAMwDngUOhpa11W3ALSKylWCb/xPt2FdUE5eL3Au+i//6T1ifcjIn7fgLe+6dSt5SG+/eGNN2zXbSJiI7CLa/N9kur6rDIhmssWjopK0jrFz0NP0/votsilia/lVGXvFH0rP6Oh3LGNNFNddJm/XO2c1UlBWz7p8/Y8qB56iSeDaOvJ6cS2/DG5/gdDRjTBfT5t45QxvPFJE/hB5f6/h4Jlw9UtKYdv1f2T9rMTsSxjFtywMU3DuJlW8+ZTd+GWPCEs6du/cSvPt2Y+hxk4j8LtLBTMsGHzeFibcvYe30v+MXDzmf3MTWe6ay7v2F9gPAGNOiVpt6RGQtMElV60LzbmCVqk7ohHyANfW0JuD3s+KVRxi05iH6UMgG73j0tJ8y7itfQ1zHPKyyMSZKtKupB0hrNJ3aMZFMR3F7PEy9+Iek376WT0bfRlbtPo5fchWf/e5k1rzzL/sNwBjzBeGc8X8TuJfg9fcCnAbcoarPRT5ekJ3xH5vqqgrWvPxnBm+aQx8OscM1hEMTZjPxvGvtS2BjYki7ruoRkb7ACQQL/6eqmt/xEZtnhb9tamuqWfPGE2SufYyhdbsoJJ2tg7/ByPN+QGafgU7HM8ZE2DEXfhHJaWmHqtppfQhY4W8fratj7XsvIJ8+xoTqZdSqm7WpZ5A47VrGTpth3wMYE6XaUvjrgA1AfdeYjW/kUlU9s8NTNsMKf8fZvXk1Bxb/ieMKX6cnlex29Wf/0EsZcda1ZPYb7HQ8Y0wHakvh/xFwCVBCsJO1hapaHtGUzbDC3/GqKspYv3guKeufYYx/EwEV1ifm4hv/DcaefjlJyfYdvjHdXZvb+EVkKMFBWC4EdgH3qOrqiKRshhX+yNr12Wr2v/d3hu1/lWyKqNR4NvY8Bc+ESxlzyoUkJPZwOqIxpg3a++XuOII9bF4F/FRV53d8xOZZ4e8cAb+fvE8XUb7iOUYffps0yqnQBPJ6ngTHXcCoky8mJbU9/fMZYzpTW5p6hhEs9hcCewg297yqqtWRDNoUK/ydz1dbw6aPXqV67UJGHH6PDEqpVTefJUykYsjZ9D9hJgNHjHc6pjGmBW39cnctwf7ySzlqpCxVfSACOZtkhd9Zfl8tm1e8Tenql+lX8C6D6vYBsFf6sK/XV/COOothJ8wgNT3T4aTGmMaaK/yeplYO+TWfF/vkiKQy3YInzsvYaTNg2gwA9m5dz74Vr5Kw823GF75G0qF/E/hQ2Bw3kqKsqSSNPJ1hU862ZiFjuijrltm0S21NNVtXvkPJxiWk5X/E8NrP8EqAgAo7PUM5lD4Z99CT6DfuNPoOGmn3DBjTiaw/ftMpqirK2LbqHco+e4+UguUMq95EkgSHiywilT1JY6nKmkjSkFwGjvsKGb37O5zYmOhlhd84wu+rZefGpRR99hGufSvILlvHgMB+XBL8e1dABvsTR1HVayzx/SeQNSKHfkPH4fa01AppjAlHW77cvUlVHxKRk1X1w4gnbIEV/uhSVnKY3Rs+oWz7UjwF68gs38zAwB7coR8GVepln2cgxT2G4es1ioS+Y+k15Hj6DjmOOG+8w+mN6T7aUvhXq+okEVmpqi322xNpVvijX3VlOXu3rObw9lXU5W8gqWQL2dU7yKaoYR2fusl3ZXM4YSBVKUOQXsNI7D2C9IGjyR440noeNeYobbmqZ5OI7ASyQoOxNOyLYF89nTYQi4l+CUnJjJh4Ckw85QvLy0oOc2DbOkr2bMBf8BnxJdtJq9zNqMpVJBbUwqbgegEV8qUXh719qUjsRyBlAK70QSRmDaZn9lCy+g+1biiMCWmxjV9E+gCLgJlHv6aquyKY6wvsjN8cTevqOJS/m8LdeZQf2ELg8E48pbtJrtxHhi+fTD3c0HRUr5QeHHb1otSbRXVCNoEe2bh69sGb1pfE9H6kZA4gvXc/+wFhokZbzvgJ9bs/UUS8wKjQ4s9U1ReBjMaETVwusvoNIavfEGDGl1731dZwcN8OivN3UFm4C9+RPbjK9uOtPEhybQF9ineQcaQYj3x5dLJKjafYlUaZO42quHRq4zMIJGYgSb1wJWfiTc4koWcmiWm9SU7LpGd6ln33YLqVVi+dEJHTgX8AOwk28wwUkWtU9f0IZzOmzeK88fQbOoZ+Q8c0u07A7+fQof2UFOylomg/1cX7CZQVIhWFeKoKia85TEptASlVm0k7UoJXAs3uq0ITKJNkKl3JVHtSqPGk4I9Lpi4uhbqEVCQ+BUlIxZ3Uk7jEVLw9UvEmpZLQI5WEHin0SEkjPiHJ7nMwnSKca+YeAM5V1c8ARGQUMA+Y0tJGIpIAvA/Eh95ngareFert8zkgA1gJXKWqtW0/BGPaxu3xkNlnEJl9BrW6rtbVUV5eQunhAiqOHKSqpJDa8iIC5UXUVRUj1cW4q4vx+Erx+stIrd5PYmUFPagkWSsbLl9tSUCFShKokkRqJIEaVyI+VwI+dyIBdwIBdyIBTyLqSUQ9CRCXiHiTkLhEJC4RlzcBd1wibm8invgkPPGJeLyJxMUnEBefhMebQFx8IvHxCXjjE3G53R3xx2i6oXAKf1x90QdQ1c0iEhfGdjXAmapaHlr/AxF5A7gFeFBVnxORR4FrgUfaEt6YziIuF8k900numQ5DRh/TtnWBAOUVpVSWFVNVepjqihJqK0vwV5XiryqjrrqMupoyqK1EfBW4fBW4/VW4/ZXEBSpJ8JcSV1tIvFbj1RritYZEappspjoWPnXjw0OtxOHHg584fBJHQNwEJI6AeILPrjjqxENd6FldHjQ0ry4PuDyh5zjU5UZCy8XlAZcb3MF1JLSuyx1cLi4PuN24Gi0Xlzv08CBuFy6XB5fL3bCey+1GxIWrYTt3aNqFuDyhZxdulwdcLlwuF66GfdbP1y8LPYs0PMfKb1zhFP7lIvIE8HRo/kpgRWsbafBb4/qBW+JCDwXOBK4ILZ8L/BIr/CaKudzuz39o9B/aYfv11dZQXVVBTVUFtdWV+Kor8NVU4q+pwl9bRaC2kjpfDYHaaupqq1B/TcODQC34a5BALRKohTofrkAtUufDVefDVVeLS/241I+7zodXK3Gr//MHgc+fCeBRP57QdEtNYkdTIAAEBAIIAaBO6p+Dy+pCr9eFpusE/I2m6wi9JvX7+3y6/jUNradHrVsXeqhIw3wwhwS3QwhI/XTj7aXp6UbracP7BvcTfJ/Q8sbrNHp/Gr8mwfmvn/s3xo05sa1/TZoUTuH/HvB94EaCbfzvA38NZ+ci4ib4Q2IE8BdgG1Csqv7QKnuBJu/ZF5HZwGyAQYNa/1XcmFgT540nzhvf0BlendZRE6ihNlBLTaDmC9P1z75ALbV1tfgCvqaf63z4Aj78df7gdF2j6YAPv/obltU/Ahr4fF79+AP1y/346wIE6gIENEBA/Y2mA9RpHXV0/Z4DnFDfMijAea6ObwlvtfCrag3Bdv5j7oZZVQPAJBFJAxYCxzW1WjPbzgHmQPByzmN9b2O6okBdgEp/JRW+ii88Kv2VVPoqqfJXff7sDz43flT7q6kJ1HxhujpQTY2/htq6jikQLnER54rD4/LgdXkbpuPccXjEg8f1+cMtbrwuL0meJNwud8Pr9dNulxu3uBvWrd+32xWc9ogHl7ga1nOJ64vPLteXltc/3BJqpsHVsD9BvvhaaN365UdP1x9v/fL6bY6erl+38fL6+cbriwgN/4Wm69+j8WvB/7+4nog0/YFEQKd0iKKqxSLyLjANSBMRT+isfwCwvzMyGNMRVJUqfxXFNcWU1JRQUltCSU0JpbWllNaUUlpbSllt2ecPXxnlteXBh6+cSn9l2O+V4E4g0ZNIoieRBE8CCZ7gfM/4nvRO6h1c5g4uj3fHE++Ox+v2kuBOIN4Tj9flJd4TT7wruLzh4fr8Oc4dh9cdLO71D7fLvvSNdhEr/CKSBfhCRT8ROBu4D3gHuJTglT3XEBzoxRjHVPurKaou4lDVIYqqijhcfbjh+Uj1EQ7XBJ+Lq4sprilu8cza4/LQ09uTnt6eJMclk+xNJjspu2G6R1yP4HRccDopLokkTxJJcUnB+dB0gjvBCrCJmBYLf6iN/l5VvbUN++4LzA3twwXMV9VXRWQj8JyI3A2sAp5ow76NaZWqcrj6MPmV+eRX5HOw4iAFlQXBR1UBhZWFFFYWUuYra3L7lLgU0hPSSU9Ip19yP8b1GkdaQhpp8cFHanwqqd5UUuNTg8U+vicJ7oRO/ZXdmLZo7c7dgIhMERHRY+y/WVXXApObWL4dmHpsMY35MlXlSM0R9pTtYW/ZXvaV72Nf+T72l+/nQMUBDpQf+NLZucflISsxi6ykLIanDefEvieSlZhFZmImmYmZZCRm0CuhFxkJGXjdXoeOzJjICqepZxXwkoj8C6ioX6iq/45YKmMaKa8tZ2fpTnaU7GBn6U52le5id+ludpftpsJX8YV1MxIy6J/cn9Hpozlj4Bn06dHn80dSH9IT0hu+mDMmVoVT+DOAIoLX39dTwAq/6VCVvkq2FG9hy5EtbCvextbirWwv3k5BVUHDOm5x0z+5P4N6DmJy78kMTBnIgJQBDEgeQL/kfiTFJTl4BMZ0D+FczvmdzghiYsuhqkNsLNpI3uG8hseesj0Nryd6EhmWOoxp/aYxNHUow1KHMSR1CAOTBxLnDufGcWNMc8LppG0UwTtrs1X1eBGZAMxU1bsjns5EhQpfBRsObWDtobWsK1zHhqINHKw82PD6oJRBjMkYw4XDL2RU+ihGpo+kX3I/a5IxJkLCaep5HLgVeAyCX9qKyLOAFX7TpILKAlYcXMHKgytZXbiazUc2U6fBfmUG9xzMlOwpjOs1jnGZ4xidPppkb7LDiY2JLeEU/iRVXXrUJWr+5lY2sedQ1SGWHljKp/mfsjx/ObvLdgPB5pqJWROZPWE2E7MmMj5zPKnxNsiJMU4Lp/AfEpHhhLpWEJFLgQMRTWW6tNpALSsLVvLB3g/46MBHbDmyBYAUbwpTsqdw+ejLye2Ty+j00XhcnXJzuDHmGITzr/L7BPvMGSMi+4AdBHvoNDHkUNUh/rv3v7yz5x0+OfAJVf4q4lxx5GTncHPOzUzrO40xGWPsblNjuoFwrurZDpwtIj0Al6o2fZujiTr7y/ezeNdiluxawprCNShKnx59mDl8Jqf2P5UT+pxgl08a0w2Fc1VPL+Au4BRAReQD4NeqWhTpcKbzFVQWsGjnIt7Y8QbrDq0DYHT6aL436XucMfAMRqePti4JjOnmwmnqeY5gH/yXhOavBJ4n2OmaiQKVvkqW7F7Cy1tfZmn+UhTluIzjuDnnZs4ZfA6Detp4CMZEk7Du3FXV3zSav1tELopUINM5VJV1h9axYPMC3tz5JlX+KgYkD+C7E7/LeUPPY1jqMKcjGmMiJJzC/46IzALmh+YvBV6LXCQTSZW+Sl7Z9grzN89n85HNJHoSOW/oecwcPpOc3jnWjGNMDGi28ItIGcFLOIXgAOnPhF5yERxL966IpzMdZk/pHp7Ne5YXt75Iua+c4zKO4xcn/YLzh55Pj7geTsczxnSiZgu/qqZ0ZhATGWsL1/LUhqdYsmsJbpebcwefyxXHXcGEzAl2dm9MjArr7ppQ/zxDGq9v3TJ3XarK0vylPLb2MZblLyMlLoVrx1/LFWOuICspy+l4xhiHhXM559+BCcAGoC602Lpl7oJUlY8PfMyjax5lVcEqshKz+EnuT7h01KXWnGOMaRDOGf80VR0b8SSmXdYUruGhlQ+xLH8Z2UnZ/PzEn3PxyBJPRf8AABBSSURBVIuJd8c7Hc0Y08WEU/g/FpGxqrox4mnMMdtVuosHlj/A23veJiMhg9un3s5loy6zYQONMc0Kp/DPJVj884Eaglf5qKpOiGgy06KSmhIeXfMoz+U9h9ft5QeTfsBVY6+yLhSMMa0Kp/D/HbgKWMfnbfzGIXVax0tbX+LBFQ9SUlvCxSMu5geTf0BmYqbT0Ywx3UQ4hX+3qr4c8SSmVZuPbObuT+5mVcEqJveezM9P/DmjM0Y7HcsY082EU/jzQiNuvUKwqQewyzk7U22gljlr5/DEuidI8abwm5N/w8zhM21oQmNMm4RT+BMJFvxzGy2zyzk7yfpD67nzwzvZWryVmcNncmvuraQlpDkdyxjTjYXTH/932rJjERkI/APoQ/C7gTmq+pCIZBDs3XMIsBO4XFWPtOU9opm/zs+ctXOYs3YOmYmZ/OWsv3DagNOcjmWMiQLh3MD1JKFhFxtT1f9pZVM/8GNVXSkiKcAKEVkMfBt4S1XvFZHbgduB2445eRTbXbqbO/57B2sPreWCYRdwx4l3kOK1HjSMMR0jnKaeVxtNJwAXA/tb20hVDxAam1dVy0RkE9AfuBCYHlptLvAuVvgbvLnjTe766C7cLjf3n3Y/M4bOcDqSMSbKhNPU80LjeRGZByw5ljcRkSHAZOBTIDv0QwFVPSAivZvZZjYwG2DQoOgfCKQ2UMsflv+BeXnzmJQ1iftPv58+Pfo4HcsYE4XC6qTtKCOBsCuxiCQDLwA3q2ppuD1CquocgoO8k5ub+6WmpmiSX5HPLe/ewrpD67h67NXcPOVm4lxxTscyxkSpcNr46/vlr5dPmE0zIhJHsOj/s9HlnwdFpG/obL8vUHCMmaPKmsI13PzOzVT6Knlw+oOcPdhGtDTGRFY4TT1t+lZRgqf2TwCbVPWBRi+9DFwD3Bt6fqkt+48GL297mV9+9Euyk7J5/JzHGZE+wulIxpgY0NIIXC0256jq7lb2fTKhrh5EZHVo2c8IFvz5InItsBu4LPy40UFV+euav/LomkeZ2mcqfzz9j3ZtvjGm07R0xv8anw+9WE+BLKA34G5px6r6wVHbNnbWMWSMKr46H7/5+Dcs3LqQi0ZcxC9O+oW15xtjOlVLQy+ObzwfujLnNuBs4J6IpopSVf4qfvzuj/nvvv/y3Qnf5fuTvm/DHxpjOl04X+6OBH4OnAj8EbhRVX2RDhZtKnwV3LDkBlYXrubOaXdy+ejLnY5kjIlRLbXxH0+w4I8Dfg9cq6qBzgoWTUpqSrhhyQ1sKNrAfafeZzdlGWMc1dIZ/xpgD8G2/qnA1MbNEqp6Y2SjRYfi6mKuW3wd24q38cD0Bzhz0JlORzLGxLiWCn9rffGYVpTWljJ78Wy2F2/n4TMf5pT+pzgdyRhjWvxyd25nBok2lb5KblhyA1uKt/DQGQ9Z0TfGdBlt6bLBtKLaX80P3/4h6w+t5/7T77fulI0xXYoV/g4WqAtw+39vZ1n+Mu459R7OGXyO05GMMeYLbOy+DqSq3LfsPt7a/RY/PeGnfG3Y15yOZIwxX9Jq4ReRUSLyloisD81PEJH/jXy07ufJDU8yL28eV4+9mm+N/ZbTcYwxpknhnPE/DtwB+ABUdS0wK5KhuqNFOxfx4IoHmTFkBj/O/bHTcYwxplnhFP4kVV161DJ/JMJ0V5uKNvG/H/wvk7Im8dtTfotLrAXNGNN1hVOhDonIcEJ98ovIpYSGVDRQVFXEje/cSGp8Kg+e8SBet9fpSMYY06Jwrur5PsGRsMaIyD5gB2AN2IAv4OOWd2/hSPUR5p43l8zETKcjGWNMq8IZiGU7cLaI9ABcqloW+VjdwwMrHmBlwUruPfVexvUa53QcY4wJSzi9c8YDlwBDAE99fz2q+uuIJuvi3tr1Fs9seoYrj7uSrw77qtNxjDEmbOE09bwElAArgJrIxuke9pbt5c4P7+T4Xsdzy5RbnI5jjDHHJJzCP0BVrR/hEF/Ax63v3QrA/affb1/mGmO6nXCu6vlIRMa3vlps+NOqP7G+aD2/Ofk3DEgZ4HQcY4w5Zi0NxLIeqAut8x0R2U6wqUcAVdUJnROx61hxcAVPbXiKy0ZdxlmDY3bYYGNMN9dSU09/YFJnBenqKnwV/PyDn9M/uT8/yf2J03GMMabNWir8O1R1V6cl6eLuX3Y/ByoO8NSMp0iKS3I6jjHGtFlLhb+3iDR7yYqqPhCBPF3S+3vf54UtL3Dt8dcyufdkp+MYY0y7tFT43UAywTb9mFXpq+TuT+5mRNoIbph0g9NxjDGm3Voq/Afac5OWiPwd+BpQoKrHh5ZlAM8TvBlsJ3C5qh5p63t0hj+v/jMHKg7w9HlP26Wbxpio0NLlnO09038KOPr6/9uBt1R1JPBWaL7L2lC0gX9u+iffGP0NJvW277mNMdGhpcLfrusVVfV94PBRiy8E6gdxnwtc1J73iCR/nZ9fffQreiX04qacm5yOY4wxHabZwq+qRxftjpCtqgdC+z8A9G5uRRGZLSLLRWR5YWFhBKK0bF7ePDYd3sQdJ95Bijel09/fGGMipcuOGKKqc1Q1V1Vzs7KyOvW9i6qKeGT1I5zc/2TOHnR2p763McZEWmcX/oMi0hcg9FzQye8flj+v/jNV/ip+esJPqe+N1BhjokVnF/6XgWtC09cQ7PmzS8k7nMcLm19g1phZDEsd5nQcY4zpcBEr/CIyD/gYGC0ie0XkWuBe4BwR2QKcE5rvMlSV+5beR1p8GtdPvN7pOMYYExHhdMvcJqr6zWZe6rK9my3ZvYTlB5dz57Q7SY1PdTqOMcZERJf9crez+ev8PLTyIUakjeCSkZc4HccYYyLGCn/IS1tfYlfpLm6cfCNul9vpOMYYEzFW+IGaQA2PrHmECVkTmD5wutNxjDEmoqzwA8/lPcfByoPcNPkmu3zTGBP1Yr7wl9eW87d1f+Okvicxte9Up+MYY0zExXzhf3rT0xTXFFt/PMaYmBHThb/CV8EzG5/hjIFnMC5znNNxjDGmU8R04Z//2XxKa0u5bvx1TkcxxphOE7OFv9pfzdwNc5nWdxrjs8Y7HccYYzpNzBb+F7e+SFF1EbMnzHY6ijHGdKqYLPy+Oh9Prn+SSVmTyM3OdTqOMcZ0qpgs/K9vf539Ffu5bsJ1dt2+MSbmxFzhV1We2vAUo9JHcWr/U52OY4wxnS7mCv/S/KVsLd7Kt477lp3tG2NiUswV/mc2PkNGQgbnDzvf6SjGGOOImCr8u0t3897e97hs1GXEu+OdjmOMMY6IqcL/bN6zuF1uvjH6G05HMcYYx8RM4S+vLefFrS8yY8gMspKynI5jjDGOiZnC/+LWF6nwVfCt477ldBRjjHFUTBR+VeX5z55nQtYE64zNGBPzYqLwryxYyc7SnVw26jKnoxhjjONiovC/sPkFkuOSOXfwuU5HMcYYx0V94S+pKeE/u/7D+UPPJykuyek4xhjjuKgv/K9tf42aQA2XjLrE6SjGGNMlOFL4RWSGiHwmIltF5PZIvY+q8sKWFzgu4zjG9hobqbcxxphupdMLv4i4gb8A5wFjgW+KSESq8oaiDWw+splLR10aid0bY0y35MQZ/1Rgq6puV9Va4Dngwki80YLNC0j0JHL+UOuXxxhj6jlR+PsDexrN7w0t+wIRmS0iy0VkeWFhYZveaGDKQK4YcwXJ3uS2JTXGmCjkceA9m+oLWb+0QHUOMAcgNzf3S6+H49rx17ZlM2OMiWpOnPHvBQY2mh8A7HcghzHGxCQnCv8yYKSIDBURLzALeNmBHMYYE5M6valHVf0i8gNgEeAG/q6qGzo7hzHGxCon2vhR1deB1514b2OMiXVRf+euMcaYL7LCb4wxMcYKvzHGxBgr/MYYE2NEtU33RnUqESkEdrVx80zgUAfGcVK0HEu0HAfYsXRV0XIs7T2Owar6pUHGu0Xhbw8RWa6quU7n6AjRcizRchxgx9JVRcuxROo4rKnHGGNijBV+Y4yJMbFQ+Oc4HaADRcuxRMtxgB1LVxUtxxKR44j6Nn5jjDFfFAtn/MYYYxqxwm+MMTEmagp/awO4i0i8iDwfev1TERnS+SlbF8ZxfFtECkVkdejx/5zIGQ4R+buIFIjI+mZeFxF5OHSsa0Ukp7MzhiOM45guIiWNPpNfdHbGcInIQBF5R0Q2icgGEbmpiXW6/OcS5nF0i89FRBJEZKmIrAkdy6+aWKdj65eqdvsHwe6dtwHDAC+wBhh71Do3AI+GpmcBzzudu43H8W3gz05nDfN4TgNygPXNvH4+8AbBUdmmAZ86nbmNxzEdeNXpnGEeS18gJzSdAmxu4u9Yl/9cwjyObvG5hP6ck0PTccCnwLSj1unQ+hUtZ/zhDOB+ITA3NL0AOEtEmhoG0kmdNhB9Z1DV94HDLaxyIfAPDfoESBORvp2TLnxhHEe3oaoHVHVlaLoM2MSXx7zu8p9LmMfRLYT+nMtDs3Ghx9FX3XRo/YqWwh/OAO4N66iqHygBenVKuvCFNRA9cEnoV/AFIjKwide7i3CPtzs4KfSr+hsiMs7pMOEINRdMJniG2Vi3+lxaOA7oJp+LiLhFZDVQACxW1WY/k46oX9FS+MMZwD2sQd4dFk7GV4AhqjoBWMLnZwHdUXf4TMKxkmCfKBOBPwEvOpynVSKSDLwA3KyqpUe/3MQmXfJzaeU4us3noqoBVZ1EcAzyqSJy/FGrdOhnEi2FP5wB3BvWEREPkErX+/W91eNQ1SJVrQnNPg5M6aRskRDO59blqWpp/a/qGhxdLk5EMh2O1SwRiSNYLP+pqv9uYpVu8bm0dhzd7XMBUNVi4F1gxlEvdWj9ipbCH84A7i8D14SmLwXe1tA3JV1Iq8dxVFvrTIJtm93Vy8DVoatIpgElqnrA6VDHSkT61Le3ishUgv+uipxN1bRQzieATar6QDOrdfnPJZzj6C6fi4hkiUhaaDoROBvIO2q1Dq1fjoy529G0mQHcReTXwHJVfZngX5KnRWQrwZ+Us5xL3LQwj+NGEZkJ+Akex7cdC9wKEZlH8MqKTBHZC9xF8IsrVPVRguMunw9sBSqB7ziTtGVhHMelwPdExA9UAbO64ElFvZOBq4B1oTZlgJ8Bg6BbfS7hHEd3+Vz6AnNFxE3wh9N8VX01kvXLumwwxpgYEy1NPcYYY8Jkhd8YY2KMFX5jjIkxVviNMSbGWOE3xpgYExWXcxrTUUSkF/BWaLYPEAAKQ/OVqvoVR4IZ04Hsck5jmiEivwTKVfUPTmcxpiNZU48xYRKR8tDzdBF5T0Tmi8hmEblXRK4M9am+TkSGh9bLEpEXRGRZ6HGys0dgTJAVfmPaZiJwEzCe4B2ko1R1KvA34IehdR4CHlTVE4BLQq8Z4zhr4zembZbV918jItuA/4SWrwPOCE2fDYxt1G16TxFJCfUfb4xjrPAb0zY1jabrGs3X8fm/KxdwkqpWdWYwY1pjTT3GRM5/gB/Uz4jIJAezGNPACr8xkXMjkBsaLW0jcL3TgYwBu5zTGGNijp3xG2NMjLHCb4wxMcYKvzHGxBgr/MYYE2Os8BtjTIyxwm+MMTHGCr8xxsSY/x92PU9ab1Gr2wAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret1 = run_simulation(duration, y0=y0, model=m)\n", "ret1.plot(**opt_kwargs)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Simulating with `gillespie` (Bars represent standard error of the mean):" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEGCAYAAABiq/5QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXhU5dn48e89kz2EJRCSsMmOshMChbpvlEqLWpfSWvX17VtotT9p7UtrXVNr1dZWi31b9yrVFqRYBcWq4FpEhCSyC8i+mJCEneyZuX9/zAQTMkkmy+RkMvfnuuaaOWfOc859OOQ+Z57znOcRVcUYY0zkcDkdgDHGmLZlid8YYyKMJX5jjIkwlviNMSbCWOI3xpgIE+V0AMHo0aOH9u/f3+kwjDEmrOTk5BSpasrp88Mi8ffv35/s7GynwzDGmLAiInsCzbeqHmOMiTCW+I0xJsJY4jfGmAgTFnX8xhjjhIqKCnbs2EFJSYnToTQoISGBQYMGERMTE9TylviNMaYeO3bsoGvXrgwbNgyXq31WkHi9Xg4ePMj27dsZPnx4UGVCuici0lVEFonIFhH5TEQmi0iyiCwTkc/9791Cse01S54kP2sw3nu7kJ81mDVLngzFZowxHVhJSQmpqantNukDuFwuUlNTKSkpIT8/P7gyIY5pLvCmqp4JjAE+A24H3lHVIcA7/ulWtWbJk4zMuYs0CnEJpFHIyJy7LPkbY5qsPSf9ai6XCxFh8eLFeDyexpcPVSAi0hk4D3gWQFUrVPUocDkwz7/YPOCK1t5239yHiZeKWvPipYK+uQ+39qaMMaaWbz/5Md9+8mNHtl1WVkZpaWmjy4XyVDYQKASeE5FPReQZEUkEUlU1D8D/3jNQYRGZKSLZIpJdWFjYpA331MDL99SiJq3HGGPag1deeQURYcuWLQ0uJyJBrS+UiT8KyAAeV9VxQDFNqNZR1adUNVNVM1NS6jxx3KACCbx8gfRo0nqMMaYpXv30AJ/uPconuw5z9kPv8uqnB1plvfPnz+ecc85hwYIFrbK+UCb+/cB+Vf3EP70I34ngoIikA/jfC1p7w/sy5lCqtZs1lWoM+zLmtPamjDEG8CX9X/5rAxUeLwAHjpbyy39taHHyP3nyJB999BHPPvts+0/8qpoP7BORYf5ZFwObgSXAjf55NwKLW3vbE6bPYuP4+8knBa9Cvrcbnwz5KROmz2rtTRljDAAPv7WV0sraN1ZLKz08/NbWFq331VdfZerUqQwdOpTk5GRyc3NbtD4Ifaue/wf8XUTWA2OBB4CHgEtF5HPgUv90q5swfRZpWdspmrmWFDlKbElwzZyMMaY5vjga+KZqffODNX/+fGbMmAHAjBkzmD9/fovWByF+gEtV1wKZAb66OJTbraln7wGsTZzEkC8WU1nxMNExsW21aWNMBOnVNZ4DAZJ8r67xzV7noUOHePfdd9m4cSMigsfjQUT43e9+F/SN3EDafwPV1jDuBnpwlI3v/9PpSIwxHdScrw0jPtpda158tJs5XxtWT4nGLVq0iBtuuIE9e/awe/du9u3bx4ABA1ixYkWLYo2IxD/ygqspIBlX7vNOh2KM6aCuGNebB781ihi3L6327hrPg98axRXjejd7nfPnz+fKK6+sNe+qq67iH//4R4tijYi+eqKiY9jR91t8Ze+z7Nu+gb6DRzkdkjGmA7piXG/mr94LwEuzJrd4fe+//36debfeemuL1xsRV/wAQ6bNpgoXX7z1R6dDMcZ0YC/NmtwqST+UIibx90jrx7quFzOy4HWOHz3kdDjGGOOYiEn8AF0vvJVEKWPzG39xOhRjjHFMRCX+IWPP5bPo4fT9/AU8VVVOh2OMMY6IqMQPUJoxk956kPXvts6jz8YYE24iLvGPvuQ68kkhJtv65jfGhMBz03yvdiziEn9UdAy7h1zPiIr1bM1+1+lwjDGmQW63m7FjxzJmzBgyMjJYuXJli9cZcYkfYNT02RylE6Xv/s7pUIwxHcn6hbB/DexZAY+O9E23UHx8PGvXrmXdunU8+OCD/PKXv2zxOiMy8ScmdeWzftcxtuRjdm78pPECxhjTmPUL4bVbwVPumz62zzfdCsm/2vHjx+nWreXDlEfEk7uBDL9iDsVz53HkrYdg5CtOh2OMCXfv3AeVp3XSVlnqmz/62mavtrS0lLFjx1JWVkZeXh7vvtvyKuqIvOIH6JKcwvpeVzP2+Hvs277B6XCMMeHu2P6mzQ9SdVXPli1bePPNN7nhhhtQ1RatM2ITP8CQy2+nkijyXn/A6VCMMeGuS5+mzW+GyZMnU1RURFPHIT9dRCf+Hml9WZcynXFH3uLAzk1Oh2OMCWcX3wPRp/W9Hx3vm99KtmzZgsfjoXv37i1aT0QnfoBB37qXKtzkvdp6B8cYE4FGXwvffAzc/sGeuvT1Tbegfh++rOMfO3Ys3/72t5k3bx5ut7vxgg2I2Ju71Xr0OoOPe83gK1+8wK5NnzBgxFecDskYE65GXws583yfb1raKqv0eDyNL9REEX/FDzD8mns4KfEcW3qv06EYY8LdTUtbLemHiiV+fC18Nve/ibElH7NlzXKnwzHGmJCyxO83+upfUERXPMuyUK+XTQ+cw6YHznE6LGOMaXWW+P0SOnVhx1k3M6JiA2uXt2w8S2OMac8s8dcw/ls/ZberHz0/vo8qrzgdjjHGhIQl/hqiomM4ccH99NaDFFTGOh2OMSYM3fTmTdz05k1Oh9GgkCZ+EdktIhtEZK2IZPvnJYvIMhH53P/e8h6HWtGo8y7n04Sv8lU2cMIT7XQ4xpgIl5+fz4wZMxg0aBDDhw/nsssuY9u2bS1aZ1tc8V+oqmNVNdM/fTvwjqoOAd7xT7crKd96mGgq8VRWOh2KMSaMLN25lPWF68k+mM2URVNYurNlzTpVlSuvvJILLriAHTt2sHnzZh544AEOHjzYovU6UdVzOeB/woF5wBUOxNCgvM0fUaRdmCwbyMsaxJolNlqXMaZhS3cuJWtlFhXeCgDyivPIWpnVouT/3nvvER0dzQ9/+MNT88aOHcu5557bolhDnfgVeFtEckRkpn9eqqrmAfjfewYqKCIzRSRbRLJb2iFRU6xZ8iQjc+6il+swLoF0ihiZc5clf2NMg+bmzqXMU1ZrXpmnjLm5c5u9zo0bNzJ+/PiWhlZHqBP/2aqaAXwduEVEzgu2oKo+paqZqpqZkpISughP0zf3YeKlota8eKmgb+7DbRaDMSb85BfnN2m+k0Ka+FX1C/97AfAKMBE4KCLpAP73glDG0FQ9NfCvi55a1MaRGGPCSVpiWpPmB2PEiBHk5OQ0u3x9Qpb4RSRRRJKqPwNTgI3AEuBG/2I3AotDFUNzFEjgXxcH6YY3BJ0lGWM6htkZs4lzx9WaF+eOY3bG7Gav86KLLqK8vJynn3761Lw1a9bwwQcfNHudENor/lRghYisA1YDS1X1TeAh4FIR+Ry41D/dbuzLmEOpxtSaV6Ix7Pd2Z80iq+4xxgQ2beA0sr6aRYzLlz/SE9PJ+moW0wZOa/Y6RYRXXnmFZcuWMWjQIEaMGEFWVha9evVqUayNdsssItcAb6rqCRG5C8gA7lfV3IbKqepOYEyA+YeAi5sZb8hNmD6LNUCfnIdI5TAFksK+jJ8R89nLDN38B/Z9/nX6DqmzW8YYw7SB01i0bREAz019rlXW2atXLxYubL0B2yG4K/67/Un/HOBr+JpgPt6qUbQzE6bP4nBsbz6LHUVa1nYmXP4jet/wLOUSQ+lL/0NVZUXjKzHGRKTnpj7Xakk/VIJJ/NUV29OAx1V1MRDTwPIdUo9eZ7Bj4n0MrdrGmhfvdjocY4xptmAS/wEReRK4FnhDRGKDLNfhjL/s+2QnXUzm7qfZkv2O0+EYY9qA1+t1OoRGNTXGYBL4tcBbwFRVPQokA3OaHlrHMOS/n6LQ1Z2ur8/i2KGWPTZtjGnfEhISyM/Pb9fJ3+v1kp+fT2UTuphp9OauqpaISAFwDvA5UOV/j0hduvXg4Defpv/ib7H52RsZ879vIK6I/AFkTIc3aNAg1q5dyxdffIFI++2qvbKykr1796KquILIR8G06rkXyASGAc8B0cCLwNktjLVdG3HHinq/G5pxAau23MakbQ+zav6vmXSdjdVrTEcUExNDfHw8y5Yto0uXLkElVaecOHGC1NRU4uPjG102mL24EpgOFMOpp3GTWhRhB/CVGXfwacLZjN82l80f/7vWdzZsozEdx8iRI7nggguIi4tDRNrly+VyMWTIEKZPnx7UL5NGr/iBClVVEVE49RRuxBOXi4E/+Bv5j51L2lszyeu1nPQzhjkdljGmlYkIGRkZZGRkOB1Kqwnmin+hv1VPVxH5AbAceLqRMhGhS7ce6HfmE62VlPxtBiUnjzkdkjHGNKrRxK+qvwcWAS/jq+e/R1X/FOrAwkW/oWPZef5jDKjaxZYnrkfb8d1/Y4yB4Kp6UNVlwLIQxxK2xlx0Lav2b2DSzsf4+Nmf0tnpgIwxpgH1XvGLyAkROV7j/XjN6bYMMhx85Xu/YnXyN5l84HkOVCQ4HY4xxtSr3sSvqkmq2rnGe+ea020ZZDgQl4uMH/2VtfGTuFhXc6DC7oEbY9qnRuv4RaRfoFdbBBduoqJjGHrLQj7XPpyruWxe9abTIRljTB3BtOpZWuP1DrAT+HeDJSLYpncX0EWKiaWCbm/8kGXPWoduxpj2JZguG0bVnBaRDGBWyCIKY9UDtVeP2ZsuR+i69wmWPxfFJTfZ073GmPahyc8f+wdgmRCCWMJefQO1j9j1HLs2r3EoKmOMqS2YvnpuqzHpwjcCV+ARySNcTy2EAE9Lp8oRji28ku1XLmDwGOvKwRjjrGCu+JNqvGLx1fVfHsqgwlW9A7VLD8qIo+cr17BlzfI2jsoYY2oLpo7/V20RSEewL2MOXWrU8QOUagz7x/+cvmMv4vhz36Tf699lQ+nTjDrPzp3GGGcE05xzmYh0rTHdTUTeCm1Y4WnC9FlsHH8/eZqMVyGfFDaOv58J02eR1m8IcT94i4PudIa9cxPZS55wOlxjTIQKpqonxT/yFgCqegToGbqQwludgdqnf9kAqkevM+h+67t8HjuCzNxfsOpvd1vfPsaYNhfUYOs1H9gSkTMADV1IHVvnrt0ZfNtb5CRdyKSdj7Hm/26gorys1jLWn78xJpSCSfx3AitE5AUReQH4EPhlaMPq2GLjEhj3k5f5uNeNTDz8Gtt/fzGHCw44HZYxJkIE0y3zm/iacL4ELATGq2rQdfwi4haRT0Xkdf/0ABH5REQ+F5GXRCSmucGHM5fbzeSZj5Gd+TADK7ZS/pfz2b7uI6fDMsZEgGBu7gowFchQ1deABBGZ2IRtzAY+qzH9W+BRVR0CHAG+34R1dTiZ35jJviv+hRsPff91OasXPYJ6rSbNGBM6wVT1/AWYDHzHP30C+HMwKxeRPsA04Bn/tAAX4RvYBWAecEUT4u2Qhow7j6ibV7AtbhQTN/6K4xVQ4W2/gzobY8JbMNnlK6p6C1AGp1r1BFs980fg50B105XuwFFVrfJP7wd6ByooIjNFJFtEsgsLw+tB4RF3rGDEHSuaVCa5Z2+Gz1nGa/FXMIADjKlcR17WINYseTJEURpjIlUwib9SRNz4W/KISApfJvJ6icg3gAJVzak5O8CiAes1VPUpVc1U1cyUlMBPxHY0uW88yyUlb5DuOoJLIJ0iRuTcxepXH3c6NGNMBxJM4n8MeAXoKSK/AVYADwRR7mxguojsBhbgq+L5I75B26ufGO4DfNHUoDuqQJ28JUgFfT/9HXu25DoUlTGmowmmVc/f8VXXPAjkAVeo6j+DKPdLVe2jqv2BGcC7qnod8B5wtX+xG4HFzYy9w+mpgau0UjlM+vxL+XjeHVRWlLdxVMaYjqahMXeTq19AATAf+Adw0D+vuX4B3CYi2/HV+T/bgnV1KA118rYx6Wwm7/oz+x6ayJbVNu69Mab5GrrizwGy/e+nv7KbshFVfV9Vv+H/vFNVJ6rqYFW9RlXtEtZvX8YcSrX2ffNSjWF/xs/J+N8l5E7+PxK9JzjzjatZPfe7HCnMcyhSY0w4E9X232Y8MzNTs7ObdK4JW2uWPEmfnIdI5TAFksK+jDm1+vspPnGUDX+/g/F5CyiVWDYP+SEZV/+CmNg4B6M2xrRHIpKjqpl15geT+EVkOnCef/J9VX29leNrUCQlfuBUPz0NNQnd81kORxf/gjFla9gv6RR85XbGTbmBzQ+d12hZY0xkqC/xB/Pk7kP4nr7d7H/NFpEHWz9E0xRnnDWeMbcvZ/0Ff6VKoshYNZvtD0wkrzLRnvw1xjQomOaclwGXqupfVfWv+LpvmBbasEywRl9wFX3vWMvqMfeTVHWUS3QVWlHMxhVLmtzls/UKakxkCLZfgK41PncJRSCm+dxRUUy88v/R7fb1vM1XSJUjjFx+PVsfPJt17/3T+vw3xtTS6NCL+Nrvfyoi7+F78vY8rFvmkGpu/fz6t19glH5OD46Sp905WObm/A/+h13/uZ+i0TMZ8/Xvh+wmcDD3JYwx7UMwD3DNByYB//K/JqvqglAHZppmzZInGZlzF+ly2NfdgxxiomsrS7teDygT1t7BsQfP5OPnfkFR/j6nwzXGOKihB7gyql9AOr4O1fYBvfzzTDsSqLuHeKlg/NE36X/XWtad/wx5cYOZvOcJOj8+huxHrmLTyjesGsiYCNRQVU82sAmo7kegZgdriq/vHdNO9NTCgF3g9dQixOVizIXXwIXXsHfbWvKW/YmzCt+g89vL2bu8N18MuJrBF0f0sAjGRJSGEv/PgKuAUnydrL2iqifbJCrTZAWSQhp1+/opkB6k1ZjuN3Qs/YY+S2nxCdYsm0fSxheZtGMuq7e+QV/XIVLlMHlZg9if8fNaD44ZYzqOeqt6VPVRVT0H+DHQF3hHRBaKyNg2i84Erb7uHvZlzAm4fHxiEhOu+DFn3rWKdwfOYaR7D+muw7W6g/73n39GWWlxW4RvTUmNaUPB3Nzdha8HzbeBicDQUAdlmm7C9FlsHH8/eZqMVyGfFDaOvz+oq/bhO58nIUB30OMOvoznoYHk/OEKct54lhPHDocqfMfYCcdEonqrekRkIL7ulC/Hd1N3AfAbVS1ro9hME02YPotNG1/gML0ZcceKWlU8Dan3/oAcIbv7Nxl8+AOSV79HxSdz2BA3huL+l9B7wnT6Dh4F+PsXKj9AKofJzxpcp38hY0z70lAd/3ZgPb6r/eNAP+Bm37C5oKqPhDw60ybqvz+QwsRbX6SqsoLNOe9yfO0SehW8z6itv4Otv2O/pLE2agwXV7x3qkVRGoV0ybmLNdBmyd+eITCmaRpK/Pfx5bCIndogFuOQfRlz6JJzV63moKUaw77xc0gDoqJjGD5pKkyaCsD+7Rs5kPM6cbvfZXzxauJddZuR9sn9HXTwq3474ZhwVW/iV9WsNozDOGjC9FmsgdrdQY+vv7qmz+CR9Bk8Ergd772Be/BI1SJ2/HocRd3G4R4wmV4jziO93xDEVfe2klUVNZ1TJx072XUMwXTZYMJIc/8gm3t/oL5qooN0pzSqM6MKXyeh6GVYA4fowr6E4ZSmjCGhfyZ9R3yVHateZ2SNXxttWVUUqSccp5L3TW/eBMBzU59r0+06ue1rn/I1glw4c22bbrcxlvhNi9RXTbR//C+YMH0WVZUVbN+8mkNbV+I6kEPqiQ302b0K154n4AMo93YLWFXUN/fhkFYVVXdx4dS9CbtybjonTxxOCOVJo6FWPbNVda6InK2qH7X6lk2H0Fg1UVR0DIPHnMPgMV82mTxx7DB7N63ixM7VTPz80YDr7amFZD9yNZXdhxKXPpzu/UeS3v8somNiay3X3Kv2+rq4CPUJp9qveviehVwY8i21H8V7P3U6BOPX0BX/TcBc4E+A9c1j6tXUaqKkLsmM+Opl8NXLyM96sZ6qomT6Hs8l9fgy2AWshEp1s8+VyuG4vpQm9aeoDC4+9kqzrtob6uIiHLTkxBGJJ52WaK/VNS3RUOL/TER2Aykisr7GfAFUVUeHNDITEeqvKrqdCdNnceLYYfJ2bODYvk1UFWwl9thOupbsZWjJpxzTxMAtinIeYlXBVlzd+hGfcgadUweQ0nsACZ2+vBEdbBcX9YnU+wPhyn5t1NZQq57viEga8BYwve1CMpGksaqipC7JJGWcDxnn1yqnXi+xv+oWcJ2pHKbnvr/i3l97CMrjJHLY1Z3jMSnsc4/hwqoPaz2xXKox7Bpza6OJ3+n7A05Zs+RJepXvJ50jdrILcw3e3FXVfGCMiMTwZVcNW1W1MuSRmTbn1I3G5rQoEpeLgw08eNb9jk0cPLCLo/m7KCncQ+WRfbhOfEFMyUE6VRQwwbOLjZ4z6OcqoqccoUC7sU97MHn93ZSsu5+jrq6ccHelNLobFbHJeOKTkYTuuDr1oO+6Pzp6f2DNkid59MDRNk3ATp/slu5cytr8bKqAKYumMDtjNtMGBjcC7NKdS/k8xtOssh1Vo616ROR84G/AbnzVPH1F5EZV/TDEsRnToAYfPIuJpdeAM+k14Mx6y3evqmLNg5eyRYeQfO4s9OgXfHyiECkuJKq0kNjywyRVFJBUuo2uR44RIx4AvEo99wcK2XnfGMqikiiPSqIquhPe6CS8cV2Q2CQkrgvuhM5Ex3eh34kulLu8HNj5GXGJSSQmdSU2LiHgcw41tTQBN/ek0dKb4S1Jvkt3LiVrZRZV/n/zvOI8slZmATS6jpaU7ciCac75CDBFVbcCiMhQYD4wvqFCIhIHfAjE+rezSFXvFZEB+Pr9SQZygetVtaL+NRkTWFMfPDudOyqKJHclSVQy4sKrG1xWvV5OnjzG8cMFyPOXkU7dm8AHSeZYXC9iqk7QpewL4kuKSaSETlqCS76sdnosfgCfdq+kMErY/O7VXH4omltLd+FRoYQ4SiWecomj3BVPpSuOSnc8HnccHnc8fY7lBEzAfXJ+S7bbhTs6HndMPFGxCUTFxhMVE090bBzRsQl8/vFiMjf/ts1vhrc0+c7NnUuZp3YXYWWeMubmzm20fEvKAjz2z9s45PIdq4ufGcHlXb7GrdcE11uNU2WDIara8AIi60+/kRtoXoByAiSq6kkRiQZWALOB24B/qeoCEXkCWKeqjze0rszMTM3Ozg5id0y4akm79rYue/pVN/h+adTXG6rX46Gk+DglJ47y/PL7+ad3FWU1ruzjvF6+eWIQl8YPRCqLcVUW464qxV1VQrSnhGhvGdHecmK1jDTvQVwBErBXCTi/pnxvN3KSypnbrSv5UW7SqjzMPnKUjBOxRLuUKqKplGg84sYj0XgkyvfuiqZ72V7SpW7vrHmazN6UC1BXFLiiUZcbcUWjrijEFQUuN7/QtzlMSZ2y3V2d+H3yf+NyRYErCpc7CnG5/a8oxO3C5YpixrafEyhLCfDapBdwud24XC7EFeV/d+F2RYHLxUVvTUMDlBaENd9ehcvlRkROvdf8xfXYP2/jhZNv1TlW13dqPAk7VbbOforkqGrm6fODueLPFpFngRf809cBOY0VUt8ZpXrglmj/q3rkru/6588DsoAGE7/p+Fpyf6GtmydW/9LolfMg6Rxp9JeGy+2mU+dudOrcjbcqV1EWXbs6p8zl4oOEHdzzP0sa3XZ+1mByEovrJu/iBCpm/JPK8hKqykupqijFU1GCt7IcT0UZ3opSDu74Lff1SD6VUPKio8jqkcw9epiecZfg8lQg3kpc3kpc3gpcWoVLq3B7K9nIILrqyTonu+3edEYWLSNKq4jCgxvPqSqxaof79wv4a+GQ5yTj1/wvHsAj4EHwAF6pfvfNS+mVSkFU3VTVvaqKkpcvwQt4EbziSzCeGp+7pfTgcJS7TtmunipWPDbQXxZUBPV/9sUhLO7WNeCxWnzsLbr8YYB/ecF7qqzgxfd5ceeYestW/PEsFPwv3/LV2wfhnUStt+ytdf8ZmyWYxP8j4BbgVnyH70PgL8GsXETc+E4Sg4E/AzuAo6pa5V9kP9C7nrIzgZkA/fr1C2ZzxjTZ7zOGANDUZ0EnTJ/FtfmPA11ZOHNt0F1cFEYFviyvb/7pFo6YXutqsDp5X9//a/x42FjKPeVUeCoo95TX+lzhqeCPR56udRUJvoTyx+Rkbpx8KRXeCiq9lVR6KqnyVvk+e7/8/NpncWySXRRGCSlVylBPH8rO6IVHB1PlrfK9tIoqj++zR6uo8nroeTyfgui69y56VnkZPSC4v22X+k4ENaeLoqK4tnd6g+VGHutESdLxOlfPvU924SepdU8INUk9tSGFUcLvezTcb2VDZed1azjtSj2VMMH+HwlGo4lfVcvx1fM3uYJJVT3AWBHpCrwCnBVosXrKPgU8Bb6qnqZu20SOxH7j2nybzb1ZmVKlFETX/QNOqVJe2/EapVWllFSW+N6rfO81X6vLVlMRIHk/XbKMp/82puGNRwW+cVwQ5eLh7IdPTbvERbQrmihXFDGuGKJd0VR4KjgWcwz1j91UEC0URufRrySatMQ0EqIScLvcREkUUa6oU5/dLjcTdz/DfSnJdZLvTw4fZv+Ue3CJC7fLjVvcvs81310u9JM1ROXN54nkTqd+5fzw8Emk9w10PvciXPjKu8SFIKfKigh9nphBbnlJgF9IZRz/3hun9re6rIic+nzTwq/Xe6xeuu59gFPbqy4rIgjCFX+bVG/ZN2/6FMRX9lQ5/zvAxc+MqLdsa2mTvnpU9aiIvA9MArqKSJT/qr8P8EVbxGBCq6X9qIRLPyyqyiufv8IDqx+odbPy7o/u5v2979MrqRcnKk58+ao8wcmKk75X5UkGlnTmeIAr0J4lnbljxR21thXnjiM+Kp74qHjiouKIi4qjwlt/O4gfjfkRMe4Y4txxxEbFEuOKITYqllhXLDHuGO5ccSeHyg7VKZeamMorl79CtCuaaFc0blfdK+Epi6ZwtOJo7X8LlApvBc987ZkG/83yl80nq+hwneQ7viSRb479UYNlAfL/8UvSOM63So7Xnv/5S6Rd95sGy3q1iGnFMK245LT5paQnD2uw7OVdvhawrv3yLl+jR3yPZiAIFtQAABeGSURBVJeNdkc3u2xrCVniF5EUoNKf9OOBS4DfAu8BV+Nr2XMjvoFejGmWpTuXsr5wPRXeima30faqlwMnD1BUWsSh0kMcLjt86v1I2REOl/vej5Yd5Wj50YDJt9JbyZt73iTKFUXnmM50julMp+hOdIrpRGpC6qnPidGJ9F/7IasrN56qMrkwfhLfuu5nJEYnkhCVQEJ0AnHuuHoTcF5xXp356Ynp3Dz25gb3c86EOWStzKrVyiXOHcdPx/+UpJikBsvmF+c3aX5N+zLmcFHOXUwr/vIaz3czfE5QVWQtaVHUkie0b73mEfjnbSw+9tapYxVs6xqnygarwcTvr6N/SFUDj9jdsHRgnn8dLmChqr4uIpuBBSJyP/Ap8Gwz1m3MqWaC1Yn49GaCqsrhssPkl+STX5zPweKDFJQU+F6lBRSWFLLn+B486mHqy1PrrD8pOolucd3oFteNXp16MaL7CLrGdeW5jYF/lQhC7vdyT/1kr9fYW7j2qbF0r2h6/y+zM2YHTN6zM2Y3Wrb6hHjXh7dTBaR3Sg/6RJmWmBbwhJOW2HjqburN8NO1JHk3NshQY2695hFWNPNYOVU2GI09uesRkfEiItpYu8+6ZdcDdSpfVXUnvkHbjWmR+tpo3/PRPTy+7nHyTubVuTqPckWREp9CSkIKCVEJp5r6dYntwpWDrmTqwKl0j+tOclwyMe6YgNt9c9eb9SbBRpN+C7UkeVeXn7f8TgAWXv120NttyQkHmn8zHFqWvFt60umogqnq+RRYLCL/BIqrZ6rqv0IWlXGEU/XsjVXXnKw4ye7ju9l1bBe7j+9mz/E97D2+N2DyBajwVjCs2zAu7HshaYlpX74S0ugW1w2XuE79WvCqF4Bj5cdYsHUBZ3Y/kxHdRzQYb0uTYEs1N3m3dJsA93x0DxXeCtITm3bCaYmWPqjXkpNORxVM4k8GDuFrf19NAUv8psUCVdfc/dHdvLbjNTzqYefRnRSUFpxa3i1uenfqTb/O/UiISqCkqu6DQemJ6fzhgj80uN2WPNHpZBJ00rSB01i0bRHQ9hcHzR0hrtq9RTZseE3BNOe8qS0CMZGlqLSIzYc285tVv6mTgCu9laz8YiXDuw9nUq9JDOgygIFdBtK/S3/6dup7qlVE9UmjOVfeLblZCc4mQWhZE1Ynmr+a9iWYTtqG4nuyNlVVR4rIaGC6qt4f8uhMWGisqqa4sphNRZtYX7SeDYUb2HRoEwdLDja63gXfWNDg9y258m7JzcrWEInJN1yvujvSACzVgqnqeRqYAzwJvpu2IvIPwBK/CVhVc+/Ke1lb4PtjWVu4lm1Htp2qSz+j8xmMTx3PiO4jGNFjBLd/eDv5JXWvsoNNwM298na6nr69P6/QHkXa+MShPOEEk/gTVHX1aa0Vqupb2DirrW/QBqorL/eUs2DrAuKj4hmTMoaZo2cyJmUMo3qMoktsl1rL/mT8TxxJwJFaTw+RedJx6qTRXn8tBJP4i0RkEP6uFUTkaiBwcwoTtpryIFSFp4LcglxW7F9Rb8saQVj5nZVEuRr+L+ZkAna6nr4lwi1eiLwr9vYsmMR/C74+c84UkQP4hr6+LqRRmTbV2INQ4LsZ+5/9/+G9fe+xKm8VpVWlRLuiiXHFBHySNS0xrdGkXy1cE3A4xdpaInGfO6JgWvXsBC4RkUTApaonQh+WaUv1NW18JPsRikqLWL5nOesK16EoaYlpTB80nXN7n8uEtAm8t+89R+vKwZKRMU0VTKue7sC9wDmAisgK4D5VrdvbkwlL9TVhLCgt4PfZv2dYt2H8aOyPuLDvhQzrNqzW06mRXFduTLgK5rf4Anx98F/ln74OeAlfp2umHWluh2WpCakBW9YkRSex4BsL6Ne54T7Tw7WqBsIvXmNaQ1BP7qrqr2tM3y8iV4QqINM8wdTT16SqbCjawKJtizhUWvfHW5w7jjsn3dlo0m8tloCNaTvBJP73RGQGX45sdzWwNHQhmeYItguCksoSXtvxGgu3LWTbkW3ER8XzzcHfpEd8D57f+LxV1xgTAepN/CJyAl8TTsE3QPqL/q9c+MbSvTfk0ZmgNdYFwb7j+/jHln/w6vZXOVl5krOSz+Keyfdw2YDLSIxOBCD3YC5gV9/GdHT1Jn5VbXhkBtOu1NcFQfe47tz2/m0s37Mct8vNlDOm8N2zvsvoHqND3oWwMaZ9Cqqhtb9/nv41l7dumduXQF0QuHBRVFbEqi9W8f1R3+e7Z36XlIQUB6M0xrQHwTTn/CswGtgEeP2zrVvmEGpOtwvVI07d/dHdVKmvR43E6ERmjZnF1UOvPlWdEypWPWRM+Ajmin+Sqg4PeSSmRdYVruNf2/9FlVYR7Yrm5xN+zpVDriTWHet0aMaYdiaYxP+xiAxX1c0hj8Y02Z7je3gk+xHe3fcuyXHJ9E3qS0p8CjPOnOF0aMaYdiqYxD8PX/LPB8rxtfJRVR0d0shMg46VH+OJdU+wYMsCYtwx/Hjsj7l++PXc8s4tzV6nVdcYExmCSfx/Ba4HNvBlHb9xiFe9LN6+mEdzHuVYxTGuHHwlPx73Y3rE9wAseRtjGhdM4t+rqktCHokBGu52YduRbdy/6n4+LfiUcT3HcedX7mRY8jCHIzbGhJtgEv8W/4hbr+Gr6gGsOWco1NftQpW3in0n9vHshmdJikni12f/mumDpuMSl8MRG2PCUTCJPx5fwp9SY5415wyB+rpdyFqZRZVWMX3QdOZkzqFrXFeHIjTGdATB9Md/U3NWLCJ9gb8BafjuDTylqnNFJBlf7579gd3Atap6pDnb6Gjq63ahSqv488V/5rw+57VxRMaYjiiYB7iewz/sYk2q+t+NFK0CfqaquSKSBOSIyDLgv4B3VPUhEbkduB34RZMj74Dq63YhNSHVkr4xptUEU0n8Or7eOJcC7wCd8XXS1iBVzVPVXP/nE8BnQG/gcnxNRPG/WxfPfrMzZhPtiq41L84dx0/H/9ShiIwxHVEwVT0v15wWkfnA8qZsRET6A+OAT4BUVc3zrztPRHrWU2YmMBOgX7+26RO+tTSny4UKTwXrCtdR6a1EEBS17pGNMSER3GjYtQ0Bgs7EItIJeBn4iaoeD7ZHSFV9Ct8g72RmZtapaupI8ovzue3929hQtIEbht/AxqKNuMRlbfKNMSHRaFWPiJwQkePVL3zNOoOqkxeRaHxJ/+81mn8eFJF0//fpQEHzQu8Y1hWu4ztLv8OOozt49IJHmTNhjjXTNMaEVDBVPc3ql198l/bPAp+p6iM1vloC3Ag85H9f3Jz1dwRLdiwha2UWqQmpPH3p0wzuNtjpkIwxEaChEbgarM5R1b2NrPts/F09iMha/7w78CX8hSLyfWAvcE3w4XYMqspf1v2FJ9Y9wcS0ifzh/D9Y23xjTJtp6Ip/KV8OvVhNgRSgJ+BuaMWquuK0sjVd3IQYO5RKbyW//vjXvLL9Fa4YfAX3TL6nTkseY4wJpYaGXhxVc9rfMucXwCXAAyGNqoMqrSrlZ+//jP8c+A+zRs/ilrG3BBz+0G7qGmNCKZgHuIYAdwJfAf4A3KqqlaEOrKMprizm5uU3s7ZwLXdPuptrh13rdEjGmAhVb/MRERnpb7P/Mr52+yNV9RlL+g2r7l0z+2A2UxZNYenOpRwrP8bMt2eyrnAdvz33t5b0jTGOEtXATeRFxAPsw1fX7zn9e1W9NbShfSkzM1Ozs7PbanPNVt27Zs2O1mLdsSTHJVNUWsTvz/89F/W7yMEIjTGRRERyVDXz9PkNVfU01hePOU2g3jXLPeXkFefx+CWPc07vcxyKzBhjvtTQzd159X1nAquvd03Akr4xpt2wR0RbUVpiWsD56YnpbRyJMcbUzxJ/K5qdMZs4d1yteXHuOGZnzHYoImOMqcsSfyu6bMBljOs57tR0emI6WV/Nst41jTHtSjDt+IcCj+PrTnmkiIwGpqvq/SGPLsw8t+k5Ps77mNSEVPom9bUHsYwx7VIwV/xPA78EKgFUdT0wI5RBhaO3dr/FozmPMrX/VPp06uN0OMYYU69gEn+Cqq4+bV5VKIIJV58d+oy7VtzF2JSx/Oac3wTshsEYY9qLYBJ/kYgMwj/urohcDdQdGDZCHSo9xK3v3UqX2C48euGjxLhjnA7JGGMaFMwIXLfgGwnrTBE5AOwCvhfSqMJEpaeS296/jSNlR5j39Xn0iO/hdEjGGNOoYAZi2QlcIiKJgMs/cLoBHsl5hNyCXB469yFGdB/hdDjGGBOUYFr1xAJXAf2BqOr6a1W9L6SRtXPv7HmHFz97kevOuq5Oc01rzWOMac+CqepZDBwDcoDy0IYTHvaf2M/dH93NyO4juW38bU6HY4wxTRJM4u+jqlNDHkk7c9ObNwF1r94rPZXM+WAOAA+f/7DdzDXGhJ1gWvWsFJFRjS8WGf706Z/YeGgjvz771/RJsvb6xpjw09Bg6xsBr3+Zm0RkJ76qHgFUVUe3TYjtR87BHJ7f9DzXDL2Gi8+I2GGDjTFhrqGqnt7A2LYKpL0rrizmzhV30rtTb/4383+dDscYY5qtocS/S1X3tFkk7dzDax4mrziP56c+T0J0gtPhGGNMszWU+HuKSL1NVlT1kRDE0y5Uj5tb4a1gyqIpTDljCi9//jLfH/n9Wr1vGmNMOGoo8buBTvjq9CNG9bi5Fd4KAPKK85i3eR4943ty89ibHY7OGGNarqHEn9eSh7RE5K/AN4ACVR3pn5cMvITvYbDdwLWqeqS52wiFQOPmAnjUY003jTEdQkPNOVt6pf88cHr7/9uBd1R1CPCOf7pdqW/c3MNlh9s4EmOMCY2GEn+L2iuq6ofA6dnycqB6EPd5wBUt2UYo1Ddubn3zjTEm3NSb+FU1FJe4qaqa519/HtCzvgVFZKaIZItIdmFhYQhCCczGzTXGdHTtdsxdVX1KVTNVNTMlJaXNtjtt4DR+lvmzU9M2bq4xpqMJpq+e1nRQRNJVNU9E0oGCNt5+ULYe2YogDO8+nAXfWOB0OMYY06ra+op/CXCj//ON+Hr+bFe2HN7Cy9teJiUhhfioeKfDMcaYVheyxC8i84GPgWEisl9Evg88BFwqIp8Dl/qn2w1V5berf0vX2K70SuzldDjGGBMSIavqUdXv1PNVu+3dbPne5WQfzObuSXfz713/djocY4wJiXZ7c7etVXmrmJs7l8FdB3PVkKucDscYY0KmrW/utluLty9mz/E9PHbhY7hdbhs+0RjTYdkVP1DuKefxdY8zOmU0F/S9wOlwjDEmpCzxAwu2LOBgyUFmj5tN9WDyxhjTUUV84j9ZcZJnNjzD5PTJTEyf6HQ4xhgTchGf+F/47AWOlh+1LhmMMREjohN/cWUxL25+kQv7XsiIHiOcDscYY9pERCf+hVsXcrziOD8Y9QOnQzHGmDYTsYm/rKqMeZvmMSl9EqNSRjkdjjHGtJmITfyvbn+VQ2WHmDl6ptOhGGNMm4rIxF/preS5jc8xNmUsmamZTodjjDFtKiIT/xs73+CL4i/4wegfWLt9Y0zEibjEr6o8v+l5hnYbyrm9z3U6HGOMaXMRl/hX569m+9HtfO+s79nVvjEmIkVc4n9x84skxyVz2cDLnA7FGGMcEVGJf+/xvXyw/wOuGXoNse5Yp8MxxhhHRFTi/8eWf+B2ufn2sG87HYoxxjgmYhL/yYqTvLr9Vab2n0pKQorT4RhjjGMiJvG/uv1ViiuL+d5Z33M6FGOMcVREJH5V5aWtLzE6ZbR1xmaMiXgRkfhzC3LZfXw31wy9xulQjDHGcRGR+F/e9jKdojsx5YwpTodijDGO6/CJ/1j5Md7e8zaXDbiMhOgEp8MxxhjHdfjEv3TnUso95Vw19CqnQzHGmHbBkcQvIlNFZKuIbBeR20OxjaU7lzJl0RQeXP0gUa4odh3bFYrNGGNM2Ilq6w2KiBv4M3ApsB9YIyJLVHVza21j6c6lZK3MosxTBkCVt4qslVkATBs4rbU2Y4wxYcmJK/6JwHZV3amqFcAC4PLW3MDc3Lmnkn61Mk8Zc3PntuZmjDEmLDmR+HsD+2pM7/fPq0VEZopItohkFxYWNmkD+cX5TZpvjDGRxInEH6gvZK0zQ/UpVc1U1cyUlKZ1sZCWmNak+cYYE0mcSPz7gb41pvsAX7TmBmZnzCbOHVdrXpw7jtkZs1tzM8YYE5ba/OYusAYYIiIDgAPADOC7rbmB6hu4c3Pnkl+cT1piGrMzZtuNXWOMwYHEr6pVIvJj4C3ADfxVVTe19namDZxmid4YYwJw4oofVX0DeMOJbRtjTKTr8E/uGmOMqc0SvzHGRBhL/MYYE2Es8RtjTIQR1TrPTrU7IlII7Glm8R5AUSuG46SOsi8dZT/A9qW96ij70tL9OENV6zwBGxaJvyVEJFtVM52OozV0lH3pKPsBti/tVUfZl1Dth1X1GGNMhLHEb4wxESYSEv9TTgfQijrKvnSU/QDbl/aqo+xLSPajw9fxG2OMqS0SrviNMcbUYInfGGMiTIdJ/I0N4C4isSLykv/7T0Skf9tH2bgg9uO/RKRQRNb6X//jRJzBEJG/ikiBiGys53sRkcf8+7peRDLaOsZgBLEfF4jIsRrH5J62jjFYItJXRN4Tkc9EZJOI1BmkIhyOS5D7ERbHRUTiRGS1iKzz78uvAizTuvlLVcP+ha975x3AQCAGWAcMP22Zm4En/J9nAC85HXcz9+O/gP9zOtYg9+c8IAPYWM/3lwH/xjcq2yTgE6djbuZ+XAC87nScQe5LOpDh/5wEbAvwf6zdH5cg9yMsjov/37mT/3M08Akw6bRlWjV/dZQr/mAGcL8cmOf/vAi4WEQCDQPppJAPRN+WVPVD4HADi1wO/E19VgFdRSS9baILXhD7ETZUNU9Vc/2fTwCfUXfM63Z/XILcj7Dg/3c+6Z+M9r9Ob3XTqvmroyT+YAZwP7WMqlYBx4DubRJd8IIaiB64yv8TfJGI9A3wfbgIdn/DwWT/T/V/i8gIp4MJhr+6YBy+K8yawuq4NLAfECbHRUTcIrIWKACWqWq9x6Q18ldHSfzBDOAe1CDvDgsmxteA/qo6GljOl1cB4SgcjkkwcvH1iTIG+BPwqsPxNEpEOgEvAz9R1eOnfx2gSLs8Lo3sR9gcF1X1qOpYfGOQTxSRkact0qrHpKMk/mAGcD+1jIhEAV1ofz/fG90PVT2kquX+yaeB8W0UWygEc9zaPVU9Xv1TXX2jy0WLSA+Hw6qXiETjS5Z/V9V/BVgkLI5LY/sRbscFQFWPAu8DU0/7qlXzV0dJ/KcGcBeRGHw3P5actswS4Eb/56uBd9V/p6QdaXQ/TqtrnY6vbjNcLQFu8LcimQQcU9U8p4NqKhFJq65vFZGJ+P6uDjkbVWD+OJ8FPlPVR+pZrN0fl2D2I1yOi4ikiEhX/+d44BJgy2mLtWr+cmTM3dam9QzgLiL3AdmqugTff5IXRGQ7vjPlDOciDizI/bhVRKYDVfj2478cC7gRIjIfX8uKHiKyH7gX340rVPUJfOMuXwZsB0qAm5yJtGFB7MfVwI9EpAooBWa0w4uKamcD1wMb/HXKAHcA/SCsjksw+xEuxyUdmCcibnwnp4Wq+noo85d12WCMMRGmo1T1GGOMCZIlfmOMiTCW+I0xJsJY4jfGmAhjid8YYyJMh2jOaUxrEZHuwDv+yTTAAxT6p0tU9auOBGZMK7LmnMbUQ0SygJOq+nunYzGmNVlVjzFBEpGT/vcLROQDEVkoIttE5CERuc7fp/oGERnkXy5FRF4WkTX+19nO7oExPpb4jWmeMcBsYBS+J0iHqupE4Bng//mXmQs8qqoTgKv83xnjOKvjN6Z51lT3XyMiO4C3/fM3ABf6P18CDK/RbXpnEUny9x9vjGMs8RvTPOU1PntrTHv58u/KBUxW1dK2DMyYxlhVjzGh8zbw4+oJERnrYCzGnGKJ35jQuRXI9I+Wthn4odMBGQPWnNMYYyKOXfEbY0yEscRvjDERxhK/McZEGEv8xhgTYSzxG2NMhLHEb4wxEcYSvzHGRJj/D4+1HGpbBv65AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret2 = ensemble_simulations(duration, ndiv=20, y0=y0, model=m, solver='gillespie', repeat=N)\n", "ret2.plot('o', ret1, '-', **opt_kwargs)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Simulating with `meso`:" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEGCAYAAABiq/5QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXhU5dn48e89kx0CARII+w5VthBAoWqLO5W+iHWjtdXyWrFqf+BrS901tVZtrbbYahW1iqgoYhUEK4JrKaKQsC8i+5aQsAayJ3P//pgJJmQymSyTyczcn+uaa+acOc859+GEe855znOeR1QVY4wxkcMR7ACMMcY0L0v8xhgTYSzxG2NMhLHEb4wxEcYSvzHGRJioYAfgj+TkZO3Vq1ewwzDGmJCSmZl5SFVTTp8fEom/V69erFq1KthhGGNMSBGR3d7mW1WPMcZEGEv8xhgTYSzxG2NMhAmJOn5jjAmG0tJStm/fTmFhYbBD8SkhIYG+ffsSExPj1/KW+I0xphbbt28nKSmJgQMH4nC0zAoSl8vFwYMH2bZtG2eeeaZfZQK6JyKSJCLzRGSLiGwWkTEi0l5ElojIN573doHY9soFz5GT0Q/Xg23JyejHygXPBWIzxpgwVlhYSKdOnVps0gdwOBx06tSJwsJCcnJy/CsT4JhmAB+o6neAYcBm4C7gI1XtD3zkmW5SKxc8x+DM+0glD4dAKnkMzrzPkr8xpt5actKv5HA4EBHmz59PRUVF3csHKhARaQN8D3gRQFVLVfUYcDkwy7PYLGBiU2+7e9bjxEtptXnxUkr3rMebelPGGFPNtc99wbXPfRGUbRcXF1NUVFTncoH8KesD5AEvichqEXlBRFoBnVQ1G8Dz3tFbYRGZIiKrRGRVXl5evTbcUb0v31EP1Ws9xhjTErzzzjuICFu2bPG5nIj4tb5AJv4oIB34h6oOBwqoR7WOqs5U1ZGqOjIlpcYTxz7livflcyW5Xusxxpj6eHf1flbvOcaXO49wzmMf8+7q/U2y3jlz5nDuuefyxhtvNMn6Apn49wH7VPVLz/Q83D8EB0WkM4DnPbepN7w3fTpFWr1ZU5HGsDd9elNvyhhjAHfSv/tf6ymtcAGw/1gRd/9rfaOT/8mTJ/nvf//Liy++2PITv6rmAHtFZKBn1oXAJmABcINn3g3A/Kbe9qgJN7NhxMPkkIJLIcfVji/7/x+jJtzc1JsyxhgAHl/8NUVl1W+sFpVV8Pjirxu13nfffZdx48YxYMAA2rdvT1ZWVqPWB4Fv1fP/gNdEZB2QBjwCPAZcLCLfABd7ppvcqAk3k5qxjUNT1pAix4gt9K+ZkzHGNMSBY95vqtY2319z5sxh0qRJAEyaNIk5c+Y0an0Q4Ae4VHUNMNLLVxcGcrtVdezamzWtRtP/wHzKSh8nOia2uTZtjIkgXZLi2e8lyXdJim/wOg8fPszHH3/Mhg0bEBEqKioQEf70pz/5fSPXm5bfQLUpDL+eZI6x4dO3gh2JMSZMTb90IPHRzmrz4qOdTL90YC0l6jZv3jyuv/56du/eza5du9i7dy+9e/dm2bJljYo1IhL/4LFXkUt7HFkvBzsUY0yYmji8K4/+aAgxTnda7ZoUz6M/GsLE4V0bvM45c+ZwxRVXVJt35ZVX8vrrrzcq1ojoqycqOobt3X/E2XteZO+29XTvNyTYIRljwtDE4V2Z89UeAN68eUyj1/fpp5/WmDd16tRGrzcizvgB+o+fRjkODiz+a7BDMcaEsTdvHtMkST+QIibxJ6f2YG3ShQzOXUj+scPBDscYY4ImYhI/QNL5U2klxWx6/5lgh2KMMUETUYm/f9p5bI4+k+7fzKaivDzY4RhjTFBEVOIHKEqfQlc9yLqPm+bRZ2OMCTURl/iHXnQdOaQQs8r65jfGBMBL492vFiziEn9UdAy7+v+MQaXr+HrVx8EOxxhjfHI6naSlpTFs2DDS09NZvnx5o9cZcYkfYMiEaRyjNUUf/ynYoRhjwsm6ubBvJexeBn8Z7J5upPj4eNasWcPatWt59NFHufvuuxu9zohM/K0Sk9jc4zrSCr9gx4Yv6y5gjDF1WTcX3psKFSXu6eN73dNNkPwr5efn065d44cpj4gnd705c+J0CmbM4ujix2DwO8EOxxgT6j56CMpO66StrMg9f+g1DV5tUVERaWlpFBcXk52dzccfN76KOiLP+AHatk9hXZerSMv/hL3b1gc7HGNMqDu+r37z/VRZ1bNlyxY++OADrr/+elS1UeuM2MQP0P/yuygjiuyFjwQ7FGNMqGvbrX7zG2DMmDEcOnSI+o5DfrqITvzJqd1ZmzKB4UcXs3/HxmrfbXzkXDY+cm6QIjPGhJwLH4Do0/rej453z28iW7ZsoaKigg4dOjRqPRGd+AH6/uhBynGS/W7THRxjTAQaeg38z1Pg9Az21La7e7oR9fvwbR1/Wloa1157LbNmzcLpdNZd0IeIvblbKblLT77oMomzD8xm58Yv6T3o7GCHZIwJVUOvgcxZ7s+TFzXJKisqKupeqJ4i/owf4MyrH+CkxHN80YPBDsUYE+omL2qypB8olvhxt/DZ1GsyaYVfsGXlUlYueI72Jfs5o2Q9ORn9WLnAuncwxoQPS/weQ6+6k0MksXvR4wzOvI/OcgSHQCp5DM68z5K/MSZsWOL3SGjdlu1n3Mow19fES2m17+KllO5ZjwcpMmOMaVqW+KsY8aP/o6Mc9fpdRz3UzNEYY0xgWOKvIio6hoOS7PW73FrmG2NMVZM/mMzkDyYHOwyfApr4RWSXiKwXkTUissozr72ILBGRbzzvje9xqAntS/8tRRpTbV6RxrA3fXqQIjLGRLKcnBwmTZpE3759OfPMM7nsssvYunVro9bZHGf856tqmqqO9EzfBXykqv2BjzzTLcaoCTfzZb//I9vVDpdCDilsGPEwoybcHOzQjDEt3KIdi1iXt45VB1dxybxLWLSjcc06VZUrrriCsWPHsn37djZt2sQjjzzCwYMHG7XeYFT1XA54nnBgFjAxCDH4NPZn97BZeuMQOPqDpy3pG2PqtGjHIjKWZ1DqcjcOyS7IJmN5RqOS/yeffEJ0dDS//OUvT81LS0vjvPPOa1SsgU78CnwoIpkiMsUzr5OqZgN43jt6KygiU0RklYisamyHRA2RFF3BAW1P/OJfU1pS3OzbN8aElhlZMyiuqJ4riiuKmZE1o8Hr3LBhAyNGjGhsaDUEOvGfo6rpwA+A20Tke/4WVNWZqjpSVUempKQELsJaxDpcbJQB9HLtJfN1e6LXGONbTkFOveYHU0ATv6oe8LznAu8AZwEHRaQzgOc9N5AxNNSge5ZxccZiMhPPJ33Xi+z9Zm2wQzLGtGCprVLrNd8fgwYNIjMzs8HlaxOwxC8irUQksfIzcAmwAVgA3OBZ7AZgfqBiaAo9r3uKEokh/61f4QpAZ0nGmPAwLX0acc64avPinHFMS5/W4HVecMEFlJSU8Pzzz5+at3LlSj777LMGrxMCe8bfCVgmImuBr4BFqvoB8BhwsYh8A1zsmW6xklN7sGXIdAaVrmPlPHt61xjj3fg+48n4bgYxDndz8M6tOpPx3QzG9xnf4HWKCO+88w5Lliyhb9++DBo0iIyMDLp06dKoWOvslllErgY+UNUTInIfkA48rKpZvsqp6g5gmJf5h4ELGxhvUIy6Yhprv1nE0E1PsPebH9C9f43dMsYYxvcZz7yt8wB4adxLTbLOLl26MHdu0w3YDv6d8d/vSfrnApfiboL5jyaNooUTh4Ou179IicRQ9OYvKC8rrbuQMSYivTTupSZL+oHiT+KvrNgeD/xDVecDMT6WD0vJXXqy/ayHGFC+lZWv3h/scIwxpsH8Sfz7ReQ54BrgfRGJ9bNc2Blx2Y2sSryQkbueZ8uqj4IdjjGmGbhcrmCHUKf6xuhPAr8GWAyMU9VjQHsgYjuu6f+/M8lzdCBp4c0cP9y4x6aNMS1bQkICOTk5LTr5u1wucnJyKCsr87tMnTd3VbVQRHKBc4FvgHLPe0Rq2y6Zg//zPL3m/4hNL97AsN+8jzgi8gLImLDXt29f1qxZw4EDBxCRYIdTq7KyMvbs2YOq4vAjH/nTqudBYCQwEHgJiAZeBc5pZKwha0D6WFZsuYPRWx9nxZzfM/q66k/2bnzkXMD9EJgxJnTFxMQQHx/PkiVLaNu2rV9JNVhOnDhBp06diI+Pr3NZf/biCmACUACnnsZNbFSEYeDsSfewOuEcRmydwaYv/h3scIwxATJ48GDGjh1LXFwcItIiXw6Hg/79+zNhwgS/rkzqPOMHSlVVRUTh1FO4EU8cDvrc9Ao5T51H6uIpZHdZSueeA4MdljGmiYkI6enppKenBzuUJuPPGf9cT6ueJBG5CVgKPF9HmYjQtl0y+uM5RGsZha9MovDkcVYueI72Jfs5o2Q9ORn9bJB2Y0yLI6pa90IiF+Pua0eAxaq6JNCBVTVy5EhdtWpVc26yXtZ+PJchn01hseN7jHV9UW2w9iKNsYFcjDFBISKZVQbBOsWfqh48ib5Zk30oGXbBNazYt57h22YT76j+VG+8lNI963GwxG+MaSFqreoRkRMikl/lPb/qdHMGGQrO/unv6ChHvX7XUQ81czTGGFO7Ws/4VTXiW+7Uhzgc5EgKnak5WliuJNPwHrmNMaZp1XlzV0R6eHs1R3ChZl/6dIq0ejdGRRrD3vSIfdDZGNMC+VPHX3Wk4DigN/A1MCggEYWwURNuZiXQLfMxOnGEg9qODT1/xsVWv2+MaUHqPONX1SFVXv1xD59oj6TWYtSEmzkS25Uvo9JxOWM4e++LbF+/IthhGWPMKfV+/tgzAMuoAMQSVto4y5Ab3qOIONq9fTU7N60MdkjGGAP40Y5fRO6oMunAPQJXB1W9NJCBVdXS2/H7snfbemJfnUA0ZRy+4g36DTu3zjLW148xpinU1o7fnzP+xCqvWNx1/pc3bXjhq3u/IZRdv5Bi4uj4ztVsWbk02CEZYyKcP90y/645AglnXfsMIud//03+S/9Dj4U/YX3R8wz5nv12GmOCw5/mnEtEJKnKdDsRWRzYsMJPao/+xN20mIPOzgz8aDKrFjwb7JCMMRHKn6qeFM/IWwCo6lGgY+BCCl/JXXrSYerHfBM7iJFZd7LilfvRFjyyjzEmPPk12HrVB7ZEpCdQd89uxqs2SR3od8diMhPPZ/SOp1j59+spLSkOdljGmAjiT+K/F1gmIrNFZDbwOXB3YMMKb7FxCQy//W2+6HIDZx15j21/vpAjufsBrFtnY0zA+dstczIwGne3zF+o+t/rmIg4gVXAflX9oYj0Bt7APWh7FvAzVS31tY5Qbs5Zl1ULZzJ45T0clSQ29ryOc3Y9Y906G2OaRIObc4p7HK9xQLqqvgckiMhZ9dj2NGBzlek/An/xPAV8FLixHusKOyN/OIW9E/+FkwoG7ZxVLelDlW6djTGmifhT1fMMMAb4sWf6BPC0PysXkW7AeOAFz7QAFwDzPIvMAibWI96w1H/494i6dRmdmqBb542PnHvqATBjjPHGn8R/tqreBhTDqVY9Mb6LnPJX4LdAZdOVDsAxVS33TO8DunorKCJTRGSViKzKy6vZ1XG4ad+xKwclxet3uZLczNEYY8KZP4m/zFNPXznYegrfJvJaicgPgVxVzaw628uiXm8yqOpMVR2pqiNTUrwnxHDjrVvnQo1hT9qvgxSRMSYc+ZP4nwLeATqKyB9w98z5iB/lzgEmiMgu3DdzL8B9BZAkIpVPDHcDDtQ36HA1asLNbBjxMNnaHpdCtrZno6snbTe8zO4tWcEOzxgTJvzplvk13NU1jwLZwERVfcuPcnerajdV7QVMAj5W1euAT4CrPIvdAMxvYOxhqbJb582xQ0h9cDty1k10Kt9H5zkX88WseygrLQl2iMaYEOdrzN32lS8gF5gDvA4c9MxrqDuBO0RkG+46/xcbsa6wJg4HI//nZsp/uYINiecwZufT7H3sLLZ85X3c+8Y+A2A3ho2JDL46acvEXf9eW718H383oqqfAp96Pu/APZiL8VNyaneSf7OArMWz6frFg/R5/yq++mI8/X/yBO1SOgPupD84875TzUFTyaNt5n2sBHsGwBhTTa1n/KraW1X7eN5Pf/md9E3TSb/0Z7T+dRYrUq9j+JEPcD6dzorXfkdpSTHdsx63ZwCMMX7xawQuEZkgIn/2vH4Y6KBM7VolJjH6l89wYNISdsYNYvQ3T5L7WBod1XuT1/o8A2CMiQz+PLn7GO6nbzd5XtNE5NFAB2Z863nGCIbdtZR1Y/9JuUSRq+28LtdczwDY/QFjQkedA7EAlwFpquoCEJFZwGqso7aAqc+Qi0PHXknFuZfzwd+ncsHRt2r087N3xHRSAxGkMSZk+TvYelKVz20DEYhpOGdUFONvf4Z1aQ9ywOV+BiDH1Y7/kEZMYnvr898YU40/if9RYLWIvOw528/Evwe4TDM7+4pfcTSuK+uih7J70C0MlR0M++wX7Hp4OCvffdpnv//WHbQxkcOfB7jm4O6S+V+e1xhVfSPQgZmGi3YoZ197J+3v3sjKtEcAZdSaezj+6Hf44qU7OZSzt9rylU1BO8sRHOJuCjo48z5L/saEKV8PcKVXvoDOuDtU2wt08cwzLVxMbByjJt5Gr/vWsPb7L5Ad148xu5+lzT+GserJK9m4/H3U5WoRTUHt5rAxzcfXzd1VwEagsp1g1Qe5FHffOyYEiMPBsPOvhvOvZs/WNWQv+Rtn5L1Pmw+XsmdpV7ppntfH9KwpqDHhyVfi/zVwJVCEu5O1d1T1ZLNEZRqsrhZBPQak0WPAixQVnGDlklkkbniV3PJCUr2MBZArydYiyJgw5OvJ3b+o6rnAr4DuwEciMldE0potOhMw8a0SGTXxV3znvhVs6vcLCr10B7065XKKiwqCFKFpiaxKLjz4c3N3J+4eND/E3cfOgEAHZZrXBdffx8bTuoNe7erHD/JeoOKxPmQ+MZHM91/kxPEjXstbi6DQYsnb1FrVIyJ9cHenfDnum7pvAH9Q1drbBJqQNWrCzWzcMJsjdGXQPctILi1h3fKFFK97h35HPqP9V59Q+uV01scNo6DXRXQdNYHu/YZY53BBUpm46/OwnzGVfNXxbwPW4T7bzwd6ALe6h80FVX0y4NGZoImOiWXo2Cth7JWUl5WyKfNj8tcsoEvupwz5+k/w9Z/YJ6l0c5XW3iLIEn/Y+V2y+zbf3CDHYRrHV+J/iG+HRWzdDLGYFioqOoYzR4+D0eMA2LdtA/szFxK362OGFH7ptUx9WgStXPAc3Ur204kj5GT0Y2/69HpdLTTm7DdYZ87BPGMP1eQ9+YPJALw07qVmLRuOak38qprRjHGYENKt32C69RsM3EV2Rj86U7Nn0IOaROHvh3Oo3XCcvcfQZdD36NyjP+KoflsplKuKrLoldBTsWd3gstfMdLdnmTtlTb3LttSTEn86aTOmVvvSp5NUJXGDu3O4Va2+R2/XbobkLSTh0NuwEg7Tlr0JZ1KUMoyEXiPpPui7vh8ea+GJP5hC8ay9MQnUNC1L/KZRRk24mZVAt8zH6MQRciWFvSOm8z+epF1eVsq2TV9x+OvlOPZn0unEerrtWoFj97PwGbhqGePN36qixlYTGROJfLXqmaaqM0TkHFX9b3MGZYKjoZeUp7cIqvrQV1R0DP2GnUu/Yd82Hzxx/Ah7Nq7gxI6v6PnNLDpTs5noQZLY/+RVlHUYQFznM+nQazCde51BdEzsqWVCuZooFM/Yg60x1TXB0pjjHMi/EV9n/JOBGcDfAOubxzSZxLbtGfTdy+C7l7FyQecaVUWFGsNaOYPh+Vl0yl8CO4HlUKZO9jo6cSSuO0WJveh5cGmjqonsasFEKl+Jf7OI7AJSRGRdlfkCqKoODWhkJiLUVlU0zpOATxw/Qvb29Rzfu5Hy3K+JPb6DpMI9DChcTSylXtfZUfNY8cIdONr1ID6lJ2069Sala28SWn87lEQoXy0Y01i+WvX8WERSgcXAhOYLyUQaX1VFiW3bk5j+fUj/frUy6nKR81B/OlPzXsBBbceovf/EuU+rzc+nFUccHciPSaFb0U6vVwvdsv7k99VCl5J9dOaoXS2YkOPz5q6q5gDDRCSGb7tq+FpVywIemTE+iMPBvvTfem1RtG/k3SSP+zkH9+/kWM5OCvN2U3Z0L44TB4gpPEjr0lw6ebmvANBJD1H4YEeOOZI44UyiKLodpbHtqYhvjyR0wNE6mcMHdjM295UGXy2sXPAcf9l/rME/Go0tH4oW7VjENzEVlAOXzLuEaenTGN9nfLDD8qklH6c6W/WIyPeBV4BduKt5uovIDar6eYBjM8an2qqJKv9zden9Hbr0/o7XsjkZ/Uj19vwBHdideglRRXnElhwhsTSXxKKtJB09ToxUuMu62hHv8HK1sOpRdqx5huKoREqiEimPbo0rOhFXXFskNhGJa8uhnN1ccGhOjR+NFeWlnD3xthrPOZyusVVUjUlGwUpki3YsImN5BuWe1l/ZBdlkLM8AaLHJP5jHyR/+NOd8ErhEVb8GEJEBwBxghK9CIhIHfA7EerYzT1UfFJHeuPv9aQ9kAT9TVe+VtcbUwVc1kS9706fT1tvVwog7Ge3lP5i6XJw8eZz8I7mkvnSW13V2kqOsjTuDmPITtC0+QHxhAa0opLUW4hB3tVNtPxq91jyBa+39FBJHkcRTInGUOOIpc8RR5oynwhlHhTOebsczvVdRZf6RVU4Hzuh4nDHxRMUmEBUbT1RMPNGxcUTHJvDNF/MZuemPDUpGjU1kT711B4cdZeRFCRe+MIjL217K1Kv96/VlRtYMiiuqdxFWXFHMjKwZfiX+xmy7oWUb83xKc9x/8ifxR1cmfQBV3Soi0X6UKwEuUNWTnuWXici/gTuAv6jqGyLyLHAj8I+GBG9MQ1VeLXTJfJTOHK1xtXA6cTho3aYdrdu0I0dSvF4t5EoKw3/77xrzXRUVnCzIp/DEMTrO9N6reUc5ypfdJiNlBTjKCnCWF+EsLyS6opC48nyiS/OI1eLaq6g4TOeVv/G90z6uVA5lPUI50ZRJNBXipEKiqZAo97sjmm7Fe/i4dRQz2iWTE+UktbyCaUePkZ75GF/uWYE6osARjTqciCMadUQhjihwOFlyfB0LWm+jONp9NZMbLcw+uZj8f1zDuB4X43BEgSMKhzMKcTg9ryjE6cDhiCKnINvr7uQUZLP76zU4nE4cDgfiiPK8O3A6osDh4JWPHmVO6Wc1tu16Yyq3THwMh8OJiJx6r3rF9dRbdzD75OIaZXnrjjqTf8dGDG7UHA81+pP4V4nIi8Bsz/R1uAdc90lVFagcuCXa86ocuesnnvmzgAws8Ye8UOy2YNSEm7km5x9AEnOnrGn01cLeEdO9rsPhdFb70chsVcCMdknVEuiIglaMuemvdW47J6Of1/LpBQmUTnqLspJCykuKKC8toqK0EFdZCRWlxbhKixi57kEWtUqoUfYHJ4+ysv0EHBWliKsMh6sMh6sUh5bj0HKcrjJWtSrmoeT2FHsSY3Z0FBnJ7XlAj3DuoSVEaTlRVOCk4lSVWKV7u3U7Va5SscPBJ9EbuXflYiqACoEKhArAJZXv7nkpXTqRG1UzVXUoL6fw7YtwAS4El7gTTEWVz4uSU04l7qrbfu/kRwx5qo+nLKgI6vnsjkOY3y7Ja9n5xxfT9onenuUF16myggv354m0Y3Wrkhr/1sNPxvL6X89AwfNyL1+5fRDu0jwWtfZ2nJpuRDx/Ev8twG3AVNy/YZ8Dz/izchFx4v6R6Ac8DWwHjqlquWeRfUDXWspOAaYA9OjRw5/NmRD25/T+ADSkC63GlG2IURNu5qmSr5l/fDF5UUJKubqrAPw4G5s7aIL7LPK0BPqzXpcy1Y9t+yr/q4FplFSUUFpRSklFSbXPpRWlvL6jEzOSo2uUPUprXBdeTKmrlDJXGWUVZZS7yt2fXe7P/9n0ltfk/df27XlryAjKXeXul5ZTXuH+XKHllLsqyCs86HVf8qKEob3r/r89+Hhr8hPzq20/zuUitSCJa7r6TmGi6nV+XpRwe6eUBpf9c7Lvfiv/E5XMtioxV/5b94ttw4a2vgcy7FaSwt+SY2uUPUY51/ks6b86E7+qluCu5693N8yqWgGkiUgS8A5whrfFaik7E5gJMHLkSO9HwJggWLRjEbNL/lO9CqDkP/TdsajOOueFssFrAp1XsYre29+jqLyIwrJC93u5+73q66viryj1Uv75wiU8/8own9uO6xBPscNVo+yMDgkUr3r81DyHOIh2RBPliCLGEUO0I5rcKO83nXOjHPR2xJAQlYDT4SRKoohyRJ367HQ4Wb75X+RG16z3SClXrhp5Gw5x4HQ4cYrT/bnqu8PBOa/fxeqSwppnzwXFbJ44Ewfu8g5xIMipsiLCXQt/Xuu2n/nR26f2t7KsiJz6PHnuD2ot++Z1nwKc2l5lWRFBEH40/0cUF1ZP8MUOB4e6tCbris9B3GVPlfO8A4yd/V2KXSdqlH0+tUPzJf6moKrHRORTYDSQJCJRnrP+bsCB5ojBBFYodnvbkCaCqspfM//q9WbjH1b8ga1Ht3Ki9MS3r7ITnCw96X6VnaSwvNDreo+WHOWeZfdUmxfnjCM+Kp74qHjiouKIi4qj1FV7O4hbht1CjDOGOGccsVGxxDhiiI2KJdYRS4wzhl8u/aXXciVOZfmPlxPtiCbaEY3T4ayxzCXzLiHbS11759adeeHSF2qNCeCpAwXVrlLAfcZ+edtLuSXtFp9lAVz6C8YXwPiCwtPmF9Glx4U+y17e9tJatz2w/cAGl02OT/ZZ9mAtVzkHCw4S7fR9i/SIy/sVQW3zGyJgiV9EUoAyT9KPBy4C/gh8AlyFu2XPDbgHejGmWXlrIvjAfx9g46GN9Enqw+GiwxwpPsLR4qMcKXG/Hys+xrGSY7Um3xNlJ3hl0yu0iWlDm5g2tI5uTeuY1nRK6HTq8zvfvMPJspr/gZPjk3l53MskRCWQEJ1AnDOufgm4VWduTbvV5z53btXZa9nUVqkkxiT6LDstfRoZyzOq/eDFOeOYlvbvJRcAABb/SURBVD7NZznAfSP0rTtqVo352bImt9ab6cl13pdpzLYbUza1VWqt/9aBLOsvn4nfU0f/mKpOb8C6OwOzPOtwAHNVdaGIbALeEJGHgdXAiw1YtzGAO4Gvy1tHqau0xlm7qnKk+Ag5hTnkFORwsOAguYW55BbmsnjX4hoJvNRVyuzNs09NJ0Yn0i6uHe3i2tGldRcGdRhEUlwS87bO40Rp9UtxgNSEVD686sNTl+zeDOowyGsC/c3I39CzTc8697cxCbgxZSv/Te/7/C7KcZ/p1+chqqlXP8mymWl0KK1/t8z1vZnelNtuaNlgHSd/1fXkboWIjBAR8bTS8ZuqrgOGe5m/A/eg7cY0ysLtC8n4IuNUAs8uyObeZffy/LrnKddysk9m10juUY4oUuJTaj1rF4QPr/qQ9nHtiXHGeF1mYLuBXv9j3j7idp9JHxqfQBtTvim2PWvpvQDMvepDv8o0hfo2vW0Jgnmc/OFPVc9qYL6IvAUUVM5U1X81WRTG+HCy9CS78nex8/hOduXvYnf+bvbk72HLkS3oaW0DKrSCPSf2cH738zm/+/mktkr99pWQSru4djjEUWuVSeWyvgQ7gTamfLCSd2M1tOltMLXk4+RP4m8PHMbd/r6SApb4TZOoWl0z9s2xfL/b90mITmDbsW3sOLaD3KLcU8s6xUnX1l3p0aZHjaRfqdxVzhNjn/C5zcZeTo/vM555W+cBoXVD2xjwrznn5OYIxESWQ0WH2HR4E/O3zWfpnqW41N3M8HDxYf617V9ESzQD2g9gdJfR9G7bmz5t+9CrbS+6t+5+qlWEr7P2ulSenT/w3wcodZXSuVXTX04b01L500nbANxP1nZS1cEiMhSYoKoPBzw606wC1SSzoKyAjYc2su7QOtbnrWfj4Y21NnerlJyQzBs/fMPnMqF81t6qR43bX6YODx7y/dBUS9SY4xzIvxF/qnqeB6YDz4H7pq2IvA5Y4jdAzZY1Px/0c9rFtSPrYBZr8taw9ejWU2f0Pdv0ZESnEQzqMIhByYOY/MFkr1U2OQU5dW7XztqNaRh/En+Cqn51WmuF8toWNpFl0Y5FPLj8wWotax796lEA4qPiGZYyjClDpzAsZRhDkofQNrZttfKNbbMcqXXtoXjFUN+mlCZw/En8h0SkL56uFUTkKsB7d3kmIpRWlJKVm8Wyfct4bctrlLtqngckxyez5KolRDl8/4k1R5vlQImkH5pI1pgf2cb8jQTy78ufxH8b7j5zviMi+3EPfd1UXUaYEHGo6BD/2fcfPtn7CSuyV1BUXkS0I9pr0gc4XHS4zqQPkVtdE8wfjVC8WoDG9QDbmKuNcPyB96dVzw7gIhFpBThUteYjiyYsHTh5gCW7l7B091LW5q1FUVJbpTKh7wTO63oeo1JHMXH+xEY/Xt7Y6pqWelZlTEvlT6ueDsCDwLmAisgy4CFVPRzo4EzzK60oZfam2fx7579Zf2g94H5S9Za0Wzi/+/kMbDew2tOpoVxVE8qC9YMVqlcLpjp/qnrewN0H/5We6euAN3F3umbCQGFZIX9e9WcyD2aiKOsOraNLqy7cnn47F/e8mB5tau8zPVKrakKZXeUYv57cVdXfV5l+WEQmBiog0zxUlfWH1jNv6zwWbl9ImZZV+/5I8RFSW6X6TPqVIrVlTSSy4xse/En8n4jIJGCuZ/oqYFHgQjKNUddDWIVlhby3/T3mbp3L1qNbiY+KJ8oZRVl59cRfn8GsjTGhpdbELyIncDfhFNwDpL/q+cqBeyzdBwMenWkye/P38vqW13l327ucLDvJGe3P4IExD3BZ78sY8/oYr2X8eYjKGBN6ak38qup7ZAYTEtblrePljS+zdPdSnA4nl/S8hJ+c8ROGJg89dZO2OQZ+qItVIRjTfPwagcvTP0+vqstbt8wtl6ryZfaXPLfuOVbmrCQxOpEbh9zIT77zE1ISag4wbS1zjIks/jTn/CcwFNgIVI7UbN0yt0CqyvGS42QXZPOLD39BSnwKvxn5G64acBWtolvVWs5a5hgTWfw54x+tqmcGPBLTKGvz1nL/svvZmb8TgDYxbZiaPpWJ/fxrgGUtc4yJHP4k/i9E5ExV3RTwaAxQv+6Rd+fv5slVT/Lx3o+rzc8vzecPK/5AtCPaztyNMdX4k/hn4U7+OUAJ7lY+qqpDAxqZ8el4yXGeXfssb2x5gxhnDInRiZwoq96bRnM2ybSrBGNChz+J/5/Az4D1fFvHb4LEpS7mb5vPXzL/wvHS41zR7wp+NfxXXDD3Aq/LW5NMY8zp/En8e1R1QcAjMXXaenQrD694mNW5qxnecTj3nn0vA9sPBFpGk0xjTGjwJ/Fv8Yy49R7uqh7AmnM2p9KKUmaum8mL618kMSaR35/zeyb0nYBDHKeWsSaZxhh/+ZP443En/EuqzLPmnM1kw6EN3P/f+9l2bBsT+k5g+sjpJMUl1ViuKZpkWj29MZHBn/74JzdkxSLSHXgFSMV9b2Cmqs4Qkfa4e/fsBewCrlHVow3ZRjhTVZ5Z8wwz180kOT6Zpy98mu91+57PMtYk0xjjD38e4HoJao6Grar/W0fRcuDXqpolIolApogsAX4OfKSqj4nIXcBdwJ31jjxMLdqxiLV5aylzlZGZm8nwlOE8fdHTJMZYDxrGmKbhqHsRFuLujXMR8BHQBncnbT6paraqZnk+nwA2A12By3E3EcXzbl08eyzasYj7/3s/Za5ve8rcfGQzn+/7PIhRGWPCjT9VPW9XnRaROcDS+mxERHoBw4EvgU6qmu1Zd7aIdKylzBRgCkCPHnX3CR/qSitKeeiLh6olfbDukY0xTc+vTtpO0x/wOxOLSGvgbeB2Vc2vOmyfL6o6E/cg74wcObJGVVM4ySnI4Y5P76CwvLDW7/1ldfvGmLr4U8df2S9/pRz8rJMXkWjcSf+1Ks0/D4pIZ8/Zfmcgt54xh5W1eWu5/ZPbKSwrpF1sO46W1LzPbW3xjTFNqc46flVNVNU2VV4DTq/+8Ubcp/YvAptV9ckqXy0AbvB8vgGY35DAw8GC7QuY/MFk4pxxvHbZa9x51p3EOeOqLWNt8Y0xTc3XCFw+q3NUdU8d6z4HT1cPIrLGM+8e4DFgrojcCOwBrvY/3PCgqjyz9hmeXfssZ6WexRPff4KkuCT6tesHWPfIxpjA8lXVs4hvh16spEAK0BFw+lqxqi47rWxVF9YjxrBS5irj91/8nne2vcPEfhN5YMwDRDuiT31vbfGNMYHma+jFIVWnPS1z7gQuAh4JaFRhqqi8iF9/+mv+s/8/3Dz0Zm5Luw1/b3YbY0xT8efmbn/gXuBs4AlgqqqW+S5lTldQVsCtS29lTd4a7h99P9cMvCbYIRljIpSvOv7BuBP+IOBPwI2qWtFcgYWT4yXHuXXprWw8vJE/nvdHxvUeF+yQjDERzNcZ/1pgL+66/rOAs6pWS6jq1MCGFh6OFR/jpiU3sf3Ydp4c+yQX9PDeb74xxjQXX4m/rr54TB3yS/OZsmQKO47t4KkLnuLcrucGOyRjjPF5c3dWbd8Z3yZ/MJkKVwUuXHxz7BtmnD+jXknfWvMYYwKpIV02mDq41MW2Y9soLC/k8e8/Xmd3ysYY05ws8Tex97a9x+rc1ShKUmwSpRWlwQ7JGGOq8adbZuOnhdsXcv/y+1FP10bHSo6RsTyDRTsWBTkyY4z5Vp2JX0QGiMhHIrLBMz1URO4LfGih55GvHqHitBavld0qG2NMS+HPGf/zwN1AGYCqrgMmBTKoULR412JOlJ7w+l19ulU2xphA8yfxJ6jqV6fNKw9EMKFq8+HN3Lfsvmp97lRl3SobY1oSfxL/IRHpi6dPfhG5CsgOaFQh5HDRYaZ+MpW2sW357ajfWrfKxpgWz59WPbfhHgnrOyKyH9gJ/DSgUYWIsooy7vj0Do4WH2XWD2YxqMMgEmMSrVtlY0yL5s+YuzuAi0SkFeDwDJxugCcznyQrN4vHznuMQR0GAdatsjGm5fOnd85Y4EqgFxBV2V+Pqj4U0MhauI92f8Srm1/lujOuszN6Y0xI8aeqZz5wHMgESgIbTmjYd2If9//3fgZ3GMwdI+6o8b2d6RtjWjJ/En83VbV+hD3KKsqY/tl0AB7//uPEOGOCHJExxtSPP616lovIkLoXiwx/W/03NhzewO/P+T3dErsFOxxjjKk3XwOxbABcnmUmi8gO3FU9AqiqDm2eEFuOzIOZvLzxZa4ecDUX9ozYYYONMSHOV1VPVyCtuQJp6QrKCrh32b10bd2V34z8TbDDMcaYBvOV+Heq6u5mi6SFe3zl42QXZPPyuJdJiE4IdjjGGNNgvhJ/RxGp2WTFQ1WfDEA8LdLn+z7n7W/e5sbBNzK84/Bgh2OMMY3iK/E7gda46/QjVmFZIQ+veJh+Sf24Ne3WYIdjjDGN5ivxZzfmIS0R+SfwQyBXVQd75rUH3sT9MNgu4BpVPdrQbTSHv6/5O9kF2cz+wWxrummMCQu+mnM29kz/ZeD09v93AR+pan/gI890izT5g8lc8941vLb5Na4deC1pHe0+tzEmPPhK/I1qr6iqnwNHTpt9OVA5iPssYGJjthFIqsru/N10iOtgvWsaY8JKrYlfVU9P2k2hk6pme9afDXSsbUERmSIiq0RkVV5eXgBC8S23MJfC8kLuPvtuEmMSm337xhgTKC12zF1VnamqI1V1ZEpKSrNu+3DRYQ4UHKBNTBsu6nFRs27bGGMCrbkT/0ER6Qzgec9t5u375bef/5YKrSC/NJ9L377UBks3xoSV5k78C4AbPJ9vwN3zZ4vy/Lrn+Srn25EmswuyyVieYcnfGBM2Apb4RWQO8AUwUET2iciNwGPAxSLyDXCxZ7rFUFWeXftsjfnFFcXMyJoRhIiMMabp+dMtc4Oo6o9r+arF9m62dM9SSl2lXr/LKchp5miMMSYwWuzN3eZW7ipnRtYMosT7b2Fqq9RmjsgYYwLDEr/H/G3z2Z2/mx9/58fEOeOqfRfnjLO2/MaYsGGJHyipKOEfa//B0JShTB81nYzvZhDjcHfP0LlVZzK+m2Hj6hpjwkbA6vhDyRtb3uBg4UEeOfcRRITxfcYzb+s8wMbPNcaEn4g/4z9ZepIX1r/AmM5jOKvzWcEOxxhjAi7iE//szbM5VnLM6vCNMREjoqt6CsoKeHXTq5zf/XwGJQ+q9p1V8RhjwlVEn/HP/Xou+aX53DTkpmCHYowxzSZiE39xeTGzNs5idOfRDEkZEuxwjDGm2URs4n9327scLj7MlKFTgh2KMcY0q4hM/GWuMl7a8BJpKWmM7DQy2OEYY0yzisjE//6O9zlQcICbht6ESESPJW+MiUARl/hVlZc3vsyAdgM4r+t5wQ7HGGOaXcQl/q9yvmLbsW389Iyf2tm+MSYiRVzif3XTq7SPa89lfS4LdijGGBMUEZX49+Tv4bN9n3H1gKuJdcYGOxxjjAmKiEr8r295HafDybUDrw12KMYYEzQRk/hPlp7k3W3vMq7XOFISUoIdjjHGBE3EJP53t71LQVkBPz3jp8EOxRhjgioiEr+q8ubXbzI0ZWiNztiMMSbSRETiz8rNYlf+Lq4ecHWwQzHGmKCLiMT/9ta3aR3dmkt6XhLsUIwxJujCPvEfLznOh7s/5LLel5EQnRDscIwxJujCPvEv2rGIkooSrhxwZbBDMcaYFiEoiV9ExonI1yKyTUTuCsQ2Fu1YxCXzLuHRrx4lyhHFzuM7A7EZY4wJOc0+9KKIOIGngYuBfcBKEVmgqpuaahuLdiwiY3kGxRXFAJS7yslYngHA+D7jm2ozxhgTkoJxxn8WsE1Vd6hqKfAGcHlTbmBG1oxTSb9ScUUxM7JmNOVmjDEmJAUj8XcF9laZ3ueZV42ITBGRVSKyKi8vr14byCnIqdd8Y4yJJMFI/N76QtYaM1RnqupIVR2ZklK/LhZSW6XWa74xxkSSYCT+fUD3KtPdgANNuYFp6dOIc8ZVmxfnjGNa+rSm3IwxxoSkZr+5C6wE+otIb2A/MAn4SVNuoPIG7oysGeQU5JDaKpVp6dPsxq4xxhCExK+q5SLyK2Ax4AT+qaobm3o74/uMt0RvjDFeBOOMH1V9H3g/GNs2xphIF/ZP7hpjjKnOEr8xxkQYS/zGGBNhLPEbY0yEEdUaz061OCKSB+xuYPFk4FAThhNM4bIv4bIfYPvSUoXLvjR2P3qqao0nYEMi8TeGiKxS1ZHBjqMphMu+hMt+gO1LSxUu+xKo/bCqHmOMiTCW+I0xJsJEQuKfGewAmlC47Eu47AfYvrRU4bIvAdmPsK/jN8YYU10knPEbY4ypwhK/McZEmLBJ/HUN4C4isSLypuf7L0WkV/NHWTc/9uPnIpInIms8r18EI05/iMg/RSRXRDbU8r2IyFOefV0nIunNHaM//NiPsSJyvMoxeaC5Y/SXiHQXkU9EZLOIbBSRGoNUhMJx8XM/QuK4iEiciHwlIms9+/I7L8s0bf5S1ZB/4e7eeTvQB4gB1gJnnrbMrcCzns+TgDeDHXcD9+PnwN+DHauf+/M9IB3YUMv3lwH/xj0q22jgy2DH3MD9GAssDHacfu5LZyDd8zkR2Orlb6zFHxc/9yMkjovn37m153M08CUw+rRlmjR/hcsZvz8DuF8OzPJ8ngdcKCLehoEMpoAPRN+cVPVz4IiPRS4HXlG3FUCSiHRunuj858d+hAxVzVbVLM/nE8Bmao553eKPi5/7ERI8/84nPZPRntfprW6aNH+FS+L3ZwD3U8uoajlwHOjQLNH5z6+B6IErPZfg80Sku5fvQ4W/+xsKxngu1f8tIoOCHYw/PNUFw3GfYVYVUsfFx35AiBwXEXGKyBogF1iiqrUek6bIX+GS+P0ZwN2vQd6DzJ8Y3wN6qepQYCnfngWEolA4Jv7Iwt0nyjDgb8C7QY6nTiLSGngbuF1V80//2kuRFnlc6tiPkDkuqlqhqmm4xyA/S0QGn7ZIkx6TcEn8/gzgfmoZEYkC2tLyLt/r3A9VPayqJZ7J54ERzRRbIPhz3Fo8Vc2vvFRX9+hy0SKSHOSwaiUi0biT5Wuq+i8vi4TEcalrP0LtuACo6jHgU2DcaV81af4Kl8R/agB3EYnBffNjwWnLLABu8Hy+CvhYPXdKWpA69+O0utYJuOs2Q9UC4HpPK5LRwHFVzQ52UPUlIqmV9a0ichbu/1eHgxuVd544XwQ2q+qTtSzW4o+LP/sRKsdFRFJEJMnzOR64CNhy2mJNmr+CMuZuU9NaBnAXkYeAVaq6APcfyWwR2Yb7l3JS8CL2zs/9mCoiE4By3Pvx86AFXAcRmYO7ZUWyiOwDHsR94wpVfRb3uMuXAduAQmBycCL1zY/9uAq4RUTKgSJgUgs8qah0DvAzYL2nThngHqAHhNRx8Wc/QuW4dAZmiYgT94/TXFVdGMj8ZV02GGNMhAmXqh5jjDF+ssRvjDERxhK/McZEGEv8xhgTYSzxG2NMhAmL5pzGNBUR6QB85JlMBSqAPM90oap+NyiBGdOErDmnMbUQkQzgpKr+OdixGNOUrKrHGD+JyEnP+1gR+UxE5orIVhF5TESu8/Spvl5E+nqWSxGRt0Vkped1TnD3wBg3S/zGNMwwYBowBPcTpANU9SzgBeD/eZaZAfxFVUcBV3q+MyborI7fmIZZWdl/jYhsBz70zF8PnO/5fBFwZpVu09uISKKn/3hjgsYSvzENU1Lls6vKtItv/185gDGqWtScgRlTF6vqMSZwPgR+VTkhImlBjMWYUyzxGxM4U4GRntHSNgG/DHZAxoA15zTGmIhjZ/zGGBNhLPEbY0yEscRvjDERxhK/McZEGEv8xhgTYSzxG2NMhLHEb4wxEeb/A/+U6uJPiWMyAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret2 = ensemble_simulations(duration, ndiv=20, y0=y0, model=m, solver=('meso', Integer3(4, 4, 4)), repeat=N)\n", "ret2.plot('o', ret1, '-', **opt_kwargs)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Make a model with intrinsic rates. This model is for microscopic (particle) simulation algorithms." ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "with species_attributes():\n", " A | B | C | {'radius': radius, 'D': D}\n", "\n", "with reaction_rules():\n", " A + B == C | (ka, kd)\n", "\n", "m = get_model()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Simulating with `spatiocyte`. `voxel_radius` is given as `radius`:" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEGCAYAAABiq/5QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXxU5b348c93JnsIBEggrLKj7ARQcKt7qVTUupTWVuv1Cq3tT3rby9VSl9Ral3qrF1ur4oobSrUqGiuCaxGXQGQHkX2RkLBD9sx8f3/MJCSQ5SQzk8nMfN9tXjPnzHnO8z2Ofs+Z5zzneURVMcYYEztc4Q7AGGNM67LEb4wxMcYSvzHGxBhL/MYYE2Ms8RtjTIyJC3cATmRkZGifPn3CHYYxxkSUZcuW7VXVzOPXR0Ti79OnD0uXLg13GMYYE1FEZFt9662pxxhjYowlfmOMiTGW+I0xJsZERBu/McaEQ0VFBZs2baKkpCTcoTQqJSWF/v37k5CQ4Gh7S/zGGNOATZs2kZ6ezuDBg3G52mYDidfrZc+ePWzcuJEhQ4Y4KhPSIxGRdBF5VUTWi8g6EZkgIp1EZKGIfON/7RiKuvPmP05BzgC8d3agIGcAefMfD0U1xpgoVlJSQteuXdts0gdwuVx07dqVkpISCgoKnJUJcUyzgHdV9WRgJLAOuBV4X1UHAu/7l4Mqb/7jDFt2G1kU4RLIoohhy26z5G+Maba2nPSruVwuRIQ333wTj8fT9PahCkRE2gNnA08BqGqFqh4ELgXm+DebA1wW7Lp75T9AslTUWZcsFfTKfyDYVRljTB0/fPwzfvj4Z2Gpu6ysjNLS0ia3C+WprB9QBDwjIl+JyJMikgp0VdXdAP7XLvUVFpGpIrJURJYWFRU1q+IuWv/2XXRvs/ZjjDFtweuvv46IsH79+ka3ExFH+wtl4o8DsoFHVXU0UEwzmnVUdbaqjlXVsZmZJzxx3KhCqX/7Qslo1n6MMaY53vhqF19tP8gXW/Zzxn0f8MZXu4Ky37lz53LmmWfy8ssvB2V/oUz8O4GdqvqFf/lVfCeCPSLSDcD/Whjsindkz6BU63ZrKtUEdmTPCHZVxhgD+JL+7/65igqPF4BdB0v53T9XBZz8jx49yqeffspTTz3V9hO/qhYAO0RksH/V+cBaYD5wnX/ddcCbwa573ORprB5zNwVk4lUo8Hbki4H/xbjJ04JdlTHGAPDAgq8prax7Y7W00sMDC74OaL9vvPEGEydOZNCgQXTq1In8/PyA9geh79Xz/4AXRWQlMAq4B7gPuFBEvgEu9C8H3bjJ08jK2cjeqcvJlIMkljjr5mSMMS3x7cH6b6o2tN6puXPnMmXKFACmTJnC3LlzA9ofhPgBLlVdDoyt56PzQ1lvbV169GV56ngGfvsmlRUPEJ+Q2FpVG2NiSPf0ZHbVk+S7pye3eJ/79u3jgw8+YPXq1YgIHo8HEeHPf/6z4xu59Wn7HVSDYfS1ZHCQ1R/9I9yRGGOi1IzvDiY53l1nXXK8mxnfHdxAiaa9+uqrXHvttWzbto2tW7eyY8cO+vbty+LFiwOKNSYS/7BzrqSQTrjynw13KMaYKHXZ6B7c+4PhJLh9abVHejL3/mA4l43u0eJ9zp07l8svv7zOuiuuuIKXXnopoFhjYqyeuPgENvX6Aadtf4odG1fRa8DwcIdkjIlCl43uwdwvtwPwyrQJAe/vo48+OmHdzTffHPB+Y+KKH2DgpOlU4eLbBf8X7lCMMVHslWkTgpL0QylmEn9GVm9WpJ/PsMK3OXxwX7jDMcaYsImZxA+Qfu7NpEoZa9/5e7hDMcaYsImpxD9w1Fmsix9Cr2+ex1NVFe5wjDEmLGIq8QOUZk+lh+5h5QfBefTZGGMiTcwl/hEXXEMBmSQstbH5jTEh8Mwk318bFnOJPy4+ga0Df8rQipV8vfSDcIdjjDGNcrvdjBo1ipEjR5Kdnc2SJUsC3mfMJX6A4ZOnc5B2lH7w53CHYoyJJivnwc482LYYHhrmWw5QcnIyy5cvZ8WKFdx777387ne/C3ifMZn4U9PSWdf7GkaVfMbm1V80XcAYY5qych68dTN4yn3Lh3b4loOQ/KsdPnyYjh0Dn6Y8Jp7crc+Qy2ZQPGsOBxbcB8NeD3c4xphI9/5dUHncIG2Vpb71I65u8W5LS0sZNWoUZWVl7N69mw8+CLyJOiav+AE6dMpkZfcrGXX4Q3ZsXBXucIwxke7Qzuatd6i6qWf9+vW8++67XHvttahqQPuM2cQPMPDSW6kkjt1v3xPuUIwxka5Dz+atb4EJEyawd+9emjsP+fFiOvFnZPViReZkRh9YwK7Na8IdjjEmkp1/B8QfN/Z+fLJvfZCsX78ej8dD586dA9pPTCd+gP4/uJMq3Ox+I3hfjjEmBo24Gi55GNz+yZ469PItB9C+D8fa+EeNGsUPf/hD5syZg9vtbrpgI2L25m61jO4n8Vn3KZz27fNsWfMFfYeeFu6QjDGRasTVsGyO7/31uUHZpcfjaXqjZor5K36AIVfdwVFJ5lDuneEOxRgT6a7PDVrSDxVL/Ph6+Kztcz2jSj5jfd4iANbccyZr7jkzzJEZY0zwWeL3G3HlLewlHc/CHNTrDXc4xhgTMpb4/VLadWDTKTcxtGIVyxcFNp+lMca0ZZb4axnzg/9iq6s3hxY/TqfyXZxSvoqCnAHkzbeRPI0x0cMSfy1x8Qms73Ixp8l6usl+XAJZFDFs2W2W/I0xjlz/7vVc/+714Q6jUSFN/CKyVURWichyEVnqX9dJRBaKyDf+18BHHAqiUQWvkSwVddYlSwW98h8IU0TGmFhWUFDAlClT6N+/P0OGDOHiiy9mw4YNAe2zNa74z1XVUao61r98K/C+qg4E3vcvtxldtP5Hobvo3laOxBgTaXI357KyaCVL9yzlolcvIndzYN06VZXLL7+cc845h02bNrF27Vruuece9uzZE9B+w/EA16XAOf73c4CPgFvCEEe9CiWTLE5M/oWSQVYY4jHGRIbczbnkLMmhwutrMdhdvJucJTkATOrXshm5PvzwQ+Lj4/n5z39es27UqFEBxxrqK34F3hORZSIy1b+uq6ruBvC/dqmvoIhMFZGlIrI00AGJmmNH9gxKNaHOulJNYEf2jFaLwRgTeWblz6LMU1ZnXZmnjFn5s1q8z9WrVzNmzJhAQztBqK/4z1DVb0WkC7BQRNY7Laiqs4HZAGPHjg1sDNJmGDd5GnlAz2X30ZX9FGpH8lK/wyWTp7VWCMaYCFRQXNCs9eEU0it+Vf3W/1oIvA6cCuwRkW4A/tfCUMbQEuMmT2N/Yg/WJQ5nV4dRXFSSy45vVoQ7LGNMG5aVWn9jcEPrnRg6dCjLli1rcfmGhCzxi0iqiKRVvwcuAlYD84Hr/JtdB7wZqhiC4aRrHqZcEjj8j1/hDcFgScaY6DA9ezpJ7qQ665LcSUzPnt7ifZ533nmUl5fzxBNP1KzLy8vj448/bvE+IbRX/F2BxSKyAvgSyFXVd4H7gAtF5BvgQv9ymzN05mKGzlxMRlZv1g+fwdCKleS9al06jTH1m9RvEjmn55Dg8t0j7JbajZzTc1p8YxdARHj99ddZuHAh/fv3Z+jQoeTk5NC9e/eAYm2yjV9ErgLeVdUjInIbkA3crar5jZVT1c3AyHrW7wPOb2G8YTHu8ums+CaXEWv/wo5vvkevgSccljHGMKnfJF7d8CoAz0x8Jij77N69O/PmBW/CdnB2xX+7P+mfCXwXXxfMR4MaRRsnLhc9rn2Kckmg9JX/pKqyoulCxpiY9MzEZ4KW9EPFSeKvbtieBDyqqm8CCY1sH5Uyup/EplPvYlDVBvJeuD3c4RhjTIs5Sfy7RORx4GrgHRFJdFgu6oy5+AaWpp3P2K1PsH7p+w1uZ2P5GxM9vBEwTHtzY3SSwK8GFgATVfUg0AmI2aeZBv7HbIpcnUl/exqH9gX22LQxpm1LSUmhoKCgTSd/r9dLQUEBlZWVjss0eXNXVUtEpBA4E/gGqPK/xqQOHTPYc8kT9HnzB6x96jpG/vc7iCsmfwAZE/X69+/P8uXL+fbbbxGRcIfToMrKSrZv346q4nKQj5z06rkTGAsMBp4B4oEXgDMCjDViDco+h8/X/4bxGx7g87l/ZPw1NlevMdEoISGB5ORkFi5cSIcOHRwl1XA5cuQIXbt2JTk5ucltnRzF5cBkoBhqnsZNCyjCKHDalJl8lXIGYzbMYu1n/wp3OMaYEBk2bBjnnHMOSUlJiEib/HO5XAwcOJDJkyc7+mXiZKyeClVVEVGoeQo35onLRb8bn6Pg4bPIWjCV3d0X0e2kweEOyxgTZCJCdnY22dnZ4Q4laJxc8c/z9+pJF5EbgUXAE02UiQkdOmagP5pLvFZS8twUSo4eCndIxhjTpCYTv6r+L/Aq8Bq+dv47VPWvoQ4sUvQeNIrN33mYvlVbWP/YT9E2fPffGGMARLXVRjxusbFjx+rSpUvDHUajPn/udsZvfpi3ki5hbOmnviGdJZMd2TMYZ0M6G2PCQESW1Zr9sEaDbfwicgTfRCrif635CFBVbR/0KCPYaT/5A7n3rOOC0gU1c/ZmUUSHZbeRB5b8jTFtRoNNPaqapqrta722r73cmkFGAnG5yK5abhO1G2PaPCf9+HvXt15Vtwc/nMjWVYt8v4eOYxO1G2PaEifdOWtPE58E9AW+BoaGJKIIFqyJ2qvH+Rk6c3GQIjPGmGOc9OoZXutvIL7pEy0j1aOhidpXn/TTMEVkjDEnavbzx/4JWMaFIJaIN27yNFaPuZvd2gmvwm4yyPMOInvbk2xZmxfu8IwxBnDWxv+bWosufDNwndieYQBf8l+z+nn204OhMxdTtXEVVS9MJn3e5Wy8/GUGjLThmo0x4eXkij+t1l8ivjb/S0MZVDTpNWA4lde+TRlJdHn9KtbnLQp3SMaYGOdkWOY/tEYg0axHv6EU/Me/OPzMJfR++8esKn2C4WfbudMYEx5NXvGLyEIRSa+13FFEFoQ2rOiT1XsgSTcuYI+7G4Pfv56l8x8Ld0jGmBjlpKkn0z/zFgCqegDoErqQoldG95PofPMHfJM4lLH5t/D5c7fb2D7GmFbnaLL12g9xichJ1B3CwTRD+/TODPjNApalncv4zQ+T97drqSgvC3dYxpgY0uQgbSIyEZgNfOxfdTYwVVVbrbknEgZpay6vx8MXT/0XE76dw9r4YWTdOI9OXXqQN/9xei67zwZ5M8YErKFB2hyNzikiGcB4fAMSfKbqfAwCEXEDS4Fdqvp9EekLvIxv0vZ84KeqWtHYPqIx8Vdb+vZshuXN5ICks+akazhj69/rjPdTqgmsHnO34+RvT/0aY6o1lPid3NwVYCKQrapvASkicmoz6p4OrKu1fD/wkP8p4APADc3YV9QZ+/2p7Ljsn7jxMHTLHBvkzRgTck7a+P8OTAB+5F8+AjziZOci0hOYBDzpXxbgPHwTuwDMAS5rRrxRaeDos4m7aTFd5UC9n9sgb8aYYHKS+E9T1V8CZVDTqyeh8SI1/g/4H6C660pn4KCqVvmXdwI96isoIlNFZKmILC0qiv4HhTt16cEeyaz3s0LJaOVojDHRzEnir/S301dPtp7JsUTeIBH5PlCoqstqr65n03pvMqjqbFUdq6pjMzPrT4jRZmc9g7yVaALbR/02TBEZY6KRk8T/MPA60EVE/oRvZM57HJQ7A5gsIlvx3cw9D98vgHQRqX5iuCfwbXODjlYnDPKmnVjjPYkOq59l2/r8cIdnjIkSToZlfhFfc829wG7gMlX9h4Nyv1PVnqraB5gCfKCq1wAfAlf6N7sOeLOFsUelcZOnsT+xB+sSh5N15ybk1BvpWrWTbnMv5LM5M6msKA93iMaYCNdg4heRTtV/QCEwF3gJ2ONf11K3AL8RkY342vyfCmBfUU1cLsZeMo2qn3/O6rQzmLDlEXbcdyrrv1wY7tCMMRGssUHalnFssvXjKdDPaSWq+hHwkf/9ZnyTuRiHMrJ6kfHf88lf8Dw9PruTfu9cyZefTWLgj/9Cx8xu4Q7PGBNhGptsva+q9vO/Hv/nOOmb4Mn+7k9p99t8Ps+6htH738X9SDafv/gHG/LBGNMsTubcRUQm4xuqAeAjVX07dCGZxqSmpTP+539n27obOPjmLYz/5kF23jeXwtNuxVNeSs/yXXRlPwU5A5o93IM99WtMbHAyA9d9+KZafNG/arqInKGqvwtpZKZRJ50yhpNOWcTKj16j/Sd34vn0bwxzbat58jeLIjosu408sLF+jDF1OOnOeTFwoao+rapP4xu+YVJowzJOjTjnCnrNXE5P134b7sEY44jTydbTa73vEIpATMu54+Loyr56P+uiRa0y5v+ae86saSoyxrRtTtr47wW+EpEP8fXwORuwZp4Qakkbe6FkksWJQ1sUakdK7x7N3hFTGfm9G0hITApGiMaYCObkAa65+IZk/qf/b4KqvhzqwEzz7KhnuIdSTWBZp+8DyrjlMzl078l89swt7C3YEZ4gjTFtQoNX/CKSfdyqnf7X7iLSXVVtDIE2ZNzkaeRB3Ulcxsxg0uRpqNfLio9fQ754nAnbHqPi0SdY2uFcksffwJDxExGX0xY/Y0xruf7d6wF4ZuIzQd93Y009S4E1UNN+UPtBLsU39o5pQ8ZNnsaa1c+znx4MnbmYLP96cbkYee5VcO5VbN+wnN0L/8opRe/Q/r1FbF/Ug2/7XkllfHsGBNAVNFDWldSY1tNY4v8tcAVQim+QtddV9WirRGVCpvegUfQe9BSlxUfIWziHtNUv4N7wDiOtK6gxQRfKq/ZANPbk7kOqeibwK6AX8L6IzBORUa0WnQmZ5NQ0xl32K06+7XN6ug/U2xW0Z/79lJUWhylCY9qG69+9viaBRwsnN3e34BtB8z18Y+wMCnVQpnV1bWCGr666D899/Vj2l8tY9s5THDm0v97t8uY/TqfyXZxSvoqCnAHkzX88lOGeIJCupNYNNTZEY/IORGM3d/vhG075UmAHvuaeP6mqDQwTZRrqCrpHOrOj0+kM2P8xnb78kIovZrAqaSTFfS6gx7jJ9BownLz5jzNs2W3WTGRMBGmsjX8jsBLf1f5hoDdwk2/aXFDVB0MenWkVO7Jn0KFW8gZfV9CdY27h1MnTqKqsYO2yDzi8fD7dCz9i+Nd/hq//zE7Joqe3ouEnhqM88cfiDelA2qwDbe9uq+3lkaixxH8Xx6ZFbNcKsZgwaagraPUVe1x8AkPGT4TxEwHYuXE1u5a9TdLWDxhe8kW9+7QJ4o1puxpM/Kqa04pxmDBrqCtofXoOGEbPAcOAW9mdM4Bu9TUTaTolfxzN3o6jcfedQPehZ9Ot98B6nxnIm/94QKOKmuaxK2fjaFhmYxqyM3sG6fU0Ey1NPZu+3m0ML3qblL2vQR7sowM7UoZQmjmSlD5j6TX0dDZ9/rbdI4ggxdu/CncIrS4aj9kSf5Rp7fbmhpqJLvEn7arKCjau/ZJ9Xy/BtWsZXY+soufWz3Fteww+hnJvR5JdsXmPIBB21R79QnnCaaxXz3RVneUfe//TkEVgIl5jzURx8QkMGHkmA0Ye6zJ55NB+tq/5nCObv+TUbx6qd59dtIilD15JZedBJHUbQuc+w+jW5xTiExLrbBepzURXz/Y9DjNv6vIWlQ8kKUTjFWxTwnXMbfWfdWNX/NcDs4C/AseP22NMi6V16MTQ0y+G0y+mIOeF+ruS0oleh/PpenghbAGWQKW62eHqyv6kXpSm9WFvGZx/6PUWNxNF6kkjVrXVJNqQ3M25fJPgoQq46NWLmJ49nUn9nE1lEkhZJxpL/OtEZCuQKSIra60XQFV1RNCiMDGr4a6ktzJu8jSOHNrP7k2rOLRjDVWFX5N4aDPpJdsZVPIVhzS13mainsvu4/PCr3F17E1y5km079qXzB59SWl3bCqJQJ8/sJNG8+RuzmV5wdKQJbK2JndzLjlLcqjyj3C2u3g3OUtyAJo87kDKOtVYr54fiUgWsACYHJTaTJsWjv7oTXUlTevQibTs70D2d+qUU6+XxD90rHefXdlPlx1P496pddYfJpX9rs4cTsikZ+mWBoap+HOT9xbsobXmCUYiC9fVc0vLzsqfRZmn7rOuZZ4yZuXParJ8IGWdavTmrqoWACNFJIFjQzV8raqVQandGJrXlbSauFzsaWjyGcmk88w17Nm1hYMFWygp2kblgR24jnxLQske2lUU0pX6h5/oqnspubMLB13pHHGnUxrfkYrETniSOyEpnXG1y6DXiv+L2YfWWiLQRBauq+dAyhYUFzRrfbDKOuVksvXvAM8BW/E18/QSketU9ZOgRWFMCzTUTLRjzAyyEhLp3vdkuvc9ud6yBTkDGri30JltWRcRV1pEYvl+0ioKSSvdQPqBQySIBwCvUneQcr8uWsTmu0ZSFpdGeVwaVfHt8Man4U3qgCSmIUkdcKe0Jz65A72PdKDc5WXX5nUkpaaRmpZOYlKKo7kRIu3qN9BEFq6r50DKZqVmsbt4d73rmxJIWaecdOd8ELhIVb8GEJFBwFxgTGOFRCQJ+ARI9NfzqqreKSJ98Y370wnIB36qqhUN78mY+lU3E3Vfdi/dOHBCM1FjGhumYnw95dXr5ejRQxzeX4g8ezHdOPHJ5D104lBSdxKqjtCh7FuSS4pJpYR2WoJLjjU7PZzcl686V1IUJ6z94Eou3RfPzaVb8KhQQhKlkky5JFHuSqbSlUSlOxmPOwmPO5lPk0t5LnlnnavQO/79ezbmLeCMdsNxxyfjTkgmLjGFuMRk4hKSiU9MIj4xhY/2fsb/rn241a9+A01k4bp6DqTs9Ozp5CzJqXPiSHInMT17ekjLOuUk8cdXJ30AVd0gIvEOypUD56nqUf/2i0XkX8BvgIdU9WUReQy4AXi0JcGbtiNc49WMmzyNqwseBdKZN3W5o2ai6nLNOWmIy0W79h1p174jedn/U+9Da9U3pI/n9Xg4WnyYkiMHeXbR3fzD+zll/iv7wnjh+S6VHDzyHS5M7odUFuOqLMZdVYq7qoR4TwlJVYeJrygiUct4q308FVL3V0EFHnKPLGT6ujmNHvPTPbtTHl/3P/kyTxkPfTCD0567iSriqZR4POLGI/F4JM736ornf7scpCzOe0LZP398OxnvvIq64sAVj7rciCsedcUhrjhwubmYXjxPIRV4asomEMclOoyluU/gcsWBKw6XOw5xuf1/cYjbhcsVR+e4dPZWHTjheDISOrHt6+W43G5cLhfiivO/unC74sDloktSJnvKCk8o2zWlK+VlJbhcbkSk5rX2L65ATljVJ8PbPrmVKqBbu26OfyEFUtYpUdXGNxB5Gt+YPc/7V10DxKmq4zFORSQFWAz8AsgFslS1SkQmADmq+t3Gyo8dO1aXLl3qtDoTgQLp1x6OsnnzH6970nDYq+f8J4dSGH9ic06XSi/v/+eaJssPnzO8wc/ePm0OleUlVJWXUlVRiqeiBG9lOZ6KMrwVpdxw5Kl6m6hQePLwKFyeCsRbictbictbgUurcGkVbm8lV/c40mDZT7YeIE6riMODG09Nk1htuakpzOqYTkGcm6wqD9MPHGRScQkKeACPgAfBA3il+tW3bmFKMv/XKZ3yWkk50etl2oFDnFlWhhfwIngF//6OvV+SlMTzHdpT4ToWfLxXmXLkCGPKyv1lQUVQ/3tfHMLyxATeap9KpRwrG6fKxCNlnFzh8W8veGvKCl58731/ghfFi4C4/PX41tfeRmvVT63PPOoBhCsnPsPQk09r8HtvjIgsU9Wxx693csX/C+CXwM34vvpPgL87rNQNLAMGAI8Am4CDqlrl32Qn0KOBslOBqQC9e/d2Up0xraalvzSK4urLng2vP16XSm+DJ45eg0dR7imnwlNBuae8zvsKTwVd5j9Rf9kqL1+ffyEV3goqvZVUeiqp8lb53nt977us/UeDZaefdh5V3irfn1ZR5fG992gVVV4PHq8Hj3o4Un4QAfYnpnBb10RmqhcvjV94Vht2qB2FKYcpihMyq5QuJe15uLOLhx2VrqvSJTzfoT3Pd2h62/rqfbuD8HYL6m0OX8tgHAJMcgW/JbzJxK+q5fja+Zs9DLOqeoBRIpIOvA6cUt9mDZSdDcwG3xV/c+s2kSW19+hwh9AsLb3RmVmlFMafmOQzq5S3Nr1FaVUpJZUlvtcq32vtv1/v389dmZ1qmooAkrxefr1/PyOfG9lo3fc0Unbm0gdq1rnERbwrnjhXHAmuBOJd8Y3W+2avBFLiUnC73MRJHHGuuJr3bpcbt7jJWrOTS3YuohsHKJDO/KvXREqyTyFO4nCJq2Y7l7jqvrpc6Bd5fH/fM6TsP5YASzSBBZ2n0v6s83DhK+8SF4LUlBXxvS/86F+MWPsk3TjAHslg5ZAb6H3+5bj8TWYucdWUrS4jCJsWzmP85vtJ3l+3SS+v50yGfO+nADX1VZcVEWr+J8J1z5wOwAv/8UWdz3z/P7Zd9Wtt1b9Gxw46q9HvtSVaZaweVT0oIh8B44F0EYnzX/X3BL5tjRhMaEXi2DEtSd6qyuvfvM49X95T50bn7Z/ezkfbP6J7WneOVBw59ld5hKMVR31/lUfpV9Kew2mHT0igXUraM3PxzDp1JbmTSI5LJjkumaS4JJLiksguTiJH9p/QbDK6JJlfjPwFCe4EktxJJMYlkuBKIDEukURXIgnuBPo9cS05e08sm12SwpIfLSHeFU+8Kx63y33CcRfkDKi37JiSVC757pON/jPLm/84w3Y+UnNPpDv7+PH2f7A6825HzWMFL/2OlOO6z6ZIBRO+eYWsa/7UZN3nrJ1VU3c39pK+5iFWJ3Vusu6EVY/V22130MpHyPjBfzUZd978x5m165CvOfDuYc16yC9v/uM8tOug70QZggcEQ5b4RSQTqPQn/WTgAuB+4EPgSnw9e67DN9GLMS3S0ivv+nqp3PHpHazZu4Z+6f3YV7qP/WX7OVB2gP3lvteDZQc5WH6QCu+JP70rvZW8u+1d4lxxtE9oT/uE9rSLb0e7hHZ0Tela8z41PpU+yz/hy8rVNc0H5yaP5wfX/JbU+FRS4lJIiU8hyRTahW8AABZ0SURBVJ1UbwLO8z7OectuY1LxseulUk1g9ZhbuGlUEw+eZf9PA2X/h3EJaY2W3ZE9o4GyM5ps5uqV/0BAzz100aIGus82PedDIHUHUm8gD/m1xgOCjSZ+fxv9fao6owX77gbM8e/DBcxT1bdFZC3wsojcDXwFPNWCfRvTZBdDVWV/2X4KSgooKC5gT/EeCksKKSwpZMHWBSck8ApvBc+ve75mOS0+jY5JHemY1JHu7boztPNQ0pPSeWZ1/b9qBCH/J/kn/GQ/wahfcvXsUXSuaP5N5UC6sIarbCAJFBqeGrRQMpo86QRSdyD1BnLCCfRE6URTT+56RGSMiIg21f3nxLIrgRMablV1M75J240JSEMP2Nzx6R08uuJRdh/dfUJyj3PFkZmcWe9VO/iS93tXvkenpE4kuBPq3ebdLe822M2vyaQfBC29sRyusoEkUGjiQb0Q1h1IvYGccAI9UTrhpKnnK+BNEfkHUFy9UlX/GbQoTJsQrnb63M25rCxaSYW3ot7mmqMVR9l6eCtbDm1h6+GtbDu8je2Ht9ebfMF35T6442DO7XUuWalZx/5SsuiY1BGXuLjo1YsaTN5N9dNujQdsokkgCRRC86Cek7oDqTeQE06gJ0onnCT+TsA+4Lxa6xSwxG8CVt1cU30FXn2j9K1Nb+FRD5sPbqaw9NgDOG5x06NdD3q3701KXAolVSUn7LNbajf+cs5fGq03kORdfVK649M7qPBW0C01+A/YRJOmBuJzuo/WeFAvWPUGcsIJ9ETphJPunI4f1DLGqb2le1m7by1/+vxPJzTXVHorWfLtEoZ0HsL47uPp26Ev/Tr0o0+HPvRq14t4t+/B8eqTRjiS96R+k3h1w6tAZPVkCpeWDMR3vDv3tmtx3S1t3mqpcN1PccrJIG2D8A2p0FVVh4nICGCyqt4dtChMVCuuLGbN3jWs3LuSVUWrWLNvDXtK9jRZ7uXvv9zo55GcvFs681a1libBQMsa56pPdgfp2eyTXahPVk6aep4AZgCPg++mrYi8BFjib4PC0U5/fBv9z4b+jI5JHcnfk8/youVsOLABr/rGeTmp/UmM6TqGoZ2HMjRjKLd+cisFJScOeuV0AK9YvfIO19hIkSoST3ahjNlJ4k9R1S+P661Q1dDGJrbkbs7lziV31mmjv/fLewFIjktmZOZIpo6YysjMkQzPGE6HxLrPyf96zK/tRmkEicQEGqhoPGYniX+viPTHP7SCiFwJ1N+dwsSECk8F+YX5LN65mBfXv0iV98TrgIzkDBZeuZA4V+P/itmN0tZnvxaMk8T/S3xj5pwsIrvwTX19TUijMq2uqS6Ve0v38u+d/+bDHR/y+e7PKa0qJd4VX2/SB9hXuq/JpF8tVptrjAkXJ716NgMXiEgq4FLVI6EPy7Sm+rpU5izJ4UDZATzqYdG2RawoWoGiZKVmMbn/ZM7qcRbjssZx2ZuXhXy2oLYqFk9SgfxasF8abYeTXj2dgTuBMwEVkcXAXaq6L9TBmdbR0BOw9+fdD8DgjoP5xahfcG6vcxnccXCdp1PtYSbTWuzEETxOfou/jG8M/iv8y9cAr+AbdM1Egcamksu9PJfe7RueD6EttNEHcuUdi1ftxjh6cldV/1hr+W4RuSxUAZnWoaqs2ruqpm29Pt1SuzWa9KtZG71p66yJqi4nif9DEZkCzPMvX4lv+kQTgUoqS3hr01vM2zCPDQc2kByXzJiuY2pu7Faz5hpjAtdWTxoNJn4ROYKvC6fgmyD9Bf9HLuAovnZ/EyF2HN7BS+tf4o2Nb3C08iindDqFOybcwcV9LyY1PpXczbkR21xjjGmeBhO/qjY+M4Npc+rrktkrrRfPrnmWRdsW4Xa5ueiki/jxKT9mRMaIOjdprbnGmLYllL8WHHW09o/P06f29jYsc9tSX5fMmf+eiRcvafFp3DD8Bn588o/JTMkMc6TGmHBz0p3zaWAEsAbw+lfbsMwh1JLxdurrkunFS/uE9rx35XukxqcGNUZjTORycsU/XlWHhDwSE5CGJiU5UnGkVZK+NQ8ZEzlcDrb5TEQs8bdR2w5vY/oHDfe+iYUnaI0xzePkin8OvuRfAJTj6+WjqjoipJGZRh0qP8RjKx7j5fUvk+BO4MLeF/LvXf+2J2iNMU1ykvifBn4KrOJYG78JE696eXPjmzy07CEOVRzi8gGX86vRvyIjOSPgLpnWXGNMbHCS+Ler6vyQR2KatOHABu7+/G6+KvyK0V1G8/vTfs/gToNrPrcumcYYJ5wk/vX+GbfewtfUA1h3ztZU4alg9srZPLXqKdIS0vjjGX9kcv/JuMTJLRpjjKnLSeJPxpfwL6q1zrpztpLVe1dz+6e3s/HgRib3n8yMsTNIT0oPd1jGmAjmZDz+61uyYxHpBTwHZOG7NzBbVWeJSCd8o3v2AbYCV6vqgZbUEc1Ulb8v/zuzV84mIzmDR85/hLN7nh3usIwxUcDJA1zP4J92sTZV/Y8milYBv1XVfBFJA5aJyELgZ8D7qnqfiNwK3Arc0uzIo1Tu5lxWFK2g0lvJssJljM4czSMXPEJago2gYYwJDieNxG/jG40zF3gfaI9vkLZGqepuVc33vz8CrAN6AJfi6yKK/9WGePbL3ZzL7Z/eTqW3smbduv3r+GTnJ2GMyhgTbZw09bxWe1lE5gKLmlOJiPQBRgNfAF1Vdbd/37tFpEsDZaYCUwF69256TPhIV+Gp4K7P7qqT9ME3E9as/Fk2+bgxJmha0i1kIOA4E4tIO+A14NeqethpOVWdrapjVXVsZmZ0DyxWUFzAz979GSVVJQ1+bowxweKkjb96XP5qBThskxeReHxJ/8Va3T/3iEg3/9V+N6CwmTFHlRVFK/j1h7+mpLKEjokdOVB+4n1uG3bBGBNMTpp6WnRXUXyDvT8FrFPVB2t9NB+4DrjP//pmS/YfDeZvmk/Okhy6pnTliQuf4OsDXwc8cbk9uGWMaUpjM3A12pyjqtub2PcZ+Id6EJHl/nUz8SX8eSJyA7AduMp5uJGhqWGVVZW/r/g7j614jFOzTuUv3/kL6UnpDOg4AAjvxOXGmOjX2BV/LsemXqymQCbQBXA3tmNVXXxc2drOb0aMUaXSW8kfP/sjr298ncsGXMYdE+4g3hVf87kNu2CMCbXGpl4cXnvZ3zPnFuAC4J6QRhWlSqtK+e1Hv+Xfu/7NtBHT+OWoX9aZ/tAYY1qDk5u7A4HfA6cBfwFuVtXKxkuZ4xVXFnPToptYXrSc28ffztWDrw53SMaYGNVYG/8wfAl/KPBn4AZV9bRWYNHkUPkhblp0E2v2reH+s+5nYt+J4Q7JGBPDGrviXwHswNfWfypwau1mCVW9ObShRYeDZQe5ceGNbDq4iQfPeZDzep8X7pCMMTGuscTf1Fg8pgmHKw4zdeFUNh/czMPnPcyZPc4Md0jGGNPozd05DX1mmubxerhp0U18c/AbZp07y5K+MabNcDIev2mG2qNrAlxzyjU2nLIxpk2xxB9EuZtzyVmSU2egtdc2vMbwjOE2960xps2wufuCaFb+rDrDLcCx0TWNMaataDLxi8ggEXlfRFb7l0eIyG2hDy3y7C7eXe96G13TGNOWOLnifwL4HVAJoKorgSmhDCoSLdi6oMHPbHRNY0xb4iTxp6jql8etqwpFMJFq3b513Lb4Nnqn9SbJnVTns+aOrmmMMaHmJPHvFZH++MfkF5ErgfrbNGLQvtJ93PzhzXRI7MCc780h5/QcElwJAHRL7UbO6Tk2uqYxpk1x0qvnl8Bs4GQR2QVsAX4S0qgiRKWnkt989BsOlB1gzvfmkJGcYaNrGmPaPCcTsWwGLhCRVMDlnzjdAA8ue5D8wnzuO+s+hnYeGu5wjDHGESejcyYCVwB9gLjq8XpU9a6QRtbGvb/tfV5Y9wLXnHKNNeUYYyKKk6aeN4FDwDKgPLThRIadR3Zy+6e3M6zzMH4z5jfhDscYY5rFSeLvqao2jrBfpaeSGR/PAOCB7zxAgjshzBEZY0zzOOnVs0REhje9WWz461d/ZfW+1fzxjD/SM61nuMMxxphma2wiltWA17/N9SKyGV9TjwCqqiNaJ8S2Y9meZTy75lmuGnQV55/U8LTB1pvHGNOWNdbU0wMY1VqBtHXFlcX8fvHv6dGuB/899r/DHY4xxrRYY4l/i6pua7VI2rgH8h5gd/Funp34LCnxKeEOxxhjWqyxxN9FRBrssqKqD4Ygnjbpk52f8No3r3HDsBsY3WV0uMMxxpiANJb43UA7fG36Mef6d68H4JHzH+Huz+9mQPoAbhp1U5ijMsaYwDWW+HcH8pCWiDwNfB8oVNVh/nWdgFfwPQy2FbhaVQ+0tI7W8Lflf2N38W6e/97z1nXTGBMVGuvOGeiV/rPA8f3/bwXeV9WBwPv+5TaruLKYF9e9yA8H/5BRXew+tzEmOjSW+Bvur+iAqn4C7D9u9aVA9STuc4DLAqkjVHI357KyaCXr9q8D4OROJ4c5ImOMCZ4GE7+qHp+0g6Grqu7273830KWhDUVkqogsFZGlRUVFIQilftXz5lZ4KwDwqpf7v7yf3M25rRaDMcaEUpudc1dVZ6vqWFUdm5mZ2Wr12ry5xpho19qJf4+IdAPwvxa2cv1Namh+XJs31xgTLVo78c8HrvO/vw7fyJ9tSufkzvWut3lzjTHRImSJX0TmAp8Bg0Vkp4jcANwHXCgi3wAX+pfbDFUlLT7thPU2b64xJpo4GZa5RVT1Rw18FFBvoVBatH0RWw5v4bL+l/HOlneo8FbQLbUb07On22QrxpioEbLEH2mqvFXMyp/FgPQB5Jyew86jOwEbadMYE33abK+e1vbmxjfZdngbN4++GbfLHe5wjDEmZOyKHyj3lPPoikcZkTmCc3qdA9iVvjEmetkVP/Dy+pfZU7KH6aOnUz2ZvDHGRKuYT/xHK47y5KonmdBtAqd2OzXc4RhjTMjFfOJ/ft3zHCw/aN01jTExI6YTf3FlMS+sfYFze53L0Iyh4Q7HGGNaRUwn/nlfz+NwxWFuHH5juEMxxphWE7OJv6yqjDlr5jC+23iGZw4PdzjGGNNqYjbxv7HxDfaV7WPqiKnhDsUYY1pVTCb+Sm8lz6x+hlGZoxjbdWy4wzHGmFYVk4n/nc3v8G3xt9w44kbrt2+MiTkxl/hVlWfXPMugjoM4q8dZ4Q7HGGNaXcwl/i8LvmTjwY385JSf2NW+MSYmxVzif2HtC3RK6sTF/S4OdyjGGBMWMZX4tx/ezsc7P+aqQVeR6E4MdzjGGBMWMZX4X1r/Em6Xmx8O/mG4QzHGmLCJmcR/tOIob2x8g4l9JpKZkhnucIwxJmxiJvG/sfENiiuL+ckpPwl3KMYYE1YxkfhVlVe+foURmSNsMDZjTMyLicSfX5jP1sNbuWrQVeEOxRhjwi4mEv9rG16jXXw7LjrponCHYowxYRf1if9Q+SHe2/YeF/e9mJT4lHCHY4wxYRf1iT93cy7lnnKuGHRFuEMxxpg2ISyJX0QmisjXIrJRRG4NRR25m3O56NWLuPfLe4lzxbHl0JZQVGOMMREnrrUrFBE38AhwIbATyBOR+aq6Nlh15G7OJWdJDmWeMgCqvFXkLMkBYFK/ScGqxhhjIlI4rvhPBTaq6mZVrQBeBi4NZgWz8mfVJP1qZZ4yZuXPCmY1xhgTkcKR+HsAO2ot7/Svq0NEporIUhFZWlRU1KwKCooLmrXeGGNiSTgSf31jIesJK1Rnq+pYVR2bmdm8IRayUrOatd4YY2JJOBL/TqBXreWewLfBrGB69nSS3El11iW5k5iePT2Y1RhjTERq9Zu7QB4wUET6AruAKcCPg1lB9Q3cWfmzKCguICs1i+nZ0+3GrjHGEIbEr6pVIvIrYAHgBp5W1TXBrmdSv0mW6I0xph7huOJHVd8B3glH3cYYE+ui/sldY4wxdVniN8aYGGOJ3xhjYowlfmOMiTGiesKzU22OiBQB21pYPAPYG8RwwilajiVajgPsWNqqaDmWQI/jJFU94QnYiEj8gRCRpao6NtxxBEO0HEu0HAfYsbRV0XIsoToOa+oxxpgYY4nfGGNiTCwk/tnhDiCIouVYouU4wI6lrYqWYwnJcUR9G78xxpi6YuGK3xhjTC2W+I0xJsZETeJvagJ3EUkUkVf8n38hIn1aP8qmOTiOn4lIkYgs9//9ZzjidEJEnhaRQhFZ3cDnIiIP+491pYhkt3aMTjg4jnNE5FCt7+SO1o7RKRHpJSIfisg6EVkjIidMUhEJ34vD44iI70VEkkTkSxFZ4T+WP9SzTXDzl6pG/B++4Z03Af2ABGAFMOS4bW4CHvO/nwK8Eu64W3gcPwP+Fu5YHR7P2UA2sLqBzy8G/oVvVrbxwBfhjrmFx3EO8Ha443R4LN2AbP/7NGBDPf+OtfnvxeFxRMT34v/n3M7/Ph74Ahh/3DZBzV/RcsXvZAL3S4E5/vevAueLSH3TQIZTyCeib02q+gmwv5FNLgWeU5/PgXQR6dY60Tnn4DgihqruVtV8//sjwDpOnPO6zX8vDo8jIvj/OR/1L8b7/47vdRPU/BUtid/JBO4126hqFXAI6Nwq0TnnaCJ64Ar/T/BXRaRXPZ9HCqfHGwkm+H+q/0tEhoY7GCf8zQWj8V1h1hZR30sjxwER8r2IiFtElgOFwEJVbfA7CUb+ipbE72QCd0eTvIeZkxjfAvqo6ghgEceuAiJRJHwnTuTjGxNlJPBX4I0wx9MkEWkHvAb8WlUPH/9xPUXa5PfSxHFEzPeiqh5VHYVvDvJTRWTYcZsE9TuJlsTvZAL3mm1EJA7oQNv7+d7kcajqPlUt9y8+AYxppdhCwcn31uap6uHqn+rqm10uXkQywhxWg0QkHl+yfFFV/1nPJhHxvTR1HJH2vQCo6kHgI2DicR8FNX9FS+KvmcBdRBLw3fyYf9w284Hr/O+vBD5Q/52SNqTJ4ziurXUyvrbNSDUfuNbfi2Q8cEhVd4c7qOYSkazq9lYRORXff1f7whtV/fxxPgWsU9UHG9iszX8vTo4jUr4XEckUkXT/+2TgAmD9cZsFNX+FZc7dYNMGJnAXkbuApao6H9+/JM+LyEZ8Z8op4Yu4fg6P42YRmQxU4TuOn4Ut4CaIyFx8PSsyRGQncCe+G1eo6mP45l2+GNgIlADXhyfSxjk4jiuBX4hIFVAKTGmDFxXVzgB+CqzytykDzAR6Q0R9L06OI1K+l27AHBFx4zs5zVPVt0OZv2zIBmOMiTHR0tRjjDHGIUv8xhgTYyzxG2NMjLHEb4wxMcYSvzHGxJio6M5pTLCISGfgff9iFuABivzLJap6elgCMyaIrDunMQ0QkRzgqKr+b7hjMSaYrKnHGIdE5Kj/9RwR+VhE5onIBhG5T0Su8Y+pvkpE+vu3yxSR10Qkz/93RniPwBgfS/zGtMxIYDowHN8TpINU9VTgSeD/+beZBTykquOAK/yfGRN21sZvTMvkVY9fIyKbgPf861cB5/rfXwAMqTVsensRSfOPH29M2FjiN6Zlymu999Za9nLsvysXMEFVS1szMGOaYk09xoTOe8CvqhdEZFQYYzGmhiV+Y0LnZmCsf7a0tcDPwx2QMWDdOY0xJubYFb8xxsQYS/zGGBNjLPEbY0yMscRvjDExxhK/McbEGEv8xhgTYyzxG2NMjPn/ZFOhBd9s+skAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret2 = ensemble_simulations(duration, ndiv=20, y0=y0, model=m, solver=('spatiocyte', radius), repeat=N)\n", "ret2.plot('o', ret1, '-', **opt_kwargs)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Simulating with `egfrd`:" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEGCAYAAABiq/5QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXxU5dnw8d81M1khEAIhYd9B2QmL4NLiWl5R1LqU1lbr4yu01lf62IdqFTW1Vq22KvaxbkVFrSDFKiiKgqLWgkKI7DvInpWwZs/M/f4xk5jAZHKSyeRkZq7v5zOfmXNy7nOuMwPXnLnPvYgxBqWUUtHDYXcASimlWpYmfqWUijKa+JVSKspo4ldKqSijiV8ppaKMy+4ArOjUqZPp3bu33WEopVRYWbt2baExJvX09WGR+Hv37k1WVpbdYSilVFgRkX3+1mtVj1JKRRlN/EopFWU08SulVJQJizp+pZSyQ0VFBbt376akpMTuUAJKTEykX79+xMbGWtpeE79SStVj9+7dJCcnM2jQIByO1llB4vF4yMvLY9euXQwePNhSmZCeiYgki8hCEdkmIltFZIKIpIjIMhHZ6XvuEIpjr1n8ArmZ/fE82J7czP6sWfxCKA6jlIpgJSUlpKWltdqkD+BwOEhLS6OkpITc3FxrZUIc02xgqTHmLGAEsBW4B/jEGDMA+MS33KzWLH6BoWtnkU4BDoF0Chi6dpYmf6VUo7XmpF/N4XAgIixatAi3293w9qEKRETaAd8D5gAYYyqMMceAq4C5vs3mAlc397F7ZD9BglTUWZcgFfTIfqK5D6WUUnX86IVV/OiFVbYcu6ysjNLS0ga3C+VXWV+gAHhFRL4Rkb+LSBsgzRiTA+B77uyvsIhME5EsEckqKCho1IE7G//bdzaFjdqPUkq1Bu+88w4iwrZt2wJuJyKW9hfKxO8CMoDnjDGjgGIaUa1jjHnRGDPGGDMmNfWMHscB5Yv/7fOlU6P2o5RSjfHuN4f4Zv8xvv62iPMe+5R3vznULPudN28e559/PvPnz2+W/YUy8R8EDhpjvvYtL8T7RZAnIl0AfM/5zX3gAxkzKTV1mzWVmlgOZMxs7kMppRTgTfq/+9dGKtweAA4dK+V3/9oYdPI/deoU//nPf5gzZ07rT/zGmFzggIgM8q26GNgCLAZu9q27GVjU3MceO2U6m0Y/TC6peAzkejrw9YD/ZuyU6c19KKWUAuCJj7ZTWln3xmpppZsnPtoe1H7fffddJk2axMCBA0lJSSE7Ozuo/UHoW/X8P+AfIrIBGAk8AjwGXCoiO4FLfcvNbuyU6aRn7qJw2jpS5RhxJdaaOSmlVFMcPub/pmp9662aN28eU6dOBWDq1KnMmzcvqP1BiDtwGWPWAWP8/OniUB63ts7d+rCuzXgGHF5EZcUTxMTGtdShlVJRpGtyAof8JPmuyQlN3ueRI0f49NNP2bRpEyKC2+1GRHj88cct38j1p/U3UG0Oo26iE8fY9Nk/7Y5EKRWhZv5gEAkxzjrrEmKczPzBoHpKNGzhwoXcdNNN7Nu3j71793LgwAH69OnDl19+GVSsUZH4h068jnxScGS/ancoSqkIdfWobjz6w2HEOr1ptVtyAo/+cBhXj+rW5H3OmzePa665ps66a6+9ljfffDOoWKNirB5XTCy7e/yQc/bP4cCujfToP8zukJRSEejqUd2Yt3o/AG9NnxD0/j777LMz1t15551B7zcqrvgBBkyeQRUODn/0tN2hKKUi2FvTJzRL0g+lqEn8ndJ7sj75Yobmv8+JY0fsDkcppWwTNYkfIPnCO2kjZWz54G92h6KUUraJqsQ/YOQFbI0ZTI+dr+OuqrI7HKWUskVUJX6A0oxpdDN5bPi0ebo+K6VUuIm6xD/8khvJJZXYLB2bXykVAq9M9j5asahL/K6YWPYO+BlDKjawPetTu8NRSqmAnE4nI0eOZMSIEWRkZLBy5cqg9xl1iR9g2JQZHKMtpZ8+bncoSqlIsmEBHFwD+76Ep4Z6l4OUkJDAunXrWL9+PY8++ii/+93vgt5nVCb+NknJbO15IyNLVrFn09cNF1BKqYZsWADv3Qnucu/y8QPe5WZI/tVOnDhBhw7BT1MeFT13/Rl89UyKZ8/l6EePwdB37A5HKRXuPnkIKk8bpK2y1Lt++A1N3m1paSkjR46krKyMnJwcPv00+CrqqLziB2ifksqGrtcx8sQKDuzaaHc4Sqlwd/xg49ZbVF3Vs23bNpYuXcpNN92EMSaofUZt4gcYcNU9VOIi5/1H7A5FKRXu2ndv3PommDBhAoWFhTR2HvLTRXXi75Teg/WpUxh19CMO7dlsdzhKqXB28QMQc9rY+zEJ3vXNZNu2bbjdbjp27BjUfqI68QP0++GDVOEk593m+3CUUlFo+A1w5TPg9E321L6HdzmI+n34ro5/5MiR/OhHP2Lu3Lk4nc6GCwYQtTd3q3Xq2otVXadyzuHX+Xbz1/QZco7dISmlwtXwG2DtXO/rW5Y0yy7dbnfDGzVS1F/xAwy+/gFOSQLHlzxodyhKqXB3y5JmS/qhookfbwufLb1vYWTJKratWQ7A5kfOZ/Mj59scmVJKNT9N/D7Dr7ubQpJxL8vEeDx2h6OUUiGjid8nsW17dp99O0MqNrJueXDzWSqlVGumib+W0T/8b/Y6etJ51UNUecTucJRSKiQ08dfiionl5MSH6WbyyK+MszscpVQYumXpLdyy9Ba7wwgopIlfRPaKyEYRWSciWb51KSKyTER2+p6DH3GoGQ373lV8k3gu57KRk+4Yu8NRSkW53Nxcpk6dSr9+/Rg8eDCXX345O3bsCGqfLXHFf6ExZqQxZoxv+R7gE2PMAOAT33KrcrzLBZwwCYyryiY3sz9rFuukLUqphi3Zs4QNBRvIysvisoWXsWRPcM06jTFcc801TJw4kd27d7NlyxYeeeQR8vLygtqvHR24rgIm+l7PBT4D7rYhDr/WLH6Bc3Y9RYKjAoB0Cmi/dhZrgLFTptsbnFKq1VqyZwmZKzOp8HhzR05xDpkrMwGY3LdpM3KtWLGCmJgYfvGLX9SsGzlyZNCxhvqK3wAfi8haEZnmW5dmjMkB8D139ldQRKaJSJaIZAU7IFFj9Mh+ggSpqLMuQSrokf1Ei8WglAo/s7NnU+Yuq7OuzF3G7OzZTd7npk2bGD16dLChnSHUV/znGWMOi0hnYJmIbLNa0BjzIvAiwJgxY4Ibg7QROpsC8NOgp7MpbKkQlFJhKLc4t1Hr7RTSK35jzGHfcz7wDjAOyBORLgC+5/xQxtBY+ZLqd32eBDcanlIqsqW3SW/UeiuGDBnC2rVrm1y+PiFL/CLSRkSSql8DlwGbgMXAzb7NbgYWhSqGpjiQMZNSE1tnXamJZQP98YRgsCSlVGSYkTGDeGd8nXXxznhmZMxo8j4vuugiysvLeemll2rWrVmzhs8//7zJ+4TQXvGnAV+KyHpgNbDEGLMUeAy4VER2Apf6lluNsVOms2n0w+SYFDwGcknl0+Tr+AFfsWah1vMrpfyb3HcymedmEuvwXjh2adOFzHMzm3xjF0BEeOedd1i2bBn9+vVjyJAhZGZm0rVr16BibbCOX0SuB5YaY06KyCwgA3jYGJMdqJwxZg8wws/6I8DFTYy3RYydMp3Nm16niG4MufdLLvd4WP/4ToZv+QsHdv4fegw447SUUorJfSezcMdCAF6Z9Eqz7LNr164sWNB8E7aDtSv++31J/3zgB3ibYD7XrFG0cuJw0O2mOZRLLKVv/V+qKisaLqSUikqvTHql2ZJ+qFhJ/NUV25OB54wxi4DYANtHpE5de7F73EMMrNrBmjfutzscpZRqMiuJ/5CIvADcAHwgInEWy0Wc0ZffSlbSxYzZ+xLbsj6xOxylVAvwhMEw7Y2N0UoCvwH4CJhkjDkGpAAzGx9aZBjwXy9S4OhI8vvTOX7Ef7dpncRFqciQmJhIbm5uq07+Ho+H3NxcKisrLZdp8OauMaZERPKB84GdQJXvOaINufdLv+vbd+hE3pUv0XvRD9ky52ZG/M8HiCMqfwApFfH69evHunXrOHz4MCKtd6j2yspK9u/fjzEGh4V8ZKVVz4PAGGAQ8AoQA7wBnBdkrGFrYMZEvtp2F+N3PMFX8/7A+Bt1rl6lIlFsbCwJCQksW7aM9u3bW0qqdjl58iRpaWkkJCQ0uK2Vs7gGmAIUQ01v3KSgIowA50y9l28Sz2P0jtlsWfWh3eEopUJk6NChTJw4kfj4eESkVT4cDgcDBgxgypQpln6ZWBmrp8IYY0TEQE0v3KgnDgd9b3uN3GcuIP2jaeR0XU6XXoPsDksp1cxEhIyMDDIyMuwOpdlYueJf4GvVkywitwHLgZcaKBMV2nfohPnxPGJMJSWvTaXk1HG7Q1JKqQY1mPiNMX8GFgJv463nf8AY89dQBxYueg4cyZ7vP0Ofqm/Z9vzPMK347r9SSoHFYZmNMcuAZSGOJWyNuOgGvjq4kfF7nmHVnP+mnd0BKaVUAPVe8YvISRE5Uev5RO3llgwyHJzz09+zOuVKXPv/Q4eyQ5xdvlGnbVRKtUr1XvEbY6K+5U5jiMMBPc9l6JGPambw0mkblVKtUYN1/CLS09+jJYILNz3XPanTNiqlWj0rdfy1p4mPB/oA24EhIYkojOm0jUqpcGBlyIZhtZdFJAPQegs/8iWVdM6cGD5POtLFhniUUsqfRvc/9k3AMjYEsYS9+qZt3O3uzLdb1tgUlVJK1WVlrJ67ai068M7AdeZlrWLslOmsAbqvfYw0isiXVLb2v4Uhu/9OzIJr2HXNfPqP0FE7lVL2snLFn1TrEYe3zv+qUAYVzsZOmU5RXDe2xg0jPXMXF/7sPipvep8y4un8zvVsW7Pc7hCVUlHOSh3/71sikEjWre8Qcv/rQ068ciU93/8JG0tfYtj36v/urB7Lv76hoZVSKhhWmnMuE5HkWssdROSj0IYVedJ7DiD+to/Ic3Zh0Ce3kLX4ebtDUkpFKStVPam+mbcAMMYcBTqHLqTI1alrLzre+Sk744YwJvtuvnrtfh3bRynV4ixNtl67w5aI9AJM6EKKbO2SO9L/ro9Ym3Qh4/c8w5r/vYmK8jK7w1JKRRErif8+4EsReV1EXge+AH4X2rAiW1x8IqN+/Tarut7MuKL32PXniynKP2R3WEqpKGFlWOaleJtwvgUsAEYbYyzX8YuIU0S+EZH3fct9RORrEdkpIm+JSGxD+4hEDqeTCdOeIWvME/St2E75377PrvX/sTsspVQUsHJzV4BJQIYx5j0gUUTGNeIYM4CttZb/BDxljBkAHAVubcS+Is6YK6Zx4Op/4cRNj39dxeqFT2I8WpOmlAodMSZwkhGR5wAPcJEx5mwR6QB8bIxpsPeuiHQH5gJ/BO4CrsTb+SvdGFMlIhOATGPMDwLtZ8yYMSYrK8vSCYWrovxDHJrzU8pKTtJdjpAm3g5gBzJm6sieSqkmEZG1xpgxp6+3Usd/jjHmV0AZ1LTqsVo98zTwW7xfHAAdgWPGmCrf8kGgWz0BTxORLBHJKiiI/I7CKZ27UXr2DQxx7KOLowiHeId1Hrp2lo7pr5RqVlYSf6WIOPG15BGRVL5L5PUSkSuAfGPM2tqr/Wzq9yeHMeZFY8wYY8yY1NRUC2GGv57r/kKin2Gdu+uwzkqpZmQl8T8DvAN0FpE/Al8Cj1godx4wRUT2AvOBi/D+AkgWkeoew92Bw40NOlJ1Nv5/2aSZAvZty7a0j82PnF/T81cppfyx0qrnH3irax4FcoCrjTH/tFDud8aY7saY3sBU4FNjzI3ACuA632Y3A4uaGHvEyRf/v2zyTAe6zLuUVXPvpbKivIWjUkpFmkBz7qZUP4B8YB7wJpDnW9dUdwN3icguvHX+c4LYV0Spb1jnXUPuYFPSeUz49lkOPDaObat13nulVNMFGqRtLd769/rq5ftaPYgx5jPgM9/rPUBjmoNGDX/DOh8YPZMLpkwH/ofsj16n26oH6fvBdaxeNZkBP/kLHVJ1ihelVOMEmmy9T0sGorzGTpnO5k2vU0Q3htz7Jem1/pbxg59RfO6VfPWPexmdM5/SZ1fw1YBfkHHd3cTGxdsWs1IqvFiagUtEpojIn32PK0IdlKpfm6Rkxv/ibxyeuoxv44cwfueT5D82kuylr+qAb0opS6z03H0Mb+/bLb7HDBF5NNSBqcB6nT2aEfcsZ8PEl6kSFxlfzWDXI+PIqWyjPX+VUgFZueK/HLjUGPOyMeZlvMM3TA5tWMqq4ROvpce961g94mEOlycw1L2dwZWbyMnsx5pFz9kdnlKqFbI62XpyrdftQxGIajqny4U4Yxnn2Ea64ygOgS4UMjT7AZY+N9NyFZD2AVAqOlhJ/I8C34jIqyIyF29rHysduFQL6pH9BAl+ev2OzPknex8exZp3nw3puP/6paFU+LDSgWseMB74l+8xwRgzP9SBqcapr9dvZzkKGMauu5fjj57FqlfupjD3QMsGp5RqVeptzikiGaetOuh77ioiXY0x1sYQUI3WlEnW8yWVdM5M/vmSSu9Z61j/+dvI1y8wYd/zVDz3ElntLyRh/K0MHj8JcVit8VNKRYJAHbiygM1Qk01qd+QyeMfeUa3EgYyZtF87q051T6mJ5cDomaQ7HIy48Hq48Hr271hHzrK/cnbBB7T7eDn7l3fjcJ/r6H9xVE+LoFRUCZT4fwNcC5TiHWTtHWPMqRaJSjVafb1+Tx/Lv+fAkfQcOIfS4pOsWTaXpE1vMH73bFZv/4AeDu88ADmZ/TiY8VudB0CpCFXvb3xjzFPGmPOBO4AewCciskBERrZYdKpRxk6ZTlFcN7bGDSM9c1fAxJ3QJomxV9/BWbO+4tO+Mxnq/G4egC4UMmTtLD589jeUlRa34BkopVqClZu73+IdQfNjvGPsDAx1UKplDd7z6hnzACRKBaPy3sb9WF/W/uVq1n4wh5PHi0IWQzCtgrRFkVKNE+jmbl+8wylfBRzAW93zR2NM6NoEKlt0NgV+h+LrLEfJ6ngl/Ys+J2X1Ciq+nsnG+BEU976EbmOn0KP/MADWLH6B7uWHSKOI3Mz+YTVdZPUXRlNuqAdTVik7Barj3wVswHu1fwLoCdzunXsdjDFPhjw61SICtQgad+cbVFVWsGXtp5xYt5iu+Z8xbPvjsP1xDko661wjuLhiRc1N5XQKaL92FmsgbJK/HYL90tAvHRWMQIn/Ib6bFrFtC8SibBKwRRDgioll8PhJMH4SAAd3beLQ2veJ3/spo4tXk+DwN13k46CJX0WAW5beAsArk16xOZLmE2hY5swWjEM1k6ZcAVptEVSte/+hdO8/FLgHz4P+R/BIM4Xs/sMoCjuMwtlnAl2HfI8uPQdon4FWIBp/LdiVvIM5bihjDnTFr6JIoHkAAqmvmiiPjpS62jGs4H0SC9+GNXCE9hxIHExp6ggSe4+hx5BzSencLazvEdjl9528LasXtHBZO0XilbddNPGroNRXTXRw9N2MnTKdqsoKdm1ZzZHtK3EcWkvayY103/sVjn3Pw+fwhXsoYx07mnyPIJgvDf3CaTxNvi2neP83Idt3oFY9M4wxs0XkPGPMf0IWgQprDVUTuWJi6T/ifPqP+K655cnjRezf/BUn96xmwM65fgeX6772UbJ2LaOy40DiuwymY++hdOl9NjGxcTXbrVn8AkNrfek05ksjmLLNIVyvukOZjFTLCXTFfwswG/grcPq4PUrVaGw1UVL7FIacezmcezmeB5/yu00aR3GcyCbtxDL4FlgJlcbJAUcaRfE9KE3qTa+85X6/NHpkP9HgjeX6RjO1UrZauCZvu+ivhdYjUOLfKiJ7gVQR2VBrvQDGGDM8pJGpqBCoKWl65i5OHi8iZ/dGjh/YTFX+duKO7yG5ZD8DS74hjgo/e/T2S/jq73fh6NCThNRetEvrQ2q3PiS2bV9nG799F0xhs52bqiuYXwtL9ixhQ8EGKjwVXLbwMmZkzGBy39DPB2XXcUMtUKueH4tIOvARMKXlQlJ2saOVR0NNSZPap5CU8X3I+H6dcsbjIfehAXThzESdZzow9sDLOA/WnYLyBG0ocnTkRGwqqaTQhTN7IudJR7pYiHvN4hd46tAxunC0SfcWmlo2Gi3Zs4TMlZlUeLz/RnKKc8hcmQlgKQk3NXk3x3HX5WZRBa3uSyPgzV1jTC4wQkRi+W6ohu3GmMqQR6aiQmObklYTh4ODGb8l2d+N5TG/o9Okn5N36FuO5X5LScE+Ko8ewHHyMLElebStyGenpyvJjlNnlvV0oP2DnTnmSOakM5nSmA5UxKXgTkhBEjviaNuJI4f3MTH/tbC8txCOZmfPpsxdd8CAMncZs7NnN5hIg0nezXHcKt+vysZ+aYRag616ROT7wGvAXrw/jnuIyM3GmC9CHJuKEk1tStrQl0bXPmfRtc9ZfssOAFa/+xw9vnmcNIrIk45kd5hESkpHNhQX4CotIK68iKSKfJJKd5B89Dix4gYg19PBf6e1rEfZs+5vlLmSKHclURXTFk9MEp749khcEhLfnp6bnvV/Mzv7ccwVt1nq52DXL4Yle5awM9bd4lewucW5jVpfWzDJ267jQujfayvNOZ8ELjPGbAcQkYHAPGB0oEIiEg98AcT5jrPQGPOgiPTBO+5PCpAN/MwY47+yVqkGNPVLA2Dc1b9k85Z/1JQN9N/KeDycOnWcE0X5pL8yzu82aXKU9fFnE1t1kvZlh0koKaYNJbQ1JTjEW+1U3+zHaaYQz+9TKCGeUkmgXOIpdyRQ6Yin0pmA2xmP25lAbglcVPl5zaB66RTQbu0sluxZT9qAUThjEnDGJuCKS8QVl4ArNoGYuHhi4hJxxcaz/d9vN+lLI9gr2GASWXqbdHKKc/yub0gwyduu47bErwUriT+mOukDGGN2iEiMhXLlwEXGmFO+7b8UkQ+Bu4CnjDHzReR54FbguaYEr1qPSO8FKg4Hbdt1oG27DuQGuCE96rcfnrHe43ZzqvgEJSeP4X7pEv/3JejI3u5TkMpiHJXFOKtKcVaVEOMuIb7qBDEVBcSZMrq5K1iR5GJ2h07kupykV7mZcfQYo4s+IH3NPwKewxr3APKT8vh9j7bkutqSXuXml9sepmLNqwxyHqaKGColBrc4cUsMbnF5nx0x/LnzMcpcdb+2ytxlPP75/XT6YCHG4QJHDMbhRBwxGIcLcbjA4eQrDvOaJ6tOIrv/i3vZs/YTJiSejcPhAocLh9OFOJy+hwtxOnA4XFzf7kKeL/knFbVqmOMcsfwk7Ur2bV+Hw+nE4XAgDpfv2YHT4QKHg87xqeSV5Z/xXqQlplFeVoLD4UREap5r/+KakTGDzJWZda7c453xzMiYEfB9huC+NIL9tWCFlcSfJSJzgNd9yzfinXA9IGOMAaonbonxPapn7vqJb/1cIBNN/CoILd2LtaEb0qdzOJ01Xxpr6rsvMfpuJli48n7v8S481CmFMl+CyolxkdkphQfMUYZPWUFleQlV5aVUVZTirijBU1mOu6IMT0UpO3Y/y+zUdnXKPpLajhkU4or/Hg53BeKpxOGpxOGpwGGqcJgqnJ5KCp3+f6sUSSUDC5fhMlW4cOPEXVMlVu2e7j2pPO1SsVI8vH18GXdsnoMbcAu4EdyAR6qfveuuAOLaJPBq+3YUOp10cru58fhRJuz/PSXi/RXlQfCIN8G4a72+Ij6e19u3o8LxXROuGI/h0pwdfPlMX19ZMCIY32tvHIIH4fK4WD5pm8Bxh9DOY/j+iSKK/vUr5sodvu0FT01ZbxkDDI51UNA2hir57rhOY+ifd4A/P302BnwP7/bVxwchp734bXHm74ukqawk/l8CvwLuxBvOF8DfrOxcRJx4vyT6A88Cu4Fjxpgq3yYHgW71lJ0GTAPo2bOnlcOpMPbnjAEAhEML7+p7C13XPkoXjlq+IV1d9pny7Sw6/hEFLiG1ynBV+x9wp8U6+qdTvkv61cocDp5O6cCyQSMpd5dT4a6g3F1e53WFu4LXjr7gt+xrKS5+OuFSKjwVVHoqqXRXUuWp8r72eF933vJP8mPOvP/QucrDjHMuospT5X2YKqrc3tduU0WVx01RSR7+MlmRyzC8T+P/bxe4XDzdsQNPN7qkV6VDeL19O173P8xUvU44hffaJfKexe2HHm9LfuKJms+5c0k7/t3+FP9uoFznSo//97qyvorCxmsw8RtjyvHW8zd6GGZjjBsYKSLJwDvA2f42q6fsi8CLAGPGjPG7jVIAbXqOavFj5g/tzrSiJKpIokvbLswY2t1Sufd2vcdr5V9Q7vuPnR8jzC37nNLVf2Jwx8GUVpVSUlnifa7yPtd+5Lv83/zNdzkY8dqIgMcWl9Pv+lyXkyeynqhZdoiDGEcMLoeLWEcsMY4Yfl1UxEOpdb904j0efl1UxKIesSS6EnE6nLjEhcvhqnntdDj5eutCcmLOTDXpVW6uGXMnDnHgdDhxitP7uvazw/tc9dXXnLfrLe+9CenEqgE/ov35F3mrafCWd4gDQWrKinhf53/2IcO3/J0uHCVPOrFh8K30vPgaHOKoOd/qstVlql/vXraAAev/WlN25/A7OGvSjTVlq49Xvb2IIAgbP3iVUXseIqGo7i+7Db1+T8YVt4F4y9aU8z2D71ed3/f6aMDPtzFaZKweY8wxEfkMGA8ki4jLd9XfHTjcEjEo5U9jvzSMMbyz8x0eWf1I3Trr/9zPZ/s/o2tSV05WnPzuUXmSUxWnvI/KU5RUlZyxzwpPBW9sfeOM9fHOeBJcCSS4Eoh3xRPviifWGVvTNLHOecS24abBNxHrjCXeGU+cK45YRyxxrjjiHHHEOmO5Z/ldHOPM43dwJPH+1I+JccQQ44jB6TjzCyI3sz+ZhUXM7pBc995CSRuu/MHfA75n/1i5gKc7ec5IZDcf9XDjyF8GLAu+JrC7Xq6pHutKIZfveIlNSb0sNZ+duGV2TdkuFJK8+Sk2xXe0VPac9Y/VLbv+ETa5Ehss23fdbL+tt3p98yQxV90esOzYkjb1vtfNJWSJX0RSgUpf0k8ALgH+BKwArsPbsudmvBO9qDAXzt3xPcbDoVOHKCwt5BVNoOoAABbTSURBVEjpEYrKimqej5Ydpajc+3ys7BjHyo/5TbyVnkqW7luKy+GiXWw72sW2o21MW9rGtiUtMa3m9etbXvcTgffq771r3iPRlUhiTCLxzni/Cbi6xcfpNxzvH39/gzf+7rngAR749ywqqKpZF4uL354/i6TYpIBlD2TM5KK1s5hc/N11WqmJZVM99zVqGzjoLu7d9jDPpbStSWS/LDpF97NmNVDSK5jhNewqG0zP8GDea6sCJn5fHf1jxpiZTdh3F2Cubx8OYIEx5n0R2QLMF5GHgW+AOU3Yt1INMsZQVFZEbkkuucW55BXnkV+S732U5rPn2B4KSr2tcya9PemM8kkxSXSI70CH+A50bduVIR2HkByfzCub/H+5CUL2T7NrfrL7s3zf8npbe/Rq16vBc6pO7rO+uIcq8FYzWWwaGUzZYO9rALxSu2zGLMt9D4JJonaVrX8okk4NJu9g3murGuq56xaR0SIivlY6lhljNgBn/I42xuzBO2m7UkDTu9QbY1iwYwHrCtZR5ani3DfPZVDKIFwOFznFOeScyjnj6tzlcJGakIrL4eJI2ZE6f4t1xDJt+DSu6n8VKfEpxDpj/R536bdL603egZI+BNdEsNrkvpOZu/w+ABZc97HlcsGWHTtlOjfkPgcks2DaukZdfQZTNpgkalfZxrb8Ol0w75cVVqp6vgEWicg/geLqlcaYfzVzLCoKWelSf6riFHtP7OXb49+y98Re9p3Yx/4T+9l9bHedxH6y8iRZeVl0b9udwR0Hc2GPC0lvk/7dIzGdDvEdcIiDyxZehsfUbSVR4ang7Z1vM31E4CurYJJ3MFfd0SqYJGpX2aYORdJSrCT+FOAI3vb31QygiT/C2FFPX19nlT+s+gPv7nqXPcf2kF/6XQccpzjp1rYbPdv1ZN+JfX7r293GzV8m/iXgcYPpWRls8g7mqjucPVjYtKm7g0midpWtLt/UXuWhZqU55y0tEYiKLoWlhWw5sqXeTinFVcWcrDjJ+K7j6dO+D33b96V3+970aNuDGKe3N9Dwuf5HBg91d3zwJu+FOxYC4XlD2w7B9O4OJonaVbY1szJI20C8PWvTjDFDRWQ4MMUY83DIo1NhoaE6+uLKYjYXbmZD4QY2Fmxk85HN5JXkBdxnlzZdmH/F/IDbBJO8m6OuPRh29D1oDk29aleti5WqnpeAmcAL4L1pKyJvApr4ld86+gdXPsi6/HUArCtYx46jO2rq03u168XotNEM6TiEIZ2GsPf4Xh5b/ViTEnBz1LU/8J8HqPBU0KWN1rWr1iWUX7JWEn+iMWb1aa0VqurbWEUXf3X05e5y5m+fT4IrgRGpI5g2fBojUkcwrNMw2sfV7Sc/Om008a74JiXgYJN3OFfXBPOLIZiydg7GF+kDAbYkK4m/UET64RtaQUSuA5pvtCDVrJp6g7YxTSor3BVk52fz5cEv662jF4SVP16Jy9HwP7FgEnA4J2+7ROP7ZNeXRjDHDWXMVhL/r/COmXOWiBzCO/X1jSGLSLU4K00qC0sL+ffBf7PiwAq+yvmK0qpSYhwxxDr8DyGQ3ibdUtJXqrWLxF8aVlr17AEuEZE2gMMYczL0YamWVF+TyieznqSwtJDl+5azvmA9BkN6m3Sm9JvCBd0uYGz6WFYcWGHrTdJwFexVdzRetavmY6VVT0fgQeB8wIjIl8BDxpgjgUuqcFFf88f80nz+nPVnBnUYxC9H/pILe1zIoA6D6vRODfebpMEkUE2+KlxZ+S0+H+8Y/Nf6lm8E3sI76JqKAGmJaeSWnJn8k2KSmH/FfHq2CzxmenPUs2sCVqrlWOq5a4z5Q63lh0Xk6lAFpFqGMYaNhRtZuGMhR0rP/PEW74znvvH3NZj0lVLhx0riXyEiU/ludrrrgCWhC0mFUkllCe/tfo8FOxaw4+gOElwJXNn/SjoldOLVTa+GZXWNUqpx6k38InISbxNOwTtBevVMEQ68c+k+GPLoVKMEapJ54MQB3tz2Ju/uepdTlac4O+VsHpjwAJf3uZw2Md4JHrLzsgGtOlEq0tWb+I0xgWdmUK1KfU0y95/Yz85jO1m+bzlOh5PLel3GT87+CcM7DW9wCGGlVGSy1NDaNz5P79rb67DModOUTlj1Ncn82/q/kRSTxK3DbuUnZ/2E1MTUZo1VKRV+rDTnfBkYDmwGqgcw12GZW5lAI1Iuu35ZTXVOqGj1kFLhw8oV/3hjzOCQR6KCkhKfcsaMUuAd5dJq0tfkrVR0cDS8CatERBN/K7XvxD5mfDrDb9LXHrRKKX+sXPHPxZv8c4FyvK18jDHG/ywYqkUcLz/O8+ufZ/62+cQ6Y7lj5B2kJqTyx6//qE0ylVIBWUn8LwM/AzbyXR2/sonHeFi0axFPrX2K4xXHuab/Ndwx6g46JXQC4L097wFabaOUqp+VxL/fGLM45JGoBu04uoOHv3qYb/K/YVTnUdx3zn0MShlkd1hKqTBjJfFv88249R7eqh5Am3OGir9OWJf2upQXN7zInI1zSIpN4g/n/YEp/abgECu3aJRSqi4riT8Bb8K/rNY6bc4ZAv46YT2w8gGezHqS/NJ8pvSbwswxM0mOT7Y5UqVUOLMyHv8tTdmxiPQAXgPS8d4beNEYM1tEUvCO7tkb2AvcYIw52pRjRBp/nbAq3BUUlhXy7MXP8r3u37MpMqVUJLHSgesVfNMu1maM+a8GilYBvzHGZItIErBWRJYBPwc+McY8JiL3APcAdzc68ghUXycsY4zlpK83dZVSDbFS1fN+rdfxwDXA4YYKGWNy8M3Na4w5KSJbgW7AVcBE32Zzgc/QxA94pyv0N4dtept0G6JRSkUqK1U9b9deFpF5wPLGHEREegOjgK+BNN+XAsaYHBHpXE+ZacA0gJ49I39M+Ap3Bb3b9T4j8WsnLKVUc2tKs5ABgOVMLCJtgbeBXxtjTlgtZ4x50RgzxhgzJjU1vAYWu2XpLTUDrVmRW5zLz5f+nFU5qzi/6/nEOGIA73ALmedmaicspVSzslLHXz0uf7VcLFbNiEgM3qT/j1rNP/NEpIvvar8LkN/ImCPK+oL1/HrFrympLOGpiU9xSa9LmjQ6p1JKWWWlqqdJ4/KLd7D3OcBWY8yTtf60GLgZeMz3vKgp+48Ei3cvJnNlJmmJabx06Uv079Df7pCUUlEg0AxcAatzjDH7G9j3efiGehCRdb519+JN+AtE5FZgP3C99XAjgzGGv63/G8+vf55x6eP4y/f/om3zlVItJtAV/xK+m3qxmgFSgc6AM9COjTFfnla2tosbEWNEqfRU8odVf+CdXe9wdf+reWDCAzV1+kop1RICTb04rPayr2XO3cAlwCMhjSpClVaV8pvPfsO/D/2b6cOn86uRv9LpD5VSLc7Kzd0BwH3AOcBfgDuNMZWhDizSFFcWc/vy21lXsI77x9/PDYNusDskpVSUClTHPxRvwh8CPA7caoxxt1RgkeR4+XFuX347m49s5k8X/IlJfSYF3F5b8yilQinQFf964ADeuv5xwLja1RLGmDtDG1pkOFZ2jNuW3cbuY7t5cuKTXNTzIrtDUkpFuUCJv6GxeJQftYdVvuSfl+AUJ4WlhTxz0TOc3+18u8NTSqmAN3fntmQgkeD0YZXzSvIAuHnwzZr0lVKths7k0Yz8DasM8PG+j22IRiml/NPE34zqG1a5vvVKKWUHTfzNqL7hk3VYZaVUa9Jg4heRgSLyiYhs8i0PF5FZoQ8t/IxIHXHGOh1WWSnV2li54n8J+B1QCWCM2QBMDWVQ4eijvR+xdO9ShncaTqwjFtBhlZVSrZOVGbgSjTGrTxtaoCpE8YSlrUe2MuvLWYxMHcmcH8xh+rLpgHbEUkq1Tlau+AtFpB++MflF5Dp8UyoqOFJ6hDtX3En7uPY8deFTxDpj7Q5JKaUCsnLF/yvgReAsETkEfAv8NKRRhYlKdyV3fXYXR8uOMvf/zKVTQie7Q1JKqQZZmYhlD3CJiLQBHMaYk6EPKzw8ufZJsvOzeeyCxxjScYjd4SillCVWRueMA64FegOu6rp+Y8xDIY2slftk3ye8sfUNbjz7Rr15q5QKK1aqehYBx4G1QHlowwkPB08e5P7/3M/QjkO5a/RddoejlFKNYiXxdzfGBB5HOALVN+F5pbuSmZ/PBOCJ7z+hN3OVUmHHSuJfKSLDjDEbQx5NGPjrN39l05FNPD3xabondfe7jTbjVEq1ZoEmYtkEeHzb3CIie/BW9QhgjDHDWybE1mNt3lpe3fwq1w+8not7Re20wUqpMBfoir8bMLKlAmntiiuLue/L++jWthv/M+Z/7A5HKaWaLFDi/9YYs6/FImnlnljzBDnFObw66VUSYxLtDkcppZosUOLvLCL1NlkxxjwZgnhapS8OfsHbO9/m1qG3MqrzKLvDUUqpoARK/E6gLd46/ahSe/rES/95KaVVpfRP7s/tI2+3OzSllApaoMSfE0wnLRF5GbgCyDfGDPWtSwHewtsZbC9wgzHmaFOPEQqnT5+YW+KdROXHZ/1Ym24qpSJCoEHagr3SfxU4vf3/PcAnxpgBwCe+5ValvukTF+1eZEM0SinV/AIl/qDaKxpjvgCKTlt9FVA9iftc4OpgjhEKOn2iUirS1Zv4jTGnJ+3mkGaMyfHtPwfoXN+GIjJNRLJEJKugoCAEofin0ycqpSJdq51z1xjzojFmjDFmTGpqaosdd0bGDOKccXXW6fSJSqlI0tKJP09EugD4nvNb+PgNmtx3cp25c3X6RKVUpLEyVk9zWgzcDDzme251d0y3FW1jTe4aOid2pmdSTx13RykVcUJ2xS8i84BVwCAROSgit+JN+JeKyE7gUt9yq2GM4U+r/0RyXDJd23S1OxyllAqJkF3xG2N+XM+fWu3oZsv3LycrL4v7x9/Ph99+aHc4SikVEq325m5Lq/JUMTt7Nv2T+3PtgGvtDkcppUJGE7/Pol2L2HdiH3eOuhOnw2l3OEopFTItfXO3VSp3l/Pc+ucYnjqciT0mAjqZilIqcukVPzB/23zySvKYMWoG1ZPJK6VUpIr6xH+q4hR/3/h3JnSZwLgu4+wORymlQi7qE//rW1/nWPkx7ZmrlIoaUZ34iyuLeWPLG1zY40KGdBpidzhKKdUiojrxL9i+gBMVJ7ht2G12h6KUUi0mahN/WVUZczfPZXyX8QxLHWZ3OEop1WKiNvG/u+tdjpQdYdrwaXaHopRSLSoqE3+lp5JXNr3CyNSRjEkbY3c4SinVoqIy8X+w5wMOFx/mtuG3abt9pVTUibrEb4zh1c2vMrDDQC7odoHd4SilVIuLusS/Onc1u47t4qdn/1Sv9pVSUSnqEv8bW94gJT6Fy/tebncoSilli6hK/PtP7Ofzg59z/cDrz5hXVymlokVUJf43t72J0+HkR4N+ZHcoSillm6hJ/KcqTvHurneZ1HsSqYmpdoejlFK2iZrE/+6udymuLOanZ//U7lCUUspWUZH4jTG8tf0thqcO18HYlFJRLyoSf3Z+NntP7OX6gdfbHYpSStkuKhL/2zvepm1MWy7rdZndoSillO0iPvEfLz/Ox/s+5vI+l5MYk2h3OEopZbuIT/xL9iyh3F3OtQOvtTsUpZRqFWxJ/CIySUS2i8guEbknFMdYsmcJly28jEdXP4rL4eLb49+G4jBKKRV2XC19QBFxAs8ClwIHgTUistgYs6W5jrFkzxIyV2ZS5i4DoMpTRebKTAAm953cXIdRSqmwZMcV/zhglzFmjzGmApgPXNWcB5idPbsm6Vcrc5cxO3t2cx5GKaXCkh2JvxtwoNbyQd+6OkRkmohkiUhWQUFBow6QW5zbqPVKKRVN7Ej8/sZCNmesMOZFY8wYY8yY1NTGDbGQ3ia9UeuVUiqa2JH4DwI9ai13Bw435wFmZMwg3hlfZ128M54ZGTOa8zBKKRWWWvzmLrAGGCAifYBDwFTgJ815gOobuLOzZ5NbnEt6m3RmZMzQG7tKKYUNid8YUyUidwAfAU7gZWPM5uY+zuS+kzXRK6WUH3Zc8WOM+QD4wI5jK6VUtIv4nrtKKaXq0sSvlFJRRhO/UkpFGU38SikVZcSYM/pOtToiUgDsa2LxTkBhM4Zjp0g5l0g5D9Bzaa0i5VyCPY9expgzesCGReIPhohkGWPG2B1Hc4iUc4mU8wA9l9YqUs4lVOehVT1KKRVlNPErpVSUiYbE/6LdATSjSDmXSDkP0HNprSLlXEJyHhFfx6+UUqquaLjiV0opVYsmfqWUijIRk/gbmsBdROJE5C3f378Wkd4tH2XDLJzHz0WkQETW+R7/1444rRCRl0UkX0Q21fN3EZFnfOe6QUQyWjpGKyycx0QROV7rM3mgpWO0SkR6iMgKEdkqIptF5IxJKsLhc7F4HmHxuYhIvIisFpH1vnP5vZ9tmjd/GWPC/oF3eOfdQF8gFlgPDD5tm9uB532vpwJv2R13E8/j58D/2h2rxfP5HpABbKrn75cDH+KdlW088LXdMTfxPCYC79sdp8Vz6QJk+F4nATv8/Btr9Z+LxfMIi8/F9z639b2OAb4Gxp+2TbPmr0i54rcygftVwFzf64XAxSLibxpIO4V8IvqWZIz5AigKsMlVwGvG6ysgWUS6tEx01lk4j7BhjMkxxmT7Xp8EtnLmnNet/nOxeB5hwfc+n/Itxvgep7e6adb8FSmJ38oE7jXbGGOqgONAxxaJzjpLE9ED1/p+gi8UkR5+/h4urJ5vOJjg+6n+oYgMsTsYK3zVBaPwXmHWFlafS4DzgDD5XETEKSLrgHxgmTGm3s+kOfJXpCR+KxO4W5rk3WZWYnwP6G2MGQ4s57urgHAUDp+JFdl4x0QZAfwVeNfmeBokIm2Bt4FfG2NOnP5nP0Va5efSwHmEzedijHEbY0binYN8nIgMPW2TZv1MIiXxW5nAvWYbEXEB7Wl9P98bPA9jzBFjTLlv8SVgdAvFFgpWPrdWzxhzovqnuvHOLhcjIp1sDqteIhKDN1n+wxjzLz+bhMXn0tB5hNvnAmCMOQZ8Bkw67U/Nmr8iJfHXTOAuIrF4b34sPm2bxcDNvtfXAZ8a352SVqTB8zitrnUK3rrNcLUYuMnXimQ8cNwYk2N3UI0lIunV9a0iMg7v/6sj9kblny/OOcBWY8yT9WzW6j8XK+cRLp+LiKSKSLLvdQJwCbDttM2aNX/ZMuduczP1TOAuIg8BWcaYxXj/kbwuIrvwflNOtS9i/yyex50iMgWownseP7ct4AaIyDy8LSs6ichB4EG8N64wxjyPd97ly4FdQAlwiz2RBmbhPK4DfikiVUApMLUVXlRUOw/4GbDRV6cMcC/QE8Lqc7FyHuHyuXQB5oqIE++X0wJjzPuhzF86ZINSSkWZSKnqUUopZZEmfqWUijKa+JVSKspo4ldKqSijiV8ppaJMRDTnVKq5iEhH4BPfYjrgBgp8yyXGmHNtCUypZqTNOZWqh4hkAqeMMX+2OxalmpNW9ShlkYic8j1PFJHPRWSBiOwQkcdE5EbfmOobRaSfb7tUEXlbRNb4HufZewZKeWniV6ppRgAzgGF4e5AONMaMA/4O/D/fNrOBp4wxY4FrfX9TynZax69U06ypHr9GRHYDH/vWbwQu9L2+BBhca9j0diKS5Bs/XinbaOJXqmnKa7321Fr28N3/KwcwwRhT2pKBKdUQrepRKnQ+Bu6oXhCRkTbGolQNTfxKhc6dwBjfbGlbgF/YHZBSoM05lVIq6ugVv1JKRRlN/EopFWU08SulVJTRxK+UUlFGE79SSkUZTfxKKRVlNPErpVSU+f/5A+2ieKMsXAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret2 = ensemble_simulations(duration, ndiv=20, y0=y0, model=m, solver=('egfrd', Integer3(4, 4, 4)), repeat=N)\n", "ret2.plot('o', ret1, '-', **opt_kwargs)" ] } ], "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 }