Skip to content

Instantly share code, notes, and snippets.

@hannorein
Created March 31, 2021 13:34
Show Gist options
  • Save hannorein/349881d3bc5ecc94a757ea56ecd22494 to your computer and use it in GitHub Desktop.
Save hannorein/349881d3bc5ecc94a757ea56ecd22494 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",
"while len(samples)<10000:\n",
" a,b = np.random.random(2)\n",
" a = x[0] + a*(x[-1]-x[0])\n",
" if pdf(a)>b:\n",
" samples.append(a)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "powered-infrared",
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD8CAYAAAB0IB+mAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAl80lEQVR4nO3deXhV1dn///d9puCIA1p9CBRERaljpVi1VqtVgghosQqOFSpqRa36U0BUUEBA/TmCWFSKWkVxZC5q1WortuBQBxRL8bHECZwVNGe6v38YH2MkmOQMe59zPq/ryh97Ze+zb9YVPllZe511zN0REZHyFwm6ABERKQ4FvohIhVDgi4hUCAW+iEiFUOCLiFQIBb6ISIWIBV3A+rRr1847deoUdBkiIiXl2Weffd/dt2rcHurA79SpE0uWLAm6DBGRkmJmb66rXVM6IiIVIpSBb2Z9zGzqJ598EnQpIiJlI5SB7+5z3H1I27Ztgy5FRKRshDLwRUQk/xT4IiIVomiBb2bbmdmtZnZfse4pIiLfyCnwzWyama0ys5cbtdeY2TIzW25mwwHcfYW7D87lfiIi0nq5jvCnAzUNG8wsCkwGegHdgIFm1i3H+5Su0W2/+RIRCVBOb7xy9yfNrFOj5h7AcndfAWBmdwP9gKW53KusNfxlMFpLUUWkMArxTtv2wMoGx7XA3ma2JTAO2NPMRrj7+HVdbGZDgCEAHTt2LEB5RdDUaL45o3yFv4gUSNG2VnD3D4DTmnHeVDN7B+iTSCT2KnxlLdRkmBcgnBX+IpJHhQj8t4AODY6r69tKV0tH5iIiIVSIwF8M7GBmnfkq6AcAx7bkBdx9DjCne/fupxSgvtLU3F8o+ktARJqQU+Cb2QzgQKCdmdUCo9z9VjMbCiwEosA0d3+lha/bB+iz/fbb51JeyzQOVAWniJSZXFfpDGyifT4wP5fXFhGR/ArlfvgFndIp9wehxXyoLCIlJZSBHwp6CCsiZSaUgR/IHH65K/e/bETke4Vyt0zthy8ikn8a4VcijfZFKlIoAz8vD201B988Cn+RihHKwC8a/VIQkQoSyjl8fYi5iEj+hTLw9dA2INq7X6SsVfaUjjRNc/siZSeUI3wREck/Bb6ISIUI5ZSO1uGHjKZ3RMpCKEf4emgbYnqwK1KyQjnClxKhkb9ISQnlCF9ERPJPI3zJD432RUKvaIFvZhsBNwJJ4Al3v7NY9xYRkRyndMxsmpmtMrOXG7XXmNkyM1tuZsPrm38F3OfupwB9c7mviIi0XK4j/OnAJOD2rxvMLApMBg4BaoHFZjYbqAZeqj8tk+N9Qyvjxio2523fklpvx4e+KZ+xIZ/7BiSJESFLBKeKFJvb57Tlc7ayj+loq6i21SSsbLtGRAKW64eYP2lmnRo19wCWu/sKADO7G+jHV+FfDbzAev6yMLMhwBCAjh075lJewWXdeN2rWZztysvemdeyHVnm1XxJ1XfObUMdCdI4RhbjSxJkiH7rnAhZqm01P7L/ZZfIG+xmK9gr8m82tLpi/ZPyo/GSTc3pi4RCIebw2wMrGxzXAnsD1wOTzKw3MKepi919qpm9A/RJJBJ7FaC+nLznm/FYZk8ez+7BP7I78wkbA7A5n7Fz5E2OjfyF7ewd2tv7VNv7tLNP2JgviDcauWfd+Jw2fOybsIrNeNN/wJvZrVnu1bziP2RBem8A4qTZ3f7DvpGXOTj6PLvZCsyK/s/OjR7oioRC0R7auvsa4ORmnpv7B6Dk0Vu+JXMz+zA381Ne8u0AaM9qaqKL+UnkNXrYMjrYqhYFccScTfmCTe0LOrKK7rxOwwH/J74RL2S7sCjbjUXZbkzKHMn1mf5swwccEn2WvtFF7GWvEzHP879WRMpVIQL/LaBDg+Pq+rZmC8PWCmu9irmZn3Jv5gAW+04A7G7LuSA2g4Mjz7Oj1RZ0pN3W1nBA9EUOiL4IwEe+MY9l9+SRzF7cl/k5d2QOpdpWcUTk7xwd/SsdI6sKV0w+abQvEphCBP5iYAcz68xXQT8AOLYA9ymI17PtuT1zKLMy+/EZG7Kdvc3/F7uHPpFF/DDAUN3cPqd/9Cn6R59ijVfxSLY7D2b248ZMPyZljmT/yIscF32UX0aeI2bZwOoUkfAy99ZPCZjZDOBAoB3wHjDK3W81s8OAa/lqkmKau49rzet3797dlyxZ0rriWrDXS9aNJ7K7My3Ti79ld6WKJL0jzzAw9jjdbVmo58zf9c25J3Mgd6cP4h22pD2rOTn2Z46JPsEm9kXQ5a2fRvgiBWFmz7p79++05xL4hdJgSueUf//73617kWYEfsqjzMnuw5R0X/7t1WzDB5wQe4Rjo4+xuX3euvsGJO0RHs3+mGnpXvzTd2Zj1nJc9FF+G5vPVvZp0OW1jH4RiOSkpAL/a4Ua4Sc9yn2ZA7gx05da35qu9l9Oj82hd+SZ76ymKUX/ym7HzenezMvuTYI0A6OPcVpsDtvYR0GX1jwKfJGclFTgF2qEn/IoD2T254bMEdT61uxh/2ZobBYHRZ4vy9UuK7LbMCXTlwczPyNCluOjj3J6bHb4R/wKfJGclFTgfy1fI3x3mJ/dm6vSR/OGb8vu9h9+H7uPAyP/CvX8fL6szG7F9ZkjeSCzPwlSnBxdyGmx2Wwa1jl+Bb5ITio28BdldmZ8+lhe9C7saCs5P3YPv4w8VxFB39gb2W24Jt2f2dn92IzPGBp7iBOij1Bl6aBLa5rCX6TFmgr8UO6Hb2Z9zGzqJ5+0/j/7f7Lb8tvkeQxMXcwHvilXxaewIDGcQ6KVGfYAnSPvcn1iMnMTI9g18gZj0yfwy+RVzM/0IMS/90UkT8pyhH/1I69z419epQ1JzojN4uTon2ljqQJUWNqeyuzCuPTxvOYd+Ym9xsXxO9gt8kbQZTVNo32RZimpEX6ufrBpFQOij/NE1bmcHpujsG/C/tGXmZcYweWxW3jDt6VfcgzDUqfwvm8adGkiUgChHOEXax2+fONT34Ab0kfyx0wNG1DHObH7OTH6cHjftavRvkiTSmqE7+5z3H1I27YK7WLZ1L5gZPwu/pwYzh6R/3BZ+kQOT45jSXbHoEsTkTwJZeBLcLaPvM3t8QncFL+GT30jjkqO5vzUED70TYIu7dtGt/3mS0SaRYEv32EGNdHFPFp1PqdHZ/Fg5mccVHcVM9MHkPUKXeIkUgYU+NKkDa2OYfF7mJ8YwQ72FhekT2VA8iKWZ/8n6NJEpBVCGfj5WIcv+bNj5C3uSYxhYmwqy7wDvZITuCbVnzov2ufniEgehHKVzteKtT2yNN/7viljUicwK7sf29nbTIjfTI/IsqDLWjet5JEKVVKrdCS82tmnXJeYzG3xCSSJcXRyFBelTuYz3yDo0kTkeyjwpVUOiL7IwsQwBkfnc1fmYA6pu4K/ZPYMuiwRWY+iBb6ZbWdmt5rZfcW6pxTWRlbHxfE/8UBiFG1tDYNT53N28ozwLOHU0k2Rb2lW4JvZNDNbZWYvN2qvMbNlZrbczIav7zXcfYW7D86lWAmnPSL/YU5iJGdH72dedm8OqbuCuZm9w7Uhm8JfpNkj/OlATcMGM4sCk4FeQDdgoJl1M7NdzWxuo6+t81q1hE7CMpwTv5+5iZG0t/cZmjqb01O/Z7X25REJjWatq3P3J82sU6PmHsByd18BYGZ3A/3cfTxweF6rlJKxU2QlDyRGcXOmN9ek+/NM3ZVcGr+NvpGnw7MtdcNRvlbySAXJZQ6/PbCywXFtfds6mdmWZnYTsKeZjVjPeUPMbImZLVm9enUO5UlQYpbl9Ngc5icupLO9y9mpoZyaOodVrukUkSAV7aGtu3/g7qe5e5f6vwKaOm8qcCnwXCKRKFZ5UgDbR97mvsRoLozdyRPZ3Tm07gpmZfYN19y+SAXJJfDfAjo0OK6ubxP5P1FzhsTmfWu0f5rm9kUCkUvgLwZ2MLPOZpYABgCz81GUtkcuP1+P9kfE7uLx7B4cWnclczN7B13W+mllj5SZ5i7LnAEsArqaWa2ZDXb3NDAUWAi8Csx091fyUZT20ilPUXNOjc1lXmIkHW0VQ1Nnc0byzGDX7TcMdQW7lDntpSOBSHuEP2T6cG26P21Zw+XxWzg0+mzQZTVNq3mkhJTUXjoa4Ze/mGU5IzaL2YmRbG0fMSR1HucmT+MT3zDo0kTKVigDX3P4lWPnyEoeSlzMWdEHmJXdj551E3kys2vQZYmUpVAGvkb4lSVhGc6N38cDiVFsZF9yYmoEF6VOZo1XBV2aSFkJZeBrhF+Zdo+sYF7iQn4bncedmYM5LDk+nB+iroe8UqJCGfga4VeuNpbiovidzEiMJUOEXycvYXxqQPCfrqWQlzKgVToSWp97G8alj2dG5iC62n+5Oj6FH0XeDLqspmklj4RESa3SEQHY2L5kfPwW/hi/go98E45IjmFyuh9p14+tSGvof46E3i+iL7CwahiHRhZzZfoYjk5ewhvZbYIu67s07SMhF8rA1xy+NLa5fc7kxA1cH7+B5d6ew5KXc0f6l+HdiE3hLyEUysDXKh1pSt/oIh6uuoCfRJZxcXoQJ6WG8a5vHnRZIiUhlIEvsj7b2EfcFp/ImNg0Fme70rNuIrMz+wRdlkjoKfClJJnBCbFHmZ+4kO3sHc5KncmZyaF87BsFXZpIaIUy8DWHL83VOfIu9yYu5bzYTBZke4RzawbN50tIhDLwNYcvLRGzLGfGHuLBxCg2sbWcmBrBJanf8IXrE9NEGgpl4Iu0xq6RN5ibGMmg6HxuzxxK7+TlvJDtEnRZ36bRvgRIgS9lpY2luCT+J+6Kj+VLT9A/OZprUv1JeTTo0kQCp8CXsrRvdCkLqobTN/I012X6c1RyNP/Jbht0WSKBKmrgm9kRZnazmd1jZocW895SedraWq5JTGFy/Dre9K3pHfY3a4kUWLMD38ymmdkqM3u5UXuNmS0zs+VmNnx9r+HuD7n7KcBpwDGtK1mkZXpH/8HCqmH0iLz2f2/Wes83C7oskaJryQh/OlDTsMHMosBkoBfQDRhoZt3MbFczm9voa+sGl15Uf51IUfzAPv6/N2v9M7sTPesmMi+zd9BliRRVi7ZHNrNOwFx336X+eB9gtLv3rD8eAeDu45u43oAJwCPu/mgT5wwBhgB07NhxrzffbOV2uFoFIU1Ykd2Gc1Jn8C/vwhGRv3FpfDptbW3QZX1D2yxLjpraHjnXT5VoD6xscFwLrG/YdCbwS6CtmW3v7jc1PsHdpwJT4av98HOsT+Q7tou8y32J0UzO9OOG9JH8o25n/v/4FPaNLg26tO9qPHDRLwPJQVEf2rr79e6+l7uftq6w/5reaSuFFrcMv489wAOJUWxgdRybuoixqeP40uNBlyZSMLkG/ltAhwbH1fVtIiXh68/RPTH6MLdketM3OZZXsj8MuiyRgsg18BcDO5hZZzNLAAOA2bmXJVI8G1iSy+LTmR6fwMe+MUckx3Bjug8Zt6BLE8mrlizLnAEsArqaWa2ZDXb3NDAUWAi8Csx091dyLUp76UgQDoy+yMKqYRwSeZYr0gM5JnkJK7NbBV2WSN40O/DdfaC7b+vucXevdvdb69vnu/uO7t7F3cfloyjN4UtQNrfPmRy/jqvjN7LMO1CTnMDM9AF6s5aUhRYtyyy27t27+5IlS1p3sZZlSo7e8i05L3Uaz2R/xKGRxYyP38KW9lmwRbV0lU7D/wda4VMxmlqWGcq9dDTClzBobx9wV/xyLozdyRPZPehZdwWPZfYIuiyRVgtl4GsOX8IiYs6Q2DxmJy6inX3MoNQFjEgNZq1XBV2aSIuFMvBFwmanyEpmJS7m1Ogc7s78gsOS43k+bHvti3yPUAa+pnQkjKoszYj4DGYkxpLyKEclR3N1GPfa14esSBNCGfia0pEw+2nkNRZUDadf5Gmu1177UkJCGfgiYbepfcHVQey1r9G75CDXzdMKwsz6AH223377oEsRWa/e0X/QPbKM81OncnF6EI9mf8yV8alsbR8XtxD9ApBm0Dp8kTxwhzsyhzAufRwbUsf4+C3URBcHXda3aR1+xSipdfgipcYMTow9wrzEhVTbak5LncN5yVP5zDcIurTvp2miiqHAF8mj7SNv80BiFGdGH+TB7P7U1E3gn9muQZclAijwRfIubhnOi9/LvYlLiVmGY5IXMyE1gDoP+JGZRvIVL5SBr3X4Ug72ivyb+YkRHBN9gpsyfTkieRmvZ9sHXZZUsFAGvtbhS7nYyOqYEL+Fm+NXsco35/DkOG5N15DVXvsSgFAGvki5OST6HH+uGsb+kZcZkz6RE1IjeMe3CLqs79K0T1lT4IsUyVb2KbfEr2J87Gaez25Pz7qJzMrsE3RZUkGKFvhmtrOZ3WRm95nZ6cW6r0iYmMHA2OPMT4ygi73N2akzOSt5Bp/4RkGXJhWgWYFvZtPMbJWZvdyovcbMlpnZcjMbvr7XcPdX3f004Ghgv9aXLFL6OkXe497EpZwbu5d52Z9SUzeBv2d+FHRZ36bpnbLT3HVi04FJwO1fN5hZFJgMHALUAovNbDYQBcY3un6Qu68ys77A6cAdOdYtUvJiluWs2IMcEPkX56R+x3GpkQzOzuf82D20sVRhb64Qr0jNGuG7+5PAh42aewDL3X2FuyeBu4F+7v6Sux/e6GtV/evMdvdewHH5/EeIlLLdIyuYl7iQE6MPc2vmMPomx7I02zHosqQM5TKH3x5Y2eC4tr5tnczsQDO73sz+AMxfz3lDzGyJmS1ZvXp1DuWJlI4NLMll8elMj0/gI9+Yfsmx3JQ+nEzYl29q2qekFO2tf+7+BPBEM86bambvAH0SicReha5LJEwOjL7IwshwLkwNZkL6WB7L7MnViSlU2/tBl/YNhXvJymWE/xbQocFxdX2biORgC/uMKfFruSo+haX+Q3rVTeD+zP6F32tfyl4ugb8Y2MHMOptZAhgAzM5HUXqnrVQ6Mzgq+hQLEiPYyf7LeanTOSN1Nh/5xkGXJiWsucsyZwCLgK5mVmtmg909DQwFFgKvAjPd/ZV8FKW9dES+0iGymrsTYxgWm8Ej2b3oWTeRv2Z2C7osKVHNXaUz0N23dfe4u1e7+6317fPdfUd37+Lu4/JVlEb4It+ImnN6bA4PJi6mra3hpNRwRqVO4gtPBF2alJhQbq2gEb7Id+0SeZM5iZEMis7ntkxPDk+O46Vs56DLkhISysDXCF9k3dpYikvif+LO+DjWeBuOTF7KpHS/8C/flFDQh5iLlKD9oq+wMDKci1Inc1X6GB7P7ME18Sl0jKwqzA21FLMsaIQvUqLa2hpuSEziuvgkXvdqeiXHMzN9QGkt39Qbt4oqlIGvOXyR5usXfZqFVcPYLbKCC9KnMiR1Lh/4JkGX9W0K9lAIZeBrhC/SMv9jH3Jn/HIuiv2Jv2Z3o2fdRB7L7BF0WRIyoQx8EWm5iDm/jc1nTuIi2tknDEpdwIWpQaz1qqBLk5BQ4IuUma6RWmYlLubU6BxmZA6id/Jyns92CbosCYFQBr7m8EVyU2VpRsRnMCMxlqTHOCo5mmtS/Ul5NOjSJEChDHzN4Yvkx08jr7Ggajj9Ik9zXaY/RyVHsSK7TdBlSUBCuQ5fRPJnU/uCqxNTODjzHBemBnNYcjwjY3dyfPRRLJ/v12q8Amd0C/9Cb3h9S6+VZgnlCF9E8q939B8srBrGTyLLuDg9iEGp81nlIf0rWss4CyKUga85fJHC2MY+4rb4REbHpvN09kfU1E1kYaZ70GVJkYQy8DWHL1I4EXN+E3uYuYmRbGsfcGrqXC5IncLn3ibo0qTAQhn4IlJ4O0Te4sHEJZwRfYj7MgfQKzmBJdkdgy5LCkiBL1LBEpbh/PhMZiYuA+Do5CVcmTqaZCks39Q8f4sp8EWE7pHXWZAYzlHRvzI5cwS/Sl7G8uz/BF2W5FlRA9/MNjKzJWZ2eDHvKyLfb2P7kiviN3NT/Gre8nb0Tl7O9PShZLXXftlo7mfaTjOzVWb2cqP2GjNbZmbLzWx4M15qGDCzNYWKSHHURJewsOoC9oksZXT6N5yUGsa7vnnQZX1F0zg5ae4br6YDk4Dbv24wsygwGTgEqAUWm9lsIAqMb3T9IGB3YCmgpQAiIbe1fcIf41dwZ+ZgxqaPp2fdRC6P30rv6D+a/yIK5dBpVuC7+5Nm1qlRcw9gubuvADCzu4F+7j4e+M6UjZkdCGwEdAO+MLP57p5dx3lDgCEAHTt2bPY/RETyywyOj/2FfSOvcE7qDM5Inc1fMk8xOj6dTe2LoMuTVshlDr89sLLBcW192zq5+0h3/z1wF3DzusK+/ryp7t7d3btvtdVWOZQnIvmwXeRd7kuM5qzo/czK7kuvugk8k90p6LKkFYq+Ssfdp7v73PWdo3faioRL3DKcG7+fexOXErcMA5MXMT41kDrXdlylJJfAfwvo0OC4ur5NRMrUjyPLmZcYwYDo4/wh04d+yTG8lu3w/RdKKOQS+IuBHcyss5klgAHA7HwUpa0VRMJrI6tjfPxWbo1fyfu+GX2TY7klfVj4l282XOFToQ+Um7sscwawCOhqZrVmNtjd08BQYCHwKjDT3V/JR1Ga0hEJv4Ojz/PnqmH8PPIvxqaP5/jUCN72LfJ/owoP6Xxq7iqdgU20zwfm57UiESkZ7exTbo5fzT2ZA7ksfSI1dRMZE59Gv+iioEuTdQjl1gqa0hEpHWYwIPYE8xMj6GJvc3bqTM5MDuUT3yjo0qSRUAa+pnRESk+nyHvcm7iU82IzWZDtQU3dBP6e+VHQZeWmzKaTQhn4GuGLlKaYZTkz9hAPJEaxgdVxXGokY1LH86XHgy5NCGngi0hp2y3yBvMSF3Ji9GFuzRxG3+RYlmaL9M75MhuV51MoA19TOiKlbwNLcll8OtPjE/jIN6Zfciw3pQ8nk/WgS6tYoQx8TemIlI8Doy+ysGo4B0eeY0L6WAbe/AwrP1wbdFkVKZSBLyLlZQv7jCnxa7kqPoWlb39Kr+ue4v5na3EN9osqlIGvKR2R8mMGR0WfYsHZ+7Pztptw3r3/4ozU2XzkGwddWsUIZeBrSkekfHXYYkPuHrIPF9R05ZHsXvSsm8hfM7sFXVZFCGXgi0h5i0aM3x24PQ8mLqatreGk1HBGpU7iC08EXVpZU+CLSGB2ibzJnMRITo4u4LZMTw5PjuOlbOegy1q3MljuqcAXkUC1sRSj4nfwp/jlrPE2HJm8lEnpfmTCvvtmCQpl4OuhrUjl+Vn0Zf5cNZyayGKuSh/D0clL+G9266DL+n4lNPIPZeDroa1IZdrM1nBD/Aaui0/ida+mV3I896QP1PLNPAll4ItI5TKDftGn+XPVcHaLrGBYeginps7hA98k6NJKngJfREKpvX3AnfHLGRn7E09k96Bn3UQey+wRdFklTYEvIqEVMeeU2HxmJS6inX3KoNQFjEwNYq1XBV1aSSpa4JvZgWb2lJndZGYHFuu+IlL6do6s5KHExQyJzuWuzEH0Tl7O89kuQZdVcpr7mbbTzGyVmb3cqL3GzJaZ2XIzG/49L+PA50AboLZ15YpIpWpjKS6M38Vd8XHUeZyjkqO5JtWftJfwREWRV/g06zNtgenAJOD2rxvMLApMBg7hqwBfbGazgSgwvtH1g4Cn3P2vZvYD4GrguNxKF5FKtE/0VRZEhjM6dRLXZfrzRHZ3rn1/DZ3blfhHKjYM/dGFWZLerF+N7v4k8GGj5h7Acndf4e5J4G6gn7u/5O6HN/pa5e7Z+us+ApqcgDOzIWa2xMyWrF69uhX/JBEpd21tLdckpnBD/Hr+17fhsOue4s5/vIlr/eZ6NXeEvy7tgZUNjmuBvZs62cx+BfQENuOrvxbWyd2nmtk7QJ9EIrFXDvWJSJnrE32G7pHXOb/D3Yx88GUee3UVE3xTtrJPgy4tlIo2+eXuD7j7qe5+jLs/8T3n6o1XItIs29qH3D6oB6P6dONvy9+nZ90VPJzRWHFdcgn8t4AODY6r69typq0VRKQlIhHj5P06M/fMn7GNfciQ1HkMS53C53XpoEsLlVwCfzGwg5l1NrMEMACYnZ+yRERabocfbMJDiYs5PTqLmZkDOOy6p3j2zcaPHytXc5dlzgAWAV3NrNbMBrt7GhgKLAReBWa6+yv5KEpTOiLSWgnLMCx+D/ckxpB159c3LeKqhctIeTTo0gLX3FU6A919W3ePu3u1u99a3z7f3Xd09y7uPi5fRWlKR0Ry1SOyjAVn70//H1cz6fHl/Cp5Kcuz/xN0WYEK5TsWNMIXkXzYpE2cK3+9Ozcdvxe1vhW9k5dzW/rQit19M5SBrxG+iORTzS7bsLBqGPtEljIq/RtOSg3jvU+/DLqsogtl4GuELyL5trV9zB/jVzAmNo1/Znei57VPsuCld4Iuq6hCGfgiIoVgBifEHmVe4kJ+uMWGnH7nc5w78wU++zIVdGlFEcrA15SOiBRSl8g73Hf6vpx18A7MeuFtaq59in++Uf7LN0MZ+JrSEZFCi0cjnHvIjtx72j7EosYxUxcxYcFrJNPZ77+4RIUy8EVEiuXHHTdn/ln7M+AnHbjpr/+h3+S/8/p7nwVdVkGEMvA1pSMixbRRVYzxv9qNW07szqpPv+TwG/7GLU+tIOsWdGl5FcrA15SOiAThl91+wMJzfs7Pd2jH2HmvcnxqBG/7FkGXlTehDHwRkaC027iKm0/szoRf7coL2e2pqZvI7Mw+QZeVFwp8EZFGzIwBPToyPzGCLvY2Z6XO5KzkGXzipf2pWgp8EZEmdIq8x72JSzk3di/zsj+lpm4CT2e6BV1Wq4Uy8PXQVkTCImZZzoo9yAOJUWxgdRybuogxc5fyZSoTdGktFsrA10NbEQmb3SMrmJe4kBOiD3Pr396g36S/s/Tt0vooxVAGvohIGG1gScbEp/PHk3/Ch2uTHDH57/whfTiZElm+qcAXEWmhX3TdmoW//zkH7bQ149PHMjB5EbXeLuiyvlfRAt/MImY2zsxuMLOTinVfEZFC2GKjBFOO/zFXxm5iqf+QXnUTeOC5WjzEm+039yMOp5nZKjN7uVF7jZktM7PlZjb8e16mH1990HkKqG1duSIi4WFm/Dr2JAsSI9jJ/su5M//F0Lue56M1yaBLW6fmjvCnAzUNG8wsCkwGegHdgIFm1s3MdjWzuY2+tga6Ak+7+7nA6fn7J4iIBKtDZDV3J8ZwQU1XHl76Lj2vfZInX18ddFnf0dzPtH0SaLx3aA9gubuvcPckcDfQz91fcvfDG32t4qtR/Uf115beeiYRkfWImvO7A7fnwd/tx6YbxDlx2j8ZPfuVUC3fzGUOvz2wssFxbX1bUx4AeprZDcCTTZ1kZkPMbImZLVm9Ony/IUVE1meX9m2Ze+bP+M2+nZj+9P/S+/qneKk2HO8pKtpDW3df6+6D3f1Md5+8nvOmApcCzyUSiWKVJyKSN23iUUb3/RF3DO7B53Vpjrzx70x+fDmZbLAPdHMJ/LeADg2Oq+vbREQE2H+HrVj4+5/Tc5dtuHLhMo75wyL++8HawOrJJfAXAzuYWWczSwADgNn5KErvtBWRcrHZhgkmDdyTa47ZnWXvfkav655k5uKVgSzfbO6yzBnAIqCrmdWa2WB3TwNDgYXAq8BMd38lH0VpLx0RKSdmxpF7VrPg9/uza3VbLrj/RU6941k+8E2KWkesOSe5+8Am2ucD8/NakYhImarefEPu+u1PueVvK7hq4ev0zEzkyvhUfhF9oSj3D+XWCprSEZFyFYkYQ37ehVlD96OdfcrJqQsYmRrEWq8q/L0LfgcREfmOnbfdlIcSF3NKdC53ZQ7i8OQ4Xsh2Keg9Qxn4msMXkUrQxlKMjN/FnfHL+dIT9E+O5rr0kaQz2YLcL5SBrykdEakk+0aXsqBqOH0ii7gm/WuOumkRb36wJu/3CWXga4QvIpWmra3l2sSN3BC/nrc//oIvU/kf5Ycy8DXCF5FK1Sf6DE8N+wVdt8n/ks1QBr6ISCWrikUL8rqhDHxN6YiI5F8oA19TOiIi+RfKwBcRkfxT4IuIVAgFvohIhQhl4OuhrYhI/oUy8PXQVkQk/yyITfiby8xWA2/WH7YFGg/5G7c1PG4HvF+g0tZVS76uWd95TX2vOX2zrrYw91dzr8tXf62rXf21/u9Ven+t7/tB99cP3X2r77S6e0l8AVO/r63hMbCkmLXk65r1ndfU95rTN6XWX829Ll/99X39U8n91dT3Kr2/1vf9sPZXKKd0mjCnGW3rOqcQWnOf5l6zvvOa+l5z+mZdbWHur+Zel6/+Wle7+mv936v0/lrf90PZX6Ge0smFmS1x9+5B11Eq1F8to/5qGfVXyxSqv0pphN9SU4MuoMSov1pG/dUy6q+WKUh/le0IX0REvq2cR/giItKAAl9EpEIo8EVEKkRFBr6ZHWFmN5vZPWZ2aND1hJ2ZbWdmt5rZfUHXElZmtpGZ3Vb/c3Vc0PWEnX6mWiZfmVVygW9m08xslZm93Ki9xsyWmdlyMxu+vtdw94fc/RTgNOCYQtYbtDz11wp3H1zYSsOnhX33K+C++p+rvkUvNgRa0l+V+jPVUAv7Ky+ZVXKBD0wHaho2mFkUmAz0AroBA82sm5ntamZzG31t3eDSi+qvK2fTyV9/VZrpNLPvgGpgZf1pmSLWGCbTaX5/Sev6K6fMirX2wqC4+5Nm1qlRcw9gubuvADCzu4F+7j4eOLzxa5iZAROABe7+XIFLDlQ++qtStaTvgFq+Cv0XKM2BVM5a2F9Li1xe6LSkv8zsVfKQWeXyg9meb0ZX8NV/vvbrOf9M4JfAUWZ2WiELC6kW9ZeZbWlmNwF7mtmIQhcXck313QNAfzObQvG2FCgF6+wv/Uw1qamfr7xkVsmN8PPB3a8Hrg+6jlLh7h/w1dyhNMHd1wAnB11HqdDPVMvkK7PKZYT/FtChwXF1fZusm/qr9dR3LaP+apmC9le5BP5iYAcz62xmCWAAMDvgmsJM/dV66ruWUX+1TEH7q+QC38xmAIuArmZWa2aD3T0NDAUWAq8CM939lSDrDAv1V+up71pG/dUyQfSXNk8TEakQJTfCFxGR1lHgi4hUCAW+iEiFUOCLiFQIBb6ISIVQ4IuIVAgFvohIhVDgi4hUCAW+iEiF+H/ZzyFCVbf1EAAAAABJRU5ErkJggg==\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