PlotVariantLocus.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib\n",
    "import numpy\n",
    "import re\n",
    "import seaborn\n",
    "import sys\n",
    "from matplotlib import pyplot as plt\n",
    "seaborn.set_context({\"figure.figsize\": (12, 6)})\n",
    "matplotlib.rcParams['axes.labelsize'] = 16\n",
    "matplotlib.rcParams['xtick.labelsize'] = 16"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_calls(fh):\n",
    "    for line in fh:\n",
    "        if line.startswith('#'):\n",
    "            continue\n",
    "        fields = line.strip().split()\n",
    "        seqid = fields[0]\n",
    "        position = int(fields[1])\n",
    "        matemapstr = re.search('MATEMAPPINGS=([^;\\n]+)', fields[7]).group(1)\n",
    "        matemaps = list()\n",
    "        for mm in matemapstr.split(','):\n",
    "            seqid, start, end = re.split(r'[:-]', mm)\n",
    "            if seqid != 'chr17':\n",
    "                continue\n",
    "            matemaps.append((int(start), int(end)))\n",
    "        matetargetstr = re.search('MATETARGET=([^;\\n]+)', fields[7]).group(1)\n",
    "        seqid, targetstart, targetend = re.split(r'[:-]', matetargetstr)\n",
    "        targetstart = int(targetstart)\n",
    "        targetend = int(targetend)\n",
    "        yield position, targetstart, targetend, matemaps"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [],
   "source": [
    "def plot_variant(variant, targetstart, targetend, matemaps, extend=800):\n",
    "    plt.plot(variant, 1.0, 'o', color='red', markersize=14)\n",
    "    plt.plot([targetstart, targetend], [1.0, 1.0], '|', color='green', markersize=14)\n",
    "    for mstart, mend in matemaps:\n",
    "        y = numpy.random.normal(loc=1.0, scale=0.005)\n",
    "        plt.plot([mstart, mend], [y, y], color='blue')\n",
    "    ticks = [variant - extend, variant, variant + extend]\n",
    "    plt.xticks(ticks, ticks)\n",
    "    plt.yticks([])\n",
    "    plt.xlim((tstart - 800, tend + 800))\n",
    "    plt.ylim((1.0 - 0.15, 1.0 + 0.15))\n",
    "    plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "10"
      ]
     },
     "execution_count": 42,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "with open('cc11.test1.vcf', 'r') as fh:\n",
    "    calls = list(get_calls(fh))\n",
    "len(calls)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArMAAAFuCAYAAACSr5B3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAE3lJREFUeJzt3X+w5Xdd3/HXOwRIMglgSLSFGViqpkIKQzsrg1NGJJjEFrWCWJkKSXRKtbZOqbUitUgQxh+gFEVRnLHFKeqUdICBCulGIGNtG3DTAkV+hBB+qSA/YkxSSQjuxz8+3+ueHO7dvbvcu7vvu4/HzJmTPedzvud7zt189nm+9/v9nhpjBAAAOjrjZK8AAAAcLzELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtnHsvgCy64YOzbt2+XVgUAAJIbb7zxs2OMC7cz9phidt++fTl48ODxrRUAAGxDVX1su2PtZgAAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtHXmyV4BAOjkOc9J3vWuk70W9/bYxyYvf/nJXgs4OWyZBQCgLVtmAeAY2AIKpxZbZgEAaEvMAgDQlpgFTklXX3/1yV6F43foUHLttckllyTnnJOccca8vuSSefuhQyd7DY9L658JsGfVGGPbg/fv3z8OHjy4i6sDMNULK+MF25+fThk33ZRcemly663JnXd+6f3nnps8+MHJgQPJRRed+PX7MrT9mQDtVNWNY4z92xlryyzATrnppuRxj0s+8YnNQzaZt3/843PcTTed2PUD2IPELMBOOHQoueyy5Pbbk6P9xmuMOe6yy9rucgBwqhCzADvhwIG5a8F2d90aY46/7rrdXS+APc55ZuEk821Ce8RLXpLcccexPeaOO5Kf/dnk8st3Z50ATgO2zALshBtuOLGPAyCJLbOcok6nrZW2gO4Rd911fI+7++6dXQ+A04wtswA74ayzju9x97//zq4HwGnGlllOSbZW0s7jH5+8/e3H9zgAjpstswA74Ud/dH4hwrE499zkuc/dnfUBOE2IWYCdcNll85u9qrY3vmqOv/TS3V0vgD1OzALshDPOmOeafcADjh60VXPcgQPzcQAcN7MowE656KLkne9MHvaw5LzzNh9z3nnz/ne+c44H4MsiZoFT0gue+IKTvQrH56KLkltuSa65JnnSk5Kzz55bX88+e/75mmvm/Q1Dtu3PBNjTamz3qxeT7N+/fxw8eHAXVwcAgNNdVd04xti/nbG2zAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKCtGmNsf3DVZ5J8bPdW54S6IMlnT/ZKAKcF8w1wIu2FOefhY4wLtzPwmGJ2L6mqg2OM/Sd7PYC9z3wDnEin25xjNwMAANoSswAAtHU6x+yvnewVAE4b5hvgRDqt5pzTdp9ZAAD6O523zAIA0NyuxWxVXV5Vb6uqT1XV3VX1R1X12qp61CZj/2FV/V5V3VlVt1fVwaq6ZOX+J1fVa6rqw1X1+eX6V6rqK7d47kdW1TVV9dll/Aer6l+tjTmjqp5XVR+tqruq6t1V9Z2bLOs/VdX7l/W6cxn3Q1V1n03GPruqPrC83g9W1Q+s3b+vqsYRLs84lvcYOKzBnPPDVfWmqvrk8v/71Vss69uq6req6qaqOlRV1x/hNV9cVQeW1/G5Zb46f23M9UeYc6498rsKbGYPzTf3qarnV9VHltfxoap6zhZjn1BV/2t5zk9V1cuq6uy1MR89wnzzq9t4a4/Lmbu14CTnJ7kxySuTfCbJw5L8WJIbqurRY4yPJUlVfX+SX1ouL8oM7McmOWdlWT+Q5NwkL05yS5KvTfLCJJdX1WPGGHduDKyq/UneluT6JP80yZ8v489dW78XJfmRJD++rOczklxTVd86xnjzyrizk7wiyYeTjCSXJ/mFJF+T5K//8lTVs5O8KslPJ/ndJE9O8sqqqjHGryzDPpnkGzZ5r16c5AlJ/vsm9wHbc6rPOc9OcnuSNyzL38p3LOtzQ5KzthpUVQ9ZnvMDSZ6e5EFJXprkv1XVE8YYh5ahP5jkAWsP/4YkL0vyxiOsB7C1vTLfvDLJVcu6vSPJk5L8XFWdO8Z48crzPibJdZmd8q1JHpE53zw0yXevLO+pSe6/9hxPS/Jvs5vzzRjjhF2S/O3MIPw3y5/3Jfl8kucc5XEXbnLbNy7L+r6V285I8r4krz/K8r4yyd1JXrh2+1uTvGcbr+O3k9yx8uczk3w6yW+sjfuPmSctvu8RlnVO5l+4a07kz8LF5XS4nCpzzsbY5frMZTlXH2nc8t+/n+T6Lcb9hyS3JXnQJuv4tKOsy68vc+D5J/tn5OKyVy7d5pvMAP/L9fsyw/vzq/NDktcn+dBqzyS5Yln23zvKurw1c2PefXbrvT/R+8x+brn+4nL9fUkOJTnipucxxmc2ufkPluuHrtz2TUkembnF4UguT3K/JK9Zu/01SR5dVY84yuM/l8OvIZlbOS7cZHn/OcmDM7e6buVpSc5L8htHeU7g2J0qc07G4S2lOzIuybcn+Z0xxm0rj/29JB9P8o+2elBVnZPku5K8aYxx6zafCzi6bvPN4zID+S1rt1+b+Vuhf5AkVXXfJN+S5LVjjHtWxr02yRdy5PnmYZlbe39zjPGX21in47LrMbvsj3G/qvrazF/Dfypzy2YyI+8DSZ6x7CPyxaq6uar+xTYW/cTl+v0rt21E41lVdUNV3VNVn66qX1zbr+PizK0SN68t8w+X63vt81LTmVX1oJr71V6Ze/9luni5fu92lrfmysytuvZdgx1wis45O2pZ9iPypXNOMuedI805T40P0LAjms83G3H5hbXb716u/85y/dWZcXuv+WaMcVfmLphHmm+elaSyy/PNidgy+47MN+amJI9JcskY49PLfQ/J3NfjpUl+Jsllmftk/NL6zsyrquq8JC/P/CG/YeWuhyzX/yXJgSSXJnlJ5n4lv7Uy7vwkt41l+/eKW1fuX/WUJPck+bMk1yR5xRjjRWvLy3L/dpa38ToemuSSzE8sX9xsDHDMTsU5Z6d9ReY/EOtzTjLnnU3nnMUVmR+g17fGAMeu83zzweX68Wu3bxzbc/7a9fHON/93jPH/jmP9tm03DwDb8KzMgw/+VuYBV9ctByd8NDOmz0ty1Rjjdcv4t1XVviTPq6pfXA/Oqjoz81PPQ5P8/bUI3Ijz14wxfmL57+trnnngZ6rqkWOM1U852/U/knx9kgdmHtj1I/O4rvHjx7GsVc9a1vnVX+ZygMP2wpyzK5aDxr45yS/4AA07ou18M8Z4X1X9bpIXVtUtOXwA2MbZDLa7y9OmqurxSS7KysHyu2XXt8yOMd4/xnjHGOO3M0Pw3Mwj/pLD+5dct/awA0m+KsnfXL2xqs7I3FT9zUm+Y4zxnrXHHWl5SfJ3l+s/S/Kgqqq1cRufLu61H9kY48/HGAfHGG8dY/y7JD+V5MeWLasby0vm1pKjLm/FFUnetcnrAI7TKTrn7LTbMg+8WJ9zkjnvbDXnPDNz3reLAeyAPTDfXJV5UNm1mS3z6iTPW+775HK9VeMkR55vrsj8rfZu/pYqyQn+0oTlQIWbM09rlRzep3Qr658KfjXzFBDPGGO8dZPx213eH2aeOuKr1+7f2O/jfUdZzsHM927jQLGN5714bdyWy6uqr8/ckds/KrBLTqE5Z0eNMf4iyUfzpXNOMuedreawK5O8e4zx7t1YLziddZxvxhh/PMb4pswtwY9O8jeSvGu5+/eX6w9n7kpxr/mmqs7K3CK9WePcP/OUp28eY3z2WNfrWJ3QmK2qr0rydZlvTDJP9ZDMswus+pYkfzTG+NTKY38+c7+Q7x1jvCGbe0vmG77Z8pIZocn8BHJPku9ZG/fMJO8dY3zkKC/liZlbRW5Z/vy/M0/Btdnybk3yPzdZxpWZRzzu+icWOF2dQnPObnhjkqdU1QM3bqiqJyR5eDY5n+NyfspHxQdo2BWd55sxxp+MMd6b5K7M3Qw+kHku24wxvpDZTf942Q1iw9MzNwxudv7Yb8vckntC5ptd22e2ql6f5P8keU/meVQvSvKvMwPu55dhb07y9iSvqqoLMuPwuzJ3kv7elWU9N8kPZ5639UPLfhgbPjPG+HCSjDE+V1U/neT5VXV75omF9yf5icxzwN68jPt0Vb0sc5+VO5b1/O7Mg7G+feV5n7Ksx5syT3dzXuapKv5ZkleNMf5kWd49VfX8zC9J+OPML024JPO0HD+0/EVYfW/ul/mJ5S0rO4oDX4ZTec5Zlrk/87yTGxsRHlVVT99Yr2Vra6rq4Zn76Cfz1H6HVsb9wVhOxp55UMkzk7xxWYcHZh4M8o4c/kd01RXLe/GbW72HwPbsofnmn2cG7Ecyt8pemXnWhCevnd7r6swvcnltVf3ysuyXJvmvY4wbN3mLrsjcLeJ3tnwTd9JuncA2yXMzvx3jtiR/kXnU3KuS7Fsb94Akv5zkTzNPD/GeJP9kbcz1mVtCN7u8em1sZf6luHlZ3seS/GTWvrggyX2S/Pvl/ruX53362pivS/K6JJ9Yxvxp5mb378nKic1Xxn9/5hGNd2eeXPgHt3hvnrqs+3fu1vvv4nK6XRrMOa8+wjL3rYy76gjjrlpb5qMz95/7/zm8v9uDN3lv7pv5LUVvOtk/JxeXvXDZQ/PNv1zW/a7M3yS/LsnFW7zmb8z8TfRdy+t5eZJzNhl3YeZvv19xon4etTwxAAC0c6K/AQwAAHaMmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0NZfActIXCxZfbaWAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x115c32320>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "pos, tstart, tend, matemaps = calls[0]\n",
    "_ = plt.plot(pos, 1.0, 'o', color='red', markersize=14)\n",
    "_ = plt.plot([tstart, tend], [1.0, 1.0], '|', color='green', markersize=14)\n",
    "for mstart, mend in matemaps:\n",
    "    y = numpy.random.normal(loc=1.0, scale=0.005)\n",
    "    _ = plt.plot([mstart, mend], [y, y], color='blue')\n",
    "ticks = [pos - 800, pos, pos + 800]\n",
    "_ = plt.xticks(ticks, ticks)\n",
    "_ = plt.yticks([])\n",
    "_ = plt.xlim((tstart - 800, tend + 800))\n",
    "_ = plt.ylim((1.0 - 0.15, 1.0 + 0.15))\n",
    "_ = plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArMAAAFuCAYAAACSr5B3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAE2hJREFUeJzt3X+w5Xdd3/HXOwRImARoSLSFGdigpkIKQzuRiVNGfkmwRa0gVKZCEp1Sra1Tai1ILZoUpmqoFEVRnGmLU9Qp6QADFWgiJLW2DbJpgSI/QgjQqiA/YkyC2RDMp398vrd7crz37t3l3t197308Zs6c7Dmf8z3fc+7ms8/zvd/v99QYIwAA0NFpJ3oFAADgWIlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANo6/WgGn3vuuePAgQN7tCoAAJDceOONXxhjnLeTsUcVswcOHMjBgwePba0AAGAHqurTOx1rNwMAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBo6/QTvQKwmRe/OHn/+0/0WtzXE56QvOY1J3otAIBVtswCANCWLbOclGwBBQB2wpZZAADaErMAALQlZgEAaEvMAgDQlgPAgJPSFddfkSuecsWJXo1jc++9yTXXJFddldxwQ3LoUHLGGcnFFycveUlyySXJaf22JbT+mewipw6Ek0u/2RTYF678L1ee6FU4NjfdlJx/fvK85yXXXZfcdVcyxry+7rp5+6MfPcc10/ZnApzSbJkF2C033ZQ88YnJ7bfPgN3MnXcmX/rSHPe7v5tccMHxXUe+araAwsnFllmA3XDvvXP3ge1CdsMYc9wll8zHAXDMxCzAbrjmmuTWW48cshvGmOOvvXZv1wvgFCdmAXbDVVcld9xxdI+5447kZ35mb9YHYJ8QswC74YYbju/jAEgiZgF2x6FDx/a4u+/e3fUA2GfELMBuOOOMY3vcAx+4u+sBsM+IWYDdcPHFx/dxACRxnlk44Xyb0CniJS9J3ve+eR7ZnTrrrOSlL927dQLYB2yZBdgNl1ySPOxhSdXOxlfN8c94xt6uF8ApzpZZOMFsAT1FnHbaPNfskb4BLJkh++AHz/Gn2aYA8NUwiwLslgsumF9R+8hHJmefvfmYs8+e9/sqW4BdIWaBk9JPPvknT/QqHJsLLkhuuSW5+urkqU9Nzjxzbn0988z556uvnvc3DNm2PxPglFZjp1+9mOSiiy4aBw8e3MPVAQBgv6uqG8cYF+1krC2zAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaKvGGDsfXPX5JJ/eu9U5rs5N8oUTvRLAvmC+AY6nU2HOedQY47ydDDyqmD2VVNXBMcZFJ3o9gFOf+QY4nvbbnGM3AwAA2hKzAAC0tZ9j9ldO9AoA+4b5Bjie9tWcs2/3mQUAoL/9vGUWAIDm9ixmq+qZVfWeqvpsVd1dVb9fVW+qqsduMvZvVtVvV9WdVXV7VR2sqqet3P/0qnpjVX2iqu5arn+pqr5mi+d+TFVdXVVfWMZ/rKr+0dqY06rqZVX1qao6VFUfqKrv3mRZ/66qPrKs153LuB+uqvttMvZFVfXR5fV+rKp+cO3+A1U1trk8/2jeY+CwBnPOj1TV26vqM8v/71dssazvqKpfr6qbqureqrp+m9d8YVVds7yOLy7z1TlrY67fZs551/bvKrCZU2i+uV9VvbyqPrm8jo9X1Yu3GPukqvrvy3N+tqpeXVVnro351DbzzS/v4K09Jqfv1YKTnJPkxiSvS/L5JI9M8mNJbqiqx40xPp0kVfUDSX5hubwiM7CfkORBK8v6wSRnJXllkluSfEOSK5M8s6oeP8a4c2NgVV2U5D1Jrk/yd5P8yTL+rLX1e0WSH03y48t6Pj/J1VX17WOMd6yMOzPJa5N8IslI8swkP5fk65P8/788VfWiJK9P8lNJfivJ05O8rqpqjPFLy7DPJPnmTd6rVyZ5UpL/vMl9wM6c7HPOi5LcnuSty/K38l3L+tyQ5IytBlXVw5fn/GiS5yZ5aJJXJflPVfWkMca9y9AfSvLgtYd/c5JXJ3nbNusBbO1UmW9el+TyZd3em+SpSf5VVZ01xnjlyvM+Psm1mZ3y7UnOz5xvHpHke1aW9+wkD1x7juck+afZy/lmjHHcLkn+cmYQ/pPlzweS3JXkxUd43Hmb3PYty7K+f+W205J8OMlbjrC8r0lyd5Ir125/d5IP7uB1/EaSO1b+fHqSzyX51bVx/zbzpMX332ZZD8r8C3f18fxZuLjsh8vJMudsjF2uT1+Wc8V245b//p0k128x7l8nuS3JQzdZx+ccYV3+zTIHnnOif0YuLqfKpdt8kxngf7Z+X2Z437U6PyR5S5KPr/ZMkkuXZf+1I6zLuzM35t1vr977473P7BeX668s19+f5N4k2256HmN8fpOb37dcP2LltqckeUzmFoftPDPJA5K8ce32NyZ5XFWdf4THfzGHX0Myt3Kct8ny/n2Sh2Vudd3Kc5KcneRXj/CcwNE7WeacjMNbSndlXJLvTPKbY4zbVh7720n+T5K/tdWDqupBSZ6X5O1jjFt3+FzAkXWbb56YGcjvXLv9XZm/FfobSVJV90/ybUneNMa4Z2Xcm5J8OdvPN4/M3Nr7a2OMP9vBOh2TPY/ZZX+MB1TVN2T+Gv6zmVs2kxl5H03y/GUfka9U1c1V9Q92sOgnL9cfWbltIxrPqKobquqeqvpcVf382n4dF2Zulbh5bZm/t1zfZ5+Xmk6vqofW3K/2stz3L9OFy/WHdrK8NZdlbtW17xrsgpN0ztlVy7LPz5+fc5I572w35zw7PkDDrmg+32zE5ZfXbr97uf4ry/XXZcbtfeabMcahzF0wt5tvXpikssfzzfHYMvvezDfmpiSPT/K0Mcbnlvsenrmvx6uS/HSSSzL3yfiF9Z2ZV1XV2Ulek/lDfuvKXQ9frv9DkmuSPCPJVZn7lfz6yrhzktw2lu3fK25duX/Vs5Lck+SPk1yd5LVjjFesLS/L/TtZ3sbreESSp2V+YvnKZmOAo3Yyzjm77S9k/gOxPuckc97ZdM5ZXJr5AXp9awxw9DrPNx9bri9eu33j2J5z1q6Pdb75X2OM/30M67dje3kA2IYXZh588OjMA66uXQ5O+FRmTJ+d5PIxxpuX8e+pqgNJXlZVP78enFV1euannkck+etrEbgR528cY/zE8t/X1zzzwE9X1WPGGKufcnbqvyb5piQPyTyw60fncV3jx49hWateuKzzG77K5QCHnQpzzp5YDhr71iQ/5wM07Iq2880Y48NV9VtJrqyqW3L4ALCNsxnsdJenTVXVxUkuyMrB8ntlz7fMjjE+MsZ47xjjNzJD8KzMI/6Sw/uXXLv2sGuSfG2Sv7R6Y1Wdlrmp+luTfNcY44Nrj9tueUnyV5frP07y0KqqtXEbny7usx/ZGONPxhgHxxjvHmP8syT/MsmPLVtWN5aXzK0lR1zeikuTvH+T1wEco5N0ztltt2UeeLE+5yRz3tlqznlB5rxvFwPYBafAfHN55kFl78psmTckedly32eW660aJ9l+vrk087fae/lbqiTH+UsTlgMVbs48rVVyeJ/Srax/KvjlzFNAPH+M8e5Nxu90eb+XeeqIr1u7f2O/jw8fYTkHM9+7jQPFNp73wrVxWy6vqr4pc0du/6jAHjmJ5pxdNcb40ySfyp+fc5I572w1h12W5ANjjA/sxXrBftZxvhlj/MEY4ymZW4Ifl+QvJnn/cvfvLNefyNyV4j7zTVWdkblFerPGeWDmKU/fMcb4wtGu19E6rjFbVV+b5Bsz35hknuohmWcXWPVtSX5/jPHZlcf+bOZ+Id83xnhrNvfOzDd8s+UlM0KT+QnkniTfuzbuBUk+NMb45BFeypMzt4rcsvz5f2Segmuz5d2a5L9tsozLMo943PNPLLBfnURzzl54W5JnVdVDNm6oqicleVQ2OZ/jcn7Kx8YHaNgTneebMcYfjjE+lORQ5m4GH808l23GGF/O7Ka/vewGseG5mRsGNzt/7Hdkbsk9LvPNnu0zW1VvSfI/k3ww8zyqFyT5x5kB97PLsHckuS7J66vq3Mw4fF7mTtLft7Kslyb5kczztn582Q9jw+fHGJ9IkjHGF6vqp5K8vKpuzzyx8EVJfiLzHLA3L+M+V1Wvztxn5Y5lPb8n82Cs71x53mct6/H2zNPdnJ15qoq/l+T1Y4w/XJZ3T1W9PPNLEv4g80sTnpZ5Wo4fXv4irL43D8j8xPLOlR3Fga/CyTznLMu8KPO8kxsbER5bVc/dWK9la2uq6lGZ++gn89R+966Me99YTsaeeVDJC5K8bVmHh2QeDPLeHP5HdNWly3vxa1u9h8DOnELzzd/PDNhPZm6VvSzzrAlPXzu91xWZX+Typqr6xWXZr0ryH8cYN27yFl2auVvEb275Ju6mvTqBbZKXZn47xm1J/jTzqLnXJzmwNu7BSX4xyR9lnh7ig0n+ztqY6zO3hG52ecPa2Mr8S3HzsrxPJ/kXWfvigiT3S/LPl/vvXp73uWtjvjHJm5P832XMH2Vudv/erJzYfGX8D2Qe0Xh35smFf2iL9+bZy7p/9169/y4u++3SYM55wzbLPLAy7vJtxl2+tszHZe4/96Uc3t/tYZu8N/fP/Jait5/on5OLy6lwOYXmm3+4rPuhzN8kvznJhVu85m/J/E30oeX1vCbJgzYZd17mb79fe7x+HrU8MQAAtHO8vwEMAAB2jZgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANDW/wN17F9UDdgqfgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1164984e0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArMAAAFuCAYAAACSr5B3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAE75JREFUeJzt3XuQZHdZx+Hvm3vihoKQiAVVsAqscomgLiFqISRCAgjGC6IUWIIICKgVr1EuIYDITVEBQ6AKvBVRiJISSy4bC1NWgBg2EpAoiREEAokGkpAAYRPc4x/nDDSd2d3ezczOvDvPU3Wqt8/8uvuc6ZnffPb0mZ4ahiEAANDRQWu9AQAAsK/ELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtHbI3g4899thh8+bNq7QpAACQXHrppZ8fhuG4RcbuVcxu3rw527dv37etAgCABVTVpxYd6zQDAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmgXXtrAvPWutNWMzOncl73pOcfHJy1FHJQQeNlyefPK7fuXOtt3CftXkOgA2phmFYePDWrVuH7du3r+LmAHyzenFleNHi89SauPLK5FGPSq6/PvnSl27/8U2bkrveNdm2LdmyZf9v3x3U4jkADihVdekwDFsXGevILMAdceWVyQknJJ/5zPIhm4zrP/3pcdyVV+7f7QM4wB2y1hsA0NbOnckppyQ33ZTs6VWuYRjHnXJK8olPjKchsKJOPz257LK13opv9uAHJ3/0R2u9FXBgE7MA+2rbtvHUgkVP1xqGcfwFFySnnrq628Z+d9VVtz84f9VVaxvYYpqNQMwC7KtXvSq5+ea9u83NNyevfKWYXQVrHW3r8cgwbARilg1nPf7AcfSkqYsv3r+3Y13zPQxrw0lbAPvqq1/dt9vt2LGy2wGwgTkyyz5Zb0c39+bIpqMnrJgjjkhuuWXvb3f44Su/LQAblCOzAPvqxBP37+0AuB1HZtknjm5Ckt/6reRDH9r1+8suZ9Om5IwzVm+bADYYR2YB9tUpp4x/2atqsfFV4/hHPWp1twtgAxGzAPvqoIPG95q90532HLRV47ht2/zBBIAVZEYFuCO2bEkuuSS55z2To49efszRR48fv+SScTwAK0bMAuvaix7+orXehD3bsmX8E7XnnZecdFJy5JHj0dcjjxyvn3fe+PGmIdviOQA2rBoW/TOMSbZu3Tps3759FTcHAICNrqouHYZh6yJjHZkFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAWzUMw+KDq65L8qnV25xVc2ySz6/1RgAbinkH2N8OpHnnXsMwHLfIwL2K2a6qavswDFvXejuAjcO8A+xvG3XecZoBAABtiVkAANraKDH7prXeAGDDMe8A+9uGnHc2xDmzAAAcmDbKkVkAAA5AKxKzVXVqVb2vqq6tqh1VdXVVvb2q7j837geraltV/W9V3VxV/1pVP7/M/R1RVa+uqmuq6paq+mBV/dAuHvseVfWWmcf+ZFW9fDfb+gNVtbOqhqo6ZO5jr6iqj1bVjdPjfryqzqyqo+bGHVxVv1pVH6uqL0/beX5Vffcyj/e4qrqoqm6YlvdX1Wl7+pwCu7fe552qunCaZ+aX0+fGPb6qzq2qK6e56cLd7PNC80lVPaCq3lFVn5vmqMur6jfm5zxg7xxA887BVfXC6T52VNV/zo+ZGftzVXVpVd1UVddV1QVV9bBlxi20z6thpSa2Y5JcmuTsJNcluWeS305ycVUdPwzDp6bQ+8ckFyd5RpKvJHlCkjdX1eHDMLxh5v7enORHkvxmkk8keW6S91bV9w/DcNnSoKranOT9ST6Z5FeS/E+SzUnus9xGVtWhSd44jfu2ZYbcKcmfJrkiyY4kP5Dk+Um+L8nsD4yXJjkjycuTvC/j+7o9P8k/VdWDhmG4enq8Ryd5Z5J3JHnZdNtnJDm/qh4/DMM/LLedwEI6zDsfTfKsuXX/PXf9x5I8eNrGI3a1s4vOJ1V19yQXJvlsktMzvufkDyd5VZLjMs5dwL45UOads5M8NWPP/EuSk5L8flVtGobhd2ce95kZu+mcaT+PSvJrSS6YtvHD07i92eeVNwzDqixJvjPJkOTXp+u/l+TWJJvmxn0wyQdnrj9out3TZtYdkjEw3zl32/ckuSTJoQtu0/OSfCzjD4IhySEL3Obl09hjZ9Z9Lsm5c+O+axr3rJl15yb5TJKDZ9YdnOTqJH+1Wp97i2WjLutp3skYlBctsM0Hzfz7oiQX7mLcQvNJkmdO+7Jl7vZ/neSatX6OLJYDbek272QM8P9Lctbc+tcnuSXJMTPrPpDkA3Pjjp727+Uz6xba59VaVvOc2S9Ml1+bLg9Lctv0iZr1xXzz6Q4/Oo1729KKYRi+lnEiPrWqDk+Sqrp3klOTvG4Yhtv2tDHT+Bckec50//u6H0v7ctPcuBuny4Pmxn15GIb/W1ox/ftLcb4yrIZ1Ne8sYhiGnQsOXXQ+OWy6XG6OMu/Ayus275wwbce759a/J+OrQ4+ZWbdc73xl2u75eWeRfV4VK/oA0zkYh1XVfTMelr42yV9NH/6z6fK1VXX3qrpzVT0j48tffzhzNw9I8slhGL4yd/eXZ/xkLR1S/8Hp8pbp/I0d0zlkf1FVd11m885Jct4wDP+8wH4cUlWbquqRGQ+nv2UYhhtnhpyd5ClVdVpV3amqvmNad3WSt8+Me1OS+1TV86vquGk5M+NLA6/f03YAe7bO553vqaovVtVtNZ6P//Q7sKuLzifnZTy14PVV9e3THPXjSX42yR/cgccHJs3nnaX/EN86t37HdPnAmXVnJ3lkVT192o97ZJxvbst4isSSRfd5dazwofbtGQ+ZD0n+M8n95j7+kIzBtzTm1iRPnxuzLcnFy9z3I6fbPGy6/tvT9ZuS/EmSkzO+vPaFaTtmX7p7SpLrk3zrdP2s7OI0g4xP4jCz/HlmXtabGfeCjF8QS+OuSHLvZcY9NskNM+NuSvLY1T7kbrFslGUdzzsvyXju2MMznnP/t9NtX7CbfdnlaQbTxxeaTzL+ELx8ZtzOJGeu9XNlsRwoS+d5J8n9p3XPnnvcM6f1b5xb//NJvjqzL9ckecgy273HfV6152OFn9z7JXlokidlPEH66iSbp4/dN8mnk7w3yeMy1vprM9b9k/fhyX3edH3+vJKfntY/Zrp+TJL/TfKLM2POyq5j9ogkW6cvhN+ZvnjeOjfm2RkPpb84ySMynuS8Pcl/Jbn7zLgTM0b0WzO+RHBqxv+5fTnJSWv9zWixHAjLepx3drOt509zx6ZdfHx358wuNJ9k/CWvj2Y8v+4npznqJdMPljPW+vmyWA6Epfu8k+SCqY1OTXLnJD8+zS9DkjfMjDttuu0fT/vxuOm21yV54My4hfZ51Z6PVXyi75zxHK1zpuvnZYy9Q+fGvTXjS2IHTdffluSKZe7vidMn+QHT9WdN1395btxdpvXPm66fneQjGaP2ztPyimnMXZN8yx7246nT2BOn68dk/B/KS+fGHZMxfP9wZt2HkvzzMvd5UZIP789vPItlIyzrZd7ZzfYt3d/37+Lju4vZheaTJK+e5qK7zI17WcagPXZ322ixWPZu6TjvJLlHxl8WWzqK+sWZxzlzGlMZ3zXhL+fu77CM77xw/sy6hfZ5tZZVOyl3GM8xvSrfOOfj+CQfGW5/8vIlGaPyW6frlyf59pp7b9eMh8Vvne5zadzuLP1Sxf2TfHfGw/E3TMvSW9N8PuMnene2T5dL+7ElyeEZf7B83TAM12d8Iu83s/r4mdvP+tDcOGAFrKN5Z4+buuC4WYvOJ8cnuWoYhhvmxl2S5NDs4q0LgX3Tcd4ZhuGzwzA8ImPUHp/x7UqX3grsounybtO2zvfOrRkPEs7PO4vs86pYtZitqrtlfLuq/5pWXZvkwVV12NzQh2Y80nn9dP3vM064PzVzX4dkPJy+bRiGpROUL57u89S5+3v0dLn0yT894/unzS5/Pn3skRnPfd2dh0+Xs/uRjL8N+HVVdUzGL+TPzqy+NuM5JPNOmBsHrIB1NO/sypMzvmT3b4vsz5xF55NrM/6i2F3mxj10ujT3wArqPO8Mw/C5YRg+Nm3X6Uk+nvGIbTIe/NuR2/fOYRnfG3t+3llkn1fHCh1iPz/JCzOeW3FSxkPVH8942H3LNOYJGf9X8N5p3CkZfyNuSPKaufv764yfxF/IeN7F30yfjO+dG/dz0+3Pme7vOdPt/ilJ7WZ7z8rcObMZj95uy3ji9A9n/EWLV2T8AnjX3O3/PuMT/JJp7BMzHjG5NcnWmXG/Pj3OudP9PXbatyHJr6zmIXeL5UBf1vO8k+RhSf4hydOn+/qJJH833e6Mufu717SdT0jyHxmPwixdv9fMuIXmk4zn1t6W8QfcE6fHf+k0P71jrZ83i6XzcgDNO89O8rSM59T/TMa36bo5yQlz41433f610+OelvGc2SHJaTPjFt7nVXleVujJPSPjCdA3Znz/sSsyvlXF5rlxj8lY/NdNn7TLpifk4LlxRyZ5TcbS/2rGv07xiF089s9m/EMIOzL+ht3rsotfrpi5zVm5fczebfoh8cmMAfuF6YfBc5McPnf7o6Yv5n/P+MsX10xfQCcs81hPnrZ/6RSHf0nypLX+hrRYui/red7J+CrNuzMeudiR8b1gP7Dc936+cV7+cstT58YuNJ9kDNp3Tdv25YyB/IIkR67182axdF4OoHnnl6ZtXzpq+o5M5+jOjTtkGnvZtB/XTft1yjJjF9rn1ViWah4AANrx12AAAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaOv/AckyY14rBfgyAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x116401278>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArMAAAFuCAYAAACSr5B3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAEtBJREFUeJzt3XuspHV9x/HPl5tAwWABjTTBbYPrhdCaukWseAECUiNqrU1FqbEaa6tt1WikaVWgjRKlpo1Sr71gq9iWRtKLFvFSarUgLLYasYpokEtFsUoBBQT36R/PszIOZ8+evZyz57vn9UqezM4zv3nmOWc2v/OemWdmahiGAABAR3vs6h0AAIDtJWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaGuvbRl8yCGHDOvWrVumXQEAgOTKK6/81jAMhy5l7DbF7Lp167Jx48bt2ysAAFiCqvraUsc6zAAAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVlg1TrzkjN39S5su02bkosuSo4/Ptl//2SPPcbT448f12/atKv3cLu0vC+ANaGGYVjy4A0bNgwbN25cxt0BuFedVRnOWPoctctdfXVy4onJt7+d3H77fS8/4IDk4IOTiy9O1q9f+f3bAe3uC6C1qrpyGIYNSxnrmVmAneHqq5Ojj06uv37hkE3G9dddN467+uqV3T+A3ZSYBdhRmzYlJ52U3HprsrVXu4ZhHHfSSW0POQBYTcQswI66+OLx0IKlHrY1DOP4j3xkefcLYA0QswA76k1vSm67bduuc9ttyRvfuDz7A7CGiFmAHXXZZSt7PQB+SMwC7Kg779y+6911187dD4A1SMwC7Kh9992+693vfjt3PwDWIDELsKOOOWZlrwfAD4lZgB316lePX4iwLQ44IDn99OXZH4A1RMwC7KiTThq/2atqaeOrxvEnnri8+wWwBohZgB21xx7jZ83e//5bD9qqcdzFF4/XA2CHmEkBdob165PLL08OPzw58MCFxxx44Hj55ZeP4wHYYWIWWLXOeOIZu3oXts369clXv5pccEFy3HHJfvuNz77ut994/oILxssbhmy7+wJYM2pY6tcvJtmwYcOwcePGZdwdAADWuqq6chiGDUsZ65lZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtFXDMCx9cNXNSb62fLuzIg5J8q1dvRPAmmHOAVba7jDvPGQYhkOXMnCbYnZ3UFUbh2HYsKv3A1gbzDnASltr847DDAAAaEvMAgDQ1lqM2Xft6h0A1hRzDrDS1tS8s+aOmQUAYPexFp+ZBQBgN7EsMVtVT66qj1fVTVV1V1XdUFV/V1WPXGDsU6rqE1V1e1XdWlUbq+r4mcsfXVUXVdWNVXXntM0PVdVjF9jWkVX1gar6n6r6blVdVVWvqqq9FtnXZ1fVUFU3zK1/cFWdPe3PLVV1c1V9rKqesMA2Tqmq86vq6qraVFWXLHJ7T62qT1bVd6blU1X19EV+ncBWrPY5p6quneaZ+eUZC2zzRVX1xenn+FJV/cYCY8w5sAvtLnNOVe1fVWdNc8kdVXV9Vf1VVa2bG3dsVZ1XVZ+vqnuq6tpFfjePq6qLq+qbVXVbVX2mql6w5F/udthi5O2gH09yZZK3Jbk5yeFJfjfJZVV11DAMX0uSqnpxknOn5Q8zxvWjkuw/s62DklyT5LwkX0/ywCSvSPJvVXXsMAyXT9s6LMklSW5M8vKMn692QpI3JTk0yenzO1lVByX5kyQ3LfAzPDrJryT5yySXJdknyUuSXFJVTxuG4Z9nxj5j2u/Lkuy7pV9KVZ2c5B+TfCDJ66fVL0pyYVWdMgzDB7d0XWBRHeacDyc5c27dl2bPVNWLkrwzydlJPjpt721VVcMwvH1mqDkHdq3dYs5J8mcZ55Mzkmycfo6zknysqn5mGIbbp3EnJHn8NGZIcuBCv5Sq+umMc9dlGeea7yV5VpI/r6r7zc1jO88wDCuyJHnY9At45XR+XZI7krx8O7Z1YJK7krx1Zt2vT9tfPzf2b5J8fQvbeVfGO/u8JDfMXXZQkr3m1u2V8T/CJ+bW7zHz708muWQLt3d+kuuT7Dmzbs8kNyR5/0rdFxbLWlhW05yT5Nok793KbeyV5JtJ3jO3/i8y/tHae2adOcdiWWVLwzln/yT3JHnD3PqTp9t58sy62TnnvUmu3cI235Dk+0kOmFt/aZJLl+t3v5LHzP7vdHrPdPqCJJuSvGM7tvXdjHfyPTPr9plOb50be0sWOJyiqh6X5LQkL13oBoZhuGUYhnvm1t2T5L+S/MTc+k1L3O99knx3GIYfzFz3B0luX2gfgR2yquacJXhsxmdX3ju3/q+THJzk2M0rzDmwKnWbc/acloW2l9ltbuOcc3fGiJ/1f9u5j0uyrJNZVe1ZVftU1UMzvnR2U5L3Txcfm+SLSZ5dVV+ZjsG4pqoWjMuq2qOq9q6qwzM+XZ8k754ZckHGZy/OraqfrKr7V9UvJvnVJG+e29beGZ+VPWcYhmu24efZJ+MfnP9e6nXmvCvJEVX1+1V16LS8LuOjt3MXvyqwNat1zpmcUlXfm46vu2z+2LUkR06nn59bf9V0ep9j8ZbAnAPLqPOcMwzDbRkfLP9OVR1XVQdU1ZFJzkny2SQf245fyXnT6Vuq6rCqOmg6fOqEJH+8HdtbmmV+yn3zsRVDki8necTMZV/M+Gjg5ozHVRyf5O3T2JctsK2/n9nWN5Icu8CYIzJO/JvHbUryugXGvSbj8Sn7TufPy9xhBos8fb4pyeMXGbPFl/ymy5+S5Dsz+3hrkqcs5/1gsayVZRXPOW9N8ryMx5w9K+Nxb0OS02bG/N60bt+56+41rX/tFn5mc47FsouWznPONG7PJH86s70h4/Guhy7yM2/xMIPp8p/LeCjT5u19P8kLl/V+WOY7+RFJHpPk1IwHSt+QZN102dXTD/nMuev8S8ZHNjW3/qemX9AzMz5auCXJhpnLD03yuSSXJ/mlJE9K8gfTL/H0uf8IdyQ5eWbdedlKzCZ5zvSf5qytjFvs+LVjknw7yfuSPHla3p/x5YTjlvO+sFjWwrIa55wt7OeeSa5Icv3Mup0es+Yci2V5l85zzrT+7IyHALwyyRMyHn755YyR/mNb2NZix8w+NMl1Gd+P9NSMz8i+JeOhB89dtvthBe/wg6Y75h3T+UunO/nAuXGvmNYftsi29pn+k1w0s+6cjI+AHjA39vXTHX3IdP5DST447c/m5fyM7w48KMl+C9zeKdM23r2En3OxPyxXZO7NYzPX+c+Vui8slrWwrJY5Z5Ftvnq63QdP539z9vzMuAdO61+6he2YcyyWVbA0nHOOnM6/cG7cQ7OFZ4+nyxeL2QuSfCUzb1id1r8v4yESeyy2j9u7rNgbAIZhuCXjS/tHTKuuWmR4Mj4LuqVtfT/jo5MjZlYfleSaYRi+Mzf88iR7z4x9ZO592W3zcmqSw6Z/nz175ao6IeOdc2GSF29ln7fmqIyPduZdkfHRHbCTrKI5Z6u7Ord/R85dvvlY2S8scXuzzDmwQhrOOUdNp1fM3faXM0b59swRRyX57DAMdy+wjwdnfHC+061YzFbVg5I8PGOxJ2McJuPLXrNOzviS/0Kf/bp5W/sn2TCzrWR8yv6IqnrA3PDHTKc3TqfPTnLc3PLhjI8YjsvMmyJq/MDif8j4dP9pw9LfzbclN2V8CWHe0TP7B+wEq2jOWWh7e2X8HOvrZm730ozz0HPnhp+W8VCBT21pe4sw58AKaTjnbD49em7s+ozPMm/PHHFTkkdNb5if38c7M85lO92yfGlCVV2Y5DMZH1XcmmR9xqfV78m977j7UJJ/TfLOqjokyVeT/HKSk5L82sy23pnxh9+YcaJ/SJLfSvLgjO/g2+wdGf8IXFxV52T8iIwnJXlVkguHYbg+SYZhuGyB/X1+kruGYbhkZt3DMx6O8K2MT+0/uqp+eJ3Z7VTVQ3LvH4yDk2yqqmdN568Ypg9PznhA9h9V1fm59+N3npfk55O8bH6/gKVZzXNOVZ2a5OnT7V+f5EEZPxLwZzO+KpQkGYbh7qp6bcYvSbgx4wePH5/x431+e3qmZvM+mnNgF9od5pwk/57xUwvePAXy5i9NeE3G42jfM7OPhyZ54nT28CT7z8w5XxiGYfMrR+dmfDX7n6rqbRnfo/S06Xb/eHYe26mW6biR0zMeCH1Lxm9/+FLGj6xYNzfu/hnfRfeNjMd7fC7Jc+bGvCDJf2S80+7M+Cjl/CRHLXC7x2S8876e8Q0OV2W8U+5zHOzc9c7Lfb804fn50Xf3/ciyDWOfPzf2uUk+nXsPcfh0klOX436wWNbKsprnnGnMx6fbvHvax49m5gPJ57b54ozHyt2V8Y0YL1lgjDnHYtmFy+4y52R8MPzmaa65I2P8/m2Sh82Ne9Iic86Zc2N/IeOnJ9yc5LaMn8//ksx8ecvOXmq6YQAAaMc3wAAA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQ1v8DAFPABMz2N9AAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1163f4fd0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArMAAAFuCAYAAACSr5B3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAEyZJREFUeJzt3XvMZHddx/HPty2lLS0CpSJWyhKgXCqhmqUWAyEUWgmKiMYE5FKjQREUUQygiAW8ICiIitzxhkCkSg2XUlq5iIBQtkJVUAqWQCFcii22FahAf/5xztphOs/us9vdfZ7vs69XcvJ0z/zmzJmdp799z5kzMzXGCAAAdHTIRu8AAADsLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFuH7cngW9/61mPbtm37aVcAACC56KKLvjTGOG49Y/coZrdt25YdO3bs3V4BAMA6VNWn1jvWaQYAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCm9Iz3/XMjd6FfeO665LzzktOOy056qjkkEOmn6edNq2/7rqN3sN12zKPCbCl1Bhj3YO3b98+duzYsR93B2BSz6qMs9Y/P21Kl1ySnH56csUVyTXX3PDyo49Ojj02Of/85MQTD/z+7aEt8ZgALVTVRWOM7esZ68gswP5wySXJKackl122OmSTaf2nPz2Nu+SSA7t/AFuEmAXY1667LjnjjOSqq5Ldvfo1xjTujDNanXIAsFmIWYB97fzzp1ML1nsa1xjT+Asu2L/7BbAFiVmAfe15z0uuvnrPrnP11clzn7t/9gdgCxOzAPva+99/YK8HcBA7bKN3ADajJz0p+fCHN3ovrnfyyckLX7jRe8G6fe1re3e9a6/dt/sBcBBwZBZgXzviiL273k1vum/3A+Ag4MgsrOAoKDfKqacm73zn3l0PgD3iyCzAvvaUp0xfiLAnjj46eepT98/+AGxhYhZgXzvjjOmbvarWN75qGn/66ft3vwC2IDELsK8dcsj0WbM3v/nug7ZqGnf++dP1ANgjZk6A/eHEE5MLL0xOOCE55pjVY445Zrr8wgun8QDsMTELbEpn3e+sjd6FG+/EE5NLL03OPju5//2TI4+cjr4eeeT057PPni5vErJb4jEBtpwa6/26xSTbt28fO3bs2I+7AwDAwa6qLhpjbF/PWEdmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0FaNMdY/uOryJJ/af7uzIW6d5EsbvRPAQcOcAxwoneeb248xjlvPwD2K2a2oqnaMMbZv9H4ABwdzDnCgHCzzjdMMAABoS8wCANCWmE1evtE7ABxUzDnAgXJQzDcH/TmzAAD05cgsAABt7feYraofqKp3VNXnq+raqvpMVb2+qu6+MOZdVTXWWM5b2t7JVXVeVV1TVVdV1Rur6k5LY7ZX1cur6j+q6itV9emqek1V3WFp3E/u4nZHVX3H0vjHztu8tqo+VlWPW7p822629/Cl8WdW1UXz/bi8qi6oqvve2L9zOJhtlTmnqo6qqmdV1SVV9dWquqyq/rKqti1tc1fbe9rS2B+qqvdU1ZXz8t6qeuiN/1uHg9dmnnPmscdW1R9W1aXzXPLJqnpRVd3gY6+q6siqemZVfXy+L1+oqjdX1eELYzbdnHPYvt7gCrdKclGSFye5PMkJSZ6W5P1VdY8xxqeSPD7JzZeud+8kL0jyxp0rqurOSf4xyb8leWSm/T8rybur6uQxxhfnoQ9PclKSP0rykSTHJ3lGkh3zuMvmcW+Zb2dRJXlTkkvHGJ9fuO3HJnlZkuck+fskD0jy4qqqMcZL5mGfW7G9JPmtJPdJ8raF7f3MvL2Xzn8fRyX55SQXVNW9xxgfWrEdYPe2xJyT5JVJfmS+vR3z/XhWkrdX1T3HGNcs7PeyJyR51NJ9edD85zck+e159WOTnFNVDxljvGXFdoDd27RzTlXVvP0Tk/xGkn9Pcvckz06yfe6NMY+9SZK3JrlDptb5aJLjkpye5NCl/V62sXPOGOOAL0nukmQkefIuxrwqybVJbrWw7pVJvpzkFgvrvivJ15I8b2HdcSu2d/sk1yV59m727b7zvj1hYd1hSb6Y5C+Wxv5ppg8jvskutndUkquSnL20/n1J3re07pgk/5vkORvxuFgsW3VpOOccleQbSX5naeyD5rE/sJttfiLJB5fWvTbJZUkOXVh3aJLPJHndRj9GFstWWjbLnJMpYkeSn1ka+7h5/V0W1j1t7pXb7cX93dA5Z6POmf2v+ec3Vl1YVUcl+fEkbxpjXLFw0alJ/mmM8eWdK8YYn8n0DOZhC+suX97mmJ4ZXZ7p2cuunJkpKF+3sO7emZ6d/NXS2FcnOTbTUde1/GimSP2LpfWHZ/qlWfSVJF+Pc5lhX+s25xw6L8tzxM79WHOOqKr7JLljVs85/zPG+ObCPn4zyTW72h6wVzbLnLPz9ID1zCWPz3Tg7bLsgc0w5xywCayqDq2qw+dD6C9L8vl86+S96GFZHYDfzDTpL7s2yR2r6ohd3P7dknx7pkPsa405MtMv15uXfrlOmn/+29JVPjL/vHvWdmamo7rnLa1/cZIHVtVPV9Utqur4JC/KFLOv2sX2gHXoPOeMMa7O9GT5iVV1/6o6uqpOSvJ7SS5O8va1tpnVcZxMH9Fzp6p6elUdNy+/kWRbprkHuBE26ZzzkSTvTvKM+Tzbo6vqlEynHLx1jPHv83VPSHK7JJdW1Svmc3W/VlVvr6qTd3PXN37OOYCH3HdkOqQ9knw8yd12MfZtSb6Q5LCl9a/PdHj6Jgvrjsn0DGMkue0a2zssyT9kispb7uJ2HzFv54eX1v/avP6IFdsdSZ6xxvaOz/SL+YI1Lv+pTC8d7Px7+VySex2ox8Ri2cpL5zlnvuzQJH+ycB9GkvdnxcuLC9c5Yt63N6xx+YOTXLmwvauSPHijHyuLZSssm3XOSXKzTOetLs4lb05y5MKYUxfmhLfPc8XDkvzLfNsnrHG7m2LOOZAvLT16/sv6ifnOXFBL78pNkqr6ziQPTPKaMcby4fk/yhSIL62q46vq9kn+LMnR8+XXrXHbL0ry/UkeNca4chf7uPMo6rnruUPr8OhMR7//fPmC+d18f5Lp2dsDkzwk05Hfc6vqu/fR7cPBrPuc81uZ3lDxK0nuN9+fY5O8taputsb2fiTJt2X1nHNqplOlzs107u2DMr0h7eyquv8u9hFYn80657xi3q/HZZpLHpdke5K/qaqdHbjz51eSPGSMce4Y45wkP5jkyExv8Fplc8w5G/Ts5RaZSv6lKy57SqZ6v+ca1318kv/O9ZV/Qa4/ifoGb8RK8ruZHvxH72afbpvp3JY/WHHZz2XFM6JMh/O/5Y0bS5d/NMmHVqyvTM/IXr20/vAklyY5ZyMeF4tlqy4N55yT5tv66aX1d57X/+Ia2zw3Uxyv2q8PJnn3ivXvWTVPWSyWvV82y5yTKUZHkgcsrT99Xv/Q+c8737D2tyu2cXGS89fY100x52zISf9jOrH5E0nutOLiM5NcPMa4eI3rvjhTRH53psPepyf5ziQfGGN8fXFsVT09yVOTPHGM8erd7NajMr2st3z+SnL9ubEnLa3fea7sR5evUFX3SnK3NbZ3m/k+fHBx5RjjfzP90txtN/sK7IGGc8495p/Lc8THM/0DeYM5oqbPqD0jyWuX92thmztWrP/gqu0Be28TzTkr55IkF84/d/6/f2mSr+7iLt3giPBmmnM2JGar6jZJ7prkP5fWb88UiKsm9/83xrh2jPGRMcZlVXWPTIfrX7I4pqqemOlluqePMdZzovFjkvzLGOPDKy77p0wfwfXIpfWPSnJFkveuuM6ZmY66vHbFZVdmeoZ1ytI+H57k5CSfXcf+AuvUcM7Z+Xmzy3PEiZmO+KyaI3YVxzu3ea8V609ZY3vAXtpEc87KuSTJ980/Pzvf3tcznQJwn8XTmOY3ht01N4zhZBPNOfv9SxOq6pwk/5zpJOKrMn3m2S9lCr3nLw1/zLz+NWts67syveT/vkwxuD3Jr2Y68fh1C+MenuSFmT5B4B3zeRs7XTXG+JYjqVX1vZmeAT151e2OMb5eVc/I9CUJn830pQmnZXoD1y/MR1QXt3d4pg80fuu4/gOOF7d3bVW9IsnPV9UVmU/ETvLzmd7l96RV+wHs3laYczJ9aPrFSZ5fVbfM9V+a8OuZXn5c9Y/HY5L861j7C1f+OMnvV9Vrc/3HDD4m03l2v7jGdYDd2ORzzs4vLPjLqvrNJP+RKU7PyvQZsOcsXO+sTEds31JVz8/05q6zMr0atCqWN8+cs7/OF1k4N+Kpmb4Z48uZTiz+WKY3PW1bGneTTJ+P9qZdbOs2mULyS5ke5I9m+sdg+d2Af55vfdfe4vKuFdv9w0wfiXWb3dyXn01yyXzbH0/y+DXGPWy+rR/bxbYOyxSvH05y9Xzf35XkjP39mFgsW3nZKnNOpjd7PX+ea76a6R+ev87Ch5wvjP2e7OYD2udxj0zygUyvDl05//cjNvoxs1g6L5t9zsn0kVuvSvLJTJ+g9MlMbwo7fsXtn5LknfP9+O8kf5fkTivGbao5p+YbAwCAdnzrCwAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABt/R+rNyLGCS989AAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1157b2ba8>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArMAAAFuCAYAAACSr5B3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAFJpJREFUeJzt3X+w5Xdd3/HXOwkkS5M0QAgdRnGThmiNYDqkMdRqEjFJ1YG2EgQtSJW2irUKFQJMO0NSqyPBSGppRaUMoqA0Kk4plW4KoaMOP7JBgoEJIQYQHBICGHZDfpCQb//4fK+cPdzdvbve3T3v3Mdj5js395zP+d7vOffms8/zvZ9zbk3TFAAA6OioI30AAABwsMQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG0dcyCDTz755Gn79u2H6FAAACC5/vrrPzdN02M2MvaAYnb79u3ZuXPnwR0VAABsQFV9cqNjLTMAAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW8cc6QMA2J/L3n1ZLjv/siN9GBv34IPJjh3JFVck731vcu+9yXHHJeeem1x6aXLRRclRPc4ltHvs1/HCFyYf/OCRPoo9nXVWctVVR/oo4KGhx2wKbGmX/7/Lj/QhbNzNNyennpo885nJtdcm99yTTNP4eO214/LTThvjGmj12ANbkjOzAJvl5puTc85Jdu0aAbueu+5KvvSlMe7970/OOOPwHuMW5AwoPLQ5MwuwGR58cCwf2FfIrpmmMe6ii8btADhoYhZgM+zYkXzhC/sP2TXTNMZfc82hPS6AhzgxC7AZrrgi2b37wG6ze3fyylcemuMB2CKsmaWNVXtFslcjs4f3vvfw3g6AJGKWTXYog/OWW8ZrZw7E8ccnp59+aI4H9nDvvQd3u/vu29zjANhixCxtHEyUOnvKYXPccePttw7Uscdu/rHAEbBqvz1L/BuwVYhZNpVJgy3r3HPH+8gezO0AOGhiFmAzXHppct11B7YW5vjjk5e+9NAdExxGTmZwpHg3A4DNcNFFyaMfnVRtbHzVGH/hhYf2uAAe4sQswGY46qjxXrMnnrj/oK0a43bsGLcD4KCZRQE2yxlnjD9R+/jHJyecsP6YE04Y1/tTtgCbQswCK+8V573iSB/Cxp1xRnLrrcnVVycXXJBs2zbOvm7bNj6/+upxfZOQbfXYA1tSTRv904tJzj777Gnnzp2H8HAAANjqqur6aZrO3shYZ2YBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQVk3TtPHBVXck+eShO5xD7uQknzvSBwFsOeYe4EjoPPd8wzRNj9nIwAOK2e6qauc0TWcf6eMAthZzD3AkbJW5xzIDAADaErMAALS11WL21470AQBbkrkHOBK2xNyzpdbMAgDw0LLVzswCAPAQsukxW1UXV9W7quq2qrqvqj5dVf+jqr55ady3V9WOqvpsVe2uqg9U1Y+us7/jqupVVfWZqrqnqt5TVd+5zrijqurlVfWJqrq3qm6oqmfs5RgfWVVXVdVfLBzjG9YZt62qLquqj83jbq+q/1VVD18Y86Kquq6qPj9/3Vuq6sqqevQ6+3tBVd007+svqupnq+phG3xogX1Y5bmnqs6vqmkf27l7uU+nVdXd85jTl64z98AKWOW5Zx73iKq6vKpunvf3qap6Y1VtXxhzdFW9eL4fty8c3/Or6qil/T2nqv6kqu6Y7+8nquq/V9Xj1/nal1TVn87Hd1tVvaaqTjiwR3gDpmna1C3JDyZ5VZJLkpyX5LlJPpxkV8Z7hiXJk5Lck+TaJP8kyYVJfjXJlOQFS/t7U5I7k/yrJE9N8vvzbc9aGvdzSe5L8uIkF8z7ezDJ9y6Ne+R8PDcm+eEk35nk2Un+y9K4hyV5d8b76v74PO4ZSV6bZNvCuJ9N8vIkT0tyfpIXJfmrJB9IctTCuJfPx3PlfH9fkuTuJK/b7O+BzbYVt1Wee5KcmOTcdbYPJ/lMkqP3cp/eMV8/JTl96Tpzj822Atsqzz3zuDfP/8+/ZB73vCSfSPLnSY6fxxw/f81fno/vqfOc8WCSVy3t76eT/Md53HlJ/nWST83bCUuPy5TkDUkuzmipzye5ZtO/B4fpG/2N8x36mfnzn0/y5bUHcWHce5K8Z+Hzb51v9yMLlx2T5KNJ/ufCZafM39DLl/b3ziQfWrrstRmBeuJ+jvll8w/i1x/E/f2x+bifPH9+XJLdSd6wNO7F8w/KmYfzfzybbatsqzT3rHNs37DePxQL1/9QktuTvDDrxOxebmPusdlWYFuVuSfJI5I8kOTnl8b94/nrXDx/fnSSR61zP16f5N4snMTby/29eN7fMxYuuyXJu5fGXTKP+9597e9At8O1Zvbz88cH5o8PT3J/xjONRV/Mnksfnj6Pe8vaBdM0PZDkd5JcXFXHzhdfPO/zt5b291tJnlhVpyZJVf2tjLOxr5umadd+jvknklw9TdOn9jNuPcv391synvX84dK4dySpJP/0IL4GsH8rMffsxXMz/v//jeUrquqRSX4pIzrv3Mc+lpl7YDWsytxz9LwtN8/avHLU/DW+Mk3TF9a5H9clOTbjL4ntyx73t6pOTvJ3s/7ckyT/bD/7OyCHLGbn9RcPr6onZJz6vi3Jb89Xv2H++MtV9biqOqmq1k6nv3phN2cm+fg0TXcv7f7DGd/E0xfG3ZfxLGB5XJKsrVt5cpJtSW6vqt+d147cVVV/sPiPzrzu4+uT3FpVv15Vu+b1Hu+sqrP2cn+PmdelnJvk8iTvnKbphvnqr8wfv7x0s/vmj9+y3j6BA7eic896fjjJB6ZpunGd665IctM0Tb+5j9snMffAqljFuWeapt1JfjPJT1XVBVV1fFWdmbEs4oaMM7n7cl5G+H5mL/f32Kp6UsaT748k+T/z1Xube+7PODO7qXPPoTwz+76MB/rmjLUi3zVN02eTZJ68z89Yb/GXGeu8/muSH5+m6XcW9vGo+bplX1i4fu3jndN8Dnsf4x43f/zFjAf66RlrPf5+kncvLEpeG/fSJKdlrKn9wSSPmcftsci5qo7P+AZ9KeNXBp/Kns86PpbxK73lF3k8Zen4gL+5VZx79lBVT0nyhKx/VvY7MkL3J/Z6D7861twDq2NV554fSfLWJO/KWHZ0Y8brgi6cpmk5Nv9aVV2c5AeS/OJ8dnjZ7RlLEG7IOFH43dM03Tvf379Kcke+du75tozfCm3q3HMoY/a5GXfihzJOb1+z9sq5+VnL72U8g3haku/OWMv62qr654fwmNbu761Jnj1N0zXTNL0545v1+CTPWRp3d5KnTdP0v6dpemuS78v4hv2bpf3eneQfJPmOJD+V5Kwkb6uqY5Jkmqa7Mtad/GRVPXt+RnZBxhqar2T8YwNsjlWce5Y9LyNC37x4YY13SvnVJK+epukjG9iPuQdWx6rOPf8po29enK++QO3RSf5wXn75NeZ3YvjtjBesvXIv+31qkn+Y5PlJTsq4vyctXP+fk1xSVT9ZVY+qqicn+ZUcirnnMC2EPinjNPVr58+vzngV3cOWxr0pyecyvxI3Y83IR9fZ3w9knKY+c/78lRnPDmpp3DnzuO9bWqB85Tr7/GKSX5v2XLj9e+uMuyHJjv3c3/Pn2z974bITM16R+OB83T0Zryy8I8nrD8f3wWbbatuqzD1L1x2bcfbkretcd2nGr/MeNx/7SRlnaKeM3yCdsJ/7a+6x2VZgW5W5J2M5wpTk+UvjnjBf/tPrfK3TMs4eX5/9vFh+6TYPJnnZwmVrT84fmL/W/UmuSrIzybs28/E+LC8Am6bpzox1HWtrPZ6Y5IZpmu5fGvr+jGcLp8yffzjJqVX1iKVx35yxDuOWhXHHZiw2Xh6XjHUca+P2Ze2Zwq352kXa643bm53zx79+X8hpmnZN0/T9SR6b8euHU5K8MWNR9R/vZ3/AQVihuWfR0zPeIvBrlhjMt/s7+eqvIdd+FZmMt9z6o3Vus8jcAytgheaeJ84fr1s6vo9lxPbfW7y8qr4uYx3trox3Otjfi+XX9ndrxpP0xbnny9M0/VjGXPOtGXPQSzJCelPnnsMSs1X12CTflPGsJBmLos+qhT8+MPu2jGcaa2s+3paxruOZC/s6JsmzMs6Orr2I4R0Zxb98qv45SW6cpunjSTJN06czJvsLq6oW9vmUjLMX183j7k/y9iT/aPEU/LxW9puy9EOxjvPmj3++fMU0TXdM0/Rn01iU/aKMZ2RX72d/wEFYlblnyfMy/r9/+zrX/ULG+0Aubmu/4ntOkn+5t/s6M/fAClihuee2+eM5S8d3RsbZ479cuOwxSf7v/OmF0zR9bkN3dtz2zIwoX2/uuXOapg9N490Snp8R4a/f6L434pjN3FmSVNVbM84gfCij7M/ImDgfyHgD3iR5TcYk+raq+m8ZZ0GfnvEiq1dP84LkaZr+tKrekuSqGn+t5uNJXpDk1Cx8A6dp+mxV/VKSl1fV7vnrPyvJd837XfSyjFfb/W5VvS7jRV0/l+Sm7Ll+7RUZz5jeXlVXZrxf4ysynsm8Zr6vfzvjB+pNGS+0mDJ+YP5dxnKE3194XJ6VseD5oxlnZb4/49cGz5j/cQH+BhrMPamqUzKWO/3KOmdoMk3TTRlz0eJtts//+b5pmm6ZLzP3wIpY8bnnjzLmhCtrvOXfzozXCP2HjOWVvzHfh20ZbbQ9yY8m+br5LO2aj6ydpa2qP854QdlNGSH+pCQ/k+TTSX594XG5MONdC27MaKiLMpZN/dtpmj6x0cd3Qw7BOpGXZqyzuDPjxQkfzVgzsX1p3Pdk/IWtOzJeXffB+U4evTRuW8ZbPtw2P2jvS3L+Ol/36IxvziczXk34oSSX7OUYvyfj7Oq9Ge+N9sYkj11n3DkZi5/vzvim/0EW3rg8X312cXPGq4m/mPFD8++ztLYt4x+PP5v3tSvJjiTfvtmPv822Vbcmc8+LsvBHDTZ4v/5Flv5ogrnHZludbdXnnowzpldmPPG9J+NdT96S5BsXxmyf55m9becvjL1ynlN2JbkrY0nDq5KcsvR1z8tord3zPPUnGS+q3/TvQc1fEAAA2jlcfwEMAAA2nZgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANDW/wfCRAVSM25adgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x115fedd30>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArMAAAFuCAYAAACSr5B3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAFDZJREFUeJzt3X/0ZHdd3/HXewkhxoTyI1CBQ1yrRAOC2G4x1BzEcLJBak8r2goUBWk5YKGH2HpItQRii2JAbFTkV9UmFfEYrVSOhRhEo6fVmG4EqYBZfpmQYAqYQAI1CzGf/nHvF8Zhvpvvbr6733nvPh7n3DM7dz5z5853Np88v7N37tQYIwAA0NGund4BAAA4XGIWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLZOOJTBp5122ti9e/cR2hUAAEiuvfbaT44xHrSVsYcUs7t3786+ffsOb68AAGALqur6rY51mAEAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErPAWrroqot2ehfuubvuSq64IjnnnOTkk5Ndu6bLc86Z1t91107v4SE5Jl4T4JhTY4wtD96zZ8/Yt2/fEdwdgEn9SGW8bOvz09rZvz8599zklluSz3zmS28/5ZTkgQ9MrrwyOeOMo79/h6H9awK0UVXXjjH2bGWsd2YBttv+/cnjHpd89KOrQzaZ1t9wwzRu//6ju38AxxAxC7Cd7ror2bs3ue225O7+5WuMadzeve0OOQBYF2IWYDtdeeV0aMFWD+EaYxr/jncc2f0COEaJWYDt9MpXJrfffmj3uf325OKLj8z+ABzjxCzAdrr66qN7P4DjnJgF2E533HF49ztwYHv3A+A4IWYBttNJJx3e/e5zn+3dD4DjhJgF2E5nnXV07wdwnBOzANvpxS+evhDhUJxySnLBBUdmfwCOcWIWYDvt3Tt9s1fV1sZXTePPPffI7hfAMUrMAmynXbumc83e9753H7RV07grr5zuB8AhM3sCbLczzkiuuSY5/fTk1FNXjzn11On2a66ZxgNwWE7Y6R0AWOVl3/Kynd6Fe+aMM5IPf3j6Zq+LL57OI3vgwHTWgrPOmo6RPffcVu/Itn9N7qHzz0/e/e6d3ou/6bGPTS65ZKf3AnaWmAXW0kVPvGind+Ge27UrOe+8aTkGHBOvCXDMEbMAsAXeAYX11OfftwAAYImYBQCgLYcZAMBRsm4fIvMBMo4F3pkFAKAt78wCwFHiXVDYft6ZBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbdUYY+uDqz6R5PojtztH3WlJPrnTOwEcN8w5wNHUec75yjHGg7Yy8JBi9lhTVfvGGHt2ej+A44M5Bziajpc5x2EGAAC0JWYBAGjreI/ZN+70DgDHFXMOcDQdF3POcX3MLAAAvR3v78wCANDYEY/Zqjqvqn6nqm6uqgNVdWNVXV5Vj1wx9ilV9ftV9Zmquq2q9lXVOQu3/72quqKqbqqqO+Ztvq2qHr+0nSdV1Zuq6kNV9Vfz5euq6sGb7OOZVfWrVfXJefx1VfWigzynp1XVqKobV9x2clX9p3kfD1TV/6mqf77Jdv5JVb1rfi7XV9VLqupeB/t5Age3znNOVV00zx2rljsWxp1RVT9VVe+Z9+0vquqtVfUNK56DOQd20DrPOfPY06vqsqq6YR67v6peXlVfvmLsc6vqz+bncV1VPX/FmHtV1YVV9ZF53Aeq6vxNfjZnV9UfzI97c1X9ZFV92VZ/tlt1wnZvcIUHJLk2yWuTfCLJ6Un+XZKrq+rRY4zrk6SqnpfkNfPyHzOF9mOTnLywrfsl+WCSS5P8RZIHJ/mBJL9XVWePMa6Zxz0/ySlJXp7kw0kekeRHkpxXVY8ZY3xmY4NVtSfJ7yS5Ksm/TPLpefwpq55MVd0vySVJbt7k+f56kscneUmS65I8NcmbqqrGGG9a2M55Sf5bkp9P8m+SfGOSH0tyapILNtk2cPfWec75uSRXLO3vl8/r3rqwbm+Sb01yWZI/nvfjxfNzOHuMce3CWHMO7Ky1nXPmYP3tJPdOcmGSG5L8/XnsI5J898YDV9Vzk7whySvm+zwpyWvnueR1C/v42iTPnp/DH2Waq36iqk4ZY7x8YXuPSfKOJL+V5NuTfFWSVyV52OLjbosxxlFfknxtkpHk387Xdyf5qyTnH8a2Tk1yIMnPLKx70IpxT5gf8zkL63YleV+StxzC471xfmEuTXLj0m1nz4/x7KX1v5nkY0nutbDuXUl+b2ncS5N8LslX7MTrYrEcq8u6zDmbbO975nH/cGHdaZk/07Cw7m8luTXJf11YZ86xWNZwWZc5J9MvxiPJ3qWxP57kziQnz9dPSPLxJJctjfuFTF+6cO/5+ulJ/jrJRUvjXjM/vwcsrHtLkg9s3Hde973z/vzd7fx579Qxs385X945Xz4nyV1JXn8Y2/psphd5Y1sZY3xixbj/PV8+bGHdE5OcmeQnt/JAVfXNSZ6Z5AWbDDlrvnz70vorkjxk4/aqenim38betDTuFzP99vRtW9kfYMvWZc5Z5VlJ/m+mX5I3tvfJMc/8C+s+nWT/0vbMObCe1mXOOXG+vG1p7KcyvaFX8/XHJ3lQVs8RD8z0i3OSPG6+36o556TMc0lV3TvJk5NcPsb4/MK4yzP9Av2PV+z/YTtqMTsfY3FiVT0i09vYNyf55fnms5P8WZKnzcd93FlVH6yqldFYVbuq6t5VdXqm3waS5D/fzS58y3z5/oV1Gy/OSVV1dVV9vqo+XlU/vXxMx/zCvDHJq8YYH9zkMf56vvzc0voD8+XXz5ePmi//dHHQGOMjSf5fki85zgY4NGs65yxv9+GZ/onul8YYd242bh77gExzyOL2zDmwJtZ0zvntTO+OXlxVj6yqU+ZjdF+U5PVjjM/O41bOEUneO19uzBFbnXO+OlPcLs85dyT5ULZ5zjma78z+UaYnuz/JY5KcM8b4+HzbQzMdu/GqTG997810nMVravUHsTbK/vok35nkKWOM9232wFV1aqbjXN+f5L8v3PTQ+fJXklyZ5Nwkr8x07OyblzZzQZL7ZDqWZDPXzZdnLa3fOHD7AUuXt67Yxq0LtwOHbx3nnGXPzDQPX7aF5/Mzmd5FuWRhnTkH1sfazTlzPJ6daZ55b5Lbk7wz06FIL1zYxGZzxC1Lt2/HnHNLtnvOOYrHj5yZ5JuSPD3TgdI3Jtk937Y/0zEUT126z9sz/WazfOzY38l0APNTM70on0qyZ5PHPSHTi3Z7kscs3fbG+XF/emn9BfP6M+frX5PpWJAnL4y5NF96zOwJmY7B/UCmF/b+Sf7FfN+R5IJ53DPm61+3Yn9vTPLzR+t1sViO1WUd55wVY9+f5I+38Fx+KCuOvzXnWCzrs6zjnJPp3dHfzRShz8x0XO0PZjrs4HUL43543r+TVmx7JLlwYd07Mh1fe16mD6x9R6ZAHRvbTPIP5utPXrG//zPJO7f1Z79DL/j95hfm9fP1P5yf9KlL435gXv/Qg2zrxPkvyRUrbtuV5JeS3JHkSStuf8W8/X+0tP4b5/XPmK+/Lcn/mPd7Y3lzkpvmP3/Zwn0fmeTd8/3H/Jf0/PnP3zuP+bb5+uNX7NNnMx3KsOP/YVosx8qyLnPO0tjHzY/1orsZ9/x53L/f5HZzjsWyZsu6zDmZPuMzknz10vrnzuu/Yb7+/fP1hyyNe/C8/gUL6x6W6QxQG3POp5M8b/7zS+cxZ87Xn75in96X5Fe38+e9Ix8AG2N8KtOpJ75mXvXegwxPpoOmN9vW55K8Z2Fbi16f6fQPTxtjvHPF7Vt93EcmeUqmt8s3lqdn+meDW7Nw6MEY431jjMdmOgXF1yd5eKbTayTJ/1p63I1jVJIkVbU70yk6Nv2nBODQrdGcs+hZST6fLz2k6Quq6nsynQbn1WOMH91kf8w5sGbWaM55dJJbxxgfWlq/cYqvM5f271FL4zaObf3CHDHGuGmM8cRMUfvoJF+R6RfqZHrXNZmOiz2wvL2qOinTu87bOufsSMxW1d9O8nWZnmwynb4hmd6yXvTkTP+Uv9k5XVNVJyfZs7CtjfWvznTs6/eNMTY7Zu3tmX7Yqx43SfbNl0/L9CGNxeW3Mp2u4lvzxYOzv2CM8edjjI2/HC9McuXGX6Yxxg1J/iTJ8onNn5npf27LnxIE7oE1mnM2xp6YaV55+1j9qeRU1Xck+S9Jfm6M8YMH215izoF1skZzzs1J7l9VyyH8TfPlTfPlH2ZqmlVzxC354i/GXzDG+NgY408zvSt8fqYPuF013/a5TGc4+GdVtfidBt+V6fNHb802OuJfmlBVb8l00u/3ZDpG44xMb6vfmeTV87C3ZTqm4w1VdVqmEwD/00wHSH/fwrbekOmHui/TD/0rM03aD8l0rsaNcRdkOin4LyT5QFUtHqj8iYUJ/i+r6hVJLqyq2zJ9ecKeTOdevGzMZy0YY1y94nk9O8mBMcZVS+t/KNMB2x/LdD62F8yX37y0iR9O8pvzc/rlTIc2vCTJTx3sLzVwcOs85yz49kwfgFj5wa+qekKmeeFPkly6tL0DY4x3LYw158AOWvM559J53Nuq6kczfWnCnkxfoHBt5kgdY3y+qi7M9CUJN2U6C8I5mU4p9q/nON147O/PFLAfyfSu7LMyfcjsSWOMxXeYL0pydZLLq+pnM51r91VJfm38zS9+ueeOwnEjF8w/sE9lOgXMdZlOWbF7adx9k/xspvMtbryl/oylMc9J8geZzt+2cXqHNyd59NK4q/LFYzmWl0uXxlamF/qD+eInB/9DFk7yu8nzujRLHwCb17880wt8YH4ulyV5+CbbeGqm/1kdyPQX7KVZOMm5xWI59GXd55x5/G/M2zxxk+dw0UG29+dLY805FssOLus+52Q6VODyJB/N9OHQ/Ul+Isn9VzyX5823H8j0wdJ/tWLMC+fneEem8P71JI/a5GfzhEzv+t4xP+9LMn9Rw3YuNT8YAAC0s1PfAAYAAPeYmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0Nb/B+Bvc1JRjJBrAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1165b8f28>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArMAAAFuCAYAAACSr5B3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAEjhJREFUeJzt3XuMZnddx/HPb7u23dJW7CXIfUVckQCpstaiSOtCdwUiahFFILRClMRqrCIUFKEElIsoNKKhkpRLEBMbJMqlWESQaKh0CytggXVbWihI6SVNL0Kh7M8/zqmdTueZ287l+c68XsnJs/M8v+ecM8+z+9v3nDlzpvXeAwAAFW1Z7x0AAIDlErMAAJQlZgEAKEvMAgBQlpgFAKAsMQsAQFliFgCAssQsAABliVkAAMoSswAAlLV1KYNPOOGEvn379lXaFQAASC6//PIbeu8nLmbskmJ2+/bt2bt37/L2CgAAFqG1ds1ixzrNAACAssQsAABliVkAAMoSswAAlCVmAQAoS8wCAFCWmAUAoCwxCwBAWWIWAICyxCwAAGWJWQAAyhKzAACUJWYBAChLzAIAUJaYBQCgLDELAEBZYhYAgLLELAAAZYlZAADKErMAAJQlZgEAKEvMAgBQlpgFAKAsMQsAQFliFgCAssQsAABliVkAAMoSswAAlCVmAQAoS8wCAFCWmAUAoCwxCwBAWWIWAICyxCwAAGWJWQAAyhKzAACUJWYBAChLzAIAUJaYBQCgLDELAEBZYhYAgLLELAAAZYlZAADKErMAAJQlZgEAKEvMAgBQlpgFAKAsMQsAQFliFgCAssQsAABliVkAAMoSswAAlCVmAQAoS8wCAFCWmAUAoCwxCwBAWWIWAICyxCwAAGWJWQAAyhKzAACUJWYBAChLzAIAUJaYBQCgLDELAEBZYhYAgLLELAAAZYlZAADKErMAAJQlZgEAKEvMAgBQlpgFAKAsMQsAQFliFgCAssQsAABliVkAAMoSswAAlCVmAQAoS8wCAFCWmAUAoCwxCwBAWWIWAICyxCwAAGWJWQAAyhKzAACUJWYBAChLzAIAUJaYBQCgLDELAEBZYhYAgLLELAAAZYlZAADKErMAAJQlZgEAKEvMAgBQlpgFAKAsMQsAQFliFgCAssQsAABliVkAAMoSswAAlCVmAQAoS8wCAFCWmAUAoCwxCwBAWWIWAICyxCwAAGWJWQAAyhKzAACUJWYBAChLzAIAUJaYBQCgLDELAEBZYhYAgLLELAAAZYlZAADKErMAAJQlZgEAKEvMAgBQlpgFAKAsMQsAQFliFgCAssQsAABliVkAAMoSswAAlCVmAQAoS8wCAFCWmAUAoCwxCwBAWWIWAICyxCwAAGWJWQAAyhKzAACUJWYBAChLzAIAUJaYBQCgLDELAEBZYhYAgLLELAAAZYlZAADKErMAAJQlZgEAKEvMAgBQlpgFAKAsMQsAQFliFgCAssQsAABliVkAAMoSswAAlCVmAQAoS8wCAFCWmAUAoCwxCwBAWWIWAICytq73DsA0OuecZN++9d6LezrppORNb1rvvQCA6eLILDCVzvvYeeu9C8t38GDyoQ8lu3YlRx2VbNky3O7aNdx/8OB67+GylH5PgA2r9d4XPXjnzp197969q7g7AIP2ypb+isXPT1Nj//7k9NOTm25Kbrvt3o8ffXRy/PHJJZckO3as/f4dgrLvCVBOa+3y3vvOxYx1ZBZgpezfn5x8cvKVr8wdsslw/5e/PIzbv39t9w9gA3LOLMBKOHgw2b07ueWWZKHvePU+jNu9O7nqquE0BNaE8+Fh4zGDAqyESy4ZTi1Y7KlbvQ/jP/zh1d0vgA3OkVmAlfD61ye33rq059x6a/K61yV79qzOPnEv03IEdOYR4n37ktNOW9fdcXSY0hyZBVgJl166ts8DIIkjswAr41vfWt7z7rhjZfeDEhwFhZXjyCzASjjyyOU974gjVnY/ADYZR2YBVsIppyQf/ejyngeZvistOI+WKsQsK8IkzKb34hcnl102+fqyczn66OTcc1dvnwA2ATFLCQcOLK0RDhxY/bgWzNzD7t3Db/a6/fbFXZ6rtWH86aev/r5RgvkElkfMsiJWexKetiO/cC9btgzXmj355IV/cUJrybHHDuP9wgSAQyJmKcERC0rYsSP55CeHo7Q33TT3dWePOSY57rghZHfsWPt9BNhgHBIAptIrTn3Feu/C8uzYMfyK2osuSn7mZ5Jt24ajr9u2DR9fdNHweMGQLfueABta64v91YtJdu7c2ffu3buKuwMAwGbXWru8975zMWMdmQUAoCwxCwBAWWIWAICyxCwAAGWJWQAAyhKzAACUJWYBAChLzAIAUJaYBQCgLDELAEBZYhYAgLLELAAAZYlZAADKErMAAJQlZgEAKEvMAgBQlpgFAKAsMQsAQFliFgCAssQsAABliVkAAMoSswAAlCVmAQAoS8wCAFCWmAUAoCwxCwBAWWIWAICyxCwAAGWJWQAAyhKzAACUJWYBAChLzAIAUJaYBQCgLDELAEBZYhYAgLLELAAAZYlZAADKErMAAJQlZgEAKEvMAgBQlpgFAKAsMQsAQFliFgCAssQsAABliVkAAMoSswAAlCVmAQAoS8wCAFCWmAUAoCwxCwBAWWIWAICyxCwAAGWJWQAAyhKzAACUJWYBAChLzAIAUJaYBQCgLDELAEBZYhYAgLLELAAAZYlZAADKErMAAJQlZgEAKEvMAgBQlpgFAKAsMQsAQFliFgCAssQsAABliVkAAMoSswAAlCVmAQAoS8wCAFCWmAUAoCwxCwBAWWIWAICyxCwAAGWJWQAAyhKzAACUJWYBAChLzAIAUJaYBQCgLDELAEBZYhYAgLLELAAAZYlZAADKErMAAJQlZgEAKEvMAgBQlpgFAKAsMQsAQFliFgCAssQsAABliVkAAMoSswAAlCVmAQAoS8wCAFCWmAUAoCwxCwBAWWIWAICyxCwAAGWJWQAAyhKzAACUJWYBAChLzAIAUJaYBQCgLDELAEBZYhYAgLLELAAAZYlZAADKErMAAJQlZgEAKEvMAgBQlpgFAKAsMQsAQFliFgCAssQsAABliVkAAMoSswAAlCVmAQAoS8wCAFCWmAUAoCwxCwBAWWIWAICyxCwAAGWJWQAAyhKzAACUJWYBAChLzAIAUJaYBQCgLDELAEBZYhYAgLLELAAAZYlZAADKErMAAJQlZgEAKEvMAgBQlpgFAKAsMQsAQFliFgCAssQsAABliVkAAMoSswAAlCVmAQAoS8wCAFCWmAUAoCwxCwBAWWIWAICyxCwAAGWJWQAAyhKzAACUJWYBAChLzAIAUJaYBQCgLDELAEBZrfe++MGtXZ/kmtXbnTVzQpIb1nsngE3DnAOslY0y3zy0937iYgYuKWY3itba3t77zvXeD2BzMOcAa2UzzjdOMwAAoCwxCwBAWZs1Zv96vXcA2FTMOcBa2XTzzaY8ZxYAgI1hsx6ZBQBgA1i1mG2t7Wmt/Utr7euttTtaa9e21v6utfbIOcY+pbX28dbaba21W1pre1tru2Y8fkxr7Q2ttY+Nj/fW2mlzrOes8bFJy/eP4+7fWnvNuJ2bW2vXt9Y+0lp7whzrfG1r7TPjuG+21r7QWnt5a+2oGWO2L7DdZ47jDmut/f74ulzXWru1tfap1trzW2u+sIBlmub5Zhz7ttba58f13dZa+8/W2m+31g6btc4zW2vvaa1dM67j7RM+3/m2+5IZ4369tfbB1tpXW2u3t9Y+11p7UWvt8OW90kCyoeaco1prbxzniDtaa59trT171pipb5ytq7ju45JcnuSvklyf5CFJXpLk0tbao3vv1yRJa+0FSd48Lq/KENgnJTlqxrqOT/K8JJ9K8uEkZ0zY5geSPG7WfS3J+5Jc1Xv/+njfY5P8SpK3Jbk0yeFJfjPJx1prT+u9v3/G848dx30xyR1JfjLJH47r+PlxzP/Msd0keXWSxyf5p/HjbUleluSdSc5PcluSpyR5a5JHJHnRhM8LmN80zzfJ8G//L5JcmaQn2ZNhDnh4kt+ZMe45SU4ct/uMeT7fueabs8fn/+OM+14+ruvCJDdmmI9eleTkBdYPzG+jzDl/P67zZRk654wk72qttd77u8Yx0984vfc1W5L8cIYX9YXjx9uTfDPJOQs8r83485PGdZy2yG3+9Dj+7Bn33TfJ1lnjtmZ4Iz++iHW+ZlznCfOMOSrJLUkumnHfYUmOm2PshUm+lWTbWr4fFstGXqZlvpln7N8muXXWfVtm/PnaJG9fwud7IMlls+47cY5xLx/38WHr/R5ZLBtpqTbnZAjRnuSsWePen+RrSQ6bZ11T1Thr/a3tG8fbO8fb5yU5mOQt8z2pj6/GMp2Z5NsZ3sS71ndz7/3OmYPGj/cleeAi1jn785jLGUmOSfKOGdv4bu/9pjnGXpbkiAy/tQNYGVMx38zjxsyaQ3rvB5ez0dba45P8YGbMN+P6rp9j+GXj7WLmOmDxqs05p4y3F88a96Ek95/x+FymqnFWPWbHcygOb639UJILknw9d7/oj0/yhSTPbK1d2Vq7s7V2oLV29gpte1uGb6W9f8ILPHPs4RkOo39+wuNbW2tHt9aelOT3klzYe795nlWemeQbGf5SLOTUJDdnOJQPLNM0zzdtsLW1dt/W2tMzzBF/vhLbztL+Qzs1w3+w+1do27BpFZ9zvjvefnvWU+8Ybx81z+anqnFW85zZu/xHhvNLk+HbYLt6798YP37AuPxpkj/IcG7HM5K8ubW2tfd+/iFu+xcynPP6joUGJjkvyYOSPHv2A621RyX57Iy73pnkNyatqLX2wCS7kpw/+wjwHGP3JPnlJH+00FhgQdM83zw1w7ltyfCtvdf23l91iNtMa+3IDJ/HB3rvNy4w9jEZzpe7sPd+3aFuGyg953xxvD0l9zw6e9f5scfNtdKpbJw1OIfkR5L8RJJfzXCy9LVJto+P7c/wAp8x6zkXZ/jqps2xvkWfT5LhK4brMuv82DnGPSvDkYpXTnj8yCQ7M3x18dIM54n8zTzre8m4j49ZYLuPTHJTko8stI8Wi2XhZZrnmyTfO84jT0zyJxmOhvzxPOtb1DmzSZ457uPTFhh3/yRfSvK5JMes93tlsWyEpfKck+GA5hVJ/jtDwH5fkudnOM+3Jzl3wnqnrnHW+k2/b4ZDzW8ZP/7E+IIcM2vc7473P2C5b/Q4cd+Z5I0LjPu58Q1+6xI+j7PGfThlwuNXJPn0Aut4WJKvjn/5j13L98Fi2QzLNM43s55zXoZv8z1wwuOLjdkPZvh23/fMM+b4MWKvnrQ9i8VyaEvFOWcMzn3jNnuGyD5n/PNzJ6xn6hpnTX8ArA/nmB7IcGmIJPmvBZ6yrB+GGD0nw0/WTTzFoLX2xCQXJXlvkhcsYd17x9uHz36gtfbjGb5Sm2+7D8rwlcotSfb03m9ZwraBRZi2+WYOezP83MIPLHej43Uldyd5d+/9OxPGHJvh0jnHJ3lS7/2ry90eMFnFOaf3fkXv/aTxvkcleXDuPrf132evYFobZ01jtrV2vwzXGrtyvOu94+2eWUN/Nsm1/Z7XTFuq5yb5TO9934R9eVySf8jwgj+nL+2niE8db6+c47EzM3y19O4J2z0xyT+PH57ee79hCdsFFmma5psJTs1w9OOqQ9juvP+hteGXu3wgw39Uu3vvBw5hW8A8Ks85vfere+93xfdvJbmk916mcVbtB8Baa+/NcAHgz2So8x0ZDq3fmeTPxmEfTPLRJBe01k7I8AI/I8ORhl+btb4nJ7lPkkePd506Puf23vvFs8b+WIavMF44Yd8ekWGCvyHDidmPba39/+O990vHcY9J8oYMR2+vynBpiSdk+AGKi3vvn5i13sMznL92cb/7BPCZj2/LcIRke4ZLdjxo/ArmLlc4SgtLN+XzzVPH9b8vyZczXM7myRl+iPSC3vvXZox9ZIZv+yXDBcgf2lr7pfHjf+33vtTWc5N8tvf+6QkvzXuS/FSGOes+rbWZl9q5co71AYuwgeaclya5JsN1ZR+S4ZevPCTDvDF7vdPbOKt47si5Gc6VuDnJ/2b4qbkLMp4YPWPcsUn+MsNJzN/O8BfjWXOs7+rcfU7HzOXqOcaen+Q7Se43Yd/OmrCunvGSb+O4+2X46uNLGU6IvjHD9dLOTnLEHOv9xXEdT5+w3e3zbTeLvEiyxWK55zLl880jMvyWna9kuOTNdUn+LcOVU7bMGnveYueHJD+aGRdon7Dt+eabs9b7fbNYqi4baM559dg4d417R5IHT1jv1DZOG3cAAADKWevfAAYAACtGzAIAUJaYBQCgLDELAEBZYhYAgLLELAAAZYlZAADKErMAAJQlZgEAKOv/AJzd9c/cnQtFAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x11635eb70>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArMAAAFuCAYAAACSr5B3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAErRJREFUeJzt3XuMpXddx/HPty29mBZqKYiXwIplRQiKYS0lEgGRNsSAqBiBloJEFC+NqNiaEKXGC4K3KoaLBq2RW6yAyiXYilZRqXSrRfDSpRILqEhLqbRAC9iffzzP0vF0ZndmZ2fnfHder+TJzjnnN8/znDPNr+/zzPOcqTFGAACgo2O2ewcAAOBQiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2jpuI4NPP/30sWvXri3aFQAASK655pqbxhj3Wc/YDcXsrl27snfv3kPbKwAAWIequmG9Y51mAABAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgreO2ewcA4Gj0/Ocn11673XsxefjDk0su2e69gK3hyCwAAG05MgsstYuvvDgXP/bi7d6N9bvzzuTyy5OXvjS56qrk9tuTE09MzjorufDC5Oyzk2N6HUdo9zNYEo6EwpFRY4x1D96zZ8/Yu3fvFu4OwP9XP1MZL1r/PLWt9u1LnvCE5Oabk9tuu/vjJ5+c3PveU+zu3n3k9+8QtfoZAEeFqrpmjLFnPWN7HR4AWFb79iVnnpl8+MOrh2wy3f+hD03j9u07svsHcJQSswCbdeed0+kDn/xkcrDfdo0xjTv77On7ANgU58zCEbJMVzbv5wrnw+Tyy6dTC9Z72tYY0/grrkjOOWdr9w3gKOfILMBmvfSlya23bux7br01eclLtmZ/AHYQR2bhCHEE9Ch21VVH9vtYld9+wM7kyCzAZt1++6F93x13HN79ANiBHJkF2KwTT0w+85mNf98JJxz+fdnBHAGFncmRWYDNOuusI/t9AHyBI7OwBuffsW4XXphcffXany+7mpNPTi66aOv2CWCHELO0tdWxef31G2uTZOqTM87Ymv1hiZ199vSXvT71qfV9PFfVNP4JT9j6fQM4yolZWMOhRKkjpzvUMcdMnzV75pkH/8MJVck97zmNP8aZXgCbJWZpSzSyVHbvTt7znuko7c03r/65s6eckpx22hSyu3cf+X0EOAo5LAAstRc95kXbvQvrt3t38sEPJpddljzucclJJ01HX086abp92WXT481CttXPANhxaqz3zy8m2bNnz9i7d+8W7g4AADtdVV0zxtiznrGOzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKCtGmOsf3DVjUlu2Lrd2XKnJ7lpu3cC2FHMO8B26D73PGCMcZ/1DNxQzHZXVXvHGHu2ez+AncO8A2yHnTT3OM0AAIC2xCwAAG3ttJj9re3eAWDHMe8A22HHzD076pxZAACOLjvtyCwAAEeRwx6zVXVOVf15VX20qu6oqo9U1R9U1UMWxn1jVV1eVR+rqlur6u+r6jmrrO8X5nEfr6pRVc9eY7u/W1X/UlWfrKrbquq9VXVBVR27ytjnVtW/zvt3XVU9b5Uxl87bW1wuWRi32pj9y0+usa/3qKr3zWO+9yAvKXAQyz7vVNWzquqNVXXDvL5L11jfj1bV1fN2b6+q66vqV6rq3gvjzDuwBJZ57qmqL62qF1fV3qq6papurKp3VtU3rbK+J1XV66pqX1XdWVVXrrHdfz/A3PPKA7xO75jH/Nzar+ahO24L1nlakmuSvDzJjUnun+Qnk1xVVQ8bY9xQVV+b5M+SXJXkuUk+neSpSV5dVSeMMV6xYn0XJLk2yVuTnH+A7Z6U5GVJ/i3JSHJOkl9PckaSH9k/qKqem+RVSV4878Pjk7y8qmphu5n3/8kL9/3Xwu1HrbIvP5TkvCR/ssa+viDT578Bh8dSzzuZ5oP7JLkiyXcd5Hm8Kcn7k9ya5OuT/HSSx1XVnjHGnfM48w4sh2Weex6R5LuT/O687eOT/GCSK6vqyWOMt65Y31OSPHwed+IBtvvtSU5YuO87kvxE1ph7qurpSb7uAOvcvDHGli9JvjrTi/3j8+1fSPLZJCcvjHt3kncv3HfM/O8Z8zqevYHtvj7JrStuH5fkY0l+b2Hc72T6YOF7rLjv0iQfOcTne32Sq9d47IFJPpXkGfPz+d4j8TOwWHbasizzzsr1zV9/JMmlG1jf98/78IiDjDPvWCxLsCzL3JPk1CTHLYw5Lsl1Sf5qte3OX/91kis3sN13ZjrQd+wqj31xko8mefr8fH5uK17zI3XO7Mfnfz8//3t8ks8l+czCuP/JwqkP464jEYe63c+vuP2oTEdHXrMw7veT3DvJozexrSRJVT06yVcl+b01hrwiyRuS/O1mtwUc0LLMO4djfVlc50rmHVgqSzH3jDFuGWMszkWfz3Tk98sPx3ar6v5JHpfktWOM/11lyEuSvH+M8fpDWf96bVnMVtWxVXV8VT0o06/1P5rpXUMyHfVMkt+oqi+rqlPnX/8/PsmvbWKbVVXHzev7ziTPSvKrK4Y8dP73/Qvf+k/zvw9ZuP++VXVTVX1+Po/kosVz4VbxrEzvwO72g6uqc5PsSXLRep4PsDFLOu8cyjqPq6ovqqqzkvxMkneOMd57gG8x78A26jL3VNXxmQ7s/cuhbnfBM5NUVnkjPb/JPj/TKVBbaivOmd3v7zKdr5FMv/765jHGx5JkjPH+qnpskjdnOn8jmd61PG+M8YZNbPNbk7xl/nok+cUxxs+uePy0+d9PLHzfzQuPJ9M7l2syhe6Jmc4TeXGSByVZ9eKJqjox0/lwbxtjfHzhsS/O9B/ZRWOMm6rq5A08L2B9lnHe2ZB5brh1xV1/mgOcZ2vegaXQZe65OMlXJDl3E9td6fwk/zDGeN/KO+doflWSXx5jXHeYtrWmrYzZZya5Z6ZztV6Q5IqqevQY49/ndy5vzBSKz8t06P3bkryyqm4fY7z2ELf5riTfkORemd7xvGC6rmu8cKMrGmNcsnDX26vqtiTPr6qXjDE+sMq3PWXe9qWrPPZLmU7UfvVG9wVYt9bzzuzT8/pOzHQB2AuTvKWqvmXxV4Yz8w5sv6Wfe6rqGZkuTvvZMca7DnGbK9d3VpLd+f8Xu+53YaaL1H5+s9tZl604EXeVE4BPTXJLklfOty/LNMHeY2HcazNdiHXMKus4lJOhL07yv0m+fL79A/M6vnRh3H3n+3/oIOs7cx739DUef3umC8wWn9cj5/14zPxanJrka+d1XTDfriPxs7BYdsqyLPPOKo9v9AKwx8778LQ1HjfvWCxLtCzj3JPkSZlORfrtdaxnXReAZfoEh88mOX3h/vtnCvZzV8w9p87P55fmr+92sdhmliNyAdgY45ZMh93PmO96WJL3jjE+tzD0PZkuxLrvYdr03kznBX/lfHv/ubEPXRi3/1zZf17neu/2Z9Oq6n5Jzk7yulWe19fM+3FlplMcPpFk//lvvzHfvtc6tw2swxLNO4djfcldz+MLzDuwfJZt7qmqx2cK6jdn+nSUTauqE5I8Lcnbxxg3LTz8wEy/WXpN7pp79p/e+YL564cdjv3YbytPM/iCqvqSJA/O9C4kmU6MfnhVHT/G+OyKoY9McnvuOod1sx6TKTw/ON9+d6Z3Qedm+sy3/c6bt/k3B1nfufP6rl7lsfOSHJvVryZ+R6ar/Va6X6aTw385yduS3HaQbQMbsETzzuFYXzId2Vlk3oEls0xzT1U9KskfZ/r4rPPG5j4tYaUnZfrYrdXmnmtz97knSf4iU+C+OlPsHzaHPWar6s1J/j7JPyb5ZKbzKX4008dF/Mo87DczvUt4S1W9PNPh6Cdn+hyyX1v5w66qx2T6OK37zXftmc9dzRjjD+cx35rkezKdCP2hJKckeWKS70vyqjHGf87jP1dVP5XpjyT8R6ag/eYkz0lywf7tVtUDMn1c1xsyveAnZLoA7Nnz+lb7n8r5Sd43xviHxQfGGB/N9B/zytdp1/zldWOMK1d7LYH1WeZ5Zx77kNz1G6CTkjygqp463/7LMcaNVXWvTAH62iQfyPQ/pTOT/FimI6pvWuWpm3dgGy3z3FNVD870pvWmTL/ef0RVfWHfxxhXrdjuAzKdf5tMR4vvXDFHXT3GuGHhqZ+f6aPA3rb4msxHpq9c5bVKkhu2ZO7ZgnNFLsr0KQC3ZLqQ4bpMV7TtWhj3xPnJ3pjpyt1rM13ld+zCuCszTep3W1aMeXCmif7DSe5I8t+Zzvk4N6ufi/L9SfbNYz+Q5AcXHj8tyR8luSHTu6ZPZ/qP9YfXWN/XZ8UHJK/zddoVH15usRyWZdnnnUznsq26viSPnceckOkPuOzL9AcO/idTxL4wySmrPGfzjsWyzcsyzz2ZDsCtNe+Mhe0eaOyzF8beJ9OnMbxsg6/Vlv3RhJo3AAAA7RypvwAGAACHnZgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANDW/wHeuBCTXiBCigAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x115f6ff98>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArMAAAFuCAYAAACSr5B3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAEbFJREFUeJzt3XuMpXddx/HPt5S2wJYUWsRLLJsIi0jBIptSgtKi7VYIigJBiSYiFxUDiEHaVE0oESFUCYiIcg2C4KVGEm9IK1Cj0lq2WC6KLIKlxaSyUNould7ozz+epzg9zHRnpju78915vZIns3vmN+d5nnMmv33vM79zpsYYAQCAjo441AcAAADrJWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaOvItQw+4YQTxvbt2zfoUAAAILn88su/NMZ4wGrGrilmt2/fnt27d6/vqAAAYBWq6vOrHWuZAQAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAto481AcAa/XiFydXXHGoj+LOTj45ed3rDvVRAMDWI2YBYBX8Rxo2JzFLOybureW8i8/Leaefd6gPY/9uvz258MLk/POTSy9NbropOeaY5NRTk7PPTnbtSo7oubKrzXMAbEk1xlj14J07d47du3dv4OEA3Fm9vDJetvp56pDYsyc588zk2muTr371mz+/bVty/PFT7O7YcfCP725q8RwAh5WqunyMsXM1Y3teJgDYLPbsSU45Jbn66uVDNpluv+qqadyePQf3+AAOc5YZQKyFY51uv31aPnDDDcn+fso1xjRu167kc59ru+QAYLMxmwKs14UXTksLVrtca4xp/EUXbexxAWwhrsxCXAFlnc4/P9m3b21fs29f8upXJ2edtTHHBLDFuDILsF6XXnpwvw6AbyJmAdbrppvW93U333xgjwNgCxOzAOt1zDHr+7qjjz6wxwGwhVkzC7Bep56afOhD6/s6Nj3vcgI9uDILsF5nnz39QoS12LYtOeecjTkegC3IlVmA9dq1a/rNXjfeuLq356qaxp955sYfG3ebK6DQg5iFDbbZflTpx5QH0BFHTO81e8op+//FCVXJfe87jfcLEwAOGDMqwN2xY0dy2WXJiScmxx67/Jhjj50+f9ll03gADhhXZmGDuQp697zstJcd6kPYvx07pl9Re9FF0y9EuPTS6e23jj56erHXOedMSwuaXpFt8RwAW1aN1f4axiQ7d+4cu3fv3sDDAQBgq6uqy8cYO1cztudlAgAAiJgFAKAxMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoK0aY6x+cNXeJJ/fuMPZMCck+dKhPghgSzHvAAfb4TTvPGiM8YDVDFxTzHZVVbvHGDsP9XEAW4d5BzjYtuq8Y5kBAABtiVkAANraKjH75kN9AMCWY94BDrYtOe9siTWzAAAcnrbKlVkAAA5DByxmq+qsqvpgVV1TVTdX1Req6s+q6nsWxj2uqi6sqi9W1b6q+mhVPXuZ+3vlPO7LVTWq6lkr7PfeVfXyqtpTVV+rqqur6p1VtX2ZsT9TVZdX1Q1VtbeqLqqqH1jneZw+H9fidt3CuGOr6rer6uJ5v6OqTl/t4wqsrMm88/yq+o/5+K6qqt+oqnsujDmpqt40z0+3VNWyPzKrqu0rzDujqo5bz7kAa3O4zDvzuKdX1b9W1U3z+byhqo5dGLPpe+dAXpm9f5LLk7wgya4k5yZ5eJJLq+pBSVJVj0zy90numeR5SZ6a5CNJ3lZVz1+4vxcmuVeSv97Pft+a5KVJ3pLkSUl+Pcnjk3ygqrbdMaiqfi7JO5JcluRpSZ6b5KgkF1XVo9ZyHgtelOSxS7YzFj5/fJJnJ7ktyUX7ORdgbTb7vHNukt9L8jdJnpzkd5O8JMnvL9zfo+f7uSrJ7lWc96ty53nnsUn2rfNcgLU5LOadqnpmkguSfCzJU5Kcl+SZSf5ihf1v3t4ZY2zYluShSUaSl8x/f2WSW5JsWxh3SZJLFm47Yv744Pk+nrXM/d97ftBeuXD7D89fc9aS2z6c5MML446dj+dVazmP+bbT59vO2M/X1pI/nzF/zekb+bjbbFt52yzzTpJjMgXmOxbG/UqS25M8fHG/859fMU3Ny57b9nkfz13F47Dfc7HZbAdmazrv/GeSixfGPX2+vyctuW3T985Gr5n98vzxtvnjUUluTfK1hXHXZ+Eq8Rjj9lXc/z3m7YaF2++49L30Po9aZtz/zsezv8dh8TxWbczPKnDQbJZ556Qk25K8b2Hc3yWpJD+2xv2uyUbcJ7CiVvNOVZ2Q5LtWGJckP76KY7qTQ9k7Bzxmq+oeVXVUVT0kyZuSXJPkj+dPv2P++Pqq+vaqOq6qnpfkh5K8dq37GmPsS/KuJC+qqidU1baqeniS38p02fwDS4a/MckZVfWceb/fkeQNmb7Z3rbG81jq3VX19Xmty3uq6sS1ngdw92zSeefr88dbFu7i5vnjSWvd9xKvqqrbqur6qvrLqnrE3bgvYB2azzsrjbs10xXV5eanTds7R27Aff5LpvVfyXQJ+wfHGF9MkjHGJ+fFwO9N8ovzmFuT/MIY40/Wub+fTfL6JB9cOIYzxxjfeJLGGG+vqmSK2rfON18zj9uzlvOYXZ/kNUn+IdP/lB6V5FeTXFJVj1oYC2yszTjvfCbTj/VOnfd9h8fOH++/jv3enOkfzQuT7E3y3ZnmnQ9X1SljjE+t4z6B9Wk774wxvlJVe+dxSz0m0xXcpfPT5u+dDVg38rBMD8YzMy2Q/kKS7fPnHpLpBQ7vz7Qo+YcyPTG3JvmpFe7vLtd7ZXohxPWZFjc/PslPZ3oydye5z5JxT8l0uf935v0+OdMC5b1JTlrLedzFuX9fph8xvGKFz1sza7NtwLaJ5523JLkxyU8mOS7JE+Zjuy3J+1a47xXXzK4w/jsz/QPzR+s5F5vNtr6t+7yT5NcyXaF9QaZ4fXSSf5vHfWo/576pemejn+jjMq3n+IP57xck+WySey6Me3eSL2XJiyBW8+RmevXgSPKchdsfMt/+S/PfK8n/JHnXwrijknwuyXvXch77GfvvSd6/GZ5cm20rbptl3plvu2+mVwbfPn/ua5lejbw3ydtXOP41xez8NX+b5NMrfE7M2mwbvHWcd+YGetMcpSNTaL8uUxx/cBXnvGl6Z0NfADbGuC7TpfcHzzc9IsnHxhi3Lgy9LNNbOnzLGndxxzqxjyzs9zOZvqkeNt/0wPm+F8fdkmmtycNyF5Y5j/3xoi84RDbRvJMxxg1jjKdmmoMeOe/rnUlOSPJPa9zv/ph34BDpOO+MMW4ZY/z8fPv3zuNfmimQVzs/bYp5Z0NjtqoemGlN12fnm65JcnJVHbUw9DFJbkpy7Rp3cc388ZSF/e7I9L+k/55v+kqmtWaL445KcvKSccta5jxWGrcz09tzXLa6wwcOtE0073zDGGPvGOMTY3oRxy9nujJzwRr3u6z5RRjfH/MOHDKd550xxnVjjI+PMa5N8pwkRyd5+10dzGbrnQP2ArCqem+Sjyb5eKb1WzsyPXi3ZVo4nEzvHnBBkr+qqjdmuvT9o5nWm7x2LHnBVlWdluQBSb51vmlnVX01ScYYfz7f9o+Zrqy+pqrul+nS+ImZ3kj4+iR/OI+/uarekuQFVXVtpjcmvlemdSLbk7x4jeeRqnp3kv+ax16XaUH0uZm+oV6/8Ng8Mcl98v//szptfluMG8cYi2+LAazSZp535vv7iUxr0T6d5H6Z3jj9GUmeNv8Dc8e4e2d6E/Rk+gcxVfX0+e9XjjF2z7e9JtNFiEsy/cjwoZnmnduT/ObCY7OacwHW6DCad87M9K4Fn8z0/rS7Mr1Y7YVjjCuXjNv8vXMA14uck2kB9HWZ3r/105nWYmxfGPfEJBdnmoj3Jbki04N3j4VxF2e6fP1N28K44zN983wm0zfL1Un+NMlDF8YdmSler5j3u3fex651nse5mb6Rr8+0zuTqJG9O8m3LPDZXrnAuVx6MtSQ22+G6NZh3npHkE/Ox3ZDpXQget8x5bF9pv1ny5ueZfrvORzL9tOnWTFdr3rO437Wci81mW9t2GM07p83zyb5MLxj75yQ/ssy4Td87Ne8cAADa2ejfAAYAABtGzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaOv/ABzBKZUy64TLAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x115f7f0b8>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArMAAAFuCAYAAACSr5B3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAE4hJREFUeJzt3XuMpXd93/HP18GAjdf1hVu4OMslTrAEGHVrLrmADXaaNBGhStIGUGoFEI2rqLRcjCgNS5IqtpUmFX9wCRDlouaCQVBSimPngtIoBXeNAjgENsHF2AkQwDXGgF2SffrH75n45PjM7OyuZ/d8d18v6dHoPOf3PPOcObO/ec8zzzlb0zQFAAA6OulYHwAAABwuMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW/c7lMEPfvCDp927d+/QoQAAQHLDDTd8cZqmh2xn7CHF7O7du7Nv377DOyoAANiGqrp5u2NdZgAAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCywlvZ+YO+xPoTDc+BAcs01yUUXJaeempx00vh40UVj/YEDx/oID1vb5wQ4rtU0TdsevGfPnmnfvn07eDgAQ72+Mr1u+/PTWti/P7n44uS225I777z3/aedlpx9dnLttcm55x794ztCLZ8ToKWqumGapj3bGevMLMB9Yf/+5IILkltuWR2yyVj/mc+Mcfv3H93jAzhOiVmAI3XgQHLJJckddyQH+2vXNI1xl1zS+pIDgHUhZgGO1LXXjksLtnvZ1jSN8dddt7PHBXACELMAR+qqq5KvfOXQtvnKV5Irr9yZ4wE4gYhZgCP1wQ8e3e0A+HtiFuBI3XXX4W1399337XEAnIDELMCReuADD2+7Bzzgvj0OgBOQmAU4Uk972tHdDoC/J2YBjtSrXjX+Q4RDcdppyeWX78zxAJxAxCzAkbrkkvE/e1Vtb3zVGH/xxTt7XAAnADELcKROOmm81+zppx88aKvGuGuvHdsBcETMpAD3hXPPTa6/PjnnnGTXrtVjdu0a919//RgPwBETs8Baet0zX3esD+HQnXtuctNNydVXJxdemJxyyjj7esop4/bVV4/7m4Zsy+cEOO7VtN3/fjHJnj17pn379u3g4QAAcKKrqhumadqznbHOzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANoSswAAtCVmAQBoS8wCANCWmAUAoC0xCwBAW2IWAIC2xCwAAG2JWQAA2hKzAAC0JWYBAGhLzAIA0JaYBQCgLTELAEBbYhYAgLbELAAAbYlZAADaErMAALQlZgEAaEvMAgDQlpgFAKCtmqZp+4OrvpDk5p07nKPmwUm+eKwPAjghmG+Ao+V4mm++ZZqmh2xn4CHF7PGiqvZN07TnWB8HcPwz3wBHy4k637jMAACAtsQsAABtnagx+0vH+gCAE4b5BjhaTsj55oS8ZhYAgOPDiXpmFgCA48COxGxVfU9V/UFVfa6q7q6qW6vqHVV13sKYD1TVtMlyzdL+zq+qa6rqzqq6o6reW1WPX/F5z6mqX62qz1TV16tqf1X9bFU9aItjfUZVHZg/7/2W7vumqvp3VXVjVX21qj5bVe+uqidtsb+Tq+pj8/5evHTfD1XVu6rq5vn4PllVP1dVu7bzdQXubd3nmy0+98tW7POUqtpbVX8xP5bPV9V/r6r7b/LYz5jnpamqnrN030uq6n9U1V/N89eNVfXKzfYFbM+6zznz2DOr6r/MYzeO8Ve2eEyPraqvzcd3r8+9MG4t55z7HXzIYTkryQ1J3pjkC0nOSfLqJB+sqidO03RzksuSnL603dOT/EKS926sqKpvTfI/k9yY5AXzMb8uyR9V1fnTNP3NPO5BSX4vyclJ/mOSzyT5J0len+Rbk/yL5YOsqpOTvCXJ55M8fMXj+Jkklyf5uSR/kPH+bf8hyR9W1ZOnabp1xTavmMet8or5uF6T5NYkT0myN8mFVfWMaZoObLIdsLkO881Hk7x0ad2nF2/M89H7kzwmY875eJKHJLk4yTdt8tiv3GR9kvxUkuuS/HKSLyX5zow57YIkP7zFdsDW1nrOqaozk/xxkinJazPmmkck+Y4tHtMbk3w5ySkHeezrOedM03RUliTflvGFffkWY96e5O4kZy2se1uS25OcsbDuUUnuSnLVwrpL5v1fsrTPK5L8bZJTV3y+12R8A/2nedv7Ld3/10l+Y2ndt89jX7pif49N8tUkz5/HvHjp/oes2ObH5rEXHa3nwmI53pd1mm+SfCDJH2/jmF+d5I4kj97mY/yOeb758flYnrN0/6r55qfmsY891s+RxXI8LWs257w54z+4On2bx/78jJN6L5s/x+M3Gbe2c87RvGb2S/PHv111Z1WdmlHuvzNN020Ldz0tyf+apun2jRXTOCN6Y5LnLYzbOI19x9Kub8+4nKKWPt/jMn5juSzJNzY55vtvsr9k9SUab0ryW0n+ZNXOpmn6worV/3v++MhNjgE4dGs132zTZUmunqbploMNXPir0hVJblo1xnwDR9VazDnzGdwfS/K2aZqWx646rjMzzha/Ivf0zapxaz3n7GjM1rjm9P7zafS3JPlckt/cZPjzkuxK8qtL6/8uyf9bMf7uJI+rqgfOt38vyV8kubKqzquq06rqoiT/Nsmbp2n66tL2b874wfFHWzyENyZ5YVU9t6pOr6rHzutuTfKOpcf6giR7Mi5LOBTPnD/++SFuByxY8/nmKVX15ar6RlV9tKpetHTs5yR5dJKbquqt83Vzd1XV71fV+SuO51UZP9yu2uTxbeaZSQ4k2X+I2wFL1nTO+ccZlwp8vqreOV9be2dVvaeqHrPi81yV5BPTNP36QR7ues85O3zafV/G6eUp40l4whZjfzfjNPfyn/rfkRGPJy+s25XxG8SU5JsX1j8049qTaWF5a5KTlvb5wiS3JXnofHtvVlxmMN/32oxvto39fTLJ45bGnDkf+4vn27uz4jKDFft+ZJK/SXLdTj4PFsuJsKzxfPPTSV6SMak/N8m75rGvXRjztHndHUl+P8n3Zfzw++j8uc9ZGPv4JF/P/Ce+JM/Kij/5rXjMT5q3e+uxfq4sluNhWcc5J8m/XJhLfjvjmvvnZ1x2cHOSXQtjvysjms+bb1+aFZcZdJhzdvqJfkKSpyb50YyLpW9NsnvFuEdkBOMvrLjvO+cv2tsz4u9bkrwz41T+lORh87gHJvnDjNh8YZLvzjhtfkeSNy3s76yMgPzXC+v2ZvU1sz8xPxGvn5+8H5q/eT+V5BEL496WcWnBxvv27s5BYjbJafO+/jrJo47mP0CL5Xhc1nG+2eJY3z3PLafNt58x7/9z+YfXvj16/mFz5cK665L81sLtg/5gSfLNSf5Pxp8udx3s+CwWy8GXdZxzcs9rdv50o0nm9U+d1//EfPv+GS8yvWJhzKVZHbNrP+cczSf9jIzfNN684r5XzV+YJ2+y7WUZr7Lb+E3kutxzIfXJ85h/M9+3fNb0JYv7zrhM4CMZUXvGvFwxjzk7yYPmcWdlXID9M0v7O2v+5vnFhW+Qv8s467KxvyfN+/vJ+XYt7eOU+ZvytiRPPBb/CC2W43lZl/lmi+P7kXnc0+fbGy8eedeKsR9Jcu3Cdl9Nct7CfPP987bPTfKPVmx/9vwD5dNJHnmsnxuL5Xhc1mXOSfI98+3/vOLzfDnJLy0c02czQntjLrls3vYpmQO0y5yzU2/NdS/TNN1eVX+Zcbp62b9K8pFpmj6yybZvrKq3z9veMU3TLVX1/iQfmqZp48VbT0zyf6dp+tTS5tfPH5+Q8UPhvIzY/FLu7YtJ/luSH0xybpIH5J6LlzeO5baq+tS8v439npTxiuVlb5iXMzNfWD1fRP3OjOtrL56m6WOrHjNw+NZovjnooc4fb8o4U7uZjbftOy/JqUn+bMWY92T8sDpjY0VVnZ7x582zk3zXNE1/tY1jAg7RGs05q+aGRYtzycOTrJoTPjzv6/w0mXOOWsxW1cMy3tbqvy6t35Pxxfr3W20/TdPdmb+YVfXEJM/JeMXehs8lObOqHj9N018urH/q/HHjC/qyLHzhZ5dmfLM9J+Oalo39JeP90RbfE+6sjG+4D8+rrkly4dL+Hp5xEfjPJ3lfkjvnbU/KePwXJfn+aZo+uNVjBg7PGs03m3lBRrx+bP5836iq9yX57qp60DS/mGN+Ydi355456Fdy71+cz0/yixl/cvzQwmM9NWP+eUySZy0dJ3AfWpc5Z5qmW6tqX5KLq6qm+VRpVT09431vN07QXZExnyz6pxkvYn9hxuUMSZc5Z4dOt7874019n5sRei9N8omMs5PnLo19Q8ZbYz10k309KuN9YP9ZxpP76ozfBH57adzujD//788I0wuTvHJety9LL8pY2nZvVl8z+zsZp/l/OsmzM06378t45eGeLfa3O6vfZ/ZN8/qfzXjBx+LiulmL5TCWdZ5vMl5g8b4kL5rnkH+e8defKcnlS/s8L+MX3w8k+YGMt/G5MeMX7Idt8fifldXv+fj+jLMwP7livrnX+0FaLJbtLes858xjn51xze27knxvRhTfkvGuSads8bguzRbvM7swbu3mnJ16oi/PuBj69iRfyyj8t2TpwuiM/8niCxnvu7bZvh6W8ZYUX8wIy48neXlWv/PAeRmvDLwl46zH/oyzo2ce5Hj3ZnXMnjp/w34845qRz2b8YLrgIPvbndUx++n8w1chLi57d/KJtliO12Wd55uMv+K8P+Osyd0ZsfonSX50k89/Qcb19F/L+IH2niP4wbLZXDMlufRYP28WS9dlneechbHfm3EW9q6Myyp/LVv8Ujxvc2mOLGaP2Zyz8ep7AABo52j+D2AAAHCfErMAALQlZgEAaEvMAgDQlpgFAKAtMQsAQFtiFgCAtsQsAABtiVkAANr6/+rptFOxwZybAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1161012b0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "for pos, tstart, tend, matemaps in calls:\n",
    "    plot_variant(pos, tstart, tend, matemaps)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}