{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f403fd5de10>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEBCAYAAABxK3LCAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XlclNX+wPHPYRMVREPBBRICF1CRRfMq4prZ1cqwRe2n\nV1OzNC2915ZbXbPFcst2LStNvWJ2Lb1mWrlA7gsguIBroo4iIrLvzJzfH8BckFWdhYHzfr2eF/M8\nc+Y5Z84M5zvPeZ7nHCGlRFEURWl4rMxdAEVRFMU8VABQFEVpoFQAUBRFaaBUAFAURWmgVABQFEVp\noFQAUBRFaaBUAFAURWmgVABQFEVpoFQAUBRFaaBszF2A6rRs2VJ6eHiYPN/s7GyaNm1q8nwtiaqj\n6qn6qZmqo+rdTf1ERUXdkFK2qildnQ4AHh4eREZGmjzfiIgIBgwYYPJ8LYmqo+qp+qlZdXW0evVq\nOnToQO/evU1bqDrkbr5DQoiLtUlXpwOAoigN0/jx4wFQY5UZlzoHoChKnbV69WpzF6FeUwFAUZQ6\na/z48Vy7ds3cxai3LK4LqLCwEI1GQ15entHycHJyIj4+3mj7rw+cnJy4cOECbm5u2Nramrs4Sj3W\npk0bCgsLsbGxuOaqzrO4GtVoNDg6OuLh4YEQwih5ZGZm4ujoaJR91xcZGRkUFBSg0Wjw9PQ0d3GU\nes7W1ladDzACi+sCysvLw9nZ2WiNv1I7QgicnZ2NeiSmNExVNfTqf97wLC4AgPoi1BXqc1CM4fLl\ny/rHWq3WjCWp/ywyACiKUn+VnlMaMmQIVlZWpKSk6J87fPiwuYpVL6kAoChKnTRy5EgA7rnnHv22\nXr16mas49ZIKAHfAwcHB3EWoVGm5EhISCAsLM3NpFOXOlHb7WFtb67ft2LFD/1h1PRqOUQKAEMJd\nCBEuhIgTQpwUQrxUsv0eIcR2IcTZkr8tjJF/Q6cCgGLJKgsAgwcPLpdGCMH58+fJzs42adnqG2Nd\nBloE/ENKGS2EcASihBDbgQnATinlfCHEa8BrwKt3msnMmTOJiYkxSIFL+fv78+6779YqrZSSV155\nhW3btiGE4M0332TUqFFVbo+IiGDOnDk4Ojpy7tw5Bg4cyNKlS7GyqhiHv/zyS86fP8+iRYsA+O67\n74iMjOTzzz9nyZIlrFixAoDJkyczc+bMcq997bXXiI+Px9/fn/HjxxMaGsq4ceP0/yyff/45ffr0\nQafTMX36dHbt2oW7uzu2trZMnDiRJ554gqioKP7+97+TlZVFy5Yt+e6772jTps3dVK2i1EpSUhIA\nJ06cKLc9Pj4eHx8f/bq3t3e55y9cuIA5Bo+0ZEY5ApBSJkopo0seZwLxQDtgBLCqJNkq4DFj5G8q\nP/30EzExMcTGxrJjxw5efvllEhMTq9wOxSexPvvsM+Li4jh//jw//fRTpft+/PHH2bhxo359/fr1\njB49mqioKFauXMmhQ4c4ePAgX3/9NUePHi332vnz5xMSEkJMTAyzZs3CxcWF7du3Ex0dzfr163nx\nxRf15U9ISCAuLo41a9Zw4MABoPhmuxkzZrBhwwaioqKYOHEib7zxhjGqUFEqKP3ef/HFF+W2d+7c\nmdmzZ1f5Ok9PT9U9dJuMfiOYEMIDCAAOAa5SysSSp64BrpWknwJMAXB1dSUiIqLc805OTmRmZgLU\n+pf67dJqtfo8qpKZmcmuXbsIDQ0lJyeHJk2a0KdPH3bv3s3u3bsr3e7o6EhQUBCtWrUiJyeH0NBQ\ndu3axdChQyvs397ennvvvZedO3fi5eVFfHw8fn5+LFu2jGHDhqHT6QAYPnw427dv1/8ayszMJCcn\nh6KiIv17SE9PZ/bs2Rw/fhxra2vOnTunL/8jjzyiH3Y2JCSE3NxcoqOjOXHihP6wW6vV4urqWq5O\nSusoLy+vwmekQFZWlqqXGlRVR2fPngWgWbNmFZ4fPnw4/v7+jB07tsr9CiH44YcfaNWqxtGQ6zRT\nfIeMGgCEEA7Aj8BMKWVG2egspZRCiAp3fEgplwPLAXr06CFvHQ41Pj7e6Hfp1uZOYEdHR+zs7LC3\nt9entbW1pXHjxlVub9KkCTY2Nvrt9vb2NGrUqMq8xo4dyy+//ELnzp0ZOXIkzZo1q/CaRo0alcvL\n0dGxQj4ffvghbm5uhIWFodPp9OlvLaeNjY2+nF26dNEfEVRXR/b29gQEBNS2ahsMNRx0zaqqo5s3\nb/Ljjz+ydu3aKuvw//7v/wDIz8/H1taWa9eu0a5dO/3zTz31FCdPnsTX19cYRTcJU3yHjHYVkBDC\nluLGf62UsrSfI0kI0abk+TbAdWPlbwohISGsX78erVZLcnIyu3fv5v77769yOxR3AV24cAGdTsf6\n9evp27dvlfsPDQ3lv//9L+vWrWP06NH6PDdt2kROTg7Z2dls3LiRkJCQcq9zdHQs92s9PT2dNm3a\nYGVlxZo1a/Qn2YKDg/nxxx/R6XQkJSXpf2106tSJ5OTkcl1CJ0+eNFi9KUp1Su8ur01/fqNGjbCy\nsqJt27ZIKTlz5oz+uS5dupCVlWWsYtYLxroKSADfAvFSyiVlntoMjC95PB74rzHyN5XQ0FD8/Pzo\n3r07gwYNYuHChbRu3brK7QA9e/Zk+vTp+Pj44OnpSWhoaJX7b9GiBT4+Ply8eFEfQAIDA5kwYQL3\n338/vXr1YvLkyRV+gfv5+WFtbU337t356KOPmDZtGqtWraJ79+6cOnVKP8vQ448/jpubG76+vowd\nO5bAwECcnJyws7Njw4YNvPrqq3Tv3h1/f3/2799vpFpUlPL27dsHcEezYXXo0IGcnBz9uhrTqwZS\nSoMvQF9AAseAmJJlGOAM7ATOAjuAe6rbT1BQkLxVXFxchW2GlpGRYZT9hoeHy+HDhxtl33cqMzNT\nSinljRs35H333ScTExNr9brSOjLF52GJwsPDzV2EOq+qOippO+TNmzfveN9FRUX6/RQ3c5bnbr5D\nQKSsRVttlHMAUsq9QFWn4wdXsV0xg4cffpi0tDQKCgr417/+pT9SURRza9KkyR2/1tramqtXr9K2\nbVug+MSwVKOJVmBxw0FbsgEDBlR6UqdXr17k5+eX27ZmzRq6detm9DKpK1WUuqpRo0Z39fo2bdoQ\nGRlJjx49ABUEKmORAUBKWa+u9z106JC5i3BH1D+TUtcFBQUxffp0Pv/8c6D4vFpqaqqZS1V3WNxY\nQPb29qSkpKjGx8yklKSkpGBvb2/uoij1jIuLC4GBgQbb32effaZ/nJaWxquv3vHgA/WOxR0BuLm5\nodFoSE5ONloeeXl5qmGrQV5eHs2bN8fNzc3cRVHqmevXr3P9umGvEC/ba7Bw4UKuX7/OypUrDZqH\nJbK4AGBra2v0KQgjIiLUzU01UHWkGNMjjzxi8H2WDQLfffcdrVq1YuHChQbPx5JYXBeQoij1m4OD\nQ4WB3gylbNfxokWLmDhxolHysRQqACiKUmfk5+eTlZV1RzeB1VbZILBy5cp6dUHJ7VIBwES++uor\nhBCVLleuXDF38RSlTpg6dSqA0f8nbr2IpKEGARUAjESn09G/f399I//8889XmdbNzU2fbu/evRXu\nCVCU+i4zM5Nx48axfft2gHLDORiLlJL+/fvr14UQ1Q6AWB+pAGBgy5cvRwiBtbU1u3fvrjTN+fPn\nyc/P1w/wVlZISAj29vZVHi3cupw8eZLCwkJ1WaxiEAMGDEAIQVFRkUnzbdasGf/+97/RaDQABr8K\nqCoREREcPHhQv96nT58GdTSgAoAB6HQ6fYP83HPPVXh+zZo16HQ6/fgb9913H3Z2dqxbtw4pZbmJ\nX25X165dsbOzw8rKqlxgiIqKupu3pDRA33//PX/88QdQfLVd6Xfpr3/9qz7N1atXuXbtmn5dq9Xy\nyCOP8MADD+hHma1ObY9uw8PDb7P0d65Xr14VppYUQnDhwgWTlcFsajNgkLmWygaDM4XaDsKUmZlZ\nbsCpssvbb799V2UoKiqSly9flsuWLZOTJk2Svr6+VeZ1u8uNGzfuqmxSmn6wM51OJ3U6XbltZQcL\nS09Pl7Gxsfp0+fn5Mjk5Wfbs2VPGx8fL8PBwGRUVJbOzs/X1sGvXLnn27FlZVFQkDxw4II8cOSJv\n3rwpAXnffffJ9PR0mZ+fL2fMmKF/TVhYmATktm3bytXpggULpKurq8zIyJDHjh2Ty5Ytk6tXr5aA\nnDRpkpSyeAC9K1euyPnz58tVq1bJvLw8WVBQID/66CMJyBkzZsizZ8/KAwcOyC+//FK+++678ty5\nc/o8bG1tJSCPHj0qU1NTZUFBwW3VYXJycrkyHz9+XGq12lp9ZyZNmlRuPS4urlxdli6DBg2qMADb\niBEj5Jo1a8ql0+l0+sezZ8+usJ/AwMC7+LbcOQcHhwpluXz5slnKYorB4MzeyFe31OUA0KRJk0r/\nUVJSUoxfwDK0Wq2Mi4uTTzzxhMEChFosZ7G2tq4xzdy5c+XVq1drvc/ly5dLQL744otVpgkICDDq\n+9q4caNJ/4/Kys/Pr7RMV69eNWk5VACogwFg06ZNlX45bv11am7p6enlfj2qRS21Xcq69bnQ0NAK\n21JTUw2S70svvaR/vGXLFjP95/xP6RHcrcvmzZtNkr8pAoA6B1BLUhbfRfjYY+XnsS/t269rJ46a\nNWuGl5dXpR966WBYwcHBzJkzh+bNm5usXI0bN67yuXvuuUf/eNSoUfz000+8/PLLhIWF8fHHH7Nr\n165y6VNSUvj222/169euXePpp5/m6aefJjc3V//ZZGVl8fPPP3Po0CHCw8PJy8tDSql/fseOHVy5\ncgWdTkd6ejoFBQVIKSkoKCAtLY0rV67o6y4nJwcpJT/88AMLFy6sULfZ2dls3ryZ3377jdTUVKSU\naLVavvnmG9LT09FqteTl5ZGYmEhmZiYbNmzggw8+QEpJly5dCAgI4Pvvv0dKqb9L9ZNPPqn811sV\nSicPqkx+fn6F1/7www/861//qjR9SkoKfn5+DBs2jEWLFvHTTz+Ve/6VV16p8P2p7FzAkSNH+PPP\nP/X1Xqrs44CAAP2FEdV9T0xl3LhxSCkJCwsrt/3RRx/Vnx8pnbzGYtUmSphrqStHAPv27avwKyAv\nL88sZasr1IQn1TNV/RQUFMjCwkIppZQfffSRPHToUIU0qamp8qOPPpJ//vlnue0ffvih/Pjjj6WU\nUs6aNUsCcvTo0bXKNysrq9x6YmKi/PXXX2+r7KV1VPo/tWLFCnnz5k25aNGiOndELaWUR48erfYI\nZvTo0bKoqMhg+akuIDMGgKtXr8rCwsIKfZ2vvvqqWcpU16gAUD1Lq5+goCAJyCeeeMJkeZbWUXp6\nupw1a1adbPSrMmXKlFp1a4WFhcn8/Pw7ykN1AZlJRkYGbdu2xdbWlqNHj+q35+fnM3/+fDOWTFGM\n48SJEwCcPHnS5Hk3a9aMJUuW1Llu1Op89dVX+kb0yy+/rDLd008/TaNGjSrcv9O3b18OHz5Mbm6u\nCUtdkckDgBDiISHEaSHEOSHEa6bOvzYyMzMrbJNSYmdnZ4bSKIrxTZ48GYApU6aYuSSW57nnntMH\nA51Ox+rVq2t8zb59++jVqxdNmjSp8ibPmTNnGr3sJg0AQghr4Avgr4AvMEYI4WvKMtSGvOUk2a3r\nilLfjBkzBoDBg9WU3XdDCKE/eVx2ycjI4OWXX76tfcXGxhqplP9j6vkA7gfOSSn/BBBCfA+MAOIM\nmYlWq+Xy5ct3/PrS29FBNf5KwxAcHKy+60bk6OjIwoULq51/4Pz58yxYsICvv/4awCTdzaYOAO2A\nsi2zBuhl6Exu3rxp9EljFEVRDMnLy4vly5ezfPlyoHicImOrczOCCSGmAKUdkVlCiNNmKEZL4EZJ\necyQvUXQ15FSKVU/NVN1VL27qZ/2tUlk6gBwBXAvs+5Wsk1PSrkcWG7KQt1KCBEppexhzjLUdaqO\nqqfqp2aqjqpnivox9VVAR4AOQghPIYQdMBrYbOIyKIqiKJj4CEBKWSSEmA78BlgDK6SUpr/wWFEU\nRTH9fQBSyq1Syo5SSi8p5TxT519LZu2CshD1so6EEM2FENNKHrcVQmy4w13Vy/oxMFVH1TN6/Qh1\n6Zei/I8QwgPYIqXsauaiKIrR1bmrgBTFzOYDXkKIGOAs4COl7CqEmAA8BjQFOgCLATtgHJAPDJNS\n3hRCeFF8s2MrIAd4Vkp5yvRvQ1FqpsYCUpTyXgPOSyn9gVtv3ewKjAR6AvOAHCllAHAA+FtJmuXA\nDCllEDAbWGqSUivKHajTXUAtW7aUHh4eJs83Ozubpk2bmjxfS6LqqHqqfmqm6qh6d1M/UVFRN6SU\nrWpKV6e7gDw8PIiMjDR5vhEREQwYMMDk+VoSVUfVs7T6yc/PJywsjPHjx2NlZZqOAUurI1O7m/oR\nQlysTTrVBaQoCh9++CETJ06sMPuVUr/V6SMARVGql5+fT25uLnZ2dqxYsQIbGxv69u3LtWvX0Ol0\nnDt3ji1btpCQkIC3tzdz586le/fuZGRkMGnSJDIzM5k5cyYFBQUAREdHM3bsWDO/K8VUVABQlDoq\nJyeHsLAwnnrqKRwdHcnNzSUxMZHz589z7tw5YmJiWLduHVlZWdja2lJYWFjpftzd3bn33nuJiIgg\nKCgIW1tbGjduTEZGBs7Ozjz88MP6tL/++iuLFi3C2toaKB6hcsWKFaSlpdGoUSNmz57NypUrWbp0\nKdbW1syYMYPZs2cTHR3Nhg0bGD58OH379iUpKYlTp07h7OyMt7c39vb2Jqkz5faoAKAoJiSlpLCw\nkIyMDJKSkvj555/5888/6d69O7GxscTHx3P16lWaNWvG+fPnyczM5Nlnn0UIUWG45saNGzNixAi6\ndetGamoqjz32GC1atODYsWM0atQIe3t7OnfujIeHB0IIUlJS2Lx5M6dPn+bq1atMnDiR4OBg1q1b\nx7Zt2zh69Cjx8fG4urrSunVr+vTpox+auNRHH30EFA8fnZGRwSuvvMLGjRs5ePAgUkrmz5+Pvb09\neXl5+tc0a9aMXr164efnxz/+8Q9atGhBREQEaWlpAJw5c4bPPvuMDz/8UE26ZGJ1+iqgHj16SHUS\nuG5SdVQ5KSXJycns3buXdu3acfXqVY4fP87evXvJyMjg4sWLXLt2rcrXe3t7o9FoCA4OxtnZmfbt\n29O4cWOklDRp0gQXFxe8vb3x9vamTZs2Bh+tdu3ataxfv56cnByOHTtGcnIyS5cuZdy4cVy7do3P\nPvuM3NxcPv74Yxo3bszXX3/N4sWLOXv2LIsXL0ar1bJgwQI6duzIK6+8wunTp9m0aROJiYlcunRJ\nP9uVTqcDiieg2blzJ1D8nerfv79B348lu8uTwFG1GUhOBYBKqMatZg29jpKSkrhy5Qq5ubmkpqYS\nFRVFdHQ0u3fv1v+yLSWEICAggFatWuHi4kLHjh1xcHDA1dUVb29vPD090Wg0tG3bFhcXFzO9o8pp\ntVp9d1B1yl6yWFBQgBACW1vbcmliY2NZvHgxrVu3pkmTJuzevZvLly/j6urK/v37mTRpEg8++CBP\nPvmkGoYd0wQA1QWkKFXQarWcOXOGy5cvk5WVxalTp9i3bx8ajYZjx45VSN+sWTP69OmDo6MjLVu2\nZNCgQfoGv3Xr1tXm1bJlS2O9jbtSm8YfKHe9elXdON27d2fNmjX69bIN3IMPPsi3337Lt99+y6hR\no0hMTKyxzpS7pwKAolD8qzUmJoY9e/aQkpLCkSNHOHr0KCkpKeXSeXt707FjR/z9/XFwcCA4OJh7\n7rmHv/zlLzRv3lyfrqEfId2uJUuW8Nxzz7F//34A/vjjD0aNGmXmUtV/KgAoDYpOp8PKyorLly/z\n888/s3fvXnbu3ElycnK5k6zOzs4MHjyYgQMH0rp1a9q2bUvbtm1xc3MzY+nrr65du7Jv3z5iYmII\nCAjg5s2b5i5Sg6ACgFKvSSk5ePAgBw4cYMeOHfz++++0bNmS5ORkdDodrVu3pm/fvvj4+NC1a1cC\nAwNp3bo1jo6Oqh/aDLp27YoQgqSkJHMXpUFQAUCpN3Q6HdHR0URERBAdHc3p06fRaDRcv34dgNat\nW9O+fXv69OmDq6sro0ePJjAw0GRDH9RFhYWFaDSacpdtmoqTkxPx8fEVtv/66680adKk0ucakqrq\npyx7e3vc3NwqnHCvLRUAFIul0+nYu3cvYWFhHD58mKtXr+p/ObZo0YJevXrRtWtXBgwYQN++ffH0\n9MTGRn3ly9JoNDg6OurvFTClzMxMHB0dK2zXarXY29vj7e1t0vLUNVXVTykpJSkpKWg0Gjw9Pe8o\nD/XfoFiMvLw8tmzZwsaNG9m1axc5OTlkZGRgb29Pv3796NKlCw888ABDhgyhbdu25i6uRcjLyzNL\n418dW1tbioqKzF2MOk8IgbOzM8nJyXe8DxUAlDorJydH350TERHBnj17KCgooGnTpvTu3ZsOHTrQ\nt29fhg4dirOzs7mLa7HqUuMPYGNjQ05OjrmLYRHu9rNTAUCpMwoLC/ntt9/YsWMHJ0+eZM+ePeTn\n5wPQpk0bnn/+eYYPH87gwYNrfX26YnnUEYDp1BgAhBDuwGrAFZDAcinlJ0KIucCzQOnxx+tSyq0l\nr/knMAnQAi9KKX8r2f4Q8AlgDXwjpZxv2LejWJqrV6+ydetWtmzZwh9//EFaWhp2dnZ07tyZF154\ngb/+9a90796dVq1qnNtCqSdsbW3RarX6S3YV46nNEUAR8A8pZbQQwhGIEkJsL3nuIynl4rKJhRC+\nwGigC9AW2CGE6Fjy9BfAEEADHBFCbJZSxhnijSiWQUrJoUOH+O9//8tvv/3G0aNHAXBzc2PkyJEM\nGzaMIUOG0KxZMzOXVDEXGxsbnnvuOT777DP69Olj7uKYXUJCAvv37+fpp582+L5rDABSykQgseRx\nphAiHmhXzUtGAN9LKfOBC0KIc8D9Jc+dk1L+CSCE+L4krQoA9dyVK1fYuXMnu3bt4tdffyUpKQlr\na2uCg4N5//33GT58ON26datzfdENzcyZM4mJiTHoPv39/fn4449v6zWllzRqtVqDlsUcbh1LqbZj\nK5WVkJBAWFiYeQJAWUIIDyAAOAQEA9OFEH8DIik+SkilODgcLPMyDf8LGJdv2d6rkjymAFMAXF1d\niYiIuJ0iGkRWVpZZ8rUkNdXR9evX2bFjB/v37ycuLg4pJY6OjvTs2ZPx48cTEhKCg4MDADdv3uSP\nP/4wUclNw1K+Q05OTmRmZgLFw2EYutEtKCjQ7/9WWq2WEydOMHLkSPz9/YmNjcXHx0cfMLKzs8nM\nzGTWrFlER0eTm5vLiBEjeOONNwB466232Lp1KzY2NgwaNIh58+bx/PPP07hxY2JjY7lx4wZffPEF\n69at4/Dhw/To0YMvv/wSoMp9ViYqKopXX32VnJwc7Ozs+Pnnn9m8eTPR0dF8+OGHADz55JO8+OKL\nhISE0KZNG5555hkiIiL48MMPefbZZxk5ciTh4eG89NJLBAYG8o9//IOUlBQaN27MZ599RseOHXn+\n+edxdHTk6NGjXL9+nblz5zJy5Ehefvllzpw5g5+fH2PGjGH69OnlypeXl3fn3zUpZa0WwAGIAkaW\nrLtS3JdvBcwDVpRs/xwYW+Z13wJPlCzflNk+Dvi8ujyDgoKkOYSHh5slX0tSWR2dPXtWfvjhh7J3\n795SCCEB2b17d/n666/L2NhYqdVqTV9QM7GU71BcXJzZ8s7IyJAXLlyQgNy7d6+UUspnnnlGfvDB\nBzIwMFBu375dSillSkqKlFLKoqIi2b9/fxkbGytv3LghO3bsKHU6nZRSytTUVCmllOPHj5ejRo2S\nOp1Obtq0STo6Ospjx45JrVYrAwMD5dGjR6vcZ2Xy8/Olp6enPHz4sJRSyvT0dFlYWChXrlwpX3jh\nBX264cOH6z9zQK5fv17/XPv27eWCBQv064MGDZJnzpyRUkp58OBBOXDgQH3Zn3jiCanVauXJkyel\np6enlLL4uzR8+PAq67GyzxCIlLVo12t1BCCEsAV+BNZKKX8qCRxJZZ7/GthSsnoFcC/zcreSbVSz\nXbFAWVlZREdHs2rVKnbt2kVCQgJQPOrjm2++yUMPPaT6cJUaubu7ExwcDMDYsWP55JNPAPRXAv3w\nww8sX76coqIiEhMTiYuLw9fXF3t7eyZNmsTDDz9cblazRx55BCEE3bp1w9XVlW7dugHQpUsXEhIS\n8Pf3r3Sffn5+Fcp2+vRp2rRpQ8+ePQFqdW7K2tqaxx9/vNy20oHtsrKy2L9/P08++aT+udIr3QAe\ne+wxrKys8PX1vavr+2urNlcBCYp/xcdLKZeU2d5GFp8fAAgFTpQ83gyECSGWUHwSuANwGBBAByGE\nJ8UN/2jA8J1ailHpdDrWrFnDp59+SnR0NABNmjRhyJAhvPTSS4wYMeKO70pUGqZbz/2UXvmj1Wq5\ncOECixcv5siRI7Ro0YIJEyaQl5eHjY0Nhw8fZufOnWzYsIHPP/+cXbt2AdCoUSP9fkofl64XFRVV\nuc/bYWNjo5/UBij3ent7+wr9/KXDZet0Opo3b17luZay5ZUmmKulNtdYBVPcXTNICBFTsgwDFgoh\njgshjgEDgVkAUsqTwA8Un9z9FXhBSqmVUhYB04HfgHjgh5K0Sh0npeTw4cPMnz8fHx8fJkyYwMWL\nF5k2bRo///wzGo2GTZs2MXPmTNX4K7ft0qVLHDhwAICwsDD69u2rnzUsIyODpk2b4uTkRFJSEtu2\nbQOKf0mnp6czbNgwPvroI2JjY2udX1X7rEynTp1ITEzkyJEjQPHwDEVFRXh4eBATE4NOp+Py5csc\nPny4Vnk3a9YMT09P/vOf/wDF/1s1ld3R0bHK8yh3qzZXAe2l+Nf7rbZW85p5FJ8XuHX71upep9Qt\n586dY+3ataxevZo///wTKJ4L9q233qJNmzYMHDjQzCVU6oNOnTrxxRdfMHHiRHx9fZk6dSrff/89\nRUVFdO/enYCAADp37lyuqygzM5MRI0aQl5eHlJIlS5bUkMv/VLXPytjZ2bF+/XpmzJhBbm4ujRs3\nZseOHQQ+zWB3AAAgAElEQVQHB+Pp6Ymvry8+Pj4EBgbWOv+1a9cydepU3nvvPQoLCxk9ejTdu3ev\nMr2fnx/W1tZ0796dCRMmMGvWrFrnVaPanCgw16JOApvetWvX5Jdffil79+4tASmEkIMGDZKrVq2S\n58+f16dryHVUG5ZSP3XhJHCXLl0qPHf27FkZGRkpc3NzzVCyuiEjI6NW6Yx+Elip36SU/Pbbb3z4\n4Yfs2LEDAF9fXxYsWMDTTz+tJkFRTM7NzY20tDTS09Oxt7c3d3HqLRUAGrAbN26wdu1ali9fTlxc\nHG3btmXOnDk88sgjBAUFqRuzFKPz8PDgxIkTFbbb29tja2tLdna2ycoSGhrKhQsXym1bsGABQ4cO\nNVkZTE0FgAbo+PHjzJ8/nx9++IGioiICAwNZs2YNTz31VJUTeiv1l5SyTgZ7BwcHMjIyTDYm0MaN\nG42eh6HJu7xSSAWABqKoqIitW7fy7bffsnnzZhwcHJgxYwYTJkyo9PpnpWGwt7cnJSUFZ2fnOhcE\nWrVqRWpqKsnJybi6upq7OHWOLJkQ5m66yFQAqOeysrJYtWoVn3zyCWfPnqVFixa8/fbbzJgxgxYt\nWpi7eIqZubm5odFoTHLT0a3y8vJqbLyysrKIiYmhbdu2DW42t9rUT+mUkHeqYdVoA3L9+nWWLVvG\nxx9/TFpaGvfffz//+c9/ePTRR1U3j6Jna2trtns3IiIiCAgIqDZNo0aN8Pf3p3Pnzqxbtw4PD48G\nMxdEbernbqnBtuuZs2fP8uyzz+Lu7s7cuXPp168f+/fv59ChQzzxxBOq8Vcsyn333cfq1as5duwY\n3t7eeHt7M23atAY/YbyhqABQDxQVFfHTTz8xZMgQOnXqxJo1a5g4cSLx8fH897//pXfv3uYuoqLc\nsccee4zIyEjeeOMNunTpwooVK+jatStLly41d9EsnuoCsmCnTp1i6dKlhIWFkZKSgru7O3PmzOH5\n55+ndevW5i6eohhM165dee+994Di7s2JEyfywgsvkJGRwauvvlrnTmBbChUALMzly5fZsGED69ev\n59ChQ9jZ2TFy5EhGjx7Nww8/3GD6R5WGy8XFhU2bNjFu3Dj++c9/kpCQwLJly1QQuAMqAFiAuLg4\n/SQU+/fvByAgIID58+czYcIEdYmc0uDY2Niwdu1a3N3dWbRoEa1ateLdd981d7EsjgoAdZCUkpiY\nGH766SfWr1/P2bNngeJBod577z2efPJJOnbsWMNeFKV+s7KyYsGCBdy8eZP33nsPKSXvvPOOmkj+\nNpg8AAghHgI+oXg2sW+klPNNXYa66tq1a6xatYqVK1dy+vRprKysGDhwIC+99BKhoaG0bdvW3EVU\nlDpFCMGXX36JlJJ58+bx559/8s0339CkSRNzF80imDQACCGsgS+AIRTPCXxECLFZSlmnJobfuHEj\nf/zxB2+99ZZJ8jty5Ajz5s3jl19+oaioiL59+zJ79mweffRRXFxcTFIGRbFUNjY2fPPNN3To0IHX\nX3+dI0eO8Mcff6gfTLVg6iOA+4FzUso/AYQQ3wMjKJ48xmCklOTk5NQ4FKpOp6t0+6effgrA3Llz\nOXDgAHZ2djRu3JgmTZrQqFEjGjVqhLW1NVZWVhX+1ubwMycnh8TERGJjY1m6dCk7d+6kefPm/P3v\nf+eZZ56hc+fOhqwORan3hBC89tpreHh4MGbMGLZt28akSZPMXaw6z9QBoB1wucy6Buhl6EwMOXbI\n7V5DX3oVjo2NjT6glG7XarVYWVmVmwO0bdu2LFq0iClTptRqvlFFUapWOjfwsmXL8Pf3JzAwUF0d\nVA1xt6PJ3VZmQjwBPCSlnFyyPg7oJaWcXibNFGBKyWon4LTJCvg/LYEbZsjXkqg6qp6qn5qpOqre\n3dRPeyllq5oSmfoI4ArgXmbdrWSbnpRyObDclIW6lRAiUkrZw5xlqOtUHVVP1U/NVB1VzxT1Y+rr\npY4AHYQQnkIIO2A0sNnEZVAURVEw8RGAlLJICDEd+I3iy0BXSClPmrIMiqIoSjGT3zEhpdwqpewo\npfSSUs4zdf61ZNYuKAtRL+tICNFcCDGt5HFbIcSGO9xVvawfA1N1VD2j149JTwIrSl0nhPAAtkgp\nu5q5KIpidGooCEUpbz7gJYSIAc4CPlLKrkKICcBjQFOgA7AYsAPGAfnAMCnlTSGEF8U3O7YCcoBn\npZSnTP82FKVmatAMRSnvNeC8lNIfePmW57oCI4GewDwgR0oZABwA/laSZjkwQ0oZBMwG1KD1Sp1l\nsC6gmsb4KfkFtYj/Xfb5uZTym+r22bJlS+nh4WGQ8t2O7OxsmjZtavJ8LYmqo+qp+qmZqqPq3U39\nREVF3TDZfQC3McbP+rI3fdXEw8ODyMhIQxTxtkRERDBgwACT52tJVB1VT9VPzcxRR1qtltzcXP1S\nUFBAUVGRfiksLCy3Xtm2mtaLior0+VX2A1sIgRACKysr/d+yQ8qULjdu3ODNN9+8o/cphLhYm3SG\nOgdgkjF+FEVpOAoKCkhPTyctLY20tLRyj8uuZ2dnl2vUc3Jyyq2X3V5YWGjut1Vrvr6+dxwAastQ\nAaC2Y/w8LoToB5wBZkkpL1eSRlEUCyelJDc3t9rG+9ixY3z//fdVNuy5ubnV5iGEwMnJiaZNm9Kk\nSRMaN26sX1xcXCpsK13Kbrezs8PGxgZbW1tsbGz0y92sW1tblxt/qOzjygakLF20Wq1+0el0HDx4\n0GifTylTXgX0M7BOSpkvhHgOWAUMujVR2bGAXF1diYiIMGERi2VlZZklX0ui6qh69al+Shvz1NRU\n0tLS9H/T0tK4efOmvtHOzs4mKytLv2i12mr3a2Njg4ODQ7nF1dUVLy8vHBwcaNq0qX57ZY8bN25s\n1MlfpJQUFhZSWFhYYzAyBhsbG6N/hwwVAGozxk9KmdVvgIWV7ajsWEA9evSQt/YRFhYWotFoyMvL\nu/tSV8HJyQl7e3uj7b8+KFtH9vb2uLm5YWtra+ZS1R2WcA5Aq9WSlJTElStXuHLlComJiSQmJnLt\n2jX9kpSUxPXr16v8f3NycsLFxYWWLVvi6upK8+bN9YuTk1O1jw8ePMjAgQNN/K4thym+Q4YKAPox\nfihu+EcDT5dNIIRoI6VMLFl9FIi/k4w0Gg2Ojo54eHgYbZjXzMxMHB0djbLv+qK0jqSUpKSkoNFo\n8PT0NHexlBI5OTn6hv3KlStoNJpy66UN/q2/0oUQuLi40Lp1a1q3bo2Pjw+urq64uLhUWFq1akWj\nRo3uuIxqmGbzM0gAqGqMHyHEO0CklHIz8KIQ4lGgCLgJTLiTvPLy8oza+Cu3RwiBs7MzycnJ5i5K\ng6HVaklMTOTixYtcunSJixcv6h9funQJjUZDWlpahdc1a9aMdu3a0a5dOwYPHoybm5t+vV27drRp\n0wYXFxdsbNT9oQ2FwT5pKeVWYOst2+aUefxP4J+GyEs1/nWL+jwMq6ioCI1Gw4ULF/RL2YZeo9GU\nu9QQwNnZmfbt2+Pl5UX//v1p165dhQbewcHBTO9IqatUqFcUM8jNzeXChQucO3eO8+fPc/78ef3j\nhISEcg28lZUV7dq1o3379gQHB3PvvffSvn172rdvz7333su9996rGnfljqgAoChGkp6eXmkDf/78\neTQaTbm0Tk5OeHt7ExgYyJNPPsl9992Hp6cnnp6euLu7qxPsilGoAHAH0tLSCAsLY9q0aQbZX58+\nfdi/f3+F7RMmTODhhx/miSeeqPK1NaXZs2cPzz//PLa2tvzyyy+89NJLbNhQ9QjHCQkJPPzww5w4\ncaLCcwMGDGDx4sX06KEmcSqVn5/P+fPnOXPmTLnlxIkTpKamlkvbunVrvLy8GDx4MF5eXnh5eeHt\n7Y2Xlxf33HOP6kpTTE4FgDuQlpbG0qVLax0AioqKqj2xVlnjbyhr167ln//8J2PHjgWotvFXqpaW\nlkZcXBzx8fGcOnVKv1y4cKHclTSurq507NiR3r17069fPzp27IiXlxf33Xef6qZR6hyLDgAzZ84k\nJibGoPv09/fn3XffrTbNa6+9xvnz5/H392fIkCEAbNu2DSEEb775JqNGjSIiIoJ//etftGjRglOn\nTnHmzBmWLFnCihUrAJg8eTIzZ84EwMHBgaysLKSUzJgxg+3bt+Pu7o6dnZ0+z3feeYeff/6Z3Nxc\n+vTpw1dffVXjL8ZvvvmGH374gd9++41t27Yxb948/a97rVbLa6+9RkREBPn5+bzwwgs899xz5V6f\nm5vLM888Q2xsLJ07dzbLzTCmlpaWxsmTJ4mLi+PkyZP6x1evXtWnadSoER07diQgIIDRo0fTqVMn\nOnfuTIcOHXBycgIs4z4ARbHoAGAu8+fP58SJE8TExPDjjz/y5ZdfEhsby40bN+jZsyf9+vUDIDo6\nmhMnTuDp6UlUVBQrV67k0KFDSCnp1asX/fv3JyAgQL/fjRs3cvr0aeLi4khKSsLX15eJEycCMH36\ndObMKb6oaty4cWzZsoVHHnmk2nJOnjyZvXv36ruIEhIS9M99++23ODk5ceTIEfLz8wkODubBBx8s\nF1SWLVtGkyZNiI+P59ixYwQGBhqqCs0uNTVV38iXbewTExP1aZo0aYKvry8PPPAAXbp0oUuXLvj4\n+NC+fXusra3NWHpFMQyLDgAff/yxUfabmZlZ67R79+5lzJgxWFtb4+rqSv/+/Tly5AjNmjXj/vvv\n198ctXfvXkJDQ/XDu44cOZI9e/aUCwC7d+/W76tt27YMGvS/kTLCw8NZuHAhOTk53Lx5ky5dutQY\nAKrz+++/c+zYMX2XUHp6OmfPnqVjx47lyvPiiy8C4Ofnh5+f3x3nZy7Z2dmcPHmS48ePc/z48Uob\n+qZNm+Lr68uDDz5Ily5d8PX1pUuXLtx7771GHWpAUczNogNAXWeosc7z8vKYNm0akZGRuLu7M3fu\n3LseCkNKyWeffcbQoUPLbS97lGBJdDodFy5c4NixY+WW8+fP64fkLW3ohw4dqm/kfX19VUOvNFjq\nW38HHB0d9UcJISEhrF+/Hq1WS3JyMrt37+b++++v8JqQkBA2bdpETk4O2dnZbNy4kZCQkHJp+vXr\np99XYmIi4eHhAPrGvmXLlmRlZRnkRO7QoUNZtmyZfnjcM2fOkJ2dXaE8YWFhAJw4cYJjx47ddb6G\nkJ+fz9GjR1mxYgUvvvgi/fr1o0WLFnh7ezNy5Ejefvttjh8/jr+/P3PnzmXjxo2cO3eOjIwMDh8+\nzMqVK3n55ZcZNmwYHh4eqvFXGix1BHAHnJ2dCQ4OpmvXrvz1r3/Fz8+P7t27I4Rg4cKFtG7dmlOn\nyk8DGxgYyIQJE/TBYfLkyeW6fwBCQ0PZtWuX/ldp7969AWjevDnPPvssXbt2pXXr1vTs2fOu38Pk\nyZNJSEggMDAQKSWtWrVi06ZN5dJMnTqVZ555Bh8fH3x8fAgKCrrrfG9XamoqsbGxHD16lJiYGGJi\nYoiLi9PfKNW0aVO6d+/O2LFj6d69O/7+/nTp0kXNNKUotWCwKSGNoUePHvLWGcHi4+Px8fExar5q\nMLia3VpHhvhckpOTiY6OJioqSv+3bJdUmzZt8Pf31y8BAQF4eXnVyV/w6iqgmqk6qt7d1I8QIkpK\nWeMNO+oIQDGLpKQkoqKi9Et0dDSXL/9vfiBvb2/uv/9+nnvuOQICAvD398fV1dWMJVaU+kcFgHoi\nNDSUCxculNu2YMGCCid5zeHmzZtERkZy5MgRIiMjiYyM1A+FIISgY8eOhISEEBgYSFBQEP7+/jRv\n3tzMpVaU+s8iA4CUUt02f4uNGzeaLe+y3Yg5OTlER0dz5MgR/XLu3Dn98x06dKBfv3706NGDoKAg\nAgICVHebopiJxQUAe3t7UlJScHZ2VkHAzEqnCrx27RpnzpxhzJgx+ruMAdzd3enZsyeTJk2iZ8+e\nBAUFqV/2ilKHWFwAcHNzQ6PRGHUCkry8PDUlZCW0Wi35+fnk5+eTl5dHYWEhWq2Wc+fO8fHHH9Ox\nY0ceffRRevbsSc+ePWndurW5i6woSjUsLgDY2toaferBiIiICpdoNjSFhYXExMSwf/9+Dhw4wIED\nB7h06RJQPFm1l5cXDzzwAL169eKhhx5ixowZ6ohMUSyMxQUAxThu3LjBgQMH2L9/P/v37+fIkSP6\nwd/c3d3p3bs3M2fOpFevXgQEBHDo0CF1CZ+iWDgVABognU5HXFycvrHfv38/Z8+eBYp/3QcGBvLc\nc8/Rp08fevfujZubm5lLrCiKMagA0ABkZGRw6NAh/S/8gwcPkp6eDkCrVq3o06cPkydPpk+fPgQF\nBdG4cWMzl1hRFFNQAaAe0mg07Nmzh71797Jv3z6OHTumv3S2W7dujBkzht69e9OnTx+8vLxU372i\nNFAqAFg4KSWnTp1iz549+uXixYtA8UQzvXv3Zs6cOQQHB9OrVy+aNWtm5hIrilJXqABgYQoLCzl6\n9Ki+sd+7dy8pKSkAuLi4EBISwqxZswgJCcHPz6/aqSgVRWnYVOtQxxUUFHDkyBHCw8OJiIjgwIED\n5OTkAODl5cUjjzxCSEgIISEheHt7q+4cRVFqTQWAOqaoqIjIyEjCw8MJDw9n3759+ga/e/fuTJo0\niZCQEPr27UubNm3MXFpFUSyZCgBmptVqOXr0qL7B37NnD1lZWQB06dKFiRMnMnDgQPr374+zs7OZ\nS6soSn2iAoCJ6XQ6YmNj9Q3+7t27ycjIAKBz586MGzdO3+C7uLiYubSKotRnKgAYmU6n4+TJk/oG\n/48//iA1NRUoHvN+1KhRDBw4kAEDBqguHUVRTEoFAAOTUhIfH1+uwb9x4wYAnp6ehIaG6ht8dYet\noijmpALAXZJScvbsWX2DHxERQVJSElA8hs6wYcMYOHAgAwcOpH379mYuraIoyv8YLAAIIR4CPgGs\ngW+klPNveb4RsBoIAlKAUVLKBEPlbypSSv788099Yx8eHs7Vq1cBaNu2LQ888AADBgxg4MCB3Hff\nfeqyTEVR6iyDBAAhhDXwBTAE0ABHhBCbpZRxZZJNAlKllN5CiNHAAmCUIfI3tpycHMLDw9m6dStb\nt27VT1Tu4uKi/3U/cOBAOnTooBp8RVEshqGOAO4Hzkkp/wQQQnwPjADKBoARwNySxxuAz4UQQpad\nT7AOSUhI4JdffuGXX34hPDycvLw8mjRpwgMPPMDs2bMZOHAgPj4+qsFXFMViGSoAtAMul1nXAL2q\nSiOlLBJCpAPOwA0DlUHvxo0bdzVpTFFREXl5eUDx3bZTpkxh+PDh9OvXT80UpihKvVHnTgILIaYA\nU0pWs4QQp81QjJaUBKbz58/z6aef8umnn5qhGHWavo6USqn6qZmqo+rdTf3U6ooTQwWAK4B7mXW3\nkm2VpdEIIWwAJ4pPBpcjpVwOLDdQue6IECJSStnDnGWo61QdVU/VT81UHVXPFPVjZaD9HAE6CCE8\nhRB2wGhg8y1pNgPjSx4/Aeyqq/3/iqIoDYFBjgBK+vSnA79RfBnoCinlSSHEO0CklHIz8C2wRghx\nDrhJcZBQFEVRzMRg5wCklFuBrbdsm1PmcR7wpKHyMzKzdkFZiHpZR0KI5sDTUsqlQoi2wKdSyifu\nYFf1sn4MTNVR9YxeP0L1wijK/wghPIAtUsquZi6KohhdnbsKSFHMbD7gJYSIAc4CPlLKrkKICcBj\nQFOgA7AYsAPGAfnAMCnlTSGEF8U3RbYCcoBnpZSnTP82FKVmhjoJrCj1xWvAeSmlP/DyLc91BUYC\nPYF5QI6UMgA4APytJM1yYIaUMgiYDSw1SakV5Q7U6S6gli1bSg8PD5Pnm52dTdOmTU2eryVRdVQ9\nVT81U3VUvbupn6ioqBtSylY1pavTXUAeHh5ERkaaPN+IiAgGDBhg8nwtiaqj6qn6qZmqo+rdTf0I\nIS7WJp3qAlIURWmgVABQFEVpoFQAUBRFaaDq9DmAyhQWFqLRaPSjdRqDk5MT8fHxRtt/fWCqOrK3\nt8fNzQ1bW1uj56UoDY3FBQCNRoOjoyMeHh5GG4s/MzMTR0dHo+y7vjBFHUkpSUlJQaPR3NXw3oqi\nVM7iuoDy8vJwdnZWE7E0AEIInJ2djXq0p9Q9qamp/Pvf/0an05m7KPWexQUAQDX+DYj6rBued999\nl3HjxpnlEvCGxiIDgKIo9VfpuaUrV26dUkQxNIMFACHEQ0KI00KIc0KI1yp5/u9CiDghxDEhxE4h\nRK1mrFEUpWGxtrYGiruCFOMySAAQQlhTPADWXwFfYIwQwveWZEeBHlJKP4onhV9oiLzrEg8PD27c\nqP0MbhEREezfv/+283FwcADg6tWrPPHE7Y1UPGDAAIMeWk+YMIENGzbcdRpFKVUaANLS0sxckvrP\nUEcA9wPnpJR/SikLgO+BEWUTSCnDpZQ5JasHKZ42skG70wBQqm3btqphVeotFQCMz1CXgbYDLpdZ\n1wC9qkk/CdhW2RNlJ4V3dXUlIiKi3PNOTk5kZmYC8Oqrr3L8+PE7LnRlunXrxvvvv6/PozIXL15k\n5MiR+Pv7Exsbi4+PD1999RVSShYvXsyvv/5KYWEhq1evpmPHjty8eZMXXniBhIQEGjduzKeffoqj\noyPLli3D2tqa1atXs2jRItq1a8cLL7xASkoKLVu2ZOnSpbi7u5OQkMCkSZPIzs5m2LBhQPFlmBcv\nXuSpp57i0KFD5OXlMWvWLI4ePYqNjQ3vv/8+/fr1Izc3l6lTp3LixAk6duxIVlYW2dnZZGZmsnPn\nTt5//30KCgrw9PRk6dKl+qOLW82fP59t27aRl5dHr169WLJkCZmZmRQWFpKbm0tmZiZdu3YlNDSU\n7du3Y29vz7fffouXlxeFhYXs2LGDRYsWcf36dd555x0ee+wxsrKyGDNmDGlpaRQWFvKvf/2L4cOH\nV8g7Ly+vwvegrsvKyrK4MptaVXV0+XJxUxIfH9+g69AU3yGT3wcghBgL9AD6V/Z82Unhe/ToIW8d\nDCk+Pl5//bmdnZ3+cNFQSvdZ3TXuDg4OnD17lpUrVxIcHMzEiRNZs2YNQgjatWtHTEwMS5cuZdmy\nZXzzzTe8/vrr9OzZky1btrBr1y6mTp1KTEwMU6dOxcHBgdmzZwPwyCOPMHHiRMaPH8+KFSt4/fXX\n2bRpE2+88QbTp0/nb3/7G1988QUAjo6OODg4YGVlhaOjI8uXL8fOzo6TJ09y6tQpHnzwQc6cOcPX\nX3+Nk5MTp0+f5tixYwQGBtK0aVPy8/NZsmQJ4eHhNG3alAULFvD1118zZ86cSt/zP/7xD+bNmwfA\nuHHj2L59O0899RS2trY0btwYR0dHhBC0atWKkydPsnr1at588022bNmCra0tKSkpHDhwgFOnTvHo\no48ybtw4GjduzObNm2nWrBk3btzgL3/5C6NGjapw5Y+9vT0BAQGG+HhNRg10VrOq6qj0f9rJyalB\n16EpvkOGCgBXAPcy624l28oRQjwAvAH0l1Lm322mH3/88d3uolLV/fov5e7uTnBwMABjx47l008/\nBWDkyJEABAUF8dNPPwGwd+9efvzxRwAGDRpESkoKGRkZFfZ54MAB/WvGjRvHK6+8AsC+ffv0rx83\nbhyvvvpqhdfu3buXGTNmANC5c2fat2/PmTNn2L17Ny+++CIAfn5++Pn5AXDw4EHi4uL076GgoIDe\nvXtX+X7Dw8NZuHAhOTk53Lx5E29v70rTjRkzRv931qxZ+u2PPfYYVlZW+Pr6kpSUBBTf6PX666+z\ne/durKysuHLlCklJSbRu3brKcij1X+n/X2X/I4phGSoAHAE6CCE8KW74RwNPl00ghAgAvgIeklJe\nN1C+ZnPrr9TS9UaNGgHFv2KKioqMlt/dklIyZMgQ1q1bV2PavLw8pk2bRmRkJO7u7sydO7fKm7PK\nlrPs49J6Kc0bYO3atSQnJxMVFYWtrS0eHh7qpi9F3/DX5oeYcncMchJYSlkETAd+A+KBH6SUJ4UQ\n7wghHi1JtghwAP4jhIgRQmw2RN7mcunSJQ4cOABAWFgYffv2rTJtSEgIa9euBYoP61q2bEmzZs1w\ndHQs9yXv06cP33//PVDcOIaEhAAQHBxcbntNeZw5c4ZLly7RqVMn+vXrR1hYGAAnTpzg2LFjAPzl\nL39h3759nDt3DiiefOLMmTOV7ru0UW7ZsiVZWVnVnnhev369/m91RxQA6enpuLi4YGtrS3h4OBcv\n1moIc6WeUwHAdAx2H4CUcquUsqOU0ktKOa9k2xwp5eaSxw9IKV2llP4ly6PV77Fu69SpE1988QU+\nPj6kpqYyderUKtPOnTuXqKgo/Pz8eO2111i1ahVQ3Oe/ceNG/P392bNnD5999hkrV67Ez8+PNWvW\n8MknnwDwySef8MUXX9CtW7cqb46ZNm0aOp2Obt26MWrUKL777jsaNWrE1KlTycrKwsfHhzlz5hAU\nFARAq1at+O677xgzZgx+fn707t2bU6cqn7q2efPmPPvss3Tt2pWhQ4fSs2fPKt9ramoqfn5+fPLJ\nJ3z00UfV1uH//d//ERkZSbdu3Vi9ejWdO3euNr1S/2m1WrKzswHVBWQSUso6uwQFBclbxcXFVdhm\naBkZGdU+f+HCBdmlSxejl6Muq6yO2rdvL5OTkw2elyk+c0MLDw83dxHqvMrqKC0tTQISkB4eHqYv\nVB1yN98hIFLWoo1VQ0EoilJnlP7qt7OzU11AJmBxw0HXBR4eHpw4ccLcxTCK0NBQLly4UG7bggUL\nGDp0aI2vTUhIMFKplIaitNFv164dGo3GzKWp/ywyAEgp1SiRRrJx40ZzF6EcWXLFkNIwlB4BtGvX\njgsXLpCfn1/uCjLFsCyuC8je3p6UlBTVMDQAsmRCGHt7e3MXRTGRskcAoE4EG5vFHQG4ubmh0WhI\nTk42Wh55eXmq0amBqeqodEpIpWEoewQAxQGhVatW5ixSvWZxAcDW1tbo0wNGRERY3NADpqbqSDGG\nylY2WEkAABbVSURBVAKAYjwW1wWkKEr9pbqATEsFAEVR6gx1BGBaKgAoilJnZGRkYG9vzz333AOo\nAGBsKgAoigJAYmIiW7duNWsZMjMz9eNkgeoCMjYVABSljIKCAq5fLz9YrU6nu61RShMSEtDpdOh0\nOnJzc28r/9TUVL777rvbfp0hTJgwgeHDhxMVFWXyvEtlZGTg6OhIs2bNAHUEYGwqACh1UmFhYa3T\n5ufn88033+jnY87MzCQqKgopJTqdjldeeYW33nqL8PBwgoOD6dmzJ6dOnSIqKorff/+d5ORklixZ\nwqFDhxg6dCgeHh78+OOPXL9+nZycHPr27YuDgwOLFy8mPj6ewYMH8/DDD/PWW28xb9484uLi9Pm9\n8cYbeHp68vnnnzNkyBCcnJzYvn07Bw8eZNq0aaxcuZLt27fz9ttvs2zZMpKSkvD39ycwMJBLly4x\nffp0nnnmGf1cEJs3b+bpp5+u1QTpqamprFq1ivT09Eqfz8vLQ6vVVvqclJLff/8dgJ07d9a67g0t\nLS2N5s2b62emUwHAyGozYJC5lsoGgzMFNZBXze6kjjQajRw+fLjcv39/ue1ZWVly48aNMisrS0op\n5ebNm6W1tbVcsWKFjI2NlS4uLnLGjBlSp9PJa9euyXfeeUe+/vrrcsaMGfLAgQNy/PjxEpDBwcHy\ngw8+kM2bN5eA7NSpk7znnnv0g4sB0tPTUzo4OEgHBwdpbW1d7rnKFmtraymEkAEBAfpt9vb2laZt\n2rRppdubNGki7ezsqsxDCCEB+eCDD+rT2draytOnT8s2bdpIQE6ZMkVKKeWZM2fklStXytXftWvX\npJRSjho1SgLyscceq1D3+/btk61atZIjRoyo9LO5ceOGvjyjR4++7c/2TlT2HerTp48cPHiwlFLK\npk2byr///e8mKUtdZIrB4MzeyFe3qABQd8TExMj8/Hz9emV19Mcff8jHH39choaGyvDwcLlu3TqZ\nlZUlf/nlF9m3b1/p7Oysb5hPnDghJ0yYID/44AM5fPhwCch7771XXrx4UXbu3FnfMJa+BpAzZsyQ\n7f6/vTOPjqpM8/DzVhKJYCAkMI0IBxNZZFFEYWwUEQE1LCIgSzPKDGqjQDPYjtIgbiNoIwjK0tiK\nCI22zKBwFJqtYUA8YhNBgtBpD0iw2RohbEHWQFK/+eMW1QlCEiGpqqS+55w65y7fvfe9v3rvfb/t\nft811xR6OZ9bvu6664LLnTt31qBBg5SYmKhu3brp7bffVocOHXTDDTdo3759mjlzpgDFxsbqlVde\nUZcuXbR06VJNmjRJy5cv1759+zRr1iwNHTpUlSpV0vjx45WXl6fhw4frxhtv1NatW7V27VrNnTtX\nGRkZmj17tsaPH6+EhASlpaUpJydHTZs2Vb9+/bRr1y7ddNNNatiwob777jt9/fXXWrBggXJycjRh\nwgQ1adJECxcu1LBhw4L2L168WD6fL7jeuHFjAapataoA1ahRQ1lZWZKk1157TYCefPJJxcXFBY+Z\nPHmyOnXqpEGDBmnnzp2qWbNmcN+5AJydna0bbrhBffv21ZdffilAPp9PjRo1ColPXciHmjZtqgce\neECSVKtWLQ0cOPCix/v9/rIyLSJwAcAFgJDi9/t18uRJnT17Vjk5OZKkw4cPB3PYgwcP1pkzZ7Rn\nzx4NHjxYDzzwgLp166YXXnhBt99+uwDVrFmzUG63bdu2Sk5OVo0aNdSmTRs99NBDF8wF9+7du9Bx\nY8aMCb7gV69erT59+ghQQkKCMjIydPr0aeXk5Ojll1/WxIkTlZ+frw8//FBz5swp9sWQn5+vSZMm\nafXq1cVqkpeXV2L9Tpw4Ebz2T/WhHTt2KD4+Xrfffrsk78XeuXNnvf766zpw4ICSk5OVlJSk4cOH\nKyEhQXXq1NG0adN05ZVXFtJx48aNSklJCZYiCpZC1q9fr+TkZHXt2lWS9NxzzwX3nys99OnTR2am\nI0eOFLLP7/dr9+7dxWp79uzZEt9zQY0yMzP1ww8/6JprrtGjjz4qSWrSpMkFSzN+v1+TJ09W1apV\nNXz4cG3ZsqVCBgMXAEIcAAYPHqyHHnoo6gJAVlaWTp8+rf79+8vn8wWrONLS0oK58Qv9GjRooEaN\nGgVz72PHjtXx48e1bds2PfvssxoyZEgw7aZNmyR5L98pU6Zo6tSpmj9/vnr27KknnnhCfr9fH3zw\ngdq2batHHnlEfr9fGRkZwbkAjh07pokTJ2rjxo3hlKrEXIoP7d69W6dOnbrgvv379wfnWtiwYUOw\nxFO9enWlp6erdu3awdzykiVLlJKSorVr12ro0KGKjY3V+++/L0l66aWXBGjt2rVKSkpS165dg/9h\nTEyMli9fLkC//e1v9fzzz2vdunXKysrS3XffLUC9evXS2bNn9Yc//EFpaWlq2bKlXnzxRUnS4sWL\nFR8frxEjRvwkjQ4cOCBA7du3V+XKlfXUU09Jknr27Cmfz6cOHTpo1KhR8vv92r9/vzp16iRA1apV\nC/rXgAEDip3Ho7zhAkAIA0BeXl7QmSpyAFi9erX69++vRx55RMOHDw/mrAvm2IcNG6YRI0YEt82a\nNUsnTpxQjx49VKdOHY0aNUrTpk0LnvPIkSMXzCmfPn1aPXv2LJQ2WihrH8rNzdXSpUu1Y8eOYtOe\na1uRpEOHDgWrkgB98cUX+uijjwTo5ptvVm5ubrANBVDt2rV1xx13CNB9990nQGPHjlXlypXl8/mC\npbTNmzerTZs2wWqkPXv2SPIyF1OmTNGBAwc0e/bsYAD3+/3q3r27ZsyYoU8++aSQD06cOFGS9OKL\nLxba/s4776hhw4aKj4/XG2+8ofz8fH399dcaMmSIzEyNGzfWtm3bykDt8OACQAgDQEZGRtDRlixZ\nErLrlhX5+fnBHNHJkyf1zDPPqF69esF7TEhICOai2rVrp6uuukrXXnttoRzoN998o3nz5hU676VW\ncUQbkazPypUr1bJlSz3++OOSvP901qxZwWCSnp6u1157TdOmTQv6y7hx4+T3+9W8efNC1U0HDhxQ\nXFycBg4cqJiYGPXo0UOAJkyYoJYtW/6o1FipUiVlZGRozpw5wW0F2z8AzZ07V5K0a9cude3aVZ9/\n/rnq168frMpas2bNj+5p1apVSk5OVmJioqZMmXLRklSksnfvXh09elSS93+cPn3aBYBQBoBJkyYF\nHfD8l155Y9++fcGHLzU1VbGxsQJ01113qVWrVvr2228lSWfOnAkes3fv3p80nWMkv+AigYqgz8mT\nJ4PPRHZ2tiQFg0KzZs2C6bp27RpM99lnnyk1NTW4PnLkSD333HO65557NG/ePNWpU0fVq1f/Ua+o\nc72dAH311Vc/smXFihW67bbbisycfffdd8FSSEpKip5//nmNGDFCBw8eLH1xzsPv92vPnj3Kz8+X\n5GXA1qxZo/T0dP3lL3/R008/rTlz5mj69OkaP368nnrqKXXp0kWDBg1Sv379gj3BUlNTlZSUpCuu\nuEJ9+vQJnu+nUtIAUO5GAy0rVq9eHVw+Nyl1eeLNN99k3rx5NG/enHnz5nH48GEGDBjA8ePH6d27\nN2lpabRr167QMXFxccHlq6++OsQWOyKdK6+8kqVLl3LmzJngkMwPP/wwu3btonfv3sF0vXv3ZtGi\nRQC0bt2aLl26MHXqVFq3bs3YsWMLnbNevXq0atUKgGnTpjF+/Hh27tzJwIEDmTJlCjk5OTRr1uxH\ntnTs2JGOHTsWaW9KSgqff/45K1euZNiwYYwZMwaADz/8kPvvv58qVapwxx13UKtWLTZu3IjP5+Ou\nu+6ibt267Nu3j+3bt1O9enXq16+PmZGZmUmDBg04e/Yshw8fZsuWLVSuXJnMzExSU1PJysoiNzeX\nnTt3sm3bNlasWMGNN97I8ePHOXLkSJHfbpgZ119/Penp6eTl5dG1a1datmxJZmYmVatW5dSpUxw6\ndAifr4w/1SpJlCjJD0gDtgJZwMgL7K8EzA3s/xK4trhzlkYJYNWqVerfv7+2bt0qSdq8ebPefPNN\njRo1Sq+88oo2b96skSNHCgh2MXzrrbcu+7qhwu/3B4vT1atXV2xsrFq0aKF169aV6XUrQg63LIkm\nfY4cOSJAP//5zyVJ27ZtU+fOnbVhw4YLpl+wYIHmz5+vTz/9VMuWLVOvXr2UnZ2tDRs2aNmyZaVi\nU15eng4ePKiVK1cqKSlJ8fHxF/zuw+fzqX79+oV6TCUlJRXqflzcLzExUXXq1FH37t1VpUoVtWjR\nQn379tU777yj0aNH66WXXtKuXbu0ePFi7dixQzk5OSWqolq5cuUl3z+hLAGYWQwwDbgb2AOsN7OF\nkr4pkOxR4Iik+mb2C2Ac0Lc0rn8+ubm5jB49mvj4eMaNG8eJEyf4+OOPadiwIRs3bkTyppSUvC83\nY2JiuPPOOxkwYAAPP/wwJ0+eLAuzSg3v/4UFCxYwaNAg9u/fz6233sqqVauIi4sjNjbWTZnpCBmJ\niYls2bKF2rVrA1C/fn0WL1580fTdunUDvFL3vffeG5xvujQnfomJiSE5OZn27duTnZ2Nz+fjxIkT\npKen88MPP9CgQQN8Ph9z585l0aJFDB48mLS0NL7//numT59OamoqaWlpbN++ncTERBITE0lNTSU7\nO5sGDRqwf/9+WrRoQaVKlahWrVrwun6//6K59rp16/6keyjz3D+UTgkAaA38ucD6M8Az56X5M9A6\nsBwLHASsqPNeaglg7969wXrvpk2bas2aNXrwwQfVrl07/eY3v9HmzZuVnZ2tzMxMTZ06VZmZmZKk\ntWvXCtCrr756Sdcta/bt26cZM2aoVq1awa9Ob7nlFs2cObPQR1qhIJpyuJeC06d4nEZFE4pGYFMg\nN3k5mFkvIE3SLwPr/YFbJQ0tkCYzkGZPYH17IM3B8871GPBYYLURXrVSqKmBF6AcF8dpVDROn+Jx\nGhXN5ehTT1KxRaqIawSWNB2YHk4bzOwrSS3DaUOk4zQqGqdP8TiNiiYU+pRWJdM/gIIVXHUC2y6Y\nxsxigWrAoVK6vsPhcDh+IqUVANYDDcwsxcyuAH4BLDwvzULgPwLLvYBVKo36J4fD4XBcEqVSBSQp\nz8yG4jX0xgAzJf3NzEbjNUYsBN4F3jezLOAwXpCIVMJaBVVOcBoVjdOneJxGRVPm+pRKI7DD4XA4\nyh9uRjCHw+GIUlwAcDgcjiglqgOAuc9lHQ5HFBOVAcDMXjSzBq4X0sUxs85m9rMC6y5YFsDM/tvM\nWofbjkjG+VDRRIIPRVUAMLNfmNlq4AbgsHPIH2Nmt5nZTmAo8LaZ9QVwwdLDzHqZ2RdAE2B7uO2J\nRJwPFU0k+VDU9AIys+bARqCLpKXhtidSMbNfAXmS3jaz3sA9wBpJs83MJ8kfZhPDhpk1BFYAoyW9\nG9gWIyk/vJZFFs6HLo6ZNQKWEyE+FHFDQZQmZlYViAMOS9pkZkuAFDOrhDdg3U7g/4A90Zo7MbMr\nJJ0psKk13pDdAMuAPOBRM/uTpMMhNzDMFNRH0rdmtgCoambX4o1ZlWNmK4C/nadj1GBmCXjf/xwN\nPEdt+OcYXs6HCuuz1cwWEyE+VGGrgMzsCWAb8HtgUmDzE3jDUK8DrgDaA6PxhrGOOszsBWCBmbUp\nsPkt4B4zS5B0DEgHvgX6hMPGcHIRfSYCA4BFQD5QE3gSeDDkBkYAZvY03vP0e2ByYLPzoQAX0WcC\nEeJDFbIEYGYpQCegGXAW+MTMfilphpkNAHZLWmdmVfACQEr4rA0PZtYHeABYBdxnZl9KOgtkAhnA\nr4ExwFEgmwqcWbgQF9AnXVKepJ1m9iywV1KGmcUBjwP1L1CaqrCYWTwwArgF71k7BWw1s3clfWZm\nm/AyXC8ThT5UhD6zJW0IZC52h9uHKswfYmY1CqweBaoDlSXl4EXc28ysraT5ktYBSDqBN6tP3I9O\nWAE5T6MlwH3A/wJVgR6B7ceA94DuZtZK0km82dyqUcEpoT5IWiQpI7B8Fk+bE9Hw8jezc0MM5+IF\nx36SdkjajzfcS/vA/llAjyj0oaL0mUFAH0kLIsGHyn0AMLNKZvY74DMzG2NmHYEcvGJXE/AeWOAA\n0DpQ/4+Z/YuZjcOr/vkyPNaHhgtodI+k45J2AX8FNgHtzKyWpPxAgHwPGGtmc4F/owJrVEJ97jrX\npfFc7zEzqxvwoT541RwVlgIarTazl4E7JX0OnDIPH978HZsAotiHitLneuCr844Lqw+V+wCAN8Lo\nNUBbvDr/d/ByGweAm82sXiDdQqCvpNxA4/B04CrgdknrQ292SDlfoxmB6i8CubN04DQQnOlb0mS8\nBqrFQEtJq0JtdAgpqT59AtsUCALjgAQ8H6rI+kBhjb4F3jOzKoEePTEFevb8/dwBAR96nOjzoRLp\nE2A8/3wPhVyfctsN1Mws8CA+BiRJejWw/X/werF8gNew8qWkmYF9y4EhkrLMrJqko+GyPxQUodH7\nePMzDwusV8IrCXUAdgMGTJaUFybTQ8Jl6OOTNCHwgJ8Il/2h4CdolAx8KKmDmd0L1DnXzbEicxn6\n1A20SYbVh8pVCcC8iWSAQh+VVAWSzhXLgf8CHgJO4hVBu5jZW+Z138vFm7SeivryL6FGw/EaNhsH\n0uXiadMTGAbsqqgv/1LSZ2dge4V8+V+KRkBHIN7M3sVr+M0Olb2hppT02R84Prw+pFKYFL6sf3i9\nlSbgdaPqeN6+a/DqaW8qsG0C8G5g+Wd4D+0T4b6PCNNoPDA7sFwJ+BvwUrjvw+lT7jT6Y2D5MbwZ\n/n4d7vtw+pT8F/ElgEBEnQLUwmvYHWFmvzrXmCvpH8B8YLiZ1QoctgzYbd4XdvslTZFXH1khuUSN\nlgN/N7M4eTncWyS9GAbzyxynT/FchkbnhjJYC6RKmkQFpKLqE/FtAIEG22XAvZKOBerPOgPrJf0x\nkMbwuqDlAhuAIcD7kt4Ik9khxWlUNE6f4rkMjf4o6fUwmR0yKqo+EV8CkPQDsAPvyzmAL/DG9LnN\nzGoH0ggYhdfboAMwKVoeXHAaFYfTp3guQ6OIfbmVJhVVn4gPAAE+Bm4ys6slHQc243XLSwIws6bA\nIXkf6PST9F4YbQ0XTqOicfoUj9OoaCqcPuUlAKwBDhKIvvK+oPtX4Eoz6wa0gqgfb9xpVDROn+Jx\nGhVNhdOnXIwFJOn7QDfOV80sC1gPnAn8/qRIb8gIAU6jonH6FI/TqGgqoj4R3whcEDPrhPe16m3A\n7yT9LswmRRxOo6Jx+hSP06hoKpI+5SoAAJg3cp5UQT9UKg2cRkXj9Ckep1HRVBR9yl0AcDgcDkfp\nUF4agR0Oh8NRyrgA4HA4HFGKCwAOh8MRpbgA4HA4HFGKCwAORwHMLNHMhgSWa5vZvHDb5HCUFa4X\nkMNRADO7FlgkqVmYTXE4ypxy8SWwwxFCXgWuM7Ov8aaHbCypmZkNALoDVYAGeOPCXwH0xxv9sbOk\nw2Z2HTANqIk3KdFASVtCfxsOR/G4KiCHozAjge2SbsKb1akgzfBmBWsFvAKclNQCb6z3fw+kmQ78\np6RbgKeBN0NitcNxCbgSgMNRcj6VdAw4ZmZHgT8Ftv8VuNHMrsIbHuCjAuOBVQq9mQ5HyXABwOEo\nObkFlv0F1v14z5IPyAmUHhyOiMdVATkchTkGJFzKgYFJQ/5uZr3BGxbYzJqXpnEOR2niAoDDUQBJ\nh4AvzCwTeO0STvEg8KiZbcKbSP7+0rTP4ShNXDdQh8PhiFJcCcDhcDiiFBcAHA6HI0pxAcDhcDii\nFBcAHA6HI0pxAcDhcDiiFBcAHA6HI0pxAcDhcDiiFBcAHA6HI0r5f3Cm2l60tM4+AAAAAElFTkSu\nQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f4010893860>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "URL = 'http://golem.fjfi.cvut.cz/utils/data/{}/{}'\n",
    "# function for reading 1D y(t) signals\n",
    "def read_signal1d(shot_number, signal_id):\n",
    "    #url = URL.format(shot_number, signal_id)\n",
    "    url = signal_id\n",
    "    return pd.read_table(url,names=['time',signal_id],\n",
    "                         index_col='time')\n",
    "# read the specified signals\n",
    "shot_no = 29395\n",
    "U_l = read_signal1d(shot_no, 'loop_voltage')\n",
    "I_p = read_signal1d(shot_no, 'plasma_current')\n",
    "P_OH = U_l*I_p  # vectorized, time-aligned operation\n",
    "B_t = read_signal1d(shot_no, 'toroidal_field')\n",
    "H_a = read_signal1d(shot_no, 'photodiode_alpha')\n",
    "# combine into a data frame table\n",
    "df = pd.concat([U_l, I_p, B_t, H_a], axis='columns')\n",
    "# plot the data table in subplots from 4 to 25 ms\n",
    "df.loc[4e-3:25e-3].plot(subplots=True, ylim=(0,None))\n",
    "plt.show()   # display the figure in a window\n",
    "plt.savefig('graph.jpg')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "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.4.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
