{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "4509bb79",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAicAAAGdCAYAAADJ6dNTAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQAARLZJREFUeJzt3XtcVHX+P/DXyGVARRDRARQQM0tCQAcyMBNMsbHIVncf5O4ilprEmBhlRbRp/syxNQmLkc1bZNumeYnUKKRMschVUbzRelsUNJBAZURzUDy/P/xy1mkGZWRgzjCv5+NxHg/POZ/5nPfMAefN53wuMkEQBBARERFJRCdrB0BERER0KyYnREREJClMToiIiEhSmJwQERGRpDA5ISIiIklhckJERESSwuSEiIiIJIXJCREREUmKo7UDMNeNGzfwyy+/wM3NDTKZzNrhEBERUQsIgoBLly7B19cXnTrdvm3EZpITrVYLrVaLhoYGnDx50trhEBER0V2oqKhAnz59bltGZmvT19fV1cHDwwMVFRXo1q2btcMhIiKiFtDpdPDz88PFixfh7u5+27I203LSpOlRTrdu3ZicEBER2ZiWdMlgh1giIiKSFJtJTrRaLYKCghAREWHtUIiIiKgN2VyfE51OB3d3d9TV1fGxDhERkY0w5/vb5vqcEBG1lCAIuH79OhobG60dClGH5+DgAEdHR4tM88HkhIg6pIaGBlRWVuLKlSvWDoXIbnTu3Bk+Pj5wdnZuVT1MToiow7lx4wbKysrg4OAAX19fODs7c9JGojYkCAIaGhrw66+/oqysDPfee+8dJ1q7HZtJTpomYWPzLBHdSUNDA27cuAE/Pz907tzZ2uEQ2QVXV1c4OTnh9OnTaGhogIuLy13XZTOjddRqNUpLS7Fnzx5rh0JENqI1f7kRkfks9TvH31wiIiKSFCYnREREJCk20+eEiMgS0jYeatfracYPatfrEXUEbDkhIpKQwsJCxMXFwdfXFzKZDLm5ua2qT6/XIywsDDKZDCUlJeLx2tpaPPbYY/D19YVcLoefnx9mzJgBnU4nltm+fTvGjRsHHx8fdOnSBWFhYfj0009bdN2vvvoKQ4cOhaurK7y8vDB+/HiD89999x2ioqLg5uYGHx8fvPrqq7h+/bpBGUEQ8O6772LAgAFijAsWLDAoo9VqMXDgQLi6uuK+++7D6tWrjWK5ePEi1Go1fHx84OLigoEDByIvL8+gzNKlSxEYGAgXFxcolUrs3Lmz2fc2ffp0yGQyZGZmmjwvCAJUKpXJ+7dv3z6MHj0aHh4e6NGjB5577jnU19eL5w8cOICJEyfCz88Prq6uGDhwIJYsWWJQx9WrVzF58mQMGjQIjo6OeOqpp0zG8emnnyI0NFQc3vvMM8+gtrZWPB8dHQ2ZTGa0Pf7442KZuXPnGp339vZu9rOxFCYnREQScvnyZYSGhiIrK8si9b3yyivw9fU1Ot6pUyeMGzcOmzZtwrFjx5CTk4Nvv/0WSUlJYpmioiKEhIRgw4YNOHjwIJ599llMmjQJmzdvvu01N2zYgISEBDzzzDM4cOAAfvzxR/z5z38Wzx88eBBjx47FY489hv3792PNmjXYtGkTXnvtNYN6UlJSsGLFCrz77rv4z3/+g82bN+PBBx8Uz2dnZyMtLQ1z587FkSNH8NZbb0GtVhvE19DQgNGjR+PUqVNYv349jh49iuXLl6N3795imbVr12LWrFlIT0/H/v37MXz4cKhUKpSXlxu9t9zcXPz73/82+Zk2yczMNDl0/ZdffsGoUaPQv39//Pvf/8Y333yDI0eOYPLkyWKZ4uJi9OzZE//85z9x5MgRpKenIy0tzeDnobGxEa6urpg5cyZGjRplMoYffvgBkyZNwpQpU3DkyBGsW7cOe/bswdSpU8UyGzduRGVlpbgdPnwYDg4O+NOf/mRQ1wMPPGBQ7tChtm99tJnHOu01lPjWJl+N0wrDk3FLQETUllQqFVQqlUXq+vrrr7F161Zs2LABX3/9tcG57t274/nnnxf3AwICkJycjEWLFonHXn/9dYPXzJw5E/n5+fjiiy8QFxdn8prXr19HSkoKFi1ahClTpojH77vvPvHfa9asQUhICN58800AQP/+/aHRaDBx4kTMmTMHbm5u+Pnnn5GdnY3Dhw8bvPZWn3zyCaZPn474+HgAQL9+/bBr1y688847YnyrVq3C+fPnUVRUBCcnJ/G93iojIwNTpkwRv7gzMzORn5+P7OxsaDQasdzZs2cxY8YM5OfnG7Qu3OrAgQPIyMjAnj174OPjY3Buy5YtcHJyglarFUe1aLVaDB48GCdOnED//v3x7LPPGrymX79++Omnn7Bx40bMmDEDANClSxdkZ2cDAH788UdcvHjRKI5du3ahb9++mDlzJgAgMDAQ06dPx9///nexjKenp8Fr1qxZg86dOxslJ46Oju3SWnIrm2k54VBiIqKWO3fuHKZNm4ZPPvmkRXO9/PLLL9i4cSNGjBhx23J1dXVGX2q32rdvH86ePYtOnTph8ODB8PHxgUqlwpEjR8Qyer3eaA4MV1dXXL16FcXFxQCAzZs3o1+/ftiyZQsCAwPRt29fTJ06FefPn79jPbt378a1a9cAAJs2bUJkZCTUajUUCgWCg4OxYMEC8Q/dhoYGFBcXIzY21qCe2NhYFBUVifs3btxAQkICZs+ejQceeMDke79y5QomTpyIrKwsk1/mer0ezs7OBsNtXV1dAdxs6WjOnT5zU6KionDmzBnk5eVBEAScO3cO69evbzapAoCVK1fi6aefRpcuXQyOHz9+HL6+vggMDMTTTz+N//73v2bFcjdsJjkhIqKWEQQBkydPRlJSEsLDw29bduLEiejcuTN69+6Nbt26YcWKFc2WXb9+Pfbs2YNnnnmm2TJNX1xz587FG2+8gS1btqB79+4YMWKEmFiMGTMGRUVF+Oyzz9DY2IizZ89i/vz5AIDKykqxntOnT2PdunVYvXo1cnJyUFxcjD/+8Y/itcaMGYMVK1aguLgYgiBg7969WLVqFa5du4aamhqxnvXr16OxsRF5eXl44403sHjxYrz99tsAgJqaGjQ2NkKhUBi8D4VCgaqqKnH/nXfegaOjo9gSYcqLL76IqKgojBs3zuT5kSNHoqqqCosWLUJDQwMuXLggtk41ve/f++mnn/D5559j+vTpzV7XlKioKHz66aeIj4+Hs7MzvL294eHhgQ8++MBk+d27d+Pw4cMGj30AYOjQoVi9ejXy8/OxfPlyVFVVISoqyqDvSltgckJE1MF88MEH0Ol0SEtLu2PZ9957D/v27UNubi5OnjyJ1NRUk+W2b9+OyZMnY/ny5c22HAA3WxgAID09HRMmTIBSqcRHH30EmUyGdevWAbjZKrFo0SIkJSVBLpdjwIAB4l/0Dg4OYj16vR6rV6/G8OHDER0djZUrV+L777/H0aNHAQB/+9vfoFKp8NBDD8HJyQnjxo0T+2/cWk+vXr2wbNkyKJVKPP3000hPTxcfizT5fR8RQRDEY8XFxViyZAlycnKaXQZh06ZN2LZtW7OdZIGbfTc+/vhjLF68GJ07d4a3tzf69esHhUIhxnurI0eOYNy4cXjzzTcxevToZus1pbS0FDNnzsSbb76J4uJifPPNNygrKzPoU3SrlStXIjg42KBPD3DzMeOECRMwaNAgjBo1Cl999RUA4OOPPzYrHnMxOSEi6mC2bduGXbt2QS6Xw9HREf379wcAhIeHIzEx0aCst7c37r//fowbNw4ffvghsrOzjf6K37FjB+Li4pCRkYFJkybd9tpN/SyCgoLEY3K5HP369TPoYJqamoqLFy+ivLwcNTU1YmtDYGCgWI+joyMGDBggvmbgwIEAINbj6uqKVatW4cqVKzh16hTKy8vRt29fuLm5wcvLS6xnwIABBl/+AwcORFVVFRoaGuDl5QUHBweDVhIAqK6uFltTdu7cierqavj7+8PR0RGOjo44ffo0XnrpJfTt21f8zE+ePAkPDw+xDABMmDAB0dHRYr1//vOfUVVVhbNnz6K2thZz587Fr7/+Kr7vJqWlpRg5ciSmTZuGN95447afuSkajQbDhg3D7NmzERISgjFjxmDp0qVYtWqV0f29cuUK1qxZY9RqYkqXLl0waNAgHD9+3OyYzMHkhIiog3n//fdx4MABlJSUoKSkRBw2u3btWvFxhimCIAC42Teiyfbt2/H4449j4cKFeO655+54baVSCblcLrZuAMC1a9dw6tQpo46oMpkMvr6+cHV1xWeffQY/Pz8MGTIEADBs2DBcv34dJ0+eFMsfO3YMgHGHVicnJ/Tp0wcODg5Ys2YNnnjiCbFfx7Bhw3DixAmxRaepnqaVc52dnaFUKlFQUGBQZ0FBAaKiogAACQkJOHjwoPh5lpSUwNfXF7Nnz0Z+fj4A4LXXXjMqA9xsmfroo4+MPieFQoGuXbti7dq1cHFxMWgZOXLkCGJiYpCYmHjb+3U7V65cMZpKvilBa7rPTT7//HPo9Xr89a9/vWO9er0eP//8s1FnX0uzymid9957DytWrIAgCBg1ahSWLFnCFUOJiADU19fjxIkT4n5ZWRlKSkrg6ekJf3//FtXx+3Jdu3YFANxzzz3o06cPACAvLw/nzp1DREQEunbtitLSUrzyyisYNmyY2BrQlJikpKRgwoQJYuuCs7Oz2EFz9+7dmDRpEr777jux30pSUhLmzJkDPz8/BAQEiCOAbh0FsmjRIjz22GPo1KkTNm7ciIULF+Lzzz8Xv0BHjRqFIUOG4Nlnn0VmZiZu3LgBtVqN0aNHi60px44dw+7duzF06FBcuHABGRkZOHz4sMEjh+effx4ffPABUlJS8MILL+D48eNYsGCBQd+R1NRUJCQkIDw8HJGRkVi2bBnKy8vFRyA9evRAjx49DD5TJycneHt7iyOJvL29TXaC9ff3N2gVycrKQlRUFLp27YqCggLMnj0bCxcuhIeHB4D/JSaxsbFITU0VP3MHBwf07NlTrKe0tBQNDQ04f/48Ll26JCZDYWFhAIC4uDhMmzYN2dnZGDNmDCorKzFr1iw8+OCDRsOgV65ciaeeesroPQLAyy+/jLi4OPj7+6O6uhrz58+HTqczaoGzOKGdVVdXC/369RN+++034fr160JUVJRQVFTU4tfX1dUJAIS6uro2ie+1DQfFTdg003AjIpvw22+/CaWlpcJvv/1m7VDM9v333wsAjLbExESxzJw5c4SAgIAW11lWViYAEPbv3y8e27ZtmxAZGSm4u7sLLi4uwr333iu8+uqrwoULF8QyiYmJJmMZMWKEUbxlZWXisYaGBuGll14SevXqJbi5uQmjRo0SDh8+bBBTTEyMeO2hQ4cKeXl5RnGfPXtWGD9+vNC1a1dBoVAIkydPFmpra8XzpaWlQlhYmODq6ip069ZNGDdunPCf//zHqJ6ioiJh6NChglwuF/r16ye8/fbbwvXr1w3KaLVaISAgQHB2dhaGDBki7Nix47afaUBAgPDee+/dtgwA4YsvvjA4lpCQIHh6egrOzs5CSEiIsHr1aoPzc+bMMfmZ//5+BwQEmCx3q/fff18ICgoSXF1dBR8fH+Evf/mLcObMGYMyR48eFQAIW7duNfke4uPjBR8fH8HJyUnw9fUVxo8fLxw5cqTZ93y73z1zvr9lgvC79p029uuvvyI8PBwHDhyAi4sLHnnkEXz22We45557WvR6nU4Hd3d31NXVoVu3bhaPj/OcENm+q1evoqysTJzxs6Np6vSZk5Nj1TiIfu92v3vmfH+b3eekJVMr324a4J49e+Lll1+Gv78/fH19MWrUqBYnJkREdLOD6v/7f//P2mEQtRmzk5M7Ta18p2mAL1y4gC1btuDUqVM4e/YsioqKUFhY2Lp3QURkR8rKyuDn52ftMIjajNnJiUqlwvz5840WcWpy6zTAAwcORGZmJvz8/MQx5d9++y369+8PT09PuLq64vHHH8euXbuavZ5er4dOpzPYiIiIqOOy6FDilkwD7Ofnh6KiIly9ehWNjY3Yvn17s+smADfHaru7u4sb/1ogIiLq2CyanLRkGuCHHnoIY8eOxeDBgxESEoJ77rkHTz75ZLN1pqWloa6uTtwqKiosGTIRERFJTJvMc3K7aYAB4O23327xxDJyuRxyubzdViUmIiIi67Joy0lLpgG+W1yVmIiIyD5YNDlpyTTAd0ur1SIoKAgRERGtqoeIiIikzezHOneaWvlO0wDfLbVaDbVaLU7iQkRERB2T2S0ne/fuxeDBgzF48GAAN9ckGDx4MN58800AQHx8PDIzMzFv3jyEhYWhsLAQeXl5Rgs1mYstJ0RkEZtT2nczk0ajQUREBNzc3NCrVy889dRTBovomUuv1yMsLAwymUxcfwUAamtr8dhjj8HX1xdyuRx+fn6YMWOGwXQN27dvx7hx4+Dj44MuXbogLCwMn376aYuu+9VXX2Ho0KFwdXWFl5eX0fQT3333HaKiouDm5gYfHx+8+uqruH79ukEZQRDw7rvvYsCAAWKMCxYsMCij1WoxcOBAuLq64r777sPq1auNYrl48SLUajV8fHzg4uKCgQMHioshNrnd5KG/N336dMhkMmRmZpo8LwgCVCqVyYlK9+3bh9GjR8PDwwM9evTAc889h/r6evH8gQMHMHHiRPj5+cHV1RUDBw7EkiWGs5NfvXoVkydPxqBBg+Do6IinnnrKZByffvopQkND0blzZ/j4+OCZZ55BbW2teD46Ohoymcxoe/zxx8Uyc+fONTpvag0hSzO75SQ6OtpoRcPfS05ORnJy8l0HZQpbTojIHuzYsQNqtRoRERG4fv060tPTERsbi9LSUnTp0sXs+l555RX4+vriwIEDBsc7deqEcePGYf78+ejZsydOnDgBtVqN8+fP41//+hcAoKioCCEhIXj11VehUCjw1VdfYdKkSejWrRvi4uKaveaGDRswbdo0LFiwACNHjoQgCDh06H9Lgxw8eBBjx45Feno6Vq9ejbNnzyIpKQmNjY149913xXIpKSnYunUr3n33XQwaNAh1dXWoqakRz2dnZyMtLQ3Lly9HREQEdu/ejWnTpqF79+5ifA0NDRg9ejR69eqF9evXo0+fPqioqICbm5tYT9PkoUuXLsWwYcPw4YcfQqVSobS01GgRxdzcXPz73/82WjzvVpmZmSYXs/3ll18watQoxMfHIysrCzqdDrNmzcLkyZOxfv16AEBxcTF69uyJf/7zn+LUG8899xwcHBwwY8YMAEBjYyNcXV0xc+ZMbNiwwWQMP/zwAyZNmoT33nsPcXFx4mc8depUfPHFFwCAjRs3oqGhQXxNbW0tQkNDDRZoBIAHHngA3377rbjftDhjW7LKqsRERGTaN998Y7D/0UcfoVevXiguLsYjjzxiVl1ff/01tm7dig0bNuDrr782ONe9e3c8//zz4n5AQACSk5PFFYQB4PXXXzd4zcyZM5Gfn48vvvii2eTk+vXrSElJwaJFizBlyhTx+K3zWa1ZswYhISFii3v//v2h0WgwceJEzJkzB25ubvj555+RnZ2Nw4cPNzsX1ieffILp06cjPj4eANCvXz/s2rUL77zzjhjfqlWrcP78eRQVFcHJyUl8r7e6dfJQ4GZykZ+fj+zsbGg0GrHc2bNnMWPGDOTn5xu0LtzqwIEDyMjIwJ49e+Dj42NwbsuWLXBycoJWq0WnTjcfXGi1WgwePBgnTpxA//798eyzzxq8pl+/fvjpp5+wceNGMTnp0qWLOLHpjz/+iIsXLxrFsWvXLvTt21dcfTkwMBDTp0/H3//+d7FM08rSTdasWYPOnTsbJSeOjo7t0lpyK4t2iG1LfKxDRPaorq4OgPEXyZ2cO3cO06ZNwyeffILOnTvfsfwvv/yCjRs3YsSIEXeM53ax7Nu3D2fPnkWnTp0wePBg+Pj4QKVS4ciRI2IZvV5vtCicq6srrl69iuLiYgDA5s2b0a9fP2zZsgWBgYHo27cvpk6divPnz9+xnt27d+PatWsAgE2bNiEyMhJqtRoKhQLBwcFYsGCBOC1FSyYPBYAbN24gISEBs2fPxgMPPGDyvV+5cgUTJ05EVlaWyS9zvV4PZ2dnMTFpihe42dLRnDt95qZERUXhzJkzyMvLgyAIOHfuHNavX99sUgUAK1euxNNPP23UQnf8+HH4+voiMDAQTz/9NP773/+aFcvdsJnkhEOJicjeCIKA1NRUPPzwwwgODjbrdZMnT0ZSUhLCw8NvW3bixIno3LkzevfujW7dumHFihXNll2/fj327NmDZ555ptkyTV9cc+fOxRtvvIEtW7age/fuGDFihJhYjBkzBkVFRfjss8/Q2NiIs2fPYv78+QCAyspKsZ7Tp09j3bp1WL16NXJyclBcXIw//vGP4rXGjBmDFStWoLi4GIIgYO/evVi1ahWuXbsmPv7573//i/Xr16OxsRF5eXl44403sHjxYnGurZZMHgoA77zzDhwdHcWWCFNefPFFREVFYdy4cSbPjxw5ElVVVVi0aBEaGhpw4cIFsXWq6X3/3k8//YTPP/8c06dPb/a6pkRFReHTTz9FfHw8nJ2d4e3tDQ8PD3zwwQcmy+/evRuHDx8WW4+aDB06FKtXr0Z+fj6WL1+OqqoqREVFGfRdaQs2k5wQEdmbGTNm4ODBg/jss8/Met0HH3wAnU6HtLS0O5Z97733sG/fPuTm5uLkyZNITU01WW779u2YPHkyli9f3mzLAXCzhQEA0tPTMWHCBCiVSnz00UeQyWRYt24dgJutEosWLUJSUhLkcjkGDBgg/kXf1J/hxo0b0Ov1WL16NYYPH47o6GisXLkS33//vdhB+G9/+xtUKhUeeughODk5Ydy4cZg8ebJRPb169cKyZcugVCrx9NNPIz09XXws0uR2k4cWFxdjyZIlyMnJMdmXBLjZQrNt27ZmO8kCN/tufPzxx1i8eDE6d+4Mb29v9OvXDwqFwmQ/jiNHjmDcuHF48803MXr06GbrNaW0tBQzZ87Em2++ieLiYnzzzTcoKytrduTsypUrERwcjAcffNDguEqlwoQJEzBo0CCMGjUKX331FQDg448/Nisec9lMcsLHOkRkT1544QVs2rQJ33//Pfr06WPWa7dt24Zdu3ZBLpfD0dER/fv3BwCEh4cjMTHRoKy3tzfuv/9+jBs3Dh9++CGys7ON/orfsWMH4uLikJGRgUmTJt322k39LIKCgsRjcrkc/fr1E1enB26O9Lx48SLKy8tRU1MjtjYEBgaK9Tg6OmLAgAHiawYOHAgAYj2urq5YtWoVrly5glOnTqG8vBx9+/aFm5sbvLy8xHoGDBhg8OU/cOBAVFVVoaGhoUWTh+7cuRPV1dXw9/eHo6MjHB0dcfr0abz00kvo27ev+JmfPHkSHh4eYhkAmDBhAqKjo8V6//znP6Oqqgpnz55FbW0t5s6di19//VV8301KS0sxcuRITJs2DW+88cZtP3NTNBoNhg0bhtmzZyMkJARjxozB0qVLsWrVKqP7e+XKFaxZs8ao1cSULl26YNCgQTh+/LjZMZnDZpITPtYhInsgCAJmzJiBjRs3Ytu2bUZfWi3x/vvv48CBAygpKUFJSYk4bHbt2rW3XTqkaSSmXq8Xj23fvh2PP/44Fi5ciOeee+6O11YqlZDL5QbDn69du4ZTp04ZdUSVyWTw9fWFq6srPvvsM/j5+WHIkCEAgGHDhuH69es4efKkWP7YsWMAjDu0Ojk5oU+fPnBwcMCaNWvwxBNPiP06hg0bhhMnTogtOk31+Pj4wNnZuUWThyYkJODgwYPi51lSUgJfX1/Mnj0b+fn5AIDXXnvNqAxws2Xqo48+MvqcFAoFunbtirVr18LFxcWgZeTIkSOIiYlBYmJii5d6+b0rV64Y9G0B/tea9PsRt59//jn0ej3++te/3rFevV6Pn3/+2aizr6VxtA4RkYSo1Wr861//wpdffgk3NzfxL3p3d3ex8+Sd/H74a9euXQEA99xzj9gKk5eXh3PnziEiIgJdu3ZFaWkpXnnlFQwbNkxsDWhKTFJSUjBhwgQxFmdnZ7GD5u7duzFp0iR89913Yr+VpKQkzJkzB35+fggICBBHAN06CmTRokV47LHH0KlTJ2zcuBELFy7E559/Ln6Bjho1CkOGDMGzzz6LzMxM3LhxA2q1GqNHjxZbU44dO4bdu3dj6NChuHDhAjIyMnD48GGDRw7PP/88PvjgA6SkpOCFF17A8ePHsWDBAoO+I3eaPLRHjx7o0aOHwWfq5OQEb29vcSSRt7e3yU6w/v7+BglmVlYWoqKi0LVrVxQUFGD27NlYuHAhPDw8APwvMYmNjUVqaqr4mTs4OKBnz55iPaWlpWhoaMD58+dx6dIlMRkKCwsDAMTFxWHatGnIzs7GmDFjUFlZiVmzZuHBBx80Gga9cuVKPPXUU0bvEQBefvllxMXFwd/fH9XV1Zg/fz50Op1RC5ylMTkhIpKQpr4Qtz4KAG4OKW7qTzF37lzk5OTg1KlTd30dV1dXLF++HC+++CL0ej38/Pwwfvx4vPbaa2KZnJwcXLlyBRqNxmBI7YgRI7B9+3YAN/9CP3r0qDg6BriZeDg6OiIhIQG//fYbhg4dim3btqF79+5ima+//hpvv/029Ho9QkND8eWXX0KlUonnO3XqhM2bN+OFF17AI488gi5dukClUmHx4sVimcbGRixevBhHjx6Fk5MTYmJiUFRUJCZXAODn54etW7fixRdfREhICHr37o2UlBS8+uqrYpn4+HjU1tZi3rx5qKysRHBwsEUmDzVl9+7dmDNnDurr63H//ffjww8/REJCgnh+3bp1+PXXX/Hpp58aTHgXEBBgcL/Hjh2L06dPi/tNE6M2tYpMnjwZly5dQlZWFl566SV4eHhg5MiReOeddwziOXbsGH744Qds3brVZLxnzpzBxIkTUVNTg549e+Khhx7Crl272uSzuZVMuNOMahJx66rEx44dQ11dHbp162bx66Rt/N9EQRqn3/Vaj1sCIpK+q1evoqysTJzxs6NpSlJycnKsGgfR793ud69pEtWWfH/bTMsJZ4glIrppx44dKCwstHYYRG3GZpITIiK6qayszNohELUpmxmtQ0RERPaByQkRERFJis0kJ5yEjYiIyD7YTHLCSdiIyFw2MhiRqMOw1O+czSQnREQt5eTkBODmHBxE1H6afueafgfvFkfrEFGH4+DgAA8PD1RXVwMAOnfu3OyCbUTUeoIg4MqVK6iuroaHh4fJhQzNweSEiDqkpqnEmxIUImp7Hh4eJqfxNxeTEyLqkGQyGXx8fNCrVy+DqdWJqG04OTm1usWkCZMTIurQHBwcLPYfJhG1D3aIJSIiIkmxmeSE85wQERHZB5tJTjjPCRERkX2wmeSEiIiI7AOTEyIiIpIUJidEREQkKUxOiIiISFLaPTk5evQowsLCxM3V1RW5ubntHQYRERFJVLtPwnbfffehpKQEAFBfX4++ffti9OjR7R0GERERSZRVH+ts2rQJjz76KLp06WLNMIiIiEhCzE5OCgsLERcXB19fX8hkMpOPZJYuXYrAwEC4uLhAqVRi586dJuv6/PPPER8fb3bQRERE1HGZnZxcvnwZoaGhyMrKMnl+7dq1mDVrFtLT07F//34MHz4cKpUK5eXlBuV0Oh1+/PFHjB079u4iJyIiog7J7D4nKpUKKpWq2fMZGRmYMmUKpk6dCgDIzMxEfn4+srOzodFoxHJffvklxowZAxcXl9teT6/XQ6/Xi/s6nc7ckImIiMiGWLTPSUNDA4qLixEbG2twPDY2FkVFRQbHWvpIR6PRwN3dXdz8/PwsGTIRERFJjEWTk5qaGjQ2NkKhUBgcVygUqKqqEvfr6uqwe/dujBkz5o51pqWloa6uTtwqKiosGTIRERFJTJsMJZbJZAb7giAYHHN3d8e5c+daVJdcLodcLodWq4VWq0VjY6NFYyUiIiJpsWjLiZeXFxwcHAxaSQCgurraqDXFXFyVmIiIyD5YNDlxdnaGUqlEQUGBwfGCggJERUW1qm6tVougoCBERES0qh4iIiKSNrMf69TX1+PEiRPifllZGUpKSuDp6Ql/f3+kpqYiISEB4eHhiIyMxLJly1BeXo6kpKRWBapWq6FWq6HT6eDu7t6quoiIiEi6zE5O9u7di5iYGHE/NTUVAJCYmIicnBzEx8ejtrYW8+bNQ2VlJYKDg5GXl4eAgIBWBco+J0RERPZBJgiCYO0gzNHUclJXV4du3bpZvP60jYfEf2ucVhiejFti8esRERHZA3O+v626tg4RERHR79lMcsIOsURERPbBZpITDiUmIiKyDzaTnBAREZF9sJnkhI91iIiI7IPNJCd8rENERGQfbCY5ISIiIvvA5ISIiIgkxWaSE/Y5ISIisg82k5ywzwkREZF9sJnkhIiIiOwDkxMiIiKSFCYnREREJCk2k5ywQywREZF9sJnkhB1iiYiI7IPNJCdERERkH5icEBERkaQwOSEiIiJJYXJCREREksLkhIiIiCTFZpITDiUmIiKyDzaTnHAoMRERkX2wmeSEiIiI7AOTEyIiIpIUJidEREQkKUxOiIiISFKskpyUlZUhJiYGQUFBGDRoEC5fvmyNMIiIiEiCHK1x0cmTJ2P+/PkYPnw4zp8/D7lcbo0wiIiISILaPTk5cuQInJycMHz4cACAp6dne4dAREREEmb2Y53CwkLExcXB19cXMpkMubm5RmWWLl2KwMBAuLi4QKlUYufOneK548ePo2vXrnjyyScxZMgQLFiwoFVvgIiIiDoWs5OTy5cvIzQ0FFlZWSbPr127FrNmzUJ6ejr279+P4cOHQ6VSoby8HABw7do17Ny5E1qtFj/99BMKCgpQUFDQundBREREHYbZyYlKpcL8+fMxfvx4k+czMjIwZcoUTJ06FQMHDkRmZib8/PyQnZ0NAOjTpw8iIiLg5+cHuVyOsWPHoqSkpNnr6fV66HQ6g42IiIg6LouO1mloaEBxcTFiY2MNjsfGxqKoqAgAEBERgXPnzuHChQu4ceMGCgsLMXDgwGbr1Gg0cHd3Fzc/Pz9LhkxEREQSY9HkpKamBo2NjVAoFAbHFQoFqqqqAACOjo5YsGABHnnkEYSEhODee+/FE0880WydaWlpqKurE7eKigpLhkxEREQS0yajdWQymcG+IAgGx1QqFVQqVYvqksvlkMvl0Gq10Gq1aGxstGisREREJC0WbTnx8vKCg4OD2ErSpLq62qg1xVxclZiIiMg+WDQ5cXZ2hlKpNBp9U1BQgKioqFbVrdVqERQUhIiIiFbVQ0RERNJm9mOd+vp6nDhxQtwvKytDSUkJPD094e/vj9TUVCQkJCA8PByRkZFYtmwZysvLkZSU1KpA1Wo11Go1dDod3N3dW1UXERERSZfZycnevXsRExMj7qempgIAEhMTkZOTg/j4eNTW1mLevHmorKxEcHAw8vLyEBAQYLmoiYiIqMMyOzmJjo6GIAi3LZOcnIzk5OS7DsoUdoglIiKyD1ZZlfhusEMsERGRfbCZ5IQdYomIiOxDu69KfLek1CE2beMh8d+a8YOsGAkREVHHYzMtJ0RERGQfbCY54WMdIiIi+2AzyQk7xBIREdkHm0lOiIiIyD4wOSEiIiJJsZnkhH1OiIiI7AOHErfW5hTD/bgl1omDiIiog7CZlhMiIiKyD0xOiIiISFKYnBAREZGk2Exywg6xRERE9sFmkhNOwkZERGQfbCY5ISIiIvvA5ISIiIgkhckJERERSQqTEyIiIpIUm5kh1qZxFlkiIqIWY8sJERERSYrNJCec54SIiMg+2MxjHcku/NeMtI2HxH9rnKwYCBERkY2xmZYTIiIisg9MToiIiEhSmJwQERGRpDA5ISIiIkmxSnLi6OiIsLAwhIWFYerUqdYIgYiIiCTKKqN1PDw8UFJSYo1LExERkcTxsQ4RERFJitnJSWFhIeLi4uDr6wuZTIbc3FyjMkuXLkVgYCBcXFygVCqxc+dOg/M6nQ5KpRIPP/wwduzYcdfBExERUcdjdnJy+fJlhIaGIisry+T5tWvXYtasWUhPT8f+/fsxfPhwqFQqlJeXi2VOnTqF4uJi/OMf/8CkSZOg0+nu/h0QERFRh2J2cqJSqTB//nyMHz/e5PmMjAxMmTIFU6dOxcCBA5GZmQk/Pz9kZ2eLZXx9fQEAwcHBCAoKwrFjx5q9nl6vh06nM9iIiIio47Jon5OGhgYUFxcjNjbW4HhsbCyKiooAABcuXIBerwcAnDlzBqWlpejXr1+zdWo0Gri7u4ubn5+fJUMmIiIiibHoaJ2amho0NjZCoVAYHFcoFKiqqgIA/Pzzz5g+fTo6deoEmUyGJUuWwNPTs9k609LSkJqaKu7rdLoOk6AYrL8zfpAVIyEiIpKONhlKLJPJDPYFQRCPRUVF4dChQ6ZeZpJcLodcLodWq4VWq0VjY6NFY5WMzSmG+3FLrBMHERGRlVn0sY6XlxccHBzEVpIm1dXVRq0p5lKr1SgtLcWePXtaVQ8RERFJm0WTE2dnZyiVShQUFBgcLygoQFRUVKvq1mq1CAoKQkRERKvqISIiImkz+7FOfX09Tpw4Ie6XlZWhpKQEnp6e8Pf3R2pqKhISEhAeHo7IyEgsW7YM5eXlSEpKalWgarUaarUaOp0O7u7uraqLiIiIpMvs5GTv3r2IiYkR95s6qyYmJiInJwfx8fGora3FvHnzUFlZieDgYOTl5SEgIKBVgXb4PidEREQE4C6Sk+joaAiCcNsyycnJSE5OvuugTGHLCRERkX3g2jpEREQkKTaTnLBDLBERkX1ok3lO2oK9PtbhRG1ERGRvbKblhIiIiOyDzbSccLQOOIssERHZBZtpOeEMsURERPbBZpITIiIisg9MToiIiEhS2OekI7q1bwr7pRARkY2xmZYT9jkhIiKyDzbTckLNu3UuFADQOFkpECIiIguwmZYTIiIisg9MToiIiEhSbCY54do6RERE9sFmkhN2iCUiIrIPNpOcEBERkX3gaB07YTSihyscExGRRDE5sVecqI2IiCSKj3WIiIhIUpicEBERkaTYzGMdrq3T9tgvhYiIpMBmWk44lJiIiMg+2EzLCVlBSzrNsmMtERFZmM20nBAREZF9YMsJmYUrIBMRUVtjywkRERFJitWSkytXriAgIAAvv/yytUIgIiIiCbJacvL2229j6NCh1ro8ERERSZRVkpPjx4/jP//5D8aOHWuNyxMREZGEmZ2cFBYWIi4uDr6+vpDJZMjNzTUqs3TpUgQGBsLFxQVKpRI7d+40OP/yyy9Do9HcddAkbWkbDxlsRERE5jA7Obl8+TJCQ0ORlZVl8vzatWsxa9YspKenY//+/Rg+fDhUKhXKy8sBAF9++SUGDBiAAQMGtC5ysh2bU/63ERER3YHZQ4lVKhVUKlWz5zMyMjBlyhRMnToVAJCZmYn8/HxkZ2dDo9Fg165dWLNmDdatW4f6+npcu3YN3bp1w5tvvmmyPr1eD71eL+7rdDpzQyYiIiIbYtE+Jw0NDSguLkZsbKzB8djYWBQVFQEANBoNKioqcOrUKbz77ruYNm1as4lJU3l3d3dx8/Pzs2TIREREJDEWnYStpqYGjY2NUCgUBscVCgWqqqruqs60tDSkpqaK+zqdjglKB8BFBomIqDltMkOsTCYz2BcEwegYAEyePPmOdcnlcsjlcq5K3NFxjR4iIvo/Fn2s4+XlBQcHB6NWkurqaqPWFHNxVWIiIiL7YNHkxNnZGUqlEgUFBQbHCwoKEBUV1aq6tVotgoKCEBER0ap6iIiISNrMfqxTX1+PEydOiPtlZWUoKSmBp6cn/P39kZqaioSEBISHhyMyMhLLli1DeXk5kpKSWhWoWq2GWq2GTqeDu7t7q+oi22C8yOCK/+3w0Q8RUYdldnKyd+9exMTEiPtNnVUTExORk5OD+Ph41NbWYt68eaisrERwcDDy8vIQEBBguaiJiIiowzI7OYmOjoYgCLctk5ycjOTk5LsOyhR2iCUiIrIPVlv4z1zsEEtERGQfbCY5YYdYIiIi+9Am85y0BXaIJVNu22kWYMdZIiIbZDMtJ0RERGQfmJwQERGRpNhMcsI+J0RERPaBfU6ow+Mig0REtsVmkhMii7l1kUGAnWaJiCSGj3WIiIhIUmwmOeEkbERERPbBZpITIiIisg9MToiIiEhS2CGWCJxplohISthyQkRERJJiM8kJR+sQERHZB5tJTjhah4iIyD7YTHJCRERE9oEdYolaiNPgExG1DyYnRHeL0+ATEbUJPtYhIiIiSWFyQkRERJJiM8kJhxITERHZB5tJTjiUmIiIyD6wQyyRhd06qofT4BMRmc9mWk6IiIjIPjA5ISIiIklhckJERESS0u7JyaVLlxAREYGwsDAMGjQIy5cvb+8QiIiISMLavUNs586dsWPHDnTu3BlXrlxBcHAwxo8fjx49erR3KERERCRB7d5y4uDggM6dOwMArl69isbGRgiC0N5hEBERkUSZnZwUFhYiLi4Ovr6+kMlkyM3NNSqzdOlSBAYGwsXFBUqlEjt37jQ4f/HiRYSGhqJPnz545ZVX4OXldddvgIiIiDoWs5OTy5cvIzQ0FFlZWSbPr127FrNmzUJ6ejr279+P4cOHQ6VSoby8XCzj4eGBAwcOoKysDP/6179w7ty5u38HRERE1KGY3edEpVJBpVI1ez4jIwNTpkzB1KlTAQCZmZnIz89HdnY2NBqNQVmFQoGQkBAUFhbiT3/6k8n69Ho99Hq9uK/T6cwNmUhyDCZqGz/IipEQEUmPRfucNDQ0oLi4GLGxsQbHY2NjUVRUBAA4d+6cmGDodDoUFhbivvvua7ZOjUYDd3d3cfPz87NkyERERCQxFk1Oampq0NjYCIVCYXBcoVCgqqoKAHDmzBk88sgjCA0NxcMPP4wZM2YgJCSk2TrT0tJQV1cnbhUVFZYMmYiIiCSmTYYSy2Qyg31BEMRjSqUSJSUlLa5LLpdDLpdDq9VCq9WisbHRkqESWd/mFMN9rr9DRHbOoi0nXl5ecHBwEFtJmlRXVxu1ppiLqxITERHZB4smJ87OzlAqlSgoKDA4XlBQgKioqFbVrdVqERQUhIiIiFbVQ0RERNJm9mOd+vp6nDhxQtwvKytDSUkJPD094e/vj9TUVCQkJCA8PByRkZFYtmwZysvLkZSU1KpA1Wo11Go1dDod3N3dW1UXERERSZfZycnevXsRExMj7qempgIAEhMTkZOTg/j4eNTW1mLevHmorKxEcHAw8vLyEBAQ0KpA2eeE7Br7pRCRHTE7OYmOjr7jdPPJyclITk6+66BMYcsJ2RuDuVCcrBgIEVE7a/e1dYiIiIhux2aSE3aIJSIisg82k5xwKDEREZF9sJnkhIiIiOyDzSQnfKxDRERkH2wmOeFjHSIiIvtgM8kJERER2QcmJ0RERCQpNpOcsM8JERGRfTB7hlhr4QyxRMYMZpEdP8iKkRARWY7NJCdEdAdcf4eIOgibeaxDRERE9oHJCREREUmKzSQn7BBLRERkH2wmOeEkbERERPbBZpITIiIisg9MToiIiEhSmJwQERGRpDA5ISIiIklhckJERESSYjMzxGq1Wmi1WjQ2Nlo7FCKbwinuicjW2ExywrV1iCyAU9wTkQ3gYx0iIiKSFCYnREREJClMToiIiEhSmJwQERGRpDA5ISIiIklp9+SkoqIC0dHRCAoKQkhICNatW9feIRAREZGEtftQYkdHR2RmZiIsLAzV1dUYMmQIxo4diy5durR3KERERCRB7Z6c+Pj4wMfHBwDQq1cveHp64vz580xOiIiICMBdPNYpLCxEXFwcfH19IZPJkJuba1Rm6dKlCAwMhIuLC5RKJXbu3Gmyrr179+LGjRvw8/MzO3AiaiObUww3IqJ2ZnZycvnyZYSGhiIrK8vk+bVr12LWrFlIT0/H/v37MXz4cKhUKpSXlxuUq62txaRJk7Bs2bLbXk+v10On0xlsRERE1HGZ/VhHpVJBpVI1ez4jIwNTpkzB1KlTAQCZmZnIz89HdnY2NBoNgJsJxx/+8AekpaUhKirqttfTaDR46623zA2TiMxgsP6OkxUDISKChUfrNDQ0oLi4GLGxsQbHY2NjUVRUBAAQBAGTJ0/GyJEjkZCQcMc609LSUFdXJ24VFRWWDJmIiIgkxqLJSU1NDRobG6FQKAyOKxQKVFVVAQB+/PFHrF27Frm5uQgLC0NYWBgOHTpkqjoAgFwuR7du3fDJJ5/goYcewqOPPmrJkImIiEhi2mS0jkwmM9gXBEE89vDDD+PGjRtm18lViYmIiOyDRVtOvLy84ODgILaSNKmurjZqTTGXVqtFUFAQIiIiWlUPERERSZtFkxNnZ2colUoUFBQYHC8oKLhjx9c7UavVKC0txZ49e1pVDxEREUmb2Y916uvrceLECXG/rKwMJSUl8PT0hL+/P1JTU5GQkIDw8HBERkZi2bJlKC8vR1JSkkUDJyIioo7J7ORk7969iImJEfdTU1MBAImJicjJyUF8fDxqa2sxb948VFZWIjg4GHl5eQgICGhVoFqtFlqtFo2Nja2qh4gs4PeTs8UtsU4cRNQhmZ2cREdHQxCE25ZJTk5GcnLyXQdlCjvEEhER2Yd2X1vnbrHlhMi6OFEbEbUXi3aIbUvsEEtERGQfbCY5ISIiIvvA5ISIiIgkxWaSE07CRkREZB9sJjlhnxMiIiL7YDPJCREREdkHJidEREQkKTaTnLDPCRERkX2wmeSEfU6IiIjsg80kJ0RERGQfbGb6eiKSvlunuAcAzfhBVoqEiGwZkxMiojtg0kXUvmwmOeHCf0REts1g8UgmeFYl9YTbZpITtVoNtVoNnU4Hd3d3a4dDRETUcWxOMdyPW2KdOP4PO8QSERGRpDA5ISIiIkmxmcc6RGSDbm0qtnIzsUV11PdFJBFsOSEiIiJJYXJCREREkmIzj3U4lJiIqAOR2OgQkhabaTnh2jpERET2wWaSEyIiIrIPTE6IiIhIUmymzwkRUbvgMGEiq2NyQkRElsXOrtRKVklO/vCHP2D79u149NFHsX79emuEQERkX5gwkA2xSp+TmTNnYvXq1da4NBHZos0p/9uIqMOzSnISExMDNzc3a1yaiIiIJM7s5KSwsBBxcXHw9fWFTCZDbm6uUZmlS5ciMDAQLi4uUCqV2LlzpyViJSIiInPZYMuj2cnJ5cuXERoaiqysLJPn165di1mzZiE9PR379+/H8OHDoVKpUF5e3upgiYiIqOMzu0OsSqWCSqVq9nxGRgamTJmCqVOnAgAyMzORn5+P7OxsaDQaswPU6/XQ6/Xivk6nM7sOIiIish0WHa3T0NCA4uJivPbaawbHY2NjUVRUdFd1ajQavPXWW5YIj4jISNrGQwb7GicrBUJEIot2iK2pqUFjYyMUCoXBcYVCgaqqKnF/zJgx+NOf/oS8vDz06dPntuvlpKWloa6uTtwqKiosGTIRERFJTJvMcyKTyQz2BUEwOJafn9/iuuRyOeRyOVclJiKSsFtboNj61AJtOBNxR2gNtGjLiZeXFxwcHAxaSQCgurraqDXFXFyVmIiIyD5YNDlxdnaGUqlEQUGBwfGCggJERUW1qm6tVougoCBERES0qh4iIiKSNrMf69TX1+PEiRPifllZGUpKSuDp6Ql/f3+kpqYiISEB4eHhiIyMxLJly1BeXo6kpKRWBapWq6FWq6HT6eDu7t6quoiIiEi6zE5O9u7di5iYGHE/NTUVAJCYmIicnBzEx8ejtrYW8+bNQ2VlJYKDg5GXl4eAgIBWBco+J0RERPbB7OQkOjoagiDctkxycjKSk5PvOihT2HJCRERkH6yytg4RERFRc2wmOWGHWCIiIvtgM8kJhxITERHZB5tJToiIiMg+tMkMsW2Bo3WI6FZGs2COH2SlSKjN3DqLKmDxmVRJumym5YSPdYiIiOyDzSQnREREZB+YnBAREZGksM8JEbUr9hWxPQYrDvN+WZW9/P7YTMsJ+5wQERHZB5tJToiIiMg+MDkhIiIiSWFyQkRERJLCDrFEJDlGnf6crBSIRNlLp0iyXzbTcsIOsURERPbBZpITIiIisg9MToiIiEhSmJwQERGRpDA5ISIiIklhckJERESSwqHERETthEOA78xgHZ8ONIScw+PNYzMtJxxKTEREZB9sJjkhIiIi+8DkhIiIiCSFyQkRERFJCpMTIiIikhQmJ0RERCQpVklOtmzZgvvuuw/33nsvVqxYYY0QiIiISKLafZ6T69evIzU1Fd9//z26deuGIUOGYPz48fD09GzvUIiIiEiC2r3lZPfu3XjggQfQu3dvuLm5YezYscjPz2/vMIiIiEiizE5OCgsLERcXB19fX8hkMuTm5hqVWbp0KQIDA+Hi4gKlUomdO3eK53755Rf07t1b3O/Tpw/Onj17d9ETERFRh2N2cnL58mWEhoYiKyvL5Pm1a9di1qxZSE9Px/79+zF8+HCoVCqUl5cDAARBMHqNTCZr9np6vR46nc5gIyIioo7L7D4nKpUKKpWq2fMZGRmYMmUKpk6dCgDIzMxEfn4+srOzodFo0Lt3b4OWkjNnzmDo0KHN1qfRaPDWW2+ZGyYREbWAwVo2NrLWj6VibslaR5YqQ+axaJ+ThoYGFBcXIzY21uB4bGwsioqKAAAPPvggDh8+jLNnz+LSpUvIy8vDmDFjmq0zLS0NdXV14lZRUWHJkImIiEhiLDpap6amBo2NjVAoFAbHFQoFqqqqbl7Q0RGLFy9GTEwMbty4gVdeeQU9evRotk65XA65XM5ViYmIiOxEmwwl/n0fEkEQDI49+eSTePLJJ82qU61WQ61WQ6fTwd3d3SJxEhERkfRY9LGOl5cXHBwcxFaSJtXV1UatKebSarUICgpCREREq+ohIiIiabNocuLs7AylUomCggKD4wUFBYiKimpV3Wq1GqWlpdizZ0+r6iEiIiJpM/uxTn19PU6cOCHul5WVoaSkBJ6envD390dqaioSEhIQHh6OyMhILFu2DOXl5UhKSrJo4ERERNQxmZ2c7N27FzExMeJ+amoqACAxMRE5OTmIj49HbW0t5s2bh8rKSgQHByMvLw8BAQGtCpQdYomIiOyD2clJdHS0yYnUbpWcnIzk5OS7DsoUdoglIiKyD1ZZlfhusEMsERGRfbCZ5IQdYomIiOyDzSQnREREZB+YnBAREZGk2Exywj4nRERE9sFmkhP2OSEiIrIPNpOcEBERkX1ok4X/2lLTHCs6na5N6tdfqRf/rXPSG578v2tKuYxROXsu83/l2rPM72Oy5zJG5ZorY+J32VL1tISlfhbv6lqWirkl9VxpwX1tph6zy7Tnz2Irvgvu6mfRUmXa8GfIUv93WlpTrHeaKw0AZEJLSklA0wyxDQ0NOHnypLXDISIiortQUVGBPn363LaMzSQnTW7cuIFffvkFbm5ukMlkFqtXp9PBz88PFRUV6Natm8XqJcvg/ZE+3iPp4z2Sto5+fwRBwKVLl+Dr64tOnW7fq8TmHut06tTpjhlXa3Tr1q1D/lB0FLw/0sd7JH28R9LWke9PS5efYYdYIiIikhQmJ0RERCQpTE7+j1wux5w5cyCXy60dCpnA+yN9vEfSx3skbbw//2NzHWKJiIioY2PLCREREUkKkxMiIiKSFCYnREREJClMToiIiEhSmJwAWLp0KQIDA+Hi4gKlUomdO3daOyS7UVhYiLi4OPj6+kImkyE3N9fgvCAImDt3Lnx9feHq6oro6GgcOXLEoIxer8cLL7wALy8vdOnSBU8++STOnDnTju+i49JoNIiIiICbmxt69eqFp556CkePHjUow3tkPdnZ2QgJCREn7YqMjMTXX38tnue9kR6NRgOZTIZZs2aJx3ifTBDs3Jo1awQnJydh+fLlQmlpqZCSkiJ06dJFOH36tLVDswt5eXlCenq6sGHDBgGA8MUXXxicX7hwoeDm5iZs2LBBOHTokBAfHy/4+PgIOp1OLJOUlCT07t1bKCgoEPbt2yfExMQIoaGhwvXr19v53XQ8Y8aMET766CPh8OHDQklJifD4448L/v7+Qn19vViG98h6Nm3aJHz11VfC0aNHhaNHjwqvv/664OTkJBw+fFgQBN4bqdm9e7fQt29fISQkREhJSRGP8z4Zs/vk5MEHHxSSkpIMjt1///3Ca6+9ZqWI7Nfvk5MbN24I3t7ewsKFC8VjV69eFdzd3YV//OMfgiAIwsWLFwUnJydhzZo1YpmzZ88KnTp1Er755pt2i91eVFdXCwCEHTt2CILAeyRF3bt3F1asWMF7IzGXLl0S7r33XqGgoEAYMWKEmJzwPplm1491GhoaUFxcjNjYWIPjsbGxKCoqslJU1KSsrAxVVVUG90cul2PEiBHi/SkuLsa1a9cMyvj6+iI4OJj3sA3U1dUBADw9PQHwHklJY2Mj1qxZg8uXLyMyMpL3RmLUajUef/xxjBo1yuA475NpNrfwnyXV1NSgsbERCoXC4LhCoUBVVZWVoqImTffA1P05ffq0WMbZ2Rndu3c3KsN7aFmCICA1NRUPP/wwgoODAfAeScGhQ4cQGRmJq1evomvXrvjiiy8QFBQkfmnx3ljfmjVrsG/fPuzZs8foHH+HTLPr5KSJTCYz2BcEwegYWc/d3B/eQ8ubMWMGDh48iB9++MHoHO+R9dx3330oKSnBxYsXsWHDBiQmJmLHjh3ied4b66qoqEBKSgq2bt0KFxeXZsvxPhmy68c6Xl5ecHBwMMo8q6urjbJYan/e3t4AcNv74+3tjYaGBly4cKHZMtR6L7zwAjZt2oTvv/8effr0EY/zHlmfs7Mz+vfvj/DwcGg0GoSGhmLJkiW8NxJRXFyM6upqKJVKODo6wtHRETt27MD7778PR0dH8XPmfTJk18mJs7MzlEolCgoKDI4XFBQgKirKSlFRk8DAQHh7exvcn4aGBuzYsUO8P0qlEk5OTgZlKisrcfjwYd5DCxAEATNmzMDGjRuxbds2BAYGGpznPZIeQRCg1+t5byTi0UcfxaFDh1BSUiJu4eHh+Mtf/oKSkhL069eP98kU6/TDlY6mocQrV64USktLhVmzZgldunQRTp06Ze3Q7MKlS5eE/fv3C/v37xcACBkZGcL+/fvFodwLFy4U3N3dhY0bNwqHDh0SJk6caHKIXZ8+fYRvv/1W2LdvnzBy5MgOPcSuPT3//POCu7u7sH37dqGyslLcrly5IpbhPbKetLQ0obCwUCgrKxMOHjwovP7660KnTp2ErVu3CoLAeyNVt47WEQTeJ1PsPjkRBEHQarVCQECA4OzsLAwZMkQcJklt7/vvvxcAGG2JiYmCINwcZjdnzhzB29tbkMvlwiOPPCIcOnTIoI7ffvtNmDFjhuDp6Sm4uroKTzzxhFBeXm6Fd9PxmLo3AISPPvpILMN7ZD3PPvus+H9Xz549hUcffVRMTASB90aqfp+c8D4ZkwmCIFinzYaIiIjImF33OSEiIiLpYXJCREREksLkhIiIiCSFyQkRERFJCpMTIiIikhQmJ0RERCQpTE6IiIhIUpicEBERkaQwOSEiIiJJYXJCREREksLkhIiIiCSFyQkRERFJyv8H/bHef9BJygcAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "{'whiskers': [<matplotlib.lines.Line2D at 0x7fcfad78cb60>,\n",
       "  <matplotlib.lines.Line2D at 0x7fcfadd20440>,\n",
       "  <matplotlib.lines.Line2D at 0x7fcfadd21dc0>,\n",
       "  <matplotlib.lines.Line2D at 0x7fcfadd221b0>],\n",
       " 'caps': [<matplotlib.lines.Line2D at 0x7fcfadd20860>,\n",
       "  <matplotlib.lines.Line2D at 0x7fcfadd20fb0>,\n",
       "  <matplotlib.lines.Line2D at 0x7fcfadd22750>,\n",
       "  <matplotlib.lines.Line2D at 0x7fcfadd22a80>],\n",
       " 'boxes': [<matplotlib.lines.Line2D at 0x7fcf4885e7b0>,\n",
       "  <matplotlib.lines.Line2D at 0x7fcfadd21ac0>],\n",
       " 'medians': [<matplotlib.lines.Line2D at 0x7fcfadd21610>,\n",
       "  <matplotlib.lines.Line2D at 0x7fcfadd22ba0>],\n",
       " 'fliers': [<matplotlib.lines.Line2D at 0x7fcfadd21700>,\n",
       "  <matplotlib.lines.Line2D at 0x7fcfadd236b0>],\n",
       " 'means': []}"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAigAAAGdCAYAAAA44ojeAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQAALhhJREFUeJzt3X9U1nWe///HxY+LgAAFk0uEkkG3doOaso5FOuYU1Bx/cOX4qXMsT+2xs02pE+qAWc2ncj1SUDHNGv2wOdNOHcf9jCG4pg3sbJJm7VHXH+CeaZTVFQRiUgRE5BJ4f//wy3u5AE0UvV7Xxf12znWG6/16Mjw58b6uh6/363q9HZZlWQIAADBIkK8bAAAA6IuAAgAAjENAAQAAxiGgAAAA4xBQAACAcQgoAADAOAQUAABgHAIKAAAwToivG7gU3d3dqqurU1RUlBwOh6/bAQAAF8GyLLW2tiohIUFBQReeI/HLgFJXV6ekpCRftwEAAC5BTU2NEhMTL1jjlwElKipK0rlfMDo62sfdAACAi9HS0qKkpCT7ffxC/DKg9FzWiY6OJqAAAOBnLmZ5BotkAQCAcQgoAADAOAQUAABgHAIKAAAwDgEFAAAYh4ACAACMQ0ABAADGIaAAAADj+OVGbRg+urq6tG3bNtXX12vMmDGaMmWKgoODfd0WgCHA+Y0LYQYFxiouLtb48eM1bdo0zZ07V9OmTdP48eNVXFzs69YAXCbOb3wfAgqMVFxcrDlz5qihocHreENDg+bMmcOLGODHes7v+vp6r+P19fWc37A5LMuyfN3EYLW0tCgmJkbNzc3ciycAdXV1KSEhQY2NjeetGT16tOrq6pgOBvwM5/fwNpj3b2ZQYJytW7de8MVLkhobG7V169ar0xCAIcP5jYtFQIFxysrKhrQOgDk+++wzr+eRkZEKDw9XZGTkBesw/BBQYBwCChC4PvnkE6/nbW1tam9vV1tb2wXrMPwQUGCco0ePDmkdAHP0Xfh+uXUIXAQUGKe1tXVI6wCYw+FwDGkdAhcBBcY5e/bskNYBMEdnZ+eQ1iFwEVAAAFeNx+MZ0joELgIKAAAwDgEFAAAYh4AC47CIDghcQUHebzshISEaM2aMQkJCLliH4Ye/ABgnLi5uSOsAmOPOO+/0et7Z2an6+vp+i2L71mH4IaDAOGlpaUNaB8AcfTdku9w6BC4CCoxzsf9y4l9YgP+Jj48f0joELgIKjDN69OghrQNgjptuusnrec8di/veubhvHYYfAgqM43K5hrQOgDkKCgok/e8i966uLq//7TneU4fhi4AC4/QOHk6n02us93MCCuB/wsPDlZWVJcuy5HQ6ddtttyk9PV233XabnE6nLMtSVlaWwsPDfd0qfCzk+0sA37j22mt16tQpr2Mej2fA4wD8R0lJidxut0pLS7Vnzx6vsaysLJWUlPimMRiFgALjNDY2SpJOnTolh8OhjIwMTZs2TZ9//rnKy8vtcNJTB8D/lJSUqL29XTk5OTp48KAmTJiggoICZk5gI6DAOCNHjpR0btFcYmKiysrKVFZWJkkaN26campq1NXVZdcB8E/h4eFavXq1r9uAoViDAuNs3LhRkpScnKyDBw/q888/19q1a/X555/rL3/5i8aNG+dVBwAIPAQUGKe6ulqSdOjQIf30pz9VWFiYZsyYobCwMP30pz+1x3v+FwAQeAgoMM6ECRMkSQ899JAqKyuVnp6u6Ohopaenq6qqSllZWV51AIDA47Asy/J1E4PV0tKimJgYNTc3Kzo62tftYIi1t7crIiJCTqdTJ0+e1H/8x3+ovr5eY8aM0aRJkzRixAh5PB6dPn2aBXUA4EcG8/59WTMoeXl5cjgcys7Oto9ZlqWXX35ZCQkJCg8P17333qsDBw54fV9HR4cWLVqkUaNGKTIyUrNmzVJtbe3ltIIA0rNPgsfj0YgRI7RlyxZNnDhRW7ZsscMJ+yQAQGC75ICyc+dOvf/++7rlllu8jufn5+vNN9/U6tWrtXPnTrlcLmVkZKi1tdWuyc7O1oYNG7Ru3Tpt375dp06d0owZM+ydBIGSkhI7pOTn5+vGG29Ufn6+HU7YJwEAAtslBZRTp07p0Ucf1Zo1a7w+6mlZln71q1/phRde0OzZs5Wamqp//ud/1unTp7V27VpJUnNzs37zm9/ojTfe0P3336/bbrtNH3/8sSorK/Vv//ZvQ/NbISCUlJTo5MmTuueee5SUlKR77rlHJ0+eJJwAAeLYsWOKjY1VaGioYmNjdezYMV+3BINcUkBZsGCBpk+frvvvv9/r+OHDh9XQ0KDMzEz7WFhYmKZOnaodO3ZIknbv3q2zZ8961SQkJCg1NdWu6aujo0MtLS1eDwS+3NxcjRo1Sl9++aVqamr05ZdfatSoUcrNzfV1awAuU1hYmBITE9XU1KTOzk41NTUpMTFRYWFhvm4Nhhh0QFm3bp3+8z//U3l5ef3GGhoaJPW/TXZ8fLw91tDQIKfT2W+Trd41feXl5SkmJsZ+JCUlDbZt+Jnc3FwVFBQoLi5Oa9asUX19vdasWaO4uDgVFBQQUgA/FhYWJo/HI0kaM2aMfve732nMmDGSzt3OgpACaZABpaamRs8++6w+/vhjXXPNNeet67kbZQ/Lsvod6+tCNcuXL1dzc7P9qKmpGUzb8DMej0eFhYWKj49XbW2tnnzySblcLj355JOqra1VfHy8CgsL7Rc4AP7j2LFj9rl7/Phx1dXVad68eaqrq9Px48clnXsN4HIPBhVQdu/ercbGRk2cOFEhISEKCQlRRUWFfv3rXyskJMSeOek7E9LY2GiPuVwueTweNTU1nbemr7CwMEVHR3s9ELiKiorU2dmplStXKiTE+24MISEhWrFihTo7O1VUVOSjDgFcqrS0NEnnZk5iY2O9xmJjY+27lPfUYfgaVEC57777VFlZqb1799qPO+64Q48++qj27t2rH/zgB3K5XCovL7e/x+PxqKKiQunp6ZKkiRMnKjQ01Kumvr5eVVVVdg2Gt54dYmfMmDHgeM9xdpIF/E/PJzpfe+21AcdXrlzpVYfha1A3C4yKilJqaqrXscjISMXFxdnHs7OztWrVKk2YMEETJkzQqlWrFBERoblz50qSYmJiNH/+fC1dulRxcXGKjY3VL37xC6WlpfVbdIvhKSUlRZK0adMm/f3f/722bdtmb9Q2ZcoUbdq0yasOgP+IiopSU1OTli1bpnnz5vUbf/HFF+06DG+XvZPsvffeqx/+8If61a9+JencWpJXXnlF7733npqamjRp0iS9/fbbXsHmzJkzysnJ0dq1a9Xe3q777rtPRUVFF734lZ1kA5vH41FkZKQiIyM1YsQI/c///I89dsMNN+jkyZNqa2tTW1ubnE6nDzsFMFjHjh1TYmKipHNrUHpf5jlx4oTi4uIkSbW1tRo7dqxPesSVM5j3b7a6h5HcbrdKS0vldDq1ePFizZ8/X7/5zW/sxbFs1gb4r96f4nG5XFq5cqVefPFFe/2i0+lUR0eHL1vEFUJAgV/r6urS+PHjFRwcrCNHjnjtMBwSEqIbbrhB3d3dOnjwoIKDg33YKYBL1Tuk9EY4CWxX7V48wJWwbds2HTlyRB9//LFOnz6twsJCLVy4UIWFhWpra9NHH32kw4cPa9u2bb5uFcAl6ujoUG1trUaOHKmQkBCNHDlStbW1hBPYBrVIFrga6uvrJUmpqalyOp1eN6PsOd67DoB/Gjt2rE6cOOHrNmAoZlBgnJ4dJauqqgYc7zneUwcACDwEFBhnypQpGjdunFatWqXu7m6vse7ubuXl5Sk5OVlTpkzxUYcAgCuNgALjBAcH64033tCmTZvkdrv11VdfqbW1VV999ZXcbrc2bdqk119/nQWyABDAWIMCI82ePVvr16/X0qVLvXYYTk5O1vr16zV79mwfdgcAuNL4mDGM5vF4VFRUpOrqaqWkpOiZZ55hczYA8FODef9mBgXGKi4u1tKlS3XkyBH72FtvvaU33niDGRQACHCsQYGRiouLNWfOHKWlpXmtQUlLS9OcOXNUXFzs6xYBAFcQl3hgnJ6dZNPS0lRSUqKgoP/N0d3d3XK73aqqqmInWQDwM+wkC7/Ws5Ps888/7xVOJCkoKEjLly9nJ1kACHAEFBin906yA2EnWQAIfAQUGIedZAEABBQYh51kAQB8zBjG6dlJds6cOcrKytKDDz6o8PBwtbe367PPPtOnn36q9evXs0AWAAIYn+KBsXJzc1VYWKjOzk77WEhIiBYvXqz8/HwfdgYAuBRs1Aa/V1xcrNdff13Tp0/XT37yE3sGZcuWLXr99dd11113sVkbAAQwZlBgHPZBAYDAxD4o8GvsgwIAIKDAOOyDAgAgoMA47IMCACCgwDjsgwIAIKDAOD37oGzatElut9vrbsZut1ubNm3S66+/zgJZAAhgfMwYRpo9e7bWr1+vJUuWKD093T4+btw4rV+/no8YA0CAYwYFRnM4HL5uAQDgAwQUGKm4uFhz5sxRWlqa1yWetLQ0zZkzR8XFxb5uEQBwBbFRG4zDRm0AEJjYqA1+jY3aAAAEFBiHjdoAAAQUGIeN2gAABBQYh43aAADsgwLj9GzUNmfOHM2aNUspKSk6c+aMrrnmGlVXV2vz5s1av349C2QBIIARUGCk2bNna9asWSotLe03lpWVxUZtQABobm7W9OnTdfToUV1//fX69NNPFRMT4+u2YAgCCoyUm5ur0tJSxcfHa968efrBD36g//7v/9ZHH32k0tJS5ebmKj8/39dtArhE48ePV3V1tf28pqZGI0aMUEpKig4dOuTDzmAK9kGBcTwejyIjIxUXF6fa2lqFhPxvju7s7FRiYqKOHz+utrY2OZ1OH3YK4FL0DicPPvigXnrpJb3yyiv67LPPJImQEsDYBwV+raioSJ2dnVq5cqVXOJGkkJAQrVixQp2dnSoqKvJRhwAuVXNzsx1O2tratGXLFt11113asmWL2traJEnV1dVqbm72ZZswAAEFxul58ZoxY4a6urq0detW/f73v9fWrVvV1dWlGTNmeNUB8B/Tp0+XdG7mJCIiwmssIiJCmZmZXnUYvliDAuOkpKRIklasWKFNmzappqbGHktKSrJfuHrqAPiPo0ePSpJeeuklnThxQlOnTlVdXZ0SEhJUUVGhF198UWVlZXYdhi/WoMA4Ho9H4eHh/fZA6S0oKEjt7e2sQQH8zOTJk/Xll1/K6XTK4/H0G+85fs8992j79u0+6BBXEmtQ4NeCg4MvGE6kcxu2sQ8K4H8+/fRTSRownPQ+3lOH4YuAAuNs3LhxSOsAmKOrq2tI6xC4CCgwzvz584e0DoA57rrrriGtQ+AioMA4TU1NQ1oHwBwHDx60v+57Gaf38951GJ5YJAvjOByOi671wz9fYFjj/B7eWCSLgBEbG6uEhASNGDFCCQkJio2N9XVLAIZQcHCwli1bxqJ39ENAgdFOnDihuro6nTx5UnV1dTpx4oSvWwJwGUaNGuX1fOLEiZo2bZomTpx4wToMP1zigXEiIyN1+vTp762LiIiwt8YG4B+4xDO8cYkHfu3HP/7xkNYBAPwPAQXGqaysHNI6AOa42BmUwcy0IDARUGCclpaWIa0DYI59+/bZX+/atUvjxo1TZGSkxo0bp127dg1Yh+GJgALjxMTEDGkdAHOkpaXZX99xxx2qqanRk08+qZqaGt1xxx0D1mF4IqDAOLm5uUNaB8AsvRe/dnV16a233vLa2p7FsZAIKDDQgQMHvJ67XC4tXLhQLpfrgnUA/IdlWdq/f7+91sThcGj//v2EE9hCfN0A0FdnZ6fX84aGBq1evfp76wD4l7S0tO+9czmGL2ZQYJxvv/3W/jozM1MjR46U0+nUyJEjlZmZOWAdACCwMIMC4/RcynE6nSorK7OPezwelZWVyel0yuPx9LvkAwAIHMygwDg33nijpHOBRJIyMjK0atUqZWRkeB3vqQMABB62uodxTp06paioKEnnbiTWe3V/7+etra269tprfdIjAGDwBvP+zSUeGOeDDz6wv46Li9PUqVMVGRmptrY2VVRUqLGx0a7Lzs72UZcAgCuJgALjVFdXS5KefvpprVmzRn/4wx/ssZCQED399NN655137DoAQOBhDQqMk5KSYn+dkJDgNZaQkGDvk9C7DgAQWFiDAuN4PB6Fh4eru7tb06dP14svvqjU1FRVVVVp5cqV+vTTTxUUFKT29nY5nU5ftwsAuEiDef9mBgXGCQ4Othe/7tq1S/v371dLS4v2799v30zs2muvVXBwsC/bBABcQQQUGGfbtm1qaWnRo48+quPHj+upp57S2LFj9dRTT+n48eOaO3euWlpatG3bNl+3CgC4QggoME59fb0k6d1331VbW5sKCwu1cOFCFRYWqq2tTe+++65XHQAg8AwqoLzzzju65ZZbFB0drejoaN19993asmWLPW5Zll5++WUlJCQoPDxc9957b78bunV0dGjRokUaNWqUIiMjNWvWLNXW1g7Nb4OAMGbMGElSVVWVnE6nsrOz9U//9E/Kzs6W0+lUVVWVVx0AIPAMKqAkJibq1Vdf1a5du7Rr1y79+Mc/VlZWlh1C8vPz9eabb2r16tXauXOnXC6XMjIy1Nraav9/ZGdna8OGDVq3bp22b9+uU6dOacaMGV6bcWF4mzJlisaNG6dVq1b1u5FYd3e38vLylJycrClTpvioQwDAFWddppEjR1offPCB1d3dbblcLuvVV1+1x86cOWPFxMRY7777rmVZlnXy5EkrNDTUWrdunV1z7NgxKygoyPrss88u+mc2Nzdbkqzm5ubLbR+G+uSTTyyHw2HNnDnT2rFjh9XS0mLt2LHDmjlzpuVwOKxPPvnE1y0CuEwHDhywgoKCLElWUFCQdeDAAV+3hCtsMO/fl7wGpaurS+vWrVNbW5vuvvtuHT58WA0NDV53mw0LC9PUqVO1Y8cOSdLu3bt19uxZr5qEhASlpqbaNYAkzZ49W+vXr1dlZaXS09MVHR2t9PR0VVVVaf369Zo9e7avWwRwGRwOh26++WZ7lrS7u1s333yzHA6HjzuDKQYdUCorK3XttdcqLCxMP/vZz7Rhwwb93d/9nRoaGiRJ8fHxXvXx8fH2WENDg5xOp0aOHHnemoF0dHSopaXF64HAN3v2bFVWVsrtdistLU1ut1v79+8nnAB+rncICQ0N1S9/+UuFhoYOOI7ha9AB5cYbb9TevXv19ddf6+mnn9bjjz+u//qv/7LH+/5hWZb1vX9s31eTl5enmJgY+5GUlDTYtuGH3G63oqKiVFJSosrKSpWUlCgqKkput9vXrQG4RL3fL2pqauTxeLRixQp5PB7V1NQMWIfhadABxel0avz48brjjjuUl5enW2+9VW+99ZZcLpck9ZsJaWxstGdVXC6XPB6PmpqazlszkOXLl6u5udl+9P4jRmByu90qLS2V0+nUc889p0OHDum5556T0+lUaWkpIQXwU2lpaZLOzZwkJiZ6jSUmJtozKT11GL4uex8Uy7LU0dGh5ORkuVwulZeX22Mej0cVFRVKT0+XJE2cOFGhoaFeNfX19aqqqrJrBhIWFmZ/tLnngcDV3t5uh5PW1lbl5eUpJSVFeXl5am1ttUNKe3u7r1sFMEg9a06ee+65AccXL17sVYfha1AB5fnnn9e2bdt05MgRVVZW6oUXXtDWrVv16KOPyuFwKDs7W6tWrdKGDRtUVVWlJ554QhEREZo7d64kKSYmRvPnz9fSpUv1pz/9SXv27NFjjz2mtLQ03X///VfkF4T/ycnJkSQtWbKk3712evZF6V0HwH8EBZ1723n11VcHHC8sLPSqw/AVMpjib7/9VvPmzVN9fb1iYmJ0yy236LPPPlNGRoYkKTc3V+3t7XrmmWfU1NSkSZMmqaysTFFRUfb/R2FhoUJCQvTwww+rvb1d9913nz788EPuqwLbwYMHJUlPPvnkgOPz589Xfn6+XQfAf1RWVurmm2/W2bNnVVtb63WZp7a2VmfPnrXrMLxxN2MYZ+HChXr77bf13HPPKS8vr9/4smXLlJ+frwULFmj16tU+6BDA5ej7KZ7FixersLDQDifSueUDCDyDef8moMA47e3tioiIsNeg9L7M4/F4FBUVJY/Ho9OnTys8PNyHnQK4VBf65KYfvi3hIg3m/ZuLfDBOeHi4srKy7DCybNky/eUvf9GyZcvscJKVlUU4AfyYZVk6cOCAvdYkKChIBw4cIJzAxgwKjNXzUeO+srKyVFJScvUbAgBclsG8fw9qkSxwNZWUlKi9vV05OTk6ePCgJkyYoIKCAmZOAGAYIKDAaOHh4SyEBYBhiDUoAADAOAQUAABgHAIKAAAwDgEFAAAYh4ACAACMQ0ABAADGIaAAAADjEFAAAIBx2KgNRvN4PCoqKlJ1dbVSUlL0zDPPeN08EAAQmAgoMFZubq4KCwvV2dlpH8vJydHixYuVn5/vw84AAFcal3hgpNzcXBUUFCguLk5r1qxRfX291qxZo7i4OBUUFCg3N9fXLQIAriDuZgzjeDweRUZGKi4uTrW1tQoJ+d+Jvs7OTiUmJur48eNqa2vjcg8A+JHBvH8zgwLjFBUVqbOzUytXrvQKJ5IUEhKiFStWqLOzU0VFRT7qEABwpRFQYJzq6mpJ0owZMwYc7zneUwcACDwEFBgnJSVFkrRp06YBx3uO99QBAAIPa1BgHNagAEBgGsz7Nx8zhnGcTqcWL16sgoICjR07Vo899phSUlJUXV2tjz/+WI2NjcrJySGcAEAAYwYFxnK73SotLe13PCsrSyUlJVe/IQDAZeFTPPB7xcXFA4YTSSotLVVxcfFV7gjAUCsvL5fD4bAf5eXlvm4JBmEGBcbp6urSNddc47WDbF8hISE6c+aMgoODr2JnAIaKw+E475gfvi3hIjGDAr+2efPmC4YT6dxi2c2bN1+ljgAMpb7hZObMmRccx/BEQIFxfv7znw9pHQBz9L6Ms2/fPlmWpY0bN8qyLO3bt2/AOgxPXOKBcQbzryc//PMFhrXe5/dA5+/3jcO/cYkHAGC0vpd1emRmZl7lTmAqAgoA4Kr713/9V+3du9frUzx79+5VWVmZr1uDIQgoAICrpncAue2227zGej8nqICAAuOEhoYOaR0Ac2RkZPQ7Nnbs2Iuqw/BCQIFxEhIShrQOgDn27t3b79ixY8cuqg7DCwEFxomJiRnSOgDm6LmMExQU1O8yTllZmYKCgrzqMHwRUGCcEydODGkdAPM8++yzysjIkGVZ9iMjI0NPPfWUr1uDIQgoMM6IESOGtA6Aed56660Bj7/33ntXuROYioAC4zz88MP21z3TvT1633undx0A/7Bnzx5JUnd3t44cOeI1duTIEXV3d3vVYfgioMA4vRfM9bxY9ejq6hqwDoB/+OEPf2h/nZycrODgYD3zzDMKDg5WcnLygHUYnkJ83QDQ1+HDh4e0DoBZLMuyt7Tv7u7WO++8028cYAYFxpkwYYIkadKkSf3uyxMUFKQ777zTqw6A/7Esq99lnD179hBOYONmgTBOe3u7IiIi5HQ6dfz4cX3wwQeqrq5WSkqKnnzyScXFxcnj8ej06dMKDw/3dbsAgIvEzQLh18LDw5WVlSWPx6O4uDjV19dr0aJFqq+vt8NJVlYW4QQAAhgzKDCW2+1WaWlpv+NZWVkqKSm5+g0BAC7LYN6/WSQLY5WUlKi9vV05OTk6ePCgJkyYoIKCAmZOAGAYIKDAaOHh4Vq9erWv2wAAXGWsQYHRmpubNXnyZF1//fWaPHmympubfd0SAOAqYAYFxho/fryqq6vt5zU1NRoxYoRSUlJ06NAhH3YGALjSmEGBkXqHkwcffFBfffWVHnzwQUlSdXW1xo8f78v2AABXGJ/igXGam5vtGwG2tbUpIiLCHjt9+rQiIyMlSSdPnlRMTIwvWgQAXAL2QYFfmz59uqRzMye9w4kkRUREKDMz06sOABB4CCgwztGjRyVJL7300oDjL774olcdACDwEFBgnOuvv16S9Morrww4vnLlSq86AEDgYQ0KjMMaFGB4+Oabb3TzzTerq6tLwcHBOnDggG688UZft4UriJ1k4ddiYmKUkpKi6upqRUZGKikpSfHx8fr2229VU1MjSUpJSSGcAH6s753Ku7q6dNNNN0kSdzSGJC7xwFCHDh2yZ05qamq0a9cuO5xERESwDwrgx3qHk7CwMK1cuVJhYWEDjmP4IqDASG63W6dPn1ZoaKjGjh0rl8ulsWPHKjQ0VKdPn5bb7fZ1iwAuwTfffGN/fezYMZ05c0YvvPCCzpw5o2PHjg1Yh+GJNSgwTnt7uyIiIuR0OtXa2iqn02mPeTweRUVFyePx6PTp09w4EPAzISEh6urqUlhYmM6cOdNv/JprrlFHR4eCg4PV2dnpgw5xJbEPCvxaTk6OJGnJkiXyeDx66KGHdMstt+ihhx6Sx+NRdna2Vx0A/9HV1SVJ+uUvfzngeG5urlcdhi9mUGCcBx54QGVlZUpLS1NlZWW/8Z7jmZmZ+uMf/+iDDgFcKmZQhjdmUODXJkyYIEmqrKyUw+HQvHnztG/fPs2bN08Oh8MOLT11APzHgQMHJEkdHR2qq6vzGqurq1NHR4dXHYYvZlBgnL/+9a8aPXq0JPX7b9zz316SGhsbdd111/mkRwCXru+neHJzc5Wfn2+HE4mPGgcqZlDg1/7hH/7B/nrUqFG6/fbbNXnyZN1+++0aNWrUgHUA/Efv8NHR0aF//Md/JJygHwIKjFNdXS1JiouL09mzZ7Vnzx59+eWX2rNnj86ePavY2FivOgD+x7Is/fnPf1ZwcLAkKTg4WH/+858JJ7ARUGCclJQUSdLx48cHHD9x4oRXHQD/dOONN6qzs1OWZamzs5Nt7uGFgALjvP/++17PR40aJZfL5XV5Z6A6AP6lvLxcDofDfpSXl/u6JRiEe/HAOM8//7zX8+++++68dWvWrLkaLQEYYgNtZ5+ZmSmJNSg4hxkUGKe0tHRI6wCYpW84mTlz5gXHMTwNKqDk5eXpzjvvVFRUlEaPHi23293vfgmWZenll19WQkKCwsPDde+99/b7PHtHR4cWLVqkUaNGKTIyUrNmzVJtbe3l/zYICGfPnh3SOgDm6H0ZZ9++fbIsSxs3bpRlWdq3b9+AdRieBhVQKioqtGDBAn399dcqLy9XZ2enMjMz1dbWZtfk5+frzTff1OrVq7Vz5065XC5lZGSotbXVrsnOztaGDRu0bt06bd++XadOndKMGTPY2hiSdNG7R7LLJOB/ei7jSNLOnTu91qDs3LlzwDoMT5e1UVvPhloVFRX60Y9+JMuylJCQoOzsbC1btkzSudmS+Ph4vfbaa3rqqafU3Nys6667Th999JEeeeQRSed2D0xKStLmzZv1wAMPfO/PZaO2wDaY6V2uVQP+hfN7eLtqG7U1NzdLkr0vxeHDh9XQ0OCVfMPCwjR16lTt2LFDkrR7926dPXvWqyYhIUGpqal2TV8dHR1qaWnxegAAAsPf/M3f+LoFGOiSA4plWVqyZIkmT56s1NRUSVJDQ4MkKT4+3qs2Pj7eHmtoaJDT6dTIkSPPW9NXXl6eYmJi7EdSUtKltg0A8KGf//zn9tcFBQWyLEvffPONLMtSQUHBgHUYni45oCxcuFD79+/X73//+35jfafwLMv63mm9C9UsX75czc3N9qOmpuZS2wYA+NCvf/1r++ucnByvNSg5OTkD1mF4uqSAsmjRIm3cuFGff/65EhMT7eMul0uS+s2ENDY22rMqLpdLHo9HTU1N563pKywsTNHR0V4PAAAQuAYVUCzL0sKFC1VcXKx///d/V3Jystd4cnKyXC6X18fDPB6PKioqlJ6eLkmaOHGiQkNDvWrq6+tVVVVl1wAAgOFtUAFlwYIF+vjjj7V27VpFRUWpoaFBDQ0Nam9vl3Tu0k52drZWrVqlDRs2qKqqSk888YQiIiI0d+5cSVJMTIzmz5+vpUuX6k9/+pP27Nmjxx57TGlpabr//vuH/jcEABjj4Ycftr9+8MEHvcZ6P+9dh+FpUB8zPt8akd/+9rd64oknJJ2bZXnllVf03nvvqampSZMmTdLbb79tL6SVpDNnzignJ0dr165Ve3u77rvvPhUVFV304lc+ZhzY+BgiELg4v4e3wbx/X9Y+KL5CQAlsvIABgYvze3i7avugAAAAXAkEFACAT40fP97XLcBABBQAgE8dOnTI1y3AQAQUAABgHAIKAAAwDgEFAAAYh4ACAACMQ0ABAADGIaAAAADjEFAAAIBxCCgAAMA4BBQAAGAcAgoAADAOAQUAABiHgAIAAIxDQAEAAMYhoAAAAOMQUAAAgHEIKAAAwDgEFAAAYBwCCgAAMA4BBQAAGIeAAgAAjENAAQAAxiGgAAAA4xBQAACAcQgoAADAOAQUAABgHAIKAAAwDgEFAAAYh4ACAACMQ0ABAADGIaAAAADjEFAAAIBxCCgAAMA4BBQAAGAcAgoAADAOAQUAABiHgAIAAIxDQAEAAMYhoAAAAOMQUAAAgHEIKAAAwDgEFAAAYBwCCgAAMA4BBQAAGIeAAgAAjENAAQAAxiGgAAAA4xBQAACAcQgoAADAOAQUAABgHAIKAAAwDgEFAAAYh4ACAACMQ0ABAADGIaAAAADjEFAAAIBxCCgAAMA4BBQAAGAcAgoAADAOAQUAABiHgAIAAIxDQAEAAMYhoAAAAOMQUAAAgHEIKAAAwDiDDihffPGFZs6cqYSEBDkcDpWUlHiNW5all19+WQkJCQoPD9e9996rAwcOeNV0dHRo0aJFGjVqlCIjIzVr1izV1tZe1i8CAAACx6ADSltbm2699VatXr16wPH8/Hy9+eabWr16tXbu3CmXy6WMjAy1trbaNdnZ2dqwYYPWrVun7du369SpU5oxY4a6urou/TcBAAABw2FZlnXJ3+xwaMOGDXK73ZLOzZ4kJCQoOztby5Ytk3RutiQ+Pl6vvfaannrqKTU3N+u6667TRx99pEceeUSSVFdXp6SkJG3evFkPPPDA9/7clpYWxcTEqLm5WdHR0ZfaPgzlcDguuvYy/nwB+ADn9/A2mPfvIV2DcvjwYTU0NCgzM9M+FhYWpqlTp2rHjh2SpN27d+vs2bNeNQkJCUpNTbVr+uro6FBLS4vXAwAABK4hDSgNDQ2SpPj4eK/j8fHx9lhDQ4OcTqdGjhx53pq+8vLyFBMTYz+SkpKGsm0AAGCYK/Ipnr5TeJZlfe+03oVqli9frubmZvtRU1MzZL0CAADzDGlAcblcktRvJqSxsdGeVXG5XPJ4PGpqajpvTV9hYWGKjo72egAAgMA1pAElOTlZLpdL5eXl9jGPx6OKigqlp6dLkiZOnKjQ0FCvmvr6elVVVdk1AABgeAsZ7DecOnVKhw4dsp8fPnxYe/fuVWxsrK6//nplZ2dr1apVmjBhgiZMmKBVq1YpIiJCc+fOlSTFxMRo/vz5Wrp0qeLi4hQbG6tf/OIXSktL0/333z90vxkAAPBbgw4ou3bt0rRp0+znS5YskSQ9/vjj+vDDD5Wbm6v29nY988wzampq0qRJk1RWVqaoqCj7ewoLCxUSEqKHH35Y7e3tuu+++/Thhx8qODh4CH4lAADg7y5rHxRfYR+UwMY+CUDg4vwe3ny2DwoAAMBQIKAAAADjEFAAAIBxCCgAAMA4BBQAAGAcAgoAADAOAQUAABiHgAIAAIxDQAEAAMYhoAAAAOMQUAAAgHEIKAAAwDgEFAAAYBwCCgAAMA4BBQAAGIeAAgAAjENAAQAAxiGgAAAA4xBQAACAcQgoAADAOAQUAABgHAIKAAAwDgEFAAAYh4ACAACMQ0ABAADGIaAAAADjEFAAAIBxCCgAAMA4BBQAAGAcAgoAADAOAQUAABiHgAIAAIxDQAEAAMYhoAAAAOMQUAAAgHEIKAAAwDgEFAAAYBwCCgAAMA4BBQAAGIeAAgAAjENAAQAAxiGgAAAA4xBQAACAcQgoAADAOAQUAABgHAIKAAAwDgEFAAAYh4ACAACMQ0ABAADGIaAAAADjEFAAAIBxCCgAAMA4BBQAAGAcAgoAADAOAQUAABiHgAIAAIxDQAEAAMYhoAAAAOMQUAAAgHEIKAAAwDgEFAAAYBwCCgAAMA4BBQAAGIeAAgAAjENAAQAAxiGgAAAA4/g0oBQVFSk5OVnXXHONJk6cqG3btvmyHQAAYAifBZR/+Zd/UXZ2tl544QXt2bNHU6ZM0U9+8hMdPXrUVy0BAABDOCzLsnzxgydNmqTbb79d77zzjn3sb//2b+V2u5WXl3fB721paVFMTIyam5sVHR19pVvFVeZwOC661kd/vgAuEef38DaY9++Qq9STF4/Ho927d+u5557zOp6ZmakdO3b0q+/o6FBHR4f9vKWl5Yr3iMv33Xff6Y+f/E4RXef/73X6dJuqq//b69htrouf2Fvx9E/7HUtJ+YEiIiLP+z2jkm/WlJ/8n4v+GQD64/zGleaTgPLdd9+pq6tL8fHxXsfj4+PV0NDQrz4vL0+vvPLK1WoPQ6SkpES1v39eL98bduFC7z8D/d+nrh3ET/m3/odO/f+P83j5/3XouuQ03XTTTYP4OQB64/zGleaTgNKj71SfZVkDTv8tX75cS5YssZ+3tLQoKSnpiveHy+N2u/XHrhZtGOS/sEpKSgb1M/r6vn9h3bfsZl68gMvE+Y0rzSdrUDwejyIiIvSHP/xBDz30kH382Wef1d69e1VRUXHB72cNSmDjGjUQ2C7mHOfcDkyDef/2yad4nE6nJk6cqPLycq/j5eXlSk9P90VLMMjFvjDxAgb4p+87dzm3IfnwEs+SJUs0b9483XHHHbr77rv1/vvv6+jRo/rZz37mq5ZgkPNd7us9DsB/ne8c59xGD58FlEceeUTHjx/XihUrVF9fr9TUVG3evFk33HCDr1qCYXgBAwIb5zIuxGf7oFwO1qAAAOB/jF+DAgAAcCEEFAAAYBwCCgAAMA4BBQAAGIeAAgAAjENAAQAAxiGgAAAA4xBQAACAcQgoAADAOD7b6v5y9Gx+29Jy/tt8AwAAs/S8b1/MJvZ+GVBaW1slSUlJST7uBAAADFZra6tiYmIuWOOX9+Lp7u5WXV2doqKiLnjHWwSGlpYWJSUlqaamhnsvAQGG83t4sSxLra2tSkhIUFDQhVeZ+OUMSlBQkBITE33dBq6y6OhoXsCAAMX5PXx838xJDxbJAgAA4xBQAACAcQgoMF5YWJheeuklhYWF+boVAEOM8xvn45eLZAEAQGBjBgUAABiHgAIAAIxDQAEAAMYhoAAAAOMQUGCsL774QjNnzlRCQoIcDodKSkp83RKAIZKXl6c777xTUVFRGj16tNxut7755htftwWDEFBgrLa2Nt16661avXq1r1sBMMQqKiq0YMECff311yovL1dnZ6cyMzPV1tbm69ZgCD5mDL/gcDi0YcMGud1uX7cC4Ar461//qtGjR6uiokI/+tGPfN0ODMAMCgDA55qbmyVJsbGxPu4EpiCgAAB8yrIsLVmyRJMnT1Zqaqqv24Eh/PJuxgCAwLFw4ULt379f27dv93UrMAgBBQDgM4sWLdLGjRv1xRdfKDEx0dftwCAEFADAVWdZlhYtWqQNGzZo69atSk5O9nVLMAwBBcY6deqUDh06ZD8/fPiw9u7dq9jYWF1//fU+7AzA5VqwYIHWrl2r0tJSRUVFqaGhQZIUExOj8PBwH3cHE/AxYxhr69atmjZtWr/jjz/+uD788MOr3xCAIeNwOAY8/tvf/lZPPPHE1W0GRiKgAAAA4/AxYwAAYBwCCgAAMA4BBQAAGIeAAgAAjENAAQAAxiGgAAAA4xBQAACAcQgoAADAOAQUAABgHAIKAAAwDgEFAAAYh4ACAACM8/8BHL7523JBYWUAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import pyBigWig\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "bw1 = pyBigWig.open(\"/mnt/10tb/home/tsoi/DashaP/data/ATACseq_Egor/inc_qnorm_genome_bychr/Sorted008a9efa5f80a61f25c0b54c9b8e8cf0.forward.incqn.bw\")\n",
    "bw2 = pyBigWig.open(\"/mnt/10tb/home/tsoi/DashaP/data/ATACseq_Egor/inc_qnorm_genome_bychr/Sorted0c5ec366b5af84de8e731d40407231a3.forward.incqn.bw\")\n",
    "\n",
    "bw1 = pyBigWig.open(\"/mnt/10tb/home/tsoi/DashaP/data/ATACseq_Egor/bw/Sorted008a9efa5f80a61f25c0b54c9b8e8cf0.forward.bw\")\n",
    "bw2 = pyBigWig.open(\"/mnt/10tb/home/tsoi/DashaP/data/ATACseq_Egor/bw/Sorted0c5ec366b5af84de8e731d40407231a3.forward.bw\")\n",
    "\n",
    "chrom = list(bw1.chroms().keys())[0]              # первая хромосома\n",
    "end = bw1.chroms()[chrom]                         # её длина\n",
    "vals1 = bw1.values(chrom, 0, end)\n",
    "vals2 = bw2.values(chrom, 0, end)\n",
    "\n",
    "plt.hist([vals1, vals2], bins=50, label=[f\"1, {max(vals1)}\", f\"2, {max(vals2)}\"], alpha=0.6)\n",
    "plt.yscale(\"log\")\n",
    "plt.legend(); plt.show()\n",
    "\n",
    "plt.boxplot([vals1, vals2])\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "a0ed7384",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "base",
   "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.12.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}