{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Understanding the birds and the foxes with eigenvalues" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import numpy as np\n", "import numpy.linalg as la\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 24, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "array([[ 0.2, -1. ],\n", " [ 0.9, 0.3]])" ] }, "execution_count": 24, "metadata": {}, "output_type": "execute_result" } ], "source": [ "A = np.array([[0.2, -1], [0.9, 0.3]])\n", "A" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": true }, "outputs": [], "source": [ "eigenvalues, eigenvectors = la.eig(A)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "array([ 0.25+0.94736477j, 0.25-0.94736477j])" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "eigenvalues" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let's make a time variable:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": true }, "outputs": [], "source": [ "t = np.linspace(0, 10, 1000)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "And compute a solution for each eigenvalue/eigenvector pair as `sol1` and `sol2`, where each is a $1000\\times 2$ array (time steps $\\times$ number of species):" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": false }, "outputs": [], "source": [ "sol1 = np.exp(eigenvalues[0] * t.reshape(-1, 1)) * eigenvectors[:, 0]" ] }, { "cell_type": "code", "execution_count": 19, "metadata": { "collapsed": false }, "outputs": [], "source": [ "sol2 = np.exp(eigenvalues[1] * t.reshape(-1, 1)) * eigenvectors[:, 1]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "And plot the time behavior of the solution for both eigenvalue/eigenvector pairs:" ] }, { "cell_type": "code", "execution_count": 29, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 29, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAD8CAYAAABjAo9vAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd4VFX6wPHvSU9Ib7RAEiCEUEJvUpQmTUUFbKCoIOqq\n2HbVXX/WbTZWRUFEsawgFjoIogiCdEKHhBoghQAhEFJIn/P74wKLSAnJzNyZyft5njyQyZ173wnh\nzZlz3/MepbVGCCGE63AzOwAhhBDWJYldCCFcjCR2IYRwMZLYhRDCxUhiF0IIFyOJXQghXIwkdiGE\ncDGS2IUQwsVIYhdCCBfjYcZFw8PDdUxMjBmXFkIIp7Vp06YTWuuIqx1nSmKPiYkhKSnJjEsLIYTT\nUkodrsxxMhUjhBAuRhK7EEK4GEnsQgjhYkyZY7+UsrIyMjIyKC4uNjsUu/Lx8SEqKgpPT0+zQxFC\nuAiHSewZGRkEBAQQExODUsrscOxCa01OTg4ZGRnExsaaHY4QwkU4zFRMcXExYWFhNSapAyilCAsL\nq3HvUoQQtuUwiR2oUUn9nJr4moUQtuVQiV0IIVxW6RlY9BycOWnzS0liv8ChQ4do2bLlHx4fM2YM\nycnJV33+r7/+yk033WSL0IQQzkxrmPsobJgCmZtsfjmHuXnqyD799NNLPl5RUYG7u7udoxFCOJ21\nEyF5LvR9DeL62fxyMmK/SHl5OSNGjCAhIYFhw4Zx5swZbrjhhvMtEPz9/Xn22Wdp3bo1a9eu5ccf\nf6RZs2a0a9eO2bNnnz/PihUraNOmDW3atKFt27bk5+eb9ZKEEGbK3ARLX4VmN0G3J+1ySauM2JVS\nTwNjAA3sAB7QWle51OO1BbtIPpJnjdDOa14vkFdubnHV4/bs2cPUqVPp1q0bDz74IJMmTfrd1wsL\nC+ncuTPjx4+nuLiYuLg4li1bRpMmTbjzzjvPH/fOO+8wceJEunXrRkFBAT4+PlZ9PUIIJ1BaCDNH\nQ0AduOUDsFOxRLVH7Eqp+sA4oIPWuiXgDtxV3fOapUGDBnTr1g2AkSNHsmrVqt993d3dnaFDhwKw\ne/duYmNjiYuLQynFyJEjzx/XrVs3nnnmGSZMmEBubi4eHjLrJUSN8+u/4dRBuG0y+IXa7bLWyjYe\ngK9SqgzwA45U52SVGVnbysXlhxd/7uPj87t59cuVK77wwgsMHjyYRYsW0a1bN5YsWUKzZs2sH7AQ\nwjEd2WrMrbcbBTHd7Xrpao/YtdaZwDtAGpAFnNZa/3TxcUqpsUqpJKVUUnZ2dnUvazNpaWmsXbsW\ngBkzZtC9++X/QZo1a8bBgwc5cODA+ePPOXDgAK1ateL555+nY8eO7N6927aBCyEch6UCFoyDWhHQ\n73W7X94aUzEhwBAgFqgH1FJKjbz4OK31FK11B611h4iIq/aJN018fDwTJ04kISGBkydP8uijj172\nWB8fH6ZMmcLgwYNp164dkZGR57/23nvv0bJlS1q3bo2npycDBw60R/hCCEewZRpkbYMB/wbfYLtf\nXmmtq3cCpYYDA7TWo89+fh/QRWv9p8s9p0OHDvrijTZSUlJISEioVizOqia/diFcTkkBfNAOQmLg\nwSVWvWGqlNqkte5wteOsUe6YBnRRSvkpY8K5D5BihfMKIYTzWfMBFByDG/9ptyqYi1ljjn09MBPY\njFHq6AZMqe55hRDC6eQfhTUToMVt0KCjaWFYpSpGa/0K8Io1ziWEEE5r9ftQXgJ9XjY1DFl5KoQQ\n1pB/DJI+g9Z3QWgjU0ORxC6EENawZgJUlEKPZ82ORBK7EEJUW0G2MVpvNRzCGpsdjST2C02YMIGE\nhARGjBhhdihCCGey9kMoK4KefzE7EkDa9v7OpEmTWLp0KVFRUWaHIoRwFiUFkPQ5NB8C4XFmRwPI\niP28Rx55hNTUVAYOHMj48eO59dZbSUxMpEuXLmzfvh2AcePG8frrxvLgJUuW0LNnTywWC9nZ2Qwd\nOpSOHTvSsWNHVq9eDUjrXiFqhK3ToeQ0dH3c7EjOc8wR++IX4OgO656zTisY+MZlvzx58mR+/PFH\nli9fzmuvvUbbtm2ZO3cuy5Yt47777mPr1q288cYbdOzYkR49ejBu3DgWLVqEm5sbTz75JE8//TTd\nu3cnLS2N/v37k5KSIq17hXB1lgpY9xFEdTK1bv1ijpnYTbZq1SpmzZoFQO/evcnJyeH06dMEBQXx\nySef0LNnT959910aNzZukixduvR3W+fl5eWRn59/vnXviBEjuP3222WKRwhXs2ex0Za3r2Mt43HM\nxH6FkbU9XKp/zrn2vDt27CAsLIwjR/7XmdhisbB27Vp8fX1/9xxp3SuEi1s3CYIaQrObzY7kd2SO\n/RJ69uzJ9OnTAWOD6vDwcAIDAzl8+DDjx49ny5YtLF68mPXr1wNw44038uGHH55//tatWwFp3SuE\nSzuyFQ6vhs4Pg7tjjZElsV/Cq6++SlJSEomJibzwwgt8+eWXaK0ZPXo077zzDvXq1WPq1KmMGTOG\n4uJiJkyYcP745s2bM3nyZEBa9wrh0jZ+Ap61oN29ZkfyB9Vu21sV0rb392ryaxfCKRWfhvHNoNUw\nYy9TO7Fn214hhKhZtn8HZWeg/QNmR3JJktiFEOJaaG0sSKrbGuq3MzuaS3KoxG7GtJDZauJrFsKp\nZWyE47scdrQODpTYfXx8yMnJqVGJTmtNTk6OLFwSwpkkfQ5e/sb8uoNymBqdqKgoMjIyyM7ONjsU\nu/Lx8ZGFS0I4i6JTsGu20XPdO8DsaC7LYRK7p6cnsbGxZochhBCXt/17KC926GkYcKCpGCGEcHhb\npxt9p+q1MTuSK5LELoQQlXE8BbK2Qut7zI7kqiSxCyFEZWybAcrd2CXJwUliF0KIq7FUGIuS4vqB\nf4TZ0VyVJHYhhLia1F8hPwta3212JJVilcSulApWSs1USu1WSqUopbpa47xCCOEQts0An2CId45G\nftYqd3wf+FFrPUwp5QX4Wem8QghhruI8SFkIbe4BD2+zo6mUaid2pVQg0BO4H0BrXQqUVve8Qgjh\nEJLnQXmRkdidhDVG7I2AbOBzpVRrYBPwpNa60ArnFkKIPyivsHD45Bnyispwd1PUC/YlrJbX+Z3O\nrGrbNxDWBOq3t/65bcQaid0DaAc8obVer5R6H3gBeOnCg5RSY4GxAA0bNrTCZYUQNUlxWQVzt2Ty\nw44s1h88SWm55XdfjwrxpW9CbUZ2aUiTSCst9z+daeyS1OtvYItfGjZijcSeAWRordef/XwmRmL/\nHa31FGAKGBttWOG6QogaoLzCwlfrDvPhsv3kFJYSG16Le7tE07xuIKG1vCirsJB+qoi1B3L4ekMa\nX6w5xE2JdXnppubUDqxmg73kuYCGlkOt8lrspdqJXWt9VCmVrpSK11rvAfoAydUPTQhR0+0+mscz\n324jOSuP7k3Cebx3EzrHhl5yymV091hOFJTwxepDTPktlRV7snlrWCIDW9WtegA7ZkLdNhDWuBqv\nwv6sVcf+BDBdKbUdaAP8y0rnFULUULM3Z3DrxNVkF5QwaUQ7vhrdiS6Nwq44jx7u782f+8fz89M9\naRzpz6PTNzP+pz1Vawd+MhWObHa60TpYqdxRa70VuOo+fEIIcTVaa95duo8Jv+yjS6NQJtzdlsiA\na5tSiQ6rxbcPd+Hlubv4YNl+TheV8erNLXBzu4Z58p2zjT9b3HZN13YEDtO2VwghLBbN6wuT+WLN\nIYa3j+Lft7fCw71qEwveHu68MbQVQX6eTFmZSi1vD54f0KzyJ9g5Cxp0geAGVbq+mSSxCyEcgtb/\nS+qju8fy4qCEaxthX4JSir8ObEZBSTkf/XqAesG+3Nsl+upPPJYMx5Nh0DvVur5ZJLELIRzCpF8P\n8MWaQzzYLZb/G5xgtZp0pRSv39KCY6eLeXX+LprVCaBjTOiVn7RrNig3aD7EKjHYmzQBE0KYbvbm\nDN5esodb29SzalI/x8PdjffuakODEF/GzdjCqcIrLI7X2qiGie0J/pFWjcNeJLELIUy1I+M0L8ze\nQddGYbw1rHW1p18uJ8DHkw/vaUdOQSl/mbnt8pUyR7bAqYNOWQ1zjiR2IYRpThaW8si0TUT4e/Ph\nPW3x8rBtSmpZP4jnBsSzNOU487cdufRBO2eBmyck3GzTWGxJErsQwhQWi2bcjC1kF5QweWR7wvzt\n0znxgW6xtGkQzGsLkskpKPn9F7WG5PnQuDf4htglHluQxC6EMMXUVQdZtf8Er9/SglZRQXa7rrub\n4s2hieQXl/HPRSm//2LWNjid5tSjdZDELoQwQfKRPN5esof+LWpzZ0f714nH1wlgdPdGzN6cybb0\n3P99IWWBUQ0TP8juMVmTJHYhhF0Vl1Xw1LdbCPLz5N+3J9qm1W4lPNarMeH+Xvx9YfL/bqTuXgjR\n3aBWmCkxWYskdiGEXU34ZR97jxXw1rBEQmt5mRZHgI8nz94YT9LhUyzacRRO7Ifs3dDsJtNishZJ\n7EIIu0nJymPKylSGtY+iV7z5NeJ3dGhAXKQ/7y7diyVlgfFgs8HmBmUFktiFEHZRYdH8dfYOAn09\neXFQgtnhAMaN1Cf7xrH/eAG5m2cbLXqdsDfMxSSxCyHsYvr6w2xNz+WlmxIIMXEK5mKDWtbluohS\nQk9tx+IC0zAgiV0IYQfH84p568c99IgL59Y29c0O53fc3BR/a5wKwHI6mRyNdUhiF0LY3Js/7qG0\n3MLfh7Q0rQrmSlqcXkm6W33e3kzVNuVwMJLYhRA2tS09l1mbM3igewwx4bXMDuePzpxEHV5FQewA\ndh8rYNX+E2ZHVG2S2IUQNqO15rUFuwj39+bxXk3MDufS9i4BSzmNe95FRIA3n/520OyIqk0SuxDC\nZuZvO8LmtFyeGxBPgI+n2eFc2u6FEFAPrwYdGNU1mhV7s9l7LN/sqKpFErsQwibOlJbz70W7aVU/\niGHtoswO59JKz8D+X4zadTc3RnSOxsfTjalOPmqXxC6EsInPVx/iaF4xL9/c3GY91qvtwC9QXgQJ\nRpljSC0vbmsbxbxtmZwuKjM5uKqTxC6EsLrTZ8r4eMUB+iZEXn0bOjOlLASfYKM/zFkjOjekuMzC\nnM0ZJgZWPZLYhRBW99GKA+SXlPPn/vFmh3J5FWWwdzHEDwT3/83/t6wfROuoIKavT3Pa0kdJ7EII\nqzp6upjPVx/ktjb1aVYn0OxwLu/QKig+fcmmX/d0bsi+4wUkHT5lQmDVZ7XErpRyV0ptUUottNY5\nhRDOZ8KyfVi05ul+Tc0O5cp2LwQPX2O3pIvc3LoeAd4eTFt32ITAqs+aI/YngZSrHiWEcFkHTxTy\n7cZ07unUkAahfmaHc3kWC+z+AeL6gtcf4/Tz8uDWtvVZvPOoU95EtUpiV0pFAYOBT61xPiGEc/rP\nz3vx9nDj8d5xZodyZZmbID8Lml1+C7yh7aMoLbeweEeWHQOzDmuN2N8DngMslztAKTVWKZWklErK\nzs620mWFEI5i37F8Fm4/wv3XxRARYJ+Nqats9wJw84CmN172kNZRQTSKqMXszZl2DMw6qp3YlVI3\nAce11puudJzWeorWuoPWukNERER1LyuEcDAfLNuPr6c7Y3o0MjuUK9PaKHOM6QG+IZc9TCnF0HZR\nbDh0kvSTZ+wYYPVZY8TeDbhFKXUI+AborZSaZoXzCiGcxP7jBSzYfoT7usaYut1dpWTvhpMHzi9K\nupJb2xothudsca5Re7UTu9b6r1rrKK11DHAXsExrPbLakQkhnMbE5fvx8XDnoR6xZodydSlnC/fi\nr74FXv1gX7o0CmX25gynqmmXOnYhRLWkZhcwb2sm93aNJszfwefWwZhfj+oIgXUrdfjt7aI4lHOG\nzWm5Ng7Meqya2LXWv2qtXWNvKSFEpUxcfgAvDzcecvS5dYDcNMjadslFSZczsGUdvDzc+GG781TH\nyIhdCFFlh3MKmbs1kxGdox2/EgaM2nWAhMuXOV4swMeTnnERLNqRhcXiHNMxktiFEFU2cfl+3N0U\nD/d0gtE6GPPrkc0hrPE1PW1wYh2O5hWzJd05WgxIYhdCVEnGqTPM3pzJPZ0aEhnoY3Y4V1d4AtLW\nXNM0zDl9E2qfnY45aoPArE8SuxCiSqauMjajeMhZRut7FoG2VKrM8WLONh0jiV0Icc1yz5TyzYZ0\nbmldj/rBvmaHUzkpCyGoIdRJrNLTnWk6RhK7EOKaTVt3mKKyCsZe7ySj9ZJ8SF1ujNZV1XZzcqbp\nGEnsQohrUlxWwRdrDnFDfIRj91u/0L6foaK0SvPr5xjTMeEs2XXU4RcrSWIXQlyTWZszOFFQylhn\nmVsHo/e6Xzg07FKt0/RrXpvM3CJSsvKtFJhtSGIXQlRahUXzycpUEqOC6NoozOxwKqe8BPb+ZGyB\n5+ZerVP1blYbpWBpyjErBWcbktiFEJX2c/JRDuWc4eGejVFVnKu2u4MroTT/mhYlXU5EgDdtGgRL\nYhdCuAatNR+tSKVhqB8DWtYxO5zKS1kAXgEQe71VTteveW22Z5zm6Oliq5zPFiSxCyEqZcPBk2xL\nz+WhHrG4uznJaN1ScXYLvH7gaZ1FVP0SagOOPR0jiV0IUSkfr0wltJYXw9o3MDuUyktfD2dOVGlR\n0uU0ifQnOsxPErsQwrntPZbPst3HGdU1Bl+v6t2AtKuUheDuBU36We2USin6JtRmzf4cCkrKrXZe\na5LELoS4qikrU/H1dOe+rtFmh1J5Whu91xvdAD7Wrbfv17w2pRUWVu51zP2bJbELIa4o63QR87Zm\ncmfHBoQ4+rZ3Fzq6w+i/Xo1FSZfTITqEQB8Plu8+bvVzW4MkdiHEFX2++hAVFs3o7k6w7d2Fdi8E\n5Qbxg6x+ag93N3rERbBib7ZDrkKVxC6EuKy84jK+Xp/G4MR6NAj1Mzuca5OyEBp0Af8Im5z++qYR\nHM8vYfdRx1uFKoldCHFZX69Po6Ck3Hk20jgn5wAc32XVapiL9Wxq/MJY4YDz7JLYhRCXVFJewWer\nDtK9STgt6weZHc612b3Q+NMG8+vn1AnyoVmdAFbskcQuhHAS87Yc4Xh+iXM1+zonZaHRdz3EtlU8\n18dHkHT4pMOVPUpiF0L8gcWi+XjlARLqBtIjLtzscK5N/lHI2GCV3jBXc33TCMoqNGsP5Nj8WtdC\nErsQ4g+W7T7OgexCHrm+kfM0+zpn9w/GnzachjmnQ3Qofl7urNjrWGWP1U7sSqkGSqnlSqkUpdQu\npdST1ghMCGGej1ceoH6wL4Na1TU7lGu3eyGENoLIBJtfysvDjesah/PrHscqe7TGiL0ceFZrnQB0\nAR5TSjW3wnmFECbYdPgkGw+dYnT3WDzdnexNfdEpo01vws1V3gLvWl0fH0HGqSJSTxTa5XqVUe1/\nNa11ltZ689m/5wMpQP3qnlcIYY6PV6QS5OvJnR2dqNnXOXuXgKUcEobY7ZI9z96DWLP/hN2ueTVW\n/XWslIoB2gLrrXleIYR9HMgu4OeUY9zXNZpa3h5mh3PtUhZAYH2o19Zul2wY6kf9YF9W73ecG6hW\nS+xKKX9gFvCU1jrvEl8fq5RKUkolZWc7Xt2nEAI+/S0VT3c3Rl0XY3Yo1660EPYvNW6autlvCkkp\nxXWNw1ibmkOFxTHm2a3y6pVSnhhJfbrWevaljtFaT9Fad9Bad4iIsM0SXyFE1R3PL2bWpkyGt48i\n3N/b7HCu3f6lUF5slzLHi3VrEs7pojKSj/xhTGsKa1TFKGAqkKK1/k/1QxJCmOGL1Ycos1gY08MJ\nFySBMQ3jFwYNu9r90tc1Njb2Xn3AMebZrTFi7wbcC/RWSm09+2H9dmpCCJspKCln2rrDDGhRh9jw\nWmaHc+3KS4wbp/GDwN3+9wYiA32Ii/RnjYMsVKr2d0BrvQpwshUMwua0NkrPzuRAaQGUnjFaqHp4\ngacf+NcG3xC7laSJK/tmQxp5xeXO2T4AjBLHkjxIuMW0ELo1CefbjemUllvw8jC3TNQJb3sLh3Pm\nJGRuhqwtkLUNTuyH0+lGQr8SD18Ibgh1Whp9PRp0gqiO4O5pn7gFAKXlFqauOkin2FDaNgwxO5yq\nSZ4HXgHQ6HrTQujaOIwv1hxiS9opOjcKMy0OkMQuqqKiHNLWwoFf4MByI5lzthogtDFENDO2Iwtu\nALUiwMsfvPyMUXx5CZQVQv4xyMuEkwchfQPsnGU83ysAYnsYN8ASbgbvAJNeZM2xYNsRsk4X86/b\nWpkdStVUlMOeRdC0P3iYd9O3S6Mw3BSsPpAjiV04CYsF0tcZCXjXXGPndzcPiOoEvf5m3LCqmwg+\nVWzveuYkHF4NB5bBvqXGf9SFzxj9tDuNNUbzwurONfuKrx3ADfFOWq2WttaY8jOhGuZCQb6etKof\nxNoDJ6BfU1NjkcQurqwgG7b8F5K+gNNpxvRJ/ABocTs07mW9EbVf6P9G6Vobo/jt38LOmbDje2OK\n5rpxdl0qXhMs33OcvccK+M8drZ2v2dc5KQvAwwea9DU7Eq5rEs4nK1MpLCk3dYGXJHZxaekbYP3H\nxtylpQxie0KflyF+IHj72/baSkHDzsZHv9dh69ewbiJ8d6+xorDvq8ZUj6i2j1ekUj/Yl5tb1zM7\nlKqxWIzE3riP7X8uK6Fb43A++vUAGw+d5Ib4SNPicLIOP8KmtDYWeXw+GKb2g30/Q8cx8NhGGLUA\nEofb/z+Ptz90HgtPbIYhk4x3EP8dAtOHw6lD9o3FxWw6fIoNh046Z7Ovc45sgfwjpk/DnNMuOhgP\nN8WGgydNjUNG7MIY9exeCL+9Y9wIDagHA96AdveBl4PUNLu5Q9sR0HIobJgCK96EiV3ghueh6+NS\nSVMFH6844LzNvs5JmW/c62na3+xIAPDz8iAxKoj1Jid2J/01LaxCa+Nm5Se9jGmOkgK45QN4cht0\nedRxkvqFPH2g2zh4bD006QNLX4WpNxqbF4tK23/caPY1ylmbfYHx85syH2J6GPdoHESn2DC2Z+RS\nVFphWgxO+i/q3IrLKjh4opAD2QUcyS0iO7+E7PwSThSUUlRWQWm5hZLyChQKb083fDzc8fN2JzLA\nm8gAH2oHetMowp+42v5E+HtX7aZXxib45VVjYUdQQ7h1MiTeYYyMnUFQFNw13ajQWfAkTO4Bg96C\nNiPk5molfLIyFS93N+5zxmZf5xzbBSdTjXdsDqRzbCiTVxxgS/oprmtszraCkthtrKCknO3puWxJ\nz2Vrei57juaTfuoMF2624u3hRmSgN2G1vKnl7U6QrydeZ+c8i8srKCmzkFNQSkpWHtn5JVzYQC7E\nz5OW9YPoFBNKp9hQWjcIxsfzCsn55EH4+WVjpOMXZky5dHjQ1PrfamlxK0R1gDmPwLzHjF9UN78P\nnr5mR+awjuUVM2dLJnd2bOCczb7OSZ5rrGY2cbXppbSPCUEpWJ96UhK7qyguqyDp0Cl+25/Nb3tP\nkHI073wSbxRei8SoIG5vV5/GEf40jvAnKtSXAG+PSo+6Kyya7PwSDmQXsOdoPnuP5bM1PZfxP+8F\nwMfTjR5xEdzYvDZ9E2oTUsvLeGJZEax6D1a9a8xJXv8CXPe4aywACoqC++bBb+Nh+T8he48xmg+K\nMjsyh/TZ6oOUWyw85KzNvsCYhtk1x5iG8Xes+vtAH0+a1w009QaqJHYrOJ5fzJJdx/g5+RjrU3Mo\nKbfg4aZoFx3CuN5xtG0YTJsGwQT7eVX7Wu5uijpBPtQJ8qFbk/+NBnLPlLLx0ClW7cvmp2QjFnc3\nRc8mYTwetZ92yW+icg8bNx9v/AcEOml52+W4ucP1z0GdVjDrIZhyA9w5DRp2MTsyh5JXXMbX69IY\n1KouDcP8zA6n6o7thJz9DjcNc07n2DCmrz9sWt8YSexVlJlbxI87j/LjziySDp9Ca4gNr8XdnRrS\ns2k4nWPD7HpTKtjPi37Na9OveW1evaUFOzJPs2bjRlpuf4H2aZs5QBQbEybS+8ZhRAb62C0uu4sf\nCA8tgxl3GWWRQ6caq1cFAF+tPUx+STmPXN/Y7FCqZ9ccUO4OU+Z4sU6xoXy2+iA7MnNpH23/G7uS\n2K9BXnEZi7ZnMWtzBhsPnQKgWZ0AnuwTx8CWdWla298hVu+psiIS904kced7aE8v9iU8zzune/Pz\n1hw8ti/n9nb1eahnIxpHmL+gwyYimsLon+HrO4xqn8HjjfsINdyZ0nKmrjrIDfERtKxfxdYPjuDc\nNExsT6hlzhz21XSMMZqprUs9KYndEVVYNKv2n2DWpgyW7DpKSbmFxhG1+Ev/eAa1qutYvau1ht0/\nwI9/NZb/txqO6vd34gLr8jFw6EQhn65K5fukDL5NSue2tvV5pl9TokKc+C355dQKg1Hz4fv7YeHT\nUJgD1//F7KhM9fX6NE4WlvJE7yZmh1I9R7cb1TDdnjI7kssK8/cmLtKfDQdP8lgv+19fEvtlZJw6\nw4wNaXyflMHx/BKC/YyFHLe3i6J1VJBDjMx/J+cALH7OWDka2Rzu/wFiuv/ukJjwWvzj1lY81bcp\nU1am8sWaQyzclsWILg15onccobWqfw/AoXjVgrtmwPzHYfk/jN3rb3ihRpZDFpdVMGVlKl0bhZky\ngrQqB5+GOadzo1DmbjlCeYUFDzuv7JXEfoEKi2bl3mymrTvMsj3HUUDvZpEMax9Fr2aReHs4YI13\naaFRDbLmA3D3hv7/hk4PXXElZri/N38blMAD3WJ47+d9fLnmEHO2ZPJc/2bc1bEBbm4ulPjcPWDI\nRCMRrHgDdAX0erHGJffvNxkDlPfubGN2KNVzbhqm0Q0OtSjpUjrFhjFtXRopWfm0irLv1JckduBE\nQQnfJaXz9fo0Mk4VERHgzeO9mnB3p4bUC3bQemitjeZHS/5mbGqReKfRMCugTqVPUTfIlzeHJfJg\n91hemreTv83ZwbdJ6fxjSEu7/yDalJu7saLWzQ1Wvg2WCqOhWQ1J7mUVFib/eoB2DYPp2tjcPuHV\nlrXV6BHU489mR3JVnWKMXzzrD+ZIYrcXrTUbD51i2rrDLN6ZRVmFpmujMP46MIEbW9R27KZIJ/bD\n4r8Y7QDs/647AAAb/klEQVQiW8ADiyH6uiqfLr5OAN+O7cLcrZn884fdDJm4irE9G/N0vzjHfJdS\nFW5ucNP7xoKWVf8x6vd7PGN2VHYxZ0smmblF/OPWlo43hXitds0x1mE0G2x2JFdVJ8iH6DA/Nhw8\nafcNwmtcYs8vLmPOlkymr0tjz7F8Anw8GNE5mpFdGtIk0sEX65QWwsp3jGkXT18Y8KbRfdEKm/cq\npbitbRS9m9Xm34tSmLziAMt2H2P88DauM3p3c4PB7xrfx19eM/Zc7fCA2VHZVIVFM2n5flrUC3Te\njTTOOT8N08vhp2HOaR8dwsq92Wit7fpLtcYk9uQjeUxbf5i5WzI5U1pBq/pBvDm0FTe3roefl4N/\nG7SGXbPhp5eM7eRa3wP9XgN/6/d7DvL15I2hifRvWYcXZm3n1kmreaJ3E57oHYe7K8y9u7nBrR9B\nUa5RLeMbYrQlcFELtx/hUM4ZJo9s5/yj9SObITfNWDXtJNpHhzB7cyZpJ88QHWa/CjoHz2jVU1Je\nweIdR5m27jBJh0/h7eHGza3rcW+XaFo3CDY7vMo5tgsWPw+HfjNWVQ77zC6rKXvFR/LTU9fzyvyd\nvLd0H+tSc3j/rrbUdoXFTe6ecMd/4avbYNYYYzu/xibUpNlYhUUzcfl+4iL9ubF55e+9OKwdM8Hd\nC5oNMjuSSmsfbdSzJx06JYm9utJPnmH6+jS+T0onp7CUmDA//m9wAsPaR1llWb9dFOXC8n/Bxk/B\nJxAG/wfa32/X7otBfp68d1dbusdF8NLcnQx6/zfevbMNPZs6+Vt6MDbXvudb+HwQfHcfPLgEajc3\nOyqrWrj9CHuPFTDh7rbOX+lkqTD224270XiX5SSaRgYQ4O3BprRTDG1vv95FLpPYKyyaFXuPM21d\nGsvPlir2TajNvV2j6dY43Hl+sC0W2PKVMQdcdAraPwC9/8/UOcVh7Y3a/ce+3syozzfweK8mPN23\nqfN8Ty/HNxhGfAef9IGv74SHfrHJ9JYZyissvL90H83qBHBTq7pmh1N9B1dCwTFoNdzsSK6Jm5ui\nbXQImw+fsut1rZLYlVIDgPcBd+BTrfUb1jhvZWScOsPMTRl8n5RBZq5RqvhErybc5cilipeTugJ+\nfsnYxahBF6O/eN3WZkcFQFztAOY91p1X5u/kg2X7ST6Sx7t3tSHQx8l3LgqKgnu+MUbuM+6CUQuN\n0byTm7Mlk9QThUwe2d75fwGDsaG5d6DD7JR0LTpEh/Du0r2cLiojyNc+/1+qndiVUu7ARKAfkAFs\nVErN11onV/fcl1NcVsFPycf4PimdVftPAMYmsn8d1Iz+Leo4dqnipRxPMXqk7/sJghrA7Z8YIxMH\nu9nl6+XOm0MTaRUVzGvzd3HbxNV8cl8HGjl7z5l6bY3v+bcjYc7DMPxL4yarkyqrsDBh2T5a1g+k\nf4vaZodTfWVFkDwfmt/ilH3220eHoDVsTc/lejtNY1pjxN4J2K+1TgVQSn0DDAGsnthTsvL4ZkMa\nc7ce4XRRGfWDfRnXO45h7aNoEOqEo6z8o0b/8C3TwCvAWGDU6WFj+zcHpZTi3i7RxEX686fpmxky\ncTUT7m5LLxN3ZLeKhJvgxr/DT/9ntB/o87LZEVXZ90kZpJ8s4vX7XaBuHWDvEijNd7ppmHNaNwjG\nTcGmQyedKrHXB9Iv+DwD6HzxQUqpscBYgIYNG1bpQtPXH+a7pAz6t6jDHR2inGvu/EIF2bDmfdg4\nFSrKoPMj0PMvTlObC9ClURjzH+/G2P9uYvQXG3n1lhbc1zXG7LCqp+vjxiYdv403psCaDzE7omtW\nUl7Bh8v20bZhsPPXrZ+z43vwr2N0c3RC/t4eJNQNZFOa/ebZrZHYL5VZ9R8e0HoKMAWgQ4cOf/h6\nZYzrHcefb4x3nsqWixXmGAl9wydQXmyMQG54AUKdcyebqBA/Zj7alSe/2crL83aRfvIMfx2Y4Jy/\nbMGY+ho83pgam/snCI+HyGZmR3VNvtmQzpHTxbw1rLVrjNaLThlTlB0fcp79eC+hfXQIszZl2K0h\nmDWukAE0uODzKOCIFc77B5GBPs6Z1POPws+vwHutYPUEYzn0n9bD7VOcNqmf4+flweSR7RnVNZpP\nfjvIEzO2UFxm3u7s1ebhDXd+BZ5+8M09UHza7IgqrbCknA+W7adzbCjdmjh5T5hzkudBRSkkOuc0\nzDnto0MoLK1g99F8u1zPGol9IxCnlIpVSnkBdwHzrXBe53c8BeY+djahvw/xA+Cx9TD0U2MzCBfh\n7qZ49ZYW/N/gBH7YkcXIT9dzqrDU7LCqLrAe3PEl5B6G2WONElQn8MlvqZwoKOGFgc1cY7QOsP17\nCIuDus7dlfLcQqXNdpqOqXZi11qXA48DS4AU4Dut9a7qntdpWSxGT/Tpw2FSF2NRRbtRMG6zsWo0\nIt7sCG1CKcWYHo2YNKId2zNPM/SjNWScOmN2WFUXfR0MeAP2/ggr3jQ7mqs6nl/MlJWpDG5Vl7YN\nnWcBzxXlpsHh1ZB4h8NViF2r+sG+1A70ZpOd6tmtUseutV4ELLLGuZxW/lGjumXzl8YPZK0I6PV/\n0HG0U90Ura5BreoSGeDNg19sZNhHa/lqdCfiajt4c7XL6TgGjmwxEnuDTtCkj9kRXdZ7S/dRWm7h\nL/1daOCw7RtAQ+u7zI6k2pRSdIgOJemQk4zYa7TSQmNEPuMeeLcFLPs7BEcbI/OndxlbsdWgpH5O\nh5hQvn24KxVac8fHa9mWnmt2SFWjFAx6ByITjCmZvCyzI7qk/cfz+XZjOiO7RBPjSFs1VofFAlun\nG5UwwVWronM07aJDyMwt4ujpYptfSxL7tSotNPYVnTka3m4CMx+EzE3Q5VF4fBPcvxBaDjVuwtVg\nCXUDmflIV/x9PLjnk3WsPruQzOl4+RkLlsqKYNZoqCg3O6I/eGPxHvw83Z1/L9MLpa0xNtRoM9Ls\nSKzm3Dz71nTbj9pdpleMzWht7Ce6fynsWwKHVhl36X1DjV2LWg415mOduBTLVqLDajHzkeu4b+oG\nHvh8IxPubsuAlk7YZTCiKdz8Hsx+CH79l0MtXlqfmsPSlGP8pX88Yf4uNJjYMt1YtOfg+5peixb1\nAln27PXE2uFdlST2i5UWGtUs6RsgbS2krYPC48bXwptCp7HQpK+xUfQV9hUVhtqBPnz7cBce+GIj\nf5q+iTeGJnJHhwZXf6KjSbzDaJ3823jjF3mTvmZHRIVF89qCZOoG+fBgt1izw7GekgKjzLHVUJfo\n23OOp7ub3dpvOFdizzkAxbnGjclaEVXvG1FRBnlHjJucuWlGWdvxFKP3+clUzq+vCo6Gxr2N/ueN\nboBQF/rPY0fBfl5MH9OZh7/axHMzt5NfXM7o7k74vRz4FmRsMubbH/4NguqbGs7XG9JIzsrjw3va\n4uvlQu8Yk+dCWaFLTcPYm3Ml9rUTIWnq/z73CoBa4eDtDx6+Ro8VTz9jJG2xgKXc+KgoNRaaFOUa\nK9lKL14koIykXbuFMTKr3QLqtzfqmYVV+Hl58OmoDjz1zVb+vjCZ/OIynuwT51z11p6+Rn37lBuM\n+fZRC62yLWFVnCosZfxPe+jaKIzBrtCW90JbpkNYE6MSSVSJcyX2Ln+CuH5QmH3244TxZ2khlJ2B\nsmIjgZeXGhveurmf/dPDSNK1WxhN+n2CIbCucbc9uCEERoGHE65odTLeHu58cHdbXpi9g/eW7iOv\nqJyXbkpwruQeHgc3vQezx8DKt6HXX00J452f9pBfXM6rt7Rwru/f1eQcMG6c9nnZ6WvXzeRciT28\nifEhnJaHuxtvDU3E39uDz1YfpKCkjH/fnuhc+6kmDocDv8DKt4wpuuiudr38zszTfL0hjfuviyG+\njpOuEbicLV+BcoPWd5sdiVOTckdhd25uildubs64PnF8l5TBuBlbKC13jmX75w1627gHM/shY3rP\nTiwWzcvzdhLq58VTfV2nLQVgvNPeMg2aDpBp0GqSxC5MoZTimX5Nz/eXeei/SRSVOlHzMO8AGDYV\n8rNgwVNGWawdTFt/mM1pubw4OMFuu/HYze6FxtRqhwfNjsTpSWIXphrToxFvDm3Fyn3Z3PfZevKK\ny8wOqfLqtzf2o02ea4w0bexIbhFvLt5Nj7hwbmtrbkWOTSR9Ztzzatzb7EicniR2Ybo7Ozbkg7vb\nsjU9l3s+WUdOQYnZIVXedU8ay94XPwcn9tnsMlprXpq7E4uGf93WyrVumILxvTv0m9EwTxb7VZsk\nduEQbkqsx5T7OrDvWAF3fLzWLv00rMLNDW77GDx8jPYS5bb5pbRwexa/7D7Oszc2dc5tIK9m0xdG\n9Vrbe82OxCVIYhcOo1d8JP99sBPH8koYNnkNh04Umh1S5QTWgyEfwtHt8MvrVj99dn4Jr8zfRWJU\nEPdfF2P185uurNho+NXsJghwgc23HYAkduFQOjcKY8ZDXSgsKWf4x2vZY6cdZ6qt2WCjze/aD2H/\nL1Y7rdaa52dtp6CknHeGt7bLtmp2lzzXqCzq8IDZkbgMF/wpEc6uVVQQ3z3cFTcFd3y8lq3O0vb3\nxn9ARALMecTYsNwKZmxIZ9nu47wwoBlNnbWv/ZVoDesmGX2YYpxzs2pHJIldOKS42gHMfOQ6gnw9\nGfHJOtYccIK2v56+Rglk8WmY96dql0AePFHI3xcm071JuGtOwYDRZC9rG3R+xLhfIaxCvpPCYTUI\n9WPmI12pH+LL/Z9vZGnyMbNDurraLaD/P2HfT7B+cpVPU1xWwRMzNuPl4cY7w1vj5kwrc6/FuklG\niw8X2CXJkUhiFw4tMtCHb8d2JaFOAI9M28S8rZlmh3R1HcdA/CD4+WXI2l6lU7y+MJmdmXmMH96a\nOkE+Vg7QQZw6bCxKan8/eLnIzk8OQhK7cHghtbyY/lAXOsSE8NS3W5m+/rDZIV2ZUnDLh+AXZpRA\nll5bdc/szRl8vT6NR29oTN/mLlwlsmEKoKDTQ2ZH4nIksQun4O/twRcPdKJ3fCQvztnJ5BUHzA7p\nymqFGfXtOfvhxxcq/bSUrDxenLOTzrGhPNvPxXrBXKgkHzZ/Bc2HQFCU2dG4HEnswmn4eLoz+d72\n3NK6Hm8s3s2/FqVQYbFPj5YqaXQ9dH8KNv8Xds296uHH84sZ/cVGAn09+ODutq5Z2njO1q+h5LTR\niltYnQv/5AhX5Onuxrt3tmFU12imrEzl4a82UVjieBtMn9frRaOnzIJxkJt+2cOKSit46MskTp0p\nY+qojkQGuui8OhhdHFdPgAadoUFHs6NxSdVK7Eqpt5VSu5VS25VSc5RSwdYKTIjLcXdTvDakJa8P\nacHyPccZ+tEaMnOLzA7r0tw9Yeinxo5esx+Cij/+EqqwaJ75bivbM08z4e62tKwfZEKgdrTje8jL\ngB5/NjsSl1XdEfvPQEutdSKwFzBnOxlRI93XNYbP7+9IZm4RQz5cxabD9uuLfk1CG8FN/zE2R1/5\n9u++pLXmxTk7WLzzKC8OSqCfK98sBbBUwKp3oXYrYzc0YRPVSuxa65+01ueGIOsAuQsi7Kpn0wjm\n/Ok6anl7cNeUtXy55hDaTr3Rr0niHZB4l7Hr0uE1gJHUX1uQzDcb03midxPG9GhkcpB2sHsh5OyD\nHs/I1nc2ZM059geBxVY8nxCV0iQygHmPdaNnXASvzN/F4zO2UOCI8+6D3zF2XZr1EPrMKf75Qwpf\nrDnEg91iecaVK2DO0Rp+Gw+hjY1qGGEzV03sSqmlSqmdl/gYcsExLwLlwPQrnGesUipJKZWUnW2d\nPhpCnBPs58Un93XguQHxLN6RxS0frCL5SJ7ZYf3e2V2XdMFRtk++n09XpTKqa7TzbehdVft+MtoH\ndH9Keq7bmKru21al1CjgEaCP1vpMZZ7ToUMHnZSUVK3rCnE561JzGDdjC6fOlPJU36Y83LORw5QO\nFpaU8+PHzzP05KcsbfxX+ox8vmYkdYsFpvQ06tcfTzJuKotrppTapLXucLXjqlsVMwB4Hrilskld\nCFvr0iiMJU/15MYWdXh7yR6GTl7L/uMFZofF4ZxCbp+0hueybiAzrCt9D41HZW01Oyz7SJ4LR3fA\nDX+TpG4H1RqxK6X2A95AztmH1mmtH7na82TELuxlwbYjvDRvJ2dKKhjbsxF/6tUYPy8Pu8exZNdR\nnpu5HaXgg7vb0qOeG3zcw5iSeHgl+IbYPSa7qSiHSV2MHZIeXS3TMNVglxG71rqJ1rqB1rrN2Y+r\nJnUh7Onm1vX46eme3JRYlw+X76fv+BUs2HYEi51WrJ4uKuOZb7fy8FebiArxZf5j3ekRF2G0HBj+\nJeRlGf3bLRa7xGOK7d8YlTC9X5SkbieOMfEohA1FBvjwnzvb8P0jXQn09eSJGVsY/MEqftp11Gal\nkRUWzbcb0+gzfgXzth1hXJ845vypGw3DLtivtEFHo8Xv3h9h9Xs2icN0ZUXw6xtQr62x9Z2wC/u/\nJxXCJB1jQvlhXA/mb8vk/aX7GPvVJprW9uferjHc1rY+/t7V/+9QVmFh0Y4sPvr1ALuP5tOuYTCf\n39+RVlGXWU3aaSykr4dlf4eoDhDrYrsIrf0QTqfDrZOkbt2Oql0VUxUyxy7MVl5hYd7WI3y+5iA7\nM/Pw9/bgxua1GdCyDj2bRuDjWfkpA601+44XMG9rJnM2Z3LkdDFNIv15sk8cNyXWvXrVS0k+fNLb\n2Pfz4d8gsG41X52DyDsCH7SHJn3gzmlmR+MSKjvHLold1Ghaa7am5/L1+jR+Sj7G6aIyvNzdaFk/\nkHYNQ2gc6U/DUD/C/L3w8XBHA/nFZZwoKCE1u5DkI3msTc0h63Qxbgq6x0VwX5doejeLvLZdj47v\nhk/7GHt/PrDI2GbP2c1+GHbNgcfWQ2is2dG4hMomdpmKETWaUoq2DUNo2zCEf1VYWHsgh9X7T5B0\n+BT/XXuY0oor39QM9/emc2woXRuH0b9FHSICvKsWSGQzuH0KfDMC5j8Bt3/i3FMX6RuNm6bdn5Gk\nbgJJ7EKc5enuRs+mEfRsGgEY0zVH84pJO3mG3DNlFJdVABDg40loLU8ahfsTUsvLegE0Gwx9XoJf\nXofIBOjxrPXObU/lpUab4sD6Rk8YYXeS2IW4DA93N6JC/IgK8bv6wdbS/Rk4nmIk94hmRrJ3Nmve\nh+PJcPe3RhsFYXdS7iiEI1EKbvnA2Jxj1kNGbxVnkr0XVrwFLW6H+AFmR1NjSWIXwtF4+sKd043V\nqNOGwcmDZkdUOZYK4/6Apx8MfNPsaGo0SexCOKLAunDvbLCUwbTbocAJOqL+Nh7S18HAt8A/0uxo\najRJ7EI4qoh4uOc7o+3A9GFGvbujSt9grDBtNRxa32l2NDWeJHYhHFmDTnDHl0ZnxBl3Q2mh2RH9\n0ZmTMGsMBNWHwePNjkYgiV0Ix9e0P9z2MRxeDV/f6VjJvaIcZj4I+Vkw9DPwcfGNuJ2EJHYhnEHi\ncMdM7r+8BqnLYfB/jKZmwiFIYhfCWSTe8b/kPm2oMQVipo2fwpoJ0HEMtLvX3FjE70hiF8KZJN4B\nwz6DzE3w2QDITTMnjp2z4Ic/Q9OBMEBKGx2NJHYhnE2L2+DeOZB/FD7tC5mb7Xv95PlGg6+GXWH4\n5+AuC9gdjSR2IZxRTHcYvQTcveCz/rDpC7BHp9bNX8H3o4yNM+6e4RpdKF2QJHYhnFVkAoxdAdHd\nYMGTMPdRKMq1zbUsFbDsnzD/cWjUC+6bC77BtrmWqDZJ7EI4s1phMHIWXP88bP8WJnWFvT9Z9xoF\nx41KnJVvQZuRcPc34FXLutcQViWJXQhn5+YOvf4GY5YadeRfD4cZ90D2nuqd12KBbd/AxE5wcKVR\n0jjkQ/CwYqtiYROS2IVwFfXbw8MroM/LRiKe1MVYPJS+4drm3yvKIWUhTOkJcx6GsCbwyG/QcbRz\nb/5Rg8jtbCFciYe3sUFHu1Gw+j3Y9F+jNDG8qdHbvVEvqJtodI680JmTkLUV9v8CO2dD/hEIjobb\npkCrYca7AuE0ZM9TIVxZSYEx9548z1jYZCk3HvcOAr8Q0BajuVjRKeNxdy+IvR7aj4KmA8Dd07zY\nxR/Ydc9TpdSfgbeBCK31CWucUwhhBd7+xhRKx9FG8s7cDEe3Q94R43M3D6N/emisUWXTsKvcGHUB\n1U7sSqkGQD/ApCVwQohK8Q2BJn2MD+HSrHHz9F3gOcD+czpCCCH+oFqJXSl1C5Cptb7qxoxKqbFK\nqSSlVFJ2thPsBiOEEE7qqlMxSqmlQJ1LfOlF4G/AjZW5kNZ6CjAFjJun1xCjEEKIa3DVxK617nup\nx5VSrYBYYJsyalujgM1KqU5a66NWjVIIIUSlVfnmqdZ6B3B+x1ql1CGgg1TFCCGEuWTlqRBCuBir\nrTzVWsdY61xCCCGqTkbsQgjhYkxpKaCUygYOV/Hp4UBNm8eX11wzyGuuGarzmqO11hFXO8iUxF4d\nSqmkyvRKcCXymmsGec01gz1es0zFCCGEi5HELoQQLsYZE/sUswMwgbzmmkFec81g89fsdHPsQggh\nrswZR+xCCCGuwKkSu1JqgFJqj1Jqv1LqBbPjsTWlVAOl1HKlVIpSapdS6kmzY7IHpZS7UmqLUmqh\n2bHYg1IqWCk1Uym1++y/dVezY7I1pdTTZ3+mdyqlZiilfMyOydqUUp8ppY4rpXZe8FioUupnpdS+\ns3+GXOkcVeU0iV0p5Q5MBAYCzYG7lVLNzY3K5sqBZ7XWCUAX4LEa8JoBngRSzA7Cjt4HftRaNwNa\n4+KvXSlVHxiH0VuqJeAO3GVuVDbxBTDgosdeAH7RWscBv5z93OqcJrEDnYD9WutUrXUp8A0wxOSY\nbEprnaW13nz27/kY/+HrmxuVbSmlooDBwKdmx2IPSqlAoCcwFUBrXaq1zjU3KrvwAHyVUh6AH3DE\n5HisTmu9Ejh50cNDgC/P/v1L4FZbXNuZEnt9IP2CzzNw8SR3IaVUDNAWWG9uJDb3HsaOXBazA7GT\nRkA28PnZ6adPlVIuvemo1joTeAdjO80s4LTW+idzo7Kb2lrrLDAGblzQIdeanCmxq0s8ViNKepRS\n/sAs4CmtdZ7Z8diKUuom4LjWepPZsdiRB9AO+Ehr3RYoxEZvzx3F2XnlIRj7OdQDaimlRpoblWtx\npsSeATS44PMoXPDt28WUUp4YSX261nq22fHYWDfglrO9/b8Beiulppkbks1lABla63PvxGZiJHpX\n1hc4qLXO1lqXAbOB60yOyV6OKaXqApz987gtLuJMiX0jEKeUilVKeWHcbJlvckw2pYytqaYCKVrr\n/5gdj61prf+qtY462wL6LmCZ1tqlR3JndxtLV0rFn32oD5BsYkj2kAZ0UUr5nf0Z74OL3zC+wHxg\n1Nm/jwLm2eIiVuvHbmta63Kl1OPAEoy76J9prXeZHJatdQPuBXYopbaefexvWutFJsYkrO8JYPrZ\nAUsq8IDJ8diU1nq9UmomsBmj8msLLrgCVSk1A7gBCFdKZQCvAG8A3ymlRmP8ghtuk2vLylMhhHAt\nzjQVI4QQohIksQshhIuRxC6EEC5GErsQQrgYSexCCOFiJLELIYSLkcQuhBAuRhK7EEK4mP8HspZi\nAl23TMMAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(t, sol1[:, 0].imag, label=\"birds\")\n", "plt.plot(t, sol1[:, 1].imag, label=\"foxes\")\n", "plt.legend()" ] }, { "cell_type": "code", "execution_count": 30, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 30, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAD8CAYAAABjAo9vAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8VFX6x/HPyaSRQihJaCGkkJDQS6ihg9JFQUUEsSFg\nb+uq61p2XcuqWBAQQewoFlBcRJHeQUIvCb2FlkYgIaTO+f1xwZ+FkmTKnRme9+s1LyCZufc7Ck/u\nnHvOc5TWGiGEEJ7Dy+wAQggh7EsKuxBCeBgp7EII4WGksAshhIeRwi6EEB5GCrsQQngYKexCCOFh\npLALIYSHkcIuhBAextuMk4aGhuqoqCgzTi2EEG5rw4YNWVrrsCs9z5TCHhUVRUpKihmnFkIIt6WU\nOlSe58lQjBBCeBi7FXallEUptUkpNddexxRCCFFx9rxifxhItePxhBBCVIJdxtiVUhHAAOAl4DF7\nHFMIIUpKSkhPT6ewsNDsKE7l7+9PREQEPj4+lXq9vW6evg38HQi20/GEEIL09HSCg4OJiopCKWV2\nHKfQWpOdnU16ejrR0dGVOobNQzFKqYFAhtZ6wxWeN0YplaKUSsnMzLT1tEKIq0BhYSE1a9a8aoo6\ngFKKmjVr2vQpxR5j7MnAdUqpg8BMoKdS6vM/P0lrPVVrnaS1TgoLu+I0TCGEALiqivoFtr5nm4di\ntNZPA0+fD9Md+JvWeqStxxVC2JnWcDodTmyDM0fh3CnwsoBPINSIhvDGUK2+2SmFHZiyQEkI4SRa\nw6FVsON72PUTnEm//PNrxEB8P2g9CsITnJPRxR08eJCBAweyffv2P3x99OjRPPbYYzRu3Piyr1+6\ndClvvPEGc+c6bya4XQu71nopsNSexxRCVEJZCWz5Eta+Bxk7wbsKNOwFnR+BOi2hWiQE1ABthcIz\nkLMfjm2CfYvg16mwdhLE9IBez0K9Nma/G5f0wQcfXPTrZWVlWCwWJ6f5I1l5KoQn0RpS/weTO8AP\nDxpDLYMnw9/3wy0zoN09UL8tBNcCiw94+0FQGES2hw7jYMQ38Fgq9Hoejm+BaT3hu3uhIMfsd2aq\n0tJSRowYQWJiIjfeeCMFBQV07979t9YoQUFBPP7447Ro0YI1a9bw888/k5CQQOvWrZk9e/Zvx1m2\nbBktW7akZcuWtGrViry8PIfklaEYITzFmeMw91HY/ROENoLhMyG+L1T0RlxQGHR5DNqOhpVvwap3\njCv5IVMhprsjkpfLv/63g53Hztj1mI3rVuX5QU2u+Lxdu3Yxffp0kpOTueuuu5g8efIfvn/27Fna\nt2/P+PHjKSwsJC4ujsWLF9OwYUOGDRv22/PeeOMNJk2aRHJyMvn5+fj7+9v1/VwgV+xCeIJt38Lk\n9rB/CVz7Ety7Ghr1q3hR/z3/qtD7eRizBKpUh89ugJVvG58KrjL169cnOTkZgJEjR7Jy5co/fN9i\nsTB06FAA0tLSiI6OJi4uDqUUI0f+/1yS5ORkHnvsMSZMmEBubi7e3o65tpYrdiHcWWkxzP8HrJ8G\nEe3ghilQM9a+56jTAkYvgjn3w8LnIfcw9H8DvJx7XVieK2tH+fP0wz//2d/f/w/j6pearvjUU08x\nYMAA5s2bR3JyMvPnzychwf43qeWKXQh3lZ8JHw8winrHB+DOefYv6hf4BcFNH0Pyw5AyHWbfY9yg\nvUocPnyYNWvWAPDll1/SuXPnSz43ISGBAwcOsG/fvt+ef8G+ffto1qwZTz75JG3btiUtLc0heaWw\nC+GOsvfB9GuMOek3fgR9XjJuhjqSUnDNv6H3C7D9W/j+PrBaHXtOF9GoUSMmTZpEYmIiOTk53Hvv\nvZd8rr+/P1OnTmXAgAG0bt2a8PDw37739ttv07RpU1q0aIGPjw/9+vVzSF6lTRgvS0pK0rLRhhCV\nlL4BvrjJ+P3wr4xZLs62/A1Y/CK0vxf6vmLbWP5lpKamkpiY6JBju7qLvXel1AatddKVXitj7EK4\nk0NrYMaNEBgKI2c7bujlSro8bkyBXDsJQupBpwfNySEuSgq7EO7iQlEPrg13/Gj8ahal4Nr/GCtZ\nFzwHtZpAbE/z8og/kDF2IdyBKxX1C7y8jMVPYYnwzZ3G6lXhEqSwC+Hqjm+BGTe5VlG/wC/IWNEK\nRnEvLTY3jwCksAvh2nIOwOc3gn8IjPrBtYr6BTWiYfBEOL4ZlrxkdhqBFHYhXFd+Jnw+BKwlcNts\n4yalq0ocBK1vN9oPHFhudpqrnhR2IVxRUb4xpn7mONz6NYQ1MjvRlfV9xZil8/39Rn4PMWHCBBIT\nExkxYoTZUcpNCrsQrsZqhe/GwomtxmrP+u3MTlQ+voFw3UQ4fRiWvGx2GruZPHkyCxYsYMaMGWZH\nKTcp7EK4mqUvQ9pc6PMyNOprdpqKadARku6Gde8ZC6nc3Lhx49i/fz/9+vVj/PjxXH/99TRv3pwO\nHTqwdetWAB566CH+/e9/AzB//ny6du2K1WolMzOToUOH0rZtW9q2bcuqVasA57TulXnsQlyG1prM\nvCJOFZRQWFKGUhDk503NQD9CAhywhH/bt7D8dWMHo/bj7H98Z+j9vLFb0w8Pwthl9mt18NNTRgsF\ne6rdDPq9eslvT5kyhZ9//pklS5bwr3/9i1atWvH999+zePFiRo0axebNm3n11Vdp27YtXbp04aGH\nHmLevHl4eXnx8MMP8+ijj9K5c2cOHz5Mnz59SE1NdUrrXinsQvzO0dxzrNqTRcqhHLYcOc2hnLMU\nlly8H0rNQF8ahgfRProGHWNDaRtVHW+LDR+Cj24wOihGdoL+4x22TN/h/EOg/+vw1QhYP93YwMMD\nrFy5klmzZgHQs2dPsrOzOX36NCEhIUybNo2uXbvy1ltvERtrrAZeuHAhO3fu/O31Z86cIS8v77fW\nvSNGjGDIkCFERETYPasUdnHVy8wrYtbGdOZtO87W9NMAhFTxoXVkNbrEhRJZM4CagX74+3ihNeQV\nlZCVV8z+rHx2HjvDxCV7mbB4L6FBflzXoi4jOkQSGxZUsRB5J2HmCAgKh2GfgbevA96pEyUMMLbW\nW/oyNLsJAmvafszLXFk7w8X6al1oz7tt2zZq1qzJsWPHfvue1WplzZo1VKlS5Q+vcUbrXins4qq1\n+UguH606wLxtxykp07SICOHJvgn0TgwnNiwIL6/yXTGfKSxh1Z4s5mw+xudrD/HR6gNck1iLh3rF\n0bReyJUPUFYK394Jhafh7gVGHxh3p5QxS+a9ZKO4DxhvdiKbde3alRkzZvDss8+ydOlSQkNDqVq1\nKocOHWL8+PFs2rSJ/v37c/3119O+fXuuvfZaJk6cyBNPPAHA5s2badmy5W+te5s1a8b69etJS0uT\nwi6ErbYcyeWthbtZuiuTYD9vRrRvwMgODWgYXsGr7POq+vvQr1kd+jWrQ1Z+EZ+uPsgnaw4xaOJK\nbmwdwRN9GxEefJlx1MX/hkOr4IapULtpJd+VCwpPhLZ3w/oPIOkuo5+MG3vhhRe48847ad68OQEB\nAXzyySdorbn77rt54403qFu3LtOnT+eOO+5g/fr1TJgwgfvvv5/mzZtTWlpK165dmTJlCm+//TZL\nlizBYrHQuHFjh7Tulba94qpxNPccL/+Yyo/bjlMtwIexXWMZ1bEBgX72v745fa6ESUv28tGqAwT4\nevPi9U0Z1LzOX3fWSfsRZt5qFL6Bb9k9h+kKcmBCK4hoCyO/rfDLpW2vtO0V4qKKSsv4YMUBJi7e\ni0bzSO847u4cTbC/4zamCKniwz/6JzKsbX3+9s0WHvpyE/O3n+CVoc2oeuG8Ofvhu3uhTkvo84rD\nspgqoAZ0ftTYUu/QGmM6pHA4mccuPNrW9FwGTljJ6/N30b1RGAsf68YjveMdWtR/LzYsiG/GduSJ\nPo2Yv+ME109axd6MfCg5B1+PMsaib/4UfByzW71LaDcGgmoZG3NchRthm8Hmwq6Uqq+UWqKUSlVK\n7VBKPWyPYELYoqTMypsLdnPD5NXkFZby0Z1teW9kGyKqBzg9i7fFi/t7NOTz0e05XVDC9ZNWcWzm\nQ8ac7CFToXoDp2dyKt8A6PqEcR9h3+IKv9yM4WKz2fqe7XHFXgo8rrVOBDoA9yulGtvhuEJUyuHs\nAoZMXs2ERXsY3KIu8x/pSo9G4Vd+oYN1iKnJDw925rag9dTd9zW74++B+D5mx3KO1rdDSGSFr9r9\n/f3Jzs6+qoq71prs7GybFi7ZPMautT4OHD//+zylVCpQD9h52RcK4QDzd5zgb99swUsppoxsTd+m\ndcyO9Af19En+XjqFXT6NGbCtK8+tOciojlFmx3I8b1/o/hTMuQ92zTPmuZdDREQE6enpZGZmOjig\na/H397dp4ZJdb54qpaKAVsC6i3xvDDAGIDIy0p6nFYKSMiuv/ZzGtBUHaB4RwqRbW1O/hvOHXS6r\nrARmjUbhRdSYL+g+L5vn5uxAKcVtHTx8OAag+TBY9l9YMR4a9S/XylofHx+io6OdEM6z2O3mqVIq\nCJgFPKK1PvPn72utp2qtk7TWSWFhYfY6rRDknC1m5AfrmLbiAKM6NuCbcR1dr6gDLH0F0tfDoLfw\nC4tm8ojW9EoI59nvt/PthnSz0zmexRs6P2K0TjiwzOw0Hs0uhV0p5YNR1GdorWfb45hClMeek3lc\nP2kVm47k8tawFvx7cFP8vC1mx/qrA8thxZvQaiQ0HQqAr7cXk0a0pnPDUP7+7RZ+3n7C5JBO0OJW\nCKptXLULh7HHrBgFTAdStdZv2h5JiPJZtjuTIZNXU1BcxldjOnBDK/s3U7KLghyYPdbYhKLfa3/4\nlr+Phamj2tA8ohqPfLWJrem5JoV0Eh9/6PSA8YMuXRYpOoo9rtiTgduAnkqpzecf/e1wXCEu6dM1\nB7nzo1+JqBHAnAeSaRVZ3exIF6c1zHkAzmbC0OnGZhR/EuDrzbRRSdQM9GP0Jykcyz1nQlAnanMH\n+FczPsEIh7C5sGutV2qtlda6uda65fnHPHuEE+LPtNa8tWA3z83ZQc+EWnw7riP1qlW58gvNkjId\ndv0IvV+Aui0v+bSwYD8+urMt54rLuPuTFM4WlTototP5BRu95nf9CJm7zU7jkWTlqXAbVqvmxbmp\nvLNoDzcnRfD+bW0c0ufFbk7uhPnPQMPe0OG+Kz49vlYwE0e0Ju3EGZ75bptnz91uOxosfrBuitlJ\nPJIUduEWyqyaJ2dt5cNVB7grOZpXhzTHUs62uqYoOQez7ga/qnD9e+BVvn9q3eLDeKx3PN9vPsaM\ndYcdHNJEQWFGn/YtXxr3IIRdSWEXLq+otIwHv9zINxvSeaR3HM8OTCx3r3TT/PJPyNgJN7xnbJ5R\nAff3aEj3RmH8+387PftmaodxUFIAGz81O4nHkcIuXNq54jLu+XQD87ad4NmBjXmkd/xfW9+6mrQf\njR7kHR8whmEqyMtL8dbNLQkL9uO+GRvJKyxxQEgXULsZRHWBX6cai7eE3UhhFy7rTGEJoz5cx8o9\nmbw2tDl3d3aDFYhnjhn7ltZuDr2eq/Rhqgf6MmF4K47lnuOFHzy4O0eH++DMUUj9n9lJPIoUduGS\nsvOLGD51LZuP5PLu8Nbc3La+2ZGuzFoGs8dAaRHc+CF4+9l0uDYNqnN/j4bM2pjOT9uO2ymki4nv\nA9WjYe17ZifxKFLYhcs5fvocN7+/hn2Z+UwblcSA5q7VyOuSVr0DB1cYi5BC4+xyyId6xdGsXghP\nf7eNjDOFdjmmS/GyQPuxkP4rHNtkdhqPIYVduJSDWWe58b01ZJwp4tO72tPdBdrtlkv6BljyEjS5\nwWgbYCc+Fi/eGtaSwpIynprtoVMgWwwH7yqQ8pHZSTyGFHbhMtJOnOHGKWs4V1LGl2M60C66htmR\nyqcoz5jaGFwHBr5drq6FFdEwPIi/XduIxWkZzN3qgUMyVapBs6Gw7Vso/Ev/QFEJUtiFS9h0+BTD\n3l+Lt5fi67EdaFovxOxI5TfvCcg9BEOmGUXKAe7oFEXziBBe+GEHp84WO+QcpmpzF5SchW1fm53E\nI0hhF6ZbvTeLER+so1qAD9+M60jD8GCzI5Xf1m+MRTZd/+7QjZq9LV68OqQ5p8+V8NK8VIedxzT1\nWhsziVI+kn1R7UAKuzDVgp0nuePj9dSvHsA3Y120j/ql5ByAuY9C/Q7Gnp4O1rhuVcZ2i+HbDems\n3JPl8PM5lVKQdCec3C5dH+1ACrswzfebjjLu8w0k1qnKV2M7EF618ns8Ol1ZCcy+B5QXDJ1mbCLh\nBA/2jCM6NJBn52ynqLTMKed0mmY3gW8QpHxodhK3J4VdmOKztYd49OvNtI2qzozR7akW4Gt2pIpZ\n9t/zuyG9DdWct9Wjv4+FF65rwoGss0xfecBp53UKv2CjuO+YDedOmZ3Grblwa7yLKMqH/JNGb+vf\nHllQnG80XSopgJJCKCs25sd6+YCXt3E15R8CVar//yO4jvEPMqh2uRs0CfuYvHQvr/28i96J4Uy8\ntTX+Pi6449HlHFwJy984vxvSEKefvlt8GNc0rsXExXsZ0iqC2iFu9EnnStrcARs+gu2zjA6QolKU\nGfNik5KSdEpKJcbRfngINn7y1697+xsPnyrGw+JrrAK0lhq/lhUZ06hKL7KBgcUXQiIgNB5qNTn/\naAo1Gxo/HITdaK3578+7mLJsH9e1qMv4m1vgY3GzH6r5mTClM/gFwZhlxq8mOJJTQK83l9G3SW0m\nDG9lSgaH0Nr472vxhTFLzE7jcpRSG7TWSVd6nntdsbe8FRp0gsBQCAyDwHAIqAne5fwYX3IOzuXC\nuRw4c9yYopZ7GE4dhMxdsGcB6PPjln4hUL8dRHYwzlkvqfznEX9RZtU8/8N2Pl97mBHtI/n34Kau\n3Xb3Yqxlxrh6YS6MnGVaUQeoXyOAcd1imbBoD7e2j6RDTE3TstiVUsa/8/n/gIxUCE80O5Fbcq8r\ndkcrLTIK/MntcORXOLwGMtOM7/kGQUx3iLsGGl4DIfXMTOpWikutPPb1ZuZuPc7YbjE81TfB9Ts0\nXsyy12HJf2DQBGhzu9lpOFdcRu83lxHs783cBzvj7W6ffi7lbBaMbwQd7oVr/2N2GpfimVfsjubt\nB3WaG4+WtxpfK8gxCvzehcYVfdpc4+v12kDTG40l5FXdpJeJCc4VlzHu8w0s253JU/0SGNct1uxI\nlXNgOSx9GZrdDK1HmZ0GgCq+Fp4ZkMh9MzYya2M6w9o67yauQwWGQnxf2PIV9HoeLD5mJ7ILq1VT\natX4ejv+B7CH/Ih3oIAakDAABr4Fj2yD+9Yaf9nKSmD+0/BmInw80NgsoPis2WldyumCEkZOX8eK\nPZm8OqSZ+xb1/AyYNRpqxBp/D1zo00a/prVpFVmNNxfspqDYg/ZJbTkCzmYYF1QeYk9GPk1fmM+S\nXRkOP5cU9opQyhjz6/IYjFsB96+H7k9B3gn44UF4oxHMfQxObDM7qeky8goZNnUN29JPM/HW1tzS\nzk2vJq1lRlEvPAM3f2LquPrFKKV4pn8iJ88U8aEnTX+Mu8a4j7Z5htlJ7CblUA7FpVZiQgMdfi4p\n7LYIizcK+wPr4c6fjSv7TZ8bd/U/6A07vjcKw1XmcHYBN01Zw+GcAj68oy39m7nxUNWy1+DAMuj/\nujFjygUlRdWgT5NaTFm2n6z8IrPj2IfFB5oPg10/w9lss9PYxYZDpwgN8iXSCaurpbDbg1JGn5Ah\n78PjadDnFeMG0De3w7utYd3Uq2aYZsuRXIa8t4rcghJmjG5P57hQsyNV3q6fYdmr0OJWu7bidYS/\n903gXEkZExbtMTuK/bQcAdYS2PaN2UnsYuOhU7SOrO6UiQN2KexKqb5KqV1Kqb1KqafscUy3FVAD\nOt4HD26Amz8zpmT+9AS81QSWvgqFp81O6DALdp7klqlrqeJrYda9nWgVWd3sSJWXtdeY2linBQx8\n06XG1S8mNiyI4e3q88W6w+zPzDc7jn3Uamw0Btv6ldlJbJaZV8TB7AKSopzzb8Lmwq6UsgCTgH5A\nY2C4Uqqxrcd1e14WaHwdjF4Ad/0CkZ1g6SvwdjPj472H9Z3+dM1Bxn6WQnytIGbfm0zDcNcai66Q\nojyYeasxHDDsc2PRmxt4uFc8ft5evPHLLrOj2E/zm+HYRsjeZ3YSm2w8bLRIaNPATQo70A7Yq7Xe\nr7UuBmYCg+1w3L84V1yG1eqGLT0j28PwL2DscmiQbOy083YzWP660SbBjVmtmpd+3Mlzc3bQMyGc\nL8d0ICzYtr0+TWW1wnfjIHsv3PSxU/vA2Cos2I+7O0czb9sJdhzzkE+GTYcCyu2HYzYcOoWvxYsm\ndZ2zz4A9Cns94Mjv/px+/mt29+7iPXR5bQlvL9xN+qkCR5zCseq0gOFfwpilxorWxf+BCa2MbnZl\n7jdVLb+olLGfb2DaigPc3rEB79+WRICvmy+NWDneWKtw7YsQ3dXsNBV2d5cYqvp789aC3WZHsY+q\ndSGqM2z92q37tG84dIpmESFO64tkj8J+scHHv/wfUEqNUUqlKKVSMjMzK3WiNg2qEx0ayNsLjQJ/\n2/R1/G/LMQpL3GzmSd1WcOtXcPdCqBlr9PR+rxPs+slt/vIezDrLkMmrWJyWwfODGvPCdU3cr0XA\nn+38wfhh2+wm6HCf2WkqJaSKD2O6xrAwNYPNR3LNjmMfzW+GnH3GkIwbKiotY1v6aacNw4B9Cns6\nUP93f44Ajv35SVrrqVrrJK11UlhYWKVO1CuxFp+Pbs+Kv/fgoZ5x7M88y4NfbqL9y4t44YcdpJ1w\ns3Hr+m3hzp9g2AyjR82XtxiLnY669l/g5bszuW7iSjLyivjsrnbcmRztni0Cfu/oBpg9BiLawnXv\nuvzN0su5Izma6gE+vOkpV+2J1xlNwba653DM9qOnKS6z0tqJkwnsUdjXA3FKqWillC9wC/CDHY57\nSfVrBPDoNfGs+HsPPru7HV3jw/hi3WH6vr2CG99bzZzNR91nEwKlIHGgsaK1/xtGb5ppPWHOA0Yn\nQRditWreX7aPOz76lbrVqvDD/Z3p1NCNpzNekHsEvrgFgsLgli/d5mbppQT5eTOuWyzLd2ey/mCO\n2XFsV6UaxPcxWvm64ZDlhkPOvXEKdijsWutS4AFgPpAKfK213mHrccvDy0vRJS6Md4e3Yt0/evHP\nAYlk5Rfx8MzNdHplMa/9nMaRHDcZi7f4QLt74KFN0OkBYx/Nd9vA2iku8Zf51NliRn+awis/pdG3\naW1m3duJyJputI3dpRSegS9uNhrA3fqNUdw9wKiOUYQG+THeU2bINLvJaDFwYJnZSSpsw6FTNKgZ\n4NRJBXaZx661nqe1jtdax2qtX7LHMSuqeqAvo7vEsPjx7nx6VztaN6jOlGX76Pr6Eu7+eD1L0jIo\nc4cZNf5VjY529642Nvj9+Ul4vwscWGFapJSDOfSfsIKVe7L413VNmHRrawL93PwmKUBpMXxzh9HR\n8+ZPIDzB7ER2U8XXwv09Ylm7P4fVez1gf9S4PkYrbTebHaO1ZsOhU069WgcPXHnq5aXoGh/GtFFJ\nrHyyJw/0aMiW9NPc+fF6ur+xhPeW7iPbHZZdhzWC274z5lEX58MnA40idDrdaRFKy6y8u2gPw6au\nxdfbi9n3deL2TlHuP54OxrTG7++FfYuM7e1ie5idyO6Gt4ukTog/by30gLF2H39oPAhS/wfFbvIp\nHDicU0BWfrEUdnuqW60Kj1/biNVP9WTira2oG1KF//6cRsdXFvPIzE2kHMzBjH705aYUJA6C+3+F\n7k8bs2YmtoUV442hAwfam5HH0PdWM37Bbvo3q8PcBzvTtJ5z5uA6nNbGauDt30LvF1ymDa+9+ftY\nGNctlvUHT7Fuvwf0W2l2k3GRs3eB2UnKLeWg88fXwcML+wW+3l4MbF6Xr8Z2ZMGjXbm1fSSLUjO4\nccoa+r2zgs/XHiK/yPxx7EvyqWI0G7v/V4jtCYv+bUyP3LvI7qcqs2o+WLGf/hNWcjingIm3tuLd\n4a0I9veMntiAsQJ4/QfQ6SHo/KjZaRxqWNv6hAb5MXHJXrOj2K5BZwgIhR3fmZ2k3DYcPkWwnzfx\n4cFOPe9VUdh/L65WMC9c14R1z/TilSHNsHgp/vn9dtq/tJB/fr/NtadMVm8At8yAEbNAW+HzIfDV\nSGNWhx1sPpLL9ZNW8Z8fU+kWH8Yvj3ZjYPO6djm2y1g9EZb912jqdc2/zU7jcP4+Fu7pEs2KPVls\nOr+s3W1ZvI02Hbvnu81wzIaDp2jVoDpeTl7jcdUV9gsCfL0Z3i6SuQ925rv7OtG3aR2+SUn/bcrk\n95tceMpkXG9jemTPZ2HPQmN4ZvkblR6eOXW2mKdnb+OGyas4eaaQCcNbMfW2Nu7dGuBiVk2AX56B\nxoNh4DtuPVe9IkZ0aEC1AB8mecJVe5MboKQA9vxidpIryi0oZtfJPNo6eRgGruLCfoFSilaR1Rl/\nc4vfpkxmny3mka820/GVxbwyL5U9J/PMjvlX3n7Q9W/wwK9GoV/8IkzuaBT6ciosKWPa8v30GL+U\nr1OOcFdyNIse78Z1Lep6xg3S31v5Fix4FpoMgaHTjau/q0SQnzd3JUezMDXD/XvINEg2NuBwg+GY\n9efH19ubsNH4VV/Yf69agDFlctFj3fj87va0i6rB9JUHuOat5Vw3cSWfrD5Iztlis2P+UbVIY+bM\nyNnGFeiMoTBzBJw6dMmXlJRZ+Xr9EXq8sZSX5qXSPKIaPz7UmWcHNvassXQwbpQufx0WvmDcfBsy\nzWP20KyI2ztFEeznzeQl7t0l0eiaOvj8cIxr73Hw64FsfL29aB7h/EkHyoxZIUlJSTolJcXp562M\nrPwi5mw+xqwN6ew8fgYfi6JHo3CGtI6ge6MwpzX1KZfSIlgzyShk2gpdHjduEPr4A1BQXMpX64/w\nwYoDHM09R4v61XiybyM6xXrA6tGLsVph/j9g3XvGbjzXv2cUhqvU6/PTmLx0Hwse7UpDJ9/Ms6uD\nK+HjAXDjR9B0iNlpLum6iSvx97Hw9diOdjumUmqD1jrpis+Twl5+qcfPMHtjOt9tOkZWfhEBvhZ6\nNAqnb9NH17rUAAAcwUlEQVTa9EgIJ8hVFu2cTof5z8DO76F6NCeS/8XHmY2Yuf4wuQUltI2qztiu\nsfRKDPe8IZcLSovgu7HGR/YO98G1L4HX1f0BNTu/iM7/XUK/prV5c1hLs+NUnrXM2ES+fnsY9pnZ\naS4qv6iU5i/M5/4eDXn82kZ2O255C7uLVCL3kFinKs8MaMyTfRNYsz+bn7af4JcdJ/hx23F8vb3o\nGhdKt0bhdI0LpUFNx29Ye0khEZwe+AFbqg4kdv0L1Js7ijbWNpxu8Cg39k6mTYMa5mVzhoIc+HoU\nHFwB17wInR68am6UXk7NID9GtI/ko9UHeaR3vPu2hLgwHLPxU2M/AxfbYByMNgJWDe2jnT++DlLY\nK8Xb4kWXuDC6xIXx4uCmpBzM4aftJ1iw8yQLUzMAqF+jCp0bhtEptiYt61cjonoVh14da605kHWW\nlXuz+GXHSdbuz6bUGkBsjbd4sdYKeqV/wDXHR8OBR6Huw27f6OqSTmyHmcMh7wTcMBVaDDM7kUu5\np2sMn645xLQV+3nx+qZmx6m8JjfAr1Nh98/Q7Eaz0/zFrwey8fZStG5QzZTzy1CMHV0oriv2ZLFi\nTxZr9mVxttiYMhka5EfL+tVoWq8qsWFBxIYFERMWWKkxeq01J88UsetkHrtP5LH5SC7rDuT8tkN9\nTGgg1zapzbVNatEyopoxh/b0Ufjln7BjNlRrAL2eM2aIeNLwxI7v4Pv7wD/EaIUc0cbsRC7pyW+3\n8v3mo6x6qiehQW46pdVqNYZjIpKMtR0u5sb3VlNq1Xx/f7JdjytDMSZQShETFkRMWBC3d4qipMxK\n2vE8Nh85xaYjuWw+ksuitJO/7aWhFIQF+REWbDxCg/wI8LXg5+2Fn7cFjaawxEpRaRn5haVk5BVx\n8kwhJ88U/WGlbN0Qf7rEhdIuugbto2sQE3aRj6Yh9eCmj6DNHfDzUzDrblj1DvR+HmJ7ufdQRfFZ\n4ybpho8hop0x7hpc2+xULuuerjF8veEIn64+yGN2HP91Ki8vaHI9pHxk7FHr5zo3gwtLytiSnstd\nydGmZZDC7kA+Fi+aRYTQLCKE287fGC8sKeNg9ln2ZuSzNyOf47mFZOYXkZlXxK4TeZwrKaO41Eph\nSRleSuHn7YW/j4UqvhZqVfUnvlYwXeLCiAkLJL5WMPG1gqkR6Fv+UDHdYNxKo0vekpfg86EQ1cXo\nmRJxxQsB13N8C3x7t7FHafLD0OOf4F2B/x5XoYbhQVyTWItP1hxibLdY9+3U2eQGWDcFdv0MzW8y\nO81vNh3OpaRM0z7GvHtZbvp/1H35+1hIqF2VhNpVzQvhZYEWtxj/MDZ8DMtegw96Ga1Ruzxm7Mfq\n6ooLYPlrsPpdY8HKqDnGDy1RLmO7xfLLzpN8tf4Id3U278rSJhHtIKg2pP7gUoX91wM5KIWpkxQ8\naIBVVJi3H7QfCw9vhp7/hKMp8GEf+Kg/7F3omvuvam1coU3uYKwmbT7M6F0vRb1C2jSo/tsCvJIy\nq9lxKsfLy9h9bO9Cl+ods+5ANom1qxJSxbyFcFLYhTE+2fUJeGQb9H0VTh00hmimdIb1040xTFdw\nZL2xMOXLYcYemLfPhesnQ4CHT990kLHdYjiae465W/+yRbH7SBxk9I7ZZ/9Op5VRXGpl4+FTpg7D\ngBR28Xu+gdDhXnhoMwyeBMoLfnwMxifCj48b49nOvorXGvYtNn7QTO8NWbuNvWHvXQ3RXZybxcP0\naBROfK0g3l+237X3JbicBslQpbqxAYcL2HY0l8ISK+2izC3sMsYu/srb12hr23IEHN1g9C7f+Jnx\na82G0HSoMVUyrJHjZtPkZxobYWz8DDJ2QGC40c2y/TiXXJDijry8FGO6xvK3b7awdFcmPRLCzY5U\ncRYfaNQfUucaWx2afON89d5slIIOJjT++j0p7OLSlDJmykQkQZ+XjZtU22cZvWiW/ddoQBbTw9hW\nrkEyBNlQGLQ2rsb3LYY9C2D/UtBlULu58emh2U3GPQFhV9e1qMv4X3YxZdk+9yzsYAzHbJ4BB5dD\nw96mRlm1L4vGdapSvSIz1RxACrson4Aaxhz4NndA3klI+x/sW2IsCtr4ifGc4DpQpwWEJUC1+sZC\nqMAw8A0C3wCjeJcWQsk5Y2XomaOQsx9ObIMTW+FspnGcmg0h+SHjxmh4olnv+Krg6+3F3Z2j+c+P\nqWw6fIpWkc7vHW6zmB7gE2gMx5hY2M8Vl7HxUC53JEeZluECKeyi4oJrQdvRxqOsFI5thPQUOL4Z\njm02tuyzlpTvWF4+EJ5gTLWMSDK2/qvewLH5xR/c0i6SCYv28P6y/Uy5zQ1X6/r4Q/y1kPYjDHjT\ntA6eKYdyKC6z0jHW3GEYkMIubGXxhvrtjMcF1jLIPwm5h6Eg21gZWpwPymIMp3j7G1f3VesaK0Sv\nwv7oriTIz5tRHaOYtHQv+zPzL75y2dUlDjI+PR5ZBw06mRJh1V6jP4zZN05BCrtwBC+LUbSreth+\nqR7s9k5RTF2xn2krDvDKkGZmx6m4uGvB4mcMx5hU2Ffvy6JVZDWXWMlr03RHpdTrSqk0pdRWpdR3\nSilzWpkJIWwSFuzHjW0imLUxnYy8QrPjVJxfsDGMl/o/UxbWnS4oYfvR0y6zaY2t89gXAE211s2B\n3cDTtkcSQpjhni4xlJRZ+WT1QbOjVE7iIDh9xLjX42RrD2Rj1ZDc0AMKu9b6F631hTaDa4EI2yMJ\nIcwQHRpIn8a1+WzNIc7+rnuo22jUz7iPY8JipdV7s6jiY6FlfdcYtLDnytO7gJ8u9U2l1BilVIpS\nKiUzM9OOpxVC2MvYbjGcKSxl5vojZkepuIAaENXZlMK+al82baNr4OvtGov5r5hCKbVQKbX9Io/B\nv3vOM0ApcMmO91rrqVrrJK11UlhYmH3SCyHsqlVkddpF12D6iv3u2RwscZCx0C1zl9NOefJMIXsz\n8kl2gWmOF1yxsGute2utm17kMQdAKXU7MBAYod224YQQ4oKxXWM4drrQPZuDJQw0fk39wWmnXLU3\nC3Cd8XWwfVZMX+BJ4Dqttev0zRRCVFqPRuHEhbtpc7CqdYw+7TudV9iX7sokNMiXxnVM3GPhT2wd\nEJoIBAMLlFKblVJT7JBJCGEiozlYDGkn8li+J8vsOBWXONBoUXHqkMNPVWbVrNiTSde4MGNvYRdh\n66yYhlrr+lrrlucf4+wVTAhhnsEt61Grqh/vL9tndpSKuzAck/ajw0+17ehpThWU0K2Ra903dI1b\nuEIIl+Lr7cVdydGs3pfN1vRcs+NUTM1YCG8MaXMdfqpluzJRCrrESWEXQriB4e0jCfbz5v3l+82O\nUnEJA+HwGqOvvwMt251B84hqFdtQ3gmksAshLqqqvw+3dojkp23HOZztZnMjEgeBtsKueQ47RW5B\nMZuP5NIt3rWu1kEKuxDiMu5KjsbipfhgpZtdtdduZmwE48DhmJV7s7BqpLALIdxLrar+XN+yHl+n\nHCE7v8jsOOWnFCQMMnbiKjzjkFMs25VJSBUfWkSEOOT4tpDCLoS4rDFdYygssfLpGsdPH7SrxIFQ\nVgx7F9j90Fprlu3OpHNcKN4W1yujrpdICOFS4moF0zsxnE/XHORccZnZccqvfnsICDU2uraz7UfP\nkJFXRHcXHIYBKexCiHIY2y2WUwUlfLPBjZqDeVkgob+xOXqpfYeRFqSexEtBTxfdAFwKuxDiipIa\nVKdVZDWmrdhPqTs1B0sYBMV5sH+ZXQ+7cOdJ2jSoTs0gP7se116ksAshrkgpxdiusRzJOcdP20+Y\nHaf8YrqBbzCk2a+Vb/qpAnYeP0PvxFp2O6a9SWEXQpTLNY1rERMayNTlbtQczNsP4q+FtHnGJut2\nsCg1AzD+e7gqKexCiHKxeCnu6RrDtqOnWbMv2+w45ZcwEAqy4PBauxxuYepJYsICiQkLssvxHEEK\nuxCi3G5oVY/QID+muFObgbhrwOJnl8VKZwpLWLs/m2tceBgGpLALISrA38fCnclRLN+dyc5jjln4\nY3d+wRDT3Zj2aOMQ0rJdmZSUaZcehgEp7EKIChrZvgEBvhamLnejlr6JA+H0YaNPuw0Wpp6kRqAv\nrSKr2ymYY0hhF0JUSEiAD8PbRfK/rcdJP+UmzcEa9QflZdNipcKSMhalZtA7MRyLC22qcTFS2IUQ\nFXZX52gU8OHKg2ZHKZ/AUIjsaNM4+4o9WeQXldK/WR07BnMMKexCiAqrV60Kg1rUZeb6w+QWFJsd\np3wSBkLGTsiu3BDSvG3HCani41KbVl+KFHYhRKWM6RpDQXEZn6x2k+ZgCQOMXytx1V5YUsaCnSfp\n06QWPi7Y9OvPXD+hEMIlJdapSu/EcD5cdYC8whKz41xZ9QZQp0WlxtndaRgGpLALIWzwUK84Tp8r\ncZ+WvgmDIP1XOHO8Qi9zp2EYkMIuhLBB84hq9GgUxrQV+8kvKjU7zpUlDjR+3fVjuV/ibsMwIIVd\nCGGjh3rFkVtQwmfucNUelgA1Yis0HLN0Vwb5RaUMaF7XgcHsyy6FXSn1N6WUVkq5x+cUIYTdtIqs\nTtd446q9oNjFr9qVMq7aD66Ac6fK9ZJZG48SFuxHcmxNB4ezH5sLu1KqPnANcNj2OEIId/Rwrzhy\nzhbz+Vo3uGpPGATWUtj9yxWfmnO2mKW7Mri+ZV2X3ALvUuyR9C3g74Cb9PEUQthbmwbV6dwwlKnL\n97v+9nn12kBQ7XL1aJ+79RglZZohrSOcEMx+bCrsSqnrgKNa6y12yiOEcFMP944jK7+YGetc/Krd\ny8uY0753EZScu+xTZ208SkLtYBLrVHVSOPu4YmFXSi1USm2/yGMw8AzwXHlOpJQao5RKUUqlZGZm\n2ppbCOFi2kbVoFNsTd5bus/1Z8gkDoSSAti3+JJP2ZuRz5YjuQx1s6t1KEdh11r31lo3/fMD2A9E\nA1uUUgeBCGCjUqr2JY4zVWudpLVOCgtzzZ29hRC2eaJPI7LPFvPhygNmR7m8qC7gH3LZ2TGzN6bj\npWBwS/eZDXNBpYditNbbtNbhWusorXUUkA601lq70YaIQgh7ahVZnT5NajF1+X5yzrpwDxmLD8T3\ng13zoOyvq2ZLyqx8syGdbvFhhFf1NyGgbdznNq8Qwi387dpGFBSXMnnJXrOjXF7iQCjMhUOr/vKt\nRaknycwrYkT7BiYEs53dCvv5K/csex1PCOGe4moFM7R1BJ+uPcSx3MvfnDRVbC/wrnLR4ZgZ6w5T\nN8SfHgnhJgSznVyxCyHs7pFr4kHDOwv3mB3l0nwDoGEvSPsRrNbfvnwo+ywr9mQxrG2ky2+ocSlS\n2IUQdlevWhVu69iAbzYcYW9GntlxLi1hIOQdg2ObfvvSl78eweKlGNa2vonBbCOFXQjhEPf3aEig\nrzf/+THV7CiXFt8HlOW3xUrnisv4av1heieGUzvE/W6aXiCFXQjhEDUCfXm4dxxLd2WyJC3D7DgX\nF1ADojr/Ns7+7cZ0ThWUMLpLjMnBbCOFXQjhMKM6RhETGsiLP+6kuNR65ReYIXEQZO/BejKND1ce\noEVECEkNqpudyiZS2IUQDuPr7cU/ByayP/Msn645aHacizu/Zd7+FTM5kHWW0V1iUMo9b5peIIVd\nCOFQPRqF0y0+jHcW7SErv8jsOH9VtS7US0KlzaVetSr0a3rRxfNuRQq7EMKhlFI8OzCRc8VlvDIv\nzew4F3UovCexpXt4KMnfrdrzXor7vwMhhMtrGB7M2G4xzNqYzqq9rrWOUWvNm+nxAAwJ2GxyGvuQ\nwi6EcIoHe8YRVTOAZ77bRmGJ6/RsX7MvmzlHAjgVGIPP7nlmx7ELKexCCKfw97Hw0g3NOJhdwLuL\nXWNFqtaatxbuplZVP4Jb3mD0jTmbbXYsm0lhF0I4TXLDUIa2juD9ZfvZfvS02XFYmJrB+oOnuL9H\nQ7ybDAJthd0/mR3LZlLYhRBO9ezARGoG+fLIV5tNHZIpLrXy8rxUYsMCGd4uEuq0hJD6l+3R7i6k\nsAshnKpagC9v3NSCvRn5vPqTebNkPl1zkANZZ/nngMb4WLxAKWNO+77FUJRvWi57kMIuhHC6LnFh\n3JkcxcerD7J0l/PbDZw4Xcg7C/fQJS6U7o1+t6Nb4iAoK4K9C5yeyZ6ksAshTPFk3wQa1Qrm0a82\ncySnwGnn1Vrz7JztFJdZeXFw0z+uMo3sCAE13X44Rgq7EMIU/j4WptzWhlKr5t4ZG5w23j5v2wkW\n7DzJ49fGExUa+MdvelmMVr67f4Zi5/2wsTcp7EII00SHBvL2sJZsP3qGf3y3Da21Q893NPcc//hu\nG80jQrgrOfriT2o6FIrzYc98h2ZxJCnsQghT9UqsxSO945i98ShvO3DHpZIyKw9+sZEyq2bCLa0u\n3TogqjME1YLtsxyWxdG8zQ4ghBAP94rjWO453lm0h9BgP27rYN9NpLXW/GfuTjYezmXira3+OgTz\ne14WaHIDpHwEhWfAv6pdsziDXLELIUynlOLlG5rRKyGc5+ZsZ8a6Q3Y9/rQV+/lkzSHu6RLNwOZ1\nr/yCpkON2TG73LPFgBR2IYRL8LZ4MWlEa7rHh/HMd9uZsmyfXY4789fDvDwvjQHN6/B0v8TyvSii\nLYREuu1wjBR2IYTL8Pex8P5tSQxqUZdXf0rjqVlbKz1bRmvN1OX7eGr2NrrFhzH+phZ4eZVzAw2l\noOkQY7FSQU6lzm8mKexCCJfi6+3F28Nacn+PWGauP8LQ91az41jF+srkF5Xy8MzNxpV6szpMG5WE\nv4+lYkGaDgVrKeycU7HXuQCbC7tS6kGl1C6l1A6l1Gv2CCWEuLpZvBRP9Engg1FJnDhdyKB3V/L8\nnO0cP33usq8rLbPy3aZ0eo1fytytx3iiTyPeHd4KX+9KlLrazaBmnFsOx9g0K0Yp1QMYDDTXWhcp\npcLtE0sIIaB341osjurOa/PT+HzdYb749TDd4sPomVCLxnWrUiPAlxKrlcM5Bazdn83cLcc5mnuO\nJnWrMmVkG1pF2rAptVLGVfuy/0LeCQh2ny3zlC0LApRSXwNTtdYLK/K6pKQknZKSUunzCiGuPkdy\nCvhk9UHmbTvOsdOFf/m+j0XRPromozo2oFdiLSzlHU+/nMzdMKkt9P0vdBhn+/FspJTaoLVOuuLz\nbCzsm4E5QF+gEPib1nr9JZ47BhgDEBkZ2ebQIftOZxJCXB201hzKLmDXyTzyCkvxUlC3WhWa1gsh\nyM8BS3OmdAYvHxizxP7HrqDyFvYr/ldQSi0ELvYZ5Jnzr68OdADaAl8rpWL0RX5aaK2nAlPBuGK/\n0nmFEOJilFJEhQZefpGRPbUYDvP/YVy9h8U755w2uuIdBa11b61104s85gDpwGxt+BWwAqGODi2E\nEE7T7CZQFtjyhdlJys3WWTHfAz0BlFLxgC/gWluQCyGELYLCoWFv2PIVWF1nE+7LsbWwfwjEKKW2\nAzOB2y82DCOEEG6t5XDIOwYHlpmdpFxsutOgtS4GRtopixBCuKb4fuAfAltmQmxPs9Nckaw8FUKI\nK/HxhyZDIPV/UJRndporksIuhBDl0fJWKClwixYDUtiFEKI8ItpCjVjY/KXZSa5ICrsQQpSHUsac\n9kMrIWe/2WkuSwq7EEKUV6sRxpz2DZ+YneSypLALIUR5Va0L8X1h8wwoLTY7zSVJYRdCiIpIuhPO\nZkLaXLOTXJIUdiGEqIjYnsa2eRs+MjvJJUlhF0KIivCyQJtRcGA5ZNtnX1Z7k8IuhBAV1XLk+Zuo\nH5ud5KKksAshREVVrQON+p2/iVpkdpq/kMIuhBCV0XY0FGTDtm/NTvIXUtiFEKIyYrpDeGNYOxlc\nrKmtFHYhhKgMpaDDfXByu8u185XCLoQQldXsJggMgzWTzU7yB1LYhRCisnz8jbH2PfONPVFdhBR2\nIYSwRdLdYPEzxtpdhBR2IYSwRVAYtLgFNn8BeSfMTgNIYRdCCNt1fgSspbBqgtlJACnsQghhuxox\n0HwYpHwI+Rlmp5HCLoQQdtH1b1BWBKvNv2qXwi6EEPZQM9aY/rh+OuRnmhpFCrsQQthL1yeM3jHL\nXzc1hk2FXSnVUim1Vim1WSmVopRqZ69gQgjhdkLjoM0dkDIdsvaaFsPWK/bXgH9prVsCz53/sxBC\nXL26Pw3eVWDh86ZFsLWwa6Dq+d+HAMdsPJ4QQri3oDBj+mPaXDi40pQIthb2R4DXlVJHgDeAp22P\nJIQQbq7Dfcb2eXMfNaVf+xULu1JqoVJq+0Ueg4F7gUe11vWBR4HplznOmPPj8CmZmebeMRZCCIfy\nDYCBb0HWbljxptNPr7QNfYSVUqeBalprrZRSwGmtddUrvS4pKUmnpKRU+rxCCOEWZo2GHd/DuJUQ\nnmDz4ZRSG7TWSVd6nq1DMceAbud/3xPYY+PxhBDCc/R5BfyC4bsxTh2SsbWw3wOMV0ptAV4Gxtge\nSQghPERQGAyeBMe3wILnnHZab1terLVeCbSxUxYhhPA8Cf2h/b2w7j2I7ABNbnD4KWXlqRBCONo1\n/4KIdvDdvXB0o8NPJ4VdCCEczdsPbvnCuGL3C3b86Rx+BiGEEMZ4+6jvnXIquWIXQggPI4VdCCE8\njBR2IYTwMFLYhRDCw0hhF0IIDyOFXQghPIwUdiGE8DBS2IUQwsPY1La30idVKhM4VMmXhwJZdozj\nDuQ9Xx3kPV8dbHnPDbTWYVd6kimF3RZKqZTy9CP2JPKerw7ynq8OznjPMhQjhBAeRgq7EEJ4GHcs\n7FPNDmACec9XB3nPVweHv2e3G2MXQghxee54xS6EEOIy3KqwK6X6KqV2KaX2KqWeMjuPoyml6iul\nliilUpVSO5RSD5udyRmUUhal1Cal1FyzsziDUqqaUupbpVTa+f/XHc3O5GhKqUfP/53erpT6Uinl\nb3Yme1NKfaiUylBKbf/d12oopRYopfac/7W6I87tNoVdKWUBJgH9gMbAcKVUY3NTOVwp8LjWOhHo\nANx/FbxngIeBVLNDONE7wM9a6wSgBR7+3pVS9YCHgCStdVPAAtxibiqH+Bjo+6evPQUs0lrHAYvO\n/9nu3KawA+2AvVrr/VrrYmAmMNjkTA6ltT6utd54/vd5GP/g65mbyrGUUhHAAOADs7M4g1KqKtAV\nmA6gtS7WWueam8opvIEqSilvIAA4ZnIeu9NaLwdy/vTlwcAn53//CXC9I87tToW9HnDkd39Ox8OL\n3O8ppaKAVsA6c5M43NvA3wGr2UGcJAbIBD46P/z0gVIq0OxQjqS1Pgq8ARwGjgOntda/mJvKaWpp\nrY+DceEGhDviJO5U2NVFvnZVTOlRSgUBs4BHtNZnzM7jKEqpgUCG1nqD2VmcyBtoDbyntW4FnMVB\nH89dxflx5cFANFAXCFRKjTQ3lWdxp8KeDtT/3Z8j8MCPb3+mlPLBKOoztNazzc7jYMnAdUqpgxhD\nbT2VUp+bG8nh0oF0rfWFT2LfYhR6T9YbOKC1ztRalwCzgU4mZ3KWk0qpOgDnf81wxEncqbCvB+KU\nUtFKKV+Mmy0/mJzJoZRSCmPsNVVr/abZeRxNa/201jpCax2F8f93sdbao6/ktNYngCNKqUbnv9QL\n2GliJGc4DHRQSgWc/zveCw+/Yfw7PwC3n//97cAcR5zE2xEHdQStdalS6gFgPsZd9A+11jtMjuVo\nycBtwDal1ObzX/uH1nqeiZmE/T0IzDh/wbIfuNPkPA6ltV6nlPoW2Igx82sTHrgCVSn1JdAdCFVK\npQPPA68CXyul7sb4AXeTQ84tK0+FEMKzuNNQjBBCiHKQwi6EEB5GCrsQQngYKexCCOFhpLALIYSH\nkcIuhBAeRgq7EEJ4GCnsQgjhYf4PvqT3GbXMdSsAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(t, sol2[:, 0].imag, label=\"birds\")\n", "plt.plot(t, sol2[:, 1].imag, label=\"foxes\")\n", "plt.legend()" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.5.2+" }, "widgets": { "state": {}, "version": "1.1.2" } }, "nbformat": 4, "nbformat_minor": 0 }