{ "cells": [ { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import math\n", "%matplotlib inline \n", "\n" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAEACAYAAABS29YJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcHXWZ7/HP0+lO0tkhO1kIayARCDsSwWhAGWVAFAS8\ngLjcO8oorhdH0en0zEtwuePMqNc7vlQMyIg6LgiCQsREQAFZEpYkQAAzhEASsnRnTzrJc/+o6uTk\npE/32eqcqt/5vl+vep2qU9VVv4cKT1c/9atfmbsjIiJha6p3A0REJHlK9iIiDUDJXkSkASjZi4g0\nACV7EZEGoGQvItIAEkv2ZjbJzOab2WIze8bMrk3qWCIi0jtLqp+9mY0Dxrn7IjMbAjwOvMvdlyZy\nQBERKSixK3t3X+Xui+L5zcBS4JCkjiciIoXVpGZvZlOAE4FHanE8ERHZX+LJPi7h/Bz4RHyFLyIi\nNdac5M7NrAX4BXCru9+et06D8oiIlMHdrdSfSfIGrQE3A+vc/VM9rPdyGpwVZjbH3efUux1JUXzZ\nFnJ8aY3NDAMmAdOA6fE0LZ42A0uAxfHnEmCxO+sP3E95uTPJK/uZwBXAU2a2MP7u8+7+uwSPmSZT\n6t2AhE2pdwMSNqXeDUjYlHo3IEFT6nnwOKlPZF9Cz0/q3Qn9YeAmYElPSb3aEkv27v4gemhLRALV\nQ1KflvOZe6X+MPADoqS+oT6tTbhm3+Dm1rsBCZtb7wYkbG69G5CwufVuQILmVnNnOeWX/IQ+DdhC\nlNAXE/U2/CE1ulIvVWI1+z4PHHjNXkSypUBSnw4cy76k3n21vpg6JfVyc6eSfULMbJa7L6h3O5Ki\n+LItLfGpV17vesqRabxBKyLSp2pf9HX/Iuuj98sW9r9Kr9uVeiHV/kWoK3sRqZtq5IE+bpTm1tT3\ndmtMU1IvpNB/G5VxRCRzSskDfdwoze39kqmkXoiSfUakpSaaFMWXbWmJr6c80MeN0j6Telpiq1S1\nk71q9iJSF3FSx4zz2P8BpNzeL91dGmv28FGodGUvIonq/UapDQW/jxR0aUwblXFEJJWKGCbggLFf\nwNaFkgfMbDgwG5jq7jdWYX8q42RBKHXDQhRftlUSXw9X6oWGCeiz/GIJpPl6nTt37zSzx4Hjan3s\nYijZi0iPSrhRmuphAiphZocQxTsbWA0scfd5NW9HuxnQCgwrdx9K9gkJ+aoQFF/W5cZX5DAB3QN6\npf5GabXPnbvfZ2bvBdqBoWb2nrxN1vV2TGu3fkRJeng85c73toy121/j5WHALmBjuXGoZi/SIEpI\n6ntr60kn9bTnAWu1YYznIFbzSS7mPxnAcCbQnygZj4in4Xs/1zCeJziC89iY8/0gor+EOvuYNu43\nP4cFzOGI+PuN3uY7QTdoU0c132zLcnzFDRNwcwe8/17q3PsliTyQf+6s3ZqJku5BfXzmz4/gXkYy\nnj2sYwcj6eQ4VgIdOVPn3s8OtnEHZ7KGYzief+FtPB6v3+RtvqeMONQbJwuynCyKofjqr5JhAtIS\nX7F5IC6FDAcOJkrI3Z/58wfxLIdxDE0537USJd0NRIl5Q9587mdH3nKnt/n2asVbCiV7kQZTwksy\nUj9MQJy0uxP0SObwZ+ZwFTAy/u6AyaNSRmUymGvU9VIkUL10aeyuqaeq94u12wCiJD0qnkbmLI8s\nMA0jqkGvB9bFu3p7PL8eeC5nfgOw3uawHujwNt9dk8BymFnJ5Zc87u79qtKYCunKPiFp+TM5KYqv\nkn3X/0ZpD3XtJqIr6VHA6LzPQtNAosS8tsBn/rQe2JCbtGtRs682MxtPdN5mA6uIumP+PoHj6Mpe\nJAsK3CitydgvcfI+CBgTT6Nz5sdwCdOs3Yi/Hx1vuwl4PZ7WxtPrwErgybzv1xLdeGzEl49YXnfM\nXckdyZqJ8nRLPJWds3VlL1KhPnq/5A8TUPaVurVbC1GyHguMiz/H5HyX+zmKKHmviafX8+bzp3Xe\n5l0lB1+hmuUBMyNKlP2BAXmfuVNP3+03rYChnTD4H2Hmt+EvBv3HR0dpibfp6TN/Pn86YJ1Bi8Me\noIvoF0oX0GUwWjdoRRJUxNgv+e8oXdpXUo+vwEcTJe/caWwPy0OJrqhX503dSXx1zufa7n7ZdRUl\n2QFEPWIGxtPeeYM/ObwzZ93AePtCn/nzucv9877PT+p7gB3Azh4+8+cLTV3vhtNnwKvLYMhUWPXF\n6K+eru71PXx2FVjudZ3BdvXGyQDVtLMrSupvuRjmb6LM3i/xVfhY4BBgfM5n7jSOKNF3EtV+X4s/\nu6fVeZ/ry+mv3UOANh7OfQ2eIHrgZxAwOGd+EFFCHlTgu9Ye5nOngTnzA4iS6PZ42hZ/7gC2GZzh\n8Luc9dtz1ufP7+hhyv1+J7DjfJjxG7if/ITutb/BWwl1vcyIkJMhhBFf728+mtcF5y4kL6kzxzYQ\n3cickDcdkvd5MFF55FWiJN79mTutAlb3egVuNoDoin4oMCRvvtA0OG8+f2r9A3S9NfrFtQXYGn9u\ny1nOnbbFU6H53Gl73vx2vPAvqCzeoK0VJXuREhXd+8V2L2XiIyuZNWcLR8wbFv/MxHiakPO5neim\n5StESXxl3uerP/0v1r93MUPJfZT+wLFPhvUw353Mu+chqr1vIkrOhT63xJ+bc5YLTVvTcpWrPFCY\nkr1IAX0OE9DU9SxDXlvJuCc3cPi87Uy9sx8HLR8NTI5/bgJRH/AVwIrm3bw6eisbDu2g84RVbJ25\ngq63/JXdEzcxhH1PbRZ67L4/+Y/T75vfmDN1j4myKee7TXuX3Xck+J+s7pQHClOyz4hQ/pQspJ7x\n9ZzUfRr4dJp27WDwmlcZuayTQx7tYvKD/Zn08AgGvT7RoGXALlYN28HasVvYOKWDbcespWvaGnjD\nGlqOXUvr4K7oyc77YMzsqFfEBvb1Ec+denvkvhPYQr3+5ypCWv59qoxTmPrZS8M4oPdL87YTaNo1\ng6aBR9K8vYvhyzc2jXlyT9P4x1v3THh8ePPIxV2jWL99wiZaD9tAy9HraDpmMQOPvp8hh29g8EHb\n6Gzad8NuIwf2Kd/7ANA1cOxzcHeaE7ZIKRK7sjez84B/A/oB33f3r+atD/rKXoq390q9qWs6Q199\nI027Traugcey7eCJ1rLF+41cunvP+CdbbPQzNmrQ0u1H9Vu8a/qWdf2OXkfr5E62Tupk9REbWDFy\nG6+wr/dKdxfE7mkt7sk9/CJlUR4oLBNlHDPrRzTGxTlEN60eBS5396U52+gkNxgzjAGdk23ks2cP\nbdow27paT9yxffTknVsmD6P/JvMxSxhw8GI/pHWxH9W0hKks7Txq+7pVkztZMXUdy45cz7J+vvem\n6KvA6tBr2qFTHigsK2Wc04AX3H05gJn9BLgQWNrbD4UklLphIb3FZ4aNPvzW40b1W3XJ9t1Dz9qy\nc+xRm7ceOtJajhjgzVtpad7E4CHP7xnX//mthzY/v3qaPfPQ8Vtfe+qUV1l0xHKeJ7pBuraeJZRG\nPn9ZF3JslUgq2U8g+h+22yvA6QkdS+pk1pQ5Zq3rzt3OoIs7do05rWPnpMkbtx421PpPbV67ejZb\nRyzxg1qXbZs4/M9rJh1064LpTU898M51y/5w5sssxr3s16uJSOmSSvYNf1MrpCsLmzXHzlq94a1d\nDLi8c/fo0zfsnDSpdfQjQ/742jH9mlq2MHjYs7tGtL64cWzr4uXHD7nnsaPGPnfn+zr/NO/MFb49\nGhgwe0I6fz0JOb6QY6tEUsl+JVHXuG6TiK7u92Nmc4Hl8WIHsKj7RJnZLNh34rRci+Umhk2+ffXE\ngYuv2rDxpQu37Rk7YadfMJjOT/d7oOleWge9vHvUiOEdYwY8/9zwnd/7752DXpr/Usfy/2Czu9lb\nZ63I359ZyuLTctqWu6WlPdVYNrPDgEP72P7dwBZ3v6ev/z7x/NXx4nLKlNQN2maiG7SziW6k/YUG\nu0Gb5rqhGUa/7ZNszFNnH9zyyvnsGnjqth3jD9m28YiBPmAzg4ct3XPwoBc2j2p++eXR/VY9PMlX\n3DZxzLz5c+bv+8eS5viqQfHVrB1B9bM3s8OB0939tj62awa+6O5zetkm/Tdo3X2XmX0MuIeo6+UP\nchO91MbeLo22ezrDXz7TmrfMbOka8AZrGT+S/puamls2MnDQyt3jW15cN5n/fujY8c/ecemaZ398\n3Epfk9Xyi0id/Z27f66vjeIceZeZXeXut9SiYXqCNgAHPFHasvkkmneczM7Bh1n/zWZjnurXMnKx\nH9G82Kb5cx3H73r2sTd2rLp95grubO3yA8prIrWShTxgZu8D/t3dR/ex3QnA2e7+rRL2fYu7X1Vg\nXfqv7CUZPQ8TsGc6tmc6zdt3M3LZtqbxjw62cQsHHNm0dOvbNj/NKRvWvzRjFfce/wS/BR7CvaOu\nQYhkz1JgQRHbnQ/cXuK+XzezI939hZJbVSIl+4RUUjfs4SUZOcPv7tnC4DWrGL+wi0PvH86kP08a\nNnTJ62evXbv1/OcZ88YVbJy2iLuanXuAP+C+oVox7d/GdNR8k6L4siuB2N4K/KGI7U4Fbihx308C\nJwNK9iHrPanHL8nov/FFjpi3g2m/WMHhvx/K4NcnN+9m48wVrPvAQrZf+CDjRuzgBeDXRC+BWKbx\nXCQk1m6l/XueA/H7dffytopKRbOAz5rZScAZRO8seIzofuQ73f2D8XaDPO//PTO7ANgNnAU8DZwH\nfNndn4032QAcXUHbiqaafQ30UH4p/Oaj0c+8zLuvGsX4hacQ/SM7HPjTsO38+cbf0/KhhZwyYDdv\nInoTzy+BO3BfW/OgRKog7XkgHvrlGXc/Nh7vayfwcXe/yMwMeNHdD4+3vc/dZ+f87GSgv7u/YGaP\nE/V6eBPwB3ffGm9zDnCqu9/Yw7FVs0+rol+SEXVF/SHRm486if78eztRX9rjgIeA+f328JFVX6d1\n1DbeD3yGqG54C3CZnkAVqYlTicb2wt1/Z2Y3Aj+K170RWJSz7X4D7bn7ywBmNhbY5NH9st/k7X84\nlP7y+XIo2ZehwEsy8pL6zR3w/nnATeS9o9TabRzRS5avJfptv5Kom2ob8IDPYTRwFXAr0VuRfgh8\nDvfVtYmwbyHXfEHxZVmVYzsLWGBmF7n7r4C3AF+J110F/MjMznf33wCrzGyIu2+O23EM0Tt4TyL6\nSxwze4e7352z//HUaMwwJfte9PHmo9zyyyPkJXWzq2e5v38BgLWbWTvTgQuIBoSbSpTc7wY+6W3+\nanzAmcBPif6B/RS4HHhMNXiRunkeOAF40cxagQ5374zXbSG6Mn8qXv4j0SCQ3Tdz30b0asnXgIFm\n9i4OHElgBvD95Jq/j2r29HijND+pd79wursMszT3Sr3HfbZbEzATuIgoybcAdxDdSL1/70umzZri\n9dcBY4H/A9xMXNMTCVma8kClzGwE8Fl3/2KR2w8EbnD3TxdYr5p9uYrq/RIl84eJrtQXu1N010Vr\nNwNOAS4DLiWqxf0CuBh40ttyfrNGJ/pK4LNEr7H7GvCrtLwIWkRK4+4dZrbWzEZ5cZ0mLgO+m3S7\nugV5ZV/gRum+F0/vu0LvvmJf0teVeq/Ha7c3EJVcLgP2ALfxM17yxT63h8Y1Ax8E5gALiZL8/Vkr\n1YRc8wXFV8N2hDY2jgEfdvfv9bHdJOAkd/91L9voyr5bCb1fDqipV3zsdhsOvA/4MDAa+AlwCbDQ\n29xtzv6j+sUNPhf4BtF7Ts/H/YlqtEVE0iHuZ99roo+3W8H+7/xIXCau7Pu4UZqb1LvLMFVL6vu1\nIyrTzCRK8O8C5hGd2N97m+/pJYBjiWrxU4H/DdyetSt5kSSEVLOvtmpf2Wcl2S8FRnDgjdJEkvoB\nx2+3YcCHgP8FGNHd81u8zdf0/oM2kqhccxlwI/BtPL4xKyJK9r1o1DLOKe5sqfVBrd0mAJ8gSvTz\niJL9g/vdaC3gOrPrvgafBH4FHBvaU65pqfkmRfFlV8ixVSITyb7Wid7a7TiiXjJ/S/TE6sneFr08\nve8ftlbgq7PhvcCluP8xqXaKiBQrE2WcWrF2mwl8keghim8B3/U2L75MZHYc8GOiUtNHkhpxUiQU\nacwDadGoZZxEWbsdQ/QI9InAPwPv8jbfUfwOrAn4ONEvis8Ct+gGrIikSVO9G1BP1m7jrd2+SzRu\nxYPAVG/z75eY6EcRDXtwOXAG7jfj7vkvVA6N4su2kOMLObZKNOSVvbXbUKIukH9P1P/+mJLKNXt3\nZJOAe4lGsvsC7l3VbKeISLU0XM3e2u0S4JtESfpL3hYNQ1r6jqx7MLNv4v6N6rVQpHE0Ss3ezGYA\nV7j7Z0v4GdXsy2HtNgr4NtEoc+/yNn+k/J3ZScBdwOfxHoZEEBGJmdmniV5a0tnXtklqiJq9tduF\nRMOQrgROrDDRzyJ6/d9He0v0odcNFV+2hRxf2mLz6C//gmPg1ErQV/bWbgcB/w6cCVzqbf5AZTu0\nC4mGR7gU9/mVt1BEGkTdS1XBJntrt3OI3vD0a+AEb/PKHswyuwr4KvA3uD/e1+ahP8Gn+LItS/GZ\nUeKNRcfyUqt7+cm2iJeGF9eoOgsy2Vu7fQj4MnCFt/nvK9+hnUOU6N9CaSdYRCpUSaKuVPzS8CXx\nS8P/ieh5nA7gZTO7Dmgt8KM3u+/31L2u7KspHpXyS0Qv7j7b2/z5yndqhxG9C/ayUhJ96ONzKL5s\nCzm+asZW4KXhd8Wrv1bKrqrRnkoEk+yt3ZqB7wAnA2d6m6+qfKc2iGggsxsI9H8MESmsyJeGF7Wr\naretVEH0s7d2G0z08pAW4BJv802V79QM+E+gC7hawx+IVF/a+9mb2bXse2n4McCfgVfc/bES9vEx\nooERJwFzgX91941F/FzjjWff637abTRwJ/Ac8GFvq9JTrGafIRoC4Szct1VlnyKyn7Qn+3qqdrJP\npJ+9mX3dzJaa2ZNm9kszG57IcdptIvAn4D7g6iom+nOIBjR7d7mJPm19fatN8WVbyPGFHFslknqo\n6l5gurufADwPfL7aB7B26w/8nOiNUdcX80KR4na894bs5XiZQymIiKRMIsne3ee5730n6yPAxAQO\n8zVgDXBD1fYY3ZD9JXBjpTdkQ+3p0E3xZVvI8YUcWyVq0Rvng8Bt1dxhPJjZBURvkCr8ou/S3UD0\n4pFvVnGfIiJ1V3ayN7N5wLgeVn3B3e+Mt7ke2OnuPy6wj7nA8nixA1jU/Vu5u+6Wv8wcXgW+w6+4\nnic5gTZ63b7Y5QvNrvg4fOAcODJ3PPpy92dmnywmngwvK75sL6civm7V3H/uvusdXzX++8TzV8eL\nyylTYr1xzOxq4H8Cs919ew/rS76jbO02CHgY+I63+X9UpaH7GnQ78Cfcv16d3YX70AoovqxLS3xJ\n9MZJS2yVykTXSzM7D/gX4M3uvrbANuUk+5uA/sCVVbshGzVmFtE4OsfSwy8mEUmGmen5lV5kIdkv\nI0rK3W9/esjdr8nbpqQGW7t9kKg75Gne5pur2Ngm4DHgK7j/rGr7FRFJQKr62bv7Ue5+qLufGE/X\n9P1ThVm7nUA0ENl7qproI1cC24H/quZOQ+/rq/iyLeT4Qo6tEqkfG8farQn4EfAJb/Ol1d25DSYa\nHfNiDYcgIiHLxHAJ1m6TvM1XJNCIfySq019e9X2LiCQgVTX7og5c7zExzA4hehHByew/7rSISGql\nqmafEf8MfC+pRB963VDxZVvI8YUcWyVSX7NPhNkM4J3A1Ho3RUSkFhqvjGNmwDzgF7j/v5ofX0Sk\nAirjFO8k4Ajge/VuiIhIrTRisr8KuAX3XUkeJPS6oeLLtpDjCzm2SjRWzd6sBbgMeFO9myIiUkuN\nVbM3eydwPe5n1vS4IiJVopp9ca4kehpXRKShNE6yj96Dex5Qk8HOQq8bKr5sCzm+kGOrROMke3gP\nMB/3dfVuiIhIrTVOzd5sPvAt3H9Zs2OKiFSZxsbp/WCTgSeACbjvqMkxRUQSoBu0vfsfwM9rmehD\nrxsqvmwLOb6QY6tE+Mk+Gh5BvXBEpKGFX8YxO5moB86RekGJiGSdyjiFXQncqkQvIo0s7GRv1gxc\nDtxa+0OHXTdUfNkWcnwhx1aJsJM9vA14Cfdl9W6IiEg9hV2zN7sNeAD37yR6HBGRGlE/+wMPMAx4\nGThCT82KSCh0g/ZAdR0eIfS6oeLLtpDjCzm2SoSc7C8CflLvRoiIpEGYZZzoQapVwMm4v5LIMURE\n6kBlnP1NAvYAK+vdEBGRNAg12Z8GPFrPB6lCrxsqvmwLOb6QY6tEYsnezD5jZnvM7OCkjtGLU4G/\n1OG4IiKplEjN3swmAd8DpgInu/v6HrZJsmY/H/gK7vcksn8RkTpJW83+G8B1Ce27d2ZNwEnAo3U5\nvohIClU92ZvZhcAr7v5UtfddpKnA6/Tw10QthV43VHzZFnJ8IcdWieZyfsjM5gHjelh1PfB5ojFp\n9m7ey37mAsvjxQ5gkbsviNfNAih5GQ4FHi3756u0DMwws7odX/EpvkaOL6TleP5qIsspU1Vr9mb2\nBuA+YGv81USi7o+nufuavG2TqdmbfRt4Efd/rfq+RUTqrNzcWdaVfSHu/gwwNqdRf6XADdoEnYae\nnBUR2U/S/exr28/drD8wHVhY0+P22JSw64aKL9tCji/k2CpR1Sv7fO5+eJL778HxwAu4b6nxcUVE\nUi2ssXHMrgFOwv3DVd2viEhKpK2ffb3oyVkRkR6EluyjMXFSIPS6oeLLtpDjCzm2SoST7M2GAlOA\nZ+rcEhGR1AmnZm/2ZuBG3M+s2j5FRFJGNfsUlXBERNImpGSfqpuzodcNFV+2hRxfyLFVIqRkryt7\nEZECwqjZm40GlgEH476nKvsUEUmhRq/Zn0r0GkIlehGRHoSS7FNXwgm9bqj4si3k+EKOrRKhJPvo\nyl5ERHqU/Zq9mQFrgBm4r6x4fyIiKdbINftDgS4lehGRwkJI9qnqX98t9Lqh4su2kOMLObZKhJDs\nU3dzVkQkbUKo2S8AbsD93or3JSKScuXmzmwne7N+wAZgCrV9z62ISF006g3aY4DVaUz0odcNFV+2\nhRxfyLFVIuvJ/mhgSb0bISKSdlkv43wUOAH3j1SnVSIi6daoZZzxwKp6N0JEJO2ynuzHAa/VuxE9\nCb1uqPiyLeT4Qo6tEllP9rqyFxEpQtZr9o8Cf4976p6gFRFJgmr2IiJSUHaTvVkTMAZYXe+m9CT0\nuqHiy7aQ4ws5tkpkN9nDSGAT7jvq3RARkbRLpGZvZh8HrgF2A3e5++d62Kaymr3ZccBtuL+h7H2I\niGRMubmzOYGGvAW4ADje3bssehl4Esaher2ISFGSKON8FLjR3bsA3P31BI4B0c3ZVPaxh/Drhoov\n20KOL+TYKpFEsj8KONvMHjazBWZ2SgLHAF3Zi4gUrawyjpnNI0q2+a6P93mQu59hZqcCPwMOL7Cf\nucDyeLEDWOTuC+J1swAKLc+FUzfCmmv37avX7Wu93P1dWtqj+BRfo8Tn7gvS1J5Kl+P5q4ksp0xV\nv0FrZr8FvuLuf4yXXwBOd/d1edtVeoP2NuBO3H9cSXtFRLIkTQ9V3Q68FcDMjgb65yf6Kkn1A1Wh\n1w0VX7aFHF/IsVWi6r1xgJuAm8zsaWAncFUCx4AUD4ImIpI22R0bx6yT6HWEG6rXKhGRdEtTGSd5\nZoOAAUQ3dUVEpA/ZTPbd3S7r9WdJEUKvGyq+bAs5vpBjq0SWk73q9SIiRcpmzd7s3cCVuF9U3VaJ\niKRbY9XsUz5UgohI2mQ12ad+qITQ64aKL9tCji/k2CqR1WSvK3sRkRJktWb/G+C7uN9Z3VaJiKRb\nI9bsU13GERFJk6wm+9R3vQy9bqj4si3k+EKOrRLZS/Zm/YDRwJp6N0VEJCuyV7M3GwMsxj2p1x2K\niKRWI9XsU9/tUkQkbbKY7DPR7TL0uqHiy7aQ4ws5tkpkMdnryl5EpERZrNn/A3AQ7p+rfqtERNJN\nNXsRESkoi8leNfsUUHzZFnJ8IcdWiSwme13Zi4iUKIs1++eBv8X9ueq3SkQk3VSzFxGRgrKV7M0G\nAy3Axno3pS+h1w0VX7aFHF/IsVUiW8m+ewC0FL9oXEQkjbJVszebCXwd9zOTaZWISLo1Ss0+E90u\nRUTSJmvJPjM3Z0OvGyq+bAs5vpBjq0TWkr2u7EVEypC1mv0PgIdw/34yrRIRSbdU1ezN7DQz+4uZ\nLTSzR83s1CrtWlf2IiJlSKqM8zXgS+5+IvCP8XI1qGafEoov20KOL+TYKpFUsn8NGB7PjwBWVmm/\nurIXESlDIjV7MzsUeBBwol8ob3T3FXnblFZ3il40vh1oxX1XFZsrIpIZ5dbsmys44Dyiskq+64Fr\ngWvd/VdmdglwE3BuuceKjQbWK9GLiJQuqSv7je4+LJ43oMPdh+dt48DNwPL4qw5gkbsviNfPAuhe\nvsLsw++Df3iH+5E9rU/h8id7iyeAZcWX7eVg48ut2aehPVWK5+o4pOVAWzlX9kkl+yeAT7n7H81s\nNvAVdz81b5tSyzjnAZ/C/e3VbW0yzGxW94kLkeLLtpDjCzk2KL+Mk1SyPwX4v8AAYBtwjbsvzNum\n1GT/AeDNuF9dxaaKiGRKzWv2vXH3x4DTq7zbzHS7FBFJmywNl5Cpbpeh9/VVfNkWcnwhx1aJLCV7\nXdmLiJQpO2PjmD0AXI/7/cm1SkQk3VI1Nk5CdGUvIlKmLCV71exTRPFlW8jxhRxbJbKR7M2GAAZs\nrndTRESyKBs1e7Mjgd8RPz0rItKoQq/Zj0f1ehGRsmUl2Wfu5mzodUPFl20hxxdybJXISrLP1M1Z\nEZG0yUrN/gZgC+5fTrZVIiLp1gg1e13Zi4iUKSvJXjX7lFF82RZyfCHHVomsJHtd2YuIVCArNfsf\nAtfh/nqyrRIRSbdUvbykqAOX2WARkUYW+g3azAm9bqj4si3k+EKOrRJK9iIiDUBlHBGRDFEZR0RE\nClKyT0jodUPFl20hxxdybJVQshcRaQCq2YuIZIhq9iIiUpCSfUJCrxsqvmwLOb6QY6uEkr2ISANQ\nzV5EJENji+LJAAAESUlEQVRUsxcRkYLKTvZmdomZLTaz3WZ2Ut66z5vZMjN71szeVnkzsyf0uqHi\ny7aQ4ws5tkpUcmX/NHARcH/ul2Y2DbgUmAacB3zHzBrxL4gZ9W5AwhRftoUcX8ixla3sJOzuz7r7\n8z2suhC4zd273H058AJwWrnHybAR9W5AwhRftoUcX8ixlS2JK+5DgFdyll8BJiRwHBERKVJzbyvN\nbB7R+1/zfcHd7yzhOPXp8lNfU+rdgIRNqXcDEjal3g1I2JR6NyBBU+rdgDTqNdm7+7ll7HMlMCln\neWL83QHMLOhfAmb2/nq3IUmKL9tCji/k2MrVa7IvQW6fzzuAH5vZN4jKN0cBf8n/AfWxFxGpnUq6\nXl5kZiuAM4C7zOy3AO6+BPgZsAT4LXCN1+vJLRERAer4BK2IiNRO4v3fzey8+OGqZWb2uQLbfDNe\n/6SZnZh0m6qpr/jMbJaZdZrZwnj6Yj3aWQ4zu8nMVpvZ071sk+Vz12t8GT93k8xsfvzg4zNmdm2B\n7TJ5/oqJL+Pnb6CZPWJmi8xsiZndWGC74s+fuyc2Af2I+tlPAVqARcCxedu8A7g7nj8deDjJNtUh\nvlnAHfVua5nxnQWcCDxdYH1mz12R8WX53I0DZsTzQ4DnAvt/r5j4Mnv+4vYPij+bgYeBN1Vy/pK+\nsj8NeMHdl7t7F/ATooeucl0A3Azg7o8AI8xsbMLtqpZi4oP9b2Bnhrs/AGzoZZMsn7ti4oPsnrtV\n7r4ont8MLCV6BiZXZs9fkfFBRs8fgLtvjWf7E11Yrs/bpKTzl3SynwCsyFnu6QGrnraZmHC7qqWY\n+Bw4M/4z6+54OIlQZPncFSOIc2dmU4j+gnkkb1UQ56+X+DJ9/sysycwWAauB+R51fslV0vmrVtfL\nQoq9+5v/2zcrd42LaecTwCR332pmfwPcDhydbLNqKqvnrhiZP3dmNgT4OfCJ+Ar4gE3yljN1/vqI\nL9Pnz933ADPMbDhwj5nNcvcFeZsVff6SvrLPf8BqEvsPpdDTNgUfwkqhPuNz903df465+2+BFjM7\nuHZNTFSWz12fsn7uzKwF+AVwq7vf3sMmmT5/fcWX9fPXzd07gbuAU/JWlXT+kk72jwFHmdkUM+tP\nNBrmHXnb3AFcBWBmZwAd7r464XZVS5/xmdlYM7N4/jSi7q75tbesyvK561OWz13c7h8AS9z93wps\nltnzV0x8GT9/o8xsRDzfCpwLLMzbrKTzl2gZx913mdnHgHuIbjD8wN2Xmtnfxeu/6+53m9k7zOwF\nYAvwgSTbVE3FxAdcDHzUzHYBW4HL6tbgEpnZbcCbgVHxA3RtRL2OMn/uoO/4yPC5A2YCVwBPmVl3\nkvgCMBmCOH99xke2z9944GaLhodvAn7k7vdVkjv1UJWISANoxJeKiIg0HCV7EZEGoGQvItIAlOxF\nRBqAkr2ISANQshcRaQBK9iIiDUDJXkSkAfx/d/R7R+xHpv8AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "x = np.linspace(0.1, 3)\n", "\n", "f1x = [(i-1)/i for i in x]\n", "f2x = [math.log(i) for i in x]\n", "f3x = x-1\n", "\n", "fig = plt.figure()\n", "ax = fig.gca()\n", "\n", "ax.plot(x,f1x, color='r', label='f1x')\n", "ax.plot(x,f2x, color='g', label='f2x')\n", "ax.plot(x,f3x, color='b', label='f3x')\n", "ax.legend([r'$\\frac{x-1}{x}$',r'$ln(x)$',r'$x-1$'])\n", "ax.grid(True)\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "From the plot above, \n", "$\\frac{x-1}{x} \\leq ln(x) \\leq x-1$ holds only for $ x \\geq 1$ and not $ x > 0$" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.6" } }, "nbformat": 4, "nbformat_minor": 0 }