Skip to content

Instantly share code, notes, and snippets.

@domitry
Created December 12, 2014 08:00
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save domitry/40c35d10fc3414fd0117 to your computer and use it in GitHub Desktop.
Save domitry/40c35d10fc3414fd0117 to your computer and use it in GitHub Desktop.
Huee
Display the source blob
Display the rendered blob
Raw
{
"metadata": {
"name": "",
"signature": "sha256:1e750f2d985be7eff52bef9401c0fc95bfd38e274be52b1698612243b5db9f16"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"%matplotlib inline"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 1
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from gillespie import *"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 2
},
{
"cell_type": "heading",
"level": 1,
"metadata": {},
"source": [
"Solution by Gillespie"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"A = Species('A')\n",
"B = Species('B')"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 3
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"sim = Simulator(A, B)"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 4
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"alpha1 = 10.0\n",
"alpha2 = 10.0\n",
"beta = 1.3\n",
"gamma = 1.3\n",
"c1 = 0.01\n",
"c2 = 0.01"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 5
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"sim.reaction(lambda a, b: alpha1/(1.0 + b**beta), lambda sim: sim.increment(A, 1))\n",
"sim.reaction(lambda a, b: c1*a, lambda sim: sim.decrement(A, 1))\n",
"sim.reaction(lambda a, b: alpha2/(1.0 + a**gamma), lambda sim: sim.increment(B, 1))\n",
"sim.reaction(lambda a, b: c2*b, lambda sim: sim.decrement(B, 1))"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 6
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"sim.increment(A, 0)\n",
"sim.increment(B, 100)"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 7
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"sim.run(30000)"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 8
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"sim.plot()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEPCAYAAACHuClZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd4VGXax/HvJKGF3ruGLl1BsUsAZS0IKip2sL/qClh2\nBRUJ6CrYVsW1N3AVxY6iKLgERRBUinQpBgxVCZ0QSDLvH8+ZyUwmZZIp58zk97muuXLmmTNz7lDm\nPk8HERERERERERERERERERERERERERERkRK9AewAlvuUPQGsBpYBHwO1fV4bDawD1gD9fcp7Wp+x\nDng2gvGKiIgDnQmcgH8yOQdIsI4nWA+ATsBSoBKQAqwHXNZri4Be1vGXwLkRi1hERMolofRTyu17\nYHehsllAvnW8EGhhHQ8CpgJHgQxMMjkZaArUxCQUgCnARRGLWEREyiWSyaQ0N2BqGgDNgEyf1zKB\n5kWUb7HKRUTEQexKJg8AR4B3bbq+iIiEUZIN1xwGnA/08ynbArT0ed4CUyPZQkFTmKd8S1Ef2qZN\nG/eGDRvCGqiISAWwAWhrdxClScG/A/5cYCXQoNB5ng74ykArzC/n6YBfiOk/cVFyB7w7FowdO9bu\nEIISC3HGQoxut+IMN8UZXoA7HF/2kayZTAV6YxLHH8BYzPDfypiOeIAFwO3AKmCa9TPXKvP8grcD\nbwHVMMlkZgRjFhGRcohkMrmyiLI3Sjj/UetR2C9A17BEJCIiEWHnaK4KKTU11e4QghILccZCjKA4\nw01xOpOr9FNihtX8JyIiwXK5XBCGXGDHaC4REcerV68eu3cXnncdu+rWrUtWVlbEPl81ExGRIrhc\nLuLpO6W43ydcNRP1mYiISMiUTEREJGRKJiIiEjIlExERCZmSiYhIDEtNTaVevXocOXLE1jiUTERE\nYlRGRgaLFi2iUaNGTJ8+3dZYlExERGLUlClTOPvss7n22muZPHmyrbFonomISBFiYZ5J27ZtGTdu\nHL169aJz585kZmbSqFGjIs/VPBMREYdyucLzKI958+axZcsWBg4cSLt27ejUqRPvvmvffoNKJiIi\n5eR2h+dRHpMnT6Z///7UrFkTgMsuu8zWpi41c4mIFMHJzVzZ2dk0adKE/Px8atSoAUBOTg579uxh\n6dKldOvWLeA9auYSERE/n376KUlJSaxevZply5axbNkyVq9ezZlnnsmUKVNsiUk1ExGRIji5ZnLe\neefRpUsXnnjiCb/yDz74gBEjRpCZmUlCgn9dIdI1EyUTEZEiODmZlIeauUREguAa52L5juV2h1Fh\nKZmISMzLyjabPnV7qRur/lzFyp0rbY6o4lEzl4jEvMZPNmbnwZ1+Ze6xoX0fqJmrbFQzEZGYVziR\nAPR/uz+ucfF0v+xsSiYiDnfwyEHeXPKm3WE4kmucyy9h/HDDD/x8888AzNo4C4C5GXO9r1/50ZW4\nxrlo9WwrAHJyczhw5EAUI45fSXYHICLFW/vXWo77z3EADD1+KAmugvu/fHe+3/OKJt+d7/f86f5P\nc1rL0wLOS52cyqKbFtHrtV7esow9Gbjdbq755Br+9/v/2PXPXRGPN95V3H+JIjHAk0gAEscnekcr\nfbXuKxLHJ/Ldpu/sCs12i7Ys8h6POn0Ud516V8A5D/d5GMAvkXgkjE/gw1UfejvvJTSqmYg41MLM\nhQFl3V7qRoIrwXtX3vut3iF3NMeiHQd28Nri16hTtQ6779sd8HpKnRQy9mSwLmtdwGvPnfscw2cO\nD/i8xjUaRyzeikA1E4lphdvM48kpr58CwJDOQ/zKCzfvbMjaELWYnOKVX17h9SWvs+fwniJf/+jy\nj/jiyi8Y3HGwt+ytQW/hHuvmzpPvpHPDzn7nX/3x1RGNtyJQMpGY9f6K9+0OISLcbrdfgnzv0vdw\nj3Xz0gUv+Z33wJkPANB2Uttiv1QLy9yXGZCMYtGu7JL7OHo07cEF7S9gYIeBzLhqBj/d/BNDjx/q\nfX3JrUsYlzqOK7pcAcC3v38b0XjDLSUlheTkZGrWrEm9evUYMGAAmZmZtsYUyWTyBrAD8J2SWg+Y\nBfwGfAPU8XltNLAOWAP09ynvaX3GOuDZCMYrMSAnN4cPVn4AwBUfmS+CtvXa2hlSWPyx9w9c41xM\nWzmNfTn7vOWTLypYUvzWE28l+4Fstt69FcBvzsC2/dtK/HzXOBejZo+i5b9bkjg+MczRFy/fnc+U\nZcUvPJiVncX6rPVl/tyv1n8V9LnntzufE5ud6FdWKbESD/V+iKmDp3JZp8vKfH27uVwuvvjiC/bv\n38+2bdto3Lgxd955p60xRTKZvAmcW6hsFCaZtAe+tZ4DdAKGWD/PBV6gYBLNi8CNQDvrUfgzpQK5\n5YtbuPzDyzl09JC3bH3W+pieXOZ2uznmmWMAGPLhEOpMNPdYHRt05Lru1/mdWzWpKk1rNuXlAS/z\nj9P/4S0//93z/c7bun+rtwnwcO5hACb+MDGSvwYnvXoSzy18zq9sYeZChn46tNiEUf/x+rSb1K7M\n1/pt12/m599/K3ughVze+XLA1NpiUZUqVRg8eDCrVq2yNY5IJpPvgcI9YwMBz63WZOAi63gQMBU4\nCmQA64GTgaZATcAzbGOKz3ukglnz1xrvXW71R6sDsOK2FQB8ue5L2+LytT9nP65xLn7Z+ktQ598x\n4w4Sxhf933DVHcV/OdzS8xbqVK3j7XzP2JPBtv3bWJ+1Htc4F82fbu49t8ajNQLe79uM9mPmj0HF\n6jFy5kiGfDgkoPznrT8zYuYIrvn4Go7kHQHgtDfMUN01f60JOP/gkYOlXuv33b+z48AOvzLPjcO0\nS6fRrn7ZE1FhjaqbbW5b/rtlyJ8VTZ4/h0OHDvH+++9z6qmn2hpPpHsuU4DPga7W891AXZ9rZ1nP\nJwE/Au9Yr70GfIVJLBOAc6zyM4F/AhcWca0KsZyK50ugIozgWfXnKjq/0Jl3LnmHlrVactZbZwWc\nkzsml6SHzaBEJ/yZ+H5JBxOP7/nLb1tO1xfNf5VXBrzCzT1vDuqabyx5gxun3xjUub2P7c3cTWYS\n37Djh9GrWS9u//J23hz0JsOOH1bq+zdkbaDtJNOs6Pv7ud3ugKToHusOGBzhec83G77hb//9m99r\n+Q/le5b28HKNc3Fs7WPJGJnBXTPv4pmFz3hfO/LgESolVio15mD4xvnaha9xY48bg1pOJVyDP8r6\nbzclJYVdu3aRlJTEwYMHadSoETNnzqRLly7FvifSy6nYOTTYbT0kSLn5uX7PJ8ybQNdGXbmg/QU2\nRRRZY+aMAUoeaZOYUND+v/rP1XRs2DHicZVVcTcAbyx5w+95l0Zd2HPfHupMrENqSmrQn3/DCTcU\nmUxa1Grh13Sz7s511K1al/Fzx/Pcoud4a+lbvLX0LQCu/+x6rul2DUkJgV8Jc36fQ98pfYv8vZIr\nJdMguQFP9X8q4PWivrjOe+c8rupyFdd9el3Aa71e68VPN//kff7qL68CsGnvpiK/tMOVSAAO3n/Q\nW9u96fObuLFHcMnZrhsYl8vFZ599Rt++fXG73Xz66af07t2bVatW0bixPUOcoz2aawfQxDpuCngW\n1NkC+NYxWwCZVnmLQuVbivvwtLQ07yM9PT1cMTvGJ6s/8R67xrkY/e1oBkwdEFRzQSz6ePXHAWWr\n71hNzoM5AFzb7VqgoN280wudmL1xtu1DhW878TbObn023278lm4vFmyfOnLmSLKys1ixcwWHcw/7\nJYCjY44CULtqbdxj3WVuvtlxr2kKmnHVDHIezKFKYhWW3rqUA6MLlgppW68t9ZPr8+x5RY9jGf7V\ncI7mHfUrc7vd/N+M/yv2uoeOHmLz3s1c9oHpxM4dk8trF74GQM9XegKw896CdbNmrp/pl0hmXTuL\nxbcsBkwzma9bvril5F86jJIrJXuPT2lxSsCNm5O5XC4uvvhiEhMT+eGHH0o9Pz093e+7MmxxhO2T\nipaCfzPX48AuYCKm872O9bMT8C7QC2gOzAbaYmouC4HhmH6TGcBzwMwirhXXzVw9Xu7Bku1Linzt\n9hNvp2vjrlzX/TqSKyXz7I/PMvzk4QFNBrHkaN5RKj9SmTpV63iHvR68/6Dff3pfRSWQDvU7sObv\naxg8bTAfr/6Y5bctp0uj4psBQlXtX9U4nHuYFbetoMuLwV+ne+PuLP2/pRGLa/uB7TSq3shv6ZUD\nRw5Q87GapKak8v6l79P4Sf+72fGp4xnTewxJ45PIc+f5vfbe4Pe8I+kK89ypF27u27p/Ky5cNHu6\nWUEMow9QvbKpDbR6thUZezK8TV2e99/c42Y2793MgswF7MvZR95Ded7RaJGoFSzZtoQer/Rg1Omj\nmHDOBMcO7GjVqhWvvfYa/fr1w+12M336dC699FJ+/fVXOnYsunYey6sGTwXmAx2AP4DrKej/+A3o\naz0HWAVMs35+BdxOQRPY7Zg+lHWYjvmiEkncKyqRzL52NgDNazXnthm3MWnhJFbsXMHIr0cWOfPX\nafLd+Tzz4zNFznuo/EhlwCyH4R7rxj3WXWwiKc7aXWuZt3met4bj6Y8IpyqPVME1zsVtX9zmHTV1\nbJ1jg37/huEbIppIAJrUaBKwhleNyjVwj3UzZ+gcGlVvRKs6rfxefyj9ITbu3uhNJKkpqZzX9jwA\nhnQZws57d+Kyvn9u7XkrAC8PeLnYGJrVbEbTmk29z91j3d5EAjD/hvkAPDn/Sb8vvBcueIGZ18xk\n76i9uMe6SXAlkP1ANocfOFzmP4dgdGjQAYAJP0wo5Uz7XXjhhdSsWZPatWszZswYpkyZUmwiiYbY\nvXUNFLc1k5zcHKr+qypgmkM88xDqVatH3Yl1vXfuXRt1ZflOM61n3vXzOP2Y0+0JOAh5+XnejvOp\ng6d6J495eO5M943aR80qNUv9vHmb53Hmm2cCpjmnuKGo4bqbLa4p7Z5T7+HJ/k/6vf7qha8ysMPA\ngLv/89udz4yrZoQlnlAdOHKAi9+/mNkbZwe89vrA17n++OtLrOm6xrn45ZZf6NG0B2BGYTVIbkBW\ndpZfcv3r0F/8tuu3IhdkLPxn+vKAl7mlZ/Saujw6/acTq/9aDWlF9/vEKu0BH7y4TSa/7viVE14+\ngbyH8gJee23xa9z8eeCon76t+vLtdc6d1fv0gqe555t7vM/3j95PjcpmCOuCPxZw2huneUfVlNX2\nA9tp+lRTv7I5Q+cwZs4Yvr/++9ACp+QRPLOvnU2/1v3IPprN0fyj1KpSy/vaq7+8SvNazWlSownT\nVk7jvtPvo261usV+lh3+2PsHo78dzTvL3ykou+sPWtRqUcK7YP4f8zm1xakhNa3e8NkNvLm0YKl9\nO0fnuca5lEzK+vmhfoCDxFUycbvd7MrexTu/vsPIr0easiL+c+W780kcn0i1pGpk52b7f4YDhsp6\nHMk7Qk5ujreWUdQX8uo7VnNM7WO8o2pCiX/q8qlc9fFVgOmgd7lctJvUDvdYN//99b9syNrAkC5D\nOK7BcaV8UoHP137OwPcGFvu6bx9ArPP8/ewdtdcvIUbSX4f+ouETDb3P7fz3+7/f/0e/1v2UTMpA\na3M51KIti2j4RENvIilOgisB91g3daqaWdPzb5jP9Cumez/DCY7mHaXKI1WoNaHgS6lhckOe6v8U\ny28rWG2n4386cvkHl4flmld2vRKASzpeQrv67Whe00zic41zce0n15I2N42O//FvX35q/lNFLv2x\nL2cfWdlZ/GPWPwJe6964O4/1ewwgbhIJmDkf7rHuqCUSgAbJDbzHZRkaHQl9WwUOhZaSaQl6h3pq\ngf+4/dKarNYPX8/PW3/m1Janeu8+XvjpBXo1D9zHIVo2793Mip0ruODdgnkwT81/intOu4c/D/1J\nl0Zd6NKoC5tHbvYuJzJjnelDePGCF0O+vu+dbbVK1Yo8Z+/hvXy/+XsGtB/AvbPuBfBbwuTAkQPU\nnlAbgEf7Psr9/7sf8B9Z5na76deqX8jxOoldIwGdVJuWslEzl0P5NgNlP5BN1aSqZXp/2+fasmH3\nBvvbnQu5pOMlTDx7Iu0mtfP7vSo9XMlvbH8k4vaNp0blGhw4coBjah/D5r2b2Tdqn7fmdFyD45h4\n9kT6tupLzccCO/9X3r6STg07hT0+cZZgZsDHEjVzVTAHjhzwfun1bNoT91h3mRMJwNsXvw0QMAnN\nLp541v611ruwn+/vteufu3j23GfJGJERsQToaVLbee9OfrzRrEe1ee9mAL8muDV/rWHQe4OKTCQd\n6ndQIhEpgmomDnPL57fw6mKzjETeQ3kh7fEdyXW89hze4+2nAdPUszdnL3Wq1qHnKz3pk9KHpxY8\nxVnHnsV3m77DPdZN38l9mZMxp+A9NjdplGemfM6DOVROrByBaMRp6tWrx+7dgbs4xqq6deuSlRW4\nRbFqJnHKk0iAkBIJwPFNjg/5Mwr7MfNHXONc1J1Y1zvZMN+dT8L4BOpOrMv+nP0s3rbY2+czd9hc\nb9K459SCocD3n3F/WOMqjwU3Lggoy7wrky+u/MKvbPPIzTx45oMASiQVSFZWFm6325YHaWa+i+f5\npIWTIA2/R7vn2vm9Z9aGWZAG36z/ht5v9g74zKISSTgpmTjI1+u/9h6H4679rlPuIt+d77cceaiG\nflqwW137Se0BuP/bgsTg21xUmO/CfP/q96+wxVRep7Q4hTlD53iTyoXtL6R5reZc0P4C3hv8HmCW\nuG9ZuyUP933Y9pqUVCyr/lxF0vgk3G43d34VuPHVuqx1fttWe7Zv7v/f/t6VoaNJycRBzn3H7PtV\n0rIUZeFZCHHr/q3lWgwy353PK7+8Qo+Xe+Aa56L1s605o+UZ3hpGw+pmTsCvO34t8v1vDPRfFbd/\nm/68f+n77B21t8yxREpqSionNDkBgEnnTfKWX3TcRXRt1JXOjToX91aRiPEM789z5zF+7viA17s2\n8l8ayDXOVeKCnNGgPhOH8G2/P3T/oWKHspZV/cfrk5WdxYbhG2hdt3WZ3nv1x1fz7vJ3A8rb1mvL\n0O5DGTNnDI/1e4zR347mko6XsPPgTuZtngfY3x8iEst8943x2Hr3VrYd2MZLP7/EKxe+wkmvnhSw\n0rJHWUaAagZ8oJhNJnfMuIMXfn4BCP+XcFZ2Fu0mtePzKz8vcj0kX/tz9nubqXw3nSrsrUFvMfGH\niWb9IstHl3/EwA4DmfHbDOon1+eMY84I3y8hUgEVt7FYSeel9U5j+MnDy7RMj5JJIEcnE9+FDXf9\ncxc1Ktdg8bbFtKnbhkZPmm1D99y3h9pVa4f92oOnDSb12FTuPDmw3dXjzi/v5Pmfng8oP67Bcay6\nfRWHcw+T/KiZpHf4gcNk7sv0u3Mqapc8EQlNXn4e3236jh5NexT73XDo6CEqJ1YmLz+PyomVy/z/\nUMkkkGOTycqdKxk8bTBrd631ll3d9Wq/xfQSXAlFLuQYDp47l8MPHKZKUpUSzynM926o8FDjsm5R\nKyLOo6HBMSL7aDZdXuzil0gAv0QCsPbv/q+H09/amL22V/25qsjX8/JNEqtdpTYrblvhLb/xBP8V\nexfetJDsBwoWkxxz1hgW3LhAiUREVDOJlINHDlLjsRoM7zWc5xY9B5jayDmtz2HYZ8MCzo/0F/KV\nH13JeW3P81t3yqNwjWPP4T1mn5QINbuJiHOEq2ZS1oUe62H2YS96LKh4eRYs9CSSYccP481BZq+G\n3Yd3k5ufS5+UPtRPrh+V9X96NOnB4m2LA5JJVnbgRKY6VeuotiEiZRJMMpkLXGid+wvwJ/ADcFcE\n44p5Qz4c4j0ecfIInjn3Ge/zkaeUvKx8JHRv0p2v1n8VUD7nd7O8ycrbV0Y7JBGJI8H0mdQG9gGX\nAFOAXsDZkQwq1hXe5W9gh+I3VIqWBskN+OvQX35lbrebSz+4FECLF4pISIJJJolAU+BywLNhtdpA\nSrD9wHYAruxiNmiye6MfMEN812Wt81vm3bN5VlKCtrURkdAEk0zGA18DG4BFQBtgXSSDige39ryV\ndwe/i3usO+yLLZZHcqVkmtdszrpd5q9uYeZCTnn9FACOjnHGMvUiErtKuyVNBFoC3XzKNgCDIxZR\njOv/dn8AXrjgBZsjCdS+fnu+XPclHRt2ZH3WekCr4IpIeJR2y5wHXBmNQOLBtv3bmLVxFhD68vGR\n8NX6r7xb0y7IXECdqnXIeTDH5qhEJB4E8403D3geOBPo4fOQQto/b5Zkt3Pf9ZJMuWgKAM8vep7/\n/PQf9hzeY3NEIhIvgpmokk7RHe59whtKyGyftNjmuTZs3L2RHffuoFH1RrbGUhzfJVAW3rTQsYlP\nRKIjmpMWU0O9SEVwNO8oG3dvBHBsIgG4ossVvLfCbPykRCIi4RLsmNABQCfAd4H8wB1bKrAH/vcA\nAIM6DLI5kpK9c8k77MvZR8taLe0ORUTiSDDJ5GWgGtAXeBW4DFgYyaBikefL+ZMhn9gcSckSXAnM\nuGpG6SeKiJRBMO1ky4GumPW4ugE1gJmA03Y/sq3PxHdTKa1pJSKxJJpL0HvWHD8ENAdygSYhXnc0\nsBKTqN4FqmAWkZwF/AZ8A9QpdP46YA3QP8Rrh1W+O9+bSNrXb29zNCIi9ggmmXwB1AWewCz0mAFM\nDeGaKcDNmOHFXTETI68ARmGSSXvgW+s5mL6aIdbPc4EXgow7Kub/Md97HMk9SUREnCzY5VR2Ax9h\nEsFxwJgQrrkPOAokY/pskoGtwEBgsnXOZOAi63gQJnkdxSSy9ZjFJm23btc6znzzTAA2j9xsczQi\nIvYpqQN+MAXzS1wEzjX5uJzXzAKeAjZjmtC+xtRIGgM7rHN2WM8BmgE/+rw/E9PcZjvPJEWAlrU1\nOkpEKq6SksmFlLw6cHmTSRtgJKaWsxf4ALim0DnuUq5d5GtpaWne49TUVFJTU8sZYunOefsc7/Hb\nF78dseuIiIRTeno66enpYf9cO7btHQKcA9xkPb8WOAUz9LgPsB2z5P0cTJOap+9kgvVzJjCWwOHJ\nUR3N5ZlJ/tBZDzGuz7ioXVdEJJyiOQN+LKYmULipq7yTFtdg+lyqAYcxG20tAg4CQ4GJ1s9PrfOn\nY0Z8PY1p3mpnnW+bu7++23s8NnWsjZGIiDhDMMnkIAVJpBpmNvyqEK65DLNj489APrAYeAWoCUwD\nbsR0tF9unb/KKl+FGZZ8OzZvzvX6ktcBzSkREfEoT9WmCmYeSO8wxxKqqDVzuca5OPOYM/nu+u+i\ncj0RkUiJ5qTFwqrjkNFUdvAkrPtOv8/mSEREnCOYZq7lPscJQCMq8CKPL/78IgBnHOO01WREROwT\nTNUmxec4FzMHxImbhkelmcszikv9JSISD6LZzNUEM9EwAzNhsBpwcqgXjmWNqzcu/SQRkQokmGTy\nEnDA5/lBq6zCydyXCcD3139vcyQiIs4SbAd8vs9xHmZxxgqn5b/Nkimt6rayORIREWcJJpn8DgwH\nKgGVgRHAxkgG5UR5+Xne46SEYDeoFBGpGIJJJv8HnA5swfSZnALcEsmgnOjPQ38C8HCfh22ORETE\neYK5xd6BWU+rQtu2fxvdG3fnwbMetDsUERHHCaZm0gGzWdVK63k3oMJ9o/Z4pQfLdiyzOwwREUcK\nJpm8CtwPHLGeLweujFhEDpSVnWV3CCIijhZMMknGf7l3N86ctBgR+e586j9eH4BNIzfZHI2IiDMF\nk0z+BNr6PL8U2BaZcJxn+Y6C1WSOqX2MjZGIiDhXMB3wf8csEX8cZq/234GrIxmUk4yYOQKAP+76\nw+ZIREScqyzrsVTH1GT2RyiWUEVkbS6txSUi8SwaOy3e43Ps+03q2XHx6VAvHiuqV6pudwgiIo5W\nUjKpSdE7Ghbevjfubbl7i90hiIg4WknJJC1aQTjVj5k/AlCrSi2bIxERcbZgRnO1BD7BjOr6E/gI\naBHJoJwiY08G4G1TFBGRYgSTTN4EpgPNrMfnVlnc+2j1R3aHICISE4K55V4GdA+izG5hH82lkVwi\nEu+iudPiLuBazB4mScA1wF+hXtjpRs8eDcDsa2fbHImIiPMFuwf8JMzS8wDzgTuBzRGKqbzCWjPx\n1EqOjjmq/UtEJG5FY56JRwZwYagXiiUHjxwE4Kxjz1IiEREJQknflJMw80mKylhuzO6LcanP5D4A\nzB021+ZIRERiQ0nJ5P+AFcA0zJpcUJBY4rpH+qetP9kdgohITCkpmTQFLgMuB/KA94EPgD1RiMs2\nby19C4Cbe9xsbyAiIjEk2E6XFsAVwN3AfcDbEYuo/MLSAa/hwCJSkURzaHBPYARmSPBXwC+hXhSo\nA3wIrAZWAScD9YBZwG/AN9Y5HqOBdcAaoH8Yri8iImFUUjJ5GJM47gLmAicBN2K+/EP1LPAl0BGz\np/waYBQmmbTH7Dk/yjq3EzDE+nku8EIpcZfbkTyzM3HumNxIfLyISNwqqWqTj9kI61ARr7kxSaA8\nagNLgNaFytcAvYEdQBMgHbMh12grlonWeTMxi1D+WDimUJu5Nu3ZRMqzKWriEpEKIxrzTAp/2YdL\nK8yCkW9ilmT5BRgJNMYkEqyfja3jZvgnjkygeSQCS5ubFomPFRGJeyUlk4wIXrMHZjvgn4BnKGjS\n8nBT8vDjIl9LS0vzHqemppKamlqmwLo16saIk0eU6T0iIrEkPT2d9PT0sH+uHWurNwEWYGooAGdg\nmrJaA32A7ZhhyXMwzVyeRDPB+jkTGAssLPS5ITdzVXmkCjf3uJnnz38+pM8REYkV0RzNFW7bgT8w\nHe0AZwMrMUvbD7XKhgKfWsfTMcOSK2MSUDtgUbiDysvP40jeEZZuXxrujxYRiXslNXN9C/QDHgf+\nGebr3gm8g0kQG4DrMasST8OMGMvATJYEM3psmvUzF7idCMzA37LfbM374gUvhvujRUTiXmkz4E8D\nBgLvEbj3++IQrrsMM9S4sLOLOf9R6xExbyx5A4CujbtG8jIiInGppHayyzC1hNOBn4t4vU9EIiq/\nkPpMNPNdRCqiaAwN/sB6PASMD/VCTtexQUdqV61tdxgiIjEpmM06xgODgLMwzVxzMZ3lcaVd/XYM\n7T609BPpBFlAAAAQHUlEQVRFRCRAMKO5JmD2LlmJWUtrOPBYJIOKtrz8PKavnU7rupGapykiEt+C\nqZlcAByPWYYe4C1gKWZuSFy495t7AWhfv30pZ4qISFGCqZm48V/Btw5xtjlW7aq1aVS9EcmVku0O\nRUQkJgVTM3kMMwx4DqbHvzeBy5/EtN3Zuxl1elz9SiIiURVMzWQqcCrwCfCRdfxeJIOKtucWPUdi\nQqLdYYiIxKxgaiZg9oD/LJKB2GVD1gYAOjXsZHMkIiKxy461uRxl6/6tAPRr1c/mSEREYleFTyZb\n9m/hsk6XeWaBiohIOZSWTJKAtdEIxC5zM+ZSNamq3WGIiMS0YG7HP8NMVNwU4VhCVa61uVzjXFRN\nqkr2A9kRCElExNmisTaXRz3M7PdFwEGrzI1ZTTguXN758tJPEhGRYgWTTMYUURYXkxb35ewDYEK/\nCaWcKSIiJQkmmaQDKUBbYDaQHOT7HO/ln18GoGnNpjZHIiIS24IZzXULZin6l63nLTATGGNeg+QG\nDOowyO4wRERiXjDJ5A7gDGCf9fw3oFHEIoqiTXs30a1xN7vDEBGJecEkkxzr4ZFEnPSZZOzJ4Nja\nx9odhohIzAsmmcwFHsD0lZyDafKKi82xMvZkkFInxe4wRERiXjDJZBTwJ7AcuBX4EngwkkFFi5KJ\niEh4BDMqKw+YDCzENG+tIQ6auXLzc9l2YBsta7e0OxQRkZgX7E6LLwEbreetKaihxKwt+7bQMLkh\nlRMr2x2KiEjMCyaZPA30AdZbz9tgEklMJ5Ol25dSvXJ1u8MQEYkLwSSTfRQkEjA1lH3FnBszLnr/\nIrtDEBGJGyUlk8HWz58xtZBp1vPLrLKYNqD9APqk9LE7DBGRuFBSMrmQgo72nZi938GM7Ir5Ndt3\nHdpFr+a97A5DRCQulJRMhkUriGjLPprNgswFtKrTyu5QRETiQjDzTFoD/8asx/W59ZgehmsnAkso\nmABZD5iFWa7lG6COz7mjgXWYYcn9Q73wzPUzAWhWs1moHyUiIgTXAf8p8BrmSz/fKgvHPJMRwCqg\npvV8FCaZPA7cZz0fBXQChlg/m2NWLm7vE0uZ7T68mz4pfbRVr4hImASTTA4Dz4X5ui2A84F/AXdb\nZQMp6JeZjFn6fhQwCJgKHAUyMCPLegE/lvfiN06/sbxvFRGRIgSTTCYBacDX+C/4uDiE6/4b+AdQ\ny6esMbDDOt5hPQdohn/iyMTUUEKS1jst1I8QERFLMMmkM3AtZuKib9NSecfVDsCMDlsCpBZzjpuS\nm9KKfC0tLc17nJqaSmpq4MdvP7AdgKu6XlV6pCIicSY9PZ309PSwf24wnQYbgI7AkTBd81FMcsrF\nDDGuBXwMnIRJLtuBpsAc4DhMUxeAZ2/dmcBYzFphvtxud+ldOa5x5ld2j4355cVEREJm9R2H3IEc\nzGiu5UDdUC/k436gJdAKuAL4Hya5TAeGWucMxXT8Y5VfAVS23tMOWBTGeEREJETBNHPVxQzJ/YmC\nPhM3psM8HDxVhAmYWfY3YjraL7fKV1nlqzC1mdsp52iynFwT/uEHDpc7WBERCRRM1Sa1mPL08IUR\nFqU2c72++HVu+vwmNXGJiFjC1cwVTM0kPdSLOMW+nH0kV0q2OwwRkbgTTDY6QEGzUmWgklVWq9h3\n2KPUmok630VE/EWzZlLD5zgB01dySqgXtsvfT/q73SGIiMSdYEZz+crHjLI6NwKxRFzruq35ey8l\nExGRcAumZjLY5zgB6AlkRyacyNm6fysbd2+kdd3WdociIhJ3gkkmvvua5GKG7Q6KVECRMm/zPAAq\nJVayORIRkfgTTDIZFukgomF/zn6GHT/M7jBEROJSSclkbDHlnlrK+DDHElEZezJIqZ1idxgiInGp\npA74g5ghwL4PN2aG+n2RDy28Hvn+EQ4ePWh3GCIicamkmsmTPse1gOHA9cB7wFORDCrcsrKzAOjf\nJuRNGkVEpAil9ZnUB+4CrgamAD2A3ZEOKtxOevUkAM5ufbbNkYiIxKfSaiYXA68A3YD9UYkoAjbu\n3mh3CCIica2kPpO7MTsaPghsxSQTz2Nf5EMLrxuOv8HuEERE4lZJNZOyzo53rARXAuP7xNTgMxGR\nmBI3CaM467PWk+/Op0mNJnaHIiISt+I+mcz/Yz4AiQmJNkciIhK/4j6ZrNu1jofOesjuMERE4lrc\nJ5M/D/1J05pN7Q5DRCSuxX0y2Zuzl5qVa9odhohIXIv7ZPLeiveoWUXJREQkkuI6mczNmAvA39r8\nzeZIRETiW1wnk7W71gJQJamKzZGIiMS3uE4mi7ct5ppu19gdhohI3HPZHUAYud1ut1+Ba5z59dxj\n3UWdLyJS4blcLghDLojrmgnA/Bvm2x2CiEjci9tkkpObQ5XEKpzY7ES7QxERiXtxm0w27N7AsXWO\npVJiJbtDERGJe6VtjhWTdhzYQecXOtsdhohIhWFHzaQlMAdYCazAbAcMUA+YBfwGfAPU8XnPaGAd\nsAYoce9dt9vNzPUzwxyyiIiUxI7RXE2sx1KgBvALcBFmf/m/gMeB+4C6wCigE/AucBJms67ZQHsg\nv9Dnuvfn7KfmY2a2+4D2A/j8ys8j/buIiMS0WB7NtR2TSAAOAKsxSWIgMNkqn4xJMACDgKnAUSAD\nWA/0KuqDV+5c6T0+eORgeKMWEZFi2d0BnwKcACwEGgM7rPId1nOAZkCmz3syMcknwK87fuW67tex\naeQmvr3u24gELCIigezsgK8BfASMwOwr78ttPYpT5Gtr/lpD54adOab2MeGJUEREgmJXMqmESSRv\nA59aZTswfSnbgabATqt8C6bT3qOFVRbg3UnvkpqSStqsNFJTU0lNTY1A6CIisSs9PZ309PSwf64d\nHfAuTJ/ILuAun/LHrbKJmI73Ovh3wPeioAO+LYG1EzdpsPCmhfRqXmSXioiIFBKuDng7aianA9cA\nvwJLrLLRwARgGnAjpqP9cuu1VVb5KiAXuJ1imrmqJVWjS6MukYpbRESKEVcLPZIG+Q/lezKtiIiU\nIpaHBkdMcqVkJRIRERvEVTLpUL+D3SGIiFRIcZVMGlZvaHcIIiIVUnwlk2QlExERO8RVMmmQ3MDu\nEEREKqS4SiaqmYiI2CO+kon6TEREbBFXyURrcomI2COukknnhtpdUUTEDvE0w8+dn6/Z7yIiZaEZ\n8EVQIhERsUdcJRMREbGHkomIiIRMyUREREKmZCIiIiFTMhERkZApmYiISMiUTEREJGRKJiIiEjIl\nExERCZmSiYiIhEzJREREQqZkIiIiIVMyERGRkCmZiIhIyJRMREQkZEomIiISMiUTEREJWSwlk3OB\nNcA64D6bYxERER+xkkwSgecxCaUTcCXQ0daIyik9Pd3uEIISC3HGQoygOMNNcTpTrCSTXsB6IAM4\nCrwHDLIzoPKKlX9gsRBnLMQIijPcFKczxUoyaQ784fM80yoTEREHSLI7gCC5gzmpd29YsQKSkyE3\nF5o1Mw8nWbsWfvkl8tfZtg0yM6FSpcDXEhKgc2dITCz+/dGKMxSxECMoznBTnMG77DK47rroXMsV\nncuE7BQgDdNnAjAayAcm+pyzHmgT3bBERGLeBqCt3UFESxLmF04BKgNLidEOeBERsdd5wFpMDWS0\nzbGIiIiIiIj4c9JkxpbAHGAlsAIYbpXXA2YBvwHfAHV83jMaE/saoH/UIjUSgSXA59ZzJ8ZZB/gQ\nWA2sAk52YJyjMX/ny4F3gSoOifENYIcVl0d54uppfcY64NkoxfkE5u98GfAxUNuhcXrcg+nHrefg\nOO/E/JmuwL+/2a44HSUR0+yVAlTC/r6UJsDx1nENTLNcR+Bx4J9W+X3ABOu4EybmSpjfYT3RHa59\nN/AOMN167sQ4JwM3WMdJmC8VJ8WZAmzEJBCA94GhDonxTOAE/L9UyhKXZ4DOIsxcL4AvKRgIE8k4\nz6Hgz2WCg+MEcxM5E/idgmTitDj7YG4iPOM7GzogTkc5FfOX6DHKejjFp8DZmIzf2CprYj0Hc0fg\nW5uaiRm5Fg0tgNmYf2SemonT4qyN+aIuzElx1sPcNNTFJLvPMV+ETokxBf8vlbLG1RRzN+txBfBS\nFOL0dTHwX+vYiXF+AHTDP5k4Lc5pQN8izgtbnLEyabE4Tp7MmIK5O1iI+c+7wyrfQcF/5maYmD2i\nGf+/gX9gquYeTouzFfAn8CawGHgVqO6wOLOAp4DNwFZgD+YO0Ekx+iprXIXLtxD9/2M3YO6MKSIe\nu+McZF3310LlTouzHXAW8COQDpxolYctzlhPJkFNZrRBDeAjYASwv9BrbkqOOxq/0wBgJ6a/pLi5\nRk6IMwnoAbxg/TxIYM3T7jjbACMxNw/NMH/31xQRg91/lsVd16n/hzweAI5g+qKcJhm4HxjrU+bU\nuXtJmNrzKZibyGnhvkCsJ5MtmPZKj5b4Z1M7VMIkkrcxzVxg7gCbWMdNMV/kEBh/C6ss0k4DBmKq\n5VMx1d+3HRhnpvX4yXr+ISapbMc5cZ4IzAd2AbmYzuJTHRajr7L8HWda5S0KlUcr3mHA+cDVPmVO\nirMN5iZiGeb/UgvgF0xtz0lxYl37Y+v4J0yLRAOcF6dtnDaZ0QVMwTQh+XqcgnbJUQR2JlbGNOls\nIPp3Nr0p6DNxYpzfAe2t4zQrRifF2R0zOqaada3JwB0OijGFwA74ssa1EDOKzkXkOmILx3kuZoRc\ng0LnOS1OX0V1wDslzluBcdZxe0yzrBPidBQnTWY8A5Pxl2KakJZg/gLqYTq7ixqOeT8m9jXA36IZ\nrKU3BaO5nBhnd8ydlO8QUafF+U8KhgZPxtROnRDjVEw/zhFM3+L15YzLM0R0PfBcFOK8ATMcdRMF\n/49ecFCcORT8efraiP/QYCfFWQnT+rAcU3tKdUCcIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiKxpD4F\n8xy2YWYGL8EskfN8GK/zJP7j/gsbDlwbxuuJiIhNxmKW7A+3mpilvkM9R8RWsb42l0g0eZaZSKVg\nCZo0zKz374AM4BJMTeNX4CvMkj9gZhOnAz9jlvn2rI81CDMj3WMCZjb9MswGUWBqQruAzmH7TUTC\nTMlEJHStMPvCDMTsuzELs79FNnABZimLScBgzMKQbwL/st57BibBgGlWuwiTNLoDj/hcYxFmCXER\nR0oq/RQRKYEbUwPJwyz4mAB8bb22HLPgXntMgvDUQBIxaycBHIPpjwHYCxwGXge+sB4eW4HWkfgF\nRMJByUQkdEesn/nAUZ/yfMz/MRem6eq0Yt7vaSHIxWyT2g+4FPi7dYz1GU7fe0QqMDVziYQmmOXj\n12L23PZsz1sJs/Q3mJVxPf0n1TGr+H6F6ezv7vMZTTF9MiKOpGQiEjy3z8+ijiGw9uDG1FYuBSZS\nsD3Bqdbr8yjYQrUWpmN/GfA9cJfP5/SyykRERALUoGAnyeLUCuIcEVsl2h2ASAV3BGiH6S/JKOac\nWzHDin+NTkgiIiIiIiIiIiIiIiIiIiIiIiIiIiIiUkH9P76rVG/OIOFQAAAAAElFTkSuQmCC\n",
"text": [
"<matplotlib.figure.Figure at 0x7efede2f60d0>"
]
}
],
"prompt_number": 9
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from scipy.integrate import odeint\n",
"from matplotlib.pyplot import *"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 10
},
{
"cell_type": "heading",
"level": 1,
"metadata": {},
"source": [
"Solution by scipy.integrate.odeint"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"def dsolve(init, t, dt=0.1):\n",
" def integrand(y, t):\n",
" du_dt = alpha1/(1.0 + y[1]**beta) - c1*y[0]\n",
" dv_dt = alpha2/(1.0 + y[0]**gamma) - c2*y[1]\n",
" return np.array([du_dt, dv_dt])\n",
" t = np.arange(0, t, dt)\n",
" y = odeint(integrand, init, t)\n",
" plot(y, xdata=t)\n",
" \n",
"dsolve([0, 100], 3000)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEACAYAAACznAEdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAF3dJREFUeJzt3XuQlPWd7/H3MNzkOiIww01BwcVhWYNsKSbHdbKHVbwc\ntRIXY8xGjSa7IZXN7tmzGzm1pexWJUetuDkaN7qGJAfXAyljjDGGEPDS0SQiQQcBBwQGxzhchvsd\ngbnsH88D0w4MzPR093OZ96uqq5/+zdPd35/P2B9+v9/z9IAkSZIkSZIkSZIkSZIkSZLUZT8AGoBV\nWW1DgCXAOmAxUJb1s9nAemAtcFVW+9TwNdYDDxewXklSgVwBTOGjgfAg8E/h9teB+8PtSmAF0AsY\nC2wASsKfLQMuDbcXAjMKVrEkqWDG8tFAWAuUh9sV4WMIRgdfz9pvETANGAGsyWr/DPB4IQqVJOWu\nRw7PKSeYRiK8Px4OI4H6rP3qgVGnaN8UtkuSYiSXQMjWEt4kSQnXM4fnNBBMFW0lmA7aFrZvAsZk\n7TeaYGSwKdzObt90qhe+4IILWmpra3MoSZK6tVpgfFdfJJcRwvPA7eH27cBzWe2fAXoD44AJBIvJ\nW4F9wGUEi8x/lfWcj6itraWlpSW1t/vuuy/yGuyb/bN/6bsBF+TwWX6SM40QFgBXAkOBD4B7Cc4q\nehq4C6gDZob71oTtNUAjMIvW6aRZwP8DziI4y2hRPoqXJOXPmQLh1nbap7fT/s3w1tabwOSOFiVJ\nKr6uLiqrE6qqqqIuoWDS3Dewf0mX9v7lS8mZdymqlnA+TJLUQSUlJZCHz3NHCJIkwECQJIUMBEkS\nYCBIkkIGgiQJMBAkSSEDQZIEGAiSpFAu33YqpVr2xZEtWd/ufrz9VG1t26UkMhC6ieaWZvYf2c+e\nD/ew98heDh07dOJ2+Njhjzw+dOwQhxsP82HjhxxrOsax5mOt99nbp7hvbG6kuaU5rzfo2IdwRz+w\nO/shXpJ1AWh4Regp29q2S0kTt99ev7qigxqbG9m8fzNb9m9h64GtbD2wlYaDDSe2dx3exZ4P95y4\n7T+6n/69+lPWt4zBfQfTr1e/k25n9TzrI9t9e/alV2kvevXo1eH70h6llJaUUtqjlB4lPbp8K6Gk\nUx/CHd33TM+XkiRfX10Rt99+AyFLY3MjtbtqWb1tNe/ufJeNuzfy3p73eG/3e2zav4lh/YYxcuBI\nygeUU9G/gooBFZQPKKe8fzlD+w2lrG/ZidugPoMo7VEadZckFYCBkDLHmo6xsmElS+uXsmzzMlY1\nrGLtjrVUDKhgcvlkJp4zkfPPPp9xZ49jXNk4zh18Ln169om6bEkxYCAkXFNzE8s3L+dXtb/i5fde\nZvnm5ZxXdh6Xj76cy0ZdxsUVF1M5rJIBvQdEXaqkmDMQEujwscMsXL+QZ9Y8w+LaxYwYMIKrL7ia\n6edP5/Ixl1PWtyzqEiUlkIGQEM0tzby48UWefPtJXlj3AlNHTmVm5Uyuu/A6Rg8aHXV5klLAQIi5\nXYd38cPqH/L4m4/Tv1d/7r7kbm6uvJmKARVRlyYpZfIVCF6HkGfbDm7jod89xNzquVw34TqevOlJ\npo2e5umMkmLPQMiTvR/u5RuvfYO5b83ls5M/y4q/XsGYwWOiLkuSOsxA6KKm5ia+X/197n3lXq6/\n8HpWz1rNyIEjoy5LkjrNQOiCNdvXcMfP7qBXj14svG0hl4y4JOqSJClnfttpDpqam/jW777FFT+8\ngtsvvp1X73zVMJCUeI4QOmnnoZ3c9uxtHDx2kGVfXMb5Z58fdUmSlBeOEDrhzc1vMvWJqUwePpmX\nP/+yYSApVRwhdNDP3/05X3j+Czx23WPcXHlz1OVIUt4ZCB3wvTe/x72Ze/nFZ3/BpaMujbocSSoI\nA+EMvv36t/nOsu/w6h2vMuGcCVGXI0kFYyCcxqPLHuU7y75D5o4M5w4+N+pyJKmgDIR2zH1rLt/6\n3bcMA0ndRty+YCcWX263aMMi7vzZnbx252uMHzI+6nIk6bT8crsCWdmwks//9PP89JafGgaSuhWv\nQ8iy49AOblhwA49c8wifOPcTUZcjSUXllFGouaWZ6+dfzx8P/2Me/IsHI6lBknKRrykjRwihB37z\nAPuO7OMbf/6NqEuRpEh0JRBmA+8Aq4D5QB9gCLAEWAcsBsra7L8eWAtc1YX3zbvXP3idh994mB/d\n/CN6lfaKuhxJikSugTAW+CJwCTAZKAU+A9xDEAgXAi+FjwEqgVvC+xnAd7vw3nl1+Nhh7vjZHfz7\ntf/u3ziW1K3l+qG8DzgG9CM4U6kfsBm4AZgX7jMPuCncvhFYED6nDtgAxOI7IP755X9mSsUUPl35\n6ahLkaRI5Xra6S7gIeAPwGHgVwQjg3KgIdynIXwMMBJYmvX8emBUju+dN69/8DrzV89n1ZdXRV2K\nJEUu10C4APg7gqmjvcCPgc+12aclvLXnlD+bM2fOie2qqiqqqqpyLPH0mpqbmLVwFg9d9RBD+w0t\nyHtIUiFkMhkymUzeXzfX05RuAf4CuDt8/FfANODPgU8CW4ERwCvARFrXEu4P7xcB9wFvtHndop12\n+vjyx1mwegGZ2zPHT9mSpESK+rTTtQQBcFZYxHSgBvg5cHu4z+3Ac+H28wSLzr2BccAEYFmO791l\nOw/t5N5X7uWRGY8YBpIUynXK6G3gSWA50Ay8BTwBDASeBu4iWDyeGe5fE7bXAI3ALE4/nVRQczJz\nmDlpJhdXXBxVCZIUO3H753HBp4zq9tQx9YmprPnKGob3H17Q95KkYoh6yiix/uXX/8KsP51lGEhS\nG93q207XbF/DC+teYP1X10ddiiTFTrcaIfzrq//KP1z+D5T1LTvzzpLUzXSbNYTaXbVcNvcyNn5t\nI4P6DCrIe0hSFFxD6KR/e/3f+NLULxkGktSObrGGsO3gNuavns+ar6yJuhRJiq1uMUJ4dNmjzKyc\nScWAiqhLkaTYSv0I4WjTUZ548wkyd2SiLkWSYi31I4Rn1zzLpOGTmDh0YtSlSFKspT4QHlv+GF/+\n0y9HXYYkxV6qA+Gdbe+wbuc6bvyjG6MuRZJiL9WB8B9v/gd3T7nbv5MsSR2Q2kXlo01Hmb9qPsu/\ntDzqUiQpEVI7Qli4fiGThk9ibNnYqEuRpERIbSA8tfIpPje57V/1lCS1J5WBsPvwbpZsXMJfTvrL\nqEuRpMRIZSD8uObHXHXBVX6rqSR1QioDYcHqBdw2+baoy5CkREldIGw/uJ3qLdXMGD8j6lIkKVFS\nFwjPv/s8V4+/mr49+0ZdiiQlSuoC4dm1z/KpiZ+KugxJSpxUBcLeD/fy2vuvce2Ea6MuRZISJ1WB\nsHD9Qq4ceyUD+wyMuhRJSpxUBcJz7z7HTX90U9RlSFIipSYQGpsbWVK7hGsmXBN1KZKUSKkJhN9v\n+j1jBo9h5MCRUZciSYmUmkD45YZfcs14RweSlKvUBMKiDYu8GE2SuiAVgbD94HbW7VzHx8d8POpS\nJCmxUhEISzYu4ZPjPknv0t5RlyJJiZWKQHhp40tMHzc96jIkKdFSEQiZ9zNUja2KugxJSrTEB8IH\nez9g/5H9VA6rjLoUSUq0xAfCr9//NVeOvZKSkpKoS5GkREt8IGTqMlx53pVRlyFJideVQCgDngHW\nADXAZcAQYAmwDlgc7nPcbGA9sBa4qgvv+xGZOtcPJCkfuhIIDwMLgYuAPyH4oL+HIBAuBF4KHwNU\nAreE9zOA73bxvQGo31fPng/3uH4gSXmQ64fyYOAK4Afh40ZgL3ADMC9smwcc/+rRG4EFwDGgDtgA\nXJrje5/wmz/8hivOu4IeJYmf+ZKkyOX6SToO2A78EHgL+B7QHygHGsJ9GsLHACOB+qzn1wOjcnzv\nE5bWL+Xy0Zd39WUkSeQeCD2BSwimfi4BDtI6PXRcS3hrz+l+1iFL65cybfS0rr6MJInggz0X9eHt\n9+HjZwgWjbcCFeH9CGBb+PNNwJis548O204yZ86cE9tVVVVUVVWdsoAjjUdYtW0VU0dMzbELkpRM\nmUyGTCaT99ftysn7rwJ3E5xRNAfoF7bvBB4gGDGUhfeVwHyCdYNRwIvAeE4eJbS0tHRs4PBG/Rv8\nzS/+huq/ru5CFyQp+cLrsLp8MVauIwSArwL/H+gN1AJ3AqXA08BdBIvHM8N9a8L2GoIF6Fl0ccpo\naf1Spo1yukiS8iVul/d2eIRw609u5eoLruaOj91R2IokKebyNUJI7PmaLihLUn4lMhC2H9zO7sO7\nufCcC6MuRZJSI5GBUL21mikjpnhBmiTlUSI/Uau3VPOx8o9FXYYkpUoiA2FFwwqmjJgSdRmSlCqJ\nDITqLdVMqTAQJCmfEhcIB44e4A97/8DEoROjLkWSUiVxgbCyYSWVwyrpVdor6lIkKVUSFwgrtq5w\nukiSCiBxgVC9pdoFZUkqgMQFwtsNb3Nx+cVRlyFJqZOoQGhuaWbNjjVMGj4p6lIkKXUSFQgf7P2A\nQX0GUda3LOpSJCl1EhUINdtrmDTM0YEkFUKiAuGd7e9QOawy6jIkKZUSFQg122sMBEkqkMQFglNG\nklQYiQmElpYWarbXcNGwi6IuRZJSKTGBUL+vnv69+zPkrCFRlyJJqZSYQHC6SJIKKzGBsGbHGr/h\nVJIKKDGBsH7nev+GsiQVUHICYdd6JgyZEHUZkpRaiQmEDbs2MH7I+KjLkKTUSkQgHG06yub9mxlb\nNjbqUiQptRIRCBt3b2TM4DH+lTRJKqBEBMKGXRtcP5CkAktEIKzfud71A0kqsGQEgmcYSVLBJScQ\nzjEQJKmQEhEIriFIUuHFPhCOn3J6Xtl5UZciSakW+0Co31fPyIEj6dmjZ9SlSFKqxT4Q6vbUcd5g\nRweSVGiJCASvUJakwot9ILy/531HCJJUBF0NhFKgGvh5+HgIsARYBywGyrL2nQ2sB9YCV3X0Der2\nOkKQpGLoaiB8DagBWsLH9xAEwoXAS+FjgErglvB+BvDdjr73+3veNxAkqQi6EgijgWuBuUBJ2HYD\nMC/cngfcFG7fCCwAjgF1wAbg0o68Sd2eOk85laQi6EogfBv4R6A5q60caAi3G8LHACOB+qz96oFR\nZ3qDxuZGthzYwuhBo7tQpiSpI3INhOuBbQTrByXt7NNC61RSez8/rU37NjG8/3B6l/bufIWSpE7J\n9WqvjxNMD10L9AUGAf9JMCqoALYCIwhCA2ATMCbr+aPDtpPMmTPnxHbZxDLXDySpjUwmQyaTyfvr\ntvev+864EvhfwP8AHgR2Ag8QLCiXhfeVwHyCdYNRwIvAeE4eJbS0tLQ2Pfn2kyyuXcxTn3oqD2VK\nUjqVlJRAHj7P8/V9EMc/xe8HngbuIlg8nhm214TtNUAjMIsOTBl5lbIkFU8+AuHX4Q1gFzC9nf2+\nGd46rH5fPZeMuKQLpUmSOirWVypv2r+JUQPPeDKSJCkP4h0I+zYxapCBIEnFEO9AcIQgSUUT20A4\n0niEfUf2Maz/sKhLkaRuIbaBsHn/ZioGVNCjJLYlSlKqxPbTdtP+TYwcODLqMiSp24hvIOxz/UCS\niim+geCCsiQVVXwDwVNOJamo4hsIjhAkqajiHQiOECSpaOIbCC4qS1JRxTIQWlpa2HJgi6edSlIR\nxTIQ9h3ZR88ePenfu3/UpUhStxHLQGg42EB5//Iz7yhJypt4BsKBBsoHGAiSVEzxDISDDQzvPzzq\nMiSpW4llIGw7uM0pI0kqslgGQsMB1xAkqdjiGQgHXUOQpGKLbyA4QpCkoopnIHiWkSQVXTwDwRGC\nJBVdPAPBEYIkFV3sAuHg0YM0tTQxsPfAqEuRpG4ldoFwfLqopKQk6lIkqVuJXSBsO7jN6SJJikAs\nA8GvrZCk4otdIOw4tIOh/YZGXYYkdTuxC4Sdh3ZyzlnnRF2GJHU78QuEwwaCJEUhdoGw49AOzuln\nIEhSscUuEBwhSFI04hcIh3a6qCxJEYhfIBze6ZSRJEUgfoHgWUaSFIlcA2EM8ArwDrAa+NuwfQiw\nBFgHLAbKsp4zG1gPrAWuau+Fdx3e5QhBkiKQayAcA/4emARMA74CXATcQxAIFwIvhY8BKoFbwvsZ\nwHfbe+8+PfvQu7R3jmVJknKVayBsBVaE2weANcAo4AZgXtg+D7gp3L4RWEAQJHXABuDSU72w00WS\nFI18rCGMBaYAbwDlQEPY3hA+BhgJ1Gc9p54gQE7iGUaSFI2eXXz+AOAnwNeA/W1+1hLe2nPKn+1e\ntJs5m+cAUFVVRVVVVRdLlKR0yWQyZDKZvL9uV/7oQC/gBeCXwP8N29YCVQRTSiMIFp4n0rqWcH94\nvwi4j2BUka3l1mduZf6n53ehLEnqXsK/H9PlPyKT65RRCfB9oIbWMAB4Hrg93L4deC6r/TNAb2Ac\nMAFYdqoXdg1BkqKR65TRJ4DPASuB6rBtNsEI4GngLoLF45nhz2rC9hqgEZhFO1NGnnIqSdHINRB+\nQ/uji+nttH8zvJ3W2X3PzrEkSVJXxO5K5bPPMhAkKQqxC4SyvmVn3kmSlHexC4TBfQZHXYIkdUux\nCwRHCJIUDQNBkgQYCJKkUJevbMuzlqbmJnqUxC6nJCm2or5SuWAMA0mKhp++kiTAQJAkhQwESRJg\nIEiSQgaCJAkwECRJIQNBkgQYCJKkkIEgSQIMBElSyECQJAEGgiQpZCBIkgADQZIUMhAkSYCBIEkK\nGQiSJMBAkCSFDARJEmAgSJJCBoIkCTAQJEkhA0GSBBgIkqSQgSBJAgwESVLIQJAkAcUPhBnAWmA9\n8PUiv7ck6TSKGQilwKMEoVAJ3ApcVMT3j1wmk4m6hIJJc9/A/iVd2vuXL8UMhEuBDUAdcAz4EXBj\nEd8/cmn+pUxz38D+JV3a+5cvxQyEUcAHWY/rwzZJUgz0LOJ7tXRkp+nTobERmpry98YlJfF4bl0d\ndOYfKnGpuyM2boTXXovmvYuhthZ++9uoqygc+yeAYv7vNw2YQ7CGADAbaAYeyNpnA3BBEWuSpDSo\nBcZHXURn9CQoeizQG1hBN1tUliS1ugZ4l2AkMDviWiRJkiTFVVouWKsDVgLVwLKwbQiwBFgHLAbK\nsvafTdDntcBVRauy434ANACrstpy6c/U8DXWAw8XsN7OOFXf5hCc/VYd3q7J+lmS+gYwBngFeAdY\nDfxt2J6W49de/+aQjmPYF3iDYGq9Bvg/YXtajl+7SgmmkMYCvUj22sJ7BAcs24PAP4XbXwfuD7cr\nCfrai6DvG4jfV4lcAUzhox+anenP8ZMWlhFchwKwkNYTC6J0qr7dB/zPU+ybtL4BVAAfC7cHEEzV\nXkR6jl97/UvTMewX3vcElgL/jQIfvzh8AKXtgrW2Z27dAMwLt+cBN4XbNwILCPpcR/Df4FLi5TVg\nd5u2zvTnMmAEMJDWEdOTWc+J0qn6Bqc+8y5pfQPYSvABAXAAWENw3U9ajl97/YP0HMND4X1vgn84\n76bAxy8OgZCmC9ZagBeB5cAXw7ZygqkJwvvycHskQV+PS0q/O9uftu2biHc/vwq8DXyf1uF40vs2\nlmA09AbpPH5jCfq3NHyclmPYgyD0GmidHivo8YtDIHTogrWE+ATBL+Y1wFcIpiWytXD6/ibtv8WZ\n+pM0jwHjCKYitgAPRVtOXgwAfgJ8Ddjf5mdpOH4DgGcI+neAdB3DZoJ+jAb+DPhkm5/n/fjFIRA2\nESwQHTeGjyZakmwJ77cDPyWYAmogmO+EYPi2Ldxu2+/RYVvcdaY/9WH76Dbtce3nNlr/J5tL6xRe\nUvvWiyAM/hN4LmxL0/E73r+naO1f2o4hwF7gFwSLw2k6fqeUlgvW+hHM1QH0B35LsNL/IK1nTt3D\nyYtAvQn+RVNLca8c76ixnLyo3Nn+vEEwn1lCvBbtxvLRvo3I2v57YH64ncS+lRDMF3+7TXtajl97\n/UvLMRxK63TXWcCrwH8nPcfvtNJwwdo4ggOyguA0uOP9GEKwrnCq08T+N0Gf1wJXF63SjlsAbAaO\nEqzz3Elu/Tl+2tsG4JGCV90xbfv2BYIPmJUE88/P0To/C8nqGwRnpDQT/D4ePwVzBuk5fqfq3zWk\n5xhOBt4i6N9K4B/D9rQcP0mSJEmSJEmSJEmSJEmSJEmSJElSGvwXHD3kMovaOC0AAAAASUVORK5C\nYII=\n",
"text": [
"<matplotlib.figure.Figure at 0x7efedbc05590>"
]
}
],
"prompt_number": 11
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"dsolve([0,100], 300)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAEACAYAAABS29YJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4FFW6x/EvBJB1DGsgLCasISAg+yLYKiKgw6KIwQ0E\nXIBB1GFU9Dpk5rkyigrecRkHGQUZARFkGySy2QxrWMOSEJKwCAkQVlmEEJL0/aOKJERCOkl3qpff\n53n6qerq6ur3WPjm9Huq64CIiIiIiIiIiIiIiIiIiIiIiIiIiLjQOGAPsNdcB6gGrAQSgBVAYK79\nJwCJQDzQq+TCFBGRomqJkejLAwEYCb4RMBl4zdzndeBdcz0ciAHKAiFAElC65MIVEZGbKSgRhwHR\nQBqQCawFHgX6ATPNfWYCA8z1/sAc4BpwGCPZd3RpxCIiUmgFJfu9QHeMsk1FoC9QDwgCUs19Us3n\nAMFAcq73JwN1XRWsiIgUTZkCXo8H3sOoy/+KUaLJzLOPw3zk51aviYhICSgo2QN8aT4A3sHoracC\ntYETQB3gpPl6ClA/13vrmdtu0KhRI8eBAweKGLKIiN86ADQuyhudGTytZS4bAI8As4ElwFBz+1Bg\nkbm+BIgAygGhQBNgy2+iPXAAh8Phs4+JEydaHoPap/b5W9v8oX0YF8gUiTM9+/lAdYxB19HAeYyr\nb+YBIzAGYgeb+8aZ2+OADHN/lXFERCzmTLLvcZNtZ4Ge+ew/yXyIiIiH0DXwbmCz2awOwa3UPu/l\ny20D329fcZSy6HMdZv1JREScVKpUKShi3lbPXkTEDyjZi4j4ASV7ERE/oGQvIuIHlOxFRPyAkr2I\niB9QshcR8QNK9iIifkDJXkTEDyjZi4j4ASV7EREvkJaRVqz3O3PXSxERKQFpGWkcPHeQxDOJJJ1N\nIvFsIolnjfXUS6kFH+AWdCM0EZES5HA4SLmYwv7T+9l/Zn/2Mv50PCcuneCOwDtoUq0Jjas1pkm1\nJjSpbqw3uL0BZQPKQhHztpK9iIgb/Jr+KwlnErITevyZePaf3k/CmQQql6tMsxrNaFbdfJjroVVD\nKVM6/4JLce566cybJgBPAVnAHuBZoBLwLXAHOTNV/ZJr/+EYE5O/hDFZeV5K9iLiE079eorYU7HE\nnowl7lSckdzP7Of05dM0rtb4hoQeViOMptWbElg+sEif5c5kHwKsAZoDVzES/A9AC+A0MBl4HagK\nvAGEY8xR2wGoC6wCmmL8ochNyV5EvMr1pB53Ko7Yk7FGgj8VS0ZWBi1qtqBFzRY0r9mcsBphNKve\njAa3NyCgdIBLYyhOsi9ogPYCxtyzFTF66hWBYxi993vMfWYCdoxk3x+YY77nMJAEdAQ2FyU4EZGS\ndvry6ZxkfjKWuNNGck/PTKdFrRbZiX1g84GE1wynTuU615OwRyso2Z8FPgSOAFeAH4GVQBBwfWg4\n1XwOEMyNiT0Zo4cvIuJR0jPTiT8dz64Tu9iduptdqcYyLSMtO6mH1wxnQNgAWtRq4TVJPT8FJftG\nwMsY5ZzzwHcY9fvcHOYjPzd9LTIyMnvdZrNp7kgRcZvUS6nZyfz6MuFMAqGBobQKakXroNaM6zSO\nVkGtqPe7eh6T1O12O3a73SXHKqhFjwMPACPN508DnYH7gHuBE0Ad4CcgDKOUA/CuuYwCJgLReY6r\nmr2IuFxGVgb7Tu0j5kTMDcn9WuY1WtduTeug1tnJPbxmOBXKVrA65EJx5wBta+AbjAHXNGAGsAXj\nKpwzwHsYCT6QGwdoO5IzQNuY3/bulexFpFjSM9OJPRnLjuM72H58OzuO72DPyT3U+1097qp9V05i\nr92aulXqekxvvTjcOUC7C/ga2IZxRc0OYBpQBZgHjCDn0kuAOHN7HJABjObWJR4RkQKlZaSx9+Re\nth/bnp3c407FEVo1lLZ12tKuTjsiWkbQpnYbfnfb76wO1yPpR1Ui4lGuZV5jz8k9bEnZwtaUrew4\nsYP9p/fTpHoT2tZpS9vabWkX3I7WQa2pVK6S1eGWKHf/qModlOxFBIfDweFfDrMlZQvRKdFEp0QT\ncyKGkMAQOtXtRIfgDrQLbsedte70uvq6OyjZi4hXOHflHFuPbSU62UjsW1K2UKZ0GTrV60TH4I50\nqteJ9sHtVYrJh5K9iHicLEcWsSdj2XB0AxuPbiQ6JZpjF4/Rrk47OtXtRMe6RnL3lcHTkqBkLyKW\nu3ztMltStrDhyAY2HN3ApuRN1KhYg271u9Gtfjc61+tMeM1wl99CwJ8o2YtIiTt+8Tgbjm7ITu6x\np2JpFdQqO7l3rd+VoMpBBR9InKZkLyJu5XA4OHjuIPbDdtb+vJb1R9Zz/up5utbvmp3c2we31yCq\nmynZi4hLORwOks4mZSd3+2E7DhzYQmzcc8c93N3gbsJqhFG6lGY2LUlK9iJSLA6Hg4QzCdmJ3X7Y\nTkDpgOzkbgux0ahqIw2kWkzJXkQK7dC5Q6w6uIo1h9dgP2ynbOmy2EJs2Y/QwFAldw+jZC8iBTp7\n5SxrDq1h1cFVrDq4ikvpl+jZsCf3h96PLcRGSGCIkruHU7IXkd9Iy0hj49GNrDywklWHVrH/9H66\n39GdnqE96dmwJy1rtVRy9zJK9iKCw+Fg78m9RCVFsfLgSjYlb6JlrZb0DO3JA40eoHO9zpQLKGd1\nmFIMSvYifurC1QusOriK5YnLWZ60nHIB5ejduDe9GvXCFmIr8sTW4pmU7EX8hMPhIPZULD8k/sDy\npOVsO7aNLvW60LdJX/o07kPT6k1VmvFhSvYiPuxS+iVWHVyVneDLlC5Dn8Z96NukL/eG3Ot3t/n1\nZ0r2Ij7m2MVjLNm/hMX7F7P+yHo61e1E3yZ96dukL82qN1Pv3U+5O9k3A+bmet4QeBv4N/AtxhSF\nhzFmq/rF3GcCMBzIBF4CVuQ5ppK9SC4Oh4M9J/dkJ/gDZw/Qp0kf+jfrz4ONHuT28rdbHaJ4gJLs\n2ZcGUjDmmB0LnAYmA68DVblxHtoO5MxD2xRjWsPrlOzF713LvMa6I+tYHL+YJQlLAOjfrD/9mvWj\ne4PulA0oa3GE4mncOQdtXj2BJOAo0A+4x9w+E7BjJPv+wBzgGkaPPwnjj8PmogQo4kvSMtJYcWAF\n8+Pm85+E/9C4WmP6NevHkogluu5d3KqwyT4CI5EDBAGp5nqq+RwgmBsTezJGD1/EL12+dpmopCjm\nx83nh8QfaFO7DYPCB/G3+/9G3d/pfw0pGYVJ9uWA32OUbPJymI/8/Oa1yMjI7HWbzYbNZitEKCKe\n7VL6JX5I/IH5cfP58cCPdAjuwKDwQUx5cAq1K9e2OjzxEna7Hbvd7pJjFeY7Y39gFNDbfB4P2IAT\nQB3gJyAMo5QD8K65jAImAtG5jqWavficS+mXWLp/KfP3zWflgZV0qd+Fx8Ifo3+z/tSsVNPq8MQH\nlNQA7VxgOUZ9HoyB2TPAexgJPpAbB2g7kjNA25gbe/dK9uIT0jPTWXFgBbP3zGZZ4jK61u/K4PDB\n9GvWj+oVq1sdnviYkkj2lYCfgVDgormtGjAPaMBvL718E+PSywxgHPBjnuMp2YvXynJksf7Iembv\nmc38uPmE1QjjiTuf4LHwx9SDF7fSj6pE3MzhcLArdRez98xmzt45VKtQjSEthxDRMoKQwBCrwxM/\noWQv4iYpF1L49+5/M3PXTK5kXOGJlk8w5M4htKzV0urQxA8p2Yu4UFpGGovjFzNj1wyik6MZFD6I\noa2H0rV+V10HL5ZSshcpJofDQXRKNDNjZjIvbh7tg9szrPUwBoQNoELZClaHJwKU7C9oRXxKyoUU\nZu2exYyYGThwMKz1MGJeiKH+7fWtDk3EpdSzF7+TmZVJVFIU03ZMY93P6xgUPohn2zxL53qdVaYR\nj6YyjogTki8k8+XOL5m+Yzp1qtTh+bbPE9EyQveDF6+hMo5IPjKzMlmetJxp26ex/sh6hrQcwpIh\nS2hTu43VoYmUKPXsxScdu3iML7Z/wfSd0wmuEswL7V7g8RaPqxcvXk09exGMK2o2JW/i4y0fE5UU\nRUSLCJYOWapevAjq2YsPSMtIY+7euXy85WPOp53nDx3/wLA2wwgsH2h1aCIupQFa8UtHzx/lH9v+\nwfQd02kX3I6xHcfSu3FvSpcqbXVoIm6hMo74lU1HNzFl8xTWHFrDU3c+xfrh62lavanVYYl4NPXs\nxStkZmWyeP9iPtz0IccvHueVzq8wrM0wqtxWxerQREqMevbisy5fu8yMmBlM2TSFGhVrML7reAaG\nDSSgdIDVoYl4FSV78Uipl1L5dOunfL7tc7o16MbMATN1IzKRYlCyF49y8NxBJm+YzLex3xLRIkL1\neBEXcfayhUBgPrAPiAM6YcxUtRJIAFaY+1w3AUjEmKe2l6uCFd+179Q+nln4DB2+6ED1CtXZ/4f9\n/OPhfyjRi7iIs9+JZwJrgS8xvg1UAt4CTmPMRfs6UJUb56DtQM4ctE2BrFzH0wCtALDz+E7eWfcO\n646s46WOLzGm4xhdHy+SD3dfZ387sBNomGd7PHAPkArUBuxAGEavPgtjInKAKCAS2JzrvUr2fm7j\n0Y28s+4dYk7EML7LeJ5v97xuZSBSAHdfjRMKnAK+AloD24GXgSCMRI+5DDLXg7kxsSdj9PBFWH9k\nPX/+6c8c+uUQr3d7nQWDF1C+THmrwxLxec4k+zJAW+APwFbgI4xyTW4O85Gf37wWGRmZvW6z2bDZ\nbE6EIt5qS8oW3v7pbRLOJPB2j7d5utXTlA0oa3VYIh7Nbrdjt9tdcixnvg7UBjZh9PAB7sYo1TQE\n7gVOAHWAnzDKONf/ELxrLqOAiUB0rmOqjOMndh7fyZ/tfybmRAxvdX+L4XcNp1xAOavDEvFKxSnj\nOHM1zgngKMYgK0BPIBZYCgw1tw0FFpnrS4AIoBzGH4gmwJaiBCfea+/JvTw671Eemv0QDzR8gMSx\nibzY/kUlehGLOHud/VjgG4wEfgB4FggA5gEjgMPAYHPfOHN7HJABjObWJR7xIQfPHeR/1vwPqw+t\n5k9d/8SsgbOoWLai1WGJ+D3dG0dc4tSvp/jf//4v3+z5hnGdxvFy55d13xoRF9O9ccQyl69dZuqm\nqUzdPJUhLYcQNyaOWpVqWR2WiOShZC9FkpGVwYyYGUy0T6Rb/W5sHrmZxtUaWx2WiORDyV4KxeFw\nsCxxGa+tfI1alWqx8PGFdKzb0eqwRKQASvbitH2n9vHKj6/w8/mf+eCBD+jbpK/uQiniJTR/mxTo\n3JVzjFs+jh4zetCncR92v7ibh5o+pEQv4kWU7CVfGVkZfLb1M8I+DeNa1jX2jdnHuM7j9MtXES+k\nMo7c1JpDaxgXNY6aFWuy8umVtApqZXVIIlIMSvZyg2MXj/Hqj68SnRLNh70+ZGDYQJVrRHyAyjgC\nGCWb/9v8f7T+vDWNqzUmdnQsjzR/RIlexEeoZy9sTt7MqGWjqFq+KuueXUdYjTCrQxIRF1Oy92Nn\nr5xlwqoJLE1Yyge9PmBIyyHqyYv4KJVx/JDD4eCb3d8Q/mk45QLKETcmjifufEKJXsSHqWfvZ46e\nP8qoZaM4cv4IS4Ys0a9fRfyEevZ+IsuRxefbPqfttLZ0qtuJbc9vU6IX8SPq2fuBxDOJjFw6kqsZ\nV7EPtdOiVgurQxKREuZsz/4wsBvYSc6sU9WAlUACsAIIzLX/BCARiAd6uSJQKbyMrAze3/A+Xf7V\nhYFhA9kwfIMSvYifcnZE7hDQDjiba9tk4LS5fB2oijH/bDgwG+gA1AVWYUxpmJXrvZq8xM2Szibx\nzMJnKF+mPNP7Tadh1YZWhyQixeTuOWizPyfP837ATHN9JjDAXO8PzAGuYXwjSAJUHC4hDoeDz7d9\nTpd/dSGiZQSrnlmlRC8iTtfsHRg99Ezgn8AXQBCQar6eaj4HCAY253pvMkYPX9ws5UIKI5aM4MyV\nM/pxlIjcwNmefTfgLqAPMAbonud1B7eeVFw1Gzebu3cubae1pWv9rmwcvlGJXkRu4GzP/ri5PAUs\nxCjLpAK1gRNAHeCkuU8KUD/Xe+uZ224QGRmZvW6z2bDZbM5HLdl+SfuFUctGEXMihmVPLKN9cHur\nQxIRF7Hb7djtdpccy5lCf0UgALgIVMK48uYvQE/gDPAexsBsIDcO0HYkZ4C2MTf27jVA6wKbjm7i\nie+foG/jvnzQ6wMqlK1gdUgi4kbFGaB1pmcfhNGbv77/NxgJfxswDxiBMRA72NwnztweB2QAo1EZ\nx6WyHFm8t/49Por+iH8+/E8GhA0o+E0i4tesuhmKevZFdPzicZ5e+DTpmel888g31L+9fsFvEhGf\nUFKXXorFlicup+20tnRv0J01Q9co0YuI03S7BC+QkZXBW6vfYs7eOXw76Ft63NHD6pBExMso2Xu4\nE5dOMGTBEMoFlGPHCzuoUbGG1SGJiBdSGceDrT+ynvbT2tOjQQ9+eOIHJXoRKTL17D2Qw+Fg6uap\nvLfhPb7q/xV9m/S1OiQR8XJK9h7mwtULDF88nMO/HCZ6ZDQhgSFWhyQiPkBlHA+SeCaRztM7U61C\nNdYPX69ELyIuo2TvIVYeWMndX93N2I5jmfb7aZQvU97qkETEh6iMYzGHw8HHWz5m0rpJzBs0j3tC\n7rE6JBHxQUr2FrqacZXRy0az9dhWNo3YRGjVUKtDEhEfpWRvkdRLqTwy7xFqV67NxhEbqVyustUh\niYgPU83eArEnY+k0vRM9Q3vy3WPfKdGLiNupZ1/C1hxaQ8T8CKY8OIWnWj1ldTgi4ieU7EvQzJiZ\nvLbqNeY9Ng9biM3qcETEjyjZlwCHw8Ff1v6Fr3d9jX2oneY1m1sdkoj4GSV7N0vPTGfkkpHEn45n\n04hNBFUOKvhNIiIu5uwAbQCwE1hqPq8GrAQSMGatCsy17wQgEYgHerkmTO90Kf0SD89+mAtXL2Af\nZleiFxHLOJvsx2FMM3h9eqk3MJJ9U2C1+RyM+WcfN5e9gc8K8Rk+5fTl09z/9f00uL0B8wfPp2LZ\nilaHJCJ+zJlEXA/oC0wnZzqsfsBMc30mcH0S1P7AHOAaxry0SRgTj/uVo+eP0v2r7twbci9f/P4L\nypRWtUxErOVMsp8K/AnIyrUtCEg111PN5wDBQHKu/ZKBusWM0avEn47n7q/uZuRdI3m357vX54wU\nEbFUQcn+YeAkRr0+v6zlIKe8k9/rfmFrylbunXkvf7H9hT92/aPV4YiIZCuovtAVo2TTFygP/A6Y\nhdGbrw2cAOpg/EEASAFyz4Jdz9z2G5GRkdnrNpsNm81W2Ng9yn9//i+D5g1ier/p9GvWz+pwRMQH\n2O127Ha7S45VmBrDPcB44PfAZOAM8B7G4GyguQwHZmPU6esCq4DG/LZ373A4fKfDv/rgaiIWRDD3\n0bnc3/B+q8MRER9lloWLVBsu7Mjh9Qz9LjAPGIExEDvY3B5nbo8DMoDR+HgZJyopimcWPsOCwQvo\ncUcPq8MREbkpq0YPfaJnv3T/UkYuHcmixxfRpX4Xq8MRER9Xkj17MS2IW8CYH8aw7IlltA9ub3U4\nIiK3pGRfBN/u/ZaXf3yZqKeiaFO7jdXhiIgUSGWcQlq4byGjlo1i5dMruTPoTqvDERE/ojJOCVmW\nsIwXl73I8ieXK9GLiFdRsnfSygMreXbxsywdspS2ddpaHY6ISKH45U3KCmvt4bU8+f2TfP/493Sq\n18nqcERECk09+wJsPLqRQd8N4ttB33J3g7utDkdEpEjUs7+FXSd2MfDbgXw94GvuC73P6nBERIpM\nyT4fB88dpO/svnzc52P6NOljdTgiIsWiZH8TqZdS6TWrF291f4vBLQYX/AYREQ+nZJ/H+bTz9P6m\nN0+1eorRHUZbHY6IiEvoR1W5pGWk0eebPjSv0ZxP+36qiUdExKMU50dVSvamzKxMBs8fTECpAOY8\nOoeA0gFWhyQicgP9gtYFxq8Yz9krZ4l6MkqJXkR8jpI98MmWT4g6EMXG4Ru5rcxtVocjIuJyfp/s\nlyUsY9K6Sawfvp6qFapaHY6IiFsUdDVOeSAaiMGYfepv5vZqwEogAViBMS3hdROARCAe6OXKYF1t\n5/GdDFs8jO8f/56GVRtaHY6IiNs4U+ivCFzG+BawHmMe2n7AaYy5aF8HqnLjHLQdyJmDtimQleeY\nlg/QJl9Ipsu/ujD1wakMCh9kaSwiIs4ozgCtM9fZXzaX5YAA4BxGsp9pbp8JDDDX+wNzgGsYc9Mm\nYUw+7lF+Tf+Vh2c/zEsdX1KiFxG/4EyyL41RxkkFfgJigSDzOeYyyFwPBpJzvTcZo4fvMRwOB8MW\nD+OuOncxvut4q8MRESkRzgzQZgFtgNuBH4F787zuMB/5uelrkZGR2es2mw2bzeZEKMU3ad0kjp4/\nytpha/WjKRHxaHa7Hbvd7pJjFTbbvQ1cAUYCNuAEUAejxx+GUbcHeNdcRgETMQZ5c7OkZr9k/xJG\nLxvNlue2EFwluMQ/X0SkONxZs69BzpU2FYAHgJ3AEmCouX0osMhcXwJEYNT3Q4EmwJaiBOZqcafi\nGLlkJAsGL1CiFxG/U1AZpw7GAGxp8zELWI2R8OcBIzAGYq/fGjLO3B4HZACjuXWJp0Scu3KO/nP7\n8/4D72umKRHxSz5/b5wsRxYPz36YptWb8lHvj0rkM0VE3MHdl156tUnrJnEx/SLvP/C+1aGIiFjG\np2+XsPrgaj7b+hnbnt9G2YCyVocjImIZn+3Zp1xI4emFTzNr4CwNyIqI3/PJZH8t8xoRCyIY3WE0\n9ze83+pwREQs55PJ/s3Vb1KlXBXe7P6m1aGIiHgEn6vZL92/lHlx89jx/A5Kl/LJv2UiIoXmU8n+\n2MVjPLf0ORYMXkD1itWtDkdExGP4TNc3y5HF0EVDGdV+FN0adLM6HBERj+IzyX7KpilcuXaFt3q8\nZXUoIiIexyfKONuPbWfyhslsfW4rZUr7RJNERFzK63v2l9IvMWTBEP7e5+/cEXiH1eGIiHgkr783\nzsglI8l0ZPJV/69ccjwREU9VnHvjeHXNY1nCMlYfWs3uF3dbHYqIiEfz2mR/9spZXvjPC8waOIsq\nt1WxOhwREY/mtWWcp75/imoVqvH3Pn93UUgiIp7N78o4C/ctJDolmpgXYqwORUTEKzhzNU59jDlm\nY4G9wEvm9mrASiABWEHO9IUAE4BEIB7o5apgAU5fPs2YH8Ywo/8MKpWr5MpDi4j4LGe+DtQ2HzFA\nZWA7MAB4FjgNTAZeB6piTDgeDswGOgB1gVVAUyAr1zGLXMYZ/N1g7rj9Dt7vpclIRMS/uHumqhMY\niR7gErAPI4n3w5ifFnM5wFzvD8wBrmHMT5sEdCxKcHktil/ErtRd/PXev7ricCIifqOwP6oKAe4C\nooEgINXcnmo+BwgGknO9Jxnjj0OxXLh6gbHLxzLt4WlUKFuhuIcTEfErhRmgrQwsAMYBF/O85jAf\n+fnNa5GRkdnrNpsNm812yw9/c/Wb9G7Um3tC7nEuWhERL2e327Hb7S45lrO1n7LAf4DlwEfmtnjA\nhlHmqYMxiBuGUbcHeNdcRgETMb4NXFeomv2mo5t4dN6jxI6OpWqFqk6/T0TEl7i7Zl8K+BcQR06i\nB1gCDDXXhwKLcm2PAMoBoUATYEtRggNIz0znuaXPMfXBqUr0IiJF5EwZpxvwFLAb2Glum4DRc58H\njMAYiB1svhZnbo8DMoDR3LrEc0vvb3ifkMAQBrcYXPDOIiJyUx79C9oDZw/QaXontj+/XXe0FBG/\n5+4yjmVe+fEVxncdr0QvIlJMHnu7hGUJy4g/Hc93j31ndSgiIl7PI5N9WkYa46LG8UnfT7itzG1W\nhyMi4vU8sowzZdMUWtZqSe/Gva0ORUTEJ3hcz/7I+SN8uOlDtj23zepQRER8hsf17MevGM/YjmMJ\nrRpqdSgiIj7Do3r2G45sYHPyZmYMmGF1KCIiPsVjevYOh4M/rvgj79z3DhXLVrQ6HBERn+Ixyf67\nuO9Iz0znyVZPWh2KiIjP8YgyztWMq7yx6g2m95tO6VIe8/dHRMRneERm/WzrZ7So1YL7Qu+zOhQR\nEZ9k+b1xzl45S9gnYawdtpbmNZtbFI6IiOfz6nvjTFo3iUeaP6JELyLiRpbW7FMupPBVzFfsHbXX\nyjBERHyepT37SesmMeKuEdSpUsfKMEREfJ5lPfvDvxxmbuxc9v9hv1UhiIj4DWd69l8CqcCeXNuq\nASuBBGAFEJjrtQlAIsYctb3yO+hf1/6VMR3GUKNijcLGLCIiheRMsv8KyHv7yTcwkn1TYDU5k4yH\nA4+by97AZ/l9xtKEpbza5dUihCwiIoXlTLJfB5zLs60fMNNcnwkMMNf7A3OAaxjz0iYBHW920Fc7\nv0pg+cCbvSQiIi5W1AHaIIzSDuYyyFwPBpJz7ZcM1L3ZAcZ2GlvEjxYRkcJyxQCtw3zc6vXf+GDS\nB9nrNpsNm83mglBERHyH3W7Hbre75FjO/hIrBFgK3Gk+jwdswAmgDvATEEZO7f5dcxkFTASi8xwv\n+xe0IiLiHCt+QbsEGGquDwUW5doeAZQDQoEmwJYifoaIiLiIM2WcOcA9QA3gKPBnjJ77PGAExkDs\nYHPfOHN7HJABjObWJR4RESkBlt8ITUREnOPVN0ITERH3U7IXEfEDSvYiIn5AyV5ExA8o2YuI+AEl\nexERP6BkLyLiB5TsRUT8gJK9iIgfULIXEfEDSvYiIn5AyV5ExA8o2YuI+AElexERP6BkLyLiB9yV\n7HtjTF2YCLzups8QEREnuSPZBwCfYCT8cGAI0NwNn+OxXDVBsKdS+7yXL7cNfL99xeGOZN8RSMKY\nrvAaMBfo74bP8Vi+/g9O7fNevtw28P32FYc7kn1djLlqr0s2t4mIiEWcmXC8sJyaXPahh8ydHUVf\nFvW97nbkCKxaVTKfZQW1zzuVKgU//wyrV1sdyY1KuXAm7J9/hjVrXHMsV8blCdzRnM5AJEbNHmAC\nkAW8l2ufJKCRGz5bRMSXHQAaWx3EdWUwAgoBygEx+NkArYiIv+gD7MfowU+wOBYREREREXEHX/vB\n1WFgN7Ddf8RfAAAC3ElEQVQT2GJuqwasBBKAFUCgJZEVzZdAKrAn17ZbtWcCxrmMB3qVUIzFcbP2\nRWJcNbbTfPTJ9Zq3ta8+8BMQC+wFXjK3+8I5zK9tkfjG+SsPRGOUvuOAv5nbvfLcBWCUdkKAsvhG\nPf8QxsnIbTLwmrn+OvBuiUZUPN2Bu7gxGebXnnCMc1gW45wm4fm34LhZ+yYCr95kX29sX22gjble\nGaOc2hzfOIf5tc2Xzl9Fc1kG2AzcjYvOXUk33Fd/cJX3qqZ+wExzfSYwoGTDKZZ1wLk82/JrT39g\nDsa5PIxxbju6P8RiuVn74OZXpnlj+05gJACAS8A+jN+5+MI5zK9t4Dvn77K5LIfROT6Hi85dSSd7\nX/zBlQNYBWwDnjO3BWGUCjCXQRbE5Ur5tScY4xxe583ncyywC/gXOV+Tvb19IRjfYqLxvXMYgtG2\nzeZzXzl/pTH+oKWSU7Jyybkr6WRfQj9rKlHdMP7R9QHGYJQJcnPgW+0uqD3e2NZ/AKEYJYLjwIe3\n2Ndb2lcZWACMAy7mec3bz2FlYD5G2y7hW+cvC6Md9YAewL15Xi/yuSvpZJ+CMchyXX1u/MvkjY6b\ny1PAQoyvUakY9UWAOsBJC+Jypfzak/d81jO3eZuT5PxPNJ2cr8Le2r6yGIl+FrDI3OYr5/B62/5N\nTtt87fwBnAeWAe3w0nPnaz+4qghUMdcrARswRsQnk3Ol0Rt41wAtGOcn7wDtzdpzfYCoHEbP6gDu\n+VW2q4VwY/vq5Fp/BZhtrntj+0oBXwNT82z3hXOYX9t85fzVIKcEVQH4L3A/XnzufOkHV6EY/7Fj\nMC4Fu96eahh1fG+89HIOcAxIxxhfeZZbt+dNjHMZDzxYopEWTd72DcdIILsxar6LuHGMxdvadzdG\nKSCGnEsRe+Mb5/BmbeuD75y/O4EdGO3bDfzJ3O4L505ERERERERERERERERERERERERERERERERE\npGT9P5XjBxTL4KCCAAAAAElFTkSuQmCC\n",
"text": [
"<matplotlib.figure.Figure at 0x7efedbc97b90>"
]
}
],
"prompt_number": 12
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"sim.reset()"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 13
},
{
"cell_type": "heading",
"level": 1,
"metadata": {},
"source": [
"initial condition: A == B"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"sim.increment(B, 100)\n",
"sim.increment(A, 100)"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 14
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"sim.run(10000)"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 15
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"sim.plot()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEPCAYAAAC+35gCAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8E3X6wPFPekKhHAVaWoqW+0bkEnDFiAqKyCGgrLqL\nwHrxQ1lPRF0ouCji4rUr64VaVBDUFdQVBNTgAnKJnAW5j5a2FFpoCz1oM78/vpM0PWjTNskk6fN+\nveaVyWQy87SEPvneIIQQQgghhBBCCCGEEEIIIYQQQgghhBBCCCGEMNBUYDewR98HiADWAAeA1UAj\nh/OnAweB/cBgz4UphBDCaF1RCaMOEIhKFG2AecDT+jnTgLn6fmdgBxAMxAGHgADPhSuEEKIy7vyj\n3BHYDOQBRcA6YDQwHEjQz0kARur7I4AlwCXgGCpp9HVjfEIIIarInUljD3AdqjoqDBgKxAJRQJp+\nTpr+HCAGSHJ4fxLQwo3xCSGEqKIgN157P/Ayqt3iAqrqqajUOZq+XU5FrwkhhPAwdyYNgA/0DWAO\nqvSQBjQHUoFo4LT+ejLQ0uG9sfqxEtq0aaMdPnzYXfEKIYS/Ogy0relF3N3QHKk/XgHcASwGvgbG\n68fHA8v1/a+BcUAI0ApoB2wpfcHDhw+jaZpXbTNnzjQ8Bl+JS2KSmGpDXN4YE6ojUo25u6TxBdAE\n1bg9GTiP6i21DJiEavC+Uz83UT+eCBTq50v1lBBCeBF3J42B5RzLAG66zPkv6psQQggvJOMgXMBs\nNhsdQrm8MS6JyTkSk/O8MS5vjMlVTEYHUA2aXj8nhBDCSSaTCVzwN9/d1VNCCOH1IiIiyMzMNDoM\nl2jcuDEZGRluu76UNIQQ1VZQVEDo30MpmlFEgMl3a7tNJhP+8nflcj+Lq0oavvuvLIQw3NI9SwFY\ne2Qtp7JP2Y9rmsZTq59i+X7Vo37+xvmczztvSIzCtSRpCCGqLTQoFIAhnwyhxastyL2UC8DniZ/z\nj1/+wailo7BqVp5c8yS3L7ndyFCFi0ibhhCi2n4/83uJ52EvhpU55/HvHwdgSJshHolJuJeUNIQQ\n1TbDMqPSc97Y/AYAdYLquDsc4QGSNIQQ1XLi/IkKX/9i7Bclni//ffllzhTOMJvNREREUFBQYGgc\nkjSEENWyZPcSACb0mADAnV3uLPH66M6jSzxff2K9ZwLzQ8eOHWPLli1ERkby9ddfGxqLdLkVQlSL\naZb686HN1EjOSqZFg7LL3+ReymXJniX0jO7J1e9cTea0TBrVaVTmPKN5e5fb2bNns23bNq655ho2\nbdrEN998c9lz3d3lVpKGEKLKNE0jYLaqqNBmVv7/MfdSrr2R3JnzPc3bk0bbtm2ZNWsWffv2pUuX\nLiQlJREZGVnuuTJOQwjhVfIL89mcvBlwPgHUDa7rzpDczmRyzVYd69evJzk5meHDh9OuXTs6d+7M\n4sWLXfsDVoEkDSFEldSZU4f+C/tX+X0//PkHAHal7XJ1SG6naa7ZqiMhIYHBgwcTHh4OwNixY0lI\nSHDhT1c1Uj0lhKgSW1sGVK2qqdBaSPALwVV+nyd4a/VUbm4uzZs3x2q1Ur9+fQDy8/M5d+4cO3bs\noHv37mXeI9VTQgi/EBQgY4mravny5QQFBbFv3z527tzJzp072bdvH9dddx2LFi0yJCZ3J43pwF5g\nN2qp11AgAlgDHABWA41KnX8Q2A8MdnNsQogaCA8Jr/J7Xh/yOoDMQ+WkRYsWMXHiRGJjY4mMjCQy\nMpKoqCimTJnC4sWLsVqtHo/JndVTccCPQCcgH1gKfAd0Ac4A84BpQGPgGaAzKrH0AVoAa4H2QOnf\nilRPCWEAx2qpg48cpHXj1lWe2TZhRwL3rbiPd4a9wwO9HnB1iNXmrdVT1eHL1VNZqLXBw1BzXIUB\np4DhgK0VJwEYqe+PAJbo7zkGHAL6ujE+IUQ1tY1oW62p0P981Z8BePDbB10dkvAQdyaNDGA+cAKV\nLM6hqqWigDT9nDT9OUAMkOTw/iRUiUMI4SdM1e13KryGO1um2gB/RVVTnQc+B+4tdY6mb5dT7mvx\n8fH2fbPZ7Nfr8QrhDYqsRbQIb0FydnKNr7X1/q30ea8P6RfSaVavmQuiE+WxWCxYLBaXX9edSaM3\nsBE4qz//D9AfSAWa64/RwGn99WSgpcP7Y/VjZTgmDSGE+w3+ZLBLEgZAt8huAJzLOydJw41Kf6Ge\nNWuWS67rzuqp/UA/oC6q8eUmIBH4BhivnzMesE19+TUwDggBWgHtgC1ujE8IUYHcS7n2RZV+PPqj\n/bhtkF512RZuWnloZY2uI4zhzpLGTmARsA3VA2o78C4QDiwDJqEavG1TYybqxxOBQmAyFVddCSHc\nZNqaaczbOA+AzX/ZbD9+4dkLhAWXXWipOqaumkq94HpM6jnJJdcTnuGLrVLS5VYIN3PsXjuq4yi+\n2v8V7w57l/t73e/y63vD6HDpcus8GREuhKjQV/u/AmB8j/GVnOm8/9z5H7pHdWdM5zEuu6bwDBnX\nL4Qoo19sPzYlbbI/H9N5DCGBIS67/qhOoygoKmDcl+O4VHSJ4MBgl11buJeUNIQQZZw8f5Ljfz1u\nb7/4IvGLSt5RdeGhahqSCSsmuPza/iIuLo6wsDDCw8OJiIhg2LBhJCUlVf5GN5KkIYQooaCogPSL\n6cSEx5D6RKrb7jOkzRAAPt39qdvu4etMJhPffvst2dnZpKSkEBUVxSOPPGJoTJI0hBB2mqZx95d3\nU1BUQFBAEOGh4TzQ8wE2TNzg8nsFBgS6/Jr+LDQ0lNGjR5OYmGhoHJI0hBB2C7Yu4Mt9X5Y49s7t\n7zCg5QC33G/HgzsA/KbnkjvYfjcXL15k6dKl9O9f9QWwXEkawoUQdr+l/ubR+7Vv0h6ARTsXubR3\nlqs5dhGuiap2L9Y0jZEjRxIUFMSFCxeIjIxk1apVLomluiRpCCHsIutFAnBu2jmP3K9OUB0AjmQe\n8cj9qsuosSQmk4kVK1YwaNAgNE1j+fLlXH/99SQmJhIVFVX5BdxAqqeEEHaapvH3G/5OwzoNPXI/\n26y3s3+e7ZH7+TKTycSoUaMIDAxkwwbXtzE5S5KGEMJu7oa5nDh/wqP3bFK3CQARL0d49L6+wtam\noWkaK1asIDMzk06dOhkWj1RPCSFKaBDawKP3S30yleAXgsnMy0TTNFlzo5Tbb7+dwMBATCYTcXFx\nLFq0SJKGEMJ4eYV5ALx444sevW9QQPGfoe8Ofsdt7W/z6P292dGjR40OoQypnhJCAHDm4hliwmMM\nmdLDOsMKQNqFtErOFEaTpCGEACAzN5OcghxD7m2rkpr09SQGJQwyJAbhHEkaQggAZq2bRVZ+ltFh\n8NOxn4wOQVRAkoYQAqDMSHBP2/PwHkPvL5zj7qTRAfjNYTsPPApEAGuAA8BqoJHDe6YDB1HLxQ52\nc3xCCAejOo4y7N5dIruwd/JeAOZtmGdYHKJinuzbFgAkA32BR4AzwDxgGtAYeAboDCwG+gAtgLVA\ne9RysTaycp8QLjbn5zk8/9PzFP6t0PCJBG1TdnhyFLas3Oc8T3a5vQk4BJwEhgPX68cTAAsqaYwA\nlgCXUOuHH0IlmU0IIdzGNqDP6IRhlMaNG/vN+JDGjRu79fqebNMYh0oIAFGArW9dmv4cIAZwXGEk\nCVXiEEK40bvb3zU6BLs5g+Z4/J4ZGRlomuYXW0ZGhlt/V54qaYQAt6OqokrT9O1yyrwWHx9v3zeb\nzZjN5ppFJ4TwGg/1fojnfnyOImtRrS35uILFYsFisbj8up5KGrcCvwLp+vM0oDmQCkQDp/XjyUBL\nh/fF6sdKcEwaQojq23ZqG7mXcgF46caXDI5GiagbQWS9SFJyUohtEGt0OD6r9BfqWbNmueS6nqqe\n+iPFVVMAXwO2yfPHA8sdjo9DlUxaAe2ALR6KUYha58ZFNzLwo4FA8fKr3uD0hdN8mWhsF2BRPk8k\njXqoRvD/OBybC9yM6nI7SH8OkAgs0x9XApOpuOpKCFEDjoP5OjfrbGAkZf31+79yNNP75l6q7Xyx\nu4B0uRXCBcwfmVl3fJ39uVELDZVn0opJfLDjA8C74vJlrupyKyPChailHBNG4uREAyMpq2lYU/t+\n2Jww+wy8wniSNISohaaunApAgCmA9KfS6dTMuPUZyvP8wOe5u9vdAOQW5rLntEwx4i0kaQhRC725\n5U0Adj60s8S3em8RHhrOJ6M+sT/3hokUhVLVpBEBdHdHIEIIz4usF2l0CJdlMpn49YFfAdXLS3gH\nZ5LGOqABKmH8CrwPvObOoIQQ7veXq//i1UkDoGd0T6NDEKU4kzQaAlnAHcAi1FxQN7kzKCGE+9gG\n870y+BWDI3HOgqELjA5BOHAmaQSiRm3fCfxXPyZ94ITwUbvSdgHQILSBwZE458bWqmpqU5LvzVs6\nccVEMnLdOxeUpzmTNGYD3wOHUaOz26DWuxBC+KB+C/sBqueUL2jTuA1QXELyFZqm8eGOD2kyr4nf\nTLsOlSeNQNRcUN2Bh/Vjh4HR7gxKCOFe111xndEhOM02aeGHOz7k9iW3GxyN8zLzMu37W5L9Zzak\nyiYsLELNG/WqB2IRQrhZkbUIgKHthhocSdV9vOtjAC4VXSI4MNjgaCr31Oqn7Pv+1KDvTPl0PfAv\n4Dqgp8MmhPAxtvr1Jwc8aXAk1ff0mqc5n3fengC9VXhoOAB7J+/1iSTnLGeSxtVAF1TbxnyHTQjh\nY55e+zQAQQGeXLSz5i4+e5Evxn4BwKbkTTR6uRFz18+t5F3GahHegsf7Pe51E0HWlExYKEQtYsT6\n265kix/U/FSnHj/ltd/iG85tSI/mPVh337rKT/YAT09YOAx4GpjhsAkhfIg/fNma0GOCff/MxTP0\nX9jfwGgqlpWfRaDJ/1YedKaM+g5QF7XuxXvAWGCzO4MSQrjegq1qkNzeyXsNjqT6Fty2gDaN23Ag\n4wCLdi7i15RfjQ6pQvf1uM/oEFzOmZLGAODPQAYwC+gHdHBnUEII15q3YR5TVk4BvG+xpaqoE1SH\n5wY+R8LIBPuxDSc2GBhR+ayaFYBWjVoZHInrOZM0bCNqLgItgELU+t7OagR8AexDrch3DWoeqzWo\nlftW6+fYTEcNHtwPDK7CfYQQlzFt7TSjQ3C5wr8VAjBrnWvWvnalMxfPANA7prfBkbieM0njW6Ax\n8ApqwsJjlFzvuzJvAN8BnVCDBPcDz6CSRnvgB/05QGfgLv3xFmCBkzEKIZwwrus4o0NwGdugvzVH\n1pRoIPcGS/csBaBucF2DI3E9Z6cRyQS+BOKAjsDfnLx+Q9T4jg/054XAeWA4YCtfJgAj9f0RqIR0\nCZWcDqEmSBRCVNP5vPP2/RvibjAwEtf75o/f2Pe9adxG24i2dGjin7X4FTWEj6Z4YkITZScp/I8T\n128FpAMfAlehSip/BaKANP2cNP05QAzgOCtZEqpKTAhRTetPrAfAOsNq63bpN4a1H0bKEylEz48m\n6IUgmoU14/RTp40Oi7O5Z/1qFLijipLG7VQ8m60zSSMINXp8CrAVeJ3iqigbrZL7lHktPj7evm82\nmzGbzU6EIkTtNGzJMAC/Sxg2zesXN7GmX0w3MBKl73t92XpqK4/0fcTQOCwWCxaLxeXXrShp3OeC\n6yfp21b9+Reohu5UVGN6KmraddtXg2TUBIk2sfqxEhyThhBC3N3tbhbvXmx0GOQV5rH1lPpz16Ru\nE0NjKf2FetYs13QYcKZNYyZqMJ/tsSqD+1KBk6gGb1CLN+0FvgHG68fGA8v1/a+BcUAIqmqrHWo6\ndiFENY3uNJrXh7xudBhu9ekdn/LzfT8DEPxCMEcyjxgSx5rDa+z7TcKMTRru4kzSuKBvOYAVGIpq\nEHfWI8CnwE5U76k5wFzgZlSX20H6c1BdcpfpjyuByciCT0LUyJf7vvTLXjylXXelmu690FpImzfb\nEDTb8/NrDf9suH3/P/ucqcH3Pc78Vv9R6vkrqLEVztoJ9Cnn+OWWjH1R34QQNXT24lkAOjbtaHAk\nnlekFdHt393Y/fBuj92zSd0mnM1Vv/PJfSZ77L6eVJ0xEPWQHk1C+ISmrzQFoEfzHgZH4hlr/7SW\n+OvjMenz8u05vYeLly567P5nc8/ycO+H2fXQLsZ0HuOx+3qSM0ljt8O2F/gdNWBPCOEjfGU98Jq6\nsfWNzDTPxDrTyjvD3gHw2BrdeYV5AETXj6ZbVDeP3NMIziSN2x22IaixFP90Z1BCCNdZOHyh0SEY\n4i89/wLA/47/z6330TQNq2YlLUcNPXv2umfdej+jOZM0mqMmKzyG6j5bFzV/lBDCi+UX5hMcEOyX\nM606I8Ck/rzd/Z+72Z22m1c2vOLye2iaRsDsAAJnB7LvzD56RfeyT2/ir5xJGm+jek7ZXNCPCSG8\nWGpOKlH1o+x/PGujsOAwALq/3d2+aqEzTp4/Sas3WnHw7MEKzwuYXfy7vfXTW4mqH1XB2f7B2U+T\n1WG/CPDvVCqEHziVfYqY8BijwzDU2afPlnienJVMlwVdyj33VPYp3tz8JqZZJp778TmOnTtG+3+1\n54V1L5Cdn11mEatLRZfKXOO7g9+5Lngv5UzSOAo8CgSjBt1NBYwZOSOEcNoLP7/AjtQdRodhqDpB\ndUo8j30tlsT0RHtXZJs1h9fQ4tUWTF01FYCPd31sf22GZQYN5jZg4EcD+S3lN/vxa95XtfTnpp1j\nzqA57voRvI4zSeMh4FrUdB5JqEWYHnBnUEKImlt5aCUFRQVGh2G49RPWlznW9JWmvLn5Tfvzd7e/\nW/l1Tqyn57tqEsJtp7bxW6pKIA3rNOTZ655l10O7yH0ut6JL+AVnkkYaao2LSH37I8VzRQkhvNiH\nIz40OgTDXXvFtWgzNU49fqrEcVupArD3fLq17a28fNPLAMwYOINt928rc721R9bS5z01Xvnm1jfb\nj3eL6lamZOOPnJn2sgNqMaTmQBfUVCDDgb+7Ma6KaKXrFoUQJV0qukS9F+uR+1yu3/fmqYrxy8fz\n12v+ai8xaDM1Nidtpt/CfiwYuoCH+zyMVbOSV5hnb0TfeHIj/WP7k1OQQ4O5Jce7HJhygHZN2nn8\n56gOfZbjGk917ExJ4z3gWcBWzt2NKm0IIbxUSk4Kzeo1k4RRSsLIBK6OvpqdD+0kJDCEQmsh/Rb2\nA2B059GA6qprSxgAA1oOwGQyER4aXmL51od7P+wzCcOVnEkaYcBmh+caamU9IYSXOpV9ihbhMtvP\n5UTUjaCgqIDgF4LtxyLrRVb6vq33byU8JByAt4a+5bb4vJkzExamA20dno8BUtwTjhDCFdJy0kos\nTiRKim0QW+J51jNZTr83a7rz5/ojZ5LGFOBd1Nrgp1BdcO9xZ1BCiJpZsmeJfZlXUbGC5wsIDgyu\n/EQBVK1RpB6qOivbTbE4SxrChaiEaZb6r63NlP8rl1NoLWRf+j6/nlzQkasawisqaTzhsO/4yTPp\nz1918h7HgCzUSPJLQF8gAlgKXKm/fidwTj9/OjBRP/9RqrZ2hxBC9+wf/HvivJoKCgiqNQnDlSpq\nCA8H6utbeDmbszTADFyNShgAzwBrUMvA/qA/B+iMGhPSGbgF1dW39k6cI0Q12EriXSLLny5DiJqo\nqKQR78L7lC4SDQeu1/cTAAsqcYwAlqBKJMeAQ6hEs8mFsQjh135N+RWAu7vdbXAkwh858y2+JfAV\nqhdVOvAlEFvhO0rSgLXANuB+/VgUaqQ5+qNtasgY1FQlNknIKoFCVMnutN00C2tmdBjCTznTe+pD\n4FNUuwOonlMfAjdf9h0lXYvqotsMVSW1v9TrGiXbTEqTljwhquCx7x/jfP55o8MQfsqZpNEMlSRs\nPgIeq8I9bGM60lEllr6o0kVzIBWIpnguq2RUycYmVj9WQnx8vH3fbDZjNpurEI4Q/m3qNVMrP0n4\nPYvFgsVicfl1nel+9SMqaSzWzx8HTABudOK9Yai1N7JRXXZXA7OAm4CzwMuotoxG+mNn/T59UdVS\na1EDCx1LG9LlVogKmGaZuK3dbXx797dGhyK8iCe63NpMRK0JbutiuxGVNJwRhSpd2O71KSpxbAOW\nAZMo7nILkKgfTwQKgclI9ZQQVdasnrRpCPeocdYxgJQ0hKhA1wVdWTJ6iYxBECV4oqTxT9S3/PJu\noqEG3gkhvEh2fjZ70/cSHR5tdCjCT1WUNB4C9qCqi2yrl9gSiHzVF8ILHT9/HIDGdRobHInwVxUl\njWhgLKq9oQg17cfnFE/3IYTwMkt2LwGQdTSE21Q0uO8M8G/gBuA+oCGqgfpP7g9LCFEd2QXZRNeX\nqinhPs6MCO8FTAXuBVYCv7o1IiFEtf1zyz+dWkxIiOqqqHrqBWAosA/4DLXkq6zYJ4SXmzNojtEh\nCD9WUdJ4DrXg0lX69pLDaxrQ3Y1xCSGqKTUn1egQhB+rKGm09lgUQogayynIAeBPV0mzo3CfipLG\nMU8FIYSouWPnjgEQEhhibCDCr8kCR0L4iW7/lhHgwv0kaQghhHBaRUnjB/1xnicCEULUnEyLLtyt\nshHhA1BLs36GmkLEcfqQ7W6MSwhRDTfE3WB0CMLPVZQ0ZgIzUOtazC/ndfl0CuElrJqV4IBgbml7\ni9GhCD9XUdL4XN9mALM9E44QojoycjOoH1Kf0KBQo0MRfs6ZRZhmAyOAgajqqXXAN+4MSghRNefy\nztG4rsxsK9zPmd5Tc1FrZ+xFTSnyKCVHh1cmEPiN4kQTAawBDqBW8WvkcO504CCwHxhchXsIUatZ\njlk4knnE6DBELeDMKk67gR6o6dFBJYEdgLOdwh9HTXoYjmpUn4eaQXceMA1oTMn1wftQvD54e8Ba\n6nqycp8QpZhmqf/K2kz5vyHK56qV+5wpaWiULA00wvlFmGJRkx6+T3Gww4EEfT8BGKnvjwCWoCZF\nPAYcAvo6eR8har06QXWMDkHUAs60abyE6l77E+oP//WokoEzXgOeAho4HIsC0vT9NP05QAywyeG8\nJFSJQwhRiV7RvXhr6FtGhyFqAWeSxhJU43cfVAnjGSDFifcNA06j2jPMlzlHo+JSi5S1hXBCYnoi\nLRu2NDoMUQs4kzRArRG+oorXtg0MHArUQZU2PkaVLpoDqagBhKf185MBx099rH6sjPj4ePu+2WzG\nbDZXMTQh/MPpC6c5cf4EuYW5RNWLqvwNotawWCxYLBaXX7fGjSJOuh54Ergd1QB+FngZVWppRMmG\n8L4UN4S3pWxpQxrChdDZGsBBGsFFxVzVEO5sScMVbJ/oucAyYBKqwftO/XiifjwRKAQmI9VTQlzW\n6sOrjQ5B1EKVJY0g1PiMDjW8zzp9A8gAbrrMeS/qmxCiEvvP7Lfvfzb6MwMjEbWJM0WVFagBfcfd\nHIuzpHpKCOCmRTfRNbIr8wfPJzAg0OhwhJfzZPVUBKq0sQW4oB/TUI3cQggDZORm8MPRH7BqVkkY\nwqOcSRp/K+eYfNUXwkBbk7cC8PGojw2ORNQ2zhZV4lA9mdYCYahkk+WmmCoj1VOi1pNpQ0RVeXIa\nkQdQU6S/oz+PBb6q6Y2FEEL4HmeSxv8Bf6C4ZHEAiHRbREKICqXmpAKw/QFZPFN4njNJI1/fbIKQ\nNg0hDLPy4EoAukU5O9G0EK7jTNJYBzyHasu4GVVVJYswCWEAq2Zl4tcTAQgK8OTYXCEUZ5LGM0A6\nal2NB4HvgOfdGZQQonyd3+oMwDvD3qnkTCHcw5mvKkWodS82o6ql9iPVU0IY4uroqwkNCuWBXg8Y\nHYqopZxJGrcBbwO2tSRbU1ziEEJ40InzJ3h9yOtGhyFqMWf67P6OShyH9OdtUAmjpvNRVZeM0xC1\nUlZ+Fg3nNiTj6Qwa121sdDjCx3hynEYWxQkDVInDqIF9QtRaSVlJhIeES8IQhqqoemq0/rgNVbJY\npj8fqx8TQnjQjJ9mkF2QbXQYoparKGncTnGD92nUQkqgelLJCvZCeNDO1J18ue9Lo8MQwmMr97mS\ntGmIWqfunLrkFeax9f6t9I7pbXQ4wgd5sk2jNfAaar6pb/TtayfeVwfVTXcHajW+l/TjEcAa1HQk\nq1HLvdpMBw6iuvUOduIeQtQKPaN7AtArupfBkYjazpkut8uB91HJwqofc+arfh5wA3BRv8961BxW\nw1FJYx4wDTV40LZG+F36o22N8PYO9xSiVtqavJWNJzfy1tC3bN8WhTCMM0kjD3izmte/qD+GAIFA\nJipp2NpHEgALKmmMAJYAl1Brhx8C+gKbqnlvIfxC3/f7ApCclWxwJEI4Vz31TyAe6A/0dNicvf4O\nIA34CbUCYJT+HP0xSt+PAZIc3puEKnEIIYC+LfoaHYIQTpU0ugB/QlU1OVYV3eDEe61AD6Ah8H05\n79GouKqr3Nfi4+Pt+2azGbPZ7EQoQvieKd9Nse8Paz/MwEiEr7FYLFgsFpdf15kK0sNAJ6Cghvf6\nG5AL/AUwA6lANKoE0hFVRQUwV39cBcxENaY7kt5TotaQFfqEq3iy99RuoDpDUJtS3DOqLmpa9d9Q\nPa/G68fHoxra0Y+PQ7V/tALaAVuqcV8h/IJ8ORLeyJnqqcaoLrBbKV6MSUM1aFckGtXQHaBvHwM/\noBLHMmASqsH7Tv38RP14IlAITEZm0xW12PHzxwG48OwFgyMRopgzRRXzZY5bXBdGlUj1lPAr8zbM\nY9raaTzY60HeHva2/bhUTQlXclX1lC92+pakIfzG94e+55ZPb7E/12ZqrD+xnkMZh5iwYoL9mBA1\n5aqk4Uz1VA7F1UQhQLB+rEFNby5EbZaUlWRPGI3qNOJc3jl76cJm2Zhl5b1VCMM40xBeHwjXt7rA\nHcACdwYlRG3Q8rWWANzY6kYyp2WWe87YLmM9GZIQlXImaTiyono73VLZiUKIyzufd96+nzAyAZBq\nKOEbnKmeGu2wHwD0Qo23EEJU09FzRwFIfyqdpmFN7ce1mRqF1kICTAGYfLLJUfg7Z5KG47oahahu\nsiPcFZAQtcGO1B0MbTe0RMKwCQpw5r+lEMbwxa8y0ntK+DzpTis8zRO9p2Ze5rjtUz67pjcXoraK\nCY/hpRt1Y+geAAAWyklEQVRfqvxEIbxMRUnjAmVHZNdDjeRuiiQNIaqlyFrEqexT9I/tb3QoQlRZ\nRUnjHw77DYBHgQnAZ8B8dwYlhD87knmE0MBQ2jVpZ3QoQlRZZS1uTYDHgHuARah1NMrvUC6EcEr7\nf7U3OgQhqq2yksYo4F2gO5DtkYiE8HNNw5oS2yDW6DCEqJaKWtKtqDU0LpXzmoZx04hI7ynhszRN\nI+6NOJaOWUq/2H5GhyNqEU/0nqrqaHEhRCVWHVrFifMnaBvR1uhQhKgWSQxCeNDJrJMMaz+s3EF9\nQvgCSRpCeNCGkxvo10KqpYTvcnfSaIlaA3wvsAfVbRcgAlgDHABWU7wsLMB04CBqtcDBbo5PCI/J\nzs9m0c5FZOVnGR2KENXm7mlEmuvbDtQU678CI1HjPc4A84BpqCVlnwE6A4uBPkALYC3QHtUobyMN\n4cIn2aYOyXsuj9CgUIOjEbWNqxrC3V3SSEUlDFALN+1DJYPhqPXD0R9H6vsjgCWoHlvHgENAXzfH\nKIRHtItoR4vwFpIwhE/zZJtGHHA1sBmIAtL042n6c4AYIMnhPUmoJCOEz2sT0abEGuBC+CJPzcFc\nH/gSmErZQYIaZee4Kv16CfHx8fZ9s9mM2WyucYBCuNuhjEO0adzG6DBELWGxWLBYLC6/riemRg8G\nvgVWAq/rx/YDZlT1VTSqsbwjql0DYK7+uAo12+5mh+tJm4bwSeEvhZP0WBIN6zQ0OhRRC/lKm4YJ\nWAgkUpwwAL4Gxuv741FLyNqOjwNCgFZAO2CLm2MUwu2Ss5IJMAXQINSoiRSEcA13V09dC9wL7AJ+\n049NR5UklqGmWT8G3Km/lqgfT0StEjiZiquuhPAJnd7qRHZBtu3bnhA+yxc/wVI9JXyOaZaJqHpR\npD6ZanQoopbyleopIWq9n47+BMBHIz8yNhAhXEBKGkK4mW1QX9GMIgJM8j1NGENKGkL4ANuUIfv+\nb58kDOEX5FMshBvFvR4HQMemHY0NRAgXkaQhhBtl5snqyMK/SNIQwk1Onj9J4zqNufDsBaNDEcJl\nPDWNiBBea+XBlWhoDG03lLzCPIqsRdQLqVeja1734XV0btqZgVcOJCw4zEWRCmE8SRqi1hu6eGiZ\nY60bt+bwo4erdb2CogLWn1jP+hPrmW2eXdPwhPAqUj0l/NrhjMMMWzyM/x74b7mvT187vcRzW4P1\nkcwjXCq6RGpO1QbjJWclE/r34qnPe8X0qmLEQng3Gach/FZ2fjYN5hbP9TShxwQ+GPFBiXO6LuhK\np2ad+Hzs5yWO28ZWlCfliRSa129e5nhWfhYN56rJCFfds4obWt1ASGBITX4EIVxGxmkIUQFN03ju\nx+cAWPOnNQB8uONDOr/VucR5py+cZv7g+WXe/8ukXy577Z7v9GT5fjXHZkp2CjkFOQD2hPGvW//F\nkLZDJGEIvyQlDeGX/r3130z+brK9bULTNAJmq+9ILw56kdCgUJbtXcbOtJ1cfPZiuRMJ7kzdSfeo\n7phMJi4VXaKgqIDBnwxm48mNAGyatIl+C/sBMGPgDGb/PLtGbSFCuJOrShqSNIRfaji3IVn5WSRO\nTqRTs04AFFmLCHqhbN8PbabznyerZmXP6T0M/XQoydnJZV4veL6A4MDg6gcuhJtI9ZQQl2GaZSIr\nP4uNEzfaEwZAYEAgRTOKuKvLXQAMvHIgn43+rErXDjAF0D2qO38b+DcAtt6/lbdve5sn+z9J3nN5\nkjCE35OShvArXyR+wdjPxwJVK0FUh6Zpsj6G8BlS0hCiHAu2LgDgq7u+cvu9JGGI2sjdSeMDIA3Y\n7XAsAlgDHABWA40cXpsOHEStIT7YzbEJH/ftgW8xzTLZt/QL6ZzKPsXW+7cysuNIo8MTwi+5+6vS\ndUAOsAjoph+bB5zRH6cBjYFngM7AYqAP0AJYC7QHrKWuKdVTgk93fcq9X91b7mvWGVYpBQhRiq9U\nT/0PKD3N53AgQd9PAGxfCUcAS4BLqHXDDwF93Ryf8FG2hHFgygFOPX6KpWOWEmAKYNU9qyRhCOFG\nRsw9FYWqskJ/jNL3Y4BNDucloUocQthpmsYdy+4gwBRA0mNJRIdHA3Bnlzu5s8udBkcnhP8zesJC\nTd8qer2M0ZMm0K3llQCYzWbMZrPrIxNeafSy0Szfv5y2EW3tCUMIUZbFYsFisbj8up4ox8cB31Dc\nprEfMAOpQDTwE9AR1a4BMFd/XAXMBDaXup5GPCTek0OntjWbvlr4lrScNJrPb86AlgNYec9KGoQ2\nqPxNQgjAd9o0yvM1MF7fHw8sdzg+DggBWgHtgC3lXSCisAudP63PgbMH3Byq8AbHzx3HNMtE8/nN\n6di0IxsmbpCEIYRB3J00lgAbgQ7ASWACqiRxM6rL7SCKSxaJwDL9cSUwmctUT/0+dQ8AHf7Vgcxc\nWU7T3634fYV9/5NRnxgYiRDCF7uZaJqmMWasxuoO7bjrmkG8d/u7Rsck3OSnoz8xaNEgbm17K6/f\n8jrtm7Q3OiQhfFKtn7AwIwN6Dkzj+NjmrLxnJbe0vcXouIQLWDUrAz8cyP097+e+FfcBcG/3e0kY\nmUCASSYwEKK6an3SAHj7bZj64Ufc/Mhyvr13eSVvE95O0zSGfDKENUfWlDh++snTNKvXzKCohPAP\nkjSAoiLocU0We4c1YdWf/svgNjLziC+zrZZ3d7e7Sb+QztIxSynSimga1tTgyITwfa5KGkaP06iR\nwEBY+98GxAz6nKGFY3l/0LfcN+g6o8MS1bDmsCpdvDr4VR7r/5jB0QghLsfnK4mjoiDVMpIxDV5l\nwurhTHnrKzIyQKan8h3/3vpvBn+iSomSMITwbj5dPVXahH8s46MLd4FlJl1yH2LT2ubUr+/h6ESV\nvPrLqzyx+gkAUp9IJap+VCXvEEJUh7RplMOqWdmRuoNe7/YCIC75KTbOfono5oGejE+UI/1COvlF\n+cQ2iAVUo3fv93qzPWU79UPqc+TRI9LYLYQbSdKoQEFRAT/8/gtDP/gj7P4jYevn88YbMH48BMtq\nnB53Lu8cXRd0JTk7mXFdx/HZnuIlVt8a+hYP9nqQwABJ7EK4kyQNJxxNS6fXRx1ok/F/bJs3GzDR\nty+cPAnz5sG95S/HIFwkpyCHsZ+PZdWhVQDMHzzfXhU1utNomZlWCA+SpOGkg2cPcuOiG7m7y3i6\nnH2Wwry6nDkD06fDDz/Avn0wbhw0alT5tQQUWYsotBYSGhRa4vie03v4LeU37u52N4XWQkICQ3hi\n9RO8tuk1Fg5fyD3d7iE0KJT8wvwy7xVCuJ8kjSrYnrKdsZ+PJTkrmQd6PcD4q8bzwZxefPAB5OVB\n27bw7bfQoQOcPw/h4RDg8/3KXGfdsXUM/2w4cY3iOHbuGFn5WVwVdRU7HtoBwLDFw/jvwf8CEBIY\nQkFRAZOunsTHuz7m6NSjxITHGBm+EAJJGlV+k1WzEm+JJ78wn7d/fZuY8BgSRibQtXFfFi6E+Hi4\n6ir46Sd1/pgx8MUXEBkJffrA3r3w/PPQoAFcfTU0b85le2bl50Ooj36ZvlR0iZScFL5M/JLHVz9O\nXKM4krKSMMeZia4fzU/HfuKm1jfx0Y6PAIgJj+FU9ik+GP4Bjeo04o5ld9ivNfP6mcSb4435QYQQ\nJUjSqIHs/Gze3vY2//jlH0zpM4XhHYajpXXjf7/k0aF1GMF189h79Cx7Egu5tf+VZGXBrFmQng5Z\nWcXXmTEDTCawWuGFF2D4cEhKgu3b1etTpqi2k9xciIgoGcOmTep9cXHQrBkUFqprff455OSoxLR+\nPUycCGfOwI8/QmYmzJkDFy7Ae+/Bn/8MISE1+lUAkJGbwcnzJ/lq/1e8tP4lCooKaBjakP4t+/Pa\nkNdo36R9mXmfMnMzWXloJb1jehNVL4qGdRraX9M0je0p2+kV06vmwQkhXEKSRs0vwnvb3+Pn4z+z\n4vcV5BTkANAnpg9bT20tce7/JvyPHs17UC+4PoWFKnG89hqsWAF79hSf98QTMGqUKp18/LFKJDZd\nu8JNN6mSyuLFcOhQ+XENHKjaWdLTy742dapKNEFB6v4A116rkkmDBqp3WFAQ3HorHCvcRGr6JbKO\ntSU/7ms6Nu3A9VdeT7wlnk7NOhEUEMSmpE1sTt7MztSdZBdkc1eXu7izy53c3PpmwkPDa/LrFUJ4\nGUkaLnSp6BLpF9MJCgjCcsxChyYdaNekHTkFOTz07UN8tf8rAIIDghnWfhg3tb4JEyY0NP7QciCb\nTm5lwJV96BrZFYDE9EQ2J22ma5NebDy6nZd+iSezII06WlPqFMQSeHgo2R0X0DuuExdyCwk4OZC7\nbujKe7tfY1/2Fm5ufTPRYVfSrG5zzhyNYdiAtoy5+uYSMV+8CL/9BtOmqTm49iQdZej1Uew7ks3u\nyGeg7SowWaH+aYLOdaSw0f4S7+8X24/ukd25o9MdXB93PZqmUSeoju2DJYTwM/6cNG4BXgcCgfeB\nl0u97vKkUZkLBRdIzUnl6LmjrDu2jjMXzwCw5sgae7XNwYyDNAtrhlWzEhoUSmyDWM7nnad9k/bc\n1eUuosOjeWXjK1x3xXWcPH+S/Wf3c23La7lUdIlfU36lXkg91hxewyN9HyGyXiSJ6YlYsZKSncL3\nh78nsl4kA68cWCKusOAwekX34oejP/Dz8Z8ptBaiaRrdI3sQ33MhZ1LqkBWaSP7eW8jJy+PTnzez\nb01fwusFM3RIMM2awYgR0KoVNGmiqtauvFJ1BBBC+Bd/TRqBwO/ATUAysBX4I7DP4RyPJ43KWCwW\nevTrQUFRAQBN6jZx6WC1A2cP8PPxn2kYqtoNMnIzOJhxkCsaXsGBswe4ouEVTLp6ElbNypHMI/Rt\n0ReTyYTFYsFsNtuvU1gIO3fCsmVw+rSqWktJgXPnoKBAVW01bAgDBkCXLtC7t2pzAdXeEhEBLVpA\nWpragoPV+TExkJwMZ8+q15s0qfh35RiTN5CYnOONMYF3xuWNMfnrLLd9gUPAMf35Z8AISiYNr+Pu\nD0j7Ju2dXrHOcSqO0nEFBUGvXmpzVFCgtrAwWLdONbx//LFqewkLU+fs3q0eGzVS57ZurZJOSAi0\naaPaYerUUV2YmzSBwYOLSy+O9u2z0KmTmeHDVRuMN3Rt9sb/4BKT87wxLm+MyVW8LWm0QK0lbpME\nXGNQLLVGSEhxL6wbblCPY8eWPMdqVW0nOTkqOdStq55brbB6tXr/zTerEsiGDbBjh+rldccdJbsm\nf/aZSiivvQbvvw9N9aUyTCbV5blePdi2TXVbtomJUSWiygqYR4+qe73xBmRnqxhsTp9WiQ3UOT17\nqnuCGqPz66/qfMfecXl50KmTqrJzlJsLa9ao3nLl6d4dGjdW17ftX3ll+SWwoiLYtavsz5aSojo4\nNG5c8c8shKd5W9LwrnonYRcQoDbHP2KBgWq77bbiY1FRKlHccUfZa4CqHps0CUaPVqUam5QUOHBA\n/bHv1Qs6dlTHNU11Ye7RQ42NqUhuLnz0ETz4oPqD3bu3SkI2t94K0dFw/DgcO1Z8PDhYJcv4eDXQ\n06agADZvVn/YHeXnqwT3hz+UjSEzUyWBCxfg1ClISFDJ5/ffyx/Xk5Ojetu1aFHy+IkT8O67lx8L\n5EjTVNIeMKA4EbrD77+r5Fode/eqUqc7xi/l58P8+a6/blXl5RV/Tkv/rhYuVP/O/sDb2jT6AfGo\nxnCA6YCVko3hh4A2ng1LCCF83mGgbaVn+Zgg1A8WB4QAO4BORgYkhBDCu92K6kF1CFXSEEIIIYQQ\nQgj3ugXYDxwEpnnwvi2Bn4C9wB7gUf14BLAGOACsBhwnWJ+OinM/MNiNsQUCvwHfeElMjYAvUN2k\nE1G934yOaTrq3243sBgINSCmD4A0PQab6sTQS7/GQeANN8X1CurfbyfwH6Chw2ueiKu8mGyeQLVz\nOs7mZmRMj6B+V3so2fZqVEx9gS2ovwlbgT4ejsmrBKKqrOKAYDzb3tEc6KHv10dVn3UC5gFP68en\nAXP1/c56fMGoeA8B7hqR8DjwKfC1/tzomBKAifp+EOoPjpExxQFHUIkCYCkw3oCYrgOupuR/8KrE\nYOu0sgX1hwHgO4o7jbgyrpsp/pnnGhBXeTGB+vK2CjhKcdIwMqYbUEnfth6obZCUkTFZgCH6/q2o\nL7uejMmr9Ed9YGye0TcjLEeNWt8PROnHmuvPQWV0x5LQKlTPMFeLBdaiPry2koaRMTVE/YEuzciY\nIlBJvjEqiX2D+qNoRExxlPwPXtUYoik50HUc8LYb4nI0CvjEgLjKi+lzoDslk4aRMS0DBpVznpEx\nLQFsy2H+ETf823nBeFynlTfwr8VlznWnOFR234z6D5+mH0+j+A9ADCo+G3fF+hrwFKq4bmNkTK2A\ndOBDYDvwHlDP4JgygPnACeAUcA717dDofzuqEUPp48lujM1mIurbp9FxjdDvsavUcSNjagcMBDah\nvuH39oKYnqH48/4KxZ2JXBaTLyUNbxj4Vx/4EpgKZJd6TaPiGF0d/zDgNKru8nLjbTwdUxDQE1ig\nP16gbGnQ0zG1Af6KSvYxqH/D0qvDezqmy93DGz7jjp4DClDtQEYKA54FZjoc84YxZkGoEmw/1Je3\nZcaGA8BCVJvrFcBjqHYPl/KlpJGMqtO0aUnJDOluwaiE8TGqegrUt0PbOOVo1B9xKBtrrH7MlQYA\nw1FF9SWoYvLHBseUpG+2BUm+QCWPVANj6g1sBM4ChaiG3f4Gx2RTlX+rJP14rIdiuw8YCtzjcMyo\nuNqgkv5O1Oc9FvgVVTIz8neVhPo8gfrMW4GmBsfUF/hK3/+C4rYKb/hMeZyRA/9MwCJUdZCjeRTX\nEz5D2QbDEFSVzWHc+83oeorbNIyO6WfANrtivB6PkTFdherZUle/dgLwfwbFFEfZhvCqxrAZ1SPN\nhOsaLUvHdQuqt1nTUud5Mq7SMTkqryHciJgeBGbp++1RVUJGx7Qd9fcA4EaKv8B5+jPlNYwa+PcH\n1LeIHajqoN9Qv9gIVEN0eV0mn9Xj3E9xbwZ3uZ7i3lNGx3QV6oPq2F3T6JieprjLbQKq1OjpmJag\n2lQKUG1zE6oZg6175CHgTTfENRHV9fI4xZ/1BR6OyxZTPsW/K0dHKNnl1qiYglGl+92oko/ZoJgc\nP1O9UUlgB/ALqv3VkzEJIYQQQgghhBBCCCGEEEIIIYQQQgghhBBCCCGEEE0oHpOQghop+xtqqph/\nufA+/6BkX/7SHgX+5ML7CSGEcLOZqOnmXS0cNQV1Tc8RwlC+NPeUEJ5im17BTPH0LPGo0eQ/A8eA\nO1Alh13AStQ0N6BG11qAbajpp21zS41Ajf62mYsapb4TNRspqJLNWaCLy34SIVxMkoYQzmuFWrtk\nOGqdgjWo9R1ygdtQ00r8ExiNms7hQ2CO/t4/oBIJqOqwkajkcBXwd4d7bEFNty2EVwqq/BQhBGra\n8pVAEWoCxADge/213aiJ49qjEoGtRBGImhsI1FTVKfr+eSAPNY31t/pmcwpo7Y4fQAhXkKQhhPMK\n9EcrcMnhuBX1f8mEqnIacJn320r2hagpq28ExgBT9H30a3jbuhpC2En1lBDOcWZ69N9R60TblocN\nRk1JDWrWWFv7Rj3UjLYrUY3uVzlcIxrVZiKEV5KkIURZmsNjeftQtjSgoUofY4CXKZ5Gv7/++nqK\nlwNtgGpg3wn8D7XCmk1f/ZgQQoharD7FC+JcTgMnzhHCUIFGByBELVEAtEO1Zxy7zDkPorrr7vJM\nSEIIIYQQQgghhBBCCCGEEEIIIYQQQgghhBBCCCFc7P8B/5kYTvfuSuMAAAAASUVORK5CYII=\n",
"text": [
"<matplotlib.figure.Figure at 0x7efedbe69e50>"
]
}
],
"prompt_number": 16
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"dsolve([100,100], 3000)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEACAYAAAC08h1NAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGMRJREFUeJzt3XmQnPV95/F3z6XRMadGGt0MkSXAIAy2BS5zuCHCBh9Y\n5dQ6pJyUyjYuJ95gnMsI78GkNpsAicvG2bKrNo7ZIXGcVUgg4MQbCZmGNV4wBAQ6LIQEQueMjtFI\no5FGc/X+8TyjGc0hafp6+ni/qrr6eZ7p5+nvTw98+unf8zz9A0mSJEmSJEmSJEmSJEmSJOWxHwAd\nwOZRyxqBDcAOYD1QP+pv9wNvAduBj+aoRklSBt0EXMu5wf8w8PVw+j7gwXD6vcAmoBJoAXYCZTmp\nUpKUUS2cG/zbgeZwel44D8HR/n2jXvd/gA9luzhJ0tSkckTeTND9Q/g8/CGwANg36nX7gIWplyZJ\nyoZ0u2KS4eN8f5ck5ZGKFNbpIOjiaQfmA4fC5fuBxaNetyhcdq56knSl8K6SVNp2Ae/JxIZSOeJ/\nClgTTq8Bnhy1/C6gCrgUWAb8YtzaXdDReZJkMlmUjwceeCDyGmyf7SvF9hVz25LJJMDSFPJ6QhcK\n/h8BPwcuA/YCnye4iuc2gss5b2Xkqp5twLrw+SfAV5ikq6f9WHe6dUuSUnShrp7fmGT5qkmW/2n4\nOK/Dx09e6CWSpCyJ5Dr7Q8eL94g/Ho9HXUJW2b7CVsztK+a2ZVosgvdMfuefn+OeO2+O4K0lqTDF\nYjHIUGZHcsR/tLt4j/glKd9FEvzHTtnHL0lRiST4u055xC9JUYkk+I+fNvglKSqRBP+JMwa/JEUl\nkuA/2WcfvyRFJZLg7xnwiF+SohJJ8J8y+CUpMpEEf++QXT2SFJVIgv9M0iN+SYpKJMHfFzP4JSkq\nkQR/v8EvSZGJJPgHy+3jl6SoRBP8FR7xS1JU0gn+e4HNwJZwGqAR2EAwOtd6oH7CNau6GRpyHHZJ\nikKqwX8VcDewEngf8EmC8SDXEgT/cmBjOD/eUCWd3adTfGtJUjpSDf7LgZeAXmAQeA74NeBOoC18\nTRuweqKVY/01dByzn1+SopBq8G8BbiLo2pkBfBxYBDQDHeFrOsL5ccoHaujosp9fkqJwocHWJ7Md\neIigH78H2ERw5D9aMnyMk/xZD9898uc83zKPeDzuWJmSNEYikSCRSGRl25kac/e/A/sITvLGgXZg\nPvAsQbfQaMmaez/Mg6se4iufvDFDby9JxS1fxtydGz4vAT4D/B3wFLAmXL4GeHKiFatiNRxx3F1J\nikSqXT0AjwOzgX7gK8Bx4EFgHfBFYDfw2YlWrI7VcKzH4JekKKQT/DdPsKwTWHWhFaeXzeKY4+5K\nUiQiuXN3RkUNJ3q9nFOSohBJ8M+qquFEr0f8khSFSIK/pqqG7j6DX5KiEEnw10+vpbvveBRvLUkl\nL5Lgb5xZx6nBE1G8tSSVvEiCv2lWHaeGPOKXpChEEvxzamvpTRr8khSFSIK/ub6O/jK7eiQpCpEE\n/4LGOvrLPeKXpChEEvzzG2sZrDD4JSkKkQT/wqZah1+UpIhEEvzVVRUwMJ12R+GSpJyLJPgByvrr\n2H/E7h5JyrXIgr9ioJYDRw1+Scq1yIK/MllHx3Ev6ZSkXIss+KuTdRw67hG/JOVaOsF/P7AV2Eww\n7OI0oBHYAOwgGIi9frKVp5fVcbjb4JekXEs1+FuALwHvB1YA5cBdwFqC4F8ObAznJzSjvJZjPXb1\nSFKupRr8JwjG2p1BMHzjDOAAcCfQFr6mDVg92QZqKuvoPO0RvyTlWqrB3wl8E9hDEPhdBEf6zUBH\n+JqOcH5CtdPqON5r8EtSrqU62PpS4GsEXT7HgX8AfnPMa5LhY5zW1lYOvvAinaePkbjlV4nH4ymW\nIUnFKZFIkEgksrLtWIrr/TpwG3B3OP9bwIeAW4FbgHZgPvAscPmYdZPJZJIv/OWjPP/uc+z8i/+V\nYgmSVDpisRikntnnSLWrZztB0E8PC1kFbAOeBtaEr1kDPDnZBppm1XFq0K4eScq1VLt6XgceA14B\nhoBXgf8J1ADrgC8Cu4HPTraBObV1nHYwFknKuVSDH+Dh8DFaJ8HR/wU119XRF/NyTknKtcju3G2u\nr6W/zCN+Scq1yIJ/4ew6BhyMRZJyLrLgXzSnjmSVXT2SlGuRBX/tjGkQG+JEz5moSpCkkhTdQCxl\nMWJ9dew9bHePJOVSZMEPUNFfz74jXVGWIEklJ9LgrxpqYN/RY1GWIEklJ9Lgr042cPCYwS9JuRRp\n8M8sb+RgV2eUJUhSyYk0+GsqGjh80iN+ScqlSIO/bloDnacMfknKpUiDv2F6A129Br8k5VKkwd80\ns4ET/Qa/JOVSpME/t6aBkwMGvyTlUqTBP6++gVNJg1+ScinS4F/Y2MCZmMEvSbmUavBfBrw26nEc\n+CrQCGwAdgDrgfrzbWRRUwP9FQa/JOVSqsH/JnBt+PgAcAp4AlhLEPzLgY3h/KQumdvAYJXBL0m5\nlImunlXATmAvcCfQFi5vA1afb8UFs2ugvJdTvf0ZKEOSdDEyEfx3AT8Kp5uBjnC6I5yf/M3LYsT6\n6nmn3aN+ScqVdAZbB6gCPgXcN8HfkuFjnNbW1rPTZW9PY8/hY1zZMjfNUiSpeCQSCRKJRFa2HUtz\n/U8DvwPcHs5vB+JAOzAfeBa4fMw6yWRy5PNg1u9dz7c/9gh33/6hNEuRpOIVi8Ug/cwG0u/q+Q1G\nunkAngLWhNNrgCcvtIFqGjjgTzNLUs6kE/wzCU7s/tOoZQ8CtxFcznlrOH/+jZQ10HHc4JekXEmn\nj78HaBqzrJPgw+Ci+dPMkpRbkd65C/40syTlWuTB708zS1JuRR78TTMbON7n8IuSlCuRB39zbSMn\nBw1+ScqVyIN/8ewmepJHoy5DkkpG5MHfMreJM+VHoi5DkkpG5MG/dH4T/ZUGvyTlSuTBf+m8BpLT\nuujrH4y6FEkqCZEHf3VVBbEzdbzb0RV1KZJUEiIPfoCK/iZ2HrS7R5JyIS+Cv3qwid2HDH5JyoW8\nCP4ZsdnsPWrwS1Iu5EXw11Y0caDL4JekXMiL4K+vaqKj2+CXpFzIi+BvmtHEkVMGvyTlQjrBXw88\nDvwS2AZcDzQCGwgGYlkfvuaCmmua6Dpj8EtSLqQT/I8A/wpcAVxNMN7uWoLgXw5sDOcvaEFDE92D\n/l6PJOVCqsFfB9wE/CCcHwCOA3cCbeGyNmD1xWxsyewmepIe8UtSLqQa/JcCh4FHgVeBvyIYg7cZ\n6Ahf0xHOX5A/1CZJuZNq8FcA7we+Gz73ML5bJxk+LuhX5s1moMrgl6RcSHWw9X3h4+Vw/nHgfqAd\nmBc+zwcOTbRya2vr2el4PM6Hb7iJZNUJevsGqK5KZ/x3SSoOiUSCRCKRlW3H0lj3eeBugit4WoEZ\n4fKjwEME3wDqmeCbQDI5/otA2do5bP7trVzZMjeNkiSpOMViMUgvs89K5/D6HuCHQBWwC/g8UA6s\nA74I7AY+e7Ebq+xrYlf7EYNfkrIsneB/HVg5wfJVqWyseqiJdzoOp1GOJOli5MWduwA1ZfPYfbjj\nwi+UJKUlb4K/saqZPZ3tUZchSUUvb4J/zoxm2k96xC9J2ZY3wT+/tpkjpw1+Scq2vAn+JY3NdA0Y\n/JKUbXkT/Eub59GDwS9J2ZY3wX/ZwmZ6Kzy5K0nZljfB/94lzQxWH2Jo6KJ+3keSlKK8Cf7G2ukw\nWMXew8ejLkWSilreBD9A5Zlmtu2xn1+Ssimvgn/64Dx2HLCfX5KyKa+Cv6as2Z9tkKQsy6vgb6hs\nZu8xg1+Ssimvgn/OjGbauw1+ScqmvAp+f7ZBkrIvr4L/ktnzODbgyV1JyqZ0BmLZDZwABoF+4Dqg\nEfjfwCWMjMDVdbEbXNo8j5McTKMkSdKFpHPEnwTiwLUEoQ/B+LobgOXARsaPt3teV7cspLdyfxol\nSZIuJN2unrED/94JtIXTbcDqqWxsxaXzGKo+wqne/jTLkiRNJt0j/meAV4Avhcua4exPbHaE8xet\nuqqCst45bH7Hfn5JypZ0+vhvAA4Ccwi6d7aP+XsyfExJdf9CNr+7n+uvWJxGaZKkyaQT/MNnYQ8D\nTxD083cA84B2YD5waKIVW1tbz07H43Hi8fjZ+ToW8eYB+/kllbZEIkEikcjKtsf20V+sGUA50A3M\nBNYDfwysAo4CDxGc2K1n/AneZDI5+ReBq9f+Lssal/OPX/9qiqVJUvGJxWKQemafI9Uj/maCo/zh\nbfyQIPxfAdYBX2Tkcs4pmT9rIftOeMQvSdmSavC/A1wzwfJOgqP+lF3SsJAdb29NZxOSpPPIqzt3\nAZbNW0jXoEf8kpQteRf8Vy5eSE+5wS9J2ZJ3wX/N0oX0V+937F1JypK8C/4Fs2sgWe7Yu5KUJXkX\n/ABVvQvZtMvuHknKhrwM/llDi9i8Z2/UZUhSUcrL4J9T2cIvD74bdRmSVJTyMvgX17bwdufuqMuQ\npKKUl8H/nqZLOHjKI35Jyoa8DP4Vi1s4Org76jIkqSjlZfCvXNbCqardUZchSUUpL4P/fb8yn6Fp\nR+k62Rt1KZJUdPIy+Ksqy6k4vYiXtu+JuhRJKjp5GfwANQMtvPaOJ3glKdPyNvibKlvYtn931GVI\nUtHJ2+BfNOsSdnktvyRlXLrBXw68BjwdzjcSDLy+g2BErvpUN/yephb29+xOszxJ0ljpBv+9wDZg\n+DeU1xIE/3JgI+PH271oVy+5lKNDb6dZniRprHSCfxHwceD7jAwAfCfQFk63AatT3fjNVy6jZ9pb\naZQnSZpIOsH/LeCPgKFRy5qBjnC6I5xPyVUtzSTL+th1oDP1CiVJ46Q62PongUME/fvxSV6TZKQL\n6Bytra1np+PxOPH4+E2UlcWYfnoZic1vsXTB9SmWKUmFKZFIkEgksrLt2IVfMqE/BX4LGACqgVrg\nn4CVBB8E7cB84Fng8jHrJpPJixtWccnv38Unln2S7/3Ob6ZYpiQVh1gsBqln9jlS7er5BrAYuBS4\nC/gpwQfBU8Ca8DVrgCfTKa6lZjnbOuznl6RMytR1/MOH8A8CtxFcznlrOJ+y9zYvY3f3jjRLkySN\nlongf47gah6ATmAVweWcHwW60tnwyqXLOTLkEb8kZVLe3rkLEF+xjFPT32Jo6OLOCUiSLiyvg3/p\ngkZiQxVse/dQ1KVIUtHI6+AHmNV7Oc+8/suoy5CkopH3wb+w6ip+vnNr1GVIUtHI++C/as4Kthza\nHHUZklQ08j74b1h2Ffv6tkRdhiQVjbwP/o+9/0q6p2/xyh5JypC8D/4rlsyhbKiaf39rf9SlSFJR\nyPvgB6g/s4J/e83uHknKhIII/kumX8WL7xj8kpQJBRH81yy4mq1HNkVdhiQVhYII/k9c+0EO8ErU\nZUhSUSiM4L/uCvqq97Hn0PGoS5GkglcQwV9dVUFNzzU8/sKrUZciSQWvIIIfYOn0D/LMtpejLkOS\nCl7BBP+HlnyQzUft55ekdKUa/NXAS8AmYBvwZ+HyRmADwQhc64H6dAscducHV9Je5hG/JKUr1eDv\nBW4BrgGuDqdvBNYSBP9yYGM4nxG3vX8ZgxUneGWHd/BKUjrS6eo5FT5XAeXAMYIhGNvC5W3A6jS2\nf46K8jKaz9xE27PPZ2qTklSS0gn+MoKung7gWWAr0BzOEz43p1XdGNc3f4Sfvm3wS1I6KtJYd4ig\nq6cO+DeC7p7RkuFjnNbW1rPT8XiceDx+UW/4H677CF/48V9NvVJJKjCJRIJEIpGVbccytJ3/ApwG\n7gbiQDswn+CbwOVjXptMJlP7ieW+/kGm/dfZbPnyDq5smZt6tZJUYGKxGGQos1Pt6mli5Iqd6cBt\nwGvAU8CacPka4Mm0qhujqrKcub038v1nEpncrCSVlFSDfz7wU4I+/peApwmu4nmQ4ENgB3BrOJ9R\ntyz6OE9t/5dMb1aSSkamunqmIuWuHoAXtr7LTY+tpPdPDlJVWZ7BsiQpf+VDV09kbrjyEqr6m3ls\nozdzSVIqCi74Aa6Z8Qke/fnTUZchSQWpIIP/yzd/hpd71jkAuySloCCDf82qlQC0PWN3jyRNVUEG\nf1lZjA/XfI5Hnv1h1KVIUsEpyOAH+E+f+hxvDP09J3rORF2KJBWUgg3+2z6wjIYz7+PrbeuiLkWS\nCkrBBj/APSvv5bG3vu1JXkmagoIO/v981x0MlJ3km0/8NOpSJKlgFHTwV5SX8eXlrfzxz77hUb8k\nXaSCDn6Ab9396wxyhj989PGoS5GkglDwwV9RXsbDt3yHR978Gu8cPBZ1OZKU9wruR9omc/Xa36W7\nv4tdf/43lJVF0SxJyp6S/pG2yTxz30O0D23hMw9/O+pSJCmvFU3wz22YyTN3/zNPd/4FX/jLR6Mu\nR5LyVtF09Qz7yctv8ql1t3FL7Zf5yTfup6K8aD7bJJWwfOjqWUwwnu5WYAvw1XB5I7CBYASu9YwM\nz5gzd6y8jBe/9P94+eh65vzBKp54YUuuS5CkvJbqp8e88LEJmAX8O7Aa+DxwBHgYuA9oANaOWTer\nR/zDevsG+Ny3v8sTR/+Ehf1xvnbjb3PPpz7iqF2SClImj/gz1dXzJPA/wsdHgA6CD4YEcPmY1+Yk\n+IftOXScP3zsb/nx/u/TO20vi/pWce3c67hx2fuIX3UZKy6dR3VVRc7qkaRU5FvwtwDPAVcBewiO\n8oe33TlqflhOg3+0l9/cx/fWb+CV/a/ybu8mTlbtYqj6MGW9c5nWP5fKZA3TYjVML6ulumwmFWWV\nVMQqqSirCKbLKqgMn8OdQCz8Jxyen2jZuPkxywHKYl6CKmlyT93/B5AnwT+LIPT/G8FR/zHODfpO\ngn7/0ZIPPPDA2Zl4PE48Hk+zjNSd6u3n9bcPsuvgYY50d3O0u5vOnm66e3voG+ynf3CAvsF+BoYG\n6B9+HuoHYPgDLMnIB9nYZePmh1876sNv9PqSBND97l669+w7O9/+sxchD4K/Evgx8BNg+OL57UAc\naAfmE5wAjrSrR5KKQT5c1RMD/hrYxkjoAzwFrAmn1xB8C5Ak5ZFUPz1uBJ4H3oCz/RT3A78A1gFL\ngN3AZ4GuMet6xC9JU5RvJ3enyuCXpCnKh64eSVKBMvglqcQY/JJUYgx+SSoxBr8klRiDX5JKjMEv\nSSXG4JekEmPwS1KJMfglqcQY/JJUYgx+SSoxBr8klRiDX5JKjMEvSSUmneD/AdABbB61rBHYAOwA\n1gP1aWxfkpQF6QT/o8DtY5atJQj+5cDGcL6kJBKJqEvIKttX2Iq5fcXctkxLJ/j/L3BszLI7gbZw\nug1Yncb2C1Kx/8dn+wpbMbevmNuWaZnu428m6P4hfG7O8PYlSWnK5sndJCMDsUuS8kS6A/e2AE8D\nK8L57UAcaAfmA88Cl49ZZyewNM33laRSswt4TyY2VJGJjYzyFLAGeCh8fnKC12SkcElS7v0IOAD0\nAXuBzxNczvkMXs4pSZIklabbCc4DvAXcF3EtqdoNvAG8BvwiXHa+G9fuJ2jvduCjOavy4k31RrzJ\n2vOBcBtvAY9ksd6pmqh9rcA+gn34GnDHqL8VWvsWE5xL2wpsAb4aLi+GfThZ21opjv1XDbwEbAK2\nAX8WLi+GfXdWOcGJ3RagkqCxV0RZUIreIdgxoz0MfD2cvg94MJx+L0E7KwnavZP8+5mMm4BrOTcY\np9Ke4QsEfgFcF07/K+Nv7ovKRO17APj9CV5biO2bB1wTTs8C3iT4/6oY9uFkbSum/TcjfK4AXgRu\nJAf7LpchdB1BobuBfuDvgU/n8P0zaezVUJPduPZpgnMh/QTt3snIzskXU7kRb6L2XE9wBVcNI9+A\nHiN/bt6bqH0w8RVthdi+doIwADgJ/BJYSHHsw8naBsWz/06Fz1UEB8fHyMG+y2XwLyQ4CTxsHyM7\nsZAkCU5gvwJ8KVw22Y1rCwjaOaxQ2jzV9oxdvp/8b+c9wOvAXzPyVbrQ29dC8O3mJYpvH7YQtO3F\ncL5Y9l8ZwYdbByPdWlnfd7kM/mK5mesGgv8A7wD+I0FXwmgXunGt0P4divFGvO8BlxJ0IxwEvhlt\nORkxC/hH4F6ge8zfCn0fzgIeJ2jbSYpr/w0RtGMRcDNwy5i/Z2Xf5TL49xOcrBm2mHM/pQrFwfD5\nMPAEQddNB0F/JARfuw6F02PbvChclu+m0p594fJFY5bnczsPMfI/1PcZ6X4r1PZVEoT+3zBy70yx\n7MPhtv0tI20rtv0HcBz4F4KTtMWy74Dg5MUugq9sVRTmyd0ZBH1pADOBFwjOrD/MyFVKaxl/MqaK\n4AhlF+nfLZ0NLYw/uTvV9rxE0N8YI79OnsH49s0fNf17wN+F04XYvhhBn+63xiwvhn04WduKZf81\nMdJNNR14HvhVimPfneMOgjPzOwkuSyo0lxL8w28iuLxsuA3nu3HtGwTt3Q58LGeVXryp3og3WXuG\nLyfbCXwn61VfvLHt+wJBmLxB0Ef8JOf+mGChte9Ggu6CTYxc3ng7xbEPJ2rbHRTP/lsBvErQvjeA\nPwqXF8O+kyRJkiRJkiRJkiRJkiRJkiRJkqTS9P8B9C16Bz5INVwAAAAASUVORK5CYII=\n",
"text": [
"<matplotlib.figure.Figure at 0x7efede2e79d0>"
]
}
],
"prompt_number": 17
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"sim.reset()"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 18
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"sim.increment(B, 3)\n",
"sim.increment(A, 5)"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 19
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"sim.run(10000)"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 20
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"sim.plot()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAZAAAAEPCAYAAABsj5JaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VPXZ//H3QNgXISJhrciiAooVFQQXIiriBrjWVvjx\nVK1UrUvbSwv1UaJdRNyeKlp3ARVcqiJSZVGIRSsiIPuuICaQoOw7JJnfH/cZZ7IyTObMOTPzeV3X\nXOfMd07m3JBk7nx3EBERERERERERERERERERERERERERkRTyMlAILIkoewRYASwC3gWOinhtBLAG\nWAn0iyg/zXmPNcA/XIxXRER84hzgVEonkAuBGs75KOcB0AVYCNQC2gFrgYDz2lygh3P+IdDftYhF\nRCRqNQ5/ScxmA9vKlM0ASpzzL4E2zvlAYCJwCFiPJZCeQEugEZZEAMYDg1yLWEREouZmAjmcG7Aa\nBUArIC/itTygdQXl+U65iIh4zKsEci9wEJjg0f1FRKSaMjy45/8AlwDnR5TlA20jnrfBah75hJu5\nQuX5Fb1phw4dgt98801cAxURSQPfAB29DqIi7Sjdid4fWAY0K3NdqBO9NnAc9g8KdaJ/ifWHBKi6\nEz2YDEaOHOl1CFFJhjiTIcZgUHHGm+KMLyAY6we8mzWQiUAfLFl8D4zEhurWxjrTAb4AbgWWA285\nxyKnLPSPuhUYC9TDEshUF2MWEZEouZlAfllB2ctVXP9351HWfODkuEQkIiJx4+UorLSUnZ3tdQhR\nSYY4kyFGUJzxpjj9I3D4S5KG05wnIiLRCgQCEGMu8GIUloiI72VmZrJtW9m50MmradOmbN26Na7v\nqRqIiEgFAoEAqfSZUtm/pzo1EPWBiIhITJRAREQkJkogIiISEyUQERGJiRKIiEgSy87OJjMzk4MH\nDyb83kogIiJJav369cydO5fmzZszefLkhN9fCUREJEmNHz+eCy64gCFDhjBu3LiE31/zQEREKpAM\n80A6duzIAw88QI8ePejatSt5eXk0b968wms1D0RExEcCgfg8YvHZZ5+Rn5/PgAED6NSpE126dGHC\nhMTu0acEIiISo2AwPo9YjBs3jn79+tGoUSMArrnmmoQ3Y6kJS0SkAn5uwtq3bx8tWrSgpKSEhg0b\nAnDgwAG2b9/OwoUL6datW7mvUROWiIgwadIkMjIyWLFiBYsWLWLRokWsWLGCc845h/HjxycsDtVA\nREQq4OcayMUXX8xJJ53EI488Uqr87bff5s477yQvL48aNUrXD9yogSiBiIhUwM8JJBZqwhIREd9Q\nAhERkZgogYiISEyUQEREJCZKICIiEhMlEBERiYkSiIiIxEQJREREYqIEIiIiMVECERFJMu3ataN+\n/fo0atSIzMxMLrvsMvLy8hIeh5sJ5GWgEFgSUZYJzABWA9OBJhGvjQDWACuBfhHlpznvsQb4h4vx\niogkhUAgwJQpU9i1axebNm0iKyuL22+/PeFxuJlAXgH6lykbjiWQ44FPnOcAXYBfOMf+wDOE12b5\nJ3Aj0Ml5lH1PEZG0VadOHa666iqWL1+e8Hu7mUBmA9vKlA0AQjuejAMGOecDgYnAIWA9sBboCbQE\nGgFznevGR3yNiEjaCi2MuHfvXt5880169eqV8BgyEny/LKxZC+eY5Zy3AuZEXJcHtMYSSmTDXr5T\nLiLiucAD8VnQPDjyyFb9DQaDDBo0iIyMDPbs2UPz5s2ZOnVqXGI5EolOIJGCzkNEJCkd6Qd/vAQC\nAd5//3369u1LMBhk0qRJ9OnTh+XLl5OVlXX4N4iTRCeQQqAFUIA1T212yvOBthHXtcFqHvnOeWR5\nfmVvnpOT89N5dnY22dnZcQhZRNLBrl0QCMDSpfDMM15HE71AIMAVV1zBsGHD+Pzzz7nyyiurvD43\nN5fc3Ny43DvRCWQyMBR42DlOiiifADyONVF1wvo9gsBOrD9kLjAEeLKyN49MICKS/HbuhBUroH17\nOOYYd+5RVATTpsFll7nz/m4J9YEEg0EmT57Mtm3b6Ny582G/ruwf1w888EDMMbiZQCYCfYBmwPfA\n/cAo4C1sVNV64Frn2uVO+XKgCLiVcPPWrcBYoB7wIZD4hj4RSbhVq+DEE8PP3doccN685EseAJdf\nfjk1a9YkEAjQrl07xo8fH1UCiSdtaSsivhQo8+nkxq/3H/4ATzwRfr56NXTqFLq/trQ9HM1EFxFf\nKS62D3aAoUPtQx1g+/b43ufgwXDyyMqCiRPDyUOiowQiIr5QUmId2BkZ4Q/2p5+2D/Xu3WH48Kq/\nviJXXAH5lQy7GT3ajnfdBR9+CNddF1vc6UxNWCLiC/fcA488En6+ahUcf7ydh5qzNmyAtm1Lf11x\nMdSsCbt3w9690Ly5NXfViPjzuF8/mD4d+vaFmTNh2TLo2hWuvhrefrvieNSEdXiqgYiIJ4JBGDYM\n5syxD/fI5FFYGE4eAGvW2PFnP7PjgQPw7beWWEI1lkaNrCnq66/LD8OdPt2OM2fasWtXO/btG/9/\nVzpRAhGRI7ZvX7hPIhi0D/yQMWPsg72kpOr3uOceeP556NULLrrIyqZMsfdr3rz0tR07hs8/+ADq\n1oUOHcJlK1aEz7t3h9/9zs6vvx7q1Kn4/qedBrfcUnWMUjUlEBGJ2nPPWZ9E/frQtKnVAmrUgBYt\nYONGm1MRWhT2wQeto7oytWuXfh4IwKWXVn79Zmfa8YABpcuvvhpeeKH89U8/Da+9Bvv3W1IKBq25\nas8e2LLFhu9K9agPREQOq2lTuOaaij+oK9KqlSWU+++HyHlq7dvDP/8JN90E9epZTeHaa+Hoo6FW\nrcO/b+TQ3meegT59YOtWOOccK1u1yvpJ9uyBCy+0RBerzMxMtm0rux5s8mratClbt24tV16dPhAl\nEJEUk5dnHc3791fefBONNWtK90NE6t0bOneGl16y5yeeCCtXhl9//nm4+WY7LymxD/78fGjTpvT7\nTJtmHdzROnjQ+jxq+Lzt5KOP4LHHYOzY8v9mv1ECMUogktYKCmDRIugfsWNOcXH0H7bFxdbZfOaZ\nMGFCuB+hrNCv2c6d8Oab8Jvf2PNQ7SD0+mefhWsGZ5wBX31V/r1++AGaNYsuvmTyzTfhfpslS+Ck\nk7yNpypKIEYJRNJa2ZnbACecAE2awMMPw5AhMHeudVDXrGmvz5oF2dnlh71GWr3a+gzmzLE5E9EK\nBm1U1A8/hMsaNIBu3eCLL9xbmsQPiopKN8nNmAEXXOBdPFVRAjFKIJK2IhPA4MHwl7/AccdVfG3X\nrjYPIiQ0j6Iib7wBv/hFfOKaPBlOOSU8FDcdTJsWrhGuXGkJ3W+UQIwSiKStyNpH6Nfg6aetGapd\nO1i/Prr32boVMjPtfPPm+KyA27evdcCn65DZq66Cd9+1cz9+RCmBGCUQSVtl+x9C9u+3ORMHD8LC\nhdCzp5UvW2ad4JHNVsuWQZcuiYk33YwaBSNGwI4d0Lix19GUpgRilEAkrRw6ZI/zzrO+je3b4aij\nKr8+GLSmrW7dYNCgcFkoiejXxz3FxTZ67JNP/Df7XQnEKIFIygsGbX7Dv/9tNYuyr8UiP9/6QFq0\nqH58Urn69W0Gv98+prQWlkiKCwatP2PWLPsrtmzyCC15HovWrZU8EmHMGDsebomXZHKkWScT25d8\nsQuxVJdqIJKyKhqiG0k/+smhfn1bsr59e68jCXO7BvIp0BhLHvOBF4EnqvwKEYmbssnjiivs2LWr\njW46//zExySx2bev9CKQyS6aPdGPAnYCNwHjgZHAEjeDEhETWbMoLAyvUqsaR3Jq2ND2LQntR5Ls\noqmB1ARaAtcC/3bK9OMr4qKiIlsiJDRCqqio9BLngUDphySHXbvseNJJVa9UnCyiSSAPAtOAb4C5\nQAdgjZtBiaSzF16wZTBefDFcVtlMcUk+06bZsToLXfrF4RJITaAt0A0IzSP9BrjKzaBE0lW/fuFV\nbENmz/YmFnFHv37hPVNatvQ2luqKpvL7FXCG24HEgUZhSdILNUedfDI89FDVGyxJcgt9r+++G0aP\n9jIOdycSPgHUAt4E9kSUL4jlhi5SApGkFlruIrR/hqS+0Pd51y7rEwmtQ5bYGNxNILlU3Gl+Xiw3\ndJESiCSV4mLrHK9TB556Cu64w8r1Y5w+9u61Je5DvPjeaykTowQiSSX012dJSXi01ejR1qQh6ePq\nq+Gdd+zci8UWE5FALgO6AJELKDwYyw1dpAQivrVxo20126NHuCyUQM4/35Ynidfy6ZJcgkH7IyLD\nmZX35pu2T3yiuJ1AngPqAX2BF4BrgC+BG2O5oYuUQMR3Skpg3rzwMuohf/sb3Htv6TL9+Ka3sv1e\nifp5cHspk97A/wO2Ag8AZwI+3FdLxH/atCmfPCCcPAoK7HjTTYmLSfypoMAGUYTs3u1dLNGKJoHs\nc457gdZAEVDdtTtHAMuwJVEmAHWwtbZmAKuB6UCTMtevAVYC/ap5b5G4y8uznedGjrTnq1fbchWb\nNoWv2brV9uIYPNiev/GG7Rk+ezY880ziYxZ/ycqCv//dBlds2mTLnvhdNNWW+4GnsCasp52yF4D7\nYrxnO2Am0Bk4gA0P/hDoCvwIjAb+BDQFhmN9LxOwuSitgY+B44GyiyKrCUs8UVRkM8dDfvwRmjUL\nPy8osFVYGzVKfGwih+N2E9aDwDbgHezD/0RiTx5gCzMeAupjiznWBzYCA4BxzjXjAGfPNAYCE52v\nWQ+sBSK6IkW89dvfln4emTzA/rJU8pBUVNVqvFcRnv8RoPxckHdjvOdW4DFgA9Y8Ng1rusoCCp1r\nCp3nAK2AORFfn4fVREQ8tXevHV96yY733w8PRoxNnDcPTjwx8XGJJEpVCeRyql51N9YE0gG4C6vN\n7ADeBgaXuSZ4mHtX+FpOTs5P59nZ2WRnZ8cYokjVduyAJhG9dLfcAg88AN27w5491h+SCovlSerJ\nzc0lNzc3Lu/lxUTCXwAXYvuLAAzBRnb1xWa3F2DLx8/CmsuGO9eNco5TsT1JvizzvuoDkYTp29e2\nlw3Rj54kK7fngYzE/uIv24wV60TCU4DXsU7x/cBYbJn4Y4EtwMNY0mhC6U70HoQ70TtSvhaiBCKu\nW7gQzj7bahkhxcXhmeQiyaY6CSSaHQn3EP6wrofNSl8ey80ci7CdDedhI6kWAM8DjYC3sAmK67EN\nrHDu9ZZzLAJuRRtaiQfGjAkvww02TLdDByUPSV+xZJ062DyNPnGOpbpUA5G4mzbN+jpyc2H48HD5\nhg3Qtq1nYYnETaIXU8zEmpw6xnJDFymBSNyVXV7i0UfhzjvD6xaJJDu3m7CWRJzXAJrjv4UUReIq\nGITevcuXDx6s5CESEs2vwuUR50XYHI1D7oQj4g/vvw9znNlHzzwDt96qkVYiZUVTbTkT68De6Txv\njC1DUnYYrdfUhCVxMX8+nH66nR84ALVrw759UK+et3GJuMHtPpCFQHfCa0/VxEZQnRrLDV2kBCLV\nkpsL55XZZ1M/UpLq3F4LC0ovXFiMJRGRlBKZPJ57zuZ3iEjlokkg64A7gFpAbeBO4Fs3gxJJtNBo\nq9COgDffrPkdIocTza/Ib4GzgHxsIcMzgZvdDEokkebNC58XFqrZSiRaXqyF5Rb1gcgRKS6GsWPD\nuwHu3Kll1yX9uN0HcgLwCbaDIEA34H9juZmIn8yaFU4e556r5CFypKLJOv8B7gaexUZeBYCl2A6C\nfqIaiEStpARqRgwFOXRIEwQlPbldA6lP6TkfQTSRUJLcpZfacc0a6/NQ8hA5ctEkkB8ove7V1cAm\nd8IRcd9338HUqdCwIXT024puIkkkmmpLB2y59d7Y3ujrgOuxJdf9RE1YclitW8PGjXa+ZQtkZnob\nj4jXErUabwOsxrIrlhslgBKIHFbk6rr6cRFxbzXeP0acR/6qhXYmfDyWG4p45XHnJ3bOHOjZ09tY\nRFJBVQmkERXv/Fd2a1sR3zvmGPjxRztX8hCJj6oSSE6ighBx0/Ll4eQxdaq3sYikkmhGYbUF3sNG\nY/0AvAO0cTMokXh68kk7btgAF13kbSwiqSSajpOPgdeB15zn1zuPC90KKkbqRJcKBQLQoAHs3u11\nJCL+4/ZEwmOAV7DJg4eAsdi2tiK+F/qbQk1XIvEXTQLZAgzB9gDJAAYDP7oZlEi8jB1rx7PO8jQM\nkZQUTbWlHfAUtow7wH+B24ENLsUUKzVhSTmheR/60RCpmFvzQELWA5fH8uYiXtq7146PPeZtHCKp\nqqqs8xQ236Oia4LYLoV+ohqIlPLZZ3DhhbBvn9eRiPiXWzWQ32LLtr8FbAzdyznqk1p87+mnYf9+\nr6MQSV1VJZCWwDXAtUAx8CbwNrA9AXGJVEtJCcycCZMnex2JSOqKttrSBrgO+APwJ+BV1yKKnZqw\nBIBVq+DEE+1c29SKVM3teSCnAXdiw3c/AubHcqMymgD/AlYAy4GeQCYwA1gNTHeuCRkBrAFWAv3i\ncH9JYSedFD5X8hBxT1VZ5y/AJdiH/BvANOK3E+E44FPgZawZrQFwLza/ZDRWy2kKDAe6ABOAM4DW\n2Mz444GSMu+pGogQDEIN58+iE06AlSu9jUfE79zaD6QE2zxqbwWvBYFusdwQOAr4Gmhfpnwl0Aco\nBFoAucCJWO2jBHjYuW4qttDjnLIxKYHIkiXQrRv88AM0a+Z1NCL+59YorLIf8PFyHLYo4yvAKViT\n2F1AFpY8cI5ZznkrSieLPKwmIlJKaNLgJZcoeYgkQlUJZL2L9+wO/A74Cvg/rKkqUpCqhwpX+FpO\nTs5P59nZ2WRnZ1cjTEkm778fPl+2zLs4RPwuNzeX3NzcuLxXTNWWamoBfIHVRADOxpqp2gPnAQXY\nEOJZWBNWKLmMco5TgZHAl2XeV01YaSxU+7j3XrjvPqhTx9t4RJKF26Ow4q0A+B7rCAe4AFgGfAAM\ndcqGApOc88nYEOLaWNLpBMxNVLDif5s22XHePPjrX5U8RBKlqiasT4DzsVFR98T5vrdje4zUBr4B\nfo2t9vsWcCPWfHatc+1yp3w5UATcimbCS4T//V87nnaat3GIpJuqqi3LgZuwoba/ovxe6AtcjCsW\nasJKU4EADBsGzz7rdSQiycetYbzXYLWBs4B5Fbx+Xiw3dJESSJoKBGDxYjj5ZK8jEUk+biWQkPuB\nB2N58wRTAklD06fbPuclJeGOdBGJntsJBGAgcC7WhPUp1uHtN0ogaWjQIBvCq2+9SGzc3lBqFLaM\nyOvOTe4AemNDb0U8078/TJsGV17pdSQi6SmarLME+Dm2pDvYaKmFgN9anFUDSSM33ggvv2zny5dD\n587exiOSrNyeBxKk9Mq4TdAwWvFYaO5Hq1bhpdtFJLGiacJ6CBuyOwvLUn0ov/SISELt2gUzZsAF\nF3gdiUj6irba0grrBwli61dtci2i2KkJK03k50ObNrBjBzRu7HU0IsktEaOwkoESSJoIDdfVt1uk\n+pJtLSyRmN10kx179/Y2DhFRDUSSyCefhPs8du+GBg28jUckFbjZhJWBrZR7QixvnmBKICku1HRV\nWAjNm3sbi0iqcLMJqwjbavbYWN5cJF7y8+1YUKDkIeIX0WSd2cCp2B4ce5yyIDDAraBipBpIClPH\nuYg73F7K5L4KyvRrLAkTqn00beptHCJSWrRZpx3QEfgYqI8lnp0uxRQr1UBSVKj2UVAAWVnexiKS\natwexnsz8DbwnPO8DfBeLDcTidVFFyl5iPhNNAnkNuBswjWO1YC6MSUh3nP+VJkyxds4RKS8aBLI\nAecRkoH6QCQBPv88vFR7RjS9dSKSUNEkkE+Be7G+jwux5iw/biglKebss+34yCPexiEiFYum46Qm\ntjd6P+f5NOBF/FcLUSd6CtmwAY49FsaMgdtu8zoakdSViMUU6wAnYkljJXAwlpu5TAkkhZx+Osyf\nr3kfIm5zex7IpcCzwLfO8/bAMODDWG4ocjjz59vj8ce9jkREqhJN1lmFJZG1zvMOWPLw2/pYqoGk\niB494KuvYM8eqF/f62hEUpvb80B2Ek4eYDURv00ilBQRDFryyM1V8hDxu6qasK5yjvOwGsdbzvNr\nnDKRuJs82Y6hEVgi4l9VVVvGEh5pFajg/NfuhRUTNWGlgEAAzj0XPv3U60hE0oNbnej/E8sbisSq\nuNiOAwd6G4eIRCearNMeuB1bUDGUcOKxnHtNrCksD7gcyATexPYeWQ9cC2x3rh0B3AAUA3cA0yt4\nP9VAktySJdCrF+zcCTW02bJIQrg9jHcSNnHwA6DEKYvHJ/WdwHKgkfN8ODADGA38yXk+HOgC/MI5\ntsZWBD4+IhZJEc89B506KXmIJItoflX3A08CM4Fc51HdFuo2wCVYYgplvgHAOOd8HDDIOR8ITAQO\nYTWTtUCPat5ffOiNN+AEvw0OF5FKRVMDeQrIwZYwiVxUcUE17vsEcDfQOKIsCyh0zgud5wCtgDkR\n1+VhNRFJIcXFsGUL3FfR9mUi4kvRJJCuwBDgPEo3G50X4z0vAzYDXwPZlVwTpOpmsgpfy8nJ+ek8\nOzub7OzK3l78ZsUKO3bu7G0cIqkuNzeX3NzcuLxXNB0n3wCdid/6V3/HElIRUBerhbwLnIEllAKg\nJTALW39ruPN1o5zjVGAk8GWZ91UnehLr1w8WLIAff/Q6EpH04vZM9CVAPHej/jPQFjgOuA7rWxkC\nTAaGOtcMxTrvccqvA2o7X9MJmBvHeMRjRUUwY4bWvhJJNtE0YTXFVuD9inAfSDyG8YaEqg2jsNnu\nNxIexgs2Uust51gE3Ir/lpKXanjtNTsOHuxtHCJyZKKptmRXUp4bvzDiQk1YSaplS/jzn+H2272O\nRCT9uD0PJDeWNxapTHExDB0KL78MJSVQUABdu3odlYgcqWiyzm7CTUa1gVpOWeNKv8IbqoEkgY8+\ngksuKV+ub52IN9yugTSMOK+B9X2cGcvNJH0Fg5XPMB83ruJyEfG3I100ogQbHdXfhVgkhf33v+Hz\nYNCasXbuhBdeUOe5SLKKptpyVcR5DeA0oA/Qy5WIYqcmLJ9ZuhROPtn6Om64wcpKSmzJdhHxh+o0\nYUXzRWMJ94EUYUNsX8Bmk/uJEoiP3HADvPJK6bLcXOjTx5NwRKQSbieQZKEE4gM7d8JRR4WfL1oE\nDzwAs2fDZr/9ySEiriWQkZWUhz6lH4zlhi5SAvHIgQNQt275cn07RPzPraVM9mDDdSMfQWym+J9i\nuZmkprJNVV9/reQhkg6izTqNsZ0Ab8SWFXkM9YEI4SaryZPh8su9jkZEjpSbiykeDfwVWIRNIOyO\n1T78ljwkgdatsxpGURFcdpmVVTQ5UERSW1UTCR8FrgCeB7oBuxISkfhWMAjnnw+zZpUuHz8eatb0\nJiYR8U5V1ZYSbA+QQxW8FkRLmaSVkhK44gprqoo0aRIMHOhNTCJSfW4tZXKks9QlRV18MUydaueP\nPQa//70mA4qI5oFIFRYvtr6N/Hx7fsUV8O673sYkIvHl9mKKkmZmzYK+fcPPu3eHefNU6xCR0tRM\nJeUMGxY+f/BBmD9fyUNEylMNRErZutUeS5bASSd5HY2I+JkSiPwkLw9eegl69VLyEJHDUxNWGvvu\nO+jRA955B3r3hrZtIScHrrvO68hEJBmkUsu2RmFFoajIJv3dcw88+mjp1+68Ex55BGrV8iY2EUk8\njcKSqFTUEV5QAJ99BqeeCu3bJz4mEUleqoGkibLJY/hweOghb2IREf9QDUSqtHp1+Fw5VkTiRTWQ\nNBCqfei/R0TKcnM5d0lyh5ylMJcu9TYOEUk9asJKUSUlVuN48UVo3Bi6dvU6IhFJNUogKaiwEFq0\nCD+/5x7vYhGR1OVFE1ZbYBawDFiKbZULkAnMAFYD04EmEV8zAlgDrAT6JSzSJLRwYenkccYZ8PDD\n3sUjIqnLi070Fs5jIdAQmA8MAn4N/AiMxrbNbQoMB7oAE4AzgNbAx8Dx2IZXkdSJTrjDfNs226tc\niyCKSFWSrRO9AEseALuBFVhiGACMc8rHYUkFYCAwEdsZcT2wFuiRoFiTQjBoOwX+7Gf2vKAAmjRR\n8hARd3ndB9IOOBX4EsgCCp3yQuc5QCtgTsTX5GEJJ+198glccEHpsrvvhqysiq8XEYknLxNIQ+Ad\n4E5gV5nXgs6jMmqrAu66q/TzBQtsSRIRkUTwKoHUwpLHq8Akp6wQ6xspAFoCm53yfKzjPaSNU1ZO\nTk7OT+fZ2dlkZ2fHMWR/WbDA5nZs2QKZmV5HIyLJIjc3l9zc3Li8lxet5AGsj2ML8PuI8tFO2cNY\n53kTSnei9yDcid6R8rWQtOlEP3DAtpxdvBh2la27iYgcgWRbC+ssYDCwGPjaKRsBjALeAm7EOsuv\ndV5b7pQvB4qAW0njJqxgEOrWtfOPPvI2FhFJb6k0TictaiC/+hVMnAiffgrnnut1NCKS7KpTA1EC\nSSKbNkGrVvDkk3D77V5HIyKpQAnEpHQCCQahZUtbpqSkRHM8RCQ+km0iocRg6lRo1sw60JU8RMQP\nlECSxCWXQI8eULu215GIiBglkCQwYYId//IXb+MQEYmkBJIE3nnH9jBvrQVcRMRHvF4LS6qwfz+8\n+iq8+y5s3Oh1NCIipSmB+NigQTBtGpxyio3AEhHxEzVh+dTQoZY8AN5/39tYREQqogTiM0uXwi9/\nCePH206CwSAce6zXUYmIlKcmLB8IBmH2bKtpPP64bQa1eDGcfLLXkYmIVE4JxAcmTIDBg+384ovh\nww+9jUdEJBpKIAlUdgmSWbPg/PPtfNQoGDbMah8iIslACSQBgkFo1Aj27CldftRRdnz11XANREQk\nWSiBuOj11+G222DHDnt+//3QrRvMmAF5eTY58KyztLaViCSnVPro8tVqvM8+C7fcAkcfbbWPyZPV\nKS4i/qPVeH1kxgzIzrbkkZMDP/4I69YpeYhI6lECiZMDB+Chh2zV3F69YN8+GDnS66hERNyjPpA4\nWLQIfv5zO581y2ogIiKpTjWQati922aLX3IJnHGGbTmr5CEi6UI1kCMQDMKVV8KkSaXLp0yBSy/1\nJiYREa/0+EEGAAAIN0lEQVQogUTh009L1ywaNIBWrWwG+WmnaRiuiKQnNWFV4Y03bCn1UPL41a9s\nj47du2H1ajj9dCUPEUlfqoFUYO1a6NsXvv8eXnvNmq3q1fM6KhERf1ENJML27baYYY8e0LGjLa1+\n/fVKHiIiFUnrBBIMQnGxzRofNAiaNoXmzWH9epg5E7p29TpCERH/SqsmrB07YMECGDsWvv4aliyx\n8mOPtYUOx42DIUPUryEiEo1U+qiscC2sYBDWrIHp0+GPf4TGjaF/f2ua+vnPISMDmjXzIFoRER+o\nzlpYyVQD6Q/8H1ATeBF4uKqLN22C++6DefNg82bo2RO++AK6d09EqCIiqS9Z+kBqAmOwJNIF+CXQ\nuaILt261meHt2kH9+jBmjI2meu89fySP3Nxcr0OISjLEmQwxguKMN8XpH8mSQHoAa4H1wCHgDWBg\n2YveeMPmZrRoYYnkH/8IctZZQWrUsKatYDDId9u/Y82WNQkMvbRk+aFKhjiTIUZQnPGmOP0jWZqw\nWgPfRzzPA3qWveiXbw6l1ZXb+bL9Who+uvyn8ga1GtC+aXs279nM7oO72XNoDx2adqBB7Qb2Zjvz\nyGqQxdZ9W7mo40V0bNrxp9eO1AlHn0C9WpWP+123bR0z180sVVa7Zm26HNOlVFmdmnVijiEa+4v2\ns/fQ3kpf99PeKiLiT8mSQKL6NBszoifrtn/L8UdfSs/WPemW1Y1AIEDh7kIKdhfQqE4j2jdtz6Hi\nQ6zdupaDxQcBqFWzFoeKD7Fg0wKCBFm9ZTXb9m874iD3F+1nyuopVV6z7rt1fP+f70uVzd4wm4a1\nGxKI6Mfatn8brRu1DnVwxV3ezjzqZtSlXkb5ZLe/aD/B/wZ5+YmXXbl3pGAwSMPaDel0dKcj/tpV\nS1Yxf+J8F6KKL8UZX36PMxgMsnTzUrZ+sZWXnngp5vdZc/sa6mbUjWNk8Zcso7DOBHKwPhCAEUAJ\npTvS1wIdEhuWiEjS+wbo6HUQbsrA/pHtgNrAQirpRBcRESnrYmAVVtMY4XEsIiIiIiKSrvoDK4E1\nwJ88juVloBBYElGWCcwAVgPTgSYRr43A4l4J9EtQjABtgVnAMmApcIdT7rdY6wJfYk2Wy4GHfBon\n2Fylr4EPfBzjemAxFudcp8yPcTYB/gWswL7vPX0Y5wnY/2PosQP7PfJbnKH7LsM+lyYAdXwaZ8LV\nxJq02gG18L5v5BzgVEonkNHAPc75n4BRznkXLN5aWPxrSdy8nBaAs4s7DbGmwc4+jbW+c8wA5gBn\n+zTOPwCvA5Od536McR32wRHJj3GOA25wzjOAo3waZ0gNYBP2h5nf4mwHfIslDYA3gaE+jNMTvYCp\nEc+HOw8vtaN0AlkJZDnnLZznYFk+ssY0FRtt5oVJwAX4O9b6wFdAV/wXZxvgY+A8wjUQv8UIlkCO\nLlPmtziPwj7wyvJbnJH6AbOdc7/FmYn9gdgUS8YfABfGK85kzywVTTBs7VEslcnCmrVwjqFvWiss\n3hCvYm+H1Zq+xJ+x1sD+Iiok3OzmtzifAO7GhpaH+C1GsPlUHwPzgN84ZX6L8zjgB+AVYAHwAtDA\nh3FGug6Y6Jz7Lc6twGPABmAjsB1ruopLnMmeQJJtunSQqmNO9L+nIfAOcCewq4JY/BBrCdbc1gY4\nF/srv2wcXsZ5GbAZawevbF6V1zGGnIX9sXAxcBvW5Fo2Dq/jzAC6A884xz2Ub1XwQ5whtYHLgbcr\nicPrODsAd2F/KLbCfucHVxBHTHEmewLJx9odQ9pSOnv6QSFWRQRoiX3YQPnY2zhliVILSx6vYk1Y\n4N9YwTop/w2chr/i7A0MwJqHJgJ9sf9TP8UYssk5/gC8h60x57c485zHV87zf2GJpAB/xRlyMTAf\n+z8F//1/ng78F9gCFAHvYk3/fv3/TCg/TjBsR/lO9FCb4nDKd1bVxqrt35C4lQECwHis6SWS32Jt\nRnh0SD3gP8D5PowzpA/hPhC/xVgfaOScNwA+x9ru/RYn2Pf5eOc8x4nRj3GCLew6NOK53+I8BRtp\nWc+53zis9um3OD3jpwmGE7F2xoNY38yvsU6sj6l4uNyfsbhXAhclMM6zsaahhYSHIfb3YawnY+3g\nC7Hhp3c75X6LM6QP4VFYfovxOOz/cSH2gRL6XfFbnGAfel8Bi7C/mI/yaZwNgB8JJ2bwZ5z3EB7G\nOw5rffBjnCIiIiIiIiIiIiIiIiIiIiIiIiIiIiIi1XE04Xkym7AZ0l9jy76MieN9HgWyq3j9DmBI\nHO8nIiIJNBJbsj3eGhHem6M614gkXLKvhSWSSKElHbIJL1mSg83u/Q+2YdOVWI1iMfARttwO2Bpe\nudhKuFMJr0M0EJsRHDIKmzW8CHjEKduFrWXUNW7/EpE4UAIRqb7jsFWCBwCvYctldwP2AZdiS0c8\nBVyFLW73CvA352vPxpIKWJPZICxRnAL8NeIec7HViEV8I+Pwl4hIFYJYTaMYW2OqBjDNeW0Jtrjm\n8VhSCNU0amJrpgH8jPAquTuA/cBLwBTnEbIRaO/GP0AkVkogItV30DmWAIciykuw37EA1izVu5Kv\nD7UEFGFLrJ8PXA38zjnHeY9k2/9GUpyasESqJ5qlrlcBxxDeGrQWtmw2wHeE+0MaYKuifoR12J8S\n8R4tsT4WEd9QAhGJXjDiWNE5lK8lBLFaydXAw4SX0O/lvP4Z1i8C0BjrnF+E7bH9+4j36UF4320R\nEREaEt59rzKNo7hGJOFqeh2ASJo7CHTC+j/WV3LNMGwI8OLEhCQiIiIiIiIiIiIiIiIiIiIiIiIi\nIiIiPvX/Abbq4A3PvuRoAAAAAElFTkSuQmCC\n",
"text": [
"<matplotlib.figure.Figure at 0x7efedbb95d90>"
]
}
],
"prompt_number": 21
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"dsolve([5,3], 2500)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEACAYAAACznAEdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAG0xJREFUeJzt3XmUFGWaqPGnWAVRaVSQTUGEBm1tleMObbmhqBeY8R7F\n3hztdu2x9fa4oTNX2u7xiL2paLeizW0U3LdRh1ZAKbVFcEFcQFZlFFRQ9kXWyvvHF0UlRRVUZWVm\nRGY+v3PiRGRkZOabQVS+fN/7RQRIkiRJkiRJkiRJkiRJkiRJjTYaWAJ8mLauHTARmAtMANqmPTcM\nmAfMBgakre8bvcc84M4cxitJypH+wBFsnxBuB66Llq8HbouWDwZmAM2BbsB8oCx67i3g6Gh5PHBG\nziKWJOVMN7ZPCLOBDtHyftFjCK2D69O2exE4FugIfJy2fihwby4ClSRlrkkGr+lA6EYimlclh07A\norTtFgGda1m/OFovSUqQTBJCulQ0SZIKXLMMXrOE0FX0FaE7aGm0fjHQNW27LoSWweJoOX394tre\nuEePHqkFCxZkEJIklbQFwEGNfZNMWgjPARdEyxcAz6atHwq0ALoDPQnF5K+A1cAxhCLzT9Jes50F\nCxaQSqWcUiluvvnm2GNIyuS+cF+4L3Y+AT0y+C3fwa5aCI8AJwL7AJ8D/5cwquhx4GfAQuDcaNtZ\n0fpZwBbgCqq7k64A/ga0IowyejEbwUuSsmdXCeH8OtafWsf6W6OppneBQ+sblCQp/xpbVFaOlJeX\nxx1CYrgvqrkvqrkvsq9s15vkVSrqD5Mk1VNZWRlk4ffcFoIkCTAhSJIiJgRJEmBCkCRFTAiSJMCE\nIEmKmBAkSYAJQZIUyeRqp5KkLKo6Hzf9vNza1uWaCUFSg2zdCmvWwNq18O239Zs2boTNm2HLljCv\nz1S1bWVlmLZurV7e2VTf7VKpHX90d/aDXN95Q7etqaxsx+WyPF1TwktXSCVq/Xr4+uswLV1avfz1\n17BsGaxaBatXh6lqedWq8Lo99oDdd4dWreo3tWwJzZpB8+b1n6q2b9oUmjQJU/ryzqZdbVdWFuaw\n44/uzn6Q6zvP9DWZytalK0wIUhFKpeCLL2DuXPif/4HPPw/TZ59VL2/aBO3bw777Vk9Vj/feG/ba\nK0x77rn9fPfdq39MlQwmBEmkUvDJJ/Dee/DRRzBnTpjmzoU2baBnT+jWDbp2hf33D/OqqW3b/HVF\nKLdMCFIJ+uILeP11mDoVpk+H998P3TdHHgmHHgq9e0OvXmFq2zbuaJUvJgSpBHz1Fbz0ElRUhESw\nYgX06wfHHQd9+8IRR8A++8QdpeJmQpCKUCoF77wDL7wA48fD/Plw6qlw8snwgx9Anz7232tHJgSp\niMycCY88EqamTWHwYDjrLDjhhDDSRtqZbCUEz0OQYrJmDYwbB/fdB998A0OHwhNPhG4gi72KQ9IO\nO1sIKnpz58Jdd8HDD0N5OVx+OZxyil1Bypy30JQKzAcfwHnnhW6g73wnPH76aTjtNJOBksHDUMqx\njz+GIUPg9NPhqKPCeQO/+Q106RJ3ZNL2TAhSjixdCldcEUYH9esXEsE114TzBqQkMiFIWVZZCffc\nAwcfDC1awOzZIRG0ahV3ZNLOOcpIyqKZM+Hii0NN4LXXQlKQCoUtBCkLKivh9tvDqKGf/tRkoMJk\nC0FqpC+/DElgwwZ4991wETmpENlCkBrh5ZfDheVOOAEmTzYZqLDZQpAykErByJFw663hchMnnRR3\nRFLjmRCkBtq0KQwnfestePNN6N497oik7DAhSA2wdi2ccw7sthtMmRJuQiMVC2sIUj0tWxYuRd2l\nCzz1lMlAxceEINXDl1+GM45PPBEeeCDcAF4qNiYEaReWLg1XIz3/fBgxwktTq3gl7dD28tdKlGXL\nwgiiIUPgllvijkaqXRIufz0MmAl8CDwMtATaAROBucAEoG2N7ecBs4EBjfhcKS9Wr4YBA2DgQPj1\nr+OORsq9TDNKN+AVoA+wEXgMGA8cAnwD3A5cD3wHuAE4mJA0jgI6A5OAXkBljfe1haBE2LwZzj4b\nunWDe++1m0jJFncLYTWwGWhNGLraGvgCGASMibYZAwyJlgcDj0SvWQjMB47O8LOlnEql4LLLQuH4\nnntMBiodmSaE5cAfgM8IiWAloauoA7Ak2mZJ9BigE7Ao7fWLCC0FKXH+8z9hxgx47DFHE6m0ZHq4\n9wCuJnQdrQKeAH5cY5tUNNWl1ueGDx++bbm8vJzy8vIMQ5Qa7tlnYdQomDbN8wyUXBUVFVRUVGT9\nfTNtDJ8HnAb8PHr8E+BY4GTgJOAroCMwGehNqCMA3BbNXwRuBqbVeF9rCIrNnDnQvz+88AIcbYem\nCkjcNYTZhATQKgriVGAW8DxwQbTNBcCz0fJzwFCgBdAd6Am8leFnS1m3Zg380z+Fi9WZDFSqGpNR\nriP86FcC0wmthT2Ax4H9CcXjcwn1BYAbgYuALcBVwEu1vKctBOVdKgXnnQd77QX33x93NFLDZauF\nkLTxEyYE5d1f/wp33RXqBrvtFnc0UsOZEKQsmDs33NymogIOOSTuaKTMxF1DkArepk3wwx/C8OEm\nAwlsIaiEDRsGH30Ezz3nyWcqbNlqIXjajUrSO+/A6NHw4YcmA6mKXUYqOZs2wUUXwR//CO3bxx2N\nlBwmBJWcESOga9dQP5BULWmNZWsIyqmZM6G8HKZPD0lBKgaOMpIaKJWCSy8NN7oxGUg7MiGoZIwd\nCxs3wiWXxB2JlEx2GakkrF4NvXvDM8/AMcfEHY2UXZ6pLDXAv/0brFwZLlMhFRsTglRPVYXkmTMd\nZqriZEKQ6mnAgHB/5F/+Mu5IpNxwlJFUDxMmwMKFcPnlcUciJZ8JQUWrshKuuw5uuw2aN487Gin5\nTAgqWuPGQevW4U5oknbNGoKK0oYN8N3vhqTQr1/c0Ui5ZQ1B2om774YjjzQZSA1hC0FFZ/VqOOgg\nePVV6NMn7mik3LOFINVh5Eg4/XSTgdRQthBUVFatCq2DN96AXr3ijkbKD1sIUi3uugsGDjQZSJmw\nhaCisXIl9OwJU6aEuVQqbCFINdx5J5x1lslAypQtBBWFVaugRw+YOjXUEKRSYgtBSvOXv4TagclA\nypwtBBW8DRvgwAPhpZfg0EPjjkbKP1sIUuTBB8NZySYDqXFsIaigbd0abo05ejT07x93NFI8bCFI\nwNNPw777es0iKRtMCCpYqRSMGAE33ABlSWvrSgXIhKCC9fLL8O234faYkhrPhKCC9bvfwbXXQhOP\nYikrktbQtqisepk1C045JdwvuWXLuKOR4mVRWSXtrrvgsstMBlI2NSYhtAWeBD4GZgHHAO2AicBc\nYEK0TZVhwDxgNjCgEZ+rErd8OTz2WEgIkrKnMQnhTmA80Ac4jPBDfwMhIfQCXo4eAxwMnBfNzwD+\n3MjPVgm7/34YNAg6dIg7Eqm4ZNrntBfwHnBgjfWzgROBJcB+QAXQm9A6qARGRNu9CAwHptZ4vTUE\n7dTmzeEids8+G85OlhR/DaE78DXw/4DpwP3A7kAHQjIgmlf9H64TsCjt9YuAzhl+tkrYM89At24m\nAykXMk0IzYAjCV0/RwLrqO4eqpKKprrYFFCD3XknXH113FFIxalZhq9bFE1vR4+fJHQLfUXoKvoK\n6AgsjZ5fDHRNe32XaN0Ohg8fvm25vLyc8vLyDENUsXn7bVi8GAYPjjsSKV4VFRVUVFRk/X0b0+f0\nGvBzwoii4UDraP0yQq3gBsIooxsIxeSHgaMJXUWTgIPYsZVgDUF1+ulP4bDD4Jpr4o5ESpZs1RAa\n8wbfBx4AWgALgAuBpsDjwP7AQuBcYGW0/Y3ARcAW4CrgpVre04SgWi1fHu55sGAB7L133NFIyZKE\nhJALJgTV6o474J13YOzYuCORkifuUUZS3qRSMGoUXHpp3JFIxc2EoMR7442QFLzngZRbJgQl3qhR\ncMkl3vNAyrWk/YlZQ9B2LCZLu2YNQSXhoYfgrLNMBlI+mBCUWBaTpfwyISixpkyBrVuhf/+4I5FK\ngwlBiWUxWcqvpP2pWVQWACtWQPfuMH8+7LNP3NFIyWZRWUXtoYfgzDNNBlI+mRCUOBaTpXiYEJQ4\nb74Z7oz2gx/EHYlUWkwISpz77rOYLMUhaX9yFpVLnMVkqeEsKqsojR1rMVmKiwlBiVFVTL7kkrgj\nkUqTCUGJMXUqbNwIJ54YdyRSaTIhKDEsJkvxStqfnkXlElVVTJ43D/bdN+5opMJiUVlFZdw4GDjQ\nZCDFyYSg2FlMlpLBhKDYTZ0K334L5eVxRyKVNhOCYudlrqVkSNqfoEXlErNyZSgmz51r/UDKlEVl\nFYVx4+D0000GUhKYEBSbVKr63ANJ8TMhKDbTpllMlpLEhKDYjBoFF18MTTwKpUSwqKxYrFoF3brB\nnDnQvn3c0UiFzaKyCtq4cTBggMlAShITgvLOYrKUTCYE5d1bb8G6dXDSSXFHIimdCUF5ZzFZSiaL\nysori8lS9llUVkEaOxZOO81kICWRCUF5U1VMvvTSuCORVJvGJoSmwHvA89HjdsBEYC4wAWibtu0w\nYB4wGxjQyM9VAZo6FTZssJgsJVVjE8JVwCygquP/BkJC6AW8HD0GOBg4L5qfAfw5C5+tAnPvvWGo\nqcVkKZka86fZBTgTeIDqYsYgYEy0PAYYEi0PBh4BNgMLgfnA0Y34bBWYFSvgv/4L/uVf4o5EUl0a\nkxD+BFwLVKat6wAsiZaXRI8BOgGL0rZbBHRuxGerwDz4IJx1FuyzT9yRSKpLpgnhbGApoX5Q11Cn\nFNVdSXU9rxKQSoXuIovJUrI1y/B1xxO6h84EdgP2BB4itAr2A74COhKSBsBioGva67tE63YwfPjw\nbcvl5eWUe23kgvf66+H2mP37xx2JVBwqKiqoqKjI+vtm48S0E4FrgP8F3A4sA0YQCspto/nBwMOE\nukFnYBJwEDu2EjwxrQj96Edw1FFw9dVxRyIVp2ydmJZpC6Gmql/x24DHgZ8RisfnRutnRetnAVuA\nK7DLqCR88w3893/DyJFxRyJpV7x0hXLq97+HDz+EMWN2va2kzGSrhWBCUM5UVkLv3vC3v8Hxx8cd\njVS8vJaREm/iRGjdGo47Lu5IJNWHCUE5M3IkXHllGGEkKfmS9qdql1GRWLAAjj0WPvsMWrWKOxqp\nuNllpES75x648EKTgVRIbCEo69auhQMOgHffDTfDkZRbthCUWGPHhrOSTQZSYTEhKKtSKbj77lBM\nllRYTAjKqoqKcP7BySfHHYmkhjIhKKtGjoR//VeHmkqFKGl/thaVC9iCBXDMMbBwIbRpE3c0Uumw\nqKzEueOOcItMk4FUmGwhKCuWLYOePWHmTOjYMe5opNJiC0GJcu+9MGSIyUAqZLYQ1GgbNkD37uFi\ndt/7XtzRSKXHFoISY9w4OPxwk4FU6GwhqFEqK0MiGDkSTjkl7mik0mQLQYkwfjy0bOmJaFIxMCEo\nY6kU/Pa3cOONnogmFQMTgjI2aRKsXg3nnBN3JJKywYSgjFW1Dpp4FElFwT9lZeS112DxYhg6NO5I\nJGWLCUEZ+e1vYdgwaNYs7kgkZYsJQQ02bRrMmQM/+UnckUjKJhOCGuw//iO0Dlq0iDsSSdlkQlCD\nTJ4Mn3wCP/tZ3JFIyjYTguotlQqjim65BZo3jzsaSdlmQlC9Pf88rFvnyCKpWDlGRPWydSvcdBPc\neqvnHUjFyj9t1cvDD4c7oZ19dtyRSMqVpF2BxqudJtC6ddC7Nzz6KJxwQtzRSKrJq50qb373O+jX\nz2QgFTtbCNqpzz8PN7+ZPh0OOCDuaCTVJlstBBOCdupHP4IDD4Tf/CbuSCTVxYSgnJsyBc49F2bP\nDgVlSclkDUE5tXkzXHop/P73JgOpVGSaELoCk4GZwEfAL6P17YCJwFxgAtA27TXDgHnAbGBAhp+r\nPPnjH6FLFzjvvLgjkZQvmTYx9oumGUAb4F1gCHAh8A1wO3A98B3gBuBg4GHgKKAzMAnoBVTWeF+7\njBLg00/hqKPg7behe/e4o5G0K3F3GX1FSAYAa4GPCT/0g4Ax0foxhCQBMBh4BNgMLATmA0dn+NnK\noVQKfvELuPZak4FUarJRQ+gGHAFMAzoAS6L1S6LHAJ2ARWmvWURIIEqYBx8Md0L71a/ijkRSvjX2\nWkZtgKeAq4A1NZ5LRVNdan1u+PDh25bLy8spLy9vVICqv88+Cy2DiRO9mqmUZBUVFVRUVGT9fRvT\n59QceAH4O3BHtG42UE7oUupIKDz3JtQRAG6L5i8CNxNaFemsIcSkshIGDIBTTgk3v5FUOOKuIZQB\nfwVmUZ0MAJ4DLoiWLwCeTVs/FGgBdAd6Am9l+NnKgXvugfXrQwtBUmnKNKP0A14DPqC662cY4Uf+\ncWB/QvH4XGBl9PyNwEXAFkIX00u1vK8thBjMmAGnnRZOROvZM+5oJDWUZyorK1avhr59w13Qzj8/\n7mgkZcKEoEZLpcKJZ+3awb33xh2NpExlKyF4x7QSdvfdMG8evPlm3JFISgJbCCVq0iT48Y/hjTeg\nR4+4o5HUGLYQlLG5c8NlrR9/3GQgqZpXOy0xy5eH+yLfeiuceGLc0UhKEruMSsi6dWF46QknhNti\nSioOjjJSg2zcCIMGQadOMHo0lCXtX15SxkwIqrctW+CHPww3vXniCWhm5UgqKhaVVS+bN4fRRCtX\nwnPPmQwk1c2fhyK2cWM48WzLlpAMdtst7ogkJZmjjIrU2rUweHBoETz9tMlA0q6ZEIrQ4sXQv3+4\nJ/Kjj0KLFnFHJKkQmBCKzPvvw3HHwdChcP/91gwk1Z8JoYiMHQunngp/+ANcf71DSyU1jP9/LAIb\nNsDVV8Mrr8DLL8Nhh8UdkaRCZAuhwH30ERx/PCxbBu+8YzKQlDkTQoHasgVuuw1OOgkuvzxcqG7P\nPeOOSlIhs8uoAE2fHpJAmzahVXDAAXFHJKkY2EIoIMuWwWWXwZlnwsUXw8SJJgNJ2WNCKADr14eR\nQ336hHMKPv4Yfv5zaOK/nqQssssowTZuhAceCPcuOPZYmDwZDjkk7qgkFSsTQgKtWAGjRsHIkWHU\n0HPPQd++cUclqdjZ6ZAgs2fDlVfCgQfCzJnw/PMwfrzJQFJ+2EKI2erVYcjo6NHw6adw4YXh3ILO\nneOOTFKpSdrFDUriBjnr1sHf/w5PPRXmJ50EF10EZ5wBzZvHHZ2kQuMd0wrMkiUwYQI88wxMmgTH\nHAPnnAP//M/Qvn3c0UkqZCaEhPv2W/jHP8K5AhMmwMKFoSUwaFCY9t477gglFQsTQsJ8/jlMmQJv\nvhnmM2fC4YfDgAFw2mlw9NFeilpSbpgQYrJ1K8yfH+47UDW99164d/Hxx4d7ERx/fBgZ1Lp13NFK\nKgUmhBxbswbmzYO5c6vnc+bArFmhz//7399+6t7d+w9IiocJoVEfAt98A4sWha6ezz+vXv7ss5AA\nVq+Ggw6CXr3C1LNnmH/ve15VVFKymBDqfINwpm/VD336j33V8qJF0KoVdO0a7jvctev2U8+e0KmT\n1wqSVBhKOiEsXx66b+bMCaN3qv5nX/Wj36zZ9j/wNX/0u3SB3XfP/ZeRpHwomYSwfj1UVMDUqWEE\nz/vvh1tGfve7YTrwwB3/h2+XjqRSUtQJobISXnoJ7rsv3Cf4yCOhX79wxc8jjgjdORZwJSko1IRw\nBnAH0BR4ABhR4/nUtGkpLr00PLjyShgyBNq1y2uMklRQCjEhNAXmAKcCi4G3gfOBj9O2SbVvn+JP\nf4KhQ1Os3rSKxasXs/zb5azfvH7bVJmqrP4CaU2F5k2a06Jpi21T86bbP65tatm0JS2btaRJWbIq\nyBUVFZSXl8cdRiK4L6q5L6q5L6plKyHk89zZo4H5wMLo8aPAYLZPCPzvkb/mvvWvcNmI9ygrK6PT\nHp3Yu9Xe7N5id1o3b02rZq1o2qQpAOn1hhQptlRuYdPWTTtMm7du3mHdxq0bw3zLRjZu3UizJs22\nJYf0RFHXfNs2u3q+Pu9Ry3OvTH7Fgz3iH34190U190X25TMhdAY+T3u8CDim5kbN9ljOTX1vom/H\nvuzdOj8X/EmlUmyu3LwtOaQnirrmO9tmxbcrdtymge+x4dUN3HrLrSFp7CSptGzasvZtdvV8fd4j\nmjdv2pwmZU1oWtaUJmVNwnKTptutK7OoIxW8fCaEep1gcOfAO3Mdxw7Kysq2dSHtwR55//za3Lzp\nZv79pn/f1prZWYKqzzZrNq7hm63f1JmgdvYeWyq3sDW1lcpUJVsro3n0uGoqo6zWRFGfx2VRS7es\nrGyH5bKyMpa+vZQn//zktqRTtX5nr6v5Humvq2u5tvfIh4Yk009nfMqrY17N7HMS+p0y9cn7n/CP\nB/+R88+pUgr/6cnnNzwWGE4oLAMMAyrZvrA8H+iRx5gkqRgsAA6KO4iGaEYIuhvQApgB9IkzIElS\nfAYSRhrNJ7QQJEmSJGlHZwCzgXnA9THHki8LgQ+A94C3onXtgInAXGAC0DZt+2GE/TMbGJC3KHNj\nNLAE+DBtXSbfvW/0HvOA/I9GyI7a9sVwwii896JpYNpzxbwvugKTgZnAR8Avo/WleGzUtS+GU+TH\nRlNCF1I3oDmlU1v4lHCgp7sduC5avh64LVo+mLBfmhP203wgWWfSNUx/4Ai2/xFsyHevGgzxFuH8\nFoDxVA9YKCS17YubgV/Vsm2x74v9gMOj5TaE7uU+lOaxUde+yOmxkYQflfQT1jZTfcJaKag5ymsQ\nMCZaHgMMiZYHA48Q9s9Cwv46msL1OrCixrqGfPdjgI7AHlS3rh5Me00hqW1fQO0jAIt9X3xF+FED\nWEs4abUzpXls1LUvIIfHRhISQm0nrHWuY9tikgImAe8AF0frOhC6D4jmHaLlToT9UqUY91FDv3vN\n9Ysprn1yJfA+8Fequ0hKaV90I7ScpuGx0Y2wL6ZGj3N2bCQhISTjnpn5dwLhH3kg8AtC10G6FDvf\nN8W833b13YvdX4DuhC6DL4E/xBtO3rUBngKuAtbUeK7Ujo02wJOEfbGWHB8bSUgIiwkFlCpd2T6j\nFasvo/nXwDOELqAlhL5DCE29pdFyzX3UJVpXTBry3RdF67vUWF8s+2Qp1T98D1DdPVgK+6I5IRk8\nBDwbrSvVY6NqX4ylel8U/bFRiiestYZt18jYHXiDMCrgdqpHWd3AjsWzFoT/HSwgefeyaKhu7FhU\nbuh3n0boJy2jMAuHVbqx/b7omLb8f4CHo+Vi3xdlhD7uP9VYX4rHRl37oiSOjVI7Ya074R9vBmFI\nWdV3bkeoK9Q2vO5Gwv6ZDZyet0hz4xHgC2AToX50IZl996rhdPOBu3IedW7U3BcXEX4IPiD0Ez9L\ndZ85FPe+6Ee4nM0MqodVnkFpHhu17YuBlO6xIUmSJEmSJEmSJEmSJEmSJEmSJElKov8Pdu9OOzHy\n0eYAAAAASUVORK5CYII=\n",
"text": [
"<matplotlib.figure.Figure at 0x7efedbb99c90>"
]
}
],
"prompt_number": 27
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": []
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment