Skip to content

Instantly share code, notes, and snippets.

@hannorein
Created March 31, 2021 13:48
Show Gist options
  • Save hannorein/9230f8d1d7ac8298c1ad085a03c44bed to your computer and use it in GitHub Desktop.
Save hannorein/9230f8d1d7ac8298c1ad085a03c44bed to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "identified-kernel",
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "opposed-loading",
"metadata": {},
"outputs": [],
"source": [
"def pdf(x):\n",
" return 0.004*(x+0.001/x)**-2 # maximum = 1"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "sublime-consortium",
"metadata": {},
"outputs": [],
"source": [
"x = np.logspace(-2,2,100)\n",
"samples = []\n",
"cur = (x[0]+x[-1])/2\n",
"while len(samples)<1000000:\n",
" prop = cur + np.random.normal()\n",
" if prop<x[0] or prop>x[-1]:\n",
" continue\n",
" if np.random.random()<pdf(prop)/pdf(cur):\n",
" cur = prop\n",
" samples.append(cur)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "powered-infrared",
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD8CAYAAAB0IB+mAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAl1ElEQVR4nO3deZxVdf3H8dfnbuOaG1IKIqhIkZommUuWmcuwuwuouZCIhaZmAoqKsg1m7iiiIJq5EJoiq5aRlvgLXFLcAjUENFHcUedun98fjjmNzHhn7nLOvff9fDzmj/ude879zPcxvOfL93zP95i7IyIilS8SdAEiIlIaCnwRkSqhwBcRqRIKfBGRKqHAFxGpEgp8EZEqEQu6gJa0a9fOO3fuHHQZIiJl5Yknnnjb3bdu2h7KwDezvkDfnXbaiSVLlgRdjohIWTGzFetrD+WUjrs/4O5DNttss6BLERGpGKEMfBERKTwFvohIlVDgi4hUiZIFvpntYGZTzWxmqT5TRES+kFfgm9k0M1tjZkubtNea2UtmttzMRgC4+yvuPjifzxMRkbbLd4Q/Haht3GBmUWAS0BPoDgw0s+55fo6IiOQpr3X47v6ImXVu0rwXsNzdXwEws7uA/sDzuZzTzIYAQwA6deqUT3nhM7rRMtPR7wdXh4hUpWLceNUBWNno9Srg+2a2FTAO2MPMRrr7hPUd7O5TgCkAPXr0KM+ns+QS7Ap/ESmxkt1p6+5rgaG5vLfxnbZlY3QzN4k1157TsfpDICKFU4zAXw1s1+h1x4a2ypNLmBfq/Ap/EclTMQJ/MdDVzLrwWdAPAAa15gTu/gDwQI8ePU4tQn35KXbIt+Vz9cdARHKQV+Cb2Z3AAUA7M1sFXOzuU81sGLAAiALT3P25Vp43XFM6QYW8iEgB5btKZ2Az7XOBufmcO3DlFPKa+hGRHIRye+SSTelUYlBW4s8kIgURysAPZEqnnEb0uVL4i0gjoQz8oo7wKzHYRURyEMrAlyLQaF+k6oUy8As+paNRvYhIOPfD1yMOi2z0Zl98iUjVCOUIvyAUZrnRVI9I1QjlCN/M+prZlPffVwCJiBRKKANfUzoB0VSPSEWr3CkdyY+mekQqTihH+CIiUngKfBGRKhHKKZ3Q7ZZZ7TS9I1IRQjnC10XbENOFXZGyFcrAFxGRwgvllI6UCU31iJQVjfBFRKpEyUb4ZrYxcD2QBBa6++9L9dlSAhrti4ReXiN8M5tmZmvMbGmT9loze8nMlpvZiIbmI4CZ7n4q0C+fzxURkdbLd4Q/HbgOuO3zBjOLApOAg4FVwGIzmwV0BJ5teFsmz88NraRHWe1bs9rbsdrbsZZN+dA34iM2JEmcCFkiODUk2cI+YnM+Ymt7j062hu3tTTa2+qB/hPw1XcGjEb9IKOT7EPNHzKxzk+a9gOXu/gqAmd0F9Oez8O8IPE0L/7MwsyHAEIBOnTrlU17RpT3Cc96ZxdlvsjTbmRe8Ey/7tqSbdGuUDJvyMQnSOOBE+IQE69jwS+dsz7vsEnmVXezf7BZ5he9FXmQz+7hEP1GRaLpHJBSKMYffAVjZ6PUq4PvANcB1ZtYbeKC5g919ipm9AfRNJBJ7FqG+vKzMtuMv2T14OLsHi7Pd/hva27CWb0VW8JPIU+wQeZ0OvE1He5t29j4bkMTsy+eq9xjvswlv+ha85u1Z4e1Znu3IUu/MwuzuZDMRImTZxV5l38hzHBR9kj1sGVHzEv/UIlIJSnbR1t3XASfn+N7iPdO2DV7LtmdWdh9mZ/bhRf/sfx1d7A2OiD7KXpEX2SvyIl+391p93hpL0573aG/vsSuv/s/3PvYanvEuLMp8m0XZ7tyc6cXkTD/a8T4HRZ+gX+Qx9o68QEThLyI5Kkbgrwa2a/S6Y0NbzsKwtcI6r+GBzD7cnfkxT3lXAHrYS4yK3c6BkSfZIfKfon7+RlbP3vYie0de5Gzu4QPfkIXZ3Xkw04MHMvtwV+ZAtmEt/aJ/Z0B0IV2KXE/BaHpHJDDFCPzFQFcz68JnQT8AGFSEzymKF7PbcWvmEGZl9mUdG7KzrWRE7A76RB+no70dWF1fs0/oF11Ev+giPvEED2X35L7Mftyc6c2NmX78IPIsx0X/xMGRJ4hZNrA6RSS8zL3tUwJmdidwANAOeBO42N2nmlkv4CogCkxz93FtOX+PHj18yZIlbSuuFXu9ZN1YmP0OUzO9+Ht2Fzagnj7RxxkYfZjv2rL1zr+HxRrfnLszB3Bn+kBepx0deIuTY/M5NrqQTe2ToMvLnUb7IgVjZk+4e48vtecT+MXSaErn1GXLlrXtJDkEftojzM7uw/XpfvzLt+MbrOWnsQcZFH2YzW1d2z43IBk3/pz9Ljene/EP/xab8jHHRx9icGwe7eyDoMv7agp8kYIpq8D/XLFG+CmPcm9mfyZl+vOaf52dbSU/j82id+Rx4lb+twg8k+3ClHQf5mS/Tw0pBkX/zNDYbNq34cJyIBT+Inkpq8Av1gg/7RHuz+7HNenDWeHf4Dv2MsNif+QnkacqcrXLy9ltuD7dj/uyPyBGhhOjDzI09gBb2odBl5Y7hb9Iq5VV4H+uUCN8d3gouyeXpY9luXfk2/Yq58RmcmDkqVDPzxfKimx7rk4fwX3ZH7Ah9fwsOpdTY3PYxD4NurSvpsAXabWqDfwnsl0ZnxrEE96NHex1fh27m9rI4qoI+qaWZ7flivTRzM1+n3a8zy9j9zAg+pfymcZS+IvkpLnAD+V++IVYh78yuzV16YHMye5Ne95lfOxmjokurOoliztFXuf6xNU8ld2RCalBXJg+hVsytYyK3c6PI09X5R9BkWpSkSP86x5exjUPPkeELKdFZ3NabDYbVcKmZAXkDn/Ofpfx6UG84tuyf+QZLoz9jp0jrbpHLjga7Ys0q6xG+PnaKBGjb+Qxfh2fwTfs3aDLCSUzOCj6JD+M/JPbMwdzVfpIeibr+Gn0Qc6K3VP+G7aJyJeEcoRfqnX48oV3fFMuTx/NnZkD2ZIPGR67i6Oij5TH6iWN9kX+R3Mj/FA+4tDdH3D3IZttptAulS3tQ8bHp/FAYhTb25uclz6No5MX8UJ2u68+WETKQigDX4KzS+TfzExcwmWxG3nVt6FPcjxjUsezzmuCLq15ozf74ktEmqXAly+JmHNM7K88XPMrjokuZGqmFwfVX86DmdA9nkBEWkGBL83a3NYxIT6VexIXs5mtY0jqV5yaPIc3fMugS2ueRvsizQpl4JtZXzOb8v77uhgXBntGlvFA4gJGxu7g0eyuHFx/GbelDybrWrgvUk5CuUrnc6XaHlly91q2PRekT+HR7G581/7FxPhNdC2HtftaySNVpKxW6Uh4dYqs4bZ4HVfEr+dV34ZeyQlcnT6cpEeDLq1lmuoRUeBL65nBEdG/8VDNr6mN/IMr00fTNzmOf2Z3CLo0EWlByQLfzHYws6lmNrNUnynF1c4+4NrEddwcv5z3fBMOT15KXWoAn3o86NJEZD1yCnwzm2Zma8xsaZP2WjN7ycyWm9mIls7h7q+4++B8ipVwOij6JA/WnMcx0YVMzvSjV3ICS7I7B11W8zS9I1Uq1710pgPXAbd93mBmUWAScDCwClhsZrP47Dm2E5ocf4q7r8m7Wgmtzexj6uI30yfyOMNTp3J08iJOic7n3NgMNrRk0OU1r2no6+KuVLCcRvju/gjwTpPmvYDlDSP3JHAX0N/dn3X3Pk2+FPZV4gfRpSyoGc7x0T8xNdOLnsk6Fme7BV2WiJDfHH4HYGWj16sa2tbLzLYys8nAHmY2soX3DTGzJWa25K233sqjPAnKJvYpY+LTuSM+ljQRjkleyJjU8XziiaBLE6lqJdse2d3XAkNzeN8UM3sD6JtIJHQvfxnbN/o8CyLDqUsPZGqmF3/J7s5v4jeyZ6SNO6CWQuMpHk3vSIXJZ4S/Gmi8lWLHhjaR/9rY6v872q/3OEcnL2ZCuazk0cVdqTD5BP5ioKuZdTGzBDAAmFWIorQ9cuXZN/o882tGcGx0ITdm+tE3OY5nsl2CLkukquS6LPNOYBHQzcxWmdlgd08Dw4AFwAvADHd/rhBFaS+dyrSpfcKE+M3cGq/jQ9+Qw5OXckXqyPDfpQsa7UtF0F46Eoj3fSMuSZ3Ivdn96W7/5or4DXwzsvKrDwwDze1LyJXVXjoa4Ve+zexjrkjcwOT4FbzpW9AvOZbJ6T5kymEHTo32pUyFMvA1h189aqNLWFAznB9HnqIuPYhjkxexIts+6LJEKlIoA18j/OrSzj5gcvwqroxP4iXvSM9kHbenf0KIZxu/oNG+lJFQBr5G+NXHDA6P/p0Ha4bz3cgyRqUHc1LqPN70zYMuTaRihDLwNcKvXtvYO9wWr+PS2C38X/ZbHFJ/GQ9k9g66rNxotC8hF8rA1wi/ukXM+WnsIeYmRtLF3uCM1JmckRzGe75x0KWJlLVQBr4IwA6R/zAzcQm/is1gXnYvDq2fyF8zuwVdVm402pcQUuBLqMUsyxmx+7gvcRFfs485MTWCi1InlddGbAp/CYlQBr7m8KWpXSL/5oHEBQyOzuW2zCH0To7n6eyOQZclUlZ0p62Unccy3Tk3NZQ32YJh0fsYFruPuGWCLqv1dMeuFElZ3Wkr0pJ9o88zr2YE/SOPcXXmSI5Kjubl7DZBlyUSegp8KUufb81wffwqVnh7eifH87v0QeVxs5ZIQBT4UtZ6Rf/Bgprh7BV5kQvTp3Bianj53Kyli7lSYqEMfF20ldb4ur3HrfGJjIlN4x/Zb3Jo/UTmZb4XdFmto/CXEghl4OvGK2ktMzgh9ifmJM5ne1vD6amzOSc5lA98w6BLEwmNUAa+SFvtGHmDmYnRnBm9h/uz+9Gzvo7Hs98MuqzW0WhfikSBLxUnbhnOid/DzMRo4pZhYHIUE1IDqfdY0KWJBKqk/wLM7DCgN/A1YKq7P1jKz5fqskfkZeYmRjI2fRw3Zvry1+xuXBW/vnyerAVfHuVr7b7kIecRvplNM7M1Zra0SXutmb1kZsvNbERL53D3+9z9VGAocGzbShbJ3UZWz/j4NKbFL+Nt35x+ybHclO5FthyerCVSYDnfaWtmPwQ+Am5z910a2qLAv4CDgVXAYmAgEAUmNDnFKe6+puG43wK/d/cnW/pM3WkrhbTWN2VE6lQeyvZgn8hzXB6fTAdbG3RZbafRvjQj7ztt3f0R4J0mzXsBy939FXdPAncB/d39WXfv0+RrjX1mIjDvq8JepNC2sg+ZEr+Cy2I38s/sjtTW13FfZr+gyxIpmXwv2nYAGk+Irmpoa84ZwEHAUWY2dH1vMLMhZrbEzJa89dZbeZYn8r/M4JjYX5mXGEFXW81ZqV9wRnIY75fjXvtazSOtVNKLtu5+DXDNV7xnipm9AfRNJBJ7lqYyqTbbR9YwI3EpN2T6cXX6CJbUd+O38RvYN/p80KWJFE2+I/zVwHaNXndsaBMJvc/32r83cTEbWj2DUqMYkzqeTz0edGkiRZFv4C8GuppZFzNLAAOAWfkWpTttpZR2i7zKnMT5HB99iKmZXvRPjuGF7HZffWCYaHpHctCaZZl3AouAbma2yswGu3saGAYsAF4AZrj7c/kWpb10pNQ2tCRj47dwS/wy1vrX6J8cy5R0by3flIqiB6CINLHWN2Vk6mc8mP0ee0ee47flvHxTSzerUlk9AEUjfAnSVvYhN8av5LLYjTyb3YHa+jruz+wbdFkiedMIX6QFK7LtOTv1c570nekX+Ttj4tPZzNYFXVbbaLRfNTTCF2mDz5dv/io2g7nZ71NbX8djme5BlyXSJqEMfK3SkTBpunzzuNT5jEsN0u6bUnZCGfga4UsY7RZ5ldmJCzgu+mduyvShf3IML5bT8k0t3ax6msMXaYOHM7tzXmoIH7AR58Xu5pTofCIW3n9LLdLcfsUpqzl8kbA7MPo0C2qG86PIM4xNn8AJqZG84VsGXZZIixT4Im30+e6bdbEpPJXdiUPrJzI78/2gyxJpVigDX3P4Ui7MYEBsIXMTI9nBXmdY6pecnTy9vB6errn9qhHKwNcqHSk3nSNvMjNxCWfFZjIruy896+v4v3J7eLpUvFAGvkg5ilmWs2L3/vfh6QOSo6hLDSDp0aBLy51G+xVNgS9SYHtEXmZOYiQDon9hcqYfhyXHsCzb0nOBREojlIGvOXwpdxtbPRPiU7kpfjlv+hb0SY5jevqQ8tp9U6P9ihPKwNccvlSKg6NPMr9mOPtFljI6fRInpobzpm8edFlSpUIZ+CKVZGv7gKnxyxkTm8bibDcOrZ/I/Mz3gi5LqpACX6QEzOCE2J+YkzifTraGoamzOTd1Gh+Wy/JNTe9UBAW+SAntGHmDexKjOSP6R+7N7E+v5ASWZHcOuqzWUfiXrZIFvpl9y8wmm9lMMzu9VJ8rEjZxy/Cr+B/4Q+ISAI5JXsTlqaNJldPyTSlLOQW+mU0zszVmtrRJe62ZvWRmy81sREvncPcX3H0ocAywX9tLFqkMe0aWMS8xgiOjj3Bd5nCOTI7m5ew2QZclFSzXEf50oLZxg5lFgUlAT6A7MNDMupvZrmY2u8lX+4Zj+gFzgLkF+wlEytgm9im/iU9hcvxKXvP29E6O53fpgwjxJrb/S9M7ZSXn7ZHNrDMw2913aXi9DzDa3Q9teD0SwN0n5HCuOe7eu5nvDQGGAHTq1GnPFStW5FTfl+gXUMrMm74556aG8mh2Nw6MPMnE+BS2tg+CLqvttO1yYIqxPXIHYGWj16sa2por4AAzu8bMbqSFEb67T3H3Hu7eY+utt86jPJHy8nV7j1vjExkdm87fsrtQWz+RhzLfDbosqSAlu2jr7gvd/Ux3P83dJ7X0Xt1pK9UqYs5JsQeZnbiAr9u7nJo6l5Gpn7HOa4IuTSpAPoG/Gmj8fLeODW0ikqedI6v5Y+IiTovO4q7MAfROTuCp7I5Bl9U6mt8PnXwCfzHQ1cy6mFkCGADMKkxZIlJjaUbG7+LOxFhSHuWo5GiuTh9O2nX7jLRNrssy7wQWAd3MbJWZDXb3NDAMWAC8AMxw9+cKUZT20hH5wt6RF5lXM4K+kUVcmT6ao5MXsyLbPuiyWkej/VAI5UPMzawv0HennXY6ddmyZW07iX6xpALNyuzDqNQppIlycew2jokuxMpoA05Aq3dKoLlVOqEM/M/16NHDlyxZ0raDFfhSoV73LflV6nQWZb/NIZHFTIjfzFb2YdBltY3CvyiKsSyzaLRKR6R529o7/D4+ngtit7MwuzuH1k/kL5nvBF2WlIFQBr7m8EVaFjHn1Nhc7k+MYiv7gJNTw7kodRKfeCLo0iTEQhn4IpKbb0VWcn/iQgZH53Jb5hD6JMexNNs56LIkpEIZ+JrSEcndBpbiwvjt3B4fzzrfgMOSl3J9ui+ZcnicolbvlFQoA19TOiKt94PoUubXjODQyBIuSw9kYHIUK7Ptgi5LQiSUgS8ibbO5reO6+DVcEb+e5317eibruDfzg/LZfVOKKpSBrykdkbYzgyOif2NeYgTfstc4J/VzhqXO4D3fOOjSJGBahy9SwTJu3JjpwxXpo2nH+/w2Ppn9ogW5Ib64tD4/L2W1Dl9ECiNqzs9jD/DHxMVsZPUcl7qAsanj+NTjQZcmAVDgi1SBXSOvMidxPidEH+TmTG/6J8fwQna7rz5QKooCX6RKbGhJxsSnc0v8Mtb61+ifHMvN6V5kw7h8U8s1iyKUga+LtiLF8+Po0yyoGc6PIv9kbPp4jk+N5HXfMuiypAR00VakSrnD3ZkDuDT9U2JkGBefSt/o40GX1TJdzM2JLtqKyP8wgwGxhcxNjGQHe50zUmdydvJ0PvANgy5NikSBL1LlOkfeZGbiEs6KzWRWdl961tfxf9lvBl2WFIECX0SIWZazYvcyMzGauGUYkBxFXWoASY8GXZoUUEkD38w2NrMlZtanlJ8rIrnZI/IycxIjOTa6kMmZfhyevJTl2W2DLusLWr2Tl1yfaTvNzNaY2dIm7bVm9pKZLTezETmcajgwoy2FikhpbGz11MVvZkr8t7zhW9E7OZ5b04doP54KEMvxfdOB64DbPm8wsygwCTgYWAUsNrNZQBSY0OT4U4DvAM8DG+RXsoiUwiHRJ9g9spzzUqdxcfokHs7uzm/iU2hv7wVd2mcaj/K1eicnOS/LNLPOwGx336Xh9T7AaHc/tOH1SAB3bxr2nx8/DtgY6A58Ahzu7tn1vG8IMASgU6dOe65YsaKVP1ID/ZdPpCDc4fbMQYxNH89G1DMhfhO10TYuly4FhX9RlmV2AFY2er2qoW293P0Cdz8LuAO4aX1h3/C+KcAlwJOJhB7XJhI0Mzgh9ifmJM6ng73N0NQ5DE+dyjqvCbo0aaWSr9Jx9+nuPvsr3qMHoIiEzE6R17k3cRG/iN7HjMyP6JWcwBPZrkGXJa2QT+CvBhrvvtSxoS1v2lpBJJwSluHX8RncnRhD2qMcnbyYK1JHkgrT8k2t5GlWPoG/GOhqZl3MLAEMAGYVpiwRCbO9Ii8xv2YEh0X+xjWZIzkqOZpXst8IuqwvU/j/j1yXZd4JLAK6mdkqMxvs7mlgGLAAeAGY4e4FebKCpnREwm9T+4QrEpOZFL+af/vX6Z0czx3pA7V8M8RCuXmamfUF+u60006nLlu2rG0n0V90kZL5j2/Buamh/C27KwdFnqAufhPt7IOgy2pZBa/mKavN0zTCFykv37B3uS1ex0Wx23gkuyu19RN5OLN70GXlr8KmhEIZ+LpoK1J+IuacEpvPA4lRtLP3OCV1HqNSJ/OJV8jy6goI/1AGvkb4IuWrW2QV9ycuZEh0Nr/P/ITeyfE8k+0SdFlC7lsriIjkrMbSnB+/gwMiT/Or1OkckbyEs2L3cHp0FlEL33XDch61t0YoR/ia0hGpDPtGn2d+zQh6Rv7B5eljOTZ5ESuzWwddVtUKZeBrSkekcmxm67g2cR1Xx6/jJe9Iz+QEZmb21/LNAIQy8EWk8vSPPsa8mpF82/7NuanT+UXql7zrmwRdVtuU6QXcUAa+pnREKlNHe5s7EuMYEbuDh7J7Ultfx6OZXYIuq2qE8sarz/Xo0cOXLGnjNqxl9pdXpNo8l92eX6Z+wXLvyMnReQyP3cUGlgq6rNYL4Q1cZXXjlYhUvm9HVjA7cQEnRedzS6Yn/ZJjeT7bKeiyKpoCX0QCs4GlGB2/jVvjdbzrm3BYcgxT0r3JugVdWkVS4ItI4H4UfYYFNSP4ceQpxqePY1DqfFb7VkGXVXFCGfi6aCtSfba0D5kcv4rLYjfybHYHauvruD+zb9BlVZRQBr7W4YtUJzM4JvZX5iZG0tVW88vUMM5M/oL3faOgS6sIoQx8Ealu20fWMCNxKefE/sCc7N70rK/jsUz3oMsqewp8EQmlmGU5M/ZH7kmMpsZSHJc6n/GpQdS7tgBrq5IFvpkdYGaPmtlkMzugVJ8rIuVt98jLzEmcz8Dow0zJ9OGw5KX8K9sh6LLKUq6POJxmZmvMbGmT9loze8nMlpvZiK84jQMfARsAq9pWrohUo42snvHxaUyN/4Y1vgV9kuOYlq4Nx/LNMtpmIdcR/nSgtnGDmUWBSUBPoDsw0My6m9muZja7yVd74FF37wkMBy4p3I8gItXiJ9GnmF8znP0jS7k0/VNOTA3nTd886LLKRk6B7+6PAO80ad4LWO7ur7h7ErgL6O/uz7p7nyZfa9w923Dcu0BNwX4CEakqW9sH3By/nHGxm1mc7cah9ROZl/le0GWVhXzm8DsAKxu9XtXQtl5mdoSZ3Qj8DriuhfcNMbMlZrbkrbfeyqM8EalUZnBc7GHmJM6nk63h9NTZnJs6jQ99w6BLC7WSXe5293uBe3N43xQzewPom0gk9ix+ZSJSrnaMvME9idFcmz6c6zKH8Xj2W1wVv54ekX8FXVoo5TPCXw1s1+h1x4Y2EZGSiVuGc+Iz+UPiEgw4JnkRv00dTcqjQZcWOvkE/mKgq5l1MbMEMACYVYiidKetiLTWnpFlzEuM4MjoI1ybOZwjk6N5ObtN0GWFSq7LMu8EFgHdzGyVmQ129zQwDFgAvADMcPfnClGU9tIRkbbYxD7lN/Ep3BC/kte8Pb2T47k9/RM9TrGBHoAiIhXpTd+cc1NDeTS7GwdGnmRifApb2wfF/dCQPAylrB6AohG+iOTr6/Yet8YncnHsVv6W3YXa+on8KfPdoMsKVCgDX3P4IlIIEXNOji1gduIC2tu7/Cx1LiNTg/nYq/NWoFAGvkb4IlJIO0dWc1/iIk6LzuKuzI/pnRzPU9kdgy6r5EIZ+Brhi0ih1ViakfG7uCM+jnqPc1RyNFenDyftoYzBogjlT6oRvogUyz7RF5hXM4I+kce5Mn00RycvZkW2fdBllUQoA18jfBEpps3sY65OTOLq+LUs923plZzAjPSPKn75ZigDX0SkFPpHFzG/ZgS7Rl7hvPRpDE2dxTu+adBlFY0CX0SqWgdbyx3x8Zwf+z1/ye7BofV1/DWzW9BlFUUoA19z+CJSShFzhsTmcF/iQrawjzgxNYKLUyfyqceDLq2gQhn4msMXkSB0j7zGrMQoTonO5dbMofRJjmNpdvugyyqYUAa+iEhQNrAUF8Vv53fx8XzoG3F4cgw3pPuSCcPjFPOkwBcRWY/9o0uZXzOCgyNPMDE9kIHJUazydkGXlZdQBr7m8EUkDLawj5gUv5rL4zfwvG9Pz/o67svs1/zyzZA/0DyUga85fBEJCzM4Kvoo8xIj6WYrOSv1C85MDeN93zjo0lotlIEvIhI220Xe4u7EpZwbu5t52b2ora/jsUz3oMtqFQW+iEiOouYMi93PvYmL2dDqGZQaxbjUIOq9ZI8Hz4sCX0SklXaLvMrsxAUcH32ImzJ96J8cw4vZ7b76wICVLPDNLGJm48zsWjM7sVSfKyJSDBtZPWPjtzAtfhlv++b0S47l5nRPsiFevpnrM22nmdkaM1vapL3WzF4ys+VmNuIrTtMf6AikgFVtK1dEJFwOjD7N/Jrh/DDyDGPTJ3BCaiRv+JZBl7VeuY7wpwO1jRvMLApMAnoC3YGBZtbdzHY1s9lNvtoD3YDH3P0c4PTC/QgiIsFqZx9wU/y3jI/dzJPZnaitr2P2M68HXdaX5BT47v4I8E6T5r2A5e7+irsngbuA/u7+rLv3afK1hs9G9e82HJtp7rPMbIiZLTGzJW+99VbrfyIRkQCYwaDYw8xJnE9n+w/D7niKc+5+mg8+TQVd2n/lM4ffAVjZ6PWqhrbm3AscambXAo809yZ3n+LuPdy9x9Zbb51HeSIipbdD5D/MTFzCmT/pyn1Pr6bnVY/yj1ebjpeDUbKLtu7+sbsPdvcz3H1SS+/VnbYiUs7iluGcg3fmD0P3JRY1jp2yiMvmv0gynQ20rnwCfzXQeB1Sx4Y2EREB9tx+C+aeuT/H9tiO6xe+zBE3/J3laz4KrJ58An8x0NXMuphZAhgAzCpEUdpaQUQqxcY1MeqO3I0bT9iT1e9+Qp9rH+V3i/6NB/A8xVyXZd4JLAK6mdkqMxvs7mlgGLAAeAGY4e7PFaIoTemISKU59NvfYMHZP2TvHbbiwvuf4+Tpi1nz4aclrSHXVToD3X0bd4+7e0d3n9rQPtfdd3b3Hd19XHFLFREpb+033YBbTvoeY/p/m0Uvr6X2qkd58Ln/lOzzQ7m1gqZ0RKRSmRkn7NOZOWf+gG0334Ahv3uCEfc8w7r6dNE/O5SBLyJS6XZqvyn3nr4fPz9gR+5espLe1zzKk6+9+9UH5iGUga85fBGpBolYhPNqv8ldp+5NKuMcPXkRVz70L9KZ4izfDGXga0pHRKrJ93fYinln7U//72zL1X9exlGTF7Fi7bqCf04oA18jfBGpNl/bIM4Vx+7OdYP24PX3PuHTVOFH+aEMfI3wRaRa9dltWx4d/mO6fWPTgp87lIEvIlLNamLRopw3lIGvKR0RkcILZeBrSkdEpPBCGfgiIlJ4CnwRkSqhwBcRqRKhDHxdtBURKbxQBr4u2oqIFJ4FsQl/rszsLWBFw8vNgKZD/qZtjV+3A94uUmnrq6VQx7T0vua+l0vfrK8tzP2V63GF6q/1tau/Wv5etfdXS98Pur+2d/cvPxTc3cviC5jyVW2NXwNLSllLoY5p6X3NfS+Xvim3/sr1uEL111f1TzX3V3Pfq/b+aun7Ye2vUE7pNOOBHNrW955iaMvn5HpMS+9r7nu59M362sLcX7keV6j+Wl+7+qvl71V7f7X0/VD2V6indPJhZkvcvUfQdZQL9VfrqL9aR/3VOsXqr3Ia4bfWlKALKDPqr9ZRf7WO+qt1itJfFTvCFxGR/1XJI3wREWlEgS8iUiUU+CIiVaIqA9/MDjOzm8zsbjM7JOh6ws7MdjCzqWY2M+hawsrMNjazWxt+r44Lup6w0+9U6xQqs8ou8M1smpmtMbOlTdprzewlM1tuZiNaOoe73+fupwJDgWOLWW/QCtRfr7j74OJWGj6t7LsjgJkNv1f9Sl5sCLSmv6r1d6qxVvZXQTKr7AIfmA7UNm4wsygwCegJdAcGmll3M9vVzGY3+Wrf6NBRDcdVsukUrr+qzXRy7DugI7Cy4W2ZEtYYJtPJvb+kbf2VV2bF2npgUNz9ETPr3KR5L2C5u78CYGZ3Af3dfQLQp+k5zMyAOmCeuz9Z5JIDVYj+qlat6TtgFZ+F/tOU50Aqb63sr+dLXF7otKa/zOwFCpBZlfKL2YEvRlfw2T++Di28/wzgIOAoMxtazMJCqlX9ZWZbmdlkYA8zG1ns4kKuub67FzjSzG6gdFsKlIP19pd+p5rV3O9XQTKr7Eb4heDu1wDXBF1HuXD3tXw2dyjNcPd1wMlB11Eu9DvVOoXKrEoZ4a8Gtmv0umNDm6yf+qvt1Heto/5qnaL2V6UE/mKgq5l1MbMEMACYFXBNYab+ajv1Xeuov1qnqP1VdoFvZncCi4BuZrbKzAa7exoYBiwAXgBmuPtzQdYZFuqvtlPftY76q3WC6C9tniYiUiXKboQvIiJto8AXEakSCnwRkSqhwBcRqRIKfBGRKqHAFxGpEgp8EZEqocAXEakSCnwRkSrx/5VELrYnxXCaAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig,ax = plt.subplots(1,1)\n",
"ax.set_xscale(\"log\")\n",
"ax.set_yscale(\"log\")\n",
"ax.plot(x,pdf(x))\n",
"ax.hist(samples,bins=x,density=True);"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "rocky-brief",
"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.9.2"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment