{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "import os\n",
    "import sys\n",
    "sys.path.append('../')\n",
    "sys.path.append('../anmodel')\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "import anmodel"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Help on package anmodel:\n",
      "\n",
      "NAME\n",
      "    anmodel\n",
      "\n",
      "PACKAGE CONTENTS\n",
      "    analysis\n",
      "    channels\n",
      "    models\n",
      "    params\n",
      "    search\n",
      "\n",
      "FILE\n",
      "    /Users/tyamada/ueda_lab/averaged_neuron/anmodel/__init__.py\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "help(anmodel)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [],
   "source": [
    "an = anmodel.models.ANmodel()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [],
   "source": [
    "an.set_sws_params()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [],
   "source": [
    "s, info = an.run_odeint()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0, 0.5, 'Membrane Potential (mV)')"
      ]
     },
     "execution_count": 33,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEGCAYAAACO8lkDAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO2dd3wc1bX4v0e9S7ZsufeCKzZgTO81NBMgBJIQEpJHCIQS3nsJhBd+kEcSAqGThJCE+uid0F3AdIxt3Kvc5SLLlq3edvf+/phZeSWttLNlyq7u9/PRR7uzsztnZu6cc885954rSik0Go1GowklzW0BNBqNRuM9tHHQaDQaTRe0cdBoNBpNF7Rx0Gg0Gk0XtHHQaDQaTRcy3BYgEfTr10+NHDnSbTE0Go0mqVi0aNEepVT/cJ+lhHEYOXIkCxcudFsMjUajSSpEZEt3n+mwkkaj0Wi6oI2DRqPRaLqgjYNGo9FouqCNg0aj0Wi6oI2DRqPRaLqgjYNGo9FouqCNg0aj0Wi6oI2DRtOJmsY2/r10h9tiaDSuoo2DRVp9ATZW1bsthsYBbnjhG6597hs27WlwWxSNxjW0cbDIb15bzsn3zGdfQ6vbomhsZsf+ZgBafH6XJdE4gT+gaGrV97oz2jhY5IsNewFoaPW5LElqcd/sdYy86W38Ab0iocYdrn/+Gybe+p7bYngObRw0rvK3jzYA4AsEOmz3+QOMvvltnl+w1Q2xNAmifHcd+xu97W2/tWyn2yL0SHObn6e/3ELA4Q6UNg42sKe+xW0Rkp6GFj8BBX94Z7Xbomji4NR7P+achz51W4yk5v456/nt6yt4a7mzRkwbhwTz/spdzLhjDp+X73FbFP78/lpueW05H67Z7bYoGhv5aO1uT3dIKvY1uS1CUhP0vBpanA1pp0TJbi+xaMs+AJZvr+Hosf1cleXhD8sBeOarrWy+82xXZdEknv2NrZx5/yfsqm1m/IACPvjlCW6LpEkhXPMcRGSYiHwoIqtFZKWIXG9u7ysis0Vkvfm/j1syano3yuM58k/W72FXrTGyamOVHnarSSxuhpV8wH8qpSYCRwLXiMgk4CZgrlJqHDDXfO8ZolEYlbXNfLK+yj5hNJokYd6ayrDb/QHF6p21jsqyckeNp8NwocxbU8nzX29z5diuGQel1E6l1GLzdR2wGhgCzAKeNHd7EjjfHQnjQwHnPvQpl/1rgduieJZWX4BWvzFKyYu9dBG3JUgdrngi/EqND88r51sPfMKK7TWOyXL2g59y2r3zHTtePPzXS8tcO7YnEtIiMhI4BPgKGKCU2gmGAQHKuvnOlSKyUEQWVlW52zufu7qSv5jx/VB9srsuOXonTjLp1ve45pnFgGE8I+FBm6FJIMsq9gOwq6bZ0ePua2xz9HiJwOkOlOvGQUQKgFeAG5RSlv1LpdSjSqkZSqkZ/fuHXR/bFsL1Jn/y5ELufn+tYzIkM42tft42h+Strazrfkfda9doXMVV4yAimRiG4Rml1Kvm5koRGWR+PghIynGYXgyTeJkuRldfv6jQl0uTaNwcrSTAv4DVSql7Qz56E7jcfH058IbTsvUWdtc1c9XTi6h3ePx0NGgHonfglnG79wPt8XeHm57DMcBlwMkissT8Owu4EzhNRNYDp5nvk45kSGY+MGc9763cxWvfbHdbFE/ide8vGdpYJNw+hwfnlbsrQBj++O5qrnl2cZftTl8r1ybBKaU+pfuO4SlOyhINPSmMkTe9zZmTBzonTArhJUXstsKySug1SxKRNRb4+/yNAPzle+7K4XpCOtWYs9oYz+0lZafReB2lH5iwuGn0tXFIMMnS60wW3FAZWk9pvEJoU+x1Q1k1iWVbdWNqrI3gkpGtrG1uH2LrRUNf29zG6ymYIxIvXuxejjYONqFc6PPu2N/EcXd9qOdcxMERf5jrtgg98quXlnHDC0u6lJxIge6AxmNo4xAlkTo44mKUsMqckf35hijLhes4StKw0yy019yWWsta6pyD99DGIUqCbdgfUPz29RVsq27s+LnZh3v8s80OS5bceHESnNf1VWpEYlLiJFISvZ5DjHyzdR9Pf7mFVZ3de1OhVCVTXSUPaJnuFLH7knWlqdVPi89PSV6W26JoUhw9WimF8ICe1djMOQ99wvTfzSYQUI6v69sdXml2W/c2MvKmt/n30h1ui6KJE20cEoybOYdUxBuqtyMbzIV1jr5zHtN/94ErMii8GfYKetLRGgcPnkqvR4eVUhAvKo2oSQIbG1yFzUmS4LJEhfa0e8bNR1l7DhpNBGqbvVP7vztlkaz9gZToyKQo2jjESZcheB7oCUXdG9NPaFh21jTxWfkeDr7tA88t9yqkVq87hU4lZdDGIUZS4cFMhXOwkyueWMgriysA+HrzPpelSQ5qmlqB6D0Z3T3xHto4xEi3ne0kauXaYYhMpQt5hWRl/roqfv3K8qi+ozsoPRN6eZyuuqCNQwwEAoqfPmUsmO7FmjBRK30PnoOmZ7xo1xdt0d5VKqGNQwzUt/rY390C5S7pWX9A8ejHG905uI20+QPUeSgh7CaBgKLNF3BbDFvQXmxknB4mr41DnHilJswbS7bz9vKdbouRcJrbAky9zZ25BF7j9n+v7DAj3yNNLyxW1Zj2Wb2LNg4pQlM8hdi8rGW8ggeu0TNfbXVbBMu4f7U08aKNQ4LxQk+ozW8t9KBTDZHxgE3oFiv378YXljDyprftF6YTs1dV0tDii7ifhy9vr0cbhwTjhca+ZlcdO/Y3uS1GauExS/qjxxZY2u9VFxcG2p1MxSc9ws2vLnNbhHa0cbCIx3RDRDbvbYi4j5d7xV7Bq/e9tjlyrzwZ8OjldY3nFmxrf71mV+cFnZx9YC3XVhKRPsBgoAnYrJRKzWET3RBUpPUtPlp8qbXQitc0oFKKE+/+0G0xNBpXOfP+TyjNd68sfI/GQUSKgWuAS4EsoArIAQaIyJfAX5VSveopPuvBT3rscXtLzXqXPfXdhxwCCvZ1N1TYYbR3FTteGcmniY1InsPLwFPAcUqp/aEfiMhhwGUiMlop9S+7BPQaSdPee5BzY1U9J98zn5mj+jonD/DwvPXtr2+Kciat6yTNjU8eappCOwDOX98V22scP2Yy0aNxUEqd1sNni4BFCZcoyVi9s85tEaLmrWXGfIgFm6qNDQ4pvj9/sK79dSqF5nbsb2JwSa7jx03mtUPeWb6Tq59ZzICibEePG+qxPvH5ZkePHS+emgQnIqtE5BYRGeOUQMlG5/kFHgvfaxJFDzf26Dvn8en6PQ4KY+B0gjKRfFZuXK/KWudGNLX5A8y4Y45jx0t2Io1WuhQoAD4Qka9E5AYRGeyAXElLc5tH8vRJZqRCnZckEx2A1Z3WEtd4D79HlnSNFU8V3lNKLVVK3ayUGgNcD4wAvhSReSLyH45I6BHi8QicSMx1cTmT+DlIYtE1SURtkzcGPXgVy/MclFJfKqV+CfwQ6AM8bJtUGlvxQs98b32r2yJYYrvFyYRO9Op0yDKxzFld6bYInsaScRCRw0XkXhHZAtwOPAoMsVMwETlTRNaKSLmI3GTnsVKBaJRTlz1d0DrrKpMjkb9lb6PxwgOjlTwgQsLQhs77RJrn8Afgu8A+4HngGKVUhd1CiUg68BfgNKAC+FpE3lRKrbL72N2RSg+m15l463tuixA1brePzqHLp7/YzGOfbXZFliA9XZLO18uN6ycilg787vKdHDe+PwXZlucMJ4y9De552JE8hxbgW0qpGUqpPzthGExmAuVKqY1KqVYMwzTLoWMnHLcVR2e6rPvgMQFbPbhmwYPzynl5kVPN3xo9DW387Rsr2bQncgmVRJLszkBTq5+KfY0dtq3dVcfPn1nMza8m2bycBBApIX27UmqdiKSLyHkicp2I3Bj8s1GuIcC2kPcVdApjiciVIrJQRBZWVXlr8Xc3sDoGuqapjXoL1TLtJhnDCv/10tJuP3PjfPY1HuhVenFFws4SvW3Or1FKebL8+I8eX8Cxf+pY8CH4rHQ2Gr0BqwnpfwM/AkqBwpA/uwjX0jt0b5VSj5oezYz+/fvbKIopkPeevR5Zs6ubmL63nIR2knnMPjjjfHVug//z+gr7DxoFkS7BNc8u5v456/h6szeXE/0qOCnUo7S0BfjXp5scG5JrNYg2VCl1sK2SdKQCGBZ6fGCHg8dPeu54u2N65qdPLuSv3z/UJWkis2F3Ax+vr+Kak8a6LUrScs8Hazlj8kCmDCl2WxQgvLG4f856ng3jNbjRNYjU3/MHFBf+7XPAG5HXhz8sp7qhleLcTC46bCh761vIzEijKCfTluNZ9RzeFZHTbZEgPF8D40RklIhkAZcAbzp4/ITiRrvqfMw5qytZVrE/7L5e4OK/f8Hd76/1RMgrFpy4xz0pKKUUD80r59yHP3VAEms0toQvkeJmkjUa4lpd0QaCtagaW41n5LA75nDE7+fadjyrnsOXwGsikga0YRhdpZQqskMopZRPRH4BvA+kA48ppVbacaxUwUpY5qJHvuCqE7xRCcXIkRyQuTmFai25QfBKeqGHG6SnHE1nkixq6yqh99hOA2bVONwDHAUsVw7V4VVKvQO848SxehP/+nRj5J1cINiqdJnn+Fm6zR0PsbOCX7c7OeayWMELOUenRbAaVloPrHDKMGiip/NoJQ+0ZU2C6UlBhT6ZFfuSb4lYrysWL2k+p9SwVc9hJ/CRiLyLMfcBAKXUvbZIlWIYN1Or61CSfXSSl/lo7W5Xjmv1jnqlj+kNKbyLVeOwyfzLMv80Hqe7hu+R57JLziFZCI7V74xXrivAq99sd1uEHrE6M9l2OdwWIEp85hBWp9YPt2QclFK32y2Ixhm8EDvtCfdVRs9c8+xi145tVZ96pWfudaJ5Frx0Re+dvY7rThln+3EiLfbzqIhM7eazfBG5QkS+b49omt7Iky7XA9Ikju5sVNgZrl7SvmHojUuKRvIc/gr81jQQK4AqIAcYBxQBjwHP2CphCuDxdu8p7pm9LvJOvRSrPV23wjZWO+LJ+Dx4baGgvfX2r6AXaQ3pJcDFIlIAzAAGAU3AaqXUWtul0yQcz/TQPB7e0sRPdzYqGcJebf4Ap907320xuuXEuz+y/RhWcw71wEf2iqLpVXhfP0SFl0ZfeUX5djfkwCsJ6Z6orG1mZ02z22J0S50DlQQsrwSnSS66jffqHrvGJjo3OW+r/+4rGS/a4u0CfE6hjYNF4ltDOnFypAzaSNlGuPLdD85d74IkycmFf/vCbRE8gTYOFvG6gu+sD7ozZl45D20b7CPctX3ss02OH9crbS0WvLg+htNEWib03/TgHSqlzku4RJqYSLYHMcnETSqS8dp6KWcD3snbuEmkhPSfHZFCk3B023YWL13vcIrNS/J5Ba8ZJK8RaSird8dyJRG6EWo0B/DKnIGeSrjosJLFoawiMg74IzAJYxIcAEqp0TbJpbEJr7R5j4iRkoQbKurEffeGyo+CHq7JC19v6/7DXoLVhPTjwN8AH3AS8BTwtF1CaaLH6sPf5k+6R1iTAHRYKTr06C7rxiFXKTUXEKXUFqXUbcDJ9onlPfRQVk2y4FbYRnuDqYXVkt3N5hKh683lO7cDZfaJpUl1vBLe6i3o6x2mk6Y7bT1i1XO4AcgDrgMOAy4DLrdLKE30JJt3kmzyJjv6emuixWptpa/Nl/XAj+0Tx7uk8sOVwqemSRJceb60N9UjkSbB3a+UuqG7yXB6EpwmVlp8AbdF6EI8E5/0pClNqhHJcwiOSNKT4TxOPDHl8t31iRMkifFyFU6Ib2lVnXPoSqsHOyheosecg1JqkflyulJqfugfMN1+8TRO8NQXW9wWwRN4ve8fz2RK7dhoosVqQjpc8vlHCZTD8+iel0aT3OhKBdERKedwKfA9YJSIvBnyUSGw107BUgnda0sOvLzyl0bjNJFyDp8DO4F+wD0h2+uAZXYJpdG4QWOrP+bv6g6AJtWIVHhvC7AFOMoZcTQaTXfEk5D2Oql5VsmNpZyDiFwgIutFpEZEakWkTkRq7RZOo9FoNO5gNSF9F3CeUqpYKVWklCpUShXZKVgqoRNhmkSQyu1owSadwvQaVo1DpVJqdaIOKiJ3i8gaEVkmIq+JSEnIZzeLSLmIrBWRMxJ1TI1G413+78utth9D54Wiw6pxWCgiL4jIpWaI6QIRuSCO484GpiilDgbWATcDiMgk4BJgMnAm8FcRSY/jOJoUweszkL0tnUYTPVarshYBjcDpIdsU8GosB1VKfRDy9kvgIvP1LOB5pVQLsElEyoGZwBexHKc3oadhpD6pnJDWeA+rhffsLLZ3BfCC+XoIhrEIUmFuS2qc6PRqlaHRaBKJ1dFK40VkroisMN8fLCL/E+E7c0RkRZi/WSH73IKxutwzwU1hfiqs3hORK0VkoYgsrKqqsnIarvGVTrbFjcejShp0FYFUw2rO4R8YeYE2AKXUMozcQLcopU5VSk0J8/cGgIhcDpwDfF8dCChXAMNCfmYosKOb339UKTVDKTWjf//+Fk/DHa54YqHbImh6OU7kbLQBTy2sGoc8pdSCTtt8sR5URM4Efo0xPLYx5KM3gUtEJFtERgHjgM7H1fRCvK53tGL0PvoWRYfVhPQeERmDeX1F5CKMshqx8jCQDcwWwxf9Uil1lVJqpYi8CKzCMD7XKKVir2mQQMTjKV9vS6dxG9ExH02UWDUO1wCPAhNEZDuwCfh+rAdVSo3t4bPfA7+P9bc1Gk1XvD4UWOM9rBoHpZQ6VUTygTSlVJ0Z9uk1eH12qrelix9DuXm39+v19qHRRIvVnMMrAEqpBqVUnbntZXtE0mg0icYJ06UjV6lFpPUcJmDMVi7uNCO6CMixUzCNRpM46ppjHj9iGR25Si0ihZUOwhhuWgKcG7K9DvgPu4TSRE+qd9q03tHEi867REek9RzeAN4QkaOUUrqEhUaj0fQSrOYctpnVU3eLSKWIvCIiQ22VzGN4fSir10n1XluKn56mF2LVODyOMUFtMEato3+b2zQaR9DKV6NxFqvGoUwp9bhSymf+PQF4u2ZFL0PrTo3b6NFKqYVV41AlIj8QkXTz7wdAr6om5/MH3BahV6PnEWg0zmLVOFwBXAzswiibcZG5rdcQ8Lhu0p02jdt4PfTncfE8h9UZ0o1KqfNslcTj7KptdlsEjYfRikeTavToOYjIuSJSBSwXkQoROdohuTQpRry9Sq/3Sp3AF9ChTY1zRAor/R44Tik1CLgQ+KP9ImlioVXnRFIer4c2t+1rjLyTJmmIZBx8Sqk1AEqpr4BC+0XSxILf65pDk/LoNphaRMo5lInIjd29V0rda49YmmjRz6XGbXToL7WIZBz+QUdvofN7jUfw63i0u2jN6PmkvL5F0RGpttLtTgmiiY9UTznoB9v7pHqJlN6G1XkOGo8T8PiD6W3pNBpNZ7RxSBECOumg0WgSiDYOKUKqmwZdPkOjcRZLxkFEBojIv0TkXfP9JBH5ib2iaaLB62GlVEdffX0NUg2rnsMTwPsYJbsB1gE32CGQV5kxoo/bIvSIDitp3MbzCWmPi+c1rBqHfkqpF4EAgFLKB/htk0oTNX6PP5jxKg6Pn56G1J8E53njl2CsGocGESnFtL0iciRQY5tUHsTrteoz03X6yE16md4Iy4CiHLdF0CQQqxrlRoyV4MaIyGfAU8C1tkmliZqZo/q6LYKml3PcOL3+VyphqWS3UmqxiJwAHISxdMBapVSbrZJpoiLV17jWHXMYWJSjS8e7iFLejyAkkmhiETOBacChwKUi8kN7RPImZdpl1rjMWVMHuS2CphdhyXMQkaeBMcASDiSiFUZ4qVcwvqyQt9npthjdkurzAHpbMjAcOq0UH6n+jCQaqyvBzQAmKf2EamIk1RuOVjypb8BT++y6YrUvsgIYaKcgGk0yk+J6UdMLsTzPAVglIu+LyJvBv3gPLiL/JSJKRPqZ70VEHhSRchFZJiKHxnsMTWrgdd3rdfk0mmixGla6LdEHFpFhwGnA1pDN3wLGmX9HAH8z/2s0Go2rGGGz3jNcyepQ1vk2HPs+4FfAGyHbZgFPmbmNL0WkREQGKaW8mwnWaNBhJU3qYbXw3pEi8rWI1ItIq4j4RaQ21oOKyHnAdqXU0k4fDQG2hbyvMLeF+40rRWShiCysqqqKVRRNkuB15asT0t7H623Ia1gNKz0MXAK8hDFy6YcYoZ9uEZE5hE9i3wL8Bjg93NfCbAt7S5VSjwKPAsyYMUPfdo9fgZR/MFP9/DS9DqvGAaVUuYikK6X8wOMi8nmE/U8Nt11EpgKjgKViTDccCiwWkZkYnsKwkN2HAjusyqhJYTyufJ0Qz+sG1uPixU2qn19nrBqHRhHJApaIyF3ATiA/lgMqpZYDZcH3IrIZmKGU2mOOgPqFiDyPkYiu8Uq+QYcNNBpNb8LqUNbLzH1/ATRg9O4vtEGed4CNQDnwD+BqG46h0SSc3jOGRdNbiOg5iEg68Hul1A+AZuD2RAqglBoZ8loB1yTy9zWpgec9N20dUh6vh/USTUTPwcwx9DfDShqNRpOU9DLdHjdWcw6bgc/MnEBDcKNS6l47hNKkHp7v+QMjSvPYsrcxpu+mesn0VCDVaz8lGqvGYYf5lwYU2ieOJla83uzjfS6deK69vppeqt9ju4lXPCc6OMP75rG1OrYOSqKxOkP6dgARKTLeqjpbpdJETcDjT6bHxYubNAccB6/fY6+bL89fPpxpR1axOkN6hogsB5YBy0VkqYgcZq9ommjw+truyRBW8nrYwePieb8Nev0C4q1raDWs9BhwtVLqEwARORZ4HDjYLsE00eH1XmW8jd7bZ+fM8pFev8d+L2m2MMQdVnLg9Lx0j60GWeuChgFAKfUp0KtCSx66Z2EJeP3B9PoFxPsGyEuKIxxavvjxkog9eg4h6yksEJG/A89hPEPfBT6yVzRv4fWG5XHbkIBem7dPMM0B18EfsP0QceHxW+R5+cBbeiZSWOmeTu//X8hr75yFA3jdZfZSowqH8rhiA+Jq0Rlp9o908rqB9Hob9Lp84C0ZezQOSqmTnBLE6/g9dNPC4XXFkQwJ6XjIzrTfOHhJcYTD4/2nJPEc3JbgAJYS0iJSglGme2Tod5RS19kjlveIp2HNmj44cYJ0Q6qHHBypehrHd/vm219AIB7FccGhYZdFSSjx5L365GUmUJLkxUudPKujld4BvgSWAx5XQ/YQT1gpJyM9gZKEJ55eZb8CJxSbdxp9d6Ry6LA4137lG498Q/rkJlCS8Og2GB1WjUOOUupGWyXxOF66aeGIp+E7kUz1+OUD4rzHDpyf13VbPJfPifIjcQ+ndmQoa3T7F+VYXpInaqwGSp8Wkf8QkUEi0jf4Z5tUHiQed8+JeHs8DTfdgWmZvkB8DmdvG2Mejnjk09cP/HG2QSeI9hrWNvtYsm2/LbJYNQ6twN3AF8Ai82+hLRJ5FK8npON5MB0xDn5vXz+Iz3NwogPgde/V6x0on8evH8RmxN9fuSvxgmA9rHQjMFYptccWKZKAeBK+TtiVeBRHhgPGoc3rGXO83/P1uHieDx3G20FxwoDF0gazbCoYafVXVwLeKBXoEt6fgey2BD0Tb6/NiQezzePejdeNV1yelwOnlgyeQ0zGIcMe42DVc/BjrB/9IdAS3NibhrLGE1Zyokm2xRFPdUS+JPAcWn1xXEMHLqLXjVe8eSW7SYacQyz32K6ktFWT8zrwe+BzDuQcFtkikUeJRbn1yctkSIn9Q/QgdsXmxBwMiN2lv+zIEcYLR5Rv9Nfwx8eMTLwg3dAag3wnHtTfkWGsEHsbPHVimTOeQ4xtMHiP7ZbRH1AxeV8lefYMRbe6nsOTIpILDFdKrbVFEo/T0hZ9wy/MyXQsidjS6cHMSk+zpEycKh8fa6+ywMaheqEEAiqmsENupjGHxYm73NLmj/o7fUzF4cTkqs5tMDNdPOXtxBpW6mOT8u1MrMbVLvmsrudwLrAEeM98P91cMrTX0OKL/sEsNBWbE72iLg0rCq3v5ZBIqQMzjyG2Xjk4pzggNhlL8jIdKScOXdugWD6wMwI2x2BcAXIcKI0CsRuHEptml1s969uAmcB+AKXUEmCULRJ5lBZfgDH986Ny0e26aeFo9QVIE7jrImOJDcuPpUOaI+jSP/vTI6L6Xm6WMz3zoOL90dEjo/pev0IHjYMvwCkTytq9FSs4abxafAEKczK44dRxgLU2GBxG7YR/0WQahxd/dhQQ/RocdsvY4jfku+P8KVF9z67SLVaNg08pVdNpm3f8RQdo9QUYWJzD8L55lr/TryAbcGakTas/QHZGensYy6neolUaW33AAWVvFSdmzsKBXtuofvlRfa9/QQ7gnHeYnZkWVagyWLPIkbCXL0B2Rlr7tbDSBvvkZTnWVptaDeWbZ7bBTIuVdJ1ug9GOPnI1rASsEJHvAekiMk5EHsJITvcamn1+stLT2pOWZ08dBMAnv+q+cG1pfrZjDb+5zU925gH5Tp5QRlFOBm9de2y33wl6QU4YrwbTOBRkH8ghnHRQf/78nWm2H9sKjS0dFQfAD48awQ+OHM6JB/Xv9nsDirJtly1Is89PdkZ6u5dz78XTGFdWwMtXHdXtd0ryshzLK7W0GfIF2+CsaUNITxPm3Hh8t98Jhg2dyIkEPYcc0/PKSBfOmzaYp66YafuxrdBs5jWzQ4zDT44dxZXHj2bqkOJuvxdth8sqVo3DtcBkjGGszwG1wA22SORR6pt9FOZktj+YN5w6jtW/O5OinAOho8uPGtHhO/0LTcXhQLetrtlHUU5m+4SYwcW5LLvtDKaENKpLZw7v8J3SfOcUR0NQ+YYYh8d/PJOLDhvK5MFFAPznaeO7/b7duqO2uQ0wBhEE+d2sKdxx/tT22lPnTes6sqvIQQNr3OOM9kKJZx88iNk3ntCetJ8wsJAjRnWsauNEtdggtc0+CnMy2nuy4wcWsuEPZzG2rLB9nwsPHdpFPqfaYNBzCCrfzPQ0Hrz0EI4f3799VOFVJ4xxSJquBNtgUUjo+rfnTOI3Z03k0OElAHznsKFhv2sHVkcrNQK3mH+9kpqmNopyM7hj1hTueHs1w0vzyKPwQA0AACAASURBVM5IJzcrnf85eyJnTB7Ilxv3AjCmfz5VdS3Mmj6YZxdscUS+2qY2CnMyOOvgQSzeuo+fn3igkb929dEU5mTQ3BbguQVb27cfP74/tU1tjsjX0GJ4DsW5mZw7bXCHEtL3f3c62/c3ceJBZdwze1379iEluY55XnXNhnxFORmU5GV2MPp3XjCVlxZVcPWJY3hz6Y4O33NKsSmlzHucySM/OIwP1+4m26z2O2FgEU9eMZMZI/rwyPwNfLWpuv17Q81qp06Eveqa2yjKzeSCQ4ewfX8TF884oMhevfpo+hdks31/E68srmjffvSYUpZv7xyxtodG0zgMLM7huHH9OhiCZ356BHXNPiYPLuKR+Rvat7d38LDfuwk+i0U5mdz8rQntYWmA/zzjIE6eOIATxvfnpUUV3f1EQom0TGiPI5KUUuclVhxvopSittlHcW4mR4/txzvXH9fh858eNxo4UOPk2LH9uH3WgaSSE/He2uY2inIMpXbXRR1DNYcM7wPACvMhzEpP49Wrj2bCwEJ+9fIyRxRHQ4sPEcjLTOehSw/p8Nm4AYWMG1DYYduCW04hLyuDf5vK2O6eeV1Ir+2b357W4ZqUFeVwzUljO+z/3g3HUZLrXK+8odVPQEFRbgYzRvZlxsiOHsIJ443Q1/b9TQDcft5kDh5azIjSfMcGHdQ2+xhSkktpQTa3nTe5w2eHmm1whynf0D65/PX7hzJlcDFX/Z8zU6ZqzA5UZnoaT/+k48CIkWauKbQSwlvXHsvgklxeXeyMMg52UIpzM/hZJw+mKCez/R4Heee642wLKUFkz+EoYBtGKOkrnOsoeYrGVj/+gOrQmwxH0NIfN+7ATXQqmVXb5GNkv56T5YNN1/msqQMPhJscuqN7Glrpm5dFmoU6TpnpQllhjgNSHSD4YBbmZCAi3XosZ0wewPsrKzloQCEiwu7aZsCBsFdIr7Inpgwu5tXF25kypLi9U+AUtU1tTBxU2OM+wTZ42qQBHDzUCJU45R1WN7RGHBodbJ/D+uZ2CMk6QbjQZjjGDyhgXWU9EwcV2mr4IxmHgcBpwKXA94C3geeUUittk8iD7GtsBSIPTZ01fTBjywraY+hOsrehlenDSnrcp29+Fu9cdxxjywo6bHfCc9hd29LBRe+Oj//7pA69Iad6I3sbjKowfSIoj7sunMZt5/kOPJQOKjaI3AYvP3okR48tZcLAA21QsL+qsFKKfY2tEb2pYX3zePu6YzvIZ3zfTukMqhtaI95fgE9/fVIHBR2813ZPaK2uN+5xpOHyL1x5FI1tfts9wh4T0kopv1LqPaXU5cCRQDnwkYhcG++BReRaEVkrIitF5K6Q7TeLSLn52RnxHicR7Nhv9A4HFfdcCkNEmDKkuMNNy0yXuGr2WKHF52dPfUt7r6wnJg0u6jBULiczPaYJftFSVW/NOAwvzQu7n90TzXfsb6YgOyNiz7w4LzNsO7BbtwXDMZHucXqadFG8uVnpNLfae49rmtpobPUzuCSyxzd5cHGHMvFZGc60wT31LZTmR26DQ/vkdVDQ+WZnpSnGSXRW2VHTTGl+Vvtoqu7ok5/lSFmeiAlpEckGzsbwHkYCDwKvxnNQETkJmAUcrJRqEZEyc/sk4BKMkVGDgTkiMl4pZX/L6YHggxnLUob52RnUm8lYu9hVYxgvKw9mZwockA+gsqaZMf1Lo/5e0IuIdXarVXbWNDGoOPrrF3yQYyltEQ1WjUM4CrIz2ocS20WwAxWbfOnUt9h7/ZRSbK1u5Kgx0bfB4Ai7Bptl3LG/KabrZxeREtJPAlOAd4HblVIrEnTcnwN3KqVaAJRSu83ts4Dnze2bRKQcY2b2Fwk6bkxsrTaqlQ+O4DmEIz8ro30CmF0E5YvJeGUZo5h8/gAZNtWFr2tuY1dtM2P6F0TeuRPBEiT1zfZewy17G9tH9kRDQZYhX63N8m2tbiInMy2mciJ5WentI3XsYmt1A0BMPVonnpGq+hYaW/2MLI1ukiMc8BwabO5Eba1u5KABPedsnCSSNrgMGA9cD3wuIrXmX52I1MZx3PHAcSLylYjMF5HDze1DMBLgQSrMbV0QkStFZKGILKyqqopDlMis3FHD6H75MY0MMDwHex/M1TuNW9E5nGCF/Gyz4duoPDZUGYpjXFn0xiHfVL52ejctPj/lu+uZOCj665eWJob3ZbNxWL2zloMGFsUUZ87PzrBdsa3cUUuawPgYlFt+dgaNrX5b10xZt6seIKYOSp7ZBu30vhpafGze2xBTG7SLHj0HpVTMXUkRmYOR0O7MLeZx+2DkMQ4HXhSR0YRP74VtMUqpR4FHAWbMmGFryHfF9loOGd5zsrc7CrLTbX8wV2yvZWBRTkwTnvJNl7mx1Wdbaeel5hq3sTT84ASvOhuV7/rKenwBxaQYBxIYoTn75osEAopVO2s5a2q4xykyeVnpVNW1RN4xDlbuqGVM/4IYO1DBDoov4kidWFmybR8AU4dGPwIpKF+jjZ281TtrUQpXBrN0h231kJVSp3b3mYj8HHhVGbNKFohIAOiH4SkMC9l1KLAjzE84xta9jWzf38RPjo2tzqDdvbZAQPH5hj0cM7ZfTN/Pb4+n2ifjFxv2MqQkl2FR1KUKUphtKAs7PYdPy43Vb2eM6Bthz/AU5Nibt1mxo4aapjaOGBV9vBwM78vOXm+rL8CCTdWcO21QTN/PD4np22UcFmzex9iygpg6QE54Dp+V70WEmDuhduBMLdquvA6cDCAi44EsYA/wJnCJiGSLyChgHLDAJRkBmLemEoBTJpbF9P18m5OBK3bUsKe+tcsEGasUmL0iu0JfzW1+Pivfw7ExGq+g52Cn8Zq3ZjcTBxUxMIaENBieg52ezYdrqhCB48bFdg3zstNtTaYu3FxNfYuPkw6K7RkJ1tuyy8DWNbfxxYY9nDwh1mfY9BxsDL3OW7ubaUNLKC1wrlZXJNwyDo8Bo0VkBfA8cLkyWAm8CKzCWDviGjdHKimleHlxBRMGFjIihkQWdEz42sHLiyrIykiL+cEM9ooabXow56yupK7Fx3kxrjiX32687JFvy94GFmyq5szJsYVswEia22UcAgHFK4srmDmyb8yKoyA7k/pmn23lH15eVEFBdkbs3mvWgdCmHby9bCdtfsUZMd7joGcTnEWfaMp317F0237OnBJ7G7QDZ5bZ6oRSqhX4QTef/R5jSVLXWbx1Hyu21/K7WZMj79wNB3q+forzEmuLa5raeG3xds6aMtDS5J5wBHttdTYoX6UU//xkE0NKcjlydGwhkeyMdLIy0mxTvk9/sYX0NOG7hw+LvHM3FGRnsNMcTpxoPlq3m63Vjfzn6d0XJYxEca5RMLK5LZDwcgtVdS28tWwnl8wc1q5EoyX4PTuS+oGA4vHPNjNxUFF78bpoKczOID1NqLGpDtmTn28hKz3N0aJ6VnDLc/A8Sinuem8tpflZXHBo7DctGOO0o2H9ff4G6lp8XHl87JUk7ZTv4/V7WLJtP1efNKbDpKdosSvhu7Omiae/3MJ50wbHHFICw3OwS7H9+f11DOuby7emxBbPB3vv8UPz1uNXKupFkkIJzvq2Q75nvtrC2so6fnb86JhnFIsIRTkZtsi3dW8jz3+9lQsPG+KpkBJo49At76/cxVebqrn+1HEd1iCIlhKbHszNexp47LNNnDdtcMyjbCDkwWxMrHzNbX5ue3Mlw/vmcVGcPSK7Yvp/encNShFXrxzMsI0NntdzX29l1c5abjxtfNQLwIRil3FYu6uOZ7/ayqUzhzE6hiGiQYLy7U+wfNuqG/nju2s4blw/ZsUY1gxSnJtJTVNi77FSij+8s5r0NOGGU+Nrg3agjUMYqupa+M1rK5g8uKjLGgjRElS++5taEyEaYPQof/XyMjLT0/jNWRPj+q0C02VOpHwAf/2wnE17Gvj9t6e0l5aOleLczISXFn9vxU5eX7KDn584hqF9oh9FFUpwtFIia+9s2dvA799ezTFjS5k1LexUH8u0t8HGxN3jVl+AX76whKLczLgV2wH5EnePAwHFf7+8lHQR/nThwXHXITKMQ2Lb4OtLtvPeyl1cf8p4BhQ5W2jSCto4dEIpxU2vLKO+xcf9351OZpyzhu1o+E98vpkFm6u59ZxJcYVDwHCZS3IzEyrfN1v38ZePNvDtQ4Z0qFAbK33ys6hOoHy765q5+dXlTB1SzC9OHhv5CxEIeoeJMmD+gOLGF5eSnibcfdE0S5Vse8IOz+GBuetYtbOWP14wtcO6A7GQm5lOVnpaQjsoT3+5hS83VvM/50xMSEmK4rwsahJoXLfvb+LW11dy+Mg+XHn86IT9biLRxqETL3y9jblrdvPrMyd0WWMgForNKpWJcpnXVdbxp/fWcPKEsrjDNUGK8zITJl99i48bXljCwKKcLjX9Y6VPXmbCer1KKW5+ZTmNrX7u++60uI0/HFhtrTpBMj4yfwOLtuzjf2dNSYxiS7BxWLSlmr99tIHvHDY05hFAoYgIxXmZCQttbtnbwJ3vruGE8f25eEbsAw1CSaTnEAgo/uvFpQSU4p7vTI8rH2cn2jiEsGVvA797axVHjynlx3Ek2EJpfzAToDia2/xc99w3FOZkJMRVDlKcm7gH8/Y3V7KtupH7vjs9YTOu++RltZesjpfnFhww/qHLV8ZDcKTYvgTIuHJHDffPWcfZUwfFHScPUpRA49DQ4uPGF5cyuCSXW8+dFPfvBUmU92qEk5aRkS7ceeHUBD4jiUtIP/75Zr7YuJdbz53E8NL4Qpp2oo2DiT+g+M+gK/+d+F35IFkZaeRnpSek4d/13lrW7Krj7oumWSp/bZWSBPWK3l6201xOcywzR8U22zgcffKyqGv2tS9cHysbq+r537dWcezYfnGNrulMsBhevAasxefnxheWUpKXxR3nT0mYYivMziBNEmMc/vetVWytbuTei6cndDZzSV5i2uCTX2xmwaZqfnvOpIgl9qOhODeT2gTMFQl6/qdOLEuYV2MX2jiY/P3jDSzcso/fzZqc8FrpJXlZcYdtPl5XxWOfbeLyo0ZwUowzPbvDkC8+xbarppmbX13G9GElXH/quARJZtA3P/68TZvfSKBmZ6Zxz8WJM/4Q4jnE6R3eO3sdayvruOvCg2OetxKOtDShKAEdgNmrKnn+62387PgxCTX+YIRf431Gtuxt4E/vreGkg/onfM5ASW4W/oCKa1Raqy/ADc8voTA7gz9ekDjP3y60ccAoenXf7HWcNXUg50+Pb2RIOIridJkbWnzc/OpyxpYVcHOco5PCUZybyf6G+B7MW99YQYsvwH0JSOJ3piTPzNvEoXz/Pn8DSytq+MO3pyZ8ZEhfU769cXgOi7ZU8+jHG7nk8GEJN/5gyBiPZ7OvoZWbXlnGpEFF3Hha4oddlsSZV1JK8ZvXlpOZlmaL4u2TAO/wwbnr25P4ifT87aLXG4dAwGhURTmZ/P78xMUoQ+mbn9m+DGUs3PPBOrbvb+LOC6ZGXCUqFvrmZ1HX4ot5Na73Vuzkg1WV/PK08YzqF1uZkZ7oG+eDuXlPAw/OK+fsqYM4a2rsk8m6IzcrnZzMtJhzDq2+AL9+ZTmDi3P5n3MSF8cPpV9hdlyVWe98dw37m9q497vT4ppz0R1987PY29Aac9jm5UUVfFa+l5vOmhD3CL5wBJV5rNdwfWUdj8zfwAWHDuH0BCTxnaDXG4eXFm3jm637ufmsiQl15UMpK8xhd21sjWpZxX6e+HwTPzhyODNGJtaVDzKgKPaGX9/i47dvrGTy4CJ+GmPl2kgEhwPvi8H7Ukrx2zdWkJWeltAEamdK87OpjtH7euLzTZTvrud3sybHNeGyJ/oXZlNVH1sbXLSlmhcWbuMnx46Kac0QK5QVZtPqC1Abw0Sz6oZW7nh7NTNH9uXSw+Obl9Qd/Qtif0aCbTA/O4NbbPD87aJXG4d9Da3c+e4aZo7sy4WHJj6cFKTM7LXF0iv64ztr6Jufxa/OnGCDZAZlhUZPa3cMDf+xTzdRVdfCHedPsW0luaDnEIv39f7KXXyyfg//dbq9E4365GfGlHOorG3mgTnrOXlCGadMHGCDZAb9C7KpiqGD4g8obnltBYOLc7j+lMTmkkIJ9sx310Vfo+qR+Ruoa27jjm9PSWguKZR2zyEGA/vGkh18ubGaX515kOdKZPRErzYO89dVUd/i43fnT7Y1OdS/MJtWfyDqhODn5Xv4YuNerj5xbMSF7+Oh/cGMUnnsb2zlHx9v5PRJAzhkeB87RAOgX0E2aQKVUcrnDyju+WAdY8sKuOyokfYIZ9KvIDsmxXb3+2tp8yv+n41eDRj3uK7FR1OUZadfXVzBml113HL2pJgL61kh1g7K7tpmnvx8M+cfMiSmVeis0jc/izSJ3nNobvNz9/trmTqkmEts8mrsolcbh/MPGcInvzrZNlc5SFlR9A1fKcU9s9cxsCiH7x1hb6Mqaw8rRafc/vHJRupbfdwYZ22iSGSmp9G/MJtdNU1Rfe+tZTtYv7ueG04dZ/tEo0HFOeyKsjLrxqp6Xl1cwWVHjYi5JLxVgh2APVH0fJvb/Nw/Zz3ThhbHvAqdVYKhzWgN7F8+LMcfULZ6NQDpaUK/gujzNs98tZXt+5u46VsTPDvZrTt6tXEAbEledaYshp75/HVVLNqyj1+cPNaWJHQopflGzzwa47W3voXHP9vM2VMH2W5cAQYW5URVFtvnD/DAnPVMGFjIWXFUNLXKoOJc9tS3RpXUf2heOVkZaVx1QuxVda1yIGxj/R4HFduvz5xg+7DL9g5UFM9Ixb5Gnl2wle/MGGa7cQUzbxPF9atrbuMvH5Zz7Nh+Ma914Sa93jg4QVmU8VSlFPfOXsfQPrmOTJQJ9oqieTAfmb+B5ja/Y9UkBxbnUFlr3Ti8vmQHG/c0cMOp422LQ4cS7GRU1li7huW763ljyXYuP2qkI8MaDyRUrV3DoGI7blw/jnZAsRVkZ5CXlR6V8Xp4XjmCcG0C6mNZoX9hdlTy/fOTTVQ3tPLfZxxko1T2oY2DAwQTobssKrfZqypZVlHDdaeMs2XYYDgGFOVYlq+ytpmnvtjCtw8Zytiy2Es1R8Og4lzLnkObP8CDc9czeXARZ0y2L8kbymBzNu5Oi6GvB+euJycz3bGia8GJnRX7rMn31BdbHFds0bTBzXsaeGlRBd87YnhC6k9ZYXBJLtv3W7t+NY1t/OvTTZw5eSDThnlnXeho0MbBAfKzMyjNz2JbdWPEfQMBw2sY1S+fCw6xbwRVZ4b3zWOrBfnAuThvKAOKcqhr9lmaofrKogq2Vjdy42njHZuFGvQcrBiwdZV1/HvZDi4/eqRjo1dK8jIpzMmwdI/rW3z845ONnDKhjIOHOqfYhvXNY+tea23wwbnryUwXrj7R/pBckOF986huaLW0XOhjn22ivsWX8GoBTqKNg0OMKM1j857IDf+dFTtZs6uOG04dZ9vQ0HAML82jYl9jxLWuK/Y18tyCrVx8+DBHi4YNLjGU7/YIPd9WX4CH5pUzbVhJzAvKx0K7fBZ6lg/MXU9eZjpXHudcqWYRsdwBePqLLexvbONaB40/wIi+eWze2xBxyPf6yjpeM0NyZQ6ugzCir9HeI13DmqY2HvtsE2dMHsDEQfbn4+xCGweHGFmaH7FR+QOK+2avY1xZAeccnJiKnFYZWZpHm19F7Pk+NNeI8/7iJGfivEHGmCuNbayq73G/FxZuY/v+Jke9BoC8rAzKCrPZtKehx/3W7Krl7WU7+fExo2ybdNkdI0oj98wbWw2v4YTx/ZnucDhkRGkedc2+iKVm7p9jGNefOZDID2WYaRwiRQCe/Hwzdc0+rj05eb0G0MbBMYaX5rGjponmtu5Hs7yxZDsbqhq48bTxjg97G97XGO2xeW/3ym3zngZeXuxsnDdIsCzHxh6Ub3Obn7/MK2fGiD4cP8750SFjywoo392z8XpgznoKszP46XH2zCbviWF986jY19TjinX/96WRa7jOYa8BaB9x1FMbXLmjhreX7+Qnx45qnxzpFEFPeUsPBrau2cg1nDpxAFOGFDslmi1o4+AQI0vzUar7XkebP8ADc9czaVBRQhZQiZaR/YyGv7kH5ftAMM57krM9NjDyNgOLctjQg+fw7Fdb2VXb7LjXEGRM/wI2VNV3GxZZsb2Gd1fs4sfHjmovJugko0rzafUHqNgXvg02tfp59OONHDeuH4eNsG9SY3eMNJVvT8bhvtnrKMrJ4CcOhuSCFOVkUpqf1WMbfOqLLdQ0tXHdKc561nagjYNDHDTQmL25ckdt2M9fWljBlr1GEtWJoZedGViUQ0leZrfyrdlVy+vBOG+hO+vdjinL77Zn3tDi468flXPk6L6ODL0Mx9iyAuqafd0Od7xv9jqKczNd8RoAJg82erLd3eNnvtrCnnp3vAaAkf3yyc5IY+X28PIt3rqPOat3c+XxoxO2kFS0TBpc1O31q2/x8c9PNnLSQf0dTeTbhTYODjGurICczDSWVuzv8llTq5/756xjxog+nDLRuSRqKCLC1CHFLKuoCfv53e+tpSA7g587ODqkM1OHlLB6Z23Y0Nzjn21iT32rrTWoIjF5sJF8XLKt6z3+Zus+5q4xFJudpVB6YvzAAjLShOXbu97jmqYDE7YOt6nAYyQy09OYNLiIZWHkU0px57tr6FeQzY+Pcce4AkwZUsy6yrqwkx0fnb+BfY1tjs39sRttHBwiIz2NyYOLWRpGcTzx+WZ217Xw62/ZPxO1J6YNLWFtZV2X+jtfb65m7prdXHXCGFfCIUEOG9GHNr9iRSflsa+hlb9/vJFTJw7gUBtrPEVi6tBisjLSWLi5usN2pRR/fGcNpflZCV2BLlqyM9IZP6CQZWE6KI+Yiu2mb7lnXAEOHlLMiu01XfIiH62rYsGmaq47ZaytNZ4iMWVwMW1+xapO3kNlbTP/+GQT504bnLTzGjqjjYODHD6yL8sqajqs17xjfxMPz1vPKRPKXOuxBZk5qi/+gOLzDXvat7X5A9z6xkoGFGVzhYs9NoBDhxsP3Rcb9nbYfue7a2hs9bs+EzU7I53pQ0v4YmNH+V5fsp0Fm6v57zMOclWxARwztpSvN+3rMF9k054GHvt0E+dPH+x6EnXmqFIaW/0s2rKvfVtzm5/f/XsVI0rzXC9ed+TovogY5W1CufPdNfgCAf779OScDR0ObRwc5PTJA/AFFLNXVwLGhLdbXltOQMFt5012WTo4cnQphdkZvLN8V/u2Rz/eyOqdtdx+3hRys+yt8RSJ0oJsZozow1vLdrZv+7x8Dy8s3MZPjx3Vntdxk9MmDWDF9tr2pGVlbTP/+9Zqpg0t9sSawSdPGECrP8C8NbsBY/j0r19ZRlZGmi2rDEbLCQf1Jys9jXeWH7jH981Zx6Y9Ddxx/hTHKgZ0R2lBNocO78N7K3a1DzyYu7qS177Zzs9PHOvo3B+70cbBQaYPLWF0/3we/dioS/Sn99fw4doqfn3mQe1jqN0kKyON86YP5t9Ld7BpTwPvLt/JPR+s5ZyDB3HmFG+sXnXBoUNZW1nH7FWVlO+u49rnvmFM/3zPzESddchgMtOFB+asp7qhlSufXkRTqz/h61bHysxRfRlZmsffPtpAU6ufW99YwYJN1dx27mRb17uwSkF2BmdOGciLC7exs6aJVxZV8Pf5G7l05jCOG9ffbfEAuHjGUNbsquO9FbtYXlHDDS8sYeKgIsfn/tiNxLosn5eYMWOGWrhwodtiWGLOqkp++tRCcjPTaWrzc+nM4fzh21M8s9j4zpomznrgExpa/LT6A0wfVsKz/3EEeVnuhkOCtPoCnPPQJ2ysakDEWJj+hZ8d2T5JzgvcN3sdD8xdD0BmuvDX7x/GaZOcqfFkhfdW7OKq/1tEZrrQ5lf87ITR3Pwt972GIBur6jnv4c9o9Qdo9QU4anQpT1xxONkZ7nquQVp9Ab79189YtdPIOwwqyuHFq45iaB/3O3jRIiKLlFIzwn7mhnEQkenAI0AO4AOuVkotEENDPgCcBTQCP1JKLY70e8lkHMAorPf+yl0cObqUCw8d4hnDEGR9ZR1Pf7mFISW5XH70SNtLhkfL7rpm/vbRBtJF+MlxoxhU7OyEvEgopXh5UQXlu+v59qFDHClpHi1zVlUyb+1ujh/XjzMdKGkeLSu21/Dcgq2M6pfPD48a6Xo4qTPVDa089ukm0gRHa2QlGi8ahw+A+5RS74rIWcCvlFInmq+vxTAORwAPKKWOiPR7yWYcNBqNxgv0ZBzcMscKCHanioEd5utZwFPK4EugRES8163RaDSaFMetQPINwPsi8mcMA3W0uX0IsC1kvwpz2046ISJXAlcCDB+eXGuzajQajdexzTiIyBwg3BCXW4BTgF8qpV4RkYuBfwGnAuGC72HjXkqpR4FHwQgrJURojUaj0QA2Ggel1KndfSYiTwHXm29fAv5pvq4AQgeDD+VAyEmj0Wg0DuFWzmEHcIL5+mRgvfn6TeCHYnAkUKOU6hJS0mg0Go29uJVz+A/gARHJAJoxcwfAOxgjlcoxhrL+2B3xNBqNpnfjinFQSn0KHBZmuwKucV4ijUaj0YTirZklGo1Go/EEKVE+Q0SqgC0xfr0fsCfiXqmFPufegT7n3kE85zxCKRW2aFVKGId4EJGF3c0QTFX0OfcO9Dn3Duw6Zx1W0mg0Gk0XtHHQaDQaTRe0cTBnWfcy9Dn3DvQ59w5sOeden3PQaDQaTVe056DRaDSaLmjjoNFoNJou9GrjICJnishaESkXkZvcliceROQxEdktIitCtvUVkdkist7838fcLiLyoHney0Tk0JDvXG7uv15ELnfjXKwgIsNE5EMRWS0iK0XkenN7Kp9zjogsEJGl5jnfbm4fJSJfmfK/ICJZ5vZs8325+fnIkN+62dy+VkTOcOeMrCMi6SLyjYi8Zb5PPlFm8gAABcpJREFU6XMWkc0islxElojIQnObs21bKdUr/4B0YAMwGsgClgKT3JYrjvM5HjgUWBGy7S7gJvP1TcCfzNdnAe9ilEg/EvjK3N4X2Gj+72O+7uP2uXVzvoOAQ83XhcA6YFKKn7MABebrTOAr81xeBC4xtz8C/Nx8fTXwiPn6EuAF8/Uks71nA6PM5yDd7fOLcO43As8Cb5nvU/qcgc1Av07bHG3bvdlzmAmUK6U2KqVagecxVqJLSpRSHwPVnTbPAp40Xz8JnB+yPdyKe2cAs5VS1UqpfcBs4Ez7pY8epdROZa4vrpSqA1ZjLAyVyueslFL15ttM809hVDZ+2dze+ZyD1+Jl4BQREXP780qpFqXUJoxClzMdOIWYEJGhwNmYpf3Nc0jpc+4GR9t2bzYO3a06l0oMUGbJc/N/mbm9u3NPymtihg4OwehJp/Q5m+GVJcBujId9A7BfKeUzdwmVv/3czM9rgFKS7JyB+4FfAQHzfSmpf84K+EBEFomx6iU43LbdKtntBSyvOpeCdHfuSXdNRKQAeAW4QSlVa3QSw+8aZlvSnbNSyg9MF5ES4DVgYrjdzP9Jf84icg6wWym1SERODG4Os2vKnLPJMUqpHSJSBswWkTU97GvLOfdmz6E3rDpXabqXmP93m9u7O/ekuiYikolhGJ5RSr1qbk7pcw6ilNoPfIQRYy4RY20U6Ch/+7mZnxdjhB6T6ZyPAc4Tkc0Yod+TMTyJVD5nlFI7zP+7MToBM3G4bfdm4/A1MM4c9ZCFkbx602WZEs2bQHCEwuXAGyHbw6249z5wuoj0MUdCnG5u8xxmHPlfwGql1L0hH6XyOfc3PQZEJBdj3fXVwIfAReZunc85eC0uAuYpI1P5JnCJObJnFDAOWODMWUSHUupmpdRQpdRIjGd0nlLq+6TwOYtIvogUBl9jtMkVON223c7Ku/mHkeVfhxG3vcVteeI8l+eAnUAbRo/hJxix1rkYy7DOBfqa+wrwF/O8lwMzQn7nCoxkXTnwY7fPq4fzPRbDRV4GLDH/zkrxcz4Y+MY85xXAreb20RiKrhxjTfZsc3uO+b7c/Hx0yG/dYl6LtcC33D43i+d/IgdGK6XsOZvnttT8WxnUTU63bV0+Q6PRaDRd6M1hJY1Go9F0gzYOGo1Go+mCNg4ajUaj6YI2DhqNRqPpgjYOGo1Go+lCb54hrenliEhwaCDAQMAPVJnvG5VSR9twzEOAa5RSP43zd34BNCilHk+MZBpNR/RQVo0GEJHbgHql1J9tPs5LwB1KqaVx/k4e8JlS6pDESKbRdESHlTSaMIhIvfn/RBGZLyIvisg6EblTRL4vxroKy0VkjLlffxF5RUS+Nv+OCfObhcDBQcMgIreJyJMi8oFZv/8CEbnL/N33zPIgmMdcZdbq/zOAUqoR2CwiyVZZVJMkaOOg0URmGnA9MBW4DBivlJqJUUL6WnOfB4D7lFKHAxean3VmBsbM5lDGYJSjngX8H/ChUmoq0AScLSJ9gW8Dk5VSBwN3hHx3IXBc/Ken0XRF5xw0msh8rcxSySKyAfjA3L4cOMl8fSowKaQqbJGIFCpjrYkggziQ0wjyrlKqTUSWYyxA9V7Ib48E3gKagX+KyNvm+yC7gQlxnptGExZtHDSayLSEvA6EvA9w4BlKA45SSjX18DtNGLV/uvy2UiogIm3qQBIwAGQopXxm6OgUjMJzv8CoTIr5Wz0dT6OJGR1W0mgSwwcYihsAEZkeZp/VwNhoftRcr6JYKfUOcAMQ+rvj6Rqm0mgSgjYOGk1iuA6YYSaNVwFXdd5BKbUGKA6WY7ZIIfCWiCwD5gO/DPnsGGBOHDJrNN2ih7JqNA4iIr8E6pRS4RLW0fzOIcCNSqnLEiOZRtMR7TloNM7yNzrmMGKlH/DbBPyORhMW7TloNBqNpgvac9BoNBpNF7Rx0Gg0Gk0XtHHQaDQaTRe0cdBoNBpNF7Rx0Gg0Gk0X/j8lSzy4SVKMEgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(s[4999:, 0])\n",
    "plt.xlabel('Time (ms)')\n",
    "plt.ylabel('Membrane Potential (mV)')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [],
   "source": [
    "als = anmodel.analysis.WaveCheck()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'SWS'"
      ]
     },
     "execution_count": 38,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "als.pattern(s[4999:, 0]).name"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "san = anmodel.models.SANmodel()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "san.set_sws_params()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [],
   "source": [
    "s, info = san.run_odeint()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x1a1cfeb9d0>]"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO2deZhcVZn/P6e36j2dfU86e0hYQww7BAi7Eh1QQUcUHSMI/hzH0SE6o+OCg6gwIghGQQZUFMQFCUsSEgghJKETsq+d7k7Se6f3rbprOb8/qqpTdKrurap7b91bVefzPP101b236j2n7jnf8573LFdIKVEoFApFepJldwIUCoVCYR1K5BUKhSKNUSKvUCgUaYwSeYVCoUhjlMgrFApFGpNjdwLCGTNmjCwvL7c7GQqFQpFSbN++/aSUcmykc44S+fLycioqKuxOhkKhUKQUQohj0c6pcI1CoVCkMUrkFQqFIo1RIq9QKBRpjBJ5hUKhSGOUyCsUCkUao0ReoVAo0hgl8gqFQpHGKJFXOJ6W7gFe29todzIUipREibxClxNtfRxr7bXN/mee3Mpdv9tO36DXFvt1Hf3MXLmaAw1dtthXKIygRD6FGPT6ae0ZSLrdyx7cwBU/eTPpdkOcaOsDwG/T823W7mvEL+Frf9qJz65EAK/uaeBwU7dt9hUBdhxvp9vtsTsZMaNEPg62VLVSUdNmm/2v/Wkn5/9wnW32pZRsrWolU58mdrCxm8c2VNpie9eJDu7+/Q6ufXijLfYBOvs9NHe5ae8d5MevHbSlwWvudlN+32rK71vNiba+pJbFtt5BHlp7mH/65WbO+u811Jy0r3cbD0rk4+C2VVu49Yl3k2pTSkn5fav51VtHWb2nIam2h/PSrno+uWoLL1TU2mLf6/PbYjecoy09tthd/tg7Q687+gZtCZ9d+sB6lvzoDe763XYef/Mo6w40JT0NP193ZOj1ZQ9u4PG3jibN9qIfrOWRN07ZX/rTN7njqW1Js58oSuRjYMPBZh5ac8jWNPzPqwdttQ+nwibH2uzxYM79/lq2H2unK4W6ylaw7KGNtoTPugcCYyJbqwO92W//dW/SG96/vV/3gfd2D8hvPNxiq/1YcNQulE7lzqffszsJGY0QYuj1LY9v5uwpI3jp3kttTJG9nLRhXKZ34PRB75M9A7xb1cplcyLucGs6W6pa6R30JcVWOqE8eYXjGR533V3baVNKIEOHI7jwf96w1X5H3yC3rdpy2vFMvR/xoEReBztnU4CzCvGbhwJdUyelyQ7K71vNnb+1Pxa79Ccbkmar223P9NUQPRF6EorYUCKvw7Pv1tidBMdQcazd7iTYzq7aDgA2HLI/FlvT2md3EhQpgBJ5HZq6kx//VHyQ8Ji83RxTwvoBMr1XlwookU9B7J6n7iDNTQpOamQylcZOt91JSFmUyOsQTU8PNXazZp/107ec6Cgp7y1AR9+g3UmwnWS1fyue3R7xuHRkDXEWSuQT5Lr/3Ri14CnMJVLPZV+9fTNsQix76C27k5AxqIHXxDEk8kKInwghDgohdgsh/iqEKAs7t1IIUSmEOCSEuM54Uu1B9dSdyQ9e3m93EjjZozz5pPXqlMOeMEY9+bXAmVLKs4HDwEoAIcQC4DZgIXA98EshRLZBW4phRFqgko6omLhCkTiGRF5KuUZKGVKaLcCU4OvlwB+llANSymqgElhixJZd2B1/jhSqkBJeqDjBwu++TmWz2pVQkQHY3M7budW2UcyMyX8eeDX4ejJwIuxcbfDYaQghVgghKoQQFS0t9s89ThXeONAMwJEmezbMUiggeeFMu/ty++tT91kCunvXCCHWARMinPq2lPLvwWu+DXiB34c+FuH6iD6xlHIVsApg8eLFjou8qUjB6bSrWSWKIMnq6UYzY3dPOxXQFXkp5TKt80KIzwIfBq6Wp2ILtcDUsMumAPWJJtJpzFi5Omm29MqwHWX8uW0n+MpVc5hUVmCD9QBbqtpwe3zk51o/1KMaekUqY3R2zfXAfwA3SynDlwK+BNwmhHAJIWYAcwD7N/tIgEiegt3eg8R+4WnodLOlqtW2R/IB/HB1cmbY2H2/FfaHa1IZozH5R4ESYK0QYqcQ4gkAKeU+4HlgP/AacI+UUu0RagF2CVBTl5vbVm3h68/vsicBQGOn2nJCkRzsdqqMYGg/eSnlbI1z9wP3G/l+J5DKN9dKQtM3DzWq2T0K67G7M5XKvTm14lUHu2+unv1MaIQyIIsKHVQZSBwl8inIz9cdxu0JRL/sboQygUxoSFMVVf71UY//S0EeWV9pdxLwB2tXMupYNBsdfYN0uz2U5OcmIRUKO7Fby1O5oVeevA6pfHOt5D9e3GN3Eqg41s7iH66z3I7yFqNj909jt/1UQIm8DnZXcKdvpZqMNlDLxoDXn4QUKOxG+VqJo0ReodBB9eaik6yfJto9ULdGHyXyKY7TPX2Fwgyi9ahV6ddHibwOyovTRlUyRTJwcjmz+3GceiiR18Hu+2e3fYVCi2QVTyf7Wn9674T+RTaiRD7FUY2AIhOI1qNOnhcdvZmpdvhe80rkFYZwsoelsJ5k3X/7nRnbE5AwSuQVCkXKoh4NqY8SeYUhkuHfdGfIs2xTEbv9W6cPejoBJfIKhQ5avuIz79YkKRWZTTSH/WBjd5KEPnV7DErkUxzlx9jLd/6+z1b73/vHPga86f+oBi0d99tcCYTDGwAl8g7H6b1RZxdvc3DyLfjtOzU8X1FrdzIyGqcvSFQirzCEs4t3ZpAJcWk1vpo4SuQVihQnAzTe0ahwjUKR4ji7CmcGTm7I0jpcI4T4gRBid/Ah3muEEJOCx4UQ4hEhRGXw/CJzkpt56BWgTOiqKxSKxDHqyf9ESnm2lPJc4GXgO8HjNwBzgn8rgMcN2lE4FOXlZjaVzT30D1o/u0crJm+3o5PW4RopZVfY2yJOjcMtB56RAbYAZUKIiUZsKZyJ6kfYj50i94OX9/OV53bYZj9ZaDYyDq8Fhp/xKoS4H7gD6ASuDB6eDIRvzVYbPNYQ4fMrCHj7TJs2zWhyFApFktl8tNXuJFhOKkdFdT15IcQ6IcTeCH/LAaSU35ZSTgV+D9wb+liEr4r4M0kpV0kpF0spF48dOzbRfKQtqVy4FMkh04uI3fl3erhG15OXUi6L8bv+AKwGvkvAc58adm4KUB936hQKhePJBEdEK1yz80R78hKSAEZn18wJe3szcDD4+iXgjuAsmwuBTinlaaEahUKhSHW2VLXR2jNgdzKiYjQm/4AQYh7gB44BdwWPvwLcCFQCfcCdBu0oFPahlltqYvfAoxN6Ev0e5+4fZEjkpZS3RDkugXuMfLcigAPKr8IJKqKBw5OnsBm14lWhUBjC7kbm2offsjcBDkeJvA52d0UVCqdjdw2pae2zOQXORol8imO3F6VQKJz9GEIl8jrYPQfW7iXbCtTAqyKlUSKvUKQ4trsBtifAelK5mVcir4OKySsU2qg64uwetxJ5h+PcoqNQJA/bw6a2WjeGEnkd7C5cetjtRTnZg1EokoUaeFUoFJZhd0ObDPN2OzPOlXB9lMjrYHfh0sNuR9rJHowiOTi7hiiUyDscu0VcoXACTg+bOhkl8jo4vXDZ7UjbHSpIBs4uAZmB03vUTkaJvA5OL1wZoLG2o35ibTKhoU9llMgrFApDJEPiHd+jtjsBGiiRdzoOd5LUwKv9KEdaoYUSeR2c7kGoCm49zi4B9pOM38fusGkqOzNK5HWwu3A5HRWPVagS4GyUyCsUBlENnfXY3aNO5XtsisgLIf5dCCGFEGOC74UQ4hEhRKUQYrcQYpEZdjIR1ZOwH6ffAbvLSDL0L4WjJbZjWOSFEFOBa4DjYYdvAOYE/1YAjxu1Yxd2exAK55PCTl7KYPdvnOkx+YeBb/JBh2c58IwMsAUoE0JMNMFW0rHbS9LD2alLD1K3eisUBkVeCHEzUCel3DXs1GTgRNj72uCxSN+xQghRIYSoaGlpMZIchSIjsdvLVTibHL0LhBDrgAkRTn0b+BZwbaSPRTgWsShKKVcBqwAWL16siusw9Cqw8jLtx+5Cm8KRBEUS0BV5KeWySMeFEGcBM4BdwXjVFGCHEGIJAc99atjlU4B6w6lVnIbdAqOwn0zw5J0eNnUyCYdrpJR7pJTjpJTlUspyAsK+SErZCLwE3BGcZXMh0CmlbDAnyQqFs0jl6XUKc3Byb0rXk0+QV4AbgUqgD7jTIjtpj558KIFRZEIJULPcEsc0kQ9686HXErjHrO9WKJxMJois3ahwTeKoFa8KhcLx2N1hTeV+hBJ5hUIHvXir3QKkUGihRN7hqJi7Qo9MKCJ2D2ym8k+sRD7FSeXClyroiaiKF1uP0xsyJw8MK5FXKBQKHfQk3MkNvRJ5hUIHFZNXpDJK5A1idcxc6YdCDyd7kQr7USKvMISSF0UycHo5c3JvTom8QZx8cxWZgSqDCi2UyKc6apGI7SiRtR5VzhJHibzDUQKiUNjuy+ji5PQpkTeIk2+uIjlYPfCpFsRhe0WzezGWEZTIKwyh5EehcDZK5BUKgyhHOwnYva2B3qpnBxcCJfJ62Hxz9UIBao609Th5yXrGoIp5wiiRVxjCwQ6Maeg3tIp0J5VXPSuR10Pv5iYnFdHt25wA1ZOwHyeHChT2o0Te6aj663iUyCqcjCGRF0L8txCiTgixM/h3Y9i5lUKISiHEISHEdcaTahMOr792T+3KBH2zOyafCb+xHqrHmDhmPOP1YSnlT8MPCCEWALcBC4FJwDohxFwppc8Ee47C7gpot32F4/0ARYZjVbhmOfBHKeWAlLIaqASWWGRLYSOqkVEonI0ZIn+vEGK3EOIpIcTI4LHJwImwa2qDx9IOy1c7WvrtCjOwu6Gz234ysDtkpoeT74GuyAsh1gkh9kb4Ww48DswCzgUagJ+FPhbhqyL+DEKIFUKICiFERUtLS4LZsA4H3zuFImOwOyZv99iXEXRj8lLKZbF8kRDi18DLwbe1wNSw01OA+ijfvwpYBbB48WKlqSlGJsws0RUYuz15e80nhQwoZpZhdHbNxLC3HwP2Bl+/BNwmhHAJIWYAc4BtRmw5FbsLnyr7CoX92N3T0MLo7JoHhRDnEtCaGuBLAFLKfUKI54H9gBe4J1Vn1tjtqdrdiChi2LdEjcsoHIwhkZdSfkbj3P3A/Ua+X+F8lABZj92OhiK1USteUxxV/+3H6nug9/WqDNiPk++BEnkdRCoPqycBJxdus8iALCrSGCXyOuh1la334rQN2N0GOXnAKVlY/QtkQkOqsA4l8jqoCqaN0xeppAOqIQW/wyuik1OnRD7FURKbBHR7c1Y/41XnvKMlxhz8dk9VTuGfWIm8QeyuYHaHaxSKTEA9/i+NsfvWObjsZAxOvwWqjCi0UCKf4qiYuP2ogdf0J5VvgRJ5g9heAW3WeBUusp8sdRMsR2/g124Z0EKJvMNxcuHJFPTjsRbbd/g02kzAdmfOAErkDZLC916RIugJjNL4ZJC6NV2JvA6p3IIng0wQGN0FcTZvUJaVlQl3wV7s7s0ZQYm8Dk5fhKGqt0JhPc5WAW2UyOugv62B1Qth9OKxSuZtx+qYvG4ZsNa+IrVRIq+D3SvtnE4mNDJ2FwE9+3ZPo82AIhBDOMbuUhIdJfI6OH3qVAbUL8dj9zx5u0PymVAGnR621UKJvA5Ov7WZ4EUptLF7nrzqzTkbJfI6+G2O16SwA5E22D6zQm8KpfLkbcfJ9VSJvA4+HZF38s1VmIPdt9juTfD0sLsnkQycvAGZHoZFXgjxFSHEISHEPiHEg2HHVwohKoPnrjNqxy70RN5q9GKBmVDB7EavN2e3CNseLsmAIpjCGm/sQd5CiCuB5cDZUsoBIcS44PEFwG3AQmASsE4IMVdK6TOa4GTjs/nu6rUxdi+EybZ71C8JeG1v6LXP230L7LafDPQacie3AUY9+buBB6SUAwBSyubg8eXAH6WUA1LKaqASWGLQli3oevI2z5G2u4LZbT8Z6M6wsrgM6Nm3+xZkQm8ylT15oyI/F7hMCLFVCPGWEOJDweOTgRNh19UGj52GEGKFEKJCCFHR0tJiMDnmY/fUKb02xu450plQwb0+u8uAsxfEpX8JcMDguwF0wzVCiHXAhAinvh38/EjgQuBDwPNCiJlEvu8RfwYp5SpgFcDixYsd91PpVXDL9y1xcunB/pkdycDn92uet3uevN33wO5GJhk4uxZqoyvyUspl0c4JIe4G/iIDSrRNCOEHxhDw3KeGXToFqDeYVluwe+DV6YVLxeStx+nhGrvtJwOnO1taGA3X/A24CkAIMRfIA04CLwG3CSFcQogZwBxgm0FbtmD/wKuzC1dudvrPwtWfRmvtPdJrY7Kz7L0Hdg/+JwO9O2z3DCstDM2uAZ4CnhJC7AUGgc8Gvfp9Qojngf2AF7gnFWfWgP3z5HUiBbYXrkwQeds9eR37OTaLrCsn/cuAgzVcF0MiL6UcBP45yrn7gfuNfL8T8Ph0VNZi7A4X6ZGXARXcq1MG7J5dY3dDnwllwO4evRHS/+4YpNhltLNjDK+eK28zudnp31UftHl2jd09CT0ywZMf9Nrb0Bsh/e+OQRZOGgHAW99YGvG81fc25Ml/7uJyzpxcerr9JBWunCzBI7efd9rxZDSCCyYG8v3SvZdYbisS5aMLAfjWjfNtsR8qA9E8ZrsFpigJZeD86SMBmDmmyHJbkSjMywbg0tljbLFvBCXyOvj8kuwswfTRRbz/X9ck3X7Ii7t24Xj+9uVLOG9aWdLTkJMl+NIVM7n5nEn8501nfOBcaX6u5fbHlLg4Z2oZZ08pY8O/L7Xc3nBCMwQ/dcF0qn50Y9Lth0KGj9x2HjUP3MQti6YkPQ2jivL4zIXTqXngJp79wgfXNY4szLPc/qSyAmaMKWL9vy9l13eutdzecEIhsx/fejY1D9zEuBJX0tOQKErkdfAGRR5gZFEeNQ/cxPeXLxw6b/Xsl5AXl5udRU52Fn/98iWs+7fLk2ZfShn8DQJF5V8um8mB718/FKZJxhRpn98/NLg4Y0wRNQ/cxE1nTQSSEy4KNbQ5WYKsLEHNAzd9oFdj9bjJqTIQyOvPPnHOB3qWyYgXe33+oXpw2Zyx7P3edUM9nCJXtuX2w8vAiMJcjtx/A3deUm653RCeYMguN5iGbd9exm/v/NDQebsXzGmhRF6H8MIV4o6Lyvn5becCcLJnwFL7IYEJn48+e1wJ2/8zsHyhpdta+74wgQtRkJfNkftvZP6EEtp6By21D4EKNPwePPbpRXzh0hlkZwnLpzD6fKf/BjefM4mng5W8tTf5ZWD66CIO/fB6AE52J+Ee+OUHGtRiVw5vfuNKlswYlRT7Hp/8QP5zs7P47kcW8rVlcwH9mLlRhupB2GyyK+eN48W7LwKsLwNGUCKvg9cvyY7grk4dFfBi6tr7LbUfWm05XORGFeVRmJdNfYfbUvuRBCbEhBH5NHZZax8CFSwngsc+cUQ+bo+frn6vpfaj/QYTRxQA0NBp8T0YamQ+WF1dOdmMKXbR2GVtGQQ+0JsLZ9KIfBqSYd/njzhdd1xpIGxitbMVCpkNLwPjSvIBaLbY2TKCEnkdBr1+ciMMeI0KxiE7+z2W24fTK7gQgpGFeZbbDxXuSGGRCaX5NHVZX7g9Pn9EgRlfGqhgTd1WN3SBUMXw5fuhuKzVv4E3isCE0pCMexAQ2Qj2S/NpToZ9v4yY/7HFgXtgtch6h4XMhuwHy4DVPWojKJHXwe3xkx9B5EMzCnoHrPUi3Z5ABc/PjZSGbHoGrBX5kP2C3NPjrsWuHMvzH0pDQYT8h2b29CThHmiVgT6L7Q94o5eB4nzr74HH58cvI0+VLMrLYcDrt3xcYsDrj2h/TFBk2ywOl4ScreG9ifzcbIrysmntsT5klShK5HVwe33kRxE4gJ4Baxfyuj2B74+UhiJXDr1Jsu+KYL8wL5t+j8/ymHi0exCa1tY/aO1v0O+JbD8vJ4vcbEGfx74yELoHVtKvYx+gb9Dihs7joyBPqwxYG5N3e3zkZouIIaOCvBzcXucu6Fcir8OAxxdR4PJzs8gSJMGTjl7Bil05lnuxA97o9gvycpDylLdvFW6Pj/ycSBU86ElbLPLuKCIPgR5OMhqZkK3hFOZlJyX/EEXkXSGRT0JDG6EMhI65k9DQRS0DeVm4Lc6/EZTI6+D2+CN2k4UQyfGkNbrqRXnWd9WHwkURusrJ8uKi3YOCpHmRke1DqDeVrJBd5IbO6kbGPahlP3kiH8mTz8/LGjpvJXoNvdX5N4ISeR0GvJE9CAjMbhiweOpWqAJHSoMrN4tBi/fW0fLiCpJUwd0eH/maXXU7vTjrK/ipcEnkhrbX4kZOuyeRvLGpiGUwN0me/KAvYv5DabC6kTGCEnkdonmREBiIsnp+rtvrIy87K+J2rnnZSbAfGniNILKhQm9lAff7JQNef8RGLmQ/GY2MVgW3WmA0G9qk2j+9HuQnSWTdUUQ2ZN/y3ozHH7UM5CuRT2263R6Koyzdz8ux3pPucXujrijMy8myvCcRGnMojNRVDhb6AQtj8iEvNdJvMGTf4t+gdzByqCCUBqsH3foGvWSJyLNbXLnZuD1+Swe/Q/cgYrgkmCYr74GUkj6Pj4K80/Ofmx0Y/LZaZHsHvVHLQEGe9Q2tEZTI69DR76GsIIrIZ2cxaHEF7+j3RN0bxJWTbbkn39EXEPlIaQh5dlaKXMh+WQT7riGBsfYedPYNatyDLMsHnjv6PJQV5kV8zJ4rCSLbqVkGrPfkewa8+PySsoLo9cDqhj5wDyLrgCsny1JHxyhK5DXw+yWd/Ro3N9d6T7qzz8OIKPbzkhAuateo4K4c6z15rUYmK0uQl50EkdUoA/m52ZY3Mnr2wVqR7+gPNbSnp8EVauiTUAai1YP83CzL70G7ZkNvfRkwghJ5DbrdXqSM7EVCcmLiWoUrFC6ysqve0T9IXk5WlHhsqIJb6Mn3BxaZaDW0Vtr3+yVdGr25/FzrG5nOPm37EJjqaxVDvakInnRorMRKkTtlP5onnW25J63lyecnwdkzghJ5Ddr7ggITLVyTBE+6Q6OCJ6urXlaQGyVUYL0XeaonoVHBLbTf5fbglzBCw4uzOh7b3jcY1dFwDc0Tt9KT1mrorbc/VA+1QmYWloFBr5+eAa+OJ5+mIi+E+JMQYmfwr0YIsTPs3EohRKUQ4pAQ4jrjSU0+LcFNj0YXa3vSVuH3S1p6BqLaD4m8lWlo7h5gTHHkvbOT4cmH9gQZVRQ9DVZ6saHNxyaOyI9u3+IK3tjpZoKGfbDWk27odDNxRL7mmICVZaBR5x7k5VhbBkI7rY4q0hqXcW64xugzXj8Zei2E+BnQGXy9ALgNWAhMAtYJIeam2sO86zsCu+tNGVkQ8bzVUyhP9g4w6PUzZWRhxPOhJwVZmYa69n6mjY5sPxmefH1HPwW52RqevLUiGyoDk8qilQFrPXm3x0dr7yCTo9jPT4InX9/Rr9HIWV8G6jr6EeLUhnTDceVa60nXdfQBGjqQCeEaEWjiPwE8Fzy0HPijlHJASlkNVAJLon3eqdS2a1fwPItnt4S2MY5WwfOyrRV5KSW17X3RBSYJnnxdez+TRxZE9CIDabBWZIdEPorIuXKtnVkRsh9NZF1JmOFU39Gv0chZXwbqO/oZV+KK+vjDQENvnf2QDkRztlw52fj8UveB73ZhVkz+MqBJSnkk+H4ycCLsfG3wWEpR297PyMLcoVV9w8nLtrYFrwtW8MlRPIg8i2Pynf0eegd9Gj2ZoBdpZQXr6IsqMIE0WHsPajv6yc0W0UNWOdkM+vz4LdqFsU6nJ2H1WgWvz09T9wCTRkS2n5UlyMvJsrQM1HX0D+3dH4l8iz35Wh1nKxljY0bQFXkhxDohxN4If8vDLrudU148QCS3K2ItEEKsEEJUCCEqWlpa4ku9xVQ2dzN7XHHU81YvRqpq6UUImDZKO1xilSd/tKUXCDyFKLL90MwOa+z7/ZKjzb3MGhv94c1We/KVTT3MGFMUccVxyD5YV8GPNPUAMDPKb2D1Bl3VJ3vx+SWzxkW/B1bPE69s7mHW2Oj10Oq1ClUtvYwrcWkuiAPnirxuTF5KuUzrvBAiB/gn4Pyww7XA1LD3U4D6KN+/ClgFsHjxYsc8KFFKycHGbm4+Z1LUawIxeesE5kBDF9NHFQ7tWz4cq2PyBxu7ADhjYknE81Z7ccfb+uj3+Jg/IbJ9CFSwbrd1+6YcbOxm0fSRGvZPhSuiiYARDjV2M6oob+jhGFHtW3QPDjZ2AzBvfGnUa6xcK9DWO0hz94BmGbA6XHOgoYv5E6PnP1mL8hLFjHDNMuCglLI27NhLwG1CCJcQYgYwB9hmgq2k0dDpptvt1S1cVs5sOdjYzfwJ0QvXkMj7rClcBxq6KMnPidpNBWu9uCGB0fgNrJzZ0OX2UNfRr1MGrA1ZHWzqZt74kqhjElZPoTzY2EV2ltD05K1cKxByNObp3AOryqDH56eyuSeqowPJWRBmBDNE/jY+GKpBSrkPeB7YD7wG3JNqM2t213YAsGCStsha5UV39nuoae3Vtp9tbSxwb10XZ0wsjSowYK0Xt7u2g5wswbzx2p68Zfmv7QS0y8CpxUjmp2HA6+NgQ1ds9i27B53MGVc81JhEwsoVn3tiuAdWzm453NTNoM/PAk1P3voFYUYwLPJSys9JKZ+IcPx+KeUsKeU8KeWrRu0km63VbbhysjhrclnUa/Kys/BLLBlVr6hpQ0pYXB49VODKtU7kewe87KnrZEn5KM3rrPTkt1W3cebkEZphECs9+a3VbQgB52uGa6zz5PfUdjLg9bNkRvR74LJwMZLX52f7sXZN+2CtJ7+tuo2ZY4uiDnxDqAxaVAaq2gBYrFEPrB6bMopa8RqF92raWDRtZNRpW2CtyG6tbiMvO4tF06ILjJVTKCuOtePzSy6YqVfBrdmF0e3xsau2gwt0BcY6T/69mjYWTCylNMoupBA+hdCaMgDwoRgExoqGbl99F32DPn2Rt2itgM8v2VbTZmsZ2FrdytRRBTohS2cPvCqRj0B77yD767t0C7eVN/fdo62cO60s6gNnpyMAACAASURBVMMqILyrbo39nCyh6cWCdZ78ezVteHwyhntgjSffP+hjx/FYvNjQFEbz07D56Enmji+OutISAvkXwpoy8M7RkwD6vTmLwiX76jvpdntjKgNWTGP1+SVbq9u4YMZozeuSserYCErkI/DGwWb8Eq4+Y5zmdVaNqjd2utlT18nlc8bo2Ldu+ty6A018qHxU1DUCIazy5NfubyI/N4uLZ+n8BkEvzuxN2jZVnsTt8XP1/PGa152a3WKuyHX2e9ha1cZVOvaFEJaFK9bub+LsKSMYF2WlaQirPPm1+5vIEnDFXL16GJxKbHLYtKKmjY4+D1fOi82+CtekEGv3NzKhNJ+zJo/QvM6qzZnW7G8E4PozJ2heZ1W4qLK5h8rmHl37YM0cZSkla/c3cdmcsbrTEq1aiLJmXyMl+Tm64SqrGto3DzXj9UuuXagt8qE0mG2/ucvN+8c7uOYMfftWrVVYsy/gaGj1ZMC6mPjr+5rIy8niinljte0nYdWxEZTID6Nv0MvGwydZtmCc5qwSsG5Z/2t7G5k1tojZ46LPKgnYtyZU8Pq+QCMTi8BYMbtmT10nDZ1url0Qm30wt4J7fH7eONjMlfPGkZutXUWsCpm9vq+RsSUuzp0SfeA/PA1m2197oAmAaxfG1tCbbb/mZC+Hmrpjs2+ByEopeX1fI5fOHkNxlHUqQ/bVwGtq8fq+Rvo9Pj58dvRFUCFcFjwVp7nLzdbqtpi8aCtWO0op+ceues6dWqa5lDyEFZ78X9+vIy87i2UxeZHmh8w2Hm6hrXeQm86eqHutFZ58Z7+HNw40c8OZE6KutA3HCk/6b+/XMXNsEXPHR19pGiL0CEJT7e+sQwj93iyE7WlvYhp2HG+nrqM/xt6sGnhNKf68vZZpowp1B5vAmlDBn3fU4vNLbj1/qu61udmCLJMH3fbWdXGwsZtbz58S0/Vme/KDXj9/31nPsgXjGKnTTQdrFgP9eXsto4vyuGq+diwWrOlNvby7ngGvn4/HUAYgFBM3L//VJ3t5r6adj58/Vbc3C+Zv9+z3S/68vZZLZo3RnNUSwmVBQ//8e7UU5mVz41mxNPRq4DVlqG3vY/PRVm5ZNCVmDwrM8+KklLxQUcuS8lHMGBN9hWGIwKCbuV7c8xUncOVk8RGN7RzCMXuO9PqDzbT1DsbRyJhbwdp7B1l3oInl507WDdWANasdX6ioZd74Es6cHH0BzvA0mBmqeHF7LVkC/mlRbHsKmj2FcUt1K7Xt/TGXAbM96d4BLy/vruemsybqhmrA+XvXKJEP48/ba5EyjsJtshf5Xk071Sd7+cSHYvPgwFyRdXt8/H1nHTecOYERUZ5GNZzAknLzBOaFihOMK3Fx+Rztwa5w+2DePfjr+3V4fDL2Rsbk1Y5HmrrZeaKDjy+eEpMXHUqDWaEKr8/PiztquXzu2Kj7tw8nNIXRZ9IUxhcqailx5XBdDPH4kH0w7x68sqeB3kFfzPUwT8XkU4NBr5/fbz3OFXPHMjXKro/DMbub+Putxyh25XDjWbEVbjA3HvvSrnq63N64GpmAF2lO4T7W2sv6Q818YvFUcmLwosFcT97vlzzzbg3nTSvTXEYfTihkZlYj8/TmGlw5WXzsvNh35jbTk1+7v4mGTje3L5kW82dOebLG09DSPcDq3Q18bNHkmDd8M3NBmpSSZ7ccY+bYIhbrrBEJkZ0lyM0WanaN01m9p56W7gHuvKQ85s+YObOjobOf1bsb+OSHpurOTR+eBjO6iVJKntpUzfwJJVw0U3vxRziu4INTzJin/vTmGrKF4DMXTY/LPphTwTccaqamtY/PXzIj5s+YGTJr7x3kxR21fOy8yYzWWMY/HJeJMfmn3qlm6qiCmAa9Q+Sb6Mn+fusxBn1+PndxecyfcZnYyFQca2d3bSd3Xlwec08KkvMw8URRIk9A4H77Tg2zxhbFHCaAU4XbjBb8/zYfwy9lXIUbzFvx+U5lKwcbu/n8pTPiKtxmTSHsdnt4oaKWD589MeYwwQftG/8Nnnqnmgml+THNqBieBjMa2ufeO47b4+fOOBqZU/aN5393bQfv1bTzuYtnkB3DmFQIl0n79wx4ffxuyzGumj+OmRr7xw/HzE3inny7mhEFudwSY7guhNXbHRtBiTyB6VK7azv53CUzYhpwDWHWwGvfoJfnth3n+jMnxBwqCuHKzTYlXPLkpirGFOdp7p8f0b5J09deqKilZ8CbgMCZ48kfauzmncpW7rh4ekwDrsPTYLQMeHx+ntl8jEtmj9bcVjeafTME7rfv1FDsyuETi+MTuHyTBp//sauBkz2DcfWkwLyB1xNtfazZ38jtS6bF1ZsOpMG5z3lVIg/85u1qSvNzuCXGAdcQZsUCX9xeS2e/hy9cGl/hhkBvwujAZ2VzDxsOtfCZC8s198qJaN+EhShen5+nN9dw/vSRnDNVf/FPOGZt0PXUpmryc7O4/UOxx6LD02C0oX1lTwONXe64BW7IvsH8N3W5eXl3PR9fPIUSjQ3ZImHG4LOUkic3VTN3fDGXzI49XAjmDbw+vbmGLCH47MWxhwtDWP0IQiNkvMhXNnfz2r5G7rioPO7WOyc7i5wsYahweX1+ntxUzTlTyzR3nIyGGZ78k5uqyMvJ4tMXJiJwxnszq/c0cLytjy9elkAjZ8L0tfqOfv7yfi0fP39qTHPzI6XBSEMrpeTxN48yc2yR7j4pUe0bLAOrNlbhl3DnxQk0MiZ48m8eauFAQxf/ctnMuMKFYI6z1dE3yJ/eO8GNZ02MaRHgcPIs3O7YKBkv8r/ccJT8nGw+n4AXDaGueuKFa/WeBmpa+7j7illxF24w7snXd/Tz5+21fHLxVM09u6PaNxiT9/slj22oZM64Yq5dEF8sHMzx5FdtrEJK+NIVMxP6vNGG9o0DzRxs7OaepbPjCheGCEyjTTz/bb2D/GHrcZafM4lpo+MLF4LxlddSSh7dUMnksoK4ZhWFMGPg9enNNfQMeLl76ayEPm9W2NQKMlrkT7T18fdd9Xzqgmm6myBFw0gF8/slj66vZN74kpj2aYls35gX96u3jiIl3JVo4TZYwdceaOJwUw/3XJmowBnz5E/2DPDH947z0fMmM2Vk/AIHgYbWqMBNGVnAzefGNx4SwpWTjdcvE354zVObqnF7fXz5ysQFDhK/B1uq2th+rJ0vXTEz7vEQMD67p9vt4bfv1HDNgvGcofEEKC2M7gT6yzcr2XmiI+HPa5HRIv/EW0fJFoIVlyfmwYGx6Wuv72vkSHMP91yVmMCBsUamudvNc++d4JZFU2JaPh7NPiRWwaUMePHTRhXy4Rj2iYlE6MEpif4GT26qZsDrT9iDg+B2xwna33y0lZ0nOrh76ayEBA6MbXfc2e/h/zbXcP3CCbob4kXDaG/qsQ2VjCl28YnFsa/PCCcnO4vsLJFwI/O7Lcfp7PfwlatmJ/R5MDbwerCxiwdfO8SGg80J29ciY0W+qcvNCxW13Lp4SlxT9oaT6EIUKSW/WF/JjDFF3BTD/hjRMDKz4zdvV+P1GRM4IzOM3j5ykt21ndy9dFbMi5+Gk5UlAvHQRASuz8Oz7x7jprMmMiuOKXvDyTdQwX+x/gjjS10xr7CNaN/A/jnPvltD94CXe65MXOCMlIH3j7ezqfIkKy6fEfegfziJDj73D/r4zdtVXDF3LGfHsONndPuJ96h/sb6SYldOQoPusWBI5IUQ5wohtgghdgohKoQQS4LHhRDiESFEpRBitxBikTnJNY9fb6zCJyV3XZ64wAEUu3LoHfDG/bn1B5vZ39DFl5fOimtO8nCKXDn0DHjjXozU1jvI77YcY/m5kymPYZ+caBjZZvXR9ZVMKM2PeRsJrTQkUsFDcVgjAgeJN7QVNW1sqWpjxeWzNB+UrWs/FDKLU2T6Br08uamaK+eN5UydZydo2jfQm3tsQyVlhbl8+oL4Z7SEk6gn/dy247T2Dhry4iHxtQpHmrp5ZU8Dn7u4nBGF8c1qihWjnvyDwPeklOcC3wm+B7gBmBP8WwE8btCOqbT2DPD7rce5OcGBpnBK8nPodscn8iEvfsrIAj6awEDTcPsen4y7gD+1qZp+j48vG/DiIfEHWW+rbmNbTRsrLp9pSOBCaYg3/70DXn67uZplZ4xLOA4bItHtlh/dUMnoojxuX5JYmGLIfoLPNfjD1uO093m417DAJdaTONDQxboDzdx58QyKYtgITItEVh0PeH38auNRLpw5SvNB3bHaT8TR+cX6SgpysxOaPh0rRkVeAqEaMgKoD75eDjwjA2wByoQQicckTOZXG6sY8PoMe3AAJa5cut2euD6zqfKk4TjsKfuByhFPQxOKw95w5gTmjE8sDhsiUU/+4bWHGVOcF9ceKdFIZKvbpzfX0NHnMaUMJLLd8q4THbx5qIXPXzoj7qm7w0lkQVr/oI8n3qriopmjOX+6UYFLbArjz9cdodiVE/cq70gUubLpG4zvHrxQUUtT1wBfuWqOYfvFrmx64uzRH23p4R+767njovKEpu7GilGR/1fgJ0KIE8BPgZXB45OBE2HX1QaPnYYQYkUw1FPR0tJiMDn6NHe7eebdGj567mRmj0s8DhsiXk9eSsnDaw8zcUS+oTjsKfuBLl48Dc2Tb1fRPeDl3iuNF+5EPPnNR0/yblUrX146O+ZNqLSINx7a5fawamMVV88fx3kJrE0YTiI7gf5s7WFGFubyWRMELpEFac+8W8PJngG+fu1cE+zHP4Vxb10nr+1r5AuXzjAlTFFakEtXHHXA7fHx2IZKzp8+kotnxbf4Kpr9brcnroeJ/++6IxTkZvMvCawPiQddkRdCrBNC7I3wtxy4G/ialHIq8DXgydDHInxVxNxLKVdJKRdLKRePHRv7vjGJ8ssNR/H4JP/vauMCB1Acp8hvONTMjuMdfOWqOYbDFBBoZCB2T76td5AnN1Vz01kTY95pUYtQBe+P0YuSUvLQmsNMKM3nUxcY9+IBCnKz6RuM/R48+XY1nf0evnaNcYGDYKjA64t5XOS9mjY2Hm7h7qWzYtqvPBb7EHu4pmfAyxNvHeXyuWMNhykAcoOzW+Jp6B5ae5gRBbl8wSSBK83Ppas/dpF/bttxGjrdfP3auQmtT4lk3y+hN8ZyeKChi3/squfOS8oTWp8SD7olTEq5LNo5IcQzwFeDb18AfhN8XQuEBxqncCqUYxsNnf38Yetxbl00xdBgYzgl+bn0DHjx+aXuAKrfL/nZmsNMG1XIx+PcH0TLPsQu8r966yj9Hh9fu8acRq7ElYMQxFzBNh45ScWxdn7w0TMNzaYIp6wwl84Y7bf3DvLUpmquXzjB0GBjOAV52UgZCFfE0jP52ZpDjC1x8ZkLy02xX5gXX0P7203VtPd5+DeTGjmAwtzsmAVux/F21h9s5hvXzaM0zi0UolFakMvxtr6Yru0b9PLYhqNcPGs0F88aY5L9gJR2ub0xbQvx8NrDlOTnsOIyY2NisWA0XFMPXBF8fRVwJPj6JeCO4CybC4FOKWWDQVuGeXR9JRJpeKApnNKgJx1LAX99XyP76rv46tVzDMfiQxQPxeT1Ra65y83/vRsKVRmLxYfIyhKMKMilvU/fvpSSn605xOSyAj6Z4JzoSIwoyKUjRpFf9XYVPYNe07x4gJGFgXhqR/+g7rWbK0+ypaqNe5bOMiVUBYFGDojpHnT2e/j121UsO2Mc58a5T5AWIwpz6YzBPsBDaw4zuijPlFh8iNL8nJgdjWfePWZaqOqU/cA9iCUNu2s7WLO/iS9eNtOyGTXhGO0rfhH4uRAiB3ATmEkD8ApwI1AJ9AF3GrRjmBNtfTxfcYJPLJ4a906PWoSHS7S8Ep9f8tDaw8waW2R4Rk00+3r88s2jeH2Sry4zx4sPMbIwLyaRXXegmd21nTx4y9lDT9Mxg7LCXDpiEJiTPQM8/U4NHzl7Utw7PWoxMiSyvR7NfU+klPx0zSEmjsjnNhMGnIfsBwftOvr0G5kn366iy21uIweBMtAeg/2tVa1sqjzJf950huEZNeGEYvJSSs3wS7fbwxNvHWXpvLGGB5yH24fYRP5nawLjMfE8u8IIhn5lKeUm4PwIxyVwj5HvNptfrD+CEMJULx6GD3xGr+D/2FXPkeYeHvvUIkPz4ocTali6dUb26zoCoaqPL57K9NHmhKpCjCjI1RUYf7CRKx9daHhe/HDKCvLo6BvE75eaK4efePMoA16f6Y1cWWFsIvvmoRZ2HO/gRx87y7RQFQRCZjlZQldk23oHeeqdwKyqhZPMCVWFKCvMpU2noQ305A4zrsTFP19obF78cErzc/H4pG7I7KlNgVlVX79mnun2IRCu0aKipo23Drew8ob5ce/2mSgZseK1qqWHF3fU8akl0xLaYU6LsoJTXlw0PD4/D687zBkTS7khzgdS6FGcn0OWCMSatfjFG4FImtFFH5GIJSb+6t5GDjR08dVlcxJe3apl3y+hRyNk1tTl5tktx/jYeVMMrW6NRGjfozYNkZUy0MhNHVVg2nhMCCEEZYX6IbNfbTxKr8mhqhAjC/N0G7lNlSfZVtPGvVfNNrWRg1Mxca1y2NE3yG/eruK6heM5a4q5jVws9iHgxY8pdnHHReWm2tciI0T+J68fwpWTZcqc6OGMLQmMjDd3u6Ne80JFLcda+/j6NXMT3qMmGtlZgjHFLk37VS09vLC9lk9dMI1JCe5Ro4VeV93j8/OT1w8yd3wxN59jrhcPpzxprZjww2sP45eSfzXZi4ewcI2G/Zd3N7CnrpOvXj3XtPGYcMoK8zQb+vqOfp5+JzAeM9fg2ohIjCzM1bTv90seePVgYDwmjmcIx8rY4AyVlu6BqNc8/uZR08djQoyJwf7Gwy28W9XKPVeaNx4TC2kv8tuPtfPq3ka+dPmsIUE2k3ElgX1vot3c3gEvD609zPnTR3L1GfHvFR5TGkpdNGsUrh+/dpB8ixo5CIh8a89g1CmEz207Tk1rHytvOMPUUFWIUUUBkW3pifwbHG7q5vmKE9xxUbmp4zEhQo1MNJEb8Pp48PWDzJ9QktBWurEwqjCPNg2RfWjtYaTE1MHGcEYW5dHl9uKJshPmS7vq2VffxTeum2fK1OHhhPafauqK7OzUtvfx2801fOy8ycyfYHzq8HCKXDmUuHKi2vf5JT965QBTRxWYNnU4VtJa5KWUPPDqAcaWuCxbcFBakENeTlZUkV21sYqTPQN868YzTJmPG4lxJfk0d0W2/15NG6/va+KuK6xp5AAmleXTN+ijq//0cEm328PP1x3hwpmjWDrPmnUQod5JQ0fkCvbjVw9S5MrhXosaubycLEYX5dHQGdn+77cc50RbPytvtKaRg0BDH01gDjR08eKOWj53SXnC2ynrMUFDZN0eHz95/RBnTi6N+/GSsTIk8lF6tD99/RAC+PdrzY3FhxNwtiLb/8uOWg42dvPN6+Zb0shpkdYiv3Z/E+/VtPOvy+aYOpIfjhCC8aUumiMU7uYuN6s2VnHjWRM4f7rxlZXRGFcS2ZOXUnL/6gOML3XxL5clvp2yHqFtius6+k87t2pjFa29g6y8wbpGLiTy9RHsv3u0lTcONnPPlbMtXTo+eWRBxPx39nt4ZP0RLp09hsvnmDMnOxJTRhZS19EfccXlA68epDQ/l3uWWtPIQSD/AHXtp/8Gz757jLqOfr51wxmmhytDjCnOQwhoiuDs7Knt5G876/n8pTMsCVeGGF+aH9F+/6CPn605zDlTRiS8pbYR0lbkB71+HnjtIDPHFpk6JzsSk0YURFyI8fC6w3j9fr553Xxr7ZcVcLJn4LRVn6v3NLDzRAdfv3aepTHAaCLb0NnPr9+u4iPnTIr72a3xUJqfS4kr5zSRDXWRJ43IN3VOdiSmjCygtv30MvDLDZV09nu474b5ljVyIfsenzytsX/7SAtvHW7h3itnWzonO9RDqB0m8u29gzy6oZKl88Zy8WzrGrmc7CzGlbhOa2SkDJSBUUV5hrbUjoUJpfkRG7mn3qmmscttaW9ei7QV+Sc3VVPV0st/3bTA9Nkcw5k9rpjK5p4PxKR3nejgj++d4J8vnG7a6tpozAnuwVPV0jt0rHfAy/2rD3DGxFJuWWTubI7hTAl6cTWtvR84/oOX9yMlfPM667rIISaPLODYMPt/fO84e+o6+eb1802fzTGcKSMLqW3vxxfmSR9p6ubJTdXcsmiKaatro9sP3INwZ2PA6+O7f9/H9NGFfOYic6csDmdSWT5CcJqz8+DrB+kZ8LLyhjMstQ+hetj9gWOr9zTwblUrX716jmmra6PaH19MY5f7A3vo1HX08+j6Sq5dMJ4LZhrfIycR0lLk6zr6eeSNI1y7YDxXzrdmsDOcOeOK6XJ7h7wor8/Pt/66h7HFLktG8k+zPz4g8oebThXwn79xhIZONz/86ELL4sAhRhe7mFCaz776rqFjGw+38MqeRu69crYlg53DWTCplL31XUMNbWvPAA++dogLZ45ieYKP1YuH+RNKGPT6qWzuAQIe5H/9fS9FrhxW3mBtTy5gPzCYuLeuc+jYb96upupkL9+7eaHljZwrJ5uZY4rYV3/K/o7j7Ty37QSfv6Tc1MVn0Zg7voTDTT1DIatut4fv/2M/Z04u5dNJGOycF5y1dKSpZ+jY917aB8B3PrLAcvvRSAuR336sjc8+tY3eAS9+v+S+F3cjkUn7YUNzbitq2gH49dvV7Kvv4jsfWWC59wAwfXQRRXnZVBwL2H//eDtPbqrm9iVTTV3Vp8WZk0cMPaOys9/Dyr/sYcaYIlYk+HDseDlr8ghaugeo73QPCWzvgJfvLz8zKV3kUDjq/eOBe/CHbcfZUtXGN66bx2iLN6ACmDAinwml+bwfvAeHGrt55I0jXL9wAkvnWe/oAJw7dSQ7T3QgpaR/0Md//Hk3E0rz+eoy6x0dgAUTS+n3+DjQGHA2fvTKAVp6Brj/o2dZ3psHhhaYbT/WBgQWQK7Z38T/u3qOZQPesWDNaGSSkRLeOtzCS7vqaejo5+0jJ/nRx85K2g97zpQyRhTk8sreBkYU5PLTNYe48awJhh7rFw+52VlcOmcMbxxo4sQVs7j3D+8zcUQ+911vfRc5xBVzx7DuQBM7T3Tw6PrKwOMV77ooaTMJLp8bmLnzj1315GQJXtnTyH03zLdkTngkZowuYtKIfFbvaWDuhBK+/4/9XDZnDJ8ycfsCPS6ePZq1+5qo6+jnnj/soCQ/lx9+7Mzk2Z81mhd31PJuVSsvVNRypLmHZ7+wxJSdNmNh6bxxCAGv7mlkX10Xz207wV1XzLJ0PCicCSPyOWNiKav3NHLZnLGs/MsezptWZvlWwnqIeB8bZyWLFy+WFRUVcX9OSslHHt3E3rpAC37r+VP4ya1nJ3WQ439ePcCv3qoCAt22F+6+KClefIi3j7TwmSe3AYFNy/7wxQsMPbMyXjr6Brn8wQ1Dy7q/d/NCU/ZKj4dP/XoLm4+2AnDdwvE8/unzLZvNEYlfvlnJg68dAmD66EL+fNfFlk1bjcSO4+3c8vhmpAxM6/y/O5dwkQl7pcdK36CXyx98k5PB9Qpfv2YuXzFpS+9YuevZ7by2rxGAy+aM4anPfciSxWfReHbLMf7rb3uBwEDsi1++eGj2mZUIIbZLKRdHPJcOIg+BxQ5PvHWU2WOL+cxF5ZbHoYfj9vh4/M2jDPr8fOnymUMLZJLJi9tr2X68nTsumm7Jgg89dhxv50/bTnDl/HFcb/L2DbHQ3O3msfWVTCor4M5LZpi6CVoseH1+nnjrKB19HlZcMXNooVwyeXl3PZuOnOT2JdOS5sGGs7euk2fereHiWWNYfu6kpM8m6egb5NH1lYwoyOWLl8+0fCxiOH6/5Kl3qqlt72fF5TMtnbIZTkaIvEKhUGQqWiKfFgOvCoVCoYiMEnmFQqFIY5TIKxQKRRqjRF6hUCjSGCXyCoVCkcYokVcoFIo0Rom8QqFQpDFK5BUKhSKNcdRiKCFEC3AswY+PAU6amJxUQOU5M1B5zgyM5Hm6lDLio9ccJfJGEEJURFvxla6oPGcGKs+ZgVV5VuEahUKhSGOUyCsUCkUak04iv8ruBNiAynNmoPKcGViS57SJySsUCoXidNLJk1coFArFMJTIKxQKRRqTFiIvhLheCHFICFEphLjP7vQYQQjxlBCiWQixN+zYKCHEWiHEkeD/kcHjQgjxSDDfu4UQi8I+89ng9UeEEJ+1Iy+xIISYKoTYIIQ4IITYJ4T4avB4Ouc5XwixTQixK5jn7wWPzxBCbA2m/09CiLzgcVfwfWXwfHnYd60MHj8khLjOnhzFjhAiWwjxvhDi5eD7tM6zEKJGCLFHCLFTCFERPJbcsi2lTOk/IBs4CswE8oBdwAK702UgP5cDi4C9YcceBO4Lvr4P+HHw9Y3Aq4AALgS2Bo+PAqqC/0cGX4+0O29R8jsRWBR8XQIcBhakeZ4FUBx8nQtsDebleeC24PEngLuDr78MPBF8fRvwp+DrBcHy7gJmBOtBtt3508n7vwF/AF4Ovk/rPAM1wJhhx5Jatm3/EUz4ES8CXg97vxJYaXe6DOapfJjIHwImBl9PBA4FX/8KuH34dcDtwK/Cjn/gOif/AX8HrsmUPAOFwA7gAgKrHXOCx4fKNfA6cFHwdU7wOjG8rIdf58Q/YArwBnAV8HIwD+me50gin9SynQ7hmsnAibD3tcFj6cR4KWUDQPD/uODxaHlPyd8k2CU/j4Bnm9Z5DoYtdgLNwFoCHmmHlNIbvCQ8/UN5C57vBEaTYnkG/hf4JuAPvh9N+udZAmuEENuFECuCx5JatnMSTLiTiPQ4+EyZFxot7yn3mwghioEXgX+VUnYJESkLgUsjHEu5PEspfcC5Qogy4K/AGZEuC/5P+TwLIT4MNEsptwshloYOR7g0bfIc5BIpZb0QYhywVghxUONaS/KcDp58xgbFOQAAAcFJREFULTA17P0UoN6mtFhFkxBiIkDwf3PweLS8p9RvIoTIJSDwv5dS/iV4OK3zHEJK2QG8SSAGWyaECDle4ekfylvw/AigjdTK8yXAzUKIGuCPBEI2/0t65xkpZX3wfzOBxnwJSS7b6SDy7wFzgqP0eQQGaV6yOU1m8xIQGlH/LIG4dej4HcFR+QuBzmD373XgWiHEyODI/bXBY45DBFz2J4EDUsqHwk6lc57HBj14hBAFwDLgALABuDV42fA8h36LW4H1MhCcfQm4LTgTZQYwB9iWnFzEh5RypZRyipSynEAdXS+l/DRpnGchRJEQoiT0mkCZ3Euyy7bdAxMmDW7cSGBWxlHg23anx2BengMaAA+BFvwLBGKRbwBHgv9HBa8VwGPBfO8BFod9z+eByuDfnXbnSyO/lxLoeu4Gdgb/bkzzPJ8NvB/M817gO8HjMwkIViXwAuAKHs8Pvq8Mnp8Z9l3fDv4Wh4Ab7M5bjPlfyqnZNWmb52DedgX/9oW0KdllW21roFAoFGlMOoRrFAqFQhEFJfIKhUKRxiiRVygUijRGibxCoVCkMUrkFQqFIo1RIq9QKBRpjBJ5hUKhSGP+Pz6mYRkfzdJaAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(s[4999:, 0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [],
   "source": [
    "an_bool = [1] * 13\n",
    "an_x = anmodel.models.Xmodel(channel_bool=an_bool)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [],
   "source": [
    "an_x.set_params(anmodel.params.TypicalParam().an_sws)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [],
   "source": [
    "s, info = an_x.run_odeint()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x1a1d2365d0>]"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO2dd3gc1bn/P0e9S5Ys2ZabXHHFxjimhF6CMQQTSIF7by435RISUrgkIRBu+CWEJCSEngLcG0K51NBDd6E3Yxv3KnfZsoolq0ur3T2/P3ZWluRdbZmdmbOz5/M8fryandl5z8w53/Oe9zQhpUSj0Wg07iTNaQM0Go1GYx1a5DUajcbFaJHXaDQaF6NFXqPRaFyMFnmNRqNxMRlOG9Cf4cOHy6qqKqfN0Gg0mqRi1apVjVLK8lDfKSXyVVVVrFy50mkzNBqNJqkQQuwJ950O12g0Go2L0SKv0Wg0LkaLvEaj0bgYLfIajUbjYrTIazQajYvRIq/RaDQuRou8RqPRuBgt8hrX0tLZyz/XHnDaDI3GUVJO5D1ePzsb2p02Q2MD1zz1GT944jN2NXY4bYpG4xgpJ/I/f349Z93+Ds0dHqdN0VjMgcPdAPR4fQ5borEDn1/S5dHvejApJ/If7TgEQIfH67Al7uLOJduouv4VfH6905jGGX705GdMv+l1p81QjpQTeY01/PXtHQB4/f4Bx70+PxNveIUnV+x1wixNgqiub+Nwp9qt35fX1TptwpB09/p49OM9+G12hLTID0Fje4/TJiQ9HT0+/BJ+++pmp03RmOCcO97lwnvfd9qMpOaupdv5xQsbeHm9vZWRFvkwvLHxIPNvWcqH1Y1Om8If39jKjc+v560t9U6borGQt7fWK+1Y1DR3OW1CUhNsCXX02BsqVmqpYZVYtacZgPX7Wzh58nBHbfnTW9UAPPbJXnbfeoGjtmgSz+FODwvveo+Drd1MHVHAm/91utMmaVyEaU9eCDFWCPGWEGKzEGKjEOJHxvFSIcQSIcR24/9h5s3VaGJHKt4X/N72Rg62BkYC7WzQwz01iSUR4Rov8GMp5XTgROBqIcQM4HpgmZRyCrDM+FsZYin4da3dvLe9wTpjNJokYfmWupDHfX7J5tpWW23ZeKBF6fBWf5ZvqePJT/c5cm/TIi+lrJVSrjY+twGbgdHAYuBh47SHgYvN3ssJJPDFe9/n639b4bQpyuLx+vH4AqNqVPSahXDaAvfwzYdC79z2p+XVnH/3e2zY32KbLRfc8z7n3vGObfczw0/+sc6xeye041UIUQUcB3wCjJBS1kKgIgAqwlxzpRBipRBiZUODs97yss11/NmIf/fXhfq25PAW7GTGTa9z9WOrgUAlGAkFtV+TQNbVHAbgYEu3rfdt7uy19X6JwG5HKGEiL4QoAJ4FrpFSRt1uk1I+IKWcL6WcX14ech9aSwjl3X3r4ZXc9sZW22xIZjo9Pl4xhoJtrWsLf6L2ojUaR0mIyAshMgkI/GNSyueMw3VCiFHG96OApBz/p2L4QWWOqjz184sJ/bg0iSYRo2sE8Ddgs5Tyjn5fvQRcYXy+AnjR7L00oalv6+aqR1fRbvP421jQDn1q4FQldcebugUejkR48p8Hvg6cJYRYY/xbBNwKnCuE2A6ca/yddCRDp93dS7fz+saDPP/ZfqdNURLVW2PJkMci4XQa7lle7awBIfjda5u5+vHVRx23+1mZngwlpXyf8I7a2WZ/3yqGKvhV17/Cwpkj7TPGRagkqE4LT7T0f2ZJYrImCu5/ZycAf/4XZ+3QyxqEYenmwHhglURLo1EdqQtMSJysvLXIhyFZvMBkwYmir/VGowr9s2LSDqHUJJZ9TZ3uWJvdocqyrrW7b2inihV2a3cvL7iwD0Wo+LBTHC3yEZAO+KAHDndx6h/e0mP2TXDCb5c5bcKQXPePdVzz1JqjlgJwQbWuUYyUFflIDodwMIrWYMyw/XBHjMsc6/hE0lBrLEjW3euu7ep0TF49Ulbkg3nR55f84oUN7GvqHPi94VP9/YPdNluW3Kg4GUp13XFHhMMViXAlKb+e/Gd7m3n04z1sGtxsNoShIZnWrVFALcIJqvOWHU2Xx0eP10dJXpbTpmhcjh5doyAK6KXGYi689z3m3rwEv1/avu9mOFTJdnsPdVJ1/Sv8c+0Bp03RmESLfBicjMm7ETUkdCA7jA06Tr51OXNvftMRGyRqhpOCLdtYRV7BpKQ8KR+uURkVC3/MJEFdGdyVyU6S4LHEhG75Do2TRVl78pqUobVbnbXHwxX6ZK3XXeGQuBQt8gZHDf1SwDOJ2TvSJS0ktS1dfFDdyLG/fFO5bRwF7vKCXZQU15DyIu+GAuaGNFjJNx9aybOrawD4dHezw9YkBy1dHiD2loV2M9Qj5UU+rPObRLlVO/CRqXMg7p6svLOtgZ89uz6ma7SjMTT9H4/ds+hTWuT9fsm3HwlsTKzimhsxi7eCadAMjYr186o9urXjJlJa5Ns9Xg6H2wjYIb30+SUPvLvTmZtbSK/PT5tCHZ9O4vdLer1+p82wBN2qjIzdw7NTWuT7o8qaGy+u2d+3Qbab6O71M/uXzoxFV41f/XPjgBnWimS9kEQrR7oNqS5a5BWjy8yCVSqrhSoo8Iwe+2Sv0yZEjfNPS2MWLfJhUMEz6fVF16TXofjIKKDtYYnm/V371Bqqrn/FemMGsWRTHR1RbBCv8ONNebTIh0GFTLvlYBsHDnc5bYa7UKxG/I8HV0R13nMObjBSn0yL9CnCDc+tc9qEPlJO5BUr4xHZfagj4jkqe6mqoOp7b+2O7CUnA4o+Xsd4YsW+vs9bDg7eGEYPobSUoCC293ipd9vYacWUTErJGbe95bQZGo2jLLzrPUfvn7ILlC26570hPWC15FJdGtvDN+X9EprDDVG1Gd3aiR9VRp5p4iPlPPkgSZNvh7BzZ0M7Vde/wvb6dvvsAf60fHvf5+tjnBnpOEnz4pOHlq7+Fbn9z3fD/hbb75lMpKzID2ZzbZvTJsTMy+sC4+lX7GoKHLBJwP745ra+zz1e9+xR6lQndzLvXfDq+lrm/OpN1tYctvW+/VuQD32429Z7m0VPhnKIwePTFQtvaxLFEC/25FuX8/72GDdPTwB2d8Qlkg+qA8+rrtW+ETi9Pj/zb1lq2/2SHS3yYejuVWTaeZJVNv0bE0lmOgCbB+31q1EPnyJbNcaLHl1jMWY8dDs6oI5qyiVxfk5i0zVJRGuXGp37qpJyIu8mVPCUD7V7nDYhKvZHGW+3w8vSocDEsnRzndMmKI3lIi+EWCiE2CqEqBZCXG/1/ZKdWETmqDMdUI9tdcnRYb3nUGfggwKjaxQwIWHoCkt9LBV5IUQ68GfgfGAGcLkQYoaV94yEmwqY6ky/6XWnTYgZp/PH4JDgox/t5sw/vu2ILUGGeiSDn5cTzy/avSBeW19LexTr8FjBoQ7nWrxWe/ILgGop5U4ppQd4Elhs8T0tw2kBGMxR684rZqBHwTXT71lezTOrapw2YwBDDan7xYsb2dUYeWmLRJLsznmXx0dNc+eAY1sPtvHdx1Zzw3NJNq8jAVgt8qOBff3+rjGO9SGEuFIIsVIIsbKhQa1Nlp0g2jG0LV29jnkl/UnG5vpP/rE27HdOpKe584iXp+IOZYMtesWYnyGlVHLZ5P/4+wpO+f3A5TSCZWWw+KcCVot8qBw7wN2UUj4gpZwvpZxfXl5usTnJJ0pbDoaJeavltPeRzGO+wZ7G0OA8+N8vbLD+pjEQ6RFc/fhq7lq6TdlN0T8JTg5UlJ5eP397f5dtQ0GtFvkaYGy/v8cAByy+p6u45ZVNA/7+9sMrlQyDBNlR38Gf36p22oyk5vY3tyo1VT+UFN21dDvff3x1VOdaTSS/zeeXXPrXDwE1Ipp/equaX7+8ieeN5aMPtffQauHWmFaL/KfAFCHEBCFEFnAZ8JLF97QMJ/LH4Hsu3VzHOpunkMfCV+//iNve2KpEKCke7HjHQwmNlJJ7l1fzxT+9b4Ml0dHZE3rpCic7E2PB1G5rFhBc66fTEygjx9+ylBN+s8yy+1m6CqWU0iuE+D7wBpAOPCil3GjlPZOdaMIdX77vI646fZIN1kQm0IdwxOZuF61l4wTBJ6mCxxlkqD6MwSRZNNRR+r9jKysiy5callK+Crxq9X1Sjb+9vzPySQ4QzLh6eVrzrN3nTIttsFBvq0+OuRDRoEKfnN0m6BmvijF4dI0CeVKTYIYSmv51Y01z8m39qHrVrpLvYZcjpEU+BrR3ejTJPppGZd7eWu/IfaN9o6qUBzWsUBct8ooTLgMrUr6Sdi304FjvwajyXMHZzbujQZUx/WpYET1eY+ikXfv7apFPUhQpX2FRSCtDcnWI4X92EW1FooqnrDqxlAWVnugdS7ZFPikBaJHXWMLDH+x22gRNgghX14Sc6aiSioZApfkHdqFFPgYUz79KcbtNXkoyEq3n6VQ4JNq7JmN5UG3DkUPt1u+opUU+SVHGY1I8bKQxT7i6JhnCSb0+P+fe8Y7TZoTljNvetvweWuQ15lC/nMeESqOFVBHRcPW4Kh2vQ1HX2k1tS7fTZoSlzYaZ4VrkFSdsPFT98qVJUgZnOTWqmvCEG+G1ao/aC5XZRcqJvLk9XhNnh2vQlY1lhPKU71m23QFLkpNL//qR0yYoQcqJvOpCPbhch4+HWm9LNGiNt45Qz/bBD3bZfl9V8lo8JENIyWpSTuRVJ9kKVJKZm1Qk47NVqU8D1OnXcBIt8oqj86i9qPS8QwmUSvapgmoVi2pokY8BnZk0miOoMuZ8qKU1dLhGi3zSokreVcQMVxJKoOx472pIdwwM8Uye+nRf+C9TBC3yihFtIe71JV1R1CQAHa6JDT0aKQVFXg+h1CQLToVDdOvMXaScyGsSiypho1RBP+8QzpZ2voZEi7xiJFtrIdnsTXb089bESsqJvJsLiYuTpkkSHClfunUzJCkn8prE0uP1O23CUZiZAKMnz2jchhZ5xTATc62ub0+cIUmMyqsOgrktE3VM/mg8CjoaKqFF3kU88tEep01QAtV9cTOT6nRDQxMrKSfy2hPSaJIbPfM8NlJO5M2gvajkQOWdgDQau9Eir3EdnR5f3NfqilzjNrTIazQ2Y6bjVXV0HakeWuQ1Go3GxWiRjwHd4aNJBG7ORyt2HXLaBM0gTIm8EOI2IcQWIcQ6IcTzQoiSft/dIISoFkJsFUKcZ95UjUajOv/38V7L76H7TWLDrCe/BJglpTwW2AbcACCEmAFcBswEFgJ/EUKkm7yXxgWoPqNUbes0mtgxJfJSyjellF7jz4+BMcbnxcCTUsoeKeUuoBpYYOZeqYJ7u+Q0Qdzc8apRj0TG5L8JvGZ8Hg3035KlxjiW1NjhhGpPUqPRJJKIIi+EWCqE2BDi3+J+59wIeIHHgodC/FRI/RJCXCmEWCmEWNnQ0BBPGmzjE92pZBrFozUa9Kxwt5ER6QQp5TlDfS+EuAK4EDhbHgm41gBj+502BjgQ5vcfAB4AmD9/vtIS8M2HVrL71gucNkOTwtjRp6ErYndhdnTNQuBnwEVSys5+X70EXCaEyBZCTACmACvM3EvjDlTXDy1w6qNfUWxE9OQj8CcgG1hi7Cz/sZTyKinlRiHE08AmAmGcq6WU8c81TyCqd3qpbZ3GaYSOpWhixJTISyknD/Hdb4DfmPl9jUYzENWHoGrUI+VmvKo+21Bt68yjukipnj80mlhJOZHXaJIZO6ogHRFyF1rkNZokoq3bG/kkkyje2NLEiBZ5xXC7E6X1Q2MW1UN+qqFFXqPRaFxMyom86kMoVcftXpTLk6dJQVJO5DXOokVUo7EXLfKKoTVQ4zR6dI27SDmR9/r8TpuQ0uhx6BqNvaScyPsV1xjtRGmcRvWQmuLmKUfKifzB1m6nTdAojBYQjdtIOZHXmMOsl6e6l2gHXr8OGWrsQ4u8Ynh0n4HrUT1kuK+5M/JJmqRBi7xi+FRXAI3r0XnQXWiRVwxdvjROo0Nq7kKLvGL4dLzWWbTCKd/5rF9RbGiRVwy3h+R1AVUfty9dkWpokVcMv+IFTG3rNBrNYLTIK4ZfB+U1Gk0C0SKvGG6XeL2sgUZjL1rkFUP1cI3b0U9fPwO3kXIiP3/8MKdNGBIdrtE4jfIdr4qbpxopJ/Kq41O8gJkVAMWTp8H9k6GUr8QSTMqJvOprZWemp9wrUYoUK/8hGVGU47QJmgSiFUUxFkwoddoETYpz6pRyp03QJBAt8orh9j1otaMMI7Wn7Cip1lpLOZGv0AVM4zCLZo9y2gRNCpFyIj+1otBpE4bE7ePIU63TKxS628Ucbi8jiUZnN01MuL14aQFxf0Xs7tQdjRZ5jaYfLtc3TQqSEJEXQvxECCGFEMONv4UQ4h4hRLUQYp0QYl4i7qNJflTXUNXt02hixbTICyHGAucCe/sdPh+YYvy7Evir2ftoNBpNInB7OGowifDk7wSuY6ATtBh4RAb4GCgRQughBRrlSbHyr0kBTIm8EOIiYL+Ucu2gr0YD+/r9XWMcC/UbVwohVgohVjY0NJgxR5MEqC6iuuNVfVTPQ6qREekEIcRSYGSIr24Efg58IdRlIY6FfDVSygeABwDmz5+vX5/iT8D1Bczt6dOkHBFFXkp5TqjjQojZwARgrQgsCDMGWC2EWEDAcx/b7/QxwAHT1mqSH8VF1A7zVK8oFTfPNG5P32DiDtdIKddLKSuklFVSyioCwj5PSnkQeAn4d2OUzYlAi5SyNjEmm0M3xzUaTSoR0ZOPk1eBRUA10Al8w6L7aDQJxd0rB2lSkYSJvOHNBz9L4OpE/bbGPSjfktIq73pUD5clGj3jVaPRJBUpptGm0SKviQnlPXFgfFle3Ne6falnN5Bqk5nMokVeMVTPvmbLlx3lU/Xdtdz+jq3GrHl2OCrjSuN3NBKN2qUhBfErXsIUN880aTY48qq/Y9WrIeUfH/bko2jRIq8Yqu+hnAzhGtWb84qbp34eVP0BotYz1CKvGKp7eWYzr9qps2ejd9XfsU8lhQqB6XCNDclT6R2nnMgr9OxD4le9gKn+AFG/IlFJAEKh7TOPSiamnMirnkEU1/gEeFFqJzDNBlfe57f8FqZQ/BUpbx+opTMpJ/KqN0VVyhyhkIoLFGCqJspIs75IqF7RqZ4HVbcP1LIx9UReoYcfCtUFIBk6Xs2QnWl9kVBJAEKhuB+UJJ680xYcIeVE3kwGWTy3MnGGhMHtTXlbVnk0cW1pflbC7AiHGQG4ZF7IbRkSipl+oWF5mQm0JHlRyVlLOZE3E67JyUhPoCWhMePlDS+wQ6DUybzhcHNIrjjXehE1Y9/oYbkJtCQ0Og/GhhZ5xTCTge3oNFT88QEm37EN6VNdo8w8PjuWhTA9jNeWIZSxnV+UY9WCwCko8maaUXbEo81kwHQbptl5/ebiSak2RjkUZuzTzw98JvOgHcT6DFu7vazZd9gSW1JO5FXveDVTwGwReZ/azw/MefJ2VOSqtyZVd4S8ij8/iK8yfmPjwcQbQiqKvAknwI76wYwAZNgg8r2q9wyjviequHnKh+TMOhp2VETx5MEsixbWSzmRV39GqdMWDI1ZL8qOAtareGtD9UrIVEvIhqQlgycfl8hnaJFPCGbCNXZkrV4T8UZb7EsCT97jNfEMbXiIqldCZvtdrCYZYvLxvGOrOl9TTuTjEalheZmMLrF+aBjEL1B2jOGH+JvKXz9xfOCDLSIa+zP8xuerEm9IGDxx2HfGMeW2DJ+E+PPgOdMr7PHk48yDwXdstY0+v4yrNVSSZ80Q6JQT+Z7e2DNwYY59Ezx6BhWwaON0di1fHa+XV2DhELH++P0yruZ8bmZgDoQdPnZPry/ma4YZAmDHJJvBeTAzXaHF0Yk/XDPMIhEdTLyVpFX2pZ7Ie2MvYIWGQNnhpRyVQWIoXyqHGspsmEkK8XnJYJ8AQHw2luRl2rIMMhydB0XUN7bHwO44KkmAHBuWrID4Rb7EotnCKSjyfiaV58fU9LXq4YfC4/WTJuAPXz4WiL7YRF8QzRFsKj/+7RNiui43yx5POSig/3FyVUzXDS+0UeS9fs6eVtHXeogGOyuhHq+fwpwMrjlnChBdHgwO37WjJdRliPzT3zkJiH0PAKtt7PEF7Lvl4lkxXWfVkhopJ/Ier5+RxTkx7cE4vCAbsGdkiMfnJzsjvS+mZ5f3Fi2dHi9wRLSjxa4NsoNe1ITh+TFdV16QA9jXWsvOTIspbhtcE8aWcJLXT3ZGWt+ziCYPDsvLsi2vdnkCIppn5MHMKFcOtTsPxjpaRodrEkS310dWelpf59wFs0cB8N51Z4a9piw/27YM3N3rIzvziH1nTaugKCeDl39wSthrgq0SOyqhDkPkC7KPxNjPPKacP35ljuX3jobOnoECAPDvJ43n304cxxnHlIe9bkRRtuW2Ben2+sjOSO9rddzx1TlMqSjgmatOCntNSV6Wbf0uPb0B+4J5cPGc0aSnCZZee1rYa4LhODv6DIKefI7REspIF1w0p5JHvrnA8ntHQ7fR75fdT+S/dcoErjxtIrNHF4e9LlbHKVpSTuTbu70U5mT2FbBrzpnC5psXUtSvc/WKk8YPuKa80BAAG9yotm4vRTmZfR2ulcW5rPvleczqlzkuXzBuwDVl+fYJQEdQRPuJ/N+/sYAvHz+GmZVFAPz43Klhr7daA1q7e4GBneU3L57FLRfP7lvb56I5R49EKrKxogy844y+BeUuOHYUS649va9zetrIQk6YUDrgGjtWxwzS2u2lMCejz7OcOrKQHb9dxOSKwr5zLp035ij77MqDQU8+KKKZ6Wncc/lxnDa1vG8U3FWnT7LJmqMJ5sGifiHhX1w4g58vms68cSUAfOX4MSGvtQJ7hjwoREtXL0W5GdyyeBa3vLKZcWV5ZGekk5uVzn9fMJ3zZo7k452HAJhUnk9DWw+L51by+Io9ttjX2tVLYU4Gi44dxeq9zXz3jCOZ9fnvnUxhTgbdvX6eWLG37/hpU8tp7eq1xb6OnoAnX5ybyRfnVA5Y+vaur81l/+EuzjimgtuXbOs7Prok17aWUFt3wL6inAxK8jIHVN63XjKbf6yq4XtnTOKltQcGXGeXQEkpjXecyX3/djxvba0n21jddNrIIh7+5gLmjx/Gfe/s4JNdTX3XjTFWd7QjnNTW3UtRbiaXzBvN/sNdfHX+EUF67nsnU16Qzf7DXTy7uqbv+MmTyli/v8V644BOQ+RHFudw6pThAwT9sW+fQFu3l5mVRdz3zo6+432OGta3NoJlsSgnkxvOn9YX7gX48XnHcNb0EZw+tZx/rKoJ9xMJJaVEXkpJa7eX4txMTp48nFd/dOqA77996kTgyBoSp0wezq8WH+k8sSMe2trdS1FOQJz+8OWBIZDjxg0DYINRmLLS03jueyczbWQh1z2zzhYB6OjxIgTkZaZz7+XHDfhuyohCpowoHHBsxY1nk5eVwT8NUbXaU27r50V99otzBzyTiqIcrj5z8oDzX7/mVEpy7fOSOzw+/BKKcjOYX1XK/KqBHvvpUwMhpf2HuwD41UUzOXZMMePL8m3rXG/t9jK6JJeygmx+edHMAd/NM/LgAcO+McNy+cu/zmNWZTFX/d8qW+xrMRyhzPQ0Hv3WwAEAVUZfTP+Z7S//4BQqS3J5brU9ohp0NIpzM/jOoBZFUU5m3zsO8uoPT7UsVAMpFq7p9Pjw+eUA7y4UwZr31ClHXoZdnTatXV6KcoeueyuNJumi2SOZNbqYjPQ021zRxg4PpXlZpEWxTk5muqCiMGdA/N5qggWsMCcDIURYO8+bOQKAY0YUMrI4p++45eGkfl7eUMyqDITnZo0u7qvc7aLVaO0ORTAPnjtjBMeOKSEtTdjWWmvq8EQckht872NLc5k1utjmcNfRIcNQTB1RAMD0UYUxDxSIhZTy5Js7PUDkIZGL51YyuaKgL8ZsJ4c6PMwdWzLkOaX5Wbz6w1OZXFEw4Lgdnnx9a8+Apm843v3pmQO8E7vCIYc6egAYFqFQ/+HSOfzyIu8R79hGgYLIefCKk6s4eXIZ00YeyYMC61dRlVLS3OmJ2LoZW5rHKz88ZYB9geuttC5AU4cn4vsFeP9nZw4Q2uC7tnoV0Kb2wDuONEz7qStPorPXZ3kLzbQnL4T4gRBiqxBioxDiD/2O3yCEqDa+O8/sfRLBgcPdAIwqHnqJAiEEs0YXD3j4menC1Joo0dDj9dHY3tPnJQ3FjMqiAUO0cjLT45roFSsN7dGJ/LiyvJDnWb221IHD3RRkZ0T0lIvzMkPmA6s1KhjmiPSO09PEUQKam5VOt8fad9zS1Uunx0dlSU7Ec2dWFg9Y3jorw5482NjeQ1l+5Dw4ZljeAKHNN5yOrjgnU0XLgZZuyvKz+kb/hGNYfpYty6WY8uSFEGcCi4FjpZQ9QogK4/gM4DJgJlAJLBVCTJVSWp8DhiBYwOLZoiw/O4N2o9PRKg62BCqhaArYYApssA+grqWbSeVlMV8X9Orjna0YLbUtXYwqjv35BQtkPEsOxEK0Ih+KguyMviGsVhF0hOKzL532Hmufn5SSvU2dnDQp9jwYHBHWYbGNBw53xfX8rMKsJ/9d4FYpZQ+AlLLeOL4YeFJK2SOl3AVUA44PYt3b1AkEhiXGSn5WRt9EIKsI2hdXJZQVGHXjtXCVyLbuXg62djOpvCDyyYMILg3R3m3tM9xzqLNvJEosFGQF7Gu12L69TV3kZKbFtcxDXlZ638gSq9jb1AEQl4dpRxlpaO+h0+Ojqiz2GHbQk++w2Bna29Rp24KG0WBW5KcCpwohPhFCvCOE+JxxfDSwr995NcaxoxBCXCmEWCmEWNnQ0GDSnKHZeKCFicPz4+rJDnjy1hawzbWtAEc106MhP9vIwBaKwI6GgABMqYhd5PMNEbWytdHj9VFd3870UbE/v7Q0EWgNWSzym2tbOWZkUVxx2PzsDMsFauOBVtIETB00Sioa8rMz6PT4LN2zYdvBdoC4HI08Iw9a2Rrq6PGy+1BHXHnQKiKKvBBiqRBiQ4h/iwmEe4YBJwI/BZ4WgdwbKgeHfPNSygeklPOllPPLy8PPSEwEG/a3MiPOztSC7Jg7+dEAABx9SURBVHTLC9iG/a2MLMqJayRAvtEUtdKTWmvsQRlPBg5O9GmzUES317Xj9UsT7ziD9h7r5hv4/ZJNta3MGBW7gII9nvzGA61MKi+I0xEKOhrWveM1+5oBmD0m/MzRcATt67TQWdtc24qUODJoIxwRY/JSynPCfSeE+C7wnAzMLlghhPADwwl47mP7nToGOBDiJ2xj76FO9h/u4lunTIjrequ9KL9f8uGORj4/eXhc1+f3xRuts/GjHYcYXZLL2BjW/QlSmB3oALPSk3+/uhGA+eNLI5wZmoIca/s1NhxooaWrlxMmxB5PhkBryEoB9Xj9rNjVxBfnjIrr+vx+MW+rludesbuZyRUFca2tb4cn/0H1IYSA48YNPULOTsyGa14AzgIQQkwFsoBG4CXgMiFEthBiAjAFWGHyXqZYvqUOgLOnV8R1fb7FnV4bDrTQ2O45aqJEtBQYXopVIaXuXh8fVDdySpyVUNCTt7ISWr6lnumjigaMe4+FguwMS1sab21pQAg4dUp8zzAvO93STsOVu5to7/Fy5jHxlZHgfAirKsq27l4+2tHIWdPiLcOGJ29ha2j51nrmjCmhrMC+tZAiYVbkHwQmCiE2AE8CV8gAG4GngU3A68DVTo6skVLyzOoapo0sZHwcHTZgfcfmM6tqyMpIi7uABb2UTosK2NLNdbT1eLkozh2o8vsqIWvs23OogxW7mlg4c2Tcv1GYY53I+/2SZ1fXsKCqNG4BKMjOpL3ba9m0/GdW1VCQnRF/azLL2pDhK+tq6fVJzovzHQdbGsFZ0Ymmur6NtfsOs3BW/HnQCkwNoZRSeoB/C/Pdb4DfmPn9RLF6bzMb9rdy8+KZkU8OwxFP1EdxXmInCrd09fL86v0smjUyqkkeoQh6UW0WiKiUkv99bxejS3I5cWJ8oYbsjHSyMtIsE9FHP9pDeprga58bG/nkMBRkZ1BrDGNNNG9vq2dvUyc//kL4xdsiUZwbWFivu9ef8GnwDW09vLyulssWjO0Tw1gJXmdF57XfL/n7B7uZPqqob5GvWCnMziA9TdBi0TpPD3+4h6z0NFsXH4sG1y9rIKXkD69vpSw/i0vmxf/wgzFAKzLI/e/soK3Hy5Wnxb9ynpX2vbu9kTX7DvO9MycNmPwSK1Z1bNa2dPHox3u4aE5l3KEaCHjyVgnUH9/YxtjSXM6fFV+8G6x9x/cu345Pypg3W+lPcBavFfY99sketta18Z3TJsY9Q1QIQVFOhiX27T3UyZOf7uXS40crFaqBFBD5NzYe5JNdTfzonCmm1lApsaiA7W7s4MEPdnHRnMq4R4VAvwLWmVj7unt9/PKljYwrzePLJj0Uq2Lev39tC1JiyksGIxxiQUvoiU/3sqm2lWvPnRrzRhL9sUrktx5s4/FP9nL5grFMjGNoYpCgfYcTbN++pk5+99oWTp0y3PSG9cW5mbR0JfYdSyn57aubSU8TXHOOuTxoBa4W+Ya2Hn7+/AZmVhYdtQZ7rARF9HCXJxGmAQEP77pn1pGZnsbPF0039VsFRlM0kfYB/OWtanY1dvCbL83qWxI3XopzMxO+JPLrG2p5Yc0BvnvGJMYMi33UT3+Co2sSubbJnkMd/OaVzXx+chmL54ScKhI1fXmwM3Hv2OP1819PraEoN9O0QB2xL3Hv2O+X/PSZtaQLwe8vPdb0Oi8BkU9sHnxhzX5e33iQH509lRFF8bckrcK1Ii+l5Ppn19He4+Wur80lM91cUq3IwA99uJsVu5u46cIZpsIMEGiKluRmJtS+z/Y28+e3d/Cl40YPWJEzXoblZ9GUQPvq27q54bn1zB5dzPfPmhz5gggEW2uJqoh8fsm1T68lPU1w25fnRLVy51BY4cnfvWwbm2pb+d0lswesex4PuZnpZKWnJdTRePTjPXy8s4n/vnB6QpYKKM7LoiWBleT+w13c9MJGPlc1jCtPm5iw300krhX5pz7dx7It9fxs4bSj1jiPh2JjVb5ENUW31bXx+9e3cNa0CtNhkCDFeZkJs6+9x8s1T61hZFHOUWuKx8uwvMyEeaFSSm54dj2dHh93fm2O6Uocjuy+1JQgG+97Zwer9jTz68WzEiNQCRb5VXua+OvbO/jK8WPiHrHSHyEExXmZCQsZ7jnUwa2vbeH0qeV8dX78Her9SaQn7/dLfvL0WvxScvtX5prqr7ISV4r8nkMd3PzyJk6eVMY3THQk9aevgCVAALp7ffzwic8ozMlISBM0SHFu4grYr17ayL6mTu782ty4Jp6EYlheVt9Su2Z5YsWRSrz/tnRmCI5sak6AjRsPtHDX0m1cMHuU6ThykKIEinxHj5drn15LZUkuN31xhunfC5Ko1mQgTLOOjHTBrZfOTmAZSVzH698/3M1HOw9x0xdnMK7MXKjQSlwn8j6/5MfBJvJXzDeRg2RlpJGflZ6QDPyH17ey5WAbt315TlTL9kZLSYK8lFfW1Rrb5E1mwYT4Zo+GYlheFm3d3r4NouNlZ0M7v355E6dMHm5qNMhggouGma2Ierw+rn1qLSV5Wdxy8ayECVRhdgZpIjEi/+uXN7G3qZM7vjo3obNTS/ISkwcf/mg3K3Y18YsLZ0RcGjwWinMzaU3AXINgS/yc6RUJa2VYhetE/v53d7ByTzM3L56Z8JXgSvKyTIdD3t3WwIMf7OKKk8ZzZpwz98IRsM+cQB1s6eaG59Yxd2wJPzpnSoIsC1Cab75fo9cX6CjMzkzj9q8mrhKHfp68ydbaHUu2sbWujT9cemzc8x5CkZYmKEpARb5kUx1PfrqP75w2KaGVOATCmmbLyJ5DHfz+9S2ceUx5wsecl+Rm4fNLU6OoPF4/1zy5hsLsDH53SeJa4lbhKpHfXNvKnUu2sWj2SC6ea24kQyiKTDZFO3q83PDceiZXFHCDydE0oSjOzeRwh7kCdtOLG+jx+rkzAZ3VgynJM/o1TIjo/e/sYG1NC7/90uyEj2QoNew7ZMKTX7WniQfe3cllnxub8EocAjaaaWk0d3i4/tl1zBhVxLXnJn64X4nJfhcpJT9/fj2ZaWmWCOiwBLTW7lm2va+zOpEtcatwjcj7/YHMUZSTyW8uTlwMrz+l+Zl928vFw+1vbmP/4S5uvWR2xF1j4qE0P4u2Hm/cu/O8vqGWNzfV8V/nTrVkz8lSkwVsd2MH9yyv5oLZo1g0O/5JReHIzUonJzMt7pi8x+vnZ8+up7I4l/++MHFx7v4ML8ymoS3+PHjra1s43NXLHV+bY2rMfjhK87M41OGJOxzyzKoaPqg+xPWLppkecRaKoCjH+wy317Vx3zs7uGTeaL6QgM5qO3CNyP9j1T4+23uYGxZNT2gTuT8VhTnUt8aXOdbVHOahD3fxbyeOY35VYpvIQUYUxZ+B23u8/OLFjcysLOLbca7UGYngMNTmOFpDUkp+8eIGstLTEtpROJiy/Gya4mwNPfThLqrr27l58UzLNi8vL8ymoT2+PLhqTxNPrdzHt06ZENeeBdFQUZiNx+unNY4JR00dHm55ZTMLqkq5/HPm5rWEo7wg/jISzIP52RncaEFL3CpcIfLNHR5ufW0LC6pKuXRe4sM0QSoMLyoeL+V3r26hND+L6xZOs8CyABWFAc+nPo4M/OD7u2ho6+GWi2eRkeAwTZCgJx9Pa+iNjQd5b3sjP/mCtRNOhuVnxhWTr2vt5u6l2zlrWgVnTx9hgWUByguyaYjD0fD5JTc+v4HK4hx+dHZi+1r6E/SU69tiXwPovnd20Nbdyy1fmpXQvpb+9HnycVSUL645wMc7m7hu4THKLV0wFK4Q+Xe2NdDe4+Xmi2da2glSXpiNx+ePuePrw+pGPtp5iO+dMTniBtNm6CtgMYrA4U4P//PuTr4wYwTHjRtmhWkADC/IJk1AXYz2+fyS29/cxuSKAr5+UpU1xhkML8iOS6Bue2MrvT7J/7OwlQGBd9zW46UrxuVyn1tdw5aDbdx4wYy4FyCLhngdjfrWbh7+cDcXHzc6rl2poqU0P4s0Ebsn393r47Y3tjJ7dDGXWdTKsApXiPzFx43mvevOsqwJGqSiKPYMLKXk9iXbGFmUw7+cYG3mqOgL18QmUv/z3k7aPV6uNbn2SyQy09MoL8zmYEtXTNe9vO4A2+vbueacKZZPOBlVnNO3oXq07Gxo57nVNXz9pPFxL2UdLcGKvDEGT7S718ddS7czZ0wxi2ZbG0cOhgxjrSj//FY1Pr+0tJUBkJ4mGF4Qe7/GY5/sZf/hLq4/f5qyk57C4QqRByzppBlMRRye8jvbGli1p5nvnzXZks7W/pTlBzzlWCqhQ+09/P2D3Vwwe5TllSTAyKKcmJbz9fr83L10O9NGFrLIxAqO0TKqOJfGdk9Mndf3Lq8mKyONq06PfxXRaDkSDon+HQcF6mcLp1k+3K/PEYqhjNQ0d/L4ir18Zf5YyytJMPo1Ynh+bd29/Pmtak6ZPDzutfadxDUibwcVMcYbpZTcsWQbY4bl2jJhIuilxFLA7ntnB929PttWzxtZnENda/Qi/8KaA+xs7OCac6ZaFqftT9BZqGuJ7hlW17fz4pr9XHFSlS3D6Y50HEb3DIMCdeqU4Zxsg0AVZGeQl5UeUyX0p+XVCAQ/SMD6Q9FQXpgdk33/+94umjo8/PS8Yyy0yjq0yMdAsMPvYJQitWRTHetqWvjh2VMsGa4WihFFOVHbV9fazSMf7eFLx41hckX8S8zGwqji3Kg9+V6fn3uWbWdmZRHnzbSuM7M/lcbsytooQ0r3LNtOTma6bYtTBSf41TRHZ98jH+2xXaBiyYO7Gzv4x6oa/uWEcQlZ3ycaKkty2X84uufX0tnL397fxcKZI5kzVp19W2NBi3wM5GdnUJafxb6mzojn+v0BL37C8HwuOc66ET+DGVeax94o7AP74qD9GVGUQ1u3N6oZh8+uqmFvUyfXnjvVtlmFQU8+mopoW10b/1x3gCtOrrJttEVJXiaFORlRveP2Hi//895Ozp5WwbFj7BOosaV57D0UXR68Z9l2MtMF3zvD+lBXkHGleTR1eKLaBvDBD3bR3uNN+OxvO9EiHyPjy/LY3Rg5A7+6oZYtB9u45pwplg1JDMW4sjxqmjsj7kVb09zJEyv28tXPjbV1caXKkoCI7o/giXq8fu5dXs2csSVxb9wcD332ReHp3b1sO3mZ6Vx5qn1LzAohoq7IH/1oD4c7e/mBjZU4wPjSPHYf6og41Hh7XRvPG6GuChvXYR9fGsjvkZ5hS1cvD36wi/NmjmD6KOv7q6xCi3yMVJXlR8wcPr/kziXbmFJRwIXHJmYFwmipKsuj1ycjeqL3LgvEQb9/pj1x0CCTjJ2Hdja0D3neUyv3sf9wl61ePAQ2RK8ozGZXY8eQ52052Mor62r5xucnWDb5LhzjyyJ7yp2egBd/+tRy5tocZhhflkdbtzfiEiB3LQ1Ukt+xocO6P2MNkY/UIn/4w920dXv5wVnJ68WDFvmYGVeWx4GWLrp7w4++eHHNfnY0dHDtuVNtH241rjQwOmH3ofAitbuxg2dW2xsHDRJcLmHnECLa3evjz8urmT9+GKdNsX80w+SKAqrrh66E7l66ncLsDL59qjWzg4dibGkeNc1dQ+5g9X8fB2LxP7TZiwf6RsgMlQc3HmjhlfW1fOuUCX2T5Owi2HLdM0RF2dYdiMWfM30Es0YX22WaJWiRj5GqsnykDO8F9Pr83L1sOzNGFSVkI4ZYqRoeyMC7hxDRu4Nx0DPt9aAg0K8xsiiHHUN48o9/speDrd22e/FBJpUXsKOhPWy4YcP+Fl7bcJBvnDKhb9E1O5lQlo/H56emOXQe7PL4eODdnZw6ZTjHj7ducls4qgwRHUrk71yyjaKcDL5lY6grSFFOJmX5WUPmwUc+2kNLVy8/PNvelq4VaJGPkWNGBmbjbTzQGvL7f6ysYc+hQGehHUP+BjOyKIeSvMyw9m052MoLwThooTP7UU6qyA/rKXf0ePnL29WcOLHUliF/oZhcUUBbtzfsMLs7l2yjODfTES8eYGZlwLMM944f+2QPje3OePEAVcPzyc5IY+P+0Pat3tvM0s31XHnaxIRtSBMrMyqLwj6/9h4v//veTs48ptzWDmur0CIfI1MqCsjJTGNtzeGjvuvy+Lhr6Tbmjx/G2dPt6yzsjxCC2aOLWVfTEvL7217fSkF2Bt+1cTTDYGaPLmFzbWvIkNffP9hFY7vH0jV+IjGzMtDJtmbf0e/4s73NLNsSECgrl6gYiqkjC8hIE6zff/Q7buk6MnHncxYthBeJzPQ0ZlQWsS6EfVJKbn1tC8MLsvnG552pJAFmjS5mW11byElvD7yzg+bOXtvmjliNFvkYyUhPY2ZlMWtDCMBDH+6mvq2Hn51v/czCoZgzpoStdW1HrW/y6e4mlm2p56rTJzkSZghy/Phh9PokGwaJQHOHh/vf3ck500cwz8I1dCIxe0wxWRlprNzdNOC4lJLfvbqFsvyshO5IFSvZGelMHVHIuhCOxn2GQF1/vnOVJMCxo4vZsL/lqH6Dt7c1sGJXEz88e7Kla+hEYlZlMb0+yaZB3nxdazf/894uvjinMmnHxQ9Gi3wcfK6qlHU1LQP2Uz1wuIs/Ld/O2dMqHPOggiyYUIrPL/lwR2PfsV6fn5te3MiIomy+6aAHBTBvXKDwfLTj0IDjt762hU6Pz/GZhdkZ6cwdU8JHOwfa98Ka/azY3cRPzzvGUYEC+PzkMj7d1TxgvsGuxg4efH8XF8+tdLyzcMGEMjo9Plbtae471t3r4+Z/bmJ8WZ7ji3ydOLEUIQLLjvTn1te24PX7+ekXknN2ayi0yMfBF2aOwOuXLNlcBwQmPt34/Hr8En550UyHrYMTJ5ZRmJ3Bq+sP9h174N2dbK5t5VcXzSI3y9o1dCJRVpDN/PHDeHldbd+xD6sbeWrlPr59yoS+fg8nOXfGCDbsb+3rnKtr7ebXL29mzphiJfb0PGvaCDw+P8u31AOBYbs/e3YdWRlpluw6FiunH1NOVnoar64/8o7vXLqNXY0d3HLxLNtmgIejrCCbeeOG8fqGg30d7Ms21/H8Z/v57hmTld6YO1a0yMfB3DElTCzP54F3A+u+/P6NLby1tYGfLTymbwyuk2RlpHHR3Er+ufYAuxo7eG19Lbe/uZULjx3Fwllq7GZzybwxbK1rY8mmOqrr2/jBE58xqTxfmZmFi4+rJDNdcPfS7TR1eLjy0VV0eXwJ31c2XhZMKKWqLI+/vr2DLo+Pm17cwIpdTfzyizMtXW8/WgqyM1g4ayRPr9xHbUsXz66q4f53dnL5grGcOqXcafMA+Or8MWw52MbrGw6yvqaFa55aw/RRRbbPHbEaYXbX8kQyf/58uXLlSqfNiIqlm+r49iMryc1Mp6vXx+ULxvHbL81SZlPf2pYuFt39Hh09Pjw+P3PHlvD4f55AXpazYYYgHq+fC+99j50NHQgR2AD6qe+c2DdZSgXuXLKNu5dtByAzXfCXfz2ec2fYs4ZONLy+4SBX/d8qMtMFvT7Jd06fyA3nO+/FB9nZ0M5Ff/oAj8+Px+vnpIllPPTNz5Gd4WxLMojH6+dLf/mATbWBuPyoohyevuokxgxz3lGLFSHEKinl/JDfmRF5IcRc4D4gB/AC35NSrhABpbsbWAR0Av8hpVwd6feSSeQhsADZGxsPcuLEMi6dN1oZgQ+yva6NRz/ew+iSXK44ucrypY5jpb6tm7++vYN0IfjWqRMYVWzvxKxISCl5ZlUN1fXtfGneaFuWYo6VpZvqWL61ntOmDGehDUsxx8qG/S08sWIvE4bn8+8nVTkephlMU4eHB9/fRZrA1jWIEo2VIv8mcKeU8jUhxCLgOinlGcbnHxAQ+ROAu6WUJ0T6vWQTeY1Go1GBoUTebLUqgaB7UwwcMD4vBh6RAT4GSoQQ6rkZGo1G43LMBmivAd4QQvyRQIVxsnF8NLCv33k1xrFaBiGEuBK4EmDcuOTaO1Gj0WhUJ6LICyGWAqGGZNwInA38l5TyWSHEV4G/AecAoYLTIeNCUsoHgAcgEK6J0m6NRqPRREFEkZdSnhPuOyHEI8CPjD//Afyv8bkG6D+YeAxHQjkajUajsQmzMfkDwOnG57OA7cbnl4B/FwFOBFqklEeFajQajUZjLWZj8v8J3C2EyAC6MWLrwKsERtZUExhC+Q2T99FoNBpNHJgSeSnl+8DxIY5L4Gozv63RaDQa86g1M0Gj0Wg0CUWpZQ2EEA3AnjgvHw40RjzLXeg0pwY6zamBmTSPl1KGXBRIKZE3gxBiZbgZX25Fpzk10GlODaxKsw7XaDQajYvRIq/RaDQuxk0i/4DTBjiATnNqoNOcGliSZtfE5DUajUZzNG7y5DUajUYzCC3yGo1G42JcIfJCiIVCiK1CiGohxPVO22MGIcSDQoh6IcSGfsdKhRBLhBDbjf+HGceFEOIeI93rhBDz+l1zhXH+diHEFU6kJRqEEGOFEG8JITYLITYKIX5kHHdzmnOEECuEEGuNNP/KOD5BCPGJYf9TQogs43i28Xe18X1Vv9+6wTi+VQhxnjMpih4hRLoQ4jMhxMvG365OsxBitxBivRBijRBipXHM3rwtpUzqf0A6sAOYCGQBa4EZTttlIj2nAfOADf2O/QG43vh8PfB74/Mi4DUCSzufCHxiHC8Fdhr/DzM+D3M6bWHSOwqYZ3wuBLYBM1yeZgEUGJ8zgU+MtDwNXGYcvw/4rvH5e8B9xufLgKeMzzOM/J4NTDDKQbrT6YuQ9muBx4GXjb9dnWZgNzB80DFb87bjDyEBD/Ek4I1+f98A3OC0XSbTVDVI5LcCo4zPo4Ctxuf7gcsHnwdcDtzf7/iA81T+B7wInJsqaQbygNUEtslsBDKM4335GngDOMn4nGGcJwbn9f7nqfiPwJLjywisWPuykQa3pzmUyNuat90Qrgm3C5WbGCGNpZqN/yuM4+HSnpTPxGiSH0fAs3V1mo2wxRqgHlhCwCM9LKX0Gqf0t78vbcb3LUAZSZZm4C7gOsBv/F2G+9MsgTeFEKuMXfDA5rxtdqlhFYh6FyoXEi7tSfdMhBAFwLPANVLKViFCJSFwaohjSZdmKaUPmCuEKAGeB6aHOs34P+nTLIS4EKiXUq4SQpwRPBziVNek2eDzUsoDQogKYIkQYssQ51qSZjd48qmwC1WdMDZCN/6vN46HS3tSPRMhRCYBgX9MSvmccdjVaQ4ipTwMvE0gBlsiAnszwED7+9JmfF8MNJFcaf48cJEQYjfwJIGQzV24O81IKQ8Y/9cTqMwXYHPedoPIfwpMMXrpswh00rzksE2J5iUg2KN+BYG4dfB4qB243gC+IIQYZvTcf8E4phwi4LL/Ddgspbyj31duTnO54cEjhMglsC/yZuAt4MvGaYPTHHwWXwaWy0Bw9iXgMmMkygRgCrDCnlTEhpTyBinlGCllFYEyulxK+a+4OM1CiHwhRGHwM4E8uQG787bTHRMJ6txYRGBUxg7gRqftMZmWJ4BaoJdADf4tArHIZQS2V1wGlBrnCuDPRrrXA/P7/c43CezMVQ18w+l0DZHeUwg0PdcBa4x/i1ye5mOBz4w0bwBuMo5PJCBY1QT2TM42jucYf1cb30/s91s3Gs9iK3C+02mLMv1ncGR0jWvTbKRtrfFvY1Cb7M7belkDjUajcTFuCNdoNBqNJgxa5DUajcbFaJHXaDQaF6NFXqPRaFyMFnmNRqNxMVrkNRqNxsVokddoNBoX8/8BzETdBYWbPu8AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(s[4999:, 0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "san_bool = [\n",
    "    1,  # leak channel\n",
    "    0,  # voltage-gated sodium channel\n",
    "    1,  # HH-type delayed rectifier potassium channel\n",
    "    0,  # fast A-type potassium channel\n",
    "    0,  # slowly inactivating potassium channel\n",
    "    1,  # voltage-gated calcium channel\n",
    "    1,  # calcium-dependent potassium channel\n",
    "    1,  # persistent sodium channel\n",
    "    0,  # inwardly rectifier potassium channel\n",
    "    0,  # AMPA receptor\n",
    "    0,  # NMDA receptor\n",
    "    0,  # GABA receptor\n",
    "    1,  # calcium pump\n",
    "]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "san_x = anmodel.models.Xmodel(channel_bool=san_bool)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "san_x.set_params(anmodel.params.TypicalParam().san_sws)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "s, info = san_x.run_odeint()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x152406abd0>]"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO2deZwcVbn3v2f2yWQmk2WyTsi+EMIWQtj3NQjEBRXkiiJXBNH3qter4vZe7xU3RH0BBaOgoqKCoCIgIZBACCSESci+TmYmmX0ms+893X3eP7p60ky6urq76nTV9Jzv5zOf6a7q7uecqnN+5znPWUpIKdFoNBpNepLhdgI0Go1Gow4t8hqNRpPGaJHXaDSaNEaLvEaj0aQxWuQ1Go0mjclyOwGRTJo0Sc6ePdvtZGg0Gs2IYuvWrceklCXRznlK5GfPnk1ZWZnbydBoNJoRhRDiiNk5Ha7RaDSaNEaLvEaj0aQxWuQ1Go0mjdEir9FoNGmMFnmNRqNJY7TIazQaTRqjRV6j0WjSGC3yGs/T3DXAS7sbXLMvpeSvW2sYDARdS4NGkyxa5DWe5+OPvc1df9hKr8/viv1/bK/jy0/v4JHXDrtiX6Oxgxb5EURH3yC7ajpSbvfiH63n7PteSbndMNWtvQAEXXq+TXuvD4CW7gF3EgCc+71Xmf21F3hsY6VraQgGJQG3boKHeHFXPRXN3W4nI260yCdAbXsf9R19rtm/7fEt3PDwxpTbPdraS3PXANWtvSz9v2uoPNaT8jR4gd9tOsJbh4+5Yruhsx+A/31+ryv2AT7wyFvM+/qL9A8GeHZbDW48VU5KyRObqpj9tRd4eU/qQ3hvV7Tw2T9u4/IHXufZbTV09Q+mPA2JokU+DnbWtPN0WTUX/GAd531/nWvp2FHd7pptgH9sr6V7wM9ft1a7loZDjV2uVqy/vONe3t0mXP6+/+I+vvTUDjaWp77B+8PbR/n2P/YAcOfvt/LEpqqU2b732V18dPXmofdfemoHF/1ofcrsJ4sW+Ti48eE3+a+/7nTFtpSSO377Duv3N7li30scaOjkqp9u4N8e2+J2Ulyl1+enuSu1oaMXd9UPvX79YDMATZ2pTUNFczff+vvu9xz769aalNn/05ajJxxr79WevMYBXt3fxO2/fcftZLiGEAKADz2yCUh9jyZs3yvc8NDGlI+RfPaP24ZeV7WExkj+8+kdQ+MlqeA/n96RMlvphBZ5jedxI/ZrZr+l20eTER93i8PN3hkTqWpJXVr6fIGU2UontMhrEsZlzXWVjeXHWPG9V91Oxqijzxdgf0PXCcdHc1mMFy3yHsdLhfjZd2vdToJn2F2b+qmsw9lT534aUkV1W+rCQumGFnkL1phM0+oZ8LvebU81FUaYINUhardj4tHsX/9Q6qeyDud9D7qfhlTxZR2PTxot8hZsNxnku+Ghja51292OUWs0qSZaqAZS73CMRLTIJ0nFKF0QBN4KIWk0mthokfc4Wk/d77m4bV9jjr411tgSeSHE/UKI/UKInUKIvwkhiiPO3SuEKBdCHBBCXGM/qZpIfP4gr+5rdDsZoxovDL5qNFbY9eTXAkullKcBB4F7AYQQS4CbgVOAa4FfCCEybdrSRPDAywe443dlbDrc4nZSlBNt4DOVK4DNBn69MPg6WnA79P7u0TaXU5A8tkReSvmylDK8/+tmoNR4vQr4s5RyQEpZCZQDK+zY0hxHSjhirDoM75A42vjlBr3t72jCLCqTqoHXho6RO5POyZj8p4B/Ga9nAJE7OdUYx05ACHGnEKJMCFHW3NzsYHKcIVrM7/mddSm0HzvoqEOSmtGMjslbYynyQohXhBC7o/ytivjMNwA/8MfwoSg/FfV2SClXSymXSymXl5SUJJOHlPO5J991OwmemDpW296HXz8taVQjXA+kaKzIsvqAlPLKWOeFEJ8ArgeukMfdzhpgZsTHSoHUub8O4gUx9SIt3QNc8IN1fOK8WXxn1VK3k6NxCan7kp7H7uyaa4GvAjdKKSPXHT8H3CyEyBVCzAEWAKN7f9gkiVaF3K5WQQntfaEtVt845M5DNDSaVDKSnT27MfmHgUJgrRBiuxDiUQAp5R7gKWAv8BJwj5RyRG4hp2N+J/Lo64dTvp/5cDZXtDLgH5FFSpMEbmvsSNYBy3BNLKSU82Ocuw+4z87va6LT6IE9c1K5j7gZP/zXAb59wxK3kzGq0TF572NL5EcDXuymnf+D448gHMkeRryY3YKjKWpovFgGvEKqYvKjoJgrQ29rYIHbIuq2fS9gdgle2dfI557cZnJWo3GOkdzQa5Ef4Yz22Q3P76y3/pBNdEPrPmYaO9rLfzxokddoNEmjY/LeR4u8BV7vprlVyVL5lHq3b4HXy4CbaE/a+2iRt8DtrrpVJXKrkt334j5X7Go0keiehDVa5DW2GO1+XM0of/ao2yLrhZ7EoMe39tAir/E87ldjcy784Xq3k6Bxmcc2VrqdhJhokbfA6/FYt8NJHr88GsV4wZNODeYlvc3j231rkbfAbRF1274X0A2Jxn1GbkXUIq+xxcgt+honSFVM3qxH7faYwEhAi7xGo/E8Zj3a0RMuSh4t8hqNJmlGj8iO3B6DFvkRjttVbOQW/fjxch6f31ln+YhIjVq8HjLSIq/xPF0DfusPKcTLEvq5J99Nyf49Znhd4FKB13szWuRHOG5XMW8X79GB16fwadxFi7zHseqJa5HVjIZojdnsmt21na6Hq7zem9Eir9FY4O0qrAm6vZbF466W3Qd5/68QYqfxfNeXhRDTjeNCCPGgEKLcOL/MmeRqhuO2F6Nxn9FQBkZBFpVh15O/X0p5mpTyDOB54NvG8ZXAAuPvTuARm3Y0HkV7uZrRTlqHa6SUnRFvCzgeIl4FPCFDbAaKhRDT7NgarXi9K+jt1DnDaMhjsvxzR11KHizv9h5Ssex7vY7ajskLIe4TQlQDt3Lck58BVEd8rMY4ptFo0oi/lFVzy682u50M5YzkcJGlyAshXhFC7I7ytwpASvkNKeVM4I/A58Jfi/JTUS+TEOJOIUSZEKKsubk52XxoNMrwdmfcfRo61HvyXsbr4Zosqw9IKa+M87eeBF4A/i8hz31mxLlSoM7k91cDqwGWL18+gttLjcYd3K40bnu5oYFnbwutm9idXbMg4u2NwH7j9XPAbcYsm3OBDimle8vyRjBuVyArdNXSjAbcHhOwg6Unb8EPhBCLgCBwBLjLOP4icB1QDvQCt9u0o/EoHm+DNBrleH0Kqy2Rl1J+yOS4BO6x89tewesj55oUMJLduBTgdh3ZeqSNc+ZOdM3+LzdU8PkrFjA2167PrAa94tXj6CbGA3jcU/N48pTz0dXuz+5p9/D+QVrkLfD6yLlG4zapaGR0PUweLfIjnNHuxWlGB26HhEYyWuQt0IVLo2PysdE1xNtokfc4Xh+517iPLiHqsWrmhYcdAS3yFuhYoEZjgY7Je9oZ0yKv0WhsMRpCmiM5h1rkLfB6AfZ6+jSa0YAO12iSRku4RuM+3pVwa7TIW+D1WKBG43Y8OBXmdY81ebTIW6ALl0YTG11DvI0WeY1G43l0jzp5tMh7HKuusNszt9wOFWg0mthokbdAexAaTWx0Q+9ttMhboGPysfHy1DFNatA1xNtokfc6Hg/XaDQab0+x1CKv0Vjg5QoM7jf0Xr8+TjCSe6xa5Ec4bpc9HY/VjIYSMJLLuRb5Ec4ILnsjBn2JNSMZR0ReCPFlIYQUQkwy3gshxINCiHIhxE4hxDIn7IxGvD7wO5K7sRpnSMmToVwuZiO5nNsWeSHETOAq4GjE4ZXAAuPvTuARu3Y00fF2E5AeeL16e90RcALdY00eJzz5nwJf4b16swp4QobYDBQLIaY5YEuj0Wg8h5cdfVsiL4S4EaiVUu4YdmoGUB3xvsY4Fu037hRClAkhypqbm+0kZ1TidtkayQNS6YK+BZpYZFl9QAjxCjA1yqlvAF8Hro72tSjHohZFKeVqYDXA8uXLdXEdhuW2BqlJhiYGUkpXY7Ze9iI17mMp8lLKK6MdF0KcCswBdhgFvBTYJoRYQchznxnx8VKgznZqNRrNCWhP3n28fA+SDtdIKXdJKSdLKWdLKWcTEvZlUsoG4DngNmOWzblAh5Sy3pkkayJxO1wykmcdaDSjAUtPPkleBK4DyoFe4HZFdjSaUY+HnUjH8PoMIi/7Oo6JvOHNh19L4B6nfns04+2i7X5PwgtI6e1KrrHPSL69esWrRqPRpDFa5DUaC7SXrhnJ/VUt8h7H6+GQ0TDw6vVprB4vIhqX0SI/wtH123283hCnA25f4pHsymiR19hCC5xG4+3HhGqRt4kWufTHKiKlS4DGy1M8tch7HO8WHY1GMxLQIq/R2MTtzpyXvUinSP8cqkOLvE3cruBul/7RMLtGo7HCdR2IgRZ5jS30mMTo8KTdRrsSyaNF3iaqq7fWUI1uSF3vsFri5fRpkddobOK2BrttPyW4HpZ0174dtMhrNBpNGqNF3iZud6V1PFgzKvC4J+22DsRCi7wVbnfF3U6AxhIP1+/0QV/jpNEir7GFrnsajbfRIm+FXtIeE+3F6t6Wxtv1QIu8RuNxLLc69rLCaFzHlsgLIf5bCFErhNhu/F0Xce5eIUS5EOKAEOIa+0l1CcsK5q59jXqsdhjUGqse3VtKHiee8fpTKeWPIw8IIZYANwOnANOBV4QQC6WUAQfsaSJwW2B05dNovI2qcM0q4M9SygEpZSVQDqxQZEvjIm43MqnAqiEbBZdAM4JxQuQ/J4TYKYR4XAgx3jg2A6iO+EyNcSztUO3JWv26FhjNaGhovfxQDq9jKfJCiFeEELuj/K0CHgHmAWcA9cAD4a9F+amoRVEIcacQokwIUdbc3JxkNtQxCuqPLUaDwFihBz7V43ZY0PLBMR4uApYxeSnllfH8kBDiV8DzxtsaYGbE6VKgzuT3VwOrAZYvX+7hS+VNtH/jPnrsXT1eFlGvY3d2zbSItx8AdhuvnwNuFkLkCiHmAAuALXZseRW3C58u++px+x5r3MdyGquHa6Ld2TU/EkKcQUhrqoDPAEgp9wghngL2An7gnpE6s8btrrgWGPexHBdRfI+syqAuIppY2BJ5KeXHY5y7D7jPzu9rvI/bjeBoQF9h99FbDWtcw22N1QKE8otgveJVrX2NNV6+B1rkRzgj2cPQaDTq0SJvgdtelNWAjtsa72UPxjEsY+Kq10qMhoscG30FkkeLvMYWuiehnlHRkFoQ9PhF8HLqtMhrNDZxW39Gg6fv+jUewZdYi7xN3K5gbnvSo8GRH8H1W+MQWuTTGLfv7UguXKMF5Ste9ewa17FeK+Hdm6BF3iZu31u9cZP7uF3B3e7NjQbcvsd20CKv0Xgc6xlWWuVVM3IlXou857EsXLp+K8d63xJ37WfoMqAct8uAHbTIW+D1m6vrt3q83lXX4ZpU4O0yEAst8ha4PXtG433UL4iLjdAqr5yRPPitRd4CL988cL+Cu21/NOD1nsRoYCTfAS3yFlhu86q4AuoK7n3cfgRkhm5olTOSq6EWeQuCI/jmapzB60XAbY13234qsG7IvVtKtMhbYLVnhtu3dhTUL9ex9OJc3mrY7TLgtv1UoD35NEZ78hrXsZxC6a7Mum0/FYxkGdAib0HQZZW39OLSv355HrcFwO0y4Lb9VGA9NpeihCSBFnkLAi7fXC3y7uP6M15db0Zio1fcehvbIi+E+LwQ4oAQYo8Q4kcRx+8VQpQb566xa8ctAm578hYV3O2ucsYocBO835tz25V313wq8LKnboWtB3kLIS4DVgGnSSkHhBCTjeNLgJuBU4DpwCtCiIVSyoDdBKcat0XeyrzrIu+2wKQAv8VNcH8KpVLzlrhtPxVY3WMvtwF2/bC7gR9IKQcApJRNxvFVwJ+llANSykqgHFhh05YrWIVrVN9dq9k9bous2/ZTQSAYdNW+VRlw+w6MhjIwkj15uyK/ELhICPG2EOJ1IcTZxvEZQHXE52qMYycghLhTCFEmhChrbm62mRznCQS83lVPTTq8aj8VWHryLjf0bodrRkERGNEibxmuEUK8AkyNcuobxvfHA+cCZwNPCSHmEv2+R71MUsrVwGqA5cuXe+5SWg68qu6qe3xUfzR4ca6H7Cw6Em7fArcbmVTg9uC7HSxFXkp5pdk5IcTdwLMypERbhBBBYBIhz31mxEdLgTqbaXUFtwfdvD5PfzTEY61j8mrxerjGbfupYCRvL2I3XPN34HIAIcRCIAc4BjwH3CyEyBVCzAEWAFts2nIFy5i8Yrz+lPqcrPSfXuN+Qx/bflamu/cgMzP9Zd560bN366mt2TXA48DjQojdgA/4hOHV7xFCPAXsBfzAPSNxZg2A3yIm73Y81u3CleOywKQC65i86pBd7POZLodLRkMZ8LCGW2JL5KWUPuDfTM7dB9xn5/e9gC/g7swKt+PBVmSPggo+aFEGVDf0VmXA7YY+Nzv9y4DbPXo7pP/dscnEghxX7Vt5kW6jRV49Xi8DuVmZbidBOT6/uw29HdK/htpkwZRCALZ+M/r4s+p7G/bivnz1Qi5ZWHKi/RQVrmnj8njiUycudSjMsxvxs+b00nEAbPvWVcptRWPupLEA/Oq25a7YD5eBxVMLo55PpcAsnVF0wrGCHPUif86cCQBcd2q0iX7qKcgNlfMPn1Xqin07aJG3IBiUZGYIJo7NZdd/X82ksan17MNjAmfNmsDjnzybm1woZJkZgg8tK+XihSU8eMuZ7zmXCpEfNyaHM2YWM6Egh833XqHc3nDC4ZDz502k8vvXpdy+35hD+cWrFlL1g/fx6YvmpDwNEwpy+Pi5s3j+8xfx17vOO+GcaiYX5TF3UgG/uPUsdv731crtDccfeO89MGtwvYgWeQv8hsgDFOZlU/bNq/jxh08fOq86Zh7+/axMQWaG4McfPp3Xvnzp0HnVs2+klAQirsGNp0/nwHevpXhMtnFeqXkgtOI0y7A/dVweVT94Hx8/dxYAeSmIB4fDJZkZAiEEVT94H0/++zkR6VN7EcINffgafON9S3jnG8d7lqkI5wwGgmQZs2iWz57AoftWctGCSQDkp8CT9weCQ2WwKC+biu9dx1evXazc7pD94HvvwUtfuJi/ffb8iPR5N16jRd6CSIEJc9NZpTz8sZBHe6x7QKn9sBeXGZGG2ZMK2P7tUOiiuUut/cCwwg2hGOz2b1/NKdOLaFacfwhVoMxh9+B/37+UT180BynVz24ZLrIA58+fxG9vDy3wVl8Gwg398epaUpjLwe+uRAj19iFUDiLzn52Zwe/vOIfz5k6kqTMFZSAo35P/jAzB3ZfO40tXLQSsY+a27RuefGQazjxpPM/cHRL6Yz3qr0GyaJG3wB+UUaeozSjOB6CuvU+p/WgiCzAuP5uCnExqFdsf8mKjzIUuHZ+vPP9gCEwU+1PH5TPgD9LRN6jUfqQnH8m0caEyUN/Rr9S+WRnIycpg0thcGhTbh1BDF20+/rRxecrzH7J/orMFMLkwF0C5s+GP6FFHtZ+Chi5ZtMhb4PMHoy74GT8mFIdULTDhmR3DZ7EIISgek5M6+1H2FC7OV28/nIasKPanFIUqWKPiChYWmOHL98MVvElxbyrsRQ5vZMJpUG1fSok/So8WoKQoNzW9ueCJvTkI9WgAmjrVNjTDwzUn2O9S39AlixZ5C/oHg+RlnxhzDI+29wz4ldrvGwytIYuehkzl9vsHg4b9E4tKQW4WPQPq17iF7kF0+wA9PvXXIFYZ6FV9D/yxykCW8jIwGJAEpYn9nCx8/uBQQ6SKAZMyMHFsSGRbe3xK7YfDQcMXfuVlZzI2N4vWHvXOTrJokbeg3x+IuthjrFHBuxWLXCyRHZsCke03GpncKBV8bG4mPT6/8ph432CA/Cj2xxjH+nyKr4E/EPX652RlkJ0p6B10rwwU5GQOOQLK7BuNTG6UHu0YY9BV+TXwm5QBw374GqmibzBATmZG1JBVfgrugR20yFswMBggL8pij7zsDDKEek8+LLLR0lCQm0W3YvsDRgWPVsEKcrOQEnpVi+xgIKoXOSbH8KRV2/cFTBf85GdnKm9kwr8fXeTUe/L9YftRZtEM3QPFzkafL3oZCF8T1SLb54vu7EFIC/q1yI9czEIFQggKctSL7HEvLnpXWXm4yGduf0yKQlZmIh8WnV7V4Rp/wHSaYEpENka4Jj8nBY1MDEdjTCrvQZT856VI5AdM7ENqGno7aJG3wExgIBTCGFA8dWsoXBKlq5ybnaF8b53jAhPFvpEm5dfAH4zqRR0XGPXhErP5+GNyM5WHKsICErU3l6PeftjRiO7Jp+Ye9PmC5EWxH07TQArugVlDn5+twzUjmn5/IKrAQkjkVM/P7fcHyMnKICPKzIKczAwGFMci+2MM/Iavi8qGJhiU+PzBmPHYVIRLogkshCp4v2L74UY0L+fEcpiXinDRoHlDHy4XqsMV/WZhU6MMpuIamHnyeVrkRzZd/X4K87KjnsvJUu9Jd/X7KcyNvnVAKjz57v5QN7wg58Q0hGcaqGzoun3m9sNxctU9iV6ffyg0dWIaMpTb7xnwk5Uhom7pG+5Nqhz8Ds8eys8+8RocF3l110BKGboHUTzprMzQ4Ldqke2N5cnnZOqY/EimvXdwaAn/cHIyM/D51d7cjt5Bxpnaz1Tek2g35sGPLzgxDTkpCNd09IbsR7sH4Z6E6grW3jfIeJN7kJetvoK394XKYLTH7IW9a5X3IFYZCNtXeQ26BvwEZfQyAKEwluqGtq3XR3F+dPu5Wep71HbQIh+DQFDS2T9oenNzUhCuae8zL1wpsW+IbHjxVyRhT1plGtp6fab2MwzvVnUFb+81LwN52ZlD4xaq6OgdZJypwIRj0gpFPtzQ5puXAaWNjDEH3fQaZGcMzQJTRVvPYNQyCKFroNq+HbTIx6CrfxApodj05qoPl8QqXOFwkcquenuvj9ysjKgx+bAnr1Lk22N48hCq4Cq9yKDR0JuLbIbyOdrtfT7TMnjck1d3Ddr7Qg1ttHuQCk/+uH1zkVV+D3pj3wPVjoYdtMjHwEpgclLQTevoMw/XpGJ2S8xw1dDAq7oKHvbkY1VwpeGivtgNfV62ei+urSdGTyJLfUy8rce8oR+KySu8Bm1WDb3icRGfP0iPL2AaslNdBu1iS+SFEH8RQmw3/qqEENsjzt0rhCgXQhwQQlxjP6mpJ7wnSHjp9HBUD7wGg5LmrgFKTOynYnZLU1c/k8zybwwEqmzowrtsmu3jn5edoXT6XF1HaAO2aePyTO2r9iLrO/qYamI/PLVUpcjWd/Sb5n/I0VB4DRqMezC1KHoaQs6WuvyHd/mcYFIGQ70574Zr7D7j9aPh10KIB4AO4/US4GbgFGA68IoQYuFIe5h3bXsvcHzHyeGEBl7VFe5j3QP4AkFmjDexn4JwSU1bH3MmFUQ9FxYYlY1MTVsfBTmZMcMlKr2o+vbQxlPTTcpAKFSgrlj3+QK09Q6a2s9LQUy+vqPf3H4KPPna9n6EIEZDp9aTDu/0Wjp+jIn9URCuEaFh/48AfzIOrQL+LKUckFJWAuXAic+O8zh1RgU3FXnFA581RuGK1ciAOpGXUlLb3mfeyGSqDxfVtvdROn5M1JkloH52S9iTNxV5xRU8bN+sDKRCZOva+4a2VR7O8RlOCq9Bex9TCvNMnyecl6V24LWmLeTslZrUg9ysTAJBqXyTtmRxKiZ/EdAopTxkvJ8BVEecrzGOjShq2nqZUJBjOj9WdSyuts2o4BaevKo0tPcO0usLmHswKepJmOUfwjFxtfcgJzPD9IHueVmhaaxBRU9nCpeBWOEiUDfwORgI0tjZz/Ti6PaFEKFwiUpPvq3P1D6EPHmVjcxQPTTtzaVm5XeyWIq8EOIVIcTuKH+rIj52C8e9eIBoblfUWiCEuFMIUSaEKGtubk4s9Yo52NjN/MljTc+rjskfbu5GCDhpgpnIqp3CeLi5G4DZE6PbV93IBIKSqmM9zDKxD+rjoQcbu5hbUhB1xTEc96RVXYODjV0AzDMph6qnUFY09xCUMK/EvB7kKZ6AcKipm7kx7KsO2VU09zC1KM90exPVZcAuljF5KeWVsc4LIbKADwJnRRyuAWZGvC8F6kx+fzWwGmD58uWeeVCilJKDDV28/0zzDkiu4gGf/fVdzJ5YMLTT33BUx+T3NYQEZvG0oqjnVTcyR1t76RsMcLKJfQhVsK5+dZtjHWjo4uw5E0zPH/fizFdE2rU/aWyO6eB3nuKB1/0NnQAsnmb+4GqVIbNj3QMc6x6I+eDsXMU9ib31nTHzH1kGvIgT4Zorgf1SypqIY88BNwshcoUQc4AFwBYHbKWM2vY+ugb8LIpRuFR78vsaOjk5RuFSPYVxf30nRXlZTDcJFShvZOpDAnPyVHORV+nJd/QNUtfRH7MMqF7Wf7Cxy1X7Bxq6yMoQzJ0Uw5NWOC5xIOxoxCwDmcp6Ej5/kMPN3TEdjaEZTh5d9eqEyN/Me0M1SCn3AE8Be4GXgHtG2syanTUdACyZHltgVAlcR+8gR1p6Ywqc6oHP3bUdLJleZDromZkhyMwQyhqZ3bUdZGUIFkyJESpQuOJ0l1EGTpk+LoZ9dTHx/sEA+xq6WBJLYBRv7bCzpoMFUwqjPgIzTJ7CGUbhehjL2VHZyBxs7GIwIC3uQThc402Jsy3yUspPSikfjXL8PinlPCnlIinlv+zaSTVbKlvJz85kaYwKnpOZQVCiZFR9S1UrACtihQoU7lvS1T/I7rpOVsw2tw9qp5G+XdnKaaXjTGOhoHbfkC2VLWQIOGvW+Bj21c1u2VnTgc8fZMWcieb2FcaDBwNBth5p45wYZRDUhmu2VLYwf/JY07UqYOxdo8j+5ooWAM6OUQ9SsVbADnrFqwlbKltZNqs4pgdzPFzi/M3dXNFCTlYGp88sNrevcApl2ZE2AkHJOXPNBQbUeVG9Pj87qts518K+UoGpauWU6eOGHvUY3b66Cr6lMiww5o2Myp7E7toO+gYDMR0NUDfwGQhKyqrarO0r9OTfrmxl1sQxpnP0IXW7oSaLFvkotHQPsK+hk3NieFCgtgXfXNHCsrcNDNoAACAASURBVJOKY3qxKncg3Hy4hexMwbKTzAUG1HnSZVVt+ONpZBQJTK/Pz7aj7ZYCo3I/9TfLW1g8tdB0SwUINfRCqHloxluHrb1YUNfQ7qrtoGvAb9mTCO8h5fQ0Vn8gyJbK1jh6Muk/8Jp2vLq/CSnh8sWTY35O1UKUuvY+9tR1cvHCkpifGwoVOFzBpJSs3dvIijkTLGeMqIqJv7KvkbzsDMtwkar91DccPIbPH+QKizIwFBN3uKFp7/WxpaqVK06ObV8Ioayhe3lvI6fPLKak0DxUAuq2dli7t4HMDMElcdYDp3vUW6pa6egb5LJFVmVA/apjO2iRj8LLexqZUZzPKTEGXUHdzIaX9zQAcO0pU2N+TlVMvrypm4pjPZb2Qc3slmBQ8vKeRi5ZWBJHI6PmGqzd20hRXlbM6ZMh+2oeP7dufxOBoOTqJdb3QIUn3dDRz47qdq5eMsXys6q22l27t5GzZ4+P2ZMJ2VcTsnp5TyO5WRlcssiikUnB/kF20CI/jF6fn43lzVy1ZIrprJIwqubHvrSngQWTx8ZcAALqBGaN0chcFafAOC2wO2raaejs59ql8TQyzntRg4Eg6/Y3cvniyaZL6cMcn6fudEPfyJSiXE6dYT7wH0bFdsdr9zUCxCfyCjz5ymM9HGzsjruRA2cb+nBv9qIFJabrVMLogdcRxr92NdA/GOS6U6dZflaFJ1/f0ceWylZWxiFweQrCNVJKnttRx7KTimMONkWmwWkP6p876snJzODyRdYCoyIe+tqBZtp6B7n+tOmWn1URMmvv9bFufxMrl04zXWkbiYrtjv+2rYYFk8fGXPGt2n6GgJWnxtebBGdFtuxIG7XtfXHVQz3wOsJ4ems1syaOiTmjIUyugpkNz26rJSjhQ2eVWn42O1OQIZwtXDtqOjjY2M1NZ820/jDOe3ED/gB/e7eGq06ZYrqP/nvsK9hP/a9bq5k0Nteymw5qQmbP7ajDFwjy4eXWZQDCDa1z9g83d7PtaDs3nVVq2ZsN23dSYINByTPbarlwQYnpxmiR5Cpo6J96p5qCnMzEGhkdrvE+1a29bK5o5aZl8RVup704KSVPlVVz7twJzJoYfXvfSEKDbs560k+VVZOXncENp1v3ZMD5Tdpe3ddEW+8gH1keXyPjtCff0j3Aq/ua+MCZ0y1DNSH7zofMni6rYcm0opiLsCLJzc5wNB781601ZGYIPhBjSw+V9jdVtFDb3seH43B0wHlPunvAzwu76rn+tOmWoRrw/t41WuQjeHprDULE50WD84N+b1e2cqSlN26BC6fBKS+uzxfgn9vruO7UaRTmWXvRYftOCtxTZdVMG5fHhfMnxfV5pz35v71biz8o4y8DDjf0++o72VXbwU1x2g+nwSlPejAQ5NltNVyysITJJg/piGZ/MCAJODSF8amyagrzsrgqjvEAcH7g9cWd9fT6Anzk7PjuQY6OyY8MBvwBnnz7CJcuLDHdO3w4Ts+R/v3mIxTmZbFyaXxedDgNTnmx/9heS9eAn48m1Mg415OoOtbD6web+fBZpWTGEYsO2XfOkw8EJU9sOsLyWeNj7pUSSXamQAjnGpknNlWRm5URtxcNznrSa/Y00Ng5wC0rTor7O07eg6bOfl7cVc+HlpXGXCMSiZNb/Uop+d2mKuZPHmu5RiRMZoYgO1Po2TVe5/kd9Rzr9vGpC+fE/Z3jXXX7hau2vY+XdjfwsRUnJbSboVMzK6SUPP5mJUumFVkuABpu36mezG/fqiIrQ/Bv585KwL5znvyr+xo52trL7RfEXwaEECFP2oEK3trj49lttXxw2QzGm+xfH408B/dTf3xjJbMmjrFcHxCJkw8O+cPmI/iDkk+ePzvu7xwPl9i/B1sqW9lT18ntF8yOK2QbRuUmaXbRIs9xgVsweWzcYQKIXAhjv3D97q0qAG5LoHCDc570G4eOcbCxmzsunJNQ4XbKfkffIE+VVXPDadPjDhOE7DvnRf7mzSqmj8vjmlPiCxNEpsEJgfvTlqMM+IN8KoFGBpzb8vrdo21sO9rO7efPjmtWTxinerT9gwH+8PZRrlg8hdkmj5yMRq6DW0s8trGS4jHZfPDM+MNloH67YztokQfeqWozWu/EBQ7sF66eAT9/2nKUa5dONX36jBm52ZmOzNF+bGMlJYW53HC69bTBSPKyMxyx/9Q71fT6Agn1pEL2nfHk99V3sqmihdvOn01WHAOukTgx+O3zB3liUxUXLZjEginmOy5Gw6m1Cr95s4rC3CxuSiBcF7YP9sMlz22vo7XHx6cumJ3Q95waeD3S0sPafY3cek5ivelQGrz7nFct8sBjGysoHpOdUBwUQk/EAfsezNNl1XT1+7kjQYELp8GuF1fe1MXrB5u57dxZMTdki0b4+ZaDNpaU+wNBfvtWFefMmcDSOBb/vNe+M/fg8Y2V5GVncPPZiQkcGIPPNiv4i7vqaewcSLiRC9u3m/+GjlAs/CNnz4y5IVs0nLgH4d704qmFnDcv9n5FKuzD8XDhbefNTvi7qh9DaYdRL/KHGrtYs6eRj587K+HWOyszg6wMewMu/kCQx96s5MyTiuMe6InECU/+Vxsqyc3K4NYEYuFhnJhh9MKuemrb+5Jr5BzwImvaevnbu7XcfPZJlkvozdJgR2CCQckjrx1m/uSxXLLAem7+cJzoSazeUIGEhGLhYZwI16w/0MT+hq6Ew4XgzMBrW4+Pp96p5vrTpjMlgXBhmBzFT4mzw6gX+V+8dpj87MyEBtsisTvw+dyOOqpb+/jspfOT+r5dT762vY9nttVwy4qTmJDAYN+QfZsVPBiU/Hx9OQsmj+XKkxOLhYMzXtzqDRUIAXdePDep7+dm2QtZvbKvkQONXdxz2byEYuFh7PYkWroHeHLLEVadMZ2ZJs8TjoXdBWFSSh5aV86M4vyYj9s0t29/4PU3b1XR4wtw1yXzkvq+U2FTFYxqkT/a0stzO+q49ZzkBA7sTWEMGAK3eGphQrMZTrSffOH65euHbQscJC+yL+9t5GBjN5+7fH6SAmfPk2/q6ufP71TzwTNL4546O5xcG568lKEycNKEMdwQxzYK0cjLysQflEk/vOaxjZUM+INJOxp2FwVuOtzCu0fbuevSeXEtQBtOns156p39g/z2zUquOWVKzEctxsLu4PeDrx5i65G2pL8fi1Et8o+8fphMIfh0kgIH9qavvbS7gcPNPUkLXMh+8vHYps6QwN10VvICZ0dkpZQ8vP4QsyeO4X1x7BUUjfCDU5K9Bo+9UYk/EOTuS5Pz4MBeQ7ux/Bg7ajq4+9J5CQ/4hsm1sUlaR+8gT2w6wnVLp8W1T0007IbsHl5fzuTC3LhXuA4nKzODzAyRtP3fbzpCZ7+fz122IKnvg72B1711nfxk7UE2HjqWtP1YjFqRb+jo55mtNXx4eWlSMbgwudkZ9CUhMKEu6iHmlhQktPhpOHbiwas3VBAISu6+JDkPDo57cX2+xNPw2sFmdtd28tnL5ictcBkZIhQPTaKCtff6+MPmI1x/2vSEpuwNx07I7KF15UwtyuODyxIPUwzZt7G1wu82VdE94Oeey5IvA3ZCdluPtPHW4RbuvHhu3IufopHslte9Pj+Pbazk0kUlnFqa2KD/e+0n39A/vP4QhblZfDLBWUXxYkvkhRBnCCE2CyG2CyHKhBArjONCCPGgEKJcCLFTCLHMmeQ6x+oNFQSkTDoGF2ZsbhY9A/6Ev/fqvtBA0z2Xzo97dWc0CnKz6B7wJ/zQjJbuAf749lFWnTGdkyYmHocNU5gXmomR6DWQUvLQq4eYUZyf8Kym4SRbwX/zZigOa0fgIPlwzTtVrWypbOUzl8wdaiyTYWhrhQRFpmfAz+NvVnLF4skxH1hvhZ1dIH++vpzxY7L52Dnxr7A1S0MyIvvk20dp7fHx+cvtlYHQuEjiZeBAQxcv7mrgkxfMZlx+fFuJJIpdT/5HwHeklGcA3zbeA6wEFhh/dwKP2LTjKHYHmiIZm5tFd38SAre+nJkT8ll1RnJx2Ej7gwGZcAF//M1K+v2BpOOwkfYBuhK8BpsOt7DNRhw2kmTCJd0Dfn77VhVXLUk+DjtkP0mBeXhdORMLcrj5bJsCl+RuqH98+wjtvYPcY1vgkhv43F3bwbr9Tdxx4Zy4NgKLRTIzjPoHA6zeUMF5cydy1qz4V3mb2U+mkXt4fTkFOZkJL4BLBLsiL4GwCzAOqDNerwKekCE2A8VCiORjEg7z6OuH8dkYaIqkMC8rYYF749AxdlS389lLkw9ThCnKS1xkO3oHeeKtI1x3avJx2DBhT75rYDCh7z20zl4cNpJkNkn7zcZKOvoG+ZxNLz5kP3GB2V7dzusHm7njojkJT90dTjIPTun1+Vm9oZIL5k9MaupuJMkuSHtoXShM8fEk5qUPpyA3k94EQ4ZPl1XT1DVg24sHGJubSXeCvdnypm6e31nHbefPTmgbi0SxK/JfAO4XQlQDPwbuNY7PAKojPldjHDsBIcSdRqinrLm52WZyrGnq7OeJTUd4/5kzbAscQGFeNl398QuclJKfrD3I9HH24rCR9oGE0vDrjRV0DfgdEbiw/UR6M5sOt7CpooXPXDLPVhw2TKLx0I6+QX71RgVXnjyZ02cWO2A/8Wm0D7x8gAkFOXzCAYE7/nSq+EXuiU1HONY9wBevXGjbfjIzrHbVdLBmTyN3XDTHkTBFUX42nQnUgf7BAA+vL+esWeMTXnxlZr+rfzChh4n/7JWD5Gdn8u9JrA9JBEuRF0K8IoTYHeVvFXA38EUp5Uzgi8Bj4a9F+amouZdSrpZSLpdSLi8pSXwhSKL84rXD+IOS/7gi+ZH0SBL15Nftb2J7dTv/54oFtuKwkfYhfk++pXuAxzdWcv1p0zh5WvJx2OH2O+O0H2rkDjC1KI9bbcZhw+RnZ9Lri/8ePPZGBZ39fr54lX2Bg+PTaOMdF9lS2cobh45x9yXzKEhwdWk0Ep3C2D3g55evH+aShSUst3hQejxkJzG75SdrDzAuPzupFb7RKMrLprMvfpH/49tHaewc4D+vXpjw4isz+0EJPXGWw331nTy/s55PXTCHiWNjPyjdLpYlTEp5pdk5IcQTwH8Yb58Gfm28rgEi14eXcjyU4xp17X08+fZRPnxWaVwP5YiHwrxsun1+gkFpOQ0yGJQ88PJBZk0cE/d+5VYkGhP/5YYK+gYDfMEBDw5CApeTmRG3/Q2HjvFOVRvfff9SR7x4gOIx2XTEWcHbenw8/mYV1506Ne6HcliRn5NJUIamEFrlSUrJj18+wOTC3IR224zFGCPc0zsQn8j/ZmMlbb2DfMmhRg5gTHZm3AK39Ugb6w8085VrF1EU53MLrCjKz+Zoa29cn+31+XnktXLOnzeR8+fFvyFhbPvHnZ14nsXwk7UHKczL4tMXJT99O17shmvqgEuM15cDh4zXzwG3GbNszgU6pJT1Nm3Z5uH15Ugkn3MgBhemKC8LKaE7jgK+Zk8De+s7+Y8rFtgebAyTSLimqbOf371VxQfOLHUkVHU8DVl0xxGTl1LywMsHKB2fn9CDUawYl59Ne5wi/8sNFfT4/I41chBqZADae63T8GZ5C1sqW7nnsvm2Y/FhxhtbMcRzDTp6B1n9RgVXLZniSKgqTHFBdlz5h5AXP9GhUFWYorysuD353711hGPdPv7zaufKQLixiicNO6rbWbu3kTsvmhvXIy7tYrev+Gng/wkhsoB+QjNpAF4ErgPKgV7gdpt2bFPd2stT71Rzy4qTKB1vb0ZNJJGedCyvJBAMxeLnlRSw6gz7sfgwxwc+rRuZn68vJ+BgqCrM2DhDVmv3NrKzpoP7bzot4Y3QYlE8Jj6Bae4a4HdvVXHj6dNZmOBOj7EIi2xbry/mw8+llDyw9gDTx+Vx8wrnGrnignAj47P87K83VtDV73fUi4fQNWiLw/6mwy28Wd7CN993siOhqjDhmLyUMmb4pbN/kEdfP8xli0psz6gZbh/iE/kH1h5k/Jhsblcciw9j6ypLKTcCZ0U5LoF77Py20zz46iEyMoTtOdHDea8nbb5q9J876jjU1M3PP7bM1rz44RQN2Y8tsjVtvTy55SgfOXumrXnx0YhnXCJoNHJzJxXYnhc/nOL8HNp7fZYhs0dfP8yAP+B4Izck8j2xRW79gSbePdrO9z94qiPjMWEKc7PIyhCWItva4+PxjZW8z6HxmEiKx+TQZtHQhsdjphQ5F6oKU5SXzWBA0j8YjNlDetyYVfWlqxY5aj/esal3qlrZcLCZe1cuTni3z2QZFSteK5q7efbdWv7tnFkxPa1kCHfV23rMC/hgIMjPXjnIydOKWLnU+unviTA2L4sMYe3FPbyuHCGEI9PFhhOPF/f8rnr2N3TxhasW2p42OpziMaFBr1ghs4aOfv6w+QgfXFbK3BLnQlUA4w1POpbIhRu5kyaMSej5rfEghKB4TDatMcoghPYp6hsM8MUrnW3kAMaPybZs5N4wxmPuuWy+Y+MxYcIx8VhjM209Ph57o5JrT5lqa3VrVPuGsxXLvpSSH685wKSxuUltZ5wso0Lk719zgLysDFv7k5gxuTA0Mt7cPWD6mafLaqhq6eVLVy1Meo8aMzIzBBPH5tLUaW6/ormbp7fWcOs5JzFtXHJ71MSixML+YCDIz9YeZNGUQq5Pco+aWIS3B+6IIbI/XXsQKXHciweYEBGuMeP5XfXsrnV2PCaS4jE5MRv6uvY+fvtWFe8/YwbzJzsXqgpj1dAHg5IfvrSfGcX5fDSJPfutKDFmqDR3mZfDR18/TLfPuVlVkUwqtLb/xqFjvF3Zyj2XzXNsPCYe0l7ktx5p41+7G/jMJfMoKXR+qtLkwlDPoKmzP+r5ngE/P1l7kOWzxnPlycntNGmdhlyauqLbB/jhS/vJz850PFQVpqQol+auAdMphH/ecpSKYz381zWLHG/kACYYnrRZQ3ugoYunt1Zz23mzbK9wjkaxRbhmwB/g/jX7OXlaUVJb6cbDhDE5tMbwpB94+SAS+JKDg42RjB+TQ1e/3/ThMc/tqGNPXSf/dc0iR0NVYcL7TzWa1MOatl5+81YVHzyz1PYK52iMzc1ibG6Wqf1AUPK9F/cxc0K+7S0cEiWtRV5Kyfdf3EdJYS7/fpGaQY6i/CxysjJoMmnBV2+o4Fj3APded7Ij83GjMaUoz9T+O1WtrNnTyF2XzGWSovm4kwvz8AWCUbuqXf2D/OyVQ5w7dwJXKGrkwjto1rdHr2A/fGk/BblZyhq5nKwMJhTkUG9Swf+w+SjVrX3cu3Kxo+MxkUwZl0eDif29dZ08+24Nt58/29FJB5FMHRcqWw0dJ6ahfzDA/WsOsHRGETcm+HjJeBkSeRNn58drDiCAL1+jppEDmFxk7mw9u62G/Q1dfOWaxUoauViktci/vLeRsiNtfPHKhbb3xjBDCBHypKNUsKbOflZvqOC6U6dy1ix7S8djEfLkTxR5KSX3vbCPqUV53HGhuvm44ZBVtDQ8+vphWnp8fOO6JcoaubDI17X3nXDurcPHWLe/iXsum6906fiM4nxq206039E3yEPrDnHRgklcvFDdYr8ZxfnUtfdFXXH5g5f2U5SX7cg2Hub2Q41HbZR78MSmKmrb+/j6ypOV9OQAJo3NQQhojBI23FXTwd+313HHhXOUhCvDTCnMi2q/zxfggZcPcvrMYq4/LfW7u6StyPv8QX740n7mlRTwkeXODnQNZ3pxPtVRKvhPXzmIPxjkK9csVm6/uWvghFWfL+yqZ3t1O1+6eqHSGGBYZI+2vHcxSn1HH79+o5JVZ0x3fKArkqK8bApzs04QmEBQ8v0X9zN9XF5Sj7VLhNLx+dS0nbgY5xfry+noG+Sr16otA6Xj8xkMyBMa2g0Hm9lwsJnPXz5f6Zzs0vGhMlAzrB609fh4eF05ly4q4fz5ziw8ikZWZgaTC3NPaGilDIVJJhTkcJeCMblIpo3Li+poPP5mJQ2d/Xx95WJljk4s0lbkH9tYSUVzD9983xLHZ3MMZ/7ksZQ3db8nJr2jup0/v1PNx8+dbWuv8nhYYCxsqmjuGTrWM+Dnvhf2cfK0Ij60TG0jF15YVd7c/Z7j331+HwBfvtrZ6WrRmDE+nyMtPe859qctR9lV28FXVy52fDbHcEIi30cgwpM+1NjFYxsr+dCy0oQfUJ6MfeA9qz4H/AH++7k9zJo4ho+f5+yUxeFMK85DCE5YdfrDl/bT4wvw9etOVmofwvWw6z3HXthVz6aKFr5w5QLHVtea2p8ylvqO/vcsTKxt7+PhdeVcvWQK58y1v0dOMqSlyNe19/Hgq4e4eskULkvysXqJsGDyWDr6BodG1v2BIF//2y5Kxubyhaucn81xgv0pIZE92Hi8gP+/Vw9R39HPd99/irI4cJhx+dlMLszlUONxkV9/oIkXdtXz+cvnKxnsHM6S6UXsruscamhbuge4f80Bzp07QVkcOJJFU4sY8Ac5bDR0Ukq+9Y/dFORmce9KtV48MDTvfXdtx9CxX22ooOJYD9+58RTlceDcrEzmlYxlT4T9bUfb+PM71dxx4RxHF5+ZsXBKIYeauodCVl39g/zPP/eydEYRt56jtpEDWGTk8WBEPfjOc3sA+PYNS5TbNyMtRH7rkVZue3wLPQOhPWS++sxOJJJvXZ+aC3uq4aWVGc9o/PXGSvbUdfLtG5Yo9x4AZk0soCAnc8j+9up2HttYyS0rZjq6qi8Wp84Yx9YjrUAoDv2tv+9mbkmBrUcrJmq/uWuAuo5+pJR88++76Rnw8z+rlqaki3yGsUXAu0dD9+BPW6rZXNHKf12zSPkGVBAaeJxalMe71e1AaEbRQ+vKWbl0KpcuUu/oAJxeWsz26naklPT5Anz1rzuZWpSnZNpqNJZMK6LXF2B/Q8jZ+d6L+2juHuC+95+q3NEBhh68ss2oh8/vrOPlvY38nysWKBvwjofULLlKARsONvPcjjrqO/p549AxvveBU1PiQUKoghflZfGv3Q0U52fz4zUHWLl0atLPLU2U7MwMLlwwiVf3NVJ9yTzu+eM2phbl8bVr1XeRw1y6qIRX9zexs6adB189RENHP3/5zHkpm0kQHtR8fkcdmRmCf+1u4GsrF6fEgwSYO6mAaePyeGFXAwunFPKdf+7hogWT+NiK1E2XO3/+RNbuaaS2vY97ntxGYV42/7NqacrsXzB/Is9sq2FTRQtPl9VwqKmb39+xwtHtC2Jx6aLJCAEv7qpnd20Hf9pSzV2XzHN0j55YTBuXz+Kphbywq56LFk7i3md2ceZJxcpm9sWLSPSxcSpZvny5LCsrS/h7UkpufPhNdhldxZvOKuX+m05L6SDH91/cxy83VACwcMpYnr7rfGWP84rGhoPN3Pb4FgAKcjJ58tPnpqxwQ2iA7eL71w9tb/A/q05J6ao+gFtWb2ZTRQsA15wyhUduPUvZbI5o/Hx9OfevOQDASRPG8Mzd5ytZm2HGtqNtfOiRt5AyNK3zd7evcGSv9Hjp9fm55P7XhsKWX756IZ+7PDVefJjP/L6MNXsaAbhw/iR+c/vZShafmfH7TVV86x+hEM2Uolyeufv8lHjxQoitUsrlUc+lg8hDaLHDo68fZn7JWD5+3uyUdM8i6R8M8Iv15QwEgtx18Tyl0/XM+OvWGrYeaeMT589i8VRn9yaJh21H2/jLlmouWzyZax3eviEemjr7eXh9OTOK87n9gjmOboIWD/5AkF+8dpiOvkE+c8ncoYVyqeSfO+p4s/wYN684aSiElEp213bw+01HOG/eRFadMT3ls0naenw8vL6c4vxsPm3z4eDJEAxKHttYSW17H5++eC4zitVN2YxkVIi8RqPRjFZiiXxaDLxqNBqNJjpa5DUajSaN0SKv0Wg0aYwWeY1Go0ljtMhrNBpNGqNFXqPRaNIYLfIajUaTxmiR12g0mjTGU4uhhBDNwJEkvz4JOOZgckYCOs+jA53n0YGdPM+SUkZ9Ko2nRN4OQogysxVf6YrO8+hA53l0oCrPOlyj0Wg0aYwWeY1Go0lj0knkV7udABfQeR4d6DyPDpTkOW1i8hqNRqM5kXTy5DUajUYzDC3yGo1Gk8akhcgLIa4VQhwQQpQLIb7mdnrsIIR4XAjRJITYHXFsghBirRDikPF/vHFcCCEeNPK9UwixLOI7nzA+f0gI8Qk38hIPQoiZQoj1Qoh9Qog9Qoj/MI6nc57zhBBbhBA7jDx/xzg+RwjxtpH+vwghcozjucb7cuP87Ijfutc4fkAIcY07OYofIUSmEOJdIcTzxvu0zrMQokoIsUsIsV0IUWYcS23ZllKO6D8gEzgMzAVygB3AErfTZSM/FwPLgN0Rx34EfM14/TXgh8br64B/AQI4F3jbOD4BqDD+jzdej3c7byb5nQYsM14XAgeBJWmeZwGMNV5nA28beXkKuNk4/ihwt/H6s8Cjxuubgb8Yr5cY5T0XmGPUg0y382eR9y8BTwLPG+/TOs9AFTBp2LGUlm3XL4IDF/E8YE3E+3uBe91Ol808zR4m8geAacbracAB4/UvgVuGfw64BfhlxPH3fM7Lf8A/gKtGS56BMcA24BxCqx2zjOND5RpYA5xnvM4yPieGl/XIz3nxDygFXgUuB5438pDueY4m8ikt2+kQrpkBVEe8rzGOpRNTpJT1AMb/ycZxs7yPyGtidMnPJOTZpnWejbDFdqAJWEvII22XUvqNj0SmfyhvxvkOYCIjLM/Az4CvAEHj/UTSP88SeFkIsVUIcadxLKVlOyvJhHuJaI+DHy3zQs3yPuKuiRBiLPAM8AUpZacQ0bIQ+miUYyMuz1LKAHCGEKIY+BtwcrSPGf9HfJ6FENcDTVLKrUKIS8OHo3w0bfJscIGUsk4IMRlYK4TYH+OzSvKcDp58DTAz4n0pUOdSWlTRKISYBmD8bzKOm+V9RF0TIUQ2IYH/IittSwAAAZ1JREFUo5TyWeNwWuc5jJSyHXiNUAy2WAgRdrwi0z+UN+P8OKCVkZXnC4AbhRBVwJ8JhWx+RnrnGSllnfG/iVBjvoIUl+10EPl3gAXGKH0OoUGa51xOk9M8B4RH1D9BKG4dPn6bMSp/LtBhdP/WAFcLIcYbI/dXG8c8hwi57I8B+6SUP4k4lc55LjE8eIQQ+cCVwD5gPXCT8bHheQ5fi5uAdTIUnH0OuNmYiTIHWABsSU0uEkNKea+UslRKOZtQHV0npbyVNM6zEKJACFEYfk2oTO4m1WXb7YEJhwY3riM0K+Mw8A2302MzL38C6oFBQi34HYRika8Ch4z/E4zPCuDnRr53AcsjfudTQLnxd7vb+YqR3wsJdT13AtuNv+vSPM+nAe8aed4NfNs4PpeQYJUDTwO5xvE84325cX5uxG99w7gWB4CVbuctzvxfyvHZNWmbZyNvO4y/PWFtSnXZ1tsaaDQaTRqTDuEajUaj0ZigRV6j0WjSGC3yGo1Gk8ZokddoNJo0Rou8RqPRpDFa5DUajSaN0SKv0Wg0acz/B3+xXcuj1tIoAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(s[4999:, 0])"
   ]
  }
 ],
 "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.7.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}