{ "cells": [ { "cell_type": "markdown", "metadata": { "collapsed": true }, "source": [ "# Birth-Death\n", "\n", "This is for an integrated test of E-Cell4. Here, we test a simple birth-death process 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` and `radius` mean a diffusion constant and a radius of molecules, 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 # 0.01\n", "radius = 0.005" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "N = 20 # a number of samples" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "y0 = {} # {'A': 60}\n", "duration = 3\n", "V = 8" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Make a model for all algorithms. No birth reaction with more than one product is accepted." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "with species_attributes():\n", " A | {'radius': radius, 'D': D}\n", "\n", "with reaction_rules():\n", " ~A > A | 45.0\n", " A > ~A | 1.5\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": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEGCAYAAACKB4k+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dd3xc9ZX38c9RlyXLsiy527hgYzAuGGNMDQRCC+CwQGLCEiCUhJCFTXbzhJQNbNhkU3jIQ3qAODEBbDqYEkoIJRCwce827pYly3JT7zrPH3MlhFEZyxqNRvN9v173dcvc0Zzrke/R/VVzd0RERAASoh2AiIj0HEoKIiLSTElBRESaKSmIiEgzJQUREWmWFO0AjkRubq6PGjUq2mGIiMSUJUuW7HX3vNZei+mkMGrUKBYvXhztMEREYoqZbW/rNRUfiYhIMyUFERFppqQgIiLNlBRERKSZkoKIiDRTUhARkWZKCiIi0iym+ymIiMQyd6e2oZHq2kaq6hqoqmugumld20B1fQPVdY1U1320rqkPrU88qj9njm+1/9kRUVIQEQmDu1NV10B5dT3lNS2W6noqauupqGmgoqaeitoGKoN1VW2wX1tPZW0DVbUNVAZL082/obFzc9rcctZYJQURkc5qbHTKa+spqazjYGUdpdV1lFTVUVoV2i6tqg/WdZRVh7bLquuDpY7ymnrCvX/3SUmkT0oSGamJpCcnkpGaRGZqEnmZqfRJSSQ9JYn05MRgO5G0pu3kRNKSE0hLTmyxJJCW1GI7OZGUxAQSEiwi/05KCiIScxobnZKqOvZV1LI/WA5Ufnz7YGUdByprKWlaV9W1e1NPTDD6piWRlZZMVnroJj4ipw9905Lom5pE37Rk+qYlkZGaFFqnJDXf7DNSE4N16GYfqRt2d1BSEJEewd05WFnHnrIaistqKC6vDq3LathbXsve8o/W+ytq2yx2SU9OJCcjhf4ZyWSnpzAsO53sPqHt7D7J9EtvsfRJJisttN0nJRGz2L2ZdxUlBRGJuNr6RopKqyksqaawpKp5u6i0mqLSGopKq9lTVkNtfeMn3puWnEBuZiq5makMy05nyvB+5GSkMCAzldzMFHIyPlr690khLTkxClfYeygpiMgRq6ipJ/9AFfkHKsk/UMWug1XsCtYFB6soLq/BD/nDPiMlkcH90hiUlcZJo3IYmJXKwL5pDOybGlqy0sjNTCEzNUl/wXcjJQUR6ZC7U1xWw7Z9lWzbV8GOfZVs31/Jjv2V7Nxfyf6K2o+dn5KUwLDsdIZlp/Op8XkMzU5naHYag/ulM7RfGoP7pdE3LTlKVyPtUVIQkWZl1XVsKa5gy97yYF3B1uIKtu2roLK2ofm8xARjWHY6I3P6cP7EQYzI6cPw/n0Y0T+dYf3Tyc1IjenK1nimpCASh0qr6/iwqIwNu8vZWFTG5uJyPiwqZ3dpdfM5iQnGiP7pjM7N4OQxOYzOzeCoARmMGtCHodnpJCdqQITeSElBpBerb2hk694K1haWsq6wjA27S9mwu4yCko9u/n1SEjl6YCanjh3A2IGZjM3L5OiBGYzMySAlSTf+eKOkINJL1NQ3sHF3Oat2lbC6oIQ1u0pYv7uMmqBFT3KiMTYvkxmjcxg/uC8TBvdl3MC+DMtOV1GPNFNSEIlBjY3O5uJylu08yIqdB1mZX8L63aXUNYSa+GSlJXH8sH5cM/MojhuaxbFDshibl6m//KVDSgoiMaC0uo6l2w+wdMdBlm4/wIqdBymrqQegb2oSk4b344bTxzB5eD+OH9qPETnpasYpnaKkINID7SmtZuHW/Szaup8Ptu1nQ1EZ7pBgcMzgLC6dOpQTRvZn6ohsxuRmqPhHuoySgkgPUFxWw3tb9vHe5n28v2UfW/dWAKFK4Gkj+3Ph8UOYPiqUBDJS9d9WIke/XSJRUFXbwPtb9/Huh3t5Z9Ne1u8uAyAzNYkZo3O4asYITh49gIlDs0hS00/pRkoKIt3APVQx/Mb6Yt7aWMyibfuprW8kJSmBk0b151vnH8NpR+dyvJKARJmSgkiE1NQ3sHDLfv62rog3Nuxh5/4qAMYNzOSamUdx5vg8ZozKIT1FA7hJz6GkINKFSqvreGP9Hl5dU8RbG4spr6knLTmB04/O5StnjuXsCQMZlp0e7TBF2qSkIHKEDlbW8uraIl5aVci7m/ZS1+DkZqZyyZQhnHvsIE47OlfDOUvMUFIQ6YSy6jpeWVPECysLeOfDvdQ3OsP7p3PdqaO44PjBnDCiv5qJSkxSUhAJU219I29s2MOzy3bx+vo91NY3Mrx/OjecMZrPThrCpGH91GFMYp6Sgkg73J2V+SU8uSSf51cWcLCyjtzMFL44YySXTBnKtJHZSgTSqygpiLRif0UtTy/N54nF+WwoKiM1KYHzJg7mX04YxhnjctVsVHotJQWRgLvz/pb9zFu0g5dX76a2oZGpI7L58WWTuHjKELI0U5jEASUFiXul1XU8vSSfv7y/nc3FFWSlJfHFk0dy1YyRHDO4b7TDE+lWSgoStzYXl/Pnd7fx1NJ8KmsbmDoim3uunMLFk4eoCanELSUFiSvuzrub9vHgO1t4c0MxKYkJXDJlKNeeehSTh2dHOzyRqItYUjCzEcBDwGCgEbjf3e8zsxzgMWAUsA34vLsfsFATjvuAi4BK4Dp3Xxqp+CS+1Dc08uKqQu5/ewtrCkrJzUzlG+eO54snjySvb2q0wxPpMSL5pFAP/Ie7LzWzvsASM3sNuA543d1/YmZ3AHcA3wYuBMYFy8nA74K1SKdV1zXw5JJ8/vD2Znbur2JsXgY/vXwSs6YOUxGRSCsilhTcvRAoDLbLzGwdMAyYBZwVnDYXeJNQUpgFPOTuDrxvZtlmNiT4OSKHpaq2gUcWbucPb2+huKyGqSOy+cHFEzlnwkD1NBZpR7fUKZjZKOAEYCEwqOlG7+6FZjYwOG0YsLPF2/KDYx9LCmZ2M3AzwMiRIyMat8SepmTw+7c2s7e8llPHDuC+2VM5ZcwAdTITCUPEk4KZZQJPAf/u7qXt/Mds7QX/xAH3+4H7AaZPn/6J1yU+1dY3Mv+DHfzq75soLqvhtKMH8NtzxjNjdE60QxOJKRFNCmaWTCghPOLuTweHi5qKhcxsCLAnOJ4PjGjx9uFAQSTjk9jX2Og8u3wX9762kfwDVcwYlcNvvjhNyUCkkyLZ+siAPwLr3P3eFi8tAK4FfhKsn2tx/OtmNp9QBXOJ6hOkPf/4sJgfv7SedYWlHD8six9dNokzx+WqmEjkCETySeE04BpglZktD459l1AyeNzMbgB2AFcGr71EqDnqJkJNUq+PYGwSwzYWlfE/L67j7Y3FjMhJ55dXncDFk4aoAlmkC0Sy9dE7tF5PAHBOK+c7cGuk4pHYd6Cill/8bSOPLNxBRkoi3//ssVxzylGkJqlpqUhXUY9m6fEaGp15i3bw81c2UFZdx9UnH8U3PjOenIyUaIcm0ut0mBTM7Erg5aCvwfeBacD/qLexdIdlOw7wX8+tZvWuUk4encMPZx2vQepEIiicJ4X/cvcnzOx04HzgHtTbWCKstLqOn7+8gYcXbmdg31R+edUJXDJ5iCqRRSIsnKTQEKw/C/zO3Z8zs7siF5LEu7+uKuTOBWvYW17DdaeO4j/OO4bMVJV0inSHcP6n7TKzPwDnAj81s1RA005Jl9tTVs0Pnl3Dy2t2M3FoFg9eO10jl4p0s3CSwueBC4B73P1g0OHsW5ENS+KJu/PMsl389/Nrqapr4I4LJ3Dj6aM15aVIFHSYFNy90sz2AKcDHxIa/fTDSAcm8WFveQ3ffXoVr64t4sSj+vOzKyYzNi8z2mGJxK1wWh/dCUwHjgH+BCQDDxPqnCbSaa+s2c13n15FWXU9371oAjecPoZEdUATiapwio8uIzTC6VIAdy8I5kcQ6ZTK2nrufmEt8xbtZOLQLObdPJXxg/QrJdIThJMUat3dzcwBzCwjwjFJL7Z6Vwm3zV/G1r0V3HLWWL5x7nhSklR3INJThJMUHg9aH2Wb2U3Al4EHIhuW9Dbuztx/buNHL61jQEYqj9x4MqeOzY12WCJyiHAqmu8xs88ApYTqFX7g7q9FPDLpNUqq6vj2kyt5ec1uzpkwkHuunEJ/DVEh0iOF1SMoSAJKBHLY1hSUcMvDSyk4WMX3LjqWG88YrV7JIj1Ym0nBzMoIzXxmfHwGNCM0qGlWhGOTGPfUkny++8wq+vdJ4bGvnMKJR/WPdkgi0oE2k4K7qzmIdEptfSN3v7CWv7y/nVPGDOBXXzyB3MzUaIclImEIp5/CyNaOu/uOrg9HYt2+8hpueWQpi7bu5ytnjuFb5x+jnskiMSScOoUXW2ynAaOBDcDEiEQkMWttQSk3PbSYveU13Dd7KrOmDot2SCJymMJpfTSp5b6ZTQO+ErGIJCa9traI2+cvIystmSe+eooGshOJUYc9HrG7LzWzkyIRjMQed+dP727j7hfXMmlYPx780nQGZqVFOywR6aRw6hS+2WI3gdDMa8URi0hiRkOj88Pn1zD3ve2cP3EQ/+8LJ5CeovmSRWJZOE8KLVsh1ROqY3gqMuFIrKiua+C2ect4dW0RN50xmu9ceCwJGsxOJOaFU6fw390RiMSOg5W13Dh3MUt2HOCuS47jutNGRzskEekiHbYVNLPXzCy7xX5/M3slsmFJT1VYUsWVv3+Plfkl/PqqaUoIIr1MOMVHee5+sGnH3Q+Y2cAIxiQ91La9FVz94EJKquqY++UZnDJ2QLRDEpEuFk5SaDCzkU2d1czsKD4+7IXEgfW7S/nXBxfR6M68m2YyaXi/aIckIhEQTlL4HvCOmb0V7J8J3By5kKSnWbHzIF+as4j05EQevvFkjh6oEVBEeqtwKppfDjqszSQ0GN433H1vxCOTHmHJ9gNcN2cR2RnJPHrjTEbk9Il2SCISQeFUNBtwATDN3Z8H+pjZjIhHJlH3wbb9fOmPCxmQmcJjN5+ihCASB8IZqey3wCnAVcF+GfCbiEUkPcIH2/Zz7ZxFDMpKY/7NpzA0Oz3aIYlINwinTuFkd59mZsugufWRps3qxZqKjAb3S2P+TTM1bIVIHAnnSaHOzBIJWhyZWR7QGNGoJGpW5h/kujmLyOubyjwlBJG4E05S+CXwDDDQzH4EvAP8OKJRSVSsKyzlmj8uol+fZB69aSaDlBBE4k44rY8eMbMlwDmEWh99zt3XRTwy6Vbb9lZwzR9DzU7n3TRTdQgicarNJwUzy2lagD3APOBRoCg41i4zm2Nme8xsdYtjd5nZLjNbHiwXtXjtO2a2ycw2mNn5R3ZZcjh2l1Rz9YMLaXTn4RtnqJWRSBxr70lhCaF6hNaGvnRgTAc/+8/Ar4GHDjn+C3e/p+UBMzsOmE1oNrehwN/MbLy7N3TwGXKEDlTUcs0fQ0NXzLtppjqmicS5NpOCux/RSGfu/raZjQrz9FnAfHevAbaa2SZgBvDekcQg7auqbeCGuR+wfX8lc6+foaErRCSsimbM7FIzuydYLj7Cz/y6ma0Mipf6B8eGATtbnJMfHGstlpvNbLGZLS4u1lw/ndXQ6Nw2fxnLdh7kvi9M1eB2IgKE16P5J8DtwNpgud3M/reTn/c7YCwwFSgE/m/Tx7RybquD7rn7/e4+3d2n5+XldTKM+Obu3LlgNa+tLeKuSyZy4aQh0Q5JRHqIcDqvXQRMdfdGADObCywDvnO4H+buRU3bZvYA8EKwmw+MaHHqcKDgcH++hOf3b23h4fd38NVPjeXaU0dFOxwR6UHCKj4Csltsd7rg2cxa/kl6GdDUMmkBMNvMUs1sNDAOWNTZz5G2vbSqkJ++vJ5Lpgzl/5x/TLTDEZEeJpwnhf8FlpnZG4SKec4kjKcEM5sHnAXkmlk+cCdwlplNJVQ0tA34CoC7rzGzxwkVT9UDt6rlUddbvvMg33hsOdNGZvPzKyZrTmUR+QRz73i+nOAv/JMIJYWF7r470oGFY/r06b548eJohxETdh2sYtav3yU9JYFnv3YaAzJTox2SiESJmS1x9+mtvdbmk0Iwh0JL+cF6qJkNdfelXRWgRFZlbT03zV1MTV0D8246WQlBRNrUXvHRYmAN0NTus2VZgwOfjlRQ0nXcnW89uZJ1u0uZc91JjBukzmki0rb2ksJ/AJcDVcB84Bl3L++WqKTL/PbNzby4spA7LpzA2ccMjHY4ItLDtdn6yN1/4e6nA18n1Fz0dTN7PKgolhjw9/VF3PPqBi6dMpSvnNnRqCQiImE0SXX3rcBzwKuEhp4YH+mg5Mht31fB7fOXc9yQLH56+WRCs6qKiLSvvYrmMYQGqZtFaAiK+cCP3L26m2KTTqqqbeCrDy8lwYzf/+uJpKckRjskEYkR7dUpbAJWEnpKKAVGAl9r+ovT3e+NeHRy2Nyd7z+7mvW7S5lz7UkaBltEDkt7SeGHfDT+UGY3xCJdYN6inTy1NJ/bzhnH2RNUsSwih6e9obPv6sY4pAusKyzlrufXcMa4XG4/Z1y0wxGRGBTu2EfSw1XU1HPro0vJTk/mF1+YSqKGsBCRTghn7CPp4ZrqEbbtreCRG2eSqx7LItJJ7c3RfHuwPq37wpHOeHJJPs8s28Vt54zTZDkickTaKz66Plj/qjsCkc7ZtreCOxesYeaYHP7t06pHEJEj017x0Toz2wbkmdnKFscNcHefHNHIpEN1DY3c/thykhMTuPfzqkcQkSPXXuujq8xsMPAKcGn3hSTh+uXrH7Ji50F+88VpDM1Oj3Y4ItILtFvRHMybMMXMUvhoeIsN7l4X8cikXR9s289v3tjElScO57OTNceyiHSNDlsfmdmngIcIzZRmwAgzu9bd345wbNKG8pp6vvn4cob378Odl06Mdjgi0ouE0yT1XuA8d98AYGbjgXnAiZEMTNr245fWkX+giie+cgqZqWpVLCJdJ5zOa8lNCQHA3TcCyZELSdrz1sZiHl24g5vPGMP0UTnRDkdEeplw/sxcbGZ/BP4S7F8NLIlcSNKWkqo6vv3kSsYNzOQbn9EI5iLS9cJJCrcAtwK3EapTeBv4bSSDktb98Pm1FJfXcP+XTiQtWcNhi0jX6zApuHsNoXoFDZUdRW9u2MNTS/P5+tlHM3l4drTDEZFeSgPixYDymnq+98xqxuZl8G/nHB3tcESkF1PTlRjws5fXU1BSxZNfPZXUJBUbiUjktPukYGaJZvbz7gpGPmnR1v089N52rjt1FCce1T/a4YhIL9duUnD3BuBE06zvUVFT38AdT69keP90/vO8Y6IdjojEgXCKj5YBz5nZE0BF00F3fzpiUQkAf3hrC1uKK/jz9SeRoU5qItINwrnT5AD7gE+3OOaAkkIEbd1bwa/f2MRnJw/hrGM017KIdI9wmqRe39E50rXcnf96djWpiQncefFx0Q5HROJIh01SzWy8mb1uZquD/clm9v3Ihxa/Fqwo4J1Ne/nWBccwMCst2uGISBwJp5/CA8B3gDoAd18JzI5kUPGsrLqOu19Yx5Th/bj65KOiHY6IxJlwkkIfd190yLH6SAQjcN/fPmRfRQ13f+54zaQmIt0unKSw18zGEqpcxsyuAAo7epOZzTGzPU3FTsGxHDN7zcw+DNb9g+NmZr80s01mttLMpnXyemLaxqIy/vTPbcw+aaSGshCRqAgnKdwK/AGYYGa7gH8HvhrG+/4MXHDIsTuA1919HPB6sA9wITAuWG4GfhfGz+9V3J27FqwhMzWJb52vPgkiEh0dJgV33+Lu5wJ5wAR3P93dt4fxvreB/YccngXMDbbnAp9rcfwhD3kfyDazuJpj8qVVu/nn5n385/nHkJOREu1wRCROhdP6aICZ/RL4B/Cmmd1nZgM6+XmD3L0QIFg3NcAfBuxscV5+cCwuVNU28KMX13LckCy+OGNktMMRkTgWTvHRfKAYuBy4Ith+rIvjaK1G1Vs90exmM1tsZouLi4u7OIzoeOAfWygoqebOS45T5bKIRFU4SSHH3e92963B8j9AZ2tBi5qKhYL1nuB4PjCixXnDgYLWfoC73+/u0919el5eXifD6DmKSqv53ZubufD4wZw8prMPYCIiXSOcpPCGmc02s4Rg+TzwYic/bwFwbbB9LfBci+NfClohzQRKmoqZerufv7KBhkbnOxceG+1QRETaHubCzMoIFeEY8E3g4eClBKAcuLO9H2xm84CzgFwzyw/O/wnwuJndAOwArgxOfwm4CNgEVAJxMbTGqvwSnlqaz81njGHkgD7RDkdEpO2k4O59j+QHu/tVbbx0TivnOqGmr3HD3bn7xbXk9Enh1k9rNjUR6RnCGo/ZzCYDo1qer6Gzj8zr6/awaOt+7p41kay05GiHIyIChJEUzGwOMBlYAzQGhzV09hGob2jkpy+vZ0xuBrPVBFVEepBwnhRmurvGb+5CTy3N58M95fzu6mkkJ4ZT1y8i0j3CuSO9Z2ZKCl2kqraBe1/byAkjs7ng+MHRDkdE5GPCeVKYSygx7AZqCLVGcnefHNHIeqk5726lqLSGX101DU19LSI9TThJYQ5wDbCKj+oUpBNKKuv4/VubOffYgcwYnRPtcEREPiGcpLDD3RdEPJI48Ie3N1NeU89/ahRUEemhwkkK683sUeB5QsVHgJqkHq7ishr+9O42Lp48lAmDs6IdjohIq8JJCumEksF5LY6pSeph+u2bm6htaOQb546LdigiIm3qMCm4e1wMORFJBQereOT9HVw+bRhj8jKjHY6ISJvC6bz2J1oZxtrdvxyRiHqhX/19E45z2zl6ShCRni2c4qMXWmynAZfRxrDW8kk791fyxOKdXDVjJMP7a9A7EenZwik+eqrlfjD66d8iFlEv89s3N5FgxtfOHhvtUEREOtSZMRbGARqwJwz5Byp5YnE+XzhpBEP6pUc7HBGRDoVTp9A0r0KT3cC3IxZRL/LbNzdjBrecpacEEYkN4RQfHdG8CvFq18Eqnli8k89PH8HQbD0liEhsaG/mtXaLiNx9R9eH03v87s1NAHztbE2gIyKxo70nhRf5aDrOJg7kAQOBxAjGFdOKSqt5/IN8rjhxBMP0lCAiMaS96Tgntdw3s1GE6hLOBX4c0ahi3IP/2EKDO7d8SnUJIhJbOmx9ZGbjzOzPwF+BJcBx7v6rSAcWqw5W1vLIwh1cMnkIIweoX4KIxJb26hSOB74HTAR+Btzg7g3dFVis+vM/t1FZ28AtZ6kuQURiT3t1CiuAnYTqFmYAM1pOCuPut0U2tNhTUVPPn/+5jXOPHcQxg9VoS0RiT3tJQWMbHaZ5i3ZwsLJOvZdFJGa1V9E8tzsDiXU19Q088I8tnDJmANNG9o92OCIindKZYS6kFQuWF1BUWsNX1XtZRGKYkkIXcHce/MdWJgzuy5njcqMdjohIpykpdIG3P9zLhqIybjpjDC0r40VEYk04/RTGm9nrZrY62J9sZt+PfGix44G3tzAoK5VLpgyNdigiIkcknCeFB4DvAHUA7r4SmB3JoGLJmoIS3tm0l+tPG01Kkh68RCS2hXMX6+Puiw45Vh+JYGLRg//YSkZKIlfN0BQTIhL7wkkKe81sLMGcCmZ2BVAY0ahiRGFJFc+vKOALJ42kX3pytMMRETli4czRfCtwPzDBzHYBW4F/jWhUMeIv722n0Z3rTxsV7VBERLpEOJPsbAHONbMMIMHdyyIfVs9XXdfAvEU7OPfYQYzI0cB3ItI7hDMdZypwOTAKSGpqcunuP+zsh5rZNqAMaADq3X26meUAjwWfsw34vLsf6OxnRNqzy3ZxoLKO608bHe1QRES6TDh1Cs8BswhVLle0WI7U2e4+1d2nB/t3AK+7+zjg9WC/R3J3/vTuNiYM7svMMTnRDkdEpMuEU6cw3N0viHgkocRzVrA9F3iT0KQ+Pc57m/exoaiMn10+WZ3VRKRXCedJ4Z9mNqnj0w6LA6+a2RIzuzk4NsjdCwGC9cDW3mhmN5vZYjNbXFxc3MVhhWfOu9vIyUjh0qnqrCYivUt7k+ysBhqDc643sy1ADaE5m93dJx/B557m7gVmNhB4zczWh/tGd7+fUGsopk+f7kcQQ6fs3F/J6+uLuPWso0lL1jTVItK7tFd8NAyYGokPdfeCYL3HzJ4hNIlPkZkNcfdCMxsC7InEZx+phxduJ8GMq2eqs5qI9D7tFR9tdfftbS2d/UAzyzCzvk3bwHnAamABcG1w2rWEKrh7lOq6Bh7/YCefOXYQQ/qlRzscEZEu196TwkAz+2ZbL7r7vZ38zEHAM0EFbRLwqLu/bGYfAI+b2Q3ADuDKTv78iHlpVSEHKuu45pSjoh2KiEhEtJcUEoFMQnUIXSboDDelleP7gHO68rO62l/e386YvAxOHTsg2qGIiEREe0mh8Eg6qPU2q3eVsGzHQX5w8XFqhioivVZ7dQq687Xw8PvbSUtO4PITh0c7FBGRiGkvKfToopzuVFJVx7PLd/G5qcM0GqqI9GptJgV339+dgfRkzy3fRXVdI1efrApmEendNFVYB9ydRxfu4PhhWUwa3i/a4YiIRJSSQgdW5JewfncZs09SZzUR6f2UFDowf9EO0pMTmaVxjkQkDigptKO8pp4FKwq4ZMoQ+qapgllEej8lhXYsWF5AZW0Ds2eo6EhE4oOSQjvmf7CDCYP7csKI7GiHIiLSLZQU2rCmoISV+SXMPmmEejCLSNxQUmjDE4vzSUlM4HMnDIt2KCIi3UZJoRW19Y08t3wXn5k4iOw+KdEOR0Sk2ygptOL1dUUcqKzjSo1zJCJxRkmhFY8v3sngrDTOGJcX7VBERLqVksIhikqreWtjMf8ybRiJCapgFpH4oqRwiKeX7qLR4QoVHYlIHFJSaMHdeWLJTqYf1Z8xeZnRDkdEpNspKbSwbOdBthRXcOV0PSWISHxSUmjhmaW7SE1K4KJJQ6IdiohIVCgpBGrrG3lhZQHnTRyswe9EJG4pKQTe2ljMgco6LjtBQ2SLSPxSUgg8syyfARkp6psgInFNSQEoqarjb+v2cMmUoSQn6p9EROKX7oDAX1cVUlvfyGUa/E5E4pySAvD0sl2Myctg8vB+0Q5FRCSq4j4p5B+oZNHW/Vw2dZjmTRCRuJGriNEAAAafSURBVBf3SeH5FYUAmjdBRAQlBRasKOCEkdmMyOkT7VBERKIurpPCpj1lrCss5dIp6psgIgJxnhQWrCgkweCzGtZCRASI46Tg7jy/ooCZYwYwMCst2uGIiPQIcZsU1hSUsnVvBZeo6EhEpFmPSwpmdoGZbTCzTWZ2R6Q+Z8GKApISjAuPHxypjxARiTk9KimYWSLwG+BC4DjgKjM7rqs/p7HReWFFAWeOzyO7T0pX/3gRkZjVo5ICMAPY5O5b3L0WmA/M6uoPWbLjAAUl1Wp1JCJyiJ6WFIYBO1vs5wfHmpnZzWa22MwWFxcXd+pDDDhzfB7nHjeo04GKiPRGPS0ptDbOhH9sx/1+d5/u7tPz8jo3zPX0UTk89OUZZKYmder9IiK9VU9LCvnAiBb7w4GCKMUiIhJ3elpS+AAYZ2ajzSwFmA0siHJMIiJxo0eVn7h7vZl9HXgFSATmuPuaKIclIhI3elRSAHD3l4CXoh2HiEg86mnFRyIiEkVKCiIi0kxJQUREmikpiIhIM3P3js/qocysGNjeybfnAnu7MJxo0rX0TL3lWnrLdYCupclR7t5q79+YTgpHwswWu/v0aMfRFXQtPVNvuZbech2gawmHio9ERKSZkoKIiDSL56Rwf7QD6EK6lp6pt1xLb7kO0LV0KG7rFERE5JPi+UlBREQOoaQgIiLNen1SMLMLzGyDmW0ysztaeT3VzB4LXl9oZqO6P8rwhHEt15lZsZktD5YboxFnR8xsjpntMbPVbbxuZvbL4DpXmtm07o4xXGFcy1lmVtLiO/lBd8cYDjMbYWZvmNk6M1tjZre3ck5MfC9hXkusfC9pZrbIzFYE1/LfrZzTtfcwd++1C6HhtzcDY4AUYAVw3CHnfA34fbA9G3gs2nEfwbVcB/w62rGGcS1nAtOA1W28fhHwV0Iz8c0EFkY75iO4lrOAF6IdZxjXMQSYFmz3BTa28vsVE99LmNcSK9+LAZnBdjKwEJh5yDldeg/r7U8KM4BN7r7F3WuB+cCsQ86ZBcwNtp8EzjGz1qYFjbZwriUmuPvbwP52TpkFPOQh7wPZZjake6I7PGFcS0xw90J3XxpslwHrOGR+dGLkewnzWmJC8G9dHuwmB8uhrYO69B7W25PCMGBni/18PvnL0XyOu9cDJcCAbonu8IRzLQCXB4/2T5rZiFZejwXhXmusOCV4/P+rmU2MdjAdCYofTiD0V2lLMfe9tHMtECPfi5klmtlyYA/wmru3+b10xT2styeF1rLloVk2nHN6gnDifB4Y5e6Tgb/x0V8PsSZWvpNwLCU0zswU4FfAs1GOp11mlgk8Bfy7u5ce+nIrb+mx30sH1xIz34u7N7j7VEJz1s8ws+MPOaVLv5fenhTygZZ/LQ8HCto6x8ySgH70zOKADq/F3fe5e02w+wBwYjfF1tXC+d5igruXNj3+e2hWwWQzy41yWK0ys2RCN9FH3P3pVk6Jme+lo2uJpe+libsfBN4ELjjkpS69h/X2pPABMM7MRptZCqFKmAWHnLMAuDbYvgL4uwc1Nj1Mh9dySPnupYTKUmPRAuBLQWuXmUCJuxdGO6jOMLPBTeW7ZjaD0P+5fdGN6pOCGP8IrHP3e9s4LSa+l3CuJYa+lzwzyw6204FzgfWHnNal97AeN0dzV3L3ejP7OvAKodY7c9x9jZn9EFjs7gsI/fL8xcw2Ecqus6MXcdvCvJbbzOxSoJ7QtVwXtYDbYWbzCLX+yDWzfOBOQhVouPvvCc3RfRGwCagEro9OpB0L41quAG4xs3qgCpjdQ//oOA24BlgVlF8DfBcYCTH3vYRzLbHyvQwB5ppZIqHE9bi7vxDJe5iGuRARkWa9vfhIREQOg5KCiIg0U1IQEZFmSgoiItJMSUFERJr16iapIl3FzAYArwe7g4EGoDjYr3T3U6MSmEgXU5NUkcNkZncB5e5+T7RjEelqKj4SOUJmVh6szzKzt8zscTPbaGY/MbOrg/HwV5nZ2OC8PDN7ysw+CJbTonsFIh9RUhDpWlOA24FJhHrVjnf3GcCDwL8F59wH/MLdTwIuD14T6RFUpyDStT5oGg/IzDYDrwbHVwFnB9vnAse1GPI+y8z6BmP/i0SVkoJI16ppsd3YYr+Rj/6/JQCnuHtVdwYmEg4VH4l0v1eBrzftmNnUKMYi8jFKCiLd7zZgejBD3lrgq9EOSKSJmqSKiEgzPSmIiEgzJQUREWmmpCAiIs2UFEREpJmSgoiINFNSEBGRZkoKIiLS7P8Dnpf16SBaFmkAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret1 = run_simulation(duration, y0=y0, volume=V, model=m)\n", "ret1" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Simulating with `gillespie` (Bars represent standard error of the mean):" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEGCAYAAACKB4k+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXxcdb3/8dcne9q0pPsSWlqgCy0FWmoBiwpWLdYFRJSiIqJXuCoXXH4oVbx6RS8orrig5YqCylJoadkrm7JIwe57oXRJky5J02Zr9pnP74+ZhLTNctrMZDKZ9/PxmEfmbHM+hynnM+e7mrsjIiICkJboAEREpOdQUhARkRZKCiIi0kJJQUREWigpiIhIi4xEB9AVgwcP9jFjxiQ6DBGRpLJixYr97j6krW1JnRTGjBnD8uXLEx2GiEhSMbOd7W1T8ZGIiLRQUhARkRZKCiIi0kJJQUREWigpiIhIi7glBTMbZWYvmNkmM9tgZjdE13/fzIrNbHX0NafVMfPMbKuZbTGz2fGKTURE2hbPJqlNwDfcfaWZ9QNWmNkz0W2/cPeftt7ZzCYBc4HJwEjgWTMb7+6hOMYoIiKtxO1Jwd33uPvK6PsqYBNQ0MEhFwMPuHu9u28HtgIz4hWfiIgcrVvqFMxsDDAVeC266jozW2tmd5vZgOi6AmBXq8OK6DiJiIgkzOV/eJXL//BqosOIubgnBTPLAxYCX3X3SuBO4BTgLGAP8LPmXds4/KgZgMzsGjNbbmbLS0tL4xS1iEhqimtSMLNMIgnhb+6+CMDd97l7yN3DwF28XURUBIxqdfiJwO4jP9Pd57v7dHefPmRIm0N3iIj0eF150ojnU0o8Wx8Z8Edgk7v/vNX6Ea12+xiwPvr+UWCumWWb2VhgHPB6vOITEZGjxbP10UzgSmCdma2Orvs2cIWZnUWkaGgHcC2Au28wswXARiItl76ilkciIt0rbknB3V+m7XqCJzs45kfAj+IVk4j0Ps3FKA9ee16CI+kd1KNZRERaKCmISMrqrc1Ku0JJQUREWigpiEjCJdsv9sWrillVWM5r2w8w87bnWbyqOObHuzuNoTBVdY3sr66n6GANb5VW89sXtrJy58HjPndnkno6ThGR7rZ4VTHzFq2jIRQGoLi8lm8tXMvuilrOGTuImoYmDtWHqG1soqYhRE19iNrGEDUNIWobmti8t5LlO8oJubcc/7UFq/nRkxvJSk+ntjFEXfQVPqr77uGKy2uZt2gdAJdMjc0AEEoKIgJ0rRVPsrUAqmsMUVnXSG1DiKaw8883Sqmqa6SqrqnV3yaq65uojv6tqm/iUH0T20qrj7pZ1zeF+cnTWzo8Z1Z6GjmZaVTXNx11vDtU1zUxZ8pQcjLTyMlMj/zNSCcnM53szDSyM9K47anNHKxpPOzY2sYQty/doqQgIqnN3Xlo+S5WFh6kMeScfcszzJkynFOH9qO8ppHy2gYqahopr22k4ohXQ1P4sM+66u7D+8mmGfTNzqB/TiZ9s9Oj7zMoyM9ha0l1uzHd8/kZ9MlKJzczckyfrHRyo8uZ6ZHS+rE3PdHmsXWNYX72yTM7vOabFq5rc/3u8toOjzsWSgoivUiy/WI/UkNTmLJD9ZRW1bO/up79VQ3sP1RPWXUDZdX1lB1qoKy6gYM1DZRU1hFq9Yu77FADf1lW2LLcLzuDE/pkkt8nkxNyMxnWP48TcjPpn5tJ0YEa/r5xH42tPiA7I415cyZy2dmj6JuVTmRQhqPNvO15itu4CRfk5/Ke8Z0PvTMyP7fN40fm58b12KCUFEQk7hqawpRU1bGvso59lfXsq6xjb2UdpZX17KuqY21RBY2hMONvfqrN4/tkpTOwbxaD8rIZfkIOk0b25+n1e6mubzpq3+En5PDSNy9s+WXelpm3PX9YQoBIEdBdL27nc+8c2+G13Dh7AvMWraO28e0BF3Iz07lx9oQOj4vF8V09dxBKCiLSJeGw09AUpqEpzBNr91BcXsPu8jp2l9eyp6KOPRV17K+uP+q4rPQ0hvTLZmj/bHIy0+iXk8Fnzj2JwXnZDOmXzeC8LAbnZTMoL4s+WUffqhauKGoznn0VdR0mBGi/uCVIMUxz2f03H15LQyhMQX4uN86eELhMvyvHd/XcQSgpiEinKmobKSyrofBA5LXrYA27DtRQdLCW4vLaljL6r9y3EoC87AxG5ucw4oRcJo/sz/ATchjeP4dhzX/75zCgT2ZLEU1zsdf1s8YFjimRxTCXTC3g/tcjRVXHU1TXleO7eu7OKCmICABNYaeuMcSS1cVsKz3EjrJD7CirYWfZIcqPaPEyoE8mowb2YdKI/owd1Id/vrmfUNgZ2i+br71/HFfMOCnweZvb7DeEwsy87fnAv3x7ejFMslJSEEkh7k5JVT1v7Ktia0k1W0uqeau0mg3FFVTVR26QNzwQGdS4ID+XsYP78qEpIzhpUB9GD+zL6IF9GD2oD3nZkVtHc5v9ULSNZUlVPT94bBO5mRmBbuxttfkP2u6+pxfDJCslBZFeqqK2kc17Ktmyr4rNe6t4Y28Vb+yrorLu7crZ/jkZDOiTSU3D4U00czLTAt0kb1+65bBf23Bs7ea7enxPLoZJVkoKIj3MsTYrdXeKy2u5++XtLN9xgJDDKfOebOkxC5Gb/8Th/fnImSMZP6wf44bmceqwPIbkZXP+j184bF+ItJkPcmPuSoVtLI6X2FNSEEkyeyvqWL2rnHXF5awtqmBdccVRZf4hdzLSjM+fP4arZ45leP+cdtvdd+XG3NUK2+5od99T9dSnEyUFkR6stiHE2qJyVu0qZ+XOg6wpKmdfZaR5Z0aaMX5YP2ZPGs7TG/ZSUXt4YmgKO0+s3cu350zq8BxduTEnss1+KotnQlFSEOlBSqsivXar6xr56G9eZuPuSpqilbhjBvXhvJMHceaofM4clc+kEf3JyUwHYMHyXW1+XpBf+125MSeyzX4sdOXm2lN/6XeVkoJIApVU1vHqtjKWbTvAa9vL2FZ6qGXbgUONXDBhCFfMGM3U0QMY2Der3c/pyq/9WNzYE9VmX2JPSUEkDtqrLK6oaeTVbWX86639vLx1f0sS6JeTwagBuWSkWcuTQUMozCtby/jwGSM7TAjQ9WIY3ZilmZKCSByFws7qXeW8+EYpL75Zyppd5YQ9MpbPjLEDmfuOUZx38mAmjezPu3/yQktCaBa0eWaii2ESSUkstpQURNrQldFGDx5qYH91PeU1jZz9w2cor2kkzeCME/O57sJTOX/cEM4alU9WxuHj83S1eaZ+7UssKCmIxMCO/Yf4+8a9PLuxhOU7DxD2SOugi88q4IIJQ3jXuMHk9+m4CCiVm2dKz6GkIHIc3J1Ne6p4ev0ent6wlzf2RSZeOW1Ef95/2lCe31JKY8hZtq0sUEKA5G6eqSeT3kNJQSQgd2fz3ioeW7Obp9bvZfv+Q6QZTB8zkO9+eBIfmDSMFTsPMm/Rupax+rtrLJ9Y0I1dQElBpFOFZTUsXl3Mo2t2s7WkmvQ047yTB/Ef7xrL7MnDGZyX3bLv3PnLEjaWj0gsKCmItKEpFKbsUAMfv/NfrNh5EIAZYwZyyyWnM+f04QxqlQha01g+x0cJsOdQUhCJCoWdl94s5ZfPvsHqXRUAFB2s5cNnDGfenEkU9JA5dDuim6t0lZKC9FpBm5UWl9ey4N+7eGj5LnZX1B22rSnsPLeplPeddoCCOE/8ItITKClISgqHnRffLOWvy3by/OYSHDj/1MHUNoY4eMSIo8daJwBdqyzWr31JJCUFSSkVtY08tHwX9766k8IDNQzOy+bLF5zK5e8YxaiBfRh70xNtHncsdQKqLJZkFrekYGajgHuB4UAYmO/uvzKzgcCDwBhgB/BJdz9okcHefwXMAWqAz7n7ynjFJ6ll+/5D/OmV7Ty8ooiahhDvGDOA/zd7AhdNHn5Yz+JE1wmIJFqnScHMPgE87e5VZnYzMA34YYAbdhPwDXdfaWb9gBVm9gzwOeA5d7/NzG4CbgK+BXwQGBd9nQPcGf0rclzcner6Jq65dznPbNpHZloaHzlzJFfPHMPpBSe0eYzqBCTVBXlS+K67P2Rm5wOzgZ8S4Ibt7nuAPdH3VWa2CSgALgYuiO52D/APIknhYuBed3dgmZnlm9mI6OdIijqeMYjCYef5zSVs3FNFdX0TuyvquO7CU7nyvJMY2i+nw2MT3YFMJNGCJIXmn0wfAu509yVm9v1jOYmZjQGmAq8Bw5pv9O6+x8yGRncrAFrPFFIUXXdYUjCza4BrAEaPHn0sYUgvFwo7j6/dze9eeIst+6pa1udmpnPKkLxOE0KzWNQJqC5BklWQpFBsZn8A3gf82MyygbROjmlhZnnAQuCr7l7Z3jyxQFsb/KgV7vOB+QDTp08/aruknlDYeWzNbn79/Ju8VXqIYf2yyUy3lqEm9lTUBR5qQiTVBbm5fxJYClzk7uXAQODGIB9uZplEEsLf3H1RdPU+MxsR3T4CKImuLwJGtTr8RGB3kPNIagpHnww+8It/8tUHV5OZnsbvPj2N9FYJoVlzs1IR6VinScHda4jcuM+PrmoC3uzsuGhroj8Cm9z95602PQpcFX1/FbCk1frPWsS5QIXqE6Qt7s7zm/fxoV+/zHX3rSLNjDs/PY0nr38Xc6aMYE95XZvHaagJkc4FaX30PWA6MAH4E5AJ/BWY2cmhM4ErgXVmtjq67tvAbcACM/sCUAh8IrrtSSLNUbcSaZJ69TFdiaSEVYUHufWpzby+/QAnDerDLy8/i4+cOZL0tLdLH9WsVOT4BalT+BiRSuKVAO6+O9rEtEPu/jJt1xMAzGpjfwe+EiAeSUG7DtRw29ObeWLtHgbnZXHLxZOZO2M0melHP+yqWanI8QuSFBrc3c3MAcysb5xjEmnRFHZ2l9cy62f/JD3NuGHWOL747pPJy27/n66GmhA5fkGSwoJo66N8M/si8HngrviGJakuHHYeWVXM2qJyGkPOx6edyI2zJzD8hO5rViqSijpNCu7+UzN7P1BJpF7hv939mbhHJilrfXEF/71kPSsLy1vKH5dtK2PZtjI1KRWJs0BjH0WTgBKBxFV1fRM/+/sW7vnXDnIz0w/ra3As01qKyPFrt0mqmVWZWWWrv5Wtl7szSOn9nl6/l/f97J/8+V87+NQ5o+mfm6m+BiIJ0O6Tgrt32sJIJIiOxi8qqarje0s28NT6vZw2oj93fmYaU0cP4G/Luj6EtYgcuyD9FNocYMjdC2MfjqQKd2fhymJueXwjtY0hvnnRBL74rpNbmpiqr4FIYgSpU2j9ky0HGAtsASbHJSLp9Uqq6vj2onU8u6mE6ScN4MeXncEpQ/IO2ycWfQ3U6kjk2AVpfTSl9bKZTQOujVtE0qs9tmY3312yntqGEN/98CSufucY0tKO7uOoIaxFEuOYZ16LTprzjngEI71XUzjM1x5czSOrijlzVD4/+8SZnDo0r8Nj1NdApPsFqVP4eqvFNCIzr5XGLSLpdarqGnmr9BCrCsu5YdY4/uu9p5LRxvAUIpJ4QZ4UWrdCaiJSx7AwPuFIbxIOO7/7x1Y27qkiOyONh/7zPM4+aWCiwxKRDgSpU/if7ghEepfSqnq+vmA1L725nzSD+qYw19+/WvUCIj1cp8/wZvaMmeW3Wh5gZkvjG5Yks9e2lTHnjpf419b9ZKYb4WgftOZeyYtXFSc2QBFpV5CC3SHRGdcAcPeDwNAO9pcU5e7c9eI2PvV/r9EvO4NBednqlSySZIIkhVDrDmxmdhJtzJ0sqa26vomv3LeSHz25ifefNowl182ktKq+zX3VK1mk5wpS0fwd4GUz+2d0+d3ANfELSZJNYVkNX7x3OW+WVDHvgxO55t0nY2Yx6ZWspqgi3SvIHM1PE2mG+iCwADjb3VWnIAC8snU/H/3ty+ytrOPez5/Dte85hcj03JFeybmZ6YftrxnQRHq2IBXNBlwETHP3x4A+ZjYj7pFJj3L5H15tGdiu2V+W7eSzd7/O0H7ZPHrdTM4fN/iw7ZdMLeDWS6eQFe2TUJCfy62XTlHrI5EeLEjx0e+AMPBe4AdAFZF+CurVnKJCYeeHT2zkT6/s4L0Th3LHFVPbnR5TvZJFkkuQpHCOu08zs1UQaX1kZllxjkt6qEP1TVx//yqe21zC52eO5TsfOo30NsYuEpHkFCQpNJpZOtEWR2Y2hMiTg6SYxlCYufOXsWF3BbdcPJkrzxuT6JBEJMaCJIU7gEeAoWb2I+Ay4Oa4RiU9Tl1jiM17qzCDuz47nVmnDUt0SCISB0GGufibma0AZgEGXOLum+IemfQYdzz3BmuKKgAYkpdNVV1TgiMSkXhpNymYWeuRy0qA+1tvc/cD8QxMeobbl27mty+81bJcWl3PvEXrAAK3IlIFs0jy6OhJYQWReoS2ahEdODkuEUmP8czGfYclhGbNQ1WoaalI79NuUnD3sd0ZiPQsS1YX8/UFa9rdrqEqRHqnQDOvmdlHiQxvAfAPd388fiFJoj28oogbH17DjDEDKTxQw56KuqP2OZahKkQkeQTp0XwbcAOwMfq6wcxujXdgkhj3v17IjQ+vYeYpg/nz1TP41kUTNVSFSAoJ8qQwBzjL3cMAZnYPsAqYF8/ApPv9ddlObl68ngsnDOHOz5xNTmZ6S73BNx9eS0MoTEF+ribKEenFgk6Um9/q/QlBDjCzu82sxMzWt1r3fTMrNrPV0decVtvmmdlWM9tiZrMDxiUx8sDrhdy8eD2zJg7l91dGEkKzS6YWMHV0PueMHcgrN71XCUGkFwvypHArsMrMXiDSEundBHtK+DPwG+DeI9b/wt1/2nqFmU0C5gKTgZHAs2Y23t1DAc4jXfTQ8l3Me2QdF0wYwu8+M43sjPTODxKRXilI57X7zewfRAbAM+Bb7r43wHEvmtmYgHFcDDzg7vXAdjPbCswAXu34MOmqJauL+ebCtZx/6mB+/5mzlRBEUly7xUdmNq35BYwAioBdwMjouuN1nZmtjRYvDYiuK4h+drOi6Lq24rrGzJab2fLS0tIuhCHPbtzH1xdEWhnNv3L6YUVGIpKaOnpSWA5sAJrvvK07sTmRobSP1Z3ALdHjbwF+Bnye9jvIHb3SfT4wH2D69OmaFvQ4/eut/Xz5vpWcPrI/f/zcO8jNUkIQkY6TwjeAjwO1wAPAI+5e3ZWTufu+5vdmdhfQ3N+hCBjVatcTgd1dOZe0b82ucr54z3LGDOrDn6+e0e5cCK1pqAqR1NBu8ZG7/8LdzweuI3LDfs7MFpjZWcd7MjMb0WrxY0Bzy6RHgblmlm1mY4FxwOvHex45WvPMadv3H+LqP/+bAX2z+MsXzmFAX02NISJvC1LRvN3MlgC5wJXAeGB1Z8eZ2f3ABcBgMysCvgdcEE0qDuwAro2eY4OZLSDSOa4J+IpaHsVeQ1OYz979GgD3fn4Gw/rnJDgiEelpOhol9WQizUQvJlIJ/ADwI3c/esyDNrj7FW2s/mMH+/8I+FGQz5Zj1xR2tuyrwh3uv+ZcTh6Sl+iQRKQH6uhJYSuwFlgCVAKjgS+bReqE3f3ncY9OYmLRil2sLDyIOwzqm8WO/Yc4a1R+5weKSMrpKCn8gLdbAOlnZZJ6ZGURNy5ch0e/ybJDDcc8H4KIpI6Ohs7+fjfGIXHy/cc2EAof3nJX8yGISHuCjn0kSejvG/ZSUdv21JmaD0FE2qKk0Ett3lvJVx9cTWZ6W/0CNR+CiLSto2Euboj+ndl94UgslFXX8x/3LKdfTgY3f+g0zYcgIoF19KRwdfTvr7sjEImNhqYwX/rbSkqq6pl/5XSueudYbr10Clnpka+6ID+XWy+dovoEEWlTR62PNpnZDmCIma1ttd4Ad/cz4hqZHJdbHt/I69sP8MvLz+LMaLPTS6YWcP/rhYCGqxCRjnXU+ugKMxsOLAU+2n0hyfFauKKIvyzbyRffNVZPAiJyXDoc5iI6b8KZZpZFZHgLgC3u3hj3yOSYbNhdwbcfWce5Jw/kWxdNPGq7nhBEJIhOxz4ys/cQmT1tB5Gio1FmdpW7vxjn2CSg8poG/vOvKxjQJ4vffGoaGelqVCYixyfIdJw/Bz7g7lsAzGw8cD9wdjwDk2DCYecbC9awt6KOB689j8F52YkOSUSSWJCflJnNCQHA3d8AMuMXkhyL/3t5G89tLuE7c05j2ugBnR8gItKBIE8Ky83sj8BfosufBlbELyQJasXOA/z46S188PThXPXOMYkOR0R6gSBJ4UvAV4DridQpvAj8Lp5BSdsu/8OrQKTS+OChBv7rvlUU5Ofy48vOoHn0WhGRrggyyU49kXoFDZXdQ7g7Nz68hv3VDSz80jvpn6PSPBGJDTVTSUJ/XbaTZzeV8K0PTmTKiSckOhwR6UWUFJJMTUMTP3xiE+8ZP4SrVY8gIjHWYVIws3Qzu727gpH2LV5VzKrCctYVV9IYCjPrtKGkpakeQURiq8Ok4O4h4GxTLWZCLV5VzLxF62gIhQEIO9z65GYWrypOcGQi0tsEaX20ClhiZg8Bh5pXuvuiuEUlh7l96RZqG0OHrdPsaSISD0GSwkCgDHhvq3UOKCl0k/ZmSdPsaSISa0GapF7d2T4SX7lZ6dQ0hI5ar9nTRCTWOm19ZGbjzew5M1sfXT7DzG6Of2gC8MzGfdQ0hMg4olJZs6eJSDwEaZJ6FzAPaARw97XA3HgGJRHlNQ18+5F1nDaiP7dp9jQR6QZB6hT6uPvrRzRAaopTPNLKLY9v4uChBv589TuYPPIEHlpRBGhuBBGJnyBPCvvN7BQilcuY2WXAnrhGJbywpYSFK4v40gWnMHmkei2LSPcI8qTwFWA+MNHMioHtREZKlTipqmvkO4vWMW5oHte999REhyMiKSRI66NtwPvMrC+Q5u5V8Q8rtd321Gb2Vtax8EvvJDsjvWW9io1EJN6CtD4aZGZ3AC8B/zCzX5nZoPiHlpqW7zjA314r5OqZY5mqSXNEpJsFqVN4ACgFPg5cFn3/YDyDSlUNTWG+/cg6CvJz+fr7xyc6HBFJQUGSwkB3v8Xdt0dfPwTyOzvIzO42s5Lm/g3RdQPN7BkzezP6d0B0vZnZHWa21czWmtm047+k5HXXS9t4Y181P7h4Mn2zg1T3iIjEVpCk8IKZzTWztOjrk8ATAY77M3DREetuAp5z93HAc9FlgA8C46Kva4A7gwTfm+zYf4g7nnuTD54+nFmnDUt0OCKSotpNCmZWZWaVwLXAfUBD9PUA8LXOPtjdXwQOHLH6YuCe6Pt7gEtarb/XI5YB+WY24lguJJm5O99dsp7M9DS+95HJiQ5HRFJYu2UU7t4vDucb5u57op+/x8yGRtcXALta7VcUXXdUfwgzu4bI0wSjR4+OQ4jd76n1e3npzf187yOTGH5CTqLDEZEUFqjg2szOAMa03j/GQ2e3NV+Dt7Wju88n0m+C6dOnt7lPMqlpaOKHj29k4vB+XHnuSYkOR0RSXKdJwczuBs4ANgDh6OrjHTp7n5mNiD4ljABKouuLgFGt9jsR2H0cn590fvP8VnZX1PGrK6aSka7ZUUUksYI8KZzr7pNidL5HgauA26J/l7Raf52ZPQCcA1Q0FzP1ZttKq7nrpW1cOrWAd4wZmOhwREQCJYVXzWySu288lg82s/uBC4DBZlYEfI9IMlhgZl8ACoFPRHd/EpgDbAVqgF47h8Plf3gVgAeuOZfvP7aRnIx0bpozMcFRiYhEBEkK9xBJDHuBeiLl/+7uZ3R0kLtf0c6mWW3s60TGWEoZL2wp4cU3Srn5Q6cxtJ8ql0WkZwiSFO4GrgTW8XadgnRB2J0fPr6Jkwf35bPnjUl0OCIiLYIkhUJ3fzTukaSQfZX1FB6o4e7PTScrQ5XLItJzBEkKm83sPuAxIsVHQMybpKaMxlCY4vJa3j1+CBdOGNr5ASIi3SjIz9RcIsngA8BHoq8PxzOo3mrxqmJW7yonFHY276lkyeqUaHUrIkkkyHwKvbYlUHdavKqYby1cSzja3a6kqp55i9YBaK5lEekxgnRe+xNt9C5298/HJaJe6valW6hvOryevrYxxO1LtygpiEiPEaRO4fFW73OAj5EivY1jaXd57TGtFxFJhCDFRwtbL0c7pT0bt4h6IXcnI91oDB09VNPI/NwERCQi0rbjaQ85Dugdw5N2k6fW76Ux5GSmHz7uX25mOjfOnpCgqEREjhZkjuYqM6tsfhFpmvqt+IfWOzSGwvzk6c1MGNaPH196BlnRQe8K8nO59dIpqk8QkR4lSPFRPOZVSBkPLS9iR1kNf7xqOrNOG8aDyyPTRjx47XkJjkxE5GjtJgUz67CIyN0LYx9O71LXGOKO597k7JMG8N6J6qgmIj1fR08KTxBpitq6INyBIcBQID2OcfUKf122k72Vdfxy7lmYtTWPkIhIz9LRdJxTWi+b2RgidQnvA/43rlH1AlV1jfz2ha28a9xgzj15UMt6FRuJSE8WpKJ5nJn9GXgKWAFMcvdfxzuwZHf3yzs4WNOo1kUiklQ6qlM4HfgOMBn4CfAFdw91V2DJrLymgbte2sbsycM448T8RIcjIhJYR3UKa4BdROoWZgAzWpeLu/v18Q0tef3x5e1U1zfxtfePT3QoIiLHpKOkoLGNjkN5TQN/emUHc6YMZ+Lw/okOR0TkmHRU0XxPdwbSWzQ/JVw/a1yiQxEROWaa9iuG9JQgIslOSSGG9JQgIslOSSFGKmoa9ZQgIkkvSD+F8Wb2nJmtjy6fYWY3xz+05HLvqzuorm/iugv1lCAiySvIk8JdwDygEcDd1wJz4xlUsjlU38Tdr2xn1sShTBqppwQRSV5BkkIfd3/9iHVN8QgmWd3/eiEHaxr58oWnJjoUEZEuCZIU9pvZKUTnaTazy4A9cY0qidQ3hbjrpW2cd/Igzj5pQKLDERHpkiBzNH8FmA9MNLNiYDvwmbhGlUQWrihmX2U9P/vEWYkORUSky4JMsrMNeJ+Z9QXS3L0q/mElh6ZQmN//8y3OHJXPzFMHdX6AiEgP12lSMLNs4OPAGCCjefwjd/9BXCNLAk+t30vhgRq+86HTNF+CiPQKQYqPlgAVRIbNro9vOMnD3Zn/4jZOHuCwtgkAAAr1SURBVNyX9582LNHhiIjERJCkcKK7XxTLk5rZDqAKCAFN7j7dzAYCDxJ5ItkBfNLdD8byvLFy+R9epaK2kc17q/jfj00hLU1PCSLSOwRpffQvM5vS+W7H7EJ3P8vdp0eXbwKec/dxwHPR5R5rT0Udg/OyuHRaQaJDERGJmY4m2VkPhKP7XG1m24gUHxng7n5GjGO5GLgg+v4e4B9Epv/scWoamqiobeQb7x9PTqamqhaR3qOj4qMCIF7tLB34u5k58Ad3nw8Mc/c9AO6+x8yGtnWgmV0DXAMwevToOIXXvsWritmwuxKA+14vZNTAPlwyVU8LItI7dJQUtrv7zjidd6a7747e+J8xs81BD4wmkPkA06dP9zjF16bFq4q5aeFawtGz7qmoY96idQBKDCLSK3SUFIaa2dfb2+juPz/ek7r77ujfEjN7hMh0n/vMbET0KWEEUHK8nx8vty/dQl1T+LB1tY0hbl+6RUlBRHqFjiqa04E8oF87r+NiZn3NrF/ze+ADwHrgUeCq6G5XEWkK26PsLq89pvUiIsmmoyeFPXHqoDYMeCTa2SsDuM/dnzazfwMLzOwLQCHwiTicu0tOyM2kvLbxqPUj83MTEI2ISOx1lBTi0vg+OmzGmW2sLwNmxeOcseDu5GSmYbXRkQGjcjPTuXH2hITFJSISSx0VH/XYG3QivLK1jL2V9cydMYqs9Mh/toL8XG69dIrqE0Sk12j3ScHdD3RnID3d3a9sZ3BeFt//6GS2lR4C4MFrz0twVCIisRVkmIuUt33/IZ7fXMINs8aRnZGuZCAivVaQYS5S3l9e3UlmuvHpc7u/s5yISHdSUuhETUMTD63YxUWnj2Bov5xEhyMiEldKCp1Ysno3VXVNfPa8kxIdiohI3CkpdMDduffVnUwc3o/pmn9ZRFKAkkIHVhYeZNOeSj573hjNrCYiKUFJoQP3vrqTfjkZXDJ1ZKJDERHpFkoK7SitqufJdXu47OwT6ZOllrsikhqUFNrx8IoiGkPOp89RBbOIpA4lhTaEw84D/y7knLEDOXVoXqLDERHpNkoKbXh1Wxk7y2r41DnqrCYiqUVJoQ33vVbIgD6ZzJ48PNGhiIh0KyWFI5RW1bN0w14+Pu1EcjLTEx2OiEi3UlI4wsKVRTSFnbkzVHQkIqlHSaGVcNh54PVCZqiCWURSlJJCK8u2l7GjrIZP6SlBRFKUkkIrDy0von9OBhedrgpmEUlNSgpRlXWNPLluDxefVaAKZhFJWUoKUY+t2U19U5hPTh+V6FBERBJGSSFqwfIiJg7vx+kF/RMdiohIwigpAFv2VrFmVzmfmD5KQ2SLSEpTUgAeWr6LzHTjkrM0RLaIpLaUTwoNTWEeWVXMrInDGJSXnehwREQSKuWTwgtbSig71MAnpp+Y6FBERBIu5ZPCIyuLGZyXxXvGD0l0KCIiCZfSSaG8poHnNu/jo2cWkJGe0v8pRESAFE8Kj6/dQ2PIuXRaQaJDERHpEVI6KSxaWcT4YXlMHqm+CSIikMJJYcf+Q6wsLOdjU09U3wQRkaiMRAdwJDO7CPgVkA78n7vfFsvPX7yqmNuXbqG4vBaA3KyUzYsiIkfpUXdEM0sHfgt8EJgEXGFmk2L1+YtXFTNv0bqWhADw46e2sHhVcaxOISKS1HpUUgBmAFvdfZu7NwAPABfH6sNvX7qF2sbQYetqG0PcvnRLrE4hIpLUelpSKAB2tVouiq5rYWbXmNlyM1teWlp6TB++u9UTQpD1IiKppqclhbZqfP2wBff57j7d3acPGXJsHc5G5uce03oRkVTT05JCEdB6QoMTgd2x+vAbZ08g94gJdHIz07lx9oRYnUJEJKn1tNZH/wbGmdlYoBiYC3wqVh9+ydRISdTtS7ewu7yWkfm53Dh7Qst6EZFU16OSgrs3mdl1wFIiTVLvdvcNsTzHJVMLlARERNrRo5ICgLs/CTyZ6DhERFJRT6tTEBGRBFJSEBGRFkoKIiLSQklBRERamLt3vlcPZWalwM7jPHwwsD+G4SSSrqVn6i3X0luuA3QtzU5y9zZ7/yZ1UugKM1vu7tMTHUcs6Fp6pt5yLb3lOkDXEoSKj0REpIWSgoiItEjlpDA/0QHEkK6lZ+ot19JbrgN0LZ1K2ToFERE5Wio/KYiIyBGUFEREpEWvTwpmdpGZbTGzrWZ2Uxvbs83swej218xsTPdHGUyAa/mcmZWa2ero6z8SEWdnzOxuMysxs/XtbDczuyN6nWvNbFp3xxhUgGu5wMwqWn0n/93dMQZhZqPM7AUz22RmG8zshjb2SYrvJeC1JMv3kmNmr5vZmui1/E8b+8T2HubuvfZFZPjtt4CTgSxgDTDpiH2+DPw++n4u8GCi4+7CtXwO+E2iYw1wLe8GpgHr29k+B3iKyEx85wKvJTrmLlzLBcDjiY4zwHWMAKZF3/cD3mjj31dSfC8BryVZvhcD8qLvM4HXgHOP2Cem97De/qQwA9jq7tvcvQF4ALj4iH0uBu6Jvn8YmGVmbU0LmmhBriUpuPuLwIEOdrkYuNcjlgH5Zjaie6I7NgGuJSm4+x53Xxl9XwVs4oj50UmS7yXgtSSF6H/r6uhiZvR1ZOugmN7DentSKAB2tVou4uh/HC37uHsTUAEM6pbojk2QawH4ePTR/mEzG9XG9mQQ9FqTxXnRx/+nzGxyooPpTLT4YSqRX6WtJd330sG1QJJ8L2aWbmargRLgGXdv93uJxT2styeFtrLlkVk2yD49QZA4HwPGuPsZwLO8/esh2STLdxLESiLjzJwJ/BpYnOB4OmRmecBC4KvuXnnk5jYO6bHfSyfXkjTfi7uH3P0sInPWzzCz04/YJabfS29PCkVA61/LJwK729vHzDKAE+iZxQGdXou7l7l7fXTxLuDsboot1oJ8b0nB3SubH/89MqtgppkNTnBYbTKzTCI30b+5+6I2dkma76Wza0mm76WZu5cD/wAuOmJTTO9hvT0p/BsYZ2ZjzSyLSCXMo0fs8yhwVfT9ZcDzHq2x6WE6vZYjync/SqQsNRk9Cnw22trlXKDC3fckOqjjYWbDm8t3zWwGkf/nyhIb1dGiMf4R2OTuP29nt6T4XoJcSxJ9L0PMLD/6Phd4H7D5iN1ieg/rcXM0x5K7N5nZdcBSIq137nb3DWb2A2C5uz9K5B/PX8xsK5HsOjdxEbcv4LVcb2YfBZqIXMvnEhZwB8zsfiKtPwabWRHwPSIVaLj774nM0T0H2ArUAFcnJtLOBbiWy4AvmVkTUAvM7aE/OmYCVwLrouXXAN8GRkPSfS9BriVZvpcRwD1mlk4kcS1w98fjeQ/TMBciItKitxcfiYjIMVBSEBGRFkoKIiLSQklBRERaKCmIiEiLXt0kVSRWzGwQ8Fx0cTgQAkqjyzXu/s6EBCYSY2qSKnKMzOz7QLW7/zTRsYjEmoqPRLrIzKqjfy8ws3+a2QIze8PMbjOzT0fHw19nZqdE9xtiZgvN7N/R18zEXoHI25QURGLrTOAGYAqRXrXj3X0G8H/Af0X3+RXwC3d/B/Dx6DaRHkF1CiKx9e/m8YDM7C3g79H164ALo+/fB0xqNeR9fzPrFx37XyShlBREYqu+1ftwq+Uwb///lgac5+613RmYSBAqPhLpfn8HrmteMLOzEhiLyGGUFES63/XA9OgMeRuB/0x0QCLN1CRVRERa6ElBRERaKCmIiEgLJQUREWmhpCAiIi2UFEREpIWSgoiItFBSEBGRFv8frJSVG5EIzj0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret2 = ensemble_simulations(duration, ndiv=20, y0=y0, volume=V, model=m, solver='gillespie', repeat=N)\n", "ret2.plot('o', ret1, '-')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Simulating with `meso`:" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEGCAYAAACKB4k+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXxU9b3/8dcnGwQCCbsQVi1gERUxApbW2rqgXUTrWltr1VZba1d/VOntvbW7vbT22s2qrRU3UCsutVZcW2tFkU1WUWQPO8iedebz+2MmIYRJcpLMZLb38/GYx8ycOWfm82XC+cz5rubuiIiIAOQkOwAREUkdSgoiIlJPSUFEROopKYiISD0lBRERqZeX7ADao3fv3j506NBkhyEiklbmz5+/w937xHotrZPC0KFDmTdvXrLDEBFJK2a2rqnXVH0kIiL1lBRERKSekoKIiNRTUhARkXpKCiIiUk9JQURE6iUsKZjZIDN72cxWmNkyM/tmdPstZlZuZouit080OGaqma0ys5VmNilRsYmISGyJHKdQC9zo7gvMrBsw38yej772a3f/ZcOdzWwUcBlwHDAAeMHMRrh7KIExiohIAwm7UnD3ze6+IPp4H7ACKG3mkMnATHevcvc1wCpgXKLiE5HMcOmdc7j0zjnJDqNDJbLMHdKmYGZDgZOAN6KbbjCzxWZ2j5n1iG4rBTY0OGwjMZKImV1rZvPMbN727dsTGLWIdJRsPLGnqoQnBTMrAh4DvuXue4E7gGOAMcBm4Fd1u8Y4/Ihl4dz9Lncvc/eyPn1iTt0hIhKIktGREpoUzCyfSEJ40N1nAbj7VncPuXsYuJtDVUQbgUENDh8IbEpkfCISH9l4cs3UMiey95EBfwZWuPttDbb3b7DbBcDS6OOngMvMrJOZDQOGA3MTFZ+ISDKlalJJZO+jicAVwBIzWxTd9j3gs2Y2hkjV0FrgOgB3X2ZmjwDLifRc+pp6Hol0nLoT1MPXnZrkSCSZEpYU3P1VYrcTPNPMMT8FfpqomEREpHka0SwiIvWUFEQySKrWU0v6UFIQkbT1xMJyFq7fzRtrdjHx1pd4YmF5skNKqNpQmANVtdSEwtSGj+ixHxdpvRyniGSvJxaWM3XWEqpDYQDKd1cwddYSAM4/qbnJEw4dv3D9bqpDYSbe+hJTJo0MdFws7k5FTYgDVSEOVtdyoCpERU0tB6tDHKwOUVEdoqIm8riyJvJ83c4DhB2+8/AiKmqi22tCVNaEqawJUVUbpqomRGWD+1CjRNDeuGNRUhCRtDRt9koqag7voFhRE2La7JUtniRjJZSbHlvMO1v3ccLAYvZW1rK/spb9VQ1ulbUcqKplX1XkftW2/YTdOe5/nuVgTQhv4w/3J9/aRK+uBfTr3pnO+Tl065xHn26d6JyfS+e8HDrn59IpL4eCvBzW7DjA88u31l8ltDYRBqGkICLt0t5f3G09ftPuipjby3dX8PjCjew+WMPugzXsqTjytnr7fhrXvlTVhvnDP9874v065+dQ1CmPok55dI3e3J2aUJiwQ0EenP3BvowZ3JOiTrl0KcijS0EuXTpF7gvzcyP30ccvLN/Kfz+5lIqaSEIKhZ19lbV87xPDWiz3xFtfOqLaKGgiDEpJQSTFpNN4gXhU4Rx5/GIOVtcyblhPtu+rZueBKnbur2bn/ip2HKhm1/5qdh2sJjfHmqxX//bDb9U/7t45j+Iu+RQXRm79uhexatv+mMcZ8PdvfIRunSNJoKhzHvm5hze91sVc99EHqkK88u5Ozj1+QKAy//qFd+sTQp2gJ/amEmFT29tCSUFE2vxrvTVVOKGws3N/FVv3VrF1byVb9lbyi2ffjnF8mO89vpTGzKBHlwJ6do3cRg3oxrJN+w6rZy/IzeH6jx3D5DGllBTm070wn9ycI4dLTbz1JcpjnEgHlBQyakD3uJU5lvac2AeUFDYZd7woKYgkQLb82m+uCucnTy9n855KNu+pYMueSrbtq2pVj5nbLxtD76JO9C7qRK+iAnp0KTjiBP/EwnK++9fFVIfClJYUBk5mUyaNZOqsJYed3Avzc5kyaWSLx7b313p7TuztiTsoJQWRLBf0l++egzWs23WA9bsOsmFXBRvfP0hBXg5VteHGbwnAA2+sY0BxIUcVd2bCMb3oX9yZo7p3pl/01rd7Jy684zU27a484tjSkkImj2n55H7+SaXMmLseaF0CritXWxJKe3+tt+fE3p64g1JSEMkQiWiw/fqMhazbeYC1Ow6wt7L2sNdLuuTTp6gTm/ZUHNZo2ykvh//+1Cg+N34wkXkxm/bdSccm/JdvU9qaUNr7a729J/a2xh2UkoJIDOlU/QOtqwIKhZ0Nuw6yatt+Vm3fT2FBLgerY889uXjjbob06sp5Y4oZ0rMrg3t1YVCPLgzqWUi3zvn1n92eExwk9pdvvMUj5kSf2NtDSUEkAzRVBfTzf6yguDCft7fs452tkduqbfsPq/Ip6pRLjnHEr/2fXTCaC08eREvae4JL5RNkU9Ix5qCUFERSSLyrgLbureKqe98EoH9xZ4b368apR/dieL8iPtC3Gx/oW0RxYX67fu1LZlFSEEkRQauAakJhVm7Zx5LyPZHbxj1HrlsbVVyYz91fKGNkv24Ud8lv8rMz+ZevtI6SgkiKaK4KKC/XWLBuN29t3M3S8j311T/dO+cxurSYj4/sw79X7aAmdCg9FObn8sPzjmPcsJ4dWg5Jb0oKInGWiCqgGx5aSKe8HEaXFvP5CUM4cVAJJ5QWM6RXl/oePqoCknhQUhCJo9YOBKsJhVm8cQ+vr97ZZJ//4sI87r9mPB/s3/2IKRcaytYqoPaUNZn/Tqn6HSkpiMRRSwPBwmFn+ea9/GfVDl5dtYN5a9+v3/+o7p3Yvr/6sGkbIlVAozlhYEmHlkOyl5KCSBy1NBDsP6t2sOtANQDD+xZxSdlAJhzdi/FH96Jn1wJVAUkgibzKUFIQiaOmpkAAeH31Tk4f2YePDO/Nh47pTb/unY/YJ12rgNIp1njJ1DIrKYjEwf6qWv61cjv9izsdkRTycoxvnzWC608/psVpH7JVpp5g05GSgkgb7TpQzfPLtzB72VZeXbWD6towPbrkUzakhEUb9lAbdlUBSdpRUhBppLkupbsOVPOPpZv5x5ItzFm9k1DYGdijkCsmDOHsUf0oG9qT3BxLu7mTIL1ilcRRUpCM1ZYTc6wupTfPWszcNTsp313Jq6t2EAo7w3p35SsfPZpzR/fnuAHdVS0kGUNJQaSBWF1KK2vCPDR3AwN7FHLtaUfzqRP6M6q/EoFkJiUFkQaa6lJqwL+/+7GUTwSqApL2UlKQrHewupan39rMjDfXNzmx3ICSwg5LCDqxSzIpKUjWemfrPh54fR2zFpSzv6qWY/p0ZfKJA5i9fAuVNYemm+iolcBEUoGSgqS0ePfiqQ2FeW75Vqa/tpY31uyiIC+HTx3fn8+OH0zZkB6YmUYVS1ZLWFIws0HAfcBRQBi4y91vN7OewMPAUGAtcIm7v2+Ra/PbgU8AB4EvuvuCRMUn2WX3wWpmzN3A/XPWsmlPJYN6FjL13GO5uGwQPbsWHLZvPEYVqwpI0lWLScHMLgaedfd9ZvZ9YCzwkwAn7FrgRndfYGbdgPlm9jzwReBFd7/VzG4GbgZuAs4Fhkdv44E7ovcirdZwrMFxP3iW6powNWHnQ8f04oeTR/PxY/uSm5PajcYiydD0PLyH/Hc0IXwYmARMJ3LCbpa7b65LHO6+D1gBlAKTo+9B9P786OPJwH0e8TpQYmb9W1UaESIJ4abHFtePNThQFSLsMGXSSB768gTOGtVPCUGkCUGSQl2n7U8Cd7j7k0BBM/sfwcyGAicBbwD93H0zRBIH0De6WymwocFhG6PbGr/XtWY2z8zmbd++vTVhSBaYt3YXNz22+Ih1CULuPPTG+iRFJZI+grQplJvZncCZwC/MrBPBkgkAZlYEPAZ8y933NtOtL9YLR/QQdPe7gLsAysrKmupBKFnmjdU7uf3Fd3ntvZ1N7tPUGAQROSTIyf0SYDZwjrvvBnoCU4K8uZnlE0kID7r7rOjmrXXVQtH7bdHtG4FBDQ4fCGwK8jmSveave5/L736dS+96nXe37ef7n/wgA4qPnJIaImMNRKR5LV4puPtBM9sGfBh4l0gD8rstHRftTfRnYIW739bgpaeAK4Fbo/dPNth+g5nNJNLAvKeumkmksWWb9vDL2St5eeV2ehcV8N+fGsXnxg+mc34uvYs6MXXWksOmq9BYA5FggvQ++gFQBowE/gLkAw8AE1s4dCJwBbDEzBZFt32PSDJ4xMyuAdYDF0dfe4ZId9RVRLqkXtWqkkjGiTVb6clDevCr51byxKJNFBfm891zRvLFDw2lS8GhP+W6MQUaayDSekHaFC4g0khc15NoU7SLabPc/VVitxMAnBFjfwe+FiAeyQKxZiu98dG3wJ283ByuP/0YrvvoMRQX5sc8Pl1XMBNJtiBJodrd3cwcwMy6JjgmyTBtGZUca7bSUNjpUpDLSzeezlFNtBuISPsEaWh+JNr7qMTMvgy8ANyd2LAk2zXVU6iiOqSEIJJAQRqaf2lmZwF7ibQr/I+7P5/wyCRrle+uoFN+zmGT0tXpqB5EqnKSbBVo7qNoElAikISqCYW559U1/N8L71IbCpOXY9SGDw1FUQ8ikcRrMimY2T4ig8eMwweRGZF24e4Jjk2yyIL17/O9WUt4e8s+zvxgP245bxTz1r6vHkQiHazJpODuLfYwEmmvA1W1TJu9kulz1nJU987cecXJTDruKAAG9ujSrh5EqgISab0g4xQGx9ru7ppIRtrl3+9u5+bHllC+u4IrTx3ClHOOpaiTlvgQSaYg/wP/3uBxZ2AYsBI4LiERScbbX1XLz55ZwUNvrOfoPl159CuncsrQnskOS0QI1vvo+IbPzWwscF3CIpKM0nhU8gUnDeDJtzax8f0Krj3taL5z1gg65+cmO0wRiWr1tXp00ZxTEhGMZJZYo5J/9/J79OqazyPX6epAJBUFaVP4ToOnOURWXtNCBtKiWKOSAQrycpUQRFJUkCuFhr2Qaom0MTyWmHAkkzQ1KnnLnsoOjkREggrSpvDDjghEMsuBqlo65+fGvFJozahkdSsV6Vgtzn1kZs+bWUmD5z3MbHZiw5JUc+mdc+ontmvJyi37+PTvXqWiJkReo7WQNSpZJLUFmRCvT3TFNQDc/X0OrasscpgnF5Vz/u//w77KWmZ8eQK/vPhECnIjf2alJYX8/DPHa1SySAoL0qYQMrPBdYPVzGwIMdZOluxWXRvmp39fzvQ56xg3tCe/u/wk+naPzGaqdQ1E0keQpPBfwKtm9q/o89OAaxMXkqSbHfuruP7BBcxds4trPjyMm889lvzcIBehIpJqgjQ0PxsdsDaByGR433b3HQmPTNLC0vI9XHf/fHbsr+L2y8YweYyqhkTSWZCGZgPOAca6+9+ALmY2LuGRScp7dulmLvrja4Td+etXPqSEIJIBglzj/wE4Ffhs9Pk+4PcJi0hSnrvzx3+9x1ceWMAH+3fnqRs+zPEDi5MdlojEQZA2hfHuPtbMFkKk95GZFSQ4LklRNaEw3398KQ/P28CnTxzAtItO0NxFIhkkSFKoMbNcoj2OzKwPcOQ6iZKxGk5qN/oHs6mqDfONj3+Ab581gkjtoohkiiBJ4TfA40BfM/spcBHw/YRGJSmj8aR2VbVh8nONo/sUBU4I6ooqkj6C9D560MzmA2cQ6X10vruvSHhkkhJiTWpXE3KmzV6pQWgiGai5NZobTmO5DZjR8DV335XIwCQ1NDWpXVPbRSS9NXelMJ9IO0KsOgIHjk5IRJIy3lzbdN5vzaR2IpI+mkwK7j6sIwOR1PLKO9u59v559OlWwJ6KWqpqD/Ut0KR2Ipkr0MprZnYekektAP7p7k8nLiRJthdXbOWrDyzgA32LuP+acfz73R1896+LqQ6FKS0pZMqkkWpPEMlQQVZeuxU4BXgwuumbZjbR3acmNDKJq7ppr1vqCfTs0i18fcYCRvXvzn1Xj6e4Sz7nn1SqSe1EskSQK4VPAGPcPQxgZtOBhYCSQoZ5Zslmvj5jIWMGlfCXq06he+f8ZIckIh0s6FSWJQ0eB5rPwMzuMbNtZra0wbZbzKzczBZFb59o8NpUM1tlZivNbFLAuCROIlcICxk7uITpV49TQhDJUkGuFH4OLDSzl4n0RDqNYFcJ9wK/A+5rtP3X7v7LhhvMbBRwGXAcMAB4wcxGuPuRazlK3D2/fCs3PLSAEwcW85erxlHUKVBTk4hkoCCD12aY2T+JtCsYcJO7bwlw3CtmNjRgHJOBme5eBawxs1XAOCDY+o/SZv9cuY3rH5zPcaXF3Hu1EoJItmuy+sjMxtbdgP7ARmADMCC6ra1uMLPF0eqlHtFtpdH3rrMxuk0SaO6aXXzlgfmM6NeN+1RlJCI0f6UwD1gGbI8+bziIzYGPt+Hz7gB+HD3+x8CvgKtpeoDcEczsWqIrvw0ePLgNIQjAko17uPreNyktKeS+q8dRXKiEICLNJ4UbgQuBCmAm8Li772/Ph7n71rrHZnY3UDfeYSMwqMGuA4FNTbzHXcBdAGVlZVorug1WbdvPF+55g+LCfB740nh6FXVq8Rh1RRXJDs2NaP418GszG0ZkgZ0XzWwd8DN3X9SWDzOz/u6+Ofr0AqCuZ9JTwENmdhuRhubhwNy2fIYcqeHU1xN+9iKVNbXk5eby4JfG079Y01WIyCFBGprXmNmTQCFwBTACaDEpmNkM4HSgt5ltBH4AnG5mY4hUDa0Frot+xjIzewRYDtQCX1PPo/hoPPX1lr2VANx49tEM7d01maGJSAoy99g1MGZ2NJFuopOJNALPBJ5298qOC695ZWVlPm/evGSHkdIm3voS5TFmNC0tKeQ/N7elWUhE0p2ZzXf3slivNXelsApYDDwJ7AUGA9fXLazi7rfFOU5JAE19LSKt0VxS+BGHegAVdUAskgADSgpjXilo6msRiaW5huZbOjAOSZCPDO/NzDc3HLZNU1+LSFOCzn0kaei1VTv46/yNjOxXREFu5KsuLSnk5585XlNfi0hMmtMgQ63beYCvPriAYb278uhXP8SXp0ca5DXeQESa09w0F9+M3k/suHAkHvZV1vCl6fMwgz9dWabpK0QksOaqj66K3v+2IwKR+AiFnW/NXMTqHQf4w+VjGdJLYxFEJLjmqo9WmNlaoI+ZLW6w3QB39xMSGpm0yW3Pr+TFt7fxo8nH8aEP9E52OCKSZprrffRZMzsKmA2c13EhSVs9t2wLv3/5PS47ZRBXTBiS7HBEJA0129AcXTfhRDMrIDK9BcBKd69JeGTSKqu37+fGR97ihIHF3HLecdQNMhQRaY0Wex+Z2UeJrJ62lkjV0SAzu9LdX0lwbNLIpXdG1hxq3IPoQFUt190/n/y8HO74/Ml0zs9NRngikgGCdEm9DTjb3VcCmNkIYAZwciIDk2DcnamzlvDe9v3cf814SpsYqayuqCISRJDBa/l1CQHA3d8B1McxRcyYu4Gn3trEjWePZKIalkWknYJcKcwzsz8D90effw6Yn7iQJKjlm/Zyy9+WcdqIPnz1o8ckOxwRyQBBksJXga8B3yDSpvAK8IdEBiUt21dZw9ceWkCPLvn8+pITyclRw7KItF+QRXaqiLQraKrsFOHufP+JpazbeYCZ154aaDlNEZEgNPdRmmi4pObYHz/P+wdr+M5ZIxg3rGeyQxORDKJZUtNA4yU13z9YQ47BQK2JICJx1mxSMLNcM5vWUcFIbNNmr6Si5vAlq8MOv3r+nSRFJCKZqtmk4O4h4GTT8Nik0pKaItJRgrQpLASeNLNHgQN1G919VsKiksNoSU0R6ShB2hR6AjuBjwOfjt4+lcig5HDfPGM4jS/VtKSmiCRCkC6pV7W0jyTWsk17cCAvx6gNO6UlhUyZNFJLaopI3AWZEG8EcAfQz91Hm9kJwHnu/pOERye8+u4Ops9Zx1UTh7J8015A8xiJSOIEqT66G5gK1AC4+2LgskQGJRH7Kmv47l/f4ug+XbnpnGOTHY6IZIEgSaGLu89ttK02EcHI4X7x7Nts3lvJtItO1HTYItIhgiSFHWZ2DOAAZnYRsDmhUQmvr97JA6+v5+qJwzh5SI9khyMiWSJIl9SvAXcBx5pZObCGyEypkiAV1SFuemwxQ3p14f+drR5GItJxgvQ+Wg2caWZdgRx335f4sLLbr55bybqdB5nx5QkUFhyqNlIDs4gkWovVR2bWy8x+A/wb+KeZ3W5mvRIfWnZasnEP9/xnDZePH8ypx+ifWUQ6VpA2hZnAduBC4KLo44cTGVS2qg2Fmfr4YnoVdVJvIxFJikAjmt39x+6+Jnr7CVDS0kFmdo+ZbTOzpQ229TSz583s3eh9j+h2M7PfmNkqM1tsZmPbXqT0NX3OOpaW7+UHnx5FcaFWPBWRjhckKbxsZpeZWU70dgnw9wDH3Quc02jbzcCL7j4ceDH6HOBcYHj0di2RwXJZpXx3Bb96biWnj+zDJ4/vn+xwRCRLNZkUzGyfme0FrgMeAqqjt5nAt1t6Y3d/BdjVaPNkYHr08XTg/Abb7/OI14ESM8uqM+MPnlxG2J0fTx6NJqUVkWRpsveRu3dLwOf1c/fN0fffbGZ9o9tLgQ0N9tsY3XbEeAgzu5bI1QSDBw9OQIgd76W3t/LCiq3cdM6xDOrZJdnhiEgWC7QcZ3S+o6EN94/z1Nmxfhp7rB3d/S4i4yYoKyuLuU8qu/TOOcCh7qWVNSFueWo5x/TpyjUfHpbM0EREAk2Idw9wArAMCEc3O9CWpLDVzPpHrxL6A9ui2zcCgxrsNxDY1Ib3Tzt3/ms163cd5IFrxlOQp9VRRSS5glwpTHD3UXH6vKeAK4Fbo/dPNth+g5nNBMYDe+qqmTLZhl0H+cM/V/HJ4/vz4eG9kx2OiEig3kdzzKzVScHMZgBzgJFmttHMriGSDM4ys3eBs6LPAZ4BVgOriMzKen1rPy8d/ejp5eSY8V+f/GCyQxERAYJdKUwnkhi2AFVE6v/d3U9o7iB3/2wTL50RY18nMsdSRntiYTkL1++mOhTm5B8/z84D1UyZNFLLaopIygiSFO4BrgCWcKhNQVrpiYXlTJ21hOpQ5J9w54FqDOjbrVNyAxMRaSBI9dF6d38qOpp5Xd0t4ZFlmGmzV1JREzpsmwP/98K7yQlIRCSGIFcKb5vZQ8DfiFQfAXHvkprxNu2uaNV2EZFkCJIUCokkg7MbbGtrl9SsNaCkkPIYCUDtCSKSSoKsp3BVRwSS6aZMGslNjy2mqvZQs0xhfi5TJmkRHRFJHUEGr/2FGKOL3f3qhESUoc4/qZR7X1vDog17ACgtKWTKpJGcf1JpkiMTETkkSPXR0w0edwYuIEtGG8fTwvXvs2jDHkpLOjOwRxetoiYiKSlI9dFjDZ9HB6W9kLCIMpC78/Nn3qZ3USf6F6sNQURSV1sm2xkOZMb0pB3khRXbmLt2F986czi5OZoWW0RSV5A2hX0c3qawBbgpYRFlmNpQmF88+zZH9+nKpacM4vMThiQ7JBGRJgWpPkrEugpZ4/GF5azatp87PjeW/FzNgioiqa3JpGBmzVYRufv6+IeTWapqQ/zfC+9yfGkx54w+KtnhiIi0qLkrhb8TqTZqWAnuQB+gL5CbwLgywsy5GyjfXcHPP3O8ltgUkbTQ3HKcxzd8bmZDibQlnAn8LKFRZYCD1bX89qVVTDi6Jx/RWgkikiZarOQ2s+Fmdi/wD2A+MMrdf5vowNLdva+tZcf+KqZMGqmrBBFJG821KYwG/gs4Dvhf4Bp3DzW1vxyyt7KGO/+1mjOO7cvJQ3omOxwRkcCaa1N4C9hApG1hHDCu4S9ed/9GYkNLX/f+Zy17Kmr49lkjkh2KiEirNJcUNLdRG+ytrOFP/17NWaP6Mbq0ONnhiIi0SnMNzdM7MpBMce9/1rK3spZvnjE82aGIiLSaRlPFka4SRCTdKSnEka4SRCTdKSnEyT5dJYhIBggyTmGEmb1oZkujz08ws+8nPrT0cv/r69hbWcs3Pq6rBBFJX0GuFO4GpgI1AO6+GLgskUGlm4rqEH/+9xpOG9GH4wfqKkFE0leQpNDF3ec22labiGDS1cNvrmfngWpu+NgHkh2KiEi7BEkKO8zsGKJrKpjZRcDmhEaVRqprw9z5ympOGdqDccM0ellE0luQNZq/BtwFHGtm5cAa4PMJjSqNPLGwnM17Kvn5Z45veWcRkRQXZJGd1cCZZtYVyHH3fYkPKz2Ews4d/3qP0aXd+eiIPskOR0Sk3YIsx9kJuBAYCuTVzX/k7j9KaGRpYPayLazZcYDfXz5WM6GKSEYIUn30JLCHyLTZVYkNJz1ceucc3J2qkDOkVxetqiYiGSNIUhjo7ufE80PNbC2wDwgBte5eZmY9gYeJXJGsBS5x9/fj+bnxtK+ylhVb9vHj80eTm6OrBBHJDEF6H71mZoloRf2Yu49x97Lo85uBF919OPBi9HnK2rynkp5dC7j45IHJDkVEJG6aTApmttTMFgMfBhaY2UozW2xmS6Lb420yUDcz63Tg/AR8Rrs9sbCcBevfZ3dFDTWhMM8u3ZLskERE4qa56qNSYEyCPteB58zMgTvd/S6gn7tvBnD3zWbWN9aBZnYtcC3A4MGDExRebE8sLGfqrCXUhByIVCFNnbUEgPNPKu3QWEREEqG5pLDG3dcl6HMnuvum6In/eTN7O+iB0QRyF0BZWZknKL6Yps1eSUXN4SuSVtSEmDZ7pZKCiGSE5pJCXzP7TlMvuvttbf1Qd98Uvd9mZo8TWe5zq5n1j14l9Ae2tfX9E2XT7opWbRcRSTfNNTTnAkVAtyZubWJmXc2sW91j4GxgKfAUcGV0tyuJdIVNKf2LO8fcPqCksIMjERFJjOauFDYnaIBaP+Dx6GCvPOAhd3/WzN4EHjGza4D1wMUJ+Ox2OW1kH2bO3XDYtsL8XKZMGpmkiERE4qu5pJCQzvfRaTNOjLF9J3BGIj4zHtydRet30797J3YeqKE6FKa0pJApk0aqPUFEMkZzSSFlT9DJMPLzr+cAAAjQSURBVGf1Tt7eso//vfAEHluwEYCHrzs1yVGJiMRXk0nB3Xd1ZCCp7p5X19KzawHnjRnAJacMSnY4IiIJoTWaA1i38wAvvr2Vz48fTOf83GSHIyKSMEoKAdw/Zx25ZnxuwpBkhyIiklBKCi2oqA7xyLwNTDruKPp1j90lVUQkUygptOCpt8rZW1nLFafqKkFEMp+SQjPcnfvmrGNEvyLGa/1lEckCSgrNWLB+N8s27eULpw7VymoikhWUFJpx/5y1dOuUxwUanCYiWUJJoQk791fxzJItXHjyQLp2CrJAnYhI+lNSaMKj8zdSHQrz+Qkdu2aDiEgyKSnEEA47M+euZ9zQnnygb5snhBURSTtKCjG8vnona3ce5PLxukoQkeyipBDDQ3PXU1yYzzmjj0p2KCIiHUpJoZEd+6uYvWwLF44dqHmORCTrKCk08tj8jdSEnMvHayZUEck+SgoNuDsz1MAsIllMSaGBOdEG5s/qKkFEspSSQgOPzttIt855nDu6f7JDERFJCiWFqL2VNTyzZDOTxwxQA7OIZC0lhai/vbWJqtowl5Sp6khEspeSQtSj8zYysl83ji8tTnYoIiJJo6QAvLN1H4s27ObisoGaIltEspqSAvDovA3k5ZimyBaRrJf1SaEmFGbWgnLO/GA/ehV1SnY4IiJJlfVJ4eW3t7HzQDUXlw1MdigiIkmX9Ulh1oJyehcV8NERfZIdiohI0mV1Uth9sJqX3t7GeSeWkpeb1f8UIiJAlieFpxdvpjoU5jNj1cAsIgJZnhQeX1jOiH5FHDege7JDERFJCVmbFNbuOMD8de9zwUkamyAiUicv2QE0ZmbnALcDucCf3P3WeL7/EwvLmTZ7JeW7KwAoLMjavCgicoSUOiOaWS7we+BcYBTwWTMbFa/3f2JhOVNnLalPCAC/+MdKnlhYHq+PEBFJaymVFIBxwCp3X+3u1cBMYHK83nza7JVU1IQO21ZRE2La7JXx+ggRkbSWakmhFNjQ4PnG6LZ6Znatmc0zs3nbt29v1ZtvanCFEGS7iEi2SbWkEKvF1w974n6Xu5e5e1mfPq0bcDagpLBV20VEsk2qJYWNQMMFDQYCm+L15lMmjaSw0QI6hfm5TJk0Ml4fISKS1lKt99GbwHAzGwaUA5cBl8frzc+PzoI6bfZKNu2uYEBJIVMmjazfLiKS7VIqKbh7rZndAMwm0iX1HndfFs/POP+kUiUBEZEmpFRSAHD3Z4Bnkh2HiEg2SrU2BRERSSIlBRERqaekICIi9ZQURESknrl7y3ulKDPbDqxr4+G9gR1xDCeZVJbUlCllyZRygMpSZ4i7xxz9m9ZJoT3MbJ67lyU7jnhQWVJTppQlU8oBKksQqj4SEZF6SgoiIlIvm5PCXckOII5UltSUKWXJlHKAytKirG1TEBGRI2XzlYKIiDSipCAiIvUyPimY2TlmttLMVpnZzTFe72RmD0dff8PMhnZ8lMEEKMsXzWy7mS2K3r6UjDhbYmb3mNk2M1vaxOtmZr+JlnOxmY3t6BiDClCW081sT4Pv5H86OsYgzGyQmb1sZivMbJmZfTPGPmnxvQQsS7p8L53NbK6ZvRUtyw9j7BPfc5i7Z+yNyPTb7wFHAwXAW8CoRvtcD/wx+vgy4OFkx92OsnwR+F2yYw1QltOAscDSJl7/BPAPIivxTQDeSHbM7SjL6cDTyY4zQDn6A2Ojj7sB78T4+0qL7yVgWdLlezGgKPo4H3gDmNBon7iewzL9SmEcsMrdV7t7NTATmNxon8nA9OjjvwJnmFmsZUGTLUhZ0oK7vwLsamaXycB9HvE6UGJm/TsmutYJUJa04O6b3X1B9PE+YAWN1kcnTb6XgGVJC9F/6/3Rp/nRW+PeQXE9h2V6UigFNjR4vpEj/zjq93H3WmAP0KtDomudIGUBuDB6af9XMxsU4/V0ELSs6eLU6OX/P8zsuGQH05Jo9cNJRH6VNpR230szZYE0+V7MLNfMFgHbgOfdvcnvJR7nsExPCrGyZeMsG2SfVBAkzr8BQ939BOAFDv16SDfp8p0EsYDIPDMnAr8FnkhyPM0ysyLgMeBb7r638csxDknZ76WFsqTN9+LuIXcfQ2TN+nFmNrrRLnH9XjI9KWwEGv5aHghsamofM8sDiknN6oAWy+LuO929Kvr0buDkDoot3oJ8b2nB3ffWXf57ZFXBfDPrneSwYjKzfCIn0QfdfVaMXdLme2mpLOn0vdRx993AP4FzGr0U13NYpieFN4HhZjbMzAqINMI81Wifp4Aro48vAl7yaItNimmxLI3qd88jUpeajp4CvhDt7TIB2OPum5MdVFuY2VF19btmNo7I/7mdyY3qSNEY/wyscPfbmtgtLb6XIGVJo++lj5mVRB8XAmcCbzfaLa7nsJRbozme3L3WzG4AZhPpvXOPuy8zsx8B89z9KSJ/PPeb2Soi2fWy5EXctIBl+YaZnQfUEinLF5MWcDPMbAaR3h+9zWwj8AMiDWi4+x+JrNH9CWAVcBC4KjmRtixAWS4CvmpmtUAFcFmK/uiYCFwBLInWXwN8DxgMafe9BClLunwv/YHpZpZLJHE94u5PJ/IcpmkuRESkXqZXH4mISCsoKYiISD0lBRERqaekICIi9ZQURESkXkZ3SRWJFzPrBbwYfXoUEAK2R58fdPcPJSUwkThTl1SRVjKzW4D97v7LZMciEm+qPhJpJzPbH70/3cz+ZWaPmNk7ZnarmX0uOh/+EjM7JrpfHzN7zMzejN4mJrcEIocoKYjE14nAN4HjiYyqHeHu44A/AV+P7nM78Gt3PwW4MPqaSEpQm4JIfL1ZNx+Qmb0HPBfdvgT4WPTxmcCoBlPedzezbtG5/0WSSklBJL6qGjwON3ge5tD/txzgVHev6MjARIJQ9ZFIx3sOuKHuiZmNSWIsIodRUhDpeN8AyqIr5C0HvpLsgETqqEuqiIjU05WCiIjUU1IQEZF6SgoiIlJPSUFEROopKYiISD0lBRERqaekICIi9f4/ddEHcDFUMP0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret2 = ensemble_simulations(\n", " duration, ndiv=20, y0=y0, volume=V, model=m, solver=('meso', Integer3(3, 3, 3), 0.25), repeat=N)\n", "ret2.plot('o', ret1, '-')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Simulating with `spatiocyte`:" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEGCAYAAACKB4k+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXxcdb3/8dcne9K0Tfd9o5uUUqCUlgIqSqGIyCKgiAsg94IiioooeLkuXFF+1ouiXhVQvIWLQIFSKiAVioIIlO6bpbS0JU26b0naZs/n98echDTNcprMZGYy7+fjMY+ZOcucz+k05zPnu5q7IyIiApAW7wBERCRxKCmIiEgDJQUREWmgpCAiIg2UFEREpEFGvAPoiL59+/rIkSPjHYaISFJZunTpHnfv19y6pE4KI0eOZMmSJfEOQ0QkqZjZey2tU/GRiIg0UFIQEZEGSgoiItJASUFERBooKYiISAMlBRERaaCkICIiDZQURESkgZKCiCS1T9/3Bp++7414h9FlxCwpmNkwM/ubma0zs7VmdnOw/AdmVmxmK4LHBY32ud3MNprZejObGavYRESkebEc5qIGuMXdl5lZd2Cpmb0YrPu5u/+s8cZmNgG4EjgBGAy8ZGbj3L02hjGKiEgjMbtTcPft7r4seF0GrAOGtLLLxcBj7l7p7puBjcDUWMUnIiJH65Q6BTMbCZwCLAoW3WRmq8zsQTPrFSwbAmxttFsRzSQRM7vezJaY2ZLdu3fHMGoR6SypWC+QqOcc86RgZvnAU8DX3b0U+C0wGjgZ2A78d/2mzezuRy1wv9/dp7j7lH79mh35VUTaoSMXqUS9wLUlWc85lseOaVIws0wiCeERd58L4O473b3W3euAB3i/iKgIGNZo96HAtljGJyIiR4pl6yMD/gCsc/d7Gi0f1GizS4E1wev5wJVmlm1mo4CxwFuxik9ERI4Wy9ZHZwKfB1ab2Ypg2XeBz5jZyUSKhrYANwC4+1ozmwP8i0jLpa+o5ZGISOeKWVJw99dovp7g+Vb2uQu4K1YxiUjXMm95McsLD1BVW8eZd7/MrTPHc8kprTVylLaoR7NIF5KMFb71F/ZFm/dx5t0vM295cej9bp+7mqraOgCKD5Rz+9zVofeX5ikpiEjcdOTCPmvBesqrjyxhLq+uZdaC9TGJNVG0N4mGpaQgInHTkQv7tgPlx7Q8kbT3wv70siJum7sqpndHsaxoFpEU0JFy/Y5c2AcX5FLczHaDC3JDHbsjcTfd91vnjeNjJw7iUGUNh6tqg0cN5fWvq2upCJYteW8/L6zZQU1dpBtW8YFybpmzkt//YxO9umVRUV1LRXUdFdW1VNYc/dxUfRKNVl2KkoJIgqmvE3j8humddsz2XiBbKv4BQu3fkQv7rTPHc/vc1UfcaeRmpnPrzPHtivs7T63i7R2lTBpawMGKGsoqazhYUcOhqhrKKmo4VBl5bNlziE17D+FB19riA+V8Y85KvjFnZZvHbUmtO2/vKOPEoT3JyUinb34GOZnpZGekNTxnZ6Zz/6ubmt0/mndHSgoiKa4jF/bWin/CJIWOXNjrP//bT0aKU/p3z+bTpw0jLyudOUu2UnK4mpLyag6UV1FSXkNJeTWlwWPL3kPUNRkvobKmjt+9cvRFNyczjfzsDPKzM+iWncHW/eUNCaGxHjkZfGvmeHIz08nLyiAvK53crHTygkdOZjq5melM+dFLRw/VANTWOU/feGar5/zcqu0dujsKQ0lBJMV15MLe0XL9phf2IQW5fGPGWKaO6s3qohL2HKxkz8FK9h2qYt+hKvYGz/sOVbH/cOS5PpntKqvkVy9vPOLz09OMnrmZ9MzNpEdOBj1yMxnaK5dNew41G48Bz9/8QfKzM+ieE0kEGelHVr2Ouu25Zvctq6jhC9NHtnnO8bo7CktJQSQG4lEE1F6dVa5/qLKGHaUV7CytYFdpZeS5rJLdZZVkZ6aRlgYHK2v41pOrmj1WVkYafbtl0atbFr27ZTGiTx698rJ48V87yUg3bjlvPL3yMinIzaIgL5OCvEzyszOIDK5wpOWFL7cY9/GDekTtnJsTzbujIQW5Ue+boaQgkuKi/cs1M92YOLgHtz6xkh2lFWwvqWBnSQVllTVH7Z+bmU7/Htm4R15/4qTB9M3PDh5Z9MnPpl9+Nr3zs+iWld7sBX7d9lIALjppcOhz7siFuaO/1jt6Yb/klCE8+lYhEJsfHUoKIs1Ipl/69dpbWRzmIldRXUvR/nK27jtM4b7DFO0/zNZ95RQdOIw1KSGvrnX+um4n/btnM7BnLmP65XPWmL4M6JHDwJ7ZDOiRw4AeOfTvnt3wS77+3/vOiydG6V+jdR25MEfj13qsL+wdoaQg0gV0pLK46UWud14mZ4zpw+vv7uHRtwp5b+9hdpRWHLFPdkYaQ3vlMrRXHpMmF/DqO7vJzkjjJ5+cxOCCyEU/Mz2xu0F15MKcyBf1jlJSEEkg7f21H7ay2N3ZXVbJxl0HeXf3QTbuOsimPYfYtPtQQ0LZd7iaZ1ftoG9+NiP75HHGmD4M753HiD55DO+dx7BeefTrnn1EUU79L/2po3pH459B4khJQSRBdOTXfmuVxbNf38L6nWW8s6OMDbsOUlJe3bA+PzuD0f26MXVUbxZv3kdOZhq/uPIURvbtRn62Lg+pSN+6SIJob9PQujqnX/dsdpVVHrXOge/PX0vP3EzGDcjn45MGMa5/PmP6d2dM/3wG9Hj/F3/9r/2JQ3pG76RC6mpFMMlMSUEkytpbBBSmaWhdnbN57yFWFR1gdVEpa4pLWLuthENVR089kpFmXPfBUXzxzFH0b1Lc05WkYkKJ5TkrKYhEUUeKgFpqGtqrWxazFrzNiq0HWFVUQllFpGlnTmYaEwb14PJTh3LC4J7sLKvgVws3xqz9uqQGJQWRKOpI7+BvnTuO255efdSgZ/sOVfG7VzZx/KDuXHTSYE4aWsCkYT0Z0y//qN62r23YA3T+r+dU/LXeUYn6b6akIBJFx9I7uLbOWbe9lDc37WXR5n0s2bLviISQk5HGeRMG8LnpIzlxSE9ys9JjFneq6siFOVEv6h2lpCDSREeGVG6td7C7s3HXQf65cQ+vv7uXNzftpTQoChrRJ48Zxw/gtFG9eeTN98jOSGPOl86I6nmJhKGkINJIR4eCbmnYhyEFOZz+k4XsLI20EBrWO5ePTRzEGWP6MG1UHwb2zGnY/qmlRdE8pdC66i9fOTZKCtJltWeoio4OBX3RSYMp3HeIe1/aQG0w+kN1rfPOroOcObovZ43ty1lj+jKsd174ExHpREoKIo20Z8TQiupa/rFhDy/9aycvr9/F7qC/QH52Ov/+wdF8eHw/ThzSk/S0rtkkVLoWJQWRRsKOGFpSXs3CdTt5Yc0O/rFhD+XVtXTPzuDD4/txzvH9efiN98hMT+PmGWM7K3SRqFBSEGmktRFDS8qr+evaHTy3ejv/3LiH6lpnYI8crpgylHMnDGDaqD5kZUSaiD721tZ2x6CyfYknJQWRRpqOGDq4Zw7nTBjAs6u2HTFU8rVnjuL8iQM5eWgBaQlULKSEIh2lpCDSxCdOGszvXnmXPQcrKSmv5uE33mNgjxw+P30EF04axMnDCrrskBEiSgoigc17DvHEkq08vbyY7SUVpKcZV5w6lItPHsK0Ub0T6o5AJFaUFKRLCtsBraK6ludXb+exxVt5a/M+0tOMD43tS7esdHrlZXH3ZZPiEL1I/LSZFMzsCuAFdy8zszuAycCP3H1ZzKMTaYcwHdA27jrInxYV8uTSrZRW1DCiTx7fPn88l08eSv8eOQ19HNpLZfuSrMLcKfynuz9hZmcBM4GfAb8FpsU0MpF2aqkD2k8XvE1uVjqzX9/C6+/uJTPdmHnCQK6aNpzpx/VRPYEI4ZJC/V/Xx4HfuvszZvaDtnYys2HAQ8BAoA64393vNbPewOPASGAL8Cl332+Rv8h7gQuAw8A1uhuR9vRKbrkDWgU3PLyUwT1zuHXmeD592jD65mdHJU6RriJMUig2s/uAGcD/M7NsIMyM3DXALe6+zMy6A0vN7EXgGmChu99tZrcBtwHfAT4GjA0e09DdiLRTSx3QstLT+OVnTmbG8QOOGnJaRCLC/GV8ClgAnO/uB4DewK1t7eTu2+t/6bt7GbAOGAJcDMwONpsNXBK8vhh4yCPeBArMbNCxnIwIwFXThpPepCgoOyONn14+ifMnDlJCEGlFm3cK7n7YzHYBZwEbiNwBbDiWg5jZSOAUYBEwwN23B5+93cz6B5sNARp3Ay0Klm1v8lnXA9cDDB8+/FjCkC5uWeF+frVwA39bv5usdMProM5p1yxkqiiWVBWm9dH3gSnAeOCPQCbwf8CZYQ5gZvnAU8DX3b20lcq85lb4UQvc7wfuB5gyZcpR6yX1vLV5H/cufId/btxLr7xMbjl3HF+YPpLrH14C6AIvcizC1ClcSuRXfn1R0LagjqBNZpZJJCE84u5zg8U7zWxQcJcwCNgVLC8ChjXafSiwLcxxJDUtfW8/P3/xHV7buIe++dl894IP8NlpI+iWre43Iu0V5q+nyt3dzBzAzLqF+eCgNdEfgHXufk+jVfOBq4G7g+dnGi2/ycweI1LBXFJfzCTS2LrtpcxasJ6X395Fn25Z3PHx4/nstBGarlIkCsIkhTlB66MCM/t34IvAAyH2OxP4PLDazFYEy75LJBnMMbPrgELgimDd80Sao24k0iT12tBnIQmtPc1K4eheyV88ayRrikuZt6KY7tkZfPv88Vw9faTuDESiKExF88/M7FyglEi9wvfc/cUQ+71G8/UEAOc0s70DX2nrcyU1NNcr+b+eXUdGGnzpw6P50odG0zMvM85RinQ9oX5iBUmgzUQgEi3N9UoG6Jufw3fO/0AcIhJJDS0mBTMrI9L6xziyFZAR+WHfI8axSQprrvMZwM7Sik6ORCS1tJgU3D1UCyORaCo+UM6df17b4vqm02K2Rk1RRY5dmH4KzfYQc/fC6Icjqaq6to4HX9vML17agONcMHEgL6/fRUV1XcM29dNiikjshKlTeK7R6xxgFLAeOCEmEUnKWVa4n9ufWs36nWWcO2EA37twAsN65zFvefERU2Aea69kETl2YVofndj4vZlNBm6IWUSSMg5W1jDrhbd56M3IdJcPfGEK504Y0LD+klOG8OhbkRtSFQWJdI5jbuAdjHp6WiyCkdTx6ju7ue2pVWwvreDq6SP51szx5Ku/gUjchalT+Gajt2lEZl7bHbOIpEsrq6jmrufW8djirYzu140nv3QGp47oFe+wRCQQ5qdZ41ZINUTqGJ6KTTjS1TTulTzlRy9SU+uUVlRzw4eP4xszxpGTqaEpRBJJmDqFH3ZGINL1NO2VvOdgFQbcPGMsX58xLr7BiUiz2pxtxMxeNLOCRu97mdmC2IYlXUFzvZIdeGJJUXwCEpE2hZmCql8w4xoA7r4f6N/K9iK4e4u9kluaQ1lE4i9MnUKtmQ2v76xmZiNoZvIbkXplFdXcNnd1i+vVK1kkcYVJCv8BvGZmrwTvP0QwHaakjrDDX6/bXsqNjyyjcN9hLpw0kIXrdlGuXskiSSNMRfMLQYe104kMhvcNd98T88gk6TyxZCt3zFtDz9xMHv3305k6qrd6JYskmTD9FAw4HzjO3e80s+FmNtXd34p9eJIMqmrq+OGf1/LIokLOGN2He688hX7dswH1ShZJNmGKj34D1AEfBe4Eyoj0U1CvZmFXaQVffmQZS9/bzw0fPo5vz/wA6Wktza0kIokuTFKY5u6TzWw5RFofmVlWjOOSJLBi6wGuf2gJZRU1/PqqU7hw0uB4hyQiHRQmKVSbWTpBiyMz60fkzkFS2J9XbuNbT6ykX/ds5t54BscP0pxLIl1BmKTwS+BpoL+Z3QVcDtwR06gkYbk7v3hpA/cu3MBpI3vxu8+dSp/87HiHJSJREqb10SNmthQ4h0jro0vcfV3MI5OEU+fON+es5OnlxVx+6lDuunQi2Rkau0ikK2ltjubejd7uAh5tvM7d98UyMEkckUHt9lNV6yzesp+PnziQWZdPItIwTUS6ktbuFJYSqUdo7i/fgeNiEpEklHnLi/nOU6uoqn2/E/vLb+/mmRXb1N9ApAtqMSm4+6jODEQS04+fX0dlzZHtCsqra5m1YH3opKD+CSLJI9RUV2Z2EZHhLQD+7u7Pxi4kSRTLC/ezq6yy2XUa1E6kawozdPbdwM3Av4LHzWb2k1gHJvH12oY9XPXAohY7oh3LoHYikjzCDJ19AXCuuz/o7g8SGfLi47ENS+Jp4bqdfHH2Ykb0yeP7nzie3Cazo2lQO5GuK+xM6QVAfWujnjGKRRLAc6u2c/NjyzlhcA9mf3EqBXlZ9MjJ0qB2IikiTFL4CbDczP5GpCXSh4DbYxqVxMW85cV8c84KTh3RiwevOY3uOZmABrUTSSVtFh+5+6NEhs2eGzymu/tjbe1nZg+a2S4zW9No2Q/MrNjMVgSPCxqtu93MNprZejOb2b7TkfZ6enkR35yzgmmj+jD7i1MbEoKIpJbWOq9NbrKofmLdwWY22N2XtfHZ/wv8GnioyfKfu/vPmhxrAnAlcAIwGHjJzMa5ey0SFa1NkvP08iJumbOSaaP68OA1p5GbpV7KIqmqteKjJcBaYHfwvnEzFCcylHaL3P1VMxsZMo6LgcfcvRLYbGYbganAGyH3l3Z6ZkWxEoKINGit+OgWoAQoB/4IfMLdPxI8Wk0IbbjJzFYFxUu9gmVDgK2NtikKlh3FzK43syVmtmT37t3NbSIhLVi7g2/OWcmUkb2VEEQEaCUpuPvP3f0s4CZgGLDQzOaY2ckdON5vgdHAycB24L+D5S0NpdFcXPe7+xR3n9KvX78OhJLaXnlnN1/903JOHNJTCUFEGoSpaN4MPAP8lUiRzrj2Hszdd7p7rbvXAQ8EnweRO4NhjTYdCmxr73GkdW9t3scNDy9hTP98Zl87lfzssC2TRaSrazEpmNlxZvZdM1sE/BBYCXzA3ee092BmNqjR20uB+pZJ84ErzSzbzEYBYwHNAR0D67aXct3sxQwuyOXh66bSM0+tjETkfa39RNwIrCJyl1AKDAdurB8u2d3vae2DzexR4Gygr5kVAd8Hzg6KnxzYAtwQfNZaM5tDZBiNGuArankUfRXVtXzhwbfolpXBw9dNO6bJcdQ/QSQ1tJYU7uT9cv38Y/1gd/9MM4v/0Mr2dwF3HetxJJzq2jre3lFGTmY6T35pOkM0dpGINKO1obN/0IlxSAzNWVzI8sIDONA9J5O120oZO6B7vMMSkQQUZkA8SWJzlxZx29zVDbd8ew5Wcvvc1cxbXhzXuEQkMSkpdGHuzn8+s4a6Jo176yfJERFpqrXWRzcHz2d2XjgSTX94bTOHqpqvr9ckOSLSnNbuFK4Nnn/VGYFIdL389k5+/Pw6cjKb/4o1SY6INKe1pLDOzLYA44NhKeofq81sVSfFJ+2wfkcZX3t0BccP6sGdF52gSXJEJLTWWh99xswGAguAizovJOmIvQcruW72YnKz0vn91VMY1DOXrIx0TZIjIqG0Or6Bu+8ATjKzLN4f3mK9u1fHPDI5ZtW1ddz4yDJ2l1Xy+A3TGdQzUkSkSXJEJKw2B70xsw8TmRNhC5GB64aZ2dXu/mqMY5Nj9OPn17Fo8z7u+dRJnDysIN7hiEgSCjMS2j3Aee6+HsDMxgGPAqfGMjA5WmsT5cxdVsQf/7mFa88cyScnD+3s0ESkiwjTTyGzPiEAuPs7gEZRSyBriku4fe5qTj+uN9+94Ph4hyMiSSzMncISM/sD8HDw/rPA0tiFJMei5HA1X35kKX26ZfHrqyaTma7+iCLSfmGSwpeBrwBfI1Kn8Crwm1gGJeG4O7c8sZIdJRU8fsN0+h7DqKciIs1pMykE8ybfEzwkgdz/6iZeWreT739iApOH92p1W7U6EpEwVNaQpBZv2cdPF6zn4ycO4pozRsY7HBHpIpQUktD+Q1V89U/LGdYrl7svO5H6iY9ERDqq1aRgZulmNquzgpG2uTu3PrmKvYcq+fVVk+meo4ZgIhI9bfVorjWzU83M3N1b21Zia97yYpYXHqCqtg6AS04ezMQhPeMclYh0NWFaHy0HnjGzJ4BD9QvdfW7MopIjzFtezO1zVzckBIAFa3cwb3mxxjASkagKU6fQG9gLfBT4RPC4MJZByZFmLVhPefWR8yKUV9dpohwRibowTVKvbWsbia2WJsTRRDkiEm1t3imY2TgzW2hma4L3k8zsjtiHJvV6d8tqdrkmyhGRaAtTfPQAcDtQDeDuq4ArYxmUvK/kcDXVtXU0bXSqiXJEJBbCJIU8d3+rybKaWAQjR/ve/DUcrqrlm+eOIysY12hIQS4/+eSJqmQWkagL0/poj5mNBhzAzC4Htsc0KgHghTU7eGbFNr4xYxxfPWcsr23cA2jIChGJnTBJ4SvA/cAHzKwY2ExkpFSJoX2Hqrhj3mpOGNyDGz8yOt7hiEiKCNP6aBMww8y6AWnuXhb7sOQH89dSUl7Nw9dN03DYItJpwrQ+6mNmvwT+AfzdzO41sz6xDy11vbBmB/NXbuOrHx3L8YN6xDscEUkhYX6CPgbsBi4DLg9ePx7LoFJZyeFq7pi3hgmDevDls1VsJCKdK0ydQm93/69G739kZpfEKqBU95O/rGP/4Sr+99rTVGwkIp0uzFXnb2Z2pZmlBY9PAc+1tZOZPWhmu+o7vQXLepvZi2a2IXjuFSw3M/ulmW00s1VmNrn9p5S83ty0l8cWb+Xfzhqlwe5EJC6spcFPzayMSDNUA7oB9aOxpQEH3b3Vwm4z+xBwEHjI3ScGy34K7HP3u83sNqCXu3/HzC4AvgpcAEwD7nX3aW0FP2XKFF+yZEmI00x8FdW1XHDvP6ipcxZ8/UPkZqXHOyQR6aLMbKm7T2luXYt3Cu7e3d17BM9p7p4RPNLaSgjB/q8C+5osvhiYHbyeDVzSaPlDHvEmUGBmg9o6Rlfy65c3smnPIe66dKISgojETZg6BcxsEjCy8fbtHDp7gLtvD/bfbmb9g+VDgK2NtisKlh3VSc7MrgeuBxg+fHg7Qkg8G3eVcd+r7/LJU4bwwbH94h2OiKSwNpOCmT0ITALW8n4RkgPRnE+hufkkmy3Xcvf7iXSmY8qUKUk38c+n73sDeL9Xsrvzn/PWkpuZznc/fnw8QxMRCXWncLq7T4jS8Xaa2aDgLmEQsCtYXgQMa7TdUGBblI6Z0Oav3MYbm/byo0sm0jc/O97hiEiKC9P66A0zi1ZSmA9cHby+Gnim0fIvBK2QTgdK6ouZurLSimp+9Nw6Thrak89M7RpFYSKS3MLcKcwmkhh2AJVEinrc3Se1tpOZPQqcDfQ1syLg+8DdwBwzuw4oBK4INn+eSMujjcBhICUm9rnnr++w52AlD159GulpzZWgiYh0rjBJ4UHg88Bq3q9TaJO7f6aFVec0s60TGXgvZby9o5SH3tjCZ6cN58Sh6pMgIokhTFIodPf5MY8khbg7P5z/L3rkZvKt8zRRjogkjjBJ4W0z+xPwZyLFR0C7m6SmtHnLi1leeICq2sgN1+WTh1CQ1/xUmyIi8RAmKeQSSQbnNVoW7SapXd685cXcPnd1Q0IAeHb1ds4a208zqIlIwggzn0JKVPrG2qwF6ymvrj1iWUV1HbMWrFdSEJGEEabz2h9ppiOZu38xJhF1UdsOlB/TchGReAhTfPRso9c5wKWkSMeyaBpckEtxMwlgcEFuHKIREWlemOKjpxq/D/ofvBSziLqoq6YNZ9aC9Ucsy81M59aZan0kIokj1IB4TYwF1P32GLg7b7y7l9zMNGrrnKpaZ0hBLrfOHK/6BBFJKGHqFOrnVai3A/hOzCLqgv7+zm5e27iH7104gQVrdwDvD4gnIpJIwhQfde+MQLqq2jrnJ8+vY2SfPD53+oiGpCAikohaTApm1moRkbsXRj+crmfusiLe2XmQ33x2MlkZmnNZRBJba3cKz/H+dJz1HOgH9Ac0PVgbKqpr+cVLG5g0tCcfmzgw3uGIiLSpxaTg7ic2fm9mI4nUJcwAfhzTqLqIRxYVUnygnJ9ePgkzjYIqIokvTEXzWOA/gGnAfwNfc/fqWAeW7Moqqvmfv23krDF9OXNM34blqmAWkUTWWp3CRCLJ4ATgp8B17l7b0vZypAf+sZl9h6rUD0FEkkprdworga1E6hamAlMbF4G4+9diG1ry2nuwkt//YxMfmziQk4YVxDscEZHQWksKGtuone5/dRPl1bXcct64eIciInJMWqtont2ZgXQVu8sqmf3GFi4+aTBj+quLh4gkFzWcj7LfvfIu1bXOzTN0lyAiyUdJIYp2llbwf2++x6WnDGFU327xDkdE5JgpKUTRb/62kdo652sfHRvvUERE2qXNpGBm48xsoZmtCd5PMrM7Yh9actlRUsGjb23l8lOHMrxPXrzDERFplzB3Cg8AtwPVAO6+CrgylkElo/tefZdad77ykTHxDkVEpN3CJIU8d3+rybKaWASTrHaXVfKnRYVcesoQhvXWXYKIJK8wSWGPmY0mmFPBzC4Htsc0qiTz+9c2UV1bx41nj453KCIiHRJm5rWvAPcDHzCzYmAz8LmYRpVE9h+q4v/eeI8LJw3muH758Q5HRKRDwkyyswmYYWbdgDR3L4t9WMnjj//czKGqWm76qOoSRCT5hRklNRu4DBgJZNSPf+Tud8Y0siRQVlHNH1/fwvknDGTcAPVeFpHkF6b46BmgBFgKVMY2nOTyyKJCyipq1OJIRLqMMElhqLufH/NIksin73uDOnfe23uYs8b05cShPeMdkohIVIRpffS6mZ3Y9mbhmdkWM1ttZivMbEmwrLeZvWhmG4LnXtE8ZrTtOVjJrrJKvvRhtTgSka6jxaRgZmvMbBVwFrDMzNab2argYr4qCsf+iLuf7O5Tgve3AQvdfSywMHifkNyd7SUVTBzSgzPH9Il3OCIiUdNa8dEQ4OTOCgS4GDg7eD0b+DuROaETyrzlxSwrPEBNnVO8v5xnVmzjklOGxDssEZGoaC0pbHb392J0XAf+amYO3Ofu9wMD3H07gLtvN4g8AN8AAAoASURBVLP+ze1oZtcD1wMMHz48RuE1b97yYm6fu4qaOgdg/+Fqbp+7GkCJQUS6hNaSQn8z+2ZLK939ng4c90x33xZc+F80s7fD7hgkkPsBpkyZ4h2I4ZjNWrCe8uq6I5aVV9cya8F6JQUR6RJaSwrpQD5grWzTLu6+LXjeZWZPE5kDeqeZDQruEgYBu6J93I7adqD8mJaLiCSb1pLC9lh0UGvcMzp4fR5wJzAfuBq4O3h+JtrH7qj+PbLZWXp0V43BBblxiEZEJPpaSwpRv0MIDACeDnpGZwB/cvcXzGwxMMfMrgMKgStidPx2G9c//6ikkJuZzq0zx8cpIhGR6GotKZwTiwMGYymd1MzyvbE6ZjSUHK5maeEBThvRi5VFJVTV1jGkIJdbZ45XfYKIdBktJgV339eZgSS6RxcXcriqlh9ePJEf/nktAI/fMD3OUYmIRFeYYS5SXnVtHf/7zy2cMboPEwb3UDIQkS4rzDAXKe+FNTvYUVrBdWeNincoIiIxpaQQwuzXtzCiTx4fGd9sfzoRkS5DSaENa4pLWPLefj5/+gjS0mLVIEtEJDEoKbRh9utbyM1M54opw+IdiohIzCkptGLfoSqeWbmNSycPoWduZrzDERGJOSWFVjy2uJCqmjqunj4y3qGIiHQKJYUW1NTW8cibhUw/rg/jB2r+ZRFJDUoKLXj57V0UHyjn6jNGxDsUEZFOo6TQgkcWFTKgRzYzjh8Q71BERDqNkkIzCvce5tUNu7nytOFkpOufSERSh654zXh0cSEGXDlVzVBFJLUoKTRRVVPHnMVbOef4AQzqqXkSRCS1KCk0sWDtDvYequKz0zp3/mcRkUSgpNDEI4veY2ivXD40tl+8QxER6XRKCo28u/sgb27ax1XThmucIxFJSUoKjTy+eCsZacblpw6NdygiInGhpBCorq1j7rIiPvqB/vTvnhPvcERE4kJJIbBw3S72HKxSM1QRSWlKCoHHF0d6MKuCWURSmZICsL2knFfe2c0Vpw5TD2YRSWm6AgJPLimizuFTmkhHRFJcyieFujpnztKtnDG6D8P75MU7HBGRuEr5pPDmpr1s3VfOp0/TXYKISMonhSeXFtE9J4OZJwyMdygiInGX0knhYGUNf1mzgwsnDSInMz3e4YiIxF1KJ4W/rN5OeXUtl01WD2YREUjxpPDUsiJG9snj1BG94h2KiEhCSNmksHXfYd7ctI/LJg/FTIPfiYgAZMQ7gKbM7HzgXiAd+L273x3Nz5+3vJhZC9ZTfKAcgG7ZqksQEamXUHcKZpYO/A/wMWAC8BkzmxCtz5+3vJjb565uSAgAsxa8w7zlxdE6hIhIUkuopABMBTa6+yZ3rwIeAy6O1ofPWrCe8uraI5aVV9cya8H6aB1CRCSpJVpSGAJsbfS+KFjWwMyuN7MlZrZk9+7dx/Th2xrdIYRZLiKSahItKTRX4+tHvHG/392nuPuUfv2ObUTTwQW5x7RcRCTVJFpSKAIajzcxFNgWrQ+/deZ4cpt0UsvNTOfWmeOjdQgRkaSWaK2PFgNjzWwUUAxcCVwVrQ+/5JRISdSsBevZdqCcwQW53DpzfMNyEZFUl1BJwd1rzOwmYAGRJqkPuvvaaB7jklOGKAmIiLQgoZICgLs/Dzwf7zhERFJRotUpiIhIHCkpiIhIAyUFERFpoKQgIiINzN3b3ipBmdlu4L127t4X2BPFcOJJ55KYusq5dJXzAJ1LvRHu3mzv36ROCh1hZkvcfUq844gGnUti6irn0lXOA3QuYaj4SEREGigpiIhIg1ROCvfHO4Ao0rkkpq5yLl3lPEDn0qaUrVMQEZGjpfKdgoiINKGkICIiDbp8UjCz881svZltNLPbmlmfbWaPB+sXmdnIzo8ynBDnco2Z7TazFcHj3+IRZ1vM7EEz22Vma1pYb2b2y+A8V5nZ5M6OMawQ53K2mZU0+k6+19kxhmFmw8zsb2a2zszWmtnNzWyTFN9LyHNJlu8lx8zeMrOVwbn8sJltonsNc/cu+yAy/Pa7wHFAFrASmNBkmxuB3wWvrwQej3fcHTiXa4BfxzvWEOfyIWAysKaF9RcAfyEyE9/pwKJ4x9yBczkbeDbecYY4j0HA5OB1d+CdZv5/JcX3EvJckuV7MSA/eJ0JLAJOb7JNVK9hXf1OYSqw0d03uXsV8BhwcZNtLgZmB6+fBM4xs+amBY23MOeSFNz9VWBfK5tcDDzkEW8CBWY2qHOiOzYhziUpuPt2d18WvC4D1tFkfnSS5HsJeS5JIfi3Phi8zQweTVsHRfUa1tWTwhBga6P3RRz9n6NhG3evAUqAPp0S3bEJcy4AlwW39k+a2bBm1ieDsOeaLKYHt/9/MbMT4h1MW4Lih1OI/CptLOm+l1bOBZLkezGzdDNbAewCXnT3Fr+XaFzDunpSaC5bNs2yYbZJBGHi/DMw0t0nAS/x/q+HZJMs30kYy4iMM3MS8CtgXpzjaZWZ5QNPAV9399Kmq5vZJWG/lzbOJWm+F3evdfeTicxZP9XMJjbZJKrfS1dPCkVA41/LQ4FtLW1jZhlATxKzOKDNc3H3ve5eGbx9ADi1k2KLtjDfW1Jw99L623+PzCqYaWZ94xxWs8wsk8hF9BF3n9vMJknzvbR1Lsn0vdRz9wPA34Hzm6yK6jWsqyeFxcBYMxtlZllEKmHmN9lmPnB18Ppy4GUPamwSTJvn0qR89yIiZanJaD7whaC1y+lAibtvj3dQ7WFmA+vLd81sKpG/ub3xjepoQYx/ANa5+z0tbJYU30uYc0mi76WfmRUEr3OBGcDbTTaL6jUs4eZojiZ3rzGzm4AFRFrvPOjua83sTmCJu88n8p/nYTPbSCS7Xhm/iFsW8ly+ZmYXATVEzuWauAXcCjN7lEjrj75mVgR8n0gFGu7+OyJzdF8AbAQOA9fGJ9K2hTiXy4Evm1kNUA5cmaA/Os4EPg+sDsqvAb4LDIek+17CnEuyfC+DgNlmlk4kcc1x92djeQ3TMBciItKgqxcfiYjIMVBSEBGRBkoKIiLSQElBREQaKCmIiEiDLt0kVSRazKwPsDB4OxCoBXYH7w+7+xlxCUwkytQkVeQYmdkPgIPu/rN4xyISbSo+EukgMzsYPJ9tZq+Y2Rwze8fM7jazzwbj4a82s9HBdv3M7CkzWxw8zozvGYi8T0lBJLpOAm4GTiTSq3acu08Ffg98NdjmXuDn7n4acFmwTiQhqE5BJLoW148HZGbvAn8Nlq8GPhK8ngFMaDTkfQ8z6x6M/S8SV0oKItFV2eh1XaP3dbz/95YGTHf38s4MTCQMFR+JdL6/AjfVvzGzk+MYi8gRlBREOt/XgCnBDHn/Ar4U74BE6qlJqoiINNCdgoiINFBSEBGRBkoKIiLSQElBREQaKCmIiEgDJQUREWmgpCAiIg3+P/wrXPBVsbB/AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret2 = ensemble_simulations(\n", " duration, ndiv=20, y0=y0, volume=V, model=m, solver=('spatiocyte', radius), repeat=N)\n", "ret2.plot('o', ret1, '-')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Simulating with `egfrd`:" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEGCAYAAACKB4k+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXyU9b33/9cneyCQCIQtbK5QFERERGmr52iLentXrbZqW6t20d5Hu53z8xR7enexG/fR2vXUIz211baKtiJatVLX2lorAkFAFomyJgHClkD2zHx+f8yVGGCSXCGZTCbzfj4e88hc11zL52L0+sz1Xc3dERERAchIdgAiItJ/KCmIiEgbJQUREWmjpCAiIm2UFEREpE1WsgPoiREjRvikSZOSHYaISEpZsWLFHncvjvdZSieFSZMmsXz58mSHISKSUsxsa0efqfhIRETaKCmIiEgbJQUREWmjpCAiIm2UFEREpI2SgoiItFFSEBGRNkoKIiLSRklBRCQJrr73Va6+99Vkh3EUJQUREWmjpCAikmIS+ZShpCAiIm0SlhTMbLyZvWhm683sTTP7YrD+m2ZWbmargtcl7fa53czKzGyjmc1LVGwiIj3VX+sEeiqRo6S2AP/m7ivNbAiwwsyeDT77obvf1X5jM5sKXAOcCowFnjOzU9w9ksAYRUSknYQ9Kbh7pbuvDN4fBNYDJZ3schmwyN0b3X0zUAbMTlR8IiJytD6pUzCzScAZwGvBqlvNbLWZ3WdmxwXrSoDt7XbbQZwkYmY3mdlyM1teVVWVwKhFJBX0pBhnoBYB9UTCk4KZFQCPAl9y9xrgHuBEYAZQCfygddM4u/tRK9wXuvssd59VXBx34iAR6WM9vbnq5tx/JDQpmFk2sYTwO3dfDODuu9w94u5R4Be8W0S0AxjfbvdxQEUi4xMRkcMlsvWRAb8E1rv73e3Wj2m32RXA2uD9E8A1ZpZrZscDJwPLEhWfiIgcLZGtj+YC1wFrzGxVsO6rwLVmNoNY0dAW4GYAd3/TzB4B1hFruXSLWh6JiPSthCUFd/8b8esJnu5kn+8C301UTCIi0jn1aBYR6WNLSssp3XaA1zbvY+6CF1hSWt4n+4ahpCAigFoA9ZUlpeXcvngNTZEoAOUH6rl98ZpQN/ee7BuWkoKISB+6c+lG6psPry6tb45w59KNCd03rERWNIuIDEitRThNkShzF7zAbfMmc/kZ7/a1jUSduqYW6psi1DcHr6bYq/xAfdxjlh+o59tPrqOhOUJDc5SGlgiNzVEa2/3taN+KDtYfCyUFkQGktfjn4ZvPSXIkqaUlEuVQYwsHG1o41Bi8Glo42NhCbfA6FPxdW17Nss37iXisb235gXq+/PAqvvnEm2BQ1xShqSV6THEsWraNvOxM8rIzyc3OIDcrk7zsDHKzMigclENeVgYNcY49tii/R9ffnpKCiKSs9r/Yz/3+89x03gmcffxwDtQ1U13fRHV9c9urpr6FmoZmauqbqWloYce+OnYdbATgpP/4U6jz5Wdn0tgSIXrEWAsONLZEuerMcQzKySQ/JzP4m8Wg7Nj7vJxM8rMzWbZ5L//14ts0tru552VnsODD0w972ujoem9fvOawIqT87Exumzc53D9YCEoKIv1Muv3aj1cUc+n0Meyva2ZvbSN7DzWxt7aJfYca2VfXzP7aJvbVNfHWzhrKdte2jYVTUd3AN59YF/ccmRnGkLwsCvOzGZqXTUNzC1WHGg/bJjvTuHLmOM47pZiCvCwG52YxJDf2d3BuFgW5WWRmGMfPfyruORqaI3z78tO6vN45JwxnwrDB/PsfVtMUiVJSlH9U8VNHWrc5ln3DUlIQkR7pqny9lbtT09BC1cEGdtU0svtgAy9u2M3Ta3bSEn23KOZLD6/iSw+vOmp/ADMoys/muEE57Nhff/TgaMCwQdn87GMzKRyUTWF+7FWQm0VskIWYuQteOOrXfnPE+eumPSy4cnqn1zu2KD9u2X53inAuP6OEh5ZtA7qf/HuybxhKCiJyzOI1kfz3P7zBy5uqGDkkj53V9VRWN7CrJpYIjmw505EheVncNm8ywwbnMHxwLiMKchg2OIeiQTlkZsRu7h39Yt9f18y5J43o9PgdVcyGqbC9bd7khBfhJJOSgoh0KRp1qg41sn1fHTv217Njf+zvY6Xlh5WNAzRFnMUry8nJzGBUYS6jh+ZxWkkhF74nj1FD8xg5NJeRQ2J/L/zBX+L+2j/U0MInz5nUaUw9+cXek337oggnmZQURASIFe80tkR5aeNutu6tY8veWrbtrWPrvjq276s76uY/oiD3qHWtDNj4nYsOK7KJpyc35578Yu/pr/1EF+Ekk5KCSJrZX9vE21WHglct71Qd4p09tWyuilXa3vCr1wEYlJPJhGGDOLF4MOefUsyE4YMYP2wQ44/Lp6RoEPk5mcxd8EKHN/WuEgL07Obck1/sA/3Xfk8oKYgMUNV1zWzYWcNbuw/x1s6DbNp9kLLdh9hzqKltm5ysDI4fPpjJo4bQ0BQhLzuT/3fVdCYOH0RxQW6XN/be+MUNx35z7s8VtqlKSUEkAfqyWWkk6mzec4h1lQfZtq+OuqYIc773PDtrGtq2GZKbxUmjCrhgyihOHlXAicUFnDSygLFF+WRmGEtKy3l+/W6aIlG+tGhVn/7i1s25f1FSEEkhLZEom3YfYs2OatZWVLO2vJr1lQfbfqkbkJ+TyQVTRjJ59BBOGT2EKaOHMHpoXoe/+jsaZA0InRh0Ux84lBRE4ugPHcjcnR376yndfoA3th9g1fYDvFlRTUNz7OZdkJvF1LFDuXb2BKaOHcqpY4fy9cfXkmHG3VfPCH2ezgZZUxl7+lFSEOknGlsirC2vpuJAPYcaWzjru8+zJ+h1m5uVwbSSQj42eyKnjy9kWkkhk4YPJiPj8F//GSEqd4/Ukzb7MvAoKYgkSW1jCyu27mfZ5n0s27yPVTsOHDaQWn52JledOY4bzp3E5NFDyM7sfKT7sD2Lj9QbPXRl4FBSEOkjDc0RVmzdz9/f3sOrb+/ljR3VRKJOZoZx2tihzD1xGK+U7aUpEuvOVd8c4anVlbz3pBGcVlLY6bF7Ui8w0HvoDkSJLNZUUhBJEHdnXUUNL2+q4m+b9vD6ln00tkTJzDCmjyvk5vefwJwThjNz4nEU5GYxd8ELbQmhVdiy/Z7UC6jNfnL010p5JQWRXlRd38xfN1XxdtUhquubueQnfwVg8qghfGLOROaeNJzZxw+nIPfo//V6Urbf03oBtSDqvoH676SkINJDW/bU8uy6XTy3fhfLt+4nEnWyMozC/GzmXzyF959SzKiheV0eJ1lj+Yi0p6Qg0k3uzpryapa+uZOlb+6ibPchAKaMHsLnzjuBrAzjnpfeYW9tEz96bhPZmRmhimKSOZaPSCslBRmwerOvQTTqlG4/wNNrKnlm7U7KD9STmWHMOWEYnzh7Ahe8ZxTjhw3qUYVvOo/l05PvKFn7DlRKCiIdcHfWllfzxBsVPPlGBRXVDeRkZvC+k0fw5Q+cwgVTRnLc4JzD9ulpRzCN5SPJ1mVSMLOPAM+4+0Ez+xowE/iOu69MeHQiSbB9Xx3l++vZU9vIpT/9G1kZxvtPKeb/mzeZC6eOYmhedof7pmtHMCWhgSPMk8L/dfffm9l7gXnAXcA9wNkJjUykD9U2tvD0mkp+v2IHyzbva1tflJ/NbfMm8/E5E0MdRxW+kurCJIXWZ+H/Bdzj7o+b2TcTF5JI33CP1RMsWraNJ1dXUtcUYURBDlkZ1jZn8IH6Zr7z1HoG52YlvLJYpD8IkxTKzexe4ELg/5lZLtB5f3uRfqymoZklpeU8+No2Nuw8yKCcTC6dPoaPzhrPFxeVtiWEVt2tE4DUrfAVCZMUPgpcBNzl7gfMbAxwW1c7mdl44AFgNBAFFrr7j81sGPAwMAnYAnzU3fdbbFzfHwOXAHXADaq3kN60YWcND7y6lSWl5dQ1RZhWUsj3PzyN/3362LbOZBUHGuLu2506AVX4SirrMim4e52Z7QbeC2wCWoK/XWkB/s3dV5rZEGCFmT0L3AA87+4LzGw+MB/4CnAxcHLwOhvVWwg9b1bq7ix9cye/emUz/3hnH7lZGXzo9LFcd85Epo8rOmr7dK4TUAITCNf66BvALGAy8CsgG/gtMLez/dy9EqgM3h80s/VACXAZcH6w2f3AS8SSwmXAA+7uwD/MrMjMxgTHEemW2sYWdlY3sLOmgZt/s4KSonxuv3gKH501/qhmpO2pTkDSXZjioyuAM4CVAO5eEfzyD83MJgXHeA0Y1Xqjd/dKMxsZbFYCbG+3245g3WFJwcxuAm4CmDBhQnfCkDRQdbCR+/++hd/8YyvV9c0U5Gbxo6tn8IGpo8jqYuhpUJ2ASJik0OTubmYOYGaDu3MCMysAHgW+5O41nUwEHu8DP2qF+0JgIcCsWbOO+lzSU/mBehb+5W0Wvb6dpkiU08YOZX1lC4caW/jOU+tpbIn2yWTwIqkuTFJ4JGh9VGRmnwU+BfwizMHNLJtYQviduy8OVu9qLRYKKq13B+t3AOPb7T4OqAhzHklf2/fV8bMXynh05Q4APjyzhBNHFvCjZze1tSLq7pzDyaZEJMkUpqL5LjP7AFBDrF7h6+7+bFf7Ba2Jfgmsd/e72330BHA9sCD4+3i79bea2SJiFczVqk+QjrRPBhkZxsfPnsBN551ISVE+cxe8oDmHRY5RqLGPgiTQZSI4wlzgOmCNma0K1n2VWDJ4xMw+DWwDPhJ89jSx5qhlxJqk3tjN80ka2F3TwM9eLOOhZdswMz4xZyKfO+9ERhe+OzR1ug41AXrKkJ7rMCmY2UFiZfrG4WX7Bri7D+3swO7+N+LXEwBcEGd7B27pKmBJTzUNzfz3S29z3yubaYk4Hz1rPJ//55MYU3h0U9H+0KxUN2dJVR0mBXfvVgsjkY4ca1+D9hPRz/jWn4k6XDZjLP/6gVOYOLzj9g5qVipy7ML0U4jb7tPdt/V+OCIxj63cwVceXd02Z3HUITcrg3+aPLLThABqVirSE2HqFJ5q9z4POB7YCJyakIgk7W3YWcP8xWuOmsS+sSXaJ/MSiKSzMK2PprVfNrOZwM0Ji0jSVk1DMz989i0eeHUrkWj8LijpUFkskkzdnnktGMvorEQEI+nJ3XmstJzvPb2BvbWNfGz2BF7YsJvK6qMHp0uHMYhEkilMncK/tlvMIDbzWlXCIpK0snlPLV9bsoZXyvYyY3wRv7rhLKaNK2yb6zhZlcUqcpJ0FeZJoX0rpBZidQyPJiYcSRfNkSj3/uVtfvJCGblZGXzn8tP42OwJZGTEWjGrslgkOcLUKXyrLwKR9LG2vJp//8Nq1lXW8L+mjeEb/3sqI4fmHbWdKotF+l6Y4qNngY+4+4Fg+ThgkbvPS3RwMrA0tUT5yfObuOcvbzNscA73Xncm804dneywRKSdMMVHxa0JASCYJW1kZzuItGrfAW3q15+hJepcOXMcX790KoWDspMdnogcIcxcy5H2HdjMbCJxhrQWOVKssjhWJwDQEnVyMjN438kjlBBE+qkwTwr/AfzNzP4SLL+fYJIbkc4s+NN66pujh61rioTvgNZTqocQ6b4wFc3PBB3W5hAb4O7L7r4n4ZFJSntm7U521jTG/Uwd0ET6ry6Lj4J5ES4CZrr7H4FBZjY74ZFJv3L1va+2DWzXmcaWCN94fC2f++0KsjPjD5KrDmgi/VeYOoWfA+cA1wbLB4H/SlhEkrK27q3lqnte5f5Xt/KZ9x7P96+YRn525mHbaLRSkf4tTJ3C2e4+08xKoa31UU6C45IU8+y6Xfzrw6swg4XXnckHg6amWZkZPeqApnoBkb4VJik0m1kmQYsjMysGop3vIukiEnV++Oxb/OzFMqaVFPLzj89k/LBBbZ+rA5pIagmTFH4CPAaMNLPvAlcBX0toVJISquua+fyiUl5+q4qPzhrHHZedRt4RxUUiklrCtD76nZmtIDaFpgGXu/v6hEcm/VrZ7oN85v7llB+o53tXTONjZ8edi0lEUkxnczQPa7e4G3io/Wfuvi+RgUn/9eKG3XzhoVJyszN46LNzmDVpWNc7iUhK6OxJYQWxeoR47QodOCEhEUm/5e7c97fNfPupdUwdM5SFn5xFiZqXigwoHSYFdz++LwOR/qv9+EXLtuxnWslQHr75HAbldHuOJhHp58L0U8DMPmRmdwWvSxMdlPQfS0rLmd9u/CKAst2H+PObu5IYlYgkSpgezQuALwLrgtcXzez7iQ5M+ocFz2yg4Yjxi+qbY+MXicjAE+b5/xJghrtHAczsfqAUuD2RgUnybdtbx8448ySDxi8SGajCFgoXAa2tjQoTFIv0I+sqavjkfcswA48zUHp3xi9SpzWR1BGmTuH7QKmZ/Tp4SlgBfC+xYUkyrdi6j6sXvkp2pjH/oskav0gkjYTpvPaQmb0EnEWseepX3H1nogOT5Pjbpj189oHljC7M47efOZuSonxGDc3v0fhFIpI6Ouu8NvOIVTuCv2PNbKy7r0xcWJIMz63bxb/8biUnFA/mN58+m+IhuYDGLxJJJ509KSwH3gSqguX2ndgc+OfODmxm9wGXArvd/bRg3TeBz7Y75lfd/engs9uBTwMR4AvuvrRbVyI98szandz64EpOLSnkgRtna7pMkTTVWVL4N+BKoB5YBDzm7oe6cexfAz8DHjhi/Q/d/a72K8xsKnANcCowFnjOzE5x90g3ziedaJ0gJ94v/adWV/KFRaWcPq6QX39qNkPzlBBE0lWHFc3u/kN3fy9wKzAeeN7MHjGzGWEO7O4v826Lpa5cBixy90Z33wyUAZrdrQ+0JoSZE4p44NNnKyGIpLkuWx8FN+nHgT8Tu1Gf0sNz3mpmq83sPjM7LlhXAmxvt82OYN1RzOwmM1tuZsurqqribSIhLX1zZ1tC+PWNsynI1bAVIumuw6RgZieY2VfN7DXgW8AbwBR3f6QH57sHOBGYAVQCP2g9XZxt47SOB3df6O6z3H1WcXFxD0JJby9s2MWtD65kWkkh991wFoOVEESEzusUyoDVxJ4SaoAJwL+Yxe7f7n53d0/m7m0D5pjZL4Ang8UdxIqoWo0DKrp7fAnnlbI9fO63K5kyeij3f2o2Q1RkJCKBzpLCHbz7a72gN05mZmPcvTJYvAJYG7x/AnjQzO4mVtF8MrCsN84phyvdtp/PPrCc44cP5oFPzaYwP1xCUFNUkfTQ2dDZ3+zJgc3sIeB8YISZ7QC+AZwfVFQ7sAW4OTjXm2b2CLEB91qAW9TyqPe0H/r6w/f8neMGZfObT8/muME5yQ5NRPqZhBUku/u1cVb/spPtvwt8N1HxpKslpeXcvnhN29DX7lDbGOHvb+9Vr2QROUqo+RQkdd25dCP1zYc/dDW2aOhrEYmvs9ZHXwz+zu27cKS3dTTEtYa+FpF4OntSuDH4+9O+CER6XyTq5GbF/4q7M/S1iKSPzuoU1pvZFqDYzFa3W2+Au/v0hEYmPXbHH9+koSVKdqbRHHm324eGvhaRjnTW+uhaMxsNLAU+1HchSW/49Subuf/VrXz2fcdz6thCDX0tIqF02voomDfhdDPL4d3hLTa6e3PCI5Nj9uKG3dzx5Do+MHUU8y9+D5kZpqGvRSSULpukmtl5xEY63UKs6Gi8mV0fDHgn/cz6yhpufXAl7xkzlB9fM4PMjHgjiIiIxBemn8LdwAfdfSOAmZ0CPAScmcjApPv21TbxmfuXU5CXxS+vP4tBORrPSES6J0w/hezWhADg7m8BGiynn2mORLnldyupOtTIvdfNYnRhXrJDEpEUFOan5HIz+yXwm2D548CKxIUkx+K7T63n1Xf28oOPnM6M8UXJDkdEUlSYpPB/gFuALxCrU3gZ+Hkig5LueWT5dn799y185r3Hc+WZ4+JuowpmEQmjy6Tg7o3E6hW6PVS2JN7a8mq+tmQtc08azvyLpyQ7HBFJcRr7KIVcfe+rbXMtAxyoa+Jzv13BiME5/OSaM8jK1NcpIj2j5ikpKhp1vvzwKnbVNPDIzecwvCA32SGJyADQ6U9LM8s0szv7KhgJ7+cvlfHixiq+fulUzphwXNc7iIiE0GlSCCa6OdNa5+CUfuG1d/Zy97NvcdmMsXxizsRkhyMiA0iY4qNS4HEz+z1Q27rS3RcnLCrpUHMkyhcWlTJx+GC+e8U0lK9FpDeFSQrDgL3AP7db54CSQh9zd96uqqW+OcJ9N5xFQa6qhESkd4VpknpjV9tI4i0pLWfFtgNEok5hfjabdh3i1LGFyQ5LRAaYLtswmtkpZva8ma0Nlqeb2dcSH5q0WlJazlceXU0kGpsTobq+mdsXr2FJaXmSIxORgSZMw/ZfALcDzQDuvhq4JpFByeH+85kNNLZED1tX3xzRPMsi0uvCJIVB7r7siHUtiQhG4quoboi/XvMsi0gvC5MU9pjZicQqlzGzq4DKhEYlbZ5fv6vDzzTPsoj0tjDNV24BFgJTzKwc2ExspFRJsH21TXzl0dWMLcxjX10TDc3vFiFpnmURSYQwrY/eAS40s8FAhrsfTHxY4u58bckaquub+ePn38uGyoOaZ1lEEi7MdJzDgW8A7wXczP4G3OHuexMdXDr74+pKnl6zk3+/aDJTRg9lyuihmmdZRBIuTJ3CIqAKuBK4Knj/cCKDSne7axr4v0vWcsaEIm563wnJDkdE0kioHs3u/u12y98xs8sTFVC6c3duX7yGhuYId33kdA2HLSJ9Kswd50Uzu8bMMoLXR4GnEh1YunrijQqe37Cb2+ZN5sTigmSHIyJppsOkYGYHzawGuBl4EGgKXouAL3d1YDO7z8x2t/aEDtYNM7NnzWxT8Pe4YL2Z2U/MrMzMVpvZzJ5eWCrae6iRb/1xHaePL+LGuccnOxwRSUMdJgV3H+LuQ4O/Ge6eFbwy3H1oiGP/GrjoiHXzgefd/WTg+WAZ4GLg5OB1E3BPdy9kILjjyXUcbGjmP6+cTmbG0aOfPnzzOapkFpGECjXMpplNBya1376robPd/WUzm3TE6suA84P39wMvAV8J1j/g7g78w8yKzGyMu6dNJ7kXNuzi8VUVfOnCk5k8ekiywxGRNBWmSep9wHTgTaC199SxDp09qvVG7+6VZjYyWF8CbG+33Y5gXVokhdrGFr722FpOGVXAv5x/UrLDEZE0FuZJYY67T01wHPFmivG4G5rdRKyIiQkTJiQypj7zo+feoqK6gUc/dg45WWptJCLJE+YO9KqZ9VZS2GVmYwCCv7uD9TuA8e22GwdUxDuAuy9091nuPqu4uLiXwkqedRU13PfKFq6dPYEzJw5LdjgikubCJIX7iSWGjUHLoDVmtvoYz/cEcH3w/nrg8XbrPxm0QpoDVKdDfUIk6nz1sTUU5Wcz/6IpyQ5HRCRU8dF9wHXAGt6tU+iSmT1ErFJ5hJntIDZUxgLgETP7NLAN+Eiw+dPAJUAZUAekxWxvDy3bxqrtB/jh1adTOCg72eGIiIRKCtvc/YnuHtjdr+3gowvibOvERmMd8K6+91UAfv7xmfznMxs498ThXD5DA9uJSP8QJilsMLMHgT8Cja0ru2qSKp27c+lG6poi3HHZqZjFq2cXEel7YZJCPrFk8MF26461SaoAhxpbeHj5dj77vhM4aaT6JIhI/xFmPoW0KN/vK+7Olr21jCjI5fP/rD4JItK/hOm89ivi9Blw908lJKIBbElpOSu2HSASdXIyM3h+/W5NlCMi/UqY4qMn273PA66ggz4E0rElpeXMX7yaSDSWX/fXNXP74jUASgwi0m+EKT56tP1y0NT0uYRFNEDduXTjYXMsA9Q3R7hz6UYlBRHpN45lTIWTgYExvkQfqjhQ3631IiLJEKZO4SCH1ynsJDayqXRDbnbGUU8KAGOL8pMQjYhIfGGKj9RmsoeWbd5HQ3OUrAyjJfpufs3PzuS2eZOTGJmIyOE6m3ltQmevvgwylUWjzneeWseYwjy+d8U0coI5l0uK8vn+h6epPkFE+pXOnhSeIlZs1L67rQPFwEggM4FxDRhPrqlk9Y5qfvCR07nyzHE8unIHgGZQE5F+qcOk4O7T2i8Hs6h9BbgQ+F5Coxogmlqi3LV0I1NGD9ETgYikhC5bH5nZyWb2a+BPwApgqrv/NNGBDQQPvraVbfvqmH/xlLhzLouI9DcdPimY2WnAfwCnAv8JfNrdI30VWKo72NDMT14o49wTh3PeKak/GZCIpIfO6hTeIDZv8lPAbGB2+9E83f0LiQ0ttS18+R321TYx/+Iph42CqroEEenPOksKGtvoGFUdbOR//rqZS6ePYfq4omSHIyISWmcVzff3ZSADyc9fKqMpEuXfPqg+CCKSWo5lmAvpRMWBen73j21cNXMcx48YnOxwRES6RUmhl/30hTIc5/MXaK4EEUk9Sgq9aOveWn6/fDsfmz2BcccNSnY4IiLdFqafwilm9ryZrQ2Wp5vZ1xIfWur50XObyMo0bvknPSWISGoK86TwC+B2oBnA3VcD1yQyqFRUtvsQS1aV88lzJjFyaF6ywxEROSZhksIgd192xLqWRASTyn72wibysjK5+f0nJDsUEZFjFiYp7DGzEwnmVDCzq4DKhEaVYt6pOsQTb1Rw3TkTGV6Qm+xwRESOWZg5mm8BFgJTzKwc2Ax8IqFRpZifvVhGTlYGn32fnhJEJLWFmWTnHeBCMxsMZLj7wcSHlTq27Knl8VUV3HDuJIqH6ClBRFJbmOk4c4ErgUlAVus4Pu5+R0IjSxH/9WIZWRnGzefpKUFEUl+Y4qPHgWpiw2Y3Jjac1LJ9Xx2LS8u5bs5ERg5RiyMRSX1hksI4d78o4ZGkoIUvv0OGoacEERkwwrQ++ruZTet6s/RSdbCRR5Zv58NnjGNMYX6ywxER6RWdTbKzFogG29xoZu8QKz4ywN19+rGe1My2AAeBCNDi7rPMbBjwMLG6iy3AR919/7GeI9Hue2UzTZGonhJEZEDprPioBJiRwHP/k7vvabc8H3je3ReY2fxg+SsJPP8xu/KeV1i1rZpLpo3hhOKCZIcjItJrOis+2uzuWzt6JSCWy4DWORzuBy5PwDl6bNEKgRIAAAosSURBVElpOau2VxNxZ9nmfSwpLU92SCIivaazJ4WRZvavHX3o7nf34LwO/NnMHLjX3RcCo9y9Mjh2pZmNjLejmd0E3AQwYcKEHoTQfUtKy5m/eDWRqAOxeoXbF68B4PIzSvo0FhGRROjsSSETKACGdPDqibnuPhO4GLjFzN4fdkd3X+jus9x9VnFxcQ/D6J47l26koTl62Lr65gh3Lt3Yp3GIiCRKZ08KlYnqoObuFcHf3Wb2GDAb2GVmY4KnhDHA7kScuycqDtR3a72ISKrp7EnBEnFCMxtsZkNa3wMfBNYCTwDXB5tdT6zTXL9y3OCcuOvHFqlJqogMDJ09KVyQoHOOAh4LhsvIAh5092fM7HXgETP7NLAN+EiCzn/MivKz2V/bFBsuNpCfnclt8yYnLSYRkd7UYVJw932JOGEwwN7pcdbvJXGJqMdKt+3nnT21XDFjLE+t2UlTJEpJUT63zZusSmYRGTDCDHMhwP/8dTND87L49hXTqKhuAODhm89JclQiIr1LSSGE7fvq+NPaSm4+70QKcrOUDERkwAoz9lHau//vW8gw4/pzJiU7FBGRhFJS6EJtYwsPL9/OxdPGMLpQw2OLyMCmpNCFxSt3cLChhRvOnZTsUEREEk5JoRPRqPOrv2/h9HGFzJxQlOxwREQSTkmhE38t28M7VbXcOPd4WqchFREZyJQUOvGrVzZTPCSXS6aNSXYoIiJ9QkmhA5v31PLSxio+cfZEcrL0zyQi6UF3uw789h9byc40rj17fLJDERHpM0oKcdQ3Rfj98u1cdNoYRg5RM1QRSR9KCnH88Y0KahpauG7OxGSHIiLSp5QUjuDuPPCPLUweNYSzJh2X7HBERPqUksIR3thRzdryGj5xzkQ1QxWRtKOkcITfvLqVwTmZXKHhsEUkDSkptLO/toknV1dwxcwSCnI1gKyIpB8lhXYeXbmDxpYoHz9bFcwikp6UFALuzqLXt3PGhCLeM2ZossMREUkKJYXAiq37Kdt9iGvPmpDsUEREkkZJIfDQsu0U5GZx6eka50hE0peSAlBd38xTayr40IyxDMpRBbOIpC8lBeDxVeU0NEdVdCQiaS/tk4K789Cy7Zw6dijTxhUmOxwRkaRK+6Swprya9ZU1XDNbTwkiImmfFB5Zvp3crAwumzE22aGIiCRdWieFhuYIT6yq4OLTRjM0LzvZ4YiIJF1aJ4Vn1+2ipqGFq87URDoiIpDmSeH3K3ZQUpTPuScOT3YoIiL9Qtomhcrqev66qYorZ5aQkaEhskVEAPpdTy0zuwj4MZAJ/I+7L+jN4y8pLefOpRspP1APwNB81SWIiLTqV08KZpYJ/BdwMTAVuNbMpvbW8ZeUlnP74jVtCQHgB39+iyWl5b11ChGRlNavkgIwGyhz93fcvQlYBFzWWwe/c+lG6psjh62rb45w59KNvXUKEZGU1t+SQgmwvd3yjmBdGzO7ycyWm9nyqqqqbh28ot0TQpj1IiLppr8lhXg1vn7YgvtCd5/l7rOKi4u7dfCxRfndWi8ikm76W1LYAbTvNDAOqOitg982bzL52ZmHrcvPzuS2eZN76xQiIimtv7U+eh042cyOB8qBa4CP9dbBLz8jVhJ159KNVByoZ2xRPrfNm9y2XkQk3fWrpODuLWZ2K7CUWJPU+9z9zd48x+VnlCgJiIh0oF8lBQB3fxp4OtlxiIiko/5WpyAiIkmkpCAiIm2UFEREpI2SgoiItDF373qrfsrMqoCtx7j7CGBPL4aTTLqW/mmgXMtAuQ7QtbSa6O5xe/+mdFLoCTNb7u6zkh1Hb9C19E8D5VoGynWAriUMFR+JiEgbJQUREWmTzklhYbID6EW6lv5poFzLQLkO0LV0KW3rFERE5Gjp/KQgIiJHUFIQEZE2Az4pmNlFZrbRzMrMbH6cz3PN7OHg89fMbFLfRxlOiGu5wcyqzGxV8PpMMuLsipndZ2a7zWxtB5+bmf0kuM7VZjazr2MMK8S1nG9m1e2+k6/3dYxhmNl4M3vRzNab2Ztm9sU426TE9xLyWlLle8kzs2Vm9kZwLd+Ks03v3sPcfcC+iA2//TZwApADvAFMPWKbfwH+O3h/DfBwsuPuwbXcAPws2bGGuJb3AzOBtR18fgnwJ2Iz8c0BXkt2zD24lvOBJ5MdZ4jrGAPMDN4PAd6K899XSnwvIa8lVb4XAwqC99nAa8CcI7bp1XvYQH9SmA2Uufs77t4ELAIuO2Kby4D7g/d/AC4ws3jTgiZbmGtJCe7+MrCvk00uAx7wmH8ARWY2pm+i654Q15IS3L3S3VcG7w8C6zlifnRS5HsJeS0pIfi3PhQsZgevI1sH9eo9bKAnhRJge7vlHRz9H0fbNu7eAlQDw/skuu4Jcy0AVwaP9n8ws/FxPk8FYa81VZwTPP7/ycxOTXYwXQmKH84g9qu0vZT7Xjq5FkiR78XMMs1sFbAbeNbdO/xeeuMeNtCTQrxseWSWDbNNfxAmzj8Ck9x9OvAc7/56SDWp8p2EsZLYODOnAz8FliQ5nk6ZWQHwKPAld6858uM4u/Tb76WLa0mZ78XdI+4+g9ic9bPN7LQjNunV72WgJ4UdQPtfy+OAio62MbMsoJD+WRzQ5bW4+153bwwWfwGc2Uex9bYw31tKcPea1sd/j80qmG1mI5IcVlxmlk3sJvo7d18cZ5OU+V66upZU+l5aufsB4CXgoiM+6tV72EBPCq8DJ5vZ8WaWQ6wS5okjtnkCuD54fxXwggc1Nv1Ml9dyRPnuh4iVpaaiJ4BPBq1d5gDV7l6Z7KCOhZmNbi3fNbPZxP6f25vcqI4WxPhLYL27393BZinxvYS5lhT6XorNrCh4nw9cCGw4YrNevYf1uzmae5O7t5jZrcBSYq137nP3N83sDmC5uz9B7D+e35hZGbHsek3yIu5YyGv5gpl9CGghdi03JC3gTpjZQ8Raf4wwsx3AN4hVoOHu/01sju5LgDKgDrgxOZF2LcS1XAX8HzNrAeqBa/rpj465wHXAmqD8GuCrwARIue8lzLWkyvcyBrjfzDKJJa5H3P3JRN7DNMyFiIi0GejFRyIi0g1KCiIi0kZJQURE2igpiIhIGyUFERFpM6CbpIr0FjMbDjwfLI4GIkBVsFzn7ucmJTCRXqYmqSLdZGbfBA65+13JjkWkt6n4SKSHzOxQ8Pd8M/uLmT1iZm+Z2QIz+3gwHv4aMzsx2K7YzB41s9eD19zkXoHIu5QURHrX6cAXgWnEetWe4u6zgf8BPh9s82Pgh+5+FnBl8JlIv6A6BZHe9XrreEBm9jbw52D9GuCfgvcXAlPbDXk/1MyGBGP/iySVkoJI72ps9z7abjnKu/+/ZQDnuHt9XwYmEoaKj0T63p+BW1sXzGxGEmMROYySgkjf+wIwK5ghbx3wuWQHJNJKTVJFRKSNnhRERKSNkoKIiLRRUhARkTZKCiIi0kZJQURE2igpiIhIGyUFERFp8/8DXkqmMIJNswgAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret2 = ensemble_simulations(\n", " duration, ndiv=20, y0=y0, volume=V, model=m, solver=('egfrd', Integer3(8, 8, 8)), repeat=N)\n", "ret2.plot('o', ret1, '-')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Simulating with `bd`:" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEGCAYAAACKB4k+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXxU9b3/8dcnKwkJBAiyLy4Bq6KAiAtdtNpi7f3VrXXpptZWe6tXu1xabXtvvb212tKr125e9dZWa6twFVHrQhW3LloEgixiBJUlCxCWhED2yef3x0xCIJPkkMxkMpP38/GYR2a+c86Zz3HwfOZ8V3N3REREANISHYCIiPQfSgoiItJGSUFERNooKYiISBslBRERaZOR6AB6o7Cw0CdPnpzoMEREksqKFSt2uvvIaO8ldVKYPHkyy5cvT3QYIiJJxcw2d/aeqo9ERKSNkoKIiLSJW1Iwswlm9pKZrTezdWZ2Y6T8FjMrM7NVkcd57fa52cw2mlmJmc2NV2wiIhJdPNsUmoFvuftKM8sHVpjZ85H37nT3n7Xf2MyOAy4DjgfGAi+Y2RR3D8UxRhERaSdudwruXuHuKyPPa4D1wLgudjkfeMTdG9z9fWAjMDte8YmISEd90qZgZpOBGcA/IkXXm9lqM7vfzIZFysYBW9vtVkqUJGJm15jZcjNbXllZGceoRUQGnrgnBTPLAx4Dvu7ue4G7gaOB6UAF8F+tm0bZvcMUru5+r7vPcvdZI0dG7WYrIiI9FNekYGaZhBPCH9x9EYC7b3f3kLu3APdxoIqoFJjQbvfxQHk84xMRkYPFs/eRAb8B1rv7He3Kx7Tb7EJgbeT5k8BlZpZtZkcCRcCyeMUnIpJIl97zGpfe81qiw+ggnr2P5gBfANaY2apI2XeBy81sOuGqoU3AtQDuvs7MFgJvEe65dJ16HolIPLVelBdce3qf7ttb8fzsuCUFd/8r0dsJnulin1uBW+MVk4iIdE0jmkWk1/prVYgcPiUFEQF0YZcwJQURSWpKZrGlpCAiCTfQLuyLi8so3lLFP97fzZzbX2RxcVmn27o7TaEWauqb2Lmvgfv/+h4rt+wJtG9PJPV6CiIiidB6UW8MtXDG7Uu54aNFfPTYI9jfGGJ/QzO1jSFqG1v/hqhr93xNaRWvbNhJqCU8Nresqo5vLlzFnc+XMDg7k/qmUPjR3NL2vKXDMF7a9r150RoALpjR1SxCwSkpiKSQRHaTTDaPryileMseGkPO7Ftf4JJZ4zlpwjBq6puoqW9mX0Nz5G8T+9q9Lt1TS3lVfdt0C+VV9dwUuTD3VIvDtr0NfKgon+zMdHIy0xmUmcagjHQGZaaTnZHGoMx0fvnSRqrrmg7at64pxPwlJUoKIiIQrl5pbnE2bK9hT20Te2obqaptZE9tE1W1TVTXNVFd19jueRM7a+qpbz7w83tHTQO/fOndDsfOykgjPzuDwdkZ5GVnkDcog537GjvOvwMMzcnklk8dR25WBrlZ6eRmZTA4O53czAxystLJzQpf7I/+7jNR929sbuF/rzily3P98TPro5aXV9V1ud/hUFIQkX6lpcWprmuicl8DO2sa2Lm/kZ01Deza38CufY3s2t/I7naP1l/OH7vz1Q7Hykw3huZkUZCbydCcTEYNGcSUUfn8ed02oOPY2JH52Sy89nSGDAongOyM9A7bHHnT01Hj3lvXxIUzxnd7fmMLciiLchEfW5AT132DUlIQkT7R0uLsrm1kW3U9O2rq2VbdwI6aerbvbaBkWw1NoRZOv20plTUNNEepRE9PM4blZlGYl8XwwVkcP3YIwwdnsXT9djLS0/jWx6cyLDeTYbnhJDAsN4vcrHTCM+4crLML+86aBo4sHNzlefT2wjxv7lRuXrSGuqYDSSknM515c6fGdd+glBRE+plkbRdobmlhfcVeyqvqKKuqo7yqnorqOiqq6qnYW8f26gYaQy0d9ivMy6Ix1EJmehofPKaQwvxsRuZlMzI/m8K8bEbmZzFicDZDczJJSzv4Ar+4uIxHlm2lMdTCT559m3lzp/Khou5nT+7Nhb23F+bWuv9vP7qaxlAL4wpymDd3aqA2gd7sG5SSgogc1Jtmzu0vRr3Q1DeFKN1Ty+ZdtWzdXcuW3XVs3VNL6Z463tleQ6jF+cRdf2nbPjPdGD10EGOG5DBz4jBGDx3E6CGDGDN0EKOGhB8j87N5enUF3350NbWNTfz93V2BL3KLi8u4edGatkRzOD1xenNhj8WF+YIZ43h42Rbg8JN/b/YNQklBZICLdnH99qNvsnT9dvIGZfD+zv1s3lXLtr31eLtanZzMdCYMzyEjzfDIG8NyM7nijMl8dvZECvOyO/yyD/LZQS/s85eUHHRRh+A9cXp7YY/3hTmRlBREBqDquiY27qhhw/Z9/OjptzpcXBtDzlOrKxgxOIvJhYM5/agRTBoxmEkjcpkwPJeJw3MpzMviiVXl3LxoTVs/+j21TdzzyntMHjE40AW2Nxf2znrcBO2Jk8oX9t5QUhBJEdGqgM49YTQbtu9j/ba9lGyr4Z3t4cf2vQ3dHs+AFf/2sS636c1FHXp3Ye+LnjgDkZKCSApYXFzGTYtWH1QN840Fq/j6ggPbZGekUTQqjznHFFJ0RD5TRuVRdEQ+l933GuVV9R2OGeTi2ttf64ls8E20/np3oqQgEgfx7kFUVdvI6tJq1pRV8+bWKpau307okF6cDuRnZ/CTT5/IsaPzmTRiMOlR6vi/PffYHl9cE9k9sy964vRX8UwoSgoi/Ui0KqBPnjiGtytqWLllD6u2VlG8ZQ+bdtW27XNU4eAOCaHVvoZmzps2JvqbEb25uCaye2br/moXiC0lBZF+IlpPnG8uXMW//t8qmiPd+4/Iz2b6hAIuOWUC08cXcML4oQwZlMmc21/s1S/2nl5cE909s7d683mpmoSUFEQSrLquiWXv7+Z7i9d0aLRtccjJSueOi07k5EnDGDt0UNQRuomsX9ev9dSipCDSx+qbQix7fzd/27iTv7+7i3Xl1Z1OjQxQ2xDiUyeN7fKYA7l+XWJLSUEkztydku01vFJSyasbKnlj0x4am1vITDdmTBzGDWcXcfpRI/jGglWUV/esFxDoF7vEhpKCSByEIjN93vTYal4uqWTb3vDFfsqoPL5w2iQ+WFTIqUcOJzfrwP+C3z63572ARGJFSUEkRrbvref5t7bzwvrtrNi8Bwe2VdfzwaJCzpp6BB+aUsiYoZ3/6lcVUM/orii2lBREogg6zmDzrv08t3Ybz63bRvGWKgCGD87EDNwhf1AGc48frS6WkjSUFEQO0+Zd+3l6TQVPr65gXfleAE4YN4RvfWwKGenGz5duaGs4Lq+uj/kauqlICbD/UFIQCWDH3nr+tLqCJ94s582t4TuCGRML+P4nP8Dc40czYXguAHNuf5G6poPXDIj1Groi8aSkINKJUIuzuLiMx1aW8reNO2lxOH7sEG7+xLF88sQxjB+W22Gf3s4FlEj6tS6gpCByEHdn5ZY9vFe5n137G1i+eQ/jh+Vw3VnHcP70sRxzRH6X+w/UmTuVUFKHkoIIsHt/I4+tKOWRN7bwbuV+0gxGDM7il5+dySmTh3e7WEyrRM/cqYuz9JaSggxY7s7yzXv4/WubeW7tNhpDLcycWMBPLz6Rhcu3kp5mnHrUiMM6prqVSrJTUpABp64xxOJVZTz42mbWV+wlf1AGl8+ewOWnTuTY0UMAeGxlaY+P39tupfq1L4kUt6RgZhOAB4HRQAtwr7vfZWbDgQXAZGATcIm777HwLF93AecBtcCV7r4yXvHJwLOtup4HXtvEw8u2UFXbxLGj87ntommcP33sQSOLgyxiL5Kq4nmn0Ax8y91Xmlk+sMLMngeuBJa6++1mdhNwE/Ad4BNAUeRxKnB35K9Ij7QOQLvlU8dz76vv8dSb5bS48/HjRnPVnMnMPnJ4hxlHe7OQvEgqiFtScPcKoCLyvMbM1gPjgPOBMyObPQC8TDgpnA886O4OvG5mBWY2JnIckcPi7uyta6K8up5P3PUXcrPS+cLpk/jSnCPbxhRE09s1h0WSXZ+0KZjZZGAG8A9gVOuF3t0rzOyIyGbjgK3tdiuNlCkpSGDuzssllfzypY2s31ZDRpoxb+5UPn/qJIbmZna7fzKPMxCJhW6Tgpl9Bngu8mv/+8BM4EdB6/vNLA94DPi6u++NtkBI66ZRyjrMMm9m1wDXAEycODFICDIAuDsvvr2D/35hA2vKqhlXkMOkEbkckZfNdWcdE/g4A3WcgUirtADb/FskIXwQmEu4yufuIAc3s0zCCeEP7r4oUrzdzMZE3h8D7IiUlwIT2u0+Hig/9Jjufq+7z3L3WSNHjgwShqSw8J3BDi741d+4+oHlVNU18pOLp/HSv57J6CGDAo8vaDVv7lRyMtMPKtP01TKQBKk+aq1g/SRwt7s/YWa3dLdTpDfRb4D17n5Hu7eeBK4Abo/8faJd+fVm9gjhBuZqtSdIV5Zv2s1Pnyth2abdjB+Ww08unsZFM8eTmR7kt050sRpnoG6lkqyCJIUyM7sHOAf4iZllE+wOYw7wBWCNma2KlH2XcDJYaGZXA1uAz0Tee4Zwd9SNhLukXhX4LCRlRZvCesP2Gm5/9m2Wvr2DkfnZ/Of5x3PpKRPJyjjwz7I33Uo1fbUMZEGSwiXAucDP3L0qUuUzr7ud3P2vRG8nADg7yvYOXBcgHhmgduyt584X3mHBG1sZnJXBt8+dypVnTD5ojAGoW6lIb3SbFNy91sx2AB8ENhAef7Ah3oGJtGppcX710kZ+/dJGGkMtXHHGZP7lo0UMH5wVdXt1KxXpuSC9j34AzAKmAr8FMoGHCFcPicSNu7N7fyNbdtfyxuY9fOy4UXzvvA8wuXBwl/upW6lIzwWpPrqQ8BiDlQDuXh4ZoSwSN5t27ueWp9axYcc+cjLTeejqU/lgUWGgfdWtVKTngjQYN0bq+x3AzLr+mSZyiEvvea2twbg79U0h7nj+HT7+36+yfNMeJg7P5YRxQwInBFC3UpHeCJIUFkZ6HxWY2VeAF4D74huWDESvvbuL8+76Cz9fuoFPnDCab8+dwrbqet7YtIc5t7/I4uKyQMe5YMY4brtoGlmRrqnjCnK47aJpak8QCSBIQ/PPzOxjwF7C7Qr/7u7Pxz0yGTCq65q49em3WLi8lInDc/n91bPZta+xVz2INH21SM8EmvsokgSUCCTmXnhrO999fA279jfy1Y8czY1nF5GTlc6c219UDyKRBOg0KZhZDeF2BOPgOYiM8LCCIXGOTVJYdW0Ttzy1jseLyzh2dD6/ueIUpo0f2va+ehCJJEanScHd1cNI4uIvGyqZ93+r2bmvgRvOLuL6s445aDQyqAeRSKIEGacQdSpSd98S+3AkldU3hbjtmfU88Npmjh45mHu/eAYnji+Iuu28uVO5edGag6qQ1INIJP6CtCk83e75IOBIoAQ4Pi4RSUppPwfRCT9YQnOLc9WcyXzn3GMZdEi30fZiNTGdiByeIL2PprV/bWYzgWvjFpGkjPAcRKvbehA1tzhZ6WmcNL6gy4TQSj2IRPreYc8xHFlc55Q4xCIp5ifPvU1dU8tBZY2hFuYvKUlQRCLSnSBtCt9s9zKN8MprlXGLSPqlaFNYd2VtWTUV1fVR31MPIpH+K0ibQvteSM2E2xgei084kgoWLt/K9xevJc2gpcOCqupBJNKfBWlT+I++CESSX0NziB8+9RZ/+McWPnhMIeeeMIpbn35bPYhEkkiQ6qPngc+4e1Xk9TDgEXefG+/gJHns2FvPtQ+toHhLFV/9yNHMmzuV9DQjLztTPYhEkkiQ6qORrQkBwN33mNkRcYxJksya0mq+8uBy9tY38evPzeS8aWPa3lMPIpHkEqT3Uaj9ADYzm8TB017IAPbUm+V85p6/k55mPPbPZxyUEEQk+QS5U/ge8FczeyXy+sPANfELSZKBu/OLFzdyx/PvcMrkYdz9+ZMpzMtOdFgi0ktBGpqfiwxYO43wZHjfcPedcY9M+q3G5hZuXrSGx1aWctGMcdx28TSyM7ofjCYi/V+QhmYDzgWOcvcfmtlEM5vt7sviH570B+2nqjj9tqUMzkpnY+V+vn5OETeeXUT4n4iIpIIg1Ue/BlqAjwI/BGoIj1PQqOYBIDxVxYHFbloHpH129gS+fs6UQMdQY7FI8giSFE5195lmVgxtvY+y4hyX9BPzl5R0WOwG4JV3VIMokoqC9D5qMrN0Ij2OzGwk4TsHGQC02I3IwBIkKfwceBw4wsxuBf4K/DiuUUm/MSIv+k2hpqoQSU1Beh/9wcxWAGcT7n10gbuvj3tkknAvl+ygqraxw3qsmqpCJHV1eqdgZsNbH8AO4GHgj8D2SJmksCXrtvGVB5czZdQQfnjB8WSlh/+pjCvI4baLpmmqCpEU1dWdwgrCPxCj9Td04Ki4RCQJ9+yaCv7l4WJOGDeUB740m6E5mfzpzQpAPYlEUl2nScHdj+zLQKR/eHp1BTc8Usz0CQX87qpTyB+UmeiQRKQPBemSipl9ivD0FgAvu/uf4heSJEprQpg5sYDfXjWbvOxA/zxEJIV02/vIzG4HbgTeijxuNLPbAux3v5ntMLO17cpuMbMyM1sVeZzX7r2bzWyjmZWYmabljrFL73mtbfW0aJ5/azs3PlLMjAkF/E4JQWTACvJ//nnAdHdvATCzB4Bi4OZu9vsd8EvgwUPK73T3n7UvMLPjgMuA44GxwAtmNsXdO46akph75Z1KrvvDSo4fN5TfXnUKg5UQRAasIOMUAAraPR8aZAd3fxXYHfD45xNeuKfB3d8HNgKzA+4rvfD6e7u45sHlHHNEHg9eNVttCCIDXJCfhLcBxWb2EuGeSB+m+7uErlxvZl8ElgPfcvc9wDjg9XbblEbKJI7WlFbz5QeWM2F4Lg99+VSG5naeENTrSGRg6PZOwd0fJjxt9qLI43R3f6SHn3c3cDQwHagA/itS3lm31w7M7BozW25myysrK3sYhrxbuY8rfruMoTmZPHT1qQwfrOmsRKSLO4XIGgrtlUb+jjWzse6+8nA/zN23tzv+fUBrL6ZSYEK7TccD5Z0c417gXoBZs2ZpBbgeqKiu44u/WYYBv796NqOHDkp0SCLST3RVfbQcWAe0/hxv/2veCU+lfVjMbIy7V0ReXgi09kx6Evijmd1BuKG5CNB6DXGwt76JK+9/g+q6Jh655jSOGpmX6JBEpB/pKil8C7gYqAMeAR53931BD2xmDwNnAoVmVgr8ADjTzKYTTiqbgGsB3H2dmS0k3OW1GbhOPY9ip/0iObN+9AJNzS089OVTOWFcoD4DIjKAdDWi+U7gTjM7ErgcWGpmm4Efu/uq7g7s7pdHKf5NF9vfCtzafchyOA5dJKexuYXMdKOypiHBkYlIfxSkofl94Angz4S7iQZbbkv6hWiL5DSFnPlLShIUkYj0Z101NB9FeEDZ+cBWwlVIt7p7fR/FJjGgRXJE5HB01aawEVhN+C5hLzAR+FrrIu3ufkfco5NeG5GXxc59jR3KtUiOiETTVVL4IQfGCqiLShLavGs/tQ3NWiRHRALrqqH5lj6MQ2JsX0MzX35gOVmZ6XzjY1OZv6SExlAL4wpymDd3qhbJEZGoNPNZCmppcb65YBXv7dzPg1+azZxjCnlhfXjcoKarEJGuBJ0QT5LIL17cyJ/f2s53z/sAc44pTHQ4IpJEulqj+cbI3zl9F4701vNvbefOF97hopnj+NKcyYkOR0SSTFd3CldF/v6iLwKR3tu0cz/fXLCKaeOG8uMLp9HaU0xEJKiu2hTWm9kmYKSZrW5XboC7+4lxjUwOS11jiK8+tIL0dOPXn5vJoMz0RIckIkmoq95Hl5vZaGAJ8Km+C0kOl7vzb0+spWR7Db+98hQmDM9NdEgikqS67H3k7tuAk8wsiwPTW5S4e1PcI5PAFi7fyqMrSrnh7CLOnHpE1G3U60hEgui2S6qZfYTwOsubCFcdTTCzKyLLbUqClWyr4d+fWMeHigq58eyiRIcjIkkuyDiFO4CPu3sJgJlNAR4GTo5nYNK92sZmrv/jSobkZHLHJdNJT1PDsoj0TpBxCpmtCQHA3d8BtLp7P3DLk+vYWLmP/750OiPzsxMdjoikgCB3CsvN7DfA7yOvPwesiF9IEsQTq8pYuLyU6886RgPURCRmgiSFfwauA24g3KbwKvDreAYl0V16z2sA/OwzJ/H9x9cya9Iwvn6O2hFEJHa6TQru3kC4XUFTZfcD7s63Fr6JA3deOp2MdM1UIiKxownxkkx5dT2le+q445KTNB5BRGJOSSFJLC4uY+WWPTSFnJzMNNTPSETiocu6BzNLN7P5fRWMRLe4uIybFq2mKRReKqeuqYXvPr6WxcVlCY5MRFJNl0nB3UPAyaaZ1RJq/pIS6ptaDiqrawoxf0lJJ3uIiPRMkOqjYuAJM/s/YH9robsviltUcpDyqrrDKhcR6akgSWE4sAv4aLsyB5QU+kBtYzNpaUaoxTu8N7YgJwERiUgqC9Il9arutpH4+elzJYRanKyMNBqbD1Qh5WSmM2/u1ARGJiKpqNtO7mY2xcyWmtnayOsTzez78Q9NVmzezQOvbeKLp0/ipxefSFZkTMK4ghxuu2gaF8wYl9gARSTlBKk+ug+YB9wD4O6rzeyPwI/iGdhA19Ac4juPrWHs0By+fe6x5GVn8PCyLYCmwRaR+AkyHDbX3ZcdUtYcj2DkgF+9uJGNO/Zx64UnkJet4SQi0jeCJIWdZnY04cZlzOzTQEVcoxrg3t62l1+//C4XzhjX6aI5IiLxEOQn6HXAvcCxZlYGvE94plSJg1CLc9NjaxiSk8m//dNxB72naiMRibcgvY/eA84xs8FAmrvXxD+sgevhZVtYtbWKOy89ieGDsxIdjogMMEF6H40ws58DfwFeNrO7zGxE/EMbeCprGvjJc29zxtEjuGC6ehaJSN8L0qbwCFAJXAx8OvJ8QXc7mdn9ZrajtStrpGy4mT1vZhsif4dFys3Mfm5mG81stZnN7NnpJLcfP7OehqYW/vOCE9DMIiKSCEGSwnB3/093fz/y+BFQEGC/3wHnHlJ2E7DU3YuApZHXAJ8AiiKPa4C7gwSfSv6+cSePF5fx1Y8cxdEj8xIdjogMUEGSwktmdpmZpUUelwBPd7eTu78K7D6k+HzggcjzB4AL2pU/6GGvAwVmNibYKSS/xuYWvv/EWiYOz+VrZx2T6HBEZADrtKHZzGoId0M14JvAQ5G30oB9wA968Hmj3L0CwN0rzKy1v+U4YGu77UojZR26vprZNYTvJpg4cWIPQuh/fvu393mvcj+/vfIUBmWmJzocERnAOr1TcPd8dx8S+Zvm7hmRR5q7D4lxHNEq0DvOABeO6153n+Xus0aOHBnjMPre9r31/HzpBs4+9gjOOlZjEkQksQINlTWzE4HJ7bfv4dTZ281sTOQuYQywI1JeCkxot914oLwHx086P3n2bZpC3mFMgohIInSbFMzsfuBEYB3QOk1nT6fOfhK4Arg98veJduXXm9kjwKlAdWs1UypbsXk3i4rLuO6so5lcODjR4YiIBLpTOM3dD/tnrJk9DJwJFJpZKeE2iNuBhWZ2NbAF+Exk82eA84CNQC2QstN1X3rPawA8/JXT+MGT6xgzdBDXqXFZRPqJIEnhNTM7zt3fOpwDu/vlnbx1dpRtnfB0GgPGYytLWVu2l7sum05ulia8E5H+IcjV6AHCiWEb0EC4Udjd/cS4RpbCQi3O/CUlTJ9QwKdOGpvocERE2gRJCvcDXwDWcKBNQXqhorqOHTUN3P35kzVyWUT6lSBJYYu7Pxn3SAaAxcVlrNyyh6aQk5OZxtbdtZw8aViiwxIRaRMkKbwdWWntKcLVR0CPu6QOWIuLy7h50RqaQuHhF3VNLdy8aA2AltUUkX4jyDQXOYSTwceB/xd5/FM8g0pF85eUUNcUOqisrinE/CUlCYpIRKSjIOsppGz30L5UXlV3WOUiIokQZPDab4ky5YS7fykuEaWoEXlZ7NzX2KF8bEFOAqIREYkuSJvCn9o9HwRcyACZgiJW3J2czPRwX9525TmZ6cybOzVRYYmIdBCk+uix9q8jI5VfiFtEKejZtdvYuqeOy2dP4LEVZTSGWhhXkMO8uVPVyCwi/UpPhtIWAakxZ3UfaA61MH9JCVNH5fOjC6bxXuV+ABZce3qCIxMR6ShIm0LrugqttgHfiVtEKWbRyjLe37mf+744i/Q0DVQTkf4tSPVRfl8EkooamkPctXQDJ00o4JwPaK0EEen/ulp5rcsqInffEvtwUsvCN7ZSVlXHbRdNa5vOQtVGItKfdXWn8DQHluNs5cBI4AhA60Z2ob4pxC9e3MjsycP5UFFhosMREQmk06Tg7tPavzazyYTbEs4BfhzXqFLAQ69vZkdNA7+4fIYmvRORpNHtNBdmVmRmvwOeBVYAx7n7L+IdWDKrbWzm7pff5UNFhZx61IhEhyMiElhXbQonAN8Djgd+Clzt7qHOtpcD/vD6Fnbtb+Tr5xQlOhQRkcPSVZvCm8BWwm0Ls4HZ7atB3P2G+IaWnOoaQ9zz6nt88JhCTp40PNHhiIgclq6SguY26oGHl21h574Gbjh7ZqJDERE5bF01ND/Ql4GkgvqmEP/zyrucdtRwZh+puwQRST5B1lOQgBa8sZUdNQ3ccLbaEkQkOSkpxEhDc/gu4ZTJwzhdPY5EJEkpKcTI4uIyKqrruf6jRRqXICJJK8g4hSlmttTM1kZen2hm349/aMkj1OLc88p7HD92CB/W6GURSWJB7hTuA24GmgDcfTVwWTyDSjZL1m3jvZ37+dqZx+guQUSSWpCkkOvuyw4pa45HMMnI3bn75Xc5snAw554wOtHhiIj0SpCksNPMjiaypoKZfRqoiGtUSeSvG3eypqyaaz98lNZLEJGkF2TlteuAe4FjzawMeB/4fFyjSiJ3v/wuo4Zkc+FMLaspIskvyCI77wHnmNlgIM3da+IfVnJYXVrF39/dxXfPO5bsDM0kLiLJL8hynNnAxcBkIKO1IdXdfxjXyJLAfX95n/zsDC6frSWrRSQ1BNEH0e0AAArxSURBVKk+egKoJjxtdkN8w0kOl97zGg3NIdaU7eVLcyaTPygz0SGJiMREkKQw3t3PjeWHmtkmoAYIAc3uPsvMhgMLCN+RbAIucfc9sfzcWNpWHc6PV845MsGRiIjETpDeR383s2ndb3bYznL36e4+K/L6JmCpuxcBSyOv+53FxWUUb9nDtr31ZKWn8cb7uxMdkohIzHS1yM5aoCWyzVVm9h7h6iMD3N1PjHEs5wNnRp4/ALxMePnPfmNxcRk3L1pDY8gBqGsKcfOiNQBcMEO9j0Qk+XVVfTQOmB6nz3Xgz2bmwD3ufi8wyt0rANy9wsyOiLajmV0DXAMwcWLfNvDOX1JCXdPBi8/VNYWYv6RESUFEUkJXSeF9d98cp8+d4+7lkQv/82b2dtAdIwnkXoBZs2Z5nOKLqryq7rDKRUSSTVdJ4Qgz+2Znb7r7HT39UHcvj/zdYWaPE17uc7uZjYncJYwBdvT0+PEytiCHsigJYGxBTgKiERGJva4amtOBPCC/k0ePmNlgM8tvfQ58HFgLPAlcEdnsCsJdYfuVS08Z36EsJzOdeXOnJiAaEZHY6+pOoSJOA9RGAY9HBsFlAH909+fM7A1goZldDWwBPhOHz+6VTbtqyUoPD95rDDnjCnKYN3eq2hNEJGV0lRTiMrtbZNqMk6KU7wLOjsdnxkJlTQN/erOCy2dP5O1t4Zk+Flx7eoKjEhGJra6SQr+9QCfCw8u20Bhq4YtnTObokXmJDkdEJC46bVNwd43KimhsbuGh1zfzkSkjlRBEJKVpjeYAlqzbxo6aBq48Y3KiQxERiSslhQAeen0zE4bn8JEpIxMdiohIXCkpdOOd7TX84/3dfHb2JNK0spqIpDglhW784fXNZKWnccmsjmMURERSjZJCF/Y3NLNoZRnnTRvNiLzsRIcjIhJ3SgpdeGJVOTUNzXz+tEmJDkVEpE8oKXTC3Xno9c0cOzqfkycNS3Q4IiJ9QkmhE8Vbq3irYi+fP20SretSi4ikOiWFTixYtpXcrHTNayQiA4qSQhT7Gpp5anU5/3TiGPKygyxjLSKSGpQUonhyVTm1jSEum923K7uJiCSakkIUC97YwpRRecyYUJDoUERE+pSSwiHeKt/Lm6XVXHbKRDUwi8iAo6RwiEfe2EJWRhoXzVQDs4gMPEoK7dQ3hXi8uIxzjx9NQW5WosMREelzSgrtPLd2GzX1zVx2yoREhyIikhBKCu08uqKU8cNyOO2oEYkORUQkIZQUIsqq6vjbuzu5eOZ4TZEtIgOWkkLEohWluMOnT9YU2SIycCkpEJ787tGVpZx21HAmDM9NdDgiIgmjpAAs37yHzbtq+fTJamAWkYFNSQF4dHkpg7PSOW/a6ESHIiKSUAM+KdQ2NvP0mgrOmzaG3CxNficiA9uATwp/XredfQ3NXKwGZhERJYXHi8sYV5DD7MnDEx2KiEjCDeikUFnTwF82VHLBjLEamyAiwgBPCk+9WU6LwwXTNfmdiAgM8KSweFUZJ4wbQtGo/ESHIiLSLwzYpLBxxz5Wl1brLkFEpJ1+1wfTzM4F7gLSgf9199tjefzFxWXMX1JCWVUdAFnpaksQEWnVr+4UzCwd+BXwCeA44HIzOy5Wx19cXMbNi9a0JQSA254tYXFxWaw+QkQkqfWrpADMBja6+3vu3gg8Apwfq4PPX1JCXVPooLK6phDzl5TE6iNERJJaf0sK44Ct7V6XRsramNk1ZrbczJZXVlYe1sHL290hBCkXERlo+ltSiFbB7we9cL/X3We5+6yRI0ce1sHHFuQcVrmIyEDT35JCKdB+qtLxQHmsDj5v7lRyMtMPKsvJTGfe3Kmx+ggRkaTW33ofvQEUmdmRQBlwGfDZWB38ghnhmqj5S0oor6pjbEEO8+ZObSsXERno+lVScPdmM7seWEK4S+r97r4ulp9xwYxxSgIiIp3oV0kBwN2fAZ5JdBwiIgNRf2tTEBGRBFJSEBGRNkoKIiLSRklBRETamLt3v1U/ZWaVwOYe7l4I7IxhOImkc+mfUuVcUuU8QOfSapK7Rx39m9RJoTfMbLm7z0p0HLGgc+mfUuVcUuU8QOcShKqPRESkjZKCiIi0GchJ4d5EBxBDOpf+KVXOJVXOA3Qu3RqwbQoiItLRQL5TEBGRQygpiIhIm5RPCmZ2rpmVmNlGM7spyvvZZrYg8v4/zGxy30cZTIBzudLMKs1sVeTx5UTE2R0zu9/MdpjZ2k7eNzP7eeQ8V5vZzL6OMagA53KmmVW3+07+va9jDMLMJpjZS2a23szWmdmNUbZJiu8l4Lkky/cyyMyWmdmbkXP5jyjbxPYa5u4p+yA8/fa7wFFAFvAmcNwh23wN+J/I88uABYmOuxfnciXwy0THGuBcPgzMBNZ28v55wLOEV+I7DfhHomPuxbmcCfwp0XEGOI8xwMzI83zgnSj/vpLiewl4LsnyvRiQF3meCfwDOO2QbWJ6DUv1O4XZwEZ3f8/dG4FHgPMP2eZ84IHI80eBs80s2rKgiRbkXJKCu78K7O5ik/OBBz3sdaDAzMb0TXSHJ8C5JAV3r3D3lZHnNcB6DlkfnST5XgKeS1KI/LfeF3mZGXkc2jsoptewVE8K44Ct7V6X0vEfR9s27t4MVAMj+iS6wxPkXAAujtzaP2pmE6K8nwyCnmuyOD1y+/+smR2f6GC6E6l+mEH4V2l7Sfe9dHEukCTfi5mlm9kqYAfwvLt3+r3E4hqW6kkhWrY8NMsG2aY/CBLnU8Bkdz8ReIEDvx6STbJ8J0GsJDzPzEnAL4DFCY6nS2aWBzwGfN3d9x76dpRd+u330s25JM334u4hd59OeM362WZ2wiGbxPR7SfWkUAq0/7U8HijvbBszywCG0j+rA7o9F3ff5e4NkZf3ASf3UWyxFuR7Swruvrf19t/DqwpmmllhgsOKyswyCV9E/+Dui6JskjTfS3fnkkzfSyt3rwJeBs495K2YXsNSPSm8ARSZ2ZFmlkW4EebJQ7Z5Ergi8vzTwIseabHpZ7o9l0Pqdz9FuC41GT0JfDHS2+U0oNrdKxIdVE+Y2ejW+l0zm034/7ldiY2qo0iMvwHWu/sdnWyWFN9LkHNJou9lpJkVRJ7nAOcAbx+yWUyvYf1ujeZYcvdmM7seWEK498797r7OzH4ILHf3Jwn/4/m9mW0knF0vS1zEnQt4LjeY2aeAZsLncmXCAu6CmT1MuPdHoZmVAj8g3ICGu/8P4TW6zwM2ArXAVYmJtHsBzuXTwD+bWTNQB1zWT390zAG+AKyJ1F8DfBeYCEn3vQQ5l2T5XsYAD5hZOuHEtdDd/xTPa5imuRARkTapXn0kIiKHQUlBRETaKCmIiEgbJQUREWmjpCAiIm1SukuqSKyY2QhgaeTlaCAEVEZe17r7GQkJTCTG1CVV5DCZ2S3APnf/WaJjEYk1VR+J9JKZ7Yv8PdPMXjGzhWb2jpndbmafi8yHv8bMjo5sN9LMHjOzNyKPOYk9A5EDlBREYusk4EZgGuFRtVPcfTbwv8C/RLa5C7jT3U8BLo68J9IvqE1BJLbeaJ0PyMzeBf4cKV8DnBV5fg5wXLsp74eYWX5k7n+RhFJSEImthnbPW9q9buHA/29pwOnuXteXgYkEoeojkb73Z+D61hdmNj2BsYgcRElBpO/dAMyKrJD3FvDVRAck0kpdUkVEpI3uFEREpI2SgoiItFFSEBGRNkoKIiLSRklBRETaKCmIiEgbJQUREWnz/wFFe6/q46GRcAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ret2 = ensemble_simulations(\n", " duration, ndiv=20, y0=y0, volume=V, model=m, solver=('bd', Integer3(8, 8, 8), 0.1), repeat=N)\n", "ret2.plot('o', ret1, '-')" ] } ], "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" }, "widgets": { "state": {}, "version": "1.1.2" } }, "nbformat": 4, "nbformat_minor": 1 }