{ "cells": [ { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "from math import erf, pi, e\n", "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 43, "metadata": {}, "outputs": [], "source": [ "inps = np.linspace(0, 5, 1000)\n", "f = lambda x: (2/pi)**0.5 * e**(-2*x*x) * erf(x*(2**0.5))" ] }, { "cell_type": "code", "execution_count": 44, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 44, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD6CAYAAACvZ4z8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3de3Sc9X3n8fdXd0uWZFmSL5JsyxgDvmAMVgyFQtKEEJMQyGmSBtL00JZTNm1o05NzuiVNT9LS012a7KbZTdhN2IZts0nqUtK0buqUkkDSUDBYDjbGGF/wVfiimy3rfv3uHzNjBiFZj6SZeUbPfF7n6Gjmuek7ifnop9/z+/0ec3dERCS68sIuQERE0ktBLyIScQp6EZGIU9CLiEScgl5EJOIU9CIiERco6M1si5kdMLPDZvbgBPs/aWZ7zWy3mT1rZmvj2xvNrD++fbeZfT3VH0BERC7NphpHb2b5wEHgvUALsBO4x91fTTqmwt0vxF/fCfyOu28xs0bgB+6+PmhBNTU13tjYOM2PISKS23bt2tXu7rUT7SsIcP5m4LC7HwEws63AXcDFoE+EfFwZMONZWI2NjTQ3N8/0dBGRnGRmxyfbF6Trph44mfS+Jb5t/A/5lJm9DnwR+L2kXSvN7CUz+6mZ3TxJgfebWbOZNbe1tQUoSUREggoS9DbBtre12N39EXdfBfwh8MfxzaeB5e5+LfAZ4LtmVjHBuY+6e5O7N9XWTviXh4iIzFCQoG8BliW9bwBOXeL4rcCHANx90N074q93Aa8DV8ysVBERmYkgQb8TWG1mK82sCLgb2JZ8gJmtTnr7AeBQfHtt/GYuZnYZsBo4korCRUQkmClvxrr7iJk9ADwJ5AOPufs+M3sIaHb3bcADZnYrMAycA+6Nn34L8JCZjQCjwCfdvTMdH0RERCY25fDKTGtqanKNuhERmR4z2+XuTRPt08xYEZGICzKOXpKMjI7xr/vOcLyjjxsuW8imFQvDLklE5JIU9NPQOzjCb/zfnbx47M3bDL9+YyNf+OBazCYahSoiEj4FfUDuzu//3W52nTjHlz6ygfetX8KX/+0gf/3cMapKi/j0raunvoiISAgU9AFt33uGp149yx+9/yo+2hSbVvCFD67lfN8Q//PpQ7xnzSLW11eGXKWIyNvpZmwAo2POf3/qAFctKec3b1p5cbuZ8ad3rmdhWREP/fOrZNsIJhERUNAH8tSrZzjS1svvvns1Bflv/Z+ssrSQ33335bx4rJOfHWoPqUIRkckp6AP4zgsnqF8wjy3rl0y4/2PvWMbiimL+z8806VdEso+Cfgqnu/p59nA7H97UQH7exCNrigvy+fjmFfzsUDtH2noyXKGIyKUp6Kewbfcp3OHD171tZea3uOf6ZRTmG9/ecSJDlYmIBKOgn8KP97eydmkFK6rLLnncovISbl2zmG173mBkdCxD1YmITE1BfwnneodoPt7JrWsWBTr+ro31tPcM8dzrHWmuTEQkOAX9JfzkYCtjDu9eszjQ8e+6spby4gK27bnUcv0iIpmloL+EZ15ro2Z+MRsCToQqKcznfeuX8OQrZxgaUfeNiGQHBf0k3J0dRzq4cVU1eZOMtpnIbWsX0z04QvMxLbsvItlBQT+J4x19tHYPsnnl9Fan/MXVNRQV5PGj/a1pqkxEZHoU9JN48WisRX7DZdML+tKiAm5cVc2PXzurJRFEJCso6CfxwtFOFpYVsap2/rTPfc+axRzv6ONIe28aKhMRmR4F/SR2HuvkHY1VM1pn/j1XxYZjPq3uGxHJAgr6CZzrHeJEZx/XLq+a0fl1C+axqraM517XImciEj4F/QT2vtEFwNWzWF/+pstrePFoJ8OaJSsiIQsU9Ga2xcwOmNlhM3twgv2fNLO9ZrbbzJ41s7VJ+z4bP++Amb0vlcWnSyLo19fNPOhvXFVN79AoL7ecT1VZIiIzMmXQm1k+8AhwO7AWuCc5yOO+6+5Xu/tG4IvAl+PnrgXuBtYBW4D/Fb9eVnvljS5WVJdSWVo442vccFk1ZvAfh7UcgoiEK0iLfjNw2N2PuPsQsBW4K/kAd7+Q9LYMSIwrvAvY6u6D7n4UOBy/XlZ7uaVr1o8FXFBaxLq6CvXTi0joggR9PXAy6X1LfNtbmNmnzOx1Yi3635vmufebWbOZNbe1tQWtPS3O9Q7xxvn+wMseXMqNq2r4+fHz9A+NpqAyEZGZCRL0E40vfNtMIHd/xN1XAX8I/PE0z33U3Zvcvam2tjZASemTihuxCb+wqpqh0TF2HT8362uJiMxUkKBvAZYlvW8ALrU841bgQzM8N3T7TsV6odbN4kZswqYVVZhB83GteyMi4QkS9DuB1Wa20syKiN1c3ZZ8gJmtTnr7AeBQ/PU24G4zKzazlcBq4MXZl50+h852s7iieFY3YhMqSgq5cnG5WvQiEqqCqQ5w9xEzewB4EsgHHnP3fWb2ENDs7tuAB8zsVmAYOAfcGz93n5k9DrwKjACfcves7rA+2NrNFYvLU3a9psYq/vGlU4yO+aTPnBURSacpgx7A3bcD28dt+3zS609f4tw/B/58pgVm0tiYc7i1h1+9fkXKrtm0YiHf3nGC185cSEl3kIjIdGlmbJKT5/oYGB7jisXTX8hsMptWxJZRUPeNiIRFQZ/kwJluAFansOumoWoeiyuKaT6moBeRcCjokxxq7QFg9aLUtejNjKYVC9WiF5HQKOiTHDzbTV1lCeUlsx9xk2zTiireON/P6a7+lF5XRCQIBX2Sg2d7Utptk9DUGOunV/eNiIRBQR83NuYcaetJabdNwpqlFRQX5LHnpFayFJHMU9DHnb4wwODIGCtry1J+7cL8PNbVVbBHSxaLSAgU9HHH4s93XVmd+qAHuGbZAva+0cWIHkQiIhmmoI87Gg/6xpr0BP3GZQsYGB7j4NmetFxfRGQyCvq44x29FBfksaSiJC3Xv6ZhAYCeOCUiGaegjzva3seK6lLy0rQezYrqUirnFaqfXkQyTkEfd6yjl8Y09c9DbOLUhoZKdp/sStvPEBGZiIIeGB1zTnT0sTJN/fMJG5ct4ODZbj1xSkQySkEPnO7qZ2h0LG03YhM2NCxgdMzZd0qtehHJHAU9cKy9DyCtXTcA1zTElinerYlTIpJBCnrgaEdiaGVpWn/OoooS6ipL2NOiFr2IZI6CnthkqZLCPBaXp2doZbINDQu0FIKIZJSCntgY+sbqsrQNrUx2dUMlJzr76OofTvvPEhEBBT0AJzr7WL4wvd02CevrY/30uiErIpmS80Hv7rSc66ehKjNBv66uAoB9b1zIyM8TEQkU9Ga2xcwOmNlhM3twgv2fMbNXzexlM/uxma1I2jdqZrvjX9tSWXwqnOsbpm9olIaqeRn5eTXzi1laWcIratGLSIYUTHWAmeUDjwDvBVqAnWa2zd1fTTrsJaDJ3fvM7LeBLwIfi+/rd/eNKa47ZVrOxYZWZiroIdZ9s/cNBb2IZEaQFv1m4LC7H3H3IWArcFfyAe7+jLv3xd/uABpSW2b6tJyLPd4vU103AOvrKjna3kvP4EjGfqaI5K4gQV8PnEx63xLfNpn7gB8mvS8xs2Yz22FmH5roBDO7P35Mc1tbW4CSUifRoq/PaIu+AnfYf1r99CKSfkGCfqIxhz7hgWafAJqALyVtXu7uTcDHga+Y2aq3Xcz9UXdvcvem2traACWlTsu5fipKCqicl9oHgl/K1fGRN3s1cUpEMiBI0LcAy5LeNwCnxh9kZrcCnwPudPfBxHZ3PxX/fgT4CXDtLOpNuUyOuElYVFFCbXmxbsiKSEYECfqdwGozW2lmRcDdwFtGz5jZtcA3iIV8a9L2KjMrjr+uAW4Ckm/ihq7lXF9Gu20S1tdVaIiliGTElEHv7iPAA8CTwH7gcXffZ2YPmdmd8cO+BMwH/n7cMMo1QLOZ7QGeAR4eN1onVG+OoQ8h6OsrOdSqJYtFJP2mHF4J4O7bge3jtn0+6fWtk5z3HHD1bApMp/MXx9BntusGYkE/5rD/zAWuW16V8Z8vIrkjp2fGvjm0MpwWPcA+jacXkTTL8aDP/GSphLrKEqpKC3lF/fQikmY5HvSZnyyVYGasr6/UyBsRSbscD/o+yoszO4Y+2fr6Sg6e7WZwRDdkRSR9cjroT3UNULcg8902CevrKhkedQ6e6QmtBhGJvpwO+jNdAyypTP9TpSaTmCGr7hsRSaecDvrTXQMsDTHoly2cR0VJgVayFJG0ytmgHxoZo71nMNQWvZmxrq6Sfac08kZE0idng/7shQEA6irD66OH2EqW+09fYHh0LNQ6RCS6cjboT3fFgj7MFj3ERt4MjYzxeptuyIpIeuRw0MfG0IfZRw+wri5+Q1YTp0QkTXI26M9kSYt+ZU0ZpUX5vKIbsiKSJjkb9Ke7BphfXEB5STiTpRLy84y1SyvYpyGWIpImORv0YY+hT7a+PjbyZmxswgd3iYjMSs4G/ekL4Y6hT7auroK+oVGOdvSGXYqIRFDuBv35/qwJ+sSSxeqnF5F0yMmgHx4do61nkCUhj6FPuHzRfIoK8jRxSkTSIieDvrV7EPfwh1YmFObnsWZJuVr0IpIWORn0Z+Jj6LPlZizAuvpKXnmjC3fdkBWR1MrJoE/Mig17+YNk6+squTAwcvFhKCIiqRIo6M1si5kdMLPDZvbgBPs/Y2avmtnLZvZjM1uRtO9eMzsU/7o3lcXPVLZMlkq2vr4C0A1ZEUm9KYPezPKBR4DbgbXAPWa2dtxhLwFN7r4BeAL4YvzchcAXgOuBzcAXzKwqdeXPzOmuAUqL8qkoKQi7lIuuWFxOQZ5pbXoRSbkgLfrNwGF3P+LuQ8BW4K7kA9z9GXfvi7/dATTEX78PeMrdO939HPAUsCU1pc/c6a5+llSWYGZhl3JRSWE+qxeXa80bEUm5IEFfD5xMet8S3zaZ+4AfTudcM7vfzJrNrLmtrS1ASbNzpmuAJRXZ022TsL6uQjdkRSTlggT9RM3eCZPIzD4BNAFfms657v6ouze5e1NtbW2AkmantXuQReXFaf8507W+vpKO3iHOXhgMuxQRiZAgQd8CLEt63wCcGn+Qmd0KfA64090Hp3NuJrl7LOizsUWvG7IikgZBgn4nsNrMVppZEXA3sC35ADO7FvgGsZBvTdr1JHCbmVXFb8LeFt8Wmgv9IwyNjGVli37N0grM9LBwEUmtKYeduPuImT1ALKDzgcfcfZ+ZPQQ0u/s2Yl0184G/j9/gPOHud7p7p5n9GbFfFgAPuXtnWj5JQK3dsaGV2diiLy0qYFXtfN2QFZGUCjS+0N23A9vHbft80utbL3HuY8BjMy0w1Vq7Y71K2diih9gN2ReOhvq7UEQiJudmxiYeCp61QV9fyemuAdp7dENWRFIj54L+Yos+C7tu4M1nyGolSxFJldwL+guDlBblM784e2bFJltbp5E3IpJauRf03QNZ220DUDmvkBXVpXqGrIikTA4GfXaOoU+2vq5SI29EJGVyLujbsnRWbLJ19RWc6Oyjq2847FJEJAJyLuhbLwywqDz7W/QA+06r+0ZEZi+ngr5ncITeoVEWVWR5iz5+Q3afum9EJAVyKuhbs3wMfUL1/GLqKku0FIKIpERuBf3FWbHZ3XUDbz5DVkRktnIy6BdnedcNxPrpj7T30js4EnYpIjLH5VbQX+y6yf4W/fr6Ctxh/2n104vI7ORU0Ld1D1JUkEfFvOycFZtsfX1s5M1edd+IyCzlVNAnniyVTc+Kncyi8mJqy4vZ26KgF5HZyamgP3shu5c/SGZmbFy2gN0t58MuRUTmuJwK+liLPvv75xM2LlvAkbZeuvo1Q1ZEZi63gv7CwJwYcZNwTcMCAHXfiMis5EzQDwyPcmFgJOsXNEt2dUPshuzuk+dCrkRE5rKcCfq2+Bj62jnSRw+xJYsvqy1j90m16EVk5nIm6FvnYNBDrJ9+98nzuHvYpYjIHJUzQZ94Bmvt/LkX9O09g5zuGgi7FBGZowIFvZltMbMDZnbYzB6cYP8tZvZzMxsxs4+M2zdqZrvjX9tSVfh0JYK+Zo4FfeKG7J6TGmYpIjMzZdCbWT7wCHA7sBa4x8zWjjvsBPDrwHcnuES/u2+Mf905y3pnrL17CICFZUVhlTAjVy0tpyg/j90KehGZoSBrAWwGDrv7EQAz2wrcBbyaOMDdj8X3jaWhxpRo7xmkcl4hRQVzq7equCCfNXUVCnoRmbEgqVcPnEx63xLfFlSJmTWb2Q4z+9BEB5jZ/fFjmtva2qZx6eDaewapmT+3WvMJGxsq2ftGF6NjuiErItMXJOgnWhhmOomz3N2bgI8DXzGzVW+7mPuj7t7k7k21tbXTuHRwHT1Dc65/PmHj8gX0DY1yuLUn7FJEZA4KEvQtwLKk9w3AqaA/wN1Pxb8fAX4CXDuN+lKmvWeQmjk2tDIhcUNWE6dEZCaCBP1OYLWZrTSzIuBuINDoGTOrMrPi+Osa4CaS+vYzqa1ncM4NrUxorC6joqRAE6dEZEamDHp3HwEeAJ4E9gOPu/s+M3vIzO4EMLN3mFkL8FHgG2a2L376GqDZzPYAzwAPu3vGg35geJTugZE520efl2dcE584JSIyXYGewOHu24Ht47Z9Pun1TmJdOuPPew64epY1zlpHb2xo5Vztowe4dtkCvvbMYXoGR5hfnP0PThGR7DG3xhrOUHv33JwslWxT40LGHHafUKteRKYnN4I+MSt2jt6MBbh2+QLMoPl4Z9iliMgck1tBP0f76AEqSgq5cnE5u45r5I2ITE+OBP3c76MHaGqs4qUT5zVxSkSmJSeCvq17kPLiAkoK88MuZVaaViykZ3CE185cCLsUEZlDciLo5/JkqWSbVlQBqPtGRKYld4J+DvfPJzRUzWNxRTHNxxT0IhJcjgT93F3nJpmZsWlFlVr0IjItORH0HT2DkQh6gE0rFvLG+X5Od/WHXYqIzBGRD/rh0THO9Q1THYGuG4CmeD+9um9EJKjIB31nBJY/SLa2roJ5hfk0H9PEKREJJvJB3xaB5Q+SFebncd2KBbxwVEEvIsFEPugTs2Jry6PRdQPwC5dV89qZ7ot/rYiIXEoOBH20um4AbrisGoAXj3aEXImIzAU5EPTR6roB2NCwgHmF+Tz/uoJeRKYW/aDvHmReYT5lEVrDvaggj6bGKnYcUT+9iEwt+kHfM0hNhPrnE264rJoDZ7vpiP/FIiIymRwI+mjMih0v0U+v0TciMpUcCProzIpNtqGhktKifHYcUT+9iFxaDgT9UCQWNBuvMD+PpsaFuiErIlMKFPRmtsXMDpjZYTN7cIL9t5jZz81sxMw+Mm7fvWZ2KP51b6oKD2JszDnXN0R1WfRa9BAbT3+otYfWCwNhlyIiWWzKoDezfOAR4HZgLXCPma0dd9gJ4NeB7447dyHwBeB6YDPwBTOrmn3ZwXT1DzM65iwsi16LHuDm1TUA/OxQe8iViEg2C9Ki3wwcdvcj7j4EbAXuSj7A3Y+5+8vA2Lhz3wc85e6d7n4OeArYkoK6A+mIzxyNyoJm461dWkF1WRH/fqgt7FJEJIsFCfp64GTS+5b4tiACnWtm95tZs5k1t7WlLrQSSwREtesmL8+4eXUNPzvUzpieIysikwgS9DbBtqCpEuhcd3/U3Zvcvam2tjbgpaeWGGMe1a4bgFuuqKWzd4h9p/QcWRGZWJCgbwGWJb1vAE4FvP5szp21qHfdANy8OvaLUd03IjKZIEG/E1htZivNrAi4G9gW8PpPAreZWVX8Juxt8W0Zkei6qSqNbtDXlhezrq6Cnx5U0IvIxKYMencfAR4gFtD7gcfdfZ+ZPWRmdwKY2TvMrAX4KPANM9sXP7cT+DNivyx2Ag/Ft2VEZ+8QFSUFFBVEe7rALVfU8vPj5+geGA67FBHJQoFW+nL37cD2cds+n/R6J7FumYnOfQx4bBY1zlh7zyDVEZwVO94tq2v53z95nedf7+C2dUvCLkdEskykm7qdvUORvhGbsGlFFfOLC3jmQGvYpYhIFlLQR0BRQR7vvKKWH+1v1TBLEXmbSAd9R28017mZyHvXLqate5A9LefDLkVEskxkg35szHOmRQ/writryc8zfrT/bNiliEiWiWzQXxhIrHMT/ZuxAAtKi9jcuJCnXlXQi8hbRTboL06WypEWPcCtaxdz8GwPxzt6wy5FRLJIZIO+MwdmxY733jWLAdSqF5G3iGzQ58I6N+Mtry7lqiXl/PCVM2GXIiJZJLpBH/GVKydzx4al7Dp+jlPn+8MuRUSyRGSDvrMnvs5NWWHIlWTWHRvqAPiXl0+HXImIZIvIBn1H7xDlJQUUF+SHXUpGNdaUcXV9Jf/8csYWCRWRLBfpoM+lETfJ7tiwlJdbujT6RkSACAd9Z+9gTt2ITfaBDUsB+IG6b0SECAd9R89QzkyWGq+hqpTrli/gn3a/gbvWvhHJdZEN+s4c7roB+PCmBg6e7eHllq6wSxGRkEUy6N1j69zk0mSp8T54TR3FBXk83nxy6oNFJNIiGfQX+kcYGfOc7aMHqCgp5Pb1S9i25xQDw6NhlyMiIYpk0Hf0xmbF5nKLHuBXmpbRPTDCk/s0U1Ykl0Uy6BPr3OTqzdiEGy6rpqFqnrpvRHJcJIM+F1eunEhenvGxpmX8x+EOXm/rCbscEQlJoKA3sy1mdsDMDpvZgxPsLzazv4vvf8HMGuPbG82s38x2x7++ntryJ9bRk3srV07mnuuXU5Sfx7eeOxZ2KSISkimD3szygUeA24G1wD1mtnbcYfcB59z9cuAvgb9I2ve6u2+Mf30yRXVfUmdv7q1cOZma+cXccc1SntjVwoWB4bDLEZEQBGnRbwYOu/sRdx8CtgJ3jTvmLuBv4q+fAN5jZpa6Mqeno3eI+cW5t87NZH7jxpX0Do3yRHNL2KWISAiCBH09kHw3ryW+bcJj3H0E6AKq4/tWmtlLZvZTM7t5lvUGkkvPig3i6oZKNq2o4m+eP8bomGbKiuSaIEE/Uct8fFpMdsxpYLm7Xwt8BviumVW87QeY3W9mzWbW3NbWFqCkS+voye3JUhO57xdXcryjj+17tf6NSK4JEvQtwLKk9w3A+DVwLx5jZgVAJdDp7oPu3gHg7ruA14Erxv8Ad3/U3Zvcvam2tnb6n2KcXF65cjJb1i3h8kXz+drThxlTq14kpwQJ+p3AajNbaWZFwN3AtnHHbAPujb/+CPC0u7uZ1cZv5mJmlwGrgSOpKX1yubxy5WTy8owHfulyDpzt5t/0TFmRnDJl0Mf73B8AngT2A4+7+z4ze8jM7owf9k2g2swOE+uiSQzBvAV42cz2ELtJ+0l370z1hxhXb7yPPrcnS03kjg1Laawu5atPH9KqliI5pCDIQe6+Hdg+btvnk14PAB+d4LzvAd+bZY3TcmFghOFRp0Z99G9TkJ/Hp37pcv7giZf54StneP/VS8MuSUQyIHIzY99c/kBBP5Ffvq6BKxeX8/APX2NoZCzsckQkAyIY9JosdSn5ecYffWANJzr7+Nbzx8IuR0QyIHJBf3H5A/XRT+qdV9Ry8+oavvr0Yc7F/wISkeiKXtD3ap2bID73gTX0DI7wX7bvD7sUEUmzyAW9+uiDuWpJBfffchl/v6uFZw+1h12OiKRR5IK+o2eIsqJ8Sgq1zs1UPv2e1TRWl/JH399L/5CeQiUSVZEL+s7eQRaq2yaQksJ8/usvb+BEZx9/vv3VsMsRkTSJXNB3aLLUtPzCqmp+6+aVfHvHCa2DIxJR0Qv6niFq1D8/LX/wvqvYuGwBf/jEyxzv6A27HBFJscgFvZYonr6igjy+es+1mMFvfauZrn49oEQkSiIV9BfXuVEf/bQtW1jK1z+xiSNtvfz2t3dp1qxIhEQq6LsHRxgaHdMSxTN04+U1PPzhDTz3egd/8MQeRkYV9iJREGhRs7miU7NiZ+0jmxpo7R7gi/96gDGHv/yVayjIj1R7QCTnRCroE7Ni1XUzO7/zrsvJM+PhH77G4PAoX7l7I6VFkfqnIpJTItVUS8yKVdfN7H3ynav4wgfX8tT+s/zKN57ndFd/2CWJyAxFKug7erRyZSr9xk0r+ea9TRxr7+ODX32Wp1/Tk6lE5qJoBX2v+uhT7d1XLeb7v3MjNfOL+c2/buaz/7CXrj4NvxSZSyIV9J29Q5QW5TOvSOvcpNLqxeX80wM38Z9uuYy/23mCd/23Z/j2juMalSMyR0Qu6NVtkx7FBfl89v1r+MHv3syVS8r54398hfd8+ad854XjDAxrQTSRbBapoG/vGdSN2DRbW1fB3/7WDTz6a5tYUFrE577/Cjc+/DR/sm0fr7zRpYeOi2ShSI2Z6+wdYlG5+ufTzcy4bd0S3rt2Mc8f6eA7O07w3RdO8NfPHaOxupR3XbmId15Zy/UrF2pYpkgWCPRfoZltAf4HkA/8lbs/PG5/MfAtYBPQAXzM3Y/F930WuA8YBX7P3Z9MWfXjdPYOsWZpRbouL+OYGTeuquHGVTV09Q3zg72n+NGrZ9m6Mxb6eQZXLC5nQ0MlV9dXsqp2Pitry1hcXkJenoVdvkjOmDLozSwfeAR4L9AC7DSzbe6evID5fcA5d7/czO4G/gL4mJmtBe4G1gF1wI/M7Ap3T3mnrrvT0TukrpuQVJYW8qvXr+BXr1/BwPAoLx7tpPn4OfacPM9Tr57l8eaWi8eWFOZRv2Aei8pLqC0vpra8mJr5xVTOK6SsOJ/ykgLmFxcyv7iAsuJ8igryKMyPfRXHX+frF4VIYEFa9JuBw+5+BMDMtgJ3AclBfxfwJ/HXTwBfMzOLb9/q7oPAUTM7HL/e86kp/029Q6MMjYzpZmwWKCnM55Yrarnlilog9kv4dNcAx9p7OdLey9H2Xk539dPWPcielvO0Xhikf5o3dPOMi+GfZ7G/LsZ/NyDPDLPYd4C8vPg23jxmVmZ5gVT8ujKb3VX0KzN7XLW0gq/ec23Krxsk6OuBk0nvW4DrJzvG3UfMrAuojm/fMe7c+vE/wMzuB+4HWL58edDa32J4ZIw7NizlKnXdZB0zo27BPOoWzOPGy2smPKZ3cISewRG6B2Lfe+OvewdHGB4dY3h0jKFRj30fSbwfY8eqrXcAAAPxSURBVHjEGXPH3XGIv4Yxj/2Cib1+cx/x92Oxl7My2xvPKbltPcuLeGqqkBRZVjUvLdcNEvQT/cIf/69jsmOCnIu7Pwo8CtDU1DSjf3lVZUV87ePXzeRUyQJlxQWUFRewWL+nRVIuyPDKFmBZ0vsG4NRkx5hZAVAJdAY8V0RE0ihI0O8EVpvZSjMrInZzddu4Y7YB98ZffwR42mN/124D7jazYjNbCawGXkxN6SIiEsSUXTfxPvcHgCeJDa98zN33mdlDQLO7bwO+Cfy/+M3WTmK/DIgf9zixG7cjwKfSMeJGREQmZ9k2k7Gpqcmbm5vDLkNEZE4xs13u3jTRvkgtgSAiIm+noBcRiTgFvYhIxCnoRUQiLutuxppZG3B8FpeoAdpTVM5ckWufOdc+L+gz54rZfOYV7l470Y6sC/rZMrPmye48R1WufeZc+7ygz5wr0vWZ1XUjIhJxCnoRkYiLYtA/GnYBIci1z5xrnxf0mXNFWj5z5ProRUTkraLYohcRkSQKehGRiItM0JvZFjM7YGaHzezBsOtJNzN7zMxazeyVsGvJFDNbZmbPmNl+M9tnZp8Ou6Z0M7MSM3vRzPbEP/Ofhl1TJphZvpm9ZGY/CLuWTDGzY2a218x2m1lKV3aMRB99/AHmB0l6gDlwz7gHmEeKmd0C9ADfcvf1YdeTCWa2FFjq7j83s3JgF/ChiP//bECZu/eYWSHwLPBpd98xxalzmpl9BmgCKtz9jrDryQQzOwY0uXvKJ4lFpUV/8QHm7j4EJB5gHlnu/u/E1v7PGe5+2t1/Hn/dDexngmcQR4nH9MTfFsa/5n7r7BLMrAH4APBXYdcSFVEJ+okeYB7pAMh1ZtYIXAu8EG4l6RfvxtgNtAJPuXvUP/NXgP8MjIVdSIY58G9mtsvM7k/lhaMS9IEeQi7RYGbzge8Bv+/uF8KuJ93cfdTdNxJ75vJmM4tsV52Z3QG0uvuusGsJwU3ufh1wO/CpePdsSkQl6PUQ8hwR76f+HvAdd/+HsOvJJHc/D/wE2BJyKel0E3BnvL96K/BuM/t2uCVlhrufin9vBb5PrEs6JaIS9EEeYC5zXPzG5DeB/e7+5bDryQQzqzWzBfHX84BbgdfCrSp93P2z7t7g7o3E/jt+2t0/EXJZaWdmZfEBBphZGXAbkLIRdZEIencfARIPMN8PPO7u+8KtKr3M7G+B54ErzazFzO4Lu6YMuAn4NWKtvN3xr/eHXVSaLQWeMbOXiTVonnL3nBlymEMWA8+a2R7gReBf3P1fU3XxSAyvFBGRyUWiRS8iIpNT0IuIRJyCXkQk4hT0IiIRp6AXEYk4Bb2ISMQp6EVEIu7/AyLeEuzUtl8YAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot(inps, [f(x) for x in inps])" ] }, { "cell_type": "code", "execution_count": 46, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.23695962509005764" ] }, "execution_count": 46, "metadata": {}, "output_type": "execute_result" } ], "source": [ ".42*pi**-0.5" ] }, { "cell_type": "code", "execution_count": 45, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.33637054899797125" ] }, "execution_count": 45, "metadata": {}, "output_type": "execute_result" } ], "source": [ "f(.43)" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [], "source": [ "g = lambda x: x - x**0.8" ] }, { "cell_type": "code", "execution_count": 39, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 39, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD8CAYAAAB0IB+mAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXyV5Z338c+PsG8BQUCWEHZBZY1QlyqKVhQQ62gFrFq1Utux7XRm+hpbbe1Mp4PzPO04bo8+DPIgVEFrXcBicResWlkSlrCGsCRsYd+SkOVczx93AmlM4Jyc7b5zvu/Xi5ee+5xznyuXxx9Xrvu6r6855xARkcavSbIbICIiiaGCLyKSIlTwRURShAq+iEiKUMEXEUkRKvgiIimiabIbcDadO3d2mZmZyW6GiEigrFy58oBz7vzax31Z8M1sEjCpf//+rFixItnNEREJFDPbUddxX07pOOcWOeemp6enJ7spIiKNhi8LvoiIxJ4KvohIilDBFxFJEQkr+GZ2i5n9j5m9ZWbfSNTnioiIJ6qCb2azzazIzNbVOj7ezDaZWZ6ZPQzgnHvTOfcA8B3gjmg+V0REIhftCH8OML7mATNLA54FbgSGAFPNbEiNlzxa9byIiNRWVgyrX4GKspifOqp1+M65pWaWWevwaCDPOZcPYGYLgMlmtgF4HHjHObeqvnOa2XRgOkBGRkY0zRMRCQbnYHc2rJoL6/4Ip45Bqw4w8IaYfkw8brzqARTUeFwIjAF+CFwHpJtZf+fc83W92Tk3E5gJkJWVpXQWEWm8ig/Bmlchex7sWwdNW8FFt8CIu6D35TH/uHgUfKvjmHPOPQU8FdYJatxpKyLSqIRCsO0Tr8hvWASVZdB9BEz4L7jkNmgZvxtO41HwC4FeNR73BHbH4XNERILjaCHkvOwV+iM7oWUHGHUvjLwLul2SkCbEo+AvBwaYWR9gFzAFmBbJCZxzi4BFWVlZD8ShfSIiiVFRBpvf8ebm8z4AHPS5GsY9BhdOhGYtE9qcqAq+mc0HxgKdzawQeMw594KZPQQsAdKA2c653AjPqykdEQmuoo3eSH71Aig+AO17wFU/hRF3QsfMc77dOYdZXbPj0THn/HtdNCsry2m3TBEJhFPHIfcNWDUPCr+EJs1g0I0w8m7ody00STvnKU6cquDV5QXM+2IHc+8bTa/zWjeoKWa20jmXVfu477dHFhHxLeegcDmsehHWvQHlJ6HzIPjGv8PQKdD2K1vS12nXkRJe/Gw78/+6k+OnKsjq3ZFjpeUxb64vC77m8EXE107shzULvNH8gU3QrA1cfKs3mu95KYQ5HZNTcIRZy/J5Z91eAG665ALuv7IPw3t1iEuzfVnwNcIXEd8JVcLWD70LsJsWQ6gCeo6Gm5+Gi74JLdqFdZrKkOPd3L3M+nQbK3ccpl3Lptx/ZR/uuTyTHh1axfVH0By+iMjZHN4O2S9BzktwbBe07gTDpno3R3W5MOzTHC8t59UVhcz5bBsFh0rodV4r7ruiD7dn9aJti9iOvQM1hy8iklTlpbDxbW80v+0TwKD/OBg/AwbeCE2bh32qwsPFvPjZdhZ8WcDxUxVcmtmRR24awvVDupLWJPYrcc7GlwVfUzoikhR713rz8mtegdIj0CEDrnkEhk+D9J4RnSp752FmfbqNP1fNz0+omp8fFqf5+XBoSkdEUlvpUVj7mrdufnc2pDWHwZO8KZs+V0OT8DcVrqgM8e76fcxals+qnUdo17Ip08ZkcM9lmXSP8/x8TZrSERGp5hzs+Is3ml//FlSUQNeLYfx/wtBvQevzIjrdyVMVvLqigBc+3Ubh4RJ6d2rNryYN4fasXrSJ8fx8NPzTEhGReDu+t2o/m9/Doa3Qoj0Mr7oA231E2MspqxUdL2XuZzuY98UOjpaUk9W7I49OSM78fDh8WfA1hy8iMVNZAVve9aZsNi8BVwm9r/C2OhgyGZpHfjdrXtEJZi3L5/VVuygPhbhhSDceuKovo3p3jMMPEDuawxeRxungVq/I57wMJ/ZB265nllN2jnww6Zxj+fbDzFy6lfc3FNGiaRNuz+rJ/Vf2pU/nNnH4ARpOc/gi0viVFXtz8tnzvDl6S4MB3/DugB1wPaQ1i/iUlSHHkty9zFyaT07BETq2bsaPxw3g7st606ltizj8EPGjgi8iwVYdD5g9z1ttc+oYnNfX24J4+DRo161Bpy0pq+S1lQXM+nQbOw4W07tTa359y8XcNrInrZqfeyM0P1LBF5FgKj4Ea//grbTZt9aLBxwy2QsU6X1FxBdgqx08cYq5n3sXYg+dLGN4rw48PP5CvnFRN19eiI2ELwu+LtqKSJ1CIdi+1CvyGxZB5Sm4YHhM4gG3HTjJrGX5vLaykFMVIa4b3JXvXd2XrN4d47I3fTLooq2I+N/peMDfw5EdXjzg0DtiEg+4tvAoz32Sxzvr9tKsSRNuHdmD7369L/27tI1R4xNPF21FJFhOxwPOg60fgAtVxQP+Mup4QOccn289yHOfbGXZlgO0a9mUH4ztxz2XZ9KlXWJjBxNJBV9E/GX/Jm/Tsup4wHbd4ev/BMPvhPP6RHXqUMjx7vp9PPfJVlYXHOH8di342Y0XMm1MBu1aRr6CJ2hU8EUk+U6dgNzXa8QDNoVBN0UUD3g2ZRUh3srZxfOfbGXr/pP07tSa//jmJdw6sgctmwVzxU1DJKzgm1lf4BEg3Tl3W6I+V0R86nQ84FwvC7bsRIPiAc+muKyCBV8WMGtZPruPljL4gvY8PXUEN17cjaZp4W+K1lhEVfDNbDYwEShyzl1c4/h44EkgDZjlnHvcOZcP3G9mr0XzmSIScCcPeNM12fNg/8aqeMBvwoi7odfoBi+nrOlIcRkvfraDOZ9t43BxOaP7nMd/3HoJVw88v9GsuGmIaEf4c4BngLnVB8wsDXgWuB4oBJab2ULn3PooP0tEgipUCVs/8sK+N70DoXIv+zXCeMBz2XO0hBeWbePlL3dSXFbJdYO78P2x/RjVO7LdLxurqAq+c26pmWXWOjwayKsa0WNmC4DJgAq+SKo5HQ/4Mhwr9OIBx3wv4njAc9l5sJjnPsnjtZWFhBzcPKw7D17dj0HdYvMXSWMRjzn8HkBBjceFwBgz6wT8BhhhZj9zzs2o681mNh2YDpCRkRGH5olIXFXHA2bPg/yPOR0PeMNvvAuxEcQDnsvW/Sd49qM83srZTVoTY8qlGUy/qi+9zot8B8xUEI+CX9cEmXPOHQQePNebnXMzzWwPMKl58+ajYt46EYmPveu8C7DV8YDpGTD2595+Nh16xfSjNu49xjMf5vGntXto0bQJ37k8k+lX9aVr+8a7hj4W4lHwC4Ga/3V7ArsjOYFzbhGwKCsr64FYNkxEYiyG8YDhWFt4lKc/3MK76/fRpnkaD17dj/uv7EPngO1amSzxKPjLgQFm1gfYBUwBpkVyAu2lI+JjzsGOz7win/umFw/Y5aIGxwOGY+WOwzzz4RY+2rSfdi2b8qNxA7jvikw6tI7d9FAqiHZZ5nxgLNDZzAqBx5xzL5jZQ8ASvGWZs51zuZGcVyN8ER86vg9WV+1nczDPiwccNsXbz6b7yJgsp6zti/yDPP3hFv6Sd5COrZvx0xsGcddlvWmfAnfFxkO0q3Sm1nN8MbC4oefVCF/EJyorIO89b26+Oh4w43L4+j83OB7wXJxz/CXvIE99sIUvtx+ic9sWPHLTYKaNyfBVIHgQqfdE5KtOxwPOhxN7vXjAy3/Y4HjAcH2+9SBPvLeZL7cfolv7lvxq0hCmjM5Iqe0P4knbI4uIp6wYNiz09rPZ8WmNeMC7vH82IB4wXF9uO8R/vbeJL/IP0bV9C/7+mv7ccWkvWjRVoW8IbY8sIl/lHOzJ8aZsascDDpsK7S+I68ev3HGIJ97bwqd5B+jctgW/nDiEaWM0oo8XXxZ8zeGLxNlX4gFbwpBboo4HDFf2zsM88f4Wlm7eT6c2zXl0wmDuHNM7sFmxQaEpHZFUUV884Mi74OLboFWHuDdhbeFRnnh/Mx9uLKJj62Z87+p+3H1Zb1o39+XYM7A0pSOSqo7uqooHnFcVD5gOo+7xLsBeMDQhTdi09zi/fXcT763fR3orb3nlPZdn0larbhLKl72tKR2RKFWUweY/e3Pzp+MBr4pJPGAkdh4s5on3N/Nmzi7aNm/KT64byH1XZqZEupQfaUpHpDGpKx5wxJ0xiQeMRNHxUp75MI/5X+6kiRnfuSKT71/dT3fGJoimdEQaq1MnvMSo7HlQ8NeqeMAbvUCR/uOijgeMxNGScmYu3crsT7dTVhnijkt78aNrB9AtXZua+YEvC76mdETOwTkoXAHZc2Hd61XxgANjGg8YiZKySl78fDvPfbyVoyXl3DysO/94/UAyO7dJaDvk7HxZ8LWXjkg9EhAPGImKyhCvrijkyQ82s+/YKa4ZdD7/fMMgLuqentB2SHh8WfBFpIbqeMDsubBx8Zl4wElPwcW3xiweMBLOOT7cWMSMdzaSV3SCrN4deXrqSEb3UZSgn6ngi/jV4R2Q85IXEVgdDzh6urduvsvgpDVr3a6j/MfiDXy29SB9O7dh5l2juH5I15QOBw8KFXwRP6k45cUDrpoL+Z94x/pdG5d4wEjtPlLCb9/dxBvZu+jYujn/Nvkipo7OoFlabENOJH5U8EX8YO86b15+zStQcrgqHvBncYkHjNTx0nKe/2Qrs5ZtwwHfu6ofP7imn/akDyBfFnyt0pGUUHoM1r3mbXWwe5UXD3jhRG/Kps/YmMcDRqoy5Jj/5U6eeG8zB0+Wccvw7vzzDYPo2VEB4UGlG69EEsk52Pm5V+Rz3zgTDzjy7rjFAzbEl9sO8djCXDbsOcboPufx6ITBDO0Z/712JDZ045VIMtWOB2zeLu7xgA2x52gJMxZvZOHq3XRPb8mz00Zy0yXddEG2kVDBF4mX0/GA87x9bU7HA/5TVTygf25KKi2vZNayfJ79aCsh5/jxuAE8eHU/bVfcyKjgi8Tawa3eSD7nZS8esE0XuPyhqnjAAclu3d9wzvHe+n38+k/rKThUwo0Xd+PnNw2m13map2+MElbwzawN8H+AMuBj59xLifpskbgrL4H1C73llDs+BWtSFQ94d9zjARuq4FAxv3xrHR9t2s/Arm15+btjuLx/52Q3S+IoqoJvZrOBiUCRc+7iGsfHA08CacAs59zjwK3Aa865RWb2CqCCL8G3u2Y84FHo2MfbgnjYtLjHAzZUeWWI/1mWz1MfbCHNjF9MHMI9l/WmqdbTN3rRjvDnAM8Ac6sPmFka8CxwPVAILDezhUBPYG3Vyyqj/FyR5Ck5DGv+4G11sLc6HnCyN5pPQDxgNJZvP8Qjb6xl874TjL+oG4/dPIQL0lslu1mSIFEVfOfcUjPLrHV4NJDnnMsHMLMFwGS84t8TyAHqHUqY2XRgOkBGRkY0zROJnVAIti/zbo5av7AqHnAYTPhdwuIBo3H4ZBmPv7ORV1YU0KNDK164J4txg7smu1mSYPGYw+8BFNR4XAiMAZ4CnjGzCcCi+t7snJsJzARvHX4c2icSvup4wJzfw+HtSYkHjIZzjrfX7OFXC3M5UlLO967uy4/HDVCGbIqKx3/1un6fdc65k8C9YZ1Ad9pKMlXHA2bPg7z3z8QDXvMoDJ4IzYIxBVJ0vJRfvLmOJbn7GNYznd9/dwyDL2if7GZJEsWj4BcCNTf/6AnsjsPniMTW/s3evHzO/DPxgFf+I4z4dkLjAaPlnOOtnN38alEuxWWVPHzjhXz3yj66KCtxKfjLgQFm1gfYBUwBpkVyAgWgSMKcOgHr3/RW2lTHAw4cDyPvSXg8YCzsO1bKI2+s5f0NRYzM6MD/um0Y/bu0TXazxCeiXZY5HxgLdDazQuAx59wLZvYQsARvWeZs51xuhOfVlI7ET33xgNf/2tvuoG2XZLewQd5es5ufv76WssoQj04YzL1X9CGtiX9XDEniafM0SR0nD8KaBd5WB/s3QLPWcNGt3nLKJMQDxsrx0nIeW5jL66t2MbxXB564Yzh9lCWb0gK1eZpG+BIzoUrI/8gr8hv/5MUD9shKajxgLK3ccYh/eCWHXYdL+NG4Afzw2v4KJJF6aYQvjdORnV40YPbvvXjAVufBsKlJjweMlYrKEE9/mMfTH26hR8dW/PcdwxnV2x9bK0vyaYQvjd/peMB5kP+xd6zftXDDv1fFA7ZIavNipeh4KT+an80X+Ye4dUQP/nXyRbRT+pSEQSN8Cb59uV6RX7PgTDzgiDth+J1JjweMtS/yD/LD+dkcLy3n32+5hNtG9Ux2k8SHAjXCFzmn0mOw7o/ecsrT8YATvAuwPogHjLVQyPH80q38dskmMju1Yd79o7mwm26iksj4suBrSkfqVDMecP2bUF7sxQOOfxyG3uGbeMBYO3Gqgp+8ksN76/cxYegF/OffDaVtC1/+rys+58tvjW68kr9xfB+snu9tdVAdDzj0W95o3kfxgPFQcKiYB+auYPO+4/xy4hDuvSJTcYPSYL4s+CJePOD7XpHf9E5VPOBlvowHjJfl2w/x4LyVlFWGmHPvaK4aeH6ymyQBp4Iv/lIdD7h6Phzf4+t4wHh6fVUh//LHNfTs2JpZ92TR73xtjyDR82XB1xx+iqmOB8ye5+05Xx0PeNNvYeANvowHjBfnHP93aT6Pv7ORy/t14rk7R5HeOnV+fokvLcuU5Nmd4xX5NX84Ew844tswfBq0757s1iVcKOT4zeINvPDpNiYOvYDffWsYLZoGa/M28QctyxR/qCsecPDNZ+IBG9lyynCVV4b46R9W82bObr5zeSa/nDiEJtr4TGJMBV/ir754wJt+C5fc7vt4wHgrrwzxDwty+NPaPfz0hkH8YGw/rcSRuFDBl/g5thtyqvazqY4HHHm3t5/NBcOS3TpfKK8M8eMF2Sxeu5dHJwzmu1/vm+wmSSPmy4Kvi7YBVlnuxQOumgd573nxgJlfD1w8YCJUqNhLgumircRGdTzg6gVwcj+0u8Dby2bEnXCeClltzjn+5Y9reHVFIb+YOIT7rwxOhKL4ny7aSuyVnYTcN7zRfMEXNeIB74Z+4yBNX6/6/O7dzby6opAfjRugYi8Jo/8jJTLOwa6V3qZl616HsuPQaQBc/2/efvMBjQdMpLmfb+eZj/KYOjqDn1yXOjeTSfKp4Et46o0HvAt6jWnU+9nE0rIt+/nVwlyuG9yVX0++SKtxJKFU8KV+9cYDPukV+5banjcSOw8W89DL2Qzs2o4npwynqaIIJcESVvDNrC/wCJDunLstUZ8rDVAdD5jzEhwt8OIBRz/g7WfTdUiyWxdIxWUVTJ/nLUCYeVcWbbS9sSRBWN86M5sNTASKnHMX1zg+HngSSANmOecer+8czrl84H4zey26JktcVJzyRvHZ82DrR96xftfCN37dqOIBk+XXb69n077jvHjvaDI6tU52cyRFhTvMmAM8A8ytPmBmacCzwPVAIbDczBbiFf8Ztd5/n3OuKOrWSuydjgd8BUoOQXovGPuwt59Nh4xkt65R+PO6vcz/soDvj+2nLY4lqcIq+M65pWaWWevwaCCvauSOmS0AJjvnZuD9NtAgZjYdmA6QkaGCExfV8YDZ87wVN9XxgCPugr5joYk27IqVfcdKefj1NVzSI52fXDcw2c2RFBfNRGIPoKDG40JgTH0vNrNOwG+AEWb2s6q/GL7COTfTzPYAk5o3bz4qivZJTc7Bzi+8Ip/7RlU84BC4YYYXD9imU7Jb2Cj966JcSsoqeXLKcJo31UVaSa5oCn5d68nqvW3XOXcQeDCcEyviMIZOFEHOy95+Nge3nIkHHHE39Gjc8YDJ9tHGIhav3ctPbxhEXwWYiA9EU/ALgV41HvcEdkfXHI/20olSzXjAzX+GUIUXD3jlT+CiW1IiHjDZSsoq+eXCdfTv0pYHtEeO+EQ0BX85MMDM+gC7gCnAtFg0SiP8BjqU743kc16uigc8H772A29u/nzNHyfS7L9so+BQCfMf+JqmcsQ3wl2WOR8YC3Q2s0LgMefcC2b2ELAEb2XObOdcbiwapRF+BMpLYMMib6uDFI8H9IsjxWU8/8lWrhvchcv66dqI+Ee4q3Sm1nN8MbA4pi1CI/yw7FntFfm1f4DSqnjAa3+RsvGAfvLcJ1s5caqCn95wYbKbIvI3fHm7n0b49Sg5DGtf8wr93jU14gHvgt5Xpmw8oJ8cKS5j7mc7mDysO4O6tUt2c0T+hi8Lvkb4NYRCsONT7+aoDQuhohS6Da2KB7wNWnVMdgulhpf+upOS8koeHNsv2U0R+QpfFnyN8PlqPGCLdO/iq+IBfetURSX/7y/buXrg+VzYTRvLif/4suCn7Ai/3njAR2DwJMUD+tx76/dx4MQp7lOgifiULwt+yjmwxZuXXz3/TDzglT+BEd9WPGCA/GFFId3TW3Jl/87JbopInXxZ8FNiSqfsJOS+6RV6xQMG3t6jpSzbsp+/v6Y/aU1097L4ky+rSqOd0nEOdq2CVS8qHrCReTNnFyEHt43qmeymiNTLlwW/0Tl50Nt+OHseFK2vigf8pncRNuNr2s+mEViSu5ehPdPp3UnbVoh/+bLgN4opnVDIiwfMrooHrCyDHqMUD9gIHThxipyCI/zDOG1fIf7my4If6CmdIwVnllMeLfDWyV/6XcUDNmIfbizCORg3WFNy4m++LPiBU2c84DXe3PyFExQP2Mh9vKmIbu1bclF3/dYm/qaCH419670iv3qB4gFTlHOOL7cd4qqB52O6FiM+p4IfqdJjkPu6t5xy10po0gwGT1Q8YIrafrCYAyfKuDTzvGQ3ReScVPDD4RwU/NUr8ooHlBqWbz8EwKWZ2tNI/M+XBd83q3ROFHl3v66aVxUP2BYuud27OarHKC2nFLJ3Hia9VTP6dlaEofifLwt+UlfpVFbA1g+80Xx1PGCvrykeUOq0ce9xBl/Qjia6u1YCwJcFPykObauKB3xJ8YASFuccm/ce5/asXud+sYgPpHbBLy+BDW97Wx1UxwP2vx5u+t/evjaKB5SzKDxcwsmySgZ2VdCJBENqFvw9q715+bWvVsUDZioeUCK2ed9xAAZ10/y9BENCC76Z3QJMALoAzzrn3k3Yh5cc8fJfq+MB01rAkMmKB5QG21RV8AdohC8BEXbBN7PZwESgyDl3cY3j44EngTRglnPu8frO4Zx7E3jTzDoCvwXiW/Cdg+2fekVe8YASYwWHiunctjntW2rqT4IhkhH+HOAZYG71ATNLA54FrgcKgeVmthCv+M+o9f77nHNFVf/+aNX74uPYnhrxgNuq4gG/7V2A7T48bh8rqWXXkVK6d1AKmQRH2AXfObfUzDJrHR4N5Dnn8gHMbAEw2Tk3A++3gb9h3r3njwPvOOdW1fU5ZjYdmA6QkRHB9gSV5bB5ibfVwZZ3a8QD/lzxgBIXuw4X64KtBEq0c/g9gIIajwuBMWd5/Q+B64B0M+vvnHu+9gucczOBmQBZWVnunC04sMUr8jnz4WQRtO3mrZkffid06hfJzyISNuccu4+UMnaQdsiU4Ii24Nd1t0m9Rdo59xTw1DlPeq47bctOwvq3vLn5nZ+fiQcccRf0v07xgBJ3h4vLKSmvpIemdCRAoq2MhUDNu056ArujPGfdquMBs+fC2j9WxQP2VzygJMX+46cA6Nq+ZZJbIhK+aAv+cmCAmfUBdgFTgGnRNupvtlYoPuTFA66aq3hA8Y3DxWUAdGytFToSHJEsy5wPjAU6m1kh8Jhz7gUzewhYgrcyZ7ZzLjfaRlVP6Qzt1R5+N+hMPODE/4aL/07xgJJ0R4rLAUhXwZcAiWSVztR6ji8GFsesRdQY4fds8QBZ3/dujup6USw/QiQqR0u8EX6H1s2T3BKR8Pny6ubpi7b9+sGN9d7HJZI0h6tG+JrSkSDx5X4CzrlFzrnp6R06JLspInU6UlxO87QmtGqmhDMJDl8WfDObZGYzjx49muymiNTpSHEZHVo3U46tBIovC/7pEX56erKbIlKnI8XldNB0jgSMLwu+iN8dKSmjQytdsJVg8WXB15SO+N2JUxW0a+nLNQ8i9fJlwdeUjvhdSVklLXXBVgLGlwVfxO9Ky0Mq+BI4KvgiDVBaXknLZvrfR4LFl99YzeGL35WWV2oNvgSOLwu+5vDFz5xzlJRrDl+Cx5cFX8TPyisdIYemdCRw9I0ViVBpRSWARvgSOCr4IhEqLVPBl2DyZcHXRVvxs/KQl+LZLE376Eiw+LLg66Kt+JlzXsG3OiOdRfzLlwVfxM+q6r3SNSVwVPBFInSm4KviS7Co4ItEyOFV/Caq9xIwKvgiEQppSkcCKmEF38wGm9nzZvaamX0/UZ8rEmvVF22bqOJLwIRV8M1stpkVmdm6WsfHm9kmM8szs4fPdg7n3Abn3IPAt4CshjdZJLmqR/giQRPuCH8OML7mATNLA54FbgSGAFPNbIiZXWJmb9f606XqPTcDnwIfxOwnEEm4qmWZGuFLwIQV2eOcW2pmmbUOjwbynHP5AGa2AJjsnJsBTKznPAuBhWb2J+Dlul5jZtOB6QAZGRnhNE8koapX6eiirQRNNBltPYCCGo8LgTH1vdjMxgK3Ai2AxfW9zjk308z2AJOaN28+Kor2icTF6Yu2uvFKAiaagl/Xt73e2U3n3MfAx+Gc2Dm3CFiUlZX1QINaJhJHWpYpQRXNKp1CoFeNxz2B3dE1x6O9dMTPQiHvn5rCl6CJpuAvBwaYWR8zaw5MARbGolHaS0f8zJ3+RVYVX4Il3GWZ84HPgUFmVmhm9zvnKoCHgCXABuBV51xuLBqlEb74mS7aSlCFu0pnaj3HF3OWC7ANpTl88TPtpSNB5cutFTTCFz/TRVsJKl8WfM3hi59pLx0JKl8WfI3wxc9OB6Co4kvA+LLga4QvfnbmxiuRYPFlwRfxN43wJZh8WfA1pSN+pmWZElS+LPia0hE/0146ElS+LPgifnYmACXJDRGJkC8LvqZ0xM9C2llBAsqXBV9TOuJn1TdeaUpHgsaXBV/Ez3TRVoJKBV8kQtpLR4JKBV8kQtpLR4JKBV8kQtpLR4LKlwVfq3TEz7SXjgSVL3mTtZkAAAZrSURBVAu+VumInzntpSMB5cuCL+JnTnvpSECp4ItESMsyJahU8EUipL10JKgSWvDNrI2ZrTSziYn8XJFYOnPRNskNEYlQWAXfzGabWZGZrat1fLyZbTKzPDN7OIxT/QvwakMaKuIXWpYpQdU0zNfNAZ4B5lYfMLM04FngeqAQWG5mC4E0YEat998HDAXWAy2ja7JIslXfeKWKL8ESVsF3zi01s8xah0cDec65fAAzWwBMds7NAL4yZWNm1wBtgCFAiZktds6F6njddGA6QEZGRvg/iUiCaIQvQRXuCL8uPYCCGo8LgTH1vdg59wiAmX0HOFBXsa963UxgJkBWVpar6zUiyeR00VYCKpqCX9e3/ZwF2jk355wnNpsETOrfv38DmiUSX9pLR4IqmlU6hUCvGo97Aruja46I/2lKR4IqmoK/HBhgZn3MrDkwBVgYi0ZpawXxM+2lI0EV7rLM+cDnwCAzKzSz+51zFcBDwBJgA/Cqcy43Fo3S5mniZ9pLR4Iq3FU6U+s5vhhYHNMWeeddBCzKysp6INbnFomW9tKRoPLl1goa4YufaS8dCSpfFnzN4YufaS8dCSpfFnyN8MXPtJeOBJUvC75G+OJnTssyJaB8WfBF/MxpLx0JKF8WfE3piJ/pxisJKl8WfE3piJ9pLx0JKl8WfBE/0146ElQq+CIRqp7S0QBfgsaXBV9z+OJrThdtJZh8WfA1hy9+FtJeOhJQviz4In6m3TIlqFTwRSIU0l46ElAq+CIROnPNVhVfgsWXBV8XbSUQVO8lYHxZ8HXRVkQk9nxZ8EVEJPZU8EVEUoQKvohIilDBFxFJEQkr+GY21syWmdnzZjY2UZ8rIiKesAq+mc02syIzW1fr+Hgz22RmeWb28DlO44ATQEugsGHNFRGRhmoa5uvmAM8Ac6sPmFka8CxwPV4BX25mC4E0YEat998HLHPOfWJmXYH/Au6MrukiIhKJsAq+c26pmWXWOjwayHPO5QOY2QJgsnNuBjDxLKc7DLSo70kzmw5MB8jIyAineSIiEoZwR/h16QEU1HhcCIyp78VmditwA9AB77eFOjnnZgIzq96z38x2VD2VDtS+9bb2sZqPOwMHzvlTNExdbYnVe872uvqeC6dv6jrm5/4K932x6q+6jp+1vzr8Z8r0V33PRdRfNL7v19meT3Z/9a7zqHMurD9AJrCuxuPbgVk1Ht8FPB3u+SL9A8w817Gaj4EViWxLrN5zttfV91w4fRO0/gr3fbHqr3P1Tyr3V33PpXp/ne15v/ZXNKt0CoFeNR73BHZHcb5zWRTGsbpeEw8N+Zxw33O219X3XDh9U9cxP/dXuO+LVX/VdVz9dfbnUr2/zva8L/vLqv42OfcLvTn8t51zF1c9bgpsBsYBu4DlwDTnXG5cWhohM1vhnMtKdjuCQv0VGfVXZNRfkYlXf4W7LHM+8DkwyMwKzex+51wF8BCwBNgAvOqXYl9lZrIbEDDqr8iovyKj/opMXPor7BG+iIgEm7ZWEBFJESr4IiIpQgVfRCRFpGTBN7NbzOx/zOwtM/tGstvjd2bW18xeMLPXkt0WvzKzNmb2YtX3StuGnIO+U5GJVc0KXMGPxUZuzrk3nXMPAN8B7ohjc5MuRv2V75y7P74t9Z8I++5W4LWq79XNCW+sD0TSX6n6naopwv6KSc0KXMHH28htfM0DNTZyuxEYAkw1syFmdomZvV3rT5cab3206n2N2Rxi11+pZg5h9h3ejYfVW41UJrCNfjKH8PtLGtZfUdWsaPbSSQoXg43czMyAx4F3nHOr4tvi5IpFf6WqSPoO787znkAOwRxIRS3C/lqf2Nb5TyT9ZWYbiEHNaixfzLo2cutxltf/ELgOuM3MHoxnw3wqov4ys05m9jwwwsx+Fu/G+Vx9ffc68Hdm9hyJ21IgCOrsL32n6lXf9ysmNStwI/x6WB3H6r2jzDn3FPBU/Jrje5H210EgFf9irEudfeecOwncm+jGBEB9/aXvVN3q66+Y1KzGMsJP9EZuQaf+ajj1XWTUX5GJa381loK/HBhgZn3MrDkwBViY5Db5mfqr4dR3kVF/RSau/RW4gh/QjdySRv3VcOq7yKi/IpOM/tLmaSIiKSJwI3wREWkYFXwRkRShgi8ikiJU8EVEUoQKvohIilDBFxFJESr4IiIpQgVfRCRFqOCLiKSI/w9mL9GUtRlFbAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "inps = np.linspace(0, 100, 10000)\n", "plt.loglog(inps, g(inps))\n", "plt.loglog(inps, inps)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Sample from any symmetric distribution" ] }, { "cell_type": "code", "execution_count": 164, "metadata": {}, "outputs": [], "source": [ "# take any symmetric distribution with std dev s\n", "# sample a, b\n", "\n", "def f(a, b, s, z):\n", " c = 1/(s*2**0.5)\n", " e1 = e**(-(c*a/z)**2)\n", " e2 = e**(-(c*b/z)**2)\n", " return c*(b*erf(c*b/z)*e1 + a*erf(c*a/z)*e2)" ] }, { "cell_type": "code", "execution_count": 165, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 340, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEGCAYAAABo25JHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAXCUlEQVR4nO3dfWxd9X3H8ffn2jht4qyEObUYCUnQ0rVRxwr2IAy1tduuStQpqTSYElQEW0PUibSs3VPYKiqxSX3Y1G5to1HDWFFX4lJWbRmKylbAWvcQlBjY2kCzeikBF1pCktLYKXHM/e6Pew03zr2OY/vc43t/n5dkcc89h3O/X+XKH5/f7zwoIjAzs3QV8i7AzMzy5SAwM0ucg8DMLHEOAjOzxDkIzMwS15p3Aeeqo6MjVq5cWXP96OgoixYtql9B84T7Tk+qvafaN8yu98HBwRcjYmm1dQ0XBCtXrmTfvn011w8MDNDT01O/guYJ952eVHtPtW+YXe+SDtVa56EhM7PEOQjMzBLnIDAzS5yDwMwscQ4CM7PEZRoEktZJOiBpSNL2Gtv8lqQnJe2XdG+W9ZiZ2ZkyO31UUguwA/h1YBjYK2lXRDxZsc1q4Fbg6og4JumNWdVjZmbVZXlEcAUwFBEHI2IM6Ac2TtrmJmBHRBwDiIgXMqzHzMyqyDIILgKerVgeLr9X6U3AmyT9h6Q9ktZlWM+cKRaDw8dP8sorRQ4fP4mf6WBmjSzLK4tV5b3JvzFbgdVAD7AM+Lakt0bET07bkbQV2ArQ2dnJwMBAzQ8dGRmZcv1cOPjiKCdOvkKhUAqFhQtauaQj30ve69H3fJRq35Bu76n2Ddn1nmUQDAPLK5aXAc9V2WZPRJwCfiDpAKVg2Fu5UUT0AX0A3d3dMdUl1llffn74+Em2fPIhxostr77XWhD/detVLF28ILPPPZtUL7tPtW9It/dU+4bses9yaGgvsFrSKkltwCZg16Rt/hHoBZDUQWmo6GCGNc1KsRhEBJdffD4tBbH4da20CLpWLOGChed5mMjMGlJmRwQRMS5pG/Ag0ALcHRH7Jd0O7IuIXeV175X0JPAK8IcRcSSrmmajWAw237mHwUPHuPziJfzn9nfRsaiNoydOccHC87jurkcZPHSMrhVL2HnTWgqFaiNjZmbzT6Z3H42I3cDuSe/dVvE6gI+Vf+a1I6NjDB46xngxeOyZYxQkWloKLF28gMPHT766bvDQMY6MjuU6TGRmdi58ZfE0dbS30bViCa0F0bViCR3tbdNaZ2Y23zXc8wjyUCwGR0bHuHfLlRw9cYqO9jak14Z+JLHzprUcGR07Y52Z2XznIDiLyrmBifH/ar/oCwV5OMjMGpKHhs6icm5gYvzfzKyZOAjOwuP/ZtbsPDQ0hbPNDZiZNQMHQQ3TnRuo9f964tjMGoWDoIZqcwPTmQyuFiC+uMzM5jPPEdQw07kBTy6bWaPxEUENM702YCJAJo4IPLlsZvOdg2AKM7k2wBeXmVmj8dBQDRMPn5nJ3UQnAsQhYGaNwEcEVXjC18xS4iOCKjzha2YpcRBU4auJzSwlHhqqwhO+ZpYSHxFU4SuDzSwlPiKYxBPFZpYaHxFM4oliM0uNg2ASTxSbWWo8NDSJJ4rNLDUOgir82EkzS4mHhiaZza0l8tivmdls+YigQlZnDPlMJDObzzI9IpC0TtIBSUOStldZf6Okw5KeKP9sybKes8nqjCGfiWRm81lmQSCpBdgBrAfWAJslramy6dci4m3ln7uyqmc6sjpjyGcimdl8luXQ0BXAUEQcBJDUD2wEnszwM2clqzOGfCaSmc1nymryUtI1wLqI2FJevh64MiK2VWxzI/BJ4DDwv8BHI+LZKvvaCmwF6Ozs7Orv76/5uSMjI7S3t8+47vFi0NqA4/ez7btRpdo3pNt7qn3D7Hrv7e0djIjuauuyPCKo9tt0cur8M7AzIk5K+hBwD/CuM/6niD6gD6C7uzt6enpqfujAwABTra+l0Sd0Z9p3o0u1b0i391T7hux6z3KyeBhYXrG8DHiucoOIOBIRJ8uLdwJdGdYzJU/omlmqsgyCvcBqSasktQGbgF2VG0i6sGJxA/BUhvVMyRO6ZpaqzIaGImJc0jbgQaAFuDsi9ku6HdgXEbuAj0jaAIwDR4Ebs6rnbDyha2apyvSCsojYDeye9N5tFa9vBW7Nsobp8jMIzCxVvrKYxp8oNjObDd9rCE8Um1naHAR4otjM0uahITxRbGZp8xEBnig2s7Qlf0TgiWIzS13yRwSeKDaz1CUfBJ4oNrPUJT805IliM0td8kcEZmapS/6IwJPFZpa65I8IPFlsZqlLPgg8WWxmqUt+aMiTxWaWuuSPCPK8qrhYDA4fP0lWz402M5uOpI8I8pwo9iS1mc0XSR8R5DlR7ElqM5svkg6CPCeKPUltZvNF0kNDeU4Ue5LazOaLpIMgb4WCWLp4Qd5lmFnikg4CT9iamSU+R+AJWzOzxIPAE7ZmZokPDUXA5zddhgRLFy/whK2ZJSnTIwJJ6yQdkDQkafsU210jKSR1Z1lPpYn5gas//TAf3vk4vrjXzFKVWRBIagF2AOuBNcBmSWuqbLcY+AjwaFa1VOP5ATOzkiyPCK4AhiLiYESMAf3Axirb/RnwGeDlDGs5g+cHzMxKlNUNzyRdA6yLiC3l5euBKyNiW8U2lwEfj4jflDQA/EFE7Kuyr63AVoDOzs6u/v7+mp87MjJCe3v7tOt8+VSR153X+HPm59p3s0i1b0i391T7htn13tvbOxgRVYffs5wsrjbz+mrqSCoAnwNuPNuOIqIP6APo7u6Onp6emtsODAww1foJzXYNwXT7bjap9g3p9p5q35Bd71n+KTwMLK9YXgY8V7G8GHgrMCDpaWAtsKteE8aeIzAzK8kyCPYCqyWtktQGbAJ2TayMiJcioiMiVkbESmAPsKHa0FAWPEdgZlaS2dBQRIxL2gY8CLQAd0fEfkm3A/siYtfUe8iWryEwMyvJ9IKyiNgN7J703m01tu3JspZK1eYHnANmlqrGP11mBjw/YGb2miSDoKO9jcsvPp8WweUXn+/5ATNLWpJBULp0QqDSj28vYWYpSzIIjoyO8dgzx3ilGDzmoSEzS1ySQeBTR83MXpPkbah96qiZ2WuSCwKfOmpmdrrkhoZ86qiZ2emSCwKfOmpmdrrkgsCnjpqZnS65IPCpo2Zmp0suCDw0ZGZ2uuSCwENDZmanSy4IPDRkZna65ILggoXncemyN9AifFWxmRmJBUGxGFx316P89/BLXLr8fO7dcqWvKjaz5CUVBBMXk71SDL4z/BJHT5zKuyQzs9wlFQQ+Y8jM7ExJBUGjnTFULAaHj58k5nuhZtbQkrrpXLUzhpYuXpB3WVVVuzleoeD5DDObe0kdETTS0JBvjmdm9ZJUEDTS0JAfnmNm9eKhoXk6NCSJnTet5cjoGB3tbT7N1cwyk+kRgaR1kg5IGpK0vcr6D0n6jqQnJP27pDVZ1tNIQ0MAhYL8BDUzy1xmQSCpBdgBrAfWAJur/KK/NyJ+OSLeBnwG+GxW9UBjDQ2ZmdVLlkcEVwBDEXEwIsaAfmBj5QYR8dOKxUVApr+aD4+cZPDQUd9nyMysgrI6R13SNcC6iNhSXr4euDIitk3a7mbgY0Ab8K6I+H6VfW0FtgJ0dnZ29ff31/zckZER2tvbq647+OIooyfHAVjU1solSxfNoLP5aaq+m1mqfUO6vafaN8yu997e3sGI6K66MiKq/gBfKf/3llrbTPUDXAvcVbF8PfCFKba/DrjnbPvt6uqKqTzyyCNV33/hpy/HJdsfiBV//ECs2v5A/Piln025n0ZTq+9ml2rfEen2nmrfEbPrHdgXNX6vTjU01CVpBfA7kpZIuqDyZxoBNAwsr1heBjw3xfb9wPunsd8ZuWDheSxcUDpJatGCVn5+0fyeKDYzq5epTh+9A/gmcAkwCFSeuhLl96eyF1gtaRXwQ2ATpb/6XyVpdbw2FPQ+4Ixhobly9MQpToy9AsCJsVc4euLUvD111MysnmoeEUTE5yPiLcDdEXFJRKyq+DlbCBAR48A24EHgKeC+iNgv6XZJG8qbbZO0X9ITlOYJbph9S9V1tLfRVT51tKsBTh01M6uXs15QFhG/O9OdR8RuYPek926reH3LTPd97rXA5FNHfXq+mVlCt5go3bundOro4NNHfeqomVlZMkFQOVm8cEErFyw8L+eKzMzmh2SC4MXRsVevIZiYLDYzs0SCoFgMPrzzcYrla+c8WWxm9pokguDI6BiPHToGQEtBfPG6y30jNzOzsiSCoPLe/t0rlvj6ATOzCkk8j8D39jczqy2JIIDX7u1vZmanS2JoyMzManMQmJklzkFgZpY4B4GZWeIcBGZmiXMQmJklzkFgZpY4B4GZWeIcBGZmiXMQmJklzkFgZpY4B4GZWeIcBGZmiXMQmJklzkFgZpY4B4GZWeIyDQJJ6yQdkDQkaXuV9R+T9KSk/5H0kKQVWdZjZmZnyiwIJLUAO4D1wBpgs6Q1kzZ7HOiOiEuB+4HPZFWPmZlVl+URwRXAUEQcjIgxoB/YWLlBRDwSESfKi3uAZRnWY2ZmVSgistmxdA2wLiK2lJevB66MiG01tv8i8KOI+PMq67YCWwE6Ozu7+vv7a37uyMgI7e3tc9BBY3Hf6Um191T7htn13tvbOxgR3dXWZfnwelV5r2rqSPoA0A28s9r6iOgD+gC6u7ujp6en5ocODAww1fpm5b7Tk2rvqfYN2fWeZRAMA8srlpcBz03eSNJ7gD8F3hkRJzOsx8zMqshyjmAvsFrSKkltwCZgV+UGki4DvgRsiIgXMqzFzMxqyCwIImIc2AY8CDwF3BcR+yXdLmlDebO/ANqBr0t6QtKuGruzKRSLwXgxyGq+x8yaW5ZDQ0TEbmD3pPduq3j9niw/PwXFYrD5zj28vf04d/btYedNaykUqk3PmJlV5yuLG9yR0TEGDx0jCAYPHePI6FjeJZlZg3EQNLiO9ja6VixBiK4VS+hob8u7JDNrMA6CBieJnTet5c0XLqZ/61okDwuZ2blxEDSBQkG0FuQQMLMZcRCYmSXOQWBmljgHgZlZ4hwEZmaJcxCYmSXOQWBmljgHgZlZ4hwEZmaJcxCYmSXOQWBmljgHgZlZ4hwEZmaJcxCYmSXOQWBmljgHgZlZ4hwEZmaJcxCYmSXOQWBmljgHgZlZ4jINAknrJB2QNCRpe5X175D0mKRxSddkWYuZmVWXWRBIagF2AOuBNcBmSWsmbfYMcCNwb1Z1mJnZ1Foz3PcVwFBEHASQ1A9sBJ6c2CAini6vK2ZYh5mZTSHLoaGLgGcrlofL75mZ2TyS5RGBqrwXM9qRtBXYCtDZ2cnAwEDNbUdGRqZc36zcd3pS7T3VviG73rMMgmFgecXyMuC5mewoIvqAPoDu7u7o6empue3AwABTrW9W7js9qfaeat+QXe9ZDg3tBVZLWiWpDdgE7Mrw88zMbAYyC4KIGAe2AQ8CTwH3RcR+SbdL2gAg6VclDQPXAl+StD+reszMrLosh4aIiN3A7knv3Vbxei+lISMzM8uJryw2M0ucg8DMLHEOAjOzxDkIzMwS5yAwM0ucg8AoFoPDx08SMaMLv82swWV6+qjNf8VisPnOPQweOkbXiiXsvGkthUK1u4OYWbPyEUHijoyOMXjoGOPFYPDQMY6MjuVdkpnVmYMgcR3tbXStWEJrQXStWEJHe1veJZlZnXloKHGS2HnTWo6MjtHR3obkYSGz1DgIjEJBLF28IO8yzCwnHhoyM0ucg8DMLHEOAjOzxDkIzMwS5yAwM0ucg8DMLHEOAjOzxDkIzMwS5yAwM0ucg8DMLHEOAjOzxDkIzMwS5yCwuvGT0Mzmp0yDQNI6SQckDUnaXmX9AklfK69/VNLKLOux/Ew8Ce2qTz7Epr49FIsOA7P5IrMgkNQC7ADWA2uAzZLWTNrsg8CxiPhF4HPAp7Oqx/LVjE9C8xGO1VOxGIwXI5PvW5ZHBFcAQxFxMCLGgH5g46RtNgL3lF/fD7xbfjJKU2q2J6H5CMfqaeL79r3nj2fyfVNWf81IugZYFxFbysvXA1dGxLaKbb5b3ma4vPx/5W1enLSvrcBWgM7Ozq7+/v6anzsyMkJ7e/tctzPvNUrf48WgtTB3WZ9X3+PF4HvPHycIhHjzhYvntK/paJR/87mWYt8T37c3vj544Wcz+7719vYORkR3tXVZPqGsWpWTU2c62xARfUAfQHd3d/T09NT80IGBAaZa36zcd31FBHf27WHw0DG6Vizh5t61dX/Mp//N0zHxfXs7h/n2yNI5/75lGQTDwPKK5WXAczW2GZbUCrwBOJphTWZzws96tnqa+L49PDCQyR8dWc4R7AVWS1olqQ3YBOyatM0u4Iby62uAh8Mzb9YgJp717BCweigURGtBmXzfMjsiiIhxSduAB4EW4O6I2C/pdmBfROwC/hb4iqQhSkcCm7Kqx8zMqstyaIiI2A3snvTebRWvXwauzbIGMzObmq8sNjNLnIPAzCxxDgIzs8Q5CMzMEucgMDNLnIPAzCxxmd1rKCuSDgOHptikA3hxivXNyn2nJ9XeU+0bZtf7iohYWm1FwwXB2UjaV+vGSs3Mfacn1d5T7Ruy691DQ2ZmiXMQmJklrhmDoC/vAnLivtOTau+p9g0Z9d50cwRmZnZumvGIwMzMzoGDwMwscU0TBJLWSTogaUjS9rzrqRdJyyU9IukpSfsl3ZJ3TfUkqUXS45IeyLuWepF0vqT7JX2v/O9+Vd411Yukj5a/59+VtFPS6/KuKQuS7pb0Qvm57hPvXSDpXyV9v/zfJXP1eU0RBJJagB3AemANsFnSmnyrqptx4Pcj4i3AWuDmhHoHuAV4Ku8i6uyvgW9GxJuBXyGR/iVdBHwE6I6It1J64FWzPszqy8C6Se9tBx6KiNXAQ+XlOdEUQQBcAQxFxMGIGAP6gY0511QXEfF8RDxWfn2c0i+Fi/Ktqj4kLQPeB9yVdy31IunngHdQerofETEWET/Jt6q6agVeX37G+ULOfA56U4iIf+PM57dvBO4pv74HeP9cfV6zBMFFwLMVy8Mk8suwkqSVwGXAo/lWUjd/BfwRUMy7kDq6BDgM/F15SOwuSYvyLqoeIuKHwF8CzwDPAy9FxL/kW1VddUbE81D6AxB441ztuFmCoNrTnJM6L1ZSO/APwO9FxE/zridrkn4DeCEiBvOupc5agcuBv4mIy4BR5nCIYD4rj4lvBFYBvwAskvSBfKtqDs0SBMPA8orlZTTpIWM1ks6jFAJfjYhv5F1PnVwNbJD0NKWhwHdJ+vt8S6qLYWA4IiaO+u6nFAwpeA/wg4g4HBGngG8Av5ZzTfX0Y0kXApT/+8Jc7bhZgmAvsFrSKkltlCaQduVcU11IEqXx4qci4rN511MvEXFrRCyLiJWU/r0fjoim/+swIn4EPCvpl8pvvRt4MseS6ukZYK2kheXv/btJZKK8bBdwQ/n1DcA/zdWOW+dqR3mKiHFJ24AHKZ1JcHdE7M+5rHq5Grge+I6kJ8rv/UlE7M6xJsvWh4Gvlv/oOQj8ds711EVEPCrpfuAxSmfLPU6T3m5C0k6gB+iQNAx8AvgUcJ+kD1IKxWvn7PN8iwkzs7Q1y9CQmZnNkIPAzCxxDgIzs8Q5CMzMEucgMDNLnIPA7BxVuzNk+f2rJN2ZV11mM+UgMDt3X+bMO0NSfu+b9S3FbPYcBGbnqMadIaF0peu3yjeCe6L8c1jSJ+pcotk5aYori83yJqkDOBURLwFbyu+toHS1+5dzLM3srHxEYDY33gu8ekvk8pOzvg5si4hDuVVlNg0OArO5sZ7T5wfuAL4REd/KqR6zaXMQmM1S+U6YlwJPlJdvBhZHxKdyLcxsmnzTObNzVHlnSODHwBeAt0TEjeX1PwBOASfK/8sdEXFH/Ss1mx4HgdksSfo4pWdm9+ddi9lMOAjMzBLnOQIzs8Q5CMzMEucgMDNLnIPAzCxxDgIzs8Q5CMzMEvf/HuVXz3HNALEAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# hrss\n", "get_hrss = lambda: 2*np.random.randint(0,2) - 1\n", "s = 1\n", "N = 100\n", "zs = np.linspace(1/10, 4, 40)\n", "arr = np.zeros((len(zs),))\n", "for _ in range(N):\n", " a = get_hrss()\n", " b = get_hrss()\n", " arr += [f(a,b,s,z) for z in zs]\n", "arr = arr/N\n", "plt.scatter(1/zs, [f(a,b,s,z) for z in zs], s=5)\n", "plt.xlabel(\"1/z\")\n", "plt.ylabel(\"f\")\n", "plt.grid()" ] }, { "cell_type": "code", "execution_count": 299, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 342, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEGCAYAAABo25JHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAScklEQVR4nO3dcWxdZ3nH8e/jBHdrna1B7jzWBCdI0SBi3Vp7bRgSONBNqZgS/ihSgqjotjZiIsBg0lQ2VKRuEohNbGNEKynrYBvEKx3SMhTRjQ5rmrRWqdtukJaOKODWa6HBzUpttLrpffaHb+itc+0kjs899n2/H8mKzz0n9z6PEvnn9z3nvCcyE0lSuXrqLkCSVC+DQJIKZxBIUuEMAkkqnEEgSYVbW3cB56u/vz83bdq04P6ZmRkuueSSzhW0Qth3eUrtvdS+4cJ6Hx8f/0FmXtZu36oLgk2bNvHAAw8suH9sbIyRkZHOFbRC2Hd5Su291L7hwnqPiImF9jk1JEmFMwgkqXAGgSQVziCQpMIZBJJUOINAkgpnEEhS4QwCSSrcqruhTC/XaCRTM7N1lyFpFXNEsIo1GsmeO+7jDR+7l+M/mKHR8CFDks6fQbCKTc3MMj5xklON5EfPv+jIQNKSGASrWH9fL0OD61nbE1x80Rr6+3rrLknSKuQ5glUsIjh48zamZmY5Ov4fRETdJUlahRwRrHI9PcFl6y6quwxJq5hBIEmFMwgkqXAGgSQVziCQpMIZBJJUOINAkgpXaRBExI6IeCwijkXELW323xgRJyLi4ebXTVXWI0k6U2U3lEXEGmA/8KvAJHAkIg5l5iPzDv37zNxXVR2SpMVVOSK4GjiWmcczcxYYBXZV+HmSpCWIzGpWrIyI64EdmXlTc/sG4JrW3/4j4kbgY8AJ4L+BD2bmE23eay+wF2BgYGBodHR0wc+dnp6mr69vGTtZHey7PKX2XmrfcGG9b9++fTwzh9vtq3KtoXYL38xPnX8CDmbm8xHxHuDzwFvO+EuZB4ADAMPDwzkyMrLgh46NjbHY/m5l3+UptfdS+4bqeq9yamgS2NiyvQF4svWAzJzKzOebm3cAQxXWI0lqo8ogOAJsiYjNEdEL7AYOtR4QEa9q2dwJPFphPZKkNiqbGsrMUxGxD7gHWAPcmZlHI+I24IHMPAS8PyJ2AqeAZ4Abq6pHktRepc8jyMzDwOF5r93a8v2HgQ9XWcNKdfpZw/19vT5HQFKtfDBNDU4/a3h84iRDg+s5ePM2enq6NwwMPWllMwhq0Pqs4fGJk0zNzHbtw2VKCz1pNXKtoRq0Pmt4aHB9Vz9ruF3oSVpZHBHUoPVZw90+XXI69E6PCLo59KTVyiCoSSnPGi4p9KTVyiBQ5UoJPWm18hyBJBXOIJCkwhkEklQ4g0CSCmcQSFLhDAJJKpxBIEmFMwgkqXAGgSQVziCQpMIZBJJUOINAkgpnEEhS4QwCSSqcQSBJhTMIJKlwBoEkFc4gkKTCGQSSVDiDYJk1GsmJ554nM+suRZLOiQ+vX0aNRrLnjvsYnzjJ0OB6Dt68jZ6eqLssSVqUI4JlNDUzy/jESU41kvGJk0zNzNZdkiSdlUGwjPr7ehkaXM/anmBocD39fb11lyRJZ+XU0DKKCA7evI2pmVn6+3qJcFpI0spnECyznp7gsnUX1V2GJJ2zSqeGImJHRDwWEcci4pZFjrs+IjIihqusR5J0psqCICLWAPuB64CtwJ6I2NrmuHXA+4H7q6pFkrSwKkcEVwPHMvN4Zs4Co8CuNsf9IfAJ4P8qrEWStIAqg+By4ImW7cnmaz8WEVcCGzPzKxXWIUlaRJUni9tdMvPj220jogf4U+DGs75RxF5gL8DAwABjY2MLHjs9Pb3o/m5l3+UptfdS+4bqeq8yCCaBjS3bG4AnW7bXAa8HxpqXWf4scCgidmbmA61vlJkHgAMAw8PDOTIysuCHjo2Nsdj+bmXf5Sm191L7hup6r3Jq6AiwJSI2R0QvsBs4dHpnZj6bmf2ZuSkzNwH3AWeEgCSpWpUFQWaeAvYB9wCPAndl5tGIuC0idlb1uZKk81PpDWWZeRg4PO+1Wxc4dqTKWiRJ7bnWkCQVziCQpMIZBJJUOINAkgpnEEhS4QwCSSqcQSBJhTMIJKlwBsESNRrJieeeJzPPfrAkrWA+qnIJGo1kzx33MT5xkqHB9Ry8eRs9PT6fWNLq5IhgCaZmZhmfOMmpRjI+cZKpmdm6S5KkJTMIlqC/r5ehwfWs7QmGBtfT39dbd0mStGRODS1BRHDw5m1MzczS39dL83kKkrQqGQRL1NMTXLbuorrLkKQL5tSQJBXOIJCkwhkEklQ4g0CSCmcQSFLhDAJJKpxBIEmFMwgkqXAGgSQVziCQpMIZBJJUOINAkgpnEJwjn0gmqVu5+ug58IlkkrqZI4Jz4BPJJHUzg+Ac+EQySd3MqaFz4BPJJHUzg+Ac+UQySd2q0qmhiNgREY9FxLGIuKXN/vdExDci4uGI+PeI2FplPZKkM1UWBBGxBtgPXAdsBfa0+UH/xcz8hcz8JeATwCerqkeS1F6VI4KrgWOZeTwzZ4FRYFfrAZn5w5bNSwAv0pekDouqbpCKiOuBHZl5U3P7BuCazNw377j3Ah8CeoG3ZOa327zXXmAvwMDAwNDo6OiCnzs9PU1fX9+y9bFa2Hd5Su291L7hwnrfvn37eGYOt9tX5cnidpfWnJE6mbkf2B8R7wQ+Ary7zTEHgAMAw8PDOTIysuCHjo2Nsdj+bmXf5Sm191L7hup6r3JqaBLY2LK9AXhykeNHgbdXWI8kqY0FgyAi/rb55weW+N5HgC0RsTkieoHdwKF5n7GlZfNtwBnTQpKkai02NTQUEYPAb0bE3zBvqiczn1nsjTPzVETsA+4B1gB3ZubRiLgNeCAzDwH7IuJa4AXgJG2mhSRJ1VosCG4Hvgq8Bhjn5UGQzdcXlZmHgcPzXru15fuljjYkSctkwamhzPxUZr6Oud/kX5OZm1u+zhoC3cClpyWV4KxXDWXmb3eikJXGpacllcLVRxfg0tOSSmEQLMClpyWVwtVHF+DS05JKYRAswqWnJZXAqSFJKpxBIEmFMwgkqXAGgSQVziCQpMIZBJJUOIOgDdcYklQS7yOYxzWGJJXGEcE8rjEkqTQGwTyuMSSpNE4NzeMaQ5JKYxC04RpDkkri1JAkFc4gkKTCGQSSVDiDoIU3kkkqkSeLm7yRTFKpHBE0eSOZpFIZBE3eSCapVE4NNXkjmaRSGQQtvJFMUomcGpKkwhkEklQ4g0CSCmcQ4I1kkspW/MlibySTVLpKRwQRsSMiHouIYxFxS5v9H4qIRyLivyLi3ogYrLKedryRTFLpKguCiFgD7AeuA7YCeyJi67zDHgKGM/MK4G7gE1XVsxBvJJNUuiqnhq4GjmXmcYCIGAV2AY+cPiAzv95y/H3Auyqspy1vJJNUuqjqBGlEXA/syMybmts3ANdk5r4Fjv808L3M/KM2+/YCewEGBgaGRkdHF/zc6elp+vr6lqGD1cW+y1Nq76X2DRfW+/bt28czc7jdvipHBO1+tW6bOhHxLmAYeHO7/Zl5ADgAMDw8nCMjIwt+6NjYGIvtn6/RyK4YDZxv392i1L6h3N5L7Ruq673KIJgENrZsbwCenH9QRFwL/AHw5sx8vsJ6zuAVQ5JU7VVDR4AtEbE5InqB3cCh1gMi4krgM8DOzHy6wlra8oohSaowCDLzFLAPuAd4FLgrM49GxG0RsbN52B8DfcCXIuLhiDi0wNtVwiuGJKniG8oy8zBweN5rt7Z8f22Vn382XjEkSd5Z7NLTkopX9FpDrjEkSQWPCLxiSJLmFDsi8IohSZpTbBB4xZAkzSl2asgrhiRpTrEjgm5ZWkKSLlSRIwJPFEvSS4ocEXiiWJJeUmQQeKJYkl5S5NSQJ4ol6SVFjgg8USxJLyluROCJYkl6ueJGBJ4olqSXKy4IPFEsSS9X3NRQJnxq95VEwGXrLvIcgaTiFRUE7c4PmAOSSlfU1JDnByTpTEUFQX9fL1e9+lLWBFz16ks9PyBJFBYEcw8iC4i5Lx9MJkmFBcHUzCwPPn6SFxvJg04NSRJQWBC88uJXcMWGn2ZN4KWjktRUTBA0Gsk7P3s//zn5LFdsvJQv3nSNl45KEgUFwekrhl5sJN+YfJZnfvRC3SVJ0opQTBB4xZAktVdMEHjFkCS1V0wQnJh+nvGJZ7xiSJLmKSIIGo3kfQcf4sXmKMCpIUl6SRFBMDUzy4MTJwFY0xN8+p1XecWQJDUVEQStS08PD67nsnUX1V2SJK0YRaw+6jOKJWlhlY4IImJHRDwWEcci4pY2+98UEQ9GxKmIuL7KWnp6wucPSFIblQVBRKwB9gPXAVuBPRGxdd5hjwM3Al+sqg5J0uKqnBq6GjiWmccBImIU2AU8cvqAzPxuc1+jwjokSYuocmrocuCJlu3J5muSpBWkyhFBu8n4Jd3PGxF7gb0AAwMDjI2NLXjs9PT0ovu7lX2Xp9TeS+0bquu9yiCYBDa2bG8AnlzKG2XmAeAAwPDwcI6MjCx47NjYGIvt71b2XZ5Sey+1b6iu9yqnho4AWyJic0T0AruBQxV+niRpCSoLgsw8BewD7gEeBe7KzKMRcVtE7ASIiF+OiEngHcBnIuJoVfVIktqr9IayzDwMHJ732q0t3x9hbspIklSTIpaYkCQtzCCQpMIZBJJUOINAkgpnEEhS4QwCSSqcQSBJhTMIJKlwBoEkFc4gkKTCGQSSVDiDQJJWgUYjOdVIMpf0WJdFGQSStMI1GsmeO+7jW089x+4D99FoLG8YGASStMJNzcwyPnGSJBmfOMnUzOyyvr9BIEkrXH9fL0OD6wmCocH19Pf1Luv7GwSStMJFBAdv3sZrX7WO0b3biGj3SPilMwgkaRXo6QnW9sSyhwAYBJJUPINAkgpnEEhS4QwCSSqcQSBJhTMIJKlwBoEkFc4gkKTCRRUr2VUpIk4AE4sc0g/8oEPlrCT2XZ5Sey+1b7iw3gcz87J2O1ZdEJxNRDyQmcN119Fp9l2eUnsvtW+ornenhiSpcAaBJBWuG4PgQN0F1MS+y1Nq76X2DRX13nXnCCRJ56cbRwSSpPNgEEhS4bomCCJiR0Q8FhHHIuKWuuvplIjYGBFfj4hHI+JoRHyg7po6KSLWRMRDEfGVumvplIi4NCLujohvNf/d31B3TZ0SER9s/j//ZkQcjIifqLumKkTEnRHxdER8s+W1V0bEv0TEt5t/rl+uz+uKIIiINcB+4DpgK7AnIrbWW1XHnAJ+NzNfB2wD3ltQ7wAfAB6tu4gO+3Pgq5n5WuAXKaT/iLgceD8wnJmvB9YAu+utqjKfA3bMe+0W4N7M3ALc29xeFl0RBMDVwLHMPJ6Zs8AosKvmmjoiM5/KzAeb3z/H3A+Fy+utqjMiYgPwNuCzddfSKRHxU8CbgL8CyMzZzPzfeqvqqLXAT0bEWuBi4Mma66lEZv4b8My8l3cBn29+/3ng7cv1ed0SBJcDT7RsT1LID8NWEbEJuBK4v95KOubPgN8DGnUX0kGvAU4Af92cEvtsRFxSd1GdkJn/A/wJ8DjwFPBsZv5zvVV11EBmPgVzvwACP7Ncb9wtQdDuac5FXRcbEX3APwC/k5k/rLueqkXErwNPZ+Z43bV02FrgKuAvM/NKYIZlnCJYyZpz4ruAzcDPAZdExLvqrao7dEsQTAIbW7Y30KVDxnYi4hXMhcAXMvPLddfTIW8EdkbEd5mbCnxLRPxdvSV1xCQwmZmnR313MxcMJbgW+E5mnsjMF4AvA79Sc02d9P2IeBVA88+nl+uNuyUIjgBbImJzRPQydwLpUM01dUREBHPzxY9m5ifrrqdTMvPDmbkhMzcx9+/9r5nZ9b8dZub3gCci4uebL70VeKTGkjrpcWBbRFzc/H//Vgo5Ud50CHh38/t3A/+4XG+8drneqE6ZeSoi9gH3MHclwZ2ZebTmsjrljcANwDci4uHma7+fmYdrrEnVeh/wheYvPceB36i5no7IzPsj4m7gQeaulnuILl1uIiIOAiNAf0RMAh8FPg7cFRG/xVwovmPZPs8lJiSpbN0yNSRJWiKDQJIKZxBIUuEMAkkqnEEgSYUzCKTz1G5lyObrb4iIO+qqS1oqg0A6f5/jzJUhab721c6WIl04g0A6TwusDAlzd7p+rbkQ3MPNrxMR8dEOlyidl664s1iqW0T0Ay9k5rPATc3XBpm72/1zNZYmnZUjAml5/Brw4yWRm0/O+hKwLzMnaqtKOgcGgbQ8ruPl5wduB76cmV+rqR7pnBkE0gVqroR5BfBwc/u9wLrM/HithUnnyEXnpPPUujIk8H3gL4DXZeaNzf3fAV4AftT8K7dn5u2dr1Q6NwaBdIEi4iPMPTN7tO5apKUwCCSpcJ4jkKTCGQSSVDiDQJIKZxBIUuEMAkkqnEEgSYX7fwM5zTE2ctrGAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# uniform\n", "get_uniform = lambda: 2*np.random.random() - 1\n", "s = 3**-0.5\n", "N = 100000\n", "zs = np.linspace(1/10, 4, 20)\n", "arr = np.zeros((len(zs),))\n", "for _ in range(N):\n", " a = get_uniform()\n", " b = get_uniform()\n", " arr += [f(a,b,s,z) for z in zs]\n", "arr = arr/N\n", "plt.scatter(1/zs, [f(a,b,s,z) for z in zs], s=5)\n", "plt.xlabel(\"1/z\")\n", "plt.ylabel(\"f\")\n", "plt.grid()" ] }, { "cell_type": "code", "execution_count": 339, "metadata": {}, "outputs": [], "source": [ "# this was a cool strategy, but I think it varies too much..." ] } ], "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.7.4" } }, "nbformat": 4, "nbformat_minor": 4 }