{ "metadata": { "name": "" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "%pylab inline\n", "import numpy as np\n", "import matplotlib.pyplot as plt" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Populating the interactive namespace from numpy and matplotlib\n" ] } ], "prompt_number": 15 }, { "cell_type": "code", "collapsed": false, "input": [ "x = np.linspace(0.,1)\n", "v = 0*x\n", "u = np.exp(-100*(x-0.5)**2)\n", "plt.plot(x,u)" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 16, "text": [ "[]" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEACAYAAABI5zaHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X1QFGeCBvBndHAlavwOkRki8hHBIBBFwc3G4LoKcTfk\nzrhX5G5vEyUs8ZbLmdzV7VayuYCVy2ruqu7ccFenW35U1LDmEmvxEkRP4+hulOBXAL84YFGHUVlB\nifEjCkPfH++CIjAMQ3e/3T3Pr4qSWZrpx1548vp299s2RVEUEBGRpQyRHYCIiNTHcicisiCWOxGR\nBbHciYgsiOVORGRBLHciIgvqt9yXLVuGsLAwTJ8+vc9tXnnlFcTGxiIpKQnHjx9XNSAREQ1cv+W+\ndOlSlJWV9fn10tJS1NXVoba2FuvWrcPy5ctVDUhERAPXb7k/+eSTGDt2bJ9f37FjB1544QUAQGpq\nKlpbW9HU1KReQiIiGrBBz7l7PB5ERER0vXY6nWhsbBzs2xIR0SCockL1/hUMbDabGm9LREQBsg/2\nDRwOB9xud9frxsZGOByOHtvFxMSgvr5+sLsjIgoq0dHRqKurG/D3Dbrcs7KyUFRUhOzsbJSXl2PM\nmDEICwvrsV19fX2PEX6wKigoQEFBgewYhiDjWOzbB3z4Ye9fs9uB118HJk3SNRIA/lzci8firkBn\nQvot9+effx779+9Hc3MzIiIiUFhYiLa2NgBAXl4eFi1ahNLSUsTExGDEiBHYuHFjQEGI9KAoQH4+\n8NxzvRf4rl3AmjXAqlX6ZyNSU7/lXlxc3O+bFBUVqRKGSGtffAG0tQGFhUBvA6L0dOC73wXefluM\n4onMineoSpCeni47gmHofSzWrweWLeu92AEgPh6YMgUoLdU1FgD+XNyLx2LwbHo9rMNms3HOnaS6\nfh2IiABOnfI9p75+PbBjB1BSol82or4E2p0cuVPQ+PBD4Mkn+z9Z+hd/ARw4AFy8qE8uIi2w3Clo\nrF8P5OT0v92oUcDixcD772ufiUgrnJahoHD6tDhRev48EBLS//YHDwIvvgjU1PQ9P0+kB07LEPmw\nYQPw4x/7V+wAMGcOMHQo8LvfaZuLSCssd7K8tjZg82b/pmQ62WzASy+JqRwiM+KVvGR5n3wCPPqo\n+BiIv/5r8T1ffQWMHq1NNiKtcOROlufvidT7PfQQMH8+8JvfqJ+JSGssd7I0j0ecHF2yJLDvz8nh\n1AyZE8udLG3TJuCHPwRGjAjs+zMygAsXgOpqVWMRaY7lTpbV0SGukglkSqbT0KHikkiO3slsWO5k\nWfv3Aw88AMyaNbj3WbYM2LoVuH1bnVxEemC5k2V1nkgd7E1IUVHA9Olca4bMhXeokiW1tgKRkUBd\nHTBhwuDfb+tWsRzBrl2Dfy+igeAdqkT32L4d+N731Cl2QKw1U1EBNDWp835EWmO5kyX9/vei3NUS\nGiqWJDh0SL33JNISy50sqbwcSEtT9z3T0sT7EpkBy50sp7UVcLuBhAR135flTmbCcifLqagAZs5U\n/xmos2cDR48C7e3qvi+RFljuZDlaTMkAwJgx4jF9J06o/95EamO5k+VoVe4Ap2bIPFjuZCmKAnzx\nBZCaqs37s9zJLFjuZCn/93/iGaj9PQQ7UGlpvBySzIHlTpZSXi6uR9fKY48BFy8CLS3a7YNIDSx3\nshQt59sBsUrkrFniihwiI2O5k6VoXe4A593JHFjuZBk3bog59+RkbffDciczYLmTZRw5AiQmAt/6\nlrb7SU0VV+R0dGi7H6LBYLmTZegxJQOIB2ePHw/U1Gi/L6JAsdzJMvQqd4BTM2R8LHeyBEVhuRPd\ni+VOlnD+vPjzkUf02R/LnYyO5U6WcOiQKNzBPi/VX0lJQH098PXX+uyPaKBY7mQJWt+Zer9hw8Ql\nl4cP67dPooFguZMl6Dnf3olTM2Rk/ZZ7WVkZ4uLiEBsbi9WrV/f4enNzMzIzM5GcnIyEhARs2rRJ\ni5xEfbp9G6iuFg/o0BPLnYzMpiiK0tcXvV4vpk6dij179sDhcGDWrFkoLi5GfHx81zYFBQW4ffs2\nfvnLX6K5uRlTp05FU1MT7Pc9Bsdms8HHrogCVl4O/M3fAMeO6bvfxkZgxgygqUm/uX4KPoF2p8+R\ne0VFBWJiYhAZGYmQkBBkZ2ejpKSk2zaTJk3CtWvXAADXrl3D+PHjexQ7kZZkTMkAgNMp5t4bGvTf\nN1F/fJa7x+NBRERE12un0wmPx9Ntm9zcXJw8eRLh4eFISkrCmjVrtElK1AdZ5Q5waoaMy+cQ2+bH\nvzXfeecdJCcnw+Vyob6+HgsWLEBlZSVGjRrVY9uCgoKuz9PT05Genj7gwET3Ky8HVq6Us+/Ocv/L\nv5Szf7Iel8sFl8s16PfxWe4OhwNut7vrtdvthtPp7LbNwYMH8cYbbwAAoqOjMWXKFNTU1CAlJaXH\n+91b7kRquHhRXGseGytn/2lpwIcfytk3WdP9A9/CwsKA3sfntExKSgpqa2tx9uxZ3LlzB9u2bUNW\nVla3beLi4rBnzx4AQFNTE2pqahAVFRVQGKKB6nxeqqwTmjNmACdPArduydk/UV98jtztdjuKioqQ\nkZEBr9eLnJwcxMfHY+3atQCAvLw8vP7661i6dCmSkpLQ0dGBd999F+PGjdMlPJHM+XYAeOABID4e\nOH4c+Pa35eUgup/PSyFV3REvhSQNPPUU8ItfAAsWyMvwt38LREYCf//38jKQdWlyKSSRkbW3A0eP\nArNny83BK2bIiFjuZFonTohVIEePlpuD5U5GxHIn0/ryS3FCU7aoKODaNaClRXYSortY7mRaVVXi\nmamy2WzA9OlifRsio2C5k2kZpdwBkaOqSnYKortY7mRKigJUVrLcifrCcidTamoSBT9pkuwkAsud\njIblTqbUOSVjlKV2ExLEnaper+wkRALLnUzJSPPtAPDgg0BYmHiuKpERsNzJlIxW7gCnZshYWO5k\nSix3It9Y7mQ6bW1ATQ0wbZrsJN2x3MlIWO5kOjU1wOTJYkVGI2G5k5Gw3Ml0jDglAwDR0eISzT89\nUphIKpY7mY5Ry33oUOCxx8SCZkSysdzJdIxa7gCnZsg4WO5kOix3ov6x3MlUWlrEA7EnT5adpHcs\ndzIKljuZSnW1WF7XKMsO3K9z6V8+UZJkY7mTqRh5SgYAxo8HRo0Czp2TnYSCHcudTMXo5Q5waoaM\ngeVOplJVJaY+jIzlTkbAcifT8HrFsroJCbKT+MZyJyNguZNp1NcDDz0EjB4tO4lvLHcyApY7mYYZ\n5tsBYOpUcUL15k3ZSSiYsdzJNKqrzVHuISGi4E+dkp2EghnLnUzDLCN3gFMzJB/LnUyD5U7kP5Y7\nmcLXXwOXLgExMbKT+IflTrKx3MkUTpwQT14aOlR2Ev90ljuXISBZWO5kCmaakgGAsDBgyBDg4kXZ\nSShYsdzJFMxW7jYbp2ZILpY7mYLZyh1guZNcLHcyPEUxx5oy92O5k0wsdzK88+eBkSOBCRNkJxkY\nljvJ1G+5l5WVIS4uDrGxsVi9enWv27hcLjz++ONISEhAenq62hkpyJlxSgYQV/fU1gJ37shOQsHI\n7uuLXq8X+fn52LNnDxwOB2bNmoWsrCzEx8d3bdPa2oqf/vSn2LVrF5xOJ5qbmzUPTcHFrOU+fDgw\nZQpw5ow585O5+Ry5V1RUICYmBpGRkQgJCUF2djZKSkq6bfPBBx/gueeeg9PpBABMMNu/ncnwzFru\nAKdmSB6f5e7xeBAREdH12ul0wuPxdNumtrYWV65cwbx585CSkoLNmzdrk5SCFsudaOB8TsvY/HgK\ncVtbG44dO4a9e/fi5s2bmDNnDtLS0hAbG6taSApet24BZ8+KVRbNKDEReO892SkoGPksd4fDAbfb\n3fXa7XZ3Tb90ioiIwIQJExAaGorQ0FDMnTsXlZWVvZZ7QUFB1+fp6ek8+Ur9OnUKePRRYNgw2UkC\nk5goliom8pfL5YLL5Rr0+9gUpe/VL9rb2zF16lTs3bsX4eHhmD17NoqLi7udUD1z5gzy8/Oxa9cu\n3L59G6mpqdi2bRumTZvWfUc2G3zsiqhXGzcCn30GmHW2T1GAsWPFVTMTJ8pOQ2YUaHf6HLnb7XYU\nFRUhIyMDXq8XOTk5iI+Px9q1awEAeXl5iIuLQ2ZmJhITEzFkyBDk5ub2KHaiQJl5vh24uwxBdTXw\n3e/KTkPBxOfIXdUdceROAZg/H/jHfwQyMmQnCVx+vliqeMUK2UnIjALtTt6hSoalKEBlpblH7gCv\nmCE5WO5kWJcuiWmNhx+WnWRwWO4kA8udDKtzvt2PK3INLSFBXPXT3i47CQUTljsZlhWmZACx6Fl4\nOFBXJzsJBROWOxmW2a+UuRenZkhvLHcyLJY7UeBY7mRId+6IG38ee0x2EnWw3ElvLHcypDNnxHK5\nw4fLTqKOxERxDoFILyx3MiQrTckAQFQU0NICtLbKTkLBguVOhmS1ch8yRFwSyUXESC8sdzIkq5U7\nwHl30hfLnQzJiuWelMRyJ/2w3MlwLl8Gbt4E7nkImCVw5E56YrmT4VRXW2PZgftNnw6cOAF0dMhO\nQsGA5U6GY8UpGQAYMwYYNw5oaJCdhIIBy50Mx6rlDnBqhvTDcifDYbkTDR7LnQylvV0sj5uQIDuJ\nNljupBeWOxlKbS3gcIhlcq2I5U56YbmToVh5SgYAYmMBjwe4fl12ErI6ljsZitXL3W4H4uOBkydl\nJyGrY7mToVi93AFOzZA+WO5kKCx3InWw3MkwWluBK1fEOu5WxnInPbDcyTCqq8UlkEMs/lPZWe6K\nIjsJWZnFf43ITIJhSgYAJk4UT5hqbJSdhKyM5U6GESzlDnBqhrTHcifDYLkTqYflTobQ0SGWw50+\nXXYSfbDcSWssdzKEhgaxHO6YMbKT6INPZSKtsdzJEIJpSgYA4uLEf9C++UZ2ErIqljsZQlWVGM0G\ni2HDgJgYsQImkRZY7mQIlZXBNXIHOO9O2mK5kyFUVQXPydROiYniP2pEWmC5k3RXrgB//CPw6KOy\nk+hrxgzg6FHZKciqWO4k3ZEjouiGDpWdRF8pKcDx4+LpU0Rq67fcy8rKEBcXh9jYWKxevbrP7Q4f\nPgy73Y7t27erGpCs7/BhYNYs2Sn0N2YMEB4OnD4tOwlZkc9y93q9yM/PR1lZGU6dOoXi4mKc7uUn\n0ev14mc/+xkyMzOhcDUkGqDDh4HZs2WnkGP2bPH3J1Kbz3KvqKhATEwMIiMjERISguzsbJSUlPTY\n7r333sOSJUswceJEzYKSdVVUBOfIHRB/74oK2SnIinyWu8fjQURERNdrp9MJj8fTY5uSkhIsX74c\nAGCz2TSISVbl8Yg558mTZSeRgyN30orPcvenqFesWIFVq1bBZrNBURROy9CAdI7ag3VMkJws5tx5\npyqpze7riw6HA263u+u12+2G0+nsts3Ro0eRnZ0NAGhubsbOnTsREhKCrKysHu9XUFDQ9Xl6ejrS\n09MHEZ2sIFhPpnYaPlwsRfDll0Bamuw0ZAQulwsul2vQ72NTfAy129vbMXXqVOzduxfh4eGYPXs2\niouLER8f3+v2S5cuxTPPPIPFixf33NGfRvZE9/re94DXXgMWLZKdRJ6XXwamTQNeeUV2EjKiQLvT\n58jdbrejqKgIGRkZ8Hq9yMnJQXx8PNauXQsAyMvLCywtEcQyv0eOBPfIHRB/fxUGakTd+By5q7oj\njtzpPjU1QGamWB0xmFVXA0uWiONBdL9Au5N3qJI0wT7f3ik+HrhwAWhtlZ2ErITlTtJUVATvzUv3\nstuBxx8XU1REamG5kzQcud81axavdyd1sdxJijt3xDK/M2fKTmIMs2fzTlVSF8udpDhxApgyBRg5\nUnYSY+DIndTGcicpON/e3ZQp4i7VCxdkJyGrYLmTFJxv785m4+id1MVyJymCeZnfvnARMVITy510\nd+MGUF8ffM9M7Q+X/yU1sdxJd8eOAQkJwLBhspMYy6xZ4lp33shNamC5k+54MrV3YWHAqFFAXZ3s\nJGQFLHfSHU+m9o0nVUktLHfSHUfufePNTKQWljvpqrkZaGkBHn1UdhJj4sid1MJyJ10dPgykpABD\n+JPXq5kzgcpKoK1NdhIyO/6Kka443+7bgw8CjzwCnDwpOwmZHcuddMX59v5x3p3UwHIn3SgKR+7+\n4Lw7qYHlTro5f17MtTudspMYG0fupAaWO+mmc9Rus8lOYmyJiUBtrVimgShQLHfSDefb/fOtb4nl\nGY4fl52EzIzlTropL2e5+ys1FTh0SHYKMjOWO+ni66/FSPQ735GdxBzmzwf+939lpyAzY7mTLvbt\nE6PRBx6QncQc5s0TI/dbt2QnIbNiuZMudu8GFi6UncI8Ro8GkpKA3/1OdhIyK5Y76YLlPnAZGeK4\nEQWC5U6aa2gArl0Tl/iR/xYuBHbtkp2CzIrlTprbvRtYsICLhQ1USgrg8QAXLshOQmbEXzfSHKdk\nAjN0KK+aocCx3ElT7e3AZ5+JkTsNHOfdKVAsd9LU4cPA5MnAww/LTmJOCxaIkXtHh+wkZDYsd9IU\np2QGZ/JkYNw44MsvZSchs2G5k6Z27WK5D9bChZyaoYFjuZNmWluB6mouOTBYnHenQLDcSTOffQY8\n8QQwfLjsJOb21FPi3MX167KTkJmw3Ekzu3eLUScNzsiR4pr3/ftlJyEz8avcy8rKEBcXh9jYWKxe\nvbrH17du3YqkpCQkJibiiSeeQFVVlepByVwUhfPtauK8Ow1Uv+Xu9XqRn5+PsrIynDp1CsXFxTh9\n+nS3baKionDgwAFUVVXhzTffxE9+8hPNApM51NUBd+4A06bJTmINnHengeq33CsqKhATE4PIyEiE\nhIQgOzsbJSUl3baZM2cORo8eDQBITU1FY2OjNmnJNDovgeQj9dSRnAw0N4vn0BL5o99y93g8iIiI\n6HrtdDrh8Xj63H79+vVYtGiROunItHh9u7qGDBE3NHH0Tv6y97eBbQBDr3379mHDhg34/PPPe/16\nQUFB1+fp6elIT0/3+73JPNraAJcLWL9edhJrWbgQKC0FXnpJdhLSksvlgsvlGvT79FvuDocDbre7\n67Xb7YbT6eyxXVVVFXJzc1FWVoaxY8f2+l73ljtZ16FDQGwsMGGC7CTWsnAh8NprgNcrFhUja7p/\n4FtYWBjQ+/Q7LZOSkoLa2lqcPXsWd+7cwbZt25CVldVtm/Pnz2Px4sXYsmULYmJiAgpC1sEpGW2E\nhwMOB3DkiOwkZAb9jtztdjuKioqQkZEBr9eLnJwcxMfHY+3atQCAvLw8rFy5ElevXsXy5csBACEh\nIaioqNA2ORnW7t3Au+/KTmFNnZdEpqbKTkJGZ1MURdFlRzYbdNoVSdTSAkRFAZcvA8OGyU5jPbt2\nAW+/zWerBpNAu5N3qJKq9uwB5s5lsWtl7lyxQuRXX8lOQkbHcidV/fa3wNNPy05hXaGhYiG2Tz+V\nnYSMjtMypJrmZiAmRjwQu48LpkgF//3fwH/8h7jclKyP0zIk3aZNwLPPsti19uyzwJkz4oOoLyx3\nUoWiAOvWAS+/LDuJ9Q0bBixdKo43UV9Y7qSKffvEuu1pabKTBIfcXGDzZuCbb2QnIaNiuZMq/uu/\ngLw8LhSml6goYMYM4KOPZCcho+IJVRq0piYgLg44exb40+KgpIPt24F/+zde8251PKFK0mzcCCxe\nzGLX2zPPAPX1wMmTspOQEbHcaVA6OoBf/1pMyZC+QkKAZct4YpV6x3KnQdmzB3jwQWDWLNlJglNu\nLrB1K3DrluwkZDQsdxqUtWt5IlWmyZOB2bOBDz+UnYSMhidUKWAXL4pnpJ47J0bvJMeOHcCqVcDB\ng7KTkBZ4QpV0t2ED8MMfsthlW7QIcLuBqirZSchIWO4UEK+XJ1KNwm4HcnLEFBlRJ5Y7BWT3bvEY\nvZkzZSchQDxXtbgYuHFDdhIyCpY7BaTzRCoZg9MplgL+zW9kJyGj4AlVGrDGRiAxETh/Hhg5UnYa\n6vTpp0BhIcAnXFoLT6iSbv75n4Ef/YjFbjSZmcCVK0BpqewkZAQcudOAuFyi2E+cAMaMkZ2G7rd3\nr1gO+MQJXsVkFYF2J8ud/HbjBpCUJBareuYZ2WmoL7m5wNChYqVOMj+WO2nutdfECpBbt8pOQr58\n9RWQkAC8/z4wb57sNDRYLHfSVHk58Od/DlRXi0sgydg++QT4u78TNzaNGCE7DQ0GT6iSZr75Rqw+\n+KtfsdjN4gc/AL79beAXv5CdhGThyJ369cYbwOnTwMcfc4EwM2lpEdMzH38sip7MidMypInjx4GM\nDKCyEpg0SXYaGqiPPgLefFP8/zh8uOw0FAhOy5Dq2trEZXX/+q8sdrNasgR47DFg5UrZSUhvHLlT\nn95+G/j8c3FTDKdjzOvSJXEJa2kp1wIyI47cSVX79wNr1og1ZFjs5vbww+JfXz/+MfDHP8pOQ3ph\nuVMP778v1mkvLgYeeUR2GlLDj34kpmjS0vhA7WDBaRnq0tEB/NM/AR98IK6TnjZNdiJS25Yt4ma0\nLVuAhQtlpyF/8GoZGpRbt4AXXxQrPv72t8DEibITkVZ+/3vxL7O33gJefll2GuoP59wpYE1N4jZ1\nu10sPMVit7bvfEcU/L//O/Dqq+KpWmQ9LPcgd+IEkJoKPP20+Kc6r4UODtHRwKFDYnmCP/sz4Pp1\n2YlIbSz3IHX2LPAP/wCkp4v12d96i1fFBJuxY4GyMnE1TXIy8J//yZK3kn7LvaysDHFxcYiNjcXq\n1at73eaVV15BbGwskpKScPz4cdVDkjoUBdi3TywAlpIiyvzIEeCv/kp2MpIlJARYtw7YsAH47DNg\n8mRxwvUPf5CdjAbLZ7l7vV7k5+ejrKwMp06dQnFxMU6fPt1tm9LSUtTV1aG2thbr1q3D8uXLNQ1s\nBS6XS9f93bwJ/PrX4tF4+fniiT3nzgH/8i9AZKSuUXrQ+1gYmaxjYbMBc+eKpQqOHROFn5oKPPus\nOAcj4zoI/lwMns9yr6ioQExMDCIjIxESEoLs7GyUlJR022bHjh144YUXAACpqalobW1FU1OTdokt\nQMsfXK9XLPK1ZYs4WfbUU2LpgP/5H/GQjRMnxIOtjbIMLH+J7zLCsZg8GVi9WvzH/wc/AFasABwO\nICtLPJ/1k0+ACxe0z2GEY2F2dl9f9Hg8iIiI6HrtdDrxxRdf9LtNY2MjwsLCVI4a3BQFaG8HWlvF\nXYaXL3f/s6lJ3Jzy5ZfAQw+J28xnzhSLRs2YAYwbJ/tvQGbywAPiiU4vvSTOzxw9Kj5+9Svx57Bh\n4ucrJkb8vE2c2PPPkSPFE6FIDp/lbvPzDNv912D29X1mfTSbr3+W3v+1e18rSs8PAKivBw4cEKPs\njg7x570fbW1iDfVvvhHXn3d+PmQIMHr03V+ee3+REhLEtcszZogTZURqsNmAKVPEx5Il4n9TFOD8\neVHyZ8+KAUZ9vfjz3kHH9eui3ENDxVVY937Y7eJrQ4eKn+vOz4cOFfv8wx+AgwfF5511cu/nna/v\nz9rX3yEY+Sx3h8MBt9vd9drtdsPpdPrcprGxEQ6Ho8d7RUdH45NPgvQo96KhoXDA39PRIdbobmkR\nUy9WUVg48GNhVVY7Fu3twNdfi4+BCuR3xIqio6MD+j6f5Z6SkoLa2lqcPXsW4eHh2LZtG4qLi7tt\nk5WVhaKiImRnZ6O8vBxjxozpdUqmrq4uoIBERDRwPsvdbrejqKgIGRkZ8Hq9yMnJQXx8PNauXQsA\nyMvLw6JFi1BaWoqYmBiMGDECGzdu1CU4ERH1Tbe1ZYiISD+q36HKm57u6u9YbN26FUlJSUhMTMQT\nTzyBqqoqCSn14c/PBQAcPnwYdrsd27dv1zGdfvw5Di6XC48//jgSEhKQnp6ub0Ad9XcsmpubkZmZ\nieTkZCQkJGDTpk36h9TJsmXLEBYWhunTp/e5zYB7U1FRe3u7Eh0drTQ0NCh37txRkpKSlFOnTnXb\n5tNPP1WefvppRVEUpby8XElNTVUzgmH4cywOHjyotLa2KoqiKDt37gzqY9G53bx585Tvf//7ykcf\nfSQhqbb8OQ5Xr15Vpk2bprjdbkVRFOXy5csyomrOn2Px1ltvKT//+c8VRRHHYdy4cUpbW5uMuJo7\ncOCAcuzYMSUhIaHXrwfSm6qO3HnT013+HIs5c+Zg9OjRAMSxaGxslBFVc/4cCwB47733sGTJEky0\n6LKU/hyHDz74AM8991zXVWkTJkyQEVVz/hyLSZMm4dq1awCAa9euYfz48bDbfZ4mNK0nn3wSY31c\nwxxIb6pa7r3d0OTxePrdxoql5s+xuNf69euxaNEiPaLpzt+fi5KSkq7lK/y9x8JM/DkOtbW1uHLl\nCubNm4eUlBRs3rxZ75i68OdY5Obm4uTJkwgPD0dSUhLWrFmjd0zDCKQ3Vf3PoNo3PZnZQP5O+/bt\nw4YNG/D5559rmEgef47FihUrsGrVqq4HE9z/M2IF/hyHtrY2HDt2DHv37sXNmzcxZ84cpKWlITY2\nVoeE+vHnWLzzzjtITk6Gy+VCfX09FixYgMrKSowaNUqHhMYz0N5UtdzVvOnJ7Pw5FgBQVVWF3Nxc\nlJWV+fxnmZn5cyyOHj2K7OxsAOJE2s6dOxESEoKsrCxds2rJn+MQERGBCRMmIDQ0FKGhoZg7dy4q\nKystV+7+HIuDBw/ijTfeACBu5JkyZQpqamqQkpKia1YjCKg3VTsjoChKW1ubEhUVpTQ0NCi3b9/u\n94TqoUOHLHsS0Z9jce7cOSU6Olo5dOiQpJT68OdY3OvFF19UPv74Yx0T6sOf43D69Gll/vz5Snt7\nu3Ljxg0lISFBOXnypKTE2vHnWLz66qtKQUGBoiiKcunSJcXhcCgtLS0y4uqioaHBrxOq/vamqiN3\n3vR0lz/HYuXKlbh69WrXPHNISAgqKipkxtaEP8ciGPhzHOLi4pCZmYnExEQMGTIEubm5mGbBJ5X7\ncyxef/11LF26FElJSejo6MC7776LcRZdAe/555/H/v370dzcjIiICBQWFqKtrQ1A4L3Jm5iIiCyI\nj9kjIrJkavSNAAAALElEQVQgljsRkQWx3ImILIjlTkRkQSx3IiILYrkTEVkQy52IyIJY7kREFvT/\nErWhBsMXFcQAAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 16 }, { "cell_type": "code", "collapsed": false, "input": [ "tfinal = .2\n", "dx = x[1]-x[0]\n", "dt = 0.5*dx\n", "nsteps = int(round(tfinal/dt))\n", "for i in range(nsteps):\n", " # Diffusion\n", " v[1:-1] = v[1:-1] + dt/dx**2 * (u[2:] - 2*u[1:-1] + u[:-2])\n", " # Reaction\n", " u = u + dt*v" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 17 }, { "cell_type": "code", "collapsed": false, "input": [ "plt.plot(x,u)" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 18, "text": [ "[]" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAEACAYAAACnJV25AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XtUlXW6B/DvVvCCmikSKqDIRYGALYqZmSNmhZfJbMYm\nmlPWZOpYLsc6XWamVZl1OlqnnJKZFjOTY3WMWpoTlkGnVDIVBS+BCioa6HYblngP4/o7fzxtBYHN\nvrzv+3svz2ctljpu937mjf3lt5/3d7EJIQQYY4yZRifZBTDGGFMWBztjjJkMBztjjJkMBztjjJkM\nBztjjJkMBztjjJlMh8Gel5eHuLg4xMbGYunSpa3+Pj8/H71790ZKSgpSUlLw0ksvqVIoY4wxzwS4\n+8vGxkbMnz8fX375JcLCwjBq1ChMmzYN8fHxLR43fvx4rFu3TtVCGWOMecbtiL2wsBAxMTGIjIxE\nYGAgMjIykJOT0+pxvMaJMcb0w22wO51OREREXP5zeHg4nE5ni8fYbDZs27YNdrsdU6ZMQWlpqTqV\nMsYY84jbVozNZuvwCUaMGAGHw4GgoCDk5uZi+vTpOHTokGIFMsYY847bYA8LC4PD4bj8Z4fDgfDw\n8BaP6dWr1+XfT548GY888ghOnz6Nvn37tnhcTEwMjhw5okTNjDFmGdHR0Th8+LB3/0i4UV9fL6Ki\nokRFRYWora0VdrtdlJaWtnhMVVWVaGpqEkIIsWPHDjF48OA2n6uDl7KU559/XrPXmjVLiMWLhaiq\nEuKZZ4To10+IGTOE2LFDsxLc0vJa6J3sa1FXJ8S77wqRlERf770nRG2tEImJQuTlaVuL7GuhJ75k\np9sRe0BAADIzM5Geno7GxkbMmjUL8fHxyMrKAgDMnTsXa9aswVtvvYWAgAAEBQXhgw8+8PHnElPa\niRPA2rVAeTkQHAy89BLwxz8Cb78N/OY3QGQk8OyzwMSJsitlMtXXA5mZwLJlQGws8OqrwO23A65O\n7JNPAq+8AqSny62Tec5tsAPUXpk8eXKL/23u3LmXf//oo4/i0UcfVb4y5rc33wTuu49C3aVnT+AP\nfwAeeQRYvRr41a+AykqgTx9pZTLJVq0C3nuPBgGpqa3/PiMDeOYZYNcuYORI7etj3uOVpxKkpaWp\n/hrnzwP//Cfw+ONt/31gIPDb3wK33ALIXIKgxbUwClnXYvVq4Ikn2g51AOjSBXjsMRrJa4W/L/xj\n+7mHo/4L2Ww8311Dr70G7NwJZGe7f9z779PXp59qUxfTl7NngUGDAKcTaDYPopULF4AhQ4DCQiAq\nSrv6mG/ZycFuQnV19OZbtw4YMcL9Yy9cAMLDgaNHgWuv1aY+ph/vvAP8+9/Axx93/Ng//Qm4eBFY\nvlz9utgVvmQnt2JMKDsbiI/vONQBGqVNmCC3HcPkWb0auPtuzx67YAH140+dUrcm5j8OdpMRgnqh\nTz7p+b+5+25gzRr1amL6dO4csHkzcMcdnj1+wAC62f7Xv6pbF/MfB7vJ5OYCAQHAbbd5/m9++Uvg\nq6/ojc6sY906+rR2zTWe/5snngD+9jegpka9upj/ONhN5pVXgKeeujIH2RO9ewPjxwOffKJeXUx/\nVq8GZszw7t/ExQFjxgArV6pSElMI3zw1kcJCWnhUXk7TGb3x3nvUjmlj805mQufOARERgMNBP9i9\nsW0bcP/9wKFDQOfO6tTHruCbpxb36qs039jbUAeAadOA/Hya/87M75NP6FOat6EOADfdBPTvTwua\nmD5xsJvE4cMUzLNm+fbve/cGfvELbsdYhTezYdry1FPA0qV0s57pDwe7SSxbBvz+97RlgK9mzODZ\nMVZw/jywaRN9SvPVHXfQnPavvlKuLqYcDnYTEIJ64zNn+vc8d94JbNxIi5aYeX36KX0682dBWqdO\n1Gfn9Q/6xMFuAseOAQ0NQEyMf89z7bXAzTfz9gJm528bxmXcOGDrVv+fhymPg90EtmyhQPZmimN7\n7r6b3vjMnC5cADZs8K8N4zJqFLBvH/Djj/4/F1MWB7sJbN1Kwa6EO++kN/7Fi8o8H9OXTz+l7xUl\ntmnu3h2w22maLdMXDnYTcI3YldCnD01nW79emedj+rJmjTJtGJebb6bvP6YvHOwGd+YMHZQxfLhy\nz8ntGHO6eBH48kv6VKYUDnZ94mA3uIIC4IYbaH8YpUyfDnzxBfdOzWb9evo0dtU583656SZg+3a6\nec/0g4Pd4JRsw7j07Uv7gXA7xlyUmg3TXL9+QFgYsHevss/L/MPBbnBqBDvAW/maTU0NfQqbPl35\n5+Z2jP5wsBtYbS2wezcwerTyz33rrfxmNZPdu2lnRiXbMC4c7PrDwW5gu3bRm9XdWZW+GjQIqK8H\nTpxQ/rmZ9nbubP+wan+5gp33jdEPDnYDU6sNA9Bip9RUoKhInedn2lIz2IcMoVCvrFTn+Zn3ONgN\nTM1gBygIdu5U7/mZdoqK1At2m43bMXrDwW5QTU204nTsWPVeY9QoDnYzOHcOcDrpgHO1cLDrCwe7\nQR04QKtEBwxQ7zVcI3bunRrb7t20gE3JtQ5X42DXFw52g1K7DQMAAwfSaUxHj6r7OkxdavbXXZKT\n6Zi906fVfR3mGQ52g9Ii2AHus5uBmv11l4AA4MYb6TxUJh8Hu0FpFezcZze+nTvpv6PauB2jHxzs\nBuR00r7aw4ap/1o8Yje26mr6io1V/7U42PWDg92AXLNhlDhYoyOuYG9qUv+1mPJ27QJGjKCj7NQ2\nejTwzTfATz+p/1rMPQ52A9KqDQMAISF0ZN6RI9q8HlNWUZE2bRgA6NEDSEjgT3h6wMFuQFoGO8Dt\nGCPTYkZMc9yO0QcOdoM5fx44dIg+XmuFtxYwLq2DfexYDnY94GA3mO3bgZEjgS5dtHtNHrEbU1UV\nbdc7ZIh2rzl2LE155HsycnGwG4zWbRiAfpDs2QM0Nmr7usw/rtG6FjfZXfr3B4KDgdJS7V6TtdZh\nsOfl5SEuLg6xsbFYunRpu48rKipCQEAA1q5dq2iBrCUZwd6nD71hDx7U9nWZf7Ruw7hwn10+t8He\n2NiI+fPnIy8vD6WlpcjOzkZZWVmbj3v66acxadIkCN5YRDX19dTrHjNG+9fmPrvxcLBbl9tgLyws\nRExMDCIjIxEYGIiMjAzk5OS0etzy5csxY8YMhISEqFYoo3ZIVBRNP9Qa99mNRQgOditzG+xOpxMR\nERGX/xweHg6n09nqMTk5OZg3bx4AwKZlQ89iZLRhXHhrAWM5fpx+DQ/X/rWHDgV+/JE2BWNyuN3I\n05OQXrhwIZYsWQKbzQYhhNtWzKJFiy7/Pi0tDWlpaR4XyoDCQmDqVDmvnZIClJRQOygwUE4NzHMy\nbpy62Gy0IVhREdBsXMg8lJ+fj/z8fL+ew22wh4WFwdHsx67D4UD4VUOAXbt2ISMjAwBw6tQp5Obm\nIjAwENOmTWv1fM2DnXmvuBh45hk5r92rFzB4MLB/P+3tzfRNVhvGZfhw+n791a/k1WBUVw96X3jh\nBa+fw20rJjU1FeXl5aisrERdXR0+/PDDVoH97bffoqKiAhUVFZgxYwbeeuutNkOd+efSJdoXPS5O\nXg3cZzcO2cFut1OwMzncBntAQAAyMzORnp6OhIQE3HPPPYiPj0dWVhaysrK0qpEB2LePdnOU2Qbh\nPrsxyLxx6sLBLpdNaDQ/0dWDZ775xz9oV8eVK+XVUFAAzJ9POwYy/fr2W2D8eLk3L5uagN69qQYZ\ns7jMxJfs5JWnBlFcTKMgmex2oKyMt2XVO9mjdYC2CU5MpBvuTHsc7AZRXCz/pmVQEB3YsHev3DqY\ne3oIduDKDVSmPQ52A2hqopGP7BE7wH12I9DqKLyOcJ9dHg52A6isBK65BujbV3YlvLWA3jU10T2Q\nkSNlV0LB/s03squwJg52A9BDf92FpzzqW3k57a4YHCy7EiApiXZ5bGiQXYn1cLAbgJ6CPSkJOHyY\n9vlm+qOX/joA9OwJhIXRwTBMWxzsBqCHG6cuXbsC11/PH7H1Si/9dRe+gSoHB7sB6GnEDnCfXc+K\nivQzYgf4BqosHOw6d+4c8P33QHS07EquGDkS2L1bdhXsak1NFKIpKbIruYJvoMrBwa5zJSW00KNz\nZ9mVXJGURFscMH2prKSZU3pa6ckjdjk42HVOb20YAEhIoBWofAaqvuzdS4MAPYmIoJXK338vuxJr\n4WDXOT3dOHXp1YvOQD1yRHYlrLm9e+nTlJ7YbHwDVQYOdp375hv9jdgBGhny1gL6sm+f/kbsAPfZ\nZeBg17GGBlrgobdRGMB9dj3S44gd4D67DBzsOlZeDgwYQK0PveERu77U1tJ2vTIPYmkPB7v2ONh1\nTI83Tl2SkjjY9eTgQWDIEFpApjcJCbRambd71g4Hu459843+bpy6DB0KHDtGR/Yx+fTaXweAbt2A\nmBhqKzJtcLDrmJ5H7F260Ju1rEx2JQzQb3/dhdsx2uJg1zE9BztAI0S+gaoPeh6xAxzsWuNg16kf\nfqA2x6BBsitpH/fZ9YNH7Kw5DnadKi4GkpNpgYde8YhdH86fp4HAkCGyK2mfK9j5PHttcLDrlJ5v\nnLrwiF0f9u+nmSd62k/oaqGhNGPH4ZBdiTVwsOuU3vvrADB4MO0+eeaM7Eqsbd8+fbdhXLgdox0O\ndp0yQrB36kSHbnA7Ri49bv7VFg527XCw61BtLa06vf562ZV0jLcWkE/vN05dONi1w8GuQ6WlQFQU\nLezQO95aQC4hjDVi583AtMHBrkN63Kq3PTxil+vkSfq1f3+5dXhi2DDA6QQuXpRdiflxsOuQEfrr\nLq4RO09jk8O1MEnP02JdAgJo9g5/wlMfB7sOGSnYr7uOthc4cUJ2JdZklP66C/fZtcHBrjNCGCvY\nAe6zy6T3rQSuxsGuDQ52nTl+nD6yGqFn6sJ9dnmMOGLnG6jq42DXGSPdOHXhEbscTU00g8oI02Jd\n7HYaBDQ1ya7E3DjYdaakhPaIMRIesctRUQEEBwO9e8uuxHPXXks180Ho6uJg15mSEmP11wGa6VBW\nBjQ2yq7EWozWX3dJTuZPeGrjYNcZI47Ye/WiewI8CtOW0frrLsnJ9H3O1NNhsOfl5SEuLg6xsbFY\nunRpq7/PycmB3W5HSkoKRo4ciY0bN6pSqBX89BN9vNbjgcQd4T679ow8YudgV5fbYG9sbMT8+fOR\nl5eH0tJSZGdno+yqs9BuvfVWFBcXY8+ePVi5ciXmzJmjasFmVloKxMbSvHCj4T679njEztrjNtgL\nCwsRExODyMhIBAYGIiMjAzk5OS0e06NHj8u/v3jxIvr166dOpRZgxDaMC4/YtVVbC3z7rTE/3cXE\nAN99x1sLqMltsDudTkRERFz+c3h4OJxOZ6vHffzxx4iPj8fkyZPx5ptvKl+lRRg52HnErq2DB+nE\npK5dZVfivYAAID6ev1/U5DbYbR5uQDF9+nSUlZXhk08+wf33369IYVZk5GAfOhQ4epTOaWXqM8qO\nju3hdoy6Atz9ZVhYGBzNzrJyOBwIDw9v9/Hjxo1DQ0MDqqurERwc3OrvFy1adPn3aWlpSEtL875i\nk3JtJWDUYO/ShT5iHzgApKTIrsb8jHJqUns42NuXn5+P/Px8/55EuFFfXy+ioqJERUWFqK2tFXa7\nXZSWlrZ4zOHDh0VTU5MQQohdu3aJqKioNp+rg5eyvO++EyI4WIifL6UhZWQI8c47squwhqlThVi7\nVnYVvtuwQYhx42RXYQy+ZKfbEXtAQAAyMzORnp6OxsZGzJo1C/Hx8cjKygIAzJ07Fx999BHeffdd\nBAYGomfPnvjggw/8+0ljUa7RuhG2X20P99m1Y/QRe1ISjdiFMPb3vF7Zfv6JoP4L2WzQ6KUM6dVX\n6RCCv/xFdiW+W7cOeOstIDdXdiXmdv48MGAA/dq5s+xqfDdwIFBQQIeis/b5kp288lQnjHzj1IVH\n7NrYv5+2cTByqAPcZ1cTB7tOmCHYBw8Gzp4FzpyRXYm5GXVh0tU42NXDwa4DdXXAoUM0CjOyTp1o\nC9n9+2VXYm5G3UrganY7B7taONh14OBBGu0GBcmuxH+um2JMPUafw+7CI3b1cLDrgBnaMC5JSby1\ngJqEMM/3y7BhQGUlL2pTAwe7DpjljQrwKExtJ07QTdPQUNmV+K9LF1qxXFoquxLz4WDXATMFu2tm\nDM9sVcfevcZf79AcDwTUwcGuA2YK9uBgOnjj6FHZlZiTmb5XAA52tXCwS3bqFG1faqZFGvxmVQ8H\nO/MEB7tkrjnJZvloDfCbVU1mDPbiYm7dKY2DXTKzvVEBDna11NUB5eXGX+/QXP/+NKipqpJdiblw\nsEvGwc48deAAHa7RrZvsSpRjs/H3ixo42CUzY7APG0Y3T2tqZFdiLmb8XgE42NXAwS5RYyPN4TXD\nKsLmAgMp3Hl+srI42JmnONglOnyYFppcc43sSpTHb1blcbAzT3GwS2TWNyrAb1Y1mPX7JSGBNsGr\nq5NdiXlwsEtk1jcqwMGutB9+oHsWERGyK1Fe9+5AZCRthseUwcEukRWCnecnK8OM6x2a44GAsjjY\nJTJzsIeG0v7s330nuxJzMPP3CsDBrjQOdknOnQO+/x6IjpZdiTp4frKyONiZNzjYJdm3j04bMvq5\nle7wm1U5HOzMGxzskpj9jQrwm1UpDQ3mXO/Q3KBBtBneqVOyKzEHDnZJONiZpw4fBgYMoO2QzcrV\nuuPTt5TBwS6JFYI9Pp42reL5yf6xwvcKwAMBJXGwS9DUdGX6mpnx/GRluE5NMjsOduVwsEtw9Cht\nIxAcLLsS9fGb1X88Ymfe4mCXYM8ewG6XXYU2+M3qP6sEe1IS3SSur5ddifFxsEuwcycwapTsKrTB\nwe6fc+doO4GoKNmVqK9nT2rd7d8vuxLj42CXYOdOIDVVdhXa4GD3jxXWOzSXmkrvD+YfDnaNCUHf\nuCNHyq5EGzw/2T8lJea/yd4cB7syONg1VlEBBAXRvGQr4PnJ/rFKf92Fg10ZHOwas1IbxoXbMb6z\nWrDb7XQDtbZWdiXGxsGuMSvdOHXhYPeNVdY7NBcUBMTG8ic8f3Gwa4xH7MxTVlrv0By3Y/zHwa6h\npiZg1y7r3Dh1SUykj9eNjbIrMRartWFcRo3iYPcXB7uGDh8G+vYF+vWTXYm2evUC+ven///Mc1YN\ndh6x+8+jYM/Ly0NcXBxiY2OxdOnSVn+/atUq2O12JCcnY+zYsSjhz91tKiqyXhvGhdsx3rNqsCcl\n0eZxNTWyKzGuDoO9sbER8+fPR15eHkpLS5GdnY2ysrIWj4mKisLmzZtRUlKCZ599FnPmzFGtYCOz\nYn/dhYPde1YN9q5daWfQ4mLZlRhXh8FeWFiImJgYREZGIjAwEBkZGcjJyWnxmDFjxqB3794AgNGj\nR+P48ePqVGtwHOyyqzCOmhrg2DFg2DDZlcjB7Rj/dBjsTqcTERERl/8cHh4Op9PZ7uPffvttTJky\nRZnqTKSxkTb/GjFCdiVycLB7Z/9+CvXAQNmVyMHB7p+Ajh5gs9k8frJNmzZhxYoV2Lp1q19FmdGB\nA7TatE8f2ZXIERVFm1mdOwf8/OGOuWHVNoxLaiqwbJnsKoyrw2APCwuDw+G4/GeHw4Hw8PBWjysp\nKcHs2bORl5eHPu2k16JFiy7/Pi0tDWlpad5XbFBWXJjUXOfOtJnV3r3AzTfLrkb/iouttTDpatdf\nD1RW0j5DPXvKrkZb+fn5yM/P9+9JRAfq6+tFVFSUqKioELW1tcJut4vS0tIWjzl69KiIjo4WBQUF\n7T6PBy9lavPnC/Haa7KrkOuRR4R4/XXZVRjDjTcKsWmT7CrkGj1aiM2bZVchny/Z2WGPPSAgAJmZ\nmUhPT0dCQgLuuecexMfHIysrC1lZWQCAxYsX48yZM5g3bx5SUlJwww03+PfTxoSsfOPU5YYbgB07\nZFehf3V11Iqx2kK2q3Gf3Xe2n38iqP9CNhs0eindqa8Hrr0WqKoy90nzHTlwAJg8mXa4ZO0rKgIe\neoj3S1m5EvjiC2DVKtmVyOVLdvLKUw2UlgKDB1s71AFg6FDg7Fng++9lV6JvO3YAo0fLrkI+HrH7\njoNdA9yGIZ060Q1kbse4x8FO4uKAEydoJhXzDge7BjjYrxg9moO9IxzsJCAAGD4c2L1bdiXGw8Gu\nASvvEXM1Dnb3Tp8GvvuOpvsxet8UFcmuwng42FVWW0s99uHDZVeiD6NH0xu1qUl2JfpUVESzYaxy\neHVHuM/uGw52le3dSyfCBAXJrkQfQkJo6+KDB2VXok/chmmJg903HOwq4/56a9yOaR8He0uxsUB1\nNX0xz3Gwq4yDvTUO9rYJwcF+tU6daOO8XbtkV2IsHOwq42BvbfRooLBQdhX68+23QLduQFiY7Er0\nhdsx3uNgV9GlS8ChQ9bepa8tw4cDZWV0fdgVPFpvGwe79zjYVVRcTCfBdO0quxJ96d4dSEjg+clX\n42BvGwe79zjYVcRtmPZxn701Dva2RUUBFy4AJ0/KrsQ4ONhVxMHePg72lmpraWqs1Xd0bIvNRu8j\nvoHqOQ52FXGwt4+DvaXiYiAmxnqHSniK2zHe4WBXyYULtD1tYqLsSvQpNpY2d+KP16SwkParZ20b\nNYpnUnmDg10lW7bQG9WqhxF3pFMnPnijOe6vuzd2LLB1K9DQILsSY+BgV8mGDcAtt8iuQt+4HXMF\nB7t7oaFAeDjPpPIUB7tKNmwAJk6UXYW+cbCT06fpdK2EBNmV6NvEifS+Yh3jYFdBdTWtIhw1SnYl\n+nbDDbzTI0C9Y97RsWO33AJs3Ci7CmPgYFfBpk3AzTdzf70jISFAv350FqqVcRvGM+PHA9u3Az/9\nJLsS/eNgV8HGjdyG8RTfQOVg91Tv3nQAyfbtsivRPw52FfCNU89Zvc8uBLViONg9c8st3Gf3BAe7\nwo4fpx47b/zlGasH+5EjtHfOwIGyKzEGvoHqGQ52hW3cCEyYQPO0WcdSUug0pZoa2ZXIwW0Y79x0\nE1BSQgsAWfs4fhTG0xy9060b9U2tOj+Z2zDe6d6d7sts3iy7En3jYFeQEBzsvrByO4ZH7N7jdkzH\nONgVVF5OO9HFxMiuxFisGuy8o6NveD57xzjYFeQardtssisxFqsGe3ExbYbWo4fsSoxl1CigshL4\n4QfZlegXB7uCNm7kaY6+iI0FLl6kGUVWsm0bcOONsqswnoAAYNw4ID9fdiX6xcGukKYmWnHK/XXv\n2WzAbbcBeXmyK9HWZ58Bt98uuwpj4vns7nGwK6S4mJbH8wnzvpk6FVi/XnYV2rl4ESgooB9ozHt8\nA9U9DnaFcBvGP5Mm0TWsrZVdiTY2bKB7C716ya7EmBITgbNngWPHZFeiTxzsCuFpjv4JCaFta7/+\nWnYl2li/HpgyRXYVxtWpE8+OcYeDXQF1dXRiUlqa7EqMbcoUa7RjhKD++tSpsisxNm7HtI+DXQFF\nRTR3PThYdiXGZpU+e0kJrbgdOlR2JcbmGrELIbsS/eFgVwC3YZSRkkI3FcvLZVeiLlcbhtc7+Cc6\nmqY+HjwouxL96TDY8/LyEBcXh9jYWCxdurTV3x84cABjxoxBt27d8Nprr6lSpN7xjVNl2GzA5MnU\npjCz9eu5DaMEm4377O1xG+yNjY2YP38+8vLyUFpaiuzsbJSVlbV4THBwMJYvX44nnnhC1UL1qqYG\n2LmTFkww/5m9HVNdTdsIjB8vuxJz4D5729wGe2FhIWJiYhAZGYnAwEBkZGQgJyenxWNCQkKQmpqK\nQIueA7dlC7UQevaUXYk53HYbze++eFF2Jer4/HO6yd6tm+xKzOGWW2gFamOj7Er0xW2wO51ORERE\nXP5zeHg4nE6n6kUZCbdhlNWrF83vNusojNswyho4ELjuOlogyK4IcPeXNoXv7ixatOjy79PS0pBm\ngvmBGzYAFr21oBrXtMc775RdibIaG2nEvmSJ7ErMxdWOGTFCdiXKyM/PR76fG+G4DfawsDA4HI7L\nf3Y4HAgPD/f5xZoHuxlUVNAXb+SkrKlTgddfp2lsZpo5smMHjTCbfQhmCrjjDuDZZ4Enn5RdiTKu\nHvS+8MILXj+H21ZMamoqysvLUVlZibq6Onz44YeYNm1am48VFpxM+t57QEYG0KWL7ErMZehQ6kGX\nlMiuRFnchlHHrbfSzqBXzeuwNLfBHhAQgMzMTKSnpyMhIQH33HMP4uPjkZWVhaysLABAVVUVIiIi\nsGzZMrz00ksYNGgQLpr1zlczQgDvvgs88IDsSszHZjPnKlRebaqOzp2B++6j9yMjNqHRUNtms5lq\nVL9lCzBnDrB/v7naBXrx+efAiy/SdTYDpxNITgZOnqRFNUxZ+/bRRnJHj1LQm4kv2ckrT330zjvA\nzJkc6moZP57me1dXy65EGa691znU1ZGYCISG8mIlFw52H1y6BHz0EX38Y+ro1o3me3/+uexKlMH9\ndfU98AC3Y1w42H2QkwOkpgJ+TBBiHpgyxRzbC9TW0ula6emyKzG3e+8FPvkEuHBBdiXycbD74J13\n+KapFqZMoePyjL6qcPNm2ms+JER2JeYWEkItvDVrZFciHwe7l06cALZvB+66S3Yl5hcRQUcN7tgh\nuxL/cBtGOw88QAMvq+Ng99KqVRTqQUGyK7EGM7RjPvuMT0vSytSpNFOtokJ2JXJxsHtBCG7DaG3q\nVODTT2VX4btDh2hDs5QU2ZVYQ9euwD33AP/7v7IrkYuD3Qt79gA//shb9GppzBjg3DnjtmOysoDf\n/panxWpp5kyaHWOiZTNe42D3gmvueie+aprp3BlYuNCYG62dOwesXAksWCC7EmsZNYrWC2zbJrsS\neXjlqYfq6mh6Y0EBHcnFtHPhAjBkCB1oEhkpuxrPvf461fz++7IrsZ4lS4BvvwX+/nfZlfiPV56q\nKC+PNqfiUNder17AQw8Bb7whuxLPNTRQvY8/LrsSa7rvPpr2eOmS7Erk4GD3EN80lWvBAvpvcPas\n7Eo8s2YNfbpITZVdiTWFh9O1X7dOdiVycLB7oLqaNvL/zW9kV2Jd4eE0ZfAf/5BdSceEoHsC//mf\nsiuxtpny0q5xAAALOUlEQVQzrTunnYPdAx98AEyeDPTuLbsSa3v8ceDNN4H6etmVuPf113Tj9Je/\nlF2Jtd11F90T++472ZVoj4O9Aw0NwN/+Bvzud7IrYSNGADExwOrVsitx7/XXgcce49lTsvXoAcyY\nQe9fq+FZMR3IzATWrqVWDM9Flu/TT4Hnn6fZJnr871FeDowdC1RW8upkPTh2jBaH7d4NDB4suxrf\n8KwYhVVXA4sX0+wGPYaIFU2ZQovEvvpKdiVtW7YMmDuXQ10vBg2iG+9PPCG7Em3xiN2NRx+lQM/M\nlF0Jay4ri0bun3wiu5KWqqupVVRWBvTvL7sa5lJTQ7trrlxJe/wbjS/ZycHejpIS4Lbb6E3at6/s\nalhzly7RVMKvvgLi4mRXc8VLL9GimBUrZFfCrrZ6Nf332bXLeKdYcStGIUIAf/gD9XI51PWne3fg\n978H/vIX2ZVcUVsL/PWvvCBJr2bMAPr0McZ0WSXwiL0Na9ZQb333buP9dLeKkydptH7okD4OsPjX\nv4APP6QVykyfiovp3FmjfQrnVowCLl0C4uPpjTphguxqmDsPP0yHcTz/vNw6hACSk2lR0u23y62F\nuTdvHhAYSOshjIKDXQEvvkg/2fl4Lf07dIimFm7YQMEqy3//N93M3bKFZ0/p3alTdCN140YgMVF2\nNZ7hYPeTwwEMH043WIy0i6CVrVoFPPcczWvv00f71/+//wMefBAoLOTDzY1i+XI6kP6LL4zxg5hv\nnvrpqadoiiOHunH8x38Ad9xBv2p96HVFBXD//UB2Noe6kcybB1RVAR9/LLsS9fCI/Wdff00n3Rw4\nQEuRmXHU1wO33kon1C9erM1r1tRQG+jBB2kGFTOWL78E5syh81G7d5ddjXvcivHRgQPAxIk0XW36\ndNnVMF+cPEnbtGZmAnfeqe5rCUE7BzY10dmaRvg4z1q7/37aBnrNGjorVa+4FeMDV6i//DKHupGF\nhtIb9OGHgYMH1X2tzExg716aE82hblwrVgBdutAc99pa2dUoy9Ij9uahzodomMPf/04Ll3bsoJOX\nlPb11xQEBQVAVJTyz8+0VV8PZGTQ0Zd6HblzK8YLBw5QX/a//otD3WxmzwbOnKFl5EqOqJ1OOih5\nxQpg0iTlnpfJpfdw51aMh1yhziN1c1q+nLZrXbAAOH9emecsLaXe/aOPcqibTWAgHaZjpraM5YK9\neajPnCm7GqaGbt1owdCFC8CwYcBbb/l+6lJVFW3Dm5ZGByT/6U+Klsp0wmzhbqlgLyzkULeK666j\nbVpzc4GPPgKSkuhgY08/0f74I02dTEykXv3Bg8DChXwqkpk1D/df/5pWqRqVJb5NDx6kg6jvuosO\nQuBQt47hw2mF4bJlwJ//TCPvoqL2H9/YCLz9NjB0KG0WVVQE/M//yFnVyrTnCvfoaNpk7sUXgYsX\nZVflPVPfPHU6gRdeAP79bzoxfsECPtnGyhoaaBT/3HP057ZurP70E3D99RTmN9ygaXlMZ44cAZ59\nFti0CXjmGVrQ1KWL9nWoMismLy8PCxcuRGNjIx5++GE8/fTTrR6zYMEC5ObmIigoCCtXrkRKSooi\nxfnq9Glg6VLgn/+kGRJPPWWsbTqZun76iU47akunTnT6Ec9PZy7ffEOf9g4coBH8vfdq25LzKTuF\nGw0NDSI6OlpUVFSIuro6YbfbRWlpaYvHrF+/XkyePFkIIcT27dvF6NGj23yuDl7Kbw6HEO+8I8TM\nmUL06yfEnDlCHD+u6kv6bNOmTbJL0A2+FlfwtbhCj9ciP1+IG28UIi5OiCeeECI3V4iLF9V/XV+y\n0+3PncLCQsTExCAyMhKBgYHIyMhATk5Oi8esW7cOD/w8Z3D06NE4e/YsTp486d1PFx9UV9NNsUce\noZkPw4fTGZhjxtBOf1lZQFiY6mX4JD8/X3YJusHX4gq+Flfo8VqMHw9s20b3YHr1ou2aQ0OBceOA\nRYto8VpdnewqidvzgZxOJyIiIi7/OTw8HDt27OjwMcePH0doaKjXxdTX0xQ119e5c9QnP3as9Vdd\nHV3QiROp95WczDMWGGPqstmAm26ir+eeo83gtmyhMwEef5zOSg4NBQYNav0VGko/EFxfPXuql1lu\ng93mYaNRXNX/ae/fTZhAgVxbe+XX2lqaWnbhAgV78//j11xDo+7Bg2n5dlralYsUEsJ9UMaYXEFB\ndGqW6+Ss+vrWg9GSEuomfP99y4FrTQ39+2uuoR0mu3alm7NX/+oTd32agoICkZ6efvnPL7/8sliy\nZEmLx8ydO1dkZ2df/vOwYcNEVVVVq+eKjo4WAPiLv/iLv/jLi6/o6Give+xuR+ypqakoLy9HZWUl\nBg4ciA8//BDZ2dktHjNt2jRkZmYiIyMD27dvx7XXXttmG+bw4cPuXooxxphC3AZ7QEAAMjMzkZ6e\njsbGRsyaNQvx8fHIysoCAMydOxdTpkzBZ599hpiYGPTo0QP/+te/NCmcMcZY2zRboMQYY0wbit+T\nzcvLQ1xcHGJjY7F06dI2H7NgwQLExsbCbrdjz549SpegGx1di1WrVsFutyM5ORljx45FSUmJhCrV\n58n3BAAUFRUhICAAa9eu1bA6bXlyLfLz85GSkoLExESkpaVpW6CGOroWp06dwqRJkzB8+HAkJiZi\n5cqV2hepkYceegihoaFISkpq9zFe5abXXXk3lFzQZHSeXItt27aJs2fPCiGEyM3NNeW18OQ6uB43\nYcIEMXXqVLFmzRoJlarPk2tx5swZkZCQIBwOhxBCiB9++EFGqarz5Fo8//zz4o9//KMQgq5D3759\nRX19vYxyVbd582axe/dukZiY2Obfe5ubio7Y9bygSWueXIsxY8agd+/eAOhaHD9+XEapqvLkOgDA\n8uXLMWPGDISEhEioUhueXIv3338fv/71rxEeHg4A6Nevn4xSVefJtRgwYADO/7yh/vnz5xEcHIyA\nALe3BQ1r3Lhx6ONmpzlvc1PRYG9rsZLT6ezwMWYMNE+uRXNvv/02pkyZokVpmvL0eyInJwfz5s0D\n4Pn6CaPx5FqUl5fj9OnTmDBhAlJTU/Hee+9pXaYmPLkWs2fPxv79+zFw4EDY7Xa88cYbWpepG97m\npqI//pRe0GRk3vx/2rRpE1asWIGtW7eqWJEcnlyHhQsXYsmSJZc3O7r6+8MsPLkW9fX12L17NzZs\n2ICamhqMGTMGN954I2JjYzWoUDueXIuXX34Zw4cPR35+Po4cOYLbbrsNxcXF6KXGYbYG4E1uKhrs\nYWFhcDgcl//scDguf6Rs7zHHjx9HmF43dfGDJ9cCAEpKSjB79mzk5eW5/ShmVJ5ch127diEjIwMA\n3TDLzc1FYGAgpk2bpmmtavPkWkRERKBfv37o3r07unfvjl/84hcoLi42XbB7ci22bduGZ555BgAQ\nHR2NIUOG4ODBg0hNTdW0Vj3wOjeVvAFQX18voqKiREVFhaitre3w5mlBQYEpbxgK4dm1OHr0qIiO\njhYFBQWSqlSfJ9ehuQcffFB89NFHGlaoHU+uRVlZmZg4caJoaGgQP/74o0hMTBT79++XVLF6PLkW\njz32mFi0aJEQQoiqqioRFhYmqqurZZSriYqKCo9unnqSm4qO2HlB0xWeXIvFixfjzJkzl3vLgYGB\nKCwslFm24jy5DlbhybWIi4vDpEmTkJycjE6dOmH27NlISEiQXLnyPLkWf/7zn/G73/0OdrsdTU1N\neOWVV9DXpAcr3Hvvvfjqq69w6tQpRERE4IUXXkD9zwf1+pKbvECJMcZMhje6ZYwxk+FgZ4wxk+Fg\nZ4wxk+FgZ4wxk+FgZ4wxk+FgZ4wxk+FgZ4wxk+FgZ4wxk/l/A3IF4kFVrxEAAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 18 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }