random.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "IPython.OutputArea.prototype._should_scroll = function(lines) {\n",
       "    return false;\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%%javascript\n",
    "IPython.OutputArea.prototype._should_scroll = function(lines) {\n",
    "    return false;\n",
    "}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib\n",
    "import matplotlib.pyplot as plt\n",
    "import glob\n",
    "import math\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "from bisect import *\n",
    "from decimal import Decimal\n",
    "from scipy.optimize import curve_fit\n",
    "import scipy.special as sps"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD+CAYAAADBCEVaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAQfUlEQVR4nO3df4hddXrH8ffTrG4l02ZLI1NJ0o5LRJqa/jBDtAhl0rIwrsYs21CStXazqMNukbZUaNNSurRQ6h9auv4oElSybV1nxRZNTMQurEEWrNVYtzGb7pKKiwmL2VU627G2ku3TP+YmXmbnZs6dufeec7/zfsHAveeeueeZk5zPfOc533tOZCaSpLL8SN0FSJJ6z3CXpAIZ7pJUIMNdkgpkuEtSgQx3SSqQ4S5JBTLcJalAPQ/3iPjZiHgwIp6IiM/1+v0lSYurFO4R8UhEnImI1+Ytn4yIb0bEyYjYC5CZJzLzs8BvANf1vmRJ0mKqjtz3A5PtCyJiFfAAcD2wCdgdEZtar90EHAIO96xSSVJlH6qyUmY+HxFj8xZvBU5m5usAETEN7AC+kZkHgAMRcQj40mLvv3bt2hwbm//21bz77rusXr16Sd/bT9bVHevqTlPrgubWVmJdR48e/V5mXrrQa5XCvYN1wJttz08B10TEBPBJ4MNcYOQeEVPAFMDo6Ch33333koqYnZ1lZGRkSd/bT9bVHevqTlPrgubWVmJd27Zt+3bHFzOz0hcwBrzW9nwn8FDb81uA+6u+X/vXli1bcqmee+65JX9vP1lXd6yrO02tK7O5tZVYF/BydsjV5cyWOQ1saHu+vrVMklSz5YT7S8AVEXF5RFwM7AIOdPMGEbE9IvbNzMwsowxJ0nxVp0I+BrwAXBkRpyLi1sw8C9wBPAucAB7PzOPdbDwzD2bm1Jo1a7qtW5J0AVVny+zusPwwTneUpMbx8gOSVKBaw92euyT1R63hbs9dkvpjOR9iapyxvYfOP37jrhtqrESS6jX04X7s9Ax72kJdkmTPXZKKZM9dkgrkVEhJKtDQ99w7GevQh/dEq6SVwJ67JBXInrskFcieuyQVyHCXpAIVe0K1k04nWsGTrZLK4chdkgpU68g9IrYD2zdu3FhnGed5bRpJpXC2jCQVyLaMJBVoxZ1QrcoWjaRhZrhXYNBLGja2ZSSpQIa7JBXIC4dJUoFq7bln5kHg4Pj4+O111tEN+++ShoFtGUkqkOEuSQUy3CWpQM5zXwb775KaypG7JBXIkXuPOIqX1CTOc5ekAjnPfYAc3UsaFHvufTC29xDHTs9c8JZ+ktRPhrskFchwl6QCGe6SVCDDXZIKZLhLUoEMd0kqkOEuSQXy8gN91mmuux9oktRPhnsDGPSSes1ry0hSgby2TMM4ipfUC55QlaQCGe6SVCBPqDZYpxaNrRtJi3HkLkkFMtwlqUC2ZYaEN/6Q1A3DfcjZf5e0ENsyklQgw12SCmS4S1KB7LkXxP67pHMcuUtSgRy5rwDtI/r9k6trrETSoDhyl6QCGe6SVKC+tGUi4hPADcCPAw9n5j/1YzuSpIVVDveIeAS4ETiTmVe1LZ8EvgCsAh7KzLsy80ngyYj4CeBuwHAfMC9XIK1s3bRl9gOT7QsiYhXwAHA9sAnYHRGb2lb5k9brkqQBqjxyz8znI2Js3uKtwMnMfB0gIqaBHRFxArgLeCYzX+lRreox58VL5YrMrL7yXLg/fa4tExE7gcnMvK31/BbgGuBbwKeBl4BXM/PBBd5rCpgCGB0d3TI9Pb2kH+DMOzO89d6SvrWvRi+hkXVdvmYVIyMjABw7/cGNyTevW1NXSQDMzs6er6tJrKt7Ta2txLq2bdt2NDPHF3qtLydUM/Ne4N5F1tkH7AMYHx/PiYmJJW3rvkef4p5jzZuuf+fms82t62vvtp59UN8bN0/UUs85R44cYan/B/rJurrX1NpWWl3LTZ/TwIa25+tbyzTE5p+MtWUjDZ/lhvtLwBURcTlzob4L+FTVb46I7cD2jRs3LrMMLZeza6SyVJ4tExGPAS8AV0bEqYi4NTPPAncAzwIngMcz83jV98zMg5k5tWZNvf1eSSpNN7NldndYfhg43LOKJEnL5uUHJKlAtYZ7RGyPiH0zMzOLryxJqqzWuXqZeRA4OD4+fnuddejC/LCTNHyaNxFbjWbQS8PBnrskFajWkbvz3MvhiF5qllpH7s5zl6T+sOeuJfNTrVJz2XOXpAI5z12SCmTPXZIKZFtGkgrkCVX1nNMipfo5cpekAhnuklQgZ8tIUoG8KqQGxl68NDi2ZSSpQIa7JBXIqZBqLNs40tIZ7uorLy4m1cPruasWjsql/vLaMpJUINsyGgqO9KXuOFtGkgpkuEtSgQx31W5s7yGOnZ5xZo3UQ4a7JBXIE6oaOp5clRbnyF2SCuQlfyWpQH6ISZIKZM9djeKMGak37LlLUoEcuWuoOXNGWpgjd0kqkOEuSQUy3CWpQIa7JBXIcJekAhnuklQgw12SCuQNslU858JrJfLaMpJUID+hqiJ5jRqtdIa7VqxOvwBs3agEnlCVpAIZ7tIFnLtxt20eDRvDXZIKZLhLUoEMd0kqkLNlpHna++t3bl54OTirRs3myF2SCuTIXcWoMqOlX7NevMSBmsaRuyQVyHCXpAIZ7pJUIMNdkgrkCVWpQTwxq17p+cg9Ij4aEQ9HxBO9fm9JUjWVRu4R8QhwI3AmM69qWz4JfAFYBTyUmXdl5uvArYa7StfttEpH5RqkqiP3/cBk+4KIWAU8AFwPbAJ2R8SmnlYnSVqSSuGemc8D78xbvBU4mZmvZ+b7wDSwo8f1SZKWIDKz2ooRY8DT59oyEbETmMzM21rPbwGuAT4P/AXwMeZaNX/Z4f2mgCmA0dHRLdPT00v6Ac68M8Nb7y3pW/tq9BKsqwsl1bV53cL3BD52eqYn6wDMzs4yMjLSXWED0tTaSqxr27ZtRzNzfKHXej5bJjPfBj5bYb19wD6A8fHxnJiYWNL27nv0Ke451rxJP3duPmtdXSiprjdunlhw+Z72nvsy1gE4cuQISz1m+q2pta20upYzW+Y0sKHt+frWMklSzZYzVHoJuCIiLmcu1HcBn+rmDSJiO7B948aNyyhDaq5+XKhsbO8h7tx8lj17DznrRh1VGrlHxGPAC8CVEXEqIm7NzLPAHcCzwAng8cw83s3GM/NgZk6tWdO5tyhJ6l6lkXtm7u6w/DBwuKcVSZKWzWvLSFKBag33iNgeEftmZmYWX1mSVFmt4W7PXZL6w7aMJBXIcJekAtX6kUDnuatE/ZrbLnXDnrskFci2jCQVyHCXpALZc5eGmHd3Uif23CWpQLZlJKlAhrskFchwl6QCGe6SVCBny0g1qPKJ035/KrXbmTbOzBkuzpaRpALZlpGkAhnuklQgw12SCmS4S1KBnC0jFcLZLEtX4r5ztowkFci2jCQVyHCXpAIZ7pJUIMNdkgpkuEtSgQx3SSqQ4S5JBfJDTNIK0ukywp0+xNPvyw6rf/wQkyQVyLaMJBXIcJekAhnuklQgw12SCmS4S1KBDHdJKpDhLkkFMtwlqUCGuyQVyHCXpAJ5bRmpcN1eH2a56++fXL3oep2uX1Pl5tTd1lfXDa/n1znoOry2jCQVyLaMJBXIcJekAhnuklQgw12SCmS4S1KBDHdJKpDhLkkFMtwlqUCGuyQVyHCXpAIZ7pJUIMNdkgpkuEtSgQx3SSqQ4S5JBer5zToiYjXwN8D7wJHMfLTX25AkXVilkXtEPBIRZyLitXnLJyPimxFxMiL2thZ/EngiM28HbupxvZKkCqq2ZfYDk+0LImIV8ABwPbAJ2B0Rm4D1wJut1X7QmzIlSd2oFO6Z+TzwzrzFW4GTmfl6Zr4PTAM7gFPMBXzl95ck9VZkZrUVI8aApzPzqtbzncBkZt7Wen4LcA3wh8D9wP8AX+vUc4+IKWAKYHR0dMv09PSSfoAz78zw1ntL+ta+Gr0E6+qCdXVnsbo2r/vgvsTHTs/0fPsXev/22rqto9P6vfh55u+zKtvqpFMNF6qz0/vOzs4yMjKy6DYXsm3btqOZOb7Qaz0/oZqZ7wKfqbDePmAfwPj4eE5MTCxpe/c9+hT3HOv5j7Fsd24+a11dsK7uLFbXGzdPnH+8Z++hnm//Qu/fXlu3dXRavxc/z/x9VmVbnXSq4UJ1dnrfI0eOsNT8u5DltE1OAxvanq9vLZMk1Ww54f4ScEVEXB4RFwO7gAPdvEFEbI+IfTMzvf+zUZJWsqpTIR8DXgCujIhTEXFrZp4F7gCeBU4Aj2fm8W42npkHM3NqzZrFe1ySpOoqNRMzc3eH5YeBwz2tSJK0bE5VlKQC1Rru9twlqT9qDXd77pLUH7ZlJKlAlT+h2tciIr4LfHuJ374W+F4Py+kV6+qOdXWnqXVBc2srsa6fycxLF3qhEeG+HBHxcqeP39bJurpjXd1pal3Q3NpWWl22ZSSpQIa7JBWohHDfV3cBHVhXd6yrO02tC5pb24qqa+h77pKkH1bCyF2SNM/QhHuH+7W2v/7hiPhy6/UXWzcXaUJdeyLiuxHxauvrtgHVteB9b9tej4i4t1X3v0XE1Q2payIiZtr2158OoKYNEfFcRHwjIo5HxO8usM7A91fFuurYXz8aEf8SEV9v1fVnC6wz8OOxYl21HI+tba+KiH+NiKcXeK33+yszG/8FrAL+A/gocDHwdWDTvHV+G3iw9XgX8OWG1LUHuL+GffYrwNXAax1e/zjwDBDAtcCLDalrgrk7fg1yX10GXN16/GPAtxb4dxz4/qpYVx37K4CR1uOLgBeBa+etU8fxWKWuWo7H1rZ/H/jSQv9e/dhfwzJy73S/1nY7gC+2Hj8B/FpERAPqqkUufN/bdjuAv805/wx8JCIua0BdA5eZ38nMV1qP/4u5S1ivm7fawPdXxboGrrUPZltPL2p9zT95N/DjsWJdtYiI9cANwEMdVun5/hqWcF8HvNn2/BQ//J/8/Do5d635GeAnG1AXwK+3/pR/IiI2LPB6HarWXodfbv1p/UxE/NwgN9z6c/iXmBv1tat1f12gLqhhf7VaDK8CZ4CvZGbH/TXA47FKXVDP8fjXwB8A/9fh9Z7vr2EJ92F2EBjLzJ8HvsIHv521sFeY+0j1LwD3AU8OasMRMQL8A/B7mfn9QW13MYvUVcv+yswfZOYvMnd7za0RcdUgtruYCnUN/HiMiBuBM5l5tN/bajcs4V7lfq3n14mIDwFrgLfrrisz387M/209fQjY0ueaqmrkPXAz8/vn/rTOuZvBXBQRa/u93Yi4iLkAfTQz/3GBVWrZX4vVVdf+atv+fwLPAZPzXqrjeFy0rpqOx+uAmyLiDeZat78aEX8/b52e769hCfcq92s9AHy69Xgn8NVsnZ2os655fdmbmOubNsEB4Ldas0CuBWYy8zt1FxURP3Wu1xgRW5n7P9rXUGht72HgRGb+VYfVBr6/qtRV0/66NCI+0np8CfAx4N/nrTbw47FKXXUcj5n5R5m5PjPHmMuIr2bmb85bref7q9Jt9uqWmWcj4tz9WlcBj2Tm8Yj4c+DlzDzA3EHwdxFxkrkTdrsaUtfvRMRNwNlWXXv6XRecv+/tBLA2Ik4Bn2fuBBOZ+SBzt0f8OHAS+G/gMw2payfwuYg4C7wH7BrAL+nrgFuAY61+LcAfAz/dVlcd+6tKXXXsr8uAL0bEKuZ+mTyemU/XfTxWrKuW43Eh/d5ffkJVkgo0LG0ZSVIXDHdJKpDhLkkFMtwlqUCGuyQVyHCXpAIZ7pJUIMNdkgr0/7skv46ofG78AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "shape, scale = 5., 10.  # mean=4, std=2*sqrt(2)\n",
    "s = np.random.exponential(.5, 10000)\n",
    "count, bins, ignored = plt.hist(s, 100, log=True,range=[0,4])\n",
    "plt.grid()\n",
    "plt.show()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAD6CAYAAACrklzBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAARyElEQVR4nO3df6xed13A8ffHwga0eB0Wr2SbtEuXYVk10JtNkei9/gjdsBRxMa1zodhRh85obIwzM2pMDJhIVGBmacYySZZe5lRYt5KBsmZRtrGVjHVzDkqpsoassOLFO5FZ8vGP53ScXe/Tnvvc5zzP6bfvV3Kz5/meX59+79mn337O95wTmYkkqTzfM+4AJEntMMFLUqFM8JJUKBO8JBXKBC9JhTLBS1KhTPCSVKiXjDsAgNWrV+eaNWsG2va5555j5cqVww1oSIxtMF2NratxgbENqquxNY3rwIEDX8/MV/ddITOH+gP8MHAzcCfwnibbbNy4MQd13333Dbxt24xtMF2NratxZRrboLoaW9O4gEfyFLm1UYkmIm6NiGMR8fiC9k0R8VREHIqIG6q/MJ7MzOuAXwJ+osn+JUnD17QGfxuwqd4QESuAm4ArgPXAtohYXy17G3APsG9okUqSlqRRgs/M+4HjC5ovAw5l5uHMfB6YBbZU69+VmVcAVw8zWElSc5ENHzYWEWuAuzPz0ur7VcCmzLy2+n4NcDm92vs7gHOBxzLzpj772wnsBJicnNw4Ozs70B9gfn6eVatWDbRt24xtMF2NratxgbENqquxNY1rZmbmQGZO9V3hVAX6fPHF0zXA47XvVwG31L5fA3yo6f7qP15kHT1jW7quxpVpbIPqamwjvcjax1Hgwtr3C6o2SVIHLCfBPwxcHBFrI+IcYCtw11J2EBGbI2L33NzcMsKQJC2m6TTJPcADwCUR8XRE7MjME8D1wL3Ak8AdmfnEUg6emXszc+fExMRS45YknUajO1kzc1uf9n2MeSrkwaNzbL/hHgCOvO+t4wxFkjplrM+isUQjSe0Z67NoMnMvsHdqaurdw9jfmmokD47mJcmnSUpSoTrxNMk2OJqXdLazBi9JhRprgneapCS1xxq8JBXKBC9JhSr2ImudF1wlnY3GmuAjYjOwed26dSM7psle0tnCi6ySVChr8JJUqLOiBt+U5RtJJXEEL0mFOusustbVR+ySVBovskpSoSzRSFKhTPCSVChn0fTRrz7v7BpJZwoT/BI5lVLSmcISjSQVyhd+SFKhnCYpSYWyBr8M1uMldZk1eEkqlAlekgplgpekQpngJalQXmQdEi+4Suoa58FLUqGcBy9JhbIGL0mFsgbfgpP1+F0bTrDd2rykMXEEL0mFMsFLUqEs0YyQUykljZIjeEkqlAlekgpliWZMfOerpLY5gpekQpngJalQPotGkgrls2gkqVCWaCSpUCZ4SSqU0yQ7xrtdJQ2LI3hJKpQj+A5zNC9pORzBS1KhHMGfIRzNS1oqR/CSVCgTvCQVyhLNGchyjaQmHMFLUqEcwZ/hfK68pH4cwUtSoUzwklQoE7wkFaqVGnxEvB14K/C9wIcz85NtHEfNOOtGOjs1TvARcSvw88CxzLy01r4J+CtgBXBLZr4vMz8GfCwizgP+HDDBj1i/i6+Szh5LKdHcBmyqN0TECuAm4ApgPbAtItbXVvmDarkkacQaJ/jMvB84vqD5MuBQZh7OzOeBWWBL9PwZ8InM/NzwwpUkNRWZ2XzliDXA3SdLNBFxFbApM6+tvl8DXA58AXgn8DDwaGbevMi+dgI7ASYnJzfOzs4O9Ac4dnyOZ7410Katm3w5nYttw/m999/Oz8+zatWqMUezuK7G1tW4wNgG1dXYmsY1MzNzIDOn+i1v5SJrZn4A+MBp1tkN7AaYmprK6enpgY71wds/zvsPdvN+rV0bTnQutiNXTwOwf/9+Bu3ztnU1tq7GBcY2qK7GNqy4ljtN8ihwYe37BVWbJGnMlpvgHwYujoi1EXEOsBW4q+nGEbE5InbPzc0tMwxJ0kKNE3xE7AEeAC6JiKcjYkdmngCuB+4FngTuyMwnmu4zM/dm5s6JiYmlxi1JOo3GBeLM3NanfR+wb2gRSZKGYqxXACNiM7B53bp14wzjrNXvDteFN0l596t0Zhrrs2gs0UhSe7o1h09j46MNpPL4NElJKtRYE7zTJCWpPWMt0WTmXmDv1NTUu8cZh07Nxw1LZyZr8GeZk8l614YT+OuXymYNXpIKZYKXpEJ5kVWSCuWNTpJUKEs0klQop1FoSfrd8er0Sal7HMFLUqG8yCpJhfJOVg2Fd7tK3WOJRpIKZYKXpEI5i0YjYxlHGi1H8JJUKGfRSFKhnEWjoWtSirFcI7XPEo0kFcoEL0mFMsFLUqFM8JJUKOfBq1X9nj4pqX2O4CWpUCZ4SSqUNzpJUqG80UlnBG+MkpbOEo0kFcoEL0mFcpqkxs7yi9QOR/CSVCgTvCQVygQvSYWyBq9OWXPDPezacILtDR9xYP1e6s8RvCQVygQvSYUywUtSoXwWjSQVymfRqEhefJUs0UhSsUzwklQoE7wkFcoEL0mFMsFLUqFM8JJUKJ9Fo2Ksafj8Guls4QhekgrlCF7F86Ynna0cwUtSoUzwklQoSzQ647RxMbVfGafeftumlUM/rtQmR/CSVCgTvCQVaugJPiIuiogPR8Sdw963JKm5Rgk+Im6NiGMR8fiC9k0R8VREHIqIGwAy83Bm7mgjWGmcDh6dY80N97zwI3Vd0xH8bcCmekNErABuAq4A1gPbImL9UKOTJA2sUYLPzPuB4wuaLwMOVSP254FZYMuQ45MkDSgys9mKEWuAuzPz0ur7VcCmzLy2+n4NcDnwR8CfAj8H3JKZ7+2zv53AToDJycmNs7OzA/0Bjh2f45lvDbRp6yZfjrENoM3YNpw/8cLng0cXfxdwv3VOFVd9m3GYn59n1apVY42hH2NbuqZxzczMHMjMqX7Lhz4PPjOfBa5rsN5uYDfA1NRUTk9PD3S8D97+cd5/sJvT+XdtOGFsA2gztiNXT7/weXufOnq/dU4VV32bcdi/fz+D/j/UNmNbumHFtZxZNEeBC2vfL6jaJEkdsJxh0sPAxRGxll5i3wr88lJ2EBGbgc3r1q1bRhhSc6OY/dLvGEt90JkPSdNyNZ0muQd4ALgkIp6OiB2ZeQK4HrgXeBK4IzOfWMrBM3NvZu6cmBhv/VKSStRoBJ+Z2/q07wP2DTUiSdJQ+KgCSSrUWBN8RGyOiN1zc4tPV5MkDW6sCd4avCS1xxKNJBXKBC9JhRrrrYzOg1cXtTVX3nntGjVr8JJUKEs0klQoE7wkFcoavDQEbT/jxvq9BmENXpIKZYlGkgplgpekQpngJalQJnhJKpSzaKQxGMWbpSRn0UhSoSzRSFKhTPCSVCgTvCQVygQvSYVyFo10hllsBs6uDSeYHn0o6jhn0UhSoSzRSFKhTPCSVCgTvCQVygQvSYUywUtSoUzwklQoE7wkFcobnaTC9Xs0cb+Xdzd5lLEv/j4zeKOTJBXKEo0kFcoEL0mFMsFLUqFM8JJUKBO8JBXKBC9JhTLBS1KhTPCSVCgTvCQVygQvSYXyWTRSgZo8T6a+zlKfLXOq/TfZ17COPcptz0Q+i0aSCmWJRpIKZYKXpEKZ4CWpUCZ4SSqUCV6SCmWCl6RCmeAlqVAmeEkqlAlekgplgpekQpngJalQJnhJKpQJXpIKZYKXpEIN/XnwEbES+GvgeWB/Zt4+7GNIkk6v0Qg+Im6NiGMR8fiC9k0R8VREHIqIG6rmdwB3Zua7gbcNOV5JUkNNSzS3AZvqDRGxArgJuAJYD2yLiPXABcBXqtW+M5wwJUlL1SjBZ+b9wPEFzZcBhzLzcGY+D8wCW4Cn6SX5xvuXJA1fZGazFSPWAHdn5qXV96uATZl5bfX9GuBy4PeADwH/A/xzvxp8ROwEdgJMTk5unJ2dHegPcOz4HM98a6BNWzf5coxtAF2NratxQXdj23D+BPPz86xatYqDR+f6rrOY+vr1dZrsp8m29diaWGr8Tfe12PZN45qZmTmQmVP9lg/9ImtmPge8q8F6u4HdAFNTUzk9PT3Q8T54+8d5/8Gxvju8r10bThjbALoaW1fjgu7GduTqafbv38/09DTb+7yo+8jV04u219evr9NkP022rcfWxFLjb7qvxbZfSlynspwSylHgwtr3C6o2SVIHLCfBPwxcHBFrI+IcYCtw11J2EBGbI2L33Nzi//SRJA2u6TTJPcADwCUR8XRE7MjME8D1wL3Ak8AdmfnEUg6emXszc+fExNJrWJKkU2tUtMvMbX3a9wH7hhqRJGkoxjqN0RKNJLVnrAneEo0ktccbkSSpUCZ4SSpU4ztZWw0i4mvAvw+4+Wrg60MMZ5iMbTBdja2rcYGxDaqrsTWN67WZ+ep+CzuR4JcjIh451a2642Rsg+lqbF2NC4xtUF2NbVhxWaKRpEKZ4CWpUCUk+N3jDuAUjG0wXY2tq3GBsQ2qq7ENJa4zvgYvSVpcCSN4SdIiOp3g+7zztb783Ij4aLX8oeqlJCeX/X7V/lREvGUMsf1ORPxrRDwWEf8UEa+tLftORDxa/SzpCZxDiGt7RHytdvxra8veGRFfrH7eOcy4Gsb2F7W4vhAR/1lb1mafLfrO4dryiIgPVHE/FhFvrC1ru89OF9vVVUwHI+IzEfGjtWVHqvZHI+KRMcQ2HRFztd/bH9aWnfJcGEFsv1uL6/Hq/HpVtay1fouICyPivio3PBERv7XIOsM73zKzkz/ACuBLwEXAOcDngfUL1vl14Obq81bgo9Xn9dX65wJrq/2sGHFsM8Arqs/vORlb9X1+jH22HfjQItu+Cjhc/fe86vN5o4xtwfq/Cdzadp9V+/5J4I3A432WXwl8Agjgx4CHRtFnDWN708lj0ns/8kO1ZUeA1WPst2l6b4Fb1rnQRmwL1t0MfHoU/Qa8Bnhj9fmVwBcW+X90aOdbl0fw/d75WrcF+Jvq853Az0REVO2zmfntzPwycKja38hiy8z7MvO/q68P8t331LapSZ/18xbgU5l5PDO/AXyKBS9aH3Fs24A9Qzx+X7n4O4frtgAfyZ4Hge+LiNfQfp+dNrbM/Ex1bBjdeXby2Kfrt36Wc562Edsoz7WvZubnqs//Re9R6+cvWG1o51uXE/z5wFdq35/m/3fEC+tk7/n0c8D3N9y27djqdtD7G/mkl0XEIxHxYES8fQxx/WL1T787I+LkW7k602dVOWst8Olac1t91kS/2Nvus6VaeJ4l8MmIOBC9dyCPw49HxOcj4hMR8fqqrTP9FhGvoJck/67WPJJ+i15J+Q3AQwsWDe18695LHAsTEb8CTAE/VWt+bWYejYiLgE9HxMHM/NKIQtoL7MnMb0fEr9H7F9BPj+jYTW0F7szM79TaxtlnnRcRM/QS/JtrzW+u+uwHgE9FxL9VI9tR+Ry939t8RFwJfAy4eITHb2Iz8C+ZWR/tt95vEbGK3l8qv52Z3xzmvuu6PIJv8s7XF9aJiJcAE8CzDbdtOzYi4meBG4G3Zea3T7Zn5tHqv4eB/fT+Fh9JXJn5bC2WW4CNTbdtO7aarSz4J3OLfdZEv9g78V7iiPgRer/LLZn57Mn2Wp8dA/6B4ZYpTyszv5mZ89XnfcBLI2I1Hem3yqnOtVb6LSJeSi+5356Zf7/IKsM739q4kDCkixEvoXcRYS3fvRDz+gXr/AYvvsh6R/X59bz4IuthhnuRtUlsb6B3IeniBe3nAedWn1cDX2RIF5gaxvWa2udfAB7M717A+XIV33nV51eNss+q9V5H7yJXjKLPasdYQ/+LhW/lxRe9PjuKPmsY2w/Ru8b0pgXtK4FX1j5/Btg04th+8OTvkV6S/I+qDxudC23GVi2foFenXzmqfqv+/B8B/vIU6wztfBtqh7bwC7qS3lXmLwE3Vm1/Qm9EDPAy4G+rE/yzwEW1bW+stnsKuGIMsf0j8AzwaPVzV9X+JuBgdVIfBHaMOK73Ak9Ux78PeF1t21+t+vIQ8K5R91n1/Y+B9y3Yru0+2wN8FfhfenXNHcB1wHXV8gBuquI+CEyNsM9OF9stwDdq59kjVftFVX99vvp93ziG2K6vnWsPUvtLaLFzYZSxVetspzcZo75dq/1Gr4SWwGO139mVbZ1v3skqSYXqcg1ekrQMJnhJKpQJXpIKZYKXpEKZ4CWpUCZ4SSqUCV6SCmWCl6RC/R8WbgXvXj3bZQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD7CAYAAACRxdTpAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAQNklEQVR4nO3df4xl5V3H8ffXRbDdiUClmejuxtlmNqsbiLbcAKbGzFSbDKXDNg0xbAiCWZjUSKyGxGI0Uf8w0qT+gHRtMwHc1jQ7rmuDu8sqKnYkJojL1oYFtrRb3Jbd0C6IXR1SQ9d+/eOewet0bvfM/cGd+9z3K9kw57nnnvM8nOHD3e99znMiM5EkleX7Bt0BSVLvGe6SVCDDXZIKZLhLUoEMd0kqkOEuSQUy3CWpQIa7JBXoon4cNCI2Av8I/E5mHr7Q/ldccUVOTEx0dK7XXnuNjRs3dvTeYeWYR4NjHg3djPnYsWOvZObbV3utVrhHxEPA+4GzmXllS/sMcB+wAXggM++tXvoIsL9uBycmJnjqqafq7v7/LC4uMjU11dF7h5VjHg2OeTR0M+aI+Gq71+qWZfYCMysOugHYA1wP7AB2RcSOiHgv8BxwtqPeSpK6VuuTe2Y+HhETK5qvAU5m5gsAEbEA7ATGgI00A/9bEXEkM7/Tsx5Lki4o6i4cVoX74eWyTETcBMxk5h3V9q3AtZl5V7V9O/BKu5p7RMwBcwDj4+NXLywsdDSApaUlxsbGOnrvsHLMo8Exj4Zuxjw9PX0sMxurvdaXL1QBMnPvBV6fB+YBGo1GdlpzskY3GhzzaHDMvdPNVMgzwJaW7c1VW20RMRsR8+fOneuiG5KklboJ96PAtojYGhEXAzcDB9dygMw8lJlzl156aRfdkCStVCvcI2If8ASwPSJOR8TuzDwP3AU8CpwA9mfms/3rqiSprrqzZXa1aT8CHOn05BExC8xOTk52eghJ0ir69oVqHZl5CDjUaDTu7MXxJu555I2fT917Qy8OKUlDybVlJKlAAw13Z8tIUn8MNNydLSNJ/WFZRpIKZFlGkgpkWUaSCmRZRpIKZLhLUoGsuUtSgay5S1KBLMtIUoEMd0kqkOEuSQXyC1VJKpBfqEpSgSzLSFKBDHdJKpDhLkkFMtwlqUDOlpGkAjlbRpIKZFlGkgpkuEtSgQx3SSqQ4S5JBTLcJalAhrskFchwl6QCeROTJBXIm5gkqUCWZSSpQIa7JBXIcJekAhnuklQgw12SCmS4S1KBLhp0B/pl4p5H3vj51L03DLAnkvTm85O7JBXIcJekAhnuklSgnod7RPx4RHwyIg5ExC/1+viSpAurFe4R8VBEnI2IZ1a0z0TE8xFxMiLuAcjME5n5IeDngXf3vsuSpAup+8l9LzDT2hARG4A9wPXADmBXROyoXrsReAQ40rOeSpJqqxXumfk48OqK5muAk5n5Qma+DiwAO6v9D2bm9cAtveysJKmeyMx6O0ZMAIcz88pq+yZgJjPvqLZvBa4FDgAfBC4Bns7MPW2ONwfMAYyPj1+9sLDQ0QCWlpYYGxsD4PiZ1deFv2pTWUsKt455VDjm0eCY12Z6evpYZjZWe63nNzFl5iKwWGO/eWAeoNFo5NTUVEfnW1xcZPm9t7fcuNTq1C2dHXu9ah3zqHDMo8Ex9043s2XOAFtatjdXbbX5JCZJ6o9uwv0osC0itkbExcDNwMG1HMAnMUlSf9SdCrkPeALYHhGnI2J3Zp4H7gIeBU4A+zPz2f51VZJUV62ae2buatN+hC6mO0bELDA7OTnZ6SEkSavwAdmSVCDXlpGkAg003J0tI0n9YVlGkgpkWUaSCmRZRpIKZFlGkgpkWUaSCmS4S1KBrLlLUoF6vuTvWmTmIeBQo9G4s5/nmWhZCvjUvTf081SStC5YlpGkAhnuklQgw12SCuQXqpJUIG9ikqQCWZaRpAIZ7pJUIMNdkgpkuEtSgQx3SSqQUyElqUBOhZSkAlmWkaQCGe6SVCDDXZIKZLhLUoEG+rCOQfDBHZJGgZ/cJalAznOXpAI5z12SCmRZRpIKZLhLUoEMd0kqkOEuSQUy3CWpQIa7JBXIcJekAhnuklQgw12SCmS4S1KB+rIqZER8ALgB+EHgwcz8236cR5K0utrhHhEPAe8HzmbmlS3tM8B9wAbggcy8NzMfBh6OiMuBjwHrMtxbl/8FlwCWVI61lGX2AjOtDRGxAdgDXA/sAHZFxI6WXX6rel2S9CaqHe6Z+Tjw6orma4CTmflCZr4OLAA7o+mjwF9n5ud7111JUh2RmfV3jpgADi+XZSLiJmAmM++otm8FrgW+BNwGHAW+kJmfXOVYc8AcwPj4+NULCwsdDWBpaYmxsTEAjp/pbl34qzYNx9LDrWMeFY55NDjmtZmenj6WmY3VXuvLF6qZeT9w/wX2mQfmARqNRk5NTXV0rsXFRZbfe/uKGvpanbqlsz682VrHPCoc82hwzL3T7VTIM8CWlu3NVVstPolJkvqj23A/CmyLiK0RcTFwM3Cw7pt9EpMk9UftcI+IfcATwPaIOB0RuzPzPHAX8ChwAtifmc/2p6uSpLpq19wzc1eb9iPAkU5OHhGzwOzk5GQnb5cktdGXL1TrysxDwKFGo3HnIPuxrPWmJm9okjTMXFtGkgo00HB3towk9cdAw93ZMpLUH5ZlJKlAlmUkqUCWZSSpQJZlJKlAhrskFciauyQVyJq7JBVooMsPrGcuRSBpmFlzl6QCGe6SVCC/UJWkAvmFqiQVyLKMJBXIcJekAhnuklQgw12SCmS4S1KBBnqHakTMArOTk5OD7MYFebeqpGHjVEhJKpBlGUkqkOEuSQUy3CWpQC75u0Z+uSppGPjJXZIK5Cf3LvgpXtJ65Sd3SSqQ67lLUoG8iUmSCmRZRpIKZLhLUoEMd0kqkOEuSQVynnuPOOdd0nriJ3dJKpDhLkkFsizTB3VKNJZxJPWTn9wlqUA9D/eIeEdEPBgRB3p9bElSPbXCPSIeioizEfHMivaZiHg+Ik5GxD0AmflCZu7uR2clSfXUrbnvBT4OfHq5ISI2AHuA9wKngaMRcTAzn+t1J4eZtXVJg1Drk3tmPg68uqL5GuBk9Un9dWAB2Nnj/kmSOhCZWW/HiAngcGZeWW3fBMxk5h3V9q3AtcBvA79H8xP9A5n5+22ONwfMAYyPj1+9sLDQ0QCWlpYYGxsD4PiZ9b108FWb/m/1y9a+trbX0TrmUeGYR4NjXpvp6eljmdlY7bWeT4XMzH8HPlRjv3lgHqDRaOTU1FRH51tcXGT5vbe3lEDWo1O3TL3xc2tfW9vraB3zqHDMo8Ex9043s2XOAFtatjdXbbX5sA5J6o9uwv0osC0itkbExcDNwMG1HMCHdUhSf9SdCrkPeALYHhGnI2J3Zp4H7gIeBU4A+zPz2f51VZJUV62ae2buatN+BDjS6ckjYhaYnZyc7PQQkqRV+AxVSSqQa8tIUoEGuirkqJVlJtb5VE1J5bAsI0kFsiwjSQUaaLh7E5Mk9YdlGUkqkGUZSSqQ4S5JBXIq5DrQbopku4d7HD9z7o1VJX0AiKTVWHOXpAJZlpGkAhnuklQgw12SCuRNTJJUIL9QlaQCWZaRpAIZ7pJUIMNdkgpkuEtSgQx3SSqQa8sMidb1Z+6+qvP3tuMaNVJZnAopSQWyLCNJBTLcJalAhrskFchwl6QCGe6SVCDDXZIKZLhLUoG8iUkX1HoT1Hq+2Wk99nM99kmjwZuYJKlAlmUkqUCGuyQVyHCXpAIZ7pJUIMNdkgpkuEtSgQx3SSqQ4S5JBTLcJalAPV9+ICI2An8CvA4sZuZnen0OSdL3VuuTe0Q8FBFnI+KZFe0zEfF8RJyMiHuq5g8CBzLzTuDGHvdXklRD3bLMXmCmtSEiNgB7gOuBHcCuiNgBbAZerHb7n950U5K0FrXCPTMfB15d0XwNcDIzX8jM14EFYCdwmmbA1z6+JKm3IjPr7RgxARzOzCur7ZuAmcy8o9q+FbgW+AjwceC/gX9qV3OPiDlgDmB8fPzqhYWFjgawtLTE2NgYAMfPnOvoGMNm/C3wjW8Nuhdw1abVV/Oscx3avbfdcbZeuuGN61xn/7Weq91x1vreXh6r9Xe7V9ban17+u6ijH2Neq2Ea8/T09LHMbKz2Ws+/UM3M14BfrLHfPDAP0Gg0cmpqqqPzLS4usvze21vWzi7Z3Ved5w+OD3QpfgBO3TK1anud69Duve2Os3dmIxf6HWl33jrnanectb63l8dq/d3ulbX2p5f/Lurox5jXqpQxd1M2OQNsadneXLXVFhGzETF/7txofOKWpDdLN+F+FNgWEVsj4mLgZuDgWg7gwzokqT/qToXcBzwBbI+I0xGxOzPPA3cBjwIngP2Z+Wz/uipJqqtW4TYzd7VpPwIc6fTkPkNVkvrDZ6hKUoGchy5JBRpouDtbRpL6w7KMJBWo9h2qfe1ExMvAVzt8+xXAKz3szjBwzKPBMY+Gbsb8o5n59tVeWBfh3o2IeKrd7belcsyjwTGPhn6N2S9UJalAhrskFaiEcJ8fdAcGwDGPBsc8Gvoy5qGvuUuSvlsJn9wlSSsMdbi3eYZrUSJiS0R8LiKei4hnI+LDVfvbIuLvIuLL1T8vH3RfeykiNkTEv0bE4Wp7a0Q8WV3rP69WIi1GRFwWEQci4osRcSIifmoErvGvVb/Tz0TEvoj4gdKu82rPn253XaPp/mrsT0fEu7o599CG+/d4hmtpzgN3Z+YO4Drgl6tx3gM8lpnbgMeq7ZJ8mOZqo8s+CvxRZk4C/wHsHkiv+uc+4G8y88eAn6A59mKvcURsAn4FaFRPd9tAc9nw0q7zXlY8f5r21/V6YFv1Zw74RDcnHtpwp/0zXIuSmS9l5uern/+L5n/0m2iO9VPVbp8CPjCYHvZeRGwGbgAeqLYDeA9woNqltPFeCvwM8CBAZr6emd+k4GtcuQh4S0RcBLwVeInCrnOb50+3u647gU9n0z8Dl0XED3d67mEO903Aiy3bp6u2YlXPsX0n8CQwnpkvVS99HRgfULf64Y+BXwe+U23/EPDN6hkCUN613gq8DPxpVYp6ICI2UvA1zswzwMeAr9EM9XPAMcq+zsvaXdeeZtowh/tIiYgx4C+BX83M/2x9LZtTnoqY9hQR7wfOZuaxQfflTXQR8C7gE5n5TuA1VpRgSrrGAFWdeSfN/7H9CLCR7y5fFK+f13WYw73rZ7gOi4j4fprB/pnM/GzV/I3lv7JV/zw7qP712LuBGyPiFM1S23to1qMvq/76DuVd69PA6cx8sto+QDPsS73GAD8H/FtmvpyZ3wY+S/Pal3ydl7W7rj3NtGEO966f4ToMqnrzg8CJzPzDlpcOArdVP98G/NWb3bd+yMzfyMzNmTlB85r+Q2beAnwOuKnarZjxAmTm14EXI2J71fSzwHMUeo0rXwOui4i3Vr/jy2Mu9jq3aHddDwK/UM2auQ4411K+WbvMHNo/wPuALwFfAX5z0P3p0xh/muZf254GvlD9eR/NOvRjwJeBvwfeNui+9mHsU8Dh6ud3AP8CnAT+Arhk0P3r8Vh/Eniqus4PA5eXfo2B3wW+CDwD/BlwSWnXGdhH8zuFb9P8G9rudtcVCJozAL8CHKc5k6jjc3uHqiQVaJjLMpKkNgx3SSqQ4S5JBTLcJalAhrskFchwl6QCGe6SVCDDXZIK9L8lzfjWrpcYXwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "s = np.power(np.random.uniform(0.,1.,10000),-0.5)\n",
    "plt.hist(np.log10(s), 100,log=True)\n",
    "plt.grid()\n",
    "plt.show()\n",
    "plt.hist(s, 100,log=True,range=[0,100])\n",
    "plt.grid()\n",
    "plt.show()"
   ]
  }
 ],
 "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.8.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}