{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Illustrates function iteration, Newton, and secant methods\n", "\n", "**Randall Romero Aguilar, PhD**\n", "\n", "This demo is based on the original Matlab demo accompanying the Computational Economics and Finance 2001 textbook by Mario Miranda and Paul Fackler.\n", "\n", "Original (Matlab) CompEcon file: **demslv06.m**\n", "\n", "Running this file requires the Python version of CompEcon. This can be installed with pip by running\n", "\n", " !pip install compecon --upgrade\n", "\n", "Last updated: 2022-Sept-04\n", "
\n" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "plt.style.use('seaborn-dark')" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "bullet_options = dict(ms=16, marker='.', ls='')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Function Iteration" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfgAAAINCAYAAADSn21IAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAABYaUlEQVR4nO3dd3hUVf4G8PdOTe+VGiCAQOhIFwgIFkSxYFlhdUUQC7DoKuyPFSu74qJSsqBiRdEFC6IsKoqCgEiVEgTpJSSBEEgyKdPv74+QIZGpmbkzd2bez/PwwMycmRy+DPPOOffccwVRFEUQERFRSFEEugNERETkewx4IiKiEMSAJyIiCkEMeCIiohDEgCciIgpBDHgiIqIQxIAnIiIKQQx4IiKiEMSAJwpCctufKhD9kVsNiOSGAU/kxIwZM9C+fXuHv1atWuX3Pq1btw7Tp0+33d66dSvat2+PrVu3+uXnt2/fHgsXLnTYH39YvHgx3n77bdvthQsXon379n7tA5HcqQLdASK5S01NRV5ent3HWrRo4efeAO+9916D2506dcLy5cuRnZ3t977Y648/zJs3D4899pjt9pgxY3DNNdf4vR9EcsaAJ3JBo9GgW7duge6GQzExMbLunz9kZGQgIyMj0N0gkhVO0RP5wLhx4zBu3LgG9/1x6vzzzz9Hx44dsWfPHtx1113o3LkzhgwZgiVLljR4XlVVFf71r39h0KBB6NatG2677Tb88MMPtp+zbds2bNu2zfba9qbo9+3bh/Hjx6NPnz7o0aMHJk2ahMOHD1/Rty1btuCBBx5A165d0b9/f8yZMwdms9mjv/cf+wMAZWVlmDVrFvr374/OnTvjzjvvxJYtWxo8t3379sjLy8Ptt9+Onj17YtGiRQCA7du3Y/z48bj66quRk5ODoUOHYuHChbBarbbnAUBeXp7tz/am6NesWYPbbrsN3bt3x4ABAzBr1iyUl5fbHl+4cCGGDx+O9evXY9SoUcjJycF1112HlStXuv33J5IzBjyRG8xm8xW/GrPIy2q14q9//StuvPFGvPnmm+jZsyfmzp2LjRs32h5/8MEHsXLlSkycOBGLFy9Gu3bt8Nhjj2Hr1q145pln0LFjR3Ts2BHLly9Hp06drvgZv/zyC+655x5YrVbMnj0bL774IoqKinD33Xfj6NGjDdr+7W9/Q8+ePfH6669j1KhReOedd/Dpp5+6/fex1x+DwYD77rsP69atw7Rp05CXl4eMjAw8+OCDV4T84sWLcd111+HVV1/FsGHDcPDgQdx///1ISEjAa6+9hsWLF6NHjx7Iy8vD//73PwDA8uXLAQB33HGH7c9/tGjRIkybNg1du3bFggUL8Oijj+Lbb7/FuHHjoNfrbe1KSkrw/PPP489//jPefPNNNGvWDDNmzLiiTkTBiFP0RC6cOXPGbpBOnToVjzzyiEevJYoiHnnkEYwZMwYA0LNnT3z33XdYv349rrnmGvz000/YtWsXFi1ahGHDhgEA+vbti5MnT+KXX37B1KlTERMTAwAOp+VfeeUVNG/eHG+99RaUSiUAYODAgRg+fDgWLlyIefPm2dqOGTMGjz76KACgX79++P7777F+/Xrcfffdbv19srOzr+jPihUrcPDgQaxYsQJdu3YFAAwaNAjjxo3D3Llz8dlnn9me36VLF0ycONF2+4svvkD//v3x73//GwpF7fhjwIABWL9+PbZv345Ro0bZfk5GRobdGpSXl2Px4sUYM2YMnnnmGdv97dq1w7333ovPP/8cf/rTnwAANTU1mD17Nvr16wcAyMrKQm5uLjZs2IA2bdq4VQMiuWLAE7mQmpqKxYsXX3F/enp6o16ve/futj9rNBokJSWhuroaALBjxw6o1Wrk5uba2giCgI8//tit166ursa+ffvw6KOP2sIdAOLi4mzB5agvQG1o1vWlsbZs2YLU1FR06tSpwXR/bm4uXn75ZZSXlyM+Ph5AbejWN3r0aIwePRoGgwGnTp3CyZMnsX//flgsFphMJrd+/u7du2E0GjFq1KgG9/fq1QtNmzbF1q1bbQEPNPyiVHcc39saEMkBA57IBY1Gg86dO/vs9SIiIhrcVigUtun+srIyJCQk2EavntLpdBBFESkpKVc8lpKSAp1O53ZfGqusrAwlJSV2Zz2A2mnxuoD/Yz/1ej1eeOEFrFq1CmazGc2aNUP37t2hUqnc7lfdcXZ3axAZGWn7c13deY49hQIGPJGPWCyWBrcbMwqMjY1FWVkZrFZrg5A/cOAAzGazyy8asbGxEAQB58+fv+KxkpISJCQkeNwnT8XGxiIrKwtz5861+3izZs0cPnf27Nn49ttvMW/ePPTv3x9RUVEAYJtCd0fdl4fz589fMc1eUlKC5s2bu/1aRMGMi+yIfCAmJgbFxcUN7tu1a5fHr9OrVy+YTKYGU+miKGLmzJm2wwTORvdRUVHIycnBmjVrGnzh0Ol0WL9+PXr27Olxn1z5Y3969+6NoqIiJCcno3PnzrZfW7ZsabAuwJ6dO3eiT58+uPbaa23hnp+fjwsXLthW0dv7mfV17doVGo0GX331VYP7d+zYgcLCQvTo0aMxf02ioMMRPJEP5Obm4ocffsDs2bNx7bXXYufOnfjiiy88fp0hQ4age/fu+Pvf/46pU6eiZcuW+Oqrr3Do0CE8/fTTAGqPp//666/YsmULOnbseMVrPPHEExg/fjwefPBBjB07FiaTCW+++SaMRmODzWF85Y/9ue222/Dhhx/iL3/5CyZNmoTMzEz8/PPPWLJkCcaOHQu1Wu3wtbp06YKvv/4aH3/8Mdq0aYODBw9i8eLFEAQBNTU1V/zM7du3o1evXg1eIyEhARMnTkReXh7UajWGDRuGgoICzJ8/H9nZ2bjtttt8XgMiOWLAE/nA7bffjlOnTmHlypVYvnw5evfujfnz5+Oee+7x6HWUSiWWLFmCV155BQsXLkR1dTWuuuoqvPXWW7YFcffeey/y8/MxYcIE/Otf/0JaWlqD1+jXrx/effddLFiwAI8//jg0Gg169eqFOXPmoG3btj77O9f5Y39GjRqFZcuW4ZVXXsG///1v6HQ6NG3aFE888QQeeOABp681Y8YMmEwmzJs3D0ajEc2aNcPDDz+MI0eO4IcffoDFYoFSqcSkSZOwaNEiTJgwAWvWrLnidSZPnoyUlBR8+OGH+OSTT5CQkIDrr78ef/3rXxsccycKZYLI1SREREQhh8fgiYiIQhADnoiIKAQx4ImIiEIQA56IiCgEMeCJiIhCEAOeiIgoBDHgiYiIQlDAN7opKdG5bkRERBTmFi6chxdemOX2xZA4giciIpK5unB//PGn3H4OA56IiEjG6of79Okz3X4eA56IiEim/hjugiC4/VwGPBERkQx5E+4AA56IiEh2vA13gAFPREQkK74Id4ABT0REJBu+CneAAU9ERCQLvgx3gAFPREQUcL4Od4ABT0REFFBShDvAgCciIgoYqcIdYMATEREFhJThDjDgiYiI/E7qcAcY8ERERH7lj3AHGPBERER+469wBxjwREREfuHPcAcY8ERERJLLy5vv13AHAJXkP4GIiCiM+XvkXocjeCIiIokEKtwBBjwREZEkAhnuAAOeiIjI5wId7gADnoiIyKfkEO4AA56IiMhn5BLuAAOeiIjIJ+QU7gADnoiIyGtyC3eAAU9EROQVOYY7wIAnIiJqNLmGO8CAJyIiahQ5hzvAgCciIvKY3MMdYMATERF5JBjCHWDAExERuS1Ywh1gwBMREbklmMIdYMATERG5FGzhDjDgiYiInArGcAcY8ERERA4Fa7gDDHgiIiK7gjncAQY8ERHRFYI93AEGPBERUQOhEO4AA56IiMgmVMIdYMATEREBCK1wBxjwREREQRHuJZUGfLSzwO32Kgn7QkREJHtyDneLVcSWExfwxd5ibDpWCosITL2+g1vPZcATEVHYkmu4F1fo8VX+WazKL8ZZnaFRr8GAJyKisCS3cDdbRWw+dgFf7CvCz8cvwCp693oMeCIiCjtyCveiCj1W7SvGl/nFKKk0+ux1GfBERBRW5BDutaP1UqzcW4yfj1+Al4N1uxjwREQUNgId7sUVenzRyNF6WowGt3TOcLs9A56IiMJCoMK97tj6yr1FHo/WlQIwsHUybu2Sib5ZiVAq3O8zA56IiEJeIML9rM6AVfuKsGpfMc55OFrPjNNidOdMjMpJR2qMtlE/nwFPREQhzZ/hbrGK+OXkRXy+pwibjpV6tBJeKQCDslNwa5cM9GmZCIWX/WTAExFRyPJXuJ+vMuKr/GKs3FuEogrPzltvEqfF6C6ZGNUpHSmNHK3bw4AnIqKQJHW4W0URO06V4fO9RVh/pBQWD4brSgG4pk0ybuua6ZPRuj0MeCIiCjlShntZjQmr95/Fyr1FOHWxxqPn1h1bvznHt6N1exjwREQUUqQId1EUsbewAp/vLcL3v5fAaHF/tK4QgGtaJ+PWrpno29KzlfDeYMATEVHIyMub79NwrzKa8c2Bc/hsTxEOl1R59Ny0GE3taL1zBtJjpR2t28OAJyKikODLkfuhc5X4fG8Rvv7tHKpNFrefJwDo3yoJt3bJxIDWSVD5abRuDwOeiIiCni/C3WC2Yt2hEny2pwh7Cys8em5SlBq3dM7A6M6ZaBIf4fHPlgIDnoiIgpq34V5QVoOVe4vwZf5ZlNWYPHpurxYJuL1LJgZnJ0OtVHj03Dq/FevwwfbT2FekQ1mNCRmxWtzYMR1jezXDHe9uR1GFAase7O3xFwcGPBERBa3GhrvFKuLn4xfw6Z5CbDl+0aPtY+MiVLipUzpu7ZKJrKSoxnX8kq/yizF77SFYxNrXbZ0chcJyPRZvPoHNxy/AYLY2+rUZ8EREFJQaE+4Xqo1Yta9xG9J0zozFbV0zcW27VESolY3tts2x0ir887vDsIjA/b2bY2L/llArFTBbRXyw/TQWbzrh1VXmGPBERBR0PAl3URSxr0iHT3YXYt2hEpg8OMUtUq3A9R3ScHuXJmifHuOLrtss+fkUzFYRw9ql4NFrWtnuVykE/KVPC5zVGfDZnqJGvz4DnoiIgoq74V5jsuDbA+fwye5CHPLwFLfWyVG4vWsT3NgxDTFa30elyWLFz8cvAABu75ppt809PZoy4ImIKDy4E+4nL1Tjsz1F+Gp/MSoN7p/iplIIGNo2Bbd3y0T3pvGSXpSmsFxvO/2ubar9mYGWSVGI1ihRZXT/71AfA56IiIKCs3C3WEVsOlaKT3YXYuvJMo9eNyNWi9u6ZuLmnAwkR2t83Gv76q/Wj3JyPJ8BT0REIc1RuJdVm/DFviJ83ohFc32zEjGmWxMMaJXkt+1j60TWC/Uqoxkalf0vFo0Nd4ABT0REMmcv3H8r1mHF7kJ8d/CcR/vCx0WoMKpTBu7ololmCZES9tq5lklRUCsFmCwiDpdUoXfLKwO+sFzPgCciouAlCAIEARDF2hXv9dUP978+8Xd8feAcVvxaiP3FOo9+Rof0GNzRrQlGtPfNKW7e0qoU6J+VhA1HS/FlfjF6t0y8os3nexu/wA5gwBMRUQAIAhARoUZkpAbKejvAWSxW1NQYodebsGBBbbhPevxpxA64Azcv2YaLHuw0p1EKGN4+FWO6NUHHjFhJF801xvh+LbDpWCm+PViCrKQo3N+7OVRKBURRxOd7i/DhjgKvXl8Q//h1yc9KSjz7FkZERMFNrVYiPv7y9Hj94K2LpDlzXsazi5ahx92Po1iZCqsHSdUkTovbuzbBzTkZSIhS+6zfUvhkdyH+ve4IRADxESo0S4hEUYUeF6pNyMmMRX5RbUauntjHdkW61NRYt15bsoAXRRGCIGDhwoWYPHmy7fYfMeCJiMJH/XC3lwnVRjMefuldfHfSCE1qlkevHchFc974taAcS7efRn6RDlVGM5rGR2BUpwwMvyoVNy/ZBgD4/pF+iI+s/bLibsBLNkX/0UcfQaVSobq6GnPnzsWgQYPQu3dvqX4cERHJnCDAYbifLK3CB1tO4oPNR2AQm0CT6t5rRmuUGJWTgTu6ZqKll/vCB0r3ZvHo3iz+ivuPnq/dnEejFBAX4XlcSxbw9957L9566y188MEHeP/999GzZ0+pfhQREQWBiIjaEWhduIuiiI2Hz+O9n0/gx9/PoXY+2b0FcG1SonBntya4vkM6ojSBXzTXGJM/24fyGhMmDchC/1ZJVzxet9Ndu7SYRq0fkCzgP/74Y8TGxmLcuHFYv349RFFEr169pPpxREQkc5GRtaeCVRrM+HxXAd7/+QSOllQiEgaI0AJwHmJKARjSNgV3dm8i+U5z/tA6OQof7TyDeRuOoUVipO20PVEUseFIKZZsOQkAGNOtSaNeP+DH4J0xGMyoqKix3U5Jcfwtxmg0o7z8ctvk5GgoFPavzWsyWVBWVm27nZQU3WAVZ31mswUXL15um5gYBZXK/rdFi8WKCxcu73eckBAFtYPTMaxWK0pLL7eNj4+ERmP/+5Yoijh/vtJ2Oy4uEloneyPXX9cQFxcBrdbxIpP6bWNjI2zfsO05f77StgAmJkZr+89qT2lpJayXVsVER2sRFeW47YULVbBYai+JGBWlQXS01mHbixerYL50+cTISDViYhxfH7msrBqmS1tBRkSoERvruG15eTWMl8431WpViItzfH5seXkNjEYzAECjUTVYLPRHFRU1MBjq2ioRH+94ClGn00Ovr10hrFYrkZDguG1lpR41l1YTq1QKJCZGO2xbVWVAdbURAKBUKpCU5LhtdbURVVW1m4UoFAKSkx1fXKOmxojKytq2giAgJcVxW73eBJ1Ob7vt7BiiwWBCRYW7bfkZAQTPZ8TJ0iq8//NJfLLjNHQGE6JgwCea59BJcRLbre0wxvgM7IV8crQGoztn4NYumbaFZqHgYrUR4z/ejdNleigEoFlCJKI1SpzVGXChuvb/913dm+BvQ7MbPE+yY/CrV6/G3//+d3z//fdIT08HAMycORN79+7FRx99hNjY2h9c959s8uTJDW4TEVH4EEURm46cx3ubT+AH2zS8iE81z6GX4pCtXRSu3IWua7N43Nc/CyO7ZKKyvMY2aAgViVEafDiuJ77ML8bag+dwplyP4go9EqM0GNo2Bbd2yUDfrCun7t3l8QheFEXccsst6NWrF2bNmoW8vDx88sknWLFihS3wPcFV9EREoafaaMGa385ixa+FOH6hboZDRCQMiIIBOyMetrXdb22Jm4yzIUIBtVLATV2a4L7+WejWPMHW5vx5HQJ7Urd8SDaCFwQB06ZNw5QpU5CWloalS5di2bJljQp3IiIKLYXleqz4tRCr8ov+cCW3K0ftANBTvxiliENabATG9m2Je3q3QGq9aXhRFGG1igz3RmjUIrvc3FxkZ2cjLy8P77zzDtq2bevrfhERUZAQRRG7Csrx311n8NPRUrub0kTCcEW4b7e2Q8vmLfDMwNa4IScDagfrHGpqjFJ0O+Q1KuA3btyIY8eOwWq1Ijk52dd9IiKiIGAwW/HtwXNYvusMDpVUoW4K3p76x9h7VC1EhroG/5o0Bp+3cHyMue4Ict3iU/KMx8fg9+/fj3HjxmHWrFlYs2YNIiIisGDBgkZ3gMfgiYiCy/lKAz7dU4TP9xTV2xve/hS8Pf9nnoIXn3/O4VkMwOVwLy+vsZ0NQ7UkOQZfUFCAhx56CBMnTsTo0aORnZ2NO+64A/n5+cjJyWlUR4mIKDgcPKvDx7vOYO3BEpj/MA9vbwrenuPmDDz5+HQIgmALcXt70QMMd285/vr0B2VlZZgwYQKGDh2KSZMmAQBycnKQm5uL1157TbIOEhFR4FisIn48fB4Tl+/BuA9/xZrfzl0R7n/UU78YOcZ3MbnFl/jp5h14wToF0f+swN914xE9eRvMFhGlpZWoqjJcceqb1SqiqsqA0tJKhruXeDU5IqJwZdZDe3Q1tMe+haAvgxiRAEPr62BocxMqLSp8mV+M5bvOoLDC/nH1+iKhx4GIBwAA/+r8HW7p0RqpMdoG13OfPn2m3T1R6qaceSqcewJ+sRkiIpIvzfG1iF03DQpDOUQoIMAKEQpoj32N6h/+gWfNj+B/xm5uv152SjRwaTO9B/u1BNTuhXt9DHffYsATEYUZzfG1iFsz3nZbgLXB7xGWSiwU/g2D4nF8b3V8oTABwODsZNzdoyl6pKuBJZcf8zTcyfcY8ERE4cSsR+y6aQAAAfaHzAoBsIoi5qpfRx/Df2BAwz3lozVK3JyTgTu7N7FdIAWmy/vxv/56Hl544UW3w53nuUuDAU9EFEa0R1dDYSh32U4hAAmowi2KzfjK2g8AkBkXgTu6NcGNHdIQrVUBEG3BLtQL+Jdecj/cAdguWkS+xUV2RERh4vTFGqi/HI+Ouo1QCtJ99P9dNx7Tpj/LaXmJcJEdERFBFEXsOVOBZTsLsOFIKT5Sl0CplC7cj5szMO2pZzwK97q2AR5vhhwGPBFRCDJbRfxwqATLdp7Bb8WXZ0rLEAOLKLg1ghehgDHrWlSMyHPa7vXX8/DSSy9iyuRpteHuZIc6e1JSYgBwRtfXGPBERCGk0mDGl/nF+O+uMyiyc/76Wksv3KDc7tZrCbDCkD0SUEc5bJOXNx/PX1pQN42r5WWFAU9EFAKKK/T4765CfLGvCFVGxzvArbH2wTPiUsShCgonWSxCgKiNg6HNSIdteCqcvDHgiYiC2MGzOny4owDf/14CixuHsA3QYLrlYSxWvQIR9k+VE1Eb1Lph8wBVhN3XYbjLHwOeiCjIiKKIn49fxIc7TmPHadenvNWJj1Dhjm5NMKZbX+jOdUTsumkQ6p0yZ9vRThsH3bB5MLYabvd1GO7BgafJEREFCaPZim8OnsOHOwpwvLTa9RMuaZEYiT/1bIqRHdMRoVZefsCsh/bQSsT9+CQAwJA1AobskbXT8n4cuded9sU8cA9PkyMiChEVehM+21OE5b8WorTK/V3fejSLx729mmFg6yQo7AWxKgKGtrcAlwK+YkSe0wV1HLkHFwY8EZFMFVXo8dHOM1i1rwg1Jqtbz1EKwLXtU/Gnns3QMcO9kZ47pAx3vd7ks9eiyxjwREQy8/u5Snyw/bTbC+eA2v3hR3fOxN09miAjzv70emNJPXLX6fQ+fT2qxYAnIpIBURSx7WQZPthxGltPlrn9vLQYDe7u0RS3dslEjNb3H+mclg9eDHgiogAyW0Ws+70ES7efxqGSKref1zY1GmN7NcPw9qlQKz3bOc5dDPfgxoAnIpKKWQ/t0dXQHvsWgr4MYkQCDK2vg6HNTagR1Vi1rxgf7Sywu+OcI31bJmJsr2bo3TLB68A1WgxYX/gddqSloEyhQOSeZzEgMxdDMobijUWL/RbuXEUvDZ4mR0QkAc3xtYhdNw0KQ/nl88sv/V6jjMV0y8P4Ut/NrddSKgSMaJ+Ksb2aoV1ajE/6t/nsRszZ8yIqzTooRBFWQYAAASJEqC1qHFl4BBOGPuyXkTsD3jM8TY6IKEA0x9cibs14220B1ga/a806zMO/Ua14HN9bezp8nSi1EqO7ZOCeHk19unBu89mNmLVzhu22te5qbpd2tTMKRrSc2hLX9BzEafkg1ugDN6tXr0bnzp1x9uxZ230zZ87EqFGjoNPxWxgRhSmzHrHrpgGwvw0sgEt7wIuYq34dWlx5XntytAaPDszC6ol9MG1IG5+Gu9FiwJw9L6K2B/b7JygECIKAOXtnw2hx//AByUujA37kyJFo1aoV3njjDQBAXl4eNm3ahLfeeguxsb4795KIKJhoj66GwlDuMNzrKAQgQajCLYrNiIQekdCjfYKAZ4c1w1f35eD+HimIVRoBU7VPf60v+AaVZp3DcK8jQkSlWYcNxT/6sjzkR42eohcEAdOmTcOUKVOQlpaGpUuXYtmyZUhPT0dlZSXuu+8+mM1mWK1WjBs3Dnfeeafd13F2LMFgMKOiosZ2OyUlxuF0kdFoRnn55bbJydFQOLgmsclkQVnZ5W0ek5KioXSwCtVstuDixcttExOjoFIp7ba1WKy4cOHyKtiEhCio1fbbWq1WlJZebhsfHwmNxv4/hyiKOH++0nY7Li4SWienw9Q/jhUXFwGtVu1W29jYCEREOG57/nwl6pZsxMRoERmpcdi2tLQSVmtt2+hoLaKiHLe9cKEKFkvt1GVUlAbR0VqHbS9erILZXNs2MlKNmBjHI5uysmqYTLVX1YqIUCM21nHb8vJqGC9dgUurVSEuLtJJ2xoYjWYAgEajQny847YVFTUwGOraKhEf73iXMJ1Ob9vwQ61WIiHBcdvKSj1qamrbqlQKJCZGO2xbVWVAdXXtKFGpVCApyXHb6mojqqpqR2wKhYDkZMfHe2tqjKisrG0rCILtmt726PWmBuc6O/9/b0JFhbtt5fUZYRVFVP/2P0RDgNJFgNZ5WbMEL2NJ7Q09gM2XfklkR1oKFFGRtml5ZwQosKn4Jwxver10HSLJeHUMPjc3F9nZ2cjLy8M777yDtm3bAgAiIyPx4YcfIjIyEjU1NbjpppswfPhwJCYm+qTTRERyYrJY8c2Bc1i6/TT+pStAhjKga5edKlMo3Ap3ABBhRYXJ/YvZkLx4FfAbN27EsWPHYLVakZycbLtfqVQiMrJ2VGMwGGC1WuFosb4nqybrj2JdqT86dqX+qNuV+qN5V+qPAFypP7Jwpf6IxXVbPWqHBa7pdHq3d5SqrDTYRm+uVFUZbKNCV6qrjbbRpis1NSbbKNYVvd7k9naYBoPZ7fel0ehJW4vbbU0m99uazVa321os7re1WkW324qi+20Bz/7fy/kzosZkwRf7irFsRwHO6mrf42XqGFhEAUrBdciLUMCYdW3tHvB+ErnnWQjnNrucogdqR/Bx6njJ+2QwcKtaKTQ64Pfv34+pU6fiueeew5o1azB//nwsWLDA9nhFRQXGjh2LkydP4qmnnkJSUpJPOkxEFGhlNSZ88mshlv96BuV6c4PH1lp64QbldrdeR4AVhuyRTi/w4msDMnOx8dwmt9qKsGJgxiCJe4QGh2TIdxq1yK6goAAPPfQQJk6ciNGjR2PKlClYu3Yt8vPzbW3i4uLw5ZdfYt26dfjqq69w/vx5n3WaiCgQzukMeG39Udy8ZCve3HLyinAHgDXWPigTo2F1MUAWIcCqja+9NKsfDckYCrVFDdFFBwUIiFHFYnBGrp96Rr7mccCXlZVhwoQJGDp0KCZNmgQAyMnJQW5uLl577bUr2qekpKB9+/bYvt29b7RERHJz+mINZq89hFve2oaPdp5xemU3AzR4wjQJEASIsH+su+5+3bB5Dq+7LpU3Fi3GkYVHIFza2MaeuvtndH0aGqXjRa8kb5LsZHf+/HlEREQgJiYGlZWVuOuuu/DKK6/gqquuuqItdy4iIrk6UlKF97adwne/l7gckddpkxKF+3o3x02a3Uj48XEoDJcXqdXtZGfVxkM3bB6MrYZL1HP76u8tf839gzBn72xUmnUQoIAIq+33GFUsZnR9Gv3TB/qlX9zJzjMB3cmuuLgYM2fOhCiKEEUR9957r91wJyKSo32FFXh36ylsPHbB7ed0aRKH+3s3x8DWSRAEARZch9KWO6E9tBJxPz4JADBmXQtD9sjaaXk/j9ztXTjm09Q+2FD8IzYV/4QKUzni1PEYmDEIgzNyOXIPAdyLnogItWcCbD9Vhne3ncaOU2VuP29AqyTc17s5ujdzsNrcVI3UN9sBAEomHvLrgro6cr8qHEfwnuFe9EREbhBFERuPXcA7v5zC/mL3AkYhANe2S8V9vZv77OIvUpF7uJN0GPBEFJYsVhHrDpXgvW2ncdjN67CrFAJu6pSOP1/dHM0THe9gKBcM9/DGgCeisGK2WPHNwXN4d+tpnLro3qZRESoFbuuaiT/1bIb02OA4Ns1wJwY8EYUFg9mK1fuLsXTbaRRWuLezYqxWhTu7N8Hd3ZsiIcrxdRrkhuFOAAOeiEKc3mTB53uL8OGOApRUurcNclKUGvf2bIbbumYixsmFneQoGMO97oJM5FvB9c4lInJTldGMT3cXYdmOAlx085oF6bFa/PnqZrg5JwMRDq4EKWfBGO6AZ9fXIPcx4IkopFQazFj+6xl8vPPKfeIdaZ4Qgft6N8eNHdOhdnBZWLnLy5sflOFO0mHAE1FIKK8x4b+7zuC/v55BpcHi1nNaJ0fhL31a4Nr2qVApgjcQg3XkTtJiwBNRULtYbcRHO8/gk92FqDK6F+wd0mPwlz4tMDg7GYogD8NQCPeUlNq9BDy53C+5xoAnoqBUWmXEhzsK8OnuQujNji/+Ul+XJnEY37cF+mUlBmUQ/lEohDuAoO233DHgiSionK8y4oPtp/HZniIY3Az2ns3j8WDflujZPD5kwiRUwp2kw4AnoqBwvtKApdsL8Ple94O9b8tEjO/bAt0c7RMfpBju5A4GPBHJWkmlAe9vO40v9hW7HewDWydhfN8WyMmMk7h3/sdwJ3cx4IlIls7p6oK9CEaLexe9HJKdjAf6tkCHdPeuthVsGO7kCQY8EclKSaUB7231LNiHtUvB+L4t0DZV3ld28wbDnTzFgCciWThfacB7205j5V73gl0AMKxdKsb3a4HslGjpOxhAoR7uRiO3qpUCA56IAup8pQHvby/ASjcXzwkAhrdPxQN9W6BNiAc7EPrhDgDl5dyqVgoMeCIKiNIqI5Z6cLqbAGDEVakY37clWiVHSd9BGQiHcCfpMOCJyK8uVBuxdFsBPt1T6FawKwRgxFVpGN+nBbLCJNgBhjt5jwFPRH5RVmPCB9sLsOLXM27tPCcAuK5DGsb3bYGspPAJdiD8wj05ufZQS2lpVYB7EloY8EQkKZ3ejGU7C/DfXWfc2iu+bir+wb4tw2rEXifcwh0AFIrgvIKf3DHgiUgSVUYz/rvrDJbtOAOdwfUqaQHAte1T8WC/FmidHPqL5+wJx3An6TDgicinakwWfPJrIZZuP+329divbZeC8f1ahvzpbs4w3MnXGPBE5BMGsxWf7SnE+9tO40K1ya3n5LZNwcR+LZGdGr7BDjDcSRoMeCLyitlixZf5xXj7l1M4V2l06zmD2iRjYv+WaJ8WujvPuYvhTlJhwBNRo1isIr45cA5LtpzEmXK9W8/pl5WIhwZkoVNGaO4V7ymGO0mJAU9E9pn10B5dDe2xbyHoyyBGJMDQ+jrUtB6JH49X4o3NJ3H8QrVbL3V1iwQ81L8lujYNrcu2umK0GLC+8DvsSEtBmUKByD3PYkBmLoZkDMUbixYz3C8xmVyfXUGeE0RRdO9qDhIpKdEF8scTkR2a42sRu24aFIZyiFBAgNX2uw7R+KtxEtZZe7p8ne5N4/DQgCz0bJ4gfadlZvPZjZiz50VUmnVQiCKsggABAkSIUFvUOLLwCCYMfTjsw508l5rq3gwYA56IGtAcX4u4NeMBAAKu/HiwirWPTDQ9ju8dhHynjFg8PCALvVsmhGV4bT67EbN2zgAAiHZqKFpFCIKAF3q+hAEZg/zdPQpy7gZ8o3cXWL16NTp37oyzZ8/a7ps5cyZGjRoFnY6hTRSUzHrErpsGwH64A7VbxwIi5qpfhxYNF9W1TY3GK6M74d0/dUOfrMSwDHejxYA5e14EYD/cAUBQCBAEAXP2zobRYvBn9yiMNPoY/MiRI/Hmm2/ijTfewKxZs5CXl4dNmzZhxYoViI3lAhqiYKQ9uhoKQ7nLdgoBSEAVblFsxlfWfmieEIkH+rZAbtsUKAQBMIfv1cHWF36HSrPrQY4IEZVmHTYU/4jhTa/3Q8/kKymp9jTJCxe4Va0vNTrgBUHAtGnTMGXKFKSlpWHp0qVYtmwZ0tPTUVRUhKeeegqlpaVQKpV45JFHcMMNN9h9HWdTDQaDGRUVlz8oUlJiHI4IjEZzg0sOJidHO9z+0GSyoKzs8uKgpKRoKJX225rNFly8eLltYmIUVCql3bYWi7XBGzQhIQpqtf22Vqu1wb7L8fGR0Gjs/3OIoojz5yttt+PiIqHVOv6nq3/YIy4uAlqt2q22sbERiIhw3Pb8+UrUHdGJidEiMlLjsG1paSWstXO5iI7WIirKcdsLF6pgsdTuTR4VpUF0tNZh24sXq2C+tI95ZKQaMTERDtuWlVXbFu9ERKgRG+u4bXl5NYyXtlHValWIi4t00rbGdv1qjUaF+HjHbSsqamAw1LVVIj7e8darOp0een3t+eNqtRIJCY7bVlbqUVNT21alUiAx0fF55FVVBlRX1460lUqF7cPU7uuu+RZWCFA4GHn+0cuaJXgZSwA9gPWXfoW5HWkpUERFwurG7IUABTYV/xT2Ae/o85e849Uq+tzcXGRnZyMvLw/vvPMO2rZtCwBQKpX4v//7P3To0AGlpaW49dZbMXjwYERFhd++0kTBoLTSgP/8eBQjjp5AX0VAl+UEvTKFwq1wBwARVlSYXM+YEDWGVwG/ceNGHDt2DFarFcnJybb709LSkJaWBgBITk5GfHw8ysvL7Qa8J4vs6o9iXfHkqkSeTAvVH827Un+WwJX6sw+u1J/VcN1Wj9rhlWs6nR46nXttKysNqKx079hhVZUBVVXuta2uNtpGm67U1Jhso1hX9HqTbXTsisFgdvt9aTR60tbidluTyf22ZrPV7bYWS8O2lQYzlu0owEc7z6DaZMHV6hhYRAFKwXXIi1DAmHUtKkbkufWzw0XknmchnNvs8Ph7fQIUiFOH16mD5D+NDvj9+/dj6tSpeO6557BmzRrMnz8fCxYsuKLdvn37IIoiMjMzveooEfmO0WzFZ3uL8M4vp1BW70vSWksv3KDc7tZrCLDCkD0SUHNmrr4BmbnYeG6TW21FWDGQq+hJIo068FFQUICHHnoIEydOxOjRozFlyhSsXbsW+fn5DdpdvHgR06dPx/PPP++TzhKRdyxWEWt+O4sx727Hqz8ebRDuALDG2gdlYvSlU+EcEyHAqo2Hoc1ICXsbnIZkDIXaoobooogCBMSoYjE4I9dPPaNw43HAl5WVYcKECRg6dCgmTZoEAMjJyUFubi5ee+01Wzuj0YjHHnsMEydORI8ePXzXYyLymCiK2HSsFGM/2IVnvv4dhRX2D5kYoMETpkmoHZ/bP44sXrpfN2weoHK8cDFcvbFoMY4sPALh0sY29tTdP6Pr09AoHS8qJfKGJBvdiKKIJ554Aq1atcLkyZOdtuVGN0TS2ltYgbyfjuHXMxVutR/aNgUzso4ie9uMBqfM1e1kZ9XGQzdsHoythkvV5aBVf2/5a+4fhDl7Z9fuZAcFrLbzE6yIUcViRten0T99YKC7LAuJibWHeTxZ4xTOArqT3Y4dOzB27Fi0b9/edt/LL7/c4HYdBjyRNI6XVmPRpuNYf6TUrfa9WyTg0WtaoWPdhWDMemgPrUTcj08CAAxZI2DIHlk7Lc+R+xXsXTjGaDFgQ/GP2FT8EypM5YhTx2NgxiAMzsjlyJ0ajVvVEoWpkkoD3vz5JL7ML3Z5LB0AOqTH4NFrWqFPy8QrHzRVI/XNdrWvO/EQF9Q5wKvCkT+5G/C8mhxRiKg0mPHBjgJ8tKMA+kubATnTPCECDw9shWHtLu0+R43CcCe5YsATBTmTxYqVe4vw1pZTuOjGvgDJ0RpM6NcCt+RkQMUdxLzibrjzGLNzrI80GPBEQUoURaw7dB7/2XQcBWWuNyiK1ihxX+/muLtHU0Q62EKZ3OfJyN3R9tZUi/WRBgOeKAjtKijDgg3Hsb/Y9RoWtVLAmG5N8Jc+LZAQ6fhaA+Q+TstTMGDAEwWRExeqkffTcWw46t7K+Bs6pGHSgCw0ieeqd19huFOwYMATBYGL1UYs2XIKn+8phMWNlfG9WyRg8qBWuCqdl272JYY7BRMGPJGM6U0W/HfXGby37TSqLl3O1pm2qdGYMqgV+mYl+aF34YXhTsGGAU8kQ1ZRxDcHzmHRphM4q3N9Jb70WC0eGZiF6zuk8ZQ3CTDcKRgx4IlkZsepMszfcAwHz7m+PHKsVoW/9GmOO7s3hVbFU96k4Itwt1hc70sQzlgfaTDgiWTixIVqLNhwDBuPXXDZVqWoXRn/QF+ujJeSr0buFy5U+bhnoYX1kQYDnijAympMeGvLSXy6pwgWN/aWHdYuBY8ObIXmiZF+6F344rQ8BTsGPFGAmCxWfLK7EG//cgoVerPL9p0z4zB1cCt0bRrvh96FN4Y7hQIGPJGfiaKIDUdKseCnYzjtxg50TeMjMHlQKwxtm8Kg8QMpwj0hoXYr1rIybsVqD+sjDQY8kR/9frYSr204ip2ny122jYtQYXzfFhjTrQnU3DPeL6Qauau5NbBTrI80GPBEfnC+yohFG49j9f6zcHWUXakQcFf3JnigTwvEcwGd31wO9yc5LU8hgQFPJCGD2YqPdxbg3a2nUW1yvVHN4DbJmDK4NVpwAZ1f8Zg7hSIGPJEERFHEj0dKMX/DMRSWuz7O3i41GtOGtEGvFgnSd44aYLhTqGLAE/nYoXOVeHW9e8fZU6I1eHhgFkZ2TIdSwWDxN4Y7hTIGPJGPXKw2YvHmE1i1rxiuTmfXqhQY16sZxl3dHFEaLjAKBIY7hToGPJGXTBYrVvxaiLd+OYlKg+vj7NddlYrHrmmFjDhewjVQ/B3uViu3YnWG9ZEGA57IC1tOXMArPxzFyYs1Ltt2zIjFE7lt0KVJnB96Ro4EYuReWsqtWJ1hfaTBgCdqhIKyGry2/hh+Olrqsm1KtAaPXdMKN3Tkld4CjdPyFE4Y8EQeqDZa8N62U/hwRwFMFucH2jVKAWN7NcN9vVvwOLsMMNwp3DDgidwgiiLWHizBgp+O4Vyl0WX7a9ul4LFBrdA0nuezy0Ggwz3+0vugvNz1oZxwxPpIgwFP5MLv5yox94cj2H2mwmXbtqnReCK3DXo2T5C+Y+SWQIc7AGg0/Kh1hvWRBqtK5EB5jQmLN5/Ayr1FLk97i49Q4eGBWRjdOZPns8uIHMKdKFAY8ER/YBVFrNpXjP9sPI5yF5dxVQjA7V2b4KH+LblvvMww3CncMeCJ6tlfrMPL647gt2Kdy7Y9msXjb0PboG1qjB96Rp5guBMx4IkAAGU1JizadBxf7C12ebW39Fgtpg5ujWvb8frscsRwJ6rFgKewZrGKWLWvCIs2nXA5Ha9RChh3dXPc17s5Inn9alliuBNdxoCnsJVfVIGX1x3BgbOVLtte0zoJj+e2QbMEnvYmV3IOd1F0NS8U3lgfaTDgKeyU15jwHzen45vGR+BvQ9tgYOtkv/SNGkfO4Q4A58+7/hIZzlgfaTDgKWxYRRGr95/Fwp+Oo6zG5LStVqXA/b2bY9zVzaFVKfzUQ2qMvLz5sg53okBhwFNYOFJShTnrDru1Wc2Q7GRMG9IGTeJ5tTe5k/vInSiQGPAU0qqMZiz5+RT+u6sALraOR/OECPxtaDb6t0ryT+fIK8EU7nFxtWs3Kiq4Fas9rI80GPAUkkRRxI+Hz+OVH4+63Dteq1JgfN8WuLdnM2g4HR8UgincAUCr5UetM6yPNFhVCjkFZTV4ed0RbDlx0WXbQW2S8UQup+ODSbCFO1GgMOApZJgsVny4owBv/3IKBrPVadvMOC2eyM3G4Gyujg8mDHci9zHgKSTsLijHP78/jOOl1U7bqRS112h/oG8LblYTZBjuRJ5hwFNQK68xYeHG41i1r9hl217N4/HUsLZolRzlh56RLzHciTzHgKegJIoivj5wDvPWH8NFF+e0J0WpMW1IG1x3VSqDIQgx3IkahwFPQefUxRq89P1hbD9V5rSdAOD2rpl4ZGArxEbwrR6MGO5EjcdPPQoaJosVS7efxju/nILRxUntbVOjMXN4W3TKjPNT78jXQincS0pcX344nLE+0mDAU1DYW1iB2WsP4ZiLRXQRKgUeGpCFu3s0hUoRvIEQ7kIp3IkChQFPslZlNGPRxhP4ZHehywvDDGydhKeGZSMzjue0BzOGO5FvMOBJtjYeLcVL3x92uRNdaowGfxuajdzsZIZBkAvVcI+79KWzokIf4J7IE+sjDQY8yc75KiNe+eEovj9U4rSdAODO7k0waUAWYrjVZdAL1XAHAK1WfelPDDB7WB9p8FORZEMURXyVfxbzNhyDzmB22rZtajRmjmiHThmxfuodSSmUw50oUBjwJAunL9bgn98dwo7T5U7baVUKTOjXEvf2bAqVkheGCQUMdyJpMOApoCxWER/vOoPXN59wuX98rxYJ+L9r26J5YqSfekdSY7gTSYcBTwFz9HwVXlx7CPlFzs+BjYtQYerg1hjVKZ0BEEIY7kTSYsCT35ktVry37TTe/uUUzFbnJ7+NaJ+Kx3PbIDla46fekT8w3Imkx4AnvzpwVocXvj2EwyVVTtulx2ox49psDGzNy7mGGoY7kX8w4MkvDGYrlmw5iQ+3n4aLXWZxZ7cmeOSaLERr+PYMNeEa7tyK1TnWRxr8BCXJ7TlTjhe+PYSTF2uctmuRGIl/jGiH7s3i/dQz8qdwDXeiQGHAk2T0JgsWbz6Bj3eecbrNrEIAxvZqhgn9WiJCrfRb/8h/GO5E/seAJ0nsLazAc9/8jlMuRu3ZKdF4+rp26MgNa0IWwx2Ija3dilWn405t9rA+0mDAk08ZzFa8sfkElu0sgLMF8iqFgAf6tMD9fZpDzQ1rQhbDvVZERO1WrAww+1gfaTDgyWfyi2pH7ScuOB+1d8yIxdPXtUN2SrSfekaBwHAnCiwGPHnNYLbizZ9P4sMdp52O2jVKAQ/1z8KfejXjtdpDHMOdKPAY8OSV/cU6PPfN7zheWu20XaeMWDxzfXu0So7yU88oUPLy5jPciWSAAU+NYrJY8fYvp/De1lNOz2tXKwVM7NcSY69uzlF7GODInUg+GPDksaPnq/Ds17/j4LlKp+06pMfgmevbow2PtYcFhjuRvDDgyW1WUcTHO89g0abjMDoZtqsUAib2b4lxHLWHDYY7kfww4MktheV6PPfN79hV4Px67Vel1Y7as1M5ag8XDHfXzp93PtsV7lgfaTDgySlRFPFV/lm88uNRVJssDtspFQIe7NsC9/duDhXPaw8bDHf3iKKLCzCEOdZHGgx4cqi0yojZaw9h47ELTtu1To7C8zdchfbpMX7qGckBw51I3hjw4c6sh/boamiPfQtBXwYxIgGG1tdhrdgPL/5wCmU1JodPFQDc26sZJg3IglbFUXuoMVoMWF/4HXakpaBMoUDknmcxIDMXQzKG4o1FixnuHoiJ0QIAKisNAe6JPLE+0hDEAM+N8DKBgaM5vhax66ZBYSiHCAUEWG2/l4nReMI0CeusPe0+t0mcFs/c0B49miX4t9PkF5vPbsScPS+i0qyDQhRhFQQIECBChNqixpGFRzBh6MMMdzelptZea4Gfd/axPp6pq5crHHaFKc3xtYhbMx6CoQIAIMDa4Pc4VGGJ+lVcq9h5xXNv6ZyBj+7ryXAPUZvPbsSsnTNQZa5d+GS9FODipWsCGgUjWk5tiWvuH8RwJ5KxRgf86tWr0blzZ5w9e9Z238yZMzFq1CjodPwWJmtmPWLXTQMACA4u5Fp7dpuIuerXoYURAJAUpcarozvhHyPaIVrDozuhyGgxYM6eFwFcDvQ/EhQCBEHAnL2zYbRwSpVIrhr9KT1y5Ei8+eabeOONNzBr1izk5eVh06ZNWLFiBWJjeelPOdMeXQ2FwfnpbkBtyCegCrcoNqOi1U34W2424qPUgMn5trQUvNYXfodKs+sv6CJEVJp12FD8I4Y3vd4PPSMiTzU64AVBwLRp0zBlyhSkpaVh6dKlWLZsGdLT0wEAjz76KLZt24Z+/fphwYIFDl/H2bEEg8GMiorLVyZLSYlxOCVoNJpRXn65bXJyNBQK+xMUJpMFZWWXQyopKRpKB6d2mc0WXLx4uW1iYhRUKqXdthaLFRcuVNluJyREQa2239ZqtaK09HLb+PhIaByMikVRbHCeaFxcJLRax/909Y9jxcVFQKtVN2zwwzqIggKCaHX4GvW9rFkCnFkCfOhWcwpiO9JSoIiKtE3LOyNAgU3FPzHgiWTKq2Pwubm5yM7ORl5eHvLy8tC2bVvbY+PGjcOcOXO87iD5nqmy1O1wp/BSplC4Fe4AIMKKCpPrmSAiCgyvDqRu3LgRx44dg9VqRXJycoPH+vbti61bt7p8DU9WTXqy21H90bEr9UfdrtQfzbtSf5bAlfqzD67Un9Vw3VYPQG+7vfN0GTQFFgwWBSgF1ydQiFDAmHUtKkbkuf0zKXhF7nkWwrnNDo+/1ydAgTh1vB96RUSN0eiA379/P6ZOnYrnnnsOa9aswfz5851OxVNgma0i3tpyEu/8cgqjFT0wVOP6yxdQu6rekD0SUPMyr+FgQGYuNp7b5FZbEVYMzBgkcY9CQ2kpt2J1hvWRRqOm6AsKCvDQQw9h4sSJGD16NKZMmYK1a9ciPz/f1/0jHyiq0GPS8j14+5dTEAGssfZBmRgNq4tBmggBVm08DG1G+qWfFHhDMoZCbVFDdPHmECAgRhWLwRm5fupZcLNaRVhd/YcLY6yPNDwO+LKyMkyYMAFDhw7FpEmTAAA5OTnIzc3Fa6+95vMOknd+OFSCe5fuwp7CCtt9BmjwhGkSAMFhyIuoPQ6rGzYPUEVI31GShTcWLcaRhUcgXNrYxp66+2d0fRoapdaf3SMiD3g8RZ+QkICvv/76ivsXL17skw6Rb+hNFry2/hg+31tk9/F11p6YaHoc87RvIEa8PD1m29FOGwfdsHkwthrury5TgNXfW/6anoMwZ+/s2p3soIAVVghQQIQV0aoYzOj6NPqnDwx0l4NGdHTtF6GqKu4bYA/rIw3JtqodP3489u/fj5qaGsTHxyMvLw9dunS5oh23JvS9I+erMHP1ARwrdbzIT6kQ8OjALNzbLQWRh79A3I9PAgAMWSNgyB5ZOy3PkXvYsHfhGKPFgA3FP2L7xZ9RZihDJGIwMGMQBmfkcuTuIW7F6hzr4xl3t6rlXvQhRBRFfLanCPM2HIPB7Pg0uGYJEXhxZAd0yrj0JjFVI/XNdgCAkomHuKAuzFwO9ycxffo/rthrgh++3mMNnWN9PONuwHO/0RBRaTDjhW8P4YfD5522u75DGqYPy0aMk41yKHzwkq9EoYuf8iHgwFkd/v7VAZwp1ztsE6lW4Klh2RjZMZ0f4gSA4U4U6hjwQUwURXyyuxDzNhyDyeL4SEv7tBjMHnkVWiZx6p1qMdyJQh8DPkhVGsx4ce0hrDvkfEr+Tz2b4tGBraBR8crAVIvhThQeGPBB6OBZHf6++gAKyhxPySdEqvHs9e0xoHWSH3tGcsdwJwofDPggIooiPt1ThNfWH3U6Jd+taRxeHNkB6bE8lYkua2y4e3KtBrKPNXSO9ZEGAz5IVBrMmL32ML4/VOK03X29m2PSgCyoFByZ0WXejNwtFl550FusoXOsjzQY8EHgcEklpn/5G047mZKPj1DhuRuvwoBWnJKnhjgtTxSeGPAyt+a3s/jnd4edblzTtUkcZt/EKXm6ki/CPSpKAwCorjb6unthgzV0jvWRBgNepoxmK15dfxSf7bG/l3yd+3o3x6T+LaFScpU8NeSrkXvdPuH88G081tA51kcaDHgZKq7QY8ZXB7C/2PG2jfERKjx3w1VcJU92cVqeiBjwMrP1xEXM/N8BlOvNDtt0zozDP2+6ChlxvBgMXYnhTkQAA142rKKI97aexuubT8DZ1X/u6t4EUwe3hppT8mQHw52I6jDgZaBCb8KzX/+OjccuOGwTqVbgHyPaYcRVaX7sGQUThjsR1ceAD7BD5yrx1Je/Ob1QTMvESMy5uSPapET7sWcUTBjuRPRHDPgAWnvwHJ7/9pDTU+CGtk3B09e14+VdySGGOxHZw9QIAItVxKJNJ7B0+2mHbZQC8Nig1ri3Z1N+YJND/gj3ixe5jai3WEPnWB9pMOD9rEJvwsz/HcQvJy46bJMUpcY/b+qAns0T/NcxCjp5efP9MnI3O5lhIvewhs6xPtJgwPvR0fNV+Nuq/U6vAte1SRz+NaoDUmO4Kx05xml5InKFAe8nPx4+j2e//h3VJovDNmO6NcHjQ1pzVzpyyt/hHhmpBgDU1Jgk/TmhjDV0jvWRBgNeYlZRxJKfT+KtX045bKNWCpg+LBu3dM70Y88oGAVi5B4TU7uhEj98G481dI71kQYDXkKVBjNmrTno9Pz2lGgNXr65Izo3ifNjzygYcVqeiDzBgJfI6Ys1ePyLfJy4UOOwTefMOLx8cwek8Hg7ucBwJyJPMeAlsPN0GaZ/+ZvT/eRv6ZyBp4ZmQ6Pi8XZyjuFORI3BgPexL/YW4aV1R2Cx2t9RXqkQ8OTQNritSyY/qMklhjsRNRYD3kcsVhELfjqGj3aecdgmKUqNl0Z1RPdm8X7sGQUrhjsReYMB7wOVBjP+8b+D2Hzc8WK6DukxePnmjrzEK7mF4U5E3mLAe+lMeQ0eX7kfx0qrHbYZ1i4Fz17fHhFqpR97RsFKbuFeVub4vU3uYQ2dY32kwYD3wu6Ccjz55W8oc3Lu5oN9W2BC/5ZQcARGbpBbuAOAycnmTOQe1tA51kcaDPhG+iq/GP/87jDMDhbTaZQCZl3XHtd14PXbyT1yDHciCl4MeA+JoojFm0/g3a2OrwSXFKXGK6M7ISeTm9eQe+Qc7hERtduI6vXcZayxWEPnWB9pMOA9YDRb8fy3v+PbgyUO27RNjcaroztxMR25Tc7hDgCxsbXvZX74Nh5r6BzrIw0GvJsq9CY8ueo37Cood9hmSHYynrvhKkRpuJiO3CP3cCei4MWAd8OZ8hr89XPn287e17s5HhmYxcV05DaGOxFJiQHvwm/FOkxbmY8L1fanjpQKAf8Y0RY3dcrwc88omDHciUhqDHgnfjpaipmrD0Bvttp9PFqjxMs3d0Tvlol+7hkFM4Y7EfkDA96BT3YXYu4PR+DgLDikx2ox77YcZKdE+7djFNQY7kTkLwz4P7CKIhb+dBwf7ihw2KZ9Wgxeu7UTUnmZV/IAw52I/IkBX4/RbMUzX/+O7w85Pg2uf6tE/POmDojWsHTkvmAO9/JybiPqLdbQOdZHGkypS6qMZvxt1W/YcarMYZtbu2TgqWFtoVIEz4czBV4whzsAGI3cRtRbrKFzrI80GPAALlQb8dfP83HgbKXDNo8OzMJ9vZsH3YczBVawhzsRBa+wD/gz5TWY8lk+Tl20f467+tKe8tdzT3nyUKiEu1Zb+zFhMJgD3JPgxRo6x/pII6wD/khJFSZ/tg/nq4x2H4/RKjH3lk7o2TzBvx2joBcq4Q4AcXGRAICSEl2AexK8WEPnWB9phG3A7y4ox+Nf7IfOwTfGlGgNFtyeg7apMX7uGQW7UAp3IgpeYRnwPx0txf+tPgCDgw1sWiRGYsHtOWgaH+nnnlGwY7gTkVyEXcB/lV+M2WsPweJgA5ur0mIw//YcJEVp/NsxCnp5efMZ7kQkG2EV8B9sP40FPx13+HivFgn4980dEaMNq7KQD3DkTkRyExZJJooi/rPpBN7fdtphm2HtUvD8DVdBo1L4sWcUChjuRCRHIR/woijilR+PYvmvhQ7b3N41E08OzYaSG9iQhxjuRCRXIR3wFquIl74/jC/2FTtsM6FfC0zo15IfzOSxcAn38nL7e0SQ+1hD51gfaYRswJutIp7/5nd8feCc3ccFAH8b2gZ3dm/q345RSAiXcAcAo5Gbj3iLNXSO9ZFGSAa8yWLFP/53ED8cPm/3caUAPHvDVdydjholnMKdiIJXyAW8wWzFjK9+w6ZjF+w+rlII+OdNHZDbNsXPPaNQEI7hrrl05USOshqPNXSO9ZFGSC0ZrzFZMG1lvsNw16oUmDu6E8OdGiUcwx0A4uMjEc9Nn7zCGjrH+kgjZEbwlQYzpq3Mx+4zFXYfj1Qr8MroTri6RaKfe0ahIFzDnYiCV0gEfHmNCVM+z8dvxfYvVBCtUWL+bTno2jTezz2jUMBwJ6JgFPQBX1ZjwqOf7MWhkiq7j8dHqLDwjs7okB7r555RKGC4E1GwCuqAr9Cb8Nin+xyGe1KUGv+5owuyU6P93DMKBQx3IgpmQRvwOr0Zj326D7+fq7T7eFqMBv8Z0wVZSVF+7hmFAoY7EQW7oAz4SoMZkz/bhwNn7Yd7kzgtFt3ZhZd7pUZhuBNRKAi6gK8ymjHls3zsd7Cgrkl8BN64swsy4iL83DMKBQz3K1VUcBtRb7GGzrE+0giqgK82WjD1s3zsK7J/KlyTOC1eZ7hTIzHc7TMYuPmIt1hD51gfaQTNRjc1Jgv+ujIfewrth3tGrBaL7+yKTIY7NQLDnYhCTVAEvN5kweMr8/FrQbndx9NiNFh8Zxc0iWe4k+cY7s5pNEpoNMpAdyOosYbOsT7SkH3A600WPPHFfuw4bT/cU2M0WHxnVzRL4II68hzD3bX4+CjEx/NsFG+whs6xPtKQdcAbzFY8ueo3bDtVZvfx5GgNFo3pghaJDHfy3OVwf5LhTkQhR7aL7MxWETNXH8AvJy/afTwpSo3FPM+dGokjdyIKdbIcwVtFEbPXHsKGo6V2H0+MVGPRmC5olcxwJ88x3IkoHMgu4EVRxPwNx7B6/1m7j8dHqLBoTBe0SeH2s+Q5hjsRhQvZBfx7207jo51n7D4Wdyncubc8NQbDnYjCiawC/vM9hVi06YTdxyJUCsy7NQft0mL82ykKCQx3Igo3sllk993vJXjp+yN2H1MpBPz7lo7o3CTOz72iUMBw945Opw90F4Iea+gc6yMNWQT8lhMXMGvNQYh2HhMAPH/jVeibleTvblEIYLh7T683BboLQY81dI71kUbAp+j3FlbgqVW/wWy1F+/AjOFtMbx9qp97RaGA4U5E4SzgAT9tZT70Zqvdxx4ZmIXbumT6uUcUChjuvqNWK6FWcxtRb7CGzrE+0gh4wFfo7V9F6N6ezXB/7+Z+7g2FAoa7byUkRCEhgXtOeIM1dI71kUbAA96eUZ3SMXVwK34wk8cY7kREtWQX8EOyk/F/I9rxg5k8xnAnIrpMVgHfq0UCXhzZASoFP5jJMwx3IqKGZBPwnTJiMfeWjtCqZNMlChIMdyKiKwX8PPjrrkpF29QY3NOjKTQMd/IQw52IyL6AB/yLIzsEugsUpBjuRESOBTzgiRqD4e4/lZXcRtRbrKFzrI80GPAUdPLy5jPc/aimhtuIeos1dI71kQYDnoIKR+5ERO7hqjYKGgz3wFCpFFBxAaxXWEPnWB9psKIUFBjugZOYGI3ExOhAdyOosYbOsT7SYMCT7DHciYg8x4AnWWO4ExE1DgOeZIvhTkTUeAx4kiWGOxGRdxjwJDsMdyIi7zHgSVYY7kREvsGNbkg2GO7yVFVlCHQXgh5r6BzrIw0GPMkCw12+qquNge5C0GMNnWN9pMEpego4hjsRke8x4CmgGO7yp1QqoFTyo8IbrKFzrI80WFEKGIZ7cEhKikZSErcR9QZr6BzrIw0GPAUEw52ISFoMePI7hjsRkfQY8ORXDHciIv9gwJPfMNyJiPyHAU9+wXAnIvIvBjxJjuFOROR/3MmOJMVwD37cZcx7rKFzrI80GPAkGYZ7aOA+4d5jDZ1jfaTBKXqSBMOdiCiwGPDkcwz30KJQCFAo+G/oDdbQOdZHGgx48imGe+hJTo5BcnJMoLsR1FhD51gfaTDgyWcY7kRE8sGAJ59guBMRyQsDnryWlzef4U5EJDM8TY68wpE7EZE8cQRPjcZwJyKSLwY8NQrDnYhI3jhFTx5juIeXmhpuI+ot1tA51kcaDHjyCMM9/FRWchtRb7GGzrE+0uAUPbmN4U5EFDwY8OQWhnv4EgSB/95eYg2dY32kwYAnlxju4S0lJQYpKdxG1BusoXOsjzQY8OQUw52IKDhxkV2YM1oMWF/4HXakpaBMoUDknmcxIDMXQzKG4o1FixnuRERBShBFUQxkB0pKdIH88WFt89mNmLPnRVSadVCIIqyCAAECRIhQW9Q4svAIJgx9mOEe5lJTYwHw/6o3WEPnWB/P1NXLFY7gw9Tmsxsxa+cM223rpQAXUft9zygY0XJqS1zTcxDDnYgoCPEYfBgyWgyYs+dFAJcD/Y8ERe2q1jl7Z8No4TmqRETBhgEfhtYX/4BKs85huNcRIaLSrMOG4h/91DMiIvIVBnwY2ly8EYKb//QCFNhU/JPEPSI50+tN0OtNge5GUGMNnWN9pMFj8GGowlQOEVa32oqwosJULnGPSM50On2guxD0WEPnWB9pcAQfhuLU8R6N4OPU8RL3iIiIfI0BH4YGZFzj0Qh+YMYgiXtERES+xoAPQ0MyhkJtUUO0Ol9kJ0BAjCoWgzNy/dQzkqPU1Fi3z7sl+1hD51gfaTDgw9AbixbjyMIjtRd4gP1z3Ovun9H1aWiUWn92j4iIfIABH2bq9pafMPRhvNDzJUSrai/woLj0Vqg7Nh+tisELPeegf/rAgPWViIgaj6vow4i9C8d8mtoHG4p/xPaLP6PMUIZIxGBgxiAMzsjlyJ2IKIgx4MOEo6vCaZRaDG96Pf7UbQwA7gVNRBQqOEUfBnjJVyKi8MOAD3EMdyKi8MQp+hDmSbgbDNwmkuzje8N7rKFzrI80eD34EMWROxFRaHJ3zwBO0YcghjsRETHgQwzDnYiIAAZ8SPEm3LlVJDnC94b3WEPnWB9pMOBDBEfuRERUHwM+BOTlzWe4ExFRAzxNLshx5E5ERPZwBB/EGO5EROQIAz5IMdyJiMgZBnwQYrgTEZErPAYfZKQKd4PB7JPXodDD94b3WEPnWB9pcKvaIMKROxERcavaEMNwJyIiTzDggwDDnYiIPMWAlzl/hXtKSgxSUmIkeW0KbnxveI81dI71kQYX2cmYP0funBUgR/je8B5r6BzrIw2O4GWK0/JEROQNBrwMMdyJiMhbDHiZYbgTEZEvMOBlhOFORES+woCXCYY7ERH5ElfRy4Acwt1o5FaRZB/fG95jDZ1jfaTBrWoDTA7hTkREwYNb1QYBhjsREUmFAR8gDHciIpISAz4A5BjuycnRSE6ODnQ3SIb43vAea+gc6yMNLrLzMzmGOwAoFPyuR/bxveE91tA51kcarKofyTXciYgo9DDg/YThTkRE/sSA9wOGOxER+RsDXmIMdyIiCgQGvIQY7kREFChcRS+RYAt3k8kS6C6QTPG94T3W0DnWRxrcqlYCwRbuREQUPLhVbYAw3ImISA4Y8D7EcCciIrlgwPtIsId7UlI0kpK4VSRdie8N77GGzrE+0uAiOx8I9nAHAKWS3/XIPr43vMcaOsf6SINV9VIohDsREYUeBrwXGO5ERCRXDPhGYrgTEZGcMeAbgeFORERyx4D3UF7efIY7ERHJHlfReyCUR+5mM7eKJPv43vAea+gc6yMNblXrplAOdyIiCh7cqtaHGO5ERBRsGPAuMNyJiCgYMeCdCKdwT0yMQmJiVKC7QTLE94b3WEPnWB9pcJGdA+EU7gCgUikD3QWSKb43vMcaOsf6SIMjeDvCLdyJiCj0MOD/gOFOREShgAFfD8OdiIhCBQP+EoY7ERGFEgY8GO5ERBR6wn4VPcO9lsViDXQXSKb43vAea+gc6yONsN6qluFORETBhlvVusBwJyKiUBaWAc9wJyKiUBd2Ac9wty8hIQoJCdwqkq7E94b3WEPnWB9phNUiu8vh/iTD/Q/Uam4VSfbxveE91tA51kcaYRPwHLkTEVE4CYsp+saE+6lTJzF8+DVYs+Yr233//OdzGDiwV4Nft956o+3xPXt2489/vgt33nkLNmz4UZK/CxERkTtCfgTfmHA3m814/vmnUVNT0+D+I0cOY9y4v+COO+6y3adQXJ5amjPnBUyaNBnR0dF48cVn0K/fAGg0Gt/9ZYiIiNwU0gHf2Gn5t99+A1FRDRd8WCwWnDhxHPff/yCSk1PsPs9kMiE7uy0iI6OgVCphNpsZ8EREFBAhG/CNDffdu3dh1arP8e67y3D77TfZ7i8oOA2j0YCsrFYOn/vAAxMxduydsFjMePDBSVd8SSAiIvKXkAz4xoa7TqfDCy/Mwl//+iTS0zMaPHb06BEIgoAVKz7G1q0/QxAE9Os3ABMmPIKYmBgAwA033IQhQ4bBYrHY7gsWViu3iiT7+N7wHmvoHOsjjZALeG9Wy8+d+y/k5HTGiBHXX/HY8eNHoVAokJmZiTlzXsXp06fxn//Mw9GjR7BgwetQKGrXK0ZGRvrs7+JPpaVVge4CyRTfG95jDZ1jfaQRUgHvTbh/883/sHfvbrz//n/tPv7AAxNx551/Qmxs7R7ArVtnIzk5BZMm/QUHDvyGTp1yfPJ3ICIi8oWQCXhvz3P/3/++xIULpbj99pEN7p8791/46KMP8OGHK2zhXqdNm2wAQEnJWQAMeCIiko+QCHhfbGIza9YLMBgMDe67++5bMX78Qxg2bASefXYmKirK8eqrebbHDxzYDwDIymrt3V9ABuLjaw8tlJfXuGhJ4YbvDe+xhs6xPtII+oD31Q51qalpdu9PTExCRkYmhg+/HjNmPI73338b1157HU6dOolXX30Zw4df73RlfbDQaIL+rUAS4XvDe6yhc6yPNIK6qv7cfnbAgGvwwgsvYenSd7B06TuIiYnF8OHXY8KEhyX7mURERI0liKIoBrIDJSW6Rj2Pe8v7Vmpq7fqCxv57UOjie8N7rKFzrI9n6urlSlDuRc9wJyIici7oAp7hTkRE5FpQBTzDnYiIyD1Bs8iO4S6tAC/FIBnje8N7rKFzrI80gmKRXV7efDz//NMMdyIiCnvuLrKT/QieI3ciIiLPyfoYPMOdiIiocWQb8Ax3/4qLi0RcXHBeCY+kxfeG91hD51gfachyip7h7n9arSzfCiQDfG94jzV0jvWRhuxG8Ax3IiIi78kq4BnuREREviGbgGe4ExER+Y4sAp7hTkRE5FsBD3iGOxERke8FPOAZ7kRERL4X8K1qiYiIyPcCPoInIiIi32PAExERhSAGPBERUQhiwBMREYUgBjwREVEIYsATERGFIAZ8GFq9ejU6d+6Ms2fP2u6bOXMmRo0aBZ1OF8CeUaDxveE91tB7dWdvL1y4sMFt8gwDPgyNHDkSrVq1whtvvAEAyMvLw6ZNm/DWW28hNjY2wL2jQOJ7w3usofc++ugjLF++HNXV1Zg7dy62b98e6C4FJ5HC0g8//CDm5OSIixcvFq+++mrx0KFDDR4bMWKEOHz4cHHFihUB7CUFgrP3xiOPPCL26tVLnDx5cgB7KH+OalhYWCiOHTtWvOGGG8SbbrpJXLNmTYB7Kl9LliwRO3XqJO7YsSPQXQla3MkujN166604fPgw3nnnHfTu3RsAYDabMXLkSCxduhTR0dG47bbbsGLFCiQkJAS2s+RX9t4bAPDLL7+guroaX3zxBRYsWBDAHsqfvRqeO3cOpaWl6NChA0pLS3Hrrbfim2++QVRUVIB7Ky8ff/wxFAoFTpw4AZVKhcGDB6NXr16B7lbQ4RR9mNq4cSOOHTsGq9WK5ORk2/179+5FdnY20tPTERMTg0GDBmHTpk0B7Cn5m6P3BgD07dsX0dHRAepZ8HBUw7S0NHTo0AEAkJycjPj4eJSXlweqm7J1991346677kJUVBSeeOIJ9OzZM9BdCkoM+DC0f/9+TJ06Fc899xwGDhyI+fPn2x47d+4c0tPTbbczMjIaLBai0ObsvUHucbeG+/btgyiKyMzM9HMPA8fdBYh1Fx6bPHlyg9vkGQZ8mCkoKMBDDz2EiRMnYvTo0ZgyZQrWrl2L/Px8APZXq/I/V3hw9d4g19yt4cWLFzF9+nQ8//zzAeppYHABon8x4MNIWVkZJkyYgKFDh2LSpEkAgJycHOTm5uK1114DAKSnpzf4dl1cXIzU1NSA9Jf8x533Bjnnbg2NRiMee+wxTJw4ET169AhUdwNCEARMmzYNn3zyCV5//XUsXboUb731VoNZQ/IdLrKjBsxmM2688UZ88MEHtkV2y5cvR2JiYqC7RjKxdetWLFu2jIvsGkEURTzxxBNo1aqVbfo5HDlaxEm+pQp0B0heVCoVpk+fjj//+c+wWq148MEHGe5kM378eOzfvx81NTUYNGgQ8vLy0KVLl0B3K2js3LkTa9asQfv27fH9998DAF5++WW0b98+wD3zH2eLOMm3OIInIiK/2L9/P8aNG4dZs2ZhzZo1iIiI4EyQhDiCJyIiyf1xAWJ2djbuuOMO5OfnIycnJ9DdC0kcwRMRkaTKyspwzz334Oqrr25w5sDDDz8Mo9GIt99+O4C9C10MeCIiohDE0+SIiIhCEAOeiIgoBDHgiYiIQhADnoiIKAQx4ImIiEIQA56IiCgEMeCJiIhCEAOeiIgoBDHgiYiIQhADnoiIKAQx4ImIiEIQA56IiCgE/T8RP3fi5v2mkQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "def g(x):\n", " return (x + 0.2)**0.5\n", "\n", "xmin, xmax = 0.0, 1.4\n", "xinit, xstar = 0.3, 0.5*(1 + np.sqrt(1.8))\n", "xx = np.linspace(xmin,xmax)\n", "yy = g(xx)\n", "\n", "n = 21\n", "z = np.zeros(n)\n", "z[0] = xinit\n", "for k in range(n-1):\n", " z[k+1] = g(z[k])\n", "x, y = z[:-1], z[1:]\n", "\n", "xticks = x[:3].tolist() + [xstar]\n", "yticks = y[:3].tolist() + [xstar]\n", "\n", "fig1, ax = plt.subplots(figsize=[6,6])\n", "ax.set(\n", " aspect=1.0,\n", " title='Function Iteration',\n", " xlim=[xmin,xmax],\n", " ylim=[xmin,xmax],\n", " xticks=xticks,\n", " xticklabels=['$x_0$', '$x_1$', '$x_2$', '$x^*$'],\n", " yticks=yticks,\n", " yticklabels=['$x_1$', '$x_2$', '$x_3$', '$x^*$']\n", " )\n", "\n", "ax.vlines(xticks, xmin, xticks, colors='white', linestyles='--')\n", "ax.hlines(yticks, xmin, yticks, colors='white', linestyles='--')\n", "\n", "\n", "ax.plot(xstar,xstar,color='white', marker='.', ms=20)\n", "ax.plot(xx,xx,'k-', linewidth=1)\n", "ax.plot(xx,yy,linewidth=4)\n", "ax.step(x,x,'C1')\n", "ax.plot(x[:4], x[:4], color='C2', **bullet_options)\n", "ax.plot(x[:4], y[:4], color='C1', **bullet_options)\n", "ax.annotate('45°', [xmin+0.1,xmin+0.05], ha='left', fontsize=11)\n", "ax.annotate('g', [xmax-0.05, g(xmax)-0.08], ha='left',fontsize=18, color='C0');" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Newton's Method" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAs0AAAINCAYAAADFp0I1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAABzlElEQVR4nO3dd3xV9eH/8de5M7n3ZgAJe++tsmUooqICzrY/tY7aaltbV6l1VAXcdW+/bq22ttW2WltFRXEAshFky15hj4y7xzm/PwJRJJAEkpyb5P18PHjkknvuOe8k9ybvnHzO52NYlmUhIiIiIiKH5bA7gIiIiIhIulNpFhERERGpgEqziIiIiEgFVJpFRERERCqg0iwiIiIiUgGVZhERERGRCqg0i4iIiIhUQKVZRERERKQCKs0iIhVoyGtANeSPXUTk+1SaRSRtXXbZZfTs2ZMlS5aUe/+oUaO49dZbazTDmjVruPjii2tk36NGjeKyyy476se/8847dOvWjW7durF+/fpyt5k2bVrZNlVRXFzMLbfcwvz588ved9lllx1T3orU9P5FRI6FSrOIpLVUKsUf//hH4vG4Lcf/8MMPWbhwoS3HriyHw8GHH35Y7n2TJ08+qn2uWLGC//znP5imeSzRRETqDZVmEUlrWVlZrF69mmeffdbuKNWuadOmNGvW7Jj3069fv3JLczwe59NPP6VHjx7HfAwRkYZOpVlE0lqPHj0477zzePnll1m6dGmF2//zn/9k7Nix9O7dm5EjR/L000+TTCYBuP/++xk0aNBBZ08nTJhAt27dWLduXdn7/va3v9G3b18eeughnnnmGQC6devG008/DUAsFuPZZ5/lzDPPpE+fPowePZoXX3zxoP1edtll3H777bz44ouMHDmSPn36cNFFF/HNN9+UbdOlS5eyYROmafLkk08yatQoevfuzahRo3jsscdIJBIVfsxjxoxh1apVrF279qD3T5s2DcMwOOmkkw55zPz587n00ks57rjjGDRoELfccgt79+4FYM6cOVx++eUAXH755QcNmbAsi5deeomRI0fSt29fLrzwwkOGzyxZsoQrr7ySwYMH069fP66++mpWr1590DZbt27l2muvpX///gwbNozXXnutwo9TRMROKs0ikvZuv/12GjduXOEwjRdeeIEJEyZw4okn8vzzz3PJJZfw0ksvMXHiRABOOeUUioqKDirfs2fPBmDevHll75s2bRonnngil19+OT/+8Y8BeOutt/jJT36CZVlcffXVvPzyy/z4xz/m+eef58wzz+SJJ55g0qRJB+X5+OOPmTp1KnfccQePPfYYu3fv5vrrryeVSgFwzz338Mtf/hKAl156iTfffJNrrrmGV199lYsvvpiXX36Z559/vsLPz7Bhw8jJyTnkbPPkyZM5/fTTcbvdB71/3rx5XHHFFWRkZPDEE09w2223MXfuXC6//HKi0Si9evUq+5xNnDjxoI9rwYIFfPLJJ0yYMIEHH3yQHTt2cPXVV5f9YjJ79mwuvvhiTNPkvvvu495772Xbtm1cdNFFZaU+HA5z6aWXsnLlSu6++24mTpzIP//5z7QfBiMiDZvL7gAiIhXJzs7mrrvu4je/+Q3PPvss48ePP2SbkpISnnvuOS688ELuuOMOAIYPH05ubi533HEHP//5zxkwYAB+v59Zs2bRt29fCgoK2LRpE7169WLu3LlceOGFxONx5syZwy233ELz5s1p3rw5AMcffzwAX375JTNnzuThhx/mnHPOAUpLa0ZGBk8++SQ/+9nP6Ny5MwDJZJJXXnmFQCAAQCgU4pZbbmHFihX07t37oPxz586lV69e/OhHPwJg0KBBZGZmlj32SFwuF6eddhoff/wx1157LQCRSITPP/+cZ599lgULFhy0/aOPPkqHDh144YUXcDqdABx33HGMHTuWf//731xyySVlH0Pnzp3LbgN4PB5efPFFcnNzAQgGg9xxxx2sWbOG7t278+ijj9KmTRtefvnlsn0PHz6c008/naeffponnniCd999l61bt/Lee++VnWnv27cvp59+eoUfq4iIXXSmWUTqhFGjRnHOOefw8ssvs2zZskPuX7hwIZFIhFGjRpFMJsv+jRo1CoCvvvoKt9vNsGHDmDlzJgCzZs2iXbt2nHXWWcydOxcoPQsbDoc55ZRTys0xd+5cnE4nY8aMOej9Bwr0nDlzyt7XuXPng0rvgfHLkUjkkP0OHjyYmTNn8tOf/pTXXnuNtWvXcumll3LeeedV6vPzwyEan3/+OT6fj8GDBx+0XSQS4ZtvvuHkk0/Gsqyyz1ObNm3o1KkTX3311RGP07lz57LCDNC6dWug9JeWcDjMkiVLGDNmTFlhhtJfek455ZSyz838+fNp06bNQTN6tGjRouwXExGRdKTSLCJ1xh133EGjRo249dZbDxmmUVhYCMCvfvUrevXqVfZv6NChAOzcuROAk08+ma+//ppoNMqsWbMYPHgwgwcPZufOnWzYsIFp06bRq1evw16gV1RURKNGjXC5Dv5DXX5+PlBaHg/IzMw8aBuHo/RbbnkzUlx11VVMnDiRaDTKgw8+yJgxYzj77LOZNWtWpT43Q4YMoVGjRmVDNCZPnsyZZ555UHmF0qnkTNPkpZdeOujz1KtXL1atWlX2eTocn8932I+ppKQEy7LIy8s75HF5eXlln5uioiIaN258yDYHPociIulIwzNEpM7Iycnhzjvv5JprruG555476L7s7GwAHnnkEdq3b3/IYw8UuZNPPplEIsGCBQvKhmH06tWLQCDAvHnzmDZt2iFnkX+YYd++fSSTyYOK84Gy2ahRo6P62BwOB5dccgmXXHIJe/bs4csvv+T555/nuuuuY+bMmXg8niM+3uVyMXr0aD766COuuOIKpk2bxp///OdDtvP7/RiGwRVXXMHYsWMPuf+HRb8qsrKyMAyD3bt3H3Lfrl27ys5QN2rUiI0bNx6yzYFffERE0pHONItInXLaaacxbtw4XnzxxbLZHqB0TK7b7WbHjh306dOn7J/b7ebRRx9ly5YtQOnZzJ49e/L3v/+dXbt2MWjQIJxOJwMHDuTdd99l3bp1Bw3NOHAm9YBBgwaRSqUOmf/4v//9LwD9+/c/qo/roosu4t577wWgSZMmXHDBBVxyySWUlJQQDAYrtY8xY8awevVqXn75ZfLy8jjhhBMO2SYQCNCzZ0/WrVt30OepS5cuPPPMM2VDKH54hroyfD4fvXv3ZvLkyWUXO0Lp2fcvvvii7HMzZMgQtmzZctCsG3v37mXRokVVPqaISG3RmWYRqXMmTJjA7NmzDzqj2ahRI6666iqefPJJgsEggwcPZseOHTz55JMYhkH37t3Lth05ciTPPvssHTp0KBuGMXjwYB544AGaNm1Kr169yrY9cAb7/fff57jjjuOkk05i8ODBTJo0iZ07d9KzZ0/mzp3LSy+9xPnnn3/QRXNVMXDgQF599dWysrtjxw5ee+01Bg0aVO5QhvIMGjSI/Px8Xn75Za644goMwyh3u9///vf86le/4sYbb+Scc84hlUrx6quv8s033/Cb3/wGKD1rDPDFF1+Qk5Nz0OfvSG688UauvPJKrrrqKi699FISiQQvvvgi8Xi87CLFc889lzfeeINrr72W8ePHEwgEeO6557SQioikNZ1pFpE6Jzc3lzvvvPOQ9//ud7/j1ltv5ZNPPuGXv/wlDz/8MP379+evf/1rWQmE0tIMpSXzgAMXzI0cOfKgsjl69Gj69OnDrbfeyiuvvIJhGLzwwgtcdNFFvPHGG/zqV7/io48+Yvz48dx3331H/THdcMMNXH311fz73//mqquu4oEHHmD48OE89dRTld6Hw+HgjDPOIJFIlDv04oDhw4fzyiuvsH37dq6//npuvvlmnE4nr732WtnFeF26dGHcuHG8+eab/OEPf6h0hhNPPJHXXnuNeDzO73//eyZMmECzZs14++236dKlC1A6A8frr7/OwIEDue+++7jttts48cQTD3vxpYhIOjAsy7LsDiEiIiIiks50pllEREREpAIqzSIiIiIiFVBpFhERERGpgEqziIiIiEgFVJpFRERERCqg0iwiIiIiUgGVZhERERGRClR6RcBdu0pqMoeIiIiIiC3y87Mq3EZnmkVEREREKqDSLCIiIiJSAZVmEREREZEKqDSLiIiIiFRApVlEREREpAIqzSIiIiIiFVBpFhERERGpgEqziIiIiEgFVJpFRERERCqg0iwiIiIiUgGVZhERERGRCqg0i4iIiIhUQKVZRERERKQCKs0iIiIiIhVQaRYRERERqYBKs4iIiIhIBVSaRUREREQqoNIsIiIiIlIBlWYRERERaZCCsST/WrS1Utu6ajiLiIiIiEi1sSyLvy0o4N3F29hWHCXT7WRk5zzuOKNrlfd1z8er+Gz1bn5zercKt1VpFhEREZE6480FBTz55ToAWmZ7yc5w0zIno8r7+aagiM9W78ZhVG57lWYRERERqTM+/XYXAKO75XPfuB5HtQ/Lsnhif/E+p3fzSj1GY5pFREREpM4ojCQA6N8m56j38emq3SzdVkKm28Gvh7ar1GNUmkVERESkzjAtCwC38+hqbDxp8sz09QBcNrANeQFvpR5nWNb+I1dg166SowomIiIiInKsfv3WN3y9peiQ97fI9vLfXw6u9H7enL+FJ75cR57fwztXDiTT7SQ/P6vCx2lMs4iIiIikvc55flKmxYodJcRTFm0bZdIo001ewFPpfRRFErwyexMAvxnWnky3s9KPVWkWERERkbR306mdATjnpTlsK45xxaA2nF3Ji/gOeHXOJkpiSTrn+Rnbq1mVHqsxzSIiIiJS720pjPD2wtKFTG44uQPOys41t59Ks4iIiIjUe89MX0/StBjSvhFD2jeu8uNVmkVERESkXvumoIipq0oXMrnh5I5HtQ+VZhERERGptyzLKltB8Ozezemc5z+q/ag0i4iIiEi9NXXVbpZsKyHD5eDqSi5kUh6VZhERERGpl76/kMnlVVjIpDwqzSIiIiJSL/3rm60UFEXJ83u4dGDrY9qXSrOIiIiI1DvfX8jk6mHtqrSQSXlUmkVERESk3nl1ziaKo0k65fkY16tqi6CUR6VZREREROqVgxcy6VjlhUzKo2W0RURERKTO+O8vB1e4zbPTN5QuZNKuEScexUIm5dGZZhERERGpNxZvLebTVbswgOtP7lBt+1VpFhEREZF6wbIsnvjiwEImzeiSH6i2fas0i4iIiEi98Nnq3SzZVkyGy8Gvh7av1n2rNIuIiIhInZdImTw9rXQhk0sHtKZp1tEvZFIelWYRERERqfP+uah0IZPGPjeXDWxT7ftXaRYRERGROq04+v2FTNrj8xzbQiblUWkWERERkTrtldnfLWRyTu9jX8ikPCrNIiIiIlJnbdgb5q39C5lcd1L1LGRSHpVmEREREamzHv9iLSnTYliHxgzrUD0LmZRHpVlERERE6qQZ6/Ywc/0+XA6D8SM71uixVJpFREREpM6JJ00e+3wtABf3a0W7xr4aPZ5Ks4iIiIjUOf/4uoDNhaVTzP1iSNsaP55Ks4iIiIjUKbuDsbIp5q47qQMBr6vGj6nSLCIiIiJ1yjPT1xNOpOjVPIsxPZvVyjFVmkVERESkzliytZgPlu8E4KZRnXAYNTPF3A+pNIuIiIhInWBaFo/sv/hvXK9m9GqRXWvHrvkBICIiIiIiRyMZxbv2fbzrPsaIFlIQz6Dbzu5s8wzlmhEdajWKYVmWVZkNd+0qqeksIiIiIiIAeNZPIWvqeByxIiwcGJikMHBiEXVmETvjKeIdTq+WY+XnZ1W4jYZniIiIiEha8ayfQvbkKzFixQAYmAA4KT3X600FyZ78Czzrp9RaJpVmEREREUkfyShZU8cDYFD+gIgD78+aOh6S0VqJpdIsIiIiImnDu/Z9HLGiwxbmAwwsHLEivGs/qJVcKs0iIiIikja86z7GqmRFtXDgXfdRDScqpdIsIiIiImnDiBaWjWGucFtMjGhhzQbaT6VZRERERNKGlZFbpTPNVkZuzQbaT6VZRERERNJGrOMZVTrTHOt4Zg0nKqXSLCIiIiJpI9ZpHKY3p4LLAMHCwPTmEOs0tlZyqTSLiIiISPpwZbCq+/UcqTVbGACUnPoEuDJqJZZKs4iIiIikDdM0CS17H8OAJE7gu5Jc9tabTfGYV6ttRcDKcNXakUREREREKrB42tucnlxI3HKx+bwPyA8uJ3PpXyARwszpQKzTWaVDMmrpDPMBKs0iIiIikhZCoSA9lj0EwMKWl9CxdS9iVg+STXoCkMrrAYY9AyVUmkVEREQkLaye/Ains5OdNKHtmFtL32k4SOX3sjcYGtMsIiIiImlg86Y1DNvxFwA2Hn8TrowsmxMdTKVZRERERGxlWRahjyfiM2Ks9PSm/dBLDt7ATJX+s5FKs4iIiIjYav7MDxkRn0bKMnCcfj8Yxnd3mik8Gz/Ds/EzW4uzSrOIiIiI2KYwFKHTovsAWNL0fBq172dzovKpNIuIiIiIbRZ98Azd2EgxAZqPnWR3nMNSaRYRERERWyxfv4nTdr4CQEHf3+H0N7E50eGpNIuIiIhIrUuaFkWf3EMjI8gWT0fyhl1ld6QjUmkWERERkVr32YwvGBv/CADr1PvAkd7Lh6g0i4iIiEit2lEcpevi+3EaFmvzTiej4wi7I1VIpVlEREREatX0ya8y0FhBDC/ZZ9175I0NA9PfFNPf9OCp6GpZep8HFxEREZF6ZfbqLZy7+3kwYEfvq8nMbnXkBxgOkk2Pq51wR6AzzSIiIiJSK6KJFHumPkJLYy973S3IHHad3ZEqTaVZRERERGrFu9Nnc3HyPQCSI+8EV4a9gapAwzNEREREpMZt2BOm+7KH8ToSbG8yBGeXMZV74P5ltAHi7UaBw1mDKQ9PZ5pFREREpEZZlsVHH77FaMd8kjhxn36/rRf1HQ2VZhERERGpUR8vL+Difc8BsLf75ZhNutqcqOpUmkVERESkxpREk+z88jk6O7YScjXCMfwmuyMdFZVmEREREakxb3yxgF+abwMQH34bljfb5kRHR6VZRERERGrEsu0l9Pj2abKNCIU5vUj2vNDuSEdNpVlEREREql3KtHjnow+40PUFANZp94NRd6unppwTERERkWr3r4Vb+HnJ8+CAok4/Itm8/9HtyDAwM5uU3baLSrOIiIiIVKvdwRhbZ/6FExxriDv9JEfcdvQ7Mxwkm/ervnBHqe6eIxcRERGRtPTcZ0sYb7wJQGzQ7zD9zWxOdOxUmkVERESk2szasJde618i3ygiHGhP9Lgr7Y5ULTQ8Q0RERESqRSie5M2PP+evzo8AiI+8B5yeY9upmcKz6YvS/bUdadsy2irNIiIiIlItnpu+nmtiL+N2pgi3PY1Eu1OqZ8eWWT37OQYaniEiIiIix+ybgiL2Ln6fk5xLSDncREZMsjtStdKZZhERERE5JvGkyUMfL+Ml118AiB3/a8zcDjanql460ywiIiIix+TVOZs4rfhftHPsJOlrRqj/dXZHqnY60ywiIiIiR231riAfzV3EJ+73AAgPuwM8fptTVT+daRYRERGRo5I0Le75eBU3O/+Gz4gRbz6QWJfz7I5VI3SmWURERESOyj++LiCwcz7nemdiYRA66Z4aWerazGhU7fusKpVmEREREamyLYURXvxqHf92vw5AtNclJPN7V/+BHE6SLQZU/36rGsPuACIiIiJSt1iWxX2frOYC61N6OjZienMIDb7Z7lg1SmeaRURERKRK/rt0O6s3beEl79sAhAb9ASuzsc2papZKs4iIiIhU2q5gjCe+XMdNrn/SyAiSbNyNaO/Lau6AZgr3lukAJFqP0DLaIiIiIpLeLMvioalraB1fz6XeqQAER9wNjpqtlEYqUaP7rwyVZhERERGplM9W7+aLNbv5h+d1nJhEO40j0XqY3bFqhS4EFBEREZEKFUUSPDR1DWMdcxjiWIHlyiA09A67Y9UalWYRERERqdATX64jEi5hovdvAIT7XYOZ3drmVLVHpVlEREREjmj2hr28v2wHv3H9j2bWblJZrQmfcLXdsWqVSrOIiIiIHFY4nuJPn6ymtbGT37jfByA4bCK4Mm1OVrt0IaCIiIiIHNZzX21ga3GMP2f+HbeVIN56OPGOZ9VqBsubXavHK49Ks4iIiIiUa8nWYt76uoDhjiWMtOZgGU6Cw+8Cw6i9EA4niZaDa+94h4thdwARERERST/xpMk9U1bhJMnD/jcBiPS5glSTbjYns4dKs4iIiIgc4pU5m1i/J8yvMz6jRWITZkZjwoN+b3cs22h4hoiIiIgcZNm2Yl6fs4kmFHGD81+QgtCJt2J5c2o/jJnCXTATgESroVpGW0RERETsF02kuPOjb0lZ8ETee3iCQRL5fYl2v9C2TEYyatuxD9DwDBEREREp89xXG9iwN8II3yaGBz8GIDjibtvO8KYLlWYRERERAWDB5kL+vqAAA5Mns/+GgUW0249IthhgdzTbqTSLiIiICKF4krs/XoUF3NN2KY0LF2O6/YRO/KPd0dKCSrOIiIiI8OSX69haFKVTVoqLSl4BIDzgBkx/c5uTpQeVZhEREZEGbub6vby7eDsAL7X9DFdkF8mcDkSOu9LmZOlDs2eIiIiINGDF0QT3TlkFwLW9UnRY/1cAQsPvBKfXxmTfsTx+uyOoNIuIiIg0ZA9/tpZdwThtczO4Jv4khpkg1u5U4u1PtTtaKYezdH5mu2PYHUBERERE7PHZql18tGInDgOePn47mZu/xHK4CQ2fZHe0tKPSLCIiItIA7Q3H+dOnawD4xYBm9F7+IACR439JKrejndHSkoZniIiIiDQwlmXxp09WUxhJ0CXfz3WZU3AWbyTla0a4//V2xzuYmcK9bQ4AiRaDbVtkRWeaRURERBqYD1fs5Is1e3A5DO4/OZfsr58GIDT0NixPwOZ0hzLiIYx4yNYMKs0iIiIiDciOkhgPf1Y6LONXQ9vRZ9UTGMkwieb9iXW9wOZ06UulWURERKSBsCyLez9eRTCWoneLLH7RejsZq97FwiA44h4wDLsjpi2NaRYRERFpIN5ZvI3ZG/fhdTmYNLozOVN/AkC058Ukm/a1OV1605lmERERkQZgS2GEJ79cB8A1IzrQfcd7uHcvxfRkExpyi83p0p9Ks4iIiEg9lzIt7vroWyIJk/5tcriohw//7NIp5sKDbsTKbGJzwvSn4RkiIiIi9dzfvy5gUUExPreTiWd0IzDvPhzRfSQbdyPS+3K741XIcmXYHUGlWURERKQ+W7s7xHMz1gMwfmRH2iTXk7n0DQCCw+8Cp9vOeBVzOEm0GWF3Cg3PEBEREamv4kmTSR9+SzxlMaxDY87t3YzA9EkYVopYpzEk2gy3O2KdodIsIiIiUk/934wNfLszSE6Gi9tHd8G7bjKegplYTi/BoRPsjlenaHiGiIiISD00Z8M+3lywBYAJZ3Qj32sS+OoeAML9fouZ3cbOeJVnpnBvnw9AovkA25bRVmkWERERqWf2heNM+uhbAH58XAtO7twE39xHcQYLSAVaET7htzYnrBojVmx3BA3PEBEREalPLMvi7o9XsScUp0MTHzec3BFH8WZ8X/8fAMFhE8CdaXPKukelWURERKQe+eeibcxYtxeP0+C+sd3JcDsJzLwHIxUj3moo8U5j7Y5YJ6k0i4iIiNQTa3aHePLLtQBcf1JHuuQHcG+egXftZCzDSXDEXWAYNqesm1SaRUREROqBaCLFHR+sKJte7v+d0BJSCQIzJgEQ6X05qSY9bE5Zd6k0i4iIiNQDz0xfz9rdYRr73Ew8syuGYZC59A1ce7/FzGhEeNCNdkes0zR7hoiIiEgdN2PdHt5auBWASWd2o7HPgxHZg2/uowCEhtyClZFrY8JjY6XBqoUqzSIiIiJ12O5gjLs+WgXAT/u3YmiHxgD4Zz+II15MIq830R4X2xnx2DicJNqOtDuFhmeIiIiI1FWmZXHXR6sojCToku/nmuEdAHDtXEzG8r8DEDzpHtsWBKlPVJpFRERE6qi/Lyhg9sZ9eF0O7hvbA4/LAZZFYPoEDCyiXc8n2WKg3THrBQ3PEBEREamDvt0R5Jnp6wH4/ciOdGjiA8C76h3c2xdguXyETrzNzojVw0zh2rEQgGSzE7SMtoiIiIhUTiSR4o7JK0iaFiM7N+H8vi0AMOJB/DPvByA04HrMQAs7Y1YbR3Sf3RE0PENERESkrnn8i7Vs2BshP+Dh9tGl08sB+BY8hTO8g1R2OyLH/9LmlPWLSrOIiIhIHfLZ6t28u3g7BnDXWd3IzSydjs1ZuI7MRS8BlK785/TamLL+UWkWERERqSN2lMS4b0rp9HKXD2rDwLaNyu7zz7gLw0wQa3sK8Xan2hWx3lJpFhEREakDUqbFpA9XUhxN0rN5FlcPbVd2n2fDVLwbp2I53ISG3wn7h2tI9VFpFhEREakD/jJvMws2F5HpdnDvmO64nPtrXCqGf8adAESOu5JUo072hazHNHuGiIiISJr7pqCI57/aAMBNozrTplFm2X2Z37yCq2g9KV9TwgNusClhDTPsP8+r0iwiIiKSxvaF49z2/gpSFpzRPZ9xvZqV3ecIbcc3/0kAQifehuXJsitmzXE4ibe3f4y2/bVdRERERMplWhYTJ3/LzmCc9o0zue3076aXA/DP+hOORIhEs37Eul1gY9L6T6VZREREJE29NmdT2TLZD5zdE5/nu9XwXNvmk/Htv7EwCI64Oy2GMNRnGp4hIiIikobmbdrHizM3AnDraZ3plOf/7k4zRWD6RACiPS4k2ex4GxLWEsvEtWMRQOnHadMvByrNIiIiImlmdzDGHR+sxLTg3N7NGder+UH3Z6x8C/euxZieLEJDbrUpZS2xLByRPWW3sWk2PZ3HFxEREUkjSdPi9g9WsjecoEu+nz+MOngKOSNWhH/WAwCEB92I5cuzI2aDo9IsIiIikkZemrmBr7cU4XM7+dO4HmS4nQfd75v7GI7oXpKNuhDp/TObUjY8Ks0iIiIiaeKr9Xt5dc5mAG4f3YV2jX0H3e/c8y2ZS/4MQHDEXeB013bEBkulWURERCQNbC+OMmnySgB+cnxLRndvevAGlkVg+kQMK0Ws45kk2pxkQ8qGS6VZRERExGbJlMlt76+kKJqkR7MAvzu54yHbeNZNxlPwFZbTS3DYRBtSNmwqzSIiIiI2e2b6BpZsKybgdfKns3vgcf2goiUiBGbcDUD4hKsxs9vakLJh05RzIiIiIjb6cs1u3lywBYBJZ3SjVU7mIdv4Fj6HM1hAKtCScL9razuivRxO4h1OtzuFzjSLiIiI2GVLYYQ7P/oWgJ/2b8XILodOH+co3oLv62cBCA2dAO5DS7XUPJVmERERERvEkya3vb+CYCxFnxbZXDeiQ7nbBWbeg5GKEW85hFjncbWcUg7Q8AwRERERGzzx5TpW7AiSk+Hi/nHdcTkPPZfp3vIV3rUfYBkOgiPuBsOm5fDsZJm4di0BIJnfR8toi4iIiDQUU1bu5J+LtgJw15juNM/OOHQjM0lgeuksGdHel5HK61mbEdOHZeEI7Sy9nadltEVEREQahI17w9w3ZTUAPx/chmEdGpe7XcbSN3Dt/RYzoxGhQX+ozYhSDpVmERERkVoSTaT44/srCCdS9Gudw6+Gti93OyOyB//cRwEIDb4FK6NRLaaU8qg0i4iIiNQCy7K4d8oqVu8K0djn5r6x3XE5yh9r4J/9EI5YEYm8XkR7XlzLSaU8Ks0iIiIiteBvCwr4eOUunA6DP53dg7yAt9ztXLuWkLH8bwAER9wDDmdtxpTDUGkWERERqWFzN+7jqWnrAPj9yI70a51b/oaWRWDaBAwsol3OI9lyUO2FlCNSaRYRERGpQVuLotz2/gpMC8b1asZPjm952G29q97FvX0+lstHaOjttZhSKqIp50RERERqSDSR4qb3llEUTdKjWYBbT+uCcZi5lo14EP/M+wAIDbgeM9CiNqOmL4eTeLtRZbftotIsIiIiUgMOXPi3av+Ffw+d0xOv6/B/5PcteBpneAep7HZEjruqFpPWAWkwrlvDM0RERERqwJs/uPCv3AVM9nMWriNz0UsABIffCa7Dbyv20JlmERERkWo2Z+M+nq7MhX/7+b+6G8OME287knj702ohYR1imTh3rwAgldfDtmW0daZZREREpBptLYpyeyUv/APwbJiKd8OnWA4XweF3wWHGPDdYloUzuBVncCtYlm0xVJpFREREqklVLvwDIBXHP+NOACJ9ryTVqFPtBJUqU2kWERERqQZVvfAPIPObV3AVrcfMzCc88He1E1SOikqziIiISDWoyoV/AI7QDnzznwAgOPQ2LE9WLaSUo6XSLCIiInKMqnrhH4B/1p9wJEIkmp1ArNuPajihHCuVZhEREZFjUFAUqdKFfwCu7QvI+PZfAARH3GPbjBBSefoKiYiIiByl0gv/llf+wj8AyyQwbQIAkR4Xkmx2fM0HlWOmeZpFREREjsKBC/9WV+HCP4CMFW/h3rUY05NFaMittZC0jnM4ibc9uey2XVSaRURERI5CVS/8AzBiRfhnPwBAeODvsXz5NR2zfnB67E6g4RkiIiIiVTV7w94qX/gH4Jv3OI7IHpKNuhDpc0XNBZRqpzPNIiIiIlWwbk+IW/9XtQv/AJx7V5G5+DUAgiPuAqe7JmPWH5aJc++3AKQad9My2iIiIiLpbl84zvh3lxGKpzi+VTZ/rMyFfwCWRWD6RAwrRazDGSTanFTzYesLy8JZvAVn8RYtoy0iIiKS7uJJk5v/u5ytRVFa5WTw8Dm98FTiwj8Az/qP8GyZgeX0Ehw2sYaTSk1QaRYRERGpgGVZ3PfJKhYVFBPwOnn8/N7k+io5vCIZITDjbgDCJ1yNmdOuBpNKTVFpFhEREanAn+duZvLynTgNeGBcTzo08VX6sb6FL+As2Uwq0IJwv2tqMKXUJJVmERERkSP49Ntd/N+MDQDcdGpnBrdvVOnHOkoK8H39DAChoRPAXfmyLelFpVlERETkMJZtL+HOj0pnbrioXyt+dFzlZso4wD/zXoxklHjLwcQ6n10TEaWWqDSLiIiIlGN7cZQb/7OMWNJkeMfG/O7kjlV6vLtgJhlr/odlOAiOuAcqM8uGpC3N0ywiIiLyA+F4it//Zxl7QnE65/m5d2x3nI4qlF4zSWB66SwZ0V6XkcrrWUNJGwDDQbz18LLbdlFpFhEREfmelGlxxwcrWL0rRGOfm8fO74XfU7XKlLHsr7j2rMT05hIa/IcaStpAGAa4M+1OoeEZIiIiIt/39LT1TF+3F4/T4NHzetEiO6NKjzcie/HPeRiA0JCbsTIqf+GgpC+daRYRERHZ793F23hzwRYAJp3Zjd4tsqu8D/+ch3HEikg26Um05yXVHbHhsUyc+9YAkGrUWctoi4iIiNhp7sZ9PDi1tJz9emg7RndvWuV9uHYtJWPZXwEInnQ3OJzVmrFBsiycRRtxFm3UMtoiIiIidtqwN8yt/1tByrQ4o3s+Vw5pW/WdWBaB6RMxsIh2OZdEyyHVH1Rso9IsIiIiDVphJMH4d5dSEkvSt2U2E87ohnEU08N5V7+He9tcLFcmoRNvr4GkYieVZhEREWmwEimTm/+7nC2FUVpme3n43J54XUdRj+Ih/DPvBSDc/zrMrKotgiLpT6VZREREGiTLsrhvyioWbinC73Hy2Pm9aezzHNW+fF8/gzO0nVR2W8LH/6qak0o6UGkWERGRBunZGRv4YPlOnAbcP64HnfL8R7UfR9EGfAtfACA4bBK4qjZFndQNKs0iIiLS4Ly9sIDX524G4PbRXRnaofFR7ysw424MM068zcnEO4yuroiSZjRPs4iIiDQon63axSOfrQXgN8Pac3bv5ke9L/fGz/FumILlcBEcfmfp6nVSvQwHiVYnlt22i0qziIiINBgLtxQxYfJKLOBHx7Xg54PbHP3OUnECM+4EINLnF6Qad6mWjPIDhoHlCdidQsMzREREpGFYuzvEjf9ZRjxlMbJzE24a1fmoppY7IHPxa7gK12Jm5hEe+LvqCyppSWeaRUREpN7bXhzl+n8voSSW5LiW2dwzpjtOx9EXZiO0E9+8xwEInvhHLG/Vl9uWSrJMnIXrAUjldrBtiIZKs4iIiNRrxdEEN7yzlJ3BOB0a+3j0vF5kuI9teevA7AdwJIIkmh5HrPtPqimplMuycBauAyCV0x5sGjau4RkiIiJSb8WSJn94bznr9oTJD3h46ke9ycl0H9M+Xdu/JmPl2wAER9xj68VpUnv0VRYREZF6KWVaTPpwZdniJU9d0Ifm2cc4h7JlEpg+AYBo9/9Hsnm/akgqdYFKs4iIiNQ7lmXx+BdrmbpqN26nwSPn9qJz/tEtXvJ93pX/xL3zG0x3gOCQW6shqdQVKs0iIiJS77wxbwtvLdwKwJ1ndmNA29xj3qcRKyYw608AhAeOx/I3PeZ9St2h0iwiIiL1yuTlO3hmeulsC+NHdmR09+opt755T+CI7CaZ24lI359Xyz6l7lBpFhERkXpj1oa93P3xKgAuHdCan/ZvXS37de5dTeaSVwFKV/5zeqplv1J3aMo5ERERqRdW7Cjhlv8uJ2VanNmjKded1KF6dmxZBGbciWEmibUfTaLdKdWzX6kcw0Gi5aCy23ZRaRYREZE6b0thhN+9s5RIwmRQ21wmntEVxzGs9vd9nvVT8Gz+EsvhITh8YrXsU6rAMLC8OXan0PAMERERqdt2B2Nc9+8l7A0n6Jrv58FzeuJ2VlPFSUYJfHUXAOETfo2Z07569it1js40i4iISJ1VGElwzb+WsKUwSsucDJ68oDcBb/XVG9+iF3EWbyLlb06437XVtl+pAsvEUbwJADO7rZbRFhEREamKYCzJ9f9eUrba37M/7kNewFtt+3eUbMW34GkAQkPvAM+xz/MsR8GycO1dDUA8q42W0RYRERGprGgixe//s4wVO4LkZrp59sd9aZ2bWa3H8M+6DyMZIdFiELEu51brvqXuUWkWERGROiWRMrn5v8tZuKWIgNfJMz/qQ4cmvmo9hnvrbDJWv4dlOCgZcQ9U00WFUnepNIuIiEidkTQt7vhgJbM27CPD5eCJ83vTrVmgeg9iJglMmwBAtOclpPJ7Ve/+pU5SaRYREZE6wbQs7p2yis9W78btNHjk3F4c16r6pyLLWP43XHtWYHpzCA2+qdr3L3WTSrOIiIikPcuyePSztXywbAdOA+4f24PB7RtV+3GM6D78sx8CIDT4JqzMxtV+DKmbVJpFREQk7T331QbeXrQVA5h4ZjdGdsmrkeP45zyCI1ZIskl3or0urZFjSN2kKedEREQkrf15ziZem7MZgFtO68yYns1q5DjO3cvJWPYXAIIj7gaHalJaMBwkWgwou20XPRtEREQkbf1z0VaenbEBgOtP6sCPjmtZMweyLALTJ2BYJtHOZ5NoNbRmjiNVZxhYGdU/FKeqNDxDRERE0tIHy3bw0NQ1AFw5pC2XDWxTY8fyrvkfnq1zsFwZpQuZiPyAzjSLiIhI2vls9W7u/vhbAC7q14pfD21XcwdLhPHPvAeAcL9rMbNa1dyxpOosE0dJAUDp10bLaIuIiIjArA17uf39FZgWnNO7GeNHdsSowcVFfAuewRncRiqrDeETfl1jx5GjZFm49qwEIB5oqWW0RURERBZuKeKm95aTNC1O65rHbad3xVGDhdlRtAHfwucBCA6fCK7qXYpb6g+VZhEREUkLS7YWM/7dpcSSJsM6NObuMd1xOmr2tGLgq3swzDjx1iOIdzizRo8ldZtKs4iIiNhu8dZirvv3EkLxFP3b5PDA2T1wO2u2prg3fYF3/cdYhpPgiLugBs9oS92n0iwiIiK2Wry1mOu/V5gfP783GW5nzR40FScwfRIAkb4/J9W4a80eT+o8lWYRERGxzTcFRWWFecD+wpxZ04UZyFzyZ1yFazEzmxAeOL7Gjyd1n0qziIiI2OKbgiJueGdpWWF+rJYKsxHaiW/uYwCEhtyK5c2p8WNK3acp50RERKTWlZ5hXko48d0Z5hofkrGff/aDOBJBEk2PI9rjwlo5phwDw0Gi2fFlt+2i0iwiIiK1ys7C7NqxkMyVbwEQHHG3rSVMKskwsHz5dqdQaRYREZHac1BhbpvL4+f1qrXCjGUSmDYBgGi3H5Ns3r92jiv1gkqziIiI1ApbCzPgXfkv3DsXYbr9hE78Y60dV46RZeIIbgfADDTXMtoiIiJSfy3aUnrRXziRYmDbXB6r5cJsxIoJzPoTAOEBv8P0N6u1Y8sxsixcu5cBEPc3s20ZbZVmERERqVGLthRx/TtLiCRMWwozgG/+kzgiu0jmdiRy3JW1emypHzT6XURERGrM9wvzIJsKs3PfGjIXvwJAaPid4PTU6vGlflBpFhERkRqx8AeF+VEbCjOWRWD6JAwzSaz9acTbjard40u9odIsIiIi1W7hliJusLswA54Nn+DZ/CWWw0Nw2KRaP77UHyrNIiIiUq3mbypMi8JMMkpgxl0ARI7/JWZuh9rPIPWGLgQUERGRajNt7R7++L/lxFOWvYUZ8C16CWfxRlL+ZoT6X29LBqk/VJpFRESkWkxZuZOJH35LyrQ4qVMT7h/XA6/Lnj9qO4Jb8S14CoDQibeDx29LDqkGhoNk0z5lt+2i0iwiIiLH7J3F23jgk9VYwJk9mjLpjK64nPYVHP/M+zCSERLNBxDrer5tOaQaGAamv7ndKVSaRURE5Nj8Zd5mnpq2HoAfHdeCm0/tjMOwaQUKwL11Dhmr38PCIHjSPWBjFqk/VJpFRETkqFiWxfMzN/Lq7E0AXD6wNdeO6IBhZ0k1UwSmTQAg2vOnJPP72JdFqodl4QjvAMD0NbPtlyCVZhEREaky07J47PO1vLVwKwC/Hd6enw9ua3MqyFj+N1x7lmN6cwgNudnuOFIdLBPXziUAxNvlg2HPhaUqzSIiIlIlSdPivimreH9Z6dm/m0Z15v+d0NLmVGBE9+Gf/SAAoUE3YmU2sTmR1CcqzSIiIlJp8aTJhMkr+Wz1bpwGTDyzG2N6NrM7FgD+uY/giBWSbNyNaO/L7Y4j9YxKs4iIiFRKNJHipv8uZ/aGfbidBveN7cEpXfLsjgWAc/dyMpb+BYDgiLvBoYoj1UvPKBEREalQMJZk/LtLWVRQTIbLwSPn9mJw+0Z2xyplWQSmT8SwTGKdxpJoPczuRFIPqTSLiIjIEe0Lx7n+30tZuTNIwOvkifN7c1yrHLtjlfGueR/P1tlYTi/BoRPsjiP1lEqziIiIHNbOkhjX/GsxG/ZGaJTp5ukf96Fb04Ddsb6TCOOfeQ8A4X7XYGa3tjmQ1FcqzSIiIlKuLYURrvnnYrYWx2ga8PDsj/vSvonP7lgH8X39LM7gVlJZrQn3+43dcaQmGAbJvF5lt+2i0iwiIiKHWLa9hPHvLGVfJEGb3Aye/UlfWmRn2B3rII7iTfgWPg9AcNgEcGXanEhqhOHAzLJ/SkOVZhERETnIV+v2cuv/lhNNmnRrGuCJC3qT5/fYHesQga/uxkjFiLcaRrzjGLvjSD2n0iwiIiJl3luyjT99spqUBUPaNeKBc3rg96RfXXBvnoZ33UdYhpPgiLts/bO91DDLwojsLr2ZmadltEVERMQ+lmXx0qyNvDRrEwBjezbljtFdcTkdNicrRypBYPokACJ9fkaqSXebA0mNskzcOxYBEG83Sstoi4iIiD2SpsUDn6zmvaXbAfjF4DZcPaw9Rpqevc1c+jqufasxMxoTHvh7u+NIA6HSLCIi0oCF4ylue38FX63fi8OAm0/tzI+Os/+iq8MxwrvxzX0UgNCQW7Aycu0NJA2GSrOIiEgDtTcc53fvLGXFjiBel4P7xvbg5M5N7I51RP7ZD+CIl5DI70O0x0V2x5EGRKVZRESkAdq8L8L17yxhS2GUnAwXj5/fmz4ts+2OdUSuHYvIWPEWAMERd4PDnrGt0jCpNIuIiDQwS7cVM/7dZRRGErTMyeCpC3rTrnF6LVpyCMskMH0iBhbRrheQbDHQ7kTSwKg0i4iINCDT1+7hj++vIJY06dEswOPn96ZJGs7B/EPeb9/BveNrTLef0NDb7I4jDZBKs4iISAPxzuJtPPjpakwLhnZoxJ/G9cTnSf8hDka8BP+s+wEID7gB09/c5kRSqwyD5IFpBbWMtoiIiNQUy7J4YeZGXpldOgfz2b2acdvpXdJzDuZy+OY/iTO8k2ROByLHXWl3HKlthgMzu43dKVSaRURE6rN40uT+T1fzwbIdAFw1pC2/Gtoubedg/iHnvrVkfvMKAKHhd4LTa28gabBUmkVEROqpfeE4N/93OYsKinEYcOtpXTi/bwu7Y1WeZRGYMQnDTBBrdyrx9qfanUjsYFkYscLSm95cLaMtIiIi1WfN7hA3vruUrcUxAl4nfxrXgyHtG9sdq0o8G6fi2fQFlsNNaPgku+OIXSwT97b5gJbRFhERkWo0Y90ebn9/JeFEita5GTx+Xm/aN0nzKeV+KBUjML20KEeO/yWp3I42B5KGTqVZRESknrAsi78tKODJL9dhAf3b5PDA2T3JzXTbHa3KMhe9hLN4IylfM8L9r7c7johKs4iISH2QSJk8OHUN7y3ZDsB5fZpz86mdcdeRGTK+zxHchn/+UwCEht6G5QnYnEhEpVlERKTOKwwnuOV/y/l6SxEOA244uSMX92tVZ2bI+CH/rPsxkmESzfsT63qB3XFEAJVmERGROm3dnhC/f3cZBUVR/B4n943rwbAOdeuCv+9zbZtHxqp3sTAIjrjH1sUsRL5PpVlERKSOmrl+L7e9v4JQPEXLnAweO68XnfL8dsc6emaKwLQJAER7XkyyaV+bA4l8R6VZRESkjrEsi7cXbuWxL9ZiWnBCq2weOqcXub66d8Hf92Ws+Dvu3UsxPdmEhtxidxxJF4ZBsnGXstt2UWkWERGpQ5Ipk4c/W8s7i7cBpUti33paFzyuunfB3/cZ0UL8sx8EIDzoRqzMJjYnkrRhODBz2tudQqVZRESkriiKJLj1/RXM31SIAVx3UgcuHdC6zl7w932+uY/iiO4j2bgbkd6X2x1H5BAqzSIiInXAmt0hbn5vGZsLo/jcTu4d250RnerH2VjnnhVkLn0DgODwu8BZt4eZSDWzLIx4celNT7aW0RYREZHyTVm5k3s+XkU0adIi28uj5/WiS349mbvYsghMn4RhpYh1GkOizXC7E0m6sUzcW+cCWkZbREREypFMmTw9fT1/W1AAwKC2udw3tkedv+Dv+zxrP8BTMBPL6SU4dILdcUQOS6VZREQkDe0Jxbnt/RV8vaUIgCsGteHqYe1xOur++OUyiQiBr+4BINzvt5jZbWwOJHJ4Ks0iIiJpZsnWYm7533J2BeP4PU4mndmNU7rk2R2r2vkW/h/OYAGpQCvCJ/zW7jgiR6TSLCIikiYsy+Lf32zj0c/XkjQtOjT28dA5PWnfxGd3tGrnKN6M7+v/AyA4bAK4M21OJHJkKs0iIiJpIJpI8cDUNXywbAcAp3bNY8IZXfF76ueP6sDMezBSMeKthhLvNNbuOCIVqp+vRBERkTqkoCjCLf9dwbc7gzgMuHZE/Zl/uTzuzTPwrp2MZTgJjrjL1lXeRCpLpVlERMRGszfs5Y4PVlIUTZKb6eb+cd0Z2LaR3bFqTipBYMYkACK9LyfVpIfNgSTtGQap3I5lt+2i0iwiImID07L485zNPP/VBiygZ/MsHjy7B82zM+yOVqMyl76Ba++3mBmNCA+60e44UhcYDlKNOtmdQqVZRESktgVjSe788Fu+XLsHgPP6NOcPozrjdTlsTlazjMgefHMfBSA05BasjFx7A4lUgUqziIhILVq9K8it/1vBpn0RPE6Dm0/tzLl9Wtgdq1b4Zz+II15MIq830R4X2x1H6grLwkiESm+6/VpGW0REpD47MJ3c41+sJZ6yaJbl5aFzetKzeZbd0WqFa+diMpb/HYDgSfeAw56lkKUOskzcBbMALaMtIiJSr5VEk9w7ZRWfrd4NwPCOjZl0Rrd6tRz2EVkWgekTMLCIdj2fZIuBdicSqTKVZhERkRq0dFsxt7+/gq3FMVwOg2tHdOCn/VvV2+nkyuNd9Q7u7QuwXD5CJ95mdxyRo6LSLCIiUgNMy+LN+Vt4dsYGUqZFy5wM7h/Xg14NZDjGAUY8iH/m/QCEBlyPGWgY47el/lFpFhERORrJKN617+Nd9zFGtBArI5dYxzOIdRrHvriDOz/6lpnr9wFwWtd8bh/dhYC34f3Y9c1/Emd4B6nsdkSO/6XdcUSOWsN79YqIiBwjz/opZE0djyNWhIUDAxMLB951H5L4ciJ/Sv2WmeG+eF0Ofj+yI+f3bdGghmMc4CxcR+Y3LwOUrvzn9NqcSOToqTSLiIhUgWf9FLInX1n2fwPzoLfOeDGP8CAZ2bdy9nlX0Dnfb0vOdOCfcSeGmSDW9hTi7U61O47IManfs6iLiIhUp2SUrKnjATCwyt3EYQAG3MP/0blRw51WzbNhKt6Nn2E53ISG32nr8sdSxxkGqZx2pHLa2fo8UmkWERGpJO/a93HEig5bmA9wYOGMF+Fd+0EtJUszqRj+GZMAiBx3ZVosgSx1mOEg1bgrqcZdwbCvuqo0i4iIVJJ33cdYlfzRWTrG+aMaTpSeMr95GVfRBlK+poQH3GB3HJFqoTHNIiIilWREC8vGLle4LSZGtLBmA6UhR2g7/nlPAhA68TYsT8OaYk9qgGVBMlp625Vh2xANnWkWERGpJDMjF5PK/cC2cGBl5NZsoDTkn3k/RjJMolk/Yt0usDuO1AeWiWfLDDxbZoBVuV9aa4JKs4iISCXs2rObDZs24KhgPPMBBiaxjmfWcKr04to2j4xV72BhEBxxt63jT0Wqm57NIiIiFVgw+1Oy/z6a7skVWBYV1mYLA9ObQ6zT2FrJlxbMFIHpEwGI9riQZLPj7c0jUs1UmkVERA6jKBRlxpsTOG3+L2ht7GS7oymbT7iZ0gnnyh+mceD9Jac+UTr+soHIWPEP3LuWYHqyCA251e44ItVOFwKKiIiU45sVK8j9/Hecby0DA5Y1Op0m5z1Opi+X4pbdyZo6HiNWhIWBgfXdW282Jac+QbzD6XZ/CLXGiBbin/0gAOFBN2L58mxOJFL9VJpFRES+J5pI8fkHb3DOlgfINUKEyWBj/0k0HXxp2VX78Q6j2XPFArxr/kfm0r9AIoSZ04FYp7NKh2Q0oDPMAL55j+GI7iXZqAuR3j+zO45IjVBpFhER2W/llp3s++BWLk9OAQM2Z3TDde6L5OWVsziHK4NY1wuwvLkAxNuNAkfDWwHQuWclmUteByA44i5wum1OJFIzVJpFRKTBS5oWkz+fykkrbmeEUYCJwdpOvyD39NvB6Tn8Aw2DVHbrstsNjmURmD4Jw0oR63gmiTYn2Z1I6qM0eZ2pNIuISIO2aW+Y+f95mCvCr+E1khQ6GxM8/SlyO42s+MGGg1STHjWeMV151k3GU/AVltNLcNhEu+NIfZUmrzOVZhERaZCSpsV/Zi+mx4I7+LVjYelwjLyTyTj7Kby+JnbHS3+JCIEZdwMQPuFqzOy2NgcSqVkqzSIi0uCs3hXkg/f/we9Cj5PvKCKOh12DbiNjwJVV//NvKl769kjDOOoh38LncAYLSAVaEu53rd1xpL5Lg9eZSrOIiDQY8aTJG7PWkL/wEe5zTgYD9vo7YY17Hk/eUfz510zh2fRl6b4b0IWAjuIt+L5+FoDQ0AngzrQ5kdRrafI6U2kWEZEGYem2Yv784efcFHqY3s4NAOztfhmpkyeCS6WvKgIz78FIxYi3HEKs8zi744jUCpVmERGp16KJFM/NWI/5zZs853oDnyNGzJ1L5PTHSHUYbXe8Ose95Su8az/AMhwER9zdMGcNkQZJpVlEROqt+ZsKefLjr7ku8gxj3XMBCLccRmT0k5j+5janq4PMJIHppbNkRHtfRiqvp82BRGqPSrOIiNQ7wViSp6atY8uSL3jF8yytnHswDRfhIbcQOeHXYDjsjlgnZSx9A9febzG9uYQG/cHuOCK1SqVZRETqlelr9/DwJyu5KPYWj3jexWlYJLLbEzzjWZJNj7M7Xp1lRPbgn/soAKEht2BlNLI5kUjtUmkWEZF6YV84zqOfr2XZt8t52v0sA1yrAIh0v7B07K3Hb3PCus0/+yEcsSKSTXoS7flTu+OI1DqVZhERqdNMy+L9pTt4ato6TopP40PPK2QZEUx3FsFTHiDW5dyaO7hhkAq0LLtdX7l2LSFj+d8ACJ50T4OZWk/SRJq8zlSaRUSkzlq9K8gDn65h7dYd3Ol6nZ94pgGQaD6A4tOfxsxuU7MBDAep/F41ewy7WRaBaRMwsIh2OZdEy8F2J5KGJk1eZyrNIiJS54TiSV6cuZG3vi6gJ+uY7H2GdsZ2LMNBuP/1hAf+Dhz6EVcdvKvexb19PpYrk9DQ2+2OI2IbfUcREZE6w7Ispq7azWNfrGV3MMqvnB9wk/ttXKRIBVpScvrTtX8m1EyVvq2HQxaMeBD/zPsACPe/HvPAn8hFalsavM5UmkVEpE7YvC/CQ5+tYfaGfeSzj7czX2CgtRiAaKdxBEc+gJWRW7uhzBSejZ8B9XMZbd+Cp3GGd5DKbkf4+F/aHUcaqjR5nak0i4hIWoslTV6fu4nX524mnrI4w/U1j3lfxp8qxHJlEhxxN9EeF9XrC/Hs4CxcR+ailwAIDp8ErgybE4nYS6VZRETS1sz1e3n4szVsKYziJc4Ljf7NGZH/QQoSeb0pGf0sqUad7I5ZL/m/uhvDjBNvezLx9qfbHUfEdirNIiKSdnaUxHj8i7VMXbUbgCH+7Tyf8X/khtYAED7+14SG3AxOr50x6y3Phql4N3yK5XARHH6XzuKLoNIsIiJpJJEyeWvhVl6cuYFIwsRpWDzWfh5n73wORyiGmZlP8WmPk2g70u6o9Vcqjn/GnQBE+l5JqlFne/OIpAmVZhERsZ1lWcxYt5cnvlzHpn0RAIY1t3ja9wqNt5ZeABRrN4qSUY9h+fLsjFrvZX7zCq6i9ZiZ+aVT94kIoNIsIiI2W7s7xONfrGXOxkIAGvvc3NtrF2esvRNn4U4sh4fQ0NuJ9P2FhgnUMEdoB775TwAQPPGPWJ4sewOJpBGVZhERsUVhOMELMzfwzuJtmBa4nQaXHt+MG5xvk7P4BQwsko26UDz6WVJ5Pe2OWz7DwPQ3Lbtd1/ln/QlHIkSi6fHEuv/Y7jgipdLkdabSLCIitSqRMvnnoq28NGsjwVjpggWjuuRx4wkOus4ej3vnNwBEel1KcNgkcGfaGffIDAfJpsfZnaJauLYvIOPbfwEQPOkeMBw2JxLZL01eZyrNIiJSKyzLYvq6vTz5vXHLXfP9/H5kR4aGPyXw4R04EiFMby4lox4m3vEsmxM3IJZJYNoEACLdLyTZ7ASbA4mkH5VmERGpcWt2h3jiB+OWfzu8PWd39pEz/TYyVr8HQLzViZSc9qSWa65lGSvewr1rMaYni9CJt9odRyQtqTSLiEiNKQwneH7mBt793rjli/u15ueD25C77xuy/3ktzpLNWIaT8KA/EO7327q1FHWaLO97LIxYEf7ZDwAQHjgey5dvcyKRH0iT15lKs4iIVLvDjVu+7qQOtM724FvwNL55j2NYKVLZbSk+/RmSzfvZnLph8s17HEdkD8lGnYn0ucLuOCJpS6VZRESqjWlZfLRiJy98tYGtxTFg/7jlUzrRv00ujpICsv5zPZ5tcwCIdj2f4Mn3a2ozmzj3riJz8WsApSv/OT02JxJJXyrNIiJyzCzLYub6fTw7Yz2rd4UAaOL3cPXQdpzduzlOh4Fn7QdkfX4zjlgRpjtA8OT7iHX7kc3JGzDLIjB9IoaVItbhDBJtT7Y7kUhaU2kWEZFjsnhrMc9MW8fCgmIAAl4nlw9sw0X9WpHpdkIiTODLO8lc/jcAEk2Pp3j0M5g57W1MLZ71H+HZMgPL6SU4bKLdcUTSnkqziIgclXV7Qvzf9A18uXYPAB6nwYUntOJng9qQk+kGwLVrKVlTrsFVuBYLg3D/awkP/D043XZGl2SEwIy7AQgf/2vMnHY2BxJJfyrNIiJSJduLo7wwcyOTl+/AtMBhwNm9mvPLoe1oluUt3cgyyfzmFfyz/oRhxkn5m1Ny2pMkWg+zN7wA4Fv4As6SzaT8zQn3v9buOCJ1gkqziIhUSmEkwWtzNvGvRVuJpywARnZuwm+Hd6BDE1/ZdkZ4F9lTx+PZ9AUAsQ5nUDLqEayMRnbErlmGgZnZpOx2XeAoKcD39TMAhIZNALevgkeI2CxNXmcqzSIickSRRIq/LyjgjXmbCcVLp4/r3yaHa0d0oHeL7IO2dW/8nOyp43FEdmO5MggOu5Nor0vqTKGsMsNR56bK88+8FyMZJd5iMLHO59gdR6RiafI6U2kWEZFyxZIm/1m8jVfnbGJvOAGUTh93zYgOnNi+Ecb3i3Aqhn/Wn/B98zIAySY9KB79LKnGXe2ILofhLphJxpr/YRkOgiPurr+/zIjUAJVmERE5yIGy/Pq8zewKxgFolZPBb4a15/Tu+Th+ULSce1eTPeUaXHuWAxDu+wtCJ94Groxazy5HYCYJTC+dJSPa61JS+b1sDiRSt6g0i4gIANFEiv8s2c7rczezO1RalptlebliUBvO7dMct9Nx8AMsi4zlbxKYcSdGMoqZ2YSSUY8Rb3+qDeltYqbKxm7H245M62W0M5b9FdeelZjeHEKDb7I7jkjlpcnrTKVZRKSBiyZSvLtkO2/8oCz/YnAbxvVqjsflOOQxRnQfWZ/fjHfdhwDE25xEyamPY/qb1Wr2tGCZdieokBHZi3/OwwCEBt9cPy/KlPotDV5nKs0iIg3UgbL8+tzN7Nlflptnefn5kLac3avZoWeW93MXzCLr0+txBrdhOdyEhtxK5PhfglH+9mI//5yHccSKSDbpUXphpohUmUqziEgDE02keGfxNt6Yt6VKZZlUAt/8J/DNfwoDi2RuR0pGP0syv08tppeqcu1aSsayvwKUXvzn0I9+kaOhV46ISANRXlluke3l54PbMu5IZRlwFG0k+5PrcO/4GoBIj4sIDr8LPP5ayS5HybIITJ+IgUW08zkkWp1odyKROkulWUSkngvFk/xn8Xb+Mv/gsvyLwW0ZW0FZBvCuepfAF3/EkQhierIJjnyQWJezayO6HCPv6vdwb5uL5cogNPQOu+OI1GkqzSIi9dTecJy3vi7gn4u2URJLAtBy/5nlypRlIx4kMO12Mr79NwCJFgMpPu1pzOzWNZ5dqkE8hH/mvQCE+1+HmdXS5kAidZtKs4hIPbOlMMJf52/h/WU7iCVLrzhv1yiTywa2ZkzPissygGvHQrKnXIuzeCOW4SA8cDzh/tdpPGw5zDSdicL39TM4Q9tJZbclfPyv7Y4jckzS4XWm734iIvXEtzuCvD5vM1NX7cK0St/Xq3kWPxvUhpM7NzlkUZJymSkyFz6Hf+4jGGaSVFZrik9/mmSLgTUbvq5yOEm2GGB3ikM4ijbgW/gCAMFhE7XQjNRtafI6U2kWEanDLMti/uZC3pi7hdkb95W9/8T2jfjZoDb0a51z8HLXR+AIbiPr0xvwFMwEINr5HIIj/4TlzamR7FJzAjPuxjDjxNucRLzDGXbHEakXVJpFROqglGnxxZrdvD53Myt2BAFwGHB6t3wuG9iGbk0DVdqfZ93HZH12I45YIZbLR8lJ9xLr/hOoZOGW9OHe+DneDVOwHK7SGU70NRSpFirNIiJ1SCxpMnn5Dv46fwub9kUA8LocnNO7OZcMaEWrnMyq7TAZIfDVPWQufQOARH5fSkY/Qyq3Y3VHr5/MFO4t0wFItB5h/zLaqTiBGXcCEOnzc1KNu9ibR6Q6pMnrTKVZRKQO2BOK88432/jXN1vZG04AkJ3h4ifHt+TCE1rSyOep8j6du5eTPeVaXPtWARA+4WpCg28GZ9X31ZAZqYTdEcpkLn4NV+FazMw8wgPH2x1HpNqkw+tMpVlEJI2t3FHCP74uYMq3u0ikSq/uaxrwcMmA1pzXpwU+z1GccbEsMpa8RmDmfRipGClfU0pOe5JEmxHVnF5qkxHaiW/e4wCEhtyK5c22OZFI/aLSLCKSZpKmxbQ1u/nH1wUsLCgue3/vFllcdEIrTu2ah6sS08aVx4jsIWvq7/FunApArP1plIx6FCuzSbVkF/sEZj+AIxEk0fQ4oj3+n91xROodlWYRkTRRHE3w3pLtvL1wK9tLYgA4HQandc3jon6t6N3i2M4cujdPI+vT3+EM78RyegkOvYNonyt0oVg94Nr+NRkr3wYgOOJuMI7ulyoROTyVZhERm63fE+athQV8sGwH0f2LkeRmurmgb3N+dFxLmmZ5j+0AqTj+2Q/iW1Q6b2+yUVeKRz9DKq/nsUaXdGCZBKZPACDa/Sckm/e3OZBI/aTSLCJiA9OymLV+H//4uuCg+ZW75Pu56IRWjO6eT4b72K8QdxauI2vKtbh3LQYg0vtnBIfdAa4qzrIhacu78p+4d36D6Q4QGnKr3XFE6i2VZhGRWlQUSfDB8h38+5ttZVPGGcDJnZtwUb9WVVqM5IgsC+/Kt8maNgEjGcb05lIy6lHiHbXQRXWz84I7I1ZMYNafAAgP/B2mv5ltWURqUjpc2KrSLCJSwyzL4puCYt5ZvI2pq3YR3z8LRsDr5Jzezfl/J7Ss+vzKR2DEigh8cSsZa/4HQLzVMEpOewIz0KLajiH7OZwkWg627fC+eU/giOwmmduJSN9f2JZDpEbZ/Do7QKVZRKSGlESTTF6+g3cWb2PdnnDZ+7s1DXBB3+ac0aMpfk/1fht2bZtH9pRrcQYLsBwuQoP+QOSE39i/6IZUO+fe1WQueRWA4PA7Nb+2SA1TaRYRqUaWZbFsewn//mYbn3y7i9j+C/u8LgdndM/ngr4t6Nk8q3qGYHyfmcQ3/yl885/AsExS2e0oHv0MyWYnVO9xJD1YFoHpEzHMJLH2o0m0O8XuRCL1nkqziEg1CMaSfLRiJ+8s3sbqXaGy93fK83FB3xac1aMZWRk18y3XUbyF7E+vx71tLgDRbj8meNK9WJ5AjRxPvsdM4S6YCUCi1dBaO6PvWf8xni3TsRwegsMn1soxRWxj0+vsh1SaRUSOwYodJbzzzTY+XrmTSKL0rLLHaXB6t3zO79uCvi2zq/+s8vd4V/+PwBe34IgXY3qyCJ58P7Gu59fY8eRQRjJauwdMRgjMuAuA8Am/xsxpX7vHF7FBrb/OyqHSLCJSRfvCcT5auYv3l25n1ffOKrdvnMn5fVswtmczcjLdNRsiHiIwYyKZK94CINGsH8Wjn8HMbluzxxXb+Ra9iLNkMyl/c8L9rrU7jkiDodIsIlIJyZTJV+v38v6yHUxft5eUWToDhttpMKpLHuf3bVF908VVwLVrCVlTrsFVuA4Lg/CA6wkP+B04a7ioi+0cJVvxLXgagNDQO8DjtzmRSMOh0iwicgSrdwV5f9kOPly+k32RRNn7ezQLMK5Xc0Z3zye3ps8qH2CZZC56Ef/sBzHMBKlAC0pOe4pEqxNr5/hiO//MezGSURItBhHrcq7dcUQaFJVmEZEfKAwn+GjlTt5ftoNvdwbL3t/Y52ZMz2aM7dWMznm1e4bPEdpB1tTxeDZPAyDW8SxKTnkIK6NRreYQ+7gLZpGx5r9YhoOSEfdALfxVQ0S+o9IsIsKB4Rf7eH/Zdmas20vye8MvTurUhHG9mjGkfWNcjtovKp4NU8maOh5HdC+WK4Pg8LuI9vypSlNDYiYJTC+dJSPa8xJS+b1sDiTS8Kg0i0iDZVkWS7eV8PHKnUxZucv+4Rc/lIzin3kfviWvAZDI60XJ6GdJNepsTx4pl1UL44ozlr2Ja88KTG8OocE31fjxRNJNbbzOKqLSLCINzprdIabsL8oFRd9NY2Tn8Isfcu5dRfaU3+LasxKA8HFXETrxj+D02ppLfsDhLJ03tgYZ0X345zwEQGjwTViZjWv0eCJppxZeZ5Wh0iwiDUJBUYQpK3fx8cqdrN393ZLWmW4HJ3fO44zu+bYNvziIZZGx7K8EZtyJkYphZuZRcupjxNuNsjeX2MY/52EcsSKSTboT7XWp3XFEGiyVZhGpG5JRvGvfx7vuY4xoIVZGLrGOZxDrNA5cGeU+ZHcozqff7mLKyp0s2VZS9n6302Bo+8aM7p7PSZ2akOG2Z3WpHzKi+8j67A94138MQLztSIpPfRzLl29zMrGLc9cyMpb9FYDgiLvBoR/bInYxLMuyKrPhrl0lFW8kIlIDPOunlF4IFyvCwoGBWfbW9OZQcuoTxDucDkBxNMHnq3fz8cpdLNhcyP7r+XAYMKBNLmd0b8opXfJqbEnro+Xe8hVZn16PM7QDy+EhNPQ2In1/AYbD7mhyJGYK97Y5ACRaDK7e5X0ti5x3f4xn2xyinc+m5Iznqm/fInVJTb7O9svPz6pwm/T6qSEi8gOe9VPInnxl2f8NzIPfxorJnvwLPuv9KK/v68WsDXtJpL47F9CnRRajuzfltG755Pk9tRu+MlIJ/HMfJfPrZzGwSOZ2omT0syTze9udTCrJiIcq3ugoeNf8F8+2OViujNKFTEQasJp6nVWFSrOIpK9klKyp4wEwKP+PYgYWpgX9l0zgt7FnSeChc56f0d3zGd09n1Y5mbWZuEocRRvInnIt7p2LAIj0/CnB4XeC22drLkkDiTD+r+4BINzvWsysVjYHEhGVZhFJW9617+OIFVW4ncOAXEI81G0tTYdcavvMF5Xh/fZfBL68HUciVDrE5JSHiHcaa3csSRO+Bc/gDG0nldWG8Am/tjuOiKDSLCJpzLvuYyyMw55l/j4LB2c651Ocl94Fw4iXEPjyNjJWvQtAvOVgSk57GjOrpc3JJF04ijbgW/g8AMHhE8GVvn8tEWlIVJpFJK1YlsWWDd+yb9G7DNn6Bd5KFGYoHeNsRAtrNtwxcm1fQPYn1+Es3oRlOAkP+j3hftfWyEUtUncFvroHw4wTbz2CeIcz7Y4jIvupNIuI7eJJk1Urvia2/L903PM5/az1Vd6HhQMrI7f6w1UHM4Xv62fxzX0Uw0qRympD8ehnSDbvb3cySTPuTV/gXf8xluEkOOIuLZUukkZUmkXEFnuCMVYumYVzzQf0LPqSU4ytZfclLQffevviadSKLjs+qNT+DExiHdPvrJwjuJWsT2/AUzALgGiXcwme/Ccsb7bNyaS6WIeZJ7zKUnEC0ycBEOn7c1KNu1bPfkXqgWp7nR0DlWYRqRWWZbFqZwnrl0wjsOEjBkS+4jzHrtI7DYjjYrVvAJEOZ9K8/7k0zcqHZBTzzzMwYsVHHNdsYWB5s4ml2YV0nrWTyfr8JhyxIky3n+BJ9xHr9iOdPaxPHE4SbUZUy64yl/wZV+FazMwmhAeOr5Z9itQL1fg6OxYqzSJSY6KJFF9v3MPWZZ+TVzCFk8w5DDcKS+90QBQvG3JPxOo6liZ9xtA0I+fgHbgyKDn1CbIn/+KwFwRalBbQklOfOOzKgLUuESEw404yl79Z+t+mx1F8+jOYuR1sDibpygjtxDf3MQBCQ27F8uZU8AgRqW0qzSJSbSzLYuPeCHPX7yC06nM67vmMUcZ8GhvB0g0MCBs+CvJOIqPXOWR0PZ0m7iPPDBDvcDrFY14ha+p4jHJWBLS82QetCGg35+7lZE/5La59a7AwiPT7LaFBN4IzDRdWkbThn/0gjkSQRNPjiPa40O44IlIOLaMtIsekOJpg3qZC5q/bhmvD5wyJf8Uox0KyjUjZNiFnNrtanIq/97nQfgQ4vVU/UDKKd83/yFj+D4x4EDO7DbFOZ5UOyUiHM8yWRebiV/DPvB/DjJPyNaPktCdJtBludzKpSWYK9/b5ACSaDziqmVBcOxbS6F9nA7DvR+/pAlGRH6qG11lFtIy2iFS7pGmxbFsxszfsY8mGLbTYNZ0zHHO5x/ENmUYc9n8vC7rzKGk3moxe55JsORi/4xi/3bgyiHW9AMubC0C83ai0marNCO8ma+p4vJs+ByDWfjQlox7BymxsczKpDUas+OgfbJkEpk0AINrtxyrMIodxTK+zaqLSLCIV2lYcZdaGfczesI/VmzZyYnIeZznmcpNjKV53smy7UGZLUl3GkuoyjmSzE3AZDpJH2G994N70BdmfjscR2YXl9BIcPolor8t0sZ9Uinflv3DvXITp9hM68Y92xxGRI1BpFpFD7A7FWbCpkPmbC1mwuZBY4TbOcM7jKsc8BjtW4HKbZdtGszuR6jKWeKcxJPN6NZyymIrhn/Ugvm9eBCDZuBvFo58l1aS7zcGkrjBixQRm/QmA8IDfYfqb2ZxIRI5EpVlEKAwnWLClkHmbSkvyhr0RWhs7OdMxj8ed8xiQseqg7RN5vYl3OotYxzGkGnexKbV9nPvWkjXlGty7lwIQ6XMFwaG3a7ljqRLf/CdxRHaRzO1I5Lgr7Y4jIhVQaRZpgIqjCb7eXLT/THIRa3aHAOhkFDDWMY8zPXPp49hw0GMSzfsT63gWsY5nYea0syF1GrAsMlb8ncD0SRjJCGZGI0pGPZY2M3dI3eHct4bMxa8AEBp+p2ZXEakDVJpFGoCiSILFW4tZsLmIBZsL+XZncP+MxxY9jY383jWXc93zaWdtKXuMZThItBxCrONZxDueiRloYVf8tGBEC8n64ha8a0tXKIy3HkHJaY9j+pvbnEzqHMsq/cXLTBJrf1rpRa0ikvZUmkXqGcuy2FYcY1FBEd8UFLOooIh1e8Jl9xuYnGCs4f/5F3KaMY+8xP7lqy2wHG7irYcT7zSGWIfRWJlNbPooDs9yumv9mO6tc8j65Dqcwa1YDhehwbcQOeHXYDhqPYukn6o+Jz0bPsGz+Ussh4fgsEk1lEqkfrHje/8PaZ5mkTouZVqs3R1iUUERiwqK+aagiJ3B+EHbOElxdvZ6Lsj8mv7RmfhjO8vus1wZxNuOJNZxDPH2p2olsu8zk/jmPY5vwdMYlkkypwMlo58h2fQ4u5NJXZWM0vjvp+Is3ki43zWaMUMkTWieZpF6KJpIsWx7SdlZ5MVbiwnFUwdt43QY9Gnq5fzsNYxIzaLt7i9xxfbC/i5tugPE259KrNMY4m1PAbfPho8kvTmKN5P9ybW4ty8AINL9QoIj7gaP3+ZkUpf5Fr2Es3gjKX8zQv2vtzuOiFSBSrNIGjMti037IizbVsLSbcUs217Cql0hUubBfyDye5z0aZnNgBYeRrmW0G3fF2Ru+gxH4XeTwZveXGIdziDeaQzxNsOPblW+BsK7+j0CX9yKI16C6ckiOPIBYl3OtTuW1HGO4FZ8C54CIHTi7foFTKSOUWkWSSOF4QTLtpewZFsxy7aVsGx7CSWxQ5cHyQ94OL5VDse3yqZfUyc9grPIXP86niWfYyS/W7465WtKvOOZxDqOIdFqCBzrqnx2M1O4diwEINnshGpfEdCIBwlMn0jGyreB0uVai09/GjO7TbUeR+qRKjwn/TPvw0hGSDQfQKzr+bWVUKTuq+Hv/ZVVx3+CitRd8aTJql1Bln7vLPKWwugh23ldDro3DdCrRRa9W2TTu0UWLdxhvBs+wbt2Mp7Z0zHM78Ywp7Jal04N12lM6ZK89exiNUd0X43s17VjEVmfXIuraAOW4SA84AbCA26o+79oSI2rzHPSvXUOGavfw8IgeNI9DWcRIJFqUlPf+6tCPw1EakE8abJmd4iVO0pYuTPIyh1B1uwOkUgdeh1uu0aZ9P5eQe6c58fldOAI7cCz/kO8n0/GXTALw/puHHMyt1Pp+OSOZ5HM76MfyFVhmWQufB7/nIcwzCSpQEtKTn+aRMvBdieT+sJMEZg2AYBoz5+WvkZFpM5RaZaGIxnFu/Z9vOs+xogWYmXkEut4BrFO48CVUW2HiSZS+wtyaTlesaOEtXvCh4xDBsjNdNO7RRa9mmfRu0UWPZtnkZ3x3bQ6juLNeJf8A++6ybi2zcfgu30km/Qk1mkMsY5nkWrcVUX5KDhC28n69Hd4tswAINppHMGRD2Bl5NobTOqVjOV/w7VnOaY3h9CQm+2OIyJHSaVZGgTP+ilkTR2PI1aEhQMDEwsH3nUfYk6fRMmpTxzVqm6RRIpVO4N8uzPIiv0lef2eEOWcQCYnw0WPZll0axagR7MA3ZoGaJWTgfGDsuvct7Z02MW6D3HvWnzQfYlmJxDrOIZYxzMxcztUOa98x7P+E7I++z2O6D4sVybBEXcT7XGRfvmQamVE9+Gf/SAAoUE3puXc5yJSOSrNUu951k8he/KVZf83MA9+Gysme/IvKB7zCvEOo8vdh2lZbCuOsnpniNW7Q6zeFWLNriBbCqOUN9F5o0w33feX4+7NsujeLEDzLO8hBRkAy8K5ZwXetR/gXfcRrr3ffneX4SDRYlDp0IsOZ2JmtTz6T4SUSkYIzLyXzCWvA5DI603J6GdJNepkczCpj/xzH8ERKyTZuBvR3pfbHUdEjoFKs9RvyShZU8cDHDS04fsMLCwMsqaOZ88VCwiZLtbuDrN6V5DVu0oL8trdoUPmQj6gid9TWo6bBui+vyQ3DXjKL8gHWCauHYvwrpuMd+2HOIs3fneXw0Wi9bDSM8odzsDy5R39xy8Hce5ZSfaUa8p+MQkf/+vSP5dr+j2pAc7dy8lY+heA0jm+dVGpSJ2mV7DUa9617+OIFVW4nYGFESviuVef5tWS8i8AczsNOjT20aVpgC55frrk++mc76exz1O5MGYK97a5eNZOxrv+I5zBbWV3WU5v6ap8nc4i3v50rcp3JEczG4hlkbH0dQJf3YORimFm5lN82hMk2p5c/fmk4SnvOWlZBKZPxLBMYp3Gkmg9rPZzidQnaTATlEpzPedcsRz/n+4m9MeJpHr0tDtOrfOu+xgL47Bnmb8vZRkMjM7kVQaTH/CUluK8AF3ySwtyu0aZuJxVfNGm4ri3fIV33Yd413+MI7Kn7C7T7Sfe/jRiHc8qXZVPCx1UzOEk3v7UKj3EiOwl67Mb8W74BIBYu1GUjHpMZ/ClehzmOeld8z6erbOxnF6CQyfYEEykHjmK7/01QaW5nvN88jHejyaTGDiESD0vzZFYnJ0FqynZsgxr1woyi9ZwXGRWpQozgNOwGNTU4pNzTiTX5674AYeTjODZ9CXedR/i2fDpQWe6TW8O8Q5nlBblNiOqddYOOZR78wyyPr0BZ3gHlsNDaOjtRPr+Qhf7Sc1KhPHPvAeAcL9rMLNb2xxIRKqDSnM955n+RdnbyHW/szFJ9UiZFtuLw+wqWEt023Ice74lq2QNTWMbaG9toa2ROOp9WzgI5ORjHkVhNuJBPBun4ln7Id6Nn2Ekw2X3mZn5xDqeSazTWSRangjOYyjkUjmpOP65j5D59XMYWCQbdaF49LOk8ur3L46SHnxfP4szuJVUVmvC/X5jdxwRqSYqzfVZIoF7ziwA3LNnQSIB7rpR2IoiCTbuDbFn61oSO1fg3reGnNBaWsY30NXYyvFG7NAHGRDFTYGzLft8HYg16kpzZwmd1v+lUsc0MIl1PLPSGY3oPjzrPyk9o7x5Gkbqu0ypQKvS8ckdzyLRfIBtS37WO/svoARINju+3DFuzsJ1ZH1yHe6d3wAQ6XUZwWETwZ1Zi0GlwfjBc9JRsgXfwucBCA6bAC4970SOWSW+99cGleZ6zPXNQoxo6bLMRjSCa/Eikv0H2pyqlGlZ7ArG2VIYoWBfhJLdGzB2f4u/ZA35kfW0tzYzxCjA/8NyvP91EsfNdncbCv2diOd2wd2sBzlteuPP70i2w0n2ge2TUcw//xcjVnzEYRoWBpY3m1insUfMbYR24l3/Md51H+IumIlhJsvuS+Z0IN5pTOny1fl9NQSgJljWd+PCLQuMg+/zfvsvsr68HSMZxvTmUjLqYeIdz7IlqjQQP3hOBr66GyMVI95qGPGOY+zNJlJfHOl7fy1Saa6n4rFitr19FzkGOCwwDdj89p006f13PN7sindQHRmSJttLYmwpjJT+2xchvHcz3sLVNA6vo6O1ma6OLQw3CggY0e8eaFD2gkjgYpenLSVZnTCbdMPbvAdZrXpDbjsyHS4qPIfjyqDk1CfInvwLoobBJ74MdrhcBA0Hm9wuTolEOD0UwWtByalPlDvG2FG8Be/6j/CunYxr27wfrMrXvXRquE5jSDXupqJcg+KpGF9u+5Sitf8gnAyzbs9UhrY4mZHNR+FNxgh8+UcyVr9Xum2rEyk57UnMgOa1lprzw+ekZ+1fuW7dVCzDSXDEXfp+IFINjvS931PL04UalmVV6iqpXbtKajqLVJO5i5/ivk1/45FHNzJkeaisNM/q6eemG9txe9tLGdT32mM+TjxpsqMkxtbiKNuKomwrjrK1OMa2wgjJ4q00iayni7GFLsYWujq20NkoINuIlLuvFE72ZbYjktMZ8rqR2bwnjqY9SOW0r5a5TQ98TkKGwcnhCKYBX2VkkHQ4yDLNQz4nzsJ1pVPDrfuw7M/8BySaHle62EjHs0jldjzmbFKxr3ZM58Fv7iWcKOakSAwLkxmZflKGxYkJJ4/vDeIP78IynIQH/YFwv99qSIzUqB8+Jx2WyfV7C+mUTPJWbmMyTvs/hjYbbndMkTrtSN/7A64sbj1uQrW9zvLzsyrcRmea6yjHtq04du085P1LV/+DN7Z/RBug36owjv2/Ejks6L8qTJuNUd7Y+DLZq7fQu8tFhzzezG+K2aL07Fw4nmJHSYwdJVG2F8e+K8X7C/KuYIw8Cunq2ELX/eV4jKOALsYWcowwlDN9sYmTUKAdycZdceV3x8zrTqpxV1I5HcDp5sB5XgsofymRqvtqx3QmbnkLHE4My8Tcf/LH3H8WKOhw8sfNf+MJby6Di3aUnlH+/qp8GCRaDCLe6SxiHc/CzGpVTcmkMr7aMZ2JC24FSv8AYR1YydFKcXVhCb8uLMIFhPxNiZ/5Msnm/ewLKw1Cec/JEeEInZJJ9jocPJXto2TBLdzd/wGGNRthb1iROupw3/tL3xqEkkEm1PLrTGea66jsi3+Ed+onR9zG3D8043D/L8/C3kO56Wf3saMkRknswHhdizyK6bK/HHc1tpTdzjVC5e7HMpzEs9thNelWWoobdyPZuGvpmVlnJRcDqQbxVIwfTz2HUDKIhYXDshgRKT3bPT0jg56JBKeFIpwWCtMu+d34ZMvhItFqGLGOZxHreAaWL7/WMst3Dvf1a5QyOa84SP94HID/+f080awNb4z+oNb/XCcNS3nPyTHBEBN278UH3NmkMf/ODmBg4HcF+Nep/9VzUqSKjvizOzOz7KRXdb7OqvVMcyhUfjkSe5iX/IzGq77FvXkTFuWPif9hQT5cYT7w+I05zXii80i2bd3Mee7ZdHMW0MO1lS6OAnIJlv9Yw0Eyqy3x3M4kcjuTyO1SejvnMOU4mgCOflq4qvpsxycEkwf/wtc0kWJEJMrtu/bSwjTL3h8zYEd+HzK7XkqkzUhMb27pHRag578tyvv69Y3GuLS4BJ8FQcPgvrzGvB/wgxVmysaPOKXZaTallYagvOfk2cEQPmC528W7WaWLFFlYBJMlek6KHIXyXmflqc7XWWVKc6XPNBu6oCHtuIFrgLuBTI5urE0SCAMTgf+jtM6+fl4Glx93cOE1LYt1+yyW7UyxbJfJsl0plu00+XaPSTR56H7TRZtr25DdLxvDUfr8dVgWUzcVkLe/LIcNg2m+TD71ZTItI4Pti4JsfmaznZHle3749ctKpZi2qQAXsNjj5pameWzZP42iZVoUf12sr5/UqB8+J/2pFDM3FeAALm/elIWZ311MrOekyNEp72d3eWeaofpeZ5WpwxrTXIclgCeAv1BanK82SodguMwjPgyApKP0zPMLFkwC9nzvvtcWJfA4DTYUmqUFeWeKlbtNImlcjg/H6XeWveigdBzzc41y6BGP82VmJrMyM4g5HAdtL+njh1+/EoeDP+dkU+xw8JecLJLf+8ZpOAx9/aTG/fA5GXI4eCE3mz1O50GFGfScFDla5f3s/tLnK3fb2nydVfpM84YN22s6S51nWRbRpMm+cJzCSJLCSKLs375wkqJofP/tA+9PEktWouF+T5bXRSOfizy/h7yAlzy/hyYBD3k+D1tm/ZRzXl9Dn3URjjTttwUs7pjJBz/vwk8v/fiYPuZ0d/+yu5i9+yusSiylbWAwJG8Yt/WaVAvJpDL09ZN0o+ekSM2z43XWvn3zCrep9Jlmv99/TGHqGsuyCCdSFEWSFEUTFEUSFO0vwqX/T5a9LYwk2Le/HFe+BDvA4cHhAb/HSWOfmyZ+D419Hpr4PTTxu7+7vf++Rj4PXtfh6/DnyTNY/8lmem2I4DhCjJQD1rX00qX/GfX+63pyq1OYtXtGpba1sBjZalS9/5zUJfr6SbrRc1Kk5qXr66zez56RTJmUxJIUR5PfvY0mKY6Vvi2KJijZf19RNFlajve/TZqV+tQcwuty0CjTTSOfm9xMN419bnIzPTTyucvef+C+Jn4Pme7q+bNCYu82mvTpjjdRce6Y22DPkpW4G7eolmOnqx9egXs4utI9PenrJ+lGz0mRmmfH66zOz9N8YLhDMJYkGEuVvo1/73YsSTCeIhg9+P3fL8fhxLHN9utxGuRmusnJdJOT4dr/1k1Opuugt419bnJ9bhplesh0O2y5cDIw5bNDCrNJ6crTB94e4E1YBD75nNiFP63FhLXP4/Ry63ETmLDgFgyMcl98xv65R249boJ+uKUZff0k3eg5KVLz0vV1Vu1nmg8U3UgiRThe+i9U9jZJqOz/379det+BbUPxFKH9hTh1lGd7f8jvcZKd4SLL6yp9m+Em2+siK8N10PtzMt3kHijDmW4yXPYU4KORc+5ZuOfMwtg/M0TSAXG3wd9PbcLFU/fgSVhlFwlaDgeJwUMpem+yjYlrz4FVhYLJEgwcWJhlb6t7VSGpfvr6SbrRc1Kk5tXm66xap5yb9M5iIvEU4USqrBAf9DZhEtl/u3pq7nccBvg9LrK8TvxeFwGvi4DHWfrWW/r+gNdVep/Hub8If1eKA14XLkfdKL5Hy7FpI00G9Cn7vwXs7d6Cp69szoYmBu33WFz38nYaf7vtoDmd9yxYitmmba3ntUM8FePL7Z8zY/s0ihNFZLtzGN78JE5uforOBtUB+vpJutFzUqTm1dbrrFpLc/tbP6hygEy3g0x3aaH1uZ34vc79b134PU78Hic+jxOf57v/+7/3/9JSXPqYunK21y6+Rx7A/9D9WA4HWBbhP9xKePxN4PreCJxkEt/jD+N75AEwDAzTJHTL7YRvvMW+4CIiIiI2q9bS/Me3F5UWXLeTDLcTn8eBz+0kc//7Mt2lBThz/+0MtwOHim7tME0a9+uFc2sBqZatKH7xzyQHDT7s5q65c8j+1c9wbt1KqlVr9n69DPS1EhERkQaqWktzXZ09o0GIRMj90dmk2rUn+OCjWNk5FT7EKC4icPPvcW7aSOG//weZmbUQVERERCT9qDQ3JKYJjiMtaVLNjxMRERGpJypTmtWW6oujLb4qzCIiIiIVUmMSEREREamASrOIiIiISAVUmkVEREREKqDSLCIiIiJSAZVmEREREZEKqDSLiIiIiFRApVlEREREpAIqzSIiIiIiFVBpFhERERGpgEqziIiIiEgFVJpFRERERCqg0iwiIiIiUgGVZhERERGRCqg0i4iIiIhUQKVZRERERKQCKs0iIiIiIhVQaRYRERERqYBhWZZldwgRERERkXSmM80iIiIiIhVQaRYRERERqYBKs4iIiIhIBVSaRUREREQqoNIsIiIiIlIBlWYRERERkQqoNDcwB2YYfPrppw/6f331/vvv06dPH3bs2FH2vttvv52zzz6bkpISG5NJZejrJ+lGz0mRmpeurzOV5gbmb3/7G2+99RbhcJhHHnmEefPm2R2pRo0dO5YOHTrwwgsvAPDMM88wY8YMXn75ZbKysmxOJxXR10/SjZ6TIjUvbV9nljQ4L730ktWrVy9r/vz5dkepFZ999pnVu3dv67nnnrMGDhxorVq1yrIsyyopKbEuuOAC65xzzrHGjRtnvfXWWzYnlfIc7uu3detW69JLL7XOOussa9y4cdbkyZNtTioNxeGek5ZlWb/97W+tAQMGWNddd52NCUXqviO9zj777DNr9OjR1umnn269/fbbtZZJKwI2MH//+99xOBxs2LABl8vFySefzIABA+yOVePOP/98Vq9ezauvvsqgQYMASKVSxONxMjMziUQijBs3jn/96180atTI5rTyQ+V9/Xbu3MmePXvo0aMHe/bs4fzzz+ejjz7C5/PZnFYagvKekwCzZ88mHA7zn//8h6eeesrGhCJ1X3mvs2QyydixY3njjTfw+/1ccMEFvP322+Tm5tZ4Hg3PaGAuuugiLrzwQnw+HzfeeCP9+/e3O1KNmz59OuvWrcM0TZo0aVL2fqfTSWZmJgCxWAzTNOv9GO+66HBfv6ZNm9KjRw8AmjRpQk5ODkVFRXbFlAbkcM9JgCFDhuD3+21KJlJ/HO51tnjxYjp37kyzZs0IBAKcdNJJzJgxo1YyqTTXE5UdNG8YBgDXXXfdQf+vr5YtW8YNN9zAXXfdxfDhw3nyyScPur+4uJhzzjmHk08+mauuuorGjRvblFTKU9HX74AlS5ZgWRYtWrSo5YTS0FT2OSkiR+9Ir7OdO3fSrFmzsv83b978oO5Tk1Sa64m0HTRvoy1btvDrX/+aX/3qV5x33nlcf/31TJkyhaVLl5Ztk52dzX//+1+mTp3K//73P3bv3m1jYvm+ynz9APbt28ctt9zC3XffbVNSaSgq+5wUkaNX0eusvL8I19YJQJXmesIwDMaPH88///lPnn/+ed544w1efvnlg34ba0gKCwv55S9/yahRo7j66qsB6N27N6eccgqPP/74Idvn5eXRrVu3ej+bSF1R2a9fPB7n2muv5Ve/+hX9+vWzK640AFX9niIiVVeZ11mzZs0OOrO8fft28vPzayWfLgSsZw53cYocavfu3WRkZBAIBAgGg1x44YU8+uijdO/e3e5oUgmWZXHjjTfSoUOHsuFGIulgzpw5vPnmm7oQUKQGJJNJxowZw1/+8peyCwHfeuutWrmI31XjR5Bac6SLU+RQ27dv5/bbb8eyLCzL4pJLLlFhrkMWLFjA5MmT6datG59++ikADz30EN26dbM5mTRkV155JcuWLSMSiXDSSSfxzDPP0LdvX7tjidQbLpeLW265hcsvvxzTNLnqqqtqbdYrnWmuJ5YtW8Zll13GxIkTmTx5MhkZGTrLISIiIlJNdKa5HvjhoPnOnTvz4x//mKVLl9K7d2+744mIiIjUeTrTXMcVFhZy8cUXM3DgwINmD/jNb35DPB7nlVdesTGdiIiISP2g0iwiIiIiUgFNOSciIiIiUgGVZhERERGRCqg0i4iIiIhUQKVZRERERKQCKs0iIiIiIhVQaRYRERERqYBKs4iIiIhIBVSaRUREREQqoNIsIiIiIlIBlWYRERERkQqoNIuIiIiIVEClWURERESkAv8fRmL12osCtJEAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "def f(x):\n", " return x**5 - 3, 5*x**4\n", "\n", "xmin, xmax = 1.0, 2.55\n", "xinit, xstar = xmax-0.05, 3**(1/5)\n", "xx = np.linspace(xmin, xmax)\n", "yy, dyy = f(xx)\n", "\n", "n = 5\n", "x, y = np.zeros(n), np.zeros(n)\n", "x[0] = xinit\n", "for k in range(n-1):\n", " y[k], dlag = f(x[k])\n", " x[k+1] = x[k] - y[k]/dlag\n", "\n", "xticks = np.array(x[:4].tolist() + [xstar])\n", "\n", "fig2, ax = plt.subplots(figsize=[9,6])\n", "ax.set(\n", " title=\"Newton's Method\",\n", " xlim=[xmin,xmax], \n", " xticks=xticks,\n", " xticklabels= ['$x_0$', '$x_1$', '$x_2$','$x_3$', '$x^*$'],\n", " yticks=[]\n", " )\n", " \n", "ax.plot(xx,yy)\n", "ax.axhline(0, color='k')\n", "ax.annotate('f',[xinit,f(xinit+0.03)[0]], fontsize=18, color='C0')\n", "ax.plot(xticks, f(xticks)[0], color='C1', **bullet_options)\n", "ax.plot(xticks, 0*xticks, color='C2', **bullet_options)\n", "\n", "ax.plot(xstar, 0, color='r', marker='*', ms=16)\n", "\n", "for xi,xinext,yi in zip(x,x[1:],y):\n", " ax.vlines(xi,0, yi,colors='C1', linestyle='--', alpha=0.3)\n", " ax.plot([xi,xinext],[yi, 0], color='C1')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Secant Method" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAs0AAAINCAYAAADFp0I1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB2UklEQVR4nO3dd5hU5fnG8e+ZurtTdikLC0jvRVRAQEFAEFQgJhoTk58pdo1iwd5b7LGLSYwtMTGJmmgSFZGmCKIoKNJ77213Z2Z3+jm/PxZWkbKU3T0zu/fnurhmdubMOc/4jsO9h/c8r2FZloWIiIiIiByQw+4CREREREQynUKziIiIiEgVFJpFRERERKqg0CwiIiIiUgWFZhERERGRKig0i4iIiIhUQaFZRERERKQKCs0iIiIiIlVQaBYRsUkmrC2VCTWIiGQDhWYRyWrLli1j7NixDBgwgB49ejBw4ECuu+46Fi1aZHdpBzVnzhwuv/zyg27z9ttv07lzZzp37szq1av3u80nn3xSuc3hCIVC3HLLLcyePbvysV/+8pf88pe/PKz9HI6a3r+ISE1SaBaRrLV8+XLOO+88du3axR133MErr7zCzTffzKZNmzjvvPOYO3eu3SUe0FtvvcWKFSsOaVuHw8EHH3yw3+fGjx9/RMdfvHgx//nPfzBN84heLyJS3yg0i0jWevXVVykoKOCll15i5MiR9O3bl7POOos///nPNGzYkN///vd2l1gtevXqtd/QnEgkmDx5Ml27drWhKhGR+kWhWUSy1o4dO4B95+Xm5eVx2223ceaZZ+71+OTJkznnnHM49thjGTBgAA888ADl5eV7bbNgwQIuueQSevfuTf/+/Rk7diybN2+ufH7JkiWMGTOG/v370717d0455RQeeOABYrFY5TadO3fm9ddf54477qBv376ccMIJXHPNNZX13nrrrbzzzjts3LiRzp078/bbbx/0fY4cOZJly5axcuXKvR7/5JNPMAyDQYMG7fOa2bNn84tf/ILjjjuOvn37csstt7Br1y4AZs2axa9+9SsAfvWrX+01ZcKyLF588UWGDBlCz549Oe+885g/f/5e+54/fz4XX3wx/fr1o1evXlxxxRUsX758r202bdrEmDFj6N27NwMGDODVV1896HsUEcl0Cs0ikrWGDBnCpk2b+NnPfsbrr7/OypUrKwP0GWecwdlnn1257bvvvstVV11Fu3bteP755xkzZgz/+9//uPLKKytfs2TJEn7+858TjUZ55JFHuP/++1m0aBEXXXQRyWSSbdu2cf7551c+/+KLL3LmmWfy17/+lT//+c971fbUU09hmiZPPvkkN998Mx9//DEPPfQQAFdeeSWDBw+msLCQN954gyFDhhz0fQ4YMID8/Px9zjaPHz+e4cOH43a793r8yy+/5IILLiAnJ4enn36a22+/nS+++IJf/epXxGIxunfvzt133w3A3XffzT333FP52jlz5jBp0iTuuusuHn30UbZu3coVV1xBKpUC4PPPP+fnP/85pmny4IMP8sADD7B582Z+9rOfVYb68vJyfvGLX7BkyRLuv/9+7r77bt566y2+/vrrQxlWEZGM5LK7ABGRI/V///d/bN++nZdffpn7778fgAYNGjBw4EB++ctfctxxxwEVZ08ff/xxTjnlFB5//PHK17dp04YLLriAadOmMWTIEH7/+9+Tn5/PK6+8gtfrBaCoqIjrrruOpUuXUlJSQteuXXnmmWfw+/0AnHzyyXz22Wd8+eWXXHHFFZX77tSpEw8//HDlz/PmzWPChAkAtGrVioYNG+LxeDj++OOrfJ8ul4vTTjuNDz/8kDFjxgAQjUb56KOPeP7555kzZ85e2z/xxBO0bduWF154AafTCcBxxx3HqFGj+Pe//835559Phw4dAOjQoUPlfQCPx8Of/vQnCgoKAIhEItx5552sWLGCLl268MQTT9CyZUteeumlyn0PHDiQ4cOH89xzz/H000/zzjvvsGnTJv773/9WXqDYs2dPhg8fXuV7FRHJVDrTLCJZ7dprr2X69Ok88cQTnHvuufj9ft59913OO+88/vKXvwCwatUqtmzZwtChQ0mlUpV/TjzxRPx+P59++ilQcZZ10KBBlYEZKsLe1KlTKztz/O1vf8Pr9bJ69Wo++ugj/vjHP7Jr1y4SicRedX0/DBcVFRGNRo/4fX5/isZHH31EXl4e/fr122u7aDTKN998w+DBg7Esq/K9tmzZkvbt21e+1wPp0KFDZWAGOOaYYwAIh8OUl5czf/58Ro4cWRmYAYLBIKeeeiqzZs0CKqaGtGzZcq+OHs2aNTukXxBERDKVzjSLSNbLz89n9OjRjB49GoBFixZx88038/jjj3PWWWdRUlICwH333cd99923z+u3bdsGQElJCY0aNTrgcfZMt3j99dcpLy+nWbNm9OzZc6+QvUdubu5ePzscjqPqidy/f38aNGjABx98wJgxYxg/fjxnnHHGXuEVKlrJmabJiy++yIsvvrjPfvZX63fl5eXtUzdUvPdwOIxlWTRu3Hif1zVu3JhwOAxAaWkpDRs23GebwsLCynndIiLZRqFZRLLS1q1b+fGPf8y1117LT37yk72e69atG9dddx1XXXUV69evJxgMAnDzzTfTt2/fffaVn58PQCAQqLxY7rumTZtGly5deOedd/jzn//Mvffey+mnn04gEADg3HPPre63tw+Xy8WIESOYMGECF1xwAZ988sk+86gBfD4fhmFwwQUXMGrUqH2e/36YPxyBQADDMPYbfLdv3155hrpBgwasXbt2n232/PIiIpKNND1DRLJS48aNcblc/P3vfycej+/z/KpVq/B6vbRu3Zp27drRqFEjNmzYwLHHHlv5p6ioiCeeeKJyIZQ+ffowffr0vaZaLF26lMsuu4z58+czZ84cOnTowLnnnlsZmLdu3cqyZcsOu9/xnjO4h2PkyJEsX76cl156icaNG3PCCSfss43f76dbt26sWrVqr/fasWNHxo0bVzmF4vtnqA9FXl4ePXr0YPz48aTT6crHw+EwH3/8Mb179wYqzopv2LBhr64bu3btyui+2SIiVVFoFpGs5HQ6uffee1m2bBk//vGP+cc//sEXX3zBtGnTeOihh3jmmWcYM2YM+fn5OJ1Oxo4dyz//+U8eeOABPv30Uz744AMuuugiFi1aRPfu3YGKrhbFxcVceumlTJ06lQkTJnDdddfRvXt3Bg0aRM+ePVm6dCl/+tOf+OKLL3jrrbc4//zzSSQShz1fORgMsmPHDqZNm1Y5PaQqffv2pbCwsLIvtWEY+93u+uuvZ8aMGdxwww1MmzaNqVOncskllzBz5szK97on9H/88ccsWbLkkOu+4YYbWLt2LZdccglTpkxhwoQJ/PrXvyaRSFRepPjDH/6QTp06MWbMGP7zn/8wefJkLr30Ui2kIiJZTdMzRCRrDRkyhDfffJOXX3658oI8j8dDt27deOqppxgxYkTltj/5yU/w+Xy89NJLvPHGG+Tl5dGrVy8ef/xxWrZsCVRM6/jrX//KE088wdixY/H5fAwePJgbb7wRj8fD5ZdfTnFxMa+99hrPP/88zZo144c//CGGYfDCCy9QWlpaOdWjKueccw7Tpk3jqquu4pprruGyyy6r8jUOh4PTTz+dv/3tb/uderHHwIEDefnllxk3bhzXXHMNbreb7t278+qrr1ZejNexY0dGjx7N66+/zvTp03nvvfcOqe6TTjqJV199lWeffZbrr78ej8dDnz59ePTRR+nYsSNQ0YHjL3/5Cw899BAPPvgghmHw05/+lJYtW7Jz585DOo6ISKYxrKO5MkVEREREpB7Q9AwRERERkSooNIuIiIiIVEGhWURERESkCgrNIiIiIiJVUGgWEREREamCQrOIiIiISBUUmkVEREREqnDIi5ts3x6uyTpERERERGxRWBiochudaRYRERERqYJCs4iIiIhIFRSaRURERESqoNAsIiIiIlIFhWYRERERkSooNIuIiIiIVEGhWURERESkCgrNIiIiIiJVUGgWEREREamCQrOIiIiISBUUmkVEREREqqDQLCIiIiJSBYVmEREREZEqKDSLiIiIiFRBoVlEREREpAoKzSIiIiIiVVBoFhERERGpgkKziIiIiEgVFJpFREREpF6KxFP8a+6mQ9rWVcO1iIiIiIhUG8uy+PucjbwzbzObQzFy3U6GdGjMnad3Oux9/fbDZUxdvoPfDO9c5bYKzSIiIiKSNV6fs5Fnpq0CoHnQSzDHTfP8nMPezzcbS5m6fAcO49C2V2gWERERkawxeel2AEZ0LuTB0V2PaB+WZfH07uB9Vo+iQ3qN5jSLiIiISNYoiSYB6N0y/4j3MXnZDhZsDpPrdnD5ya0P6TUKzSIiIiKSNUzLAsDtPLIYm0iZjJu+GoBfntiSxn7vIb3OsKzdR67C9u3hIypMRERERORoXf7GN3y1oXSfx5sFvfzv0n6HvJ/XZ2/g6WmraOzz8PbFJ5LrdlJYGKjydZrTLCIiIiIZr0NjH2nTYvHWMIm0RasGuTTIddPY7znkfZRGk7z8+ToAfjOgDblu5yG/VqFZRERERDLeTcM6AHDWi7PYHIpzQd+W/OAQL+Lb45VZ6wjHU3Ro7GNU96aH9VrNaRYRERGROm9DSZQ3v65YyOTawW1xHmqvud0UmkVERESkzhs3fTUp06J/mwb0b9PwsF+v0CwiIiIiddo3G0uZsqxiIZNrB7c7on0oNIuIiIhInWVZVuUKgj/oUUSHxr4j2o9Cs4iIiIjUWVOW7WD+5jA5LgdXHOJCJvuj0CwiIiIiddJ3FzL51WEsZLI/Cs0iIiIiUif965tNbCyN0djn4RcnHnNU+1JoFhEREZE657sLmVwxoPVhLWSyPwrNIiIiIlLnvDJrHaFYivaN8xjd/fAWQdkfhWYRERERqVP2Xsik3WEvZLI/WkZbRERERLLG/y7tV+U2z09fU7GQSesGnHQEC5nsj840i4iIiEidMW9TiMnLtmMA1wxuW237VWgWERERkTrBsiye/njPQiZN6Vjor7Z9KzSLiIiISJ0wdfkO5m8OkeNycPnJbap13wrNIiIiIpL1kmmT5z6pWMjkF32OoUngyBcy2R+FZhERERHJem/NrVjIpGGem1+e2LLa96/QLCIiIiJZLRT77kImbcjzHN1CJvuj0CwiIiIiWe3lz79dyOSsHke/kMn+KDSLiIiISNZas6ucN3YvZHL1oOpZyGR/FJpFREREJGs99fFK0qbFgLYNGdC2ehYy2R+FZhERERHJSjNW7WTm6mJcDoOxQ9rV6LEUmkVEREQk6yRSJk9+tBKAn/dqQeuGeTV6PIVmEREREck6//xqI+tLKlrMXdS/VY0fT6FZRERERLLKjki8ssXc1YPa4ve6avyYCs0iIiIiklXGTV9NeTJN96IAI7s1rZVjKjSLiIiISNaYvynE+4u2AXDT0PY4jJppMfd9Cs0iIiIikhVMy+Lx3Rf/je7elO7NgrV27JqfACIiIiIiciRSMbwr38O76kOMWAkbEzl03taFzZ6TueqUtrVaimFZlnUoG27fHq7pWkREREREAPCsnkhgylgc8VIsHBiYpDFwYhFzBoif/iyJtsOr5ViFhYEqt9H0DBERERHJKJ7VEwmOvxgjHgLAwATAScW5Xm86QnD8RXhWT6y1mhSaRURERCRzpGIEpowFwGD/EyL2PB6YMhZSsVopS6FZRERERDKGd+V7OOKlBwzMexhYOOKleFe+Xyt1KTSLiIiISMbwrvoQ6xAjqoUD76oJNVxRBYVmEREREckYRqykcg5zldtiYsRKarag3RSaRURERCRjWDkFh3Wm2copqNmCdlNoFhEREZGMEW93+mGdaY63O6OGK6qg0CwiIiIiGSPefjSm21/FZYBgYWB684m3H1UrdSk0i4iIiEgGsUi4/BjAgZbgszAACA97Glw5tVKVQrOIiIiIZAzfpw+QE91CqZVH1JEHfBuSK2+9QUIjX6m2FQEPhavWjiQiIiIichCeNZPJW/AXAG62rmXs+b+i+eYJ5C74KyTLMPPbEm9/ZsWUjFo6w7yHQrOIiIiI2M4o24Z/yvUAvJw6k84nn0Wj/ADx4I9JNeoGQLpxVzDsmSih0CwiIiIi9rJMglPH4oztYrHZin/4L+C1Xi0qnjMcpAu721sfmtMsIiIiIjbLnfcKnnXTiFlurk6OYczQrridmRVTM6saEREREalXnNsX4pv5EAAPpH5B83bHMaBtw703MtMVf2yk6RkiIiIiYo9klOCkMRhmgknpXvzLGMGbQ9vvvY2ZxrN2KgCJ1kPB4bShUJ1pFhERERGb+Gf+FlfxcrZTwC3Jy7j0pDY0C9ZuV4xDpdAsIiIiIrXOs3oiuQteA2Bs4jcUNCri/3q3sLmqA9P0DBERERGpVY6yrQSm3gjAi6mRzDCP5U+ndcSVYRf/fVfmViYiIiIidY9lEpgyFkdsFyscbfld6jxGd2/KCcfk213ZQSk0i4iIiEityZ37Ip71n5B0eLk8eiW5OblcM6it3WVVSaFZRERERGqFa/sCfJ8/AsBDqV+y0mrBmFPa0iDPY3NlVdOcZhERERGpeckogUljMMwkX+WexKvFp3JssyBnHVt08NcZBqavSeV9uyg0i4iIiEiN8396H67iFUS9hVxc/GuchsFtwzvgqCoIGw5STY6rnSIPQtMzRERERKRGeVZNIHfh3wC41bySYoKc16sFHQv9Nld26BSaRURERKTGOMq2EPjoJgBmFv6c/4Y708Tv4bKTW9tc2eHR9AwRERERqRmWSWDydThixZQ16MbFG0cCcMPQDvg8hxhDtYy2iIiIiNRluV+/gGfDDCxXDrdxDVHTyYC2DTm1QyO7SztsCs0iIiIiUu1c2+fjm/UYAJ+3v4H/bQ7idTm4cWh7DBu7YBwphWYRERERqV7JcgITr8Iwk5S1Pp2rllV0v7i4fyuOKci1ubgjo9AsIiIiItXKP+NeXCWrSPuKeNT9G3ZFU7RpmMsv+hxjd2lHTKFZRERERKqNZ9UH5C76OxYG83s9zF8XlANwy7COuJ3ZGz2zt3IRERERySiOyCYCUyvay5UdfwW3ftMICzizaxP6tCqwtbajpZZzIiIiInL0zHRFe7l4CcnCnrzmPZ8l29YR8Lq4dnC7I9+vYWDmNqq8bxeFZhERERE5arlz/4hn40wsVy7rBjzJ829vBOCqU9rQyOc58h0bDlJFvaqpyiOn6RkiIiIiclRc277BN+t3AEROuZ9Hv7YoS6TpXhTg7J7NbK6ueig0i4iIiMiRS5QRmDgGw0wRbz+Sj3JHMHHpdhwG3HZaRxxZ2JN5fxSaRUREROSI+Wfcg6t0NWl/M7ae/BAPTVoBwHkntKBzU//RH8BM41kzBc+aKWCmj35/R0hzmkVERETkiHhWvEfu4n9iYRA+7Rmen13ClnCc5kEvvxnYpvoOZJnVt68jpDPNIiIiInLYHOFNBD6+BYBoryuZTXfe/HoTALcP70Su22lnedVOoVlEREREDo+ZJjD5GhzxUpJNjqO411genLgcCxjdvSn92jSwu8Jqp9AsIiIiIocl9+s/4Nn0OZYrj/Dw53hl9hZW7yqnYZ6b646mJ3MGU2gWERERkUPm2joX3xePAxAe9FuWJJvw5y/WA3DzsA7k57rtLK/GKDSLiIiIyKFJlBGYVNFeLtZ+NGWdfsJvP1xG2rQY0qERQzs2trvCGqPuGSIiIiJySPzT78ZVuoa0vzmRIY/wz683sXhrBL/Xyc3DOmDUUE9mM8f+OdIKzSIiIiJSJe/yd8ld8kZle7n1MS9//HQBANcNbkeh31szB3Y4STXrUzP7Ppwy7C5ARERERDKbI7wR/+72cuW9x5Bo3p8HJy0nnjLp06qAs3oU2VxhzVNoFhEREZEDM9MEJl2DIxEi2eR4yk+8nv8t2MLsdSV4XQ7uGN6xxqZlZBJNzxARERGRA8r76nk8m2dhun2ERoxje9Tk6WmrALhiQBuOKcit2QLMNO4N0wFIHnMKOOxZNEVnmkVERERkv1xbviLviycAiAx6gHSwNY9NWUEknqZbUYCf9WpRK3UY6SRGOlkrxzoQhWYRERER2YeRiBCcdDWGlSbW4Szinc9l6vIdfLxiJ06HwZ0jOuJy1P1pGXsoNIuIiIjIPvzT78IZWkva34LIkIcpjaV4bMoKAC7o25KOhX6bK6xdCs0iIiIishfv8v+Ss+QtLMNBePizWN58np62il3lSdo2zOOifq3sLrHWKTSLiIiISCVHaAP+j28DoLz31SSb9+PzNbt4b+FWDOCOER3xuOpfhKx/71hERERE9s9MEZy8u71c016UnziW8kSahyctB+CnJzTnuBb5NhdpD7WcExEREREA8uaMw735C0y3n9Dw58Dh4g+frmRTKE5RwMuVA9vaUpflDdpy3O9SaBYRERERXFvmkPflUwBEBj+Amd+a+ZtCvPHVRgBuH9GRPI8NPZIdTpLN+9X+cb9fht0FiIiIiIi9jET42/ZyHX9IvNOPSaRMfjtxGRYwqlsTTmrT0O4ybaXQLCIiIlLP+afdgTO0jnTgGCKDHwbD4OVZ61i9s5wGuW6uG9Le7hJtp+kZIiIiIvWYd9k75Cx7G8twEBr+HJY3yMLNIf4yax0ANw/rQEGu274CzTTujTMBSLY42bZltBWaRUREROopR2gd/mm3A1De51pSzU4klkxz74SlpC0Y0bmQ0zoX2lwlGKmY3SVoeoaIiIhIvWSmCE66BkciTLKoD+V9rgXgD5+uYc2uKI18Hm4e1sHmIjOHQrOIiIhIPZQ3+1ncW2ZjegKEhj8LDhdz1pfwjzkV3TLuHNGRfDunZWQYhWYRERGResa1eTZ5s58GIDLoQcxgK8oSKe7/sKJbxg97FDGwXSNba8w0Cs0iIiIi9YgRD+1uL2cS63Q28c7nAPDMtFVsKo3RLOjluiHtbK4y8yg0i4iIiNQj/mm34wyvJx1sRWTQgwDMXL2Ld+ZtAeDu0zvj96pXxPfpv4iIiIhIPeFd+m9ylv8Hy3BWtpcLxZI8MHEZAD/r1YI+rQrsLXI/LI/P7hIUmkVERETqA0fpWvzT7gCg/MTrSBX1BuB3U1eyPZKgVYNcrhrYxsYKD8DhrOjPbHcZdhcgIiIiIjXMTBGcfA2OZIRksxMp7301AFOXbWfC4m04DLj3jM7kuO1ZOCQbKDSLiIiI1HF5Xz6Ne8ucivZyp1W0l9tVnuDhySsA+HXflhzbPGhzlZlN0zNERERE6jDXpi/Im/MsAJHBD2MGW2JZFg9PWk5JNEnHQh+X9G9tc5UHYaZxb54FQLJZPy2jLSIiIiLVy4iXfttervO5xDv9CIAPFm/j4xU7cTkM7j2jMx5XZk8+MBJldpeg6RkiIiIidZJlVbSXi2wkHWxNZNADAGwNx/nd1IppGZed3JpOTfx2Vpk1FJpFRERE6qCK9nL//ba9nMePZVk88OEyIvE0PZoF+OWJLe0uM2soNIuIiIjUMY7SNfg/2d1eru/1pIp6AfD2vM18vrYYr8vBPWd0xuUw7Cwzqyg0i4iIiNQl6STBSVfjSJaRaNaP8l5jANhQEuWZaasAuOqUtrRpmGdnlVlHoVlERESkDsmb/TTurV9jeoKEhz8LDidp0+K+CUuJJk16t8znvBOa211m1lH3DBEREZE6wr3pc/Jm724vN+RRzEALAP7x1UbmbgyR53Zy9+mdcRjZNS3DcuXYXYJCs4iIiEhdYMRKCEy6BgOLWJefEu/4AwBW7ijjDzNWAzB2SDua59sfQA+Lw0my5Sl2V6HpGSIiIiJZz7LwT7sNZ2RTRXu5U+4HIJEyueeDpSTSFgPaNuSHxxbZXGj2UmgWERERyXLeJW+Rs+JdLIeL0IhxWJ6K3su/n7GGpdsi5Oe4uGNER4wsm5aRSTQ9Q0RERCSLOUpWE/jkTgDKT7yBVNMTAJi1ppjX52wA4K7TO1Po99pW41Ex07i3zAYgWdRHy2iLiIiIyGFKJwlOGoORKifRvD/lva4EoLg8wT0TlgJw7nHNGNyhkZ1VHjUjHrK7BE3PEBEREclWeV8+iXvbN5jefMKnVbSXsyyL+z9cxs6yBG0b5XHt4HZ2l1knKDSLiIiIZCH3xpnkzRkHQHjIo5iBit7Lb83dzIxVu/A4DR4c1YUctz3TGeoahWYRERGRLGPEiglMvhYDi2jX80h0GA3Aih1lPDNtJQDXDGpHx0K/nWXWKQrNIiIiItnEsgh8fCvOyGZS+W2JDKxoLxdLprnz/cWV7eV+qlX/qpVCs4iIiEgWyVn8Bt6V72M5XISHPwceHwDjpq9m5Y5yGua5ufuMTmovV83UPUNEREQkSzhLVuGffjcAZX1vJNX0eABmrNrJG19vAuCeMzrTMM9jV4k1wnK67S5BoVlEREQkK6QTBCZdXdFersVJRE/4DQA7InHum7AMgP/r3YKT2za0s8rq53CSbDXE7io0PUNEREQkG/i+ePw77eWeAYcT07K4b8IySqJJOhb6uGpgW7vLrLMUmkVEREQynHvDp+R+9QcAwqf+DtNfcZHfP+Zs5PO1xXhdDh4c1RWPS9Gupmh6hoiIiEgG26u9XLefk2g/EoClWyOMm74agOuHtKNtozw7y6w5ZhrX1q8BKpYI1zLaIiIiIrIXyyLw0c04y7aQKmhHZOB9AESTae4cv5iUaTGkQyPO7tnM5kJrliNWbHcJmp4hIiIikqlyFv8D76oPsBxuwsPHgbvibPJTH69kza4ohX4Pd4xQe7naoNAsIiIikoGcxSvxT78HgLJ+N5Fq0hOAqct38M68LRjAfWd2piDX/nZs9YFCs4iIiEimSScITBqDkYqSaDGA6AlXALA1HOfBiRXt5X7VtyUntmpgZ5X1ikKziIiISIbxzXoM9/b5mN4Cwqc9DYaDtGlxzwdLCMVSdCsKcMXJre0us15RaBYRERHJIO71M8j7+o8AhIc+jumvuMjvr1+uZ876UnLdDh4Y2QWXUzGuNql7hoiIiEiGMGLFBKZcC0C02/kk2p0BwDcbS/njp2sAuGloB1o2yLWrRHsY9v+CoNAsIiIikgksi8DUG3GWbSVV0J7IwIqLAIvLE9z+3mLSFpzepZDR3ZvaXGgtczhJtBlmdxWaniEiIiKSCXIWvY539YcV7eVGPA/uPEzL4u7xS9kWSdCmYS63D1d7ObsoNIuIiIjYzFm8Av+MewEo638rqcIeALw6a13lMtmP/KAbeR57VsMTTc8QERERsVc6TmDiVRipGImWg4gefykAX64r5k8z1wJw62kdaN/YZ2eV9rFMXFvnApBqerxt85sVmkVERERs5Pv8Mdw7FmLmNCA87EkwHOyIxLnz/SWYFvywRxGjuxfZXaZ9LAtHdGflfWyanaLpGSIiIiI2ca//hLy5LwAQHvoEpq+IlGlxx/tL2FWepGOhjxuHtre5SgGFZhERERFbGNFdBCaPBSDa/Zck2o4A4MWZa/hqQyl5bicPj+5KjlvzmDOBQrOIiIhIbdvTXq58K6kGHYkMuBuAT1fv4pVZ6wG4Y0RHWjfMs7NK+Q6FZhEREZFalrPwb3jXTMRyeAgNHwfuXLaEYtwzfgkAPzm+OSO6NLG5SvkuhWYRERGRWuTctezb9nIn3Ua6sDuptMnt7y2hNJaia1M/1w1uZ2+Rsg+FZhEREZHako4TnDgGIx0n0XIw0eMuBmDc9DXM3xzC73Xy8A+64nEpomUatZwTERERqSW+zx7BtXMRZk7DyvZy01bs4PU5GwC45/TOtMjPtbnKDONwkmg73O4qdKZZREREpDa4100j75sXAQgPexLT15QNJVHunbAUgP/r3YIhHRvbWaIchEKziIiISA0zojsJTNndXq7Hr0m0OY1EyuT29xYTiac5tlmQq09pa3OVcjCaniEiIiJSkyyLwNQbcJZvI9WgE5EBdwLw9LRVLN4aIT/HxUOju+By6lzmflkmru3zAUgVHqtltEVERETqopwFf8G7ZjKW00toxDhw5TJxyTbemrsJgPtGdqEomGNzlRnMsnCUbau431jLaIuIiIjUOc6dS/F/+ltgd3u5xt1Yu6ucBycuB+DCfi0Z0LahnSXKIVJoFhEREakJqRjBSVdVtJdrNYRoz4uJJdPc9t5iypNpeh2Tz2Unt7G7SjlECs0iIiIiNcD32cO4di7BzG1EaNhTWMADE5exfHsZDfPcPDiqCy6HTXMN5LApNIuIiIhUM/faj8ib9zIA4aFPYuUV8vc5G/lwyXacDoOHf9CVxn6vzVXK4VBoFhEREalGRvkOglOuB6D82AtJtBnGF2uLefaTVQBcP6QdvY4psLFCORIKzSIiIiLVxbIITL0eR3Q7qYadKTv5DjaVxrj9vcWYFozu3pSfHN/c7irlCKjlnIiIiEg1yZn/Kt61U3e3l3uemOXmpv/OpTSWomtTP7ee1hHD0Dzmw+Jwkmg9tPK+bWXYdmQRERGROsS5czH+mQ8CEDn5DlINO/PAxGUs233h32NndcPrUvQ6Ig6nrYEZFJpFREREjl4qSnDiGIx0nHjrocSOvZDXv3fhnxYwyW6aniEiIiJylHwzH8K1aylmbmPCQ59k1roSntOFf9XDMnHuWAxAunFX25bR1plmERERkaPgWTOFvPmvAhAe9iQbk37u0IV/1ceycEY24YxsAsuyrQyFZhEREZEjZJRvJzD1BgDKe15MqPlgbvrvQl34VwcpNIuIiIgcCcskOGUsjugOUo26Eul/qy78q8M0kiIiIiJHIHfeK3jWfVzRXm74OF7/Zqcu/KvDFJpFREREDpNzxyJ8Mx8CIDLgLmZGmujCvzpOoVlERETkcOxpL2cmiLc5jZUtf6oL/+oBhWYRERGRw+Cf+QCu4mWYuYXsOOVRbvrfYl34Vw+oT7OIiIjIIfKsmUzu/L8AEBr2JPdP38VyXfhXsxxOEq0GV963rQzbjiwiIiKSRRxlWwlMuR6A8uMu5dXtHXXhX21xeir+2EihWURERKQqlklgyvU4YrtINerGR80u14V/9YymZ4iIiIhUIfebl/Gsn4bl9LKo7++4+f2VuvCvtlgmzl1LAUg37KxltEVEREQykXP7QnyfPQzAtr53cMXUOGWJNMe3CHKbLvyreZaFM7QBZ2iDltEWERERyUjJKMFJFe3loq2Hc/nSE9hUGqNFfg6/O6s7Hl34V29opEVEREQOwD/zt7iKl5POa8LdXMHcTWH8XidPnd2Dgjy33eVJLVJoFhEREdkPz+qJ5C54DYB/tbidt5bGcRrwyOhutG2UZ3N1UtsUmkVERES+x1G2hcDUGwBY3OqX3DK/CQA3DetAvzYN7CxNbKLQLCIiIvJdlklg8lgcsWLC+V35yaoRAPysVwt+fJw6ZdRXCs0iIiIi35E790U8G6ZjOnO4IHwFkZSTge0act3gdnaXJjZSn2YRERGR3VzbF+D7/BEAnnVdyJzSQjo09vHAqC44HWotZwvDQeKYgZX37aLQLCIiIgKQjBKYNAbDTDI752SeLhlIwzw3T57dHZ9Hkck2hgHuXLur0PQMEREREQD/p/fhKl5ByNWYS0p+jcfp4IkfdadZMMfu0iQD6NcmERERqfc8qz4gd+HfsDC4ovwySgjw4Bmd6dEsaHdpYpk4i1cAkG7QQctoi4iIiNjBEdlMYOpNALyYHsVMsweXn9yaEV2a2FyZABXLaJeuxVm6Vstoi4iIiNjCMglMGYsjXsIiqy2/S/6U07sUcnH/VnZXJhlGoVlERETqrdyvX8CzYQZRvIxJXEXX5g256/TOGIY6ZcjeNKdZRERE6iXX9vn4Zj0GwL3JXxELtOUPP+yG16VzirIvhWYRERGpf5LlBCZehWEm+SB9Iu85h/Hy2T1omOexuzLJUPpVSkREROod/4x7cJWsYrPVkDtSl/LQ6G60b+yzuyzJYArNIiIiUq94Vr5P7qJ/YFoG1yd/w5gRvTi5bUO7y5IMp+kZIiIiUnelYnhXvod31YcYsRIslxfnxlkAvJAezfEnjeIHPYpsLlIOynCQbHFS5X27KDSLiIhIneRZPXF3O7lSLBwYmFiAAaQsBzmt+vDjfi3tLlOqYhhYHr/dVWh6hoiIiNQ9ntUTCY6/GCMeAsDA3H1bwWGYXLDpXrxrJtlUoWQbhWYRERGpW1IxAlPGAmCw/xXk9gSgwJSxkIrVUmFyRCwTZ/FKnMUrwTJtK0OhWUREROoU78r3cMRLDxiY9zCwcMRL8a58v5YqkyNiWThLVuEsWaVltEVERESqi3fVh1iHGHEsHHhXTajhiqQuUGgWERGROsWIlVTOYa5yW0yMWEnNFiR1gkKziIiI1ClWTgFW5SV/VWyLAyunoGYLkjpBoVlERETqDsvCMhxVzmfew8Ak3u6MGi5K6gKFZhEREakbUlECk8aQc4gX9lkYmN584u1H1XBhUhdocRMRERHJeo6yLQTHX4x72zekrIozzU7DqlzM5Pv2TN8ID3saXDm1WapkKYVmERERyWqurXMJfnAxzrKtlFg+YpaHIkcxiaITce1aipEIYWFgYH176w0SHvY0ibbD7S5fqmI4SDbvW3nfLgrNIiIikrW8y/5DYOoNGOk4KzmGpelmjHR9ScrfnNCoV7FcOXhXvo931QSMWAlWTgHxdmdUTMnQGebsYBhY3ny7q8CwrEPrEr19e7imaxERERE5NJZJ3hdP4Jv9DACfOvrwdqw3T3hewDIclP7oTZLN+9tcpGSLwsJAldvoTLOIiIhkl0QZwSnXVi5K8g/32TwfHswHObcDUN5rzN6B2TJxhNYBYAZb2fpP/HIEMmT8FJpFREQkazhCG8gffyGunYuxHB6e8F7J74v78K+cBwlQTrLpCZSfOHbvF1kWrl3LAUgEWu7/ykDJXBkyfgrNIiIikhVcm78k/4NLcER3ks4t5A7PLfxza3NuzPkfvViM6fYRGv4cON12lyp1kEKziIiIZDzv4jcIfHwrhpkk2ag7Yx038956Nyd7V3KV8RZYEBn0IGZ+G7tLlTpKoVlEREQyl5nGN/NB8r75EwDRdiMZG7+cCSvLaOSK8ZLvTxjlaWIdf0i8849tLlbqMoVmERERyUhGPERg4lV4130EQKTPWG7dNZIJK7fjdhr8p/V/yNu4nnTgGCKDHwJDk5Wl5ig0i4iISMZxlKyuuOCveAWWK4fQ0Kd4cF1X3lu0CacBf+21lpbz/4dlOAid9mxG9PGVuk2hWURERDKKe8OnBCdchiNeStpXRGjUqzy71M+bc9djAI8OCdJ3zkMAlPe+htSe1eJEapBCs4iIiGSMnPl/wT/9bgwrTbLpCYTOfIlXFsR5ddYaAG4d1oYfrroORyJEsqg35SdeV/VODQfJZn0q70uWyZDxU2gWERER+6WT+GfcQ+6C1wCIdTqH8KmP8daCXTw/Yw0A1wxqyy+S/8a9+QtMt7+ivZzjEKKMYWDlNKjB4qVGZcj4KTSLiIiIrYxYMcEJl+PZOBMLg7KTbiV6wpW8v2gbj01ZAcDF/VtxYctt5L39FACRwQ9WrA4nUksUmkVERMQ2zl3LyX//ApyhtZhuH+Hh40i0Hc7U5Tu4/8OlAPysVwuu6NOQ4Jv/h2GliXX80eG1l7NMHOGNAJiBFpqikW0yZPwUmkVERMQWnjVTCEwagyMRJh1oSemoV0k36sJna3Zxx3uLMS04q0dTxg5pR2DytThD60gHWla0lzscloVr5xIAEv7mWkY722TI+Ck0i4iISO2yLHLn/gnfzAcwsEg070fojD9h5Tbi6w2l3PTfRaRMi9M6Neb24Z3IXf4fcpa9XdFebvhzWN6g3e9A6iGFZhEREak96TiBj28jZ8mbAES7/ZzIoAfB6WH+phBj31lAPGUyoG1D7h/ZBXdkPf5ptwNQ3uc6Unu6KIjUMoVmERERqRVG+XbyP7gU95bZWIaDsgH3EO15ERgG8zaFuObf8ylLpOndMp9HftAVt2ESnHQNjkSYZLMTKe9zjd1vQeoxhWYRERGpcc4di8h//0KckY2YniCh0/9AstVggH0C81Nn9yDH7STviydwb5mN6QkQOu3ZQ2svJ1JD9OkTERGRGuVZ9QHBSddgpKKkCtoRGvkq6QbtAfhmYynXvr2AskSaPi3zefLsHuS6nbg2f0ne7GcAiAx+CDPY0s63IKLQLCIiIjXEssib8yy+Wb8DINFyEKERv8fKKQAOHJiNeIjgpKsxLJNYp3OIdzrbxjchUkGhWURERKpfMkrgoxvJWf5fAMp7XkTZgLsrp1h8s7GUa/69gPJkRWDeMyUDwD/tdpzhDaSDrYgMfvDoazEcJJseX3lfskyGjJ9Cs4iIiFQrR2QzwfEX494+D8vhIjLoAWLdf1H5/MECs3fpv8lZ/h8sw1nRXs4TOPqCDAMrr/Do9yP2yJDxU2gWERGRauPa+jXB8ZfgLN+KmdOA0Bl/ItnipMrn9wrMrQp46kfdKwOzo3Qt/ml3AFB+4lhSRb1teQ8i+6PQLCIiItXCu+wdAlNvxEjHSTXsTOmoVzGDrSqfP1hgJp0kOOlqHMkIyWZ9Ke99dfUVZpk4IlsAMP1FmqKRbTJk/BSaRURE5OhYJnmzfodvznMAxNucRvh7Uyvmbqi46K88mebEVgU8+d3ADOTNfhr31q8q2tGd9iw4nPsc5sjrs3DtWAhAwtdUy2hnmwwZP4VmEREROWJGIkJg8rV4V38IQHmvKynrd8teoXfuhlKueXs+0aS538Ds3jSLvN2BOzLkYczgMbX7JkQOgUKziIiIHBFHaD354y/EtXMJlsNDeOhjxDufu9c23w3MfVsV8MT3ArMRLyUw6ZqK9nJdfkK84w9r+22IHBKFZhERETlsrk1fkP/BJThiuzBzCykd+dI+F+59vaGUaw8SmLGsivZykY2kg62JnPLbWn4XIodOoVlEREQOS86if+CfdjuGmSTZuAehka9gBprvtU2VgZk97eX++532cv7afBsih0WhWURERA6NmcI380HyvnkRgHj7UYSGPQXuvL02m72uhOv/s+CggdlRugb/J7vby/W9gVRRr9p5DyJHSKFZREREqmTESwlOvBLPumkAlJ14PeUnXrdP+69PVu7ktncXkUhbBwzMpJMEJ47BkSwj0bwf5b2uqqV3IXLkFJpFRETkoJwlqwi+fyGukpVYrhxCw54m0WH0PttNXLKNuz9YStq0GNS+EQ+N7orXtW9P3bwvn8K9bS6mN59wdbeX2x/DQarJsZX3JctkyPgpNIuIiMgBuddPJ/jhFTjipaT9zQiNfIVU4bH7bPf2vM08Mmk5FnBG1ybcc3onXM59A4570+fftpcb/AhmoEVNvwUwDExfUc0fR2pGhoyfQrOIiIjsy7LIWfAX/NPvwbDSJJv2ovTMl7B8TfbZ9K9frufZT1YD8OPjmnHzsA44jH1XoDBiJRXt5bCIdjmPeMcf1PjbEKkuCs0iIiKyt3QS//S7yV34VwBinX9MeMij4MrZazPLsvjjzLW88vk6AH514jGMOaUtxn4Cc0V7udtwRjaRym9D5JT7a/xtfPfYjvKtAJh5TWF/9UnmypDxU2gWERGRSkasmOCEy/Bs/AwLg7KTbiN6wm/2CSqmZfHkRyt54+tNAFw5sA0X9mt1wP16l7xFzop3sRwuwsPHgcdXo+9jL5aJa9t8ABKtC8Go4TnUUr0yZPwUmkVERAQA586l5I+/CGdoLabbR3jE8yTanLbPdinT4sGJy3hvYcXZv5uGduCnJzTfZ7vK/ZasIvDJnQCU9b2RVNPja6R+kZqk0CwiIiJ41kwhMPEqHMkI6WArSke+QrpRl322S6RM7hq/hKnLd+A04O4zOjOyW9MD7zidJDDpaoxUOYnm/SvOWotkIYVmERGR+syyyJ37Ar6ZD2JgkWjen9AZf8LKbbjPprFkmpv+t4jP1xTjdho8OKorp3ZsfNDd+754Ave2b2qvvZxIDVFoFhERqa/ScQIf30rOkrcAiHb7PyKDHgCnZ59NI/EUY99ZwNyNIXJcDh7/YXf6tWlw0N27N84k96vnAQgPeXSfpbZFsolCs4iISD1klG0jf8KluLfMwTKcRAbeQ+zYC/fbmaC4PME1/17Akm0R/F4nT5/dg+Na5B98/7FiApOvrWgv1/Vn+10MRSSbKDSLiIjUM67tCwiOvwhnZBOmN5/Q6X8g2XLQfrfdFo5z1b/msWZXlAa5bp4791g6N/Ef/ACWReDjW3FGNpPKb0tk4H018C5EapdCs4iISD3iWfk+wcnXYaSipAraERr1Z9IF7fa77YaSKFe9NY9NoThN/B6eP7cnbRrlVXmMnMX/xLvy/Yr2ciNqub3c/hgGqcbdK+9LlsmQ8VNoFhERqQ8si7zZz+D74nEAEi0HEzr991je/U+zWLglzNi3F1AcTdKyIIfnf9KTZsGc/W77Xc6SVfin3w1AWb+bSDU5rvrew5EyHJpPnc0yZPwUmkVEROq6ZJTA1OvJWfEuAOU9L6ZswF3g2H8M+HTVLm59dxGxlEnnJn6ePqcHjX37Xhy4j3SCwMQxGKkoiRYnq72c1CkKzSIiInWYI7KJ4PhLcG+fh+VwExn8ILFu/3fA7f87fzMPT1pO2oL+rRvwyFld8XkOLS74vngc9/Z5mN4Cwqc9DYajmt7FUbIsjOiOiru5jTVFI9tkyPgpNIuIiNRRri1fEfzgEpzl2zBzGhI6808km/ff77aWZfHiZ2t58bN1AIzq1oQ7R3TC5Ty04Ove8Cm5X/0BgPDQ32H67f/n9EqWiXvrXAASrYdqGe1skyHjp9AsIiJSB3mXvk3go5sw0nFSDTtTOupVzGCr/W6bMi0embSc/y7YAsBF/VpyxYA2GId4Rm+v9nLd/o9EuzOr7X2IZAqFZhERkbrEMvF9/ih5uxcVibcZQXj4s1ie/beJK0+kuf29xXy6ehcOA24e1oEfH3cYZ4kti8BHN+Ms20KqoD2RgfdWw5sQyTwKzSIiInWEkYgQmHQN3jUTASjvdRVl/W854NziXeUJrnt7AYu3RvC6HDw4qiuDOzQ6rGPmLPo73lUfYDncFe3l3FW3pBPJRgrNIiIidYAjtI789y/EtWspltNL+NTHiHf+8QG3X18c5Zq357OhJEZ+jounzu7Bsc2Dh3VMZ/FK/DPuBaCs382kCo89mrcgktEUmkVERLKce9PnBD+4DEdsF+m8JoTOfIlUUa8Dbr9gc4ix7yykJJqkeX4Oz57Tg9YND/MMcTpBYNLu9nLHDCR6wuVH+S5EMptCs4iISBbLWfQP/NNuxzCTJAuPJTTy5YN2rpi+cie3vbeYeMqka1M/T53dg0aH0oP5e3yzHsO9fT5mToPMai8nUkMUmkVERLKRmcL36W/Jm/cyALEOPyA89Elw5x7wJW/P28yjk5djWnBy2wY8PLobeZ7Db9/lXj+DvK//CED41McxfUVH9h5qi2GQatSl8r5kmQwZP4VmERGRLGPESwl+eCWe9dMAKOt7I+V9rj1goLAsixdmruXlzyt6MP+ge1NuH97xkHsw73Xs6C4CU64FINr9FyTanX6E76IWGQ7MYEu7q5AjlSHjp9AsIiKSRZwlqwi+fyGukpVYrlxCpz1Nov2oA26fSJk8NHk57y/cCsAl/Vtx2cmtD7kH814si8BHN+Es20qqQQciA+450rchknUUmkVERLKEe/0nBD/8DY54KWl/c0IjXyFV2OOA2xeXJ7j5f4uYuzGEw4BbT+vI2T2bHfHxcxa+jnf1hxXt5YaPO+hUkIxiWRjxkoq73gJN0cg2GTJ+Cs0iIiKZzrLImf9n/DPuxbDSJIt6U3rGi1i+Jgd8yYodZdzwzgI2heL4vU4eHt2V/m0aHnEJzl3L8X96LwBlJ9120LCecSwT9+bZgJbRzkoZMn4KzSIiIpksncD/yV3kLnodgFjncwmf+ig4vQd8yYxVO7njvSWUJ9McU5DDUz/qQZtGR7HoSDq+u71cjETLwUSPu+TI9yWSpRSaRUREMpQR3UVwwmV4Nn2OhUHZyXcQPf7yg17w9/c5G3lm2iosoHfLfB75QTcKct1HVYfv88dw71iImdOQ8LAn1V5O6iWFZhERkQzk3LmU/PEX4gytw3T7CY8YR6LNaQfcPpk2eXTKCv47fwsAPzq2iJuHdcB9BB0yvsu9/hPy5r4AQHjo45i+pke1P5FspdAsIiKSYTxrJhOYeBWOZBnpYGtKR75CulHnA25fUp7klncX8dWGUhwGXDu4HT/v1eLIOmR8hxHdSWDydQBEe/yKRNsRR7U/kWym0CwiIpIpLIvcr/+I77OHMLBItDiJ0OkvYOUe+AK+VTvLuP6dhWwsjeHzOHlwdFcGtD3yC/6+W0tg6o04y7eRatCRyMl3Hf0+RbKYQrOIiEgmSMUIfHwLOUv/DVQsHBI55bfgPPB85Jmrd3H7e4spS6Rpnp/Dkz/qTvvGvmopJ2fhX/GumYTl8BAa8Xz2tJcTqSEKzSIiIjYzyraR/8EluLd+hWU4iZxyH7Eevz7oBX9vfr2JJz9eiWnBCS2CPHZWdwryju6Cvz2cu5bhn3EfAGUn3066cbdq2a9tDINUw46V9yXLZMj4KTSLiIjYyLV9AcHxF+KMbMb05hM6/Y8kW55ywO1TaZPfTV3J2/M2AxVLYt96Wkc8rmrqaJGOE5w4BiMdJ9FqMNGeF1XPfu1kODDz29hdhRypDBk/hWYRERGbeFa8R3DKdRipGKmC9oRGvUq6oN0Bty+NJrn1vcXMXleCAVw9qC2/6HPMUV/w912+zx7BtXMRZm4jQkOfUns5kd0UmkVERGqbZZE3+2l8XzwBQKLVYEIjfo/lzT/gS1bsKOPm/y5kfUmMPLeTB0Z14ZT2jaq1LPe6j8n75kUAwkOfOOiKg1nFsjASoYq7nqCmaGSbDBk/hWYREZHalIwSmDKWnJXvAVB+3KWUnXwHOA78V/LEJdv47YfLiKVMmgW9PPGj7nQs9FdrWUb5DoKTxwIQPfbXB+0JnXUsE/emLwAto52VMmT8FJpFRERqiSOyieD4i3Fvn4/lcBMZ/BCxbj8/4PaptMlz01fz9zkbAejbqoAHR3Wttgv+KlkWgY9uxBHdTqphZyIn31m9+xepAxSaRUREaoFryxyCH1yKs3wbZk5DQme+SLJ5vwNuv7Mswe3vLearDaUAXNC3JVcMaIPTUf3/NJ2z4C9410zGcnoJjRgHLrWXE/k+hWYREZEa5l36LwIf3YKRjpNq1IXSka9iBlsecPv5m0Lc8u4itkcS+DxO7jmjM6d2bFwjtTl3LsX/6W8BKDvpdtKNutbIcUSynUKziIhITbFMfJ8/Qt5Xvwcg3vZ0wqc9g+XZ/3xky7L49zebeeKjlaRMi7YN83jsrG60aZRXM/WlYgQnXYWRjhNvdWrdaC8nUkMUmkVERGqAkYgQmHQ13jWTACjrfTXl/W46YAu3WDLNI1NW8P7CrQAM69SYu07vhM9Tc39V+z57GNfOJZi5jQkPe1JdJUQOQqFZRESkmjlC68h//0Jcu5ZiOb2Ehz5OvNPZB9x+Y2mUW/63mKXbIjgMGHNK9fdf/j7P2qnkzXsZ2N1eLq+wxo4lUhcoNIuIiFQj98bPCE64DEesmHReU0IjXyLV9IQDbv/5ml3c+f4SSmMpCnLdPDS6Cye2alCjNRrl2wlMuR6A8p4XkWgzrEaPZzvD+HbRGJ1Nzz4ZMn4KzSIiItUkZ+Hr+D+5A8NMkSzsSWjky5j+Zvvd1rQs/jxrPX/8dA0W0K0owKM/6EpRMKdmi7QsAlNvwBHdQapRF8pOur1mj5cJDAfpBu3trkKOVIaMn0KziIjI0TJT+D69n7x5rwAQ63AW4aFPgHv/rdsi8RT3frCUaSt3AvCjY4u4cWgHvK6aX7I6Z/6reNdOrWgvN3wcuGo4pIvUEQrNIiIiR8GIlRCceCWe9Z8AUNbvJsp7X3PAf0Zevj3Cre8uZl1xFI/T4OZhHfjhsfs/G13dnDsX45/5IACRk+8k3ahLrRzXdpaFkSyruOv2aYpGtsmQ8VNoFhEROULOklUE378AV8kqLFcuodOeIdF+5H633dNO7qmPV5JIWzQNeHnsrG50KwrUTrGpKMGJYyray7UeRuzYC2rnuJnAMnFv/AzQMtpZKUPGT6FZRETkCLjXTSM48Uoc8VLS/uaUjnyVdGH3/W4bjqV4YOIypi7fAcDAdg255/TO1b8c9kH4Zz6Ia9dSzNzCiqkjOtsqclgUmkVERA6HZZE77xV8n96HYZkki/pQeuaLB2zZtmBziDveW8ymUByXw2DMKW35v94tarSd3Pd51kwhd/6fAQgNexIrr2ZWFxSpyxSaRUREDlU6gf+TO8ld9HcAYl1+SnjIw+D07rOpaVm8PnsDz89YQ9q0aJ6fw0Oju9K9tqZj7GaUbSMwdXd7ueMuIdn61Fo9vkhdodAsIiJyCIzoLoITLsWzaRYWBmUn30n0+Mv2O82huDzBvROWMnN1MQCndSrkjhEd8Xtr+a9dyyQ49Xoc0Z2kGnWlrP+ttXt8kTpEoVlERKQKzp1LyH//Qpzh9ZieAOHh4w64IMic9SXcNX4J2yMJvC4H1w9px9k9m9XqdIw9cue9gmfdxxXt5UY8r/ZyIkdBoVlEROQgPKsnEph0NY5kGelga0pHvUq6Yad9tkubFq98vo6XPl+LaUGbhrk8PLobHQp9NlQNzh2L8M18CIDIgLv3W7OIHDqFZhERkf2xLHK//j2+zx7BwCLR4mRCZ7yAlbPvEtfbI3HuGr+EOetLARjdvSk3D+tArtum1mZ72suZCeJthhPr8St76sgUhkE6v3XlfckyGTJ+Cs0iIiLfl4oR+Ohmcpa9DUC0x6+IDLwPnPu2iPt09S7u/WApJdEkuW4Ht57WkZHdmtZ2xXvxf/oAruJlpPOaEB76uIKi4dCZ9myWIeOn0CwiIvIdjrKtBD+4BPfWr7EMJ5FT7id27K/32S6VNvn9jDX8dfYGADoV+nhodFdaN8yr7ZL34lk9idwFfwEgPOwprNxGttYjUlcoNIuIiOzm2j6f4PiLcEY2Y3rzCZ3+AsmWA/fZbvXOcu4ev4Ql2yIA/OT45lw7uB1el6O2S96Lo2wrgak3AFB+3GUkWw22tZ6MYVmQilXcd+XozHu2yZDxU2gWEREBPCveIzjlOoxUjFSDDpSOfBWzoO1e25iWxZtfb2Lc9NXEUybBHBd3jOjE0I4ZsFiIZRKYcj2O2C6SjbtTdtItdleUOSwTz4YZgJbRzkoZMn4KzSIiUr9ZJnlfPoXvy6cASLQaQmjE77G8wb022xqOc/+EpXyxrgSA/m0acPfpnSj077uwiR1yv3kZz/ppWK4cwsPH7XfBFRE5cgrNIiJSfyXLCU4Zi3fl+0DFlIayk+8Ax95nsj5cvI1Hp6wgHE/hdTm4dnA7zj3Ont7L++PcvhDfZw8DEBlwD+mGHW2uSKTuUWgWEZF6yRHeRHD8Rbh3LMByuAkPeYR41/P22qY0muSxKSuYuHQ7AN2KAtx3Zmfa2Hyx316SUYKTrqpoL9f2dGLdf2F3RSJ1kkKziIjUO64tc8gffwmO6HbM3EaUnvEiqeZ999pm1ppi7vtwKdsjCZwGXNy/NRf2a4nLae/Fft/n//R+XMUrSOc1JXzq73SRm0gNUWgWEZF6xbvkXwQ+uhnDTJBq1LXigr/gMZXPx5Jpxk1fzRtfbwKgVYNc7j+zM92bBQ+0S9t4Vn1I7sK/AhA+7Wms3IY2VyRSdyk0i4hI/WCm8X3+CHlf/wGAeNvTCZ32LHi+XeZ60ZYw93ywhDW7okBFK7lrBrUlx66V/Q7CUbaFwEc3AlB+/OUkW55ic0UidZtCs4iI1HlGIkxg4hi8a6cAUNb7Gsr73QhGxVSLlGnx51nreOnzdaRNi8Y+D3ef0YmT2mTomVvLJDB5LI5YMcnGPSjrr/ZyB2UYpPf8a4Kmr2SfDBk/hWYREanTHKVryX//QlzFy7CcXsJDnyDe6UeVz68rjnLPB0tYsDkMwGmdGnPLaR0pyN13yexMkTv3RTwbple0lxsxDpweu0vKbIaDdKOudlchRypDxk+hWURE6iz3xpkEJ1yOI1ZMOq8poZEvk2p6PFBxdvkfczbwwsy1xFMmfq+Tm4d14IwuTTKmldz+uLbPx/f5IwBEBt5LukEHmysSqR8UmkVEpE7KWfg3/J/ciWGmSDY5jtCZL2H6mwGwfHuE3364jMVbK5bB7te6gDtHdKIomGNnyVVLlhOYOAbDTBJvdwaxbufbXVH2SCcqbnVWPjtlwPgpNIuISN1ipvDNuI+8+a8CEOv4Q8JDHwdXLomUyauz1vHqF+tJmxYBr4vrhrTjB92bZvTZ5T38M+7DVbKStK9I7eUOh5nGs24asHsZZkfmXdgpB5Eh46fQLCIidYYRKyH44W/wbJgOQFm/mynvfTUYBgs2h/jth8tYtbMcgCEdGnHLsA40zpBlsKviWfUBuYtex8IgfNozWDkN7C5JpF5RaBYRkTrBWbyC4PsX4ipdjeXKIzT8GRLtziSWTPOHT9fwz682YlrQMM/NTUM7MKxT46w4uwzgiGwmMPUmAKInXEHymAE2VyRS/yg0i4hI1nOv+5jgh1fiSIRI+1tQOupV0o27MXtdCQ9MXMbG0hgAZ3ZtwvWnts/ozhj7sEwCU8biiJeQLOxJWb+b7K5IpF5SaBYRkexlWeTOexnfp/djWCbJZidSesaLhJ0FPDtpGe/M2wJAE7+H24d3YkC7DO27fBC5X7+AZ8MMLFcu4eHP6UI2EZsoNIuISHZKJ/B/cge5i/4BQLTLeUSGPMT0NREemTybbZGKq+1/fFwzxpzSFr83+/7Kc22bh2/WowBETrmfdIP2NlckUn9l3zeIiIjUe0Z0J8EPLsOzeRaW4aDs5DvZ1PHXPDFhFR8u2Q7AMQU53DmiE71bFthb7JFKlhOYNAbDTBFvP5JY15/ZXZFIvabQLCIiWcW5czH571+IM7wB0xOgdPjzvB3uxrN/nkNpLIXDgP/rfQyXn9yaHHf2thbzz7gHV8kq0v5mhIc8qvZyR8MwSPubV96XLJMh46fQLCIiWcOzeiKBSVfjSJaRDrZm3sm/557PLOZtWgZA+8Z53DWiE92bBW2u9Oh4Vr5P7qJ/qL1cdTEcpAu7212FHKkMGT+FZhERyXyWRe5Xz+P7/FEMLKLNT+bxwG38+b+lpC3IdTu49KTW/LxXC1xOh93VHhVHZBOBj24GINrrSpItTra5IhEBhWYREcl0qRiBj24iZ9k7ACw75qf8avM5bFkVBmBox8aMHdIu85fAPhRmmsDk63DES0k2OY6yvjfYXVHdYaYrbrUaYHbKgPFTaBYRkYzlKNtKcPzFuLfNxTKcvBK4gt+uGACYtMjP4aZhHRjQNvvayB1I7td/wLNxJpYrT+3lqpOZxrN2KqBltLNShoyfQrOIiGQk17Z5BMdfhLNsC1FngMvj1/DJtu64nQa/PrElv+7bMqsv9Ps+19a5+L54HIDwoN+SLmhnc0Ui8l0KzSIiYo9UDO/K9/Cu+hAjVoKVU0C83enE24/Gu3oSgaljMVIxVhvHcEH59ay1iujfugE3DetAqwa5dldfvRJlBCZdjWGmiLUfTbzLT+2uSES+R6FZRERqnWf1xN1LQ5di4cDAxMKBd9UHWB/dgpGOAzA1fTzXJseQ4y/g4SHtGdapMUYdbBnmn3E3rtLVpP3NiQx5RG3RRDKQQrOIiNQqz+qJBMdfXPmzgbn37e7A/GG6N2NSY/lJ75ZcdnJrfJ66+VeWZ8V75C5+4zvt5QrsLklE9qNufgOJiEhmSsUITBkLgIF1wM1MC052LeW1nx5Lh6JGtVVdrXOENxH4+BYAynuPIdniJJsrEpEDye5mliIiklW8K9/DES89aGAGcBgQsCJ0L/24dgqzg5kmMPnq3e3ljqf8xOvtrkhEDkKhWUREao131YdYh/hXT8Uc5wk1XJF98r76PZ5NszDdPkIjxoHTbXdJdZdhYPqaYPqaaL54NsqQ8dP0DBERqTVGrKRy7nKV22JixEpqtiCbuLZ+Td6XTwAQGfQAZn4bewuq6wwHqSbH2V2FHKkMGT+daRYRkVphWRbb03mkD3V7HHXyojgjESE4cUxFe7kOZxHvfK7dJYnIIVBoFhGRGrdiRxlX/3s+r65vyqEuR2JgEm93Ro3WZQf/9LtxhtaS9rcgMuRhTRcQyRKaniEiIjWmpDzJH2eu4Z15m+lnLGSM+78AWMDBoqKFgeUNEm8/qlbqrC3e5f8jZ8mbWIaD8PBnsbz5dpdUP2TIMsxyhDJk/BSaRUSk2iXTJm/N3cSLn60lEk/zf84p3O/+My7SpPLb4ixdjYWx3y4a1u44HR72NLhyarnymuMIbcD/8a0AlPe+mmTzfjZXJCKHQ6FZRESqjWlZTFi8jRc+XcOmUBwXKZ4J/JMfJscDEOv4I8JDf4dn/XQCU8ZixEsrw3PlrTdIeNjTJNoOt/ndVCMzTWDytTgSIZJNe1He5zq7KxKRw6TQLCIiR82yLGauLub5GatZvr0MgLZ5cf4W/D0tSr4EINL/VqK9rgLDINF2BDsvmIN3xbvkLvgrJMsw89sSb39mxZSMOnSGGSDvq3F4Ns/CdPsJDX9O7eVEspBCs4iIHJV5m0KM+2QVX28MAeD3OrmuR5pfr78Pd8kaLFceoeHPkvj+RX2uHOKdzsHyFgB1d66pa8sc8r54EoDI4Acw81vbXJGIHAmFZhEROSKrdpbx++lrmLZyJwAep8F5J7TgN81XUfTx1TgSYdKBYygd+Qrpxt1srtYeRiJMcNI1GFaaWMcfEu/0Y7tLEpEjpNAsIiKHZUsoxgsz1zJ+0VZMq2LJ6x90L+LSk1rRZvXf8E26H8MySTbrS+kZf8LKa2x3ybbxf3JXRXu5wDFEBqu9nEg2U2gWEZFDUhJN8uqsdfxr7iYS6YquF0M6NOLKgW1pW+DEP+12che/AUC063kVIdHpOfhODQMzt1Hl/brEu+w/5Cz9F5bhIDT8OSxv0O6S6q86/DmrFzJk/AzLsvbt97Mf27eHa7oWERHJQNFkmn/M2chrX66nLFGxnl/vlvmMOaUtPZoFMcp3kD/hMtybv8AyHJQNuJtoz4vrdThxhNbT4I0ROBJhyk4cS3nfG+wuSUQOorAwUOU2OtMsIiL7FU+Z/GfeZl6ZtY5d5UkAOhX6uOqUtpzUpgGGYeDcsYj88RfhDG/A9AQIjfg9ydan2ly5zcwUwcnX4EiESRb1obzPtXZXJCLVQKFZRET2sics/+XL9WyPJABokZ/Dbwa0YXiXQhy7zyB7Vk2ouMgtVU4qvw2hUX8m3aCDnaVnhLw5z+He/GXFLxHDnwWH/qoVqQv0f7KIiAAQS6b5z/wt/OWL9ewoqwjLTQNeLujbkh8eW4Tb6ajY0LLImzMO36xHAUgcM5DQ6X/Aymlw+Ac103jWfVyxn1ZDsr7lnGvLHPK+fBqAyKAHMYOt7C1IKtSxz1m9kyHjp9AsIlLPxZJp3pm/hde+F5Yv6teS0d2L8Lgc326cihKYehM5y/8DQPTYC4gMuOfoFuuwzKOoPnNUtJe7uqK9XKeziXc+x+6S5LvqyOes3sqA8VNoFhGpp/aE5b98sZ6du8NyUcDLhf1b8YPuTb89s7ybo2wLwfEX4972DZbDReSU3xLr8Us7Ss9I/ml34AytIx1oSWTQg3aXIyLVTKFZRKSeiSXTvD1vM699ueGQwjKAa9s3BMdfhLNsK6a3gNAZL5A8ZkBtl56xvEvfJmfZ21iGk9CIcWovJ1IHKTSLiNQT+wvLzYJeLuzXitEHCMsA3uX/JTDleox0nFSDTpSOegUzv00tVp7ZHKF1+KfdDkD5ideRKuptc0UiUhMUmkVE6riyRIr/zNvCX2fvHZYv6teKUQcJy1gmeV88gW/2MwDEWw8jPGIclqfqfqb1hpkiOOlqHMkIyWYnUt77arsrEpEaotAsIlJH7SpP8MZXG3lr7mbC8RQAzXefWT5oWAZIlBGcci3eVRMAKD/hCsr636auA9+TN/sZ3FvmVLSXO03t5UTqMv3fLSJSx2woifK32Rt4b+FW4qmKK85bN8jllycew8huVYRlwBHaQP74i3DtXITl8BA+9VHiXX5SY/WaR9KqLgO4Nn9J3u6z8JHBD2MGW9pckRxMtn7OpEImjJ+W0RYRqSOWbo3wly/XM2XZdszd3+zdiwL8um9LBndoVLkoycG4Nn9J/geX4ojuwMxtTOmZL5Fq1qeGK88+RjxEgzdG4AxvINb5x4RPe8bukkTkKGgZbRGROs6yLGavL+G1Lzbw+driysdPatOAX/dtSa9j8jEOISwDeBe/SeDjWzDMJMnG3QmNfAUz0KKmSs9eloV/2m04wxtIB1sTGfSA3RWJSC1QaBYRyUJp0+LjFTv4yxfrWbw1AoDDgOGdC/nliS3p3MR/6Dsz0/g+e4i8uS8AEG8/ktCwp8GdVwOVZz/vsn+Ts/y/Fe3lhj+nCyNF6gmFZhGRLBJPmYxftJW/zd7AuuIoAF6Xg7N6FHF+nxa0yM89rP0Z8RCBSWPwrp0KQFmf6yjvez0YB5/3XG3MNO4N0wFIHnNKxl9o6Chdi3/anQCU972eVFEvmyuSQ5JlnzP5ngwZP4VmEZEssLMswdvfbOZf32xiV3kSgGCOi58c35zzTmhOgzzPYe/TUbK64oK/4uVYTi/hYU8R73hWdZdeJSOdrPVjHpF0srK9XKJZP8p7jbG7IjkMWfM5k/3KhPFTaBYRyWBLtob551cbmbh0O8l0xdV9Tfwezu9zDD86thl5niM74+Le8CnBCZfjiJeQ9hURGvkKqSY9q7P0Oidv9tO4t36F6QkSHv6szlaK1DMKzSIiGSZlWnyyYgf//GojX28MVT7eo1mAn53QgmGdGuOqom3cweQseA3/J3dhWGmSTY4nNPJlTF/T6ii9znJvmkXenOcAiAx5RBdIitRDCs0iIhkiFEvy3/lbePPrTWwJxwFwOgxO69SYn/VqQY9mwaM7QDqJf8a95C74CwCxTmcTPvUxcB3ePOj6xoiXEph0DYZlEuvyU1umsIiI/RSaRURstnpnOW98vZH3F24ltnsxkoJcN+f0LOLHxzWnScB71McwYsUEJ1yBZ+OnWBiU9b+FaK+r4BDb0dVbloX/49twRjZWtJc75X67KxIRmyg0i4jYwLQsPltdzD+/2rhXf+WOhT5+dkILRnQpJMddPXNmnbuWk//+BThDazHdPsLDnyPRdkS17Luu8y79Fzkr/oflcBEaMQ7Lcxit/ESkTlFoFhGpRaXRJO8v2sq/v9lc2TLOAAZ3aMTPerU4rMVIDoVn7VQCE6/CkQiTDrSkdNQrpBt1rbb9VwfLe5TTTmqIo2Q1/k92t5c78QZSTU+wuSI5Gpn6OZNDkwnjp2W0RURqmGVZfLMxxNvzNjNl2XYSu7tg+L1OzupRxE9PaH7Y/ZUP4aDkfvMivpkPYFgmiWb9CJ35J6zcRtV7nLoqnaTg7bNxb5tLonl/Sn/4hrpliNRhWkZbRMRG4ViK8Yu28va8zazaWV75eOcmfs7pWcTpXZvg89TA13A6jv/j28ld8gYA0W4/JzLoQXAefi/n+irvy6dwb5uL6c0nfJray4mIQrOISLWyLIuFW8L8+5vNTFq6nfjuC/u8LgendynknJ7N6FYUqNYpGN9llO8gf8KluDd/iWU4KBtwN9GeF+uCv8Pg3vhZZXu58JBHMQPNba5IRDKBQrOISDWIxFNMWLyNt+dtZvn2ssrH2zfO45yezTiza1MCOTX7levcsYj89y/EGdmI6QkSOv33JFsNqdFjHjUzjXvjTACSLU62/YyuESshMPkaDCyiXc8j0WG0rfVINcmwz5kcpgwZP4VmEZGjsHhrmLe/2cyHS7YRTVacVfY4DYZ3LuTsns3o2TxYY2eVv8uz6gOCk67FSJWTym9LaNSfSTdoX+PHrQ5GKmZ3CRUsC//Ht+KMbCaV35bIQLWXq0sy5nMmRyQTxk+hWUTkMBWXJ5iwZDvvLdjCsu+cVW7TMJezezZjVLem5Oe6a6cYyyJvznP4Zj0GQOKYUwid/gesnILaOX4d4l3yJjkr38NyuAgPfw48PrtLEpEMotAsInIIUmmTT1fv4r2FW5m+ahdps6IDhttpMLRjY87u2aza28VVXVSUwNQbyVn+XwDKj72QsoH3gENf7YfLWbKKwCd3AVDW90ZSTY+3tyARyTj6ZhUROYjl2yO8t3ArHyzaRnE0Wfl416Z+RncvYkSXQgpq66zydzgimwl+cAnubd9gOVxEBj1ArPsvar2OOiGdJDDpaoxUOYkWJxE94Td2VyQiGUihWUTke0rKk0xYso33Fm5l6bZI5eMN89yM7NaUUd2b0qGxff9079o6l+D4i3GWb8XMaUDojBcqLo6RI+L74gnc277Z3V7uGV0kJiL7pdAsInVPKoZ35Xt4V32IESvByikg3u504u1Hgytn/y9Jm3y6upj3Fm5hxqpdpL4z/WJQ+0aM7t6U/m0a4nLY27rNu+w/BKbegJGOk2rYmdKRr2Dmt7a1pmzm3vApuV89D0D41N9h+tVeTkT2T6FZROoUz+qJBKaMxREvxcKBgYmFA++qDzCn30N42NMk2g4HKnoqL9gc5sMl25i4ZHtGTb/Yh2WSN+txfHOeBSDe5jTCw5/D8lS9ilWms2y64M6IFROYfG1Fe7luPyfRfqQtdUjtsOtzJtUjE8ZPy2iLSJ3hWT2R4PiLATDY96vNouIs8dKBz/NG+FgmLtnOxtJv2xhlyvSLfSTKCE6+Bu/qDwEoP+E3lPW/VdMIjoZlEfzwcrwrx5MqaEfxTyeAO8/uqkTEJlpGW0Tqj1SMwJSxwP4D857HTQyKpt/E3+PPE8dDrtvB4A6NOb1LYUZMv/g+R2gD+eMvxLVzMZbDQ3joY8Q7n2t3WVkvZ/E/8a4cj+VwEx4+ToFZRKqk0CwidYJ35Xs44qVVbufAosAo4/pmi/D3+hmD2jcix52ZZ2xdm74gf8KlOKI7MXMLKR35Eqmi3naXlfWcJavwT78bgLJ+N5Fq0tPmikQkGyg0i0id4F31YeUc5qpYOLigwXxCXa6phcqOTM6if+KfdhuGmSTZuAehka9gBurgRWpmGvfmWQAkm/Wr+Skn6QSBiWMwUlESLQYQPeGKmj2eZIba/pxJ9cqQ8VNoFpGsF4olSe3cRrNDCMwABiZGrKRmizpSZhrfzAfJ++ZPAMTbjyQ07Ok6PX3ASJRVvVE18X3xOO7t8zC9BYRPexoMR60dW+xVm58zqX6ZMH4KzSKSlXaWJZi2YgcfLd/Jl+tLGOd00MRh4DSqvrbZwpGRy0wb8RDBiVfiWfcxAGUnjqX8xLEKdtWkor3cHwAID30c09/M5opEJJsoNItI1tgSivHRip18tHwHczeU7nW53zeBAZwZ//KQ9mNgEm93Rs0UeYScJasIjr8IV/EKLFcOoWFPk+gw2u6y6oyK9nLX7G4vdz6JDBt/Ecl8Cs0iktHWF0eZunwHHy3fwcIte7e+7FYUYGjHxgzp0IjWwb6Yf34VIx46YPcMqGg7Z3mDxNuPqunSD5l7/QyCH16OI15K2t+M0MhXSBUea3dZdYdlEfjoJpxlW0kVtCcy8B67KxKRLKTQLCIZxbIsVu4o56PlO5i6fAcrdnw7j80Ajm8R5NROhZzaoRFFwb1X9wsPe5rg+IuwMA7apzk87OkDrgxY23Lm/wX/9LsxrDTJpicQOvMlTF9Tu8uqU3IW/R3vqgkV7eVGPF+n54eLSM1RaBYR2yVSJl9tKGH6yl3MWLWTTaF45XNOh0GflvkM7diYwR0a08jnOfB+2g4nNPJlAlPGYsRLK8Nz5a03uNeKgLZKJ/HPuIfcBa8BEOt0DuFTH8uYMF9XOItX4p9xLwBl/W8hVdjD3oJEJGspNIuILXaWJfh09S5mrNrFrDXFlCfTlc95nAb9WjdgaKfGnNKuEfmHsYx1ou0Idl4wB++Kd8lZ+HeMZAQz2Jp4+zMrpmRkQCg1YsUEJ1yOZ+NMLAzKTrqV6AlXgpFZC6vUFqumxiQdJzDxqor2csecQvT4y2rmOJIVauxzJrUiE8ZPy2iLSK2wLItl28uYsWonM1btYuHm8F4TKBr7PAxs15CB7RrRt3UBuRm64MjRcu5aRv77F+IMrcV0+wgPH5cZZ77rIN+nvyVv7guYOQ0o/tkkTF+R3SWJSIbSMtoiYqtYMs2c9aVMX7WT6St3si2S2Ov5rk39nNKuEQPbN6RzEz+OOn6m1bNmCoGJV+FIRkgHWlI66lXSjbrYXVad5F4/nby5LwAQHvqEArOIHDWFZhGpNpZlsXZXlM/WFjNrTTGz15cQT3274IjX5aBf6wa7zyg3pNDvtbHaWmRZ5M59Ad/MBzGwSDTvR+iMF7FyG9pdWZ1kRHcRmHwdANHuvyTRdoS9BYlInaDQLCJHJRRL8uW6Ej5bUxGUt4Tjez3fxO/hlPaNOKVdI3q3zCentqZdmGncW2YDkCzqY9+yuek4gY9vJWfJWwBEu/2cyKAHwXngCxrrleoepz3t5cq3kmrQkciAu6uhSMl6mfJ9IEcmQ8ZPoVlEDkvKtFi4OcTna4qZtbaYhVvCmN+ZnOx2GpzQIp/+bRrQv00DOjT2Ydg07cKIh2w5buXxy7eT/8GluLfMxjIclA24h2jPi+rtBX8HUp3jlLPwdbyrP8RyeAgNHwfu3Grbt2Q3u78P5OhkwvgpNItIlTaHYny2ppjP1xTz5bpiIvH0Xs+3bZhH/zYN6NemAb2PqcWzyRnMuX0h+eMvxBnZhOkJEjr9DyRbDba7rDrNuWs5/k/vBaDspNtIF3a3tyARqVMUmkVkHzvKEsxZV8Ls9SXMWV/C+pLYXs8Hc1z0bVVQEZRbN9hnkZH6zrNyPMHJ12KkoqQK2hEa+SrpBu3tLqtuS8cJTrwKIxUj0XIw0eMutrsiEaljFJpFhJLyJHM2lPDluoqQvGZXdK/nnQZ0bxakf5sGnNSmAV2bBnA6NMVgH5ZF3uxn8H3xOACJloMJjXgeK6fA3rrqAd9nj+LauQgzpyHhYU+C4bC7JBGpYxSaReqhUCzJV+tLd59JLt1rqWqoWK66Y6GPPq0K6N2ygF7H5OP36uvioJJRAlNvIGfF/wAo73kRZQPuBof+u9U097pp5H3zJwDCw57UMuQiUiP0bS5SD5RGk8zbFGLO+lLmrC9h6bYI31/VqH3jPPq0/DYkH84qfPWdI7KZ4PiLcW+fh+VwERn0ILHu59tdVr1gRHcSmDIWgGiPX5Noc5rNFYlIXaXQLFLHWJbF5lCcuRtL+WZjiLkbS1m1s3yf7do0zKV3ywL6tCygV8t8GubVvRZolrPmg79r69cEx1+Cs3wrZk4DQmf8iWSLk2r8uHXJEY+TZRGYeiPO8m2kGnQiMuDO6i1M6pTa+D6QmpMJ46dltEWyXNq0WLmjjLkbS5m7McQ3G0v3WXkPoHWDXI4/Jn/32eT8+rOwSA3yLnuHwNQbMdJxUg07UzrqVcxgK7vLqjdyFrxGYNrtWE4vxee+S7pxN7tLEpEspWW0ReqgWDLNwi3hyrPI8zaFKEvs3QLO6TDo2tTPcc3zOb5FkONaBGlQB88k28Yy8X3+GHlfjQMg3mY44eHPYXn8NhdWfzh3LcM/4z5gd3s5BWYRqWEKzSIZzLQs1hVHWbg5zILNIRZuCbNsexlpc+9/IPJ5nBzbPMjxLYIc3yKf7kUB9UquIUYiQmDytXhXfwhAea+rKOt3s1YYq02pWEV7uXScRKshRHuqvZyI1DyFZpEMUlKeZOGWMPM3h1i4OczCLWHC8dQ+2xX6PRzfYs9Z5Hw6NPapBdz3mWlcW78GINX0hGoJtY7QevLHX4hr5xIsp5fwqY8R7/zjo95vvXYE4+T7/BFcOxdj5jYiNOwprbAoVauB7wOpRRkyfgrNIjZJpEyWbY+w4DtnkTd8bxERAK/LQZcmfro3C9CjWZAezQIUBby2LU2dTRyx4mrbl3vTLIIfXIojtgszt5DSkS+RKupdbfuvzw5nnNxrPyLvm5cACA99EiuvsKbKkjqmOr8PpPZlwvgpNIvUgkTKZMWOMpZsDbNkW4QlWyOs2FFGMr3vdbitG+TS4zsBuUNjHy6nFmqwU86if+CfdjuGmSRZeCyhM1/GDDS3u6x6xyjfQXDK9QBEj72ARJthNlckIvWJQrNINYsl07sDckU4Xrw1zMqd5fvMQwYoyHXTo1mA7kUBejQL0K0oQDDH/rY6spuZwjfzgcozm7H2owkPewrcuTYXVg9ZFoGPbsQR3U6qYWciJ99hd0UiUs8oNIschWgyzbJtEZZui7B4d0hevbOM/ZxAJj/HRdemATo39dO1qZ/OTfy0yM/RNIsMZcRLCU68Es+6aQCU9b2B8j7Xaf6sTXIW/AXvmslYTi+hEc+DS7+4iEjtUmgWOQSmZbE5FGP5tjKW7yhj+fYyVmyPsKEkts/KegANct102R2OuzQN0KWpX/OQs4izZBXB9y/EVbISy5VDaNjTJDqMtrusesu5cwn+T38LQOTkO0g36mJzRSJSHyk0i3xPWSLFyh3lLN8eYfn2ioC8ckfZPr2Q92jk81SE4yZ+uuwOyU38HgXkLOVeP53gh1fgiJeS9jcjNPJVUoU97C6r/vpOe7l466HEjr3Q7opEpJ5SaJZ6K5EyWVcSZdWOMlbvLGfFjjJW7CjbbwcLALfToG3DPDo28dOxsY+OhT46FPrq5PLTdYZxGBdQWhY5C/6Cf/o9GFaaZNNelJ75EpavSc3VJxUOMk6+zx7CtWspZm5jwkOf1PQYOXKH830gmScDxk/LaNdxzsWL8D18P2W33U26a/1cMSuRMllbXM7qneWs3Flxu3pnGeuLo/udewwVfZA7Fvro0NhPx8KKgNy6Qa66WNRV6ST+6XeRu/BvAMQ6n0t4yCPgyrG5sPrNs3Yq+e/9CoDS0a+RaD3U5opEpK7SMtqCZ9KHeCeMJ3lif6J1PDRHk2nWFUcrQ/GqneWs2lnOhpIo+2lcAVSspNeukY92jfJo1zivIiA39lOQpw4W9YURKyY44TI8Gz/DwqDspNuJnnCFzmjazCjfTmB3e7nynhcrMIuI7RSa6zjP9I8rb6NXX2djJdUjbVpsCcdYuyvK2uIo63aVs7Y4ytpd5WyLJA74Or/323DctlFeRUhu5KNQc4/rNefOpeSPvwhnaC2m20d4xPMk2pxmd1liWQSmXI8juoNUo66UnXSb3RWJiCg012nJJO5ZnwHg/vwzSCbBnR1nUEujycowvLY4yrrd9zeUREkcaE4FFW3d2jSsOGvcds8Z5EZ5NPYpHNc7lolr61wAUk2P32c+nGfNZAITx+BIRkgHW1E68lXSjTrXfp313X7GKXfeK3jXfVTRXm74OE2TkaNXxfeBZLgMGT+F5jrM9c3XGLGKi9qMWBTXvLmkep9oc1UVTMtieyTBhpIoG0tibCiNsr44xsbSKBtKYoTjqQO+1u00aFmQS6sGubRumEfrBt/eL8jNjl8KpBZYFo7ozsr7GN8+nvv1H/F99hAGFonm/Qmd8Ses3Ia2lVqvfW+cnLsW4/vsIQAiA+7SLzJSPQ70fSDZIUPGT6G5DnPPnIHlcGKYaSyHA/enM2o1NCdSJlvCcTaURHf/iVXclsbYVBojnjIP+vomfg+tG+btE46bBXNwOvSNJweWSMeZtnkypSv/SXmqnFU7p3Bys8EMaTyARtPvJmfpvwCIdjufyKDfglMdUOzw/XFav/1D7l3ySUV7uTanEevxa7tLlDrggN8HRUPxOL12lydVyKTxU/eMOiz/Jz/E/cnHGJaFZRgkB59K6Zv/qbb9J1ImW8NxNoVibC6NsTkUY1MoXnl/eySx34U/9nAa0Cw/h2Pyc2lRkMMxBbm0LMihRUEux+TnkON2VlutUn98unU6j37zAOXJEIOicSxMZuT6KEineG77Lo6NRbEMJ5GB9xI79gJd8GeT/Y3TgPI4Pw+H2el08cWZz9On9Si7y5Qsd6Dvg7Rh4XcFuPW4uzi56UC7y5QDqM3xO5TuGQrNWcqxeROO7dsOvIFlUfCD0yunZwBYOTmUvPvhQUOCWdgEs1lzAMoTabaG42wNx9gSih92KAbwuhwcU/BtMG5ZkFvxc0EuRQGvWrhJtfp063TunnMrAIZlcko0CsB2h5Ont+2gWTpNyOFg/oBb6NTzKjtLrdf2N07d4wl+UxIC4IqmTZiZl8v9vR9hQNNT7CxVstiBvg+m5+ZiGgbG7n/j1+csM9X2+Ck012HBn/8Y75RJB93GMgyM7wzv93/en697nMxNv36QreH4QecV7+F1OWgezKFZvpdmwZzd93NoHvTSLD+HBrluXYAntSKRjnPulLMoS0WwsHBYFqdEoxwfi3N+KEKuZbHa7eKapk3YkduAfw37n/5p1gb7G6eRkTJu31lMwLJ4LRjgd40aYGDgc/k1TnJEDvR9AN+GLkCfswxlx/hVa5/msrKyoypGqpd5/q9puGwp7vXrsNj/nPjvB+QDBeY9r1+b35SnOwxh2aadlc/5PE6aBDw08VeE4qKgl6JgDs0CXoryvRTkHCQUW0nKy5NH9P5EDtfUrZOIpL7zy71lcUaknNFl5QDMyM3h5sLGhJ0OSIWZuHYCpzZVe7na9v1xMiyLX4TCBCyLJW43TzcsAMDCIqJxkiO0z/fBAehzlpnsGL9qPdOss4WZxw1cBdwP5HJkV3WmgHLgbuD3gCKuZKuWY1oS7BXE2H2RaLdYnDc2bwXgr0E/TzRsQHr395hlWoS+CrF+3Hrb6q2vvj9Ox0Vj/G3LNhLAec2LWOH99qJMjZMcqe9/zg50phL0OctEdozfocRhdc/IYkngaeCvVATnK4A0FWH6UF7rBF4A7gF2HnxzkYzn9Dkrv2ABVnvc/DPg5+scL+P9vr22NRwGTp8uNLXD98dpqdfDvwI+publ7RWYQeMkR+77nzPTMJiWl7ffbfU5yzyZOn6HfKZ5zZotNV1L1rMsi1jKpLg8QUk0RUk0WfmnuDxFaSyx+/6ex1NVtl37voDXRYM8F419Hhr7vTT2eWjk99A4z0PbrWvo/djd+BfMPejcZQuIn9CLnQ8/QbJzl6N81yKZ4aGF9/H5jk+xqrw8tWIeXP/GA7i9+z21UJl8l8ZJaoM+Z9nNjvFr06aoym0O+Uyzz+ereqM6xLIsypNpSqMpSmNJSqNJSncH4YqfU5W3JdEkxbvD8aGHYAc4PDg8FfOGG+a5aeTz0DDPQyOfh0Y+97f3dz/XIM+D13WwbhOtcb/fDRbPh9RBLuJzubA6d8PTqzfqTit1xeAWp/LZjhmHtK2FxZAWQ+vd91om0DhJbdDnLLtl6vjV+e4ZqbRJOJ4iFEt9extLEYpX3JbGkoR3P1caS1WE4923KfOQ/tPsw+ty0CDXTYM8NwW5bhrmuSnI9dAgz135+J7nGvk85FZXP+JIhMbd2u3VZu5ArJwcdixaBX5/9RxbxGbfv9r6QHS1vL00TlIb9DnLbnaMX7V2z7DDnukOkXiKSDxdcZv4zv14ikgiTSS29+PfDcflyfRR1eBxGhTkusnPdZOf49p96yY/17XXbcM8NwV5bhrkesh1O2y5cNL73n/3Ccx72sx9v92cEYvhff9/xM/7v9ouU6RGeJxebj3uLu6acwsGxn6/aPf09bz1uLv0F6RNNE5SG/Q5y26ZOn7VfqZ5T9CNJtOUJyr+lFXepiir/Pm79yue27NtWSJN2e5AnD7Cs73f5/M4Cea4CHhdFbc5boJeF4Ec116P5+e6KdgThnPd5LjsCcBHIv+HZ+Ke9RmGWTFFxHK6wOshevHl5L78AsTjGOmKXyIsh4Nkv5Mp/e94O0sWqXZ7VpCKpMIYOLAwK2+1Aljm0DhJbdDnLLvV5vhVa8u5e96eRzSRpjyZrgzEe90mTaK771dPzP2WwwCfx0XA68TndeH3uvB7nBW33orH/V5XxXMe5+4g/G0o9ntduBzZEXyPlGPdWhr1ObbyZwuDVO8+hF54BbNVaxxr1xC8/CJcX83B+M4I7ZyzALNlKztKFqkxiXScaVs+YsaWTwglSwm68xlYNIjBRafqjFIG0ThJbdDnLLvV1vhVa2huc+v7h11ArttBrrsi0Oa5nfi8zt23LnweJz6PkzyPkzzPtz/7vvNzRSiueE22nO21S97jj+B77CEshwMsi/Ibb6V87E3g+s4MnFSKvKd+R97jj4BhYJgmZbfcQfkNt9hXuIiIiIjNqjU03/bm3IqA63aS43aS53GQ53aSu/uxXHdFAM7dfT/H7cChoFs7TJOGvbrj3LSRdPMWhP70Z1J9+x1wc9cXswhe9mucmzaRbnEMu75aCBorERERqaeqNTRna/eMeiEapeDHPyDdug2RR5/ACuZX+RIjVIr/5utxrltLyb/fhdzcWihUREREJPMoNNcnpgmOg/VwrubXiYiIiNQRhxKalZbqiiMNvgrMIiIiIlVSYhIRERERqYJCs4iIiIhIFRSaRURERESqoNAsIiIiIlIFhWYRERERkSooNIuIiIiIVEGhWURERESkCgrNIiIiIiJVUGgWEREREamCQrOIiIiISBUUmkVEREREqqDQLCIiIiJSBYVmEREREZEqKDSLiIiIiFRBoVlEREREpAoKzSIiIiIiVVBoFhERERGpgmFZlmV3ESIiIiIimUxnmkVEREREqqDQLCIiIiJSBYVmEREREZEqKDSLiIiIiFRBoVlEREREpAoKzSIiIiIiVVBormf2dBh87rnn9vpZJNu99957HHvssWzdurXysTvuuIMf/OAHhMNhGyuT79I4SW3Q5yy7Zer4KTTXM3//+9954403KC8v5/HHH+fLL7+0uySRajFq1Cjatm3LCy+8AMC4ceOYMWMGL730EoFAwObqZA+Nk9QGfc6yW8aOnyX1zosvvmh1797dmj17tt2liFSrqVOnWj169LD+8Ic/WCeeeKK1bNkyy7IsKxwOW+ecc4511llnWaNHj7beeOMNmyut3w40Tps2bbJ+8YtfWGeeeaY1evRoa/z48TZXKtnsQJ8zy7KsK6+80urTp4919dVX21ihHMzBxm/q1KnWiBEjrOHDh1tvvvlmrdWkFQHrmX/84x84HA7WrFmDy+Vi8ODB9OnTx+6yRKrN2WefzfLly3nllVfo27cvAOl0mkQiQW5uLtFolNGjR/Ovf/2LBg0a2Fxt/bW/cdq2bRs7d+6ka9eu7Ny5k7PPPpsJEyaQl5dnc7WSrfb3OQP4/PPPKS8v5z//+Q/PPvusjRXKwexv/FKpFKNGjeK1117D5/Nxzjnn8Oabb1JQUFDj9Wh6Rj3zs5/9jPPOO4+8vDxuuOEGevfubXdJItVm+vTprFq1CtM0adSoUeXjTqeT3NxcAOLxOKZpaj6/jQ40Tk2aNKFr164ANGrUiPz8fEpLS+0qU7LcgT5nAP3798fn89lUmRyKA43fvHnz6NChA02bNsXv9zNo0CBmzJhRKzUpNNcRhzpp3jAMAK6++uq9fhbJdgsXLuTaa6/lvvvuY+DAgTzzzDN7PR8KhTjrrLMYPHgwl1xyCQ0bNrSp0vqtqnHaY/78+ViWRbNmzWq5QqkLDvVzJpnpYOO3bds2mjZtWvlzUVHRXtmnJik01xEZO2lepBZs2LCByy+/nMsuu4wf/ehHXHPNNUycOJEFCxZUbhMMBvnf//7HlClTePfdd9mxY4eNFddPhzJOAMXFxdxyyy3cf//9NlUq2exQP2eSmaoav/39K2FtnQBUaK4jDMNg7NixvPXWW/zxj3/ktdde46WXXtrrtzGRuqikpIRLL72UoUOHcsUVVwDQo0cPTj31VJ566ql9tm/cuDGdO3dW55hadqjjlEgkGDNmDJdddhm9evWyq1zJUof7fSCZ5VDGr2nTpnudWd6yZQuFhYW1Up8uBKxjDnTRg0h9tmPHDnJycvD7/UQiEc477zyeeOIJunTpYndp8h2WZXHDDTfQtm3byilkIjVl1qxZvP7667oQMMukUilGjhzJX//618oLAd94441aubDbVeNHkFpzsIseROqzLVu2cMcdd2BZFpZlcf755yswZ6A5c+Ywfvx4OnfuzOTJkwF47LHH6Ny5s82VSV1z8cUXs3DhQqLRKIMGDWLcuHH07NnT7rLkELhcLm655RZ+9atfYZoml1xySa11QtKZ5jpi4cKF/PKXv+Tuu+9m/Pjx5OTk6LdnERERkWqiM811wPcnzXfo0IFzzz2XBQsW0KNHD7vLExEREcl6OtOc5UpKSvj5z3/OiSeeuNeV5r/5zW9IJBK8/PLLNlYnIiIiUjcoNIuIiIiIVEEt50REREREqqDQLCIiIiJSBYVmEREREZEqKDSLiIiIiFRBoVlEREREpAoKzSIiIiIiVVBoFhERERGpgkKziIiIiEgVFJpFRERERKqg0CwiIiIiUgWFZhERERGRKig0i4iIiIhU4f8BUUEYbX6cQcMAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from turtle import title\n", "\n", "\n", "def f(x):\n", " return x**5 - 3\n", "\n", "xmin, xmax = 1.0, 2.55\n", "xinit, xstar = xmax-0.05, 3**(1/5)\n", "xx = np.linspace(xmin, xmax)\n", "yy = f(xx)\n", "\n", "n = 4\n", "x = np.zeros(n)\n", "x[:2] = xinit, xinit-0.25\n", "y = f(x)\n", "for i in range(2,n):\n", " x[i] = x[i-1] - y[i-1]*(x[i-1]-x[i-2]) / (y[i-1]-y[i-2])\n", " y[i] = f(x[i])\n", " \n", "\n", "xticks=np.array(x[:4].tolist() + [xstar])\n", "\n", "fig3, ax = plt.subplots(figsize=[9,6])\n", "ax.set(\n", " title=\"Secant Method\",\n", " xlim=[xmin,xmax],\n", " xticks=xticks,\n", " xticklabels=['$x_0$', '$x_1$', '$x_2$','$x_3$', '$x^*$'],\n", " yticks=[]\n", " )\n", "\n", "ax.plot(xx,yy)\n", "ax.axhline(0, color='k')\n", "ax.annotate('f', [xinit,f(xinit+0.03)], fontsize=18,color='C0')\n", "ax.plot(x, y, color='C1', **bullet_options)\n", "ax.plot(x, 0*y, color='C2', **bullet_options)\n", "\n", "ax.vlines(x, 0, y, colors='C1', linestyles='--', alpha=0.3)\n", "\n", "for xi,xinext,yi in zip(x,x[2:],y):\n", " ax.plot([xi,xinext],[yi, 0], color='C1')\n", "\n", "ax.plot(xstar, 0, color='red', marker='*', ms=16)" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3.9.7 ('base')", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.7" }, "vscode": { "interpreter": { "hash": "ad2bdc8ecc057115af97d19610ffacc2b4e99fae6737bb82f5d7fb13d2f2c186" } } }, "nbformat": 4, "nbformat_minor": 2 }