notebooks/04.05-Histograms-and-Binnings.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Histograms, Binnings, and Density"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "A simple histogram can be a great first step in understanding a dataset.\n",
    "Earlier, we saw a preview of Matplotlib's histogram function (discussed in [Comparisons, Masks, and Boolean Logic](02.06-Boolean-Arrays-and-Masks.ipynb)), which creates a basic histogram in one line, once the normal boilerplate imports are done (see the following figure):"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "plt.style.use('seaborn-white')\n",
    "\n",
    "rng = np.random.default_rng(1701)\n",
    "data = rng.normal(size=1000)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXEAAAD0CAYAAABtjRZ7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAQEUlEQVR4nO3dfUiV9//H8dfxWLN5lJBiTezGvtsYKW2TQ9sf1gJrJxytYoY1ScrYRmw2oVl21qxWWaPZYJoVDmKUUTkhhMb4NUnCajmEFRotgtoqbTMk9Jy1tHOu3x+xM6zl3Y5effT5+Ktz5+d9gp58znWuKx2WZVkCABgpwu4BAAADR8QBwGBEHAAMRsQBwGBEHAAMFjlUC/31119qbGzU+PHj5XQ6h2pZADBaIBBQa2urkpOTFRUV9cjjQxbxxsZGZWVlDdVyADCsVFRUyO12P3L/kEV8/PjxoUEmTJgwVMsCgNFu3bqlrKysUEMfNmQR//sQyoQJE5SQkDBUywLAsPC4w9B8sQkAButxJ97V1SWv16ubN2+qs7NTq1at0rPPPqv3339fU6ZMkSQtXbpU6enpKi0tVW1trSIjI+X1ejV9+vShmB8ARrQeI15dXa2xY8dq586dunPnjhYuXKgPPvhAK1asUE5OTuh5TU1Nqq+vV2VlpVpaWpSbm6uqqqpBHx4ARroeIz5v3jx5PB5JkmVZcjqdamxs1NWrV1VTU6PJkyfL6/WqoaFBqampcjgcio+PVyAQUFtbm+Li4obkTQDASNVjxKOjoyVJPp9Pq1evVl5enjo7O7V48WIlJydrz5492r17t2JiYjR27Nhur+vo6CDiADDIev1is6WlRdnZ2VqwYIHmz5+vuXPnKjk5WZI0d+5cXbx4US6XS36/P/Qav9+vmJiYwZsaACCpl4jfvn1bOTk5ys/PV0ZGhiRp5cqVunDhgiTp7NmzSkpKUkpKiurq6hQMBtXc3KxgMMguHACGQI+HU/bu3av29naVlZWprKxMklRQUKCioiKNGjVK48aN05YtW+RyueR2u5WZmalgMKjCwsIhGR4IpykFx21Z99qON21ZF8ODY6h+s8+NGzeUlpammpoaLvbBY9kVUjsRcfSkt3ZysQ8AGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBInt6sKurS16vVzdv3lRnZ6dWrVql5557TgUFBXI4HHr++ee1ceNGRUREqLS0VLW1tYqMjJTX69X06dOH6j0AwIjVY8Srq6s1duxY7dy5U3fu3NHChQv14osvKi8vT6+++qoKCwtVU1Oj+Ph41dfXq7KyUi0tLcrNzVVVVdVQvQcAGLF6jPi8efPk8XgkSZZlyel0qqmpSTNmzJAkzZo1S6dPn1ZiYqJSU1PlcDgUHx+vQCCgtrY2xcXFDf47AIARrMdj4tHR0XK5XPL5fFq9erXy8vJkWZYcDkfo8Y6ODvl8Prlcrm6v6+joGNzJAQC9f7HZ0tKi7OxsLViwQPPnz1dExD8v8fv9io2Nlcvlkt/v73Z/TEzM4EwMAAjpMeK3b99WTk6O8vPzlZGRIUmaNm2azp07J0k6deqU3G63UlJSVFdXp2AwqObmZgWDQQ6lAMAQ6PGY+N69e9Xe3q6ysjKVlZVJkj755BNt3bpVu3bt0tSpU+XxeOR0OuV2u5WZmalgMKjCwsIhGR4ARjqHZVnWUCx048YNpaWlqaamRgkJCUOxJAw0peC43SMMuWs73rR7BDzBemsnF/sAgMGIOAAYjIgDgMGIOAAYjIgDgMGIOAAYjIgDgMGIOAAYjIgDgMGIOAAYjIgDgMGIOAAYjIgDgMGIOAAYjIgDgMGIOAAYjIgDgMGIOAAYjIgDgMGIOAAYrMffdo+RayT+wmLAROzEAcBgRBwADEbEAcBgRBwADEbEAcBgRBwADEbEAcBgRBwADEbEAcBgRBwADMZl94DN7PwvDq7teNO2tREe7MQBwGBEHAAM1qeInz9/XsuWLZMkXbx4UTNnztSyZcu0bNkyfffdd5Kk0tJSZWRkaMmSJbpw4cLgTQwACOn1mHh5ebmqq6s1ZswYSVJTU5NWrFihnJyc0HOamppUX1+vyspKtbS0KDc3V1VVVYM3NQBAUh924pMmTVJJSUnodmNjo2pra5WVlSWv1yufz6eGhgalpqbK4XAoPj5egUBAbW1tgzo4AKAPEfd4PIqM/GfDPn36dK1du1YVFRWaOHGidu/eLZ/PJ5fLFXpOdHS0Ojo6BmdiAEBIv7/YnDt3rpKTk0N/vnjxolwul/x+f+g5fr9fMTEx4ZsSAPCv+h3xlStXhr64PHv2rJKSkpSSkqK6ujoFg0E1NzcrGAwqLi4u7MMCALrr98U+mzZt0pYtWzRq1CiNGzdOW7ZskcvlktvtVmZmpoLBoAoLCwdjVgDAQ/oU8YSEBB09elSSlJSUpMOHDz/ynNzcXOXm5oZ3OgBAj7jYBwAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGBEHAAMRsQBwGB9ivj58+e1bNkySdKvv/6qpUuX6p133tHGjRsVDAYlSaWlpcrIyNCSJUt04cKFwZsYABDSa8TLy8u1YcMG3bt3T5K0fft25eXl6dChQ7IsSzU1NWpqalJ9fb0qKyu1a9cubd68edAHBwD0IeKTJk1SSUlJ6HZTU5NmzJghSZo1a5bOnDmjhoYGpaamyuFwKD4+XoFAQG1tbYM3NQBAUh8i7vF4FBkZGbptWZYcDockKTo6Wh0dHfL5fHK5XKHn/H0/AGBw9fuLzYiIf17i9/sVGxsrl8slv9/f7f6YmJjwTAgAeKx+R3zatGk6d+6cJOnUqVNyu91KSUlRXV2dgsGgmpubFQwGFRcXF/ZhAQDdRfb+lO7WrVunTz/9VLt27dLUqVPl8XjkdDrldruVmZmpYDCowsLCwZgVAPCQPkU8ISFBR48elSQlJibq4MGDjzwnNzdXubm54Z0OANAjLvYBAIMRcQAwGBEHAIMRcQAwGBEHAIMRcQAwGBEHAIMRcQAwGBEHAIMRcQAwGBEHAIMRcQAwGBEHAIMRcQAwGBEHAIMRcQAwGBEHAIMRcQAwGBEHAIMRcQAwWL9/2z2GzpSC43aPAOAJx04cAAxGxAHAYBxOAUYwuw7ZXdvxpi3rDkfsxAHAYEQcAAxGxAHAYEQcAAxGxAHAYEQcAAxGxAHAYEQcAAxGxAHAYEQcAAw24MvuFy1aJJfLJUlKSEhQZmamtm3bJqfTqdTUVH344YdhGxIA8O8GFPF79+7JsiwdOHAgdN+CBQtUUlKiiRMn6r333tPFixc1bdq0sA0KAHjUgA6nXLp0SXfv3lVOTo6ys7P1008/qbOzU5MmTZLD4VBqaqrOnDkT7lkBAA8Z0E48KipKK1eu1OLFi3Xt2jW9++67io2NDT0eHR2t69evh21IAMC/G1DEExMTNXnyZDkcDiUmJiomJkZ37twJPe73+7tFHQAwOAZ0OOXbb7/Vjh07JEm///677t69q6efflq//fabLMtSXV2d3G53WAcFADxqQDvxjIwMrV+/XkuXLpXD4VBRUZEiIiL08ccfKxAIKDU1VS+99FK4ZwUAPGRAER89erSKi4sfuf/o0aP/eSAAQN9xsQ8AGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGIyIA4DBiDgAGCzS7gFMMKXguN0jAMOKXf+mru1405Z1BxM7cQAwGBEHAIMRcQAwGBEHAIMRcQAwGBEHAIMRcQAwGBEHAIMZdbEPF90AQHfsxAHAYGHdiQeDQW3atEm//PKLRo8era1bt2ry5MnhXAIABszOT/ODdcl/WHfiP/zwgzo7O3XkyBGtWbNGO3bsCOePBwA8JKw78YaGBs2cOVOS9PLLL6uxsTH0WCAQkCTdunVr4Av42/7TfABglxs3bgzodX838++GPiysEff5fHK5XKHbTqdT9+/fV2RkpFpbWyVJWVlZA/75T/3nCQHAHmn/t/U/vb61tfVfD0+HNeIul0t+vz90OxgMKjLywRLJycmqqKjQ+PHj5XQ6w7ksAAxbgUBAra2tSk5O/tfHwxrxlJQUnTx5Uunp6fr555/1wgsvhB6LioqS2+0O53IAMCL0dIKIw7IsK1wL/X12yuXLl2VZloqKivS///0vXD8eAPCQsEb8SfPnn39qzZo1am9v16hRo/T555/rmWeesXusXnV0dCg/P18+n09dXV0qKCjQK6+8YvdYfXbixAl9//33Ki4utnuUHpl+Suz58+f1xRdf6MCBA3aP0iddXV3yer26efOmOjs7tWrVKqWlpdk9Vq8CgYA2bNigq1evyuFwaPPmzd2OMthtWF/sc/ToUSUlJamiokJvvfWWysvL7R6pT/bv36/XXntNBw8e1Pbt2/XZZ5/ZPVKfbd26VcXFxQoGg3aP0iuTT4ktLy/Xhg0bdO/ePbtH6bPq6mqNHTtWhw4d0tdff60tW7bYPVKfnDx5UpJ0+PBh5eXl6csvv7R5ou6Muuy+v5YvXx46Lae5uVmxsbE2T9Q3y5cv1+jRoyU92AU89ZQ55+WkpKRozpw5OnLkiN2j9KqnU2KfdJMmTVJJSYnWrl1r9yh9Nm/ePHk8HkmSZVnGnOAwZ84czZ49W9KT2ZFhE/HKykp988033e4rKirS9OnTlZ2drcuXL2v//v02Tfd4Pc3d2tqq/Px8eb1em6Z7vMfNnZ6ernPnztk0Vf/0dErsk87j8Qz4vGO7REdHS3rw97569Wrl5eXZO1A/REZGat26dTpx4oS++uoru8fpzhohrly5YqWlpdk9Rp9dunTJSk9Pt2pra+0epd9+/PFHKy8vz+4xelVUVGQdP348dHvmzJk2TtN/169ftxYvXmz3GP3S3NxsLVq0yKqsrLR7lAH5448/rNmzZ1t+v9/uUUKG9THxffv26dixY5Ie7AJM+fh25coVffTRRyouLtbrr79u9zjDVkpKik6dOiVJj5wSi/C7ffu2cnJylJ+fr4yMDLvH6bNjx45p3759kqQxY8bI4XAoIuLJSeeT/7nxP3j77be1bt06VVVVKRAIqKioyO6R+qS4uFidnZ3atm2bpAcXUe3Zs8fmqYafuXPn6vTp01qyZEnolFgMnr1796q9vV1lZWUqKyuT9OAL2qioKJsn69kbb7yh9evXKysrS/fv35fX632iZh7WpxgCwHD35HwmAAD0GxEHAIMRcQAwGBEHAIMRcQAwGBEHAIMRcQAwGBEHAIP9P4W73E5pca9GAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(data);"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The `hist` function has many options to tune both the calculation and the display; \n",
    "here's an example of a more customized histogram, shown in the following figure:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD0CAYAAACLpN0/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAXXklEQVR4nO3df0xV9/3H8df9wUXGxToyNs02aEdybSozeDVdEwO2QebKGuO6rYDKaszU2KBzpZTNuJYxhpCMrbEOrTVBg4po7IyZSZeSOi4junYkdy02YKxLa9tUsaSRexe5vd7z/cOvt71qe7Byd/XD8/HXPedzPpf3x8RXPvdzfjksy7IEALjjOVNdAABgYhDoAGAIAh0ADEGgA4AhCHQAMIQ7FX/00qVLGhgYUE5OjlwuVypKAIA7zuXLlzU8PKyCggJNmTLluvaUBPrAwICWLVuWij8NAHe8vXv3at68edftT0mg5+TkSLpS1PTp01NRAgDccT788EMtW7YsnqHXSkmgX11mmT59ur71rW+logQAuGN93lI1J0UBwBAEOgAYgkAHAEMQ6ABgCAIdAAxBoAOAIQh0ADBESq5DByabjp5T4zquaoEvyZXAZMzQAcAQtoEei8X0zDPPqLy8XFVVVXrnnXdueMzPf/5zdXZ2Srry8K1169Zp6dKlWrVqlUZGRia+cgBAAttA7+7uViQSUVdXl2pqatTc3HzdMc8995wuXrwY3+7s7JTP59O+ffu0ZMkStbW1TWzVAIDr2AZ6f3+/ioqKJEmFhYUaGBhIaH/55ZflcDjix1zbp7i4WMePH5/ImgEAN2B7UjQUCsnr9ca3XS6XotGo3G63Tp06pb/+9a/asmWL/vznPyf0ycrKkiRlZmZqdHQ0CaUDN4cTkzCdbaB7vV6Fw+H4diwWk9t9pdvhw4d17tw5Pf7443r//feVlpamb37zmwl9wuGwpk6dmqTyAQBX2Qa63+/XsWPHVFZWpmAwKJ/v09nL008/Hf/8/PPP62tf+5qKi4t1+vRp9fT0aPbs2QoEApo7d25yqgcAxNkGemlpqfr6+lRRUSHLstTU1KT29nbl5uaqpKTkhn0qKytVV1enyspKpaWlqbW1dcILBwAksg10p9OphoaGhH35+fnXHbdu3br454yMDG3ZsmUCygPsjXdtHDAdNxYBgCEIdAAwBM9yAa5xM0s4XOKI2wkzdAAwBIEOAIYg0AHAEAQ6ABiCQAcAQxDoAGAIAh0ADEGgA4AhCHQAMASBDgCGINABwBAEOgAYgkAHAEMQ6ABgCNvH58ZiMdXX12toaEgej0eNjY3Ky8uLt+/du1cvvfSSHA6HVq5cqbKyMlmWpeLiYt19992SpMLCQtXU1CRtEACAcQR6d3e3IpGIurq6FAwG1dzcrG3btkmSRkZG1NnZqb/85S8aGxvTD3/4Qz388MN69913NWvWLG3fvj3pAwAAXGG75NLf36+ioiJJV2baAwMD8bbs7GwdPnxYaWlpunDhgtLT0+VwOHTy5EmdO3dOVVVVWrVqlc6cOZO8EQAAJI0j0EOhkLxeb3zb5XIpGo3Gt91ut/bs2aPy8nItXrxYkpSTk6PVq1ero6NDa9asUW1tbRJKBwB8lm2ge71ehcPh+HYsFpPbnbhSs3z5cvX29ur111/XiRMnVFBQoJKSEknSvHnzdP78eVmWNcGlAwA+yzbQ/X6/AoGAJCkYDMrn+/QdimfOnFF1dbUsy1JaWpo8Ho+cTqe2bt2q3bt3S5IGBwc1Y8YMORyOJA0BACCN46RoaWmp+vr6VFFRIcuy1NTUpPb2duXm5qqkpET33nuvysvL5XA4VFRUpPvvv18zZ85UbW2tenp65HK5tHnz5v/FWABgUrMNdKfTqYaGhoR9+fn58c/V1dWqrq5OaL/rrru0Y8eOCSoRADAe3FgEAIYg0AHAEAQ6ABiCQAcAQxDoAGAIAh0ADEGgA4AhCHQAMASBDgCGsL1TFEiVjp5TqS7B1p1QIyYPZugAYAgCHQAMQaADgCEIdAAwBIEOAIYg0AHAEAQ6ABjC9jr0WCym+vp6DQ0NyePxqLGxUXl5efH2vXv36qWXXpLD4dDKlStVVlamS5cuqba2Vh999JEyMzPV0tKi7OzspA4EACY72xl6d3e3IpGIurq6VFNTo+bm5njbyMiIOjs7tX//fu3atUstLS2yLEudnZ3y+Xzat2+flixZora2tqQOAgAwjkDv7+9XUVGRJKmwsFADAwPxtuzsbB0+fFhpaWm6cOGC0tPT5XA4EvoUFxfr+PHjSSofAHCVbaCHQiF5vd74tsvlUjQajW+73W7t2bNH5eXlWrx4cbxPVlaWJCkzM1Ojo6MTXTcA4Bq2ge71ehUOh+PbsVhMbnfi0vvy5cvV29ur119/XSdOnEjoEw6HNXXq1AkuGwBwLdtA9/v9CgQCkqRgMCifzxdvO3PmjKqrq2VZltLS0uTxeOR0OuX3+9XT0yNJCgQCmjt3bpLKBwBcZXuVS2lpqfr6+lRRUSHLstTU1KT29nbl5uaqpKRE9957r8rLy+VwOFRUVKT7779f3/3ud1VXV6fKykqlpaWptbX1fzEWAJjUbAPd6XSqoaEhYV9+fn78c3V1taqrqxPaMzIytGXLlgkqEQAwHtxYBACGINABwBC8sQj/c7zl5/ON99+maoHP/iBMOszQAcAQBDoAGIJABwBDEOgAYAhOimLCcLITSC1m6ABgCAIdAAxBoAOAIQh0ADAEgQ4AhiDQAcAQBDoAGIJABwBDEOgAYAgCHQAMYXvrfywWU319vYaGhuTxeNTY2Ki8vLx4+65du3T06FFJ0oIFC+IvjS4uLtbdd98tSSosLFRNTU1yRgAAkDSOQO/u7lYkElFXV5eCwaCam5u1bds2SdLZs2d15MgRHTx4UE6nU5WVlVq4cKEyMjI0a9Ysbd++PekDAABcYbvk0t/fr6KiIklXZtoDAwPxtunTp2vnzp1yuVxyOByKRqNKT0/XyZMnde7cOVVVVWnVqlU6c+ZM8kYAAJA0jkAPhULyer3xbZfLpWg0KklKS0tTdna2LMtSS0uL7rvvPt1zzz3KycnR6tWr1dHRoTVr1qi2tjZ5IwAASBrHkovX61U4HI5vx2Ixud2fdhsbG9PGjRuVmZmpZ599VpJUUFAgl8slSZo3b57Onz8vy7LkcDgmun4AwP+znaH7/X4FAgFJUjAYlM/36ctpLcvSE088oZkzZ6qhoSEe4lu3btXu3bslSYODg5oxYwZhDgBJZjtDLy0tVV9fnyoqKmRZlpqamtTe3q7c3FzFYjG99tprikQi6u3tlSQ9+eSTWr16tWpra9XT0yOXy6XNmzcnfSAAMNnZBrrT6VRDQ0PCvvz8/PjnN99884b9duzYcYulAfg84307VNUCn/1BMAY3FgGAIQh0ADAEgQ4AhiDQAcAQBDoAGIJABwBDEOgAYAgCHQAMQaADgCEIdAAwBIEOAIYg0AHAEAQ6ABiCQAcAQxDoAGAI2+ehA+N99jaA1GKGDgCGsJ2hx2Ix1dfXa2hoSB6PR42NjcrLy4u379q1S0ePHpUkLViwQNXV1bp06ZJqa2v10UcfKTMzUy0tLcrOzk7eKAAA9jP07u5uRSIRdXV1qaamRs3NzfG2s2fP6siRI9q/f78OHDigf/zjHxocHFRnZ6d8Pp/27dunJUuWqK2tLamDAACMI9D7+/tVVFQkSSosLNTAwEC8bfr06dq5c6dcLpccDoei0ajS09MT+hQXF+v48eNJKh8AcJVtoIdCIXm93vi2y+VSNBqVJKWlpSk7O1uWZamlpUX33Xef7rnnHoVCIWVlZUmSMjMzNTo6mqTyAQBX2Qa61+tVOByOb8diMbndny69j42N6amnnlI4HNazzz57XZ9wOKypU6dOdN0AgGvYBrrf71cgEJAkBYNB+Xy+eJtlWXriiSc0c+ZMNTQ0yOVyxfv09PRIkgKBgObOnZuM2gEAn2F7lUtpaan6+vpUUVEhy7LU1NSk9vZ25ebmKhaL6bXXXlMkElFvb68k6cknn1RlZaXq6upUWVmptLQ0tba2Jn0gADDZ2Qa60+lUQ0NDwr78/Pz45zfffPOG/bZs2XKLpQEAbgY3FgGAIQh0ADAEz3IBDDbe5/BULfDZH4TbHjN0ADAEgQ4AhiDQAcAQBDoAGIJABwBDEOgAYAgCHQAMQaADgCEIdAAwBIEOAIYg0AHAEAQ6ABiCQAcAQxDoAGAI28fnxmIx1dfXa2hoSB6PR42NjcrLy0s4ZmRkRJWVlTpy5IjS09NlWZaKi4t19913S5IKCwtVU1OTlAEAAK6wDfTu7m5FIhF1dXUpGAyqublZ27Zti7f39vaqtbVVw8PD8X3vvvuuZs2ape3btyenagDAdWyXXPr7+1VUVCTpykx7YGAg8QucTrW3t2vatGnxfSdPntS5c+dUVVWlVatW6cyZMxNbNQDgOrYz9FAoJK/XG992uVyKRqNyu690nT9//nV9cnJytHr1aj388MP617/+pdraWh06dGgCywYAXMs20L1er8LhcHw7FovFw/zzFBQUyOVySZLmzZun8+fPy7IsORyOWywXAPB5bJdc/H6/AoGAJCkYDMrns3/34NatW7V7925J0uDgoGbMmEGYA0CS2c7QS0tL1dfXp4qKClmWpaamJrW3tys3N1clJSU37LN69WrV1taqp6dHLpdLmzdvnvDCAQCJbAPd6XSqoaEhYV9+fv51x7366qvxz3fddZd27NgxAeUBAMbLNtBhpo6eU6kuAcAE405RADAEgQ4AhiDQAcAQBDoAGIKTogDG7WZOplctsL9nBROLGToAGIJABwBDEOgAYAgCHQAMQaADgCEIdAAwBIEOAIYg0AHAENxYZBieoghMXszQAcAQzNAB8MvOEMzQAcAQtoEei8X0zDPPqLy8XFVVVXrnnXeuO2ZkZESLFi3S2NiYJOnSpUtat26dli5dqlWrVmlkZGTiKwcAJLAN9O7ubkUiEXV1dammpkbNzc0J7b29vVq5cqWGh4fj+zo7O+Xz+bRv3z4tWbJEbW1tE185ACCBbaD39/erqKhIklRYWKiBgYHEL3A61d7ermnTpt2wT3FxsY4fPz6BJQMAbsT2pGgoFJLX641vu1wuRaNRud1Xus6fP/+GfbKysiRJmZmZGh0dnah6AQCfw3aG7vV6FQ6H49uxWCwe5uPpEw6HNXXq1FssEwBgxzbQ/X6/AoGAJCkYDMrns38Lid/vV09PjyQpEAho7ty5t1gmAMCO7ZJLaWmp+vr6VFFRIcuy1NTUpPb2duXm5qqkpOSGfSorK1VXV6fKykqlpaWptbV1wgsHACSyDXSn06mGhoaEffn5+dcd9+qrr8Y/Z2RkaMuWLRNQHgBgvLixCAAMQaADgCEIdAAwBIEOAIYg0AHAEAQ6ABiCQAcAQxDoAGAIAh0ADEGgA4AhCHQAMAQvib5D8BJfAHaYoQOAIZihA0iK8f6qrFpg/44FjA8zdAAwBDP0FGNtHMBEYYYOAIawnaHHYjHV19draGhIHo9HjY2NysvLi7cfOHBA+/fvl9vt1tq1a/XQQw/p448/1qJFi+LvH124cKEef/zx5I0CAGAf6N3d3YpEIurq6lIwGFRzc7O2bdsmSRoeHlZHR4cOHTqksbExLV26VPPnz9dbb72lRx55RL/5zW+SPgAAwBW2Sy79/f0qKiqSJBUWFmpgYCDe9sYbb2jOnDnyeDzKyspSbm6uBgcHNTAwoJMnT2r58uVav369zp8/n7wRAAAkjSPQQ6GQvF5vfNvlcikajcbbsrKy4m2ZmZkKhUL6zne+o/Xr12vPnj1auHChGhsbk1A6AOCzbAPd6/UqHA7Ht2OxmNxu9w3bwuGwsrKy9MADD+h73/ueJKm0tFRvvfXWRNcNALiGbaD7/X4FAgFJUjAYjJ/olKTZs2erv79fY2NjGh0d1dtvvy2fz6dNmzbpb3/7myTp+PHjmjVrVpLKBwBcZXtStLS0VH19faqoqJBlWWpqalJ7e7tyc3NVUlKiqqoqLV26VJZl6Ze//KXS09NVU1OjjRs3qrOzUxkZGSy5AMD/gG2gO51ONTQ0JOzLz8+Pf37sscf02GOPJbR/+9vfVkdHxwSVCAAYD24sAgBDEOgAYAgCHQAMQaADgCF42iKAlOK56ROHGToAGIIZOgCjTOYZPzN0ADAEgQ4AhiDQAcAQBDoAGIKTogDuCLxQ3R4zdAAwBDP0mzSZL4kCcHtjhg4AhiDQAcAQd+SSCydHANwqE5dPbQM9Foupvr5eQ0ND8ng8amxsVF5eXrz9wIED2r9/v9xut9auXauHHnpIIyMjeuqpp3Tp0iV9/etf1+bNm5WRkZHUgQDAZGcb6N3d3YpEIurq6lIwGFRzc7O2bdsmSRoeHlZHR4cOHTqksbExLV26VPPnz1dbW5seeeQRPfroo9qxY4e6urq0YsWKZI/ltsKvCMAMyfi/nKxZv22g9/f3q6ioSJJUWFiogYGBeNsbb7yhOXPmyOPxyOPxKDc3V4ODg+rv79eaNWskScXFxfrjH/+YEOiXL1+WJH344YdfquiPL5z7Uv0A4Hbw3ntf+VL9rmbm1Qy9lm2gh0Iheb3e+LbL5VI0GpXb7VYoFFJWVla8LTMzU6FQKGF/ZmamRkdHE75zeHhYkrRs2bKbHA4A3Pm23mL/4eHhhKXvq2wD3ev1KhwOx7djsZjcbvcN28LhsLKysuL7p0yZonA4rKlTpyZ8Z0FBgfbu3aucnBy5XK4vPSgAmEwuX76s4eFhFRQU3LDdNtD9fr+OHTumsrIyBYNB+Xyfrv3Mnj1bzz33nMbGxhSJRPT222/L5/PJ7/erp6dHjz76qAKBgObOnZvwnVOmTNG8efNucWgAMPncaGZ+lcOyLOuLOl+9yuXUqVOyLEtNTU0KBALKzc1VSUmJDhw4oK6uLlmWpTVr1mjRokW6cOGC6urqFA6H9dWvflWtra36yle+3JoRAGB8bAPdBP/9739VU1OjixcvKi0tTS0tLfrGN76R6rLGZXR0VLW1tQqFQvrkk0/0q1/9SnPmzEl1WeP2yiuv6OWXX1Zra2uqS/lCdpfn3s7+/e9/6w9/+IM6OjpSXcq4ffLJJ9q4caPef/99RSIRrV27ViUlJakuy9bly5e1adMm/ec//5HD4dBvf/vbhFWLVJsUd4oeOHBAs2bN0t69e7V48WK9+OKLqS5p3Nrb2/XAAw9oz5492rx5sxoaGlJd0rg1NjaqtbVVsVgs1aXY+uzluTU1NWpubk51SePy4osvatOmTRobG0t1KTflyJEjmjZtmvbt26edO3fqd7/7XapLGpdjx45Jkvbv368NGzboT3/6U4orSnRH3il6s1asWBG/zOeDDz647iTt7WzFihXyeDySrswO0tPTU1zR+Pn9fi1cuFBdXV2pLsXWF12eezvLzc3V888/r6effjrVpdyUH/zgB1q0aJEkybKsO+biiIULF+rBBx+UdHtmiXGBfvDgQe3evTthX1NTk2bPnq2f/exnOnXqlNrb21NU3Rf7otqHh4dVW1urjRs3pqi6z/d5dZeVlemf//xniqq6OV90ee7tbNGiRXrvvfdSXcZNy8zMlHTl3339+vXasGFDagu6CW63W3V1dXrllVe0ZcuWVJeTyJpkTp8+bZWUlKS6jJsyODholZWVWX//+99TXcpNO3HihLVhw4ZUl2GrqanJOnr0aHy7qKgohdXcnLNnz1o//elPU13GTfvggw+sH/3oR9bBgwdTXcqXcv78eevBBx+0wuFwqkuJmxRr6C+88IIOHz4s6crM4E75eSdJp0+f1i9+8Qu1trZqwYIFqS7HWH6/X4FAQJKuuzwXE+/ChQtauXKlamtr9ZOf/CTV5Yzb4cOH9cILL0iSMjIy5HA45HTePjF6e/+enCA//vGPVVdXp0OHDuny5ctqampKdUnj1traqkgkot///veSrtzMdfVZOpg4paWl6uvrU0VFRfzyXCTP9u3bdfHiRbW1tamtrU3SlRO8U6ZMSXFlX+z73/++fv3rX2vZsmWKRqPauHHjbVXzpLhsEQAmg9vntwIA4JYQ6ABgCAIdAAxBoAOAIQh0ADAEgQ4AhiDQAcAQBDoAGOL/AI2anw//lbLzAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(data, bins=30, density=True, alpha=0.5,\n",
    "         histtype='stepfilled', color='steelblue',\n",
    "         edgecolor='none');"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The `plt.hist` docstring has more information on other available customization options.\n",
    "I find this combination of `histtype='stepfilled'` along with some transparency `alpha` to be helpful when comparing histograms of several distributions (see the following figure):"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAD0CAYAAAC7KMweAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAS50lEQVR4nO3dfWxUdb7H8U870wIyusiVLD6VaxoqwYaUYgxKRldrg9QYo7BQocWHkN2guD6gREmsWHEsUROiEdBF666yUFwf4wOJDc/17i4pTtxZU7xqtkEI2AosziidTufcP7gcWjrTY2tPz/w679dfc85v2vm2hE+/8zu/85scy7IsAQCMket1AQCA/iG4AcAwBDcAGIbgBgDDENwAYBi/2y9w4sQJRSIRjRs3Tj6fz+2XAwDjdXV1qa2tTcXFxRo5cmSvcdeDOxKJaMGCBW6/DAAMOxs2bNDll1/e67zrwT1u3Di7gPHjx7v9cgBgvEOHDmnBggV2fp7J9eA+NT0yfvx4XXTRRW6/HAAMG+mml7k4CQCGIbgBwDAENwAYhuAGAMMQ3ABgGMdVJclkUitWrNC+ffuUn5+vlStXasKECfb4ypUrtXfvXo0ePVqStGbNGp199tnuVQwAWc4xuBsbGxWPx9XQ0KBwOKy6ujqtXbvWHv/Xv/6l9evXa+zYsa4WCgA4yTG4m5ubFQwGJUklJSWKRCL2WDKZVGtrq2pqatTe3q45c+Zozpw57lUL4zR+cdh+fP3kX3tYCTB8OAZ3NBpVIBCwj30+nxKJhPx+v3788UdVVVXpzjvvVFdXlxYuXKji4mJNmjTJ1aIBIJs5XpwMBAKKxWL2cTKZlN9/Mu9HjRqlhQsXatSoUQoEApo+fbpaWlrcqxYA4BzcpaWl2rlzpyQpHA6rqKjIHvv3v/+t2267TV1dXers7NTevXt12WWXuVctAMB5qqS8vFxNTU2qrKyUZVkKhUKqr69XQUGBysrKdPPNN2vu3LnKy8vTzTffrIkTJw5F3QCQtRyDOzc3V7W1tT3OFRYW2o8XLVqkRYsWDX5lAICUuAEHAAzj+rauQF9YLgj0Hx03ABiG4AYAwxDcAGAYghsADENwA4BhWFWCX6Q/q0JYQQIMDjpuADAMwQ0AhmGqBBmDqRTg56HjBgDDENwAYBimSpCRmDYB0qPjBgDDENwAYBiCGwAMQ3ADgGEIbgAwDMENAIYhuAHAMAQ3ABiG4AYAw3DnJAZd97seAQw+Om4AMAzBDQCGIbgBwDAENwAYhuAGAMMQ3ABgGIIbAAxDcAOAYQhuADCMY3Ank0nV1NRo3rx5qq6uVmtra8rnLFq0SBs3bnSlSADAaY7B3djYqHg8roaGBi1dulR1dXW9nrN69WodP37clQIBAD05Bndzc7OCwaAkqaSkRJFIpMf4li1blJOTYz8HAOAux02motGoAoGAfezz+ZRIJOT3+/Xll1/qgw8+0PPPP68XX3zR1UKRwr6Pex5fOsubOv4fm0sBQ8MxuAOBgGKxmH2cTCbl95/8snfffVeHDx/W7bffrgMHDigvL08XXnihrr76avcqBoAs5xjcpaWl2rZtmyoqKhQOh1VUVGSPLVu2zH78wgsv6LzzziO0AcBljsFdXl6upqYmVVZWyrIshUIh1dfXq6CgQGVlZUNRIwCgG8fgzs3NVW1tbY9zhYWFvZ537733Dl5VAIC0+ASc4aT7xUqPL1QCcA93TiLjNX5xmBUrQDcENwAYhuAGAMMQ3ABgGIIbAAxDcAOAYVgOmIncWNaXYfuaABg4Om4AMAzBDQCGYaok0zHFAeAMdNwAYBiCGwAMQ3ADgGGY44Yx0m00df3kXw9xJYC3CO5sxRawgLGYKgEAwxDcAGAYghsADENwA4BhCG4AMAzBDQCGIbgBwDAENwAYhhtwhqszdxXsz3O5IQfIaHTcAGAYghsADENwA4BhCG4AMAzBDQCGIbgBwDAENwAYhuAGAMM4BncymVRNTY3mzZun6upqtba29hjfsGGDZs+erTlz5uijjz5yrVAAwEmOd042NjYqHo+roaFB4XBYdXV1Wrt2rSTpyJEj2rhxo9555x11dHToxhtv1KxZs5STk+N64VmrP3dEAhiWHDvu5uZmBYNBSVJJSYkikYg9NnbsWL377rvKy8tTe3u7RowYQWgDgMscO+5oNKpAIGAf+3w+JRIJ+f0nv9Tv9+uNN97QCy+8oOrqavcqRcZI92nrQ/09gGzl2HEHAgHFYjH7OJlM2qF9SlVVlXbt2qU9e/bob3/72+BXCQCwOQZ3aWmpdu7cKUkKh8MqKiqyx7755hstWbJElmUpLy9P+fn5ys1loQoAuMlxqqS8vFxNTU2qrKyUZVkKhUKqr69XQUGBysrKNGnSJM2bN085OTkKBoO64oorhqJuAMhajsGdm5ur2traHucKCwvtx0uWLNGSJUsGvzIAQErMawCAYQhuADAMwQ0AhiG4AcAwBDcAGIZPec8E7D8CoB/ouAHAMAQ3ABiG4AYAwxDcAGAYghsADMOqEvxs7KENZAY6bgAwDMENAIYhuGG8xi8OM42DrMIcN3rrfifnpbO8qwNASnTcAGAYghsADENwA4BhCG4AMAzBDQCGIbgBwDAsB/SCSR+c0K3W8w4eU/sF13lYDACJjhsAjEPHjT6F9x/zugQAZyC4AcNs37895fnfXPyboSwDHiK4gSFA2GIwMccNAIYhuAHAMEyVuOXMJX/ssgdgkBDcGDa678l9/eRfe1gJ4C6mSgDAMHTc6JfzDm61H3MXJeANOm4AMIxjx51MJrVixQrt27dP+fn5WrlypSZMmGCPv/baa/rwww8lSddcc42WLFniXrXAALGOGsOJY8fd2NioeDyuhoYGLV26VHV1dfbY/v379f7772vTpk3avHmzdu/erZaWFlcLBoBs59hxNzc3KxgMSpJKSkoUiUTssfHjx2v9+vXy+XySpEQioREjRrhUKtzUfU+SkovHeFYHAGeOHXc0GlUgELCPfT6fEomEJCkvL09jx46VZVlatWqVJk+erEsuucS9agEAzsEdCAQUi8Xs42QyKb//dKPe0dGhhx56SLFYTI8//rg7VQIAbI7BXVpaqp07d0qSwuGwioqK7DHLsnT33Xfr0ksvVW1trT1lAgBwj+Mcd3l5uZqamlRZWSnLshQKhVRfX6+CggIlk0n94x//UDwe165duyRJDz74oKZOnep64QCQrRyDOzc3V7W1tT3OFRYW2o//+c9/Dn5VgIHSLTkEBhs34ACAYQhuADAMwQ0AhiG4AcAw7A6IAeu+U6CUWbsFsjc3hjM6bgAwDB03kIFYWoi+ENxAPxGq8BpTJQBgGDpuYJjr6x0CHyRhJjpuADAMHTfgIebLMRAE91DZ97HXFQAYJpgqAQDD0HEDaTCNgUxFxw0AhqHjRlajq4aJCG70Et5/zOsSBlXjF4cVOXZMklRy8RhPazFFuj9orPvODEyVAIBhCG4AMAxTJcAwwXx99qDjBgDD0HFjWIkc+x+vSzAKXbqZ6LgBwDAENwAYhuAGAMMQ3ABgGIIbAAxDcAOAYQhuADAM67iz3HDbUArIBgQ3XHPewa324/YLrvOwEmB4IbhhJO6QRDZznONOJpOqqanRvHnzVF1drdbW1l7POXLkiGbOnKmOjg5XigQAnOYY3I2NjYrH42poaNDSpUtVV1fXY3zXrl2666671NbW5lqRAIDTHKdKmpubFQwGJUklJSWKRCI9xnNzc1VfX6/Zs2e7UyGAjNHXplR8Os7QcQzuaDSqQCBgH/t8PiUSCfn9J790xowZ7lUHDLJTq2j4CDOYzHGqJBAIKBaL2cfJZNIObQDA0HNM4NLSUm3btk0VFRUKh8MqKioairrMsO/jnseXzvKmjgzRffkfAPc4Bnd5ebmamppUWVkpy7IUCoVUX1+vgoIClZWVDUWNAIBuHIM7NzdXtbW1Pc4VFhb2et7WrXRbADAUmKweTGdOnQCAC9hkapgL7z/GfiTAMEPHjSHR14XLdPuYcFu7WdKt8WZ99+Cj4wYAw9BxZyGmTgCzEdwAPMHt8wNHcGPI7P3pf1OeLxB7dQP9QXDDc1yEBPqH4EZWSjfPz+ZTMAGrSgDAMHTc/cXdkQA8RscNAIYhuAHAMEyVAHBVX+u1MTAEd5bgbklg+GCqBAAMQ8eNQZXu7kgAg4fgxoAQ0IB3mCoBAMNkZ8fNp7NnlHO+/zzt2PH/mjKElQBmoOMGAMNkZ8eNHvqary4dNXEIKwHwcxDcZ2IaBUCGI7gBGIMPJD6J4HbCboDAkOM2+b5xcRIADENwA4BhmCqRhsV0SPdNpPj4rYFL9Xvkd2u24TgvTnADMF62zYkzVQIAhqHjBtJgD3NkqswP7u7zz043w/TnuTBCX/uYSOxlguyU+cGNftt+9Av78Tc/xezHA7l9ne1bAWdDfQHUMbiTyaRWrFihffv2KT8/XytXrtSECRPs8c2bN2vTpk3y+/1avHixrr32WlcKHe5SvS1Pt4LB6S38N22xPscxcKwwQSZwDO7GxkbF43E1NDQoHA6rrq5Oa9eulSS1tbXp9ddf11tvvaWOjg7Nnz9fM2bMUH5+vjvV9mfZXgbuOdK9Ez7TGF0whJUMH2wJi4EyeZmgY3A3NzcrGAxKkkpKShSJROyxzz//XFOnTlV+fr7y8/NVUFCglpYWTZly+j9MV1eXJOnQoUMDq/Dw9wP7ujON/nbwv2c/tf/nh7Rjnx35rNe5r346K+Vz9x/5UZeN/O+03+s/J6Ipz29T79cYzn6wjrj6/dvzEq5+f3jjr4f+mnZs+vnTU55vP9Se8vy3Od+mPO/kVF6eys8zOQZ3NBpVIBCwj30+nxKJhPx+v6LRqM4++2x7bPTo0YpGe4ZGW1ubJGnBggX9r35QPenx62Poved1AcAv0tbW1mNq+hTH4A4EAorFTs+ZJpNJ+f3+lGOxWKxHkEtScXGxNmzYoHHjxsnn8w34BwCAbNHV1aW2tjYVFxenHHcM7tLSUm3btk0VFRUKh8MqKiqyx6ZMmaLVq1ero6ND8XhcX3/9dY9xSRo5cqQuv/zyX/hjAEB2SdVpn5JjWZbV1xefWlXy5ZdfyrIshUIh7dy5UwUFBSorK9PmzZvV0NAgy7L0+9//XjNnzhz0HwAAcJpjcLutq6tLTz/9tCKRiOLxuO69996MXlL49ddfa+7cufr00081YsQIr8vp5YcfftDDDz+saDSqzs5OPfLII5o6darXZdmclpdmis7OTi1fvlwHDhxQPB7X4sWLVVZW5nVZaX3//fe69dZb9eqrr6qwsNDrclJ66aWXtHXrVnV2duq2227Tb3/7W69L6uXU/5kDBw4oNzdXTz75ZEb+Pj3fq+S9995TIpHQpk2btHbtWrW2tnpdUlrRaFSrVq1yb7njIKivr9f06dP1xhtv6Omnn1Ztba3XJfXQfXnp0qVLVVdX53VJKb3//vsaM2aM/vKXv2j9+vV68snMvbjd2dmpmpoajRw50utS0vr73/+uzz77TBs3btTrr78+8FVmLtuxY4edR/fcc49Wr17tdUkpeX7n5O7duzVx4kT97ne/k2VZeuyxx7wuKaVTtT344IO6++67vS4nrTvuuMP+w9LV1ZVx7wr6Wl6aSW644QZ72s+yrIy+sL5q1SpVVlbq5Zdf9rqUtHbv3q2ioiLdc889ikajWrZsmdclpXTJJZeoq6tLyWRS0WjUXoiRaYa0qjfffFN/+tOfepw799xzNWLECL300kvas2ePHn30UW3YsGEoy+olVZ0XXHCBKioqNGnSJI+q6i1VnaFQSFOmTFFbW5sefvhhLV++3KPqUutreWkmGT16tKST9f7hD3/Q/fff721Babz99tsaO3asgsFgRgf30aNHdfDgQa1bt07ffvutFi9erC1btignJ8fr0no466yzdODAAc2aNUtHjx7VunXrvC4pNctj999/v7Vlyxb7+KqrrvKwmvSuv/56q6qqyqqqqrKKi4ut+fPne11SWi0tLVZFRYW1fft2r0vpJRQKWR9++KF9HAwGPaymbwcPHrRuueUW68033/S6lLTmz59vLViwwKqqqrKmTZtmzZ492/ruu++8LquXZ555xnrllVfs45tuuslqb2/3sKLUQqGQ9eyzz1qWdfLfv7y83Dpx4oTHVfXmeZszbdo07dixQzNnzlRLS4vOP/98r0tK6ZNPPrEfX3fddXr11Vc9rCa9r776Svfdd59Wr16dUe8OTulreWkmaW9v11133aWamhpdeeWVXpeTVvd3p9XV1VqxYoXGjRvnYUWpTZs2TX/+859155136rvvvtNPP/2kMWPGeF1WL+ecc47y8vIkSb/61a+USCTS3r3oJc+De+7cuXr88cc1d+5cWZalJ554wuuSjPbcc88pHo/rqaeeknTyJqlTe8tkgvLycjU1NamystJeXpqJ1q1bp+PHj2vNmjVas2aNJOmPf/xjRl8AzGTXXnut9uzZozlz5siyLNXU1GTkdYM77rhDy5cv1/z589XZ2akHHnhAZ52VeusJL3m+HBAA0D+eLwcEAPQPwQ0AhiG4AcAwBDcAGIbgBgDDENwAYBiCGwAMQ3ADgGH+D1u83Fv+SXrbAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "x1 = rng.normal(0, 0.8, 1000)\n",
    "x2 = rng.normal(-2, 1, 1000)\n",
    "x3 = rng.normal(3, 2, 1000)\n",
    "\n",
    "kwargs = dict(histtype='stepfilled', alpha=0.3, density=True, bins=40)\n",
    "\n",
    "plt.hist(x1, **kwargs)\n",
    "plt.hist(x2, **kwargs)\n",
    "plt.hist(x3, **kwargs);"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If you are interested in computing, but not displaying, the histogram (that is, counting the number of points in a given bin), you can use the `np.histogram` function:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[ 23 241 491 224  21]\n"
     ]
    }
   ],
   "source": [
    "counts, bin_edges = np.histogram(data, bins=5)\n",
    "print(counts)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Two-Dimensional Histograms and Binnings\n",
    "\n",
    "Just as we create histograms in one dimension by dividing the number line into bins, we can also create histograms in two dimensions by dividing points among two-dimensional bins.\n",
    "We'll take a brief look at several ways to do this here.\n",
    "We'll start by defining some data—an `x` and `y` array drawn from a multivariate Gaussian distribution:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [],
   "source": [
    "mean = [0, 0]\n",
    "cov = [[1, 1], [1, 2]]\n",
    "x, y = rng.multivariate_normal(mean, cov, 10000).T"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### plt.hist2d: Two-dimensional histogram\n",
    "\n",
    "One straightforward way to plot a two-dimensional histogram is to use Matplotlib's `plt.hist2d` function (see the following figure):"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWYAAAD3CAYAAAAua/5EAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAgVklEQVR4nO3dfVRUdf4H8PcwOEAzPKyGmSElpqlLnFRWbEszyx2XPa1uq4gYxpr7oK2GehRFHM1U5LiSra3gsicssBQXtyh394gPyVkfqDQ1yYekJ58oNE1mqmEY7u8PD/MDU+6Xy52Zy+X9OodzmOF77/1e0Ddfvvdzv9cgSZIEIiLSjAB/d4CIiFpiMBMRaQyDmYhIYxjMREQaw2AmItIYBjMRkcYwmImIABw9ehSpqakAgMuXL2P69OmYPHkykpOT8eWXXwIASkpK8OSTTyIpKQl79uzxWl8CvbZnIqIOoqCgAGVlZQgJCQEArF69Gk888QQSExNx8OBBfPrppwgJCUFRURFKS0vhdDqRkpKChx56CCaTSfX+eCWYf/jhBxw/fhyRkZEwGo3eOAQR6YTb7UZtbS1iY2MRHByseD9Xr16F3W4XamuxWBAREeF5HR0djXXr1mH+/PkAgMOHD+O+++5DWloa7rrrLixatAgHDhzAoEGDYDKZYDKZEB0djZMnTyIuLk5xn2/FK8F8/PhxTJ482Ru7JiKd2rRpE+Lj4xVte/XqVcTHxwsPBMPDw7Fjxw5POFutVpw7d87z9fPnzyMsLAwbN27Eyy+/jIKCAtxzzz0IDQ31tDGbzcK/CNrKK8EcGRkJ4Po3ukePHt44BBHpRE1NDSZPnuzJDSXsdjuMRiPOnj2LhoaGVtsGBgZ6tmk+am4uIiICo0aNAgCMGjUKL774ImJjY+FwODxtHA5Hi6BWk1eCuem3Vo8ePRAVFeWNQxCRzqgx7dnY2IjGxkbZNnKGDBmCvXv3Yty4cXj//fdx7733Ii4uDmvXroXT6UR9fT2qq6vRr1+/dvf5Znjxj4h0IyAgAAEBrRebyX0dADIyMpCVlYXNmzfDYrFgzZo1CA8PR2pqKlJSUiBJEmbPno2goCC1ut4Cg5mIdMNgMMBgMMi2uZmoqCiUlJQAAO666y4UFhb+qE1SUhKSkpLa31EZDGYi0hW5YO4IGMxEpBvtGTFrCYOZiHSDwUxEpDEilR0d4aY3BjMR6QZHzEREGsNgJiLSoI4QvHIYzESkGxwxExFpjNFolA1ekTv//I3BTOQjkiTJtpELFTX2AYitFyGyH62NPjliJiLSGAYzEZHGMJiJiDSoIwSvHAYzEelGR5wXvxkGMxHphtFolL1AymAmIvIhjpiJiDSGwUxEqhOpL5YjEjwd4SYLJRjMREQa1BGCV067fm1evnwZjzzyCKqrq9XqDxGRYkajUejjZo4ePYrU1NQW77399tuYOHGi53VJSQmefPJJJCUlYc+ePV47D8UjZpfLBZvNhuDgYDX7Q0SkmNLRckFBAcrKyhASEuJ57+OPP8Y///lPT5VHbW0tioqKUFpaCqfTiZSUFDz00EMwmUyq9L05xSPmnJwcJCcno3v37mr2h4hIsaY7/+Q+bhQdHY1169Z5Xl+5cgW5ubnIzMz0vHfs2DEMGjQIJpMJoaGhiI6OxsmTJ71yHoqCedu2bejatSuGDx+udn+IiBRTGsxWqxWBgdcnENxuNxYtWoSFCxfCbDZ72tjtdoSGhnpem81m2O12r5yHoqmM0tJSGAwGHDhwACdOnEBGRgby8vIQGRmpdv+IiISJrJUhdwNKVVUVvvjiCyxduhROpxNnzpzBihUrMGzYMDgcDk87h8PRIqjVpCiYN23a5Pk8NTUVS5cuZShTh+R2u2XbiCy1KVLmpkaJmlrLfqq1H61VQIgEM9D6+cfFxWH79u0AgHPnzmHOnDlYtGgRamtrsXbtWjidTtTX16O6uhr9+vVTre/NsVyOiHRDZKF8SZLQ0NDQ5n1HRkYiNTUVKSkpkCQJs2fPRlBQkNKutqrdwVxUVKRGP4iI2k10xHwzUVFRKCkpafW9pKQkJCUltauPIjhiJiLdMBgMslNGatxd6W0MZiLSDS6UT0SkMQxmIiKNCQgI0MUCTQxm6pBE5glFSsJcLpdsm6YbD1ojMgoTqQSQO5bIOYnoCKNGJYxGo2wwd4RzZzATka50hOCVw2AmIt0QmcpQ668Ob2IwE5Fu8OIfEZHGMJiJiDSmtYXwOxIGMxHpBu/8I/IjJYvQ3IxIKZzT6ZRtI1I7K9JnudBQ60GrIufdEXEqg4hIYxjMREQaI1Iu1xHuDGQwE5FucMRMRKQxAQEBslUZvPhHRORDHDETEWkM55iJiDSGI2aiG4gsDqPWUpxq7UfkKdkic5L19fWybW677TbZNnJERnsibdRaNlVrd9npZcSs/R4SEbVB06j5Vh+3cvToUaSmpgIATpw4gZSUFKSmpuKZZ57BpUuXAAAlJSV48sknkZSUhD179njtHDhiJiLdEFkr42ZfLygoQFlZGUJCQgAAK1aswOLFizFgwABs3rwZBQUFmDZtGoqKilBaWgqn04mUlBQ89NBDMJlMqp8HR8xEpBtyo+VbjZqjo6Oxbt06z+vc3FwMGDAAwPXprqCgIBw7dgyDBg2CyWRCaGgooqOjcfLkSa+cB4OZiHRDaTBbrdYW64d0794dAHD48GEUFxcjLS0NdrsdoaGhnjZmsxl2u90r58GpDCLSDTWrMv79738jLy8Pf//739G1a1dYLBY4HA7P1x0OR4ugVhNHzESkG0pHzDd66623UFxcjKKiIvTq1QsAEBcXh0OHDsHpdKKurg7V1dXo16+fV86DI2YSIlI6JbKspch+1CpPu3btmir9+eGHH2TbhIWFybYRWT5U7nsoUnIn0t+mi1yt6QhlZTcSWY9ZLpjdbjdWrFiBO++8EzNnzgQA/OxnP8OsWbOQmpqKlJQUSJKE2bNnIygoSLW+N8dgJiLdaM9URlRUFEpKSgAA77333k3bJCUlISkpqX2dFMBgJiLd0MsNJgxmItIN3pJNRKQxDGYiIo1hMBMRaYwaVRlaoCiYXS4XMjMzcf78edTX12P69Ol47LHH1O4bqUCt8jSRVcRE/sFfuXJFts23334r20atp1J//vnnsm1E1kIQOa/evXvLtlFjtbbg4GDZNh0hnJTo1CPmsrIyREREYPXq1bh69SrGjRvHYCYiv+vUwTxmzBhYrVYA10dkWluTlYg6p04dzGazGQBgt9sxa9YspKenq9knIiJF9BLMiiutL168iClTpmDs2LF44okn1OwTEZEiaq2V4W+KRsyXLl3C1KlTYbPZ8OCDD6rdJyIiRfRSlaFoxJyfn49r165h/fr1SE1NRWpqqtDCKURE3tSpR8xZWVnIyspSuy/kBb5czU3kwaYi/Tl37pxsm4sXL8q2ESlha7pe0ppu3brJtrn99ttl24is6CZHrXUeRH4OHSHAbqSXOWbeYEJEusFgJiLSGAYzEZHGcNlPIiKN4YiZiEiDOkLwymEwE5Fu6GXErP3JFiIiQe2pYz569ChSU1MBAF988QUmTZqElJQULFmyxFNS+vLLL2P8+PFITk7GsWPHvHYeHDFrmFytqUgtqkhtsUiNssPhkG0jUjf8ySefyLY5ePCgbJvq6mrZNiJPlL777rtl24h8D2NiYmTbiDy1Ozw8vNWvi1y46ggXt7xF6Yi5oKAAZWVlnlrz7OxspKenIyEhATabDbt27ULPnj3x3nvvYevWrbh48SJmzpyJ0tJSr5xH5/0JEpHuNFVlyH3cKDo6GuvWrfO8rqqqwtChQwEAI0aMwP79+3Ho0CE8/PDDMBgM6NmzJ9xuN7755hvvnIdX9kpE5AdKg9lqtSIw8P8nECRJ8oyszWYz6urqYLfbYbFYPG2a3vcGTmUQka6ocXGveXg7HA6EhYXBYrG0mNJzOBwIDQ295T7efPNNbNiwAfX19Z6g37Vrl9DxGcxEpBtqVWUMHDgQlZWVSEhIQEVFBYYNG4bo6GisXr0azzzzDGpqatDY2IiuXbvech8FBQXIz8/HnXfe2ebzYDATkW6oFcwZGRlYvHgxcnNzERMTA6vVCqPRiPj4eEycOBGNjY2w2Wyt7qNXr15CF5dvhsFMRLrRnluyo6KiUFJSAuD6g3OLi4t/1GbmzJmYOXOmUF+Cg4Mxbdo0DBgwwPPLYM6cOULbMpg7MJHlOkXauFwu2Ta1tbWybd577z3ZNmfOnJFt89FHH8m2EXnOpMh52e122TZOp1O2TfMLR7cisuynGou863VJTxFausHkkUceUbwtg5mIdEMLwfzRRx/h/vvvR2RkpOJ9MJiJSDe0EMwHDhzA/fffj+3bt//oaw8//LDQPhjMRKQr/p6m+cMf/gDg+t2DJ0+exOeff46+ffuiT58+wvtgMBORbmhhxNwkLy8PFRUVuP/++7Fx40aMGTMGaWlpQtsymIlIN4xGo+yFYZELx2p499138cYbbyAgIAANDQ1ISUkRDmbekk1EuqGlp2R369YN33//PYDrFUKt3YxyI46Y/USkpElOQ0ODbJtvv/1Wts3nn38u2+bChQuybT744APZNiIrrL399tuybSIiImTbTJw4UbbNHXfcIdvmnnvukW0THBws26ZLly6ybYKCglr9uq9Gex2VFqYyJk6cCIPBgMuXL8NqteK+++5DdXW10L/ZJgxmItINLQRzbm5uu/fBYCYiXfF3VcZdd93V7n0wmIlIN7QwYlYDg5mIdKM9a2VoCYOZiHRDSyPmffv2obCwsMWj21577TWhbRnMRKQbWhoxZ2dnIzMzEz169GjztgxmPxH5rS33ENCmGsnWiJTUiaz49tlnn8m2OXHihGybs2fPyrYRKWETaSNSWiayiLlIuVxrT7JoYjKZZNvIlVGKrBbYEf5U9xYtjZjvvPNO/PznP1e0LYOZiHRFKxf3unXrBpvNhoEDB3r6JFJbDzCYiUhHtDSVERUVBQC4dOlSm7dlMBORbmhhKqOmpgY9evTAr371K8X7UBzMjY2NWLp0KU6dOgWTyYTly5crfr4VEZEatBDMhYWFWLhwIWw2m+dYTU/J9npVxs6dO1FfX48tW7bgyJEjWLVqFfLy8pTujoio3ZQGs8vlwoIFC3D+/HkEBATghRdeQGBgIBYsWACDwYC+fftiyZIlQtMgCxcuBAAUFRUpOwm0I5gPHTqE4cOHAwAeeOABHD9+XHEniIjUoDSY9+7di4aGBmzevBn79u3D2rVr4XK5kJ6ejoSEBNhsNuzatQujR4/2VtdbUDwLbrfbYbFYPK+NRqNQaRYRkbcoXfazd+/ecLvdaGxshN1uR2BgIKqqqjB06FAAwIgRI7B//36fnYfiEbPFYoHD4fC8bmxsFHpSMF3X/G6gW5H7Rffdd9/J7qOmpka2jUiN8m233SbbRuSJ03LLWgJiT7fu1q2bbJuBAwfKtunbt69sG5FFaUT+xBX5/yE32uvMNcoiAgICZOvXb/Y9vO2223D+/Hn88pe/xJUrV5Cfn4/333/f8/Mwm82oq6trc3/sdjsMBgPKy8vx6KOPIjw8XGg7xUk6ePBg7NmzB4mJiThy5Aj69eundFdERKpQOpWxceNGPPzww5g7dy4uXryIp59+usUAweFwICwsrE19mT17NkaOHIkPP/wQjY2NKC8vx9/+9jehbRX/+h09ejRMJhOSk5ORnZ3tmfAmIvIXpVMZYWFhnrs3w8PD0dDQgIEDB6KyshIAUFFRgfj4+Db15euvv8bYsWNRXV2NZcuWtZhhkKN4xBwQEIBly5Yp3ZyISHVKbzBJS0tDZmYmUlJS4HK5MHv2bMTGxmLx4sXIzc1FTEwMrFZrm/ricrmwY8cO3Hvvvfjmm298E8xERFqjdCrDbDbjpZde+tH7xcXFivvy+9//Hu+88w4WLlyIoqIiPPvss8Lb8koCEemKFh7ECgBXr17FSy+9hB49euC5557Dl19+KbwtR8xEpBtaWCvjnXfewe7du1FZWYmDBw8CuF61dvr0aUyZMkVoHwzmNhJZdlGk3EtkSUqn09nq1y9fviy7D5GnW4vUn5eXl8u2ERmNiJTLJScny7YRKTsaNWqUbBuRMkCR5TpFnoDNJ1x7nxaCefjw4YiMjMTVq1c9q8kFBASgV69ewvtgMBORbmhhrYzw8HAkJCQgISEBly9f9gyw5NZXb47BTES6oYVgbvL8889j79696N69u2cRo82bNwtty2AmIt3QUjAfPXoUO3fuVDR1wmAmIt1Qeku2N9x9991wOp0ICQlp87YMZiLSDS2NmC9evIhHH33Us049pzKIqFPSUjCvWbNG8bYM5mZESuFEfqgiZVEit2fKrdYmstqVyJOrRcrlgoODZdv0799fto1I6VnXrl1l2/Tp00e2jciqcGqVwnHVN23QUjD/61//+tF7f/7zn4W2ZTATkW5ooY65ye233w7g+mOlPv74Y6GBXxMGMxHpii9vu27NjTdKTZs2TXhbBjMR6YaWRszNH0BRW1srdBduEwYzEemGloLZZrN5Pg8KCkJGRobwtgxmItINLV38KyoqwpUrV3D27FlERUUJXdRuwkvJRKQbSp9g4g3/+c9/kJycjPz8fEycOBFvvfWW8LYcMTcj8ieOSGmZyH4kSWr3sURWRhN5aOkPP/wg26bpacGtaboK3Zqf/OQnsm1EVuG65557ZNuILBoj8rMSWS1QpJyQvE9LI+aNGzdi27ZtMJvNsNvtePrppzF27FihbRnMRKQrWqnKMBgMMJvNAACLxSK05G0TBjMR6YbRaJS9wctX62L36tULq1atQnx8PD744ANER0cLb8tgJiLdaM9UxoYNG7B79264XC5MmjQJQ4cOxYIFC2AwGNC3b18sWbKkTRUd2dnZ2LJlC/bv348+ffpg7ty5wtvy4h8R6YbSi3+VlZX48MMP8cYbb6CoqAg1NTXIzs5Geno6Xn/9dUiShF27drWpL1VVVXC73bDZbDh8+DA++eQT4W0ZzESkG011zHIfN/rf//6Hfv364dlnn8Wf/vQnjBw5ElVVVZ6L3iNGjMD+/fvb1Jdly5Zh5MiRAID09HSsWLFCeFtOZRCRrii5+HflyhVcuHAB+fn5OHfuHKZPn+556ggAmM1moUXDmuvSpYtnXrlXr15tmgZhMBORbiidY46IiEBMTAxMJhNiYmIQFBSEmpoaz9cdDgfCwsLa1JeePXsiNzcXDzzwAI4dO4bu3bsLb8tgbiOROmYRIleG5VajEqnT7d27t2wbuadxA0BsbKxsGxEiddUi/4BDQ0Nl2wQGyv/z5pOr9UVpMA8ZMgSvvfYafve73+Hrr7/G999/jwcffBCVlZVISEhARUUFhg0b1qa+ZGdn44033sDevXvRp08fzJgxQ3hbBjMR6YbSYH700Ufx/vvvY/z48ZAkCTabDVFRUVi8eDFyc3MRExMDq9Xapr4EBQUhLS2tTds0YTATkW60p1xu/vz5P3qvuLhYlX61FYOZiHRDS7dktweDmYh0g8FMRKQxDGYiIg3qCMErR1Ew19XVYd68ebDb7XC5XFiwYAEGDRqkdt+EiSyhKfLDEik/EymvEimp++6772TbyD3BuWnlqtZ88803sm1ElvQUeZq0SH8iIiJk24gsZyryMxf5WenhPzH9v049Yi4sLMSwYcOQlpaGTz/9FHPnzr3po7qJiHypUwdzWlqaZwTldrvbtM4oEZG3dJpg3rp1K1599dUW761cuRJxcXGora3FvHnzkJmZ6bUOEhGJ6jTBPGHCBEyYMOFH7586dQpz5szB/PnzheYoiYi8TUtPyW4PRVMZZ86cwXPPPYe1a9eif//+aveJiKhTUxTMa9asQX19vWd9UYvFgry8PFU7RkTUVp1mKuNmtBbCan2jRf7EUevJyyLlZ3JPXpZbfQ4A4uLiZNuInFN4eLhsG5ESNpHvjciqcCJtqHPqCMErh/+6iUg3OvWImYhIizr1xT8iIi3iiJmISGMYzEREGqOXYNb+ZAsRUSfDEXMzIuVnal04CAkJkW3jcDha/brFYpHdh8iKby6XS7aNiC5dusi2EVnxTa3vcUcYGZG6OGImItIYg8Hgqcy41cetgvny5ct45JFHUF1djS+++AKTJk1CSkoKlixZIjRoUxODmYh0o2nELPdxI5fLBZvN5rmpKzs7G+np6Xj99dchSRJ27drl0/NgMBORbigN5pycHCQnJ6N79+4AgKqqKs/ibCNGjMD+/ft9eh4MZiLSDSXBvG3bNnTt2hXDhw/3vCdJkqed2WxGXV2dT8+DF/+IqFMrLS2FwWDAgQMHcOLECWRkZLR4JJvD4UBYWJhP+8RgJiLdUHJL9qZNmzyfp6amYunSpVi9ejUqKyuRkJCAiooKDBs2zCv9vWUffXo0IqIOICMjA+vWrcPEiRPhcrlgtVp9enyOmJsRKYlR62nbIseSe5aiWk+BFllCU6RGWaT+WK0a0o5Qi0q+19465qKiIs/nxcXFqvWrrRjMRKQrevilzWAmIt3gnX9EROQVHDETkW7oZaF87feQiKiT4YiZiHRDL3PMDOZmRErCRErhRH7wIm3knpKtFpFyOZa5UUfAYCYi0hi9BDPnmImINIYjZiLSjaaF8uXaaB2DmYh0g1MZRETkFRwxE5Fu6GXE3GmCWa0nYKtVCie3cpzIfiRJkt2HCLXOicjfGMxERBrDYCYi0hgGMxGRxuglmNtVlVFdXY0hQ4bA6XSq1R8iok5P8YjZbrcjJycHJpNJzf4QESmmlxGzomCWJAmLFy/GnDlzMGPGDLX7RESkWFuD1+VyITMzE+fPn0d9fT2mT5+Oe++9FwsWLIDBYEDfvn2xZMkSn67jLBvMW7duxauvvtrivZ49eyIxMRH9+/f3WsfUJvJNFSk/E/mhi6xSpwaR/naERcGJ1KJkxFxWVoaIiAisXr0aV69exbhx49C/f3+kp6cjISEBNpsNu3btwujRo73Z9RZkg3nChAmYMGFCi/dGjx6N0tJSlJaWora2FlOnTsWmTZu81kkiIhFKgnnMmDGwWq0Arg92jEYjqqqqMHToUADAiBEjsG/fPm0F882Ul5d7Ph81ahReeeUV1TpERORLZrMZwPXrZrNmzUJ6ejpycnI8AW42m1FXV+fTPvHvXCLSjaYRs9zHjS5evIgpU6Zg7NixeOKJJ1pMATocDoSFhfnyNNofzLt37xa6vZiISIsuXbqEqVOnYt68eRg/fjwAYODAgaisrAQAVFRUID4+3qd94g0mRKQbSuaY8/Pzce3aNaxfvx7r168HACxatAjLly9Hbm4uYmJiPHPQvsJgJiLdUBLMWVlZyMrK+lG74uJiVfvWFpxjJiLSGI6Ym+kIdwQ1xxplopb0cucf/2cTEWkMR8xEpBt6GTEzmIlIVzpC8MrhVAYRkcZwxExEuqGXqQyOmImINIbBTESkMZzKICLd0MtUBoOZiHRDL8HMqQwiIo3hiJmIdIMjZiIi8gqOmIlINzhiJiIir+CImYh0Qy8jZgYzEemGXoKZUxlERBrDETMR6UpHGBHLYTATUafX2NiIpUuX4tSpUzCZTFi+fDnuvvtuv/XHK8HsdrsBADU1Nd7YPRHpSFNONOVGe3z11VeyI+avvvrqR+/t3LkT9fX12LJlC44cOYJVq1YhLy+v3f1RyivBXFtbCwCYPHmyN3ZPRDpUW1ureJRqsVgQHh4unDnh4eGwWCye14cOHcLw4cMBAA888ACOHz+uqB9q8Uowx8bGYtOmTYiMjITRaPTGIYhIJ9xuN2praxEbG6t4HxEREdixYwfsdrtQe4vFgoiICM9ru93eIqiNRiMaGhoQGOif2V6vHDU4OBjx8fHe2DUR6ZAa87kREREtwrYtLBYLHA6H53VjY6PfQhlguRwREQYPHoyKigoAwJEjR9CvXz+/9scgSZLk1x4QEflZU1XG6dOnIUkSVq5ciT59+vitP34J5urqaiQlJWH//v0ICgry9eHx3XffYe7cubh27Rq6dOmCnJwc3HHHHT7vR11dHebNmwe73Q6Xy4UFCxZg0KBBPu9Hc+Xl5fjvf/+LNWvW+PS4WipXOnr0KP7yl7+gqKjIL8cHAJfLhczMTJw/fx719fWYPn06HnvsMZ/3w+12IysrC5999hkMBgOef/55v48mOwOfT2XY7Xbk5OTAZDL5+tAeJSUl+OlPf4pNmzbh17/+NQoKCvzSj8LCQgwbNgzFxcXIzs7GsmXL/NKPJsuXL8eaNWvQ2Njo82M3L1eaO3cuVq1a5fM+AEBBQQGysrLgdDr9cvwmZWVliIiIwOuvv45//OMfeOGFF/zSjz179gAANm/ejPT0dLz44ot+6Udn49PZbUmSsHjxYsyZMwczZszw5aFbSEtL89RMXrhwAWFhYX7rR9MvKLfb7Ze/HpobPHgwHn/8cWzZssXnx9ZKuVJ0dDTWrVuH+fPn++X4TcaMGQOr1Qrg+v8bf1U3Pf744xg5ciQA//5f6Wy8Fsxbt27Fq6++2uK9nj17IjExEf379/fWYYX6sXLlSsTFxWHKlCk4ffo0CgsL/dqP2tpazJs3D5mZmV7vR2t9SUxMRGVlpU/6cCOtlCtZrVacO3fOp8e8GbPZDOD692XWrFlIT0/3W18CAwORkZGB8vJy/PWvf/VbPzoVyYcef/xx6amnnpKeeuopKTY2VkpJSfHl4W/qzJkz0mOPPea34588eVJKTEyU3n33Xb/1obmDBw9K6enpPj/uypUrpe3bt3teDx8+3Od9aHL27FlpwoQJfjt+kwsXLki/+c1vpK1bt/q7K5IkSdLXX38tjRw5UnI4HP7uiu75dDhSXl7u+XzUqFF45ZVXfHl4jw0bNuCOO+7AuHHjYDab/fZn4pkzZ/Dcc89h7dq1Pv0rQosGDx6MPXv2IDExURPlSv526dIlTJ06FTabDQ8++KDf+vHmm2/iq6++wh//+EeEhITAYDAgIIBVtt7WKRcx+u1vf4uMjAyUlpbC7XZj5cqVfunHmjVrUF9fjxUrVgC4XuTuz/vz/Wn06NHYt28fkpOTPeVKnVl+fj6uXbuG9evXY/369QCuX5gMDg72aT9+8YtfYOHChZg8eTIaGhqQmZnp8z50RqxjJiLSGP5NQkSkMQxmIiKNYTATEWkMg5mISGMYzEREGsNgJiLSGAYzEZHG/B8mza1L3E/2dgAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist2d(x, y, bins=30)\n",
    "cb = plt.colorbar()\n",
    "cb.set_label('counts in bin')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Just like `plt.hist`, `plt.hist2d` has a number of extra options to fine-tune the plot and the binning, which are nicely outlined in the function docstring.\n",
    "Further, just as `plt.hist` has a counterpart in `np.histogram`, `plt.hist2d` has a counterpart in `np.histogram2d`:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(30, 30)\n"
     ]
    }
   ],
   "source": [
    "counts, xedges, yedges = np.histogram2d(x, y, bins=30)\n",
    "print(counts.shape)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "For the generalization of this histogram binning when there are more than two dimensions, see the `np.histogramdd` function."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### plt.hexbin: Hexagonal binnings\n",
    "\n",
    "The two-dimensional histogram creates a tesselation of squares across the axes.\n",
    "Another natural shape for such a tesselation is the regular hexagon.\n",
    "For this purpose, Matplotlib provides the `plt.hexbin` routine, which represents a two-dimensional dataset binned within a grid of hexagons (see the following figure):"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWYAAAD0CAYAAACo/4zqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABA9klEQVR4nO2de5BUVX7Hv/1+Ts97mBkGREBkicHHEkGzoGaLwphsaVmKiDtK3PxjdlXEUpCV0d1E0TJSpkypkYquAlExuLXmWRtxSyIk7K4KykODYI0jj6Hn0TP9ft78MfW7nr59X9PTPdPd8/tUTc1Mn3vuPfd23+89/XsdiyRJEhiGYZiKwTrVA2AYhmHyYWFmGIapMFiYGYZhKgwWZoZhmAqDhZlhGKbCsJdjp4lEAkeOHEFraytsNls5DsEwTI2QzWYRDAZxySWXwO12F72fUCiESCRialu/34+Ghoaij1VuyiLMR44cwR133FGOXTMMU6Ps2rULS5YsKapvKBTCkiVLTE8E6+vr8etf/7pixbkswtza2gpg7EK3t7eX4xAMw9QI586dwx133CHrRjFEIhHYbDb09fUhk8nobmu32+U+ojAfPnwYf/u3f4sdO3bIr/3Lv/wLdu7cibfeegsAsHv3brz55puw2+245557cN111xU9Zt0xlmOn9NRqb29HV1dXOQ7BMEyNUQqzZy6XQy6XM9xGyfbt2/Huu+/C4/HIrx07dgz//M//DMrBCwaD2LFjB/bs2YNkMom1a9fij//4j+F0Oic8biXs/GMYpmawWq2mfpTMnj0bzz//vPz/8PAwtm3bhs2bN8uvffrpp7j88svhdDpRV1eH2bNn4/PPPy/PeZRlrwzDMFOAxWIx9aNk1apVsokjm83ipz/9KR555BH4fD55m0gkgrq6Ovl/n89n2tk4XspiymAYhpkq1IR3PBw9ehS9vb14/PHHkUwm8eWXX+KJJ57AsmXLEI1G5e2i0WieUJcSFmaGYWoGrRmxchs9Fi9ejH/7t38DAHzzzTfYsGEDfvrTnyIYDOK5555DMplEKpXCyZMnsWDBgpKNXYSFmWGYmqEUwqxFa2sruru7sXbtWkiShAceeAAul6uofRnBwswwTM1gJrJDa5uuri7s3r1b97XVq1dj9erVExukCViYGYapGco5Y55MWJgZhqkZWJgZhmEqkGoQXiNYmBmGqRl4xswwDFNh2Gw2Q+FVy/yrNFiYGYapGXjGzDAMU2GwMDMMw1QYLMwMwzAVSDUIrxEszAzD1AxmRLkahJuFmWGYmsFms8mF7bVgYWYYhplEeMbMMAxTYbAwMwxjCH2tLlYMJqP/RI9RSbAwMwyjidLOOV7xK1d/2oeaHbZWBLraxw+wMDNMSTFyPBmJ30T7j2cfeu3VKm6lWGm7EmBhZpgawkh0x7OfahTnahyzGizMDMMUUK0CV63jVsLCzDBMzcDCzDAMU2GYqZVRKnNPOWFhZphJQnTcqdlwxyMY5bYBV7ON2cy4K12cWZgZpoSohaLR/+LvUoieMgRO/K0nPFrhcsp9VSNmCuVLkoRMJjNJIyqOCZXyHxwcxDXXXIOTJ0+WajwMU/WIs7ZcLgdJklTFOpfLye3lHIPaa0bt1Qqdg9FPpVO0MKfTafT09MDtdpdyPAxTM1SCABiJUTWJlRksFgusVqvuj9a5Hj58GN3d3QCA48ePY+3ateju7saPfvQjDAwMAAB2796Nm2++GatXr8ZvfvObsp1H0cL89NNPY82aNWhrayvleBiGYYqm2Bnz9u3b8eijjyKZTAIAnnjiCWzZsgU7duzAypUrsX37dgSDQezYsQNvvvkm/vEf/xHbtm1DKpUqy3kUJczvvPMOmpqasHz58lKPh2EmDTInaJkSJElCNpvV3MZMu/ij1a7XP5fLIZvNFm3uMDrHWqNYYZ49ezaef/55+f9t27bhO9/5DgAgm83C5XLh008/xeWXXw6n04m6ujrMnj0bn3/+eVnOoyhh3rNnDw4cOIDu7m4cP34cGzduRDAYLPXYGKYsKMVK/J9+stmsLIgkjuI22WwWmUxGbqO/lf2Vx9USatHeTMfMZDLycTOZDDKZjGmB1TvHWsbIjEE/SlatWgW7/dtYCLIEfPzxx9i5cyfWrVuHSCSCuro6eRufz4dIJFKW8ygqKmPXrl3y393d3Xj88cfR2tpaskExTLnI5XKabSReen31hI2EeCKV4Ixm8JlMBg6HQ/M4WrNzsZ2oFbuyiM1mUxVeEbPn/e///u948cUX8fLLL6OpqQl+vx/RaFRuj0ajeUJdSiYUlcEwzOSiDIub6H5qkVJEZPzqV7/Czp07sWPHDsyaNQsAsHjxYnz00UdIJpMIh8M4efIkFixYUJZzmHAc844dO0oxDoZhmAmjZaoQMTLnZLNZPPHEE+jo6MC9994LAPijP/oj3Hfffeju7sbatWshSRIeeOABuFyuko1dhBNMGIapGczMirXau7q6sHv3bgDAb3/7W9VtVq9ejdWrV09skCZgYWaqhmop5F6t6cy1wESEuZJgYWYqnomu5iH2yeVyef3MpDCL7XrOM7V9jHfFEdF5p+Xcy+Vy8td1tew9rbFojbMahMosNputJorlszAzFUspV/MQY4Hpt9VqVRVctf56SzGptYuzZvotZp0pQ+e00ra1+ms9YJR/G52XWDejFgSaMv/00Iu8qRRYmJmaRq9YDd2g5VqmSZIkzTRg+spNccpa0HZqYkNirtY2nm8CtZaSzaYMhmEMmcrQNqNKcrUGCzPDMEyFYSZczqi9EmBhZqoSMUtOzVxAWXKiuUJrpqRmYxXLcgJQ3b/4W60/MFaFMZ1Ow+FwFOxDrQ6Gch+is0/LJKJ3DmaplUgSnjEzTJlRs5OqpS2TgJItVs1uSyJnJNCiGIqYtUdrOQtTqRQsFgvsdrtsWzZyKCqPLwq01qxPKzJELxKjGoTKLFar1TAqg51/DFMCxEgEZWEgEaN2Mxj1N5pZGjkD0+n0hL5Ki86+8Qiq+JArVVp3JcIzZoaZZKrhhposir0WtSzKANuYGYZhKg6eMTNMEZgxBVTDjcNUJrUyY678ETI1gZ5jTNzGTLveMcxsV2x/s/s3sjNPpJ0cV+U6x1qgFGU/pxqeMTNlYzz1GrTaxYgCZTSGMmpDGZ0hbkN/049aOrUYHif2ofZMJiM775xOZ8HMS3Q+Umib1hhtNluevVdsV45V2T8ej8PhcMgrbpgRG7308WoRKzOYqZVRDbU0WJiZioVmj1qhZSRiSkGlPhaLRRZArdWRKbROLbyO9p3NZpFOp/P6xONxWK1WuR6vMppDjH8W9yfug9rVwtmU4X3K8Ynx0RSCp8Z0K2bENmaGKTMksHro1cKQJClvZipCr+n1B4BkMqkbW5xKpXRnYMpiQ2rtejZPI9NGJpPJW6uuWKpBrMzAwswwDFNh1Iows/NvmqJVanIqxqDXboRRhEc6ndZ15plZYNXIGTeRTDIz74GRkEzkfayEz0EpMXL8VYsDkGfM0wy9NOBSfWDNOPaUgig6y/TsyiKU4gwgLw2b6mSkUil5O4/HA7fbLe8/nU4jHo/Lfex2u2yrpfZYLCabUux2OxwOR14Woij6drsdbrfbdCiWeA1Ee7jYX/xf7ZqQo0tZ/lMr3VrNASiKsp4tvlrQKpGq3KbSYWGeRownEqJc+9eaodLrZstUKsdIIhWPx5FMJguOG4vF5GgGNdHPZDLIZDKQJAnJZLJgFkztWs66TCaDSCQCh8MBj8ejOka9a0APE4vFApfLJdutxcgNMUJDTYTVCi5pOTzVEB8S1SrQtWLKYGFm8ijnh9bMV+ZizRcklkpRVu47nU7rnmMikZiQeUUpqGr99fZBRXj0HJZG4qPVbva9rVZRBmonwYSFmWFKTLWKWi1QKzPmyn90MAzDmGQizr/Dhw+ju7sbANDb24vbb78da9euxWOPPSabif7+7/8et9xyC9asWYNPP/20bOfBM2bGNFzHgql0ip0xb9++He+++67sH9i6dSvWr1+PpUuXoqenB3v37kVnZyd++9vf4u2338bZs2dx7733Ys+ePWU5D54xM3khU2o2UGWdC73IDqN2vTEA3zoBtcZAjjutdq1kD2pPpVKazjcz52AUXkerphQbfqZ1/uK4jELcjNrLGX431ZDjUu9H7fxnz56N559/Xv7/6NGjuPLKKwEAK1aswIEDB/DRRx/he9/7HiwWCzo7O5HNZjE0NFSW8+AZ8zRCq0aEiNkbUrmdWvibsl1vDJlMRo4WUAsdS6VScpaexWKB0+nMy3hLJBKqWXr0fzQaRSQSkTPx/H4/6urq5P0lEglEo1F5DA6HQ063Vp4D7UNZ74LC6aieBrWL5yyGcykjP2j5Ka3lrMSQQPEaKZ2CaiFy4msUYqdWe6SaHX9A8TPmVatW4ZtvvpH/F78d+nw+hMNhRCIRNDQ0yNvQ601NTaUZvAAL8zRDDLnSw0zonN7MUC/pgvaZzWbz1uUTj00iTLNkZXsymUQymYTFYpHjldUIhUIF4W+SJMk3ms/nQzqd1qxF4fF4NMPeMpkMXC4X7Ha76oMol8upttH+xAgOtUgM8foahdeNN3GCjkmz42pJvDCiVM4/8T2LRqMIBALw+/2IRqN5r9PDvdSwKWOaUoqbsBRfd/UEXE2Ule16ogwgL4lErX8ikTAcgx5GCQ1m+muFx41nH8W+n7UwSxYpVebfokWLcPDgQQDAvn37sGTJElxxxRX48MMPkcvlcObMGeRyubLMloEiZ8zpdBqbN2/G6dOnkUqlcM899+D73/9+qcfGMAwzLko1Y964cSO2bNmCbdu2Ye7cuVi1ahVsNhuWLFmC2267DblcDj09PaUadgFFCfO7776LhoYGPPPMMwiFQrjppptYmKsQM1l2WiaN8WTpae03kUjIpSuVs05qTyaTcLvdBRXUJElCJBLB6OgoAoEA3G53QTslk2iNMZ1OIxwOw+12q5oscrkcwuEwnE6nnM4tkslk0N/fD5fLhebmZtVzSCaTmvWbyRyRyWRU2wm9cxBTto1m3bUyK9ZjIsLc1dWF3bt3AwAuvPBC7Ny5s2Cbe++9F/fee+/EB2pAUcJ8/fXXY9WqVQC+La3IVCfjtTmr9dXbRkkul0MikUA8Hs9z/NlsNjidTlgsFsTjcUSjUbk9lUrBbrfD6/XCZrMhHA5jYGBATqGORqNwOp1oamqCy+WSBZfSi2mcZLJIpVKIRCJyjeV4PI6RkREEAgF4vV65DgaZUdLpNKLRKLxeLzweDzKZDEKhEGKxmCx458+fR0tLC1paWmC1WvP6A2PlQ51Op2wPFp2ddAy73Q6Xy5Un0FrX2Gaz5dX2UDoY1VA6+moRI9MSbVPpFCXMPp8PABCJRHDfffdh/fr1pRwTM0kUI6xa+zDTN51OY2RkRLUtm80iGo1qpkSTGA4PDxeEc9HM9Ny5c3IhIbV9WK1WDA4OFlScI1s2OQqV0Ri0bSwWQygUKrBr03iCwSDi8Tja2tpUb/5UKoV0Og2n06l6DcgR6vV6C/qL/9P49ATYTEhcNQjUeJn2mX9nz57FnXfeiRtvvBE/+MEPSjkmZgqY6IfVzA1hVPTeKH5WL8aX+ou/1TAqA6pXdN6MM1LN5CEy0ToOasWLiqEaxKkYSuX8m2qKmjEPDAzg7rvvRk9PD6666qpSj4lhGKYopvWM+aWXXsLo6CheeOEFdHd3o7u7G4lEotRjYxiGGRfTesb86KOP4tFHHy31WJgKQM/2aMa5NBHbpZl4XTO2bK0xmO07kXMQ60pr7d/o+JNBrdqYa6XsZ+WPkCk7ytoKas41sq+q1aqg/8Vi82rtWhE81JZMJuWED7X+yWRSTktWtudyOQwPD+Ozzz7D8PBwXjF8GvvIyAjOnDmD0dFR1WzDdDqNvr4+hEIhzaxGt9sNh8MBIP8BRSnYaqF/hMUylkpO/dUwk2wivgdmomW0xlKLTOsZM1Mb0E1tpjCP2E4rO1NKrzI0LJvNyiuKAGOhYhTJIK4AIoa7UegZMOago/A5m82GaDSKcDicJ6Zixtzg4CD6+vrkIvkDAwPw+XyYO3cumpubMTQ0hBMnTiASiQAABgcH4XA40NbWhoaGBqTTaYyOjspjjMVi6O/vR0tLC5qbm/NuZovFArfbDZfLJUdZ2O12tLa2IhAI5EWpUFq2xTK2tBWFygGQV1JJpVKQJAkulyuvXc2RKWbo0fVWq5WhfH9FqkGUJkotnCML8zTGKEqC4nzVUK55p7ZvmuEqIQHJZDIYGhpS3Uc2m0UkEskTbBFJkhCPx3HkyBHV9mg0is8++0zzHNPpNE6fPo1oNKo6w81ms+jv70cgEMgTTPEcXC4XOjo6VJNT6MFBIq7WTusMqn39FrfXm+XRQ1LtHESRrwWxMkOtOP9YmBnGAL0b2UgIyt1uND4z7bVErQgz25irmIk6ksy0621T7nazY5zq9nJfA8Y8bGNmpgyl4wsoXKZetFGq2R7pK75ami61i3WBle3k6AOQV3eY2uPxOMLhsJx0ISZeSJKEkZERuQiWx+OBx+PJq/sQDAZx5MgRhEIhzJo1C3PmzMnLyDt//jwOHjyIb775Bq2trejs7Myrl5FIJHD+/HmEQiEEAgG0trYW1NMAgJGRETgcjjynHmGz2fD111/D5XKhtbW1ICPPbrfLq3J7vV45pZzOIZfLIZVKIRaLwev1qtqQAcjOPLUqb6KJQingZC5hvqVWojJYmKuIYmdeYptRDWVRkMXXCYrMEBEL3FPBefGY8XgciUQCLpcL8XgcZ8+ezavDHI/HEYvF4Ha7MTo6iqNHj+YVre/r60NfXx86Ojrg9XrxySefYGBgQBa0YDCIYDCIpqYmtLa2IhQKIRwOy+MLh8MIh8Pw+Xxoa2uTlw8iqPYy1eOghwiJZDKZxJkzZ2QnX0NDg1yngq5PJBKB1WqFx+OB3W4vqP9MDk6Px6P6gCBHITlV1R6Gomgri/QzY7AwM5NKKb7uGjn7lNEXau16tYtJfNSQJAmDg4MYGBgoOAb9f/78eRw+fFjzwdDb24ve3l7N/kNDQ4jFYgW1MujvSCSCmTNn6p6fVrgbzX4dDkeeKCu3odVLtK4BPUyM0rbVnIWEUUjddKcWrg0LM2OaiT4c9EQdGHtwiEsrqfXXSzJRCzFTMtHZlJYoE7XgeKpmasX5x8LMMEzNUEnCvH//frz66qt51Qhff/11U31ZmBmGqRkqyca8detWbN68Ge3t7ePuy8JcIyijMLTajfrrtVssFtmcoNYu2naV20iSJDvFtGYsTqdT1w5O+9cyZ5ANV+vmtFqtiMVi8Pl8mmOgFVO0EjbIiajV38iGTJEsWhi9j+J21fCVfLKppBlzR0cHrr766qL6sjBXCXpiJP5W/q33mrK/2rL29Fu54obYnyIjzp49i1wuh4aGBtTV1eVFLQwMDKC3txfJZBJ1dXV56cu5XA6Dg4P44osvEIlEYLfb80LPcrkcYrEY+vr6EIlE4HQ688LryOl29uxZRKNRNDQ0YNasWfJyTRRW5vP5EI/Hkc1m4fP58pxoVIsjHA7D5XKhpaVFFmg6TkNDA7xeL2KxGJxOZ4G9mcatFoIIjBW411p1m7a12+2qYXGEVvgci/QYlSTMzc3N6OnpwaJFi+Rj3nbbbab6sjBXEWp1FEqd3EAiJkljK1BTCJ2yH/BtGFp/f3/eNoODgxgeHobb7UYmk8E333yTZ2cbHR1FOByG3W5HOp3GqVOn8qI5KEY6l8vJK5PE43G5PZVKIZVKwWazycWLxP6hUAihUAh1dXW49NJL4fF48oQ+nU4jFArJ4XHZbDZvJptMJnH69Gk4nU4sWLBA9UGSSCRgsVjg9/vhcDjgcDhU61RQYSPxQaKEBFntWwa9VksrWZeTShLmrq4uAGO1W8YLC3MVQrPncma1WSwW1ZhmkaGhIQSDQdW2XC6H/v5+jI6Oas7gBwYGcObMGU3zRTKZxNdff63ZnkqlcP78ec3zCIfDCAQCmoKYyWRkgdXafyAQQCAQUG2nbxlaS0UBY+YZr9er2U5xy+MNn2O0meprde7cObS3t+PP/uzPit4HCzPDMDVDJcyYX331VTzyyCPo6ekpSBDiqIwqZ6L2QzOzZSqvqTajpLRqq9WqWj+YbNE2m01zRhuNRhGPx1WrqwFjGX/xeLzADCD2HxoaQl1dnWYFODoHtf42mw2pVEq2NauhtXArQVmBWjWUzazITBl9Wu1TPcOrJcRys3rbKEmn09i0aRNOnz4Nq9WKv/7rv4bdbsemTZtgsVhw0UUX4bHHHjMV0fHII48AAHbs2IGhoSGcPn0aF1xwgeY3LzVYmCsMPced0gGoZtM0I8ipVArJZFLelpxWZFseGRlBMBiUBdfpdKK+vh4ulwuSJCEWi2FwcBC5XA51dXXIZrOIx+Ny5mAoFEJvby8ikYgcjdHU1CTXmhgZGcHJkyflFa+tViv8fr/sGBsZGcGnn36Kvr4++Tw7OjrQ0tICq9Uqr7Y9OjoKALJwkkDb7XZ0dXWhra0N0WgU0WhUToWmG0usf0xmISpjarFYUF9fL/endPH6+nrZbOHxeFBXV5d3k4vX3m63y8ejNHaxpojVapVLfor9xffaTOgXk0+xM+YPPvgAmUwGb775Jvbv34/nnnsO6XQa69evx9KlS9HT04O9e/di5cqVpseyZ88ebN++HfPmzcOpU6dw77334oYbbjDVl4W5QjCTVaflrRdD2dTaaRtaIUQJOdsSiYQsluK+kskkgsEgrFar7BAUBcRut8Pv92NkZAQfffQREolE3iw6nU7LQj8wMIBIJFJQfyMcDmN4eBjHjh1DMBgsGMPZs2dx5swZNDU1IZPJFDykaPZ88cUXo6Ojo0DQkskk4vE4mpubZbuzeA4Wy9jqIj6fD62trXl2Xfr2kEgk0NLSgjlz5mimTVssFtnRp2zPZrOwWCzwer2aq12TIKvVwVD7jPBsO59ihfnCCy+UfSoUGXTo0CFceeWVAIAVK1Zg//794xLmN954A7/61a/gcrkQi8Vw1113sTBPJ8zcnBaLRV7hQ4uRkRFNswQJu5Yz0GKx4Ny5c4hGo6rtuVwO4XC4QJTF9lAoJBcnUpLJZOQHg1roIJkLOjs7NetYAJDD+LRErrm5WfOrsCRJaGlp0f2qTN88tKBaHGpjFMVYr51jmLUpVpi9Xi9Onz6NP/3TP8Xw8DBeeukl/O53v5O39fl8eYWxzEDFroCxJcnYlMFULEY3jRmbbTmZaseRGSphDJVMMdfnF7/4Bb73ve/hwQcfxNmzZ3HXXXflreATjUZNC+uGDRtgsVgwNDSEm2++GZdeeimOHTumWlVQCxZmhmFqhmJnzIFAQHbw1tfXI5PJYNGiRTh48CCWLl2Kffv2YdmyZabGsGbNmoLX/vzP/9xUX4KFmTFNKUqPmqkwZzSGcs4YJxobzkwtxdbKWLduHTZv3oy1a9cinU7jgQcewCWXXIItW7Zg27ZtmDt3LlatWmVqDGSXnggszDWAWjq2mtPI6XSq2pmp3ePxqNZkpmiBSCQCi8VSUEuCHHXk1FITX3LORSIRzSw4ioTo7++XE1wI+j8Sicjp1OI4yZn39ddfo6urq8DOS+2UdKIUd3E8drtdrnkhHh8AhoeH4ff7dZM+tGzYdB3E8Dm1fXCqdfEUO2P2+Xz4u7/7u4LXd+7cWbKxjQcW5gpBmW6ttY2aCNPqF1qk02k55djhcECSpDwBpnhiSZLgdDqRy+XklGjK0Ovr65NF3efzoaurS85oO3fuHD777DM5fM1ischfC3O5HIaGhnD48GEMDg4CGHtAdHV1oaGhQRaoVCoFSZIwf/58zJ49G6dPn8aZM2fkc4tGo3Jat9VqRUNDA+rr6+UQNI/HA5fLhVAohJGRETQ3N2PGjBlyqjMtX2WxWGSvu8vlkuO4u7q6MHv2bPn6xONxDA8Py9etsbERra2tcDgccnw39QfGwu98Pp/s7KH6HeLyXOT4o3OivqJQKB18WkWhGHUqqbrcRGBhrjDUBFrrxlXO6kTEWhfK/ZNwDg4OyrG7YjsF6R8/fhwDAwN5ThBgzBHyxRdfIJFIYGRkpCASg47d29uLr776CqFQKK89lUrh1KlTaGxsxKxZswrOwel04sILL0RdXR1+//vfF4T4kdiHw2EsWrQorw4GHX9gYAADAwNYsWKFagILhQheffXVaGxslAWVroHX65WF3O1257XTGOLxuLwWoNoM3eVyFaybKEKV8LQSZFiUx08lZP4Rx48fx1tvvZX3LXXr1q2m+rIwVyhmPjxG9k69WTSZBvT2kUgkCkRZhJI3tEin07ohRkbHpyJFWkiSpJk1SBi1BwKBAtEllDHHapCw6rUzk0ulPMg2bdqEH/7wh1yPmWGY6U0lmTJaWlpw6623FtW36BHmcjn09PTgtttuQ3d3N3p7e4vdFaOAnG3ij7JdL3qBkkGi0Whe6rXYTjUmtGZ7sVhMzuhTm4GQ7VWrvnAmk4HL5cLcuXMLVqUGxj4/IyMjOHbsmOqsOp1O48yZM3ImnBKr1YoLL7wQ7e3tqvu3WCyGM5W2trYCe7uI2+1GY2OjZmF8mo0rzUHKcejN4Iwqx1XK7K9aoOtt9DMZzJw5Ey+//DL++7//Gx9++CE+/PBD032LnjG/9957SKVSeOutt3Do0CE89dRTePHFF4vdHQN904Qo0vT1XmmPFgWZtqHaxU6nEw6HA5lMBrFYTI4MoJoNFDURjUbR398vl8MkW2sikZAdcLFYTDZx0H6pXgYdn4oTNTU1oaGhAdFoFL29vQiHw8hkMkgmk7IjbHR0FD6fD7NmzYLH48HXX3+N3t5eWZRdLhcAyH3mz5+Piy++WHbsud1uZLNZhEIhxONxdHZ2Yu7cuXK70onW3t6Oiy++WC5uRM44Mkv4fD40NzfL/cmxF4vF5AcWxb1SKnwqlZKvpZoA0BjE+iR6MzcW5OKoJBtzOp3GV199ha+++kp+7Xvf+56pvkUL80cffYTly5cDAC677DIcOXKk2F0xMBc/qzVLpg/a0NBQXh0LERJUrVRku92Ovr4+DA8PqzoMPR4PIpGIan1l6k8rkSjPyWq1oq6uDh0dHXIoHPCtDZzStY8dO4ZIJKJZsa6urg7XXnutalF5m82G5uZmXHjhhQXtoihfddVV8Pv9qt8UcrkcOjs74fP5VEXT6/XKoXK0P2X/XC6nadem19SiMZTjVGtnjKkEYc5kMrDb7fjZz35W9D6KFuZIJAK/3y//b7PZ5AExU4OZ5Aw94vG4rrONQtq00PtKD4zNeGnVEjWoHobWunh6dSYA7ZVACEmS5BhoLVwul2GJTiPTgxnThJ5wa7UzxlSCMG/cuBHPPvssrr/++oLwx71795raR9Eq6vf78zzyuVyORZlhmCmlEoT52WefBQC8//77Re+jaOffFVdcgX379gEADh06hAULFhQ9CIZhmFJAfgK9n2oIYSx6irty5Urs378fa9asgSRJePLJJ0s5rprDTJqtXkKBGRu0XiqwmfEZfeMxmmmQmUFrDORk1IJsy3r7UDrzlG16tTQogkLvPDOZjGHssx7k4GNTxNRQCTPmUlC0MFutVvz85z8v5VhqCj2hVFuRRIwMIMR2ZdicMiIjk8nA6/UilUoV2HDJcZhOp+XMPrE/1UIeGhpCMpksqGVBhex7e3uRTCblmsXKfdhsNvj9fsTj8bzkEQpJi8fjCAQCclQHtZPNt6mpCTNnzpQXcaXr4nQ6IUkSFi1ahJaWFiSTSXlVbFoBRZIkBAIBZDIZeZURun50Lm1tbUgkErDb7aq2aLfbnZduLV5nIH/hVDVbvBiVoQUvrFpeKkmY33777bw45tdffx133nmnqb5sFC4jZma5yjoXYiicGMqlte9MJiMLMYV2OZ1OOUyO2kWRpAcACfKJEyfy/AXJZBJOpxN2ux2RSASff/55Xlo1FRJqbm6GxWJBNBqVx0DHp+WfUqkUzp49Kwut0+mU20nAm5qa8goLkbifP38esVgMl112Gb7zne/I9Td8Ph8CgQDC4TDi8Tjq6+vluhniNaT6FG1tbZg9e7YstpSOTUtqURFzmknTtXM4HPJDihyPBD0MKK1arJshvj/iw0cpGrwiSemphASTf/3Xf8X777+PgwcP4n//938BjE1uTpw4wcJc6RgVHjJqlyRJdZko4FuBHhkZ0Xw45HI5HDx4EPF4XLU9lUrh8OHDGB4eVm2nWGmanSqP73Q6kUgkcOrUKVXzhcPhkMPP1Mbo8XiwaNEiLFy4UPVGstvtaGxs1E0iyeVyuPrqqzXFLpPJoKurS9O0kU6nUVdXp7oQK4ms1ooldEw9oWARLj2VMGNevnw5WltbEQqFcNtttwEY+xzMmjXL9D5YmBmGqSmm+oFXX1+PpUuXYunSpRgcHJSLGBmFs4qwME8RRk4iM5XFjJxMRsegZBQ9Z5keanZxtf3rYdRuNAatZBSx3aiY00RhZ1/lUAmmDOJnP/sZPvjgAzn132Kx4M033zTVl4V5klHalZUCTGIqio3S0ZbNZuVkD1EU6HcymcTIyIhs91XaN0dHR3HixAkMDw/DYhkrASoeI5vNIplMYsaMGWhtbUUwGJRXzwYgpyCPjIzAYrGgvr4+LyMunU6jt7dXNmNQlIMoXuSAS6VSeWnMdB4zZszA/Pnz5frQyuSX+vp6zJw5E06nE/F4HP39/XlmmZkzZ+Lqq69Ga2srkskkzp8/n2dH93g86OjokE0pqVSq4JrX1dXJphq1B4jyfRHfA/rbKFKEKS2VJMyHDx/Ge++9V9TxWJjLhPKGVXP0iVCbmgCQYFDxIbWiRMC3gqwsykP9qY6yKLK0z3Q6LWfdic5Em82G9vZ2tLa2or+/H+fOnZNrVlD/0dFRjIyMwOPxIBgMoq+vDxbLtyuQ0KodDocDLpdLrhNBY6AIDpvNhtmzZ2PevHl58aZUvyKXy8HlcqGjoyOv1oTX68UFF1yAVCoFm82GJUuW5NmG3W43urq6kMlkEAqF0NjYKBdfoh+XyyW/Rz6fLy+yg37Tg1Acm1patRh5wYI8uVSCjZm44IIL5Cin8cLCXGbEuhBG690ZFTGiFTy0GBgY0G3/+OOPVR2GdNxEIqE5BpvNhlgsJo9BLZKkt7cXwWCwYB9ipAmJpbgN/R0IBHDRRRdpOtNcLhdmz56temNZrVZ4PB6sWLFCsxqd0+nEzJkz5f0p92+xWFBXV6cZ0iaGGmrd3NSXBXlqqCRhPnv2LK677jpccMEF8nHZlFFhVMKNasaea6bCXbHtADTXBKQ2M+PTs4mbuc4TvXEr5cZnCqkkYabU7GJgYWYYpmYwWlGGtpkMfvnLXxa89pOf/MRUXxbmElLtTp5i07kr7RjM9KWSZswtLS0Axj7zx44dMzRlirAwlwClvVTNITTR/YvOOq39U8afGrlcDi0tLThz5ozumMQUZyUNDQ0YGhrKc+yJeL1e+cZQfggpbI0yGpX9bTYbotGoZio7HVOvHRjLSqyrq1M9Nzo/vRszk8nozrjMhP9V+wO6mqkkYV6zZk3e/3/5l39pui8Lc5EoxdioAJFYs0HtyUmOJYoMEI+RSqXkcDFyYolfx2jlESpERAJGPxQuRmFo4XC4oGRrMpmUxVIZvkYhevX19Vi8eDHOnz+P/v5+eay0ekg4HIbb7ZajQkiEJUnCzJkz0d7eDpvNhtHRUQSDwbzrMHfuXMyfP19eWZoiOWgMLS0t6OzshNPplFdhyWQy8jipSL7P51N9v+x2u7zadS6XQzqdLgiP83g8chag2mIBYiSIWgQNnS+L8tQxEWH+h3/4B7z//vtIp9O4/fbbceWVV2LTpk2wWCy46KKL8Nhjj43LDCKuXBIMBuVJkRlYmItEWQtBfF2cMakVHiKxUrux6eamoj+pVKogAoKWfSLRFQWGBESSJJw7dw7hcDgvEoNSmQOBAPr6+goECvh2PT+1cdpsNnR0dGDGjBn45JNPEAwGEYlE8o7vcDhgt9vh9XoRCATQ1taWNwutr6+XCw65XC7MmTMnL7WbwtEcDofcX2x3OByor6+HJEnyyihaIUlUu0M8PtW2IHGluiAi9JAUxyPe0OJDVAy7Y6aWYuOYDx48iE8++QRvvPEG4vE4XnnlFWzduhXr16/H0qVL0dPTg71792LlypWmx9LT0yP/7XK5sHHjRtN9WZgngFGEghZqoq2GVi0MAHLShd4xBgcHNe1aZFowWtRVa4x0A4iirDy+3+9HR0eHZntjYyO6uro02x0OBzo7OzVvNIfDgTlz5ujeiG63Wze0zel06kZ56NXvpXamsijmAfnhhx9iwYIF+PGPf4xIJIKHH34Yu3fvxpVXXgkAWLFiBfbv3z8uYd6xYweGh4fR19eHrq4uNDU1me7LwswwTM1Q7Ix5eHgYZ86cwUsvvYRvvvkG99xzT943X5/Pp7qaux7/8R//geeeew7z5s3DiRMn8JOf/AQ33nijqb4szBNAK+5XyzGlRC+ml75ia62jRzZrvUzC+vp6eWVrtf709V/LYWg0s29qakIikcDp06dV+86ZMwdtbW04c+aM6j5aWlrQ1NSUl4ko4vf74fV6kUwmVdtFm7HW+MkuX8wsSkyNZ1NFdVCsMDc0NGDu3LlwOp2YO3cuXC4Xzp07J7dHo1EEAoFxjeUXv/gF3nnnHfh8PkQiEdx1110szOVGLVVXL9JB2Ue0N4tZgblcTnaeud1uuN1upFIpWZxILMTUYFFAqH8mk0FzczMaGxuRSCTkKlfUnk6nEQgEUFdXh3Q6jVAoJFfBUtpUJUkqSPMGxlJOu7q6kEqlcOTIEXz99dewWCz4gz/4A1x77bVwu92wWq1YtGgRTpw4gb6+PmSzWcycORN/+Id/KNc6bm9vx+DgoOwQbGpqwvz58+XVquvq6hCPxxGLxeSU6VmzZsHv9+fV9xCjNsiuTNdWzNoDkOfIU3u/RLGnVVWoDwt05VKs8++73/0uXn/9dfzFX/wFzp8/j3g8jquuugoHDx7E0qVLsW/fPixbtmzcYyFntN/vl+uBm4GFuUi0oi/M9FH2pZs9Ho+rOgudTiccDgdGRkYKbMKi44kiFZT79nq9cLvdOH36tBzuJvZ3Op1obW1FJBJBLBYr+CZA9t5MJiOLnOgY83g8uOKKK7Bs2TLMnDlTLlBE2Gw2LFy4EAsXLpQFU+lsa21tRVtbG+rr6+WiSiJerxderxetra15dS4Iu90uC6+aeNLDha6l2vsAaFfEo5oiNHYW58qkWGG+7rrr8Lvf/Q633HILJElCT08Purq6sGXLFmzbtg1z587FqlWrxjWWWbNm4amnnsKSJUvw+9//HrNnzzbdl4W5BBilMiu3VXtNr0Qm7V/PUWexWApEWcRqtcqriKgdx2KxyMWLzMTqKrHZbKivr4fb7dY8vt1u15w10APC5XJp3lh2u71g2SvlPoxmtEbr+ZlJKWdRrlwmEi738MMPF7y2c+fOoseydetWvPXWWzhw4ADmzZuHBx980HTfyl8ulqkYWJCYakD8Fqn2M1kcPXoU2WwWPT09+Pjjj3HixAnTfVmYGYapGWw2m6mfyeDnP/85rr32WgDA+vXr8cQTT5juy6aMElENNSAmOlso9zma2X8lzNo55bpyqaSUbIfDIduVZ82aNa6sQRbmCSJGXOiFzolF1ml7ZbvePiRJgsPhkO3AasdwuVzyyiZq/WfMmIFgMFiwDY2lra0NIyMjcro2bUORHw0NDUin03JSidgOjH0QaYxq50IZgZRKrRyDy+WCx+ORo0OU7bSMlNbiqRQ5oWWvJ1u+3oxJTGlXQ28pLWbqqSRh7uzsxLZt23DZZZfh008/RVtbm+m+LMxFoqyVARSmMNPfytAr8etUNptFOp1WjcUlh18ikcirYyGGzNHfFC1BSzFRBIUkSYjH43JEQXt7O5LJJEKhkJza7fV65UgHCq8LBoOys7C+vl4uIA+M1eYYGhrC6OgoLJax+hsdHR1ytAPV3qDMRVoNmwSVUskpfjoQCKC9vV12HEqShFgsJq/E7fF40NTUlNcuXmfK4KOoDLrmFANusVjkmGfxgSC+hxT/Koo7XUMAHI1RJVTS0lJbt27FG2+8gQ8++ADz5s3DX/3VX5nuy8I8AbRmZRaLRY6rVYPEWKteBm0TiUQ062kAQCwWUz0+Cf/58+dVBd/lcmHGjBmIRCKqMwy3241Zs2YhHA7nFV8iHA4HZsyYga6uLlnwRGg1Ea/XW1BwicbodrvR2NiIhoaGgkgNiv/0+/3yLFztGkiSBJfLVTADFlOp1epcKOPJ1a4Bib1ecSqmMqmU98rlcmHdunVF9WVhnkKKCUsjzHz4jPYvFupRg77W67VrmRUA49kLxUfrtevtn8LjjNr1rpXRV99KuckZc1SSKWMisDAzDFMzsDBXKVrmh4n0V9tGz5lHx6SUYeU2uVxOXvFZrZ32T7WPlYklZHcOBAJIp9MFGYVaKdYidrsdzc3NyGQyGB0dLTgGZeEBUM04dDqdaGhogNVqRSwWK3BaOp1ONDY2wuFwqJ4jZSw6HA5Ns5DebBr4NnRKK11edMBWw83KGMPCXGWYKdE53owwpfAqs/PE/WnZksXC69lsFqlUShY5Em9yZikFjOyglC6dSqVkRyDwrTBRtAatck1RE2o4HA45+45MCdR/dHQULpcLLS0teeUwSTyj0SisVisaGhryalEEAgHkcjlEo1FYLGNF7akcp/IcyT5N9ZfFcZAzjtKvtd4ztXZ6r0iExdeZ2mFaC3M4HMZDDz2ESCSCdDqNTZs24fLLLy/12CadYmZOorNOa2US2rdef0DdmUfbqC3HJLY7HA7EYjHNB4jL5UIoFNJNu6YVPJTXgPp3dnaqpjSTcDY1NWm222w2NDY2wu/3F5y3eA5q7fS3ODat94kciWr9xb+r4eZkxs+0FuZXX30Vy5Ytw7p163Dq1Ck8+OCDqivCVhvlfMP0zBpAaZI3jPah58gDjOtATNSRpjfLFV8vtl25ndnXmdphWgvzunXr5K+a2Wx2XOXsGIZhysW0Eea3334br732Wt5rTz75JBYvXoxgMIiHHnoImzdvLtsAGYZhxkM1CK8RhsJ866234tZbby14/YsvvsCGDRvw8MMPy+tiVTtaNmZlpphym1KYIUrxYdJa5JWw2+2aK6KYSVfWWimEMLoOE+1vFo6ymL5MmxmzGl9++SXuv/9+PPfcc1i4cGGpxzRl6EVRqKVgi5D9Va9Og95+JWlsqad0Ol0QemaxjKVa+3w+JJNJ1UVYbTYbWlpaEI/HC2pZ0Pg6OzsRjUYxPDyct9oHMJZ23dzcjGw2i9HR0bx6ERaLBR6PB4FAABaLekF+p9MptycSibx2cuwFAgE4HA45ekR5/bRqOSuvoRYU3SHa86vhJmRKx7QW5meffRapVEouY+f3+/Hiiy+WdGClRi86YrxhcnrHIHFQc7SJqcQUHqYUQAp/o/oZVPiHxujxeOB2u5FIJJBIJOSQOGr3+Xzwer1IJBJyLQuxva6uDn6/H7FYDCMjI/B6vWhubs6LCablrEZHR+FwOFBfX5/XXl9fj0wmg3g8DqvVKqdOE36/H9lsFslkEhaLBfX19bJPgs6BYrUpHlsco9p1E38r3xe9JZ9YoKcX01qYK12E9TD7pkzkazWFtxktEqq14ggJtNaHzGIZqzWhF4Ln8XhUZ9bU7vP50NjYqLl/qqehdb3sdjsaGho0TR+0ogmdhxKaIRcb6UGvGa1IIm7L1D7TWpgZhmEqkUqqLjcRKn+EFYqek81MOzCxJ7eZGb3RSg1G7UYfYDOzEyNK5fCb6mMwTCnhGbMGWvZMpeAqU3zFWsliX6WAiWUlxbq/tK0oisrjUXo22X21CrsHAgG5HjMVnwfGbMjkiKN2sZ4GrXxNwp1Op/OK69tsNrmOBQBkMhkkk8m8c3W5XHl2Zz1xpFrSygw9rWtA1288SwSxrXl6wKaMGkcUWnLYadmMqU1vlqwUBvE3FfLJZrO6Hyw1AQe+nfmqRXNYLGOF8L1eLzKZDLxeb16dC7GdCtsrHXFUEzmTycDpdMoPBDFVWlyZRNmudx3E1yVJKhibEnJmarUbwaF0tU8tvL9syjCAxM1MDG6xX5lpdqjlCKPX9JyJZhxodXV1cDgcqoXrSVC1amWQw5HalbNbZX+1fRiNla6DVn8xwqQWbj6m9NBnw+hHjcHBQVxzzTU4efIkent7cfvtt2Pt2rV47LHHDO//UsPCbJJSCMFUi0m5z4EFk5lqaIJj9KMknU6jp6dHjqXfunUr1q9fj3/6p3+CJEnYu3fv5J7HpB6NYRimjBQ7Y3766aexZs0aecHUo0ePyhnNK1aswIEDByb1PGpGmCfqeZ8OnvvpcI7M9KYYYX7nnXfQ1NSE5cuXy6+Jvgifz4dwODyp51H1zj+1VOnxfJ0221+vHnKpBI+cjGr7N3MMvaQWAIbORaP+5Rb2yXpwsAOwdikmKmPPnj2wWCz4n//5Hxw/fhwbN27E0NCQ3B6NRhEIBMoyXi2qVpjN1jbWepPG258cT2J0htZxlPsWV2pWW4lEdGpRRqC4jd6K22J/Cn9T1tugNjMLo1K4n3IlFjWnoLJduRqLeI7Ka6Bc2kpMPRejXAjROar33plpZxiRXbt2yX93d3fj8ccfxzPPPIODBw9i6dKl2LdvH5YtWzapY6paYTbDRG9CrcgAi8WiWQtD+VsZaUECQ8KjfMJTfDPVotCbwapFclgs39bbyGazplaqVvanpZkovlitvyRJciyxcgz0Go1d6xxp/8rwPPEhSMdRttMYlK+ZaWdql2JmzGps3LgRW7ZswbZt2zB37lysWrWqVEM0RU0Lc7mgN99oZqZXB8JMiJtRiI7R/o0WKzUKXdNL4KDz09pGS9BF9EIExfHptevBYjz9MPO50/tc7NixQ/57586dJRvXeGFhZhimZijVjHmqqZmojHJgNotPr7+WMy+bzcqlPbX66ZXCNGOi0OtvZvxGH3IqzF9sf457ZkpNseFylUZNz5hL7X0XhVa0IysFWOkIFMeRyWTyHHNkBybHVy6Xk1eyttvtcn9yllmtVtVSl6LNWhRkEm/R2SZm1ylretD4xXblOZINWgu1D7/oNDXzdZNhiqFWZsxVK8x6URBmLvx4+6vNbJUOJ72ZNM2Q1aCi8VrjpIgJvRoS4gxarV0UeTVHmSjCau0A5KL8ascQRV9rDKKYa/XXameY6UTVCrOIKCqV2n+iMbrKiI9i2o3GaORoM3LUmRmf1hiM2hnGDFyPucKY6M3MYjA51MLXTIYpNzUxY2YYhgHYxlxzlLuQejV8GIxQs1GrbVML58pUL7Xw+Zv2wqy0/WoJtFFCiRGU8KG1AGspMPpAUjaf1nmQ7U1vkVczIYJmx8MwpYZnzFWOkciqhcWphcYRSlFTizIQU5VFgRY/THrha2qoOeWUYxTD5yjVOZvNypEmyuL21C6Gx4nHMCqqxBEWDDMxpq0wm0EreoDijdW2UQqsWn8KG6Nqb2r9zURRaHmgRZHVaqfwOa1zEJeHUmsXr4EWbNZgJptaicpgYS4SM1+XJpIVZ/YY5W5nYWWYyYeFmWGYmoFtzNOAYhxZ47FdaznjRHu0Wi1malfWh1ZuU4oPoN41MGO2YZjJhIW5yjESFSVqzkDxdbX9iOYMrXbRaShJUl5BelGAqfi8WJBe6dCj+sZKO3Wp0LoGZtsZptywMNcI4xXo8fYfz6xSdLgp28Xax2rtwLcCrdY+nvMb7wdXaybNMExxTHthHg9GdSZKseTRRGpR6EVBmJ09T2S2wYLMTDXit1C9bSqdCcWNnDx5Et/97neRTCZLNR6GYZiiIVOG0U+lU7QwRyIRPP3003A6naUcD8MwzLSnKGGWJAlbtmzBhg0b4PF4Sj2mikVvRZNSOtmKQc/JWI7jMEwlUiszZkMb89tvv43XXnst77XOzk7ccMMNWLhwYdkGNtmYscEq057H05/sXlq1KvRWBNFDTZD1Vo82Y+fWoxo+1Mz0ZdpEZdx666249dZb815buXIl9uzZgz179iAYDOLuu+/Grl27yjbIyUIrgsHsG6nWX/lBUcYeqxWfN3P8YmfIZsZodgwMU2lMG2FW47/+67/kv//kT/4Er7zySskGVAlM9I3TCmcTXzNaM88ME612p3esavjwMoySaS3MDMMwlUgxwpxOp7F582acPn0aqVQK99xzD+bPn49NmzbBYrHgoosuwmOPPTapxY8mLMzvv/9+KcbBMAwzJbz77rtoaGjAM888g1AohJtuugkLFy7E+vXrsXTpUvT09GDv3r1YuXLlpI2p8uvfMZrozQ7Mep85yoKpJYqJyrj++utx//33Axi7H2w2G44ePYorr7wSALBixQocOHBgUs+DhblIjEJvyh2WY/YDaDS+arC3Mcx4GO/94PP54Pf7EYlEcN9992H9+vV5WbQ+nw/hcHhSz4GFuQSIb/hUiJ3R8ad6fAwzWRQ7WTl79izuvPNO3HjjjfjBD36QZ0+ORqMIBAKTeRoszKVkqgWvFrzRDDMRihHmgYEB3H333XjooYdwyy23AAAWLVqEgwcPAgD27duHJUuWTOp5sDAzDDOteemllzA6OooXXngB3d3d6O7uxvr16/H888/jtttuQzqdxqpVqyZ1TBwuxzBMzVBMuNyjjz6KRx99tGC7nTt3lnRs44FnzAzDMBUGz5gZhqkZOPOPYRimwqgVYWZTBsMwTIXBM2aGYWoGnjEzDMMwZYFnzAzD1Ay1MmNmYWYYpqaoBuE1gk0ZDMMwFQbPmBmGqRlqxZTBM2aGYZgKg4WZYRimwmBTBsMwNUOtmDJYmBmGqRlqRZjZlMEwDFNh8IyZYZiagWfMDMMwTFngGTPDMDUDz5gZhmGYssAzZoZhaoZamTGzMDMMUzPUijCzKYNhGKbC4BkzwzA1RTEz4lwuh8cffxxffPEFnE4n/uZv/gYXXHBBGUZnDp4xMwwz7XnvvfeQSqXw1ltv4cEHH8RTTz01peMpy4w5m80CAM6dO1eO3TMMU0OQTpBuTIT+/n7DGXN/f3/Bax999BGWL18OALjssstw5MiRCY9lIpRFmIPBIADgjjvuKMfuGYapQYLBYNHmA7/fj/r6etOaU19fD7/fL/8fiUTy/rfZbMhkMrDbp8baW5ajXnLJJdi1axdaW1ths9nKcQiGYWqEbDaLYDCISy65pOh9NDQ04Ne//jUikYip7f1+PxoaGvL+j0aj8v+5XG7KRBkokzC73W4sWbKkHLtmGKYGKYWjraGhIU9sx8MVV1yB3/zmN7jhhhtw6NAhLFiwYMLjmQgWSZKkKR0BwzDMFENRGf/3f/8HSZLw5JNPYt68eVM2nqoQ5pMnT2L16tU4cOAAXC7XVA+n7ITDYTz00EOIRCJIp9PYtGkTLr/88qkeVtmotFClySKdTmPz5s04ffo0UqkU7rnnHnz/+9+f6mFNGoODg7j55pvxyiuvTKkIViIVH8cciUTw9NNPw+l0TvVQJo1XX30Vy5Ytw7p163Dq1Ck8+OCD+OUvfznVwyobYqjSoUOH8NRTT+HFF1+c6mGVnXfffRcNDQ145plnEAqFcNNNN00bYU6n0+jp6YHb7Z7qoVQkFR3HLEkStmzZgg0bNsDj8Uz1cCaNdevWYc2aNQDGHCO1/i2h0kKVJovrr78e999/P4Cxz/p0cpQ//fTTWLNmDdra2qZ6KBVJxcyY3377bbz22mt5r3V2duKGG27AwoULp2hU5UftvJ988kksXrwYwWAQDz30EDZv3jxFo5scKi1UabLw+XwAxs7/vvvuw/r166d2QJPEO++8g6amJixfvhwvv/zyVA+nIqloG/PKlSvR3t4OADh06BAWL16MXbt2TfGoJocvvvgCGzZswMMPP4xrrrlmqodTVrZu3YpLL70UN9xwAwBgxYoV2Ldv3xSPanI4e/YsfvzjH2Pt2rW45ZZbpno4k8Idd9whFxs6fvw45syZgxdffBGtra1TPbTKQaoSrrvuOimRSEz1MCaFEydOSKtWrZKOHz8+1UOZFP7zP/9T2rhxoyRJkvTJJ59IP/rRj6Z4RJNDMBiUrr/+eunAgQNTPZQp44c//KH05ZdfTvUwKo7a/q5YpTz77LNIpVJ44oknAIwFv9eyM2zlypXYv38/1qxZI4cqTQdeeukljI6O4oUXXsALL7wAANi+fTs7xJjKNmUwDMNMRyo6KoNhGGY6wsLMMAxTYbAwMwzDVBgszAzDMBUGCzPDMEyFwcLMMAxTYbAwMwzDVBgszAzDMBXG/wPusp0P1vDpYgAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hexbin(x, y, gridsize=30)\n",
    "cb = plt.colorbar(label='count in bin')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`plt.hexbin` has a number of additional options, including the ability to specify weights for each point and to change the output in each bin to any NumPy aggregate (mean of weights, standard deviation of weights, etc.)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Kernel density estimation\n",
    "\n",
    "Another common method for estimating and representing densities in multiple dimensions is *kernel density estimation* (KDE).\n",
    "This will be discussed more fully in [In-Depth: Kernel Density Estimation](05.13-Kernel-Density-Estimation.ipynb), but for now I'll simply mention that KDE can be thought of as a way to \"smear out\" the points in space and add up the result to obtain a smooth function.\n",
    "One extremely quick and simple KDE implementation exists in the `scipy.stats` package.\n",
    "Here is a quick example of using KDE (see the following figure):"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWgAAAD3CAYAAAAwos73AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAjLUlEQVR4nO3de1Bc5f0/8De7myWEJWA0MW3t2obIpBoZLqm1QzEhiNiM1nrhEhCH0tpIx1yQBhQTRJoQoqKp1ZB42yqiQEq1mVptwiSGyKi1NDRDK9igZhKjyUaMYddkF9jz+4Nv9pcNew7Lwwl79uT9mmEmex6es5+D8uHZ5xomSZIEIiLSHEOwAyAiIv+YoImINIoJmohIo5igiYg0igmaiEijmKCJiDTKJFpx69at2LVrFwYHB7F06VJkZWWpGRcR0QVPKEG///772LdvH1599VWcOnUKL7zwgtpxERFd8MJEFqrU1dUhLCwM//vf/+BwOFBWVoarr77aW3769Gl0d3dj5syZMBqNqgZMRPoyPDwMu92O+fPnY+rUqcL3OXHiBBwOR0Dfa7FYEBMTI/xek0WoBf3VV1/hyJEj2LJlCw4fPozi4mK89dZbCAsLAwB0d3cjPz9f1UCJSN8aGxuxYMECobonTpzAggULAm4QRkdHY8eOHZpP0kIJOiYmBnPmzIHZbMacOXMQHh6O/v5+XHzxxQCAmTNnAhj5gc+ePVu9aIlId7744gvk5+d784YIh8MBo9GIQ4cOYWhoSPF7TSaTt44uE3RycjJeeukl/OIXv8CxY8dw6tQpnwc981ds9uzZuOyyy1QJlIj0TY3uUI/HA4/HM+b3hAqhBJ2WloYPPvgAd9xxByRJQmVlJfuaiSjoDAYDDAbl2cNjlWuJ8DS7srIyNeMgIpqwsLAw71iY0veECuEETUSkRaGUgMfCBE1EusEWNBGRRjFBExFpVCCTFUJpQgMTNBHpBlvQREQaxQRNRKRhoZSAx8IETUS6wRY0EZFGGY3GMRPwBbGSkIhIa9iCJiLSKCZoIiKN0luCDp3OGCKiAJxJ0nJf/ng8HlRWViInJwcFBQU4ePDgqO/p7+9HZmYmXC6Xz/W+vj4kJyePuq4GJmgi0o2xkrNckm5ra4Pb7UZzczNKS0tRW1vrU753714UFRXBbrf7XHc4HNi4cSPMZvN5eR4maCLSDaPRGNDXuTo7O5GamgoASEhIQHd3t0+5wWCAzWbzOZhEkiSsXbsW9913HyIiIs7L87APmoh0I5D+ZX/f43A4YLFYvK+NRiOGhoa8x2OlpKSMqvPUU09h4cKFmDdv3gQiVsYETRRCJEkSqieXuETvJ/Jek0E0QVssFjidTu9rj8fjTc5ytm/fjtmzZ6O1tRV2ux1FRUVobGwcf9AKmKCJSDdEE3RSUhJ2796NJUuWoKurC3FxcWPeZ+fOnd5/L168GC+88ML4gg0AEzQR6YpICz4jIwMdHR3Izc2FJEmoqamBzWaD1WpFenr6eYgyMEzQRKQbons9GwwGVFdX+1yLjY0d9X27du3yW1/u+kQxQRORboTSIpRAMEETnUVp0Ezkl/98DMKJ0Eoc5xsTNBGRRjFBExFpVCB7cYTSpwkmaCLSjUASNBA6SZoJmoh0I5AN+yVJwtDQ0CRFNDFM0ESkG4G2oEMFEzTpltofY5XuJ/JeInVEk4/aSUurSTAsLGzMI608Hs8kRTNxE9rN7ssvv8TChQvR19enVjxERMJEtxvVKuEW9ODgICorKzF16lQ14yEiEsYTVf7Pxo0bkZubi1mzZqkZDxGRMIPBENBXqBCK9M9//jNmzJjh3eCaiEgLRDfs1yqhLo7W1laEhYXh3XffxYcffojy8nLU19dj5syZasdHpEjtwTnRgUW5epO5f7PIR3fRj/tqx6GmYL+/moQS9NmbUhcUFKCqqorJmYiCLpAujFBZpAJwmh0R6YjeBgknnKAbGhrUiIOIaMKYoImINCrUBgHHwgRNRLqht5WETNCkGWrPyJD7RVR7poboPZU+aovcbzLn94rEPhldC+ziICLSKL0l6NBZUkNENAbRlYQejweVlZXIyclBQUEBDh48OOp7+vv7kZmZCZfLBQAYGBjAPffcgzvvvBM5OTnYt2+f+s+j+h2JiIJEdLOktrY2uN1uNDc3o7S0FLW1tT7le/fuRVFREex2u/eazWbDtddei5dffhkbNmwYdSq4GtjFQUS6YTAYxpzF4W9sorOz07t1RUJCArq7u0fd12az4fbbb/deKywshNlsBgAMDw8jPDx8ouGPwgRNqpvM5dJKI/Jqj9Yr3U/tgTG1+0nlBhCV3ieU+mrPEO2DdjgcsFgs3tdGoxFDQ0MwmUZSZEpKyqg606dPBwDY7XasXr0aFRUVEwndL3ZxEJFuiPZBWywWOJ1O72uPx+NNzkp6e3tRWFiIkpISXHPNNao+C8AETUQ6ItoHnZSUhPb2dgBAV1cX4uLixnyvAwcOYOXKlairq8PChQtVfxaAXRxEpCOBbJbkrzwjIwMdHR3Izc2FJEmoqamBzWaD1WpFenq63/vU1dXB7XZj/fr1AEZa4fX19RN/iLMwQRORroj0nRsMhlGzMGJjY0d9365du7z/VjsZ+8METUS6EcheHKG0VwcTNAlTe5N6uVkSw8PD445B5H3Gei+lMrU3y5f7mK708V3tWRda3pRfjt5WEjJBE5FuMEETEWkUEzQRkUYxQRMRaVQg+0EzQVPQTOaBmCKDhEoDbXKDd6LLuYeGhvxeP7MbmT9K8SmR+6VXmjEwZcqUcb+P2oOE52MgMJgJkC1oIiKNEl2oolVM0ESkG2xBExFpFBM0EZFGMUFT0Kk9ECh6GKrIyj+RssHBQdk633zzjWyZ3GCg3OAhoPy8SgN+cltTKm3irtQXKrIcWWnAVC4phVKyCgRncRARaRRb0EREGsUETUSkUUzQREQaxQRNRKRRTNAYGVWvqKjAZ599BrfbjeLiYtljYWhyiSy/VvvUbLfbLVvn9OnTsmWnTp3ye93hcMjWUZrFIReH0qwQpRkASjMyIiIixn0/paXeIv8dRU7o1uOp3hf8LI7t27cjJiYGjz76KE6cOIGf//znTNBEFHRsQQO48cYbkZmZCWDkL3ooHSFDRPrFBA0gMjISwMhHzxUrVmDVqlVqxkREJERvCVp4W6fPP/8cd911F2655RbcfPPNasZERCTkTIIe6+tcHo8HlZWVyMnJQUFBAQ4ePDjqe/r7+5GZmeldoXr69GksX74ceXl5uPvuu9Hf36/68wi1oI8fP46ioiJUVlbixz/+sdoxEcQH9eTKRA9KVaonNwinNKj31VdfyZYdP37c73Wl//GV3ktu0FGpS+7Mp0N/ZsyYIVsmsrezyACd6H7Qai/11morVHS70ba2NrjdbjQ3N6Orqwu1tbWor6/3lu/duxd1dXWw2+3ea6+++iri4uKwfPlyvPHGG9i8eTPWrFmj3sNAsAW9ZcsWnDx5Eps3b0ZBQQEKCgoUR+eJiCaDaAu6s7MTqampAICEhAR0d3f7lBsMBthsNsTExPitc9111+Hdd99V/XmEWtBr1qxR/S8FEZEaRFr3DocDFovF+9poNGJoaMi7CVZKSorfOlFRUQBGPnkNDAwIRiyPC1WISDdEBwktFgucTqf3tcfjkd2h0F8dp9OJ6dOnC0SsLHTOfiEiGoNoF0dSUhLa29sBAF1dXYiLixvzvZKSkrBnzx4AQHt7O5KTk9V9GLAFTUQ6ItqCzsjIQEdHB3JzcyFJEmpqamCz2WC1WmUX4S1duhTl5eVYunQppkyZgrq6OlWe4WxM0EEmsqRXZEaG6Cb6SoO/X3/9td/rZ490n+uTTz4Zd9mXX34pW0dp83252RrR0dGydWbOnClbdvYA0bnkZg6ILqUWuZ/aJ35rdaaGEtFZHAaDAdXV1T7XYmNjR33frl27vP+OiIjAk08+KRhpYJigiUg3eKo3EZGGhWLLXw4TNBHpht6WejNBE5FuMEHTuE3m0my5MqU9mpUm2MstvwaAzz77zO/1AwcOyNZRKvO3/wGgvOez0rJtuQE/pT5IteeyKiUDpXm2cjEqxS4ySBhKySoQ7IMmItIotqCJiDSKCZqISKOYoImINCyUEvBYmKCJSDfYgia/1D41W2kWh9LyZrml2SdOnJCt88UXX8iW9fX1yZZ9/PHH467zn//8R7bs2LFjfq8rPe+sWbNkyy666CLZMjlTp05VtUzpJHClGShyZUp1QinxnC9Go3HMM1JD6QxVJmgi0g22oImINEpvCTp0ZmwTEY1BdD/o823ZsmVoa2tTXGTmDxM0EemK1pIzAJSVleFf//oXbrvtNjz66KP49NNPA6rHLg4/lAb11L6f3GDg4OCgbB2lpc9yp2YfPnxYtk5PT49smcj+zfv375et8/nnn8uWyZ3QffZZceeaNm3auMuU9nW+9NJLZcuUBh3lTgNXGlhUOglcbiArlD6eB4NWuzhiY2NRVlaG/v5+rF+/HjfddBN++MMfYsWKFUhMTJStxwRNRLqh1b049uzZg9deew19fX245ZZbUFFRgaGhIdx9993Yvn27bD0maCLSDa22oLdv346lS5fiRz/6kc/15cuXK9ZjHzQR6caZFvRYX5MtOjraJzmXlZUBGDkLUQlb0ESkG1prQTc2NqK+vh5ff/01duzYAWBkTGru3LkB1WeCJiJdEUnAHo8HVVVV6O3thdlsxrp163D55Zd7y1taWtDU1ASTyYTi4mKkpaXhyJEjKCsrgyRJiI6ORl1dHSIiInzum5+fj/z8fGzZsgX33HPPuONiglaJ2qdwu1wu2Tpyp2kDwKFDh/xe//DDD8ddB1Ce/dHb2+v3utySbUB+pgYg/4s1Y8YM2ToiZd/73vdk61x22WWyZUqzOORmmijN1FD7FG4lF8rsD9FBwra2NrjdbjQ3N6Orqwu1tbWor68HMHJKfUNDA1pbW+FyuZCXl4eUlBT88Y9/xE9/+lPk5+fjiSeewJ/+9CcUFBT43Hf37t1IS0tDTEwMmpubfcpycnLGfB4maCLSDdEujs7OTqSmpgIAEhIS0N3d7S3bv38/EhMTYTabYTabYbVa0dPTgx/84AfefWwcDgdmz5496r5n9sBROplICRM0EemGaIJ2OBw+n4KMRiOGhoZgMpngcDgQFRXlLYuMjPQm5Lq6Ovz1r3+F2+3GvffeO+q+t956KwDg3nvvxcDAAMLCwtDW1oa0tLSAnkc4QY/VZ0NENNlEE7TFYoHT6fS+9ng83jMjzy1zOp2IiopCZWUlNmzYgNTUVLz99tsoLy/HM8884/c9S0pKsGjRIuzbtw8ejwc7d+7E008/PebzCM83ObvPprS0FLW1taK3IiJSheheHElJSWhvbwcAdHV1IS4uzlsWHx+Pzs5OuFwuDAwMoK+vD3FxcZg+fbq3ZT1r1iycPHlSNq5jx47hlltuQV9fH6qrq30SvhLhFrRSn00oEF3OLbJ/s9IGKXKDgUoDgUePHpUtkxvw+/LLL2XrKJUp7RXd39/v97rS/s3njnKf7ZJLLvF7XWlKUmxsrGzZVVddNe77yZ0EDsDnY+65zGaz3+siJ3cDYqdwXygDgUpEW9AZGRno6OhAbm4uJElCTU0NbDYbrFYr0tPTUVBQgLy8PEiShJKSEoSHh2Pt2rWorq6Gx+OBJEmorKyUfc/BwUHs2LEDc+fORX9///lP0Ep9NkREwWAwGMbckN/fH0aDwYDq6mqfa2f/8c/OzkZ2drZP+dy5c/HSSy8FFNevfvUrvPHGG3jggQfQ0NCA3/zmNwHVE86mSn02RETBoLWFKmfccMMNuOGGGwAAK1euDLiecB+0Up8NEVEwaHU/6C1btmDBggX4yU9+4v0KhHCT11+fDRFRMGl1N7u//e1v2Lt3r+IYjD/CCdpfn43WiB7kqlQmN+CnNEiotCpwYGDA73W5fZ2BkZVNcuRW8SnVURokFJlgrzSYprTy78orr/R7XWnln9xAIADMmzfP73Wlg2anT58uW6b0yyXXvafUHyqSKDgQqEyrXRyXXXaZ4t7gcthpTES6osU/YoODg7j55psRFxfnja+urm7MekzQRKQbWu3iuPvuu4XqcT9oItINre4HfeWVV6KjowOvvfYaTpw4oXi02tmYoIlIN7Q6i6OiogLf/e53cfDgQVxyySV48MEHA6rHBE1EuqHVBH3ixAnccccdMJlMSEpKUpxUcDZd9EGLLNsW3b9Zrp7SKdxKyzrlyuRmdwBiy6+VTgJXmmWi9D+z3IwMpVkSc+bMkS2T22xLaTm30vx7uTjkTuAGlE/hVlqIJTdbYzL3fCbtzuIAgL6+PgAjv79jrXY8QxcJmogIEF/qfb6tWbMGDz74IPr6+rBixQpUVVUFVI8Jmoh0Q2st6MWLF3vfT5IkzJgxA8ePH0dpaSnefPPNMeszQRORbmgtQb/11luQJAkPP/wwcnNzER8fj//+97945ZVXAqrPBE1EuqG1BH1mG9pDhw4hPj4ewMiUu08++SSg+rpO0KJ7Poss9VYaJFQqkxu8O3OWmT9ut1u2TO6QUqXDS8PDw2XLrrjiCtkyuUFCpeXc3/nOd2TL5AYQrVarbB25PaQBYNq0aX6vy+3dDIjv3yxXxoHAyaXVhSpRUVHYtGkT4uPjsW/fPsV9x8/GaXZEpCtam2IHAI899himT5+Ot99+G5dccgkeeeSRgOrpugVNRBcWrbagp02bhqKionHXY4ImIt3QaoIWxQRNRLqhtUHCiWKCJiLdYIIOEtEZGSL3U1rqLVemdHK3yAwPpdVQSkuVT5065fe60nJupVkXSmVyy6KV6ijNyJg9e7bf69HR0eOOAZCfuaL0s1V7g32aXKIJ2uPxoKqqCr29vTCbzVi3bp3P1gMtLS1oamqCyWRCcXEx0tLS8M0336CqqgqHDx/G4OAg1q5d651Kp5aQSdBERIEQaSG3tbXB7XajubkZXV1dqK2tRX19PYCR04gaGhrQ2toKl8uFvLw8pKSk4Pnnn8cVV1yBRx55BD09Pejp6VE9QbNJQES6YTQaA/o6V2dnJ1JTUwEACQkJ6O7u9pbt378fiYmJMJvNiIqKgtVqRU9PD9555x1MmTIFv/zlL7F582ZvfTUxQRORbohuN+pwOGCxWLyvjUYjhoaGvGVnn7MZGRkJh8OBr776CidPnsTzzz+PxYsXY+PGjao/DxM0EemGaIK2WCw+W/96PB7vqtJzy5xOJ6KiohATE4PFixcDANLS0nxa3WrRVB+0VgYCReopDRIqvZfc0mKlgUClk6fl4rj44otl68TExMiWKS2LPrvFEeh7KQ0gyt1PaZm60tJsuTLRPZpDafT/QiU6DzopKQm7d+/GkiVL0NXV5bPPeHx8PDZt2gSXywW3242+vj7ExcUhOTkZe/bswfz58/HBBx9g7ty5qj+PphI0EdFEifwhzcjIQEdHB3JzcyFJEmpqamCz2WC1WpGeno6CggLk5eVBkiSUlJQgPDwcy5Ytw5o1a5CTkwOTyXReujiYoIlIN0Sn2RkMBlRXV/tcO/skn+zsbGRnZ/uUx8TE4KmnnppAtGNjgiYi3eBCFSIijWKCJiLSKCZoFUzWbA3RmRpqz+JQmmkgt1T5zBxMf5Q2qZebJaG0Kb/STI2IiAjZMrmZJmqfmi1ymjYgtol+KP3y0mhM0AAGBgawevVqOBwODA4O4v7770diYqLasRERjQsTNACbzYZrr70WhYWF+Pjjj1FaWorXXntN7diIiMaFCRpAYWGh92Px8PCw4sdnIqLJFEoJeCxjJuht27bhxRdf9LlWU1OD+Ph42O12rF69GhUVFectQCKiQF1wLeisrCxkZWWNut7b24v77rsPZWVluOaaa/zWlSRJ9QFBufcJ9v2UPkWInAatNDCmtNRb5H5Kg4RKZXLPLLrfslyZ6B7Ncr+IofQLSuNzwSVofw4cOICVK1di06ZNmDdvntoxEREJYYIGUFdXB7fbjfXr1wMY2eTmzObWRETBwgQNMBkTkSbp7VTv0ImUiOgCw6XeRKQb7OJQgdwPSHT5tdqUPgLJzWoQ/dgk97NQe1aI2jMrlN5LJD6l9xJdmh1Kv4ikHj39d2cLmoh0gy1oIiKN0tsgIRM0EekGW9BERBrFBD0OgfywxnMvOSIDiKInOyvt+yxHaZm13InVan8MU3ompfiU6onstyz6c1ezDumX3hJ06HTGEBGdJx6PB5WVlcjJyUFBQQEOHjzoU97S0oLbbrsN2dnZ2L17t0/ZP/7xDyxcuPC8xMUuDiLSDdEWdFtbG9xuN5qbm9HV1YXa2lrvimm73Y6Ghga0trbC5XIhLy8PKSkpMJvN+Pzzz2Gz2RRPQJoItqCJSDfCwsK8Mznkvvwl6M7OTqSmpgIAEhIS0N3d7S3bv38/EhMTYTabERUVBavVip6eHrhcLjz00EOoqqo6b8/DFjQR6YZoC9rhcMBisXhfG41GDA0NwWQyweFwICoqylsWGRkJh8OB6upqFBUV4dJLL1XvAc7BFjQR6caZBD3W17ksFgucTqf3tcfj8Q6en1vmdDoxZcoU/POf/8TTTz+NgoICfP311ygpKVH9eUKmBa32TA3RZeVyZUrvJTpLQs06osuvRZZSqz1KHkqj7hRcoi3opKQk7N69G0uWLEFXVxfi4uK8ZfHx8di0aRNcLhfcbjf6+voQHx+Pv//9797vSUlJwRNPPKHeg/yfkEnQRETnS0ZGBjo6OpCbmwtJklBTUwObzQar1Yr09HQUFBQgLy8PkiShpKRk0s5hZYImIt0QXeptMBhQXV3tcy02Ntb77+zsbGRnZ8ves6OjY5yRBoZ90EREGsUWNBHpht5WEuoiQYv8wEUHCUX2W57MPYvl7qf2ftWiQumXg0KTnv4f00WCJiIC9NeCZh80EZFGsQVNRLqhtw37QydSIqILTMi0oNUeCFQisrJOK/1ak9k60MozE52htz7okEnQRERjYYImItIovSVo9kETEWkUW9BEpBtnNuwf63tCxYRa0H19fUhOTobL5VIrHiIiYaL7QWuVcAva4XBg48aNMJvN464rOrtCTVpZ+izys+Dya6ILg1CWkiQJa9euxX333YeIiAi1YyIiEnLBtaC3bduGF1980efat7/9bSxZsgTz5s07b4EREY2X3mZxjJmgs7KykJWV5XMtIyMDra2taG1thd1uR1FRERobG89bkEREgbjgErQ/O3fu9P578eLFeOGFF1QLiIhIFBO0CiZroC0U/kNofek4UShhgj7Hrl271IiDiIjOwYUqRKQboi1oj8eDqqoq9Pb2wmw2Y926dbj88su95S0tLWhqaoLJZEJxcTHS0tJw5MgRVFRUYHh4GJIkobq6GnPmzFH1ebjUm4h0RWSKXVtbG9xuN5qbm1FaWora2lpvmd1uR0NDA5qamvD888/j8ccfh9vtxu9//3vceeedaGhowLJly/D444+r/ixsQRORboi2oDs7O5GamgoASEhIQHd3t7ds//79SExMhNlshtlshtVqRU9PD8rLyxEVFQUAGB4eRnh4uIpPMoIJmoh0QzRBOxwOWCwW72uj0YihoSGYTCY4HA5vIgaAyMhIOBwOzJgxAwDw8ccfY+PGjXj66adVeor/TxcJOpRGZYlIeywWC5xOp/e1x+OByWTyW+Z0Or0J+7333sPDDz+MRx55RPX+Z4B90ESkI6JLvZOSktDe3g4A6OrqQlxcnLcsPj4enZ2dcLlcGBgYQF9fH+Li4vDee+9h/fr1eO6553D11Vefl+fRRQuaiGgiMjIy0NHRgdzcXEiShJqaGthsNlitVqSnp6OgoAB5eXmQJAklJSUIDw9HTU0NBgcHcf/99wMAvv/976O6ulrVuJigiUg3RPugDQbDqOQaGxvr/Xd2djays7N9yrdv3z6BSAPDBE1EuqG3lYTsgyYi0ii2oIlIN9iCJiKiScEWNBHpht5a0EzQRKQroZSAx8IuDiIijWILmoh0Q29dHGxBExFpFBM0EZFGsYuDiHRDb10cTNBEpBt6S9Ds4iAi0ii2oIlIN9iCJiKiScEWNBHpBlvQREQ0KdiCJiLd0FsLmgmaiHRDbwmaXRxERBrFFjQR6UootZDHIpSgh4eHsWHDBnR3d8PtdmP58uVIS0tTOzYioknh8XhQVVWF3t5emM1mrFu3Dpdffrm3vKWlBU1NTTCZTCguLkZaWhr6+/vx29/+FqdPn8asWbOwYcMGREREqBqXUIL+y1/+gqGhITQ1NeHo0aN48803fcqHh4cBAF988cXEIyQiXTuTJ87kjYk4evTomC3oo0ePjrrW1tYGt9uN5uZmdHV1oba2FvX19QAAu92OhoYGtLa2wuVyIS8vDykpKdi8eTNuuukm3HbbbXjmmWfQ3NyMwsLCCT/D2YQS9DvvvIMrrrgCv/71ryFJEtauXetTbrfbAQD5+fkTj5CILgh2u92n1ToeFosF0dHRAeec6OhoWCwW7+vOzk6kpqYCABISEtDd3e0t279/PxITE2E2m2E2m2G1WtHT04POzk4sW7YMAHDdddfh8ccfn/wEvW3bNrz44os+1y666CKEh4dj69at+OCDD/DAAw+gsbHRWz5//nw0NjZi5syZMBqNqgZMRPoyPDwMu92O+fPnC98jJiYGO3bsgMPhCOj7LRYLYmJivK8dDodPwjYajRgaGoLJZILD4UBUVJS3LDIyEg6Hw+d6ZGQkBgYGhOOXM2aCzsrKQlZWls+1kpISLFq0CGFhYbjmmmvw6aef+pRPnToVCxYsUDVQItIv0Zbz2WJiYnyS7nhYLBY4nU7va4/HA5PJ5LfM6XQiKirKe33q1KlwOp2YPn36hOL3R2iaXXJyMvbs2QMA6Onpwbe+9S1VgyIimkxJSUlob28HAHR1dSEuLs5bFh8fj87OTrhcLgwMDKCvrw9xcXFISkry5sH29nYkJyerHleYJEnSeCu53W489NBD6OvrgyRJqKqqwlVXXaV6cEREk+HMLI6PPvoIkiShpqYG7e3tsFqtSE9PR0tLC5qbmyFJEpYtW4bMzEwcP34c5eXlcDqduOiii1BXV4dp06apGpdQgp5s33zzDUpLS3Hy5ElMmTIFGzduxKWXXhrssDAwMIDVq1fD4XBgcHAQ999/PxITE4Mdlo+dO3firbfeQl1dXVDjGGsaU7D9+9//xmOPPYaGhoZghwIAGBwcREVFBT777DO43W4UFxcjPT092GEBGOkzXrNmDT755BOEhYXh4Ycf9mlxknpCYiVhS0sLrrrqKjQ2NuJnP/sZnn322WCHBACw2Wy49tpr8fLLL2PDhg2orq4Odkg+1q1bh7q6Ong8nmCH4jONqbS0FLW1tcEOyevZZ5/FmjVr4HK5gh2K1/bt2xETE4NXXnkFzz33HH73u98FOySv3bt3AwCampqwatUqPPHEE0GOSL9CYiVhYWGhd47kkSNHzktnvIjCwkKYzWYAI62K8PDwIEfkKykpCddffz2am5uDHYriNKZgs1qt+MMf/oCysrJgh+J14403IjMzEwAgSZKmZkNdf/31WLRoEQBt/T7qkeYStL9pfTU1NYiPj8ddd92Fjz76CDabTVNx2e12rF69GhUVFZMel1JsS5Yswfvvvx+UmM6lNI0p2DIzM3H48OFgh+EjMjISwMjPbcWKFVi1alVwAzqHyWRCeXk5du7ciSeffDLY4eiXFGIOHDggpaenBzsMr56eHmnJkiXS22+/HexQ/HrvvfekVatWBTsMqaamRnrjjTe8r1NTU4MYzWiHDh2SsrKygh2GjyNHjki33nqrtG3btmCHIuvYsWPSokWLJKfTGexQdCkk+qC3bt2K119/HcBIy0IrH/cOHDiAlStXoq6uDgsXLgx2OJqmNI2JRjt+/DiKioqwevVq3HHHHcEOx8frr7+OrVu3AgAiIiIQFhYGgyEkUknICf7nywDcfvvtKC8vR2trK4aHh1FTUxPskAAAdXV1cLvdWL9+PYCRCe1n1u+Tr4yMDHR0dCA3N9c7jYnkbdmyBSdPnsTmzZuxefNmACODmVOnTg1yZMANN9yABx54APn5+RgaGkJFRYUm4tKjkJhmR0R0IeLnEiIijWKCJiLSKCZoIiKNYoImItIoJmgiIo1igiYi0igmaCIijfp/OqsZgbxC49kAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from scipy.stats import gaussian_kde\n",
    "\n",
    "# fit an array of size [Ndim, Nsamples]\n",
    "data = np.vstack([x, y])\n",
    "kde = gaussian_kde(data)\n",
    "\n",
    "# evaluate on a regular grid\n",
    "xgrid = np.linspace(-3.5, 3.5, 40)\n",
    "ygrid = np.linspace(-6, 6, 40)\n",
    "Xgrid, Ygrid = np.meshgrid(xgrid, ygrid)\n",
    "Z = kde.evaluate(np.vstack([Xgrid.ravel(), Ygrid.ravel()]))\n",
    "\n",
    "# Plot the result as an image\n",
    "plt.imshow(Z.reshape(Xgrid.shape),\n",
    "           origin='lower', aspect='auto',\n",
    "           extent=[-3.5, 3.5, -6, 6])\n",
    "cb = plt.colorbar()\n",
    "cb.set_label(\"density\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "KDE has a smoothing length that effectively slides the knob between detail and smoothness (one example of the ubiquitous bias–variance trade-off).\n",
    "The literature on choosing an appropriate smoothing length is vast; `gaussian_kde` uses a rule of thumb to attempt to find a nearly optimal smoothing length for the input data.\n",
    "\n",
    "Other KDE implementations are available within the SciPy ecosystem, each with its own strengths and weaknesses; see, for example, `sklearn.neighbors.KernelDensity` and `statsmodels.nonparametric.KDEMultivariate`.\n",
    "For visualizations based on KDE, using Matplotlib tends to be overly verbose.\n",
    "The Seaborn library, discussed in [Visualization With Seaborn](04.14-Visualization-With-Seaborn.ipynb), provides a much more compact API for creating KDE-based visualizations."
   ]
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "jupytext": {
   "encoding": "# -*- coding: utf-8 -*-",
   "formats": "ipynb,md"
  },
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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": 4
}