{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Hodgkin-Huxley Model\n", "\n", "- A.L. Hodgkin, A.F. Huxley, \"A quantitative description of membrane current and its application to conduction and excitation in nerve\", J. Physiol., 117, 500-544, 1952." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [], "source": [ "%matplotlib inline\n", "import numpy as np\n", "from ecell4 import *" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [], "source": [ "Q10 = 3.0\n", "GNa = 120.0 # mS/cm^2\n", "GK = 36.0 # mS/cm^2\n", "gL = 0.3 # mS/cm^2\n", "EL = -64.387 # mV\n", "ENa = 40.0 # mV\n", "EK = -87.0 # mV\n", "Cm = 1.0 # uF/cm^2\n", "\n", "T = 6.3 # degrees C\n", "Iext = 10.0 # nA\n", "\n", "with reaction_rules():\n", " Q = Q10 ** ((T - 6.3) / 10)\n", "\n", " alpha_m = -0.1 * (Vm + 50) / (exp(-(Vm + 50) / 10) - 1)\n", " beta_m = 4 * exp(-(Vm + 75) / 18)\n", " ~m > m | Q * (alpha_m * (1 - m) - beta_m * m)\n", "\n", " alpha_h = 0.07 * exp(-(Vm + 75) / 20)\n", " beta_h = 1.0 / (exp(-(Vm + 45) / 10) + 1)\n", " ~h > h | Q * (alpha_h * (1 - h) - beta_h * h)\n", "\n", " alpha_n = -0.01 * (Vm + 65) / (exp(-(Vm + 65) / 10) - 1)\n", " beta_n = 0.125 * exp(-(Vm + 75) / 80) \n", " ~n > n | Q * (alpha_n * (1 - n) - beta_n * n)\n", "\n", " gNa = (m ** 3) * h * GNa\n", " INa = gNa * (Vm - ENa)\n", " gK = (n ** 4) * GK\n", " IK = gK * (Vm - EK)\n", " IL = gL * (Vm - EL)\n", " ~Vm > Vm | (Iext - (IL + INa + IK)) / Cm\n", "\n", "hhm = get_model()" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Vm>m+Vm|(1.0*((((-0.1*(Vm+50))/(exp((-(Vm+50)/10))-1))*(1-m))-(4*exp((-(Vm+75)/18))*m)))\n", "Vm>h+Vm|(1.0*((0.07*exp((-(Vm+75)/20))*(1-h))-((1.0/(exp((-(Vm+45)/10))+1))*h)))\n", "Vm>n+Vm|(1.0*((((-0.01*(Vm+65))/(exp((-(Vm+65)/10))-1))*(1-n))-(0.125*exp((-(Vm+75)/80))*n)))\n", "m+h+n>Vm+m+h+n|((10.0-((0.3*(Vm--64.387))+((m**3)*h*120.0*(Vm-40.0))+((n**4)*36.0*(Vm--87.0))))/1.0)\n" ] } ], "source": [ "for rr in hhm.reaction_rules():\n", " print(rr.as_string())" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEPCAYAAABoekJnAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXm8VXW5/9+PIDMIiAgqTgWGJjmhpplHM+OaVxvM4WZq\neuumddXqV0lWQqVNTmV5bw4ZjqVWDhUqogfRSlDxiiICCioos4yHMz+/P757cTabPaw9rL3Wd53n\n/Xqd1zl77X3Wfp71HT7f5zuKqmIYhmEYlbJD3AYYhmEYfmNCYhiGYVSFCYlhGIZRFSYkhmEYRlWY\nkBiGYRhVYUJiGIZhVEVihUREeojIHBF5OPN6qIhME5EFIvKYiAyO20bDMAwjwUICXALMA4KFLpcB\n01R1DDA989owDMOImUQKiYjsAZwE3AJI5vIpwJTM31OAT8VgmmEYhpFDIoUEuA74FtCZdW1XVV2R\n+XsFsGvdrTIMwzC2I3FCIiInAytVdQ5d0cg2qNvXxfZ2MQzDSAA94zYgD0cBp4jISUAfYJCI3AGs\nEJERqrpcREYCK/P9s4iYwBiGYVSAquZtvJcicRGJqn5XVUep6j7AmcATqvoF4CHg3MzHzgUeKHIP\n+1HliiuuiN2GpPzYs7BnYc+i+E81JE5I8hB4+FPg4yKyADg+89owDMOImSR2bW1FVWcAMzJ/rwVO\niNciwzAMIxcfIhKjQhoaGuI2ITHYs+jCnkUX9ixqg1TbN5Y0RETT5pNhGEbUiAha4WB7oru2DMMw\nokKkojozFdS6sW1CYhhGt6U79l5EIaA2RmIYhtHNuP7667n55ptrdj8TEsMwjG7G8OHD2bx5c83u\nZ0JiGIZhVIUJiWEYhlEVJiSGYRgJYsKECVxxxRXbXX/wwQcZOXIknZ2def4rXkxIDMMwEsR5553H\nnXfeud31O+64g7PPPpsddkhetZ08iwzDMLoxp556KmvWrGHmzJlbr7333nv87W9/45xzzuG8887j\noosu4qSTTmLgwIEcc8wxLF++nEsuuYQhQ4YwduxYXnzxxbrabEJiGIaRIPr27cvpp5/O7bffvvXa\nvffey9ixYznwwAMBuO+++7jyyitZvXo1vXr14sgjj2T8+PGsXbuW0047jW984xt1tdmExDAMIw8i\ntfmphHPPPZf777+f1tZWAG6//XbOPffcre9/5jOf4eCDD6Z37958+tOfpn///px99tmICKeffjpz\n5sypxSMIja1sNwzDyEOci96PPvpohg0bxl/+8hcOO+wwZs+ezQMPuCOYRIThw4dv/WyfPn22ed23\nb182bdpUV3tNSAzDMBLIOeecw+233878+fOZMGECu+yyS9wmFcS6tgzDMBLIOeecw7Rp07jlllu2\n6dZK4v5gJiSGYRgJZK+99uLoo4+mqamJU045Zet1Edlm48Xc18G1emLnkRiG0S3JnL8Rtxl1R0S4\n6667WLlyJZdeeuk21ys9j8QiEsMwDKMqTEgMwzCMqjAhMQzDMKrChMQwDMOoChMSwzAMoyoSKSQi\nMkpEnhSRV0TkZRG5OHN9qIhME5EFIvKYiAyO21bDMIzuTlJXtrcBX1fVF0VkAPC8iEwDvghMU9Wf\ni8h3gMsyP4ZhGGVT7/UWaSWRQqKqy4Hlmb83icirwO7AKcCxmY9NARoxITEMowKCNSTXX3/9NntV\nGeWTSCHJRkT2Bg4GngV2VdUVmbdWALvGZJZhbOXWW2HcOBg/Pm5LKmfZMhg5EhJ4ZlJoVq6EHXeE\nIUPK+7/+/fuzcuXKaIwqk7lzYdddoR661r9//5rdK9Er2zPdWjOAH6nqAyLynqoOyXp/raoOzfkf\nzT6msqGhgYaGhnqZbJTJL38J/frBl74UtyWVIwJHHw1PPx23JZUjAj/4AUyeHLclldOnDxx6KDzz\nTNyWVI4IHH88TJ8e/Xc1NjbS2Ni49fXkyZMrXtmeWCERkR2BvwJTVfX6zLX5QIOqLheRkcCTqvqB\nnP+zLVI8QgQGD4b33ovbksoRgYMOgjofAVFTRFxENWtW3JZUjggceCC89FLcllROnHkpdVukiBsB\nuxWYF4hIhoeAYBvMc4EH6m2bUXsGDYrbgurZvDluC6qnuTluC6pnwIC4LaiedevitqB8kjpGcjRw\nNvCSiATaPBH4KXCviFwALAFOj8c8o5b06hW3BdXT1ha3BdWTOYzPazZujNuC6unsjNuC8kmkkKjq\n0xSOlk6opy1G9KxfH7cF1dOjR9wWVE8ahMTH1nwuPs5ITmTXltG9SEMr0sfCn0sahKSpKW4Luicm\nJEbspKFvvqUlbguqJw1CsmVL3BZUj49zhUxIjFjp3TtuC2pDGlrCaRASnxslwdiIj3nJhMSIlT59\n3G8fW2Hgd+EPCJ6971FVz57+5iOAjg73e9OmeO2oBBMSI1aCldS+toY7OpwPLS3Q3h63NZXR2uoq\n4dbWrsrMR3pmpg75nJd23NHZ71teMiHxnBdeiNuC6ggqYh9bYeDs79nTrV/wdS1JczP07et88Hni\nQ0eH88PXdOjsdLP/BgzwrzyYkHjMihVuSwifW5EdHbDTTv4W/o6OrsLvayXc0uLGqgYN8tcH6MpL\nvlXCAUFeGjjQv3QwIfGYoMAsWxavHdXQ3u62SPG98PvYigwIoiofK7AAVdeiHzQoHY0S3/KSCYnH\nbNjgfi9eHK8d1ZCmVqTPPuywg6uEgzzlG4EPvqeDRSRG3QkKvc8bHnZ0pKcV6VvhD0ibDyYk9ceE\nxGMCIfG14Ki6n0GD/PUhTRVY//7+TmPO9iENjRLf8pIJiccEe1T5lukCfC44AWlozQezhawSjhef\n85IJiccErUdfC057e9fUWV99CPrmfZ7+G/jgs5CkJS/5GhmWFBIROV1EBmX+/r6I/EVEDoneNKMU\nwaKlNBQcXyuwoDXfr5+/PqQhHdIUkfiYDmEiku+r6gYR+QjwMdyBU/8TrVlGGNIiJGkp/L61IgN8\nrsAC0uBDZ6eLDH1slIQRkmC528nAzar6VyAFRxH5T3u7y3RpqIR9KzgBafLBxwosIE3p4GOjJIyQ\nLBORm4AzgL+JSJ+Q/2dETHu7W4Ph28BcQHt7VwXmW8EJSEMFFrSEffYhWwwtL9WfMIJwOvAIcKKq\nrgOGAN+K1CojFB0d/q8K79nTz4ITYK35ZJCmvOSjDyWFRFU3A6uAj2QutQOLojTKCEdathfxseAE\n+NwdEZCGdAiiW5998DmqCjNraxLwbWBi5lIv4I4IbTJCkhYh8bHgBKShEk6TD5aX4iFM19angVOB\nzQCqugwYGKVRRjjSIiQ+FpyANPiQvSDRKuH48NmHMELSoqqdwQsR6R+hPUYZ+D5Gkj3Y7lvBCUhL\nSzhNg+3mQ/0JIyT3ichvgcEi8mVgOnBLtGZFx8aN8PzzcVtRG4KIxNdZW9kDpD5Xwr62IgPS4oPv\necnnyDDMYPsvgD9lfsbgFij+KmrDCiEiE0RkvogsFJHvlPv/N90Ehx0WhWX1J5j+63PhT0MFZj7E\nT1oG232NDHuG+ZCqPgY8FrEtJRGRHsCvgROAZcBsEXlIVV8Ne4+BmdGdjRu7/vaVYEFie3vXXkM+\n4XMoH5AGH9K0aaPP6eCzoBesekRkE6AF3lZVHRSNSUU5HFikqksAROQPuIkAoYVk3Tr3e/16/4Uk\nN5wfFEeKVEF2K9K3UD4gLdN/fd2aIyAt6eDreFtBIVHVAfU0JCS7A29nvV4KHFHODdascb99PQku\nmyAKCTKeb0ISCGFQgamCSNxWlYfPrciA3ErY0iEecn3wKR1KdoaIyJ75rqvqW7U3pySFIqRtmDRp\n0ta/GxoaaGho2Pp67Vr329cB6mx87xcOCs6OO7oWcVsb9PJsF7fAh8Du1lZ/fQjSorkZ+vaN26ry\nyM5LnZ0uL+24Y9xWlUdueWhthd69o/u+xsZGGhsba3KvML3qf6erAu8D7AO8BhxQEwvKYxkwKuv1\nKFxUsg3ZQpJLEDKmMSLxjaDgQJcY+loJQ1dk5ZsPwRgJdKWDb0ISlAWRrshqp53itqo8cstDU1O0\nQpLbyJ48eXLF9woza+uDqnpg5mc0bpziXxV/Y3U8B4wWkb1FpBduI8mHyrlBS4v7nYaIxPf9hfJV\nwr6Rr/D7RjBGApaX4sRnH8qe56OqL4hIWeMStUJV20Xka8CjQA/g1nJmbIETkiFD0iEkvkckQdcc\n+F0J57bmfSONPlheqi9hxki+mfVyB+AQXBdTLKjqVGBqpf/f0uIW8W3ZUkOjYiINYyTBlGWfffC1\nFRngcwUWkAYfgu38wT8fwkQkA+kaI2kH/opbnOglLS2u77S5OW5Lqsf3iCQtlXB24fc9HdLgg895\nKdsHn9KhpJCo6qQ62FE30iQkaRoj8bUCyzdQ7Rtp8CF7Qa6vecnnqCrMNvLTRGRw1uuhIvJotGZF\nR5qExPeIJHuMJA2tSN8Kf0DaBtvNh/oTZtPGXTInIwKgqmuBXaMzKVrSJiS+j5H4HpGksUvFfIgH\nn8tDGCHpEJG9ghcisjfQWfDTCSdtQuJzRJK2wXbfCn+Azy3hgLSlg29iGGaw/XJgpog8lXn9UeDL\n0ZkULYGQrF4dtyXVkz1Gsnx53NaUj88FJyANlXAaxkgsL8VLmMH2R0TkULr2tLpUVb2thgMhWbrd\nenj/SENEkqZWpG+FPyANYyS5g+0++uBzXgoz2L4DMAE4RFX/CvQTkcMjtywi0ta15fMYSdoG29Mi\n6L6ng6+NktzI0CcfwoyR3Ah8GDgr83pT5pqXBEKSlgWJaYpI0lCBmQ/xkJaurSAy9M2HMGMkR6jq\nwSIyB9ysLRHxbF9Nh6rbUXPQoHREJGlYR2Jz/+MnbWMklpfqT5iIpDVzMiEAIrILns7aamvrasGn\nQUjSFJH41gILSIsPvo+R+FwJB/gshmGE5AbgL8BwEbkKeAb4SaRWRURLi9uWuU+f9AhJWsZIfCs4\nAT4X/oA0+JA92O6zoPvaKAkza+tOEXke+Fjm0qnl7ribFNImJNknDPpY+NPYijQf4iENYuhzOhQ7\ns31o1ssVwD2Zv1VEhmZWuHtF2oQkaIX5lukCfG6BBfhc+AM6OrpOE0xDOqTBB9/yUrGI5AUKH22r\nwL61Nyda0iokPXv62wJL02C7rxVYGgfb0+CDT+WhoJCo6t51tKMupFFIfB4jSUMl3NHRdRyqb4U/\nII2D7b6mQ5CXfCsPoU5IFJFTcVujKDBDVR+O1KqISJuQBC36vn3duhhVd2a1L+QOtvtUcALS2BL2\n1YdevdzfvlXCAT5HhmFWtv8UuBh4BXgVuFhEbNZWAgi6tnbYwfnk2yLLtLQi0xBV+VqBBaSlUeLr\nIWlhIpJPAgepageAiPweeBGYGKFdkRAISa9ebk1J9tGWPpI75bGpyf32hewxkrRUwk1N/kWG2S3h\n3r1dvsrOWz6Q1kaJL3kpTDWqwOCs14MpPAifaAIhEXG/W1ritqg6fG+FZRecoFuirS0+eyoh24cd\nd3QNk9bWeG0ql+yWsEg68lIghj6R7UMwicaXOiqMkPwEeEFEpojIFOB54KpozYqGQEjAz66gXHJb\n9L61wrILDvgZlWS35sH/Shj89yENYgh++VBSSFT1HtymjX8G/gQcqap/iNqwKMgVEl/UPh+q/vdt\nZ0dU4KcP2a158FMMfa7AArIbVZCOhpVPXXTFFiQeknMpOMFjNxHZTVVfiM6saMgWkt69/R5wDzJd\n0H9qBSce0uBDGqKqtDRKfI3Qiw2nPQe8DKwp8P5xtTcHROQXwMlAK/A68EVVXZ95byJwPtABXKyq\nj5Vz79yIxGchSUvByW1F+uhDGtIhO6ryUQzTIOg+56ViXVvfADYCTcBtwCmqelzwE6FNjwEHqOqH\ngAVkZoeJyP7AGcD+uIO2bswcuhWaNAlJWkN5XwpOQBp9SIOgp8EHn8SwYEWsqter6tG4NSR7ANNF\n5D4ROShKg1R1mqoG29Q/m/lugFOBe1S1TVWXAIuAsk5qTNMYSe70TJ8yXUBuVJUGMUyDD2kQQx99\nyF2O4JMPYQbbXwcexEUK44H9ojYqi/OBv2f+3o2ucRoyf+9ezs3SFJHkCklaWmDmQ/1Jiw++N6x8\njqqKDba/DzgTFwm8BfwRuFJVq540KyLTgBF53vpusP2KiFwOtKrq3UVulXc9y6RJk7b+3dDQQEND\nA5A+IfE1DA7wOZQPyFeB+VL4A3xuCQfki25986Hegt7Y2EhjY2NN7lVssH0hMBd4ANgA7AlcKCIC\nqKpeW+mXqurHi70vIucBJ9F1BgrAMmBU1us9Mte2I1tIskmTkKRloNp3H6wCSwZp9SHKhlV2Ixtg\n8uTJFd+rmJD8kK4W/4CKv6FMRGQC8C3gWFXNruofAu4WkWtxXVqjgVnl3LulxZ3XDv5P/803RrJ+\nfXz2VEJaCn/aulR8TYc0jFX5Gt0W20Z+Uh3tyOYGoBcwzQU//FNVL1LVeSJyLzAPaAcuUtWytmpp\nbU1PRJJvjOSdd+KzpxLS0D2XlmnYuT5s3BifPZWQBjH0ObpN3LZsqjq6yHtXUcX2LGnq2kpDJezz\n4GJAvgrMt8gw34LE5cvjs6cS0iAk+SKSdevis6ccPN77tnzSNP03DeMLaSz8vqZDGgbbfV9X5XN0\nW1BIROSSzO+P1M+caElbRJKGvvk0Fn7ffEiLoKfBB18bJcUikvMzv2+ohyH1IM1C4lOmC/C5BRbg\n8wBpQFq7GH3zwefyUGyMZJ6ILAR2F5G5Oe+pqo6L0K5IyBWSDRvitaca0toS9s0Hnwt/QL4xEt/S\nIa2ztnzxodisrbNEZARuRfu/Ax6c01WctO3+63tEksaWsK8VmO9jJBaRxEvRWVuquhwYJyK9gDGZ\ny6+pqmfn2DnS3LXlU+slIC3dQr77kC8vpcEH38uDTw2rktN/RaQBmAK8mbm0p4icq6ozojQsCtI0\na8vnOecBaWjN+9yKDEiDkKQhuvW5uzrMOpJrgRNV9TUAERkD/AHIPfgq8VhEkizSUAmnpQJLm5D4\n6oOv6RBmHUnPQEQAVHUBCVzIGIY0CYnPYXBAGgp/WgTd1woswKLbeAkjCM+LyC3AnbgB98/jTk/0\njjQJSW7h33FHd+xuW5v72wfSUPjTIobZeSZIB9Wuo5yTThrSwefGYZiI5ELgVdwBV/8NvJK55h1p\nmrWV23oBvzIe+B3KBxRaVFneLnDxktso6dHDlY8tVR8YUT/SEFXlG/fcssWPvFQyIsnswHtN5sdr\n0hyRQFe3yuDB8dhULrkFp1cvVzH7FFXl+tCzp/tpaXF5zAcK5aXNm11l5gO5EUnv3s6vfL4lldxG\nSY8erkw0N0PfvvHZFYZuvdeWz0KSm+nAz4gku/CL+DfGkOsD+NcazlfZpiEv+dZV6nMvQ7cSkubm\nrlZiGqb/5iv8PhWcNFTC+QTdNx+KRSS+kOa85EOZLiokItJDRK6ulzFR0pZZQhl0mfgekeRrvfhY\ncPK1wHwoOAE+tyIDinWT+kLuhAHwMx18LdNFhURVO4CPZI7X9ZotW7btZ0yDkPjeHZFvLKR/f9i0\nKR57KsHnVmRAGiKStrbtfRgwwC8ffI5uwwxDvQg8KCL3AUHxUFX9c3Rm1Z6mpm0HDn0XkjSE8vmE\nxLfC73MrMiAtQpKvUeKTDz5Ht2GEpA+wFjg+57pXQpIbkfg+/betzc3oyMa3SjgNhT8Ngp5WIRkw\nwKLbehFm+u95dbAjcgpFJD4tusqmtXV7IbHCX3/SMHsurULimw8+R7clZ22JyH4iMl1EXsm8Hici\n34vetNqSKyQ9ezoBaW+Pz6ZqyCckvlXCVviTQVqFxLfy4PMYSZjpvzcD3wVaM6/nAmdFZlFE5HZt\ngd/jJBaRJINCYuhDdwS4Q61g2/NIwK+81Nm5/eFc4JcPkL9R4kt3dRgh6aeqzwYvVFUB784jyY1I\nwL+pptn4HpGo+h+RdHS4qDZf4fclHQqt/PYpHQIfcruofUqHYmLogw9hhGSViLw/eCEipwHvRmdS\nNDQ1bR+R+FRYcmlt9bsS7ux0BT+3JexT4c+XBuCXD2kQkkJb6vjSmocuH3wVwzBC8jXgt8AHROQd\n4OvUYdNGEfmmiHSKyNCsaxNFZKGIzBeRE8u535Yt20ckPhWWXHyPSIoVfp98yE0D8CtfFRISnyYM\nFMpLvrTmwf+8FGbW1uvAx0SkP7CDqm6M2igRGQV8nK5TGRGR/YEzgP2B3YHHRWSMqnaGuWfaIpJ8\nGc8nf4oVfl98KBaRvP56/e2phLRHJIsX19+eSsjXMAR/GlZhZm0NE5EbgKeBGSLySxHZOWK7rgW+\nnXPtVOAeVW1T1SXAIuDwsDfsDhGJby0w3yOSQoXfp3yVZiHxyYdiecmH8hCma+sPwErgM8BpwCrg\nj1EZJCKnAktV9aWct3YDlma9XoqLTEKRtoikUNeWL/6kofAX6o7wSQzzbS0Cfs08KySGPqVDsejW\nh/IQZmX7CFX9UdbrH4vIGdV8qYhMA0bkeetyYCKQPf5RbLlg3iNfJk2atPXvhoYGGhoa8s7a8qnS\nysUikvhJw2B79tEK2fhUNtLQKImja6uxsZHGxsaa3CuMkDwmImfRFYV8Dnismi9V1Y/nuy4iHwT2\nAf4vs0/kHrijfo8AlgGjsj6+R+badmQLScCWLTB8+LbXfMpouVhEEj++D5BCOrrn0tAoiSMvBY3s\ngMmTJ1d8r4JCIiKb6GrxXwrckfl7B2Az8M2Kv7UAqvoysGuWDYuBQ1V1rYg8BNwtItfiurRGA7PC\n3rs7dG355E8aCr9FJMkgDY0S3wfbCwqJqg6opyEF2Np1parzROReYB7QDlyUWRwZijQKic+rwtMw\nyJuG1nxra3qFxKfy4Ptge6jTjEVkHLB39ufrsY28qu6b8/oq4KpK7rVxIwwatO01nwpLLsUiEh82\nokxD4U/DYHtLi/9imIYFicUiEh98KCkkInIbcCDwCpC9ZsOrbeTXrYOddtr2Wv/+sHp1PPZUS769\nw4KVsYVamUmipaXr2ONsfKrAfJ9pA+no2iqWl3wRdN+7ScNEJEcAB5TTjZRE1q+HwYO3veZTYckl\nn5BAVyWWdCFpbi4uJL5EVflakcEebp2d228BkzQKtYR79XL2F2rtJ4k05KVC6dC3r/Mv37k3SSJM\nNp+NW03uNYUikrQJiS+tsEKFv0cPJ4JbttTfpnIp1IrcYQeXNj6swygUkYj4Uz7SkJcKNUqCdEh6\nXgojJLcB/xSRBSIyN/OTu1gw8aQtIsm3Lgb86VYpVPjBn3C+UCUM/uStYt2gvuy3lYa8VCgiAT8a\nh2G6tm4FzgZeZtsxEm9QdUJiEUly2LKlcOH3JV0KpQF0VWC77pr//aRQaLAd/EmHYkLikw+F0sGH\nxmEYIVmpqg9FbkmENDe7EDE3s/mSyfJRaowk6aShFVlKSHxIhzREVWnJS/l6GMAPH8IIyRwRuRt4\nmK5TErUe039rRb5oBPwpKLmo+h+RpKEVWUxIfEmHYl1bvqRDWvJSISHxIS+FEZJ+OAHJPf/DGyFZ\nt2778RHwJ5Pl0tbmBnR9XtDX3Fy6Wyjp+N6KhHR0MaYhIsm3YDrAh+g2zHkk59XBjkhJW0RSKtP5\nUHBKFX4f0qVUROKDD5s3O1vz4YsPaRAS36PbsAsSs1EAVT0/EosiIG1CsnkzDByY/z1ffCrVHZH0\nggOu8A8dmv89XyqwzZthyJD87/mUlwYU2NDJFx+amgqngw95Kcz0378Bf838TAd2wm3a6A2rVsGw\nYdtfD+Zn+7bUcsOG7bd7CfCpNe97K7JYZOhLBbZ5c/G+eR98SENe8n3iRpiurfuzX2cG3p+JzKII\nWLECRuQ5/aRHD7egrFh/fRLZsKF4ROJDwdm4EfbYI/97PlVgvncxpqFra+NG/yP0QuvCwI8yXckG\nDmOAXWptSJSsWFF4Pv+AAS4j+kQaIpJihd+XSjgNg+1NTf4LyaZNhbu2fEkH3yOSMGe2bxKRjZmf\nDbhpwN+J3rTasXx5YSEZNChdQuJD6wXS0YoslQ4++GARSTLYtKl4OiS9TIfp2krCuSRVUahrC1wG\n9E1I8m2JH+BTwSkWkbz9dn3tqYRSkeHixfW1pxJKCcmyvGeQJotS0a0PeanQhCDwI6oqdkLinsX+\nUVXfqr050VAqItmwob72VMt77+VfFwN+ZDpwhd/3mTZpiEhKjbf54EOxRokvPmzYUFxIku5DsYjk\n72SdUJjFLpmfBG9qvC3Fxkh8FJJVq2CXAqNUvhScNIyRrF9fPCLxwYd16/yf/lusUeJLOhSLSLzu\n2lLVD2a/FpG9gcuAE4ArI7WqhrS1wZo1MHx4/vd97NpavRr22iv/e74UnHzb+gf4UoGlYayq0K4P\n4E86+N4tBMUbJcH5NkkmzGD7GBH5PfAI8DwwVlVviNqwWvH22zByZOHDeXyMSFavzr8uBvwo/KpO\n3Av54EPhD3aULlb4k34ORmdn/uMVAnw4B0PVdfUWWhjqS3koJoY+CEmxMZIDgcuBA4CfAxeoake9\nDKsVS5bAPvsUft9HIVmxonDXlg+VcFNT1+FP+fCh8G/Y4BbBFTvLI+mFf9Mmlwb59mwDP9Jh/Xr3\nrAs1FH0oDxs3ui3kfc5LxcZIXgSW4la0Hw4cLl3nVaqqXhyxbTVh8WLYe+/C7/s4/XfZMr8X861Z\nAzvvXPh9H7qFio1TgR+Ff82awi15SEde8iEdVq0q3PUOfkS3xYTkgszv3AF3yXMtsZSKSAYOhDff\nrJs5VdPZCe+8A7vtlv99Hyrh1atLF/6kF5wwQpJ0H5YvLzwtHvwRklJi6IOQFMtLPhzbXGyw/fd1\ntCMy3ngDPvGJwu/71rX1zjuu4BTaW6hPH3fGRGen6z5KIsWEEFzBSXolHKYVmfTCX2x9FfghJO++\nWzwv+ZAOK1f6H90msqoRkf8WkVdF5GUR+VnW9YkislBE5otI7vkoeXn5ZfjgBwu/71vX1rx5MHZs\n4fd32MH1tTY318+mckmDkLz5JuxZZKWVD4X/3XeLHwXsgw/LlvkvJEuXwqhRhd/3wYcwB1vVFRE5\nDjgFGKf2ykUnAAAbEUlEQVSqbSKyS+b6/sAZwP7A7sDjIjJGVQueI9/aCgsXwv77F/6+gQP9ikjm\nzYMDDij+mSDjFdoHKm6WLYPddy/8fiAkqu6I5CTyxhvFu0yD7ogk+7BkSfHxQx/SoZSQ9OnjjhNO\ncoQeJi8lPR2S+GgvBH6iqm0Aqroqc/1U4B5VbVPVJcAi3CSAgsyf7wpKoW4g8K9r65VXwglJklv0\nCxfC+99f+P1gV+aWlvrZVC6LFxcv/D74kIZ0eP11eN/7Cr8vkvwI9403YN99C7/vQzqEWUeyn4hM\nF5FXMq/Hicj3IrRpNPBREfmXiDSKyGGZ67vhZpEFLMVFJgV56SUYN674l/nWtTVnDhx4YPHPJH1w\n7tVX4QMfKP6ZpBf+UkICyRf0RYtg9Ojin0l6t8r8+aXzUtJ9WLy4uJBA8n0I07V1M/At4H8zr+cC\n9wA/rvRLRWQakG+Y7/KMTUNU9UgRGQ/cCxR6zHlnj02aNAmAadNg7NgGoKGgLT5FJGvXwoIFcHjR\nOCzZma6z0/kQVkgKbd8RJ6rhhaTYyXdx0tnpWvPFIhJItqAHeWm//Yp/LsnlQbV0RAJdPhSboVYu\njY2NNDY21uReYYSkn6o+G6whUVUVkbZqvlRVP17oPRG5EPhz5nOzRaRTRIYBy4DsIak9Mte2IxCS\nf/0LPvWp4rb4NEbyxBPwkY+4xUvFSHLhf+stV7EW2mcrIMmFf/Fi1wApJRBJ9uHtt50PYdIhyXlp\n550L77MVkOR0eOstZ1+YvFTrdGhoaKChoWHr68mTJ1d8rzBjJKtEZGu7RUROA96t+BtL8wBwfOa7\nxgC9VHU18BBwpoj0EpF9cF1gs4rdKEzXVrDXlg/H7d51V2lhhGQXnNmz4eCDS38uyWL47LNwxBGl\nP5fkLsZ//ct/H55/vnT5hmSXhzSkA4SLSL4G3AR8QETeARYDn4/Qpt8BvxORuUArcA6Aqs4TkXuB\neUA7cJFq4ep/1Sr34ItNqwO3PUSfPsVPiksCb78NjY0wZUrpzya5Em5shKxGUEGS7MOzz5buXoRk\nV2DPPANHH136c0lOhxkz4NhjS3+uX7/krocJ2yhJcl6CEBGJqr6uqh8DhgH7qerRmVlTkZCZlfUF\nVT1QVQ9V1cas965S1fer6gdU9dFi95k717VWwkyX86F76wc/gAsvLLxJYDZJznRPPBGu8Ce5Aps2\nLZwYJrlbqLHRdZOWIsl56cknwwtJUn2YNs1/HyBERCIifYDPAnsDPcQNlqiq/jBi26oiTLdWQDDg\nPnJktDZVysMPu0Lzf/8X7vNJrcDmznUtw0MOKf3ZpBacBQvcbrPjx5f+bJJ9WL06fJdKEvPSwoWu\n1+HQQ0t/NqnpsGiRS4cjjyz92aSW6YAwXVsPAutwW8gneL30trz0UrgEgmRPAX7+efjP/4Q//anw\nNtO5JLU/dcoUOOuscAvDklqB3X47nHZaOB+SWoHddRd89rN++3DnnXD66YV3Ls4mqT7ccUf4dEhq\nmQ4IIyS7q2qR3aqSyUsvwZe/HO6zSZ0CPHMmfO5zcNNN4bohApJYcNatg9tuc8IYhiQKyebNcPPN\nrm8+DElMhy1b4Le/dV1bYUhiOgQ+PPFEuM8ncePGLVtcuZ4+Pdznk5iXsgkza+sfIhKykygZtLe7\nrUSK7bGVTdLGSDo74YYbXMv39tvh1FPL+/8kFv7Jk+Eznym+JUc2SQzlr7kGjjuu9BqYgCS2Iq+5\nxjVKwvqQxArs+uvhqKOKb32UTRJ9uPZaN9nBZx+yKXaw1StAJ+5s9i+KyGIgWKSvqppYcVm0yO2/\nU2p+eUCSurYWLHCRVHOzm1lTasFYPvr1S44/4Kb83nWX294lLEmrhBctgl/9yvkSlqSJ4eLFcN11\n8Nxz4f8naY2SxYudGJabDknKS2+84dLh2WfD/0/S8lIuxbq2dgMOwp0/4hXlDLRDMrq22tpcAbn6\navje9+C//9vtsVMJffu6ramTwOrVrnvuf/+3+FbZuSSpAmtpgf/4D7jiitKr2bNJUgXW2urGpy6/\nvDwfkpQOzc0uL5XrQ5LSoaUFPv9550OxPcJySVrDKpdiQrJEVT068qmLcoUk7q6tF15wA+rDhrmW\nVjmFJB9JKTibNrkFlGee6bq1yiEpFVhnJ5x3nuuS+9rXyvvffv3ib6CA8+H8892Oy1//enn/m5S8\n1NnpIvV99oFLLy3vf/v1c9sLxY0qXHCBmx16ySXl/W/SehlyKSYku4jIN8gfkaiqXhuRTVXz0kuu\n8Iclrq6tLVvc2MFtt8HPfw7nnFObbaKTEAY3Nbmxnf32g6uuKv//k1D429rgK19xi0Eff7z8tOnX\nz51CGCcdHU4A33jDDeyW60PfvvGLoaqL0BcvhkceqSwd4hbDjg746lfd/mbTp5e/pX2/fu4gsqRS\nTEh6ACV24kkmlXRtvfNOdPbk46mnXBRy8MHO3mIHDJVL3GHwypVwyiluQPemmyo7ByJuHzZvdtNL\nVV3lVewogkLE7cOWLa5LbsMGmDrV2VMucVfC7e2uAp4zxy3eq2T3ibh9aGpy6bB5Mzz6aGXnBCWh\ncViMYkKyXFUr38UrRlavLr2bZjb1HCPZuBEuuwweeABuvLH8GVlhiLPgvPYanHSS6weePLnyCCvO\nrq2VK+Hkk92sv9/+1p0FUQlxFv61a+Hf/911yf3xj6U3+ixEnOnQ1OS6RVtaXCu+1AaThYizPKxZ\n49Jh333h3nurS4e4o6piJPFgq6oZO7a8VnCwcWPUPPqoq5yam90RwFGICMRX+GfOhI9+1A0k/vCH\n1XXTxeXD66+7aZmf+ATcemvlIgLxVWBvvul8OPpot+it0soL4vNh9Wo4/ni3K+7DD1cuIhCfD0uW\nuDQ45hg3jd/HdAhLser2hLpZUWPCrh8JGDQI1q+PxhZwFeJXvwr/9V9uQdutt0Z7RkUcme6ee9wq\n3TvvdAO71RJHa372bFfov/lN+NGPqh+viiMdXnzRVV5f+Yobd6v2eNk4BP2NN5wPH/sY/P731VXA\nEG86XHQR/Oxn1adD0oWkYNeWqq6ppyG1pFwhGTrU7Z8UBS+/7KZdHnCAy1yDB0fzPdnUs/Crwk9/\n6qb3Tp9e+vTGsNS7Aps61U12uOWW2kWK9S7806e7vPbrX7vxnVpQbx9mzXIz/b73PVcJ1wJLh+gJ\ns0WKd5QrJEOGRDND6JZb3HjIL37hZpHVYkZWGOqV6drbXWGfPRv++U+3CLRW1LNP+LbbYOJEeOgh\n+PCHa3ffehb+u+9202Lvuy/cbrJhqaegP/AAfOlLLj1OPrl2963nNvJ33eWmWPucDpVgQoKLSGop\nJK2tbp54YyM8/XT47ShqRT26hTZuhDPOcBHJU09V14edj3p1bV11Vdf+WaWObC2XehX+3/zGdZ88\n8UT5eb8U9fLhV79yPkydCocdVtt710vQf/lLt6g4inSwiCQGym0ZDxzoBsDb2qobXAU31/u001yU\n8+yz4c4PqTVRt+bfece1GA891M08q/aZ5SPqCkzVdZ88+CD84x/RHCFQj8J/9dUuDWbMqH4haz6i\n9qGzE/7f/3NTrJ95JvxebOVQj00bf/ITN/Y5cybstVft7590IUnlrK1yu5BE3NhFteMkCxe6Mx6O\nO86F6XGICESb6YJZTaed5taIRCEiEK2QqLoB9b//3UWNUZ1DE2U6qLqZcTff7CLCKEQEok2H9nb4\n4hfduEhUIgLRp8MVV7hZWTNmRCMiEK0P7e3w7W9Xd49URiSVEHRvDR9e2f+//LKbMjp5sltoGCe9\nernM0d4e7ryGsLz2GpxwAnz/++G36K+UqKKqzk630vu551wXhI+z51Thu9+Fv/7ViUgtF7PmEpUP\nwf5lmzbBY49VtkgvLEFeUq3tOKUqfOc7LpqaMaPyuiMMUQl6W1tXOlSDCUmGasZJnnvOdfVcd52b\nrRE3Il1jDLUau5g3Dz7+cfjxj10rMmqiGCPp6HCDuQsWuC1Poo4Yo6jAVN2g+syZ7tTMYcNqc99C\nRFGBNTW5vdf693cTHHr3ru39c+nZ0/20ttbuuzo73TjoP//p0mHnnWtz30IEgl7LvNTa6hZ8trW5\nHpRKdm8ISGXXViUMGVJZ19bTT7uV3L/9bTJEJKCWFcDcuS4S+dnP6iMiUPsKrL3dTe9dssS1IOvR\n7dizp9vBua2tNvfr7HTrQ2bNctFU1CIC24phLdiwASZMcK33P/4xehEJqGVk1dHh1oQ9/7yb6hu1\niIDrQhapXV5qaXHd06ru9NVq08GEJEMlEcm0aa5ldddd0a1Sr5RaFZw5c1wkct11cPbZ1d8vLL17\nuxZTZ2f19wpaXmvXwt/+Fv6cmlpQq3Rob3dTyOfPd11B9ViPBNu25qtl/Xo48US3pur3v69tt2sp\nap0OCxe6nSrCHn9dC2rlQ3MzfPrTroxVs21LNiYkGYYPL2+n1gcecPtJ/fnPrqJNGrVo0c+e7VqP\nN97opvrWExEXalfrQ3OzW3Hf3u7SrJKNC6uhFoU/6MdevtxNj631VOtS1MKH9evdGOJhh7n8VO1K\n73KphQ/BmS4rV7qJGvVOh1qU6aYmt6HqTju53ShqNVnGhCTDnnvCW2+F++zdd7suhqlTyztLvZ5U\nW3D++U/45CfdospyzxKpFdWOkwSFpl8/t0CsXt0o2VSbDoEQNje78YQoB6ULUW0Ftm6di0TGj3dH\nSNdrYW42tUiH005zXUJxpUO1Pmze7MZyhw93e7DVMiJMnJCIyOEiMktE5ojIbBEZn/XeRBFZKCLz\nReTEWn7vXnu5/vNS3HwzfOtbbrD20ENraUFtqSbTzZzpuuqmTHE7l8ZFNRXYxo1u7GrECNf1GNU0\n5VJUkw7BmS69e8P991c3GFoN1cygW7fORSJHHOEWHcYhIlDbdIijQQLV+bBpkysPe+3lynWtuxUT\nJyTAz4Hvq+rBwA8yrxGR/YEzgP2BCcCNIlIz+/fay+2aWozrr4crr3RrD2q9crXWVFoJP/mki0Du\nvhv+7d9qb1c5VOpD0I0yenT9++JzqdSHQAiHD3ddELXox66USiPDIBI58ki36jsuEYHKK+GNG11k\nnpR0qMSHYILDmDFu0WSlR3gXI4lC8i4QDGENBpZl/j4VuEdV21R1CbAIOLxWX1pMSFTdtNff/MbN\n2x89ulbfGh0DBpS/Nf4jj7ixkPvuc7O04qZ///Lnt69d63aNPfRQN5Ou3n3xuVRS+INW/Jgx0bQe\ny6USMQxE5KijXAMsThGBytIhaJC8//3+pkPgw7hx0ZaHJArJZcA1IvIW8AtgYub6bsDSrM8tBXav\n1ZcOHepmCOXO3OrshG98w81ueOopN5biA0OGuMIclr/8Bc49120Z0tAQmVllMXhweT68+67bKO/4\n4103StwiAuVXYCtXOvsPOywZQgiV+fCxjzkRue66+EUEKvPhuOP8TodVq1xeGj/eNYKj9CEWjRWR\nacCIPG9dDlwMXKyqfxGRzwG/AwrNi8o7u33SpElb/25oaKAhRM0o4lqxs2a5MBDcbJnzz3dnRc+Y\nEe0q6FpTzo7G99zjdiydOhUOOSRau8qhHDF84w03e+6CC9xOvkmovKC8wv/mm86Hs86CSZOS40M5\n0e2bb7pI5Iwzqjshs9aUkw5vveXS4cwzk5UO5fpw4onwuc8VPmSusbGRxsbGmtgWi5CoasEJsyJy\np6oGHSv3A7dk/l4GjMr66B50dXttQ7aQlMMJJ7jT2CZMcBsTnn66W/QV9RYOURBmgaUqXHutazU+\n/njyxn3C7n82d64bz7n8crjwwujtKoewhf/VV10XxDe/6VZMJ4lhw9yRsaWYN8+VnST6EHbjxvnz\nXTp8/etuB4EkEXY7/LA+5DayJ0+u/GT1BARs27FIRIKd/I8HFmT+fgg4U0R6icg+wGhgVi2/+Itf\ndF1Y553n+hQnTHDrRHwTESgtJO3trtKdMsXtfps0EYFwYvjwwy58v/rq5IkIuPn6paKqRx91XXJX\nXpm8ChickKxeXfwzjz7qukST6kOYSvjRR91R0T/6UfJEBMIJeuDDD39YXx+SuNfWl4HfiEhvYEvm\nNao6T0TuBeYB7cBFqrXauMGx++5u1tKMGfCDH8C++9by7vVl+HC3pX0+li93q9R33NFt8RLXLsWl\nGDq0cMHp7HRHyd5wg9u88Igj6mtbWEaOdGM3+Qgiwquvdg2WpK5JKlaBqbqI9he/cFttHHNMfW0L\ny6BBhQVd1U0I+PnPk50OI0YUXjQddzokTkhU9Tkgb7WgqlcBV0X5/R/8YDJb5+Wyzz5ubCeXxx93\ne05dcIHb/jrumSjF2HNPN5Msl3fecVHj5s3wr3/BqFHbfyYpjBzpuq1yWb7cjb+tWuV8iGr78Vqw\nyy5uo8tcVqxw+Wj58uT7MGqU29IolxUrXDqsXOkW4Ua1lX0tGDECXnll++vZ6RCXD0ns2jJqwL77\nugHoYK+q995zW7+fd55b1fqjHyVbRMAViGwx7Ox060IOOcSdiTJjRrJFBGCPPbadVq7q1ugcfLDz\n4x//SHYFDG4acraQqLoJGgcd5H588GHvvbddcJydDh/6kPMhySICLq9nlwdV+MMfXBqMGxezD6qa\nqh/nkqGqOnas6p//rHrNNaojRqh+9auq69bFbVV41qxRHTRItalJtbFR9aijVMePV501K27LwrN2\nrerAgaotLaozZ6oec4zqwQer/uMfcVsWnlWrVHfaSbW5WfXpp50PBx2k+swzcVsWntWrXV4KfDj2\nWP98eO891f79u3z46Edr60Om7qyo3k14m9SohkmT3MDnhz/suog+9KG4LSqPoUPdTLpRo9zA++WX\nu265JMzpD8uQIW4gfY89XD/9ZZe5SR1RrC6OimHDunzYaSc/fdh5Z+fDXnu5gfeJE12Xlk8+DB7s\nJpaMGuU2jLzssuT4IFrb8erYERFNm0/dmdZW1yXx/vf7JSDZtLS4bsYxY5JR6CuhpaUrHXz24fXX\nYb/9/PYhqrwkIqhqRatmTEgMwzCMqoTE0zaeYRiGkRRMSAzDMIyqMCExDMMwqsKExDAMw6gKExLD\nMAyjKkxIDMMwjKowITEMwzCqwoTEMAzDqAoTEsMwDKMqTEgMwzCMqjAhMQzDMKrChMQwDMOoChMS\nwzAMoypMSAzDMIyqMCExDMMwqsKExDAMw6gKExLDMAyjKmIREhH5nIi8IiIdInJIznsTRWShiMwX\nkROzrh8qInMz7/2y/lYbhmEY+YgrIpkLfBp4KvuiiOwPnAHsD0wAbhSR4OjH/wEuUNXRwGgRmVBH\ne72ksbExbhMSgz2LLuxZdGHPojbEIiSqOl9VF+R561TgHlVtU9UlwCLgCBEZCQxU1VmZz90OfKo+\n1vqLFZIu7Fl0Yc+iC3sWtSFpYyS7AUuzXi8Fds9zfVnmumEYhhEzPaO6sYhMA0bkeeu7qvpwVN9r\nGIZh1BdR1fi+XORJ4Juq+kLm9WUAqvrTzOtHgCuAN4EnVXVs5vpZwLGq+pU894zPIcMwDI9RVSn9\nqe2JLCIpg2zDHwLuFpFrcV1Xo4FZqqoiskFEjgBmAV8AfpXvZpU+CMMwDKMy4pr++2kReRs4Evib\niEwFUNV5wL3APGAqcJF2hUwXAbcAC4FFqvpI/S03DMMwcom1a8swDMPwn6TN2qoYEZmQWcS4UES+\nE7c99URERonIk5lFni+LyMWZ60NFZJqILBCRx0RkcNy21gsR6SEic0Tk4czrbvksRGSwiNwvIq+K\nyDwROaIbP4uJmTIyV0TuFpHe3eVZiMjvRGSFiMzNulbQ90ILwwuRCiERkR7Ar3GLGPcHzhKRsfFa\nVVfagK+r6gG47sKvZvy/DJimqmOA6ZnX3YVLcF2kQcjdXZ/FL4G/ZyaqjAPm0w2fhYjsDXwJOERV\nDwR6AGfSfZ7Fbbj6MZu8vhdYGF5UK1IhJMDhuHGTJaraBvwBt7ixW6Cqy1X1xczfm4BXcZMVTgGm\nZD42hW6yiFNE9gBOwo2pBZMvut2zEJGdgGNU9XcAqtququvphs8C2IBrcPUTkZ5AP+AdusmzUNWZ\nwHs5lwv5nm9h+OHF7p8WIdkdeDvrdbCQsduRaXkdDDwL7KqqKzJvrQB2jcmsenMd8C2gM+tad3wW\n+wCrROQ2EXlBRG4Wkf50w2ehqmuBa4C3cAKyTlWn0Q2fRRaFfC+0MLwgaRESmzEAiMgA4E/AJaq6\nMfu9zOy31D8nETkZWKmqc9h2avlWusuzwE3vPwS4UVUPATaT03XTXZ6FiLwPuBTYG1dRDhCRs7M/\n012eRT5C+F70uaRFSJYBo7Jej2JbRU09IrIjTkTuUNUHMpdXiMiIzPsjgZVx2VdHjgJOEZHFwD3A\n8SJyB93zWSwFlqrq7Mzr+3HCsrwbPovDgH+o6hpVbQf+DHyY7vksAgqVidz6dI/MtYKkRUiew+0I\nvLeI9MINFD0Us011I7ND8q3APFW9Puuth4BzM3+fCzyQ+79pQ1W/q6qjVHUf3GDqE6r6Bbrns1gO\nvC0iYzKXTgBeAR6mmz0L3CSDI0Wkb6a8nICbjNEdn0VAoTLxEHCmiPQSkX3ILAwvdqPUrCMRkX8D\nrsfNxrhVVX8Ss0l1Q0Q+gtuS/yW6QtCJuMS/F9gTWAKcrqrr4rAxDkTkWNwWPKeIyFC64bMQkQ/h\nJh30Al4HvogrI93xWXwbV2F2Ai8A/wkMpBs8CxG5BzgWGIYbD/kB8CAFfBeR7wLnA+24rvJHi94/\nLUJiGIZhxENaurYMwzCMmDAhMQzDMKrChMQwDMOoChMSwzAMoypMSAzDMIyqMCExDMMwqsKExDAq\nRER2zmxVP0dE3hWRpZm/N4rIr+O2zzDqha0jMYwaICJXABtV9dq4bTGMemMRiWHUDgEQkYasA7Um\nicgUEXlKRJaIyGdE5GoReUlEpma2NEdEDhWRRhF5TkQeCfZAMgwfMCExjOjZBzgOd/7DnbjDhMYB\nW4BPZjbcvAH4rKoehjuE6Mq4jDWMcukZtwGGkXIUmKqqHSLyMrBD1r5Fc3Hbmo8BDgAed/sJ0gN3\nZoZheIEJiWFETyuAqnaKSFvW9U5cGRTgFVU9Kg7jDKNarGvLMKIl7+FaObwG7CIiR4I7WyZzbrZh\neIEJiWHUDs36ne9v2P6kOVXVNuA04Gci8iIwB3fokmF4gU3/NQzDMKrCIhLDMAyjKkxIDMMwjKow\nITEMwzCqwoTEMAzDqAoTEsMwDKMqTEgMwzCMqjAhMQzDMKrChMQwDMOoiv8P1okj/KFxt8EAAAAA\nSUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "run_simulation(np.linspace(0, 100, 1001), model=hhm, y0={'Vm': -75}, species_list=['Vm'])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# FitzHugh–Nagumo Model\n", "\n", "- R. FitzHugh, \"Mathematical models of threshold phenomena in the nerve membrane.\", Bull. Math. Biophysics, 17:257—278, 1955." ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [], "source": [ "a = 0.7\n", "b = 0.8\n", "c = 12.5\n", "Iext = 0.5\n", "\n", "with reaction_rules():\n", " ~u > u | -v + u - (u ** 3) / 3 + Iext\n", " ~v > v | (u - b * v + a) / c\n", "\n", "fnm = get_model()" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "v>u+v|(((-v+u)-((u**3)/3))+0.5)\n", "u>v+u|(((u-(0.8*v))+0.7)/12.5)\n" ] } ], "source": [ "for rr in fnm.reaction_rules():\n", " print(rr.as_string())" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false, "scrolled": true }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAAEPCAYAAAC6Kkg/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXd81dX5+N8nIYEMRghZhBH2JmFPIYgiIM6CiqtYq/Zr\nHW21PztshdZqbbXVukq17lncCjIVmTKFEAgrzAwyCCMDMs/vj5MLN8kdn3vvZxHu+/X6vExuPvec\nw/E85znPc57zHCGlJEiQIEGCBPGFEKsbECRIkCBBLjyCyiNIkCBBgvhMUHkECRIkSBCfCSqPIEGC\nBAniM0HlESRIkCBBfCaoPIIECRIkiM9YpjyEEJ2FEN8KIXYKITKFEA+4ee9fQoh9QojtQoghZrcz\nSJAgQYI0pYWFdVcDv5RSbhNCRANbhBDLpJRZjheEENOBnlLKXkKIUcDLwGiL2hskSJAgQeqxzPKQ\nUh6TUm6r/7kMyAI6NnrtauDN+nc2AO2EEAmmNjRIkCBBgjTBFnseQogUYAiwodGfkoGjTr/nAJ3M\naVWQIEGCBHGH5cqj3mX1EfBgvQXS5JVGvwfzqQQJEiSIxVi554EQIgz4GHhHSvmZi1dygc5Ov3eq\n/6xxOUGFEiRIkCB+IKVsvEDXhJXRVgL4L7BLSvmsm9e+AG6vf380cFJKWeDqRSll8NHheeyxxyxv\nQ3N6gv0Z7E87P4FgpeUxDrgVyBBC/FD/2e+ALgBSyvlSykVCiOlCiP1AOXCHNU0NEiRIkCDOWKY8\npJRr0GD5SCnvM6E5QYIECRLEByzfMA9iL9LT061uQrMi2J/6EuxP+yAC9XvZASGEbA7/jiBBggQx\nEyEE0s8Nc0ujrYIECRLEzqi4nuaB3gvsoNuqGXLsGPzpT1BTY3VLrOXYMXjpJShzdXroIqKgAP7y\nFygpsbol1lJTA3/4AxQV+fY9qyOi9I6qqqpScrF3b2D9GbQ8miH33AMbN0JtLcybZ3VrrOOxx+Cz\nz2DnTnjxRatbYx1vvgmPPgrbt8P//md1a6zjhRfgP/+BtWvhm2+sbo35PPvss0RFRREefhd/+pPq\ni0AIWh7NkLVr4ZNP1IRZWmp1a6xj82Z44w14/33IybG6Ndbx0UfwxRewbh1s22Z1a6zjm2+UTBw6\nBOvXW90a84mPj6e8vJytW+Ghh+BXvwqsvKDyaGacOgVnz8Lo0XDZZYGvLi5UKishKwvS0+G229Sq\n82Jl507VD/feC88/b3VrrOPQIejZEx5+GJ56yurWWMf27ZCWBrffHlg5zU555OfDyJHw6adWt8Qa\nDh6E7t1BCPj1r+G555T76mJj927o1g0iIuCBB+DVV6G83OpWmU9ZGUgJ0dFw113KIj1+3OpWmY+U\ncPgwdO0Kd9yhrLD9+61ulflICRkZkJoaeFnNTnm8+y5ERsLdd0NxsdWtMZ8DB5TyABg2DBITYelS\na9tkBcXFkFCfvL9HD2WJLVhgbZusoLAQ4uPVYiIuDq69Fl57zepWmc/Jk+q/7dqpBcVtt12c/VBb\nCxUVakwESrNTHp9/Dr//vRKSl1+2ujXmc+CAWnE7+OlP1ar7YuPECYiJOf/7nXfCf/9rXXusorDw\nvBIFNR5ee02tQC8mDh2ClBSlREGNhzfeuPgiEqurlQLVg2anPI4cgV69lKvi5ZdVWNrFREmJWmE6\nuOkmWLFCTSIXEydPNhSSGTNg3z7Ys8e6NllBQUHDVebYsVBXB99/b12brCAvDzo6XTXXv79aZC1a\nZF2bAiEkJIQDBw6c+33OnDn84Q9/8Pq9qipo21anNuhTjD2oq1Ox/UlJMGgQ9O2rIk0uJkpLlX/b\nQZs2cN118NZb1rXJChorj7AwtUF4sbkqHG4rB0LAT35y8fVDaamSBWeakzUqhNB0oLGqKmh5uKS4\nWA2Qli3V7/fee/FFG5WWQuvWDT9zuK4uJldFY7cVqEnzrbeU6X6x0NhtBcrf/9FHF1cAQVlZw0UV\nwA03wKpVyirxFyH0efRAywnyoPJwQ2PT9KqrVJiik3XX7HElJGPHqgG6bp01bbKCxpYHKEu0Rw9Y\nuNCaNllBYWFDNyYoGRk3Dj7+2Jo2WYEruYiOhlmz1CFKf5FSn8csgm4rNzRWHi1bwuzZF5fLxpXl\nIYQy0S8mV4UrywOal6tCC2VlTccDXHyuK3f94JCLC80qj4yMpKKi4tzv+fn5QbdVIDRWHgBz5qiV\nRV2dJU0yHXdCcsst6uzLmTPmt8kKXFkeoFaaq1dfPAEEZ86o0NTGzJgBu3ZdPGcdXFkeoM6EhYTA\nhg3mtykQ0tLSePfdd6mtrWXx4sWsWrVK0/eClocbSkogNrbhZ0OGqMn0u++saZPZNN4wd5CUBCNG\nqFDmiwF3yiM6Wk2cH3xgfpusoKJCnXtqTHi4WlC88YbpTbIEd8pDCLUHdKF5J5577jm+/PJLYmJi\neO+997juuus0fU9Py6NZJUZ0NUCEUCdK33gDJk2ypFmm4spt5eC22+Dtt1X4bnPnxAn3QnLbbSpR\n4AMPmNsmK3CnPEDJxYwZKgNzSLNaRjbFnfIAuPVWGD4cnn1WKdULgWHDhpGZmenz94JuKze4GyC3\n3KJW3BdDkkBPQnLddSppYkGBuW2ygvJy9/0webJKlLh7t7ltsoIzZ9wrj8GDlQtjzRpz22QFnuQi\nJQUGDLhwz3z4QnV1M3FbCSFeE0IUCCF2uPl7uhDilBDih/rnUU/luRsg8fEwYYLK69Pc8WR5REXB\nNdeoLLPNncrK8yHbjWnRAm6+Gd55x9w2WUFFhes9Dwc33wzvvWdee6zCk/KAC9N15Q81Ne4XE75i\nteXxOjDVyzvfSSmH1D+Pe3rR0wCZPbv5+7mrq9XgaNXK/TsO11Vzx5PyANUP77zT/AMpPLmtQMnF\nRx81/0wM3pTHrFkqE0NzvzCrrs7z/OALlioPKeVq4ISX1zQfofE0QK66Sp1zaM7JEh3/fk8Re5Mm\nKbfVrl3mtcsKzp71rDxSU1VfNXeXjSe3FSiXTd++sGSJaU2yBG/Ko21bmDq1+V+WVVvrWS58wWrL\nwxsSGCuE2C6EWCSE6O/p5bIy5ZpxRXQ0TJvWvA9GeXJZOQgNVSdrm7OQSKlW0p6ERAi1UdrcXTbe\n3Fag9gSbez94Ux6gAkmas1zAxaU8tgKdpZSpwPPAZ55e9jZAbrqpefv7y8vdK09nZs1q3unJq6uV\nkgwN9fzerFnq7Etzvu/Em9sKVD98/XXzvuvdXQi7M1OnwtatzfsMUF3dRaI8pJSlUsqK+p+/BsKE\nEO1dvTt37lwOHpzLm2/OZeXKlS7LmzpV3aKVm2tYky3l7Flt/sxRo+D0aZW6pTnibb/DQY8ekJys\n8hs1R6TUZnl06ADjx6v73psrWhZWEREwfXrzDqyprDzGa6/NZe5c9QSCrZWHECJB1J+5F0KMBISU\n0uWW1ty5c4mMnMsjj8wlPT3dZXmtWqloo+a66q6q0hanHhLSvK0PrcoDVD8018zLVVUqm3ALDae5\nbryx+fYDaF9YzZrVvF1XISGJ/OIXzUB5CCHeB9YBfYQQR4UQPxFC3COEuKf+lZnADiHENuBZwOPx\nNi1+zdmzm6/rytdJM6g8YOZMtQ/WHF1XWqwOBzNmwDffNE/XVW2terQoUYfrqrmehWo2ex5SytlS\nyo5SynApZWcp5WtSyvlSyvn1f39RSjlQSpkmpRwrpfR4hY0W5XHppSqfz9Gj+v077IIvk2Zzdl1p\nXWWCujgsMVEdnmxuaNnvcBATo7Ivf/21sW2yAkfwhJbU583ddVVb20xCdfWkqkp1jDe3TVgYXHkl\nfPGF5/eklFTVXljB774oD62uKykl1bUX1gUYvvQDKOtDixVWWVPpf6MswFuYbmOuv15bNGJzlgto\n3lZ5s7E89MSRjkLL6uKaazwnCFy8fzEpz6XQ9q9teXrd0/o10mC8hac2xpuQrD2ylj4v9CHqiSh+\ns/w3gTfQJPyZLD7+2P2BwW3HtpH671RaP9man331M+rkhXGy0Be3FcC118Lixcpyc8X+kv2MfnU0\nEX+J4MaPbqS27sLw9fk6Hi6EqKunnnqKWbNmNfjswQcf5MEHH/T4vYsm2soXtIapAlxxhbrD+eTJ\npn9blr2MH3/2Y9689k3237+ff37/TzbkXBj5misrfUvsNmoUHD/u+rKsjbkbue7D63h6ytPkP5TP\nx1kfszR7qX6NNRBfJ4s+faB9e9dpufce38sV71zBw2Me5vj/O87mvM28m/Gufo01EF/cVqDS+KSm\nwrJlTf9WUFZA+hvp3Dr4Vsp+W0ZxRTHPbXhOv8YaiNZAEgcRESr/mZ1deLNnz2bRokWU1W9S1dbW\nsmDBAm655RaP39PT8mg2WXV98XNHR8PEiWpwzJ59/vPiimJu/+x2FsxawISuEwB4eMzDvLT5JUZ1\nGmVAq/XF10kzJET5dxcuhPvvP/95eVU5N398My9f+TJX97kagHnp83hm/TNM6TFF51brj6/9AGrD\neOFCGDPm/GdVtVXc/PHN/HHCH7kt9TYAHr/0cX6z/Dfnfrcz7u7y8MS11yqX7lVXnf+stq6Wmz+5\nmTvS7uC+kfcB8NRlT3HDghv4xehfECLsvQb1dzx89ZX398Q8fe6QlY/5dhtVly5dGDp0KJ9++im3\n3XYb33zzDZGRkYwcOdLtd+rqlGfG2/knzUgpL/gHkDt2SNm/v9TMK69IecMNDT+79ZNb5S8X/7LB\nZwVlBbLNk23k2eqz2gu3iH//W8q77vLtOx99JOWUKQ0/u3/R/fK2T25r8FlZZZls82QbWVJREmAr\njWfJEikvu8y376xaJeWQIQ0/m7dynpzx3gxZV1d37rPauloZ//d4eaDkgA4tNZbFi5v+v/XG7t1S\nJidL6fRPlvM3z5fj/jtO1tTWnPusrq5ODnhxgFx7ZK1OrTWOzEwp+/Xz7Tv5+VK2bSulmiLtyUsv\nvSSnTZsmpZRyzpw58o9//KPbdwH5+uvvymHD/tnkc+nnvGvvJYMPVFb6FkVw1VUqn091/V7w9znf\n892h7/jzpD83eC8+Kp7esb3ZmLtRx9Yagz8rrMsvVzm/HCGa+47v470d7/Hs1GcbvBcVHsWErhNY\ndsCFT8Nm+GKFOhgzBg4dUrdRgnLTPLfhOZ6f9nyD6z1DRAhX9LiCJdn2TwblqxsToHdv9R3HVRGn\nK0/z2MrH+Ne0fxEacn7JKoRgSo8prDy0Ur8GG4Q/cpGYqPrCzsycOZOVK1eSm5vLZ599xs033+zx\n/ZoaHa0OmtGeh7dEeI1JSIDu3c/7uR9b+RiPTniUqPCmGycTukxg1WH7H0P2dcMcoE0btfexfLn6\n/bGVj/HL0b+kfUTTg/zjOo9j/dH1OrTUWPyZLFq0gClTzvu5/7zqz9w++HZS2qU0effSbpdeEOPB\nn34QQuWAc9xt8cy6Z7iixxUMTRra5N30lHS+PfStDi01Fl/3PBzMmKF/W/QkLi6O9PR05syZQ/fu\n3enTp4/H92trg8rDJf6sNqdMgaVLYd3Rdew9vpc5aXNcvje+y3jW5awLvJEG489KE877+3cX72bF\nwRU8ONp1xMaYTmNYn9M8lQeo/Z9Fi5TV8e6Od/ndJb9z+d6QxCFsO7YtwFYaj7+T5vTp9bmuqsp4\nafNL/P6S37t8b2znsWzM3ehwHdsWf8fDlVfq3xa9ufnmm1mxYoVXqwOClodbfHVbgXLZLFsG//z+\nnzw85mHCQ11LWlpiGtuPbdehlcYSiJAsXAj/2vA89wy7h+hw1ycth3ccTkZBhu3j/P3th6lT1Z0O\nL2x4mRsH3EhcVJzL9/rH9efQyUNUVFcE2FJj8bcfJk2CLVvgxXWvMaHrBHrF9nL5XofIDrQOb82h\nk4cCa6jB+NsPQ4bo3xa9ufXWW6mrq+Ohhx7y+m5trb7XDTcb5eGr2wpg3DjIOJTD8uwV3J56u9v3\nUtqlUFpVyvGK4wG20lj8FZKePaFF9Eney3if/xv+f27fiwqPomu7ruwutvf9rf72Q3w89Ohzlhc3\nvMwvRv/C7XthoWH0i+tHRkFGAK00Hn/7ITISxoyt45/rnuPhMQ97fDc1MZXtBfZeWPnbD83tXveg\n5eEGf9xWrVpB4pXzGRt9K61bur8IQwjBoPhBtp8s/NnzAOXnTpr2Br3ENJJaJ3l8NzUh1fZWmD9W\nqIOOkz8lpiqVvh36enwvLSGNH/J/8K8Sk/DXbQXQffK31FS0ZnSn0R7fuxDGQyD90JwI7nm4wZ8J\no7auluIurxG95x6v76YmXBgrLH+F5FjiG0Rk/dTre2mJabb39/tjhTrIjXuN8J13en1vSJL99z38\nXXEDHIp5lRYZP20QaeaK1IRUMgrtvagKpB+aEzU1QbeVS/yZML499C3JbRPZvGiA13ebs3m+7dg2\n6sJPsnPhRK93eqclprGtoHlOmodPHuZw1Q8cWXoNZ854fndI4hB+OGZvy8Pffig5U8L3x7/mzIZb\nvGaXTU20v+Xhbz9caLnMvBF0W7nBH7fVW9vf4u5Rt3PypPcLoi4E89xfIXlz25vMGXIb8XEhZHhZ\nRDr6wc4RNv72w9sZb3PTwBtJHdCKdV6C6wYnDCazMNPW+Z38ddd8kvUJl/e4nAkjY7xelNWrfS+O\nlR2jtLLUv0aagL/jYeG+hfo3xkKCbis3+Oq2qqiu4Is9XzB74E2MHYvXyWJg/EB2F++2dYZZf/Y8\n6mQdH+78kFsH38r48d5TkydGJxIaEkpuqX2vY/R3svho10fcOPBGJk0CN5dRnqN1y9YktU5if8l+\nv9poBoH0w6z+s5g4Eb77zvO7oSGh9I/rz47CHf410gT8VaLvZzavi3+Cbis3+Oq2Wpq9lOEdh5MQ\nncC4cd6VR1R4FJ3adGLv8b2BNdRA/JksNuRsoH1Ee/p06MOYMbDeyzEOIYTt9z386YfskmyOlR1j\nXOdxjB2rEmd6Y3DCYFsHUfjTDyVnSlifs57pvaZrUh5Qv+/RzPrhTPUZlmWrbApCiAv+Af1DdZtV\nYsS2bbW//9nuz7i277WAugTnYc8RicD5yWJAvPc9EivwZ8P8092fcl3f6wAYPRr+8hfv33G4rmb0\ntucRXH+CJz7O+phr+15LaEgoI0fCxo0qkZwnYRscr8bDrAGz3L9kIf6suD/f/TmXdb+M6PBo0tJU\nxmVvl6w1RyW6/MBy0hLTWClXIiWMGfMsd9wRT2v3QZnsPb6XTbmbuGWw58y2VvHDD/q6rbwqDyHE\nDcBiKeVpIcQfgKHAn6WUW/VrRuD4MmHU1NXw1d6vzuWxGjFC3ajnLYW1Y4U1e9Bs9y9ZiK9CIqXk\nk6xPWDBLXerRty+UlKh7DOLj3X8vLTGNT3d/GmBrjcOf4ImPsz7m8UmPAxAXp57du6F/f/ffGZww\nmLcy3gqgpcbiz6S5YNeCc2eewsJg8GB1t8WECe6/MzhhMB/s/CCAlhqLP/3wye5Pzi2qhIC4uCiy\nsgpJSXH/nbCqME4cP0GhTS8CKS2FFi003luhAS2Wxx+klP8TQowHJgNPAy8DtspR7suEsebIGlLa\npdC5bWdApa0eOBA2bVKp2t0xOGEw87fM16G1xuCrkGQWZlJTV0NaYhqgVtmjRinX1TXXuP9eWmIa\nj618LMDWGoev/XD01FGyS7JJT0k/99moUSrvmTflYecVt697YCfOnGDt0bV8OPPDc585rDBPymNQ\nwiB2FOygTtbZMj17VZX2u35ALS6/3PMlcyfOPffZ2LF3UVQEv3B/dhQpJf/4+z+46ac3kRid6H+D\nDeLZZ+HgQf3K0/J/2hFOMgN4RUr5FaDLkRshxGtCiAIhhNvdNiHEv4QQ+4QQ24UQbhMG+GJ5OLus\nHIwYocw6TzS3yeLT3Z9yfb/rG8Tya9n36B3bm9zTubaNsPFVeXy6+1Ou6nMVYaFh5z4bPdr7vkf3\nmO4Ulhdy6uwpP1tqLL66MZdmL+WSLpc0ODDrUB6eaB/Rnrat2nL45GE/W2osvo6H1YdX07VdV7q2\n63ruszFjvI8Hx36gXaMyAzn/5AotyiNXCPEf4EZgoRCilcbvaeF1YKq7PwohpgM9pZS9gLtRFo9L\ntHaMlNKl8khLg21e9oDtnqbEVyH5fM/nTfph2DDv/dAipAUD4gfYNsLG135YvH8xV/ZqmAXPYXl4\nIjQklAFxA8gszPSjlcbjaz8syV7C1J4NxVGL8gC1sLLrOSh/5MLhsnKQlgYZGe6vKnZg58PEgWRe\ncIUWJXADsBiYIqU8CcQAv9ajcinlauCEh1euBt6sf3cD0E4IkeDqRa3nPHYV7UIIwYC4hpveWpSH\n3dOU+LLSLCovYn/JfsZ0GtPg88GDYbuGsZ+WYN+IK1+EpLKmkjVH1jC52+QGnw8aBHv2KGvOE3a2\nRn2xRKWULMlewhU9rmjwec+ecPo03g8L2jjiylflsTR7aRMl2q4dxMa6vrLZGTtHIup90t6r8pBS\nlgNFwPj6j2oAs4Lbk4GjTr/nAJ1cvah1wlh+YDmXd7+8SdqFAQNg794Le7LwZXCsOLiC9JT0Bq4a\ngM6dlSL2tudn55PFvvTD2qNrGRA/gJiImAafR0RA165qTHjC7uNB62IiszCTVi1a0bN9zwafC6Hu\nNd/hxci0cz/4EnWWczqHwvJChiQ29ZCnpnpfWAWVhxNCiLnA/wN+W/9ROPC2fk3wSuPkOi6PNu/a\nNZcPP5zL3LlzWenhhNeKgyu4rPtlTT6PiFCXQ+3a5bkxdl5h+bLSXJa9jMu7X97kcyGU9eFtsrBz\nmhJffLtLs5cypbvre9kHDjx/o547BicMtm1uJ18mi8X7F3NFjytc5rLq319FI3rCzsrDl35YfmA5\nl3a7tMGtiQ60KI++Hfpy6OQhzlR7yW9jAZWVcOTISubOnXvuCQQtbqvrgGuAcgApZS7gIdpZV3KB\nzk6/d6r/rAkdO87l7rtVh6Snp7ssrLq2mlWHV3Fpt0td/l2L66o5TBZSSpYdcK08QJvrys7pOXyZ\nLJZmL2VKD9fKY9Ag70p0UPz5SCO74ctiwtV+h4MBA7wrj96xvck5nUN5VbmPrTQeX8aDN7nwlr4n\nPDScPh362HIfrLISBg5MN1V5VEp5XjKEEPoFCnvnC+D2+npHAyellC69r1rcVpvyNtE9pjsdIju4\n/HtamveIq4HxA9lVtIuauhqvjTcbrUKy5/geQAm8K7QISZuWbUiKTmJfyT5fm2k4WvuhsLyQAycO\nMDJ5pMu/a7E8YiNjadOyjS0jjbS6rcqrytmQu4FJKZNc/l2L8mgR0oK+HfradtLUMh7qZJ1ya/dw\nrTy0WB6gvBN2dF2Z7rYCFggh5qM2q+8GVgCv6lG5EOJ9YB3QRwhxVAjxEyHEPUKIewCklIuAA0KI\n/cB84F53ZWlxVSw/sNyly8rBwIHe3VatW7YmKdqeOY20ThYOl5W7dNupqd6VB6h9jwtZSJZlL2NS\nt0lN9n0cDBrkXXmAssLsGHmmtR/WHFnDkMQhbu+06d9fyYW3XJipifZ06Wrd88gszKRNyzYu760H\n6NEDiorglJfIbLvue5geqiul/Dvwcf3TG3Vo8F96VC6lnC2l7CilDJdSdpZSvialnC+lnO/0zn1S\nyp5SylRPp9q1RFt5Ux59+qgIG2/Y0b9bW6vCCFtoOPa57MAyt6srUCvNrCzvYYl2jbjSGjzx7aFv\nm0RZOdOjBxw7ptJzeMKO4wG0T5qrDq9iYlf3p2Pj4tRp8/x8z+UMjrdnuK4viwl3LitQqT369FGZ\nBzyRlphm234wO1QXKeVSKeXD9c8y/arXD28dU1ZVxg/HfmB8l/Fu3+nSBYqLtU0Wdos0cvi3vdzd\nQ52sY82RNQ1OUzcmKkqFJR496vYVwN5ComWyWH1kNZd0ucTt3x2ThZYgiq35tsrWA/jYD13d9wNo\nc13ZVYlq7YeVh1e6dd056N3bewSeI6jGbvtgprmthBBlQohSN89p/ZqgD95Msg05G0hNSCUyzH3y\nqtBQ6NVL2+CwW6SR1oGxq2gXsZGxXtMnaOoHm7qttJjnheWFFJQVMDB+oMf3eveG/V48lMM7DmdL\n/hYfW2k8WsbE2ZqzbM3f2uS8T2O0RFw5Lky7ECfNOlnHuqPrPC4uQS0mvMlFTEQMMRExHDjh5VCI\nyZimPKSU0VLK1m6eNvo1QR+8ua3WHl3LuM7jvJajZWUxvONwNudtttWFSL74t7X2wz4ve+Gd23Sm\nsqaSgjIvJ8hMpK4Oqqu9u2vWHFnD2M5jXYZkOtOjB2Rney6rR/senDp7iqLyIh9bayxa3FabcjfR\nL66f2/0OB336eFeiHSI70D6iPfuO2yuIQks/7C7eTbtW7UhqneTxPS3zA2DLNCVWnPPo4urRrwn6\n4M1ttfboWsZ2Huu1nG7dvCcP69SmEwJBzukcH1tpHFo3y9ceXet1dQXahORcLh8bua4cE4U3993q\nw55dVg569vQ+aYaIEIYmDbWd9aFlslh1eBUTunjIelhP9+7eT1fD+YWVndDSD2uPaFtc9urlfVEF\n9oy4siLaahGwsP5ZARwAvtavCfpQU6M29VxRW1fL9znfa1Ie3bt7Vx5CCNsJidaY/jVH1mhSHlqU\nKNhPSLRuCq45usarnx+U5eFNeYD9Js26Os8y4UDLfgf4oDyShrMpb5PGVpqDlklzzVHtcnHokPc6\n7XiI1or0JAOllIPqn17ASEDDPWvm4mmzOLMwk8ToROKi4ryW062b9hWWnYREy8BwZMLtE9vHa3kp\nKXBYw9EFu4UlaumHsqoysoqyGN5xuNfytLitwH77HlossJq6GtbnrNc0aaakqEnTWwTeiOQRtlKi\noFF5aHTnxsaq8k572fW1o9vKiqy6DagPl7XVXR6gz34HaLM8AEZ0tJeQaDLNj65lXJdxbs93ONO1\nq/YVlp3cVlr6Yf3R9QxJGkKrFt5NlMREOHFCleuJYUnDLrjxkFGQQac2ndwemnUmMhJiYiDXy9X1\nQ5OGsu3YNlsdovW255Ffms+JMyfoF9fPa1lCaFtYpbRL4eTZk7bKwG16qK4Q4iGn59f1B/u8DCHz\n8SQovijKt7rhAAAgAElEQVSPLl1UiGqtl6wbdts017q6Gt/Z+yoToH171QcnT3p+r19cPw6cOGCb\nXD5aVldrj67V3A+hodCxI+R42d7qHtOdsqoy2wQPaBkPG3I2MDp5tOYyU1K8h2+3a9WOjq07klWU\npblco/HWF45FldaLrBxWmCdCRMi56DO7YMWeR2sguv4JB75C5bqyFR4tjyNqcGihZUu1wvKWVTYh\nOoHo8GjbhONpiShxCIkWhFDWh7cVVnhoOH1i7ZPLR4uAbMzdyKhO2o1nx4LCE459MLu4rrSMh415\nG92mZnFFcrJ3JQr2cl1pOTyrdbPcgRblAfY7RGvFnsdcKeW8+ucvUsp3pZRn9WuCPrhTHvml+Zyu\nPO02j5MrkpO9m+dgr30PbwPjbM1ZsoqyXKaadkeXLnDkiPf3hiUNu2D6QUrJprxNjOg4QnOZnTtr\n64fhSfbZNNcyUWzK3cSIZO390KmTNuVhx37w5KndlLeJUcnaFxOdOmmfHzbmarhJyySsCNVdJoRo\n5/R7eyHEEv2aoA/uOmVz3mZGJI/w6W5lzSssG+17eBsYGQUZ9I7tTURYhOYyk5K8p6QAGNN5DOtz\nvNxdaxLe/LqHTx0mLCSM5DbJmsvUqkRHJo9kQ66X6wdNwtt4KK0s5eDJgwyKH6S5TF8mTbssJrxZ\nYDV1NWw7to2hSUM1l9mxI+TleX/PTnIhpe/XVHtDy4waV3+DYH0jZAng8jY/K3E3YWzO28zwJO9R\nNc5oFZILabLYnLfZp9U2KCHRpDw6jWH9UXsIibd+8HW1DdrHw5jOqh/scMLa26S5JX8LqQmpbpNC\nukLrompo0lB2Fu20xT6Yt/Gwu3g3yW2SaduqreYytSqPXu17UV5VTu5p67eIq6vV/l2IXheIo015\n1Aohzt0EL4RIAayXjka4GyBb8rcwrOMwn8rS6rYamTySH/J/oLLGSyiOCWhRHlpCU53Rann0i+tH\ncUUxheVeNopMwKvy8NFlBZCQ4P0aVoDE6ERiImLYU6whu6bBaFGivux3gHa3VVR4FP069LPF/o9R\ncqFFeQghGN1ptC2sD71dVqBNefweWC2EeEcI8Q6wCvidvs0IHFeWh5TSr8GhdYXVumVresf2tkVS\nPG8rzU15m3zuB60rrBARooTEBtaHt2grf5RHYqLKrquFsZ3Hsu7oOp/KNwJvLoqNeRv9skS1jAeA\ncZ3HsfbIWp/KNwItStRXz4Qv/TC281jbyIWeYbqgbcN8MTAM+KD+GVr/ma1w1TF5pXnUyTo6t+nc\n9I8eSEhQefu1MK7zONYetbeQlFeVk12SzaAE7f5t0G55gHJd2WHS9NQPdbKOLXlbfFaiWi0PgLGd\n7KE8vKWr2ZjrW6QVQHy8D3LRxR5y4W1RtTl/s89uzJgYNRlXVHh/d0wne+x7WGJ5CCFCgKkopfEV\nECmE8G3UmYCrjtmct5lhHYdpOhTnTFyc91BdB2M7j7WFkHgaHNuObWNg/EDCQzUkv3LCJ+Vhk81B\nTxvme4r30CGyA7GRsT6V6bPlkWMP5eFuPBSUFVBaWUrP9j19KjM6WoW+apk0HRaY1eegPPVDVW0V\nmYWZpCWm+VSmENplY0TyCLYXbLfctW2V2+olYAwwu/73svrPbIWrCWNL/hafTVJQysOXFZbdhcQf\nlxWolWZhofcb5ABGJY9ia/5Wqmurfa5HT7z1g6+rTFCTJni/5wXUNcV5pXmWnyz25Lbakr+FoUlD\nfV5UCaFdNjq16URUeBR7j2tIQWsgnsbDzsKddGvXjejwaJ/LjYtTd/94Izo8mt6xvfnhmJf7rQ3G\nKuUxSkp5L3AWzkVbaQ/RMAlXysOf/Q44b3lomTS7tO1CeGg42Sc0JEAyEE/muT+RVqD6NDxc26TZ\ntlVbusV0s/xQlCch2ZLn32ICtFsfoSGhjEweyfc51qZ/8+S22n5su8+rbQeOBYUW7ODS9TQe/J0f\nQOW40qI8wB4uXauUR5UQ4tylB0KIOC6AaCsppV+RVqBu0gsN1TZpgj02Bz1tiDlWmv7QoYN2IRnf\neTyrj6z2qx698CQk2wu2k5qY6le58fG+7XvYedLcVrDNb+Xhi0t3XOdxlk+a3hZVw5J8nx9AycVx\njcblJV0uYdXhVX7VoxdWKY/ngU+BeCHEE8Ba4Ek9KhdCTBVC7BZC7BNCPOLi7+lCiFNCiB/qn0fd\nldV44swtVbG2ya21HwZzxifXlU1WWK6Ux5nqMxw6eUhT0jdX+KI80lPS+fbQt37Voxfuoq2klEp5\nJPinPGJjoaRE27uXdL2E7w5/51c9euFp0tx+zP9+8HXTfM2RNX7VoxeeJs2Mwgy/laivcrHq8CpL\nz/9YojyklO8Aj6AURh5wjZTyf4FWXG/NvIDajO8PzBZCuJrhvpNSDql/HndXXkSjg9M7CnYwOGGw\nz35dB76Y5xO6TmDloZV+1aMX7iyPXUW76B3b2+fNcgexsdpXWOkp6aw+vJraOi9ZJQ3EnRLNOZ1D\ny9CWJET7d761fXuVXVcLYzuPZfux7ZRVaTRdDcDdZFFeVc6RU0fo26GvX+X6IheD4gdRWF5IfqnG\nqAsDOHu26dwAKvIuszCTwQmD/SrXF+WR1DqJ+Kh4S1O0652OHTzfYd7e8QAFwPv1T0H9Z4EyEtgv\npTwkpaxGhQG7SrioafaPbHQ1+Y7CHT6lXmiML5bHoIRBnDh7gqOnvGTPMxB3yiOjIMNvAQHfhCQh\nOoHkNsmWbg66mzQDcVmBCs/UanlEhkUyNGmopa5Md/2QWZhJ3w59fTpZ7owv4yE0JJSJKRMttUbd\nycXBEweJjYj16WS5M770A6iFlZULTL3TsYNny2MrsMXp2Vz/OH4OlGTAebbNqf/MGQmMFUJsF0Is\nEkL0d1dYE8sjQOURE+M9HbmDEBHCpJRJfHPwG7/rCxR3K4uMggwGx/uvPHyxPADSu1ovJC6Vx7Ht\nAfVD+/balQfApd0utXQ8uHNbBapE27WDU6e0vz+522TL5cLqRRXApJRJlipRI9xWbhMVSylT9K2q\naRUa3tkKdJZSVgghpgGfAS7T4y5fPvfcZJ+ens6Ogh08OOpBvxvXrp125QH1k8Whb/hx2o/9rjMQ\n3AnJ9oLtTO051e9yfRaSbpN4fdvrPDz2Yb/rDAR3QpJRmMFVva/yu9z27b3f6e7MpJRJ/HrZr/2u\nL1Dc9cO2Y9tIS/DPzw/+ycUz65/xu75Acbeo2l6w3VTlMTFlIj9b+DNq62oJDQn1/gWdcYyHlStX\nsnLlSl3K9JDl/jxCiGuACagJ/zsp5Zc61J0LOB/97oyyPs4hpSx1+vlrIcRLQoj29eHCDbjhhrnc\nfrv6ubq2mr1r99I/zq2h4hVfhWRyt8k8vupxpJR+77MEgivlIaUMeIUVGws7d2p/f2LXidz5xZ3U\n1NXQIkTT8NIVT5bHo5e4jbfwiq+Wx+hOo8kqzuLU2VN+u0YCobISWrdu+vn2gu3cMOAGv8tt1077\n3g9Avw79OFtzloMnDtItppvf9fqLJ8vjpoE3+V2uL25MUHnPkqKT2HZsm18RoIHikIv09HTS09PP\nfT5v3jy/y9RywvyvwAPATiALeEAIoUe01WaglxAiRQgRDtwIfNGo7gRRPxPXn2oXrhQHNHRb7SvZ\nR6c2nYgMi3T1qiZ8VR492/dECMG+kn1+1xkIrnya+WX5CCFIjE70u1xf+yEuKo4ubbuwJc+apHiu\nVpoV1RUBbRKD78qjZYuWjEoeZVmIpiu3VZ2sY0fBDr8jrcD38SCEsNSFZ5Tbqm1b7/eYN8bKaESr\nQnWvBKZIKV+TUv4XFR01I9CKpZQ1wH3AEmAX8KGUMksIcY8Q4p7612YCO4QQ24BnAbdLBecN8x0F\nO3zO49QYX327Qggmd5vMigMrAqrXX1xNmhkFGaQmpAZkCfkjJJd1u4xlB5b5XWcguFKiOwt30qdD\nH783icH3lSZYu+/harI4cOIAMRExxETE+F2ur8oD4NIU5dK1AlfjoayqjPyyfJ/TszjTpo1v8wPA\n5d0vZ0m2NVchWaU8JNDO6fd2aNuv8F6wlF9LKftIKXtKKZ+s/2y+lHJ+/c8vSikHSinTpJRjpZRu\nj+06Wx6BbpaDn0LS7VKWH1weUL3+4mqFFejqCpSQ+Ko8pvWaxtf7vw6oXn9xJSSB+rfBt1BdB1N6\nTGFxtjU5RF2Nh8zCTEvkYnJ3taiy4pyDu37o16FfQG7VNm3UIeI6H/5Jk7tP5vuc7y0J4TY1VNeJ\nJ4GtQog3hRBvoqKtntC3GYGjt/Jo29Z3IZnSYworDqygqrYqoLr9wZWQ6DFp+qM8JnSdwI6CHZSc\n8XGprgOulMeOAmvGw9CkoZw4c8KSe+7PnGkagbiraBcD4gYEVK5DefiSyi2lXQoxETGWXF1g1KIq\nNFT1r9YsFABtWrZheMfhlkQjmh2qC4CU8n1UYsRPgI+B0VLKD/RtRuA0UB46ua18nSwSoxPpFdvL\nklO1roREj0nTH/O8VYtWTOg6gaXZSwOq2x9cKY+s4qyAJ83WraG01Pt7zoSIEGWF7TPfCnOnPAIJ\nIgHVty1aqPJ94cpeV7Jw78KA6vYHd8ojULkA/1y603paMx5MdVsJIYY6HiARFQmVC3Ss/8xWOPY8\nSitLKSgvoEdMj4DK80d5AMzoNcMWQlJTV8O+kn0BbRKDf5YHwPRe0y1xXbkSkl1Fu/xOz+IgIkJd\n5VlT49v3pveczqL9iwKq2x/OnHGx91O0M2DlAb5HXEG98thnvVxAvQUWH9hiAvxbWE3rqVy6Zmfh\nNnvPYzPwBvBM/fN0/eP43VY4VlkOf2agsdS+bpg7uLK3NULSeHAcPHGQhKgEosKjAirXoTx8HevT\nek5j8f7Fpvu5G08WpytPc+LsCbq07RJQuUKo1Oy+uCkALu9xOasPrzb9Pu/GaTlq62rZU7wnYCUK\n/i2sxncZz76SfRSUacwuqROulEdWcRb9OgTeD/5YHgPjB1JVW2V6qnpXyuPUWT8mOCc8KY9fAaVA\nBfA6cLWUcpLjCahWA3AIyo7CwF1W4J+PG5Sf++TZk+wv2R9wG3yhsZBkFWfpMlGEh/vnpugW042Y\nVjH8kG9uqpLGQrK7eDd9O/QlRGjZ3vNMdLTvrqt2rdoxJGmI6SGajd1Wh04eIj4q3q+7KxrjjxIN\nCw3jsu6XmW6NNpaLk2dPUlZVRqc2nQIu2x/LQwjBtJ7TWLTPXGu0sVxU1VaR8LR/ed4cuJUoKeWz\nUspxqDMenYAVQogFQgj/j6cayDnloYOfH9SkKYSKl/eFEBFiiX+3ifIo0md1Bf67rmb0nsHnez7X\npQ1aaSwku4p26dYP/ux7gDX+/sbKQ4/9Dgf+KA9Q/fDV3q90aYNWXMlF3w59dTnI64/lAXB1n6v5\nbM9nAdfvC43lIrskO2BrXMuGeTbwObAUGAH0CahGg2hgeeigPEDd6+GPkMzoPYMv9n7h/UWdqKtT\n/njnQ2F6mebgv/KY2X8mH+36SJc2aKXxpJlVlKXbpOmv8ri277V8uvtTU114jftBr/0OUHJRXu77\n96b3ms7yA8s5W3NWl3ZowaVFrqNc+OPavrzH5Ww/tt1UF17jUN2s4qyA90M9bZj3EEL8XgixEZgH\nbAf6SSk/DKhGgwgJUek49HJbgf9CckXPK9ict5mico1peQPEcWuc82JKL7cV+L/CGpk8ktOVp9lV\ntEuXdmihoqLhgdFdxdZbHr1jexMXFWdqlt3GG+Z6hOk68NfyiI+KZ0jSEJbsN++gXONJU09LtG1b\n/5RHqxatmNZrGp/tNs/6OHu2oVzo4ZnwZHnsA24AvgbWA12A/xNCPCSE+FVAtRpEflk+ISKEhKjA\nfHkO/FUekWGRTOs5jU+yPtGlHd5ovLqSUrK7eLflK6wQEcKP+v2Ij3d9rEs7tFBR0dTy0EuJ+jtp\nAszsZ64V1njDXE+3lb9yAXBD/xv4366ArwPSjFF7gaDGg7/98KN+P+LjLHPlooHy0KEfPCmPP6HO\ndtQB0fVPa6f/2g7HfodeiQkDEZIbB9xompA0PgCUV5pHqxatiI2M1aX8qCg1+PzhR/3NE5LqamWB\nhtVnITlTfYbc0tyAw7Yd+Gt5QL0LL+sj01xXzm6rOlmnFhM6TZqByMX1/a5n4d6FpkWfGbkXGBnp\nv1xM6zmNDbkbTDtI21h56LG49LRhPldKOc/dE1CtBqHX4R8Hgawspvacypa8Lab4NY3060JgQjKu\n8ziOlR0zJfqssdWx9/heusd0DyinlTOBKI9+cf2IaRXD+qPrdWmLJ6RsqDwOnzxM+4j2tGnZRpfy\nA7HAEqITGJo01LQcT86ycab6DHmlefRor89iIpBFVVR4FJO7TeaLPebsjTorD8diwrA9jwuRH479\nwNAk/c4v+rthDhARFsGVva80xXXVZDNMx9UVBKY8QkNCmdl/Ju/veF+39rjDpWmuYz8EojxAWR//\n22m8NeqwwFrUp27S45CkM4FYHgA3DLjBlH6Ahlb5nuN76NG+h25XBQQiF6D64b0d7+nSFm84L6xy\nTufQpmWbgK8KCCoPDwQsJP1v4L1M4weHEf5MZyIiAhOS21Nv583tbxp+qrZJP+gYaQX+nfNw5uZB\nN/PBzg+orq3WrU2uaLxZrrcSDcQiB+W6WrRvEaWVAXSmRpwXVnZaVAFc0+caNudtNuX6amfZ0Gsf\n0FO01YP1/x0fcC0mUFZVxpFTR2y1wprWaxp7j+81/DRpeblqqwM7ua0ARnQcQXhoOGuPGhttZGSk\nFSjLw19LFFTUVe/Y3oZnIGi8Wa63Eg3EIgcVdTUxZaLh1oeUjSZNA+QikPkhIiyCWf1n8U7GO7q1\nyR2N+6FvbGAuK/Bsefyk/r/PB1yLCWw/tp2B8QN1vb0uUOURHhrObYNv4/UfXtetTa5oojx0jDCC\nwJWHEII5aXN4Y9sburXJFUae8QBVtq8n7RtzR9odvL7N2PHQpB9sZnkA3DnkTl7b9po+DXLD2bMq\neMLhvtPbIg9ULgAlF9vfMNUq1yt4wpPy2CWE2Af0EULsaPRkBFyzzmzN38rQRH3zNQaqPEAJyZvb\n36SmzseMej7grDxOnDlBRXUFya2TdStfDyG5dfCtfJz1MRXVARbkAWcBqamrIftENr1jXV557xd6\nKI9Z/Wfx3aHvDA2kcFYeUkrb7XmAijY6cOIAWUVZ+jTKBS4XVTayyEFdVwywIXeDDi1yTU2N2gc7\n577TaTHhKdpqNnAJsB91c+BVTs/VAdesM1uPbWVI0hBdy9RDSPrF9SOlXYqhaZjLyxuZpDqlX3Cg\nh5B0bN2R0Z1GG3rWwVl5ZJdk07F1RyLCIjx/yQf0UB6tW7bmun7X8XbG2/o0ygXOyuNY2THCQsPo\nENlBt/IDibZyEBYaxu2Db+e1H4yzPpyVh2Mx0aeDfgky9JALIQRzUufw6tZX9WmUC86cUW11TAmG\n73kASCmPSSkHA/mcP+uRK6U8FHDNOrMhZwPDOw7XtcxAfbsO7hxyJ69sfSXwgtzgLCR6u2pADbxA\nJ02An4/4Of/a8C/DTHTniBI9T1Q7aNVKuUIC5a6hd/Hvzf+mtq428MJc4LxhrrfLCvRZVAH8ZMhP\neCvjLcPSlTjLxYETB0iMTiQyLNLzl3wgkFBdZ+4ceicfZ33M8YrjgRfmAudFVXFFMZW1lSRFJwVc\nrtdoKyFEOuq0+Uv1zz4hxMSAa9aZgvICUhNSdS1TD98uwOxBs/k+53v2Hd8XeGEucBaSnUU7dZ8s\n9FhhgcptdPLsSdYdXRd4YS5wFhI9czk50MPyABjTaQztI9obtnHubHno7aoB/RZVfTr0IS0xzbAw\nbme50POEvQO95CI+Kp5r+15r2ALTYXkA7CzcyYC4Abp4JrSE6v4DmCKlnCClnABMAf4ZcM2AEGKq\nEGK3EGKfEOIRN+/8q/7v24UQbv1Sl3a7NOA7PBqj1worMiySe4bdwz+/16XbmlBRcV5IMgszdcvt\n5UAvIQkRIdw/8n6e2/Bc4IW5oIGQFO3U3fLQS3kIIfjl6F8aNh7KylRkGCjLw4hJUw8LDOChMQ/x\nj+//YYg1WlZ2Xi52FOxgYNxAXcsPNNrKmQdGPsCLm140JIy78aJKL7nQojxaSCn3OH6RUu4FAg5p\nEkKEAi8AU4H+wGwhRL9G70wHekopewF3Ay+7K+/q3vpvw+g1WQD8fOTPeT/zfUNMU+cVlp5ZhR3o\npTwA7hhyB8sPLDcktr2BkBTu1OW2OGf0HA8z+89k3/F9bDu2TZ8CnSgtPa889N4sB7XxqpfyuLz7\n5QgEyw4s06dAJ5rIhU0XVQBDkobQrV03Pt39qT4FOtFYLgbG66NEtSiPLUKIV4UQ6UKISUKIV1G3\nDAbKSGC/lPKQlLIa+AC4ptE7VwNvAkgpNwDthBAusx5e2/daHZrUkFat1AlVPUiMTuS6vtfx8ma3\n+s9vHEJSXFHMmeozulx040yghwSdadOyDXek3cHT657Wp0AnHHseel3B25iICP0mzbDQMO4beZ8h\n/VBWplyuYMyeh157P6CssF+N+RXPrNf/ctLy8vP9kFmYqfuiqmVL/64mdsdDYx7iyTVP6m6FOSuP\nzKJM3RZVWpTH/wFZqEuh7gd21n8WKMmA8/Izp/4zb++4nBlbt9Q/V6OeKyyAh8c+zPMbn+d0pR/5\nzT3gUB6ZhZkMjB+oa6QV6LvCAvj1uF/zdsbb5JzO0a9QzgvJ/pL9JLdO1nVzFNSkqZflAXDviHtZ\nmr1U93BVh/LQ89Y8Z/RUHgCzB84mqyhL97xfDrmorKnk4MmDui8mhNAvmATUJVECoXuqdodcSCnP\n7XnogVf3k5TyLMbcW65VvTaeCV1+b+7cued+Tk9PJz093a9GOaOn5QHQP64/U3tO5Zl1zzBvkn65\nJR2hunrdotgYvZVHYnQiPx36U55c/SQvXvmibuWWl0NsrDEuK9DXbQXKCnt47MM8tvIx/jdLv9PW\nDrfVrqJduodtg1pUVVaqE9x6FN2yRUsem/gYv/vmd3xz+ze6tdehPLKKs+ge052WLVp6/5KPOCKu\nWuuwdhVCMC99Ho9++yjX9L1Gl6uT4bzyKCwvpHJ/Jf9++t+mbZgbRS7Q2en3zijLwtM7neo/a8Lc\nuXPPPXooDtDf8gCYlz6PFza9oOtFUY0tD73RW3kA/Hrsr/lg5wccPnlYtzJPn1YX9BixWQ76Kw9Q\n4ctrjqzR9a53h+WxNX8rQxL1PfsE59PeV+u4t/vjtB+TV5qn696HQy6MWlSBvpvmoG4hDQ8N1/UO\nnLIy1c6dRTtJG5PGvHnzzs2VgWCl8tgM9BJCpAghwoEbgcb5ib8AbgcQQowGTkopTbu7UW/zHCCl\nXQq3DLqFx1c9rluZ54SkcIdhykPvSTMuKo77RtzHI8tdBtn5RWmpurhqU94mhiUN061cB3rueTiI\nCo/it+N/yyPLH9HN1+2wPLbkb9H97JMDvRdWLUJa8OdJf+a3K36r250njmgrI4JIHOg9RwgheHLy\nkzyy/BHd7jxxLKq25G0hLSFNlzLBi/IQQoQKIfTf0QOklDXAfcASYBfwoZQySwhxjxDinvp3FgEH\nhBD7gfnAvUa0xR16u60cPDrhUd7PfF+3SJvycgiLqCSjIEPXrMIO9Pb1O3hk/CNsyN3AigMrdCnv\n9Glo3VqyIWcDozqN0qVMZ8LD1eZorc5n+342/Gfkl+WzYNcCXcpzWB6b8zYbokTBmIXVzP4ziWgR\nwfzN83Upz7Go2pi7kWEdjemH8HD954jLul/GsI7DeGrtU7qU51Ae3+d+z5jOY3QpE7yfMK8Fxgu9\nnabny/9aStlHStlTSvlk/WfzpZTznd65r/7vqVLKrUa0wx1GuK1AHQr662V/5a4v79LllHF5ORyu\n2kKfDn0MCxwwQolGhkXy3NTn+Pmin1NZE3gFp09DWdhBwkLDdN8kBuXfN2LSDAsN499X/ptfLvkl\np876cd9vI8rKQLYq4fDJw7qHpzowoh9CRAjzZ8znjyv/SF5pXsDllZdDy6hKNudtZmznsTq0sClG\nycY/pvyDFza+QHZJdsBlnTqlFlXf53x/LpeWHmhxW20DPhdC3CaE+FH9c71uLbAxRgiIgzvS7iA6\nPJoXNr4QcFmnT0Nm6SrGdzYme35YmFpt1xlwg+rVfa6md2xv/rzqzwGXdfo0ZFYsJz0lPfCGucEo\nK2xcl3FM7zldFzdeaSnsr/mWcV3GER4arkPrmmKUbAyIH8A9w+7h/q/vD9iNV1YGRS020zu2t263\nKDamZUuoqtK/3M5tO/PIuEf46Zc/DXiBefo0nInaQ4gIoVu7bjq1UJvyaAWUAJeiEiQ6kiQ2e4xy\nW4Hybc6fMZ/HVz8e8GbpyZPwbd4XTO81XafWNUQIY8xzB69c9Qqv/fBawO6r06dhffEirux1pU4t\na4oR+x4O/j7l7yzNXhpw8siyMthWtoTLul2mU8uaYtSKG5Rbd0/xnoDTdZw6BTtrvmRaz2k6tawp\nRvbDr8b8CoAn1zwZUDmnT8P+kIVM7zld38g7KeUF/6h/hv5UVEjZsqUhRZ/j/R3vy+7PdZclFSV+\nfb+6WsqQdjmy3V/bycqaSp1bd542baQ8ccKw4uXS/Utlx2c6ymOlx/wuo3XcCdnmibayuLxYx5Y1\npGdPKffuNax4uTFno4z7W5zcf3y/32X06lsp2z0RKw+dOKRjyxoyfLiUGzcaVrzcXbRbdvhbB7k5\nd7PfZYwZWyc7PdUzoDK8MX26lF9+aVjxMudUjkz4e4JcdWiV32XMmCFln7+Pll/t+arJ3+rnTr/m\nXS2JEfsIIVYIIXbW/z5YCPGofurLvjhMUiPvablp4E3M6DWDWz65xa+8NqdOQfjoV7hxwI2GuShA\nWR5GmOcOLu9xOXcOuZNrP7zWrzs/6uqgrNt7TOkxhdjIWANaqDAiXNeZEckj+OPEP3Lth9dy4swJ\nv6YHsq4AACAASURBVMoojv2c3jH96dquq86tO4+RLl1QSRNfmv4SMxfMJL80368y8sK/pWWLcEOC\nSBwYaXkAJLdJ5o1r3+CGj25gf8l+v8rIr93B8eqjXNHzCl3bpsVt9QrwO8AxdewAZuvaCpsSEqJu\nITNy0gR4esrThIgQ5nw+x2f/5tGik1SlvswDox4wqHUKo4UE1BmYnu17MmvBLJ/TdJecqoRxf+MX\nox80qHUKo5UHqLMfl3W7jKvev8rne75r62o5NfgJfjnqYYNapzBaeQDMGjCLnw75KVPemeJzTjgp\nJfm9/sI9qb/U/ZCkM2bIxdSeU5k7cS5T35nKkVNHfP5+dvLj3NTj57resgralEekVHmlgHobB/RP\n/WhTzBCSsNAwFsxaQGF5ITMXzPQpvvvx739Du4Krdc+c2hgzhEQIwWtXv0br8NbMeG8GJ8+e1Pzd\nx1c/TnhJGuO6jDOwhcbugzkQQvDMFc8wIG4Al751KYXlhZq/+9z6F6mrjGZW6gwDW2iOXAD87pLf\ncVXvqxj/+nifDpS+n/k+1eEF3DVijnGNwxy5ALhn+D3cO+JeJrw+wad0Nkv2L6G03XruTtV/calF\neRQJIXo6fhFCzERdDnVRYMZkARARFsHCmxcSHR7N6P+OJrMw0+t3Xtj4AusLltP/qCFHcRpglpCE\nhYbx7vXvMiBuAMP/M5yNuRu9fuf9He/zwe7X6bRN/6STjTHafecgRITw7xn/ZnrP6Qz7zzC+Pfit\n1+8s3r+YJ9f+hZjvXiM0xNjzv2aNByEET0x+gp8N+xmjXh3F57s/9/qdtUfW8uDXDyI+e5u2rfVd\nbTfGrH4AtYE+L30eE9+YyDsZ73iNRssszOTHn/2YqKVvktg+yuO7fuFtUwToAawAzgB5wFogxd9N\nFiMeDNowl1LKTp2kPHzYsOKbUFdXJ1/d8qqMfSpW3rfwPrnv+L4m7+ScypG3fXKb7PN8H/nS+9ny\nmmuMb1dampRbthhfjzMf7PhAJj6dKO/47A65s3Bnk78XlxfL+xbeJ5OfSZZvL9suR4wwvk3Tp0u5\ncKHx9Tjz9b6vZed/dJY3LrhRbslr+j/h5JmT8jfLfiPj/x4v31u9VvbubXybbrlFyrffNr4eZ9Yc\nXiN7PNdDTn93ulxzeI2sq6tr8PfSylL5xKonZIe/dZAfbP5atm9vfJvuv1/KZ581vh5ntuRtkakv\np8pJb0ySy7OXN+mHyppK+dz3z8kOf+sg38t4X4aFSXn2rOuyCGDDXEtixGxgshAiCgiRUvrmhL3A\nMcvycCCE4M6hd3J1n6v529q/Mfa/Y0mITmBg/EDCQ8PJLslmV9Eu5qTNYcvdW1jwXhTt2hnfLjNX\nWA5uHHgjU3pM4bkNzzH5rcm0bdmWwQmDiQyL5PCpw2zJ28LNg25m+8+2s2FlLLHG7ZOfwyzLw5mp\nPaeS9fMsnt/4PNd/eD2hIaEMTRpKdHg0Oadz2Ji7kat6X0XGzzLIzkigfXvj22SW28qZcV3GsfPe\nnczfMp87v7iTiuoKhnccTttWbckvzWdj7kbSU9LZdNcmaopTmq1cDE0ayua7N/Pmtjf5xZJfUFRe\nxKhOo4iNiKWooog1R9YwKnkU3835jpSo/oSGqnbqjVflIYToADwGjAekEGI18CcppTEX7toMo06Z\neyMuKo6/T/k7T0x+goyCDLKKs6isqWRO6hxGJI8gOlxdVHDypEo9YDRGnvPwRExEDHPT5/LHiX9k\n+7Ht7D2+l7KqMm5qcxPju4w/1w/Hj2Oa8rCiH6LCo/jN+N/w/8b9P3YX7yajIIPyqnKSWicxutNo\n2kcojbH5BMTEGN8eK5QHqAy8D4x6gPtH3s++EnWZVmllKQnRCQzvOJzE6EQANu+n2SoPULnA7hx6\nJ3cOvZMDJw6w7dg2Ss6UEBsRy/wZ8+nYuiMAR45AXJxBbdDwzgfAd8D1qPToNwMfAsadQLIRVgmJ\ng7DQMIZ1HOY2N8/Jk+YJidkrbmdCRAhDkoYwJMl1ltjiYujQwfh2WGF5OBMiQugf199tgERJCc3W\n8nBGCEHv2N70ju3t8u9mLaqsUh7OdI/pTveY7i7/Vlxs3KJKy65aopTyz1LKg1LKA1LKxwGXt/k1\nR8x2W/mKWZOmHYTEE2ZaHlYqD2+cOGGO8rDKItdKSYk5Fpjd5cLI+UGL8lgqhJgthAipf24ElhrT\nHPthdyEpLIT4eOPruRCExAzlYbUF5g2zJk2rLQ9vFBZCgglL3AtBLoxSHm7dVkKIMs7f2vcL4O36\nn0OAcuAhY5pkL+wuJEVFxvk0nbG7kBw/fnG4rbxRUgLdXXswdCU8XCVgtCtFRcFFFVhkeUgpo6WU\nreufEClli/onREqpf95vm2L3wRG0PBRmWR5WbZhrxaxJMzxc35sE9SYoFwpLLA9nhBCDgRTn96WU\nnxjTJHtxIQhJ0PIIWh4Ojh2DxETj6wkLs3c/BJWHorgYBuh/IzOgLVT3dWAQsBNwvtHholAedhaS\n2loVVRLcKFYr7qDyME95XAiLKrMsMLsrD6PmBy2WxyhgQP1pxIsOOwvJ8eMqHLGFsRkYAHuvsGpq\nlJCYsUEaHg4Vvif9NY2g5aEIWuSK/HxISjKmbC3RVpsAY7Pu2Rg7C4lZqyuwt5AUFCirwwwlamfL\n4+xZpdjMiLay86IKzHVb2XU8AOTmQseOxpStRdxeB9YLIY4BjulDSikH+1upEKI96qBhV+AQcIOU\nskkKVSHEIeA0UAtUSylH+lunv9hZSHJyIDnZnLrsrDzy8sztB7tOFgUFyvoyMAP5Oey8qKqoUGnz\nL/ZzHlIq2bBSefwXuBXIpOGeRyD8BlgmpfybEOKR+t9/4+I9CaRLKUt0qtdn7CwkOTnQubM5dbVs\nqa6ztCNGCkhj7OzjNstlBRfGosrgxMKAvZVHSYk6ahBlQEJd0KY8CqWUX+hc79XAxPqf3wRW4lp5\ngEqJYhl2FpKjR6FTJ3PqsrOQGGmaN8bObiszLVE7L6qOHjV3UWVnuTByPGhRHj8IId4DvuT8bYIy\nwFDdBCllQf3PBbhPdyKB5UKIWmC+lPKVAOr0CzsLSU4OjB5tTl12njTNtjzs2g+HDkFX426ebYDd\nF1VmKQ87jwej3blalEckSmlMafS5R+UhhFgGuDKif+/8i5RSCiHcRXKNk1LmCyHigGVCiN1SytWu\nXpw7d+65n9PT00lPT/fUPM2Eh9v3hPnRozBrljl12XmFdfgwXGZSmk47TxaHD5tzuhzsvagyW3nY\nWYk29kysXLmSlStX6lK+lvs85vhTsJTycnd/E0IUCCESpZTHhBBJgMt7NqWU+fX/LRJCfAqMBLwq\nDz0JC7Ovr99Mt5WdheTAAfMmTTsrj0OHYNIkc+qy83g4ehSGuE6+rDt2VqLZ2U3lovHCet68eX6X\nr/WQoDOOq/t+4net8AXwY+Cp+v9+5qLeSCBUSllafxHVFMD/f6mf2FVIamvVZBFcaZqrPOxugZnl\ntrL7eLj+enPqsuv8AKofrrvOuPK1uK0Wcj5BYgRwHeo62kD4K/A/IcSd1IfqAgghOgKvSCmvRLm8\nPhEq7rAF8K6U0vRsvnYVktxcdXI0MtKc+uwqJBUV6pS9UQehGmNXy0NKtZhISTGnPrv2A8C+fdCr\nlzl12XV+AOMXVVrcVh85/16/eb42kErrQ2+beKmllHnAlfU/HwDSAqlHD+w6ae7bBz17mlefXSeL\ngwfVhGlGWCbYtx8KClTbzLjLA9SkaUe5OHtW9UUwcMC120pP/BG53oAJB//tgV1XFvv3m7e6Avv2\nw5495vaDXZXH7t3Qt6959dm1H7KzleIwI9sA2Fcuiouhrs7YfG9a9jyc7/WQqNDaR4xrkr2wq5BY\nMWnacYW1a5dxWUNdYdfxsHs39OljXn12tTz27g0uJgB27lRyYWS2AS1uq2jjqrc/dhWSHTvgcrfx\nbPpj1xXWrl0wbZp59dn1hPmePeYqD7tOmpmZMHCgefWFhqr9ptpa9bNdMGNR5dZtJYTo4ukxtln2\nISgkCjtbHv1NTNtp1/GQkQGDBplXn13Hw44d5vaDEPZcYDosDyPxZHks4ry7ypm4+sdGetY47Dgw\niovVxqBZZzzAnpZHVZVyU/TrZ16ddgzVlRK2b4c0E8NL7DgeQCmPRx81t06HIm3Vytx6PbF9u7Fh\nuuBBeUgpG6xrhRApqPxTlwF/MbRVNsKOK80tWyA11ZzsqQ7s2A+ZmSqaxKxwZbDnYiI3V20Qm5UU\nEexpeZSVwZEj5gYOgP0UaXU1/PADDB9ubD1eo62EEL2FEG8Ai4EtQD8p5fPGNss+2HGyWL8exowx\nt047ThZbt8KwYebWacfxsGWLeSeqHdhtwgTYtEktqsLDza3XbrKRmQlduqiL4ozE057HICHEB8DH\nwArUbYKvSClt1E3GY8cVtxXKw46TxcaNQeUBsGYNjBtnbp2hocryra01t15PrFsHY8eaX6/dZGPD\nBhg1yvh6PFke24AxwCpUTqlnhRDP1z//Mr5p9sBuq4q6OjU4zMqm68Bu/QDw3XcwYYK5ddpVeYwf\nb369dps016+3RnnYbYFplvLwtGF+Z/1/G2+aCxefNVvsJiBZWergj1nXzzqwWz/k50NRkbmRNXBe\neUhp7p6TOyoq1CbxSNPv2Dy/oIiIML/uxtTVKeXxiumXNthvYbVhAzzwgPH1eNowf8P46u2P3QaG\nFS4rsN/qatUqtdo2O7Y+JETVWVOjFInVbNyoFKiZQQMO7LSg2LsX2rQxL8eZM3bqh5ISlVXYjDB+\nkzICXbjYaWAAfPstXHKJ+fXaTYmuWgUTJ3p/zwjs5LpavdoalxXYa0x8912wHwCWLVNyYcbCJqg8\nvGCngVFbC0uWmHui2oHdlOjKlUHlAbB8Oeh075nP2GlMLFpkjVyAvfphyRKYOtWcuoLKwwt2Ghib\nNimz3Kxb0pyxU3TNoUNqv8Ps8FQHdlEeJSUqnv/SS62p3y4Lq8pKZZFPaXzXqUnYxaUrpVIeV1xh\nTn1aznn0EUKsEELsrP99sBDC5DOc1mEXAQH4+mvrVldgH0X6xRdw5ZXW5RKyy2SxeLGyOqzasLbL\neFi9WqXiMDKDrCfsMkfs3KlOuZt1VYMWy+MV4Heoe8wBdgCzDWuRzbCLgIC1pjnYZ9L8/HO45hrr\n6reL5fHVVzBjhnX122XSXLgQpk+3rn67zBFffqnmB7OiALUoj0gp5QbHL1JKCdhgyJiDXQTk4EF1\nzahVm4Jgj744cUK578zMKNwYOyiPqipleVx5pXVtsMOkWVcHn3xi7WLCDnIB8L//wQ03mFefFuVR\nJIQ4ZwgJIWYC+cY1yV7YQUBADYwf/cja8FA79MUXX8CkSRAVZV0b7KA8lixRrprkZOvaYIdJ8/vv\nITra/PM+zthBLvbuVTcomplpQMt9W/cB/wH6CiHygIPALYa2ykY4bxRbma//gw/g2Wetqx/sMVm8\n8w7cfbe1bbBDP7z3HtxisRTaYdL84AO46SZrD2zaYTwsWAAzZ5o7R3m1PKSU2VLKyUAHoI+UcpyU\n8lAglQohZgkhdgohaoUQQz28N1UIsVsIsU8IYdnthVYLSVaWWlVY6bIC6/shL08lAbTSzw/W90NZ\n2f9v70yjpaqOBfwViKCIqFEwOAQXwSmJkUTRPFf0uvJiRF6SpzHPh0N8GsckL4nGgcHh4oAMYlAU\nZbqIERGMgPOA5l0jGjVEFAFxRgbFIREDEYHLrfejurW9dPc9p/ucs8/pu7+1WPTE2cVee++q2rt2\nlQVPHH+8OxnA/aLZ1GSL5gknuJMB3I8HVTOqku6HIGVoOwE/AXoC7UVEsKOPK6po9yXgWGB8mXbb\nAzdiKeBXAX8VkXtV9eUq2q0I12kYJk2CU091X6nM9WIxfTocd5z7dBiut63uvtsMCVfRRXlcL5oP\nPgg9e8Lee7uTAdwHkjz1lHleSef1CrJtdQ+wBkvH/mkUjarqUgAp72v2BV7Pezm5DL8/BhJXHi4n\nyaefwm232d6ua1xOElWYMAEmT3bTfiGulcctt8Dgwe7az+N60Rw/Hs4+2137eVwbVRMnwhlnJL91\nF0R57KaqCV07+WK7wIqC9yuBBHJFbonLwTF7tlWI69XLTfuFuFSif/qTVfFLOvV4MVwqjxdesOJP\nLkNT87jsh+XLzaC66y437Rficl6sWWOh69dem3zbQZTH0yJygKouDPNgEZkLFKttNlhV7wvwiFCZ\ne+vr6z97XVdXR12EORtcDo4JE+AXv3DTdktcKtFx46wf0pDJ1qXFPX48nHmm+y1McNsPkybBiSe6\nSQjZEpfzYto0u1G+yy7Bft/Y2EhjY2MkbZdUHrkb5c1YrfLTROQtIF+9WVX1gHIPVtVqI/FXAYWJ\nOPbAvI+iFCqPqHE1OJYuhSVL3MawF+JKia5aZeknbr01+baL4criXrsWZsywSnFpwFU/bNoEDQ12\nzyUNuJoXqnDzzXD99cH/TUvDeujQoRW3X87z6AEciNXviJNSz58P9M7VTn8HOAFHN9tdDY7rroNz\nzkm+rGYpXCnRCRNgwADo0iX5tovhatGcOtXyWPXokXzbxXDlecycCb17J5N2PAiu5sXDD1vtele5\nzcopj2Wq+nYcjYrIscANWPjvAyKyQFX7iUgPYKKq9lfVJhH5FfAI5v1MdhFpBW4Gx3vv2X7uq68m\n2245XCjRTz4x6+rPf0623XK4UB5NTWZMTJuWbLvlcDEvVGHkSBg+PNl2y9Ghg43TpLn2WrjgAndb\nueWUxy4icj7FPQNV1esqbVRVZwOzi3z+DtC/4P1DwEOVthMVLhbNsWPt8lPQvcwkcGFp3nqrhSDu\nu2+y7ZbDhfKYNcs8DheFwErhYl7MnWspSZJKOx4EF/Pi+efNsHR5x6Wc8mgPpGSjwC1JW1jr1tnB\n6NNPJ9dmEJLuh6YmGD0a/vCH5NoMQtKLhSqMGgWXpCyXtQvPY9Qot9Z2MVz0w+jR8JvfuE1XVE55\nrFbVyk9TaoikLayGBjj8cNvXTRNJ98OsWVa/JOnLT62RtOfR2GiH5T/8YXJtBiHp8fD885ZtYUDK\ncnon3Q/Ll9t5x7hxybVZjCChum2eJC3NjRttb/vOO5NpLwxJWljNzXD11XDVVcm0F4aklceVV8LA\ngVY/PU0kbXGPHGnWdloCSPIk3Q+jRsHpp0PXrsm1WYxyyuPfE5Mi5SS5WEydah7HoYcm014YkrSw\n5syx9lznsSpGkuPhySctFf/JJyfTXhiSPChevNjCtSdOTKa9MCQ5L1autKCJpUuTaa8cJZWHqv49\nSUHSTFKex8aNZm3fcUf8bVVCUv3Q3Az19TBsWLr2tvMkaWkOHQpDhlhIZtpIsh+uuALOPz894dqF\nJNkPw4fDz38O3bol0145Ujgk00dSluaUKbDPPunb48+T1CSZPdtSkbgsdFSOpCzNp56CN96AU06J\nv61KSKofFi2yc5805DUrRlL9sGKFGZZp8DrAK49AJGFxb9hgXsfMmfG2Uw1JTJLmZrO2r7kmnV4H\nJGdMXHGFJUB0GVFTjqSMiaFDLcJqu+3ib6sSkvLIhw+3BIhp8DrAK49AJLFYNDRYZbg0nnXkSWKx\nmDULOnVKR+K/UnToYNmO4+SJJyyO/9RT422nGpIwJhYuhHnz0pOaphhJzIsVK6wkwSuvxNtOGLzy\nCEDclsWGDba/f/fd8bURBR06mKxxsXmzWZkjRqTX64D4jQlVuPhiizRLW2RRIUksmvX1cOGFbssO\nt0YSSnTYMEuImaZLw155BCDuxeKWW+Cb34S+feNrIwq23truG8TFHXfYgWi/fvG1EQVxGxOzZ5uS\nTtt9hpbEvWjOn29p12+/Pb42oiBuJfraa5aqKC1nHXm88ghAnIvFxx+bVfHYY/E8P0rinCSffgqX\nXmoLRZq9Doh3PGzaBIMGWXqatN3raEmc4yHvfV1+eTrSrpcjbiU6ZIhFmrmuHNkSrzwCEOckGTXK\nLO1vfCOe50dJnJNk3DjzvlzXaQ9Cx47xbd81NMAee8D3qy1okABxjoe5c+1Ow+mnx/P8KIlzfXju\nOYu6mzIlnudXg1ceAYhrkrz7rmWMXbAg+mfHQVyTZM0aiySJqEZN7MSlPP71Lzvzuffe9HtfEN94\naG42r2PYsPRGmhUS1/qQ977q69N55pNyxzgdxDVJhg41y2rPPaN/dhzENUlGjIAf/Qj23z/6Z8dB\nXMpjzBjLaXbQQdE/Ow7iGg933mlz7rjjon92HMS1jfnww7B6NZx2WvTPjgLveQQgjjQMr7xi0VVp\nCr1rjTgmycqVVuxpYagix26JQ3l88AH8/vfw7LPRPjdO4jCqNm607MFTpmTD+4J4+mHzZvM6rrkm\nndkFwHsegYhjcAwZYhefdtop2ufGSRz9UF8PZ50Fu+0W7XPjJA7lcfXVFl3Vq1e0z42TODyP8eOt\ndssRR0T73DiJox+mTbPIw7SUoC5GSnVauoh6cDzzjP1JW52K1oi6H5Yssf39NFVLDMLWW0erPN56\ny8bCy07qZFZO1J7oP/9pSvSRR6J7ZhJEbVTlIw+nTUu39+U9jwBEOTjyh2BDh8I220TzzKSIepIM\nGmSpxnfYIbpnJkHHjtEumpdeCr/+dXrSTgQl6vtPo0fDUUdZ1F2WiNqoGjcODjww/ZGHTjwPEfkp\nUA/sCxysqs+X+N0y4J/AZmCTqjq5Rhfl4HjoIfjww3SnnShFlP0wbx68+GK6c3mVIsptqwUL4PHH\n7aJo1ojS81i9Gm68Ef72t2ielyRR9kOWIg9dbVu9BBwLjG/ldwrUqeo/4hepNFFZ3Js3m6Wd5kOw\nckQ1SVThoousyFHHjtU/L2miVB4DB5rnkdakf+WIsh+uvNIMqp49o3lekkTpkWcp8tDJEqaqSwEk\n2Iae812/qPa4b78dtt8+feVEgxLVJJkzx6LXTjqp+me5IKpF87HH4M03LWdRFunUKZoEka+9BjNm\npC/9RlDatzeDaPNme10p+cjDF1+MTrY4Sbv9q8BjIrIZGK+qTuqIRTFJ8odg06en+xCsHFFsWzU1\n2VnHmDHpT79RiiiUR3OzeR1XX52Ni3DF2GabaJTHJZfAeeelL/1GGPKGVTXKIx95uPvukYkVK7Ep\nDxGZC+xa5KvBqnpfwMccpqrvisguwFwRWaqqT0YnZTCiUB433QTf+hYcdlg0MrkgCs+jocHCcn/w\ng2hkckEUyuOuu8yIOP74aGRyQadOsH59dc/461/t/KuhIRqZXJE3rDp1quzfZzHyMDbloapVZ+dR\n1Xdzf38gIrOBvkBR5VFfX//Z67q6Ourq6qpt/jO22aa6SbJmje1lZuEQrBzVeh759Bv33JNd7wuq\nVx4bN9o9n4kTs+t9gY2HzZvNm6zkDK8w+WEa02+EoVrDKqnIw8bGRhojWojSsG1VdBkRkW2B9qq6\nVkQ6A0cBQ0s9pFB5RE21ymPECLvsk4VDsHJUO0HGjIHvfjc76TdKUe0Z2IQJ0Ls3HHlkdDK5QMTm\nxoYNlSmPRx+FVauykfywNaoxrPKRhzNmRCtTMVoa1kOHllxSW8VVqO6xwA3AzsADIrJAVfuJSA9g\noqr2x7a8ZuUO1bcCpqnqoy7krWZvd9Wq7KXfKEU1E+TDDy39xjPPRCuTC/L3PFTDe1Br11qRp6xd\nhCtFfusqrOeQT36Y1cjDllQaiVgYeVjplpcrXEVbzQZmF/n8HaB/7vWbwIEJi1aUavZ26+stmiZL\n6TdKUU2o7lVXWfqNr341Wplc0K6dLXibNoWv9Dd6tKVbz9pFuFJUeh44fboZZcceG71MLqh0btxz\nj23nnnhi9DLFTQ3o/PipdNvq5ZctLDVLh2DlqHTbKqvpN8qRP/cIozzee8+KPGXxIlwpKvHKN2yw\nCKupU7N99lVIJUq0qQkGDzaDopooLVdk+LguOSrdtho82FzzHXeMXiYXdOxYWT9cdlk202+Uo5JD\n8yxfhCtFJV75zTfD179u6edrhW23DZ95e+pU6N4djj46HpnixnseAahkgjz9tFmY06fHI5ML8kq0\nuTl4lNCLL1pVuHHj4pUtacIqj9deszoVWb0IV4qwFvfHH9s5x+OPxyeTCzp3tu2noKxfb1vaf/xj\ndr0v73kEoFMnWyhUg/0+H4J4xRXZOwQrR7t2pkDCWFiDBtkWRZcu8cnlgrDK45JL0lmHulrCeuWj\nRsExx5jnUUuEVR5jx8Ihh9ifrOI9jwC0a/d5eGYQZXD//fDRR3DKKfHLljT5SRIkF1NjoxW7mjMn\ndrESJ0y47vz5tXERrhhhvPKslV0OQ5htq3/8w5TovHnxyhQ33vMISNBD83zyw+HDs3kI1hpBLay8\n93XVVeEjkrJAmLTsAwfauU/WL8IVI8y2VdbKLochjOcxfLiV2N1nn3hlihvveQQkb2G1dvh9223w\npS9B//7JyJU0220XbJLMmmWRWSecEL9MLgi6bTV3LixfXhsX4YoRVHlksexyGIIqjxUrYPJkeOml\n+GWKG688AhJkb3f9erMwZ87M7iFYa3TuDOvWlf9NPgRx7Nhsp98oR5BFsxaSH7ZGUI98yBC48MJs\nlV0OQ+fOwbat6uvhnHOgR4/YRYodrzwCEmSS3HgjHHwwfOc7ycjkgiAWVkODZQb9ftXZzdJLly5W\nNrUcd91lyjPLyQ9bI4gSfeYZePbZ7JVdDsO227Y+L5Ysgfvus8i7WsArj4C0djD40UcwciQ8mXjO\n32RpTXl88ontbc+ZU7veF0DXrhZ2WopNm8zaHj++tvuhNeWR5bLLYejc2QICyjF4sHmiXbsmI1Pc\neOURkNY8j+HDLdXCvvsmJ5MLWjvzuP56Szt/8MHJyeSCHXYorzwmToReveB730tOJhe0Ni8efNDy\nmv3sZ8nJ5ILWtq2eesqizO68MzmZ4sYrj4CUO/NYsQImTaqN5IetUe7M4+9/t1QLTz+drEwuReZL\nWQAACL1JREFUKOd5rFtnUWYPPJCsTC4o53lkvexyGMptW9XsvS/XAmSFchZWfT2cfXZtJD9sjXLb\nVtdcAz/9Key9d7IyuaCc8hgzBo44Avr0SVYmF5RbNG+/3fopq2WXw1BuXtx/v42Vk09OVqa4qXF7\nIDpKHZAuXmyHYLWS/LA1Sm1bLV8OU6bAokXJy+SCrl3h7be3/PzDD0151ELq+SDsvHPxlCu1UHY5\nDKWUR977Gjmy9u59ec8jIN26WVbUluQPweKuAJYWSk2Syy+Hc8+FL385eZlcUOrMY9gwu9tSC6nn\ng9CtG3zwwZaf10LZ5TCUOvPI3/s65pjkZYob73kEpHv3LZXHvHnwwgvJVABLC507wzvvfPGzRYts\nf79WQhCDUGzb6u23LVPq4sVuZHJBt27w/vtf/CxfdvmJJ9zI5IJi23fr15tRNWNGbXpf3vMIyK67\nflF55A/BslgBrBq2337L7bvBgy0BYq2EIAahmPK47DL45S9trLQVdtllS+WRL7u8335uZHJBsX64\n6Sb49rdr996X9zwC0r07rF79+ft777WSoied5E4mF+y5Jyxb9vn7fP3lmTOdieSErl3Nws6zcCE8\n/HDb8r7gc88jX5J35craKbscht12s37YuNFyua1ZY+cctex9ec8jIIXbVk1NZmnXavLDcvTqBW+8\nYa/bqvcFduZRqDwuvtjSrm+/vTuZXJBP9pjfsqmvh7POahuRh4VstZWlHFmxwt63Be/LiechIqOA\n/wA2Am8Ap6nqFsePInI0MAZoD0xS1RGJClpAoecxaZJZXP36uZLGHbvvbhFF69dbCOK6dW3P+wJb\nKNautfTazz9vHsfZZ7uWKnlEPj80f/NN88jbSuRhS/bay0out2/fNrwvV57Ho8DXVPWbwKvAoJY/\nEJH2wI3A0cD+wAARcabHu3WzPe7XX7e97RtuqM1DsMbGxrLft28PX/mK5em58EK7Ud7WvC+w//NB\nB8Ff/gIXXWR3XIqlnm+tP2uBnj0tS+y555oXGmfkYZr7s2dP29K94AI477za976cKA9Vnauqzbm3\nzwK7F/lZX+B1VV2mqpuAO4EfJyVjSzp0sPQjhx9ulvYBB7iSJF6CTM4+fSwc9ZBDoK4udpFSy6GH\n2oK5446lkx+mebGLigEDzOtqaoIzzoi3rTT3Z9++Fqq9YAH87neupYmfNByYnw4Uq/S9G7Ci4P1K\nwGnRxksvtfKZ55/vUgr3jBtnF+EGDnQtiVt++1uLsjn55Nr0QoNy4ol2UH7OOW3TC81z1lm2pTtg\nQG0ngcwTm/IQkblAsaDFwap6X+43Q4CNqnpHkd8FrBieHPvtV9sHYEHZaSfL09PW6dbNtifaOl26\n2HZVW6ddOwuaaCuIqps1WkT+BzgT+J6qbpFaTUQOBepV9ejc+0FAc7FDcxFJnaLxeDyeLKCqFfnN\nrqKtjgYuBI4opjhyzAd6i0hP4B3gBGBAsR9W+p/3eDweT2W4irYaC2wHzBWRBSIyDkBEeojIAwCq\n2gT8CngEWALMUNWXHcnr8Xg8ngKcbVt5PB6PJ7tk+oa5iBwtIktF5DURudi1PFlERJaJyMKcB/hc\n7rOdRGSuiLwqIo+KSBvJGRwOEWkQkfdE5KWCz0r2nYgMyo3VpSJylBup00uJ/qwXkZW58blARPoV\nfOf7swwisoeI/J+ILBaRRSLy69znkYzRzCqPtF0izDAK1KlqH1Xtm/tsIDBXVfcGHs+992zJFGz8\nFVK070Rkf+zcbv/cvxknIpmdfzFRrD8VuC43Pvuo6kPg+zMgm4DzVPVrwKHAL3NrZCRjNMudnapL\nhBmnZcDBj4CpuddTgf9MVpxsoKpPAh+1+LhU3/0YmK6qm1R1GfA6NoY9OUr0J2w5PsH3Z6uo6mpV\nfSH3eh3wMnZ/LpIxmmXlUewSYY0nBIgFBR4Tkfkicmbus+6qmk9A/x7Q3Y1omaRU3/XAxmgeP16D\n878i8qKITC7YYvH9GYJc1GofLKNHJGM0y8rDn/RHw2Gq2gfoh7m13y38Ui2iwvd1BQToO9+vrXMz\nsBdwIPAuMLrMb31/FkFEtgPuBn6jqmsLv6tmjGZZeawC9ih4vwdf1JqeAKjqu7m/PwBmY27qeyKy\nK4CIfBl4v/QTPC0o1Xctx+vuuc88ZVDV9zUHMInPt1F8fwZARDpgiuMPqjon93EkYzTLyuOzS4Qi\nsjV20HOvY5kyhYhsKyJdcq87A0cBL2H9eGruZ6cCc4o/wVOEUn13L/DfIrK1iOwF9AaecyBfpsgt\nbnmOxcYn+P5sFRERYDKwRFXHFHwVyRhNQ2LEilDVJhHJXyJsD0z2lwhD0x2YbWOMrYBpqvqoiMwH\nZorIz4FlwH+5EzG9iMh04AhgZxFZAVwGDKdI36nqEhGZiV14bQJ+of6S1Rco0p+XA3UiciC2ffIW\ncDb4/gzIYcDJwEIRWZD7bBARjVF/SdDj8Xg8ocnytpXH4/F4HOGVh8fj8XhC45WHx+PxeELjlYfH\n4/F4QuOVh8fj8XhC45WHx+PxeELjlYfHEwIR+VJBevB3C9KFrxWRG13L5/Ekhb/n4fFUiIhcDqxV\n1etcy+LxJI33PDye6hAAEakTkftyr+tFZKqI/DlXbOs4Ebk2V3TrIRHZKve7b4tIYy6j8cP5fEMe\nTxbwysPjiYe9gCOx2gm3Y8V3DgDWA/1zCevGAj9R1YOwQkhXuxLW4wlLZnNbeTwpRoGHVHWziCwC\n2qnqI7nvXgJ6AnsDX8NqqYDlZ3vHgaweT0V45eHxxMNGAFVtFpFNBZ83Y/NOgMWq+m8uhPN4qsVv\nW3k80VOsbGpLXgF2EZFDweou5GpIezyZwCsPj6c6tODvYq9hy2psqqqbgOOBESLyArAA+E6cgno8\nUeJDdT0ej8cTGu95eDwejyc0Xnl4PB6PJzReeXg8Ho8nNF55eDwejyc0Xnl4PB6PJzReeXg8Ho8n\nNF55eDwejyc0Xnl4PB6PJzT/D2OUMcasn02eAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "run_simulation(np.linspace(0, 200, 501), model=fnm)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": false, "scrolled": true }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEPCAYAAABV6CMBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcE/X5wPHPw4IioCKggAiiVQ6RU8WzulZRvK0iHuBR\nbKW2KlrF82fd2tZq64HWetejKoJ4IOCJ6NajgoogN4KKIqeCIJfA7j6/P74TN5tNdiebTGaSfd6v\nV14zSSYzz46SJ99bVBVjjDHGrwZhB2CMMSa/WOIwxhiTFkscxhhj0mKJwxhjTFoscRhjjEmLJQ5j\njDFpCTVxiMijIrJCRGbWctwBIlImIqflKjZjjDHJhV3ieAzoX9MBIlIE3Aa8BkgugjLGGJNaqIlD\nVd8Fvq/lsEuB54Bvg4/IGGNMbcIucdRIRNoBpwD3ey/ZMHdjjAlZpBMHMAK4Vt28KIJVVRljTOga\nhh1ALfYDRokIQCvgOBHZqqrj4g8SESuJGGNMHahq2j/II13iUNU9VXUPVd0D185xcWLSiDs28o+b\nbrop9BgKJc5cxbhpkzJtmvLHPyrduyuutjSdx011+EzqR9u2yjHHKFdcodx9t/L888qUKcpXXymr\nVys//qhUVPj/+7ZuVX74Qbn88pt4913l8ceVyy5TevbMTrwDByobNuTff/f6EmddhVriEJFngCOA\nViKyGPevrBGAqj4YZmymcKjCypUwfjw8+SS8807w1+zTB44+GpYtg3/9C5o1A4lgRWvDhrD99rDj\njnDYYe5x/vnpnaO8HMaNg4svhhUrqr737LPuEa97d5gyBbbbLrPYTXhCTRyqenYax/4qyFhMflCF\nVatg1iwoLob//jeY63TqBGeeCaecAr16QVFR3c5TUuK+mAtZURH88pfukUgVXn4ZTjqp8rWZM6FJ\nk6rHrVvnkqvJD1Fv4ygoxcXFYYfgSy7jVIWNG+HLL+Hhh+Gee/x+stjXUX37wpAhLgG0bp37X/31\n/b+5CJx4ovvvHO/OO+HKKyufxyfXWbOgW7fk56vv9zMqJJN6rqgQES2Ev6OQlJfD3LnwyCNw993Z\nOWdREQwbBuedB/vuW/dSgIkWVRg0CJ55pvp7S5dC27a5j6m+EBG0Do3jljiMb6rw7bfw6qtw++3u\nl2GmDj0ULroITj0Vdtgh8/OZ/Hf//fC731V9bcgQ+Pe/cx+LRLFhqo6SfUda4iiAvyNMW7bAZ59B\naSncdBOsXl33c+2xBwwfDqed5qqHjKmrgw+GyZOrvlZRkbsqR++LNTcXC1Cqv8MSRwH8HUHasgX+\n9z+46y7XA6YuevSA66+Hfv1gp52i2UvIFKbJk10SiZeLBFJIieOuu+6iadOm/OY3v6nyel0ShzWO\nF4gNG2DJEnj3Xfj1r9P//MEHw29/6xqRd9wx+/EZk4mDDnJVpWvXQvPm7rUG3ii0Avhez4lddtmF\nlStXZuVcljjyyMaNMHUqPPggPP10ep896ii4+mrXTz+xK6Qx+WLHHV2iWL0aWrZ0r4m4NpF//Svc\n2OqTSI8cr4+2boUZM+BPf3L/IOIfTZvC4YenThrXXefGOKhWf7z5JhxzjCUNUxhatHD/X0+d6p7f\nd5/7N/Ldd+HGVV9YiSMEFRUwZw6MHOn+h1+71v9n77oLzj238teWMfVZnz4ugZx9NowaBTvv7Lrv\nLl0admSFzRJHDpSVuWkXBg3yd/zZZ7vBUX36WAO0MX488ww89ZSbQmXZMvfv5ocfCn/UflisqiqL\nYgnirLOqVjE1apQ8aZx7rquWqqioWq00ciTst58lDWPSUVTk/v0MG+ae77AD3HxzuDEF5bbbbuOM\nM86o8tqwYcMYFvvjA2bdcTP01VeuF9Obb9Z83DPPwJFH2rgGY3Jh3bqqA0rr+vUQ1e64X3/9NV27\ndmXFihU0a9aM8vJy2rdvz9ixY+nbt2+140WEp59+mpUrV3L55ZdXeb0u3XGtxJGmDRtgwoTK0kTH\njtWTRnExfPFF1VLEWWdZ0jAmV7bfvmqyEAmm225iB5a6PtLVoUMH+vTpw4svvgjAW2+9RZMmTZIm\njSBY4qjFwoVwxx1w3HHuP3CzZlVn+gQ3/mH58sok8fbbbvS0MSZcsYZzcOM+vv02++fPxqMuzjnn\nHJ7xJvgaOXIkg/w2omaBNY6nMHOmGymdTM+e8NJLsPvuuY3JGJO+kSPdXFf9+sEuu8DHH7s2xHw3\nYMAArrzySpYsWcLYsWOZnDg3S4CsxOHZvBmee66y6JiYNO65p/LXwfTpljSMySdHH+1qBQD23x/G\njAk3nmzYeeedKS4u5oILLmDPPfekc+fOObt2qIlDRB4VkRUiMjPF+4NE5FMRmSEi74tIijJA3Y0d\n6xJF48aQ0EmB55+vTBaXXprtKxtjcql1a/cDEWDgwOxN9x+mc845h0mTJnHOOefk9LphV1U9BvwT\n+E+K978ADlfVtSLSH3gIOCgbF966FbbZpvrrd94Jl19uXWGNKUTbbOO6vzdo4P6dl5VVXVAq3wwe\nPJjBgwfn/LphLx37roh0rOH9D+KeTgF2y/Sa770HP/951ddatoQFC9yMr8aYwiZSmTyuusr9ux8y\nJOyo8ks+tXFcCLxS1w9XVLiGsfik8bvfwfr1bn4bSxrG1B+x5AFw4YXw2mvhxpNvwq6q8kVEjgSG\nAIemOqakpOSn/eLi4mpr/iYuM7p+vZs00BhTP4m4JY6Lilx3+3nzIIfty6FYvnx5le/Kugp95LhX\nVTVeVbuneL8H8ALQX1UXpjgm5cjxqVNdL4qYF190y5QaYwzApk2Vs0Zv2FB1BumojhxPV7ZHjke6\nxCEiHXBJY3CqpFGT8ePh5JMrn69da+taG2Oq2m47N9PDnnu6WogCyBOBC7s77jPA/4DOIrJYRIaI\nyFARGeod8kdgJ+B+EZkmIh/6PfemTVWThqolDWNMcnvs4WbXBetR6UfYvarOruX9XwNpL4Q6bx50\n7er2GzRw9ZjGGFOTQYPgoYfgnXeqzm0llkmqiXRVVV3FkgZY0jDG+Pff/1aWOIYP56f2jREjRrDL\nLruEGFm0FFzimD27cj/W3c4YY/xSdcnj9tvhiitg112hadOmrFy5MuzQMtY0S11JQ+9VlQ3xvapi\nvxY+/TT1JIXGGFOTLVtg223dfgF8RaZk63HgkkWMJQ1jTF1ts41bfA2ssTyZgkocjzzitvffH24c\nxpj8d9ZZlfv5PJ9VEAqqqir2y6Ciwn4lGGOyI/ZdsnGjG/NRSKyqKo4lDWNMtqxd67bxI8rru4JM\nHMYYky077ODW7wDo0yfcWKKiYBLHV1+57SmnhBuHMabwjB7tttOmwapV4cYSBQWTOD77zG1zvBCW\nMaae2LLFbVu1CjeOKCiYxLHjjm578MHhxmGMKUyNGsG557r9E04IN5awFUzi+O47t21YcGPhjTFR\n8R9vketXXqnf0xkVTOL42c/cdsmScOMwxhS2WBtHff6RWjCJo1Mnt/3f/8KNwxhT2Fq0qNwfMya8\nOMJUcAMAe/VyPR+MMSZIsfFi+fwVagMAPdOnhx2BMaY+iM1l1bJluHGEIdQSh4g8CpwArKxhzfF7\ngOOAjcAFqlqtPJE45Ug+/wIwxuSP2HdOWRkUFYUbS13ka4njMaB/qjdF5HhgL1XdG7gIsOkLjTGR\nsWaN29a3hvJQE4eqvgt8X8MhJwNPeMdOAZqLSOtcxGaMMbWJjR8DWLAgvDhyLewSR23aAYvjnn8D\n7BZSLMYYU02sajzWs7M+qLWAJSIDgddU9QcRuRHoA/xZVT8JPDovhITnSVswSkpKGDQInn4aJkwo\n5sQTi4OPzBhjgEMPhfffh5Ejoz3tUWlpKaWlpRmfp9bGcRGZqardReQw4C/A7cCNqnpgxld35+8I\njE/WOC4iDwClqjrKez4POEJVVyQcV2Xp2IMOgg8+yEZ0xhjjTz52zgmycTw2sP5E4GFVnQBsk+6F\n6mgccB6AiBwErElMGslMnhx0WMYYU1VJidtedVWoYeSEnxLHy8ASoB/QG/gRmKKqPTO+uMgzwBFA\nK2AFcBPQCEBVH/SOuRfX82oD8KtkVWSJJQ73+UyjM8aY9OTb909dSxx+EkdT4FhgpqouEJG2QHdV\nfaNuoWafJQ5jTBSMGeMWferfH159NexoahdYVZWqbgC+BQ7zXioDFqZ7oVx58EG33bQp3DiMMfXP\nGWe47WuvhRtH0GpNHCJSAlwNXOe9tA3wZIAxZeSii9z25z8PNw5jTP00aZLbHnJIuHEEyU9V1ae4\nto2pqtrbe22GqvbIQXy+xFdVuedua9VVxpgw5Mt3UJC9qjarakXchZqmexFjjKlP3nzTbY8+Otw4\nguKnxDEc2As4BvgbMAQYqar3BB+eP4kljsaNYfPm6Gd7Y0zhyodSR5CN4/8AnvcenXCD/yKTNJJZ\nutRtR40KNw5jTP310ktue9554cYRhIJayKnqa25bAH+eMSZPRf17KOslDhFZLyLrUjx+yCxcY4wp\nfHff7bZ33hluHNlmJQ5jjAlQlL+LAmvjEJEOyR51CzN3YkvIfvRRuHEYY+q3005z23ffDTeObPLT\nq2oWlVOZNwb2AOarareAY/MtWYnDve62Ucz0xpj6I6rfRXUtcdS6Hoeq7ptwoT7A79O9kDHG1Her\nVkHLlmFHkbk6tXGIyKzEhBImK3EYY6Iuit9HgZU4ROTKuKcNcCsALkn3QmGYPRu6dYPSUiguDjsa\nY4wpDH7aOEqobOMoAxYBz6vqj4FGloZUJQ73nttGKcsbY+qftWuheXO3H5Xvo8DW48gHljiMMfkg\nat9HQXbHnSgizeOetxCR19O9UIpz9xeReSKyQESuSfJ+KxF5TUSmi8gsEbkgG9c1xpgwvPCC2956\na7hxZMpPVdV0Ve1V22tpX1ikCJgPHI1rM/kIOFtV58YdUwJsq6rXiUgr7/jWqlqWcK6UJY5vvoH2\n7eG55+D00zOJ2BhjMhelUkeQ06qXi8jucRfqCFSkPNq/vsBCVV2kqluBUcApCccsA3bw9ncAViUm\njdrstpvbDhiQUazGGJMVe+3ltsuXhxtHJmrtVQXcALwrIu94zw8HLsrCtdsBi+OefwMcmHDMw8Bb\nIrIU2B4YmIXrGmNMaBYscKWOtm2jUeqoCz8DAF8Tkf2o/FK/XFW/y8K1/dyy64HpqlosIj8DJopI\nT1Vdl3hgSUnJT/vFxcUUW/9bY4yporS0lNLS0ozP46eNowEwCNhDVW/25qlqo6ofZnRhkYOAElXt\n7z2/DqhQ1dvijnkF+Kuqvu89nwRco6ofJ5wrZRsHwOrVbrTmiBEwbFgmURtjTOYWL4YOHeDQQ+G9\n98KLI7DuuCLyAK5N40hV7SoiLYA3VHX/uoX603kb4hq7jwKWAh9SvXH8TmCtqv5JRFoDU4Eeqro6\n4Vw1Jg53jNvma9HQGFNYovCdFGTj+IGq+jvgRwDvS7tRuhdK5DVyXwK8DswBRqvqXBEZKiJDvcNu\nAfYXkU+BN4GrE5OGMcbkozPPdNv588ONoy78lDimAIcAH6tqbxHZGVfi6J2LAP3wU+Lo3dtNtV5R\nUZnpjTEmTGGXOoIscfwTeBHYRURuAd4H/pbuhcI2darb9usXbhzGGJPvfE05IiJdcW0RAJPi2yGi\nwE+Jwx3nttbOYYyJglgj+bHHwmuv5f76WW8c9xrBq7zkbRV+auuIBEscxph8Feb3UhBVVZ/gejHF\nHh97j9h+3vnPf9x20aJQwzDGmJ907eq2338fbhzpKPjZcasf67YF8GcbYwqAKjRoULmfS4Et5OSd\n/BTcVCMK/FdVx6d7IWOMMdXlYy9PP9Oq3wpcBswG5gKXiUje9aoyxpioik23PnJkuHH45Wccx0yg\nl6qWe8+LcPNHdc9BfL6kU1W1ahW0agW33QZXXx1wYMYY41MY1ehBTjkyAzfdyCrveUvgbVXtUadI\nA5BO4nDHu621cxhjoiKfEoefAYB/Az4RkSdE5Alcr6pb0r2QMcaY1L791m2PPTbcOPzwOwBwV+AA\nXOP4h6oaqSVI0i1xXHwxPPAAbNgATZoEGJgxxqQh16WOIAYA9kl8ydvGBgB+ku7FgpJu4nCfcVur\nrjLGREXse2nLFmiU8VSyfq6X/cRRAcwCViV7X1WPTPdiQbHEYYwpBGVllQkjF99NQYzj+ANwBrAR\nGA28mGzlPWOMMdnR0NfIuvClbBxX1RGqeihuDMduwCQRGSMivXIWXYC++spt77473DiMMSbecce5\n7XfZWKA7IH4bx7sBZwODcUu3js7KxUX6AyOAIuCR+GVj444pBu7CLR71naoWJzkm7aoq9zm3teoq\nY0xU5HIKkiDaOH4GnAWcAnyNq66aoKqbMgk07vxFuKVjjwaWAB9RfenY5rj1P45V1W9EpJWqVsvD\nljiMMYUkV99NQYzjWAAMBF4FPgA6ABeLyJUi8oe6hVlFX2Chqi5S1a3AKFySincO8LyqfgOQLGlk\n4t573TZWbWWMMVFw8cVu+8UX4caRSk0ljhK8rrfJqOqfMrqwyABcSeI33vPBuPXNL407JlZF1Q3Y\nHrhbVZ9Mcq46lTjcZ93WSh3GmKioqICiIrcf5HdT1ntVqWpJRhHVzs/taAT0wa0+2AT4QEQmq+qC\nQCMzxpgQNfAzp0eIwuz8tQRoH/e8PfBNwjGLcQ3im4BNIvIO0BNXjVZFSUnJT/vFxcUUFxf7CqJt\nW1i2zGX1fJze2BhTmIYNc70+Z8+Gbt2yc87S0lJKS0szPk9oCzmJSENc4/hRwFLgQ6o3jncB7gWO\nBbYFpgBnquqchHPVuaqqvNz1ne7ZE6ZPr9MpjDEm63IxGDDrjeMiMszbHpZJYKmoahlwCfA6MAcY\nrapzRWSoiAz1jpkHvAbMwCWNhxOTRqZi9YiffprNsxpjTGaiPBiwpsbxT1W1p4hMU9XeOY4rLZmU\nONzn3dYayI0xUTJkCDz2GMyfD506Zf/8QYzjeAbYH2gHfJ7wtubzehyJFi6EvfeG+++H3/42i4EZ\nY0wGNmyAZs3cfhA/bANZyElE2gBvACdROTsuAKq6KN2LBSXTxOHO4bZW6jDGREmQ302BLOSkqsu9\nksUyoJn3WBKlpGGMMYVs773ddnmEVkGqtbewN1fUAuA+77FARI4IOK6cGzvWbd97L9w4jDEm3vPP\nu23btuHGEc/PmuOf4LrJzveedwJGqWriQk+hyUZVlTuP21p1lTEmKoKc9DDINccbxpIGgKp+RrgD\nB40xpt6IH5i8eXN4ccTzkzimisgjIlIsIkeKyCPAx0EHFoYRI9x2gU1oYoyJkFGj3LZr13DjiPFT\nVdUY+D1wqPfSu8B9qhqR3Je9qip3Lre16ipjTFSsWgWtWrn9bH43BdIdN19Y4jDGFLogvpuCbOOo\nV2Lz4Ed52UZjTP1zxhlue8MN4cYBVuJIcT63LYBbY4wpEPfc42bMhex9NwVS4hCRIhG5ve5hGWOM\nyYb99w87gkq1jRwvBw4TqV8rVRx/vNuuWxduHMYYE9OrV+X+tGnhxQH+elU9AOwKjAE2ei+rqr4Q\ncGy+Zbuqyp3Tba26yhgTFT17wowZbj8b301ZXzo2TmNgNfCLhNcjkziMMaY+2H//ysQRploTh6pe\nkIM4IueCC+Dxx2HFCmjdOuxojDHGJY5HH3X7YS537WeSw84iMklEZnvPe4jI/2Xj4iLSX0TmicgC\nEbmmhuMOEJEyETktG9f147HH3LZNm1xd0RhjataxIxxzjNu/8MLw4vAzjuNh4Hpgi/d8JnB2phcW\nkSLceuL9gX2As0Wk2oB677jbcEvI1qtGemOMidesGWz0WppjP27D4CdxNFHVKbEnXiv01ixcuy+w\nUFUXqepWYBRwSpLjLgWeA77NwjXTcsstbjsnq6ucG2NM3TRrFo3enn4Sx7cislfsiYgMwC3slKl2\nwOK45994r/1ERNrhksn93ks57eN03XVu261bLq9qjDHJNWsG69fDuHHueVgzXPjpVXUJ8BDQRUSW\nAl8Cg7JwbT9JYARwraqqN5bEqqqMMfXW9tu7xHHSSe55ly7hJA8/vao+B44SkaZAA1XNVkFpCdA+\n7nl7XKkj3n7AKG/8YSvgOBHZqqrjEk9WUlLy035xcTHFxcVZCXLUKDjrLJg4Efr1y8opjTGmTmIl\njphVq9L7fGlpKaWlpRnH4WcAYCvgJuAwXCnhXeBmVU0z5GrnbQjMB44ClgIf4lYanJvi+MeA8ckG\nHgYxALDq+d3WBgMaY8JUUQENG0JZGRQVudcy+V4KcnbcUcBK4DRgAK6RenS6F0qkqmW4arDXgTnA\naFWdKyJDRWRopuc3xphC06ABNGkCGzbAe++511asyH0cfkocs1R134TXZqpq90AjS0PQJY7Zs2Hf\nfV1jeaynlTHG5FpZGey0E3z9tduKQIcO8NVXdTtfkCWON0TkbBFp4D3OBN5IP8T8FetV9be/hRuH\nMaZ+e/996NTJJY2Yr7/OfRwpSxwisp7Knk9NgQpvvwGwQVW3Dz48f4IucbhruG15uSsuGmNMrl1x\nBbRoATfe6J5n2v6a9RKHqjZT1e29RwNVbeg9GkQpaeRKWZnbxhqkjDEml1Rh7Fg49dTK16691m1j\n30+54msFQBHpAXQkrvtuoU+rnvw6bmu9q4wxuTZ9Opx+OixcWLUGpGFDt+T1ffelf87A2ji8brCP\nAqcDJ8U96p3LLnPbefPCjcMYU//EShvxM+LGakDuvz/5Z4Lip1fVHKBbTn7S11GuShzuWm4b3bth\njClEvXrBvffCYYdVfT2T76Qge1V9hJu91hhjTAi+/BKWLoWDDw47EsdP4ngM+EBEPhORmd4jAmtQ\nhSO2+taVV4YbhzGm/hg7Fk4+OXnnnCZNch+Pn6qqz4ErgFlUdslFVRcFGlkacllV5a7ntlZdZYwJ\n2ttvw5lnwosvwqGHJn//F79w81a1aJHeuYOsqlqpquNU9Qtv7YxFUUoaYdhuO7eNwrz4xpjCNXq0\nSxqjRydPGgCx+Vz/8pecheWrxHEf0BwYT+UqgFofu+PGqFYOArRShzEmCHfeCXfdBS+/DD161Hxs\nXWtB6lri8LMeRxNcwjgm4fXIJI5cC2uBeGNM4auogKuugtdfd1OMdOgQdkTV+RoAGHW5LnEATJsG\nffrARRfBgw/m9NLGmAK1eTOcf77rQTV2rP82i1yXOPxUVSUuia4Aqjok3YsFJYzE4a7rtgWQe40x\nIVuzBn75S2jZEp56Cho39v/ZXCcOP43jLwMTvMckYEdgQ7oXKkR9+rjt0qXhxmGMyW9LlsDhh0P3\n7q4hPJ2kEYa0q6pEpAHwvqpGZChKeCUOd223tVKHMaYuZs+G44+HSy5xbRt1aUONYokjUSdg5zp8\nrhoR6S8i80RkgYhck+T9QSLyqYjMEJH3vckWjTEm761fD48+6sZg3HILDB+ePx1v/ExyuF5E1nmP\nH3Ddcqt9yadLRIqAe4H+uClNzhaRrgmHfQEcrqo9gD8DD2V63WxbssRt99wz3DiMMdGnCu+8A0OG\nQPv2rgH8+edh0KCwI0tPaL2qRORg4CZV7e89vxZAVW9NcfxOwExV3S3Je6HOwWjVVcaYmnz9Nfzn\nP/D447DttvCrX8HgwdCmTXbOH5mqKhHpUNMj3Qsl0Q5YHPf8G++1VC4EXsnCdbPugQfc9t57w43D\nGBMdGzfC009Dv37Qu7frRPPMMzBrlmvLyFbSCENNAwBfoXLp2Hg7e49M18LznRtF5EhgCJBi0D2U\nlJT8tF9cXExxbBx+DgwdCr/9LVx6qWvgMsbUT6owebIrWYwZA337wq9/DaecElxPqblz3Ta2XlBN\nSktLKS0tzfiavquqRKQjcC1wNHC3qv4zowuLHASUxFVVXQdUqOptCcf1wI1S76+qC1OcK/TlQvbf\nH6ZOhfnz3WLyxpj6Y+lSePJJlzDKy11V1Lnnwm7VKtaz7xe/cBMdPvYYXHBBep8NcgBgJ+B64CDg\nDuBxVd2a7oWSnLchMB84ClgKfAicrapz447pALwFDFbVyTWcK/TE4eJw2wiEYowJ2ObNMG6c+8L+\n4AMYMMAljIMPzm3vqNi11q+Hpk3T/WyW56oSke7ADUA34O/Ahapanu4FUlHVMhG5BHgdV+31b1Wd\nKyJDvfcfBP4I7ATcL+7ubFXVvtmKIShbt0KjRmFHYYzJNlX45BOXLEaNgp49XbIYMyb9L+1s2Lix\ncj+X109Z4hCRclyD9QTi1uHwqKr6qFHLjaiUOGILx4OVOowpFGVlMHMmvPWWq4rasMFVCZ13HnTs\nGG5ssdLGgAEueaX/+SxXVYnIBd5u4gGCSxxPpHuxoEQlcYBVVxmTz1Rh8WKYMqXyMW2am6H20EPd\neIvDD69cViFM339fOQliXRZxggDbOPJBlBLHypXQurXbj0hIxpgU1q2Djz6qmigqKuDAAysfBxwA\nO+4YdqTVxX6kPvqoqy6r2zkscYQdxk+s1GFM9JSXu3mhYgli8mT48kvo1atqoth99+hP/XHrrXDd\ndW4/k+8ZSxwR+js+/xz22ssVZ8uz1p3AGJOOJUuqliSmToV27aomiR498q8jy8KFsPfebn/MGNe+\nUVeWOCL2d1ipw5jcWb/eJYb4RPHjj1WTRN++sNNOYUeamY0bK3tP3Xsv/P73mZ0vyHEcnYH7gDaq\n2s0bkHeyquZwafSaRTFxzJzpfs3stptrbDPGZEd5uRstHZ8kFi50a1nEJ4o994x+lVM6Nm2CJk0q\nn2fjKy/IxPEOMBx4QFV7ixtQMUtVu9Ut1OyLYuIAK3UYk6nycvjiCze/04cfuiTx8ceuA0p8kujZ\n000eWKiWLKk6Cn3JEth118zPm/UBgHGaqOoUbwAeqqoikvHI8fpgyhT3P3WnTvDZZ2FHY0x0bdkC\nCxbAnDnuMXeu2y5Y4CYD3Gcf17tp+HBX5dSyZdgR586zz8KZZ1Y+X706/Co3P4njWxHZK/ZERAYA\ny4ILqXD09ca4L1gQbhzGRMXGjW4+t8QEsWiR6820zz7QtSuceCJcfTV07hzOiOwo2LoVttmm8vlR\nR8Gbb4aIe4VVAAAUa0lEQVQXTzw/VVU/wy2gdAjwPfAlMEhVFwUenU9RraqCyrYOsCorU3+sXeuS\nQiwxxJLEsmWux+E++1QmiX32cb2ECrmqKV233AI33FD5/IQTYMKE7F8n8F5VItIUaKCq69K9SNCi\nnDigsq2joqKwGuuM+e67qiWH2P6aNdClS2ViiCWJPfesnJbHVPfWW65kESPi7mmXLsFcL8jG8cbA\n6UBH3GSEsSlHbq5DnIGIeuJYs6ayTjLCYRqTlKorKSRWL82Z46pTEksPXbu6KTqiMC1Hvvj4Y9eG\nE2/NmuBHrAfZOP4SsAaYCvyY7gUMNG9eub9uHWy/fXixGJPK5s2uW+v8+a4zx/z5MG+eSxCNG1cm\nhm7d3KCzffZxDddWiq67F1+E006r+tojj8CFF4YTj19+ShyzVHXfHMVTJ1EvcYCrpiry1kyMeKim\ngKm6RYfmz6+aIObPd108d9/dNUh37ux6A8aqm1q1CjvywlFR4VYMffjhqq8//rhbKTD+h2bQgixx\n/E9EeqjqjDrEZTzxxfaPPqpeLDUmm9avr5oUYvuffeYGkcWSQ+fObgW5Tp1c+0O+Tb+RT2JTESVa\nuhTats19PJmoaVr12bh1OIqAvXG9qTZ7b6uq9shJhD7kQ4kjxgYFmmwpL4evvqpMDvEJYvVq9yUV\nnyA6dXLbXP6ire/WrnVjTxJLFxCNZaaDWI/je6AXrjG8mmx0xxWR/sAIXHJ6JHG9ce+Ye4DjgI3A\nBao6LckxeZM4rroK7rgDTj4ZXnop7GhMPvjhB9fWMG+ea5iOJYjPP4edd66aHGIJwhqnw7N2Lfzx\nj3DPPdXfO/ZYl0Tat899XMkEkTimqWrvjCNLdWGRItya40cDS4CPqL7m+PHAJap6vIgcCNytqgcl\nOVfeJA6w7rmmOlVYvrxy7EMsScS6tnbu7NobYo/OnV2Jor4Ojoua2bNd9fOmTdXfGzHCTUYYxW7I\nQbRx7CwifyB5iUNV9c50L5agL7AwVnIRkVHAKcDcuGNOBp7wLjhFRJqLSGtVXZHhtUP1zTdu3pkG\nDazKqr4pK3NzL8USQ/y2USPXEN21q0sOJ5zg9tu3t9JD1Kxf70oVd92V/P2BA90CS4Wa2GtKHEVA\nkB1H2wHx88Z+Axzo45jdgLxOHO3aVe5/8gn06RNeLCY4334LEydWLUV8/rlrCI31VjrkENf1sksX\n67kUZeXlLkkMH576mL/+FX7zG1d9WOhqShzLVfVPAV7b72/txBJP0s+VlJT8tF9cXExxcXGdgsoV\nVVdNtd9+VuooVIsXw9ixLkEMGOCSQ6dOVafGNtFUUeGm/Lj11tTHjBsHxx9f2c0+H5SWllJaWprx\necJs4zgIKFHV/t7z64CK+AZyEXkAKFXVUd7zecARiVVV+dbGEfOPf7iJ3Lbbzk3+ZowJx8qVcPHF\n8MILqY/55z/hvPNghx1yF1fQgmgcb6mqqzKOLNWFRRriGsePApYCH1Jz4/hBwIhCaByPF2sc//JL\n6Ngx1FCMqRcqKtyI7dqWXH35ZVcj0Lp1buIKQ9Ybx4NMGt75y0TkEuB1XHvKv1V1rogM9d5/UFVf\nEZHjRWQhsAH4VZAxhaGiwjV87rGHVVkZE4RPP4XBg91iUKkcfTQ89VRhJ4lssjXHI+COO9z4DrDk\nYUxdVVS4brFHHeU6JtTkqqvgiiuys4pePgt8WvUoy/fEAZVVVtOnu2UwjTGpbdwI991Xcy+nmDvv\nhHPPtV5ryVjiyPO/Q7Wyr36e/ynGZE1FheuZdvrp/o6fMMH1dLKBtf7UNXHYsKKIEIFXX63cN6Y+\nUXVT8IhUfRQVJU8agwfDjBkusahWPk44wf795IKVOCKmRQv4/ns3fcGHH4YdjTHZVVEBDzzgpuDw\nY+BA1w12l12Cjau+sqqqAvg7YmK/mD74AA6q1vnYmOj74QfXAJ1sVthkTjvNtUV06GAlhlyyxFEA\nf0dMfHtHWVl+jUw19cukSW7G1/Jyf8ffcgv84Q+w7bbBxmX8sTaOAiLiuhVCNGfUNPXL8uXQr1/1\n9gcRN/4hWdJ46SX3enz7gypcd50ljUJgiSOi9tkHbrzR7VvR3QTthx9cl9VkyaFtW3jzzeqfufRS\nN0tsYnJQdevN2Iy+hcuqqiKuSxe3cA9YN12TmU2b3Nxo997r/zN9+sDo0cmXPDX5z9o4CuDvSCVW\n4mjfHr7+OtxYTLRt3gzXXusWD/Jrzz3d5H428LT+sTaOAhbLiYsXu+kUTP22ZQtcc03yaqXGjZMn\njV12cb30klUrff65JQ2THitx5JFYyaNnTzc1iSlcZWXw5z/DzTf7/0zz5jB+PBx2WHBxmcJiVVUF\n8Hf4EUsebdrAsmXhxmIyowr33+9/MBy4pUjHj4cjjwwuLlN/WOIogL/Dr/heVvXoz85br77q5k9K\nx4QJbvoMY4JkbRz1SHyysK660bB0qVs/PFm7Q6qk8cQT1edaip9zyZiossSRpxKTR0VFeLHUF6rw\n5JPJk0O7dq7xOdGNN8LWrcmTw3nnWeI3+Sm0xCEiLURkooh8JiJviEjzJMe0F5G3RWS2iMwSkcvC\niDWqVOGSS9x+UVHyQVomfevXu3aHxOTQoIH7sk/UvbvrmZQsOdx8s43+N4UntDYOEfk78J2q/l1E\nrgF2UtVrE45pA7RR1eki0gyYCpwavy65d1y9auNINH++GygYU49vRVpWrIChQ930GH7ccANcfz00\naRJsXMbkSt41jovIPOAIVV3hJYhSVe1Sy2fGAv9U1UkJr9frxBETX+0xY4b7JWxgzRqXIJ591t/x\no0bBqafanEqm8OVj43hrVV3h7a8AalwmXkQ6Ar2BKcGGlb9U4a673H6PHi6R1Kd8umEDnH9+9Sqm\nnXZKnjReftkNpkusXjrzTEsaxtQk0NpXEZkItEny1g3xT1RVRSTlV5xXTfUcMExV1yc7pqSk5Kf9\n4uJiiouL6xBx/rv8cveIlT4KcTnazZth2DB48EF/x48fD8cdZ9PTG1NaWkppaWnG5wm7qqpYVZeL\nSFvg7WRVVSLSCJgAvKqqSWfgsaqq5Fatglatqr62ZQs0ahROPOkqK4Phw/3Pu/TiizYrqzHpyMeq\nqnHA+d7++cDYxANERIB/A3NSJQ2TWsuWrqSxeHHla9ts40oj116b+nO59uOPrg0isYqpUaPkSePZ\nZ5OPfzj1VEsaxuRCmCWOFsCzQAdgETBQVdeIyK7Aw6p6gogcBrwDzABigV6nqq8lnMtKHD5UVKSu\nrlmzBnbcMbhrq8Jzz7k1pP166ik45xwb62BMUPKuV1U2WeJI39dfw+6713zMwIGu+2ltM6eqwkcf\nwaOP+m93iPfGG26FOWNMblniKIC/I0xPPeVWgAvKbbfBlVdaA7UxUWKJowD+jihShcmTYeJENzZk\nyRJX5dW8OeyxBxxwgFt3urbSizEmeixxFMDfYYwxuZSPvaqMMcbkIUscxhhj0mKJwxhjTFoscRhj\njEmLJQ5jjDFpscRhjDEmLZY4jDHGpMUShzHGmLRY4jDGGJMWSxzGGGPSYonDGGNMWixxGGOMSUso\niUNEWojIRBH5TETeEJHmNRxbJCLTRGR8LmM0xhiTXFgljmuBiaraCZjkPU9lGDCHyhUA81Y2FonP\nhXyIMx9iBIsz2yzOaAgrcZwMPOHtPwGcmuwgEdkNOB54BMj7BUTz5X+mfIgzH2IEizPbLM5oCCtx\ntFbVFd7+CqB1iuPuAoYDFTmJyhhjTK0aBnViEZkItEny1g3xT1RVRaRaNZSInAisVNVpIlIcTJTG\nGGPSFcoKgCIyDyhW1eUi0hZ4W1W7JBxzC3AuUAY0BnYAnlfV85KcL+/bP4wxJgx5s3SsiPwdWKWq\nt4nItUBzVU3ZQC4iRwBXqepJOQvSGGNMUmG1cdwK9BORz4BfeM8RkV1F5OUUn7FShTHGREAoJQ5j\njDH5Ky9HjovIP0Rkroh8KiIviMiOKY7rLyLzRGSBiFyT4xjPEJHZIlIuIn1qOG6RiMzwBjl+mMsY\nvev7jTO0e+ld39eg0bDup5/7IyL3eO9/KiK9cxVbQgw1xikixSKy1rt/00Tk/0KI8VERWSEiM2s4\nJgr3ssY4I3Iv24vI296/8VkiclmK49K7n6qadw+gH9DA278VuDXJMUXAQqAj0AiYDnTNYYxdgE7A\n20CfGo77EmgR4r2sNc6w76UXw9+Bq739a5L9Nw/rfvq5P7jxSK94+wcCk0P4b+0nzmJgXK5jS4jh\n50BvYGaK90O/lz7jjMK9bAP08vabAfOz8f9mXpY4VHWiqsbGdkwBdktyWF9goaouUtWtwCjglBzG\nOE9VP/N5eGiDG33GGeq99PgaNOrJ9f30c39+il9VpwDNRSTV+KWg+P3vGOpgW1V9F/i+hkOicC/9\nxAnh38vlqjrd218PzAV2TTgs7fuZl4kjwRDglSSvtwMWxz3/xnstahR4U0Q+FpHfhB1MClG4l34H\njYZxP/3cn2THJPvBEyQ/cSpwiFdl8YqI7JOz6PyLwr30I1L3UkQ64kpIUxLeSvt+BjYAMFM1DCC8\nXlXHe8fcAGxR1ZFJjgu81d9PjD4cqqrLRGRnYKKIzPN+yWRNFuLMSQ+KTAeNegK/n0n4vT+Jvz5z\n3TPFz/U+Adqr6kYROQ4Yi6vKjJqw76UfkbmXItIMeA4Y5pU8qh2S8LzG+xnZxKGq/Wp6X0QuwNXN\nHZXikCVA+7jn7XGZNGtqi9HnOZZ5229F5EVcdUJWv+iyEGfg9xJqjtNrhGyjlYNGV6Y4R+D3Mwk/\n9yfxmN2813Kp1jhVdV3c/qsicp+ItFDV1TmK0Y8o3MtaReVeikgj4HngKVUdm+SQtO9nXlZViUh/\n3BxWp6jqjykO+xjYW0Q6isg2wJnAuFzFmCBpPaeINBGR7b39psAxQMqeJDmQqj42CvdyHHC+t38+\n7tdbFSHeTz/3ZxxwnhfbQcCauKq3XKk1ThFpLSLi7ffFddmPUtKAaNzLWkXhXnrX/zcwR1VHpDgs\n/fsZZot/Bj0FFgBfAdO8x33e67sCL8cddxyuF8FC4Locx/hLXL3hJmA58GpijMCeuJ4t04FZuY7R\nb5xh30vv+i2AN4HPgDdwsw1E5n4muz/AUGBo3DH3eu9/Sg097cKME/i9d++mA/8DDgohxmeApcAW\n7//NIRG9lzXGGZF7eRhuktjpcd+Xx2V6P20AoDHGmLTkZVWVMcaY8FjiMMYYkxZLHMYYY9JiicMY\nY0xaLHEYY4xJiyUOY4wxabHEYYwxJi2WOExGRKRl3HoDy0TkG2//exGZncF5LxC3Rkj3uNdmiUiH\nLMWdbL6erBORZ7xJ7oYlvF4iIhUi8rO41y73Xku5Lop33CIRaZHlODO+HyJS6q31cWItxx0mInNq\nWm/DRJslDpMRVV2lqr1VtTfwAHCnt98LN2I1E99QdYLDbI5WrfO5RKTI53FtgP1Vtaeq3p3k+jOB\ns+JeOwM30rg2QYzazcY5FThHVSfUeJDqe7jRyyZPWeIw2SZx2yIRecgrKbwuIo0BRORnIvKqN/X5\nOyLSOcl5FJgAdBORajOKxv9CFpEBIvKYt/+4N5ncByLyubhV2J7wfuE+lnCOO73Y3hSRVjXF5p33\nARGZDNyWcJ7GIvKYuJUHPxGRYu+tN4B2XgnssCR/41i89TC8kscaYFXcec/2zjlTRG5NerNFBovI\nFO8aD4hIA+/1/iIyVUSmi5t1OFbKuTLus0lLcCIyXEQ+9EpKJd5rTUXkZe98M0VkYLJ4iJvvzCuB\n7OfttxKRL5MdZ/KPJQ4TpL2Be1V1X9yX4une6w8Bl6rq/rjJKu9L8fkK3Mp/1yd5T1Psg5vH6mDg\nCtwEbn8HugHdRaSHd0xT4CMvtv8CN/mIbVfgYFW9KuF6vwfKVbUHcDbwhDeJ4EnA516J7L0kf8MP\nwNci0g034eDo2BsisitudcsjcaW3A0SkyqJLItIVGAgc4pXyKoBB4qaUfwg4TVV74Uoyye5TtVKG\niBwD7KWqfXFrN+wnIj8HjgWWqGovVe0OvJbk70k8pya7hsl/kZ1W3RSEL1V1hrc/FegobtbaQ4Ax\n3sShANvUcI6RwA3iFqHxQ4HYGiOzgOWqOhvAa3PpCMzAfcnGvqifAl6oJTYFxmjyyd0OBe4BUNX5\nIvIVbt0FP+0Go3HJ5hjcEgG/wv0aPwAoVdVVXuxPA4cDL3mfE+/4/YCPvXgb4yaqPBB4R1W/8mJa\n4yOOmGOAY0Rkmve8KbAX8B5wh1fymZAiEZp6whKHCdLmuP1y3BdbA+B77xdyrVS1XETuAK5NfCtu\nf7uE97Z424qEGCpI/v+8eOerLbaNNYRal6qXWHXcP3Cln3VxCSsxQUmS1wCeUNUqJbIaGqfLqFrL\n0DjFcX9T1YcSXxSR3sAJwF9EZJKq/jnF55NdL9W1TB6yqiqTS6JucZsvRWQAuPUC4qqPqhwbt/84\ncDSwc9xrK0Ski1en/0vSrxJpQGUVzjnAu2nEluhdYJD3mU5AB9zU5bURVd0EXAP8Ne51BT4EjhDX\na60I14j+34RjJgEDvKopRKSF12YxGTg8VkqTyh5Yi4A+3mt9gD2SxPQ6MMQrfSEi7URkZ3GLZ/2o\nqk8Dt8fOU4tFwP7e/gAfx5s8YYnDZFtNbQ+x54OAC0Uktm7GySnOowCquhW4m6qJ41rcr/X3cWsi\n+I0hZgPQV1yX0GLgZh+xpTrXfUADEZkBjALO92Ku6TM/vaeqo1V1epU3VJfj/sa3cWspfKyVy/zG\nPjcX+D/gDRH5FNcY30ZVvwMuwlW/TcetGwFuFbgWIjIL1y4Tn9xi55yIqx78wPt7ngW2B7oDU7wq\nrBuB2kob4BLMxSLyCdCylnth8oitx2GMyQoReRu4SlWn+ji2IzDea2g3ecZKHMaYbFkNPF5DGwsA\nXi+tccC3OYnKZJ2VOIwxxqTFShzGGGPSYonDGGNMWixxGGOMSYslDmOMMWmxxGGMMSYt/w+aFKEe\n8ItU/wAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "run_simulation(np.linspace(0, 200, 501), model=fnm, # return_type='nyaplot',\n", " opt_kwargs={'x': 'u', 'y': ['v']})" ] } ], "metadata": { "kernelspec": { "display_name": "Python [conda root]", "language": "python", "name": "conda-root-py" }, "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.5.2" } }, "nbformat": 4, "nbformat_minor": 0 }