Skip to content

Instantly share code, notes, and snippets.

@abikoushi
Created December 11, 2018 17:19
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 abikoushi/9eff2cb7525ee15b501a0efea9ae5996 to your computer and use it in GitHub Desktop.
Save abikoushi/9eff2cb7525ee15b501a0efea9ae5996 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"import numpy as np\n",
"import tensorflow as tf\n",
"import matplotlib.pyplot as plt\n",
"%matplotlib inline\n",
"a = 1\n",
"b = 0.3\n",
"def SIR(state, t):\n",
" S, I, R = tf.unstack(state) \n",
" dS = -a * S * I \n",
" dI = a * S * I - b * I\n",
" dR = b * I\n",
" return tf.stack([dS, dI, dR])\n",
"\n",
"init_state = tf.constant([1.0, 0.01, 0.0], dtype=tf.float64,)\n",
"\n",
"t = np.linspace(0, 25, num=50)\n",
"tensor_state = tf.contrib.integrate.odeint(SIR, init_state, t)\n",
"sess = tf.Session()\n",
"state = sess.run(tensor_state)"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x1393a7a58>]"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgkAAAFkCAYAAACq4KjhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3XmcjeX/x/HXNWPft0SLNSJLmSFLWSp7kSU0kTVLCBOJ\nJPFtkSJKJFkjWxGSJbv4UmYqJcquxZZQtgZz/f64+EY/28ycM/c5Z97Px+M8ZubMOff9mWFm3ue6\nr+tzGWstIiIiIv8W5nUBIiIiEpgUEkREROSyFBJERETkshQSRERE5LIUEkREROSyFBJERETkshQS\nRERE5LIUEkREROSyFBJERETkshQSRERE5LISHBKMMZWNMfOMMb8aY+KNMfWv4znVjDExxpjTxpif\njDGtEleuiIiIJJfEjCRkBL4BOgPX3PjBGFMA+BRYBtwJjADeN8bUSMS5RUREJJmYpGzwZIyJBxpY\na+dd5TGvAXWstaUvum8akNVaWzfRJxcRERG/So45CRWApf+6bzFQMRnOLSIiIomUKhnOkQc48K/7\nDgBZjDFprbV///sJxpicQC1gN3Da7xWKiIiEjnRAAWCxtfZwUg6UHCEhMWoBU70uQkREJIg1Bz5M\nygGSIyTsB2781303An9ebhThvN0A9etP4dCh4mzZAkePuk/kyAHFi8Mdd0CpUhARAenT+6nyFCY6\nOpo333zT6zJSFH3Pk5++58lP3/PktWXLFlq0aAHn/5YmRXKEhP8Cdf51X83z91/JaYABA4oTERGB\ntfDLL7BxI8TEuNsnn8DYsZA2LVSpArVrQ506UKwYGOOvLyW0Zc2alYiICK/LSFH0PU9++p4nP33P\nPZPky/UJDgnGmIzAbcCFP8WFjDF3An9Ya382xrwK3GStvdAL4V2gy/lVDuOBB4BHgOte2WAM3Hqr\nuzVs6O6zFn78ERYvhoUL4bnnoGdPyJfPBYbateGBByBLloR+hSIiIgKJW91QFvgaiMH1SRgKxAID\nz38+D3DrhQdba3cDDwLVcf0VooF21tp/r3hIEGPcqEH37rBoEfzxhwsLDRvCypXQqBHkzAlVq8J7\n78FffyXlbCIiIilPgkOCtXaVtTbMWhv+r1vb859vY629/1/PWW2tjbTWprfWFrHWfuCrL+CCDBnc\n6MHw4W6EYccOeOstd/+TT0LevPDEE/Dll24UQkRERK4uZPduKFTIhYOFC2HPHujdG5YsgfLl4a67\n4J13/pkMKU5UVJTXJaQ4+p4nP33Pk5++58ErSR0X/cUYEwHExMTE+HSyy7lzLiiMHQvz5kGaNNCk\nCXToAJUqacKjiIgEv9jYWCIjIwEirbWxSTlWyI4kXE54uFsBMXs2/Pwz9O8PX3wB997rllIuXux1\nhSIiIoEjRYWEi+XNC337wrZtbnQhUyY3p6FmTfj2W6+rExER8V6KDQkXhIVBjRqwerXrvbB3L5Qp\nA61audEGERGRlCpQ2zInO2Pg4Yehbl14/3148UWYORN69IA+fSBrVq8rFBERAeLj4eRJOHHC3S5+\n/8QJ2LzZZ6dSSPiX1KndqogWLWDIEBg61E10fOEF6NTJTXYUERG5bta6P+R//nn1219/wfHj/7y9\n0vunk2/fwxS1uiExfvvNBYQJE6BgQRg3zjVoEhGRFOTsWThyxHXuu/h25IhbT3+127Fj7tX/laRL\nB5kzuxbBmTK5W+bMl38/UybImNHdMmT45/2LPo7dto3IKlXAB6sbNJJwDTfd5C4/dO8OXbvC/ffD\noEFu0mNYip/RISIShOLj3R/4Q4fc7fff/3n/wu3wYXe7EAaOHbv8sdKnh2zZLr3lyeNaAl/4OEsW\nd806S5b/f8uc2fdD1L/95rNDKSRcp1KlYNkyGDjQLZ1cvRqmTIEbbvC6MhER+d8f/v373e3AgSu/\n//vv//+VfVgY5MrlfqnfcIN7v2BBt/XwhVvOnJd+nD27GwUIYQoJCZAqFfznP1C5spuzcNddMH26\n+1hERPzk77/h11+vfvvtN4iLu/R5mTPDjTe6V/Z58sDtt7uPc+f+JwhcCAXZs2t4+DIUEhKhZk34\n5huIioL77nPB4dln9f9LRCRRjh1z/fMvvu3d+8/7+/df+viMGeHmm+GWW1wP/sqV3cc33fRPILjx\nRvc4SRKFhES66SZ3+WHAALdN9erV8MEHLpiKiMhFzp1zr/Z37IDt293bC7edOy+93p86Ndx6K+TP\nD8WLuy53+fK5+26+2d2yZFEf/WSikJAEqVLByy+7EPv44+7yw4wZcM89XlcmIpLMrHXX+7du/ee2\nbZsLArt2/XMpICzM/cG/7TYoVw6aNXOB4MItTx4NywYQhQQfqF0bvv7aXX6oWhVefx2io72uSkTE\nD86edX/4t2y5NBBs3frPiEB4OBQuDEWLul+QhQu7UFC4MBQooIYzQUQhwUduuQVWrHBLI59+Gk6d\ncpchRESC1oEDsGmTu333nXv7ww9uIiG4Yf/ixd1yvwYN3Ntixdw8AQWBkKCQ4EOpUrlRhMyZoV8/\nF6affdbrqkREruHcOTcysHHjP6Fg0ybXLwBck56SJd12ua1bu/eLF3eXBjQ3IKQpJPjBCy+4n7k+\nfdyltWee8boiEZHz4uPdXIGNG93tq6/c9dKTJ93nCxeG0qWhc2fXIKZ0aTcyEB7ubd3iCYUEP3nx\nRfez2Lu3Cwo9e3pdkYikSPv3w9q1sH69CwUxMW4PAHCBoGxZd6mgbFm3Ba52s5OLKCT4iTGufXN8\nPPTq5YKCJjOKiF/Fx7vLBmvXwhdfuLc7d7rP5csHd9/tJkuVLQuRka6BkMhVKCT4kTHw0kvu5/bp\np11Q6N7d66pEJGTExcGGDbBmjQsE//2v23AoPNytyX7oIbcm+557XH8BkQRSSPAzY+CVV9wchR49\nXFB46imvqxKRoBQf7yYULlsGS5e6cHDihJstXbGi+yVzzz1QvrzbLVAkiRQSkoEx8Npr7ue7WzcX\nFLp08boqEQkKO3e6QLBsGSxf7jYnSpfOdXHr3x8eeMDNJdDEQvEDhYRkYoxbHhkf77acDguDJ5/0\nuioRCTinT7umK59+Cp99Brt3uwBQrhx07OhCQcWKIb/7oAQGhYRkZAwMHeqCQufObtfRpk29rkpE\nPLdvHyxY4ILB55+75YgFC7o5BTVquFauWnUgHlBISGbGwJtvukZm7drBnXe63UtFJAWx1vUmmD/f\nBYONG93wYqVKrtFKvXquWZEaFYnHFBI8YAy8954bPWzSxC1fzpDB66pExK+shS+/dLvAzZoFv/zi\n2hrXru2WPdWurW1kJeAoJHgkc2b46CO3bPmpp2DcOK8rEhGfuzBiMGMGzJzp5hfceCM88gg0auQm\nH6ZO7XWVIlekkOChkiVh9GjXCr1yZfdWRIKctfD99y4YzJgB27e7CUiNG7ttkatW1UoECRoKCR5r\n1QpWr3YTGSMjXat0EQlC+/bBpEkwebLrepgtGzRsCCNHwv33a8RAgpJCQgAYOdLNW2rSxO21kjmz\n1xWJyHU5c8atShg3DhYudEGgUSMYMgRq1tR2yRL0wrwuQCB9ejeP6bffoEMHN1opIgFs61a3vest\nt7jRgv374e233WjC1Klu6aICgoQAjSQEiKJF4f333SXLKlXUaEkk4Bw/7uYYjB8P69ZBjhzQogW0\nbevWMouEIIWEANK0qWvF3qOHWx5ZtqzXFYkIO3e6UYLx490WyzVquLDw8MOQNq3X1Yn4lUJCgHnj\nDbepW5MmEBurnVxFPGGtm1E8fDjMnet+ELt0gU6d3JbLIimE5iQEmLRp3XLqo0ehTRvNTxBJVqdP\nw8SJbsOkatXgp59gzBj4+We3nasCgqQwCgkBqEABt4pq7lwYNszrakRSgP374cUXIX9+l85vuQWW\nLHH9Dtq3V0tUSbEUEgJUvXrQsyf06+cuiYqIH+ze7S4h5M/vrvU1aQI//uj2U6hRQ3snSIqnkBDA\nBg50HVy7d/e6EpEQs3OnGyEoUgRmz4ZBg9xeCiNHuqVGIgIoJAS0jBndvKlPP3WbxYlIEm3f7pYs\nFi0K8+bB4MGwaxc8+6zrkCgil1BICHANGrjN4bp1g1OnvK5GJEj99JPrgV6smOuM+MYbLhz07OnS\nuIhclkJCgDMG3nrLdWMcPNjrakSCzNatruFR8eKwdCm8+aa71NCjhyYjilwHhYQgUKSI6wD72mtu\ntFREruHgQbdrWsmSsGqVa4a0Y4fblz19eq+rEwkaCglB4rnn3CTGbt3UO0Hkik6dckNut90G06b9\nk6w7d4Z06byuTiToKCQEiQwZYMQIdzl13jyvqxEJMPHx8OGHbs5B//6u18H27W7OgVoniySaQkIQ\nefhhqFPHLYk8edLrakQCxBdfQIUK0Lw5REbC5s0uUefM6XVlIkFPISGIXJjEuG8fvPqq19WIeGz7\ndmjcGCpXdiMJq1a5ngfqcyDiMwoJQea229yS7iFDYNs2r6sR8cCpU64V6R13wFdfwZQp8OWXbo91\nEfEphYQg1KcP3HSTm6itSYySoixbBqVLuz4H/fq5FsrNm0OYfpWJ+IN+soLQhUmMixfDJ594XY1I\nMvj9d9cMqXp1uPlm2LQJBgzQckYRP1NICFL16sGDD7qeMCdOeF2NiJ9Y67ZELVbM9SYfNw5WrIDb\nb/e6MpEUIVEhwRjTxRizyxhzyhiz3hhT7hqPb26M+cYYc8IY85sxZpwxJkfiShZwkxhHjIADB9w2\n9yIhZ9s2txNjq1ZQq5brnti2rXZmFElGCQ4JxphmwFBgAFAG+BZYbIzJdYXH3wNMAsYCdwCPAHcD\n7yWyZjmvcGE3P+H11zWJUUJIXJxLvqVKuRbKixbB1KmQO7fXlYmkOIkZSYgGxlhrJ1trtwKdgJNA\n2ys8vgKwy1r7jrV2j7V2HTAGFxQkiZ591v3u/M9/vK5ExAc2b4Zy5eCFF1xDkO+/d6MIIuKJBIUE\nY0xqIBJYduE+a60FlgIVr/C0/wK3GmPqnD/GjUATYEFiCpZLpU/vRhOmTnUb3YkEJWvd/gqRkXD2\nrFva+Npr2oRJxGMJHUnIBYQDB/51/wEgz+WecH7koAUwwxgTB+wDjgBdE3huuYInnoA8eeDll72u\nRCQR9u+HunXdxiQdOsDGjVCmjNdViQiQyt8nMMbcAYwAXgSWAHmBN3CXHJ642nOjo6PJmjXrJfdF\nRUURFRXll1qDVbp0bjQhOtq1rb/tNq8rErlO8+e7yYjh4fDZZ67vuIhct2nTpjFt2rRL7jt27JjP\njm9sArrxnL/ccBJobK2dd9H9E4Gs1tqGl3nOZCCdtbbpRffdA6wB8lpr/z0qgTEmAoiJiYkhIiIi\nAV9OynX6NBQqBDVrwsSJXlcjcg0nT0KvXjB6tFvP+/77mpgo4iOxsbFERkYCRFprY5NyrARdbrDW\nngFigAcu3GeMMec/XneFp2UAzv7rvnjAAlrL5CPp0rlJjFOmuJb2IgErNtbNPZg40YWEuXMVEEQC\nVGJWNwwD2htjWhpjigHv4oLARABjzKvGmEkXPX4+0NgY08kYU/D8KMIIYIO1dn/SypeLdegAN9yg\nuQkSoOLj3XrdChXcjNuYGOjUSX0PRAJYgkOCtXYm0AsYBHwNlAZqWWsPnX9IHuDWix4/CXga6AJ8\nB8wAtgCNk1S5/D/p07vRhA8+gB07vK5G5CJHj7q9zp991k2eWb8eihf3uioRuYYEzUlILpqTkHin\nTrm5CXXrug62Ip774Qdo0AAOHYIPP9TkRBE/82xOggS+9Omhd2+YNMk1qxPx1OzZUL48pE3reh8o\nIIgEFYWEENSxI+TKpT0dxEPnzrmtnBs3dsHgv//V2lyRIKSQEIIyZPhnNGHXLq+rkRTnyBF46CEY\nPNh1TZwxAzJl8roqEUkEhYQQ1akT5Mih0QRJZt9/7/Ze2LABFi50aVWrF0SClkJCiMqQAZ55xi1F\n373b62okRZg5080/yJjRtVauWdPrikQkiRQSQtiTT0L27BpNED+Lj4fnnoNmzdwyx3Xr3BIbEQl6\nCgkhLGNGN5owYQLs2eN1NRKS4uKgZUt49VUYMsRtR5oxo9dViYiPKCSEuM6dIVs2jSaIH/z5Jzz4\nIMyaBdOnu0Sq+QciIUUhIcRdPJqwd6/X1UjI2LcPqlZ1vQ8WL3aXGkQk5CgkpACdO0PWrG5EWCTJ\nfvwRKlVyHRTXrIFq1byuSET8RCEhBciUCXr0cCsdfv/d62okqP33vy4gZMjg3i9VyuuKRMSPFBJS\niI4dwVoYO9brSiRozZsHDzwAJUrAF1/Arbde+zkiEtQUElKIXLmgeXN45x04c8braiTovPceNGzo\ndg5bssStrRWRkKeQkIJ07w6//ur23BG5LtbCgAFuKKpzZ9diOV06r6sSkWSikJCClC7t5piNGOF1\nJRIUrIXoaBg0yO3D8NZbEB7udVUikowUElKY7t3dfLOvvvK6Eglo1rrZriNGwLvvwrPPqgeCSAqk\nkJDC1KsHBQu6F4Uil2UtdOvm/pOMGeMuNYhIiqSQkMKEh0PXru7S8r59XlcjAcdaeOopGDnSBYQO\nHbyuSEQ8pJCQArVtC2nSuFFkkf+x1iXId95xa2UVEERSPIWEFChbNmjVyoWEv//2uhoJCPHx0KUL\njB4N778PTzzhdUUiEgAUElKobt3g4EF32UFSuPh4t7zx3XddQGjXzuuKRCRAKCSkULffDrVru8nr\n1npdjXgmPh6efNI1Sxo3zl2LEhE5TyEhBeveHWJjYe1arysRT8THu5ULY8fC+PHQpo3XFYlIgFFI\nSMFq1oSiRdVcKUWy1o0gjBvn9hFv3drrikQkACkkpGBhYW5uwpw5sHev19VIsnrhhX8uMbRq5XU1\nIhKgFBJSuFat3FbSo0Z5XYkkm5Ej4aWXYMgQXWIQkatSSEjhMmVyk9nfew9OnvS6GvG7mTPd8NHT\nT0OvXl5XIyIBTiFB6NoVjh6FKVO8rkT8avlyePxxeOwxeP117cUgItekkCAULAj167tW/VoOGaK+\n/hoaNID77nMrGcL0oy8i16bfFAK45ZCbN7sXmxJiduyAOnWgWDH46CPXk1tE5DooJAgA1apB6dJa\nDhlyDhyAWrUga1ZYsMBNQhERuU4KCQK4y9PdusGnn7oXnhIC/vzTjSCcPAmLF8MNN3hdkYgEGYUE\n+Z/HHnObP40d63UlkmR//w2NGrnEt2gRFCjgdUUiEoQUEuR/0qeHli1dA764OK+rkUSLj3cNML74\nAubNc9eRREQSQSFBLtG+vdsdct48ryuRRHv+edcP4cMPoWpVr6sRkSCmkCCXKFEC7rnHNVeSIDRt\nGrz6quum2KiR19WISJBTSJD/p0MH+Pxz2LnT60okQTZudFs9P/449OzpdTUiEgIUEuT/adLETWB8\n/32vK5Hrtm+fa5ZUurQbBlI3RRHxAYUE+X/Sp3cvRsePhzNnvK5Grun0aWjY0LXL/OQTSJfO64pE\nJEQoJMhltW/v+vDMn+91JXJV1kKnTvDtty4g5M3rdUUiEkIUEuSySpWCChU0gTHgDRsGkybBuHFQ\nrpzX1YhIiFFIkCvq0AGWLIFdu7yuRC5r4ULo3Rv69HGdsEREfEwhQa6oaVPInNm9SJUAs3UrPPoo\n1K0LL73kdTUiEqIUEuSKMmaEFi00gTHgHDkCDz8MN98MU6dCeLjXFYlIiFJIkKvq0MGtrluwwOtK\nBICzZyEqCg4dcm0xs2TxuiIRCWEKCXJVd94Jd9+tCYwBo08fWLoUZs2C227zuhoRCXEKCXJNHTq4\njQT37PG6khRu7lwYOhTeeAMeeMDrakQkBVBIkGtq1gwyZdIERk/t3g2tW7umSd27e12NiKQQCgly\nTZkyQfPmLiScPet1NSlQXJxbyZAtm/tHUMtlEUkmCglyXTp0gN9+g88+87qSFKhvX4iNhRkzIHt2\nr6sRkRREIUGuS5kyULasJjAmu7lzXVfF1193M0hFRJJRokKCMaaLMWaXMeaUMWa9Meaq/WCNMWmM\nMS8bY3YbY04bY3YaY1onqmLxTIcOrsnf3r1eV5JC7Nnj5iE0aADdunldjYikQAkOCcaYZsBQYABQ\nBvgWWGyMyXWVp80C7gPaAEWBKODHBFcrnnr0UciQwTVXEj+Li3MzRrNlc99wzUMQEQ8kZiQhGhhj\nrZ1srd0KdAJOAm0v92BjTG2gMlDXWrvCWrvXWrvBWvvfRFctnsic2W0RoAmMyeC55zQPQUQ8l6CQ\nYIxJDUQCyy7cZ621wFKg4hWeVg/YCDxrjPnFGPOjMeZ1Y4w2vQ9CHTrAL7+4vgniJ/PmuX4IQ4Zo\nHoKIeCqhIwm5gHDgwL/uPwDkucJzCuFGEkoADYDuwCPAOwk8twSAyEiIiNAERr+5MA/h4YfVD0FE\nPJcqGc4RBsQDj1lrjwMYY54GZhljOltr/77SE6Ojo8maNesl90VFRREVFeXPeuUa2reHLl3cksib\nbvK6mhBy5oyb+JElC0yYoHkIInJN06ZNY9q0aZfcd+zYMZ8d37irBdf5YHe54STQ2Fo776L7JwJZ\nrbUNL/OciUAla23Ri+4rBmwGilprd1zmORFATExMDBEREdf/1UiyOHYM8uaF/v3dEn7xkWeegeHD\n4YsvoHx5r6sRkSAVGxtLZGQkQKS1NjYpx0rQ5QZr7RkgBvhf43hjjDn/8borPG0tcJMxJsNF992O\nG134JUHVSkDImhWaNHETGOPjva4mRCxa5PZkGDJEAUFEAkZiVjcMA9obY1qeHxF4F8gATAQwxrxq\njJl00eM/BA4DE4wxxY0xVYAhwLirXWqQwNauHezYAatXe11JCDh8GNq2hVq1oEcPr6sREfmfBIcE\na+1MoBcwCPgaKA3UstYeOv+QPMCtFz3+BFADyAZ8BXwAzMVNYJQgVbkyFCkC77/vdSVBzlro3BlO\nn1Y/BBEJOImauGitHQWMusLn2lzmvp+AWok5lwQmY+CJJ2DAAHj7bS3lT7Tp02HmTPdWs0BFJMBo\n7wZJtJYt3YT8Dz/0upIg9csvbhQhKsp1VxQRCTAKCZJoefJAvXq65JAo8fFuHkKGDDBypNfViIhc\nlkKCJEm7dvDNN66DsCTA6NHw+eeuH0KOHF5XIyJyWQoJkiS1a7tL6RpNSIAff3Q9Ebp0gZo1va5G\nROSKFBIkSVKlgjZt3LyEkye9riYInD3rJnPccgu89prX1YiIXJVCgiRZ27auC+PHH3tdSRB49VWI\niYEPPoCMGb2uRkTkqhQSJMkKFYL773cdGOUqYmJg0CC3DbS6KopIEFBIEJ9o1w5WrYKffvK6kgB1\n6hS0aAGlS7tNL0REgoBCgvhEo0auodL48V5XEqCeew527XKXGVKn9roaEZHropAgPpEunXuhPGmS\na7AkF1m+3O3uOHgw3HGH19WIiFw3hQTxmXbtYP9++OwzrysJIMePu5md1apBt25eVyMikiAKCeIz\nd94JZctqAuMlnn8eDh5035Qw/biJSHDRby3xqXbtYMEC+O03rysJABs2wFtvwX/+45aAiIgEGYUE\n8amoKEib1s1NSNHi4tw2mRER0F27ootIcFJIEJ/KmhWaNnWj6/HxXlfjocGDYetW941Ilagd2UVE\nPKeQID7Xrh3s2OH6JqRIP/wAL70EvXu7iRoiIkFKIUF87t57oWjRFDqBMT4e2reHggXVNElEgp5C\ngvicMW404aOP4MgRr6tJZqNGwbp1MHasax4hIhLEFBLEL1q1gnPnXIPBFGPvXujbFzp2hCpVvK5G\nRCTJFBLEL268ERo0gDFjwFqvq0kG1sKTT0KWLNoCWkRChkKC+E3Hjm4O39q1XleSDKZNc60mR492\nSzxEREKAQoL4zf33Q+HCbjQhpP3+u+uF0KQJ1K/vdTUiIj6jkCB+ExYGHTrArFlw+LDX1fhRdLSb\ngPH2215XIiLiUwoJ4letW7tVgZMne12JnyxaBFOmwLBhbiKGiEgIUUgQv8qdGxo1CtEJjMePu4kX\n1au75RwiIiFGIUH8rmNH+PFHWL3a60p8bNAgt8PjmDGuOYSISIhRSBC/q1bNdWAMqQmMmzfDm2+6\nraC1w6OIhCiFBPE7Y9wExo8/dgsBgp610LmzCwe9enldjYiI3ygkSLK4cMl+4kRPy/CNKVPctZOR\nI92+2CIiIUohQZJFrlzwyCPw3ntBPoHx6FE3etC0KdSo4XU1IiJ+pZAgyaZjR9i2DVas8LqSJHj+\neTh50i15FBEJcQoJkmwqV4bixeHdd72uJJFiYlzb5YED4eabva5GRMTvFBIk2VyYwDhnDhw44HU1\nCRQf7yYrligBTz3ldTUiIslCIUGSVcuWkCoVTJjgdSUJ9P778OWXMGoUpE7tdTUiIslCIUGSVY4c\nbs7f2LHuxXlQOHQI+vRxPabvvdfrakREko1CgiS7jh1h505YutTrSq5Tnz5uScZrr3ldiYhIslJI\nkGRXsSKULBkkHRjXroXx4+HVV91GFCIiKYhCgiQ7Y9xowty5sG+f19VcxdmzbrJi2bLQvr3X1YiI\nJDuFBPFEixaQJo17kR6wRo6E775zyx7Dw72uRkQk2SkkiCeyZYNHH3UTGM+d87qay/jtN3jhBejU\nyY0kiIikQAoJ4pmOHWHPHliyxOtKLqNXL0iXDl5+2etKREQ8o5Agnrn7brjzTjeaH1DWrIFp09xq\nhuzZva5GRMQzCgniGWOga1f49FPYvt3ras47dw66dYNy5f7ZulJEJIVSSBBPNW8OOXPC2297Xcl5\n778P33zjCgrTj4eIpGz6LSieSp/ezQ0cPx6OHfO4mCNHoF8/N4JQvrzHxYiIeE8hQTzXuTP8/TeM\nG+dxIQMGQFyca5wkIiIKCeK9vHndcsi33nL9izzx/fdu86b+/V1BIiKikCCBoUcPtxzyk088OLm1\n0L07FC7s3oqICKCQIAEiIgKqVIHhwz04+ezZsHy5O3maNB4UICISmBQSJGBER7v9lL76KhlPeuoU\n9OwJDz4Ideok44lFRAKfQoIEjHr1oFChZB5NeP1114L5zTeT8aQiIsEhUSHBGNPFGLPLGHPKGLPe\nGFPuOp8vWNNWAAAgAElEQVR3jzHmjDEmNjHnldAWHu76GM2cCb/+mgwn3LsXBg+Gp5+GIkWS4YQi\nIsElwSHBGNMMGAoMAMoA3wKLjTG5rvG8rMAkYGki6pQUok0b1zvhnXeS4WTPPON2murXLxlOJiIS\nfBIzkhANjLHWTrbWbgU6ASeBttd43rvAVGB9Is4pKUSWLPDEEzBmDJw86ccTrVzphixeew0yZ/bj\niUREgleCQoIxJjUQCSy7cJ+11uJGBype5XltgILAwMSVKSnJU0/B0aPwwQd+OsHZs26pY8WKri+0\niIhcVkJHEnIB4cCBf91/AMhzuScYY4oArwDNrbXxCa5QUpyCBaFBAzeBMd4f/2Peew+++851b9L+\nDCIiV5TKnwc3xoThLjEMsNbuuHD39T4/OjqarFmzXnJfVFQUUVFRvitSAlJ0NFSuDIsX+3hl4uHD\n8Pzz0LYtlC3rwwOLiCS/adOmMW3atEvuO+bDjXCMu1pwnQ92lxtOAo2ttfMuun8ikNVa2/Bfj88K\nHAHO8k84CDv//lmgprV25WXOEwHExMTEEBERkZCvR0KEtW635pw5XVDwmaeegsmT4aef4MYbfXhg\nEZHAEBsbS2RkJECktTZJqwkTNNZqrT0DxAAPXLjPGGPOf7zuMk/5EygJ3AXcef72LrD1/PsbElW1\nhDxj3GjCkiWwebOPDrp5M4we7fZnUEAQEbmmxFyQHQa0N8a0NMYUw/3RzwBMBDDGvGqMmQRuUqO1\n9oeLb8BB4LS1dou19pRvvgwJRU2auL2WRozwwcGsdf0QChZ0zRhEROSaEhwSrLUzgV7AIOBroDRQ\ny1p76PxD8gC3+qxCSbHSpIGuXd0qh99/T+LBFixwwxJDh2p/BhGR65Soqd3W2lHW2gLW2vTW2orW\n2o0Xfa6Ntfb+qzx3oLVWEw3kunTo4N6OGZOEg8TFuVGE6tVd72cREbkuWv8lAS1XLmjZ0nVgjItL\n5EFGjoQdO9z+DOa6F9eIiKR4CgkS8Hr0gH373KKEBDt0CAYNgk6doGRJn9cmIhLKFBIk4BUv7iYx\nvvxyIkYT+vd3owcD1exTRCShFBIkKLzwAuzZk8DRhE2bYOxYePFFd91CREQSRCFBgkLJkm404aWX\nrnM0wVp3naJoUejc2e/1iYiEIoUECRovvAB798KkSdfx4E8+gRUrYNgwSJ3a77WJiIQihQQJGiVK\nQNOm1zGacPo09OzpNn3w6cYPIiIpi0KCBJUXXoCff4aJE6/yoOHD3YOGDUuuskREQpJCggSVO+6A\nZs2ustJh3z73yS5doFixZK9PRCSUKCRI0Onf3w0UTJhw6f27j+7mYI8OnE2Tin1Pd+B43HESssup\niIhcKpXXBYgk1IXRhEFvHCZLhRWs+nkpn+/8nGybd/LVTOhaF0ZPKAFAmAkjS9osZE2b1b1Nl5Wc\n6XNSvVB1GhZryM1Zbvb4qxERCVwKCRI0Tp89zdq9a1m6cynfVfyc326P5bFPLLfnvJ06hWszcNIK\n4m6Po/nb43jo3AmOnT7Gn3//ybG/j13y/i9//kL04mieWvgUFW+pSKPijWhcvDEFsxf0+ksUEQko\nCgkS8L769SteWPkCK3ev5PTZ0+TOmJvqhaqT6Ycu7FlZnW+/vZW0c6ZD7ChYupR7ClW95jGPnDrC\n/J/m8/GWj3l++fM88/kzlMlThsbFG9P4jsYUy6X5DCIiJhCv2RpjIoCYmJgYIiK0YWRKdfT0Ufot\n68fojaMpdWMpWt3ZiuqFqlMqdymMMWzd6i49jHnzJO3fuB3KlYPZsxN8nuNxx/ls22d8vOVjFvy0\ngBNnTlA8V3Fa39WaLuW6kDFNRj98dSIi/hEbG0tkZCRApLU2NinH0kiCBBxrLdO/n0704mhOnDnB\nsFrD6Hp3V1KFXfrftVgxiIqC4/1fw/59EPPGG4k6X6Y0mWhaoilNSzTl1JlTfL7zc2b9MIt+y/sx\n7L/D6Fe5Hx0iO5A2VVpffHkiIkFDqxskoGw7vI2aU2ry2OzHqJy/Mlu6bKFHhR7/LyBcMLDtHjr9\nNYTYaj2hUKEknz996vTUv70+HzT8gJ+6/kTt22rTY3EPio4syoSvJ3A2/mySzyEiEiwUEiQgnD57\nmhdXvkjJ0SXZ8ccOFjy2gFlNZnFLlluu+rzb3uvN6fTZifruOU6f9m1NBbMXZGKDiXz/5PfcffPd\ntJ3XlpKjSjJz80zibbxvTyYiEoAUEsRzn+/4nFKjS/HKmlfoVbEX33f+nrpF6l77iatWwcyZnH7x\nNXYcyMS4cf6pr/gNxZnVZBYxHWIolL0QzT5qRsSYCBb8tEB9GEQkpCkkiGfOxZ+j84LO1JxSk5sz\n38y3nb7l5QdeJkPqDNfx5HPQvTuUL0/eXs157DF45RV8PppwsYi8EXzW/DPWtFlD1nRZeWjaQ9w7\n4V6+2f+N/04qIuIhhQTxxJlzZ2g+uzljYsYw+sHRrGi1guI3FL/+A4wbB99+CyNGQFgY/fvD/v0w\ndqz/ar7g3nz3srLVSha3WMyff/9JubHleH758/x99m//n1xEJBkpJEiyO3XmFI1mNmL2ltnMajKL\nTmU7YYy5/gMcOQL9+kGrVlC+PABFi8Ljj8N//uM+7W/GGGoWrklMhxier/w8Q9YOocyYMqz/Zb3/\nTy4ikkwUEiRZ/fX3Xzz44YMs27mM+VHzaVS8UcIPMmiQu67w6quX3P3KK3DqFAwY4KNir0Oa8DQM\nqDaA2I6xZEqTiUrjKhG9KJoTcSeSrwgRET9RSJBkc+TUEWp8UIOYfTEsbrGYWrfVSvhBtmyBkSPd\nSELevJd86qabXEB45x3YtMlHRV+nkrlLsq7dOobUGMK7Me9S+t3SrNi1InmLEBHxMYUESRYHjh+g\n2qRqbP9jO8tbLqdy/soJP4i1EB0N+fO7t5fRrZu79NC1q3t4ckoVlopelXqxqdMmbs58M/dPvp+O\n8zty7PSx5C1ERMRHFBLE7/Ye20uViVU4dOIQq1qvIvKmyMQdaMECWLwYhg6FtJfvfpgmDbz1FqxZ\nA9OmJaHoJCiSswgrW69kVN1RfPj9h5QYVYJF2xd5U4yISBIoJIhfbf9jO5UnVCbuXBxr2qyhRO4S\niTtQXJwbPaheHerXv+pDa9SAxo2hVy/466/EnS6pwkwYT5Z7ks2dN1Myd0nqTK3D04uf1goIEQkq\nCgniN98f/J7KEyqTIXUG1rRZQ+EchRN/sLfegl27YPhwuI6VEEOHwtGjbrWDl/JlzcfC5gt5s9ab\njPxyJBXHVeSnwz95W5SIyHVSSBC/2HJoC1UnViVPpjysar3qmu2Vr2r/freioXNnKHF9IxH588Nz\nz8Gbb8LWrYk/tS8YY+hRoQfrn1jP8bjjRIyJYNI3k9StUUQCnkKC+Nzhk4epN60eeTLlYXnL5eTO\nmDtpB3zmGTcH4cUXE/S0Xr0gXz546qnkn8R4ORF5I4jtGEuTEk1oPbc1Lea04M+///S6LBGRK1JI\nEJ+KOxdH45mNOfb3MeZHzSd7+uxJO+Dq1TBlCrz2GuTIkaCnpkvnGjIuXQpz5iStDF/JlCYTEx6e\nwNRGU5n/43zKjCnDl79+6XVZIiKXpZAgPmOtpetnXVn38zpmN51NoexJ3Lr5zBl3iaFiRWjdOlGH\neOghePBBN+fx5MmkleNLj5V6jG86fUOuDLm4Z/w9DFk7RDtLikjAUUgQn3lrw1uMjR3LmIfGJK4P\nwr+9/bZrnvTOOxCW+P+qI0a4aQ2DBye9JF8qlL0QX7T5gp4Ve/Ls0mepNaUWB44f8LosEZH/UUgQ\nn1i4bSFPL3maXhV70aZMm6Qf8LffXPvEzp2hTJkkHapwYejdG4YMgR07kl6aL6UOT83g6oNZ0mIJ\n3x34jrvG3MWq3au8LktEBFBIEB/44dAPPPrxo9QtUpfB1X30cr1nT8iQwWdrGPv2hRtvvGKjRs/V\nKFyDrzt+ze05b+f+yffz6ppXdflBRDynkCBJ8vvJ36k3rR75subjw0YfEh4WnvSDLl8O06fD669D\ntmxJPx4ubwwbBvPnu8aNgShv5rwsbbmUvvf25bnlz1FvWj0OnzzsdVkikoIpJEiiXVjJ8NfffzE/\naj6Z02b2wUHjoEsXuPdet/ezDzVq5Bo2du/uNpEMRKnCUvHS/S/x2WOfseGXDdp+WkQ8pZAgiWKt\npfOCzqz/ZT2zm82mQLYCvjnw8OGwbZubrHgdnRUTwhjXuPHnn93lh0BWp0gdvu74NbdkuYXKEyoz\nfP1wNV8SkWSnkCCJMnz9cMZ9PY73HnqPe/Pd65uD/vyz66z41FNQurRvjvkvxYu7lgvDh8PChX45\nhc/cmvVWVrVeRffy3YleHM0jsx7RjpIikqwUEiTBFm5bSK/Pe9G7Um9a3dXKdwd++mnInBkGDvTd\nMS+je3eoU8e1XjgQ4CsOU4en5o2abzCn2RyW7VxGxHsRxO6L9bosEUkhFBIkQfYc3UPz2c2pW6Qu\nrzzwiu8OvGQJfPSR25kpSxbfHfcyjIEJE9zb1q0hPggWETQo1oDYjrFkS5eNSuMq8e7Gd3X5QUT8\nTiFBrlvcuTiafdSMLGmzMKnBJN+sZAD4+2/o2hWqVYOoKN8c8xpuvBEmToRFi1yzpWBQKHsh1rZd\nS7sy7XhywZM0n92cv/72aC9sEUkRFBLkuj237Dli9sUw45EZ5EifsH0UrmroULcN9MiRPp+seDW1\na7u+Cc8+C19/nWynTZJ0qdLxzoPvML3xdOb/NJ9yY8vx3YHvvC5LREKUQoJcl3k/zmPof4cypPoQ\nyt9S3ncH3rMHXnoJevS47m2gfenVV91po6LgxIlkP32iNSvZjJgOMaQJT0P598sz4esJXpckIiFI\nIUGuac/RPbT+pDUP3/4wPSr08N2BrXUrGXLkgBde8N1xEyBtWpg2zS2sCNRujFdSNGdRNjyxgeal\nmtN2XlvazG3DyTMBtIuViAQ9hQS5qovnIUx4eALGl5cDZsxwLRDfftutavBIsWJuXsLYsW7uZDBJ\nnzo9Y+uPZVKDSczcPJO7x97NlkNbvC5LREKEQoJc1XPLniN2Xywzm8wke/rsvjvwoUNuFKFJE2jY\n0HfHTaR27aBxY2jf3o0qBJuWd7bkq/ZfEW/jKTe2HFM3TfW6JBEJAQoJckX/m4dQYwh333y3bw/e\no4dbe/j22749biIZ40YSMmeGFi3g3DmvK0q4O264g6/af0Wj4o1oMacF7ea240RcEE20EJGAo5Ag\nl3XxPITu5bv79uCffgoffujaHt54o2+PnQTZs8OUKfDFF25CYzDKmCYjkxpMYnz98UzfPJ2yY8uy\n6cAmr8sSkSClkCD/j1/nIRw7Bp06uZaHLVr47rg+UqUK9OsHL74I69Z5XU3iGGNoU6bN/1Y/3D32\nbkZ9NUrNl0QkwRQS5P/x2zwEgN69XVB4991k7YmQEC+8AOXLwyOPuBWawapYrmJseGIDT0Q8QZfP\nutB4ZmP+OPWH12WJSBBRSJBL+HUewooV8N57MGQI5Mvn22P7UKpU8PHHkC6dG/D4I4j/rqZLlY6R\ndUcyp9kcVu5eyV3v3sUXe7/wuiwRCRKJCgnGmC7GmF3GmFPGmPXGmHJXeWxDY8wSY8xBY8wxY8w6\nY0zNxJcs/rL32F7/zUM4edItHahSBTp29O2x/SBPHtey+eBBqF8fTp3yuqKkaVCsAd90+ob82fJT\ndWJVXlr9Eufig3B2pogkqwSHBGNMM2AoMAAoA3wLLDbG5LrCU6oAS4A6QASwAphvjLkzURWLX5w5\nd8Z/8xDAjeH/+qtbQhAWHANYRYu6OZaxscG74uFi+bLmY0WrFfSr3I8XVrxAjQ9q8Ntfv3ldlogE\nsMT8to4GxlhrJ1trtwKdgJNA28s92Fobba19w1obY63dYa3tB2wD6iW6avG5fsv7sfG3jcx4ZIbv\n5yF8+SW8+abbArpoUd8e288qVIDp0+GTT9yqzWCf+5cqLBWD7hvEspbL+PHwj5QaXYpZm2d5XZaI\nBKgEhQRjTGogElh24T7rpkwvBSpe5zEMkBkI4iu9oWXBTwt4fd3rDH5gsG/3ZQCIi4O2baFMGXj6\nad8eO5nUrw+jRrn9p15/3etqfOO+gvexqdMm7itwH00/asrjcx7n6OmjXpclIgEmoSMJuYBw4MC/\n7j8A5LnOYzwDZARmJvDc4gc/H/uZlp+0pF7Rejxd0Q9/xF95BX78EcaNczMCg1THjvD8827HyKkh\n0swwZ4aczGoyi8kNJjPvx3mUHl2a5buWe12WiASQZP2tbYx5DOgP1LfW/n6tx0dHR5M1a9ZL7ouK\niiIqKspPFaYsZ86dIerjKDKmzsjEBhN9Pw/hu+9cSOjbF+4M/ikogwbBL79AmzauB1T16l5XlHTG\nGB6/83GqFqhKq09a8cDkB4iuEM0rD7xCulTpvC5PRK5h2rRpTJs27ZL7jh075rPjm4Q0WDl/ueEk\n0NhaO++i+ycCWa21V2zCb4x5FHgfeMRau+ga54kAYmJiYoiIiLju+iRh+i7ty+vrXmd1m9VUurWS\nbw9+9ixUqgTHj8PXX7vtFkPAmTNQr55rtLR6Ndx1l9cV+U68jWfE+hH0XdaXwjkKM6XhFMrkLeN1\nWSKSQLGxsURGRgJEWmtjk3KsBF1usNaeAWKABy7cd36OwQPAFfvTGWOigHHAo9cKCJI8Fm1fxOC1\ng3nlgVd8HxAA/vMfiIlxlxlCJCAApE7tdoosUsT1UAjmZkv/FmbCiK4Y/b9OjeXfL8+ra17VUkmR\nFCwxqxuGAe2NMS2NMcWAd4EMwEQAY8yrxphJFx58/hLDJKAn8JUx5sbztyxJrl4S5dc/f+XxOY9T\n57Y69KrUy/cnWLUKXnoJBgyAitc1nzWoZMoECxZA+vRQuzYc+PcMnSBXIncJNjyxgV6VevH8iuep\nMrEKP/7+o9dliYgHEhwSrLUzgV7AIOBroDRQy1p76PxD8gC3XvSU9rjJju8Av110G574siWxzsaf\n5bHZj5E2PC2TG04mzPi4Z8Hhw66pwL33uk0QQtSFZktHj8I998COHV5X5FtpwtPwygOvsLr1ag6e\nOMid797J4C8Gc+bcGa9LE5FklKi/ENbaUdbaAtba9NbaitbajRd9ro219v6LPr7PWht+mdtl+yqI\nfw1cOZC1e9cyrfE0cmW4Uv+rRLIW2rVz3RWnToXwcN8eP8AULermJoSFuekXsUm68heY7sl3D5s6\nbaJ7+e70W96Pu9+/m9h9IfiFishlBUfrO/GJz3d8zstrXmbQfYOonL+y708wejTMnQvjx8Mtt/j+\n+AGoYEFYuxby54eqVWHpUq8r8r30qdPzWo3X2PDEBuJtPHePvZu+S/ty6kyQ96oWkWtSSEgh9v21\nj+azm1OjcA363NvH9yfYtMk1S+rSBR5+2PfHD2A33ADLl7srLHXrug6NoajsTWXZ2H4jA6sNZNj6\nYdz57p2s3rPa67JExI8UElKAuHNxNPuoGanCUvFBww98Pw/h5El49FE3/v7GG749dpDIlAnmzXPf\nhqgoGB6iM25Sh6emX5V+fNvpW27IeANVJ1al84LO/Pn3n16XJiJ+oJCQAvRY1IP1v6zno6YfkTtj\nbt+fIDoadu+GGTPc/sopVOrUMHEiPPOM+5b06RP8ez1cSbFcxVjTZg1v13mbyd9OpsSoEsz7cR4J\n6bsiIoFPISHEjY0Zy+iNo3mn7jv+6Ycwaxa89x689RYUL+774weZsDAYMgSGDYPXXoPWrV0DplAU\nZsLoendXNnfeTKncpXh4+sM8NO0htv+x3evSRMRHFBJC2Nq9a+nyWReeLPsk7SPb+/4Ee/ZA+/bQ\npIlb1SD/Ex3tFnhMm+amaJw44XVF/pM/W34WPLaA2U1n8/3B7ykxqgTPL3+ek2dOel2aiCSRQkKI\n+uXPX2g8szEVbqnA8Np+uEB+9iw89hhky+ZGEny970MIeOwx13RpzRq3RHLLFq8r8h9jDA2LN2RL\nly08e8+zvLHuDYq/U5yPf/hYlyBEgphCQgg6ffY0jWY0Ik14Gj5q+hFpwtP4/iQDB8KGDe6lcrZs\nvj9+iKhRwy2RPHMGIiNdngrlv5kZUmdg0H2D/ncJ4pFZj1BrSi11bBQJUgoJIcZaS8dPO/Ldwe+Y\n02yOfyYqrlgBL7/stkUMwbbLvla6NGzcCI8/7racfuQR+OMPr6vyr8I5CvPpY58yP2o+O47soNTo\nUvRZ2ofjcce9Lk1EEkAhIcSM2DCCyd9OZlz9cUTeFOn7E+zcCU2bwv33w7PP+v74ISpDBhgzBj7+\n2GWsO+90W1yEuoeKPsTmzpvpX6U/IzaM4PaRtzPh6wnaNEokSCgkhJBlO5fRa0kvelXsxWOlHvP9\nCY4ehQcfhOzZ3XLHEG+77A+NGsG330LhwnDffdC/f+iufrggXap09K/any1dtnBvvntpO68td425\niwU/LdB8BZEAp5AQInYe2UnTj5ryQKEHGFx9sO9PcOaMGyc/cMDNxsuZ0/fnSCFuvRWWLXO7ab/6\nKlSp4gZoQl2BbAWY8cgMNjyxgZzpc/LQtIe4b9J9fPnrl16XJiJXoJAQAo7HHafB9AZkT5ed6Y2n\nEx7m41f41rp2y6tXw+zZUKSIb4+fAoWHu00y16yB/fvhrrvcksmU8ML67pvvZkWrFXwa9Sm/n/yd\n8u+Xp9lHzdRfQSQAKSQEuXgbT5u5bdh1dBdzH51L9vTZfX+SYcNg7Fh3Ub1aNd8fPwWrWBG++Qbq\n1XM7bNeuDT/84HVV/meM4cGiD/Jtp28ZX388a/eupfg7xXnqs6c4eOKg1+WJyHkKCUHMWkuPRT34\n+IeP+aDhB5TIXcL3J5k71/UZ7tMH2rTx/fGFrFlhyhT45BPYscOthujWLfRXQACEh4XTpkwbtj21\njf/c9x8mb5pM4bcKM2DFAI6cOuJ1eSIpnkJCEHtx5Yu8/eXbjH5wNA2KNfD9CWJjXUegRo3ckkfx\nG2NcZ8bNm+GVV9weEEWKwDvvuL5VoS596vT0ubcPO7vtpENEB4asG0KBEQXov7w/h08e9ro8kRRL\nISFIDfvvMAatHsRr1V+jY9mOvj/Br7+6MfA77oDJk92mBOJ3adNC797w00/QoAE89ZSbr7B0qdeV\nJY+cGXIytNZQdnXfRYeIDgxbP4wCIwrQd2lfDp045HV5IimOfvMHoXGx4+i5pCd97+1L73t6+/4E\nx4+7gBAe7vY/zpDB9+eQq8qTB8aNg6++cg0ta9RwoWF7CpnblydTHl6v+Tq7uu+iS7kuvP3l2xQc\nUZDen/fWnAWRZKSQEGRmbp5J+/nt6Vy2My/f74dLAOfOuRl027bB/PmQN6/vzyHXLTLSrYCYPt1d\n/SlRAnr0gL17va4seeTOmJvB1Qezu8duelTowbsb36XA8AL0XNyT/cf3e12eSMhTSAgiC7ctpMXs\nFjQv3Zy3676N8cemSn36uHAwfbprCyieMwaaNYOtW+H5593Vn0KFXJvnTZu8ri555MqQi5fuf4nd\nPXbzTKVnGPf1OAoML0D7ee3ZfHCz1+WJhCyFhCCxZs8aGs9sTJ0idRhffzxhxsf/dNa6vRjeeMMt\neXzwQd8eX5IsQwbXoXHvXhg61LWtuPNOqFsXVq5MGT0WcqTPwcD7BrK7x24GVB3AZ9s/o+ToktSe\nUpslO5aog6OIjykkBIHYfbE8NO0hKtxSgRmPzCB1eGrfnsBa9xJ1wAC3iqF7d98eX3wqUyb3T7R9\nu1s6+euvrsVz+fLw0UfuilGoy5YuG30r92VX91180PADDp44SK0ptSg1uhTjYsdx+uxpr0sUCQkK\nCQFuy6Et1JpSi2K5ijH30bmkS5XOtyewFnr1cuvuhg6F557z7fHFb1KnhubNXTOmhQtdeGjSBG6/\nHUaPhj//9LpC/0sTnoYWpVsQ0yGGla1WUjhHYdrPb0++N/MxcOVATXIUSSKFhAC2++huanxQgzyZ\n8rCw+UIyp83s2xPEx7s1dsOGwdtvw9NP+/b4kiyMcZ0aly+HL7+EiAjXRTtPHhciPv889EcXjDFU\nLVCVuY/OZWvXrTS5owlD1g0h35v5aP1Ja9b9vE6XIkQSQSEhQH2z/xvuGX8P6VKlY0mLJeRIn8O3\nJ4iPh44dYdQo1265a1ffHl88Ua4czJwJe/a4+QsxMVCzJuTPD337usmPoa5ozqK88+A7/Bz9My9W\ne5HVe1Zzz/h7KDW6FG9teEudHEUSQCEhAC3ctpDKEyqTN1Nevmj7BXkz+3gZ4rlzrsXy+PEwYQJ0\n6ODb44vnbr3VhYItW2D9eqhf32XB4sXd3IXRo0O/7XOO9Dnoc28ftnfbzpIWSyh+Q3F6LunJTcNu\nouWclnyx9wuNLohcg0JCgHkv5j3qTatHtQLVWNl6JXky5fHtCc6ccX0Qpk51s95atfLt8SWgGONC\nwahRsG8fzJoFuXO7q0x587pW0O+/7z4XqsJMGDUK12BWk1n8Ev0LL1Z9kXU/r6PyhMqUGFWC4euH\nq/WzyBWYQEzSxpgIICYmJoaIiAivy0kW8Taefsv6MXjtYLqU68KI2iN8v+VzXBw8+ug/fRAaN/bt\n8SVoHDgAH37odv5et85dfSpXzjXarFfPLa30RxuOQBFv41mxawXvxb7HnC1zMMZQt0hdmpdqzkNF\nH/L9BGGRZBQbG0tkZCRApLU2NinHUkgIAKfPnqb1J62ZuXkmb9R8g+gK0b5vlHT6NDzyiJvF9tFH\n7i+BCHD4MHz2mcuOixbBX3+5yxUPPeT+m9x3H6QL4b+ZB08cZOqmqUz9biox+2LImjYrjYs3pnnp\n5lTNX9X3YV3EzxQSQsjhk4dpMKMBG3/byJSGU2h8hx9e3e/bB02bwsaNbj/iWrV8fw4JCXFxrknT\n/PnutmuXa+JUqRJUqQJVq8Ldd4duaNj6+1Y+/O5Dpn43lZ1HdnJT5puIKhlF81LNuSvPXf7pciri\nY+UvzNsAABRfSURBVAoJIWLHHzuoM7UOR04fYd6j86h4a0Xfn2TVKtfTNzzcXZCuVMn355CQZC38\n8IPrwbBqldtD4tgxt1Nl+fIuMFStChUqQMaMXlfrW9ZaNvy6gambpjJj8wwOnTxE8VzFaVaiGQ2L\nN6RU7lIKDBKwFBJCwPpf1lNvWj2yp8vOwuYLKZyjsG9PYK1rsdy3r3sJOH26m7Emkkjnzrm9Ilav\ndqFh9Wp3qSJVKjefoVIlKFvWbUpVuHDo7C5+5twZlu5cyofff8i8H+fx599/Uih7IRoWa0jDYg2p\neGtF37dJF0kChYQgdi7+HCM2jOC5Zc9R9qayzH10Ljkz5PTtSY4dg9at3aWFvn3dngypUvn2HJLi\nxce7JZarVrnbl1/C7t3uc1mzuqZOkZGXBodgf/Eddy6O5buWM2fLHOb+OJcDJw5wY8Ybefj2h2lY\nvCH3F7yfNOFpvC5TUjiFhCC1448dtJnbhi/2fkGPCj145YFXfD+LetMmt2rh0CG3XWD9+r49vshV\n/P6729I6JsZNgYmJcY2dALJlgzJloGRJuOMOt+31HXdATh9n5ORyLv4c639Zz5ytc5izdQ47j+wk\nS9os1C1Sl7q31aVm4ZrcmOlGr8uUFEghIchYaxkTM4ZeS3qRO2NuJjaYSJX8VXx/og8+cF0UixaF\njz92L91EPPb77y4sxMS4APHDD7BtG5w96z6fO/c/geFCeCha1LWVDpaRB2st3x38jjlb5jDvp3nE\n7nO/lyPyRlC7cG1q31abCrdU8P3mbCKXoZAQRH4+9jPt5rXj852f0ymyE6/XfJ1MaTL59iSnT0OP\nHq6lXuvWrnNO+vS+PYeID8XFuaDwww+webN7+8MP8NNPrt8XuMmQhQq5rHvx7bbbIF++wL6CduD4\nAT7f+TkLty9k8fbFHD51mCxps1C9UHVqF65NrdtqkS9rPq/LlBClkBAErLVM/nYy3RZ1I3OazIyr\nP45at/lh6eE338ATT8D338PIkdCuXfC8/BL5lzNnXHjYtg127Pjntn27u2xxYfQhVSoXFPLlcz0d\nLry9+P2sWb39Wi44F3+O2H2xLNq+iEU7FrH+l/XE23iK5ixKtfzVqFagGlULVOWmzDd5XaqECIWE\nALf/+H46ftqReT/Oo+WdLRlRewTZ0mXz7UmOHHE7+Iwe7fYG/uADNztMJESdPQt7914aHn7++Z/b\nr7+6yZQXZM7sAsPNN7tLF3nz/v+3efO6LbaTM1cfOXWEpTuXsmL3ClbuXsmW37cAbmOqqvmrUq2A\nCw4KDZJYCgkB6mz8WSZ/O5nen/cmPCycMQ+NoUGxBr49SXw8TJwIffq4ywwvvuga8afWtU5J2c6e\ndX3DLoSGvXvd299+g/373ef27YNTpy59XoYMLjTccIO75cr1z/v//jhnThc+fLm888DxA6zas4pV\nu1excs9Kfjj0AwBFchShav6qVLy1IhVvqcjtuW7XUku5LgoJAeZc/Dmmfz+dgasGsu2PbTxa8lHe\nqv0WN2S8wbcn2rjRbem8YYPbpGnIEPdSSESui7X/196dB8dV3Akc//50X6Pb1mVbtuWDQwbbJLap\ntb02ECeEI+RaClyBFHtUSHaXorIVasMGQzbsxkkRkmw4ttjCLGdBQraWJQRsXKFsgsEuZEzkQ47k\nQ7KOkaWxrtFodPX+0TOjmfGM7JHGOn8f16/eez1vZt60W5qf3uvXbYed9icM/uShpcXeENTWZpf+\n9c7O818jIQHy8mzk50deZmfbyx2RIj199DMXzh4ne07v4b1T77G3fi/VrdUYDLlpuawtW8u6eeu4\ndt61rClbQ1563qWrLDVtxTNJmMJdf6a+YTPMb4/+lm3vbePI2SPcsuwWXv3aq6wqWRXfN2pvh+9/\nH555xt4/tmcPbNgQ3/dQahYQsV/g2dn2Kt2F9PeHJg4ul73SF75sbrYdL/3bPT3RXzMpaeQYHA57\nucPhCF4vwuH4OmVZX+dvsiChtIuWhAOcGtpH7bkP+Xn9r3hk4BEAluZdztrSdaxb8Fk+U7qaq4qu\nIj1ZOy2r+NEkYQyMMbx5/E1+8IcfcMh5iC0VW3j21mdZO29tfN9oaMgmBg8+aNd/8Qu4996p3a1b\nqRkkJQVKS23EYmgIurpsdHZGju5um0x0d49Ec3Podk8PeL3ZwPW+ADCQXwvzPuTP8/fx53n7eHHu\nC5A4CMOJJJ27nLSOVWR1rybXs5r8/pU4UrJJT7dnMdLSRo/UVBspKSPr4eF/LCUlNBJ1LqwZR79t\nYmCMYWfdTh567yH2N+5nY/lG9nxzDxvK4/xXvcdjB0J6/HGoqbG3Nf74x1CkA7MoNR0kJo5ckhiv\nwUHo7QW3278U3O6lvvgGbjd09Xqp667mhKeK+oIqGvOqaJXXaBEvAJnepTh6VpLetYJkVyWJ7ZUY\n12K8nkT6+ghEeH+NWCUkhCYNyckXH0lJF14GR2Ji5DJ/efh6cNnFRkJC6DJSWUJCaEQqS0iwZ7Gm\n441nmiRchN6BXl4/8jpPf/w0HzR8wLp569j1jV1cv+j6+E7y4nTCE0/YOxba2+HLX7bJwpo18XsP\npdS0Enx5IrpU4BpfWANDAxxrO8bBloNUNVdxsOUg1a27cXlcAKQlpXHFnCuonFtJ5ZxKKudWcuWc\nSorS59HfL3i99nKL1xs5Bgbs4/39oev+8O8THv79g2Nw0C49npH18OXQkF2/UAwN2b4nU5FI5AQi\nPJmItBztsfCl1xu/Y9YkIQpjDPvO7GPHwR28evhVuvu72bRwE7+783fcuOTG+CYHhw/bswYvvmh/\nI9xzD9x3n46YqJQas+TEZFYUrWBF0QruuvouwP5ec7qdVLdWh8TrR17HPeAGICsli2UFy1hesDyw\nXF5o1+M+ENwlMjxsk4Xg8CcQsYT/daIt/QmJvyw4gsv8+4XvE1wWvE9wefi6f7/gfcOf53TasUXi\nQZOEME3dTTx/6Hme++Q5atprWJCzgPvX3c/dK+9mcd7i+L2RMbB7Nzz2GLz9tr3o+fDDdljleJyj\nVEqpMCJCcVYxxVnF3LD4hkD5sBmmvrOe6tZqjrUdo6athpr2Gnaf3E2ruzWwX5mjjOWFy1mav5SK\nvAoq8itYnLeYxXmLyU4d9VTHhPL/5T1b7wyvqrIT/8aDJgmAZ8DDm8ffZMcnO3in7h1SElP46uVf\n5YkvPsHmRZvje29ybS38+tfw8st2lMSrr7aXFG6/3V7EU0qpCZYgCSzMXcjC3IXcvOzmkMc6+jo4\n3n48kDjUtNew78w+XvrTS/T0j9zGUZBeEEgYKvJs8lCeW055Tjnzc+bHfzI7NSFmZZJgjKG6tZqd\ndTvZeWIne07voW+wj3Xz1vHUTU9x+5W3k5MWxzFda2psYvCb38ChQ3b0lptvtncrbN48PXuzKKVm\nhdy0XNaUrWFNWWjfKGMM7Z526lx1nDh3IhB15+r4oOEDznSdwTDSOWBu5lwW5Cywkb0gsD4/Zz5l\njjKKsopISpiVX0lT2qz5HznrPsuuE7tsYlC3k+aeZtKS0thYvpFHr3uUm5bexPLCi7hx+mIdOTKS\nGFRX2xugb7kFHnoIvvAFmygopdQ0JSIUZhRSmFEY8fZv76CXM11nqO+sD42uet6ue5v6znp6B3pH\nXg+hKKuIMkcZpY7SwLLUUUpZdhnFWcUUZRYxJ3OOJhMTaEbWtDGGkx0nqWqu4kDjAd49+W5g6tYV\nc1ewdcVWtlRsYf2C9fEbeKS9Hd5/H/buhbfegqNH7egot94KP/oRbNmiMzMqpWaN1KRUKvJtv4VI\njDG4PC4auhpo6m6isauRpu4mu97dyEeNH9HY3RjSJwJsMlGYUUhRVhFFmUWB5MG/PSdzTiB5mZMx\nh6yUrPh2NJ9lpn2SMGyGqXXVUtVcxcdNH1PVUkVVcxUdfR0AlDpK2bRwE/etvY8bFt8Qv0lTGhps\nQuCPw4dt+bx5cN11sH07fO5zdnSSaeKVV17hjjvumOzDmFW0ziee1vnEi1TnIkJBRgEFGQWsLF4Z\n9bn9Q/209LTQ0tOCs8eJ0+3E2eO0224nDV0NHGg6gLPHSaf3/HG0UxJTAgmDP3koSC8gPz0/auSl\n55GSqH3EYIxJgoh8B/gnoBg4BPyDMebAKPtvAh4DrgTqgUeNMf8dy3v2DvRy4twJal211LnqqHXV\ncqTtCAebD9Ld3w1AeU45q0tW891rv8s1JdewqmQVxVnFY/mIobq6bBLw6afwxz/apODUKfvYZZfZ\nIZIfeMAuy8unbR8D/eU58bTOJ57W+cQbT52nJKYE+i9ciHfQS7unnbPus7T1tnG217f0bbd52mh1\nt3Ks7RgujwuXxxW49TNcZnImuWm55KTlkJuWa9dTR9b92zlpOThSHGSnZoeEI9UxIxKNmJMEEbkd\n+4X/d8B+4H7gHRFZZoxpi7D/QuBN4EngTuAG4L9EpMkYs2u093rkvUfo+LSDWlctTd1NgfLM5Ewq\n8itYXrCcBzc8yOqS1awqWUVhRmGsHydUb6+9TFBdbZOC6mobDQ328cREWLkSbrvNJgTr18PcueN7\nT6WUUnGRmpQa6MdwsbyDXs71ncPlcXHOcy6QPLg8Ljq9nXT0ddDR10Gnt5PmnmaOth2ls68zUDZs\nhqO+dmpiaiBhcKQ4yErJihiOFAeZKZlkJmeSmZJJRnIGmcm+pa/cv56RnEFyQvKEXUIZy5mE+4H/\nNMY8DyAi3wJuAu4BfhJh/3uBE8aY7/m2a0Rkve91Rk0Sal21XFV+FRsXbGRJ/hIq8itYkr+Eosyi\nsVXQ0JAdHP306fOjttZOUO8fqmvRIrjySti61U6qVFlpZ4SZRpcPlFJKjS41KTUwdkSsjDG4B9x0\nebvo8nbR7e0OrHd5u+juH9nu6e8JiZaelpDt7v7ukI6co0mQBNKT0slIziA9Of289f4z/TF/lmhi\nShJEJBk77ue/+cuMMUZE3gWujfK0dcC7YWXvAI9f6P1e+MoLFzdVtNdrp2hrbT0/nE57JuD0absc\nGBh5Xl6evTxQXm7vPPAnA1dcYe9GUEoppaIQkcDZgHj0dxs2w/QN9uHud9M70It7wLcM2/YMeOxy\n0BO6PmjXm2mOw6ezYj2TUAgkAs6wcicQ7f7B4ij7Z4tIqjEm0ijTaQBHn3nGDljunw7NP2Waf72z\nM/q8rJmZIzOsFBfbSwMlJXa9pMRGtETg+PEoH2Vm6+zspKpqXFOPqxhpnU88rfOJp3U+dgkk4PD9\nCyq003WkRn7OUe9R9rIXfN+l4yEmhpkwRKQEaASuNcZ8FFS+HdhojDnvbIKI1ADPGmO2B5XdiO2n\nkBEpSRCRO4GXYvkgSimllAqx1Rjz8nheINYzCW3AEBA+Z3ER0BLlOS1R9u+KchYB7OWIrcApoC/G\nY1RKKaVmszRgIfa7dFxiShKMMQMi8jFwPfAGgNgehNcDv4zytH3AjWFlW3zl0d6nHRhX9qOUUkrN\nYh/E40XGMnPRz4C/FZG7ROQy4GkgA3gOQET+XUSCx0B4GlgsIttFZLmIfBv4mu91lFJKKTVFxXwL\npDHmNREpBH6IvWzwCfB5Y8xZ3y7FwPyg/U+JyE3Yuxn+ETgD/LUxJvyOB6WUUkpNITF1XFRKKaXU\n7DGWyw1KKaWUmgU0SVBKKaVURFMuSRCR74jISRHxiMiHIvLZyT6mmUpEtonIcFgcmezjmklEZIOI\nvCEijb76vTXCPj8UkSYR6RWRXSKyZDKOdaa4UJ2LyI4I7f6tyTremUBE/llE9otIl4g4ReR/RGRZ\nhP20rcfJxdR5PNr6lEoSgiaP2gasws4w+Y6vo6S6NKqxHVCLfbF+cg9nxsnEdu79NnBeByAReQD4\ne+yEaWsAN7bNT//p4ybPqHXu83tC271OCzk+G4D/ANZiJ/FLBnaKSLp/B23rcXfBOvcZV1ufUh0X\nReRD4CNjzH2+bQEagF8aYyJNHqXGQUS2AV8yxlzEBBlqvERkGLjNGPNGUFkT8FNjzOO+7WzssOV3\nG2Nem5wjnTmi1PkOIMcY85XJO7KZzfeHXSt2JN73fWXa1i+hKHU+7rY+Zc4kBE0etdtfZmwGM9rk\nUWr8lvpOy9aJyIsiMv/CT1HxICKLsJl9cJvvAj5C2/yltsl3ivaYiDwpIvmTfUAzTC72LI4LtK1P\nkJA6DzKutj5lkgRGnzwq9jk81cX4EPgm8HngW8AiYI+IZE7mQc0ixdgfam3zE+v3wF3AdcD3gL8E\n3pIxzT+vwvnq8efA+8YYfx8nbeuXUJQ6hzi09ZgHU1IzhzEmeFzvahHZD5wG/grYMTlHpdSlFXZq\n+7CI/AmoAzYBf5iUg5pZngSuAP5isg9kFolY5/Fo61PpTMJYJo9ScWSM6QSOA9rjeGK0AIK2+Ull\njDmJ/f2j7X6cRORXwBeBTcaY5qCHtK1fIqPU+XnG0tanTJJgjBkA/JNHASGTR8Vlogo1OhHJwjae\nURuaig/fD2wLoW0+G9tbWdv8BBGReUAB2u7Hxfdl9SVgszGmPvgxbeuXxmh1HmX/mNv6VLvc8DPg\nOd9Mk/uB+wmaPErFl4j8FPg/7CWGMuARYAB4ZTKPaybx9e9Ygv0rCuxkZ1cDLmNMA/Y64r+ISC12\navR/xc5v8r+TcLgzwmh17ottwOvYL60lwHbsGbRxT6s7W4nIk9hb624F3CLiP2PQaYzp861rW4+j\nC9W57+dg/G3dGDOlAntv8ynAg51O+jOTfUwzNbDJwBlfXddjp+deNNnHNZMC21FoGHspLTieDdrn\nYaAJ6PX98C6Z7OOezjFanQNpwNu+X5p9wAngKWDOZB/3dI4o9T0E3BW2n7b1CarzeLX1KTVOglJK\nKaWmjinTJ0EppZRSU4smCUoppZSKSJMEpZRSSkWkSYJSSimlItIkQSmllFIRaZKglFJKqYg0SVBK\nKaVURJokKKWUUioiTRKUUkopFZEmCUoppZSKSJMEpZRSSkX0/3Nb04YSsttAAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x139323f28>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.plot(t, state[:,0])\n",
"plt.plot(t, state[:,1])\n",
"plt.plot(t, state[:,2])"
]
}
],
"metadata": {
"anaconda-cloud": {},
"kernelspec": {
"display_name": "Python [default]",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.5.2"
}
},
"nbformat": 4,
"nbformat_minor": 1
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment