{ "metadata": { "name": "", "signature": "sha256:1a4ed5ca21c34718c5cd2b9e011dbd4e678c2e641c90c164f101adebf604e4bf" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "\u9019\u500bnote\u662f\u5728\u7c21\u4ecb\u5982\u4f55\u4ee5scipy\u89e3\u7c21\u55ae\u7684\u5fae\u5206\u65b9\u7a0b\u3002\u8003\u616e\u7684\u65b9\u7a0b\u5f0f\u662fLotka\u2013Volterra equations\uff1a\n", "\\begin{eqnarray}\n", "\\frac{d}{dt} x &=& x(\\alpha -\\beta y)\\\\\n", "\\frac{d}{dt} y &=& y(\\gamma -\\delta y)\\\\\n", "\\end{eqnarray}\n", "\u6211\u5011\u4e3b\u8981\u662f\u4ee5odeint\u5c0d\u4e00\u968e\u5fae\u5206\u65b9\u7a0b\u7cfb\u7d71\u7a4d\u5206\uff0c\u5728\u9032\u5165\u6b63\u984c\u4e4b\u524d\u6211\u5011\u5148\u5f15\u5165\u9700\u8981\u7684\u5e7e\u500b\u6a21\u7d44\uff0c\u5176\u4e2dmatplotlib\u662f\u5c07\u6700\u5f8c\u5f97\u5230\u7684\u7d50\u679c\u756b\u51fa\u4f86\u7684\u6a21\u7d44\u6700\u5f8c\u5f97\u5230\u7684\u7d50\u679c\u756b\u51fa\u4f86\u7684\u6a21\u7d44\uff0c\u540c\u6642 \"%matplotlib inline\" \u53ef\u4ee5\u5c07\u5716\u5d4c\u5165ipython notebbok\u4e2d\uff1bnumpy\u5247\u662f\u6578\u503c\u8a08\u7b97\u9700\u8981\u7684\u6a21\u7d44\u3002" ] }, { "cell_type": "code", "collapsed": false, "input": [ "from scipy.integrate import odeint\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "%matplotlib inline" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "markdown", "metadata": {}, "source": [ "\u5e95\u4e0b\u662f\u4e3b\u8981\u7a0b\u5f0f\uff0c\u4e00\u958b\u59cb\u6211\u5011\u5148\u8ce6\u4e88\u53c3\u6578\u503c\uff0c\u63a5\u8457\u5b9a\u7fa9\u4e00\u968e\u5fae\u5206\u65b9\u7a0b\u7cfb\u7d71\n", "$$\\frac{d}{dt} u =f(u,t)$$\n", "\u5b9a\u7fa9\u65b9\u5f0f\u662f\u4ee5 def \u5b9a\u7fa9\u51fd\u6578 $f(u,t)$\uff0c\u6ce8\u610f\u51fd\u6578\u50b3\u56de\u503c\u7684\u90e8\u4efd\"return [dxdt,dydt]\"\n", "\n", "\u5728\u6578\u503c\u6c42\u89e3\u6642\uff0codeint\u9700\u7d66\u4e88\u5fae\u5206\u65b9\u7a0b\u7cfb\u7d71\u7684\u51fd\u6578\u4e88\u5fae\u5206\u65b9\u7a0b\u7cfb\u7d71\u7684\u51fd\u6578(system)\u3001\u521d\u59cb\u503c\u3001\u521d\u59cb\u503c([x0,y0]))\u3001\u4ee5\u53ca\u6b32\u89e3\u7684\u6642\u9593\u9ede(t)\u3002" ] }, { "cell_type": "code", "collapsed": false, "input": [ "alpha=0.1\n", "beta=0.015\n", "gamma=0.0225\n", "delta=0.02\n", "def system(z,t):\n", " x,y=z[0],z[1]\n", " dxdt= x*(alpha-beta*y)\n", " dydt=-y*(gamma-delta*x)\n", " return [dxdt,dydt]\n", "\n", "t=np.linspace(0,300.,1000)\n", "x0,y0=1.0,1.0\n", "sol=odeint(system,[x0,y0],t)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "markdown", "metadata": {}, "source": [ "\u5e95\u4e0b\u662f\u4ee5matplotlib\u5c07\u5f97\u5230\u7684\u6578\u503c\u89e3\u532f\u51fa\uff0c\u9996\u5148\u662f\u5728\u76f8\u5e73\u9762\u4f5c\u5716\uff1a" ] }, { "cell_type": "code", "collapsed": false, "input": [ "X,Y=sol[:,0],sol[:,1]\n", "plt.plot(X,Y)\n", "plt.show() " ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAD9CAYAAACcJ53WAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XtYVWW+B/Dv4iKGV1Buggh5Se7sQGlKFBrUrDRMc9JT\nlmBnOuVpzJ4aS1Ork2lNZ1K6HE85ZSdHK0vTUtNGt5GNMiZk3i+BIgKmYHJTbuv88RO2pCKXvfZa\nC76f51nP2nvtzVq/Ur/75d3vel9FVVUVRERkGk56F0BERM3D4CYiMhkGNxGRyTC4iYhMhsFNRGQy\nDG4iIpNpNLhzc3ORmJiIsLAwhIeHY/HixQCAefPmISAgABaLBRaLBRs3bnRIsUREBCiNjeMuKChA\nQUEBoqOjUVpaipiYGKxZswaffPIJunTpghkzZjiyViIiAuDS2Iu+vr7w9fUFAHTu3BkhISHIy8sD\nAPC+HSIinahNlJ2drQYGBqolJSXqvHnz1D59+qiRkZFqSkqKWlxc3OC9ALhx48aNWwu2pmjSu0pK\nStSYmBh19erVqqqqamFhoVpbW6vW1taqs2bNUlNSUq4IbiOaO3eu3iVcgTU1DWtqOiPWxZqapqnZ\ned1RJVVVVRg3bhweeOABJCcnAwC8vb2hKAoURcHUqVORkZFxvdMQEZGdNBrcqqoiNTUVoaGhmD59\nev3x/Pz8+serV69GRESEdhUSEVEDjX45uX37dnz00UeIjIyExWIBAMyfPx8rVqxAVlYWFEVBcHAw\nlixZ4pBiWyshIUHvEq7AmpqGNTWdEetiTfbV6HDAFp9UUTjqhIiomZqanbxzkojIZBjcREQmw+Am\nIjIZBjcRkckwuImITIbBTURkMgxuIiKTYXATEZkMg5uIyGQY3EREJsPgJiIyGQY3EZHJMLiJiEyG\nwU1EZDIMbiIik2FwExGZDIObiMhkGNxERCbD4CYiMhkGNxGRyTC4iYhMhsFNRGQyDG4iIpNhcBMR\nmQyDm4jIZFz0LsCRqqqA3Fzg+HGgvBxQFMDJybZ3cwP8/YHevYEOHfSulojo6tpUcBcXAx98AKxd\nC1it2lzD1xf43e+AW26RfWwscMMN2lyLiOhqFFVVVbufVFGgwWmvcPYs8PTTwPvvX/31ESNku/lm\nIDhYWtOurtc/b00NkJ8vLfPjx4HsbCArC9ixAzh5svGfDQkB7rlHtsGDpSVPRNQUTc1O0wV3WRnw\n8MPAqlW2Y927A089BTz6KNCzpyaXvaaLF4Ht26WVv2aNBP3V+PkBU6cCKSlAUJBDSyQik2hzwb1/\nP2CxAJWV8nzIEODdd4GBA+16Gbs6flwC/eOPJdx/y91dwjw1FYiMdHx9RGQsbSa4z54FAgKACxfk\n+d//DkycaJdT66KqCtiwAXjvPWDduitf79VLun+mTgU6d3Z8fUSknzYR3NOnA4sWyePNm4GkpFaf\n0pBUFdi6FfjrX4Evv7zy9dRU6QoKCXF8bUTkOKYO7uJiwNNTHs+aBfzXf9mpMBPJzwfS0oDXXgOq\nqxu+Fh8v/0+GDtWnNiLShmmDe98+IDxcHmdn84u8OlVV8oXswoXAjz82fG3MGOCll9hPTmR2pgzu\nrVuB22+XxzU1HErXmNpa4NNPgdmzgaNHG7720EPA3LkyBJKIzKOp2dloNObm5iIxMRFhYWEIDw/H\n4sWLAQBFRUUYPnw4BgwYgBEjRuDcuXOtLjgjQ0L7ppukz5eh3TgnJ+APfwCOHJH/XxcvAm+/LV1M\ny5YBN94od4QqCjBnjgyjJKK2odEWd0FBAQoKChAdHY3S0lLExMRgzZo1eP/999GzZ08888wzWLhw\nIYqLi7FgwQLbSZvZ4s7Pl9EUPXoAZ8607j+IRGmpfLE7e3bD44MGAYsXy52fRGQsmnSVJCcnY9q0\naZg2bRq2bdsGHx8fFBQUICEhAQcPHmz2xQH5ld/Z2fZYUZpaDTXHvn0yMuXrrxsenzdPhh+6u+tS\nFhFdxu7BnZOTg2HDhmHv3r0IDAxEcXExAEBVVXh6etY/r7v43Llz658nJCQgISHhqucND5dQOX8e\n6NKlKZVQa1VWSrfKk082PB4XJ63xwYP1qYuovbFarbBeNrHSCy+8YL/gLi0txbBhw/D8888jOTkZ\nHh4eDYLa09MTRUVFtpM28VNj0yZg5Ejgiy9kZATpY+9eYMYMGSt/uaVLgSlT+FsQkaPY5ctJAKiq\nqsK4cePw4IMPIjk5GQDqu0gAID8/H97e3i0qcuRIwMODoa238HD5EFVVuUN1/nw5nppqm/Z29myg\nokLfOolINBrcqqoiNTUVoaGhmD59ev3xMWPGYNmyZQCAZcuW1Qd6c7zwguyzs5v9o6QhNzfg2Wcl\nxGtrZZ4Vd3fg5ZdlryjApElAYaHelRK1X412lXz33XcYOnQoIiMjoVz6ffmVV17B4MGDMWHCBJw4\ncQJBQUH45JNP0L17d9tJr9PcrxvuN2aMdJOQOezcKa3wfftsx+LiZLKviAj96iJqKwx9A86zzwIL\nFsit3HUjSshcjh8H/vM/G06U1auXtNCHDNGvLiIzM3RwK4qMI/7nP+19ZdLD+fPSB56WZjvm7y8h\nfttt+tVFZDZ2+3LS3upaaBs2OPrKpJWuXWUYoarKHZpPPgnk5UnLW1GAwEB+SBPZk8Nb3HVDyxyw\nshnprLwceO4529S8ANCnD7ByJe/cJLoaQ7a4a2pkf2lACrVx7u7AG2/Ih3RpqfSJHz8uiywrisyn\nkpGhd5VE5uPQ4H79ddk/+KAjr0pG0KmTrTulpAR4/HEZChoXJyEeFgYcOKB3lUTm4NCuEnaT0G+V\nlMhcKUuW2I6NHi13bXp56VcXkR4M2VUCAK++6ugrkpF16QL8z//Ih3l+vtxNu24d4O0tH/QzZsiU\ntURk47AW9w8/ALGxcku1m5u9r0htzZ49wNixwM8/24699RbwH//BuVOo7TLcOO7ERMBqZTcJNd9X\nXwF332177uwMrF8PjBihX01EWjBcV4nVCgwc6KirUVty113ygV9TI0MLa2qkS0VR5O/U/v16V0jk\nWA7t4/7t/M9EzeHkBDzxhIR4eTkwbRpw6JCMSFEUWWuTS7RRe+CQrpLcXLl7rrwcuOEGe1+N2rtT\np2T9ze++sx1bvlxmMSQyE0N1lbzzjuwZ2qSFXr2A9HRpiddNpfBv/yat8JAQ4OhRfesjsjeHBPfy\n5Y64ChFwxx0S4JWVssbmwYNA//4S4n/6kxwnMjuHdJUoChAZCfz4o72vRHR92dlAcrIMMayzbl3D\nkSpERmCorhJA+iCJ9BAcLI0GVZUJrgC5O1NRZN6UvDx96yNqLs2Du7ZW9hMmaH0louv7wx9so1Km\nTgV27AACAiTEX32V9xmQOWjeVbJzp0zhyX8QZFT79gFJScCl9a/Rr58snhwcrG9d1P4YpqvkH//Q\n+gpErRMWJvOk1NQAc+fKKJQbb5RW+H//NxsdZDyat7hHjAA2b+ZffjKX/ftlmobTp+X5wIHAxo2y\nEASRVgzT4rZatb4Ckf2FhgKFhbKg9axZMqwwKEha4XXzihPpRfMWt6IAMTHArl32vgqRY+3dCyQk\nAGfPyvPwcJnsqndvXcuiNsQwLW5A/rITmV14OHDmjLTCZ86UIA8MlMbJ22+zFU6O45Dgjo93xFWI\nHMPZGXjlFQnqH38EuneXpdicnIDbbgOKivSukNo6hwR3dLQjrkLkeJGRQHExUFUlMxd+/z3Qo4e0\nwrdt07s6aqs0De7iYtkHBmp5FSL9ubjIXOGqCnz9tRxLSJAAnzXLdiMakT1oGtw//SR7LjVF7cmI\nERLghYWAxQLMny/dK1FRcoyotRwS3ETtkbc3sHu3tLbnzJFJrnx9pSGzcaPe1ZGZMbiJNKYowAsv\nSCs8PV2OjRolx594QkapEDUHg5vIgYYMkQAvKgKGDgXS0gBXV7nF/sQJvasjs9A0uA8f1vLsRObl\n4SGjTmprgddekznD+/SRVvi6dXpXR0an6Z2TigJ07AhUVNj7CkRtzw8/ALGxtucvvgjMns0v99sT\nw9w5yUl5iJomJka6Uc6eldEoc+bITT333gtcuKB3dWQkmgc3x3ATNY+np4xGqaoCUlOB1atloe0+\nfWRFeyK2uIkMysUFeO89aYW/+aZ8eenvL10nO3fqXR3picFNZAKPPy4BvnWrPL/lFgnwZcv0rYv0\n0Whwp6SkwMfHBxEREfXH5s2bh4CAAFgsFlgsFmy8zp0EDG4i+0lIkADPzpY5UR5+WAJ8+nTeVt+e\nNBrcU6ZMuSKYFUXBjBkzkJmZiczMTNxxxx2NXoDBTWR/QUEyxWxpKTB8uMyT4uwsY8PPn9e7OtJa\no8EdHx8PDw+PK443ZwShv3/ziyKipunUSRY2rq0FnnlG7szs1k3GifOGnrbLpSU/lJaWhg8//BCx\nsbF4/fXX0b179yveM2/ePADy5crIkQlI4GoKRJpRFGDhQtmWLwceeMD22+6ePcBlvZ1kIFarFdYW\nrO943RtwcnJyMHr0aPx06f7106dPw8vLCwDw/PPPIz8/H0uXLm14UkXBxYsq3NykJcAbCIgcb+tW\n4Pbbbc+tVmDYMN3KoSbQ7AYcb29vKIoCRVEwdepUZGRkXPV9devyMbSJ9JGYaFulB7DND/7pp7qW\nRXbQ7ODOz8+vf7x69eoGI04ud+ZMy4siIvuJjJQAz8kBunQBJkyQAF+yRO/KqKUa7SqZOHEitm3b\nhjNnzsDHxwcvvPACrFYrsrKyoCgKgoODsWTJEvj4+DQ8qaJg61a1/hOfiIzj7FkZB370qDx//XVg\nxgx9ayLR1K4SzSaZWrVKxfjxDG4iozp/XrpPMjPl+YsvAs8/r2tJ7Z7uk0yxq4TI2Lp2lTlRysqA\n+HiZ1EpRgD//mQ0uo9MsuOu+nCQiY3N3B779VmYgHDUKePVVmZXwyScZ4EbF4CYiAICbG7B+PVBZ\nCYwdC7zxhgQ4W+DGw+AmogZcXYHPP5cAHz3a1gK/dE8dGYBmwV1SotWZicgRXF2BtWulC2X4cFnw\nWFGA+fP1row0C+6yMq3OTESO5OYm86FUVMgkVrNmSYAvXqx3Ze0Xg5uImqRjR1nguKwMGDQI+NOf\nJMBXrtS7svaHwU1EzeLuDmRkAMXFQEAAMHGiBPg//qF3Ze0Hg5uIWqR7dyA3Fzh5Ur68TEqSAN+9\nW+/K2j4GNxG1ir8/UFMDHDggz2NiJMBzcnQtq01jcBORXQwcKOO9//lPeR4cDAQGcoSZFjQL7vJy\nrc5MREZ2yy0S4J98Il0pXbsC99wjrXKyD82Cu7JS1sAjovbpvvskwOfOlfHgLi4yFpxaT7PZAQEV\nXbsCv/5q77MTkdnU1ADjxgFffCHPP/sMuPdefWsyIt2ndQVU+PkBp07Z++xEZFYlJUBoqIxEAYDD\nh4H+/fWtyUh0n9YVkBWoiYjqdOki/d4HD8rzAQOAfv3krkxqOgY3ETncTTdJ//fHHwPHjslNPU88\noXdV5qFpcLu7a3l2IjK7CROA2lrg0UeBtDQZ/13XD07Xpmlwu7lpeXYiagsUBXjnHaC0FOjdG0hO\nlmMFBXpXZlyaBrerq5ZnJ6K2pFMn4MQJICtLnvv5gevWXgODm4gMJSpKwnr+fBk26OQErFmjd1XG\nwuAmIkN69lkZbRIcLEupKQpw+rTeVRkDg5uIDKtjR+Dnn4HMTHnu4wM8/ri+NRkBg5uIDC86WrpP\nZs4E3n5bWt8//qh3VfphcBORabzyim0h8uhomdCqPU5exeAmIlPx9JTW9//9H7Bzp0xetW6d3lU5\nFoObiEzpgQdkFtJ+/YAxYwBvb3neHjC4ici0XF2BI0eALVuAX36Rm/7WrtW7Ku0xuInI9BITgepq\nICREFm0ICACqqvSuSjsMbiJqE5ydgf37gU2bgLw8oEMHYMMGvavSBoObiNqU4cOltR0cDNx5p20p\ntbaEwU1EbY6Li9y48+mnMvLEyQk4elTvquyHwU1Ebdb48cC5c/K4f39gzhx967EXBjcRtWnduklX\nyfTpwEsvyV2XZl9xh8FNRO3CX/8K7N0rj93dgYwMfetpDQY3EbUbYWEybLBzZyAuDpg1S++KWqbR\n4E5JSYGPjw8iIiLqjxUVFWH48OEYMGAARowYgXN1HUhX4eJiv0KJiOzB2VlWm3/uOZnz293dfPOd\nNBrcU6ZMwcaNGxscW7BgAYYPH47Dhw/j97//PRYsWHDNn2eLm4iM6uWXpbukokIamTk5elfUdI0G\nd3x8PDw8PBocW7t2LR566CEAwEMPPYQ1jSxNwRY3ERnZoEFAWZk8Dg4G1q/Xt56mana0FhYWwsfH\nBwDg4+ODwsLCa7xzHtasAY4dAxISEpCQkNCKMomItOHuLqNOYmOBu+4CZs+W0SeOYLVaYbVam/1z\niqo2fk9RTk4ORo8ejZ9++gkA4OHhgeLi4vrXPT09UVRU1PCkigJAxd//Dkyc2OyaiIh0MWeOhLbF\nAuze7fjrK4qC60QygBaMKvHx8UFBQQEAID8/H97e3tc+uaZjVoiI7OvFF2V+k8xMGe9t1C8tmx2t\nY8aMwbJlywAAy5YtQ3Jy8jXfqygtL4yISA933CFdvIB8T1derm89V9NocE+cOBG33norDh06hN69\ne+P999/HzJkzsXnzZgwYMABbtmzBzJkzr/nzDG4iMqMbb7Qtkdapk8z1bSTX7eNu0Ukv9XGvWgWM\nG2fvsxMROUZFhXx5CcgkVX37ans9zfq4m1eElmcnItLWDTfInZaALJFmlBkGNQ1ufjlJRGbn7AzU\n1srj/v2NcaMOW9xERNehKLbwDg4GTp7Utx4GNxFRE1w+PLB3b9uXl3pgVwkRURM5Odn6vHv2tD12\neB1anpwtbiJqa5ydgfPn5bFeE+kxuImImqlLF+DIEXl8882Ovz67SoiIWqBfP9vt8X/5i2OvzRY3\nEVEL3XEHkJoKPP00kJ3tuOsyuImIWuG992R/440yPawjsKuEiKiV6lZwjI93zPXY4iYiaqVu3YAP\nPgC2bwcOHdL+egxuIiI7uLSiIwYO1P5aDG4iIjs5flz2S5Zoex32cRMR2UlgIBASAjz6qLbXYXAT\nEdnRd9/J/rPPtLsGu0qIiOzI0xNwcwPGj9fuGgxuIiI7+9e/ZF83TNDe2FVCRGRnERGyf+YZbc7P\nFjcRkQb69gXefVebczO4iYg0MHeududmVwkRkQbuuUf2JSX2Pzdb3EREGujaVfYbN9r/3AxuIiIN\nnTpl/3NqFtyKAnh7a3V2IiJzqKy0/zk1C24fH1mbjYioPXNzs/852VVCRKSBukUVhg+3/7k1DW5H\nrQZBRGQ0e/bIXotpXjULbicnBjcRtV9z5shei54HTb+crK3V6uxERMa2di0wdKg259Y0uNniJqL2\naO9e2S9dqs35GdxERHZWN8lUv37anJ993EREdrR9u+x37dLuGuzjJiKyk+pqYMgQoFMnICZGu+uw\nq4SIyE4CA2VfVKTtdVxa+oNBQUHo2rUrnJ2d4erqioyMjAavs6uEiNqTf/93ID8f2LIF6NBB22u1\nOLgVRYHVaoWnp+c1XmdXCRG1D3/5iyya8NZbQGKi9tdrVVeJ2kiTml0lRNQezJwJPP00MH068Nhj\njrlmq1rcSUlJcHZ2xh//+Ec88sgjDV4vKpqHtDSgZ08gISEBCQkJra2ViMhQ7rlHbrR5+WXgueea\n//NWqxVWq7XZP6eojTWbG5Gfnw8/Pz/88ssvGD58ONLS0hAfHy8nVRQMHKjis8+A0NCWnJ2IyLjK\ny2XkCAAsXw5MmmSf8yqK0mhPRp0Wd5X4+fkBALy8vDB27NgrvpxkVwkRtUXp6bbQPnzYfqHdHC0K\n7vLycpRcWkitrKwMmzZtQkTdrUJ1J+aoEiJqQ2pq5I7IoUMBT08Zs92/vz61tCi4CwsLER8fj+jo\naMTFxeHuu+/GiBEjGryHo0qIqK1YsQJwcZE5SD77DDh7Vt+FYlr05WRwcDCysrIafQ+7SojI7Pbt\nA8LD5XFoKJCVBbi66lsT0IpRJdfDrhIiMqucHCA42Pb82DHgxht1K+cKnKuEiOiSY8cku+pCe9Mm\naYAaKbQBzlVCRIQNGySz6qZh/fpryS8t1ou0B07rSkTtUlUV8MQTEth33inHdu2S3PrNWAvD0ayP\nm10lRGREO3cCCQnAhQvyfPBg4MsvAS8vXctqFnaVEFGbd+qUTP6kKMAtt0hov/uuZNTOneYKbUDD\n4HZ2lgHrRER6KCwExo2TsPb3B6xWYMIEoLhYAnvqVL0rbDnNuko6dgQuXtTq7EREV9qzR2bp27rV\ndux3vwM++AAYMEC3suxOs+B2c7P1IRERaaG6GnjnHfmS8XIJCcB77wF9++pSluY06ypxc2OLm4js\n7+BBIDlZukBcXW2hPW8eUFYm3SBbt7bd0AY07iphi5uIWmvfPuCll4CPP254PDISePNN4NJs0u2K\npl0lbHETUXNlZUlQf/55w+NeXsDcuUBqqjQM2zO2uIlIN+fOAUuXAmlpwPHjDV/z9wfmzAEeflj7\nxXfNhi1uInKI6mpg/XoJ6W++ufL1xETprx49Wt8pU81As+B2d5cvCoio/SkuBlauBN5/H/jXv658\nPSBAQjo1VRYloObRLLh79ABOn9bq7ERkBLW1wI4dwCefyGIDV/s336cPMGUKMHlyw6lSqeU0C+6e\nPYEDB7Q6OxE50uUB/emncgv51QwdKn3S48cDXbo4tMR2RdPgPnNGq7MTkRays4GvvpK+6A0brv2+\nm28G7rtPArpuKlRyHM2C28tL5gogImPJz5dgXr9eQrqxQQRJSTLl6ZgxbfuGFrPRLLj79QOOHJG7\nmBRFq6sQ0W/9+iuQni6TKlmtwA8/NP7+IUMknO+8U25q4b9X49O0q8TJCfjlF8DbW6urELU/hYXS\n31wXztcLZjc3mbsjIUGG3MXGcrid2Wm6kEJ0NJCRAdx9t1ZXIWo7VFX6mHfulGDesUP+/VxPhw4S\nyHXhHBNjjJXISTuaBTcg/WObNzO4qX2rqQEOH5Zbueu2zEz5bbQpYmJk8v+6rW9fdme0d4qq2n+d\nGkVRoKoqDh6UFsDx4/LrGlFboapAbi6wfz/w448SxFlZwKFDTT+HtzdgschvpnFxsvXqpV3NZHx1\n2Xnd92kZ3AAwahQwcqRMbk5kVGVlEroHDsi0oQcOSCg3914EV1cJ4su3iAiOaaamMUxwHz4M3Hab\nDDsaPNjeVyJqqKJC+omPHbNtP/8s+6NHW7acXr9+QEgIEBoq+5AQYOBAoGtX+9dP7ZthghuQFZRT\nUoC337atAUd0NVVVQF6edEPk5gInTsj++HEgJ0e21syB4+8vS1jddJPs+/eXIO7TB3DR9Bsfousz\nVHAD8u345MmAnx/w1FPSfcJvvtsGVQVKSoCCAtuWny+3ReflyXbypIRwZWXrr+fhIXNe9O175ebv\nz6FuZF6GC25ApnX86CPgf/9Xfm1NSpJhTIMGya+e7X1ydK1VVQFFRcDZs7KdOSMjG06fbrgvLJTw\n1WLKgp49gd69gcDAhvu6rVcvtnyp/TJkcF8uJwfYskXWhsvMlD7IPn2kJVX3jzggQGYZ9PSUzcND\n9mYdoaKq0uIsKQHOn2+4lZTIHW91j8+fv/J5cbFs587Jh6DWFAXw9ZXNz8/2+PJjPj7yuHNndoER\ntZbhg/u3Kivli8ycHFv/Zl6etBCLiiSw6va1tTLf9w03XLm5uUmLzdm54ebiIjObtRWdOskHWc+e\n8uFW9wF3+fMePWzPPT2Bbt3kblYiMibTBXdTqar8yl9RYdvKy22PL16UkQPV1bK/fFu1SoK/pkZa\nh3Wbk1PDx5c/d3WVzcVFPhg6drz2/mrHOnWybZ07yzG2TInoatpscBMRtVVNzU7+4kxEZDIMbiIi\nk2lXwW21WvUu4QqsqWlYU9MZsS7WZF8tDu6NGzdi4MCB6N+/PxYuXGjPmjRjxD8o1tQ0rKnpjFgX\na7KvFgV3TU0Npk2bho0bN2L//v1YsWIFDnBlYCIih2hRcGdkZKBfv34ICgqCq6sr7r//fnzxxRf2\nro2IiK6iRcMBV61aha+//hrvvvsuAOCjjz7Czp07kZaWJiflQGUiohZpSiS3aFaI6wUzx3ATEWmn\nRV0l/v7+yM3NrX+em5uLgIAAuxVFRETX1qLgjo2NxZEjR5CTk4PKykp8/PHHGDNmjL1rIyKiq2hR\nV4mLiwvefPNNjBw5EjU1NUhNTUVISIi9ayMioqto8TjuUaNG4dChQzh69CieffbZ+uNGHN+dkpIC\nHx8fRERE6F1KvdzcXCQmJiIsLAzh4eFYvHix3iXhwoULiIuLQ3R0NEJDQxv8ueqtpqYGFosFo0eP\n1ruUekFBQYiMjITFYsFgg6zLd+7cOYwfPx4hISEIDQ3Fjh07dK3n0KFDsFgs9Vu3bt0M8Xf9lVde\nQVhYGCIiIjBp0iRcvHhR75KwaNEiREREIDw8HIsWLWr8zaodVVdXq3379lWzs7PVyspKNSoqSt2/\nf789L9Ei3377rbp79241PDxc71Lq5efnq5mZmaqqqmpJSYk6YMAAQ/y/KisrU1VVVauqqtS4uDg1\nPT1d54rE66+/rk6aNEkdPXq03qXUCwoKUs+ePat3GQ1MnjxZXbp0qaqq8md47tw5nSuyqampUX19\nfdUTJ07oWkd2drYaHBysXrhwQVVVVZ0wYYL6wQcf6FrTTz/9pIaHh6sVFRVqdXW1mpSUpB49evSa\n77frLe9GHd8dHx8PDw8PvctowNfXF9HR0QCAzp07IyQkBKdOndK5KsDd3R0AUFlZiZqaGnh6eupc\nEXDy5EmsX78eU6dONdyIJSPV8+uvvyI9PR0pKSkApEuzW7duOldl880336Bv377o3bu3rnV07doV\nrq6uKC8vR3V1NcrLy+Hv769rTQcPHkRcXBw6duwIZ2dnDBs2DJ9//vk132/X4M7Ly2vwhxIQEIC8\nvDx7XqJNysnJQWZmJuLi4vQuBbW1tYiOjoaPjw8SExMRGhqqd0l48skn8dprr8HJYKtAKIqCpKQk\nxMbG1t8EZ840AAADKElEQVTToKfs7Gx4eXlhypQpuPnmm/HII4+gvLxc77LqrVy5EpMmTdK7DHh6\neuKpp55CYGAgevXqhe7duyMpKUnXmsLDw5Geno6ioiKUl5fjq6++wsmTJ6/5frv+S+CNN81XWlqK\n8ePHY9GiRejcubPe5cDJyQlZWVk4efIkvv32W93nc/jyyy/h7e0Ni8ViqNYtAGzfvh2ZmZnYsGED\n3nrrLaSnp+taT3V1NXbv3o3HHnsMu3fvRqdOnbBgwQJda6pTWVmJdevW4b777tO7FBw7dgxvvPEG\ncnJycOrUKZSWlmL58uW61jRw4ED8+c9/xogRIzBq1ChYLJZGGyp2DW6O726eqqoqjBs3Dg888ACS\nk5P1LqeBbt264a677sKuXbt0reP777/H2rVrERwcjIkTJ2LLli2YPHmyrjXV8fPzAwB4eXlh7Nix\nyMjI0LWegIAABAQEYNCgQQCA8ePHY/fu3brWVGfDhg2IiYmBl5eX3qVg165duPXWW9GjRw+4uLjg\n3nvvxffff693WUhJScGuXbuwbds2dO/eHTfddNM132vX4Ob47qZTVRWpqakIDQ3F9OnT9S4HAHDm\nzBmcO3cOAFBRUYHNmzfDYrHoWtP8+fORm5uL7OxsrFy5Erfffjs+/PBDXWsCgPLycpSUlAAAysrK\nsGnTJt1HLfn6+qJ37944fPgwAOlTDgsL07WmOitWrMDEiRP1LgOAtG537NiBiooKqKqKb775xhBd\ngqdPnwYAnDhxAqtXr268W8ne346uX79eHTBggNq3b191/vz59j59i9x///2qn5+f2qFDBzUgIED9\n29/+pndJanp6uqooihoVFaVGR0er0dHR6oYNG3Stac+eParFYlGjoqLUiIgI9dVXX9W1nt+yWq2G\nGVXy888/q1FRUWpUVJQaFhZmmL/rWVlZamxsrBoZGamOHTvWEKNKSktL1R49eqjnz5/Xu5R6Cxcu\nVENDQ9Xw8HB18uTJamVlpd4lqfHx8WpoaKgaFRWlbtmypdH3arLmJBERacdYX9MTEdF1MbiJiEyG\nwU1EZDIMbiIik2FwExGZDIObiMhk/h8naJmxaD2+WgAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 3 }, { "cell_type": "markdown", "metadata": {}, "source": [ "\u5e95\u4e0b\u662fx\u3001y\u5c0dt\u4f5c\u5716\uff0c\u5176\u4e2dx\u7684\u5716\u5f62\u4ee5\u7d05\u8272(r)\u8868\u793a\uff0cy\u7684\u5716\u5f62\u4ee5\u85cd\u8272(b)\u8868\u793a\u3002" ] }, { "cell_type": "code", "collapsed": false, "input": [ "px, = plt.plot(t,X,'r-')\n", "py, = plt.plot(t,Y,'b-')\n", "plt.xlabel(\"t\") \n", "plt.ylabel(\"x,y\") \n", "plt.legend([px,py],[\"x\",\"y\"])\n", "plt.show() " ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEKCAYAAAAfGVI8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlcVPX3P/DXgJgLLrgAChZuCCrgpLmkFn4Qy1xLMy2V\nBFv0Y+Unf5bVp8KvWS5puWVlaZiFmrlViphGmrsJH0nFUMEQARdwQVFZ7u+P46AIA7Pc/Z7n4+Ej\nm+2ecebOue/tvE2CIAhgjDFmWC5KB8AYY0xZnAgYY8zgOBEwxpjBcSJgjDGD40TAGGMGx4mAMcYM\nTvREkJGRgV69eqFdu3Zo37495s+fDwCIjo6Gr68vzGYzzGYz4uLixD40Y4wxB5jEXkeQnZ2N7Oxs\ndOjQAfn5+ejYsSPWr1+P1atXo06dOnj99dfFPBxjjDEnVRP7Bb29veHt7Q0AcHd3R2BgIDIzMwEA\nvHaNMcbUR/QWwd3S09Px6KOP4siRI5gzZw6WLVuGevXqoVOnTpgzZw7q169/JxCTSaowGGNM15z9\nGZdssDg/Px9Dhw7FvHnz4O7ujnHjxiEtLQ1JSUlo0qQJJk2aVO45giDo9s/777+veAz8/vj9GfH9\n6fm9CYI41/GSJILCwkIMGTIEI0eOxODBgwEAnp6eMJlMMJlMGDt2LPbv3y/FoRljjNlJ9EQgCAKi\noqLQtm1bTJw4sfT2rKys0r+vW7cOQUFBYh+aMcaYA0QfLN61axdWrFiB4OBgmM1mAMCHH36I2NhY\nJCUlwWQyoXnz5vjiiy/EPrSqhYaGKh2CpPj9aZue35+e35tYJB0stofJZBKtv4sxxoxCjN9O0VsE\njDGmNg0aNEBeXp7SYTjFw8MDubm5krw2twgYY7qnh98Xa+9BjPfGtYYYY8zgOBEwxpjBcSJgjDGD\n40TAGGMGx4mAMcYMjhMBY4wZHCcCxhhT0MmTJ9GwYUMkJiYCAM6ePYvGjRtjx44dssXAiYAxxhTU\nsmVLzJw5EyNHjkRBQQHGjBmDMWPG4JFHHpEtBl5QxhjTPS38vgwaNAinTp2Cq6srDhw4ADc3tzL3\n84IyxhiTkskkzh8njB07FkeOHMErr7xSLglIjVsEjDHdU/vvS35+PkJCQhAWFoZNmzYhOTkZHh4e\nZR4jZYuAEwFjTPfU/vsSFRWF69evIzY2Fi+99BIuXbqEVatWlXkMdw0xxphObdiwAfHx8Vi8eDEA\nYO7cuTh06BBiY2Nli4FbBIwx3dPD7wu3CBhjjEmGEwFjjBkcJwLGGDM4TgSMMWZwnAgYY8zgOBEw\nxpjBcSJgjDGD40TAGGMGx4mAMcYMjhMBY4wZHCcCkQkCsHw5MHo08MMPSkfDmP4VFgIzZwKRkcCu\nXUpHo02cCET26afAjBlA9+7Ae+8BU6cqHRFj+hYVBWzfDoSEAEOGAGvWKB2RfWbPno2hQ4eWue3V\nV1/FxIkTZYuBi86JKCsLaNcO+PNPoHlz4Nw5oGNH4IsvgCeeUDo6xvRn2zbgpZeAw4eBWrWApCQg\nPBzYv5/OQQs1/75kZ2ejVatWyMzMRL169VBUVAQfHx/ExcXBbDaXPo6LzmnE/PnAqFF3voCenkBM\nDDB+PFBQoGxsjOnR7NnAO+9QEgCADh2AN9+k5KAV3t7e6NmzJ3643ZccFxeHxo0bl0kCUuNEIJLi\nYuDbb4EXXih7+7/+Ra2CTz9VJi7G9Cozk678R4woe/vEicDp08Cvv9r+WkrvVBkREYEVK1YAAFas\nWIFRo0Y5/mIO4K4hkezaBYwbR03Ue6WkAI8+Cpw6BdSuLX9sjOnRwoWUCJYvL3/f6tXA3LnAnj2W\nH2l1/74UFBTAx8cHO3bsQLdu3XDs2DH4+vqWeQx3DWnAli3WxwECAmjweOlSeWNiTM82bwYGDqz4\nviFDgIsXgd275Y3JUTVr1sSQIUPw7LPPokuXLuWSgNQ4EYgkPh547DHr90+ZAsyZQ11IjDHnFBdT\nK7xnz4rvd3WlLqK5c+WNyxkRERH466+/ZO8WAjgRiOLmTeoS6tzZ+mM6dwa8vIC4OPniYkyvkpPp\nfPLysv6Y558HduwA0tPliso5DzzwQGnLQG6iJ4KMjAz06tUL7dq1Q/v27TF//nwAQG5uLsLDw+Hv\n748+ffrg0qVLYh9aMUlJgL9/1f3/L70EfP65PDExpmc7d1pvDVjUrg08+yywbJk8MTmjpKQEc+bM\nwYgRI+Du7i778UVPBG5ubvjkk09w5MgR7N27F4sWLcKxY8cwY8YMhIeH4++//0ZYWBhmzJgh9qEV\ns38/0KVL1Y975hnqs8zIkD4mxvRs714ad6tKVJT6E8G1a9dQt25dbNu2DVMVWoEqeiLw9vZGhw4d\nAADu7u4IDAxEZmYmNm7ciIiICADUF7Z+/XqxD62Y/fsr7xayqF2bprp9/bX0MTGmZ4cP05qBqgQH\nA97e0sfjjNq1ayM/Px/Jycnw8fFRJAZJxwjS09ORmJiILl26ICcnB163O/S8vLyQk5Mj5aFl9eef\ntFbAFpGRNN1NxTPZGFO1mzeBEyeAwEDbHh8VJW08elBNqhfOz8/HkCFDMG/ePNSpU6fMfSaTCaYK\nVl9ER0eX/j00NBShoaFShSeawkJaHxAQYNvjzWagZk2a8dCjh7SxMaZHKSm0er9GDdseP2IE8PLL\n0sYkp4SEBCQkJIj6mpIsKCssLET//v3Rt2/f0sJJAQEBSEhIgLe3N7KystCrVy+kpKTcCUTlCz6s\nOXaM5jKnptr+nBkzaCYDDxwzZr8VK4CffwZWrrT9OVr9fbmbphaUCYKAqKgotG3btkz1vIEDByIm\nJgYAEBMTg8GDB4t9aEUcO2Z7E9Xi2WepRPXNm9LExJieHT4MBAUpHYW+iJ4Idu3ahRUrVuC3336D\n2WyG2WxGXFwcpkyZgq1bt8Lf3x/bt2/HlClTxD60IhxJBPffT4NYv/wiTUyM6VlKiv3nnIeHR2mX\ntFb/eHh4SPMPCgnGCHr06IGSkpIK7/vVnipQGnH0KNCnj/3PGzWKitQ99ZT4MTGmZydPAq1b2/ec\n3NxcvPwy4OdHq/xZWbyy2EkpKbYPFN9t6FDaTENH6+oYk1xJCU3OaNHC/uc++ywQGyt+THrAicBJ\np04BLVva/7y6dYGwMEBHyykYk1xmJuDh4VgV3x49gNxc4K+/xI9L6zgROOHSJSp+1bChY89/5hlg\n1SpxY2JMz06cAFq1cuy5Li40lZRbBeVxInBCWhrNZ3Z0Q4r+/ankxIUL4sbFmF6dPOl4IgCoe+j7\n73lB5704ETjh1Kmy+6Laq3Zt4PHHgbVrxYuJMT07ccKxrliLkBBaiLZ3r3gx6QEnAidYWgTO4O4h\nxmzn6ECxhclE3UP2LEYzAk4EThAjEfTtCxw6BOio9BJjkjlzBmjWzLnXeOYZ2sqSN4m6gxOBE8RI\nBDVrAv36AWvWiBMTY3qWkeF8ImjTBmjShDatYYQTgRNOn6YFKs7i7iHGqlZURC3npk2df63hw7l7\n6G6SFJ1zhBaLQtWrR8XjnF35ffMmXaEkJwMKlSNnTPUyMoCuXWktgbNOn6bS8VlZgJub86+nJFUW\nnTOKq1epBHX9+s6/1n33AYMGUSE6xljFMjIAX19xXuuBB6iLaOtWcV5P6zgROCgzk76Ujq4huBd3\nDzFWOTHGB+7G3UN3cCJwUGamuN04YWE0R/r0afFekzE9ETsRPP008NNPQEGBeK+pVZwIHCR2InBz\nA558kqa1McbKEzsReHsDDz4IbN4s3mtqFScCB4mdCADqHuKmKmMVE2MNwb24e4hwInCQZYxATKGh\nwNmz9m17yZhRiN0iAGg/kC1baPKHkXEicJAULQJXV+q35CsUxsqTIhE0bAj07Als3Cju62oNJwIH\nnTkjzZx/S5lcjS2pYExSRUXA+fOAl5f4r83dQ5wIHCZFiwCgBTPXr9MG3YwxcvEiLdysJvrmusDA\ngVRuIjdX/NfWCk4EDigqoj0EvL3Ff22Tia5QePMMxu44dw7w9JTmtevWBcLDgXXrpHl9LeBE4IDs\nbKBRI2muToA7ZXK5e4gxImUiAPjiixOBA7KyxCl8ZU1wMFCrFm+ewZjFuXPSjA9Y9OsHHDxIF3lG\nxInAAVJfnVg2zzDyFQpjd5P6nKtZExgwwLjl4DkROCAnR9qrE4CaqqtX03gEY0YndSIAjD17iBOB\nA+RIBK1b04K1hARpj8OYFsiRCMLDgWPHgH/+kfY4asSJwAE5OdJ/KQHuHmLMQo5EUL06rTQ2Yr0v\nTgQOkHrgyuKZZ4D162njGsaMTI5EABh39hAnAgfI0TUEUNdQ+/ZUC4UxI5MrEVjqff39t/THUhNO\nBA6QKxEAxr1CYexuciUCS70vo20SxYnAAXJ9KQH6Um7eDOTny3M8xtTm+nXaFrZOHXmOZ7n4MtKC\nTk4EdioqAvLyaGWxHBo1Ah55BFi7Vp7jMaY258/ThZdY28JWpWtX4No14K+/5DmeGnAisNOFC0CD\nBtSElMvo0cDy5fIdjzE1kbMFDgAuLjRR4/vv5Tum0jgR2EmuqaN3698fSEykeuyMGY3ciQAARo0C\nvv0WKC6W97hK4URgJ7mmjt6tRg0aK1ixQt7jMqYGSiSCoCA6z7dvl/e4SuFEYCc5ZwzdzdI9ZKQB\nLMYAZRIBADz/PBATI/9xlcCJwE5KdA0BQLduNFB98KD8x2ZMSUolghEjgJ9/Bq5ckf/YchM9EURG\nRsLLywtBQUGlt0VHR8PX1xdmsxlmsxlxcXFiH1Y2SrUITCYeNGbGpFQiaNQI6NUL+OEH+Y8tN9ET\nwZgxY8r90JtMJrz++utITExEYmIiHn/8cbEPKxslxggsRo6k6oi3bilzfMaUoFQiAIzTPSR6IujZ\nsyc8PDzK3S7opHNbqa4hAGjeHGjbFti0SZnjM6YEJRNB375ASgpw8qQyx5eLRJstlrdgwQIsX74c\nnTp1wpw5c1C/fv1yj4mOji79e2hoKEJDQ+UKz2ZKdQ1ZWLqHBg9WLgbG5KRkIqhencYKli8Hpk5V\nJoZ7JSQkIEHk+vQmQYJL9fT0dAwYMADJyckAgHPnzqFx48YAgHfffRdZWVn4+uuvywZiMmmi1eDr\nC+zZAzRrpszxL18GHngASE0Fbv+TMqZbggDcdx9w9Sr9VwmJiVSe+uRJWmymNmL8dsrytjw9PWEy\nmWAymTB27Fjs379fjsOKThCUvToBgHr1gIEDabELY3p36RLt361UEgCADh2oztHOncrFIDVZEkFW\nVlbp39etW1dmRpGW5OUp/6UEgBdeAL76itcUMP1T+sILoBl7zz8PLF2qbBxSEn2MYMSIEfj9999x\n4cIFNGvWDFOnTkVCQgKSkpJgMpnQvHlzfPHFF2IfVhZq+FICQI8eQEkJsHs30L270tEwJh21nHOj\nRwOtWtHFYAVzYTRP9EQQW0Hx/MjISLEPowi1fClNJmDsWGoVcCJgeqaWc65RI5pB9N13wIQJSkcj\nPhUOfaiXpRyuGoweTdtYXr6sdCSMSUctiQCgLtkvv9RnlywnAjucP6+emTqenkDv3rx7GdM3NSWC\n0FCgoADQ6FyXSnEisMO5c+pJBABdoSxZonQUjElHTYnAxeVOq0BvOBHYQU0tAoBaBBcvAocOKR0J\nY9JQUyIAgIgI2i1Qb4XoOBHYQW2JwMUFiIqiQWPG9EhticDLiy7A9LZ7GScCO6hpsNgiMpIK0ent\nCoUxQH2JAABefBH44gt9DRpzIrCD2sYIAMDHBwgL45XGTJ/UmAjCwmi23oEDSkciHk4EdlBb15DF\nhAnAokX6ukJhrLCQWroNGigdSVkuLsDLLwMLFyodiXg4EdiopIQGZhs1UjqS8h55hL6cv/2mdCSM\niefCBaBhQ3UWeouKAn76iVoseqDCf2J1unQJqF2bytKqjckE/Pvf+rpCYUyN3UIWDRsCQ4boZ/o2\nJwIbqflLCQCjRgG//w5kZCgdCWPiUPs598orwOLF1IWldVYTwfz585GXlydnLKqm1vEBC3d34Lnn\naDYDY3qg9kQQEgK0aEGlXrTOaiLIycnBQw89hGHDhiEuLk4Tm8ZISe2JAADGj6c1BTdvKh0JY85T\nejdAW0yYACxYoHQUzrOaCKZPn46///4bkZGR+Oabb9C6dWu8/fbbOKn3zTut0EIiCAgAgoKA1auV\njoQx56m9RQAATz4JnDoF/O9/SkfinErHCFxcXODt7Q0vLy+4uroiLy8PQ4cOxeTJk+WKTzXUuIag\nIhMnAnPn8lRSpn1aSARubjSVVOutAquJYN68eejYsSPeeOMNdO/eHX/99RcWL16MP//8E2vXrpUz\nRlVQ46riivTtC9y4wVNJmfZpIREAtNL4xx/pN0KrrCaC3NxcrF27FvHx8Rg2bBjc3NzoCS4u+Omn\nn2QLUC200DUE0Jzr118H5sxROhLGnKOVRODpCTz9NC3q1CqToJJRYJPJpOoB6d69gTfeAPr0UTqS\nqt24Afj5UasgMFDpaBhzjOU73Ly50pFULSWFFnamp9O+5nIS47fT5nUEAQEBCAgIwEKDrlrSyhgB\nANSoQTOI5s5VOhLGHCMI2mkRADRRo1s3ICZG6UgcY7VFcPToUbRt27bMbT/88ANq1aqFfv36iR+I\nylsETZoABw9SkTctOH8e8PenKxW1T8Fj7F75+fS9vXZN6Uhs98cfwPPPA8ePA66u8h1X0hbBsGHD\nMHPmTAiCgOvXr+OVV17BnDlzJEkCaicIVPdEjXWGrGncGHjmGeCzz5SOhDH7aak1YNG9O513Wlxg\nZjUR7Nu3DxkZGejWrRs6d+6MJk2aYPfu3XLGphqXLlG/3333KR2Jff7zH1oCf/260pEwZh8tJgKT\nCZg8GZg9W3vTt60mgmrVqqFmzZooKCjAjRs30KJFC7iosQygDLQyY+hebdoAPXrwDmZMe7SYCABg\n0CCqUrxrl9KR2MfqL3vnzp1Ro0YNHDx4EDt37sT333+Pp59+Ws7YVENLA8X3eucdYNYsLjvBtEWr\nicDVFZg0CZgxQ+lI7GM1EXz11VeYNm0a3Nzc0KRJE2zcuBEDBgyQMzbV0Mpisop07EhlJ7Q6m4EZ\nk1YTAUADxklJwKFDSkdiO6uJ4KGHHip32+jRoyUNRq202jVk8d//0hVKUZHSkTBmGy0ngho1aKzg\ngw+UjsR2xuz0t5PWE0H37sD99wOxsUpHwphttJwIAOCFF4A9e4DDh5WOxDacCGyg5TECi//+F/jw\nQ9pykzG103oiqFWLxgqmT1c6EttwIrCB1lsEABAWBtStCxiwXiDTIK0nAoCqkiYkAMeOKR1J1TgR\n2EDLg8UWJhO1CqZN41YBUz89JAJ3d+C117TRKuBEYAM9tAgAoH9/WhS3Zo3SkTBmXUkJzcXX0kp+\nayZMALZsAVJTlY6kcpwIbKCHMQKAWgUffAC89x7PIGLqdfEidWPernyvaXXrUqvg/feVjqRynAiq\nYKkzpIdEAADh4VTMa8UKpSNhrGJa2KvYHhMnAtu3q3s7S04EVbh8mbpTatRQOhJxmEzUZzl1KnDr\nltLRMFae3hKBuzvw1ls0RqdWnAiqoIeB4nv16EH107kGEVMjvSUCgGYQHT6s3hpEnAiqkJOjv0QA\n0FjB9OlcmZSpjx4TwX33AdHRwNtvq7MyqeiJIDIyEl5eXggKCiq9LTc3F+Hh4fD390efPn1w6dIl\nsQ8rmZwcwNtb6SjE17Ej0LUrYNAN55iK6TERAMCoUTTxJD5e6UjKEz0RjBkzBnFxcWVumzFjBsLD\nw/H3338jLCwMMzRUmi87W5+JAKCVxrNn02A4Y2qh10RQrRqt43n7bfWt5RE9EfTs2RMeHh5lbtu4\ncSMiIiIAABEREVivoS189PqlBGi/guHDaeCYMbXQ8zk3ZAglhO++UzqSsqrJcZCcnBx43f5kvby8\nkJOTU+HjoqOjS/8eGhqK0NBQGaKrXHY28OCDSkchnfffp4HjCRMoMTCmND0nApMJ+OQT2kb2qaeA\n2rXtf42EhAQkJCSIG5e1zeudkZ6ejgEDBiA5ORkA4OHhgby8vNL7GzRogNzc3LKBqHTz+oEDgchI\nYPBgpSORzqxZNJthwwalI2EM8PUFdu+mirl6NXw4EBgozkIzSTevF5OXlxeys7MBAFlZWfDU0DQc\nPV+dWLz6Kk1tE/kigzG7CYI+6gxVZeZMYP58IDNT6UiILIlg4MCBiLm9RVZMTAwGa+jyWs+DxRY1\nagAffURlc9U2iMWMJS+PSjjrZQGnNQ88ALz0Eg0cq4HoiWDEiBF4+OGHcfz4cTRr1gzLli3DlClT\nsHXrVvj7+2P79u2YMmWK2IeVhCAYo0UAUJ9l9erAsmVKR8KMzCjnG0CrjbduBQ4eVDoSicYIHKHG\nMYJLlyhzX76sdCTyOHQI6NuX6qc3aKB0NMyIEhKoKOKOHUpHIo+vvgKWLgX++ANwcfCyXDNjBFqV\nnW2cqxOAZkcNHarumihM34zUIgBoIkpxMfDNN8rGwYmgEnpdVVyZDz6gXcz+/FPpSJgRGS0RuLgA\nixdTN9HFiwrGodyh1c8IA8X38vCgFcfjx/PAMZOf0RIBQC3xZ56hZKAUTgSVMOKXEgCef56uVJYu\nVToSZjRGPeemTQN+/hnYu1eZ43MiqIQRWwQAJYHPPqOpbVYWgTMmCaONy1nUq0d1v8aPV2b3QE4E\nlTBqIgAAs5laBq+9pnQkzEjOngV8fJSOQhnPPgvUrw8sWCD/sTkR3Ov8eVpdtWoVcnIEQ16dWERH\n0xznjRuVjoTpWkoKfdl27MDZs0DTpkoHpAyTCfjyS9on5ORJeY/NieBuOTlA5870KUybhuykbMO2\nCABa4blkCTVXjbKWgsns4EGgZ0/g8mUUDnsOuRdLDH3x1aoVdcmOHSvvZA1OBHcbPx4YMYJWecTH\nIyvLhCaF/ygdlaJ69QKeeAJ4802lI2G6U1gIPPccDUh98gmyF/yAxsI5uBbdVDoyRb32GlBQQK0D\nuXAisPjjD7o6efddAEBh46a4YGoErxVzFA5MebNm0YwGLkrHRPX554CfH61iBHD2/q5oWvuK+or1\ny8zVlWbsvfsu8I9M16GcCCw++IBqwtasCYAGij09gWqrvgNu3FA4OGXVr0/nbGQkcOWK0tEwXSgs\npBKcM2ZQ5zhuDxQH1FF+ma0KtG0L/Oc/wAsvyLPHMScCgAarEhNp2P62M2cAn/urUdFwvhRG//5A\n797AxIlKR8J04ccfqUPcbC69KTMTaGr2BP76C8jKUjA4dZg8mVYbf/659MfiRADQv/TYsWVq32Zm\n0gYZ6NsXuGcPZqOaO5eKga1bp3QkTPM++wx45ZUyN509C/g0cwX+9S917vAuMzc3YMUK6iJKSZH2\nWJwICguB2FhgzJgyN585c3s+c9++wObNysSmMu7uwLffAuPGUdcZYw5JT6cStwMGlLm5dOro44/z\nxddtAQHUa/3cc8CtW9IdhxPB1q1AixbUTL1LaYsgJITmTqanKxKe2nTrRv2WY8fK03fJdOi774Bh\nw2gDjLtkZt5OBI89RuclF7sCQBvYNG0qzraW1nAiWLECGDmy3M2lLQIXF6BHD9pElQGgevE5OcDC\nhUpHwjRHEKhZOWpUubtKVxU3a0Y1F44flz8+FTKZgK+/BmJipNunwdiJoKAA+OUXujq5R2bmXUvd\nu3YF9uyRNzYVc3MDVq6kQllcrprZ5fBh6uPo0qXcXWVWFXftqlwFNhXy9KTlTSNHAhcuiP/6xk4E\nv/5KNWAbNy53V2nXEED9IZwIymjZEli0iMrn8qpjZrN164AnnyydMmpx7Rpw8yaVQQfAiaACTzwB\nDB9OyUDsXjNjJ4J164DBg8vdLAj3tAgefJAGt65flzc+lXv6aRrX4/ECZrP16ys8506fBu6//678\nwImgQtOnU9L86CNxX9e4iaCoCPjpJ2DQoHJ3XbwI1K5duraM/hIYCPzvf/LGqAEff0ylmRYtUjoS\npnppabQ+4OGHy911+jTtD14qJARITaVfPVbK0i27aBHw22/iva5xE8Hu3dT34+dX7q6MjApK4YaE\ncCKoQI0awOrVwP/9H4+nsyqsWwcMHEg1FO5RLhFUr05zJ48ckS8+jfDxAZYvpymlYq27M24isNJE\nBehLWS4/cCKwqlUrqgrw9NPUpcZYhSzjAxUolwgAIDiYzzkrevem9Ty3yzQ5zZiJQBAqTQRpaZwI\n7PXEE8CECcBTTxm+NBOrSE4OkJxMq4YrUGEi4HOuUu+8Q7OJxGDMRJCcTMkgOLjCu9PTgebN77kx\nOJiex4tcrJoyhU7mceN48JjdY+NGWih2VxmXu1kGi8vgRFApFxfqIhLltcR5GY2xtAbumcJmkZZW\nQSLw8AAaNABOnZI+Po0ymYBly2htwbx5SkfDVGX9eqvdQkAlXUOHD/NVRSXq1BHndYyZCKxMG7Wo\nsGsI4D5LG9SuTRd/s2YBGzYoHQ1ThatXaUls374V3n3rFnDuXAUTNBo1ogJXp09LH6PBGS8RpKdT\n/Yju3Su8WxDoIRUmgpAQukJhlfLzoyQwdixw4IDS0TDFbdlCU0br1avw7jNngCZNgGrVKriTu4dk\nYbxEsGEDFdev8FtHawiqVaPNWMoJCqJa6axKDz1ES+IHDaIWFjOwDRsqXK9jUWG3kAWfc7IwZiKo\npFuowoFii/bt+Utph0GDgLfeAvr1A3JzlY6GKaKwENi0idYPWGG1BQ7wOScTYyWCixdpX+LwcKsP\nqXCg2MLfnzYRLSiQJj4deuUVmlrarx+Qn690NEx2f/xBJ1Rp4a7yUlOB1q2t3MmJQBbGSgS//AKE\nhQG1all9SKVXJ25u9I09dkyK6HRr9mygXTtqiPEaA4OpolsIAP7+u5JEEBAAnDhBLQsmGWMlgkoW\nkVmcOEGVNa3iKxS7mUzAF1/Q7NsRI6jMEzMAQbApEaSmUmO7QjVr0v4Eqanix8dKGScRFBQA27bR\nQHEljh8H2rSp5AGcCBzi6kp7AN24AURG8ro8Q0hOpv8GBVl9iCDQxZfVFgHA55wMjJMILHsPNGxY\n6cNSUqhfRkWcAAAV0klEQVQ1ahV/KR1WvTrw44/U/fbii5wMdM/SGrCycBOgzWjq1KliYRSfc5Iz\nTiJYv77KJuqlS1T1tnSXpIrwl9IptWrRJJLUVGoZFBcrHRGTjI3jA1a7hSzateNzTmKyJgI/Pz8E\nBwfDbDajc+fO8h24uNjq3gN3s3QLVXIBQyPJeXm8LZcT3N0pGWRkABERPGagS//8Q1Pwevas9GGV\nzhiyaN+ey1FLTNZEYDKZkJCQgMTEROzfv1++A+/eTZf5VueFkpSUKsYHAKr01LYtfzGdVLs25eZz\n52gfc04GOrNmDU3MsLJw08KmRNC6NU/blpjsXUOCEgWkVq2iYvlVOH68ivEBi/bt7wyEMYfVqkV1\nia5cofLVvBOojqxeDQwbVuXDbOoaql6dNr1ISREnNlaO7C2C3r17o1OnTliyZIk8By0upquTZ56p\n8qE2tQgAHicQUY0aNHxTrx7Qpw+vQNaF9HSaCmRl74G7JSfT6VSldu24FS6hytttItu1axeaNGmC\n8+fPIzw8HAEBAeh5Vx9idHR06d9DQ0MRGhrq/EF//53KGrZqVeVD7UoEGzc6HxsDQOv0YmKAyZOp\nS3nLlkoXojK1W7OGSk67uVX6sPx8IDu7inU7FnzxVSohIQEJCQmivqZJUKSvBpg6dSrc3d0xadIk\nCsRkkqbb6KWX6Jv2xhuVPuzGDdpyIC/P6t4Zd2Rl0dzo8+erGFlm9po9G1i4kAaT27VTOhrmkC5d\ngA8+qLSUCwDs2weMH0/7V1Rp/XqqYvjzz+LEqCNi/HbK1jV0/fp1XL16FQBw7do1xMfHI6iShSai\nKCwE1q61qa/y6FFqNFSZBADA25tWwpw753yMrIzJk4Hp04FevYDNm5WOhtktPZ02b+rVq8qHJidX\nutasLJ5CKinZuoZycnLw5O0dioqKivDcc8+hT58+0h502zb6dbdaPOiOpCSgQwcbX9dkulMe18vL\nqRBZeSNHAi1a0MbckycDEydyw0szVq2ibqEqZgsBdowPAPSFOHeONrkRa1suVkq2RNC8eXMkJSXJ\ndTgSG2vTIDFAe1+EhNjx2paZQ2FhjsXGKvXww8CePcCAAdRaW7SIJo8wFRMEGuz58kubHv6//1FV\nWpu4ugKBgfRl6NLF8RhZhfS7svjqVVrZ+OyzNj3crhYBwINXMnjgAWDXLroQDA2lnayYih08CNy8\naXX3v7uVlACHDgGdOtnx+nzOSUa/ieCHH4BHHwU8Pat8aEkJXZ1wIlCfOnVoi+kBA2jXs61blY6I\nWRUTA4webVM/3vHjQOPGVJHWZjxOIBn9JoJly4AxY2x66LFj9KVs1MiO17fMa+bKaZJzcaGdzr7/\nnkpSTJvG/+yqc/MmsHIlJQIbHDhAid0uXGpCMvpMBCdO0JJFGzsg9+4Funa18xgeHrQK6p9/7I+P\nOaRXL+p9iI8H+valWbxMJX75hS6OqijjYnHwoJ3dQgC3wiWkz0TwzTfAc89VuaDFwqFEAPDG2gpo\n2hTYvp0+L7OZuo2YCnzxBRAVZfPD9+51YMy3WTNahXbxop1PZFXRXyK4eZMWnowda/NT9uxxcCIC\n1xxShJsbMHUqJYH/9//o9+f2EhWmhNRUIDHRpvU6AH1WR4860DVkMnGpCYnoLxGsWUM/0G3b2vTw\n8+epHLJdA8UW3FRVVLduNNvLZKLP77fflI7IoD7/nMbjbFqNSRdeHTva/PCyeJxAEvpLBAsWAK+8\nYvPDExKovo0N61/K40SguDp1qAE4bx4NJL/wAm0wxGRy/TrNFnrpJZufsmNHldsUWMfnnCT0lQgO\nHKAqVlXsS3y37dttKpJYscBAGpQuLHTwBZhY+ven3wc3N/qt4LEDmcTGUr9qixY2P2X7dprZ7RCe\nQioJfSWCefOAceNoFaKNtm2zqSxKxWrVojKZJ044+AJMTHXrAp99Rr9Nb70FDBwInDypdFQ6VlIC\nfPwx8PrrNj/l4kX6HXe6RaBMrUzd0k8iOHWKqpTZ0UQ9fpwmIdhVWuJe3FRVnZ49aYHgww/Txeo7\n79Be1ExkGzfSvqN2NKnj46k14ND4AEC1vUwmICfHwRdgFdFPIpg1C3j5ZaB+fZufsnEjXTW6OPOv\nEBTEM4dU6L77gClTKCGcPk07z61cyReSohEE4KOP6B/ZjoqAmzbRGhCHWWYO8cWXqPSRCDIzaWu8\niRPtetrGjVXuZ181bhGomo8PsGIFdRfNnk0thO3blY5KBxISaFR+8GCbn3LzJq07c/qcCwmh6apM\nNPpIBLNn05SRxo1tfkpGBs1Cc3h8wKJ9e+DwYSdfhEmtRw+aSzBpEvUe9ulj44YorDxBAN59F3j7\nbbvG4+Lj6XTx8XHy+F270oo0JhrtJ4K0NODbb4E337TrabGxVO/e4b5KizZtgAsXeJMaDXBxoark\nR48CTz1FheyefprzuN1+/plaAyNH2vW0lSttrgpfuW7daDEC9/OJRvuJ4J13gFdfpV3DbCQIwPLl\nwKhRIhzf1ZWuUPbsEeHFmBzc3Gg4KTWVuooef5y6K/bvVzoyDSgupilZH35oV2vg4kXqFho+XIQY\n/PwojowMEV6MAVpPBAcPUl/l7X2PbZWQQDPfbCibbpvu3alwPtOU2rWpRMXJk9RVNHQobbP72298\nsWnV8uU0IWPAALuetmwZJduGDUWIwWSiVgF3D4lGu4mgpIQGh6OjaQqbHT79lJ7q1Gyhu3Ei0LSa\nNYF//5uWg4wYQUtRHnyQahfeuKF0dCqSm0utgblz7ZopVFJCVSjGjxcxFkv3EBOFdhPBV18BRUV2\nVTwEaCHw7t12d29WrnNnKnrDvxqaVr06EBlJYwgffUR92n5+wHvvcclrADQ4/NRT9H23w+rVNI/D\nzqdV7uGHgZ07RXxBYzMJgjoawSaTCTaHkp1N8/e3bQOCg+06zrBhNPvsnXccCLIyDz0EzJkDPPKI\nyC/MlHTsGJWvio2lGWZjxwKPPWZX97g+7N1LSeDoUbvW6hQV0UyhhQuB3r1FjOfWLcoup06J1N+k\nXXb9dlqhvRaBIFA7PirK7iSwbx+1Bv7zHwni6tMH2LJFghdmSgoMpLIVp0/TQqipU++0EtLSlI5O\nJteu0fTsTz+1KwkA1L3m5QWEhYkcU/XqdNG1bZvIL2xM2ksES5bQ6F50tF1PKyykPsoPPqASQaJ7\n/HEqccF0qW5dqmy6bx+tjr1yhRqBjzwCLF5M5cx1a9Ik6texcb8Bi+xs6k2aP9+uIQXbhYfT4gTm\nNG11DR05AoSGUt9gQIBdrz91KrVuN22S6EtZVERN1WPH7JrKyrTr5k1qBK5cSd+rrl1psHnwYNrF\nVBc2bKCZFUlJdr0pQQCGDKHT9MMPJYrt2DG6AEtPl+ik1gZjdQ1dukTz+2bOtDsJbNlCsxa++krC\n70u1atT+5e4hw7jvPqpV9f33VOXk+eep/HWzZtSNtHgxcOaM0lE64ehRagbFxtqd2T7+mKb5v/ee\nRLEB9DtgMnGtLxFoIxEUFVGzNDycpnXY4fBhWji2erUIS9urMmAA8OOPEh+EqVHt2rRYav16+vGP\njKTxqJAQ2o1r6lQqaVFcrHSkNsrNpSw3e7bdG3r/8APwySd0Kji9cr8yJhMtDV+9WsKDGIP6u4YE\ngYrDZGQAP/1k11ZiBw7Qb/OCBfR9kdyVK3Q5yDMZ2G1FRbTEZMMGIC6Oqif/6190TdO7t137ucjn\n2jWaGtW1K13a2yE2liZjbNniZHl3Wx08SBk4NdWw3UNidA2pOxEIAvVP7ttHg0J169r8ekuXUvmh\nr7+mCxvZDB9O4xgvvyzjQZlWZGbSRJdff6U/NWrQtNTu3emPv7/Cv2c3btAJ4+NDJ4+Nqy6LiqgO\nXWwsVfW1c0Kf4wQBaN0aWLWKml4GpO9EUFJCsxV27qQzxsZpa2lpdEWSmkr72AcGShSwNT//TKNj\nu3fLfGCmNYJA3fAJCfR12bWLNkp6+GH607Ur0KGD3TM2HZefT+Nw9erRwIeNiyUOHABefJHmSnz3\nnV1FgMXx/vtU+HHRIpkPrA76TQQ3btDI25kzdHnRoEGVzz96lFa+r11L9WMmTaLBPNkVFgItW9Ko\noUGvUJjjMjPvJIX9+2mMy8sLMJvv/AkJAZo2FbnlcP480K8fvfjixVV2wQoCXaN9/DH1zsyYQWNx\nirRmMjNpgWl6ul29Bnqhz0SQkUG1an19qcCVldGmwkKa0bZ1K135nz1L6wTGjwcaNZI5+Ht9/DFw\n6BBdVTHmhOJiat0mJtJXKjGRJskUFFAF9DZtaPJMQAD1kDRrBnh42PmDfOAADaKNHk2j2laeXFBA\nCcpyztWsST2gUVH0d0UNH04LO+wsQKkHuksEJ5dsg/DmFJSMiUJJ1AsoKnHBpUs0gSE3F/jnHzop\nUlPpZGjRgrrjhwyhjUdUs+z/8mU6KxMSgLZtlY6G6VBeHu25ffw4kJJC/01NpeuooiJKCJY/3t7U\nqL77T+3agJtrCar/8B3cFs+Hy/RpKHj0cRQU0A9+bi5daJ89e+d8O3mSWiTh4TQJo2NHFY3PJidT\nYCdO2F2EUut0lwj83DLg0rghXGrVhIsL/bDXr09fXA8P4P77gVat6E9wsMoX7cydS3si/vSTis4W\nZgRXrlBCyMigi6dz5+5cTOXm0t4A1/NuoPD0WRSWuOJWw6YocXVDzZq06r5mTTrfmjalPy1aUM9L\nYKDE00Gd9dxzQPPmVD7AQHSXCITiYhFrQyvs5k26ZJoyReRSp4w5IS+PJjMsXUqrvSZMUFFT2kln\nz9Loenw8/dcg9LeyWC9JAKCR6m+/pSlMJ08qHQ0zuvPnqf+/TRvquvzrL+C11/STBABqvsyaRRde\nV64oHY2m6OiXV4XMZmDaNFqck5OjdDTMaASBNm958UVaoJCZCezYAXz5JdCkidLRSSMiAujZkwa/\nb95UOhrN4EQgtZdfBiIikNCxI43q6VRCQoLSIUhKM++vpITmc0ZH04//88/T4FpKCiUAK3W6NPP+\nqmIyUSmBOnWoIF1enn7em4RkTQRxcXEICAhA69atMXPmTDkPrax336VE0LMn7dBRVKR0RKLT+8mm\n2vdXXEwzZpYsoYn83t40DfTqVZq+nJIC/Pe/tBihEqp9f46oVo1WGj/4IBAcjIQvv+RNqKsgWyIo\nLi7GhAkTEBcXh6NHjyI2NhbHjh2T6/DKM5upWb52LV2pzZ9PBdsZs0VJCXXt/PorXUz8+980d9rD\ng1YD//EHzaE+cIBWV86ZQ/PqjTpjzdWV/g2++Yb+zbp0oWXPV68qHZkq2V7BzUn79+9Hq1at4Ofn\nBwAYPnw4NmzYgEDZa0AoKDCQppTu3k0n83vv0TqD7t3ppA0MpOlvBpsHbViFhTSL5+65nZY/WVm0\nsv7MGZoHmplJc6kDA+lPQABtfNCxo00r7w0rLAwYN47mmy9bRitOu3QBunWjc65VKzrnFClDoB6y\nTR9ds2YNtmzZgiVLlgAAVqxYgX379mHBggUUiFGvXBhjzEnO/ozL1iKo6odeJcsZGGPMcGQbI/Dx\n8UFGRkbp/2dkZMDX11euwzPGGLNCtkTQqVMnpKamIj09Hbdu3cKqVaswUNaNAhhjjFVEtq6hatWq\nYeHChXjsscdQXFyMqKgoYw0UM8aYSsm6jqBv3744fvw4Tpw4gbfeeqv0dj2uL/Dz80NwcDDMZjM6\nd+4MAMjNzUV4eDj8/f3Rp08fXLp0SeEobRcZGQkvLy8EBQWV3lbZ+/noo4/QunVrBAQEID4+XomQ\n7VLR+4uOjoavry/MZjPMZjM2b95cep+W3l9GRgZ69eqFdu3aoX379pg/fz4A/Xx+1t6fXj6/Gzdu\noEuXLujQoQPatm1b+tsp6ucnKKyoqEho2bKlkJaWJty6dUsICQkRjh49qnRYTvPz8xMuXrxY5rbJ\nkycLM2fOFARBEGbMmCG8+eabSoTmkB07dgiHDh0S2rdvX3qbtfdz5MgRISQkRLh165aQlpYmtGzZ\nUiguLlYkbltV9P6io6OFOXPmlHus1t5fVlaWkJiYKAiCIFy9elXw9/cXjh49qpvPz9r708vnJwiC\ncO3aNUEQBKGwsFDo0qWLsHPnTlE/P8VLTNy9vsDNza10fYEeCPfMhNq4cSMiIiIAABEREVi/fr0S\nYTmkZ8+e8PDwKHObtfezYcMGjBgxAm5ubvDz80OrVq2wf/9+2WO2R0XvD6h4NpvW3p+3tzc63K7G\n6e7ujsDAQGRmZurm87P2/gB9fH4AUKtWLQDArVu3UFxcDA8PD1E/P8UTQWZmJpo1a1b6/76+vqUf\nopaZTCb07t0bnTp1Kl07kZOTA6/bS/29vLyQo/FCdNbez9mzZ8vMCNPyZ7pgwQKEhIQgKiqqtOmt\n5feXnp6OxMREdOnSRZefn+X9de3aFYB+Pr+SkhJ06NABXl5epd1gYn5+iicCvS4k27VrFxITE7F5\n82YsWrQIO3fuLHO/yWTS1Xuv6v1o8b2OGzcOaWlpSEpKQpMmTTCpkm0QtfD+8vPzMWTIEMybNw91\n6tQpc58ePr/8/HwMHToU8+bNg7u7u64+PxcXFyQlJeHMmTPYsWMHfvvttzL3O/v5KZ4I9Lq+oMnt\nMr+NGzfGk08+if3798PLywvZt+sLZWVlwdPTU8kQnWbt/dz7mZ45cwY+Pj6KxOgMT0/P0hNs7Nix\npc1rLb6/wsJCDBkyBKNGjcLgwYMB6Ovzs7y/kSNHlr4/PX1+FvXq1UO/fv3w559/ivr5KZ4I9Li+\n4Pr167h6u7jVtWvXEB8fj6CgIAwcOBAxMTEAgJiYmNIvrFZZez8DBw7EypUrcevWLaSlpSE1NbV0\n5pSWZGVllf593bp1pTOKtPb+BEFAVFQU2rZti4kTJ5berpfPz9r708vnd+HChdJurYKCAmzduhVm\ns1ncz0+yYW47bNq0SfD39xdatmwpfPjhh0qH47RTp04JISEhQkhIiNCuXbvS93Tx4kUhLCxMaN26\ntRAeHi7k5eUpHKnthg8fLjRp0kRwc3MTfH19haVLl1b6fqZPny60bNlSaNOmjRAXF6dg5La59/19\n/fXXwqhRo4SgoCAhODhYGDRokJCdnV36eC29v507dwomk0kICQkROnToIHTo0EHYvHmzbj6/it7f\npk2bdPP5HT58WDCbzUJISIgQFBQkzJo1SxCEyn9P7H1/qtmzmDHGmDIU7xpijDGmLE4EjDFmcJwI\nGGPM4DgRMFaFy5cvY/HixUqHwZhkOBEwVoW8vDx89tlnSofBmGQ4ETBWhSlTpuDkyZMwm8148803\nlQ6HMdHx9FHGqnD69Gn0798fycnJSofCmCS4RcBYFfhaiekdJwLGGDM4TgSMVaFOnTqltaMY0yNO\nBIxVoWHDhujevTuCgoJ4sJjpEg8WM8aYwXGLgDHGDI4TAWOMGRwnAsYYMzhOBIwxZnCcCBhjzOA4\nETDGmMFxImCMMYP7/xS1VsIEXY0VAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 4 }, { "cell_type": "markdown", "metadata": {}, "source": [ "\u5e95\u4e0b\u662f\u53e6\u4e00\u500b\u7c21\u55ae\u7684\u4f8b\u5b50" ] }, { "cell_type": "code", "collapsed": false, "input": [ "from scipy.integrate import odeint\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "%matplotlib inline\n", "pi = np.pi\n", "sqrt = np.sqrt\n", "cos = np.cos\n", "sin = np.sin\n", "\n", "def deriv_z(z, phi):\n", " u, udot = z\n", " return [udot, -u + sqrt(u)]\n", "\n", "phi = np.linspace(0, 7.0*pi, 2000)\n", "zinit = [1.49907, 0]\n", "z = odeint(deriv_z, zinit, phi)\n", "u, udot = z.T\n", "# plt.plot(phi, u)\n", "fig, ax = plt.subplots()\n", "ax.plot(1/u*cos(phi), 1/u*sin(phi))\n", "ax.set_aspect('equal')\n", "plt.grid(True)\n", "plt.show()" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAALsAAAD9CAYAAADkkQYbAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXlcVHX3xz+AlJrmLploFIgIKriiiDZmuOCeVmqaZvaU\n5Vaatpc9aWiZmkumP/fnydTUtEcxKZnE3QLcUMSFBARcSRCRZc7vj+MgywzcucvcO3Dfr9e8dObe\nOffc4dzvPfd8z/ccJyIi6OhUApzVVkBHx17oxq5TadCNXafSoBu7TqVBN3adSoNu7DqVBknGnpOT\ng8DAQAQEBMDX1xfvv/++XHrp6MiOk9Q4e3Z2NqpXr478/HwEBwfj66+/RnBwsFz66ejIhmQ3pnr1\n6gCA3NxcFBQUoG7dupKV0tFRgipSBZhMJrRt2xYXLlzA+PHj4evrW7jNyclJqngdnVKIdUYkj+zO\nzs6IjY1FcnIy9u3bB6PRWEoxOV6ffvqppuRUdJ20en6SbFXSt4tQq1Yt9O3bF3/++adcIouRmJio\nKTlyytKiTnLKklMnKUgy9uvXryMjIwMAcPfuXURERKBNmzayKKajIzeSfPbU1FSMHj0aJpMJJpMJ\no0aNQo8ePeTSrRhjxozRlBw5ZWlRJzllyamTFCSHHssU7uQk2c/S0SmKFJtymBnUkg++asuRU5YW\ndZJTlpw6ScFhjF1HRyq6G6PjUFQKN0ZHRyoOY+xa9B8rsk5yytJ9dh0dO6P77DoOhe6z6+gIwGGM\nXYv+Y0XWSU5Zus+uo2NndJ9dx6HQfXYdHQE4jLFr0X+syDrJKUv32XV07Izus+s4FLrPrqMjAIcx\ndi36jxVZJzll6T67jo6d0X12HYdC99l1dATgMMauRf+xIuskpyzdZ9fRsTO6z67jUOg+u46OABzG\n2LXoP1ZkneSUpfvsOjp2RvfZNcadO0BcHHD6NPD330ByMpCSAqSlAVlZvP3OHcBkAlxdgSpVgIcf\nBurVAxo0ABo2BBo3Bry9gWbN+N/HHlP7rORDik3pxq4iBQXAqVPA/v38OnIESE0FfHwAPz/gqafY\ncN3dATc3oGZN4JFH+OXsDOTlAfn5QE4OcOMGcO0av5KSgIQE4Nw54OxZoGpVoGNHIDAQ6NIF6NSJ\nLxRHRJJNkQQuX75MBoOBfH19yc/PjxYuXFhsu0TxxYiMjNSUHLGybt0i+vFHopEjierVI2renCg0\nNJJWrSKKiyPKy5NXJ5OJ6OJFog0biKZMIWrThqh2baIhQ4hWrCC6dk24LLl0koIUm5JUstrV1RXz\n589HQEAAsrKy0K5dO4SEhKBFixZSxFY47t0Ddu4E1q0D9u4FgoOB/v2B2bOBJk0AoxEwGJQ5tpMT\n8OST/Bo2jD9LSwP27GGdpk4FunYFXnoJGDgQuN8iq2Ii2yVHRAMHDqTffvut8L3M4h2OCxeI3nmH\nR3CDgWjlSqKMDLW1Ks7t20Tr1hH16kVUty7R1KlEly6prZV1pNiU5AZiZhITExETE4PAwMBin48Z\nMwYeHh4AgNq1ayMgIACG+8OYOSRV0d4/9JABc+cCkZFG9OkDHDtmwJNP8vaYGPX1K/l+1CgDRo0C\nNmww4uefgXbtDHj6aaBPHyOaNVNXv9jY2MLuLpLb1chxtWVmZlK7du1o27ZtxT6XSTwRadN/LCnr\n0CGinj2JnniCaMkSoqws9XUSQ2Ym0YIFRHXrRtKQIUSnTqmvkxkpNiU5zp6Xl4chQ4Zg5MiRGDRo\nkFRxDsn58+zvvvgiMGQIR0HefJOjJo5IjRrA5MnAf//LkZtnngH+9S/g+nW1NZOIlKvMZDLRqFGj\naMqUKRa3SxSveW7fJnr3XfbJ58whyslRWyNluHWLaNIkogYNiJYuJcrPV08XKTYlyRqjoqLIycmJ\n/P39KSAggAICAig8PFwWxbTOrl1ETZoQjR5NlJqqtjb24fhxoi5diIKD+eFbDVQz9nKFV0Cf/dYt\nNnAPD6KvvpImy4zWfPayZBUUEM2bR1S/PtHy5RzHt6dOUmxKz42xgaNHgbZtgWrVgJMngfbt1dbI\n/jg7A++8w3MDS5fyc0pWltpaCUNPFxAAETB/PhAWBixbBjz3nNoaaYOcHH4QP3oU2LaNc3GURs+N\nUZB794DXXuPkrJ9+Au5PGejchwhYvhz4+GPgxx85cqMklWLxhhq51Vev8h/v7l1g377Shq7FfG97\ny3JyAl5/Hdi8mdMRNm5UVicpOIyx25uEBM4S7NGD/4AVOmdEBp5+GoiI4FybRYvU1sYyuhtjgdOn\ngZ49gZkzgXHj1NbGsUhM5AFi8mRg0iT55UuxKdlyYyoK0dFAaCjwzTfAiBFqa+N4eHgAv//OWZwP\nP8wujlZwGDfGHv7xiRNAnz7Ad98JM3TdZ7eM2eC/+IJTDnSfXWNcusQj+qJFwODBamvj+Hh6Art3\nA2+/zYOIFtB9dgDp6byg4u23OW6sIx8REcCoUUBUlDxxeD3OLoGcHI4k9OoFfP652tpUTL7/Hliw\nADh2jDMqpaDH2UXKIeKRvGlTjrxoQaeKKKt5cyM6deLfWs2xz2GMXQmWLgX+/BNYvZonR3SUY/Fi\n4K+/+LdWi0rrxhw7BvTrBxw8yA9TOsoTFwd068a5NE89JU5GpXBj5OTOHV5Nv3ixbuj2xNcXeO89\nnqhTYwx0GGOX0z+eNg3o3Bl4/nnt6CQXWpRVVM7bb/Ngs2KFLKJtotLNoB47xvHf2Fi1NamcuLgA\nK1cC3bvzfEaDBvY7dqXy2e/eBVq2BJYsAXr3Vlubys3kyVy+b+lS276nx9kF8sknwJkznI6qoy43\nb3JNy717eQASSqV4QJXqPyYk8CjywgvS5BRF99nFy6lbF/jwQ35gtRcOY+xS+egjfjiyp4+oUzav\nvw7ExHCmqT2oFG5MdDTH1BMSHLdwUUVl4ULgjz+ArVuF7a/77OXQuzcwYICe5KVFsrN5ruP33zkO\nXx66z14Gx47xzJ15xZGj+LQVSVZZcqpX5wXttkZlxOAwxi6Wr79mX/2hh9TWRMcar78O/PADkJmp\n7HEqtBtz6RLQoQP/W7OmamroCGDoUF67On582ftVCjdGDEuWAGPH6obuCIwdC6xfr+wxHGZkNxqN\nhUXqhZCbyy1cDhwAvLzEywH4IeryZV45n5wMZGTwKy7OiEaNHsh65BGgTh1+1a/PrV28vIBatcqW\nL0YnR5IlRE5eHvD44/yMVVYhKtWqC4wdOxY7d+5Ew4YNcfLkSSmiZOd//+MZuqKGLoSbN3kJWXQ0\nx4BjYrgDXdOm/Edwd+cJkdq1uZOdjw9/j4gTnG7d4ovi6lV2n86f5251rVtzx7qOHYGgoIrVrlEO\nXF25tv2GDcD77ytzDEkje1RUFGrUqIGXX37ZorGr6bP3789+4OjRZe+Xn89x3vBwnro+f54zIjt0\nAAICgDZteIR2FunwEfEa19hYzuM+epTvNh4eXMmgXz8u+C9WfkXi9995VvXwYev7qBpnT0xMRP/+\n/TVl7P/8wy5McjLw6KOltxcU8A+7eTPw889seAMGcKm7Dh2Uj9zk53PP0/BwPr45v37UKKB5c2WP\nrWXu3eOmxRcvchNjS2i6SJJcDcSKxmrL2z811YCuXYHo6OLbt2wxYtWqWJw+PQX16wMdOxqxaBEw\nbNiD7x88KLzB1YIFC0SfT5cuQF6eET16ABcuAGfPGtCpkxHe3sDnnxvQsyfwxx/C5Znfx8bGYsqU\nKTbrI/f5FX1v/qy8/Q8dMsLPD4iIMGDYMA02ELt06RK1bNnS4jYZxBdiS0H7558n+r//e/D+5Emi\n4cO52W2/fpH011/210mInLt3uX1kq1ZE/v5E//ufsGL/Sugkpyxb5CxeTPTKK9a3S7EphzF2oeTn\ns1GnphIdO0Y0aBCRmxv3PPrnH7urIwqTiWjbNiI/P6KgIKIjR9TWyH7ExBC1aGF9uxSbqnCPRbGx\n7PtNngwMGsR++MWLwPTplv13LeLkxLofP85d6gYO5MmWW7fU1kx5WrZ8EN6VG0nGPnz4cAQFBeHc\nuXNo0qQJVitYJ0FInsadOxzaM69IOncOmDixeLlpR8odcXHhaFJcHF8ArVoBv/1mH53klGWLnCpV\ngHbt+AFebiQ9oG7YsEEuPSSzZw/wxhuAycTVpyZPVv6YubncW+ncOX6lp3NM/vp1vuDMPPQQRxfq\n1QMaNQK8vfnl48P9mcqjTh1OlBo8mI1/xAhg1qyKm+/j78+/a69e8sp1mBlUa2Rmch3wyEiuvvva\nazwp9OST8h8rI4Mvqv37gUOHeMR94okHxvv44zxzWq9e8bz5nBzgxg1+paRwXn18PLtXfn48ydSj\nB/Dss+Ub//XrwJgxfN5btvDxKhrff8/zEStXlt4myaZEe/sCUFg8HT1K5OlJNG4ctyBPS+OHU1sj\nGGVx5Qq3Qnz6aaKaNYlCQ4nmziWKiiLKzpYmOzubaN8+oi+/ZPmPPsoP1Fu3EuXmWv9eQQHRjBlE\nTz1FFBcnTQct8scfRJ07W94mxaYcxtiLhq9MJjbABg2INm9+sM+vvxIZDMLlWCMvj2jTJqI+ffji\nGTOGw4B37tguSwiRkZFUUMCGHxpKxPOu/PLzI2rblqhdO6JOnfiieP55ogkTiHx8eJ/Nm/kCkFMn\nOWXZKufKFf7bWkKKTTlc3ZicHHZVTp/mpKEnnniw7fx5aTOQmZl861ywgHNhzI2xlFrKd+4csH07\nsGkTuzb16rG/Om0aPw9s3Mjn2bUrL0yuU4fP/+pVfj54+GHg7NkHxZ769OGZ45o1uV+ro9avbNiQ\nI0+5uTI/l4i+TAQgt/grV4g6diR64YXSoywR0dSpHE+3lexsdk3q12fZSsa109NZxxYtiBo1Inrj\nDaKff+bPLZGcTDRxIuu2cCHfdUqyYweRkxPR55/zvs2acQfuadOI/vxTuXNRksaNif7+u/TnUmzK\nYYw9IYHoySeJZs607pMPHlzcrSmPggKi1auJmjQheu45ojNnZFHVImfPEo0cSVSrFrtFBw48cD2E\ncPo0UY8eRK1bE504UXr7unVETZsSpaTw73P8ONHHH7PRt2/P5yn1GcOetGtHdPhw6c8rvLEfP05U\nr14kLVtW9n5t2/JDa1mY/ce4OKLgYKLAQKKDB8XpJcQXvXyZ6KWX2AedNYvo1i1xcojYiFev5lF+\n/vzSF/2sWUTe3pF09+6Dz/Lz+XmjTx+eSZ47lygrS9DhVPPZiYieeYYoIqL051JsSvMzqNHRQEgI\n8NZb5Xdeu3nTeracmfx8bjzQrRs3qT1wgFN65SYvD5g7l1OEPT35eeKDDzgPXixOThx2PHyY876H\nDy8ez3//fc6Tnzr1wWcuLkDfvsCuXdzy5ehR1mfJEv4ttMojj/AkoayIvkwEIFX8yZM8Gm3bJmz/\nWrWIbt60vv3iRR7J+/QhSkqSpFqZJCSw69CrF9H588oc4+5dTm7r2LG4v5+RweHYTZusfzc2lqNW\nrVtzBEiLDB9O9J//lP5cik1p1tjPnSN6/HGiH34Qtn9BAZGzM9+2LfHTT+xKfPONbb6yrWzaxG7G\nt9/KG++3hMlE9NFH/LCbmvrg80OHiB57jOj69bK/u3Ejkbs70VtvCXdt7MXIkURr15b+vMIZ+9Wr\nPGGyYsWDz8rz+7KyiKpVK/25yUT02Wf8EHrsmHJxaJOJfeYmTcjmFGKpOs2cyTH39PQHsqZMIRo1\nqvzv3rzJhtWsWekHQjV9diVGds357PfucQ7IsGG2t1IvGVfOyWE54eHsq7ZvL5+eRSHilIWffuI0\ngrZtlTmONT75hNdvDhjAvx/ADXcjI8tPqKpTh1f1z57NSxm//17dJl9mCgr4eUNWRF8mArBVvMnE\no8zQoba7GllZRNWrP3ifmUnUvTvHzYtGJ+TGZOLYdseO7C+rhclENGIEz66a3adVq4i6dhXuTsXH\nE/n6Er36KlFOjnK6CmHQIKItW0p/LsVkNTWyL1nCM4Zr14pbgGwekW7d4giOpydXmqpaVV49i/LF\nF9yE7Ndfyy+ZoSROTsCqVVzRYPFi/uzllzlCtWePMBne3hzpuXGDF4IrXaGrLG7dkha5soiEi69c\nbBH/11/8YJeQYHl7eX5ffj4/oN68yRMSb79teUST02f/5JNIatq0+MOhGOTU6T//iaT69Xlugoho\nzRqikBDbZOTnE/3rX0TNm0fS1avSdRJzfq1bE0VHl/5cislqYmTPzARefBFYtMj2Oi9mXFw48f/p\np4HAQGDePGVzQ86d43LL27drqwZM48ZAWBg/7xQUcCz+9GngxAnhMlxcgGXLuNLC009zjr69uXmT\n6/PIiujLRABCxb/2GvuJUsjL4wzAzp2VDS0SEd27x3ePJUuUPY5YCgqIunShwhnnzz4jmjRJnKwP\nPyRq08byzK9S5OYSPfSQ5TRnKSarurHv3cuxXqmLod95h419/35pcoQwezZPTCkdR5fC8eNEDRvy\ng/q5czw5ZymJrDxMJqLJk3kQsZR8pwQXLnCejyWkGLuqbkx29oPa3OUthi5rHeOaNcAvv3C7wZQU\n8XKEcPkyu0hLljyo6yIVJdaNtm7Ni82XLAGaNeOU5chI22U5OQHz5/PD/tix4sKStp5fYqIyK81U\nNfYvvmC/sH9/8TKOHePKAdu380Ldixfl088SH33EHTyU+GPIzccf84V59y7H4HfvFifHyYmb9CYm\n8tpXpblwQaHfV/Q9QQBlib90iahuXc7XFsvt25wHYk7rXbKE/X+lOH+eqF49dePpttKzJ9H69ZzZ\n6e8vTdaVK+xy7tghj27WePNNzuq0hBSTVW1k/+ADYMIEjh6IZdIkwGDgAqYAL14+dUoW9SwyZw6P\n6mrG023l9dd5VrRDB77rSak906gRr5567bXy3UUpREcrNAst+jIRgDXxR45wkldmpnBZJWO1P/3E\n+RxFZWRkED3yiPVkMEtyhPLPP5xVWTTDUM3cEaGycnOJ6tThO2inTryYWapen3/O+eZCo162nF9e\nHs+EWwtYSDFZVUb2zz5jf7JGDXHf/+cfHtVXry4uo1YtHn3i42VRsxibN/NdpGFD+WUriasrEBoK\n7NjBD61ylNH/4APOwVmyRLqsksTF8d1ekeptoi8TAVgSHxPDo7qU3Iu33rLum7/0UvFsSbno3p1L\nXDgiGzYQ9e/PFRnExttLcuYMP7/IvS5g/nyevbWGFJO1+8geFga88w6vjBfDsWNcHCgszPL2Z57h\n2utykpXFx+3ZU1659iI4mLMxH3uMqxLIgY8Prx6bOFEeeWZ++42LRSmC6MtEACXFnz/P+S+3b9su\nKzIykkwmom7dipejLsmlSzyZYm3CR4x/vGMH+6hyyLKEPWq9NGlCtHw5152RKsvM3bu87uC336TJ\nMZOby4Wirl2zvo8Uk7XryP7998Arr4jvXvfrr1wzpazWMR4e7LtHR4s7hiX272d/3ZHx8wOSkngi\nTy6qVuU8+OnTucamVPbv5wkwxUr6ib5MiCg8PJyaN29OXl5eFBYWVuZVmJPDI661rMbyKCggCgiw\nnONckhkziN5/X9xxLNGrl/KxZaV5801eudS2rbxyCwqIOnQQvnyyLMaP51KAZSHFZEV/Mz8/nzw9\nPenSpUuUm5tL/v7+FFei8GBRxX78keueiGXLFl7ELCQf5dgxDkvKlbvSqJHlgj2OxNy57ALKbexE\nXPLCx0daAl5eHg+G5S1Ql2Lsot2Yo0ePwsvLCx4eHnB1dcWwYcOwfft2q/v/5z9cBkIMRMBHHxnx\n3nvC0nbbteP01r/+Kr3N1jyN3FyunOvuLl2WNexRU71WLc7rsWVBhFC9evTg0he//CJezr59/Bt7\negrXz1ZE13pMSUlBkyZNCt+7u7vjiIUFj2PGjIGbmwciIoAuXWrDaLS9IVWVKgbcvg3UqWOE0Sis\ngdW4ccBnnxkxbVrpBlu2HP/qVaBhQwOcnUtvj42NFaS/kAZbUr4v5PweeQS4fNl4vzamMHm2nN/0\n6cAHHxjx6KNA9+62n9+qVUCnTqX/vppoIPbTTz/RuHHjCt+vX7+eJkyYYPGWExvLfrRYhgzhxlK2\nYC5fLTWPJSaGG3o5OgcOcC6SXHH2kuTnc3nC8iqyWeLaNZ6dvnGj/H0lmKx4N6Zx48ZISkoqfJ+U\nlAR3S/d6cGVaa3Hx8rh6lWOvI0fa9j03N16HunatuOOayc+vGB0ugoK4AoFSfVZdXIBXX7XcQKA8\n1q7lvlGyr0wqidirJC8vj5566im6dOkS3bt3r9wHVLF8/TXR6NHiYtFHjnB8+d69B5/ZKufQIa4c\nYAlHirMTcW3L33+XR5YlkpM5D6dkwaWy5OTlceaq0HqbUmxK9MhepUoVLF68GL169YKvry9efPFF\ntGjRQr6r8D6rV/OiATF07Ai0aCFtdK9eXYGagypw7x4QE6Nc7RyAc1o6dwa2bhX+nU2bOJ9JiXqb\npRB9mQhAqvi4OK7TLSWkFRXFvqTYXJzLlzmXx9E5cIDXkirN2rVEAwcK27eggDuLhIcLly/FpjRR\nXcAaW7dydTAxNWTMBAfz7OGiReK+X78+11EpKBCvgxbYtUvBnJMi9OsH7N0rbKZ2+3aehZW7K55V\nRF8mApAqvm1bIrO7J8WvjY/nDL20NHFymjThCsAlcSSf3cfHcnF/MbLKo0eP4hmiluTk5hI1b060\nc6dtsqXYlGZH9qtXeS1icLB0Wd7enJMzfbr47589K10PtYiJ4ZG2Qwf7HK9fv/LXu65Ywf2f+vSx\nj04AtDuyb9zIOdhykZnJvvsvv9j+3enTuVKuozJuHFcYthfR0TxqWyMjg0t7xMbaLluKTWl2ZN+7\nl3PT5aJGDY7svP46V5uyhS5dOCPPEbl+nasLv/qq/Y7ZujWQlmY9d/799zmu7u9vP50AlUtplMW+\nfVx6zYwc+SNPPw0EBRkxdqxtKanBwVz6uWhLF7l0klOOJVlhYdz+3c1NuiyhuLjwb7ZvX2k5+/fz\ng+mcOaJES0KTxp6VxTVKWraUX/a//sUjzty5wr9Tty6vdo+IkF8fJUlJ4bvZhx/a/9gdOgD3U2sK\nMfew/fZbBSr0CkG0AyQAseL37+ccaaVISuI2LLt3C//OokXCOlloieeflzev3xa2bCn9zDVxYvH6\n8WKQYrKaNPZFi8pedCsHUVHcY0loS5j0dE4sK6tBmZbYvp3Iy0u93qfx8dyD1cyOHVy/UervJ8XY\nNenGnDnDE0FFkds/Dg7mssz9+3OIszwaNuQWi2vWKKeTXLLS0oDx44Hly4Fq1dTRy9OTw8dZWcDm\nzUa89hrw3/9yWxu10KSxX7yobBK/meeeAz79lAuiJiSUv//EiVyTPTdXed3Ekp/Pte5fe43PSy1c\nXDiOHh/P9THfflueORNJSLuplI1Y8d7enBdjL1as4BwcIcfs3Zto6VLldRKDycTrOHv3Vr5GvRCe\neYZnrkeMkG+JpBST1Zyxm0xciN7evua6dbwGsryyEMeO2V66z1589BGnWEitdS8XvKBS3r+lFGPX\nnBtz+zYXUCrpayrtH48axUU7R4xgX9ca7dtzQtWnn2rHZyfikoKbNwMff2yUrXScFL2WLuV/334b\nOHJEvBw50ZyxK9JLRyAGA096zJ/PneZu37a839dfc+9QIX6+0uTnA2+8wYud//hDpfh1CX74gevJ\njB0r7QFZduS7wZRGjPg//7RP3nVZZGVxLcmyVtCsX8+ZhGq2QU9L4wzDnj3FVVlTgjVruPTIiRNE\nCxcSlViWLBkpJqu5kf2ff9Svf/7II+zKzJnD6zbfeKN0XfORI7kr34QJ6ui4ezfP6gYFca662Cpr\ncrJsGUde9u4FWrXi31FLq7w0Z+wFBdzisSRq+MdDhnAJ5SpVeHnft9/ylLeZF1804uhRDkfaS6f0\ndOCll7gpwrp1wOefF297bo8aNCUhYrdlzhzAaOSip2LkKI3mjN1kUrZ/qa3Urs0do3fv5ioHzZrx\nw9edO+yP7tzJeTa2rLsUQ0YGj5q+vrzW8+RJLk6kNvfucfGrrVv5eafo/AiRtv6WmvPZw8O5tqJW\nOXKE11jWrUs0ZQrRqVOcctCggW25NkKJj+eckjp1iF55hSgxUf5jiCUtjfutDhliuW3k8uXS+9uW\nRIrJam5kd3bW9nrPjh2Bn3/mKsFVqwK9e3Pkpk0b/r/UEZ6Iu2d//TU/E3Ttyv748ePAqlW4X9FL\nfX77jZ8ZnnmGKwRUr156nzt32G/XCpoz9po1ub17SbQS0zbzxBNAr15G/P03P8ya05GHDOFb97Rp\nXMLj8GHg0qXSD2p5eex/nzoF/PvfRsyezdP8jRuze3LuHLfOTE7mdoxFKg3a5fysycrL4zYzo0c/\neGawtiA+I4PdQK347KJrPSpFrVockXEUnJ05IhIUxD1HDxzgHJB584AXXuC+Q1evsmHn5fH+zs4c\nH69bl6sX1K7N3+nblxdbeHhozNe9z7FjnHPz+OO8rrW8/lIZGdq5EwHQns+enMxxWkfm5k1+7ujc\nmVuTE3EaRG4uVye7e1cbuStCycwkeucdXje6fr3wPJeRIznuLidSTFZzbkzdujyLKkcnB7WoU4dj\n388/z373ihXsi7u6ct3IqlWl1cKxF/n57KJ5e/Na1pMneX5B6F3n77+1NbJr7ievVo3bApZcrKs1\nn708Wc7OnBeydy8vjQsKslwv3p46CaWggBdpe3oasWEDt5VcuxZo0MA2OYmJbOxa8dk1Z+wA+6x/\n/622FvLQqhXHn19/neupDB/OE1Va5N49rsLr68vPHOPH88Uqpj6k+QHcSmFndZDRnSqFWPFDh3Jb\nmopGZiZRWBinEj//PJHRKF+etxTOniWaOpXnCnr14ipsUvU6fZpzi+RGismKHtk3b94MPz8/uLi4\nIFrO1nQAvLyU6VKtNjVqADNmAOfPc/x8/HgOWc6bZ/87WWIiHzcoiEuMVKkCHDzIM8UGg/Ro0PHj\nQECAHJrKh2hjb9WqFbZt24Zu3brJqQ8A/pFiYop/5mg+e1nUrMlL/E6f5tBkZKQR7duzuzBzJueX\nlKxRI1XvK15MAAAPY0lEQVSnW7d4MmzyZC5i1KEDl/T79FPutRQWxoOMEFlCiI19YOxa8dlFx9l9\nSmb7yEibNjwCVnScnB70VzUXFfr1V66YdfIkj/otW/Lic29vjm83asQPikWTvwCOXt29C1y5wnnt\nSUmcb3/iBI+y167xKN69O0eH2rWznHAnFzEx8ne/lorik0pjxoyBh4cHAKB27doICCi/gVi3bgbc\nuAHs2MGrbgwGAwwGgywNtooiVZ75M6n6mN87OxvRpw8wZ44BWVnAypVGJCYCFy8a8OuvQEKCETdu\nALdvG1ClClClihEuLkBBgQE5OYCrqxF163KyWuPGvJ1b/Bjg5QVERfHxAgOVPb8uXQw4fBiYMEF4\nwzdr7+VsIOZ03+m3SEhICNLS0kp9Pnv2bPTv3x8A0L17d8ybNw9t27YtLdzJCWWIL5Pu3bnqrl2r\nvDoIRBztyMnhf6tV01bs/vBhjj4dPy6/bCk2VebPExERgZMnT5Z6mQ1dSQwGIDLywfuK5LNLlePk\nxJNTjz4K1KvHSVhmQ9fC+fFoLl2O3MgyFoi90sqie/fixq7jOOzZI28FZrko040pi23btmHSpEm4\nfv06atWqhTZt2iA8PLy4cAm3nHv3ONHowgVOltJxDG7cAJ58kktWW0r7lYoUmxJt7IKES1AM4HTZ\n/v3Ft4HXsT9r1nClgy1blJGvmM+uNoMHc2wY0H12R5G1dSuXFZQqRwk0bex9+7Lfbq1+i462uHqV\n5wr69VNbE8to2o0BeHTv18++bVJ0xDFvHk+GFa10LDcV1o0B2F9X8sfTkQci4P/+Dxg3Tm1NrKN5\nYw8N5fWY69cbZZGntk+rpBw1ZUVFscF36aKsTlLQvLG7uvK6R6We7nXk4auveLGKFtfOmtG8zw4A\nqam8oODCBfWKnupY5+RJoGdPrqJQtaqyx6rQPjvAmX4DBgDff6+2JjqWmDuXU4eVNnSpOISxA0D3\n7kYsWGC5powtaNE/1qJOQmWdPcsLPsaPt49OUnAYY/fwAEJCpBcR1ZGXGTP4pXblZSE4hM9uJiGB\nFyDEx+u+uxb44w8ODZ85Yz8XpsLmxljirbf4iX/xYlnF6thIQQHQqRPwzjtcMcFeVPgHVOCB3/f5\n59w76MQJaXLk1Ekrcuwp67vveDR/8UX76iQFhzF2M/Xq8aLkCRMcu2qYI5OUxA3Lli/XzuooITic\nGwPwLbRrV+5sp1abl8oKEaddt2/PBm9vKpXPbiY+nqemjxyxTzdsHea777hq2MGDvDTQ3lQqn91M\n8+YP6oTn54uXI6dOastRWlZcHPDJJ8B//2uboes+uwxMmcJVtj76SG1NKj6ZmVyVOCyMBxpHxGHd\nGDPXrnHBn6VLtbtowNEhAoYO5bmN5cvVTfaqFG6MNRo0AH78kbspnz2rtjYVk1mzuNLY4sXazmos\nD4cx9rL8vqAg7sHZty+P9GLlyKmTGnKUkLV2LS/K2LIFePhh9XWSguZ6KonllVc4BXjAAK4pXq2a\n2ho5PkeOAN98w0WPHn9cbW2k4/A+e1FMJo7OXL/OVQnEjkQ6nPcydCiwfTvfObVCpfbZi+LszC1d\natTgyEFurtoaOSaRkWzoGzdqy9Cl4jDGLtTvq1IF+OEHfpB64QUu/ilGjpw62UuOHLIiIvh327yZ\nqwprQSe5cBhjtwVXV+66bO5A7Uh9VdVk/Xruhrd1a/HCpBWFCuWzl8Rk4uViUVHAzp1cs1ynNETA\nl19yDH3XLl7vq1V0n90Kzs7At98Cw4YBHTtyPodOcbKy+PfZupV/Hy0bulQcxtjF+n1OTsB773Fr\nlcGDgalTjZDrZuPoPnt8PDclrlGD21eWDC9q8fykINrY3333XbRo0QL+/v547rnn8I/GHePQUP6D\nbtvGCw5u3lRbI/UgYpelSxfOL1q5UvuVAWRBbE/JPXv2UEFBARERzZgxg2bMmFFqHwniFSM7m2jy\nZCJ3d6KICLW1sT+pqUR9+xK1bcu9Sh0NKTYlemQPCQmB8/1lKoGBgUhOTpbp8lOWatWABQuAVat4\n1nXs2PJTDCoCBQXchrJVK+5GeOhQxfbPLSFLusCqVasw3MqqWzHd8iy9L+r3Se2+NmXKFISEAMuW\nGbFmDeDnZ8AXXwBeXkY4OwuXt2DBAtHnU/R9yXOU4/yKbn/4YQMmTADy84346itgzBjHOT85u+WV\neU949tlnqWXLlqVeO3bsKNzniy++oOeee072W05JIiMjFZMTG0sUHEzk60u0davwVuZK6iSHrBMn\niAYMYJdt3TrbW7Rr8fyk2JQka1y9ejUFBQXR3bt3LQvXoM9uDZOJaOdOIn9/og4diHbsILr/SOJw\nREcTDR9O1LAh0TffEFn58zgkqhh7eHg4+fr60rVr16wLdyBjN1NQQLRxI1G7dkTe3kRLlxJlZamt\nVfnk5xP98gtR9+48ks+dS3T7ttpayY8qxu7l5UVNmzalgIAACggIoPHjx8uqWEnsfUs1mYj27SMa\nOJCoTh2iV1/l90VdAS3c5s+eJXrvPaLGjfmO9OGHkZSbK4tamji/kkixKdEPqAkJCdIeFjSOkxOX\n6+jaFUhJ4UXGb7wBZGcDgwbxEkA16tYQAdHRwI4d/EpNBUaNAn79FfDz49xzV1f76+UIVOjcGLkh\n4kpkv/wC/O9/vAywWzeenAkK4loqci8aycsDTp3iCbGoKH7VrAkMHMgLVTp35kzPykKlrBujBdLT\neZHDwYP8OnUKaNoUaNGCX56eXFvezY1fNWtyCYqHH+Y7R14e3ymys7lZbkoKr/VMTuZioadOAefP\ncwXj4GC+ywQHA089pfaZq0elMHaj0VgYf9WCHEuycnLYOM+c4VH/wgW+INLTuePznTvcudu8qKRK\nFe4C7eJihJubAY0bo/Dl48NuiY+PbZ2ilTw/teUA0myqEt0AladqVaBlS36VBRH7+y4u/N5orJj5\n41rDYUZ2HR1Az2fX0RGEwxi7FnOrK7JOcspy+Hx2HR1HQ/fZdRwK3WfX0RGAwxi7Fv3HiqyTnLJ0\nn11Hx87oPruOQ6H77Do6AnAYY9ei/1iRdZJTlu6z6+jYGd1n13EodJ9dR0cADmPsWvQfK7JOcsrS\nfXYdHTuj++w6DoXus+voCMBhjF2L/mNF1klOWbrPrqNjZ3SfXceh0H12HR0BOIyxa9F/rMg6ySlL\n99ltJDY2VlNy5JSlRZ3klCWnTlIQbewff/wx/P39ERAQgB49eiApKUlOvUph7r6gFTlyytKiTnLK\nklMnKYg29unTp+P48eOIjY3FoEGDMHPmTDn10tGRHdHGXrNmzcL/Z2VloX79+rIoZA3J/XRkliOn\nLC3qJKcsOXWSgqTQ44cffoj169ejevXqOHz4MGrXrl1cuJOTZAV1dEqiSBXfkJAQpKWllfp89uzZ\n6N+/f+H7sLAwxMfHY/Xq1aKU0NGxB7JMKl2+fBmhoaE4deqUHDrp6CiCaJ+9aJuZ7du3o02bNrIo\npKOjFKJH9qFDhyI+Ph4uLi7w9PTEd999h4YNG8qtn46OfIhuPWaBadOmkY+PD7Vu3ZoGDx5MGRkZ\nFvcLDw+n5s2bk5eXF4WFhVncZ9OmTeTr60vOzs70119/WT3mE088Qa1ataKAgADq0KGDaDlCdLpx\n4wY9++yz1KxZMwoJCaFbt27ZrJOQ40ycOJG8vLyodevWFB0dLUrfyMhIevTRRwu7Gf773/+2KOeV\nV16hhg0bUsuWLS1uF6qPEFlCdbp8+TIZDAby9fUlPz8/WrhwoSS9zMhq7Hv27KGC+51yZ8yYQTNm\nzCi1T35+Pnl6etKlS5coNzeX/P39KS4urtR+Z86cofj4eDIYDGUaqYeHB924ccPqdiFyhOr07rvv\n0pw5c4iIKCwszOL5laWTkOPs3LmT+vTpQ0REhw8fpsDAQFFyIiMjqX///hb1K8q+ffsoOjraqoEK\n0UeoLKE6paamUkxMDBERZWZmkre3t6jfqSSypguEhITA2ZlFBgYGIjk5udQ+R48ehZeXFzw8PODq\n6ophw4Zh+/btpfbz8fGBt7e3oONSGZ6YEDlCddqxYwdGjx4NABg9ejR+/vlnm3QScpyixwgMDERG\nRgbS09NF6VvW72Kma9euqFOnjtXtQvQRKkuoTo899hgCAgIAADVq1ECLFi1w5coV0XqZUSw3ZtWq\nVQgNDS31eUpKCpo0aVL43t3dHSkpKaKP4+TkhGeffRbt27fHihUrRMkQqlN6ejrc3NwAAG5ublZ/\nXGs6CTmOpX1KDhpC5Dg5OeHgwYPw9/dHaGgo4uLirJ5/WQjRRyhidEpMTERMTAwCAwMl62VzAzEh\nsfdZs2bhoYcewogRI0rtV3SiKSQkBGfPnkV2dnaxzLiScfyyOHDgAF5++WWkpKRg0qRJCAsLQ/X7\n7eWEyhGi06xZs0p9x9qk2YEDB9CoUSNcu3YNISEh8PHxQdeuXQVPspUc/Up+T4ictm3bIikpCdWr\nV0d4eDgGDRqEc+fOCTq+rfoIxVadsrKyMHToUCxcuBA1atSQrJfNxh4REVHm9jVr1mDXrl34/fff\nLW5v3LhxYdJYREQEvvzySzg7O2PGjBm2qgIAaNSoUaFOM2fORI0aNTB16lSbZAjVyc3NDWlpaXjs\nsceQmppqNfrUqFEjAECDBg0wePBgHD16FF27di12HABISkqCu7u7VV0AIDk5GY0bNy5zH0tyiqZz\n9OnTB2+++SZu3ryJunXrlvt72KqPUGzRKS8vD0OGDMHIkSMxaNAgWfSS1Y3ZvXs3vvrqK2zfvh1V\nq1a1uE/79u2RkJCAxMRE5ObmYuPGjRgwYECZcq35ednZ2cjMzAQA3LlzB3v27EGrVq1sliNUpwED\nBmDt2rUAgLVr11r8I5Slk5DjDBgwAOvWrQOAwhQMs+tki77p6emF53v06FEQkc2GLlQfoQjViYjw\n6quvwtfXF1OmTJFPr3IfYW3Ay8uLmjZtWhhaGj9+PBERpaSkUGhoaOF+u3btIm9vb/L09KTZs2db\nlLV161Zyd3enqlWrkpubG/Xu3buUrAsXLpC/vz/5+/uTn5+fRVlC5AjV6caNG9SjR49SoUdbdLJ0\nnGXLltGyZcsK93nrrbfI09OTWrdubTWCVJ6cxYsXk5+fH/n7+1Pnzp3p0KFDFuUMGzaMGjVqRK6u\nruTu7k4rV64UpY8QWUJ1ioqKIicnJ/L39y+0pV27donWy4yia1B1dLSEw6xU0tGRim7sOpUG3dh1\nKg26setUGnRj16k06MauU2n4fy+vyqIST5m4AAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 5 }, { "cell_type": "markdown", "metadata": {}, "source": [ "
  • SciPy-\u6578\u503c\u8a08\u7b97\u5eab\uff1ahttp://myshare.dscloud.me/scipydoc/scipy_intro.html
  • \n", "
  • \u7e6a\u5716\u6307\u4ee4\u8acb\u53c3\u8003\uff1a http://onionys.blogspot.tw/2012/01/python-matplotlib.html
  • \n", "
  • \u66f4\u591a\u7684matplotlib\uff1a http://nbviewer.ipython.org/github/c3h3/NCCU-PyData-Courses-2013Spring/blob/master/Lecture4/visualization/BasicMatplotlib.ipynb
  • " ] }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 5 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 5 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 5 } ], "metadata": {} } ] }