{ "cells": [ { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAEICAYAAAB25L6yAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABE6klEQVR4nO3dd3hUVf7H8fdJIYH0TkghCb0kgSTSpKMUwYqu+sOKir33RUXX1cW1rAV3FVFRVOwoKnZBmpQEAoRQDCG9907a+f1xhxCQMoEkM0m+r+eZZzJz79z53mTyycm5556rtNYIIYSwXjaWLkAIIcTJSVALIYSVk6AWQggrJ0EthBBWToJaCCGsnAS1EEJYOQlq0WaUUk8qpT5opW1dp5Raf4bbGKeU2neS5SFKKa2UsjuT9xGitUlQC5RSKUqpWqWU9zHPbzcFV4gZ25iolMposyJbgdZ6ndZ6wOHHpv0+53S3p5Raavr+jGj2XF+lVItPTlBK/U0ptVEpVaWUWnOc5cOUUnGm5XFKqWHNliml1HNKqULT7TmllDrd/RLWR4JaHHYQuPLwA6VUONDDcuV0GEXAP1tpOy8DC49doJTqBnwNfAB4AO8BX5ueB5gHXAREAhHA+cDNrVCTsBIS1OKwZcA1zR5fC7zffAWllINS6gWlVJpSKlcp9YZSqrtSygn4HuillKow3XqZXtZNKfW+UqpcKbVbKRXTbHuDlFJrlFIlpmUXNFvmpZRaqZQqU0ptAfqcqHCl1HtKqftNXweYWrm3mx73UUoVKaVsmrf6lVLLgGDgG1O9DzXb5BzTPhYopeaf4vv2HhChlJpwivVOSmv9i9b6UyDrOIsnAnbAy1rrQ1rrVwEFTDYtvxZ4UWudobXOBF4ErjuTeoR1kaAWh20CXE3haQtcgdGCa24h0B8YBvQFAoAntNaVwAwgS2vtbLodDpwLgI8Bd2AlsAhAKWUPfAP8BPgCdwIfKqUOd028DtQA/sBc0+1EfscIM4AJQDIwvtnjdVrrxuYv0FpfDaQB55vq/XezxWOBAcAU4Aml1KCTvHcV8CzwzPEWKqX+a/pDdLzbzpNst7khwE599HwPO03PH16+o9myHc2WiU5Aglo0d7hVfS6wB8g8vMDU5zkPuFdrXaS1LscIqCtOsc31WutVWusG0/YjTc+PApyBhVrrWq31b8C3wJWmPxSzMf0R0FonYLRcT+R3YKxSygYjoP8NnG1aNsG0vCWe0lpXa613YIRe5CnWfxMIVkrNOHaB1vo2rbX7CW4RZtbjDJQe81wp4HKC5aWAs/RTdx4S1KK5ZcD/Yfzb/P4xy3ww+qzjDrcIgR9Mz59MTrOvqwBH06iKXkD6MS3dVIxWug/Gv/rpxyw7Lq31AaASo6U/DiPws0yt89MJ6mNrdj7ZylrrQ8DTpltbqABcj3nOFSg/wXJXoOKYFrjowCSoRROtdSrGQcXzgC+PWVwAVANDmrUI3bTWh0OspaGQBQSZWsGHBWO04vOBeiDomGUn8ztwKdDN1E/7O0bfrQcQf4LXtGaQvYvRvXNJ8ydN/fgVJ7jtNnPbuzH6wZu3kCNMzx9e3rzVH9lsmegEJKjFsW4AJpv6nZuYWr5vAf9RSvlC04G7aaZVcgEvpZSbme+zGaO1+pBSyl4pNRFjtMLHpm6SL4EnlVI9lFKDMUL3ZH4H7gDWmh6vMT1eb9re8eQCYWbWe1Ja63pgAfDwMc/f0qzf/thbUz+yUspWKeWI8Z+EjVLK0dSPf3hfGoC7TAd07zA9/5vp/n3gPtPPoxdwP7C0NfZLWAcJanEUrfUBrXXsCRY/DCQBm5RSZcAvGAfd0FrvBZYDyaaukV4n2Mbh96nFCOYZGK31/wLXmLYDRsg6Y3RDLMVosZ7M7xh9toeDej1GV83aE74C/gU8Zqr3gVNs3xzLgezTfO3VGP+x/A+j+6Ya4w/j4e/VRRjHD0owDqxeZHoejD7yb4BdQALwnek50Uko6cYSQgjrJi1qIYSwchLUQghh5SSohRDCyklQCyGElWuT6Ry9vb11SEhIW2xaCCE6pbi4uAKt9XFPIGuToA4JCSE29kQjvIQQQhxLKXXCs2+l60MIIaycBLUQQlg5CWohhLBy7XZtuLq6OjIyMqipqWmvtxSiXTg6OhIYGIi9vf2pVxbiNJwyqE1TRX7S7KkwjHmCX27JG2VkZODi4kJISAgyTa7oLLTWFBYWkpGRQWhoqKXLEZ3UKYNaa70PY55fTBO6ZwIrWvpGNTU1EtKi01FK4eXlRX5+vqVLEZ1YS/uopwAHTPMWt5iEtOiM5HMt2lpLg/oKjKkc/0IpNU8pFauUipXWhRCiq9maUsSSdcm0xYykZge16dL0FwCfHW+51nqx1jpGax3j43OqqzN1PGvWrEEpxTfffNP03KxZs1izZk271vHGG2/w/vvHXiXLPM7Ox7+ilFKKq666qulxfX09Pj4+zJo166Tbi4+PZ9WqVU2Pn3zySV544YXTqq01Xn8iL7/8MlVVVU2PT/R9MNeZvl50PjmlNdz6wTY+3JxGVe2JrlNx+lrSop4BbNNa57Z6FR1EYGAgzzxz3ItNt4v6+npuueUWrrnmmlbdrpOTEwkJCVRXVwPw888/ExAQcMrXHRvU1urYoBaiNR2qb+CWD+Korq1n8dXRODm0/mC6lgT1lZyg26MjeOKJJ3j55ZebHs+fP59XXnmlRduIjIzEzc2Nn3/++S/LQkJCKCgoACA2NpaJEycCRivx2muvZdy4cfTu3Zsvv/yShx56iPDwcKZPn05dXR0AcXFxTJgwgejoaKZNm0Z2tnGhkIkTJ3LPPfcQExPDK6+8clSrMykpiXPOOYfIyEiioqI4cOAAFRUVTJkyhaioKMLDw/n666/N2rfzzjuP7777DoDly5dz5ZVXNi2rrKxk7ty5jBgxguHDh/P1119TW1vLE088wSeffMKwYcP45BNjYFBiYiITJ04kLCyMV199tWkbL730EkOHDmXo0KFH/RyeeeYZ+vfvz9ixY9m3b99xa7vuuuu49dZbGTVqFGFhYaxZs4a5c+cyaNAgrrvuuqb1fvrpJ0aPHk1UVBSXXXYZFRUVvPrqq2RlZTFp0iQmTZrUtO78+fOJjIxk1KhR5OYabY+UlBQmT55MREQEU6ZMIS0tDYCDBw8yevRowsPDeeyxx8z6foqu48mVu4lPL+HFv0XSz8/l1C84HVrrU94AJ6AQcDNn/ejoaH2sxMTEpq+fXJmg//bGxla9Pbky4S/v2dzBgwf18OHDtdZaNzQ06LCwMF1QUHDS1zS3evVqPXPmTP3777/r8ePHa621njlzpl69erXWWuvevXvr/Px8rbXWW7du1RMmTNBaa71gwQJ99tln69raWh0fH6+7d++uV61apbXW+qKLLtIrVqzQtbW1evTo0TovL09rrfXHH3+sr7/+eq211hMmTNC33nprUx0LFizQzz//vNZa6xEjRugvv/xSa611dXW1rqys1HV1dbq0tFRrrXV+fr7u06ePbmxs1Fpr7eTkdNx9c3Jy0jt27NCzZ8/W1dXVOjIysml/tdb60Ucf1cuWLdNaa11cXKz79eunKyoq9Lvvvqtvv/32o2obPXq0rqmp0fn5+drT01PX1tbq2NhYPXToUF1RUaHLy8v14MGD9bZt25qer6ys1KWlpbpPnz5N+9bctddeqy+//HLd2Niov/rqK+3i4qJ37typGxoadFRUlN6+fbvOz8/X48aN0xUVFVprrRcuXKifeuqpv/xstNYa0CtXrtRaa/3ggw/qp59+Wmut9axZs/TSpUu11lq//fbb+sILL9Raa33++efr9957T2ut9aJFi477fWz++RZdxwebUnTvh7/V//5hzxlvC4jVJ8hUs9ro2rjQqVfb/KloHyEhIXh5ebF9+3Zyc3MZPnw4Xl4t36Xx48cDsH79erNfM2PGDOzt7QkPD6ehoYHp06cDEB4eTkpKCvv27SMhIYFzzz0XgIaGBvz9/Ztef/nll/9lm+Xl5WRmZnLxxRcDxkkXYJxY9Pe//521a9diY2NDZmYmubm59OzZ86Q1RkREkJKSwvLlyznvvPOOWvbTTz+xcuXKppZ8TU1NU2vzWDNnzsTBwQEHBwd8fX3Jzc1l/fr1XHzxxTg5OQFwySWXsG7dOhobG7n44ovp0aMHABdccMEJ6zv//PNRShEeHo6fnx/h4eEADBkyhJSUFDIyMkhMTOTss88GoLa2ltGjRx93W926dWvqf4+Ojm76D+mPP/7gyy+Ni69fffXVPPTQQwBs2LCBL774oun5hx9++DhbFV3NloNFLPh6NxP6+3DfuQPa9L3a7czE5hacP+TUK7WBG2+8kaVLl5KTk8PcuXP/svz111/nrbfeAmDVqlX06nX867POnz+ff/7zn9jZHfn22dnZ0djYCPCXsy8dHBwAsLGxwd7evmk4l42NDfX19WitGTJkCH/88cdx3+9wwJnjww8/JD8/n7i4OOzt7QkJCTH7bNALLriABx54gDVr1lBYWNj0vNaaL774ggEDjv4wbt68+S/bOLyvALa2ttTX15td+8k0/x42f4/D30NbW1vOPfdcli8/de9c85+BuTXKEDzRXGZJNbd+EEeQZw9evXI4tjZt+/noUnN9XHzxxfzwww9s3bqVadOm/WX57bffTnx8PPHx8ScMaYCpU6dSXFzMzp07m54LCQkhLi4OoKn1Za4BAwaQn5/fFNR1dXXs3r37pK9xcXEhMDCQr776CoBDhw5RVVVFaWkpvr6+2Nvbs3r1alJTzR/yPnfuXBYsWNDUWj1s2rRpvPbaa03DjrZv395UQ3l5+Sm3O27cOL766iuqqqqorKxkxYoVjBs3jvHjx/PVV19RXV1NeXn5USNqWmrUqFFs2LCBpKQkwOhX379/f4vqHDNmDB9//DFg/MEbN24cAGefffZRz4uurbq2gZuXxVJb38hb18Tg1r3tpw7oUkHdrVs3Jk2axN/+9jdsbW3PaFvz588nPT296fGCBQu4++67iYmJafG2u3Xrxueff87DDz9MZGQkw4YNY+PGjad83bJly3j11VeJiIhgzJgx5OTkMGfOHGJjYwkPD+f9999n4MCBZtcRGBjIXXfd9ZfnH3/8cerq6oiIiGDIkCE8/vjjAEyaNInExMSjDiYeT1RUFNdddx0jRoxg5MiR3HjjjQwfPpyoqCguv/xyIiMjmTFjBmeddZbZtR7Lx8eHpUuXcuWVVxIREcHo0aPZu3cvAPPmzWP69OlHHUw8ntdee413332XiIgIli1b1nSw+ZVXXuH1118nPDyczMzM065RdHxaax76Yie7s8p45cph9PVtn6GaSrfB4OyYmBh97IUD9uzZw6BBg1r9vVqisbGRqKgoPvvsM/r162fRWkTnYg2fb9H2Fv32Jy/8tJ+Hpg/gtol9W3XbSqk4rXXM8ZZ1mRZ1YmIiffv2ZcqUKRLSQogW+yEhhxd+2s/FwwO4dUKfdn1vixxMtITBgweTnJxs6TKEEB3Q7qxS7v0knmFB7vzrkvB2P7jcZVrUQghxOvLKa7jpvVjce9iz+JpoHO3P7PjW6egyLWohhGipmroGbno/juKqOj67ZTS+Lo4WqUOCWgghjqOxUXP/pzvYmVHCG1dFMzTAzWK1SNeHEEIcx8u/7Oe7Xdk8Mn0g04ac/MzetiZBfZq++uorlFJNY3VP5kxnb1u6dCl33HHHab9eCNEyn8dl8OpvSVweE8S88WGWLkeC+nQtX76csWPHmnXKskyzKUTHsfFAAY9+uZOz+3rxz4uHWsX0AV0mqJ9//vmmaTfvvfdeJk+eDMBvv/3GnDlzWrStiooK1q9fz9tvv910ajEYkyk98MADDB06lIiICF577bXjTrPZfOL5zz//vGmqzm+++YaRI0cyfPhwzjnnnKbpN4UQ7SMpr5ybl8UR4uXEf+dEY29rHRFpmYOJ3z8CObtad5s9w2HGwhMuHjduHC+++CJ33XUXsbGxHDp0iLq6OtatW9c0I565vv76a6ZPn07//v3x8vIiLi6O6OhoFi9eTEpKCvHx8djZ2VFUVISnpycvvfQSq1evxtvb+6TbHTt2LJs2bUIpxZIlS/j3v//Niy++2KLahBCnJ6+8huve3YqDnS3vXn9Wu8zhYa4uM+ojOjqauLg4ysrKcHBwICoqitjYWNatW3fUBPfmWL58OXfffTcAV1xxBcuXLyc6OppffvmFW265pWlWPU9PzxZtNyMjg8svv5zs7Gxqa2sJDQ1t0euFEKen8lA9NyyNpbCilo/njSLQo4elSzqKZYL6JC3ftmJvb09oaChLly5lzJgxREREsHr1apKSkv4yR8PmzZu5+eabAfjHP/5x1DzJRUVF/Pbbb+zatQulFA0NDSileP75582upXmfV/MpSO+8807uu+8+LrjgAtasWcOTTz55mnsrhDBXfUMjd3y0jd1Zpbx1TQyRQe6WLukvrKMDpp2MGzeOF154gfHjxzNu3DjeeOMNhg8f/peDBSNHjmya7vTYyew///xzrr76alJTU0lJSSE9PZ3Q0FDWrVvHueeey5tvvtk0v3FRURHw12k2/fz82LNnD42NjaxYsaLp+dLS0qZrFb733ntt8j0QQhyhtebxrxNYvS+fpy8aypRBfpYu6bi6XFBnZ2czevRo/Pz8cHR0bJpz2FzLly9vuqrKYbNnz2b58uXceOONBAcHExERQWRkJB999BHw12k2Fy5cyKxZsxgzZsxRV3J58sknueyyy4iOjj5lf7YQ4sy9+msSy7ekc/ukPswZ2dvS5ZxQl5rmVIi2Ip/vjufjLWk88uUuZkcF8sJlERYfhifTnAohRDO/JOby9xW7mNDfh4Wz2382vJYyK6iVUu5Kqc+VUnuVUnuUUse/aqgQQli5uNQi7li+jaEBbvx3TpTVjJU+GXNHfbwC/KC1vlQp1Q2wrrErQghhhn055Vz/7lb83brzznVn4eTQMUYon7JKpZQbMB64DkBrXQvUtm1ZQgjRutKLqrjmnc042tvy/twReDs7nPpFVsKcNn8okA+8q5TarpRaopRyOnYlpdQ8pVSsUio2Pz+/1QsVQojTVVBxiGvf2UJ1bQPv3zCCIM+O1SlgTlDbAVHA/7TWw4FK4JFjV9JaL9Zax2itY3x8fFq5TCGEOD1lNXVc+84Wskqreee6sxjY09XSJbWYOUGdAWRorTebHn+OEdxdiq2tLcOGDWu6LVzY/mdXCiFaprq2gRuWbmVfTjn/uyqamJCWTetgLU7ZR621zlFKpSulBmit9wFTgMS2L826dO/enfj4+JOu09DQgK2t7Qkfm/s6IcSZq61v5NYP44hNLebVK4YzaYCvpUs6beaOS7kT+FAptRMYBjzbZhW1kTfeeKOpNRwaGtp0luCZCgkJ4eGHHyYqKorPPvvsL4+XL19OeHg4Q4cO5eGHH256nbOzM/fffz+RkZH88ccfrVKLEMJQ39DIvZ/Es2ZfPs9cFM75kb0sXdIZMWtsitY6HjjuGTOn47ktz7G36NRXRmmJgZ4DeXjEwydcfsstt3DLLbdQV1fH5MmTue+++1q0/erqaoYNG9b0+NFHH+Xyyy8HwMvLi23btgHwyCOPND3Oyspi1KhRxMXF4eHhwdSpU/nqq6+46KKLqKysZOTIkTKNqRCtrLFR88iXu/huVzaPzRzE/40MtnRJZ6xjDCJsRXfffTeTJ0/m/PPPb9HrTtb1cTiwj328detWJk6cyOGDq3PmzGHt2rVcdNFF2NraMnv27JbvgBDihLTWPPnNbj6Py+Cec/px4zjLX0arNVgkqE/W8m1LS5cuJTU1lUWLFv1l2cmmNj0VJyenkz4+HkdHR+mXFqIVaa1Z+P1e3v8jlXnjw7h7Sj9Ll9RqukyLOi4ujhdeeIF169ZhY/PXrvnDU5u2phEjRnDXXXdRUFCAh4cHy5cv584772zV9xBCGP7z837eXJvM1aN68+iMgVY/f0dLdJmgXrRoEUVFRU0HEWNiYliyZInZrz+2j3r69OmnHKLn7+/PwoULmTRpElprZs6cyYUXXnha9QshTmzRb382XTX8qQuGdKqQBpnmVIhWIZ9vy3nz9wP86/u9XDI8gOcvi8TWpmOGtExzKoTolJasS+Zf3+9lVoQ//740osOG9KlIUAshOqS31x/kn9/tYWa4Py9fPgy7DjBd6elq1z1ri24WISxNPtft790NB3n620RmDO3Jy1d07pCGdgxqR0dHCgsL5UMtOhWtNYWFhTg6Olq6lC5jybpknvomkelDevLqlcM7xMT/Z6rdRn0EBgaSkZGBTIEqOhtHR0cCAwMtXUaX8NbaZJ5ZtYfzwnvyyhVdI6ShHYPa3t6e0NDQ9no7IUQn88bvB1j4/V5mRhh90l0lpKELjaMWQnRMWmte/TWJ//yyn/Mje/Gfv0V2+j7pY0lQCyGsltaaF37ax+urDzA7KrBTD8E7GQlqIYRV0lrzzHd7WLL+IFeOCOKZi8Kx6YIhDRLUQggr1Niomf9VAsu3pHHt6N4sOH9Ilw1pkKAWQliZ+oZGHvx8Jyu2Z3LbxD48OG1Ap5u7o6UkqIUQVuNQfQN3Ld/Oj7tzeXDaAG6f1NfSJVkFCWohhFWoPFTPvGWxbEgqZMH5g7n+bBnOe5gEtRDC4kqqarnu3a3syizlxcsimR0tJxA1Z1ZQK6VSgHKgAag/0VR8QgjRUjmlNVz7zhYOFlTy3zlRTBvS09IlWZ2WtKgnaa0L2qwSIUSXk5xfwdVvb6Gkqpal15/FmL7eli7JKknXhxDCInZllHLtu1tQwMfzRhMe6GbpkqyWuedhauAnpVScUmre8VZQSs1TSsUqpWJl4iUhxMms+zOfKxb/QY9utnx+6xgJ6VMwN6jHaq2jgBnA7Uqp8ceuoLVerLWO0VrH+Pj4tGqRQojO46vtmVz/7laCPHvwxa1jCPV2snRJVs+soNZaZ5ru84AVwIi2LEoI0florXlrbTL3fBJPTIgHn94yGj9XmcfbHKcMaqWUk1LK5fDXwFQgoa0LE0J0Hg2Nmqe+SeSZVcals96bOwJXR3tLl9VhmHMw0Q9YYTqF0w74SGv9Q5tWJYToNGrqGrj7Y+NswxvGhjL/vEFdet6O03HKoNZaJwOR7VCLEKKTKaw4xE3vx7I9vYQnZg1m7lg52/B0yPA8IUSbOJBfwfXvbiW3rIb//l8UM8L9LV1ShyVBLYRodZuTC5m3LA47G8XyeaOICvawdEkdmgS1EKJVfRGXwaNf7iLIszvvXjeCYK8eli6pw5OgFkK0isZGzUs/72fR6iTG9PHif3OiceshIztagwS1EOKMVdc28MBnO/huVzZXjgjiHxcO7VJXCW9rEtRCiDOSXVrNTe/HsjurjPnnDeLGcaFd/oosrU2CWghx2ranFTNvWRxVh+pZck0MUwb5WbqkTkmCWghxWlZsz+DhL3bh5+rABzeczYCeLpYuqdOSoBZCtEhDo+a5H/ayeG0yI0M9+d9V0Xg6dbN0WZ2aBLUQwmylVXXc+fF21u7P59rRvXls1mA5aNgOJKiFEGbZl1POvGWxZJVU869LwrlyRLClS+oyJKiFEKf03c5sHvx8B04Odnw8bxTRvT0tXVKXIkEthDih+oZGnv9pH2/+nkxUsDv/uypa5pC2AAlqIcRxFVQc4s6PtvNHciFzRgaz4PwhdLOT/mhLkKAWQvzFtrRibvtgG8VVtbxwWSSXRgdauiSrl1Kawp6iPcwIndHq25agFkI00VqzdGMKz3y3B393R764dQxDA+TCs6eSUprC3B/notGMDxyPk33rXgdSgloIAUB5TR2PfLGL73Zlc84gP168LFImVTLD4ZBu0A28PfXtVg9pkKAWQgC7s0q546PtpBVV8ciMgdw8Pkzm6zBDalkqN/x4Q1NI9/Xo2ybvI0EtRBemteajLWk89U0iHj3s+ejGkYwM87J0WR1CcmkyN/54Iw26gSVTl7RZSEMLglopZQvEApla61ltVpEQol2U19Qxf0UCK3dkMa6fN/+5fBjezg6WLqtDSCpO4safbgRo05b0YS1pUd8N7AFc26gWIUQ72ZlRwp3Lt5NRXM0DU/tz28S+cmVwM+0v3s9NP92ErbJlybQlhLmFtfl7mjUoUikVCMwElrRtOUKIttTYqFmyLpnZ/9tIXX0jH88bxR2T+0lIm2l34W7m/jgXOxs73p3+bruENJjfon4ZeAg44TyGSql5wDyA4GCZA0AIa5NffogHPtvB7/vzOXewH89fGoF7D5n1zlzxefHc+sutuDm48dbUtwhyCWq39z5lUCulZgF5Wus4pdTEE62ntV4MLAaIiYnRrVWgEOLMrdmXxwOf7aC8pp6nLxzCVaN6y6iOFtiSvYU7frsD3x6+LJm6hJ5OPdv1/c1pUZ8NXKCUOg9wBFyVUh9ora9q29KEEGeqpq6Bhd/vZenGFPr7OfPhjaNkgv8WWpO+hvvX3E+wazBvTX0L7+7e7V7DKYNaa/0o8CiAqUX9gIS0ENZvb04Zdy+PZ19uOdeNCeGRGQNxtLe1dFkdynfJ3zF//XwGeg7kjXPewN3R3SJ1yDhqITqZhkbN2+uTeeHH/bh2t2fp9WcxcYCvpcvqcD7d9yn/3PRPov2ieW3yazh3c7ZYLS0Kaq31GmBNm1QihDhjGcVV3P/pDjYfLGLaED+evTgcLxkb3SJaaxbvXMyi+EVMCJzACxNewNHOslO7SotaiE5Aa81nsRn849tEAJ6/NIJLowPlgGELNepGntvyHB/t/Yjzw87nqbOfwt7G8vOdSFAL0cHlldXw6Je7+HVvHiNDPXnhskiCPHtYuqwOp66hjsc2PMaqg6u4evDVPBDzADbKOubflqAWooPSWrNyRxYLVu6muraBJ2YN5roxIXLyymmorKvkntX3sCl7E3dH3c0NQ2+wqv9GJKiF6IDyymt4bEUCPyXmMjzYnecvjaSvr+UOdnVkBdUF3PbLbewv3s8/xvyDi/tdbOmS/kKCWogORGvNV/GZPPVNIlW1DTw6YyA3jgvDVlrRpyWlNIVbf7mVwppCXp38KuMDx1u6pOOSoBaig8gurWb+igR+25tnakVH0NdXTl45XdvztnPnb3diq2x5e+rbhPuEW7qkE5KgFsLKNTZqlm9NY+GqvdQ1NvK4qS9aWtGn7+fUn3lk7SP4O/vzvyn/I8i1/ebtOB0S1EJYsQP5FTz65S62HCxiTB8v/nVJOL29Wv9ST12F1pr3dr/HS3EvEeETwWuTX8PD0cPSZZ2SBLUQVqi2vpE3fz/Aa6uTcLSz4d+zI7gsRsZFn4n6xnqe3fwsn+3/jHN7n8uzY5+1+Iks5pKgFsLKxKYU8eiXu/gzr4KZ4f4suGAwvi4dI1CsVXltOQ/+/iAbsjYwd+hc7o6622rGSJtDgloIK1FSVctzP+xl+ZZ0Aty78851MUwe6Gfpsjq8jPIM7vj1DlLLUnly9JPM7j/b0iW1mAS1EBamtebLbZk8s2oPpdV13Dg2lHvP7Y+Tg/x6nqltudu4Z/U9NOgG3jj3DUb6j7R0SadFPglCWNC+nHIe/zqBLQeLiAp255mLwxnkL5clbQ0r/lzB05ueJsA5gEVTFtHbtbelSzptEtRCWEDFoXpe+WU/72xIwcXRjn9dEs7lMUFy+ncrqG+s56W4l1iWuIxR/qN4YcILuDm4WbqsMyJBLUQ70lrzdXwWz67aQ175Ia4cEcRD0wbi4STXLmwNpYdKeWjtQ2zM2shVg67i/pj7sbPp+DHX8fdAiA4iMauMJ1fuZktKERGBbrx5dTTDg61/DG9H8Wfxn9y9+m6yK7N5asxTXNLvEkuX1GokqIVoY0WVtbz40z6Wb0nDrbs9Cy8J52/SzdGqfk79mfnr5+Nk78S7095lmO8wS5fUqiSohWgjdQ2NfLAplf/8vJ/K2gauGR3Cvef0x62H5Sei7ywaGhtYFL+IJbuWEOETwX8m/gffHp3vsmMS1EK0Mq01a/bl8/R3iSTnVzK2rzdPnD+Y/n4ygVJrKq4p5qG1D7EpexOz+83m7yP/TjfbztnXL0EtRCvam1PGM9/tYd2fBYR6O/H2tTFMHugrp363st0Fu7l3zb0UVhd2uv7o4zllUCulHIG1gINp/c+11gvaujAhOpK88hr+8/N+PtmajrODHY/NHMQ1o0PoZtdxTlPuCLTWfLrvU57b+hze3b15f8b7DPEeYumy2pw5LepDwGStdYVSyh5Yr5T6Xmu9qY1rE8LqVR6q5611ySxem0xtfSPXjgnhrsn9ZLhdG6iqq+Ifm/7Bd8nfMTZgLP8a+y/cHd0tXVa7OGVQa601UGF6aG+66bYsSghrV9/QyKexGfznl/3klx/ivPCePDhtIKHeMgVpW0gqTuL+3+/nYOlB7hh2BzdF3NShJlU6U2b1USulbIE4oC/wutZ683HWmQfMAwgODm7NGoWwGlprfkjI4fkf95FcUElMbw/evDqaKBkP3Wa+SvqKZzY9g5O9E4unLmaU/yhLl9TuzApqrXUDMEwp5Q6sUEoN1VonHLPOYmAxQExMjLS4RaezIamAf/+4jx3pJfTzdeata2I4Z5AcKGwrVXVVPLP5GVYeWMnIniNZOH4h3t29LV2WRbRo1IfWukQptRqYDiScan0hOoP49BKe/3EvG5IK6eXmyL9nR3BJVAB2tl3nX+/2tqdwDw+ufZD08nRui7yNeRHzsLWxtXRZFmPOqA8foM4U0t2Bc4Hn2rwyISxsT3YZL/60n1/25OLp1I3HZw1mzshgHO27bmC0Na01H+75kJfiXsLD0YMlU5dwVs+zLF2WxZnTovYH3jP1U9sAn2qtv23bsoSwnKS8cl7+5U++3ZmNi6Md95/bn+vHhuIs80O3qYLqAh7f8DjrM9czIXACT5/9dIe4nmF7MGfUx05geDvUIoRFJedX8Oqvf/L1jiy629ty+6Q+zBvXR075bgdrM9by+IbHqayr5O8j/84VA66Qvv9mpIkgurwD+RUs+i2Jr+MzcbCzZd74MOaNC8PL2cHSpXV61fXVvBj7Ip/s+4S+7n1ZMnUJ/Tz6Wbqs01d/COxa/3MjQS26rKS8chb9lsTKHVl0s7PhhrGhzBvfBx8XCej2sLtgN4+se4SUshSuGXwNd0XdhYNtB/veVxdD6kY4uA5S1kN9DdwZ2+pvI0Etupw92WUs+i2JVQnZONrZctO4MG4aH4a3tKDbRV1jHUt2LWHxjsV4dvfkralvdZyx0VVFkPYHpGyAlHWQswvQYOcIQSMgZDw0NkArj1CRoBZdxra0Yv67Oolf9uTh7GDHbRP7cMPYMDzldO92k1yazPx180koTOC80PP4+8i/W/dlsioLIHWDEcypGyB3N03BHHgWTHwUQsZCYEybdHkcJkEtOjWtNeuTCvjv6gP8kVyIew977jmnH9ePCZWDhO2oobGBD/Z8wGvbX6O7XXdenPAiU0OmWrqsvyrNNLoy0jYa9/l7jeftuhst5knzIeRsCIhu02A+lgS16JQaGo1Tvf/3exIJmWX4ujjw2MxBXDkiGCcZZteuUkpTeHzD48TnxzMpaBJPjH7COs4w1BoKk0zB/IfRYi5JM5Y5uELQSIi43Ggx+w8DO8v95yWfWNGp1NQ18FlcBkvWJZNaWEWYtxPPzQ7nouEBONjJiSrtqaGxgWWJy1gUvwgHWweeHfsss8JmWW7YXUMd5OyE1D+MYE7bBFUFxrIe3tB7DIy6zbj3G9rq/cxnQoJadAqFFYf4YFMa7/+RQmFlLZFB7jwyfSBTh/TEVq5N2O6SipN4YuMT7CrYxaSgSTw+6nF8evi0bxE1pZCxFdI2G8GcGQd1VcYy997Q71wIHm0Es1dfsOJx2xLUokM7kF/B2+sP8kVcBofqG5ky0Jd548MYEeopJ0xYQG1DLUt2LeGtXW/hYu/C8+OfZ1rItLb/WWgNJamQvsVoKadvPnLgT9lAz3AYfjUEjzLC2dW/betpZRLUosPRWrMhqZC31yezel8+3exsmB0VwA1jw+jr62zp8rqs+Lx4FmxcQHJpMjPDZvLQWQ/h6ejZNm9WfwiydxiBnL7ZCOiKXGNZNxdjFMbER4x+5sAYcOjY16uUoBYdRnVtA1/FZ7J0Qwr7csvxdu7Gvef0Z86oYBkDbUFltWW8EvcKn+3/jJ5OPfnvlP8yLnBc675JaSZkbIH0rcZ99g5oqDWWeYRA6AQIHmkEs+9gq+pfbg0S1MLqZRRX8cGmND7emkZJVR2D/V15/tIILhjWSw4QWpDWmh9TfuS5rc9RVFPEnEFzuHP4nfSw73FmG66rMYI4Y4vRx5y+FcqzjGV2jtBrOIy8xRguFzgCXPzOfGesnAS1sEpaazYeKOS9jSn8ssf4l/bcwX7MPTtU+p+tQGpZKs9ufpaNWRsZ7DWY16e8zmCvwS3fkNZQeAAyYyEj1rjP2QWN9cZy997Qe7QRyEFngV+4RYfJWYoEtbAqZTV1fBmXwbJNqRzIr8TTqRu3TOjDnFG9CXDvbunyurya+hreSXiHt3e9TTfbbjwy4hGuGHCF+ZP6VxYYoy8yYo37zDioKTGWdXM2Wstj7jTO+gs8C5x922xfOhIJamEVEjJL+WBTKl/HZ1Fd18CwIHdevCySmRH+MlG/lViTvoaFWxaSWZHJjNAZPBjz4MmH3B2qgOx4yNxmBHLWtiMnlCgb8B0Cgy80DvYFxIDPgE7Xt9xaJKiFxVQequfbnVl8tDmNHRmlONrbcGFkAHNGBRMR6G7p8oRJWlkaz219jrUZawlzCzv+JEr1hyAnwQjjrO1GOBfsA91oLHcPNk67PusmI5j9I6GbXLHdXBLUot0lZJayfEsaX8dnUXGonn6+ziw4fzCXRAXi1l3m37AWVXVVLN65mPcT38fexp77o+9nzqA52GsNWfFGIGeb7nMTobHOeGEPbwiIMlrLAdFGd4ZzO5/s0slIUIt2UVpdx8odWXy8JY3dWWU42NkwM9yf/xsZTHRvDzk4aEUadSPfJX/Hy3Evk1edxwX+Y7nHdQg+B+Nh41LjRJLDQ+Mc3Y/0K/caBr2iwC3Qqs/y64gkqEWbaWzUbDpYyKdb0/k+IYdD9Y0M8nflHxcO4cJhAdJ6tjZ11cTv+5rndi8hoSaXIY22vJRbQOTBj4zlDm7QK9IYGtdruHHzCJFQbgfmXIU8CHgf8AM0sFhr/UpbFyY6rvSiKj6Py+CLbRlkFFfj4mjHZTGB/C0miPAAN2k9W4OaUmMYXPZOyNlJRs52XtEF/ODUA9/6ep6tqGOm+yBsoi8w+pP9h4FnmISyhZjToq4H7tdab1NKuQBxSqmftdaJbVyb6EDKa+r4flcOn2/LYMvBIpSCsX29eWDqAKYP7SkjNyxFayjPMUI5Z0dTMFOcAkCZjeItH38+7GGPrXLhZv9xzB12Oz28B0goWxFzrkKeDWSbvi5XSu0BAgAJ6i6uvqGRdUkFrNiWyU+JOdTUNRLq7cQDU/tzcVSgjHtubw11UPAn5CYYYZyTYAT04ak8ATxCwT+S2mFXspxy3spaTVltBRf0uYA7ht9BT6eelqtfnFCL+qiVUiHAcGDzcZbNA+YBBAcHt0ZtwgpprdmRUcrX8Zl8syObgopDuPew59LoQC6JCmR4kLt0bbSHqiIjkHN3G4Gcuwvy9hw5yGfbzZjzYsB06BlpzB7nN4RGB2dWHVzFou2LyKzIZLT/aO6NvpdBXoMsuz/ipMwOaqWUM/AFcI/WuuzY5VrrxcBigJiYGN1qFQqrkJRXwcodWXyzI4uDBZV0s7Nh8gBfLo4KYNIAX7rZ2Vi6xM6poc64Cknu7iPBnLsbyjKPrNPD2wjikTcbp1j3DAfvfmB75GCt1pp1met4Zdsr7C/ezwCPAbx5zpuMCRhjgZ0SLWVWUCul7DFC+kOt9ZdtW5KwFhnFVXy7M5tvdmSxO6sMpWB0mBe3TujDtKE9ZdRGa9IayrIgL9EI4rxEY2xywb4jrWQbO/AZaFwayncw9BxqBPMpJiWKy43j1W2vsi1vG4HOgTw37jmmh07HRskf147CnFEfCngb2KO1fqntSxKWlF1azapdOXy3M4ttaSUARAa58/iswcyK8MfP1dGyBXYGVUVGN0Veoune9PXhOS8AXHqB3xDoO9m4LJTfEPDq16IJiXYX7ua17a+xIXMD3t29mT9yPrP7zcbeVv7AdjTmtKjPBq4Gdiml4k3P/V1rvarNqhLtKqukmh8Scli1K5vY1GIABvm78tD0AcwK70Ww1xlOW9lV1ZRB/j7I39MskPdARc6RdRzcwHcQDLnYCGPfwcbjHqc/4f7eor38N/6/rE5fjZuDG/dF38cVA6+gu50c3O2ozBn1sR6Qo0OdTFphFT/szub7hBy2m1rOA3u68MDU/pwX7k+Yj1wpxWyHyiF//5FAzt8H+XuhNP3IOnbdjUmH+kw6Esa+g8A1oNWGwe0r2sebO9/k59SfcbF34fZht3PVoKtw7iY/y45OzkzsIrTW7M+t4MfdOfyQkENitnE8eEgvVx6cNoAZQ3tKOJ9KVREU7DcFsSmM8/dBWcaRdewcjQN5waPA5zojjH0GGmfwtdHMcHuL9vLGjjf4Ne1XnO2duTniZq4Zcg2u3Vzb5P1E+5Og7sQaGjXb04r5KTGXH3fnkFpoXIE5Ktid+ecNYvrQngR5SrfGUQ4f1CvYfySUD99X5h1Zz647+PQ3rmDtOxB8Bhkt5jYM5GPtyt/F4p2LWZOxBhd7F26NvJU5g+bg5uDWLu8v2o8EdSdTVVvP+j8L+GVPLr/uyaOwshZ7W8WoMC9uGhfG1MF++MoBQWNazqKDRwK56fYn1FYcWc/RDbz7Q/+p4D3ACGPv/saVR2zaf9SE1prY3Fje2vkWf2T/gZuDG7cPu53/G/R/0oLuxCSoO4Hs0mp+3ZPHb3vz2JBUwKH6Rlwc7Zg4wJepg/2YMMAHV8cueKRfa6jMN8K38E/TfZIRyMWpoBuOrOsaYHRZDJtj3PsMMILZ2dcqTqVu1I2szVjLW7veYmf+TjwdPbk3+l4uH3A5TvYyr3NnJ0HdATU0auLTi1m9N5/f9uY19TcHeXbn/0YGc+4gP84K9cTetouMk62tNK67V5jU7N4UyjWlR9azcwTPPsYJIUNnGy1jr75GMDu4WK7+k6hrqGPVwVW8m/AuB0oPEOAcwGMjH+PCvhfiaCf/GXUVEtQdRGHFIdb+mc+affms3Z9PcVUdtjaK6GAPHpkxkCkDfenr69x5T9+urzUmEio6YIRxUbNgbn6WHoBrIHj1gfDLjLHH3n2Ne7cgi3RXnI7y2nK+2P8FH+z5gNyqXPp79Odf4/7F9JDp2NnIr21XIz9xK1Xf0Eh8egm/7zeCeWdmKVqDl1M3Jg3wZdJAX8b388GtRyfq0mioM66p1xTEze5L049c1gmgu4cRvqHjjVayd1+jdewZ1qEv8ZRVkcWHez7kiz+/oLKukrN6nsWC0QsYGzC28/4RFqckQW1F0gqrWJdkBPPGpELKD9VjoyAq2IN7z+nPxAE+DO3lho1NB/6FrasxWsbFB6Eo+ehbSfrR/cYOrkbwBsZAxOVGK/lwGJ/BCSHWKD4vnmWJy/g17VcApoZM5doh1zLEa4iFKxPWQILagkqqavnjQCHrkgpY/2cBaUXG8LkA9+7MivRnXD8fzu7j3fFazdXFxoiKpkA23YoPGkPfaDZnl4MbeIYal3AKv8wIYc8wo5Xs5G0VB/LaSl1DHT+m/shHez5iV8EuXOxduGbwNVw58Er8nf0tXZ6wIhLU7aiqtp7YlGI2Hihk44ECdpm6M5y62TK6jxdzzw5hbD8f+vg4Wfe/uQ31xkkexSlH3w6Hc/M5KwCcfIzwDRlnCuFQY15krz5GF4Y172sbyK/K5/P9n/Pp/k8pqC4gxDWER0c8ykV9L6KHvYxrF38lQd2Gauoa2J5Wwh/JhWw6UMj29GLqGjT2torhQR7cM6U/Y/t5ERHobl0jNLSGqkJjCFtJinFfnAIlh++P6aKwsTMO1HmGGled9gg5EsYeva12REV70lqzPW87H+/9mJ9Tf6Ze13N2wNk8PehpxvQaIzPZiZOSoG5FNXUNbEsrZnNyEZuSC9meXkJtfSM2CoYGuDH37FDG9PXmrBAPenSz8Le+usQ4cFeSZgRwSZopmFON+7rKo9fv4W0EcEC0MbTNI8Q46cMz1BiD3E5n43U0lXWVfHvgWz7e9zFJJUm42Ltw5aAruWLAFQS7ygU2hHkkqM9AeU0dcanFbE0pYsvBInakl1LbYATz4F6uXDOqN6PCvDgr1LN9527W2ugnLkkzRkuUpBmt4KZgToNDpUe/ppuzEbweIRA6wWgJu/c23QdLq7iFdhfu5rN9n7Hq4Cqq66sZ5DmIJ0c/yYzQGdK9IVpMgroF8spqiE0tZsvBImJTi0jMKqNRg62NYmiAG9efHcLIME+ie7dxMDc2GlNllmY0C+P0o++bnwYNpiAONroogkceCWG3ICOcu2BfcWsrry3n+4Pf88WfX5BYmIijrSPTQ6dzWf/LCPcOt+7jDsKqSVCfQEOjZn9uOXGpxcSlFhObWkR6UTUAjvY2DA/y4I5JfRkR6sXwYHecHFrxW3mowgjhsgzjvvmtJM0YOdFYd/RrHN3BPcg4WBc20fjaLcgIZ/dgCeI2orUmLjeOFUkr+Dn1Z6rrq+nv0Z9HRzzKrD6zZP4N0SokqE1Kq+rYnl7MtrQStqcVsz2thIpD9QB4OzsQ09uDa0eHEBPiyZBerqd/8K/+kHEmXWmm6T7jmMfpR5/2DKBswcXfCN+gEeAWaISwW5DxtXuQdE20s5zKHL5N/pYVf64grTwNJ3snzgs9j0v7X8oQryHSehatqksGdV1DI/tyytmeXkJ8Wgnb04tJzjcOntkoGNDTlYuG9yK6twdRwR4Ee/Yw7xevrgbKs43ALcs6cn84hMsyjUmCjtXd40j4Bo86EsCuAUYIO/cE2y75o7Iq1fXV/Jb2G18nfc2m7E1oNNF+0dwceTPnBJ8jfc+izXT6336tNSmFVezMKCE+vYSdGaUkZJZyqN44HdnbuRvDgjyYHRXI8GB3IgOP042htdHKLcuG8ixTCDf/OtN4XFXw1wIc3MAtAFx7gX+EMQ+FW4ARwm6BxvMd+JTnzq6hsYHY3Fi+OfANP6f+TFV9Fb2cenFL5C2c3+d8glyCLF2i6AI6VVBrrckormZXZim7MkvZmVHCroxSymqMLgxHexuG9nLjqlG9GRbkzrAgdwJdbVEVeaaWcAJszzEFcLbxXHm28fWxw9UAengZFyF19YeAGCN8XU2PXQONe+mS6HC01uwr3se3B77l+4Pfk1edh5O9E9NDpzMrbBbRftEy7lm0K3OuQv4OMAvI01oPbfuSzKO1Jq2oioTMMhKyjFZyQmYpxVXGQbbutg2M9Gnk5n61hLtW0697BX6qGJvKXCjNgYxs+PEErWDbbuDS0+gX9hsK/aYaX7v2OvreXqaZ7ExSy1JZdXAV3x/8noOlB7FTdowNGMuDYQ8yIWiCXBxWWIw5LeqlwCLg/bYt5cTqGhr5M7eCxOwy9mXkkZ2RQkleBs51hfiqYvxsSrmqeyW9e5Tj61SMS10hdjWFqBINJc02pGyM05ldehqt34Boo0V8OJRd/Y3HPTxlhEQXkVmRyY8pP/LDwR/YU7QHhSLaL5qrBl3F1N5TcXd0t3SJQph1FfK1SqmQti5ENzaSvHsrlUWZHCrOpqEsh0OlOVCeg0NNAd6UMFWVcKkyJi5CAd1Mr1U2KAc/cPYDl1BwGW0cgHPxMwLY2XTv5CMH5QSZFZn8lPITP6X8REJhAgDh3uE8EPMA00Km0dOpp4UrFOJoVpVaAZ/PxFEdGR9ciSPldp7Uuflg4zKURq8AGr0DsXHxM1rBzkY4KydvOYVZnFRKaQq/pP3Cz6k/k1iYCMBgr8HcE3UP00KmEegSaOEKhTixVgtqpdQ8YB5AcHDL5zBQNjbsm7CIbk6euHgH4O4bgLOLOzIeQpwOrTWJRYn8mvorv6X9xoHSAwBEeEdwX/R9nNP7HBmxIToMpbU+9UpG18e35h5MjImJ0bGxsWdYmhAtU9dQx9bcraxOW82ajDXkVOZgo2yI9otmSvAUpgRPkW4NYbWUUnFa65jjLbOqrg8hWqq4ppj1metZk76GjVkbqairwNHWkdG9RnNb5G1MDJqIh6OHpcsU4oyYMzxvOTAR8FZKZQALtNZvt3VhQhzP4THOazPWsjZjLTvzd6LR+HT3YVrINCYETmBUr1EylE50KuaM+riyPQoR4kRKD5WyOXsz6zPXsz5zPfnVxmn4Q7yGcEvkLYwPHM9gr8FyEorotKTrQ1id+sZ6dhfuZmPWRjZmbmRnwU4adSMu9i6MCRjD2ICxjA0Yi3d3b0uXKkS7kKAWFqe1Jr08nU3Zm/gj6w8252ymvLYchWKw12BuDL+RsQFjCfcOx85GPrKi65FPvbCIvKo8NmdvZkvOFjZnbya7MhuAnk49OSf4HMb0GsNI/5FyIFAIJKhFOymoLiA2N5at2VvZkrOFlLIUAFy7uTKi5wiuH3o9o/1H09u1t8zlLMQxJKhFm8ipzCE2N5a43Dhic2KbgtnJ3olov2gu7X8pI3qOYIDnADkIKMQpSFCLM9aoG0kuSWZb3ja2521nW+42siqzAHC2dybKL4pL+l1CjF8Mg7wGST+zEC0kvzGixarqqkgoSCA+P574vHh25O+grLYMAC9HL6L8orhmyDVE+UbR36M/tjIPixBnRIJanFSjbiS1LJWd+TuNW8FO9hfvp1EbV8jp49aHc3ufyzDfYUT5RhHkEiR9zEK0MglqcZT8qnx2F+5mV8EuEgoS2FWwi/LacsDoxhjqPZSbwm8i0ieSCJ8I3BzcLFyxEJ2fBHUXVlRTRGJhIrsLdpNYmEhCYQJ5VXkA2Cpb+nn0Y1rINMK9w4n0iSTULVQO/AlhARLUXYDWmtyqXPYU7mFv0V72FO0hsTCR3KrcpnVCXEOI8YthqPdQhnoPZaDnQJkvQwgrIUHdydQ11JFcmsz+4v3sK9rH3uK97C/aT/GhYgAUit6uvYn2i2aw12AGew1moOdAXLrJRXiFsFYS1B3U4Vbyn8V/8mfJn+wv3s+fxX+SXJpMfaNx1fVuNt3o59GPycGTGeA5gEGeg+jv0Z8e9j0sXL0QoiUkqK2c1pr86nwOlBzgQMkBkkqSmu4r6iqa1vPr4Uc/j36MCxhHf4/+9PfoT4hbiIxZFqITkN9iK9HQ2EBWRRbJpckcLD1IcmmycStJpryuvGk9Nwc3+rj1YWbYTPq596OPex/6efST0RdCdGIS1O1Ia03JoRJSy1JJLUslpSyFlNIUUspSSC1Lpa7xyIV9PR09CXML47yw8whzC6OPex/6uPfBy9FLxikL0cVIULcyrTXFh4pJL08nrSyN9PJ0UstSSStLI608rekMPgA7ZUegSyAhriGMCxhHqFsoIW4hhLmFSQtZCNFEgvo01DbUkl2ZTWZ5JhkVGcatPIP08nQyyjOO6jtWKPyd/Al2DWZG6AyCXYIJcQsh2CWYAJcA7G3sLbgnQoiOQIL6OGobasmpzCGzIpPsymyyKrLIqsgisyKTzIpM8qry0By5eru9jT0BzgEEuQQx3Hc4gc6B9HbtTZBrEAHOATjYOlhwb4QQHV2XC+qGxgYKawrJqcw5cqs68nV2ZTYF1QVHvcZG2eDXww9/J39G+o8kwDmAAOcAAl0CCXQOxKeHj5yxJ4RoM2YFtVJqOvAKYAss0VovbNOqTlNVXRUF1QXkVeU13XKrco+6z6/Kp17XH/U6R1tHejr1xN/Jn/GB4+np1JNeTr3o5dwLfyd//Jz8pItCCGExpwxqpZQt8DpwLpABbFVKrdRaJ7Z1cWDM3lZ2qIyC6gIKagqM+yrjPr86v+k+vyr/qL7hw7rbdcevhx9+Pfw4q+dZ+PbwpWePnvg5+dHTqSc9e/TEzcFNRlIIIayWOS3qEUCS1joZQCn1MXAh0OpB/bdv/kZNQw1gjJ6oqKuguKaYBt3wl3UdbR3x7u6Nd3dv+rr3ZbT/aHx6+ODbwxef7sa9bw9fnO2dJYSFEB2aOUEdAKQ3e5wBjDx2JaXUPGAeQHBw8GkVE+YeRl1DXVOwOts74+noiVd3LzwdPZuC2bu7twSwEKLLaLWDiVrrxcBigJiYGH2K1Y9r4Tir7PoWQgiLMmeoQiYQ1OxxoOk5IYQQ7cCcoN4K9FNKhSqlugFXACvbtiwhhBCHnbLrQ2tdr5S6A/gRY3jeO1rr3W1emRBCCMDMPmqt9SpgVRvXIoQQ4jjkdDohhLByEtRCCGHlJKiFEMLKSVALIYSVU1qf1rkpJ9+oUvlA6mm+3BsoOOVanYvsc+fX1fYXZJ9bqrfW2ud4C9okqM+EUipWax1j6Trak+xz59fV9hdkn1uTdH0IIYSVk6AWQggrZ41BvdjSBViA7HPn19X2F2SfW43V9VELIYQ4mjW2qIUQQjQjQS2EEFbOYkGtlJqulNqnlEpSSj1ynOUOSqlPTMs3K6VCLFBmqzFjf+9TSiUqpXYqpX5VSvW2RJ2t6VT73Gy92UoprZTq8EO5zNlnpdTfTD/r3Uqpj9q7xtZmxmc7WCm1Wim13fT5Ps8SdbYWpdQ7Sqk8pVTCCZYrpdSrpu/HTqVU1Bm/qda63W8Y06UeAMKAbsAOYPAx69wGvGH6+grgE0vU2o77OwnoYfr61o68v+bus2k9F2AtsAmIsXTd7fBz7gdsBzxMj30tXXc77PNi4FbT14OBFEvXfYb7PB6IAhJOsPw84HtAAaOAzWf6npZqUTddMFdrXQscvmBucxcC75m+/hyYojruRRJPub9a69Va6yrTw00YV9LpyMz5GQM8DTwH1LRncW3EnH2+CXhda10MoLXOa+caW5s5+6wBV9PXbkBWO9bX6rTWa4Gik6xyIfC+NmwC3JVS/mfynpYK6uNdMDfgROtoreuBUsCrXaprfebsb3M3YPxF7shOuc+mfwmDtNbftWdhbcicn3N/oL9SaoNSapNSanq7Vdc2zNnnJ4GrlFIZGPPa39k+pVlMS3/fT6nVLm4rWodS6iogBphg6VraklLKBngJuM7CpbQ3O4zuj4kY/zWtVUqFa61LLFlUG7sSWKq1flEpNRpYppQaqrVutHRhHYWlWtTmXDC3aR2llB3Gv0yF7VJd6zPrAsFKqXOA+cAFWutD7VRbWznVPrsAQ4E1SqkUjL68lR38gKI5P+cMYKXWuk5rfRDYjxHcHZU5+3wD8CmA1voPwBFj8qLOqtUvCG6poDbngrkrgWtNX18K/KZNPfUd0Cn3Vyk1HHgTI6Q7er8lnGKftdalWmtvrXWI1joEo1/+Aq11rGXKbRXmfK6/wmhNo5TyxugKSW7HGlubOfucBkwBUEoNwgjq/Hatsn2tBK4xjf4YBZRqrbPPaIsWPHJ6HkZr4gAw3/TcPzB+WcH4YX4GJAFbgDBLH+1t4/39BcgF4k23lZauua33+Zh119DBR32Y+XNWGF0+icAu4ApL19wO+zwY2IAxIiQemGrpms9wf5cD2UAdxn9INwC3ALc0+xm/bvp+7GqNz7WcQi6EEFZOzkwUQggrJ0EthBBWToJaCCGsnAS1EEJYOQlqIYSwchLUQghh5SSohRDCyv0/OVVTPuz0O8wAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "t0 = 0\n", "tf = 1\n", "N = 100\n", "\n", "#Every variable is the same as the one given in Forward_Euler.py\n", "\n", "t = np.linspace(t0,tf,N)\n", "y = np.zeros(N)\n", "z = np.zeros(N)\n", "w = np.zeros(N)\n", "h = (tf-t0)/N\n", "\n", "y[0] = 1\n", "y[1] = y[0] + h*y[0]\n", "\n", "for i in range (0,N-2):\n", " y[i+2] = y[i+1] + h*(1.5*y[i+1]+0.5*y[i])\n", "\n", "w = np.exp(t) #we vectorized the subtraction instead of iterating through w\n", "z = abs(y-w) #same here\n", "\n", "plt.plot(t,y, label = \"y - Numerical Method method\") #very simple plot of y,z,w\n", "plt.plot(t,w,label = \"w - Actual\")\n", "plt.plot(t,z,label = \"z - Error\")\n", "plt.title(\"Method with N=%d\"%N)\n", "plt.legend()\n", "plt.show()\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "interpreter": { "hash": "b0fa6594d8f4cbf19f97940f81e996739fb7646882a419484c72d19e05852a7e" }, "kernelspec": { "display_name": "Python 3.9.10 64-bit", "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.9.10" }, "orig_nbformat": 4 }, "nbformat": 4, "nbformat_minor": 2 }