Skip to content

Instantly share code, notes, and snippets.

@ceptreee
Created January 19, 2019 21:53
Show Gist options
  • Save ceptreee/068dc2356bce7c6a3eddca99fc221128 to your computer and use it in GitHub Desktop.
Save ceptreee/068dc2356bce7c6a3eddca99fc221128 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np \n",
"import matplotlib.pyplot as plt \n",
"import sympy as sympy\n",
"\n",
"# おまじない\n",
"plt.rcParams['text.usetex'] = True\n",
"\n",
"x = sympy.Symbol('x')\n",
"xx = np.linspace(0,1,1000)\n",
"\n",
"f = sympy.sin(sympy.pi*x)\n",
"ff = sympy.lambdify(x,f) "
]
},
{
"cell_type": "code",
"execution_count": 23,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xd8U+XiBvDnTbpLB13QSUlLadmUIiCCikUFcWHFPS5X6nVeJw4Ur9eN2+t14M+9FUG8AgIVNzJaNqVAW0YphRZKB91J3t8fCcMCTUeSNzl5vp+PH2lP2jw9tE/DOe8QUkoQEZG26FQHICIi+2O5ExFpEMudiEiDWO5ERBrEcici0iCWOxGRBrHciYg0iOVORKRBLHciIg3ycuaTRUREyMTERGc+JRGR28vLyzsgpYzsyMc4tdwTExORm5vrzKckInJ7QohdHf0YXpYhItIgljsRkQax3ImINIjlTkSkQSx3IiINYrkTEWmQzXIXQmQKIZa2cTzL+pjp9o1GRESdZXOcu5QyRwjxwMmOCSGyjnuMQQiRKaXMsXdIos6SUuJgXTP2HGpA6aEG7DlUj4raJtQ2GlHb1ILaRiNMZgmdEBAC0OsEgvy8EervjdAAb4QH+iAhPAAJYQGI6x4AP2+96i+JqF26OolpOIAvrX8uBpAO4C/lLoTIBpANAAkJCV18OqK27a9pxModldhcWo38shrk763Bwbpmu33+hLAADIwNwYDYEAyMDcHQhFAE+jp1LiBRu3T1uzK01dvhrR8gpZwNYDYAZGRkcDdusqvDTUb8srUCvxdWYEVxJXYcqDvhMUG+XogPC0Bsd3/EhvqjR7Afgv29EOznjW5+XvDW6SAhYZaA0WRGbaMRVfXNqGpoQXltE0oq67G7sh6lhxqw2/rnBRvLAABeOoEh8aE4PTkCY/pEID2hO/Q64ezTQHSCrpZ7FYAwewQhaq9Ddc34YfM+LNm8D38UHkSzyXz0WKCPHhmJYRiaEIp+0cHoFxOM2FB/CNH1wm0xmVFYfhibSquxqbQa60qqsLG0Grm7DiF31yG89uN2RHTzwbn9e2LCgJ4YaQiHt55jFkiNrpb7ahx79W4AcMobr0RdYTSZ8dv2A/g6rwQ5+eVHC10IIKNXd5yT1gOjksIxICYYXg4qVG+9DmnRwUiLDsblGfEAgJrGFqwsrsQfhQewrKAcuyvr8dnK3fhs5W5EdPPBZelxuDwjHslR3RySiehUhJRtXymx3jR9B8A0KeUc6/uWSinHW/88HcAaAAbrJZhTysjIkFw4jDri4OEmfLZyNz5ZuQv7a5oAWAr9jOQITBwYjcy0HogM8lWc0kJKifyyGvywaR8WbCxDccWxS0QZvbrjulG9MHFgNF/NU4cJIfKklBkd+hhb5W5PLHdqr+37a/Hu7zswb20pmoyWV+mGiEBcNiwOk9NjER3irzhh26SUWFtSha9Wl+B/6/eirtkEAIgN9cffRifiiuHxCPLzVpyS3AXLndze9v21eOXH7Vi4sQxHvjXPSY3C1DN64/SkcLtcO3e2uiYjvl1Xind/33H01XyQrxemntEbfx/TG8EsebKB5U5uq7D8MF79cTu+37AXUgI+eh2mDI/D1NG9YYjUxvVqs1liWUE5Zv9WjFU7KgEAIf7eyB5rwI2nJ3JIJZ0Sy53czqG6ZrySsw2frNwNk1nCWy9w5fAE3Hp2kstfeumKlcUH8eLSbUdLPqKbL6af1xdZw+Kg41BKaoXlTm6jxWTGJyt24ZWc7ahuaIFOAFcMj8ft4/ogNlS7pX48KSWWFx3EC0u2Yu3uKgDAwNgQzLywH4YncoQxHcNyJ7ewemclHpq7EYXlhwFYRr48MikNqT2DFSdTQ0qJ79bvxbOLClBW3QgAuGhwDB6d1M9lRgKRWix3cmk1jS14blEBPl25GwDQOyIQMyam4Zy0KLe8UWpv9c1GvPVLMd7+pQhNRjNC/L3x6KR+uCw9lufHw7HcyWUt3rwPM+dvwv6aJnjrBW45Mwm3np3MhbhOoqSyHg/P24jfth8AAIzpE4GnLx2I+LAAxclIFZY7uZzaxhY8/r98zMnbAwAYmhCKZycPQt+eQYqTuTYpJeauKcUTC/JRVd+CAB89HruwH6ZkxPNVvAdiuZNLyd1Zibu/WoeSygb4eunw4IRUXD8qkQtrdcCBw0147LvNWLDBslDZ+f174pnJA9E90EdxMnImlju5BKPJjNd+3I7XfyqEWQL9Y4Lx6pVDkBzFV+udIaXEvLWlmDl/Mw43GREV5IsXpwzGmD6RqqORk3Sm3LnIBdlVRW0Trnt3FV5bVggJ4JazkjDv1tEs9i4QQmByehwW/XMMMnp1R3ltE65/bxVeXroNJjNX0aaTY7mT3azaUYkLXvsNfxYfREQ3X3x200g8cH4qfLz4bWYP8WEB+CJ7JO7OTAEAvPrjdtz4/ipU2nEzEtIO/tRRl0kp8c6vxbjqnRUor23CaYlhWHjnGRiVdMLeLdRFXnod/pnZBx9NPQ1hgT74bfsBTHrtN6zdfUh1NHIxLHfqksYWE/75xTo8tXALTGaJm8ca8Nm0EYgK9lMdTdPG9InE93ecgaEJodhb3Ygr3l6BuWv2qI5FLoTlTp1WXtuIK2evwHfr9yLQR4+3rh2GhyamOWyzDPqrmFB/fJk9CteN7IVmkxn3fLUes34ogJnX4Qksd+qk/L01uOT1P7CupAqxof6Yc8vpOH9AT9WxPI6Plw5PXDIAT1zcH3qdwBs/F+GWT/NQ32xUHY0UY7lTh+Xk70fWW8uxt7oR6Qmh+Pa20UiL9sx1YVzFdaMS8cHfhiPIzwuLN+9H1pt/Yp91nRryTCx36pDPV+1G9se5qG824ZIhMfhs2kgubuUixvSJxLxbR6NXeADyy2pw2ZvLjy7ORp6H5U7tIqXEaz9ux0NzN8IsgTvP6YOXrxjCtWFcTHJUN3x762ikJ4SitKoBWW8txxqOpPFILHeyyWSWmDl/M15aug1CAE9eMgD3jE/hGicuqnugDz69aSTGpUahqr4FV7+zAj8VlKuORU7Gcqc2NRvNuPPztfh4xS746HV44+p0XDuyl+pYZIO/jx5vXzcMWcPi0Nhixk0f5eKbPA6V9CTctJFOqbHFhFs/XYNlBeUI8vXCOzdkYKSBE5Pchbdeh+ezBiEqyBdv/FyE++asR7PJjKtOS1AdjZyAr9zppBqaTZj2US6WFZSje4A3Ps8eyWJ3Q0IITD8/FQ9NSIWUwENzN+LD5TtVxyInYLnTCeqajPjbB6vw2/YDiOjmg8+zR2JAbIjqWNQFN5+ZhJmT+gEAHvtuM975tVhxInI0ljv9RW1jC254bxVWFFciKsgXX2SP9Ni9TbVm6hm98eQlAwAATy3cgv/+VKg4ETkSy52OOtxkxA3vrULurkOIDvHDlzeP4lK9GnPtyF6YddkgCAE8v3grXl+2XXUkchCWOwGwXGP/+wersWa3ZTmBr24ehd4RgapjkQNMGR6PFy8fDCGAF5Zsw//9xks0WsRyJzQZTcj+OBcrd1guxXw2bQQ3Y9a4yelxeG7yIADAkwu24JMVuxQnIntjuXu4FpMZt326Fr9tP4DwQB98Nm0EeoXzFbsnmDI8Hv++uD8A4JFvNx3dxJy0geXuwUxmibu+XIecLfsR4u+Nj/8+gtfYPcz1oxLx8MRUAMD0Oevx/Ya9ihORvbDcPZSUEjPmbcSCDWXo5uuFj6aehn4xHBXjibLHJuHuzBSYJXDXF+vw81YuVaAFNstdCJElhMgUQky3cTzb/vHIUV5eug1frC6Br5cO7/9tOAbHh6qORArdeU4ybh5rgNEsccsna7CupEp1JOqiNstdCJEFAFLKHABVQojMVsczARRbjxcLIdIdlpTs5uMVu/DaskLoBPD61ekYnhimOhIpJoTAgxNScVl6HBpaTJj6wWoUV3C5YHdm65X7cABHxkkVA2hd3rkAvraWukFKucbO+cjOfthUhpnzNwEAnr50IMb366E4EbkKIQSevWwgzu4bicq6Zlz37irsr+GGH+7KVrm3/rf6XxYXkVJWAXgbwNcAkk72CYQQ2UKIXCFEbkVFRaeDUtetKD6IO79YBymBe8en4EouIEWteOt1+O816RgSb1kP/ob3VqGmsUV1LOoEW+VeBeCU/2a3XrbJkVImHff2X0gpZ0spM6SUGZGRkV0KS523dV8tpn2Ui2ajGdeN7IXbxyWrjkQuKsDHC+/dOByGyEAU7KtFtvX7htyLrXJfjWOv3g0AlrY6fvylmGfQxi8CUqe8thF/e38VahuNOL9/T/zrov7caIPaFBbog4+mnobIIF+sKK7EjHkbIaVUHYs6oM1yl1LOAWA4ciPVeuMUQogjJT/betklE8AUKeVsh6alDmtsMWHaR3nYW92IoQmheOXKIdDrWOxkW1z3ALx7Qwb8vHX4Om8P3vylSHUk6gDhzN/GGRkZMjc312nP5+nMZok7Pl+LBRvLEBvqj29vG83NrKnDFm/eh398kgcpgTeuScfEgdGqI3kcIUSelDKjIx/DSUwa9tLSbViw0TJJ6b0bh7PYqVPO698TD02wzGK9+8t1HAPvJljuGvVN3h68/tORsexD0bcnlxWgzps2xoCrTotHk9GMmz7MxZ5D9aojkQ0sdw1avbMSD87dAAD410X9cVbfKMWJyN0JIfDviwdgdHI4Dhxuwk0f5qKuyag6FrWB5a4xe6sacMsneWgxSdx4eiKuH5WoOhJphLdehzeuGXZ0iOR9X6/nCBoXxnLXkMYWE/7xSR4OHG7G6ORwPHJBmupIpDEh/t545/oMBPl6YdGmfXjjZ46gcVUsd42wrPK4CRv2VCOuuz9evyodXnr+9ZL9JUV2w6tXDbHu5LQVywr2q45EJ8Gffo34cPlOfLNmD/y8dZh9XQa6B/qojkQaNi61B+4dnwIpgX9+vg5FXGTM5bDcNeDPooN4YsEWAMCsrMFcl52c4razkzFhQE/UNhkx7aNcrkHjYljubq60qgG3fbYGJrPEzWcacNHgGNWRyEMIIfDC5YPRt0cQiivqcM+X62A28warq2C5u7HGFhP+8XEeKuuaMaZPBKafl6o6EnmYQF8vzL5+GEL8vZGzpRxv/cobrK6C5e7GnlyQj42l1YgP88d/rhrKNWNIiV7hgXj5isEAgBcWb8WfRQcVJyKA5e625q8rxScrdsNHr8MbVw9DaABvoJI641J74Lazk2CWwB2fr0U5N/lQjuXuhgrLD+OhuRsBADMv7IeBcSGKExEBd2emYJTBMoP19s/XwmjiGvAqsdzdTH2zEbd+mof6ZhMuGhyDa0ZwNyVyDV56HV69agiignyxakclXliyTXUkj8ZydyNSSjzy7SZs238YhshAPD15IDfdIJcSFeR39P7PW78UYWk+JzipwnJ3I1/llmDumlL4eevw5jXD0M3XS3UkohOMMIRj+nl9AQD3fLUOJZVcQVIFlrubyN9bg5nzNwMAnrxkIJfwJZeWPdaAzLQeqG004s4v1qKF19+djuXuBuqbjbj98zVoMpoxJSMOWcPiVEciapNlgtMgRIf4Ye3uKrySw+vvzsZydwP//l8+iivq0CeqGx6/aIDqOETtEhrgg1euGAKdAN74uQjLCw+ojuRRWO4ubuHGMnyxugQ+Xjr85+qh8PfRq45E1G4jDOG4Y1wfSAnc9eU6HDzcpDqSx2C5u7DSqgY8+I1lR6UZE9OQ2pMLgpH7uWNcMk5LDEN5bRPun7OBG3w4CcvdRRlNZtz1xVrUNBqRmRaF60f1Uh2JqFO89Dq8cuUQhPh7Y1lBOd7/Y6fqSB6B5e6iXv+pEKt3HkJUkC9mZQ3meHZyazGh/njuskEAgGcXFWBTabXiRNrHcndBuTsr8dqP2yEE8PIVQxDGjTdIA84f0BPXjkxAs8mMOz9fi/pmbrDtSCx3F1Pd0IJ/frEOZmkZKzw6OUJ1JCK7eeSCfkjp0Q3FB+rwzMIC1XE0jeXuYv713WaUVjVgUFwI7h3fV3UcIrvy89bj5SuGwFsv8PGKXfh5a7nqSJrFcnchCzeWYd5ay/ICr1wxBD5e/Osh7ekfE4J7rC9cps/ZgEN1zYoTaRPbw0WU1zRixjzLMr4zJqbBENlNcSIix8kea8DwxO4or23CjG83cnikA7DcXYCUEg98swGH6lswpk8Erh3JYY+kbXqdwEtThiDQR4+FG/fh23WlqiNpDsvdBXyxugQ/ba1AsJ8XnuewR/IQ8WEBeOyi/gCAmd9a7jWR/bDcFdt1sA5PfJ8PAHjikgHoGeKnOBGR81w+LA7n9uuB2iYj7vtqPcxmXp6xF5a7QiazxL1frUd9swmTBkXj4iGxqiMROZUQAs9MHoiIbj74s/gg3vtjh+pImmGz3IUQWUKITCHE9FMcT7c+Jsv+8bRt9q/FyN1lmYX65CVc7ZE8U3g336OzV2ct3orC8lrFibShzXI/UthSyhwAVUKIzJM87GYp5RwABiGEwQEZNalgXw1eWroVADAraxBCAzgLlTzXOWk9cEVGPJqNZtz39QaYeHmmy2y9ch8OoNj652IA6ccfFEJkA8gTQhiklLOklMWtP4EQIlsIkSuEyK2oqLBLaHfXYjLjvq/Xo8UkcfWIBJzVN0p1JCLlZkxKQ3SIH9aVVOHd30+oEuogW+Ue2urt8FZvJ1n/qxRCvC2EaP14SClnSykzpJQZkZGRXYiqHbN/Lcam0hrEhvrj4YlpquMQuYRgP288PXkgAOCFJdtQWH5YcSL3ZqvcqwCE2XhMkZSyCkAegGy7pNKwbftr8WrOdgDAs5cN5CbXRMc5u28ULh8Wh2ajGdPnrOflmS6wVe6rcezVuwHA0pMcPyIUll8GdApGkxn3z9mAZpMZVw6Px5g+/JcMUWuPTOqHHsG+WLO7Cu9z9EyntVnux90ozbS+nQMAQoilxx0PPe74bMfGdW/v/r4D60uqEB3ih4cv4OUYopMJ8ffGM9bLM88v3oqiCl6e6QybQyGtN0pzji9uKeX4kxyf5aiQWlBUcRgvLrXsAP/05IEI9vNWnIjIdY1L7YHL0uPQZDRj+hyOnukMTmJyApNZYvqcDWg2mnFZehzO5ugYIptmWi/P5O06xMszncByd4IPlu9EnnWy0sxJ/VTHIXILIQF/vTyz40Cd4kTuheXuYDsP1OH5xZYdZ566dCBCAng5hqi9xqX2wOShsWgymvHQ3A1cGrgDWO4OJKXEg3M3oLHFjIuHxGB8vx6qIxG5nUcn9UN4oA9WFFfiq9wS1XHcBsvdgb7O3YMVxZUIC/TBYxf2Vx2HyC11D/TBzAstlzOfWrAF5bWNihO5B5a7g1TUNuGphVsAWG4MhQVy7RiizrpocAzO6huJmkYjHv8uX3Uct8Byd5Anvs9HdUMLxqZE4uIhMarjELk1IQSevGQAAnz0WLCxDDn5+1VHcnksdwf4aWs5vlu/F37eOjx1yQDurERkB3HdA3DfuZaNtR+dvwm1jS2KE7k2lrud1Tcb8ci8TQCAe8anID4sQHEiIu244fREDI4PRVl1I55fvFV1HJfGcrezl5ZsQ2lVA/rHBGPq6N6q4xBpil4n8OzkgfDSCXy8Yhfydh1SHcllsdztaOOearz3xw7oBPDs5EHw0vP0EtlbWnQwbj7TACmBh+ZaZn7Tidg+dmI0mfHg3A0wS2Dq6N4YGBeiOhKRZt0xrg96RwRi2/7DeOuXItVxXBLL3U7e/2MnNu+1bMBx9/gU1XGINM3PW4+nL7UsTfD6skKuHHkSLHc7KKmsx0vWFR+fvHQAArkBB5HDjUoKx5SMODSbzJg5fxOXJmiF5d5FUkrMnL8JDS0mXDg4his+EjnRgxPSEBrgjT8KD+K79XtVx3EpLPcuWpK/Hz9trUCQnxcencQNOIicKSzQBw9NSAUAPLlgC2o49v0olnsX1Dcb8fh3mwEA95/XF1FBfooTEXmey4fFY1iv7qiobcKLHPt+FMu9C179cTv2VjdiQGwwrhnRS3UcIo+k01mWJtBbx75v2MOtnAGWe6dt21+Ld3/bASGAJy8ZCL2OSwwQqZIWHYypoxNhlsCMeZu4LR9Y7p0ipcQj326C0Sxx9WkJGBIfqjoSkce7KzMF0SF+2FhajU9X7lIdRzmWeyfMXVOKVTsqER7og+nnpaqOQ0QAAn298Jh13ffnf9jq8eu+s9w7qLq+BU9b12l/eGIat80jciHn9e+Js/tGorbJiKcWbFEdRymWewc9v6QAB+uacVpiGCanx6qOQ0THEULg3xcPgK+XDvPX7cUfhQdUR1KG5d4B60uq8OnK3fDSCTzBddqJXFJ8WADuPKcPAODRbzehyWhSnEgNlns7mcyWm6hSAn8/ozf69gxSHYmITmHaGAOSIgNRfKAO7/6+Q3UcJVju7fTpyl3YWFqNmBC/o68KiMg1+Xjp8K+LLJvSv76sEGXVDYoTOR/LvR0OHm7CC9aZbzMv7MeFwYjcwJg+kZgwoCfqm014emGB6jhOx3Jvh+cXb0VNoxFjUyJxXv+equMQUTvNuCANft46/G/9XvxZdFB1HKdiuduwYU8VvswtgZdO4LEL+/EmKpEbiesegFvPSgYA/Ou7zWgxec6uTSz3NpjNEjPnb4aUwNQzeiMpspvqSETUQdljDUgIC8DW/bX4+E/PmbnKcm/DN2v2YF1JFSKDfHHHuGTVcYioE/y89Zg5yTJz9eWl21BR26Q4kXPYLHchRJYQIlMIMd3G456zXyz1ahpb8NwPlpuoD01IRZAfZ6ISuatz0qKOzlx97gfPuLnaZrkLIbIAQEqZA6BKCJF5isdlAjDYP546r+Vsx4HDTRjWqzsuHcqZqETuTAiBmRf2h49ehzl5e7Bm9yHVkRzO1iv34QCKrX8uBpDe+gFCCMNxj9GEwvJafLB8J4QAHr+oP2+iEmlA74hA3DSmNwDgsfmbNb8ssK1yb72WbfhJHmOQUp6y3IUQ2UKIXCFEbkVFRYcDOpuUEv/6Lh9Gs8RVpyVgQGyI6khEZCe3j0s+uizwl6tLVMdxKFvlXgUg7FQHhRCZ1ks2pySlnC2lzJBSZkRGRnYmo1Mt3rwPvxceQIi/N+47t6/qOERkRwE+XphxgWWv41mLC3CorllxIsexVe6rcezVuwHA0lbHK603W7MAGIQQJ1y2cScNzSY88b1lmdB7z01BWKCP4kREZG8XDIzGKEM4qupb8OJS7e652ma5SynnwFLamda3cwBACLHU+vYa6/vCcOIlHLfz9q9FKK1qQGrPIFx9WoLqOETkAEIIPH5xf+h1Ap+t3I0tZTWqIzmEzaGQUspZUsocKeXs4943vtVjZkspk6SUaxwR0hlKKuvx5s9FACw3Ub30nAJApFUpPYJw7YgEmCXw7//lQ0rt3Vxlg1k9vXALmoxmXDQ4BiMMJ7tvTERacvf4FIQGeOPP4oNYvHm/6jh2x3IH8GfRQSzatA/+3no8NJF7ohJ5gtAAH9ydmQLgyIs7bW3q4fHlbjJL/Pv7fADArWclITrEX3EiInKWa0YkoE9UN+yurMd7v+9UHceuPL7cv84twZayGsSE+GHaWE1NsiUiG7z0OjxqXXfm9WXbUV7bqDiR/Xh0udc2tuCFJZahUA9OTIOft15xIiJytrEpkchMi0Jds+nopjxa4NHl/t+finDgcDPSE0Jx4aBo1XGISJEZF/SDt17g67w92FRarTqOXXhsue8+WI/3rBvnPnYh148h8mS9IwJx4+mJkBJ4/H+bNTE00mPL/emFW9BsMmNyeiwGx7v9/Csi6qI7zumD8EAfrN55CAs2lqmO02UeWe5/Fh3ED5stQx+nn8ehj0QEBPt5417relLPLCxAY4t7D430uHI3mSWesA59vOWsJPQM8VOciIhcxRXD45EWHYzSqgbM/tW9VzL3uHKfk1eC/CNDH8dw6CMRHaPXCTx2oWVo5Js/F2FftfsOjfSocq9tbMHzi7cBAB6YkAp/Hw59JKK/GmkIx4QBPdHQYnLrLfk8qtwtQx+bkJ4QiosGx6iOQ0Qu6uGJafDx0mHe2lKsddMt+Tym3I8f+jiTQx+JqA3xYQG46QzLlnxPLdjilkMjPabcn1lkHfo4NBZDOPSRiGy45awkRHTzQe6uQ1i0aZ/qOB3mEeW+ovjYqo/3n8+t84jItiA/b9xlXTXy2UUFbrdqpObL3Xzc0Md/nMlVH4mo/a4cHn901ciP/9ylOk6HaL7c560txea9NYgO8UM2V30kog7w0uvw8ETLhtqv/bjdrTbU1nS5NzSb8Lx1lbf7zu3LoY9E1GFn9Y3EmD4RqGk04rVl21XHaTdNl/u7vxdjX00j+scE49KhsarjEJEbEkLg4YlpEAL4+M9dKK44rDpSu2i23Ctqm45ueD1jYhp0Og59JKLOSYsOxpRh8TCaJZ5d5B4TmzRb7q/kbENdswnnpEbh9OQI1XGIyM3de24KAnz0WJK/HyuKD6qOY5Mmy337/lp8vmo39DrBDa+JyC6igv1w89gkAJaJTWaza09s0mS5P7OoAGZpGcaUHBWkOg4RacS0sb3RI9gXG0urMX99qeo4bdJcuf9ReADLCsrRzdcLd49PUR2HiDQkwMcL91v3gJj1w1Y0NLvuxCZNlbvZLPHUgi0Ajkwd9lWciIi0ZvLQWPSPCUZZdSPe/d1113zXVLnPW1uK/DLLhKWpo3urjkNEGqTTCcy4wDKx6c2fi1Be65prvmum3DlhiYic5fSkCGSmRaGu2YSXl7rmxCbNlDsnLBGRMz04IQ16ncCXq3dj675a1XFOoIly54QlInK25KhuuGZEAswSeGrhFtVxTqCJcueEJSJS4a7MFAT5euHXbRX4bXuF6jh/4fblzglLRKRKWKAPbjnbMrHpmYUFLjWxyWa5CyGyhBCZQojpJzkWaj2eJYR4zjER28YJS0Sk0tTRvREd4of8shqXmtjUZrkLIbIAQEqZA6BKCJHZ6iFTAIRJKedYH5/tkJSnwAlLRKSan7ce955r2eHthcXb0NjiGhObbL1yHw7gyCj9YgDpxx+UUs6WUs62vmkAkNP6EwghsoUQuUKI3IoK+12T4oQlInIVlw6NRWqhf3CnAAAJaklEQVTPIJRWNeCjP3eqjgPAdrm33kk6/GQPEkIYAFRKKU+YrmX9BZAhpcyIjIzsZMwTzeWEJSJyEZZ7fpaJTa8vK0RVvfodm2yVexWAsHZ8niwp5c12yNMujS0mvLiEE5aIyHWM7ROBM5ItOzb996dC1XFslvtqHHv1bgCwtPUDhBBZUspZ1j+3vibvEB8u34my6kakRXPCEhG5BiEEHpxgGbH34fJdKKmsV5qnzXK33ig1HClt641VCCGWWv+fCeA5IUSeECLP0WEBoKq++ehvxQcnpHLCEhG5jAGxIbh0aCyaTeajVxdU8bL1gCOvylu9b7z1/zkAkhyQ65Te+LkINY1GnJ4UjrF9OGGJiFzLveemYMGGMny7bi9uGmPAgNgQJTncahJTaVUDPli+EwDw0IQ0CMFX7UTkWuK6B+DG0YkAgKcXboGUaiY2uVW5v7hkK5qNZlw4OAYD49T8NiQisuW2s5IR4u+N5UUH8cs2NcsSuE255++twby1pfDWC9xvnTBAROSKQgK8cfvZyQCAZxcVwKRgWQK3KffnfiiAlMA1I3ohITxAdRwiojZdN6oXYkP9UbCvFnPX7HH687tFuS8vPIBftlWgm68X7hiXrDoOEZFNft563H+e5SrDi0ucvyyBy5e72Szx7A8FAIB/nGlAOJcZICI3cdHgGPSPCca+mka898cOpz63y5f7go1l2LCnGlFBvph6BpcZICL3odMJPGxdluDNn4pw8HCT857bac/UCc1G89F9Ue/KTEGAj81h+URELmV0cgTOTIlEbZMR/1nmvGUJXLrcP1u5C7sr65EUGYgpGXGq4xARdcqDE1IhBPDJil3YeaDOKc/psuVe29iC16y/5aafnwovvctGJSJqU1p0MC5Lj4PRLPG8k5YlcNnGnP1rMSrrmjGsV3ec26+H6jhERF1yz/gU+HrpsGBDGdbuPuTw53PJci+vacT//Wa5s/zwxFQuM0BEbi8m1P/ooJBnFhU4fFkClyz3V37cjoYWE87t1wPDerVnOXkiItd3y1lJ6B7gjVU7KrGsoNyhz+Vy5V5UcRhfri6BTgDTz+cyA0SkHcF+3rh9XB8Alln3jlyWwOXKfZb1C75ieDySo4JUxyEisqtrRyYgNtQf2/YfduiyBC5V7nm7KrF48374eetwV2aK6jhERHbn66XHfedZ+u2lpY5blsBlyl1KiWcWWpYZuOkMA3oE+ylORETkGBcPjkVadDDKqhvxoXWPCntzmXLP2VKO3F2HEBbog5vPNKiOQ0TkMDrdsf1W//tTIarrW+z/HHb/jJ1gNJnxnHVxsDvGJSPIz1txIiIixxrbJwKnJ4WjptGIN36x/7IELlHuc/L2oLD8MBLCAnDNiF6q4xAROZwQx169v//HTuytarDr51de7g3NJrycsw0AcN95feHjpTwSEZFTDIoLxaRB0Wg2mvHy0m12/dzKm/S9P3Zgf00TBsaGYNLAaNVxiIic6r5z+8JLJ/DNmj3Yuq/Wbp9XablX1jXjrZ+LAFhWTdPpuMwAEXmWxIhAXD0iAWZpmedjL0rL/T/LtqO2yYixKZEYnRyhMgoRkTJ3jOuDQB89fiwox8rig3b5nMrKvaSyHp+s2AUhgAfPT1UVg4hIucggX0wbaxkC/uwP9llUTFm5v7BkK1pMEpcOiUW/mGBVMYiIXMJNYwyI6OaDtbursHjz/i5/PiXlvqm0GvPX7YWPXod7zuUyA0RE3Xy9cOc5lkXFZi0ugNFk7tLnU1Luzy6y3DS44fReiOseoCICEZHLueq0BCSGB6C4og5f5XZtUTGnl/uv2yrwe+EBBPl54dazkp399ERELstbr8N951mWOn8lZxvqm42d/lxOL/cjr9pvPSsZ3QN9nP30REQubeKAaAyKC0F5bRPe+31Hpz+PU8u9qr4F+WU1iA7xw99GJzrzqYmI3IJOJ46OIHzrF8te0p36PLYeIITIEkJkCiGmd+b48fbVNAIA7h6fAj9vfYfDEhF5gtOTI3BmSiQONxnx+rLOLSrWZrkLIbIAQEqZA6BKCJHZkeOttZjM6NsjCJelx3UqLBGRp3jg/FQIAXy8YmenPt7WK/fhAIqtfy4GkN7B4yd4YEJf6LnMABFRm/rFBOOSIbFoMXVuQpOtcg9t9XZ4B49DCJEthMgVQuT6CImz+0Z1NCMRkUe6Z3wKfDu5Uq6tj6oCENaF45BSzpZSZkgpM5J6hkAIvmonImqP+LAA/P7AuE59rK1yX41jr84NAJZ28PhfePFyDBFRh0QG+Xbq49osdynlHACGIzdKrTdOIYRY2tZxIiJSy8vWA6SUs07yvvFtHSciIrWU78RERET2x3InItIgljsRkQax3ImINIjlTkSkQcIee/W1+8mEqAWw1WlP6NoiABxQHcJF8Fwcw3NxDM/FMX2llEEd+QCbQyHtbKuUMsPJz+mShBC5PBcWPBfH8Fwcw3NxjBAit6Mfw8syREQaxHInItIgZ5f7bCc/nyvjuTiG5+IYnotjeC6O6fC5cOoNVSIicg5eliEi0iCWOxGRBjms3O25sba7a+trFUKEWo9nCSGeU5HPWdr7d6718wC06+cj/cj3hbOzOVsHuiLb2dmczfp1nnJfjI70pkPK3d4ba7uzdnytUwCEWdfGh1a/gdv7d259v8GZ2Zytnefi5uP2S9Ds+WhHV2QCKLYeLxZC2Nyn2Z21tSdGR3vTUa/c7b6xthtr82u1bkN45E64AYBWNzyx+XduLbHi1u/XoDbPhfUXfJ4QwiClnCWl1PI5sfV9kQvga2upG6SUa5wZzsV0qDcdVe5d3lhbQ9r1tVqLrVLDP8jtOQ8GDX/9x7N1LpKs/1UKId4WQrR+vJa0eS6klFUA3gbwNSznxJN1qDcdVe5d3lhbQ9r7tWZJKW92dBiF2jwPQohMD9qmsT3fE0XWYssDoMlLdVa2vi+yAORIKZOOe9tTdag3HVXudt1Y283Z/FqFEFlHtivU8P0HW+eh0nqjKAuW68xavlTXnp+PI0Jh+aHWKlvn4vhLMc/Ac14UnkyHetMh5c6NtY+xdS6s739OCJEnhMhTl9Sx2vE9scb6vjCc+M9PTWnnz0foccc1O1PT1rkAMFsIkW09PkXL5wI4+i+TjOP/hdLZ3uQMVSIiDeIkJiIiDWK5ExFpEMudiEiDWO5ERBrEcici0iCWOxGRBrHciYg0iOVORKRB/w85j6/KPI3hNAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.figure() \n",
"plt.plot(xx,ff(xx))\n",
"plt.savefig('fig0.png')"
]
},
{
"cell_type": "code",
"execution_count": 24,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEICAYAAACktLTqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xd4VFX+BvD3zKSHFNIgFZiEkNAJQUAEFYMKdozYy7ISVlddK4q4uK4du+ta8GfviqCuoEjErpSEXgIkoYRQEggpkDoz5/fH3EAMkEmZmTNz5/08jw9J7s3Md8bw5nLuOd8jpJQgIiJ9MagugIiIHI/hTkSkQwx3IiIdYrgTEekQw52ISIcY7kREOsRwJyLSIYY7KSWEyBBCZAshMlp9Pbsd3xsuhMhyXnUd15662/het3s95LkY7qTadO1PU/MXhBA5Usp59r5RSlmpnW+yd25nCSGKOnBuu+o+GVe8HvIeDHdSRguxSinlvOZQFEKEd+QxpJS5OPYLwhmGt+ekjtZ9Mi54PeQlfFQXQN5JC/YnARQLITKklKu0QzkA5rY6NwtAc3hWAsiQUs5xQj0ZAIoBXC6lvFcbKvoMQLJWw2cAzgIQAeAyKWXLED5R3a/B9i+SSu3PYgCvaYed+nqIGO6khJSyWAgBAI83D0doklt9DgDhUsp5QoglAC7D8Ve2Ra1+QTRfSU9po4RcKWVxi8+nN9fSPCwipVwlhCjWPs4VQuS1+HiCECJLu9I+rm7t2HQhRLZW+9FzW3yt3a+HqKMY7qSS6QRBHtHyEyFEeMtxbO38y1p9TwVsV8arWp03F+33CYB8IUQujl1dn0jzL4SDbdXdIvSbx8/Dgc6/HqKO4pg7qVRh74QWNxmzACzRPj7R+HbrXxIdVQzb+Ho+bMMvbdbTHq1ujDb/a8BVr4e8HK/cSQkt+Irbcd4MALkAJkALQwBZAFrOSolo/VhaYOa08dDzWg3LzJRS3gtgrjYl0dTqeGdkA8jVahmh1dWp10PUUQx3UsUE21Vyayu1oYvmK9di2MLucQAztdBtPdyS3Ppr2vd35CblQW2OeiVsM3iKtRuqJiFEDoA87eMZsAXxBNhutOZpz9W6bgAY0XyjVLu/0OnXQ9RRgpt1kArNIdn66li7os/qSLgJIZ7UrrqV6UzdbTyW8tdDno9X7uRS2lizCbar0+OGHlrMounI47V1A9QlOlr3ybjL6yHPxxuq5GrFsN1Ibeum5dz2LMNvvhHpgLFxh2hv3Sfjbq+HPBuHZYiIdIhX7kREOuTSMfeoqCjZu3dvVz4lEZHHy8/PPyCljO7I97g03Hv37o28vDxXPiURkccTQuzs6PdwWIaISIcY7kREOsRwJyLSIYY7EZEOMdyJiHSI4U5EpEN2w10IkaXtGHOy49naOTMcWxoREXWW3Xnu2pZiJ+xQp7VIbT7H1GrbMSLlpJQ4eKQRuw/VofRQHXYfqkV5TQNq6s2oaWhCTb0ZFquEQQgIARgNAiEBvggP9EV4kC8ig/2QFBmEpIggJHQPQoCvUfVLImqXri5iGgHb9mSArSFUBmwbERyl9cLOAYCkpKQuPh1R2/ZX12P59gpsLK3Cpr3V2LSnGgePNDrs8ZMigjAoPgwD48MwKD4Mw5LCEezP5qrkfrr6U9l6e7DI1ido/a3nAkBmZia7lJFDHW4w46ct5fi1sBzLiiuw/cCR484J8fdBYkQQ4rsHIj48ED1CAxAa6IPQAF90C/CBr8EACQmrBMwWK2rqzaisbURlXRPKahpQUlGLXRW1KD1Uh13axwvX7wUA+BgEhiaG49SUKIztG4WMpO4wGrre+peoq7oa7pVotTEwkbMdOtKIbzfuw3cb9+G3woNotFiPHgv2MyKzdwSGJYWjf2wo+seFIj48EI7otd5ksaKw7DA2lFZhQ2kV1pRUYn1pFfJ2HkLezkN48fttiOrmh7MH9MTEgT0xyhQJXyPnLJAaXQ33lTh29W7CsT0hiRzKbLHil20H8Fl+CXI3lR0NdCGAzF7dcVZ6D4xOjsTAuFD4OClQfY0GpMeGIj02FJdlJgIAquubsLy4Ar8VHsDSgjLsqqjFh8t34cPluxDVzQ+XZiTgssxEpMR0c0pNRCdjt5+7dtP0dQDTpJTztK8tkVJO0D6eAWAVgBPtBfknmZmZko3DqCMOHm7Ah8t34f3lO7G/ugGALdBPS4nCpEGxyErvgegQf8VV2kgpsWlvNb7dsA8L1+9FcfmxIaLMXt1x7ehemDQollfz1GFCiHwpZWaHvseVm3Uw3Km9tu2vwRu/bseC1aVoMNuu0k1Rwbh0eAImZ8QjNixQcYVtk1JidUklPl1Zgv+t3YMjjRYAQHx4IP4ypjcuH5GIkABfxVWSp2C4k8fbtr8Gz3+/DYvW70Xzj+ZZaTGYelofnJoc6ZCxc1c70mDGF2tK8cav249ezYf4+2DqaX3w17F9EMqQJzsY7uSxCssO44Xvt+HrdXsgJeBnNGDKiARMHdMHpmh9jFdbrRJLC8ow95dirNheAQAIC/RFzjgTbji1N6dU0kkx3MnjHDrSiOdzt+L95btgsUr4GgWuGJGEm89Mdvuhl65YXnwQzyzZejTko7r5Y8Y5/ZA9PAEGTqWkVhju5DGaLFa8v2wnns/dhqq6JhgEcPmIRNwyvi/iw/Ub6i1JKfF70UE8/d0WrN5VCQAYFB+G2Rf0x4jenGFMxzDcySOs3FGBmfPXo7DsMADbzJcHzk9HWs9QxZWpIaXEV2v34IlvCrC3qh4AcOGQOPzz/P5uMxOI1GK4k1urrm/Ck98U4IPluwAAfaKCMWtSOs5Kj/HIG6WOVttoxqs/FeO1n4rQYLYiLNAX/zy/Py7NiOf74+UY7uS2Fm/ch9lfbsD+6gb4GgVuOj0ZN5+ZwkZcJ1BSUYv7F6zHL9sOAADG9o3CY5cMQmJEkOLKSBWGO7mdmvomPPS/TZiXvxsAMCwpHE9MHox+PUMUV+bepJSYv6oUDy/chMraJgT5GfHgBf0xJTORV/FeiOFObiVvRwXu+HQNSirq4O9jwH0T03Dd6N5srNUBBw434MGvNmLhOlujsnMH9MTjkwehe7Cf4srIlRju5BbMFite/H4bXvqhEFYJDIgLxQtXDEVKDK/WO0NKiQWrSzH7y4043GBGTIg/npkyBGP7RqsujVykM+HOJhfkUOU1Dbj2jRV4cWkhJICbzkjGgpvHMNi7QAiByRkJ+OYfY5HZqzvKahpw3Zsr8NySrbBY2UWbTozhTg6zYnsFznvxF/xRfBBR3fzx4Y2jcO+5afDz4Y+ZIyRGBOHjnFG4IysVAPDC99tww1srUOHAzUhIP/i3jrpMSonXfy7Gla8vQ1lNA07pHYFFt52G0cnH7d1CXeRjNOAfWX3x7tRTEBHsh1+2HcD5L/6C1bsOqS6N3AzDnbqkvsmCf3y8Bo8u2gyLVWL6OBM+nDYSMaEBqkvTtbF9o/H1radhWFI49lTV4/LXlmH+qt2qyyI3wnCnTiurqccVc5fhq7V7EOxnxKvXDMfMSelO2yyD/iwuPBCf5IzGtaN6odFixZ2frsWcbwtg5Tg8geFOnbRpTzUufuk3rCmpRHx4IObddCrOHdhTdVlex8/HgIcvHoiHLxoAo0Hg5R+LcNMH+ahtNKsujRRjuFOH5W7aj+xXf8eeqnpkJIXji7+PQXqsd/aFcRfXju6Nt/8yAiEBPli8cT+yX/kD+7Q+NeSdGO7UIR+t2IWc9/JQ22jBxUPj8OG0UWxu5SbG9o3GgpvHoFdkEDbtrcalr/x+tDkbeR+GO7WLlBIvfr8NM+evh1UCt53VF89dPpS9YdxMSkw3fHHzGGQkhaO0sg7Zr/6OVZxJ45UY7mSXxSox+8uNeHbJVggBPHLxQNw5IZU9TtxU92A/fHDjKIxPi0FlbROuen0ZfigoU10WuRjDndrUaLbito9W471lO+FnNODlqzJwzaheqssiOwL9jHjt2uHIHp6A+iYrbnw3D5/nc6qkN+GmjXRS9U0W3PzBKiwtKEOIvw9evz4To0xcmOQpfI0GPJU9GDEh/nj5xyLcPW8tGi1WXHlKkurSyAV45U4nVNdowbR387C0oAzdg3zxUc4oBrsHEkJgxrlpmDkxDVICM+evxzu/71BdFrkAw52Oc6TBjL+8vQK/bDuAqG5++ChnFAbGh6kui7pg+unJmH1+fwDAg19txOs/FyuuiJyN4U5/UlPfhOvfXIFlxRWICfHHxzmjvHZvU72ZelofPHLxQADAo4s2478/FCquiJyJ4U5HHW4w4/o3VyBv5yHEhgXgk+mj2apXZ64Z1QtzLh0MIYCnFm/BS0u3qS6JnIThTgBsY+x/fXslVu2ytRP4dPpo9IkKVl0WOcGUEYl45rIhEAJ4+rut+L9fOESjRwx3QoPZgpz38rB8u20o5sNpI7kZs85NzkjAk5MHAwAeWbgZ7y/bqbgicjSGu5drsljx9w9W45dtBxAZ7IcPp41Er0hesXuDKSMS8e+LBgAAHvhiw9FNzEkfGO5ezGKVuP2TNcjdvB9hgb54768jOcbuZa4b3Rv3T0oDAMyYtxZfr9ujuCJyFIa7l5JSYtaC9Vi4bi+6+fvg3amnoH8cZ8V4o5xxybgjKxVWCdz+8Rr8uIWtCvSA4e6lnluyFR+vLIG/jwFv/WUEhiSGqy6JFLrtrBRMH2eC2Spx0/ursKakUnVJ1EV2w10IkS2EyBJCzLBzPMfx5ZEzvLdsJ15cWgiDAF66KgMjekeoLokUE0LgvolpuDQjAXVNFkx9eyWKy9ku2JO1Ge5CiGwAkFLmAqgUQmS1Op4FoFg7XiyEyHBapeQQ327Yi9lfbgAAPHbJIEzo30NxReQuhBB44tJBOLNfNCqONOLaN1ZgfzU3/PBU9q7cRwBongRbDKB1eOcB+EwLdZOUclXrBxBC5Agh8oQQeeXl5V0umDpvWfFB3PbxGkgJ3DUhFVewgRS14ms04L9XZ2Booq0f/PVvrkB1fZPqsqgT7IV764HYP3WOklJWAngNwGcAkk/0AFLKuVLKTCllZnR0dKcLpa7Zsq8G097NQ6PZimtH9cIt41NUl0RuKsjPB2/eMAKm6GAU7KtBjvZzQ57FXrhXAjjpgKw2bJMrpUxu8Tm5mbKaevzlrRWoqTfj3AE98a8LB3CjDWpTRLAf3p16CqJD/LGsuAKzFqyHlFJ1WdQB9sJ9JY5dvZsALGl1vOVQzONo4xcBqVHfZMG0d/Oxp6oew5LC8fwVQ2E0MNjJvoTuQXjj+kwE+BrwWf5uvPJTkeqSqAPaDHcp5TwApuYbqdqNUwghmkN+rjamngVgipRyrlOrpQ6xWiXu+nQt1pbY+sXMvTaTe55ShwxOCMcLVwyDEMCcb7dg0fq9qkuidrK7E5OUcs4JvjZB+7MSAAPdTT27ZCsWrrctUnrzhhGIDvFXXRJ5oHMG9MTMiWl4bFEB7vhkDeLCAzGU6yLcHhcx6dTn+bvx0g/Nc9mHoV9PthWgzps21oQrT0lEg9mKG9/Jw+5DtapLIjsY7jq0ckcF7pu/DgDwrwsH4Ix+MYorIk8nhMC/LxqIMSmROHC4ATe+k4cjDWbVZVEbGO46s6eyDje9n48mi8QNp/bGdaN7qy6JdMLXaMDLVw8/OkXy7s/WcgaNG2O460h9kwV/ez8fBw43YkxKJB44L111SaQzYYG+eP26TIT4++CbDfvw8o+cQeOuGO46YevyuAHrdlchoXsgXroyAz5G/u8lx0uO7oYXrhyq7eS0BUsL9qsuiU6Af/t14p3fd+DzVbsR4GvA3Gsz0T3YT3VJpGPj03rgrgmpkBL4x0drUMQmY26H4a4DfxQdxMMLNwMA5mQPYV92com/n5mCiQN7oqbBjGnv5rEHjZthuHu40so6/P3DVbBYJaafbsKFQ+JUl0ReQgiBpy8bgn49QlBcfgR3frIGVitvsLoLhrsHq2+y4G/v5aPiSCPG9o3CjHPSVJdEXibY3wdzrxuOsEBf5G4uw6s/8waru2C4e7BHFm7C+tIqJEYE4j9XDmPPGFKiV2Qwnrt8CADg6cVb8EfRQcUVEcBw91hfrinF+8t2wc9owMtXDUd4EG+gkjrj03rg72cmwyqBWz9ajTJu8qEcw90DFZYdxsz56wEAsy/oj0EJYYorIgLuyErFaJNtBestH62G2cIe8Cox3D1MbaMZN3+Qj9pGCy4cEoerR3I3JXIPPkYDXrhyKGJC/LFiewWe/m6r6pK8GsPdg0gp8cAXG7B1/2GYooPx2ORB3HSD3EpMSMDR+z+v/lSEJZu4wEkVhrsH+TSvBPNXlSLA14BXrh6Obv52OzYTudxIUyRmnNMPAHDnp2tQUsEOkiow3D3Epj3VmP3lRgDAIxcPYgtfcms540zISu+Bmnozbvt4NZo4/u5yDHcPUNtoxi0frUKD2YopmQnIHp6guiSiNtkWOA1GbFgAVu+qxPO5HH93NYa7B/j3/zahuPwI+sZ0w0MXDlRdDlG7hAf54fnLh8IggJd/LMLvhQdUl+RVGO5ubtH6vfh4ZQn8fAz4z1XDEOjHPVDJc4w0ReLW8X0hJXD7J2tw8HCD6pK8BsPdjZVW1uG+z207Ks2alI60nmwIRp7n1vEpOKV3BMpqGnDPvHXc4MNFGO5uymyx4vaPV6O63oys9BhcN7qX6pKIOsXHaMDzVwxFWKAvlhaU4a3fdqguySsw3N3USz8UYuWOQ4gJ8cec7CGcz04eLS48EE9eOhgA8MQ3BdhQWqW4Iv1juLuhvB0VePH7bRACeO7yoYjgxhukA+cO7IlrRiWh0WLFbR+tRm0jN9h2Joa7m6mqa8I/Pl4Dq7TNFR6TEqW6JCKHeeC8/kjt0Q3FB47g8UUFqsvRNYa7m/nXVxtRWlmHwQlhuGtCP9XlEDlUgK8Rz10+FL5GgfeW7cSPW8pUl6RbDHc3smj9XixYbWsv8PzlQ+Hnw/89pD8D4sJwp3bhMmPeOhw60qi4In1ieriJsup6zFpga+M7a1I6TNHdFFdE5Dw540wY0bs7ymoaMOuL9Zwe6QQMdzcgpcS9n6/DodomjO0bhWtGcdoj6ZvRIPDslKEI9jNi0fp9+GJNqeqSdIfh7gY+XlmCH7aUIzTAB09x2iN5icSIIDx44QAAwOwvbPeayHEY7ortPHgED3+9CQDw8MUD0TMsQHFFRK5z2fAEnN2/B2oazLj707WwWjk84ygMd4UsVom7Pl2L2kYLzh8ci4uGxqsuicilhBB4fPIgRHXzwx/FB/Hmb9tVl6QbdsNdCJEthMgSQsw4yfEM7Zxsx5enb3N/LkbeTtsq1EcuZrdH8k6R3fyPrl6ds3gLCstqFFekD22Ge3NgSylzAVQKIbJOcNp0KeU8ACYhhMkJNepSwb5qPLtkCwBgTvZghAdxFSp5r7PSe+DyzEQ0mq24+7N1sHB4psvsXbmPAFCsfVwMIKPlQSFEDoB8IYRJSjlHSlnc+gGEEDlCiDwhRF55eblDivZ0TRYr7v5sLZosEleNTMIZ/WJUl0Sk3Kzz0xEbFoA1JZV449fjooQ6yF64h7f6PLLV58nafxVCiNeEEK3Ph5RyrpQyU0qZGR0d3YVS9WPuz8XYUFqN+PBA3D8pXXU5RG4hNMAXj00eBAB4+rutKCw7rLgiz2Yv3CsBRNg5p0hKWQkgH0COQ6rSsa37a/BC7jYAwBOXDuIm10QtnNkvBpcNT0Cj2YoZ89ZyeKYL7IX7Shy7ejcBWHKC483CYftlQCdhtlhxz7x1aLRYccWIRIzty3/JELX2wPn90SPUH6t2VeItzp7ptDbDvcWN0izt81wAEEIsaXE8vMXxuc4t17O98et2rC2pRGxYAO4/j8MxRCcSFuiLx7XhmacWb0FROYdnOsPuVEjtRmluy+CWUk44wfE5zipSD4rKD+OZJbYd4B+bPAihAb6KKyJyX+PTeuDSjAQ0mK2YMY+zZzqDi5hcwGKVmDFvHRrNVlyakYAzOTuGyK7Z2vBM/s5DHJ7pBIa7C7z9+w7ka4uVZp/fX3U5RB4hLOjPwzPbDxxRXJFnYbg72Y4DR/DUYtuOM49eMghhQRyOIWqv8Wk9MHlYPBrMVsycv46tgTuA4e5EUkrcN38d6pusuGhoHCb076G6JCKP88/z+yMy2A/LiivwaV6J6nI8BsPdiT7L241lxRWICPbDgxcMUF0OkUfqHuyH2RfYhjMfXbgZZTX1iivyDAx3JymvacCjizYDsN0Yighm7xiizrpwSBzO6BeN6nozHvpqk+pyPALD3Uke/noTquqaMC41GhcNjVNdDpFHE0LgkYsHIsjPiIXr9yJ3037VJbk9hrsT/LClDF+t3YMAXwMevXggd1YicoCE7kG4+2zbxtr//HIDauqbFFfk3hjuDlbbaMYDCzYAAO6ckIrEiCDFFRHpx/Wn9saQxHDsrarHU4u3qC7HrTHcHezZ77aitLIOA+JCMXVMH9XlEOmK0SDwxORB8DEIvLdsJ/J3HlJdkttiuDvQ+t1VePO37TAI4InJg+Fj5NtL5GjpsaGYfroJUgIz59tWftPxmD4OYrZYcd/8dbBKYOqYPhiUEKa6JCLdunV8X/SJCsbW/Yfx6k9FqstxSwx3B3nrtx3YuMe2AccdE1JVl0OkawG+Rjx2ia01wUtLC9k58gQY7g5QUlGLZ7WOj49cMhDB3ICDyOlGJ0diSmYCGi1WzP5yA1sTtMJw7yIpJWZ/uQF1TRZcMCSOHR+JXOi+iekID/LFb4UH8dXaParLcSsM9y76btN+/LClHCEBPvjn+dyAg8iVIoL9MHNiGgDgkYWbUc2570cx3LugttGMh77aCAC455x+iAkJUFwRkfe5bHgihvfqjvKaBjzDue9HMdy74IXvt2FPVT0Gxofi6pG9VJdD5JUMBltrAqM2933dbm7lDDDcO23r/hq88ct2CAE8cvEgGA1sMUCkSnpsKKaO6Q2rBGYt2MBt+cBw7xQpJR74YgPMVomrTknC0MRw1SUReb3bs1IRGxaA9aVV+GD5TtXlKMdw74T5q0qxYnsFIoP9MOOcNNXlEBGAYH8fPKj1fX/q2y1e3/ed4d5BVbVNeEzr037/pHRum0fkRs4Z0BNn9otGTYMZjy7crLocpRjuHfTUdwU4eKQRp/SOwOSMeNXlEFELQgj8+6KB8Pcx4Ms1e/Bb4QHVJSnDcO+AtSWV+GD5LvgYBB5mn3Yit5QYEYTbzuoLAPjnFxvQYLYorkgNhns7Way2m6hSAn89rQ/69QxRXRIRncS0sSYkRwej+MARvPHrdtXlKMFwb6cPlu/E+tIqxIUFHL0qICL35OdjwL8utG1K/9LSQuytqlNckesx3Nvh4OEGPK2tfJt9QX82BiPyAGP7RmPiwJ6obbTgsUUFqstxOYZ7Ozy1eAuq680YlxqNcwb0VF0OEbXTrPPSEeBrwP/W7sEfRQdVl+NSDHc71u2uxCd5JfAxCDx4QX/eRCXyIAndg3DzGSkAgH99tRFNFu/ZtYnh3garVWL2lxshJTD1tD5Iju6muiQi6qCccSYkRQRhy/4avPeH96xcZbi34fNVu7GmpBLRIf64dXyK6nKIqBMCfI2Yfb5t5epzS7aivKZBcUWuwXA/ier6Jjz5re0m6syJaQgJ4EpUIk91VnrM0ZWrT37rHTdX7Ya7ECJbCJElhJhh57wnHVeWei/mbsOBww0Y3qs7LhnGlahEnkwIgdkXDICf0YB5+buxatch1SU5XZvhLoTIBgApZS6ASiFE1knOywJgcnx5ahSW1eDt33dACOChCwfwJiqRDvSJCsaNY/sAAB78cqPu2wLbu3IfAaBY+7gYQEbrE4QQphbnHEcIkSOEyBNC5JWXl3e6UFeRUuJfX22C2Spx5SlJGBgfprokInKQW8anHG0L/MnKEtXlOJW9cG/dqDzyBOeYpJQnDXcp5VwpZaaUMjM6OrrDBbra4o378GvhAYQF+uLus/upLoeIHCjIzwezzrPtdTxncQEOHWlUXJHz2Av3SgARJzsohMjShmx0oa7Rgoe/trUJvevsVEQE+ymuiIgc7bxBsRhtikRlbROeWaLfPVfthftKHLt6NwFY0up4hXazNRuASQhx3LCNJ3nt5yKUVtYhrWcIrjolSXU5ROQEQgg8dNEAGA0CHy7fhc17q1WX5BRthruUch5soZ2lfZ4LAEKIJdrnq7SvReD4IRyPUlJRi1d+LAJgu4nqY+QsUSK9Su0RgmtGJsEqgX//bxOk1N/NVbsJJqWcI6XMlVLObfG1Ca3OmSulTJZSrnJGka7w2KLNaDBbceGQOIw0nejWAhHpyR0TUhEe5Is/ig9i8cb9qstxOF6eAvij6CC+2bAPgb5GzJzEPVGJvEF4kB/uyEoF0Hxxp69NPbw+3C1WiX9/vQkAcPMZyYgNC1RcERG5ytUjk9A3pht2VdTizV93qC7Hobw+3D/LK8HmvdWICwvAtHG6WYdFRO3gYzTgn1rfmZeWbkNZTb3iihzHq8O9pr4JT39nmwp136R0BPgaFVdERK42LjUaWekxONJoObopjx54dbj/94ciHDjciIykcFwwOFZ1OUSkyKzz+sPXKPBZ/m5sKK1SXY5DeG247zpYize1jXMfvID9Y4i8WZ+oYNxwam9ICTz0v426mBrpteH+2KLNaLRYMTkjHkMSPXqKPhE5wK1n9UVksB9W7jiEhev3qi6ny7wy3P8oOohvN9qmPs44h1MfiQgIDfDFXVo/qccXFaC+ybOnRnpduFusEg9rUx9vOiMZPcMCFFdERO7i8hGJSI8NRWllHeb+fNJ+iB7B68J9Xn4JNjVPfRzLqY9EdIzRIPDgBbapka/8WIR9VZ47NdKrwr2mvglPLd4KALh3YhoC/Tj1kYj+bJQpEhMH9kRdk8Wjt+TzqnC3TX1sQEZSOC4cEqe6HCJyU/dPSoefjwELVpditYduyec14d5y6uNsTn0kojYkRgThxtNsW/I9unCzR06N9Jpwf/wbberjsHgM5dRHIrLjpjOSEdXND3k7D+GbDftUl9NhXhHuy4qPdX2851xunUdE9oUE+OJ2rWvkE98UeFzXSN2Hu7XF1Md8q3EtAAAMCUlEQVS/nc6uj0TUfleMSDzaNfK9P3aqLqdDdB/uC1aXYuOeasSGBSCHXR+JqAN8jAbcP8m2ofaL32/zqA21dR3udY0WPKV1ebv77H6c+khEHXZGv2iM7RuF6nozXly6TXU57abrcH/j12Lsq67HgLhQXDIsXnU5ROSBhBC4f1I6hADe+2MnissPqy6pXXQb7uU1DUc3vJ41KR0GA6c+ElHnpMeGYsrwRJitEk984xkLm3Qb7s/nbsWRRgvOSovBqSlRqsshIg9319mpCPIz4rtN+7Gs+KDqcuzSZbhv21+Dj1bsgtEguOE1ETlETGgApo9LBmBb2GS1uvfCJl2G++PfFMAqbdOYUmJCVJdDRDoxbVwf9Aj1x/rSKny5tlR1OW3SXbj/VngASwvK0M3fB3dMSFVdDhHpSJCfD+7R9oCY8+0W1DW678ImXYW71Srx6MLNAJqXDvsrroiI9GbysHgMiAvF3qp6vPGr+/Z811W4L1hdik17bQuWpo7po7ocItIhg0Fg1nm2hU2v/FiEshr37Pmum3DngiUicpVTk6OQlR6DI40WPLfEPRc26SbcuWCJiFzpvonpMBoEPlm5C1v21agu5zi6CHcuWCIiV0uJ6YarRybBKoFHF21WXc5xdBHuXLBERCrcnpWKEH8f/Ly1HL9sK1ddzp94fLhzwRIRqRIR7IebzrQtbHp8UYFbLWyyG+5CiGwhRJYQYsYJjoVrx7OFEE86p8S2ccESEak0dUwfxIYFYNPearda2NRmuAshsgFASpkLoFIIkdXqlCkAIqSU87Tzc5xS5UlwwRIRqRbga8RdZ9t2eHt68VbUN7nHwiZ7V+4jADTP0i8GkNHyoJRyrpRyrvapCUBu6wcQQuQIIfKEEHnl5Y4bk+KCJSJyF5cMi0dazxCUVtbh3T92qC4HgP1wb72TdOSJThJCmABUSCmPW66l/QLIlFJmRkdHd7LM483ngiUichO2e362hU0vLS1EZa36HZvshXslgIh2PE62lHK6A+ppl/omC575jguWiMh9jOsbhdNSbDs2/feHQtXl2A33lTh29W4CsKT1CUKIbCnlHO3j1mPyTvHO7zuwt6oe6bFcsERE7kEIgfsm2mbsvfP7TpRU1Cqtp81w126UmppDW7uxCiHEEu3PLABPCiHyhRD5zi4WACprG4/+VrxvYhoXLBGR2xgYH4ZLhsWj0WI9Orqgio+9E5qvylt9bYL2Zy6AZCfUdVIv/1iE6nozTk2OxLi+XLBERO7lrrNTsXDdXnyxZg9uHGvCwPgwJXV41CKm0so6vP37DgDAzInpEIJX7UTkXhK6B+GGMb0BAI8t2gwp1Sxs8qhwf+a7LWg0W3HBkDgMSlDz25CIyJ6/n5GCsEBf/F50ED9tVdOWwGPCfdOeaixYXQpfo8A92oIBIiJ3FBbki1vOTAEAPPFNASwK2hJ4TLg/+W0BpASuHtkLSZFBqsshImrTtaN7IT48EAX7ajB/1W6XP79HhPvvhQfw09ZydPP3wa3jU1SXQ0RkV4CvEfecYxtleOY717clcPtwt1olnvi2AADwt9NNiGSbASLyEBcOicOAuFDsq67Hm79td+lzu324L1y/F+t2VyEmxB9TT2ObASLyHAaDwP1aW4JXfijCwcMNrntulz1TJzSarUf3Rb09KxVBfnan5RMRuZUxKVE4PTUaNQ1m/Gep69oSuHW4f7h8J3ZV1CI5OhhTMhNUl0NE1Cn3TUyDEMD7y3Zix4EjLnlOtw33mvomvKj9lptxbhp8jG5bKhFRm9JjQ3FpRgLMVomnXNSWwG0Tc+7Pxag40ojhvbrj7P49VJdDRNQld05Ihb+PAQvX7cXqXYec/nxuGe5l1fX4v19sd5bvn5TGNgNE5PHiwgOPTgp5/JsCp7clcMtwf/77bahrsuDs/j0wvFd72skTEbm/m85IRvcgX6zYXoGlBWVOfS63C/ei8sP4ZGUJDAKYcS7bDBCRfoQG+OKW8X0B2FbdO7MtgduF+xztBV8+IhEpMSGqyyEicqhrRiUhPjwQW/cfdmpbArcK9/ydFVi8cT8CfA24PStVdTlERA7n72PE3efY8u3ZJc5rS+A24S6lxOOLbG0GbjzNhB6hAYorIiJyjouGxCM9NhR7q+rxjrZHhaO5Tbjnbi5D3s5DiAj2w/TTTarLISJyGoPh2H6r//2hEFW1TY5/Doc/YieYLVY8qTUHu3V8CkICfBVXRETkXOP6RuHU5EhU15vx8k+Ob0vgFuE+L383CssOIykiCFeP7KW6HCIipxPi2NX7W7/twJ7KOoc+vvJwr2u04LncrQCAu8/pBz8f5SUREbnE4IRwnD84Fo1mK55bstWhj608Sd/8bTv2VzdgUHwYzh8Uq7ocIiKXuvvsfvAxCHy+aje27Ktx2OMqDfeKI4149cciALauaQYD2wwQkXfpHRWMq0YmwSpt63wcRWm4/2fpNtQ0mDEuNRpjUqJUlkJEpMyt4/si2M+I7wvKsLz4oEMeU1m4l1TU4v1lOyEEcN+5aarKICJSLjrEH9PG2aaAP/GtY5qKKQv3p7/bgiaLxCVD49E/LlRVGUREbuHGsSZEdfPD6l2VWLxxf5cfT0m4byitwpdr9sDPaMCdZ7PNABFRN38f3HaWranYnMUFMFusXXo8JeH+xDe2mwbXn9oLCd2DVJRAROR2rjwlCb0jg1BcfgSf5nWtqZjLw/3nreX4tfAAQgJ8cPMZKa5+eiIit+VrNODuc2ytzp/P3YraRnOnH8vl4d581X7zGSnoHuzn6qcnInJrkwbGYnBCGMpqGvDmr9s7/TguDffK2iZs2luN2LAA/GVMb1c+NRGRRzAYxNEZhK/+ZNtLulOP48ii7NlXXQ8AuGNCKgJ8ja58aiIij3FqShROT43G4QYzXlrauaZidsNdCJEthMgSQszozPGWmixW9OsRgkszEjpTKxGR17j33DQIAby3bEenvr/NcBdCZAOAlDIXQKUQIqsjx09Y8MR+MLLNABFRm/rHheLiofFosnRuQZO9K/cRAIq1j4sBZHTwOIQQOUKIPCFEnp+QOLNfTKcKJSLyNndOSIV/Jzvl2vuu8FafR3bwOKSUc6WUmVLKzOSeYRCCV+1ERO2RGBGEX+8d36nvtRfulQAiunD8T3w4HENE1CHRIf6d+j574b4Sx67OTQCWdPA4EREp0Ga4SynnATA13yjVbpxCCLGkreNERKSWj70TpJRzTvC1CW0dJyIitZRvs0dERI7HcCci0iGGOxGRDjHciYh0SDhir752P5kQNQC2uOwJ3VsUgAOqi3ATfC+O4XtxDN+LY/pJKUM68g12Z8s42BYpZaaLn9MtCSHy+F7Y8L04hu/FMXwvjhFC5HX0ezgsQ0SkQwx3IiIdcnW4z3Xx87kzvhfH8L04hu/FMXwvjunwe+HSG6pEROQaHJYhItIhhjsRkQ4x3ImIdMhp4e7IjbU9XVuvVQgRrh3PFkI8qaI+V2nv/3O9vw9Au/5+ZDT/XLi6NlfrQFbkuLo2V9Ne50n3xehIbjol3J2xsbanasdrnQIgQuuND73+ALf3/7n2dZMra3O1dr4X01vsl6Db96MdWZEFoFg7XiyEOG6fZj1pa0+Mjuams67cu7yxto60+Vq1PWabpzmZAOh1w5P2bKZuanGOnrX5Xmi/4POFECYp5RwppZ7fE3s/F3kAPtNC3SSlXOXK4txMh3LTWeHe5Y21daRdr1ULtgod/0Vuz/tg0vHrb8nee5Gs/VchhHhNCNH6fD1p872QUlYCeA3AZ7C9J96sQ7nprHB36MbaHq69rzVbSjnd2cUo1Ob7IITI8qJtGtvzM1GkBVs+AF0O1Wns/VxkA8iVUia3+NxbdSg3nRXu3Fj7GLuvVQiR3bxdoY7vP9h7Hyq0G0XZsI0z63morj1/P5qFw/aXWq/svRcth2Ieh/dcFJ5Ih3LTKeHOjbWPsfdeaF9/UgiRL4TIV1epc7XjZ2KV9rUIHP/PT11p59+P8BbHdbsM3957AWCuECJHOz5Fz+8FcPRfJpkt/4XS2dxk+wEiIh3iIiYiIh1iuBMR6RDDnYhIhxjuREQ6xHAnItIhhjsRkQ4x3ImIdOj/AepI/ngLDXkEAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.figure() \n",
"plt.plot(xx,ff(xx))\n",
"plt.title('$f(x) = \\sin(\\pi x)$')\n",
"plt.savefig('fig1.png')"
]
},
{
"cell_type": "code",
"execution_count": 28,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEICAYAAACktLTqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xd4VFX+BvD3THpCCmmQCkxCSOiEICCCikEFsWHEXpaVsLrqWlHExXUtCHbXteDP3gVBXAGBiB0pCb0ESEIJoSQQUkifmfP7YyYhBsikzMyZufN+noeHJPfOzHcuyZvLued+j5BSgoiItEWnugAiIrI9hjsRkQYx3ImINIjhTkSkQQx3IiINYrgTEWkQw52ISIMY7qSMECJVCJEhhEht8fWMNjw2RAiRbr/qOq4t9bfyWKd9X+RaGO6k0jTL3/rGLwghMqWUC6w9UEpZZtlfb23fjhJC5HfgMW2q/2wc8b7IPTDcSQlLeJVJKRc0hqEQIqQ9zyGlzMKpXxD2MLQ9O7e3/rNxwPsiN+CpugByP5ZgnwOgQAiRKqXcYNmUCWBei33TATSGZhmAVCnlXDvUkwqgAMB1UspHLENF8wEkWGqYD+AiAKEArpVSnil8/1S/EOJtmP9XUmb5uwDA21LKLEe8L3JvDHdyOCllgRACAGY3DkNYJLT4HABCpJQLhBArAVyL089o81v8gmg8g57cSglZUsqCZp9Pa6ylcThESrlBCFFg+ThLCJHd7ONxQoh0yxl2cwnNhlXSpZTThBAZlvpb7t/u90XUHgx3UkV/hiAPbf6JECKk+fi1Zf9rWzymFOaz4g0t9puHtvsSQI4QIgvA263s1/gL4fhZtjfV3yzIG8fOm4ZsOvq+iNqDY+6kSqm1HZqfBQNYafn4TOPaLX9JtFcBzOPrOTAPv7RaT1u1uCja9LED3xe5MZ65k8NZQq+gDftNB5AFYBwsIQggHUDz2SihLZ/LEpSZrTz1ghbDMjOklI8AmGeZiqhvsb2jMgBkWeoZ1qy+Dr0vovZguJMKepjPkltabxmyaDxjLYA55GYDmGEJ3ZbDLQktv2Z5fHsuTh63zE0vg3kGT4HlgqpeCJEJINvy8XSYA3gczBdas1uczbesf1jjRVLLNYZGHXpfRO0huFgHOVpjSLY8O7ac0ae3J9SEEHMsZ93KdaT+Vp7Lad4XuSaeuZPDWMaY9TCflZ425NBsFk17nq+1C6AO1d76z8bZ3he5Jl5QJUcqgPlCamsXLee15fb7xguQNhobt5m21n82zvq+yPVwWIaISIN45k5EpEEOHXMPDw+XPXv2dORLEhG5vJycnGNSyoj2PMah4d6zZ09kZ2c78iWJiFyeEGJ/ex/DYRkiIg1iuBMRaRDDnYhIgxjuREQaxHAnItIghjsRkQZZDXchRLpltZizbc+w7DPdtqUREVFHWZ3nbllW7Izd6SxtUhv30Z9l6TEiZaSUOF5Vj4MnalB0ogYHT1SjpLIOlbUGVNY1oLLWAKNJQicEhAA8dAKBvl4I8fNCiL8XwgK8ER/mj/hQf8R29Yevl4fqt0TUJp29iWkYzEuUAeamUKkwL0LQxNIPOxMA4uPjO/lyRK07WlGLtXtLsb2oHDsOV2DHoQocr6q32fPHh/pjQEww+scEY0BMMIbEhyDAh81Vyfl09ruy5dJgYS13sPS2ngcAaWlp7FJGNnWyzoCfd5Xgt7wSrCkoxd5jVaftE+jjibhQf8R09UNMiB+6BfkiyM8TQb5e6OLrCS+dDhISJgkYjCZU1hpQVl2PspoGFFfWobC0GgdKq1F0ogYHLB8v2XoYAOCpExgcF4JzE8Mxunc4UuO7wkPX+ba/RJ3V2XAvQ4tFjYns7URVPb7ffgQrth/B73nHUW80NW0L8PZAWs9QDIkPQd+oIPSNDkJMiF/LlZA6pMFoQl7xSWwrKse2onJsKizD1qJyZO8/gez9J/DaD3sQ3sUbF/frjvH9u2OEPgxeHpyzQGp0NtzX49TZux6n1oMksimD0YRf9xzD/JxCZO0obgp0IYC0Hl1xUUo3jEwIQ//oIHjaKVC9PHRIiQpCSlQQrk2LAwBU1DZgbUEpfs87hlW5xThQWo3P1h7AZ2sPILyLN65JjcW1aXFIjOxil5qIzsZqP3fLRdN3AEyVUi6wfG2llHKc5ePpADYAONM6kH+SlpYm2TiM2uP4yTp8tvYAPlm7H0cr6gCYA/28xHBMGBCF9JRuiAj0UVylmZQSOw5X4PttR7Bk62EUlJwaIkrr0RW3jOyBCQOieDZP7SaEyJFSprXrMY5crIPhTm2152gl3v1tLxZtLEKdwXyWrg8PwDVDYzEpNQZRwX6KK2ydlBIbC8vw1fpC/G/zIVTVGwEAMSF++MuonrhuWBwCfb0UV0muguFOLm/P0Uq88sMeLN16GI3fmhclR2LKeb1wbkKYTcbOHa2qzoBvNhXh3d/2Np3NB/p4Ysp5vfDX0b0QxJAnKxju5LLyik/i1R/24LsthyAl4O2hw+RhsZgyqhf0EdoYrzaZJFblFmPerwVYt7cUABDs54XMMXrcfm5PTqmks2K4k8s5UVWPV7J245O1B2A0SXh5CFw/LB53XZjg9EMvnbG24DheXLm7KeTDu/hg+iV9kDE0FjpOpaQWGO7kMhqMJnyyZj9eydqD8poG6ARw3bA43D22N2JCtBvqzUkpsTr/OF5YsQsbD5QBAAbEBGPW5X0xrCdnGNMpDHdyCev3lWLGwq3IKz4JwDzz5fGJKUjuHqS4MjWklPh28yE8tywXh8trAQBXDIrGPyf2dZqZQKQWw52cWkVtA+Ysy8Wnaw8AAHqFB2DmhBRclBLpkhdKba263oC3fi7A2z/no85gQrCfF/45sS+uSY3h8XFzDHdyWsu3H8GsxdtwtKIOXh4Cd56fgLsuTGQjrjMoLK3GY4u24tc9xwAAo3uH49mrByAu1F9xZaQKw52cTmVtA5783w4syDkIABgSH4LnJg1En+6BiitzblJKLNxQhKeW7EBZdQP8vT3wxOV9MTktjmfxbojhTk4le18p7v9qEwpLa+DjqcOj45Nx68iebKzVDsdO1uGJb7djyRZzo7JL+3XH7EkD0DXAW3Fl5EgMd3IKBqMJr/2wB6//mAeTBPpFB+HV6wcjMZJn6x0hpcSijUWYtXg7TtYZEBnogxcnD8Lo3hGqSyMH6Ui4s8kF2VRJZR1ueXcdXluVBwngzgsSsOiuUQz2ThBCYFJqLJb9YzTSenRFcWUdbn1vHV5euRtGE7to05kx3Mlm1u0txWWv/Yo/Co4jvIsPPrtjBB65NBnenvw2s4W4UH98kTkC96cnAQBe/WEPbn9/HUptuBgJaQd/6qjTpJR455cC3PDOGhRX1uGcnqFYeu95GJlw2tot1EmeHjr8I703PppyDkIDvPHrnmOY+Nqv2HjghOrSyMkw3KlTahuM+McXm/DM0p0wmiSmjdHjs6nDERnkq7o0TRvdOwLf3XMehsSH4FB5La57ew0WbjiouixyIgx36rDiylpcP28Nvt18CAHeHnjr5qGYMSHFbotl0J9Fh/jhy8yRuGVED9QbTXjgq82Y+30uTByHJzDcqYN2HKrAVa//jk2FZYgJ8cOCO8/Fpf27qy7L7Xh76vDUVf3x1JX94KETeOOnfNz5aQ6q6w2qSyPFGO7Ublk7jiLjrdU4VF6L1PgQfPP3UUiJcs++MM7ilpE98cFfhiHQ1xPLtx9Fxpt/4IilTw25J4Y7tcvn6w4g8+NsVNcbcdXgaHw2dQSbWzmJ0b0jsOiuUegR5o8dhytwzZurm5qzkfthuFObSCnx2g97MGPhVpgkcO9FvfHydYPZG8bJJEZ2wTd3jUJqfAiKymqQ8dZqbOBMGrfEcCerjCaJWYu346WVuyEE8PRV/fHAuCT2OHFSXQO88ekdIzA2ORJl1Q248Z01+DG3WHVZ5GAMd2pVvcGEez/fiI/X7Ie3hw5v3JiKm0f0UF0WWeHn7YG3bxmKjKGxqG0w4Y6PsvF1DqdKuhMu2khnVdtgxF2fbsCq3GIE+njindvSMELPG5NchZeHDs9nDERkoA/e+CkfDy3YjHqjCTecE6+6NHIAnrnTGdXUGzH1o2ysyi1GV38vfJ45gsHugoQQmH5pMmaMT4aUwIyFW/Hh6n2qyyIHYLjTaarqDPjLB+vw655jCO/ijc8zR6B/TLDqsqgTpp2fgFkT+wIAnvh2O975pUBxRWRvDHf6k8raBtz23jqsKShFZKAPvsgc4bZrm2rNlPN64emr+gMAnlm6E//9MU9xRWRPDHdqcrLOgNveW4fs/ScQFeyLL6eNZKtejbl5RA/MvWYghACeX74Lr6/ao7okshOGOwEwj7H/9YP12HDA3E7gq2kj0Ss8QHVZZAeTh8XhxWsHQQjghRW78X+/cohGixjuhDqDEZkfZ2PtXvNQzGdTh3MxZo2blBqLOZMGAgCeXrITn6zZr7gisjWGu5trMJrw90834tc9xxAW4I3Ppg5HjzCesbuDycPi8O8r+wEAHv9mW9Mi5qQNDHc3ZjRJ3PflJmTtPIpgPy98/NfhHGN3M7eO7InHJiQDAKYv2IzvthxSXBHZCsPdTUkpMXPRVizZchhdfDzx0ZRz0Deas2LcUeaYBNyfngSTBO77YhN+2sVWBVrAcHdTL6/cjS/WF8LHU4f3/zIMg+JCVJdECt17USKmjdHDYJK485MN2FRYprok6iSr4S6EyBBCpAshplvZnmn78sgePl6zH6+tyoNOAK/fmIphPUNVl0SKCSHw6PhkXJMai5oGI6Z8sB4FJWwX7MpaDXchRAYASCmzAJQJIdJbbE8HUGDZXiCESLVbpWQT3287jFmLtwEAnr16AMb17aa4InIWQgg8d80AXNgnAqVV9bjl3XU4WsEFP1yVtTP3YQAaJ8EWAGgZ3tkA5ltCXS+l3NDyCYQQmUKIbCFEdklJSacLpo5bU3Ac936xCVICD45LwvVsIEUteHno8N+bUjE4ztwP/rb31qGitkF1WdQB1sK95UDsnzpHSSnLALwNYD6AhDM9gZRynpQyTUqZFhER0eFCqXN2HanE1I+yUW8w4ZYRPXD32ETVJZGT8vf2xHu3D4M+IgC5RyqRafm+IddiLdzLAJx1QNYybJMlpUxo9jk5meLKWvzl/XWorDXg0n7d8a8r+nGhDWpVaIA3PppyDiICfbCmoBQzF22FlFJ1WdQO1sJ9PU6dvesBrGyxvflQzGy08ouA1KhtMGLqRzk4VF6LIfEheOX6wfDQMdjJutiu/nj3tjT4eukwP+cg3vw5X3VJ1A6thruUcgEAfeOFVMuFUwghGkN+nmVMPR3AZCnlPLtWS+1iMkk8+NVmbC4094uZd0sa1zyldhkYG4JXrx8CIYC53+/C0q2HVZdEbWR1JSYp5dwzfG2c5e8yAAx0J/XSyt1YstV8k9J7tw9DRKCP6pLIBV3SrztmjE/Gs0tzcf+XmxAd4ofBvC/C6fEmJo36OucgXv+xcS77EPTpzrYC1HFTR+txwzlxqDOYcMeH2Th4olp1SWQFw12D1u8rxaMLtwAA/nVFP1zQJ1JxReTqhBD495X9MSoxDMdO1uGOD7NRVWdQXRa1guGuMYfKanDnJzloMErcfm5P3Dqyp+qSSCO8PHR446ahTVMkH5q/mTNonBjDXUNqG4z42yc5OHayHqMSw/D4ZSmqSyKNCfbzwju3piHQxxPLth3BGz9xBo2zYrhrhLnL4zZsOViO2K5+eP2GVHh68J+XbC8hogtevWGwZSWnXViVe1R1SXQG/OnXiA9X78PXGw7C10uHebekoWuAt+qSSMPGJnfDg+OSICXwj883IZ9NxpwOw10D/sg/jqeW7AQAzM0YxL7s5BB/vzAR4/t3R2WdAVM/ymYPGifDcHdxRWU1+PtnG2A0SUw7X48rBkWrLonchBACL1w7CH26BaKgpAoPfLkJJhMvsDoLhrsLq20w4m8f56C0qh6je4dj+iXJqksiNxPg44l5tw5FsJ8XsnYW461feIHVWTDcXdjTS3Zga1E54kL98J8bhrBnDCnRIywAL183CADwwvJd+CP/uOKKCGC4u6zFm4rwyZoD8PbQ4Y0bhyLEnxdQSZ2xyd3w9wsTYJLAPZ9vRDEX+VCO4e6C8opPYsbCrQCAWZf3xYDYYMUVEQH3pydhpN58B+vdn2+Ewcge8Cox3F1Mdb0Bd32ag+p6I64YFI2bhnM1JXIOnh46vHrDYEQG+mDd3lK8sGK36pLcGsPdhUgp8fg327D76EnoIwLw7KQBXHSDnEpkoG/T9Z+3fs7Hyh28wUkVhrsL+Sq7EAs3FMHXS4c3bxqKLj5WOzYTOdxwfRimX9IHAPDAV5tQWMoOkiow3F3EjkMVmLV4OwDg6asGsIUvObXMMXqkp3RDZa0B936xEQ0cf3c4hrsLqK434O7PN6DOYMLktFhkDI1VXRJRq8w3OA1EVLAvNh4owytZHH93NIa7C/j3/3agoKQKvSO74Mkr+qsuh6hNQvy98cp1g6ETwBs/5WN13jHVJbkVhruTW7r1ML5YXwhvTx3+c+MQ+HlzDVRyHcP1YbhnbG9ICdz35SYcP1mnuiS3wXB3YkVlNXj0a/OKSjMnpCC5OxuCkeu5Z2wizukZiuLKOjy8YAsX+HAQhruTMhhNuO+LjaioNSA9JRK3juyhuiSiDvH00OGV6wcj2M8Lq3KL8f7v+1SX5BYY7k7q9R/zsH7fCUQG+mBuxiDOZyeXFh3ihznXDAQAPLcsF9uKyhVXpH0MdyeUva8Ur/2wB0IAL183GKFceIM04NL+3XHziHjUG0249/ONqK7nAtv2xHB3MuU1DfjHF5tgkua5wqMSw1WXRGQzj1/WF0nduqDgWBVmL81VXY6mMdydzL++3Y6ishoMjA3Gg+P6qC6HyKZ8vTzw8nWD4eUh8PGa/fhpV7HqkjSL4e5Elm49jEUbze0FXrluMLw9+c9D2tMvOhgPWE5cpi/YghNV9Yor0iamh5MorqjFzEXmNr4zJ6RAH9FFcUVE9pM5Ro9hPbuiuLIOM7/ZyumRdsBwdwJSSjzy9RacqG7A6N7huHkEpz2StnnoBF6aPBgB3h5YuvUIvtlUpLokzWG4O4Ev1hfix10lCPL1xPOc9khuIi7UH09c0Q8AMOsb87Umsh2Gu2L7j1fhqe92AACeuqo/ugf7Kq6IyHGuHRqLi/t2Q2WdAQ99tRkmE4dnbIXhrpDRJPHgV5tRXW/ExIFRuHJwjOqSiBxKCIHZkwYgvIs3/ig4jvd+36u6JM2wGu5CiAwhRLoQYvpZtqda9smwfXnaNu+XAmTvN9+F+vRV7PZI7imsi0/T3atzl+9CXnGl4oq0odVwbwxsKWUWgDIhRPoZdpsmpVwAQC+E0NuhRk3KPVKBl1buAgDMzRiIEH/ehUru66KUbrguLQ71BhMemr8FRg7PdJq1M/dhAAosHxcASG2+UQiRCSBHCKGXUs6VUha0fAIhRKYQIlsIkV1SUmKTol1dg9GEh+ZvRoNR4sbh8bigT6TqkoiUmzkxBVHBvthUWIZ3fzstSqidrIV7SIvPw1p8nmD5UyqEeFsI0XJ/SCnnSSnTpJRpERERnShVO+b9UoBtRRWICfHDYxNSVJdD5BSCfL3w7KQBAIAXVuxGXvFJxRW5NmvhXgYg1Mo++VLKMgA5ADJtUpWG7T5aiVez9gAAnrtmABe5Jmrmwj6RuHZoLOoNJkxfsJnDM51gLdzX49TZux7AyjNsbxQC8y8DOguD0YSHF2xBvdGE64fFYXRv/k+GqKXHJ/ZFtyAfbDhQhvc5e6bDWg33ZhdK0y2fZwGAEGJls+0hzbbPs2+5ru3d3/Zic2EZooJ98dhlHI4hOpNgPy/MtgzPPL98F/JLODzTEVanQloulGY1D24p5bgzbJ9rryK1IL/kJF5caV4B/tlJAxDk66W4IiLnNTa5G65JjUWdwYTpCzh7piN4E5MDGE0S0xdsQb3BhGtSY3EhZ8cQWTXLMjyTs/8Eh2c6gOHuAB+s3occy81Ksyb2VV0OkUsI9v/z8MzeY1WKK3ItDHc723esCs8vN68488zVAxDsz+EYorYam9wNk4bEoM5gwoyFW9gauB0Y7nYkpcSjC7egtsGEKwdHY1zfbqpLInI5/5zYF2EB3lhTUIqvsgtVl+MyGO52ND/7INYUlCI0wBtPXN5PdTlELqlrgDdmXW4eznxmyU4UV9Yqrsg1MNztpKSyDs8s3QnAfGEoNIC9Y4g66opB0bigTwQqag148tsdqstxCQx3O3nqux0or2nAmKQIXDk4WnU5RC5NCIGnr+oPf28PLNl6GFk7jqouyekx3O3gx13F+HbzIfh66fDMVf25shKRDcR29cdDF5sX1v7n4m2orG1QXJFzY7jbWHW9AY8v2gYAeGBcEuJC/RVXRKQdt53bE4PiQnC4vBbPL9+luhynxnC3sZdW7EZRWQ36RQdhyqheqssh0hQPncBzkwbAUyfw8Zr9yNl/QnVJTovhbkNbD5bjvd/3QieA5yYNhKcHDy+RraVEBWHa+XpICcxYaL7zm07H9LERg9GERxdugUkCU0b1woDYYNUlEWnWPWN7o1d4AHYfPYm3fs5XXY5TYrjbyPu/78P2Q+YFOO4fl6S6HCJN8/XywLNXm1sTvL4qj50jz4DhbgOFpdV4ydLx8emr+yOAC3AQ2d3IhDBMTotFvdGEWYu3sTVBCwz3TpJSYtbibahpMOLyQdHs+EjkQI+OT0GIvxd+zzuObzcfUl2OU2G4d9KKHUfx464SBPp64p8TuQAHkSOFBnhjxvhkAMDTS3aignPfmzDcO6G63oAnv90OAHj4kj6IDPRVXBGR+7l2aByG9uiKkso6vMi5700Y7p3w6g97cKi8Fv1jgnDT8B6qyyFySzqduTWBh2Xu+5aDXMoZYLh32O6jlXj3170QAnj6qgHw0LHFAJEqKVFBmDKqJ0wSmLloG5flA8O9Q6SUePybbTCYJG48Jx6D40JUl0Tk9u5LT0JUsC+2FpXj07X7VZejHMO9AxZuKMK6vaUIC/DG9EuSVZdDRAACfDzxhKXv+/Pf73L7vu8M93Yqr27As5Y+7Y9NSOGyeURO5JJ+3XFhnwhU1hnwzJKdqstRiuHeTs+vyMXxqnqc0zMUk1JjVJdDRM0IIfDvK/vDx1OHxZsO4fe8Y6pLUobh3g6bC8vw6doD8NQJPMU+7UROKS7UH/de1BsA8M9vtqHOYFRckRoM9zYymswXUaUE/npeL/TpHqi6JCI6i6mj9UiICEDBsSq8+9te1eUowXBvo0/X7sfWonJEB/s2nRUQkXPy9tThX1eYF6V/fVUeDpfXKK7I8RjubXD8ZB1esNz5NuvyvmwMRuQCRveOwPj+3VFdb8SzS3NVl+NwDPc2eH75LlTUGjAmKQKX9OuuuhwiaqOZl6XA10uH/20+hD/yj6sux6EY7lZsOViGL7ML4akTeOLyvryISuRCYrv6464LEgEA//p2OxqM7rNqE8O9FSaTxKzF2yElMOW8XkiI6KK6JCJqp8wxesSH+mPX0Up8/If73LnKcG/F1xsOYlNhGSICfXDP2ETV5RBRB/h6eWDWRPOdqy+v3I2SyjrFFTkGw/0sKmobMOd780XUGeOTEejLO1GJXNVFKZFNd67O+d49Lq5aDXchRIYQIl0IMd3KfnNsV5Z6r2XtwbGTdRjaoyuuHsI7UYlcmRACsy7vB28PHRbkHMSGAydUl2R3rYa7ECIDAKSUWQDKhBDpZ9kvHYDe9uWpkVdciQ9W74MQwJNX9ONFVCIN6BUegDtG9wIAPLF4u+bbAls7cx8GoMDycQGA1JY7CCH0zfY5jRAiUwiRLYTILikp6XChjiKlxL++3QGDSeKGc+LRPyZYdUlEZCN3j01sagv85fpC1eXYlbVwb9moPOwM++illGcNdynlPCllmpQyLSIiot0FOtry7UfwW94xBPt54aGL+6guh4hsyN/bEzMvM691PHd5Lk5U1SuuyH6shXsZgNCzbRRCpFuGbDShpt6Ip74ztwl98OIkhAZ4K66IiGztsgFRGKkPQ1l1A15cqd01V62F+3qcOnvXA1jZYnup5WJrBgC9EOK0YRtX8vYv+Sgqq0Fy90DceE686nKIyA6EEHjyyn7w0Al8tvYAdh6uUF2SXbQa7lLKBTCHdrrl8ywAEEKstHy+wfK1UJw+hONSCkur8eZP+QDMF1E9PThLlEirkroF4ubh8TBJ4N//2wEptXdx1WqCSSnnSimzpJTzmn1tXIt95kkpE6SUG+xRpCM8u3Qn6gwmXDEoGsP1Z7q0QERacv+4JIT4e+GPguNYvv2o6nJsjqenAP7IP45l247Az8sDMyZwTVQidxDi743705MANJ7caWtRD7cPd6NJ4t/f7QAA3HVBAqKC/RRXRESOctPwePSO7IIDpdV477d9qsuxKbcP9/nZhdh5uALRwb6YOkYz92ERURt4eujwT0vfmddX7UFxZa3iimzHrcO9srYBL6wwT4V6dEIKfL08FFdERI42JikC6SmRqKo3Ni3KowVuHe7//TEfx07WIzU+BJcPjFJdDhEpMvOyvvDyEJifcxDbispVl2MTbhvuB45X4z3LwrlPXM7+MUTurFd4AG4/tyekBJ7833ZNTI1023B/dulO1BtNmJQag0FxLj1Fn4hs4J6LeiMswBvr953Akq2HVZfTaW4Z7n/kH8f3281TH6dfwqmPRAQE+XrhQUs/qdlLc1Hb4NpTI90u3I0miacsUx/vvCAB3YN9FVdERM7iumFxSIkKQlFZDeb9ctZ+iC7B7cJ9QU4hdjROfRzNqY9EdIqHTuCJy81TI9/8KR9Hyl13aqRbhXtlbQOeX74bAPDI+GT4eXPqIxH92Qh9GMb3746aBqNLL8nnVuFunvpYh9T4EFwxKFp1OUTkpB6bkAJvTx0WbSzCRhddks9twr351MdZnPpIRK2IC/XHHeeZl+R7ZslOl5wa6TbhPnuZZerjkBgM5tRHIrLizgsSEN7FG9n7T2DZtiOqy2k3twj3NQWnuj4+fCmXziMi6wJ9vXCfpWvkc8vfWFVyAAAMFElEQVRyXa5rpObD3dRs6uPfzmfXRyJqu+uHxTV1jfz4j/2qy2kXzYf7oo1F2H6oAlHBvshk10ciagdPDx0em2BeUPu1H/a41ILamg73mnojnrd0eXvo4j6c+khE7XZBnwiM7h2OiloDXlu1R3U5babpcH/3twIcqahFv+ggXD0kRnU5ROSChBB4bEIKhAA+/mM/CkpOqi6pTTQb7iWVdU0LXs+ckAKdjlMfiahjUqKCMHloHAwmieeWucaNTZoN91eydqOq3oiLkiNxbmK46nKIyMU9eHES/L09sGLHUawpOK66HKs0Ge57jlbi83UH4KETXPCaiGwiMsgX08YkADDf2GQyOfeNTZoM99nLcmGS5mlMiZGBqsshIo2YOqYXugX5YGtRORZvLlJdTqs0F+6/5x3DqtxidPHxxP3jklSXQ0Qa4u/tiYcta0DM/X4Xauqd98YmTYW7ySTxzJKdABpvHfZRXBERac2kITHoFx2Ew+W1ePc35+35rqlwX7SxCDsOm29YmjKql+pyiEiDdDqBmZeZb2x686d8FFc6Z893zYQ7b1giIkc5NyEc6SmRqKo34uWVznljk2bCnTcsEZEjPTo+BR46gS/XH8CuI5WqyzmNJsKdNywRkaMlRnbBTcPjYZLAM0t3qi7nNJoId96wREQq3JeehEAfT/yyuwS/7ilRXc6fuHy484YlIlIlNMAbd15ovrFp9tJcp7qxyWq4CyEyhBDpQojpZ9gWYtmeIYSYY58SW8cblohIpSmjeiEq2Bc7Dlc41Y1NrYa7ECIDAKSUWQDKhBDpLXaZDCBUSrnAsn+mXao8C96wRESq+Xp54MGLzSu8vbB8N2obnOPGJmtn7sMANM7SLwCQ2nyjlHKelHKe5VM9gKyWTyCEyBRCZAshsktKbDcmxRuWiMhZXD0kBsndA1FUVoOP/tinuhwA1sO95UrSYWfaSQihB1AqpTztdi3LL4A0KWVaREREB8s83ULesERETsJ8zc98Y9Prq/JQVq1+xSZr4V4GILQNz5MhpZxmg3rapLbBiBdX8IYlInIeY3qH47xE84pN//0xT3U5VsN9PU6dvesBrGy5gxAiQ0o51/JxyzF5u/hw9T4cLq9FShRvWCIi5yCEwKPjzTP2Ply9H4Wl1UrraTXcLRdK9Y2hbbmwCiHESsvf6QDmCCFyhBA59i4WAMqq65t+Kz46Ppk3LBGR0+gfE4yrh8Sg3mhqGl1QxdPaDo1n5S2+Ns7ydxaABDvUdVZv/JSPiloDzk0Iw5jevGGJiJzLgxcnYcmWw/hm0yHcMVqP/jHBSupwqZuYispq8MHqfQCAGeNTIATP2onIucR29cfto3oCAJ5duhNSqrmxyaXC/cUVu1BvMOHyQdEYEKvmtyERkTV/vyARwX5eWJ1/HD/vVtOWwGXCfcehCizaWAQvD4GHLTcMEBE5o2B/L9x9YSIA4LlluTAqaEvgMuE+5/tcSAncNLwH4sP8VZdDRNSqW0b2QEyIH3KPVGLhhoMOf32XCPfVecfw8+4SdPHxxD1jE1WXQ0Rkla+XBx6+xDzK8OIKx7clcPpwN5kknvs+FwDwt/P1CGObASJyEVcMika/6CAcqajFe7/vdehrO324L9l6GFsOliMy0AdTzmObASJyHTqdwGOWtgRv/piP4yfrHPfaDnulDqg3mJrWRb0vPQn+3lan5RMROZVRieE4PykClXUG/GeV49oSOHW4f7Z2Pw6UViMhIgCT02JVl0NE1CGPjk+GEMAna/Zj37Eqh7ym04Z7ZW0DXrP8lpt+aTI8PZy2VCKiVqVEBeGa1FgYTBLPO6gtgdMm5rxfClBaVY+hPbri4r7dVJdDRNQpD4xLgo+nDku2HMbGAyfs/npOGe7FFbX4v1/NV5Yfm5DMNgNE5PKiQ/yaJoXMXpZr97YEThnur/ywBzUNRlzctxuG9mhLO3kiIud35wUJ6OrvhXV7S7Eqt9iur+V04Z5fchJfri+ETgDTL2WbASLSjiBfL9w9tjcA81339mxL4HThPtfyhq8bFofEyEDV5RAR2dTNI+IRE+KH3UdP2rUtgVOFe87+UizffhS+Xjrcl56kuhwiIpvz8fTAQ5eY8+2llfZrS+A04S6lxOyl5jYDd5ynR7cgX8UVERHZx5WDYpASFYTD5bX40LJGha05Tbhn7SxG9v4TCA3wxrTz9arLISKyG53u1Hqr//0xD+XVDbZ/DZs/YwcYjCbMsTQHu2dsIgJ9vRRXRERkX2N6h+PchDBU1Brwxs+2b0vgFOG+IOcg8opPIj7UHzcN76G6HCIiuxPi1Nn7+7/vw6GyGps+v/Jwr6k34uWs3QCAhy7pA29P5SURETnEwNgQTBwYhXqDCS+v3G3T51aepO/9vhdHK+owICYYEwdEqS6HiMihHrq4Dzx1Al9vOIhdRypt9rxKw720qh5v/ZQPwNw1TadjmwEici89wwNw4/B4mKT5Ph9bURru/1m1B5V1BoxJisCoxHCVpRARKXPP2N4I8PbAD7nFWFtw3CbPqSzcC0ur8cma/RACePTSZFVlEBEpFxHog6ljzFPAn/veNk3FlIX7Cyt2ocEocfXgGPSNDlJVBhGRU7hjtB7hXbyx8UAZlm8/2unnUxLu24rKsXjTIXh76PDAxWwzQETUxccT915kbio2d3kuDEZTp55PSbg/t8x80eC2c3sgtqu/ihKIiJzODefEo2eYPwpKqvBVdueaijk83H/ZXYLf8o4h0NcTd12Q6OiXJyJyWl4eOjx0ibnV+StZu1Fdb+jwczk83BvP2u+6IBFdA7wd/fJERE5tQv8oDIwNRnFlHd77bW+Hn8eh4V5W3YAdhysQFeyLv4zq6ciXJiJyCTqdaJpB+NbP5rWkO/Q8tizKmiMVtQCA+8clwdfLw5EvTUTkMs5NDMf5SRE4WWfA66s61lTMargLITKEEOlCiOkd2d5cg9GEPt0CcU1qbEdqJSJyG49cmgwhgI/X7OvQ41sNdyFEBgBIKbMAlAkh0tuz/YwFj+8DD7YZICJqVd/oIFw1OAYNxo7d0GTtzH0YgALLxwUAUtu5HUKITCFEthAi21tIXNgnskOFEhG5mwfGJcGng51yrT0qpMXnYe3cDinlPCllmpQyLaF7MITgWTsRUVvEhfrjt0fGduix1sK9DEBoJ7b/iSeHY4iI2iUi0KdDj7MW7utx6uxcD2BlO7cTEZECrYa7lHIBAH3jhVLLhVMIIVa2tp2IiNTytLaDlHLuGb42rrXtRESklvJl9oiIyPYY7kREGsRwJyLSIIY7EZEGCVus1dfmFxOiEsAuh72gcwsHcEx1EU6Cx+IUHotTeCxO6SOlDGzPA6zOlrGxXVLKNAe/plMSQmTzWJjxWJzCY3EKj8UpQojs9j6GwzJERBrEcCci0iBHh/s8B7+eM+OxOIXH4hQei1N4LE5p97Fw6AVVIiJyDA7LEBFpEMOdiEiDGO5ERBpkt3C35cLarq619yqECLFszxBCzFFRn6O09d9c68cBaNPPR2rj94Wja3O0dmRFpqNrczTL+zzruhjtyU27hLs9FtZ2VW14r5MBhFp640Or38Bt/Te3fF3vyNocrY3HYlqz9RI0ezzakBXpAAos2wuEEKet06wlra2J0d7ctNeZe6cX1taQVt+rZY3ZxmlOegBaXfCkLYup65vto2WtHgvLL/gcIYReSjlXSqnlY2Lt+yIbwHxLqOullBscWZyTaVdu2ivcO72wtoa06b1agq1Uwz/IbTkOeg2//+asHYsEy59SIcTbQoiW+2tJq8dCSlkG4G0A82E+Ju6sXblpr3C36cLaLq6t7zVDSjnN3sUo1OpxEEKku9EyjW35nsi3BFsOAE0O1VlY+77IAJAlpUxo9rm7aldu2ivcubD2KVbfqxAio3G5Qg1ff7B2HEotF4oyYB5n1vJQXVt+PhqFwPxDrVXWjkXzoZjZcJ+TwjNpV27aJdy5sPYp1o6F5etzhBA5QogcdZXaVxu+JzZYvhaK0//7qSlt/PkIabZds7fhWzsWAOYJITIt2ydr+VgATf8zSWv+P5SO5ibbDxARaRBvYiIi0iCGOxGRBjHciYg0iOFORKRBDHciIg1iuBMRaRDDnYhIg/4fLOsCQJqDCsUAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.figure() \n",
"plt.plot(xx,ff(xx))\n",
"plt.title('$f(x) = ' + sympy.latex(f)+'$')\n",
"plt.savefig('fig2.png')"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"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.6.7"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment