ellipticity-exploration.ipynb

{
  "cells": [
    {
      "cell_type": "code",
      "source": [
        "import os\n",
        "import json\n",
        "import numpy as np\n",
        "import matplotlib.pyplot as plt\n",
        "from astropy.io import fits\n",
        "import lib.galaxy_utilities as gu"
      ],
      "outputs": [],
      "execution_count": 1,
      "metadata": {
        "collapsed": false,
        "outputHidden": false,
        "inputHidden": false
      }
    },
    {
      "cell_type": "code",
      "source": [
        "catalog = fits.open(os.path.join('lib', 'NSA_GalaxyZoo.fits'))\n",
        "hdu = catalog[1]\n",
        "smooth = hdu.data['t01_smooth_or_features_a01_smooth_debiased'] > 0.5\n",
        "featured = hdu.data['t01_smooth_or_features_a02_features_or_disk_debiased'] > 0.5\n",
        "spiral = featured & (hdu.data['t04_spiral_a08_spiral_debiased'] > 0.5)\n",
        "no_spiral = featured & (hdu.data['t04_spiral_a09_no_spiral_debiased'] > 0.5)"
      ],
      "outputs": [],
      "execution_count": 2,
      "metadata": {
        "collapsed": false,
        "outputHidden": false,
        "inputHidden": false
      }
    },
    {
      "cell_type": "code",
      "source": [
        "axRatios = hdu.data['SERSIC_BA']\n",
        "plt.hist(axRatios[smooth], alpha=0.6, label='Flagged as smooth', bins='scott', density=True)\n",
        "plt.hist(axRatios[featured], alpha=0.6, label='Flagged as featured', bins='scott', density=True)\n",
        "plt.legend()\n",
        "plt.xlabel('r-band Sersic BA')\n",
        "plt.ylabel('Density');\n",
        "\n",
        "plt.figure()\n",
        "axRatios = hdu.data['SERSIC_BA']\n",
        "plt.hist(axRatios[spiral & (axRatios > 0.0)], alpha=0.6, label='Flagged as spiral', bins='scott', density=True)\n",
        "plt.hist(axRatios[no_spiral & (axRatios > 0.0)], alpha=0.6, label='Flagged as no spiral', bins='scott', density=True)\n",
        "plt.legend()\n",
        "plt.xlabel('r-band Sersic BA')\n",
        "plt.ylabel('Density');"
      ],
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ],
            "image/png": [
              "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEICAYAAABS0fM3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAePUlEQVR4nO3deXRV9bn/8ffDVCpwKUi8ZQE22CtSInMCRaygFFCgIIPTVQtWoT/rVFtrQa1D1dV7r1yrotZLhVJbWgeQFCy2iBVRqyJgGAQpUqmNZdVANZACNsHn98feiSGck3MC2eck2Z/XWlk5Z49PNiHP+Q772ebuiIhIfDXLdgAiIpJdSgQiIjGnRCAiEnNKBCIiMadEICISc0oEIiIxF1kiMLNuZvaCmW0xs7fM7LoE2ww3s1IzKwq/bo0qHhERSaxFhMeuAL7r7uvNrB2wzsyec/ctNbZ7yd3HpXvQTp06eW5ubn3GKSLS5K1bt263u+ckWhdZInD3XcCu8PU+M9sKdAFqJoI6yc3NZe3atfUQoYhIfJjZX5Kty8gYgZnlAv2B1xOsHmJmG8zsWTPLy0Q8IiLyqSi7hgAws7bAYuDb7r63xur1wBfcvczMxgCFwMkJjjEDmAFw4oknRhyxiEi8RNoiMLOWBElgobs/XXO9u+9197Lw9XKgpZl1SrDdXHfPd/f8nJyEXVwiInKUImsRmJkB84Ct7n5vkm0+D/zd3d3MBhEkpj1RxSTSVJSXl1NcXMzBgwezHYo0MK1bt6Zr1660bNky7X2i7BoaClwKbDKzonDZTcCJAO7+CDAFuNLMKoADwIWucqgiKRUXF9OuXTtyc3MJPnOJgLuzZ88eiouL6d69e9r7RTlr6GWg1t9Qd38QeDCqGESaqoMHDyoJyBHMjOOPP56SkpI67ac7i0UaKSUBSeRofi+UCEREYi7y6aMiEr1ZT2+q1+P9aFLvlNs0b96c3r0/3a6wsJCdO3cye/ZsnnnmmXqNJ12VN5x26nTE5MMGr7CwkB49etCrVy8Ahg8fzuzZs8nPz4/83PFKBMuOKHd0uK/dn5k4RJqAz372sxQVFR22bOfOndkJpgkoLCxk3LhxVYkgk9Q1JCKRWLNmDUOGDKF///6cdtppbNu2DYD9+/dz/vnn06tXLyZOnMjgwYOrysbMmzePHj16MGjQIKZPn87VV18NQElJCZMnT6agoICCggJeeeUVAPbs2cOoUaPIy8vjiiuuINmkwyuvvJL8/Hzy8vK47bbbqpbPnDmTXr160adPH2644YYj9nvxxRfp168f/fr1o3///uzbt49Vq1YxbNgwJkyYwEknncTMmTNZuHAhgwYNonfv3uzYsQMIkuJZZ51Fnz59GDFiBO+9917S5X/84x9ZunQp3/ve9+jXr1/VMZ566ikGDRpEjx49eOmll+rjnyWheLUIRKTeHDhwgH79+gHQvXt3lixZctj6nj178tJLL9GiRQtWrlzJTTfdxOLFi3n44Yfp0KEDW7ZsYfPmzVXH+Nvf/sadd97J+vXradeuHWeddRZ9+/YF4LrrruP666/n9NNP57333mP06NFs3bqVO+64g9NPP51bb72V3/72t8ybNy9hrHfffTcdO3bk0KFDjBgxgo0bN9KlSxeWLFnC22+/jZnx0UcfHbHf7Nmzeeihhxg6dChlZWW0bt0agA0bNrB161Y6duzISSedxBVXXMGaNWu4//77mTNnDvfddx/XXHMNU6dOZerUqcyfP59rr72WwsLCpMvHjx/PuHHjmDJlStX5KyoqWLNmDcuXL+eOO+5g5cqVx/4Pl4ASgYgclURdQ9WVlpYydepUtm/fjplRXl4OwMsvv8x11wXdtKeeeip9+vQBghbEsGHD6NixIwDnnXcef/rTnwBYuXIlW7Z8Wq9y7969lJWVsXr1ap5+OihaMHbsWDp06JAwlieffJK5c+dSUVHBrl272LJlC7169aJ169ZcfvnljBs3jnHjjiyCPHToUL7zne9w8cUXM2nSJLp27QpAQUEBnTt3BuCLX/wio0aNAqB379688MILALz66qtVsV166aXceOONtS5PZNKkSQAMHDgw0m43dQ2JSCR+8IMfcOaZZ7J582aWLVt2THdBf/LJJ7z22msUFRVRVFTE+++/T9u2bdPa991332X27Nk8//zzbNy4kbFjx3Lw4EFatGjBmjVrmDJlCs888wxnn332EfvOnDmTRx99lAMHDjB06FDefvttAD7zmc9UbdOsWbOq982aNaOiouKof86aKo/bvHnzej1uTUoEIhKJ0tJSunTpAsCCBQuqlg8dOpQnn3wSgC1btrBpUzDjqaCggBdffJEPP/yQiooKFi9eXLXPqFGjmDNnTtX7ypbIGWecwa9+9SsAnn32WT788MMj4ti7dy9t2rShffv2/P3vf+fZZ58FoKysjNLSUsaMGcOPf/xjNmzYcMS+O3bsoHfv3nz/+9+noKCgKhGk47TTTuPxxx8HYOHChXzlK1+pdXm7du3Yt29f2sevT+oaEmkC0pnumWk33ngjU6dO5a677mLs2LFVy7/1rW8xdepUevXqRc+ePcnLy6N9+/Z06dKFm266iUGDBtGxY0d69uxJ+/btAXjggQe46qqr6NOnDxUVFZxxxhk88sgj3HbbbVx00UXk5eVx2mmnJaxO3LdvX/r370/Pnj3p1q0bQ4cOBWDfvn1MmDCBgwcP4u7ce++RJdHuu+8+XnjhBZo1a0ZeXh7nnHMOr776alo//5w5c7jsssu45557yMnJ4Wc/+1mtyy+88EKmT5/OAw88wKJFi+p2sY+RNbbSPvn5+X7UD6bR9FFpIrZu3cqXvvSlbIdxVA4dOkR5eTmtW7dmx44dfPWrX2Xbtm20atWKsrIy2rZtS0VFBRMnTuQb3/gGEydOzHbIjU6i3w8zW+fuCW9KUItARDJq//79nHnmmZSXl+PuPPzww7Rq1QqA22+/nZUrV3Lw4EFGjRrFueeem+Vo40GJQEQyql27dkkfNzt79uwMRyOgwWIRkdhTIhARiTklAhGRmFMiEBGJOQ0WizQFqaZG11UaU6njUIb6448/ZuzYsezevZtZs2ZxwQUX1Gn/mqWl69vtt99O27ZtExbMqwslAhE5KnEoQ/3mm28C1FpTqTZHU1q6oqKCFi0y+6dZXUMiEonGXob6gw8+4JJLLuGNN96oKg29bt06hg0bxsCBAxk9ejS7du0C4Kc//SkFBQX07duXyZMns3///oSlpYcPH171s+7evZvc3FwgKMExfvx4zjrrLEaMGAHAPffcQ0FBAX369Dks5rvvvpsePXpw+umnV13TY6UWgYgclaZehvqEE07g0UcfrerqKi8v59JLL+U3v/kNOTk5PPHEE9x8883Mnz+fSZMmMX36dABuueUW5s2bxzXXXJOwtHQy69evZ+PGjXTs2JEVK1awfft21qxZg7szfvx4Vq9eTZs2bXj88ccpKiqioqKCAQMGMHDgwLr9wyWgRCAiRyUOZair27ZtG5s3b2bkyJFAUCqjshT15s2bueWWW/joo48oKytj9OjRtV+8BEaOHFn1s69YsYIVK1bQv39/ICiQt337dvbt28fEiRM57rjjABg/fnydz5OIEoGIRKKyDPWSJUvYuXMnw4cPP+pjVZahrnwwTF1UlqF+44036NChA9OmTTusDPXzzz/PokWLePDBB/nDH/6Q9DjuTl5eXsKic9OmTaOwsJC+ffuyYMECVq1alfAYLVq04JNPPgE4oix3mzZtDjvXrFmz+OY3v3nYNvfdd1+6P3adaIxARCLRFMpQV3fKKadQUlJSlQjKy8t56623gKCSaefOnSkvL2fhwoVV+9QsLZ2bm8u6desAaq0wOnr0aObPn09ZWRkA77//Ph988AFnnHEGhYWFHDhwgH379rFs2bJaY06XWgQiTUEDrJzbFMpQV9eqVSsWLVrEtddeS2lpKRUVFXz7298mLy+PO++8k8GDB5OTk8PgwYOr/vjXLC19ww03cP755zN37tzDrklNo0aNYuvWrQwZMgSAtm3b8stf/pIBAwZwwQUX0LdvX0444QQKCgrq9o+ShMpQV9cA/zOJJKIy1FIblaEWkQZNZagbHiUCEckolaFueDRYLNJINbZuXcmMo/m9UCIQaYRat27Nnj17lAzkMO7Onj176jzNVl1DIo1Q165dKS4upqSkJNuhSAPTunVrunbtWqd9lAhEGqGWLVvSvXv3bIchTYS6hkREYk6JQEQk5iJLBGbWzcxeMLMtZvaWmR1xN5cFHjCzd8xso5kNiCoeERFJLMoxggrgu+6+3szaAevM7Dl331Jtm3OAk8OvwcBPwu8iIpIhkbUI3H2Xu68PX+8DtgJdamw2AXjMA68BnzOzzlHFJCIiR8rIGIGZ5QL9gddrrOoC/LXa+2KOTBYiIhKhyBOBmbUFFgPfdve9R3mMGWa21szWat60iEj9ijQRmFlLgiSw0N2fTrDJ+0C3au+7hssO4+5z3T3f3fNzcnKiCVZEJKainDVkwDxgq7snK/S9FPh6OHvoy0Cpu++KKiYRETlSlLOGhgKXApvMrPLBpjcBJwK4+yPAcmAM8A6wH7gswnhERCSByBKBu78MWIptHLgqqhhERCQ13VksIhJzSgQiIjGnRCAiEnNKBCIiMadEICISc0oEIiIxp0QgIhJzSgQiIjGnRCAiEnNKBCIiMadEICISc0oEIiIxp0QgIhJzSgQiIjGnRCAiEnNKBCIiMRflE8pEJAZmPb0p5TY/mtQ7A5HI0VKLQEQk5pQIRERiTolARCTmNEYgIrVKZwxAGjclAhHJulTJRoPN0VIiEIm5THziV6uiYdMYgYhIzKlFINKIaQ6/1Ae1CEREYk6JQEQk5pQIRERiTolARCTmNFgsIg2e7jOIlhKBSBOnOfySirqGRERiTolARCTm1DUk0oCpW0cyIbIWgZnNN7MPzGxzkvXDzazUzIrCr1ujikVERJKLskWwAHgQeKyWbV5y93ERxiAiIilElgjcfbWZ5UZ1fJGmQF0/0hBke7B4iJltMLNnzSwvy7GIiMRSNgeL1wNfcPcyMxsDFAInJ9rQzGYAMwBOPPHEzEUoIo2Cbjg7NllrEbj7XncvC18vB1qaWack285193x3z8/JyclonCIiTV3WEoGZfd7MLHw9KIxlT7biERGJq8i6hszs18BwoJOZFQO3AS0B3P0RYApwpZlVAAeAC93do4pHROJLD/CpXVqJwMyeBuYBz7r7J+ns4+4XpVj/IMH0UhERyaJ0u4YeBv4T2G5m/2Vmp0QYk4iIZFBaicDdV7r7xcAAYCew0sz+aGaXmVnLKAMUEZFopT1YbGbHA9OAK4A3gfsJEsNzkUQmIiIZke4YwRLgFOAXwNfcfVe46gkzWxtVcCJROta55xqAlKYi3VlDPw3n+lcxs8+4+8funh9BXCJNgkpISGOQbtfQXQmWvVqfgYiISHbU2iIws88DXYDPmll/wMJV/wYcF3FsIo3eucX/k3RdYdcbMxiJSHKpuoZGEwwQdwXurbZ8H3BTRDGJNAjpdOvU9odepLGoNRG4+8+Bn5vZZHdfnKGYRGIhVRJJ1WI41v3lcHEuXJeqa+gSd/8lkGtm36m53t3vTbCbSKOhT/QiqbuG2oTf20YdiGTZsuuSr/va/ZmLQ+qNWgySrlRdQ/8Xfr8jM+FIZGr7Q98AaY5+9K0VDWTXTVPuOkr3hrL/IZhCegD4HdAHuD7sNpJ06VN3Ygmuy7nF/6h6fUx/lGq55tXPIYdTayJe0r2hbJS732hmEwlqDU0CVgNKBJmS6hN9DBNJU/6EJpJJ6SaCyu3GAk+5e2n4TBmpL42s60ZEmo50E8EzZvY2QdfQlWaWAxyMLiyps5gmklq7MJZ1zFwgIo1YWonA3WeG4wSl7n7IzP4JTIg2tEaoqf4xzvLPpSmeItGqy6MqexLcT1B9n8fqOR6RevP6uxoMFklHurOGfgF8ESgCDoWLHSUCEZFGL90WQT7QSw+XFxFpetItQ70Z+HyUgYiISHak2yLoBGwxszXAx5UL3X18JFGJSKOmG9Ial3QTwe1RBiHxo4Hcxk0zuZqWdKePvmhmXwBOdveVZnYc0Dza0EREJBPSnTU0HZgBdCSYPdQFeAQYEV1o0pjpE79I45HuYPFVwFBgL4C7bwdOiCooERHJnHTHCD52939V1hcKbyrTVFIROSoaTG5Y0m0RvGhmNxE8xH4k8BSwLLqwREQkU9JNBDOBEmAT8E1gOXBLVEGJiEjmpDtr6BMzKwQK3b0k4phERCSDUj283oDbgKsJWw9mdgiY4+4/jD68BqapVhcVkVhL1SK4nmC2UIG7vwtgZicBPzGz6939x1EHKCLxo+cpZ1aqRHApMNLdd1cucPc/m9klwApAiSCmdJ+ASNORarC4ZfUkUCkcJ2gZTUgiIpJJqVoE/zrKdZjZfGAc8IG7n5pgvQH3A2OA/cA0d1+fIp5oaQxARGIoVYugr5ntTfC1D+idYt8FwNm1rD8HODn8mgH8JN2gRUSk/tTaInD3oy4s5+6rzSy3lk0mAI+FD7t5zcw+Z2ad3X3X0Z5T6of6/0XqbtbTm2pd/6NJqT47Z09dnllc37oAf632vjhcpkQgIkmpPEX9S/fO4qwysxlmttbM1paU6H42EZH6lM1E8D7Qrdr7ruGyI7j7XHfPd/f8nJycjAQnIhIX2UwES4GvW+DLQKnGB0REMi+yMQIz+zUwHOhkZsUEpSpaArj7IwSF68YA7xBMH70sqlhERCS5yBKBu1+UYr0TPPBGRKTeaDC57hrFYLGIiEQnm9NHRUQyTi2GI6lFICISc2oRxJDuHBaR6tQiEBGJObUImiB94heRulCLQEQk5pQIRERiTolARCTmlAhERGJOiUBEJOaUCEREYk6JQEQk5nQfgYhINVHVIkr1TGPI3nON1SIQEYk5JQIRkZhT11AjpBISIlKf1CIQEYk5JQIRkZhTIhARiTmNEYiI1EFt00sb62Mu1SIQEYk5tQgaIM0KEpFMUotARCTmlAhERGJOiUBEJOaUCEREYk6JQEQk5pQIRERiLlbTR9OZljm4e8cMRCIi0nDEKhGIiEQpqofaRE1dQyIiMadEICISc+oayjCVjxCRhibSFoGZnW1m28zsHTObmWD9NDMrMbOi8OuKKOMREZEjRdYiMLPmwEPASKAYeMPMlrr7lhqbPuHuV0cVh4iI1C7KFsEg4B13/7O7/wt4HJgQ4flEROQoRJkIugB/rfa+OFxW02Qz22hmi8ysW4TxiIhIAtmeNbQMyHX3PsBzwM8TbWRmM8xsrZmtLSkpyWiAIiJNXZSJ4H2g+if8ruGyKu6+x90/Dt8+CgxMdCB3n+vu+e6en5OTE0mwIiJxFWUieAM42cy6m1kr4EJgafUNzKxztbfjga0RxiMiIglENmvI3SvM7Grg90BzYL67v2VmPwTWuvtS4FozGw9UAP8ApkUVT6boPgERaWwivaHM3ZcDy2ssu7Xa61nArChjEBGR2mV7sFhERLJMiUBEJOaUCEREYk5F50REGohZT2+qdf2PJvWO5LxKBCIiGdJQH1yjriERkZhTIhARiTklAhGRmFMiEBGJOQ0W15FKSIhIU6MWgYhIzCkRiIjEnBKBiEjMKRGIiMScEoGISMwpEYiIxJymj9ag6aEiEjdqEYiIxJwSgYhIzCkRiIjEnBKBiEjMKRGIiMScEoGISMxp+qiISAOR6lGW8ItIzqsWgYhIzCkRiIjEnBKBiEjMKRGIiMScEoGISMwpEYiIxJwSgYhIzCkRiIjEnBKBiEjMKRGIiMRcpInAzM42s21m9o6ZzUyw/jNm9kS4/nUzy40yHhEROVJkicDMmgMPAecAvYCLzKxXjc0uBz509/8Afgz8d1TxiIhIYlG2CAYB77j7n939X8DjwIQa20wAfh6+XgSMMDOLMCYREakhykTQBfhrtffF4bKE27h7BVAKHB9hTCIiUkOjKENtZjOAGeHbMjPbls14gE7A7izH0JDp+tRO1yc5XZvaXPfLY7k+X0i2IspE8D7Qrdr7ruGyRNsUm1kLoD2wp+aB3H0uMDeiOOvMzNa6e36242iodH1qp+uTnK5N7aK6PlF2Db0BnGxm3c2sFXAhsLTGNkuBqeHrKcAf3N0jjElERGqIrEXg7hVmdjXwe6A5MN/d3zKzHwJr3X0pMA/4hZm9A/yDIFmIiEgGRTpG4O7LgeU1lt1a7fVB4LwoY4hIg+mmaqB0fWqn65Ocrk3tIrk+pp4YEZF4U4kJEZGYUyKoRRolMr5jZlvMbKOZPW9mSadnNUWprk+17SabmZtZbGaDpHNtzOz88PfnLTP7VaZjzKY0/m+daGYvmNmb4f+vMdmIMxvMbL6ZfWBmm5OsNzN7ILx2G81swDGf1N31leCLYIB7B3AS0ArYAPSqsc2ZwHHh6yuBJ7Idd0O6PuF27YDVwGtAfrbjbijXBjgZeBPoEL4/IdtxN7DrMxe4MnzdC9iZ7bgzeH3OAAYAm5OsHwM8CxjwZeD1Yz2nWgTJpSyR4e4vuPv+8O1rBPdKxEU6JUQA7iSoIXUwk8FlWTrXZjrwkLt/CODuH2Q4xmxK5/o48G/h6/bA3zIYX1a5+2qCWZTJTAAe88BrwOfMrPOxnFOJILl0SmRUdzlBlo6LlNcnbLJ2c/ffZjKwBiCd350eQA8ze8XMXjOzszMWXfalc31uBy4xs2KCmYfXZCa0RqGuf5tSahQlJho6M7sEyAeGZTuWhsLMmgH3AtOyHEpD1YKge2g4QUtytZn1dvePshpVw3ERsMDd/9fMhhDcb3Squ3+S7cCaIrUIkkunRAZm9lXgZmC8u3+codgaglTXpx1wKrDKzHYS9GUujcmAcTq/O8XAUncvd/d3gT8RJIY4SOf6XA48CeDurwKtCeoQSZp/m+pCiSC5lCUyzKw/8H8ESSBOfbyQ4vq4e6m7d3L3XHfPJRhDGe/ua7MTbkalU16lkKA1gJl1Iugq+nMmg8yidK7Pe8AIADP7EkEiKMlolA3XUuDr4eyhLwOl7r7rWA6orqEkPL0SGfcAbYGnwscovOfu47MWdAaleX1iKc1r83tglJltAQ4B33P3IwouNkVpXp/vAj81s+sJBo6neThlpqkzs18TfEjoFI6R3Aa0BHD3RwjGTMYA7wD7gcuO+ZwxubYiIpKEuoZERGJOiUBEJOaUCEREYk6JQEQk5pQIRERiTolAmgQzm2ZmD0Z07LIky28OK4duNLMiMxtcD+fKN7MH6rD9qrCKZ5GZbTWzGTXW9wsrv8aphIXUke4jkEbHgps2LJvlBsKyB+OAAe7+cXhTWKs67N/C3StqLg9vuKvrTXcXu/taM+sI7DCzBWExNwhKNbwcfv9dHY8rMaEWgTQKZpYbfvJ9DNjM4bfYV+oWfkLebma3Vdu30MzWhZ/eZ1RbXmZmd5vZhrDw27+Hy7ub2atmtsnM7koSUmdgd2VZEXff7e5/C/cfaGYvhuf8fWVlyDC2+8xsLXCdmZ1nZpvD868OtxluZs+Er9ua2c/CODaa2eQUl6kt8E+CG9QqE+Z5BPWeRppZ6xT7S0wpEUhjcjLwsLvnuftfEqwfBEwG+gDnVatr9A13H0hQGPBaMzs+XN4GeM3d+xI8M2F6uPx+4Cfu3htIduv+CoLE8ycze9jMhgGYWUtgDjAlPOd84O5q+7Vy93x3/1/gVmB0eP5Ed6T/gKB8QG937wP8IUksC81sI7ANuNPdD4XLTwPedfcdwCpgbJL9JeaUCKQx+UtYfz2Z59x9j7sfAJ4GTg+XX2tmGwjqHXXj0+Ju/wKeCV+vA3LD10OBX4evf5HoRO5eBgwEZhDUwHnCzKYBpxAU23vOzIqAWzj8ORVPVHv9CrDAzKYTlFqo6avAQ9XO+WGSn/viMFGcCNxgnz4p7yKCWv+E3y9Ksr/EnMYIpDH5Z+ULM5tIUIMF4Irwe816KW5mwwn+oA5x9/1mtoqggBlAebX6NYc4/P9Dytor4SfvVQQVVjcBUwkSylvuPiTVz+Du/y8cYB4LrDOzganOmSKeEjNbDwwOa9RMBiaY2c0ET7M63szaufu+YzmPND1qEUij5O5L3L1f+FU5uDrSzDqa2WeBcwk+cbcHPgyTQE+CctipvEJQERPg4kQbmNkpZla9bHQ/4C8E3TM54WAyZtbSzPKSHOOL7v66u99K0KqoOe7xHHBVte071Ba0mR0H9Cd4DOQIYKO7dwsrwH4BWAxMrO0YEk9KBNKUrCH4Y7cRWBwmiN8BLcxsK/BfBN1DqVwHXBV+yk/25Ke2wM8tePj8RoLn6t4eztaZAvx32B1VRNBXn8g94UDwZuCPBM/ure4uoEPlgDLBM7ITWRh2Q60jeJjLOoJuoCU1tluMuockAVUfFRGJObUIRERiTolARCTmlAhERGJOiUBEJOaUCEREYk6JQEQk5pQIRERiTolARCTm/j+H7DHilsij5wAAAABJRU5ErkJggg==\n"
            ]
          },
          "metadata": {
            "needs_background": "light"
          }
        },
        {
          "output_type": "display_data",
          "data": {
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ],
            "image/png": [
              "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEGCAYAAAB/+QKOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3de5RU5Znv8e9Pbp0oQ0DITBYXwRkDchNMA1GIQqJIxIB4yZExE4giR0czLHMyDpoTMWqWmYRRY6KjjBJNZrwlAgcNxmgi4h3RIGIrikBikywFNFxEDA3P+aN2k6LZ3V0NvbuK6t9nrVpd9e79Vj1s23p6v++7n62IwMzMrK5Dih2AmZmVJicIMzNL5QRhZmapnCDMzCyVE4SZmaVqW+wAmlPXrl2jd+/exQ7DzOyg8eKLL26MiG5p28oqQfTu3Ztly5YVOwwzs4OGpN/Xt81DTGZmlsoJwszMUjlBmJlZqrKagzBr7Xbu3El1dTU7duwodihWYioqKujRowft2rUruI8ThFkZqa6upmPHjvTu3RtJxQ7HSkREsGnTJqqrq+nTp0/B/TzEZFZGduzYweGHH+7kYHuRxOGHH97kM0snCLMy4+Rgafbn98IJwszMUnkOwqyMXT7vlWZ9v+vOGNToPm3atGHQoL/ut2DBAtatW8fs2bN56KGHmjWeQtVeRNu1a9eifP6pp57K3XffzSc+8YmC+4wePZrZs2dTWVmZYWQNc4Ko9eCMwvb70g+zjcPsIPexj32M5cuX79W2bt264gRTIhYtWrRPW0QQERxySOkO5JRuZGZWlpYuXcpxxx3H0KFDOf7441m1ahUA27dv58tf/jL9+/dn0qRJjBgxYk/pnDvuuINPf/rTDB8+nAsuuIBLLrkEgA0bNnDmmWcybNgwhg0bxtNPPw3Apk2bGDt2LAMGDGDatGnUd+fMiy66iMrKSgYMGMCsWbP2tM+cOZP+/fszePBgvvnNb+7T74knnmDIkCEMGTKEoUOHsnXrVhYvXswJJ5zA+PHj6du3LxdeeCG7d+8GcmcwGzduZN26dfTt25evfvWrDBw4kLfffrveGEqBzyDMrFl9+OGHDBkyBIA+ffowf/78vbb369ePJ598krZt2/LYY49xxRVX8MADD3DLLbfQuXNnqqqqWLly5Z73+OMf/8g111zDSy+9RMeOHfn85z/PMcccA8CMGTO49NJLGTVqFH/4wx845ZRTeO211/jOd77DqFGjuPLKK/nlL3/JHXfckRrrd7/7Xbp06cKuXbv4whe+wIoVK+jevTvz58/n9ddfRxJ//vOf9+k3e/Zsbr75ZkaOHMm2bduoqKgAcsmvqqqKI444gnHjxjFv3jzOOuusvfq++eab3HXXXXz2s5+tN4bBgwcfwH+B5uMEYWbNKm2IKd/mzZuZMmUKb775JpLYuXMnAE899RQzZuSGegcOHLjnS3Lp0qWceOKJdOnSBYCzzz6bN954A4DHHnuMqqqqPe+9ZcsWtm3bxpIlS5g3bx4A48ePp3Pnzqmx3H///cyZM4eamhr+9Kc/UVVVRf/+/amoqOD888/ntNNO47TTTtun38iRI/nGN77BueeeyxlnnEGPHj0AGD58OEceeSQAkydP5qmnntonQRxxxBF7kkN9MZRKgvAQk5m1qG9/+9uMGTOGlStX8uCDDx7QVd+7d+/mueeeY/ny5Sxfvpz169dz2GGHFdR37dq1zJ49m9/85jesWLGC8ePHs2PHDtq2bcvSpUs566yzeOihhxg3btw+fWfOnMntt9/Ohx9+yMiRI3n99deBfZeSpi0tPfTQQxuNoVQ4QZhZi9q8eTPdu3cH4M4779zTPnLkSO6//34AqqqqeOWV3AqsYcOG8cQTT/D+++9TU1PDAw88sKfP2LFj+dGPfrTnde2ZywknnMDdd98NwMMPP8z777+/Txxbtmzh0EMPpVOnTrzzzjs8/PDDAGzbto3Nmzdz6qmncsMNN/Dyyy/v0/ett95i0KBB/Nu//RvDhg3bkyCWLl3K2rVr2b17N/fddx+jRo1q8FjUF0Op8BCTWRkrZFlqS7vsssuYMmUK1157LePHj9/T/s///M9MmTKF/v37069fPwYMGECnTp3o3r07V1xxBcOHD6dLly7069ePTp06AXDTTTdx8cUXM3jwYGpqajjhhBO49dZbmTVrFpMnT2bAgAEcf/zx9OrVa584jjnmGIYOHUq/fv3o2bMnI0eOBGDr1q1MnDiRHTt2EBFcf/31+/S98cYbefzxxznkkEMYMGAAX/ziF3n22WcZNmwYl1xyCatXr2bMmDFMmjSpwWNRXwylQvXN7h+MKisrY79vGORlrlYGXnvtNY4++uhih7Ffdu3axc6dO6moqOCtt97ipJNOYtWqVbRv355t27Zx2GGHUVNTw6RJkzjvvPMa/fJtaYsXLy7qtR6FSPv9kPRiRKRebJHZGYSkucBpwLsRMTBl+78C5+bFcTTQLSLek7QO2ArsAmrqC97Mysf27dsZM2YMO3fuJCK45ZZbaN++PQBXXXUVjz32GDt27GDs2LGcfvrpRY62dchyiOlO4MfAT9M2RsQPgB8ASPoScGlEvJe3y5iI2JhhfGZWQjp27FjvLYNnz57dwtE03ejRoxk9enSxw2hWmU1SR8QS4L1Gd8yZDNyTVSxmZtZ0RV/FJOnjwDjggbzmAH4t6UVJ0xvpP13SMknLNmzYkGWoZmatStETBPAl4Ok6w0ujIuJY4IvAxZJOqK9zRMyJiMqIqOzWrVvWsZqZtRqlkCDOoc7wUkSsT36+C8wHhhchLjOzVq2o10FI6gScCHwlr+1Q4JCI2Jo8HwtcXaQQzQ5uhS7fLlQBy7xd7rt5LVy4kKqqKmbOnFlwn+ZacpvlMtd7gNFAV0nVwCygHUBE3JrsNgn4dUR8kNf1b4H5ySXqbYG7I+JXWcVpZs3L5b6b14QJE5gwYcI+7TU1NbRtm+3f+FmuYpocEZ+KiHYR0SMi7oiIW/OSAxFxZ0ScU6ffmog4JnkMiIjvZhWjmbW8cij3fdVVV3HeeecxevRojjzySG666aY9266//noGDhzIwIEDufHGG/fpu2vXLqZOncrAgQMZNGgQN9xwA5BbJjtjxgyGDBnCwIEDWbp0KZArR1L77506dSoXXnghI0aM4LLLLqv3WDYXl9ows2bVGsp9A7z++us8/vjjbN26lb59+3LRRRexYsUKfvKTn/D8888TEYwYMYITTzyRoUOH7ulXW1Rw5cqVAHu9//bt21m+fDlLlizhvPPO27NPvurqap555hnatGnDli1bUo9lc3GCMLNm1RrKfde+b4cOHejQoQOf/OQneeedd3jqqaeYNGnSnoqtZ5xxBk8++eReCeLII49kzZo1fP3rX2f8+PGMHTt2z7bJkycDuWKDW7ZsSU1OZ599Nm3atGnwWDaXUljFZGatSDmU+wbo0KHDnudt2rShpqamoM/t3LkzL7/8MqNHj+bWW29l2rRpe7Y1tVx4cx7LNE4QZtaiyqHcd30+97nPsWDBArZv384HH3zA/Pnz+dznPrfXPhs3bmT37t2ceeaZXHvttbz00kt7tt13331A7myqU6dOe6rW1qe+Y9lcPMRkVs5KsPpwOZT7rs+xxx7L1KlTGT48d+nWtGnT9hpeAli/fj1f+9rX9tyv+rrrrtuzraKigqFDh7Jz507mzp2738eyubjcdy2X+7Yy4HLfB6/Ro0cze/ZsKiuzK15dMuW+zcyawuW+S48ThJmVhIO93PeBWrx4cbFD2Icnqc3KTDkNG1vz2Z/fCycIszJSUVHBpk2bnCRsLxHBpk2bqKioaFI/DzGZtYDL573S6D7XnTGo0X0a06NHD6qrq/G9UayuiooKevTo0aQ+ThBmJaI5kki7du3o06dPc4VkrZyHmMzMLJXPIMqJr+Uoey01VGUGThBmrZITjRXCCcLMUjmJmOcgzMwslROEmZml8hCTWZkpZGjIrBCZJQhJc4HTgHcjYmDK9tHA/wPWJk3zIuLqZNs44IdAG+D2iPheVnGaHSh/IVu5yvIM4k7gx8BPG9jnyYjY635+ktoANwMnA9XAC5IWRkRV2hu0CoUuXzUza0aZJYiIWCKp9350HQ6sjog1AJLuBSYCrTdBWNH47MBas2LPQRwn6WXgj8A3I+JVoDvwdt4+1cCI+t5A0nRgOpB61yhL4QvqzKwAxUwQLwFHRMQ2SacCC4CjmvomETEHmAO5O8o1b4itXCGJxEnErGwVbZlrRGyJiG3J80VAO0ldgfVAz7xdeyRtZmbWgop2BiHp74B3IiIkDSeXrDYBfwaOktSHXGI4B/jHYsWZKU8+WyvgK7IPXlkuc70HGA10lVQNzALaAUTErcBZwEWSaoAPgXMid5eTGkmXAI+QW+Y6N5mbMDOzFpTlKqbJjWz/MbllsGnbFgGLsojLrJZXKB04H8PyVuxVTHaw84oos7LlWkxmZpbKZxBmVnSeyC5NThDWMjwUZXbQcYKwstSUydPTq79f0H4Lely2v+GYHZQ8B2FmZqmcIMzMLJWHmKy0eK7CrGT4DMLMzFI5QZiZWSoPMdlB5/m177HAJR7MMucEYdbMvGzWyoUThJWU59e+V+wQzCzhOQgzM0vlM4gs+EZAZs2usavjXaup+TlBWIvx8NHePFdhpc4JwspaoV/C5aCQf2s5JxtXhG1+ThB2UCrGF385JBuftVhTeJLazMxSZZYgJM2V9K6klfVsP1fSCkmvSHpG0jF529Yl7cslLcsqRjMzq1+WZxB3AuMa2L4WODEiBgHXAHPqbB8TEUMiojKj+MzMrAGZzUFExBJJvRvY/kzey+eAHlnFYmZmTVcqk9TnAw/nvQ7g15ICuC0i6p5d7CFpOjAdoFevXpkGadZaeDLboAQShKQx5BLEqLzmURGxXtIngUclvR4RS9L6J8ljDkBlZWVkHrCZWStR1FVMkgYDtwMTI2JTbXtErE9+vgvMB4YXJ0Izs9araAlCUi9gHvBPEfFGXvuhkjrWPgfGAqkroczMLDuZDTFJugcYDXSVVA3MAtoBRMStwJXA4cAtkgBqkhVLfwvMT9raAndHxK+yitMa5xIZxVUOF+iVCl9t3TRZrmKa3Mj2acC0lPY1wDH79jAzs5ZU9ElqMzt4ebVTeXOpDTMzS+UEYWZmqZwgzMwslecgWjmvUDLbm1c6/ZXPIMzMLJUThJmZpXKCMDOzVE4QZmaWqqAEIWmepPGSnFDMzFqJQr/wbwH+EXhT0vck9c0wJjMzKwEFLXONiMeAxyR1AiYnz98G/gv474jYmWGMZnaQc0mOg1PBQ0aSDgemkiuw9zvgh8CxwKOZRGZmZkVV0BmEpPlAX+BnwJci4k/JpvskLcsqODMzK55Cr6T+r4hYlN8gqUNEfJTcw8HMzMpMoUNM16a0PducgZiZWWlp8AxC0t8B3YGPSRoKKNn0N8DHM46tND04o9gRNIlrLZnZ/mpsiOkUchPTPYDr89q3AldkFJOZmZWABhNERNwF3CXpzIh4oIViMjOzEtDgHISkryRPe0v6Rt1HY28uaa6kdyWtrGe7JN0kabWkFZKOzds2RdKbyWNKk/5VZmZ2wBqbpD40+XkY0DHl0Zg7gXENbP8icFTymA78J4CkLsAsYAQwHJglqXMBn2dmZs2ksSGm25Kf39mfN4+IJZJ6N7DLROCnERHAc5I+IelTwGjg0Yh4D0DSo+QSzT37E4eZmTVdoRfKfZ/cUtcPgV8Bg4FLI+K/D/DzuwNv572uTtrqa0+LbTq5sw969ep1gOEcXLxCycpNISU5XI6j5RR6HcTYiNgCnAasA/4B+NesgmqKiJgTEZURUdmtW7dih2NmVjYKTRC1ZxrjgZ9HxOZm+vz1QM+81z2StvrazcyshRSaIB6S9DrwGeA3kroBO5rh8xcCX01WM30W2JzUeXoEGCupczI5PTZpMzOzFlJoue+ZyTzE5ojYJekDchPMDZJ0D7kJ566SqsmtTGqXvOetwCLgVGA1sB34WrLtPUnXAC8kb3V17YS1mVmxXT7vlUb3ue6MQS0QSbYKLdYH0I/c9RD5fX7aUIeImNzI9gAurmfbXGBuE+IzM7NmVOgqpp8Bfw8sB3YlzUEjCcLMzA5ehZ5BVAL9k7/4zcysFSh0knol8HdZBmJmZqWl0DOIrkCVpKXAR7WNETEhk6jMzOrh+1u3nEITxFVZBmFmZqWn0GWuT0g6AjgqIh6T9HGgTbahmZlZMRU0ByHpAuAXwG1JU3dgQVZBmZlZ8RU6SX0xMBLYAhARbwKfzCooMzMrvkITxEcR8ZfaF8nFcl7yamZWxgpNEE9IugL4mKSTgZ8DD2YXlpmZFVuhq5hmAucDrwD/m1wNpduzCsrM7GBXDvWaCl3FtFvSAmBBRGzIOCYzswPm6yUOXINDTEkZ7qskbQRWAaskbZB0ZcuEZ2ZmxdLYHMSl5FYvDYuILhHRBRgBjJR0aebRmZlZ0TSWIP4JmBwRa2sbImIN8BXgq1kGZmZmxdVYgmgXERvrNibzEO2yCcnMzEpBY5PUf9nPbXaAnl/rG+iZWXE1liCOkbQlpV1ARQbxmJlZiWgwQUSEC/KZmbVShV5JvV8kjZO0StJqSTNTtt8gaXnyeEPSn/O27crbtjDLOM3MbF+FXkndZJLaADcDJwPVwAuSFkZEVe0+EXFp3v5fB4bmvcWHETEkq/jMzKxhWZ5BDAdWR8SapNDfvcDEBvafDNyTYTxmZtYEmZ1BkLtnxNt5r6vJXWS3j+RmRH2A3+Y1V0haBtQA34uI1PtPSJoOTAfo1atXM4SdPa9QMrODQaZzEE1wDvCLiNiV13ZERFQC/wjcKOnv0zpGxJyIqIyIym7durVErGZmrUKWCWI90DPvdY+kLc051Bleioj1yc81wGL2np8wM7OMZZkgXgCOktRHUntySWCf1UiS+gGdgWfz2jpL6pA870quHlRV3b5mZpadzOYgIqJG0iXAI0AbYG5EvCrpamBZRNQmi3OAeyMi/w51RwO3SdpNLol9L3/1k5lZc3FZ8PplOUlNRCwid3Oh/LYr67y+KqXfM0Bp30nDzKzMlcoktZmZlRgnCDMzS+UEYWZmqZwgzMwslROEmZmlynQVk5lZuWiNy2F9BmFmZqmcIMzMLJWHmMzMiuTyea80us91ZxTvmmGfQZiZWSonCDMzS+UEYWZmqZwgzMwslSepE4XcBnREny4tEImZWWnwGYSZmaVygjAzs1ROEGZmlspzEM2skLkMMytf5VSzyWcQZmaWKtMEIWmcpFWSVkuambJ9qqQNkpYnj2l526ZIejN5TMkyTjMz21dmQ0yS2gA3AycD1cALkhZGRFWdXe+LiEvq9O0CzAIqgQBeTPq+n1W8Zma2tyzPIIYDqyNiTUT8BbgXmFhg31OARyPivSQpPAqMyyhOMzNLkWWC6A68nfe6Ommr60xJKyT9QlLPJvZF0nRJyyQt27BhQ3PEbWZmFH+S+kGgd0QMJneWcFdT3yAi5kREZURUduvWrdkDNDNrrbJMEOuBnnmveyRte0TEpoj4KHl5O/CZQvuamVm2skwQLwBHSeojqT1wDrAwfwdJn8p7OQF4LXn+CDBWUmdJnYGxSZuZmbWQzFYxRUSNpEvIfbG3AeZGxKuSrgaWRcRC4F8kTQBqgPeAqUnf9yRdQy7JAFwdEb4CzcysBWV6JXVELAIW1Wm7Mu/55cDl9fSdC8zNMj4zM6tfsSepzcysRDlBmJlZKicIMzNL5QRhZmapnCDMzCyVE4SZmaVygjAzs1ROEGZmlsoJwszMUjlBmJlZqkxLbZSb59e6HJSZNY/Tq7/f6D4LelzWApHUz2cQZmaWygnCzMxSOUGYmVkqJwgzM0vlBGFmZqmcIMzMLJUThJmZpXKCMDOzVJkmCEnjJK2StFrSzJTt35BUJWmFpN9IOiJv2y5Jy5PHwizjNDOzfWV2JbWkNsDNwMlANfCCpIURUZW32++AyojYLuki4PvA/0q2fRgRQ7KKz8zMGpblGcRwYHVErImIvwD3AhPzd4iIxyNie/LyOaBHhvGYmVkTZJkgugNv572uTtrqcz7wcN7rCknLJD0n6fT6Okmanuy3bMOGDQcWsZmZ7VESxfokfQWoBE7Maz4iItZLOhL4raRXIuKtun0jYg4wB6CysjJaJGAzsxZy+bxXGt3nujMGZfLZWZ5BrAd65r3ukbTtRdJJwLeACRHxUW17RKxPfq4BFgNDM4zVzMzqyDJBvAAcJamPpPbAOcBeq5EkDQVuI5cc3s1r7yypQ/K8KzASyJ/cNjOzjGU2xBQRNZIuAR4B2gBzI+JVSVcDyyJiIfAD4DDg55IA/hARE4Cjgdsk7SaXxL5XZ/WTmZllLNM5iIhYBCyq03Zl3vOT6un3DJDNoJqZmRXEV1KbmVkqJwgzM0vlBGFmZqmcIMzMLJUThJmZpSqJK6nNzGxfp1d/v8A9f5bJ5/sMwszMUjlBmJlZKicIMzNL5QRhZmapnCDMzCyVE4SZmaVygjAzs1ROEGZmlsoJwszMUjlBmJlZKicIMzNL5QRhZmapnCDMzCxVpglC0jhJqyStljQzZXsHSfcl25+X1Dtv2+VJ+ypJp2QZp5mZ7SuzBCGpDXAz8EWgPzBZUv86u50PvB8R/wDcAPx70rc/cA4wABgH3JK8n5mZtZAszyCGA6sjYk1E/AW4F5hYZ5+JwF3J818AX5CkpP3eiPgoItYCq5P3MzOzFpLlDYO6A2/nva4GRtS3T0TUSNoMHJ60P1enb/e0D5E0HZievNwmadWBh35AugIbixxDqfKxaZiPT8N8fOoz478P5NgcUd+Gg/6OchExB5hT7DhqSVoWEZXFjqMU+dg0zMenYT4+9cvq2GQ5xLQe6Jn3ukfSlrqPpLZAJ2BTgX3NzCxDWSaIF4CjJPWR1J7cpPPCOvssBKYkz88CfhsRkbSfk6xy6gMcBSzNMFYzM6sjsyGmZE7hEuARoA0wNyJelXQ1sCwiFgJ3AD+TtBp4j1wSIdnvfqAKqAEujohdWcXazEpmuKsE+dg0zMenYT4+9cvk2Cj3B7uZmdnefCW1mZmlcoIwM7NUThD7oYASIt+QVCVphaTfSKp3nXE5auz45O13pqSQ1KqWLhZyfCR9OfkdelXS3S0dY7EU8P9WL0mPS/pd8v/XqcWIsxgkzZX0rqSV9WyXpJuSY7dC0rEH/KER4UcTHuQm3N8CjgTaAy8D/evsMwb4ePL8IuC+YsddSscn2a8jsITcBZGVxY67lI4PuVV7vwM6J68/Wey4S+jYzAEuSp73B9YVO+4WPD4nAMcCK+vZfirwMCDgs8DzB/qZPoNoukZLiETE4xGxPXn5HLnrOFqLQkqsAFxDrvbWjpYMrgQUcnwuAG6OiPcBIuLdFo6xWAo5NgH8TfK8E/DHFoyvqCJiCbnVnvWZCPw0cp4DPiHpUwfymU4QTZdWQiS1DEjifHJZvbVo9Pgkp749I+KXLRlYiSjk9+fTwKclPS3pOUnjWiy64irk2FwFfEVSNbAI+HrLhHZQaOp3U6MO+lIbpUzSV4BK4MRix1IqJB0CXA9MLXIopawtuWGm0eTOPpdIGhQRfy5qVKVhMnBnRPyHpOPIXUc1MCJ2FzuwcuQziKYrqAyIpJOAbwETIuKjFoqtFDR2fDoCA4HFktaRGytd2Iomqgv5/akGFkbEzshVM36DXMIod4Ucm/OB+wEi4lmgglwRP8ugRJETRNM1WkJE0lDgNnLJobWMH9dq8PhExOaI6BoRvSOiN7k5mgkRsaw44ba4QkrQLCB39oCkruSGnNa0ZJBFUsix+QPwBQBJR5NLEBtaNMrStRD4arKa6bPA5oj404G8oYeYmigKKyHyA+Aw4Oe521vwh4iYULSgW1CBx6fVKvD4PAKMlVQF7AL+NSI2FS/qllHgsfk/wH9JupTchPXUSJbwlDtJ95D7w6FrMgczC2gHEBG3kpuTOZXc/XO2A1874M9sJcfWzMyayENMZmaWygnCzMxSOUGYmVkqJwgzM0vlBGFmZqmcIKysSZoq6ccZvfe2etq/lVRhXSFpuaQRzfBZlZJuasL+i5OqqMslvSZpep3tQ5JKuq2ljIftB18HYWVDuYtOVMyyC0n5h9OAYyPio+RCt/ZN6N82ImrqticXEjb1YsJzI2KZpC7AW5LuTIrgQa5kxVPJz1818X2tlfAZhB3UJPVO/lL+KbCSvUsN1OqZ/EX9pqRZeX0XSHox+Wt/el77NknflfRyUizvb5P2PpKelfSKpGvrCelTwMba8ioRsTEi/pj0/4ykJ5LPfKS20mYS242SlgEzJJ0taWXy+UuSfUZLeih5fpiknyRxrJB0ZiOH6TDgA3IX3dUm0rPJ1cM6WVJFI/2tlXKCsHJwFHBLRAyIiN+nbB8OnAkMBs7Oq/t0XkR8hlxBxX+RdHjSfijwXEQcQ+6eFRck7T8E/jMiBgH1lTD4NbmE9IakWySdCCCpHfAj4KzkM+cC383r1z4iKiPiP4ArgVOSz0+7Av/b5MooDIqIwcBv64nlfyStAFYB10TErqT9eGBtRLwFLAbG19PfWjknCCsHv0/q39fn0YjYFBEfAvOAUUn7v0h6mVw9qJ78tSDeX4CHkucvAr2T5yOBe5LnP0v7oIjYBnwGmE6uRtB9kqYCfckVKXxU0nLg/7L3fULuy3v+NHCnpAvIlZyo6yTg5rzPfL+ef/e5SQLpBXxTf72z4WRy91og+Tm5nv7WynkOwsrBB7VPJE0iV6MGYFrys249mZA0mtwX7XERsV3SYnKF3wB25tX32cXe/580Wpsm+Ut9MbmKta8AU8glmlcj4rjG/g0RcWEysT0eeFHSZxr7zEbi2SDpJWBEUsPnTGCipG+Ru/vY4ZI6RsTWA/kcKz8+g7CyEhHzI2JI8qid1D1ZUhdJHwNOJ/cXeifg/SQ59CNXdrwxT5OrMApwbtoOkvpKyi/NPQT4Pblhnm7JJDaS2kkaUM97/H1EPB8RV5I7C6k7r/IocHHe/p0bClrSx4Gh5G7n+QVgRUT0TCrqHgE8AExq6D2sdXKCsNZgKbkvwRXAA9bQOZMAAACTSURBVEni+BXQVtJrwPfIDTM1ZgZwcXJWUN+dug4D7pJUlYz/9weuSlYPnQX8ezKstZzcXECaHyQT0CuBZ8jdmznftUDn2olscvdAT/M/yXDWi+RusvMiueGk+XX2ewAPM1kKV3M1M7NUPoMwM7NUThBmZpbKCcLMzFI5QZiZWSonCDMzS+UEYWZmqZwgzMws1f8HQOJ1oJKiYK0AAAAASUVORK5CYII=\n"
            ]
          },
          "metadata": {
            "needs_background": "light"
          }
        }
      ],
      "execution_count": 3,
      "metadata": {
        "collapsed": false,
        "outputHidden": false,
        "inputHidden": false
      }
    }
  ],
  "metadata": {
    "kernel_info": {
      "name": "python3"
    },
    "language_info": {
      "name": "python",
      "version": "3.6.7",
      "mimetype": "text/x-python",
      "codemirror_mode": {
        "name": "ipython",
        "version": 3
      },
      "pygments_lexer": "ipython3",
      "nbconvert_exporter": "python",
      "file_extension": ".py"
    },
    "kernelspec": {
      "name": "python3",
      "language": "python",
      "display_name": "Python 3"
    },
    "nteract": {
      "version": "0.14.3"
    }
  },
  "nbformat": 4,
  "nbformat_minor": 0
}