1_Gaussian_Noise.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 90,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import random"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 91,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[0.64681121, 0.08899201, 0.5958028 ],\n",
       "       [0.81147962, 0.97623058, 0.49107451],\n",
       "       [0.62989361, 0.42443954, 0.95286122]])"
      ]
     },
     "execution_count": 91,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.random.random((3, 3))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 93,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[-0.89755835,  0.16213858,  0.90207152,  0.11642965,  0.76682339,\n",
       "         0.08232798,  0.50264859, -0.09631406, -0.94233744,  0.51212907],\n",
       "       [-1.1171383 , -0.23404022, -1.77253463,  0.54829055, -2.43962403,\n",
       "        -2.82092561, -0.55017031,  0.35646454, -0.19262477,  0.82110461],\n",
       "       [ 0.16323457, -1.35640292, -0.10409185,  0.18359682,  0.10844746,\n",
       "         0.12510907,  1.17450958, -0.93016321,  0.62160015, -0.33566053],\n",
       "       [ 0.29481277,  0.51218464,  0.79148149,  0.17905693,  1.21957955,\n",
       "        -0.59076701,  0.17927976, -1.97943024, -1.9442332 ,  0.47337426],\n",
       "       [ 0.35721529,  1.12738405,  2.06904748, -0.17774801, -0.59866904,\n",
       "         1.09447541,  0.07008963,  0.36743782,  1.95050528,  0.52485861],\n",
       "       [ 0.72246961,  0.29416237, -0.78847629,  1.91187326, -0.40041557,\n",
       "         0.83246766, -0.01893989,  2.11403677, -0.62141978, -1.36965489],\n",
       "       [ 0.86220807,  0.00560358,  0.71601784, -0.29949426, -0.08626726,\n",
       "        -0.15588242,  0.52149294, -0.74635169,  0.44768234, -0.36039751],\n",
       "       [ 2.23151098, -0.38171053, -0.29391018, -0.37563169, -0.05882703,\n",
       "         0.36766105,  0.36740065,  0.27381007, -0.5869817 ,  0.82409417],\n",
       "       [-0.6411837 ,  1.37754418,  0.18064629,  1.09125909, -2.43548823,\n",
       "        -0.07379421,  0.84011404,  0.63777825,  1.45174715, -1.46765984],\n",
       "       [-1.09805879,  1.56294336,  0.5259828 ,  0.72704649,  0.58942197,\n",
       "         1.55221329,  0.77333532, -0.72517139,  2.07619302, -0.65788277]])"
      ]
     },
     "execution_count": 93,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.random.normal(loc=0, scale=1, size=(10,10))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 128,
   "metadata": {},
   "outputs": [],
   "source": [
    "noise=np.random.normal(loc=0, scale=1, size=10000)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 129,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[  13  206 1363 3358 3419 1402  230    8]\n",
      "[-4 -3 -2 -1  0  1  2  3  4]\n"
     ]
    }
   ],
   "source": [
    "hist, bin_edges=np.histogram(noise, bins=[-4, -3, -2, -1, 0, 1, 2, 3, 4])\n",
    "print(hist)\n",
    "print(bin_edges)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 130,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAD8CAYAAAB+UHOxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAEwdJREFUeJzt3X2MXfWd3/H3Zw150CYtUIaUtR0Zbb3bkHTjoKmh4p+UZMGQKJBqkUBtYqVI3kogJVLaLiRS2SRFItpuqKJmqbzFDWmzoWiTCAu8JV5CFEUqDyZrHhyHMiVsmNjF3pqQRKhUkG//uD93LzAz9854PHemv/dLurrnfM/v3Ps9tud+5jzc41QVkqT+/MqkG5AkTYYBIEmdMgAkqVMGgCR1ygCQpE4ZAJLUKQNAkjplAEhSpwwASerUKZNuYCFnnnlmbdq0adJtSNKa8sgjj/xVVU2NGreqA2DTpk3s27dv0m1I0pqS5C/HGTfyEFCSNyV5KMmjSQ4k+UyrfznJj5Lsb48trZ4kX0wyk+SxJOcNvdb2JE+1x/albpwk6cSNswfwEnBRVf0iyanA95L8WVv2L6rqT18z/lJgc3ucD9wKnJ/kDOBGYBoo4JEku6vq+eXYEEnS4ozcA6iBX7TZU9tjoVuIXg58pa33AHBakrOBS4C9VXWsfejvBbadWPuSpKUa6yqgJOuS7AeOMPgQf7Atuqkd5rklyRtbbT3w7NDqs602X12SNAFjBUBVvVJVW4ANwNYk7wJuAP4u8PeBM4Dfa8Mz10ssUH+VJDuS7Euy7+jRo+O0J0lagkV9D6Cqfgp8B9hWVYfbYZ6XgP8IbG3DZoGNQ6ttAA4tUH/te+ysqumqmp6aGnkVkyRpica5CmgqyWlt+s3A+4EftuP6JAlwBfBEW2U38NF2NdAFwAtVdRi4F7g4yelJTgcubjVJ0gSMcxXQ2cDtSdYxCIw7q+ruJN9OMsXg0M5+4J+18XuAy4AZ4EXgYwBVdSzJ54CH27jPVtWx5dsUSdJiZDX/n8DT09PlF8EkaXGSPFJV06PGrepvAkvLZdP190y6hbE8c/MHJt2COuLN4CSpUwaAJHXKQ0A6IWvl0Iqk13MPQJI6ZQBIUqcMAEnqlAEgSZ0yACSpUwaAJHXKAJCkThkAktQpA0CSOmUASFKnDABJ6pQBIEmdMgAkqVMGgCR1ygCQpE4ZAJLUKQNAkjo1MgCSvCnJQ0keTXIgyWda/ZwkDyZ5Ksl/SfKGVn9jm59pyzcNvdYNrf5kkktO1kZJkkYbZw/gJeCiqno3sAXYluQC4PPALVW1GXgeuKaNvwZ4vqr+DnBLG0eSc4GrgHcC24A/SrJuOTdGkjS+kQFQA79os6e2RwEXAX/a6rcDV7Tpy9s8bfn7kqTV76iql6rqR8AMsHVZtkKStGhjnQNIsi7JfuAIsBf4H8BPq+rlNmQWWN+m1wPPArTlLwB/a7g+xzqSpBU2VgBU1StVtQXYwOC39nfMNaw9Z55l89VfJcmOJPuS7Dt69Og47UmSlmBRVwFV1U+B7wAXAKclOaUt2gAcatOzwEaAtvxvAseG63OsM/weO6tquqqmp6amFtOeJGkRxrkKaCrJaW36zcD7gYPA/cDvtGHbgbva9O42T1v+7aqqVr+qXSV0DrAZeGi5NkSStDinjB7C2cDt7YqdXwHurKq7k/wAuCPJvwb+Aritjb8N+E9JZhj85n8VQFUdSHIn8APgZeDaqnpleTdHkjSukQFQVY8B75mj/jRzXMVTVf8buHKe17oJuGnxbUqSlpvfBJakThkAktQpA0CSOmUASFKnDABJ6pQBIEmdMgAkqVMGgCR1ygCQpE4ZAJLUKQNAkjplAEhSpwwASeqUASBJnTIAJKlTBoAkdcoAkKROGQCS1CkDQJI6ZQBIUqcMAEnqlAEgSZ0aGQBJNia5P8nBJAeSfLzVfz/JT5Lsb4/Lhta5IclMkieTXDJU39ZqM0muPzmbJEkaxyljjHkZ+GRVfT/JW4FHkuxty26pqn8zPDjJucBVwDuBXwP+PMlvtMVfAn4bmAUeTrK7qn6wHBsiSVqckQFQVYeBw23650kOAusXWOVy4I6qegn4UZIZYGtbNlNVTwMkuaONNQAkaQIWdQ4gySbgPcCDrXRdkseS7EpyequtB54dWm221earv/Y9diTZl2Tf0aNHF9OeJGkRxg6AJG8Bvg58oqp+BtwK/DqwhcEewh8eHzrH6rVA/dWFqp1VNV1V01NTU+O2J0lapHHOAZDkVAYf/l+tqm8AVNVzQ8v/GLi7zc4CG4dW3wAcatPz1SVJK2ycq4AC3AYcrKovDNXPHhr2YeCJNr0buCrJG5OcA2wGHgIeBjYnOSfJGxicKN69PJshSVqscfYALgQ+AjyeZH+rfQq4OskWBodxngF+F6CqDiS5k8HJ3ZeBa6vqFYAk1wH3AuuAXVV1YBm3RZK0CONcBfQ95j5+v2eBdW4Cbpqjvmeh9SRJK8dvAktSpwwASeqUASBJnTIAJKlTBoAkdcoAkKROGQCS1CkDQJI6ZQBIUqcMAEnqlAEgSZ0yACSpUwaAJHXKAJCkThkAktQpA0CSOmUASFKnDABJ6pQBIEmdMgAkqVMGgCR1amQAJNmY5P4kB5McSPLxVj8jyd4kT7Xn01s9Sb6YZCbJY0nOG3qt7W38U0m2n7zNkiSNMs4ewMvAJ6vqHcAFwLVJzgWuB+6rqs3AfW0e4FJgc3vsAG6FQWAANwLnA1uBG4+HhiRp5Y0MgKo6XFXfb9M/Bw4C64HLgdvbsNuBK9r05cBXauAB4LQkZwOXAHur6lhVPQ/sBbYt69ZIksa2qHMASTYB7wEeBN5WVYdhEBLAWW3YeuDZodVmW22+uiRpAsYOgCRvAb4OfKKqfrbQ0DlqtUD9te+zI8m+JPuOHj06bnuSpEUaKwCSnMrgw/+rVfWNVn6uHdqhPR9p9Vlg49DqG4BDC9Rfpap2VtV0VU1PTU0tZlskSYswzlVAAW4DDlbVF4YW7QaOX8mzHbhrqP7RdjXQBcAL7RDRvcDFSU5vJ38vbjVJ0gScMsaYC4GPAI8n2d9qnwJuBu5Mcg3wY+DKtmwPcBkwA7wIfAygqo4l+RzwcBv32ao6tixbIUlatJEBUFXfY+7j9wDvm2N8AdfO81q7gF2LaVCSdHL4TWBJ6pQBIEmdMgAkqVMGgCR1ygCQpE4ZAJLUKQNAkjplAEhSpwwASeqUASBJnTIAJKlTBoAkdcoAkKROGQCS1CkDQJI6ZQBIUqcMAEnqlAEgSZ0yACSpUwaAJHXKAJCkTo0MgCS7khxJ8sRQ7feT/CTJ/va4bGjZDUlmkjyZ5JKh+rZWm0ly/fJviiRpMcbZA/gysG2O+i1VtaU99gAkORe4CnhnW+ePkqxLsg74EnApcC5wdRsrSZqQU0YNqKrvJtk05utdDtxRVS8BP0oyA2xty2aq6mmAJHe0sT9YdMeSpGVxIucArkvyWDtEdHqrrQeeHRoz22rz1SVJE7LUALgV+HVgC3AY+MNWzxxja4H66yTZkWRfkn1Hjx5dYnuSpFGWFABV9VxVvVJVvwT+mL8+zDMLbBwaugE4tEB9rtfeWVXTVTU9NTW1lPYkSWNYUgAkOXto9sPA8SuEdgNXJXljknOAzcBDwMPA5iTnJHkDgxPFu5fetiTpRI08CZzka8B7gTOTzAI3Au9NsoXBYZxngN8FqKoDSe5kcHL3ZeDaqnqlvc51wL3AOmBXVR1Y9q2RJI1tnKuArp6jfNsC428CbpqjvgfYs6juJEknjd8ElqROjdwDkLRyNl1/z6RbGMszN39g0i1oGbgHIEmdcg9glVorvwlKWrvcA5CkThkAktQpA0CSOmUASFKnDABJ6pQBIEmdMgAkqVMGgCR1ygCQpE4ZAJLUKQNAkjplAEhSpwwASeqUASBJnTIAJKlTBoAkdcoAkKROGQCS1KmRAZBkV5IjSZ4Yqp2RZG+Sp9rz6a2eJF9MMpPksSTnDa2zvY1/Ksn2k7M5kqRxjbMH8GVg22tq1wP3VdVm4L42D3ApsLk9dgC3wiAwgBuB84GtwI3HQ0OSNBkjA6Cqvgsce035cuD2Nn07cMVQ/Ss18ABwWpKzgUuAvVV1rKqeB/by+lCRJK2gpZ4DeFtVHQZoz2e1+nrg2aFxs602X/11kuxIsi/JvqNHjy6xPUnSKMt9Ejhz1GqB+uuLVTurarqqpqemppa1OUnSX1tqADzXDu3Qno+0+iywcWjcBuDQAnVJ0oQsNQB2A8ev5NkO3DVU/2i7GugC4IV2iOhe4OIkp7eTvxe3miRpQk4ZNSDJ14D3AmcmmWVwNc/NwJ1JrgF+DFzZhu8BLgNmgBeBjwFU1bEknwMebuM+W1WvPbEsSVpBIwOgqq6eZ9H75hhbwLXzvM4uYNeiupMknTR+E1iSOmUASFKnDABJ6pQBIEmdMgAkqVMGgCR1ygCQpE4ZAJLUKQNAkjplAEhSpwwASeqUASBJnTIAJKlTBoAkdcoAkKROGQCS1CkDQJI6ZQBIUqcMAEnqlAEgSZ0yACSpUycUAEmeSfJ4kv1J9rXaGUn2JnmqPZ/e6knyxSQzSR5Lct5ybIAkaWmWYw/gH1bVlqqabvPXA/dV1WbgvjYPcCmwuT12ALcuw3tLkpboZBwCuhy4vU3fDlwxVP9KDTwAnJbk7JPw/pKkMZxoABTwrSSPJNnRam+rqsMA7fmsVl8PPDu07myrSZIm4JQTXP/CqjqU5Cxgb5IfLjA2c9TqdYMGQbID4O1vf/sJtidJms8J7QFU1aH2fAT4JrAVeO74oZ32fKQNnwU2Dq2+ATg0x2vurKrpqpqempo6kfYkSQtYcgAk+dUkbz0+DVwMPAHsBra3YduBu9r0buCj7WqgC4AXjh8qkiStvBM5BPQ24JtJjr/On1TVf03yMHBnkmuAHwNXtvF7gMuAGeBF4GMn8N6SpBO05ACoqqeBd89R/1/A++aoF3DtUt9PkrS8/CawJHXKAJCkThkAktQpA0CSOmUASFKnDABJ6pQBIEmdMgAkqVMGgCR1ygCQpE4ZAJLUqRP9/wAkdWjT9fdMuoWRnrn5A5NuYdXrLgDWwj9cSVoJHgKSpE4ZAJLUKQNAkjplAEhSpwwASeqUASBJnTIAJKlTBoAkdWrFAyDJtiRPJplJcv1Kv78kaWBFAyDJOuBLwKXAucDVSc5dyR4kSQMrvQewFZipqqer6v8AdwCXr3APkiRW/l5A64Fnh+ZngfPnG/z4T17w3j2SlmStfHZM8qZ1Kx0AmaNWrxqQ7AB2tNmX/vLzH3zipHd14s4E/mrSTYzBPpeXfS6vtdDnsveYzy/nq/0/vznOoJUOgFlg49D8BuDQ8ICq2gnsBEiyr6qmV669pbHP5WWfy8s+l89a6BEGfY4zbqXPATwMbE5yTpI3AFcBu1e4B0kSK7wHUFUvJ7kOuBdYB+yqqgMr2YMkaWDF/0OYqtoD7Blz+M6T2csyss/lZZ/Lyz6Xz1roEcbsM1U1epQk6f873gpCkjq1ZgIgyT9PUknOnHQvc0nyuSSPJdmf5FtJfm3SPc0lyR8k+WHr9ZtJTpt0T3NJcmWSA0l+mWRVXXWxVm5nkmRXkiNJVu2l1Ek2Jrk/ycH29/3xSfc0lyRvSvJQkkdbn5+ZdE8LSbIuyV8kuXuhcWsiAJJsBH4b+PGke1nAH1TVb1XVFuBu4F9NuqF57AXeVVW/Bfx34IYJ9zOfJ4B/BHx30o0MW2O3M/kysG3STYzwMvDJqnoHcAFw7Sr983wJuKiq3g1sAbYluWDCPS3k48DBUYPWRAAAtwD/ktd8aWw1qaqfDc3+Kqu016r6VlW93GYfYPBdjFWnqg5W1ZOT7mMOa+Z2JlX1XeDYpPtYSFUdrqrvt+mfM/jQWj/Zrl6vBn7RZk9tj1X5M55kA/AB4D+MGrvqAyDJh4CfVNWjk+5llCQ3JXkW+Mes3j2AYf8U+LNJN7HGzHU7k1X3gbUWJdkEvAd4cLKdzK0dVtkPHAH2VtWq7BP4twx+Yf7lqIErfhnoXJL8OfC351j0aeBTwMUr29HcFuqzqu6qqk8Dn05yA3AdcOOKNtiM6rON+TSD3e+vrmRvw8bpcxUaeTsTLV6StwBfBz7xmr3pVaOqXgG2tPNm30zyrqpaVedXknwQOFJVjyR576jxqyIAqur9c9WT/D3gHODRJDA4XPH9JFur6n+uYIvA/H3O4U+Ae5hQAIzqM8l24IPA+2qC1wEv4s9zNRl5OxMtTpJTGXz4f7WqvjHpfkapqp8m+Q6D8yurKgCAC4EPJbkMeBPwN5L856r6J3MNXtWHgKrq8ao6q6o2VdUmBj98503iw3+UJJuHZj8E/HBSvSwkyTbg94APVdWLk+5nDfJ2Jssog9/sbgMOVtUXJt3PfJJMHb9iLsmbgfezCn/Gq+qGqtrQPi+vAr4934c/rPIAWGNuTvJEkscYHLJalZezAf8OeCuwt12y+u8n3dBcknw4ySzwD4B7ktw76Z5gcDsTBof37mVwwvLO1Xo7kyRfA/4b8JtJZpNcM+me5nAh8BHgovbvcX/77XW1ORu4v/18P8zgHMCCl1iuBX4TWJI65R6AJHXKAJCkThkAktQpA0CSOmUASFKnDABJ6pQBIEmdMgAkqVP/F3UUdiNd4L1WAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.bar(bin_edges[:-1], hist, width = 1)\n",
    "plt.xlim(min(bin_edges), max(bin_edges))\n",
    "plt.show() "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 131,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[  1   3   5   5   5   4   8   7  15  14  12  19  17  17  23  31  44  44\n",
      "  32  44  55  67  57  69  86  93  88  94 105 125 126 134 120 151 168 155\n",
      " 180 169 172 205 180 236 216 242 233 223 234 236 227 237 239 257 228 239\n",
      " 238 228 196 191 250 203 193 185 179 175 155 157 125 139 113 129 104 124\n",
      "  81  98  91  54  66  79  49  45  47  47  21  29  25  24  15  14  17  25\n",
      "  17   9  14   6   6   6   5   3   5]\n",
      "[-3.         -2.93939394 -2.87878788 -2.81818182 -2.75757576 -2.6969697\n",
      " -2.63636364 -2.57575758 -2.51515152 -2.45454545 -2.39393939 -2.33333333\n",
      " -2.27272727 -2.21212121 -2.15151515 -2.09090909 -2.03030303 -1.96969697\n",
      " -1.90909091 -1.84848485 -1.78787879 -1.72727273 -1.66666667 -1.60606061\n",
      " -1.54545455 -1.48484848 -1.42424242 -1.36363636 -1.3030303  -1.24242424\n",
      " -1.18181818 -1.12121212 -1.06060606 -1.         -0.93939394 -0.87878788\n",
      " -0.81818182 -0.75757576 -0.6969697  -0.63636364 -0.57575758 -0.51515152\n",
      " -0.45454545 -0.39393939 -0.33333333 -0.27272727 -0.21212121 -0.15151515\n",
      " -0.09090909 -0.03030303  0.03030303  0.09090909  0.15151515  0.21212121\n",
      "  0.27272727  0.33333333  0.39393939  0.45454545  0.51515152  0.57575758\n",
      "  0.63636364  0.6969697   0.75757576  0.81818182  0.87878788  0.93939394\n",
      "  1.          1.06060606  1.12121212  1.18181818  1.24242424  1.3030303\n",
      "  1.36363636  1.42424242  1.48484848  1.54545455  1.60606061  1.66666667\n",
      "  1.72727273  1.78787879  1.84848485  1.90909091  1.96969697  2.03030303\n",
      "  2.09090909  2.15151515  2.21212121  2.27272727  2.33333333  2.39393939\n",
      "  2.45454545  2.51515152  2.57575758  2.63636364  2.6969697   2.75757576\n",
      "  2.81818182  2.87878788  2.93939394  3.        ]\n"
     ]
    }
   ],
   "source": [
    "bins=np.linspace(-3,3,100)\n",
    "\n",
    "hist, bin_edges=np.histogram(noise, bins)\n",
    "print(hist)\n",
    "print(bin_edges)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 132,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD8CAYAAAB5Pm/hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAADrVJREFUeJzt3X+o3Xd9x/Hny9q5YR2u5LaLacItEsfq2OK4dELH6KjT2opp/6i0DBudEIWWKTgwKqzOEYhs1rmxlUVaTKH+CLRioNnW2jmcf0RNQmZbY7eg0aYJSdSqFcGR5L0/7jf0LN7cc+495+Sc+8nzAYdzzud+vue8vyS87ofP9/P93FQVkqR2vWTSBUiSxsugl6TGGfSS1DiDXpIaZ9BLUuMMeklqnEEvSY0z6CWpcQa9JDXupZMuAGDVqlU1Ozs76TIkaUXZt2/fD6pqpl+/qQj62dlZ9u7dO+kyJGlFSfK9Qfo5dSNJjTPoJalxBr0kNc6gl6TGGfSS1DiDXpIaZ9BLUuMMeklqnEEvSY2bijtjpWHNbnl00iWM1eFtN0+6BK1gjuglqXF9gz7J2iRfTnIwydNJ3tu1fyTJc0kOdI+beo75YJJDSZ5J8qZxnoAkaXGDTN2cAt5fVfuTvALYl+Tx7mefqKq/7e2c5BrgduC1wKuALyV5TVWdHmXhkqTB9B3RV9WxqtrfvX4BOAisWeSQjcDnquoXVfVd4BBw7SiKlSQt3ZIuxiaZBV4HfA24Drg7yZ3AXuZH/c8z/0tgT89hR1jgF0OSzcBmgHXr1i2jdOniMc6LzV7obd/AF2OTXAY8DLyvqn4K3Ae8GtgAHAM+frbrAofXLzVUba+quaqam5npu2++JGmZBgr6JJcyH/IPVdUjAFV1vKpOV9UZ4FO8OD1zBFjbc/hVwNHRlSxJWopBVt0EuB84WFX39rSv7ul2K/BU93oXcHuSlyW5GlgPfH10JUuSlmKQOfrrgLcDTyY50LV9CLgjyQbmp2UOA+8GqKqnk+wEvsX8ip27XHEjSZPTN+ir6qssPO++e5FjtgJbh6hLkjQi3hkrSY1zrxuNTev7z0grhSN6SWqcQS9JjTPoJalxBr0kNc6LsRqKF1yl6eeIXpIaZ9BLUuMMeklqnEEvSY0z6CWpcQa9JDXOoJekxhn0ktQ4g16SGmfQS1LjDHpJapxBL0mNM+glqXEGvSQ1zqCXpMYZ9JLUOINekhpn0EtS4wx6SWqcQS9JjTPoJalxBr0kNc6gl6TGGfSS1LiX9uuQZC3wIPCbwBlge1V9MsnlwOeBWeAw8Laqej5JgE8CNwE/B95RVfvHU74uhNktj066BElD6Bv0wCng/VW1P8krgH1JHgfeATxRVduSbAG2AB8A3gys7x5/ANzXPUuaQuf+Ij+87eYJVaJx6Tt1U1XHzo7Iq+oF4CCwBtgI7Oi67QBu6V5vBB6seXuAVyZZPfLKJUkDWdIcfZJZ4HXA14Arq+oYzP8yAK7ouq0Bnu057EjXJkmagIGDPsllwMPA+6rqp4t1XaCtFvi8zUn2Jtl78uTJQcuQJC3RQEGf5FLmQ/6hqnqkaz5+dkqmez7RtR8B1vYcfhVw9NzPrKrtVTVXVXMzMzPLrV+S1EffoO9W0dwPHKyqe3t+tAvY1L3eBHyxp/3OzHs98JOzUzySpAtvkFU31wFvB55McqBr+xCwDdiZ5F3A94Hbup/tZn5p5SHml1e+c6QVS5KWpG/QV9VXWXjeHeCGBfoXcNeQdUmSRmSQEb2ki0jvunrX1LfBLRAkqXEGvSQ1zqCXpMYZ9JLUOINekhpn0EtS41xeKem8XGrZBkf0ktQ4g16SGmfQS1LjnKO/iPm3YKWLg0EvaSD+bdmVy6kbSWqcQS9JjXPqRtKyuMZ+5XBEL0mNM+glqXEGvSQ1zqCXpMZ5MbZx3hQlyRG9JDXOoJekxhn0ktQ4g16SGmfQS1LjDHpJapxBL0mNM+glqXHeMCVpaIvdmOfOlpPniF6SGtc36JM8kOREkqd62j6S5LkkB7rHTT0/+2CSQ0meSfKmcRUuSRrMICP6TwM3LtD+iara0D12AyS5BrgdeG13zD8luWRUxUqSlq5v0FfVV4AfDfh5G4HPVdUvquq7wCHg2iHqkyQNaZiLsXcnuRPYC7y/qp4H1gB7evoc6dp0AbljpaRey70Yex/wamADcAz4eNeeBfrWQh+QZHOSvUn2njx5cpllSJL6WVbQV9XxqjpdVWeAT/Hi9MwRYG1P16uAo+f5jO1VNVdVczMzM8spQ5I0gGUFfZLVPW9vBc6uyNkF3J7kZUmuBtYDXx+uREnSMPrO0Sf5LHA9sCrJEeAe4PokG5ifljkMvBugqp5OshP4FnAKuKuqTo+ndEnSIPoGfVXdsUDz/Yv03wpsHaYoSdLoeGesJDXOoJekxhn0ktQ4g16SGuc2xZLGyi2MJ88RvSQ1zqCXpMYZ9JLUOOfoVyh3qJQ0KEf0ktQ4g16SGmfQS1LjnKOXNDHnu9bk+vrRckQvSY1zRL+CuNJG0nIY9FPMYJc0Cga9pKnj/jij5Ry9JDXOoJekxhn0ktQ4g16SGmfQS1LjDHpJapxBL0mNM+glqXEGvSQ1zqCXpMYZ9JLUOINekhpn0EtS4wx6SWqc2xRPGfeglzRqfUf0SR5IciLJUz1tlyd5PMn/dM+/0bUnyd8nOZTkm0l+f5zFS5L6G2Tq5tPAjee0bQGeqKr1wBPde4A3A+u7x2bgvtGUKUlarr5BX1VfAX50TvNGYEf3egdwS0/7gzVvD/DKJKtHVawkaemWezH2yqo6BtA9X9G1rwGe7el3pGv7JUk2J9mbZO/JkyeXWYYkqZ9Rr7rJAm21UMeq2l5Vc1U1NzMzM+IyJElnLXfVzfEkq6vqWDc1c6JrPwKs7el3FXB0mAIlqZd/OHzpljui3wVs6l5vAr7Y035nt/rm9cBPzk7xSJImo++IPslngeuBVUmOAPcA24CdSd4FfB+4reu+G7gJOAT8HHjnGGqWJC1B36CvqjvO86MbFuhbwF3DFnWx8SYpSePknbEXiGEuaVIMeknN6B1QeWH2RW5qJkmNM+glqXEGvSQ1zqCXpMZ5MXZMXGUjaVo4opekxhn0ktQ4g16SGmfQS1LjvBgrqUlLWRDR+l20juglqXEGvSQ1zqCXpMYZ9JLUOINekhpn0EtS4wx6SWqc6+hHyI3MJE0jR/SS1DhH9EvkqF3SSuOIXpIaZ9BLUuMMeklqnEEvSY0z6CWpcQa9JDXOoJekxhn0ktQ4b5gagDdJSVrJhgr6JIeBF4DTwKmqmktyOfB5YBY4DLytqp4frkxJ0nKNYurmj6tqQ1XNde+3AE9U1Xrgie69JGlCxjFHvxHY0b3eAdwyhu+QJA1o2KAv4LEk+5Js7tqurKpjAN3zFUN+hyRpCMNejL2uqo4muQJ4PMm3Bz2w+8WwGWDdunVDliFJy9e74OLwtpsnWMl4DDWir6qj3fMJ4AvAtcDxJKsBuucT5zl2e1XNVdXczMzMMGVIkhax7BF9kpcDL6mqF7rXbwQ+CuwCNgHbuucvjqJQSboQBl1OvZJG/sNM3VwJfCHJ2c/5TFX9a5JvADuTvAv4PnDb8GVKkpZr2UFfVd8Bfm+B9h8CNwxT1Lh445Oki5FbIEhS4wx6SWqcQS9JjTPoJalxBr0kNc6gl6TGGfSS1DiDXpIaZ9BLUuMMeklqXPN/M9ZtDyRd7JoPekkah8UGkdO2s6VTN5LUOINekhpn0EtS45yjl6QRG/UikGHn/JsMelfaSNKLnLqRpMYZ9JLUuBU1deOUjCQtnSN6SWrcihrRS9LFqHc2YzkrcBzRS1LjpmJE/+RzP3H+XZLGZCqCXpI0mOUMip26kaTGGfSS1DiDXpIaZ9BLUuMMeklqnEEvSY0bW9AnuTHJM0kOJdkyru+RJC1uLEGf5BLgH4E3A9cAdyS5ZhzfJUla3LhG9NcCh6rqO1X1v8DngI1j+i5J0iLGFfRrgGd73h/p2iRJF9i4tkDIAm31/zokm4HN3dtffO9jb3lqTLVMg1XADyZdxBh5fitXy+cG7Z/fbw3SaVxBfwRY2/P+KuBob4eq2g5sB0iyt6rmxlTLxHl+K1vL59fyucHFcX6D9BvX1M03gPVJrk7yK8DtwK4xfZckaRFjGdFX1akkdwP/BlwCPFBVT4/juyRJixvbNsVVtRvYPWD37eOqY0p4fitby+fX8rmB5wdAqqp/L0nSiuUWCJLUuKkJ+iR/neSbSQ4keSzJqyZd0ygl+Zsk3+7O8QtJXjnpmkYpyW1Jnk5yJkkTqxxa3sYjyQNJTiRpcllzkrVJvpzkYPf/8r2TrmmUkvxqkq8n+a/u/P5q0f7TMnWT5Ner6qfd6z8Hrqmq90y4rJFJ8kbg37sL1R8DqKoPTLiskUny28AZ4J+Bv6iqgZZ9TatuG4//Bv6E+eXC3wDuqKpvTbSwEUnyR8DPgAer6ncmXc+oJVkNrK6q/UleAewDbmno3y/Ay6vqZ0kuBb4KvLeq9izUf2pG9GdDvvNyzrnBaqWrqseq6lT3dg/z9xY0o6oOVtUzk65jhJrexqOqvgL8aNJ1jEtVHauq/d3rF4CDNHR3fs37Wff20u5x3sycmqAHSLI1ybPAnwJ/Oel6xujPgH+ZdBFalNt4NCLJLPA64GuTrWS0klyS5ABwAni8qs57fhc06JN8KclTCzw2AlTVh6tqLfAQcPeFrG0U+p1f1+fDwCnmz3FFGeT8GtJ3Gw9NvySXAQ8D7ztn1mDFq6rTVbWB+dmBa5OcdwpubOvoF1JVbxiw62eAR4F7xljOyPU7vySbgLcAN9S0XBxZgiX8+7Wg7zYemm7d3PXDwENV9cik6xmXqvpxkv8AbgQWvLg+NVM3Sdb3vH0r8O1J1TIOSW4EPgC8tap+Pul61JfbeKxg3cXK+4GDVXXvpOsZtSQzZ1fuJfk14A0skpnTtOrmYeZ3YjsDfA94T1U9N9mqRifJIeBlwA+7pj2NrSq6FfgHYAb4MXCgqt402aqGk+Qm4O94cRuPrRMuaWSSfBa4nvndHY8D91TV/RMtaoSS/CHwn8CTzGcKwIe6O/ZXvCS/C+xg/v/mS4CdVfXR8/aflqCXJI3H1EzdSJLGw6CXpMYZ9JLUOINekhpn0EtS4wx6SWqcQS9JjTPoJalx/wdCjuuc0DmNfwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.bar(bin_edges[:-1], hist, width = 1)\n",
    "plt.xlim(min(bin_edges), max(bin_edges))\n",
    "plt.show() "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "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.6.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}