{ "metadata": { "name": "", "signature": "sha256:f409417a38ec225c6094658c4ccb0e4ddf8ba0010090283e913b52f8c5cfa7a8" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Periodic SuperSmoother\n", "\n", "One feature of the SuperSmoother is its ability to handle periodic signals. We'll demonstrate this here:" ] }, { "cell_type": "code", "collapsed": false, "input": [ "%matplotlib inline\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "# Use seaborn for plotting defaults.\n", "# This can be safely commented-out if seaborn is not installed\n", "import seaborn; seaborn.set()" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "# Create a periodic function with interesting structure\n", "\n", "def func(t, period):\n", " phase = t / period\n", " modphase = phase + 0.1 * np.sin(2 * np.pi * phase)\n", " return np.sin(2 * np.pi * modphase)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "rng = np.random.RandomState(42)\n", "t = 50 * rng.rand(100)\n", "period = 10\n", "dy = 0.1\n", "y = func(t, period) + dy * rng.randn(len(t))\n", "\n", "phase = t % period" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "def plot_data():\n", " fig, ax = plt.subplots(2)\n", " fig.subplots_adjust(hspace=0.25)\n", "\n", " # Plot the data first\n", " ax[0].errorbar(t, y, dy, fmt='.')\n", " ax[0].set_xlabel('t')\n", "\n", " # Plot the phased data below\n", " ax[1].errorbar(phase, y, dy, fmt='.')\n", " ax[1].set_xlabel('phase');\n", " return ax\n", "\n", "plot_data();" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAesAAAFkCAYAAAAAFROsAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X2QXFd55/HvzOgNhtYgmwbCkoVs4jwF2TWxrQpGdtnI\nCXKAaLFZKlmHKAvIL7Gy8dqmJBzLQIWywXjWAmuD2GCLKHgpbfHmsMYb4wJPWbbkhZhVYnsDj2Nv\nwhYEyCDLM8PYM9LM9P5xu0dXrX69fW/fl/59qqaqp2/3vWfOdN/nnnOec+5QpVJBREREsms47QKI\niIhIawrWIiIiGadgLSIiknEK1iIiIhmnYC0iIpJxCtYiIiIZt6KXN5vZG4Fb3X1j3fPXAVuByepT\nV7n7U70cS0REZFBFDtZmtgP4PeBnDTafDWxx98NR9y8iIiKBXrrBnwbeCQw12HYOcKOZPWxmN/Rw\nDBERkYEXOVi7+1eAhSab9wNXARcB55vZ26MeR0REZND1NGbdwh3uPg1gZvcBZwH3NXtxpVKpDA01\naqCLiIgUUldBL/ZgbWZjwONm9nrgeYLW9d5W7xkaGmJycibuokhIuVxSHfeB6jl5quPkqY6TVy6X\nunp9HMG6AmBmlwEvcfc7q+PUE8A88A13vz+G44iIiAykoYzcdauiq7hk6Uq5P1TPyVMdJ091nLxy\nudRVN7gWRREREck4BWsREZGMU7AWERHJOAVrERGRjFOwFhERyTgFaxERkYxTsBYREck4BWsREZGM\nU7AWERHJuFwH6+17DrF9z6G0iyEiIpKoXAdrERGRQaBgLSIiPVEvZ/KSup91X0zNzqddBBERkcTl\ntmV9y92PsbBYYWGxwi13P5Z2cURERBKT65a1SJ7Vug3Ht21IuSQivVEvZ/J6almb2RvNbKLB85vN\n7NtmdsjMLu/lGGHhcZGdW9azYmSIFSND7NyyPq5DiKRO43+SJ+rl7I/IwdrMdgB3Aqvrnl8J7ALe\nAlwIXGlmL++lkM2Mja5mbHR1+xeKZNDU7DxTs/MKzglQnUrR9NKyfhp4JzBU9/zrgKfdfcrdjwOP\nABf0cBx98aRwwq2R+i7EWhAXyYNwL+dzM8d0rk5I5GDt7l8BFhpsWgtMhX6fAcaiHies/iQ2vm1D\nw/E+BXfJk3DAVpdiPHTB01/q5UxeEglmU0Ap9HsJONruTeVyqem2qdljLCwuLf9+2/7DjF/TvLE+\nMjLUdp+DSPXRH53U8yev38ilO+5d/lwvLFa4bf9hVq4YWX7NyhUj+p810apetu8+wMJiBWh8rth6\n8wMA7L1pU3IFLIBuPnu1c+6z0/Ndv1c6k0Sw/h5whpmtA2YJusDH271pcnKm4fNBS2PppOeOLyw2\nfT3AYvWL2uo1g6ZcLqk++qCbeh4bXcXU7PxyYDm+sMjOLeu5cjzI2dxx2Vn6nzXQro6PLyye9Lj+\ntTo/tNft+eLWq97ELXc/xj8ffQGAa3dNKPG3jW4vaOKYZ10BMLPLzOyK6jj19cDXgUPAXnf/UQzH\nAego+1tdYJIXY6OrT5nVoC7F3rSbKaLzg+RRT8Ha3f/R3TdUH+939zurj7/m7r/m7uvd/dO9HCP8\nxTt97Zq2JzGN+UneKDjHr1md6vyQHE2nTVYuFkXp5kT2/R837rrRAhSSNbXPopIhpSh00ZmcXATr\nmk4C7djoao5Mz+nqTnJLF5TJCecE6PwgeZKrYN2p4SFd4Ul+KDjHr1Wd6twgeZSLYN3NyWxqdp6l\nSuPnRUQkunbDiVEuPDVE2Znc3nWrkVryCMBLS6tOeV5JJcnTgjQiIvHLRctaRCQuasFJHhWqZd1s\n6oCmFIiI9C6JOeqa996ZQgVraD6/UnNZ+0NfPJFiSmI4UUOUnStcsJbOxT2+rC+eiEgyCjdm3S5L\nUZmHIsWi73R/JDFHXfPeO6eW9QCLu8tauQEixZbEcKKGKDtTuJZ1OxpPDYSnud1y92OxBVd96aRf\nNEVQBslABetwgLpyfILPbN+YcolEuqMu31PpArx/kvjc6bPcGXWDD6ikuqzHt23Ql0/6Ymp2niPT\nc0pqlIEwUMG6FqBA3bWgsaI80tS4QLiXTGQQROoGN7NhYA9wJjAPXO7uz4S2XwdsBSarT13l7k/1\nWNZY1O7KpROe5E1SeQZFkFRSo4YdJCuitqwvAVa5+wbgBuD2uu1nA1vcfWP1JxOBGk6Mb6nbTCS/\nwsM4p69dox4iKbyowfo84H4Ad/8WUH9Jew5wo5k9bGY39FA+SZDGl/OlXZ7BoN1ERcM4xTJon99u\nRc0GXwtMh35fNLNhd1+q/r4f+BQwA9xjZm939/ta7bBcLkUsSnfuue3fcumOewH45PWDlQ3erzoe\ndEnW82lr1zQ9xkg1H2MQ/s/lcmn57w0/F7ep2WOJ7Tvran/z1psfAGDvTZsSPd4gfX6jiBqsp4Fw\njYYDNcAd7j4NYGb3AWcBLYP15ORMxKJ0b2x0Vd+PmbZyuTRQf29akq7nxeqYdaNjtNpWJLU6rv29\ntd6huP/uIEcgOK1du2tioHIEwp/jfn2unp2e68txsqLbi5KowfogsBn4opmdCzxe22BmY8DjZvZ6\n4HngImBvxOP0xeUffxCAuz5wUcolKSYl6fTHoCVN6vPUH/34XCl5sr2oY9b3AHNmdpAguew6M7vM\nzK5w9ymCpLMJ4ADwpLvfH09x46GxWsmrZp9d3UQlflo+V5+rLInUsnb3CnB13dNPhbbvJxi37kpa\nLbAlTdcUkQaUwNYfuqFHewO1KEojV+96qOFjiY8W8kieWoGShH5+rpTd31qm1gbPwwld46/d0VhU\n/+hEFz99z/W5yorMtKzTGhv59PUXNnwsydF8ShGJYpDPHZlqWadleKj9ayQajUX1j1qBkmf6/LaW\nmWCdl5N6Hrrqs6a+G011KJIfCqLZkJlucMh+goGmMfROdSgiUQ1ysmpmWtZp0mIoydKVuRSdEk+T\nN+jJqplqWWedpsf0TnUoItK9TLWs83BVmuVu+rxQHUrRFK1rduvND7C4WMnUOTkveU1JUctaMmGQ\np2RIvikPo3+ynteUpEy1rPMgS1eaIiJJeHZ6HtA6zFmiYC19pwseKZKidc+Gbw2atUSuQT53KFi3\noSzP/ijamJ8Mlm67ZnVekW5pzFpSpzE/kewIZmwM52LGxiDlukRqWZvZMLAHOBOYBy5392dC2zcD\nHwQWgM+6+10xlDWTdIXcP6rreKgepZ3T1q5mcVFj1lkStWV9CbDK3TcANwC31zaY2UpgF/AW4ELg\nSjN7ea8FTcsgr5jTL+Gr96xfyUu6tu85xNabH0i7GKcY37ahq4sfnVeSVcQWd9RgfR5wP4C7fwsI\nn2FfBzzt7lPufhx4BLigp1KmpCjds1k9wYUNDzW+oUr4S6cTXDyOTM9xZHou7WIMpO17DnHl+EQh\nzitZMEjnhKjBei0wHfp9sdo1Xts2Fdo2A4xFPE7mDdKHJU1FuXBKW7ju8liPU7Pz1WlFkqS9N23K\n/DDJoJ0TomaDTwOl0O/D7r5UfTxVt60EHG23w3K51O4lfffJ6zdy6Y57lx/X2777wPJatbftP8z4\nNdnsQJiaPQZks45rhoeDa736Mi6Xfd2Llp9buWIk039Llsu2csXISY+zXNZ6J75vlUx/31oZGRni\ntLVrli84Gp1XsiLrn41Wn+U8nPO6FTVYHwQ2A180s3OBx0PbvgecYWbrgFmCLvDxdjucnJyJWJRk\njY2uAhqX7/jC4kmPs/g3hOdMXrtrIpNjws3KGH7++MIiK0aCfvIdl52VybqG4OSQ1bJBUHfvu/XB\n5cdZLmu9PHzf2nm2OvxQm+qV1b8h659jCD6/tfnt4c9yHs550P2FRNRu8HuAOTM7SJBcdp2ZXWZm\nV1THqa8Hvg4cAva6+48iHifTdFOK/hrkpQbjkudu8BPft+Fcft/C3bYaOovHIJ0TIrWs3b0CXF33\n9FOh7V8DvtZDuTKj3bhNow9KlqbGPDcTdAdl+QTXbAWonVvWL7cCd25ZX7jsTune2OhqRkYaZCLm\nzGteWcrs9zHviraiXI0WRRkAw0PBvMks6+UKuYjTNJKiaXLpUd3Hr9mUuSK2uLXcaI+y0HpuZWp2\nnqUK/HQqf1N16rtss17XeXH62jVpF0FEuqRgnYCsjEfVxsgAlpbyvxpRloYXJB3j2zbkIvmpkfqL\nT7Wuk1PEc4S6wWM2aHP/4tKoO6vTBD4t8tGdTlbbyuLQQhbLJNIvCtYx+/6Ps3PFXx/g8njx0G7s\nKc/ZzdKdozNzHJ3J70WZZo+kL88XfArWMasFlhUjQzw3cyz1D8aKAmTO1mi1uHTk+QSXNZ0mPqnO\npZ6CdcxqwSQrmYgnLh6yO3WrlVqXbbPhBWXYJiOLF0ZLleBHZBApWMconNBVO9ll4YSXh6lbUakb\nPH5ZzLu4etdDy49/+8b8LuHQ7d25JF5ZOSdHoWCdkJGR4Uyc8Ma3bWBdKf9TdTTel56p2XmOTM+p\nW7aPlDQZvyxehHZDwTpG4YDy6vJo2sUpnEbjfQri8QvXKXBSb1Eatu85xEvWrFz+/Qsf/a1UytEv\n6i2SRhSsY1YLKFkKIuPbNrD3pk2pliFJRVytKG21Os3S7IZm9zwX6USWzslRaFGUBCmAxEtjff03\nNrpa3bF9Vr8mvsQnz+dktawl95S0E79anYa7vtMa68tzUlAU4Tqu3ZBCRC3rmKUZNLQcp8QpPLsh\nC2X4xX+xVi1N6Umez41qWSdILT4pkryO9eVNbWxVJKzrYG1mLzKzL5vZATO7z8xe1uA1d5jZY2Y2\nYWYPmtnaeIorIv0STsg5fe2aVMb78p4U1Ku8TjOS+EXpBr8a+Ft3/4iZ/Q5wE3Bt3WvOBja5+7O9\nFlA6N0jjetIfnQboJIdg8pwUJBLXdyNKN/h5wP3Vx/cDvxHeaGbDwBnAnWb2iJm9t6cSSkfyPuE/\nLlpTORka0umvz2zfONA9CnKqli1rM9vKqa3mnwDT1cczwFjd9hcDu4Fd1f1PmNlj7v5E78UtDiWD\nSZGoVyd+6lHor6yfk1sGa3ffC+wNP2dmXwZK1V9LwHN1b3se2O3uc9XXPwi8AWgZrMvlUqvNhTNS\nTSCJ6+/+5PUbuXTHvcuPGxmEOp6aPQak+7cWqZ73ffjitq/ZvvvAcsb2bfsPM37NBYmXoUh13Ezc\n54huDUIdhyVV33Gdk6KMWR8E3gb8NfBW4EDddgP2m9nZwAhwPrCv3U4nJ7OzUlI/PFtdaCLOv3ts\ndFXTfZbLpcLXcTAUsATAtbsmEu0+bHYVPgj1XO/4wuJJj5P++weljherF0Bp/K2DUsdhSdR3q3NS\nt8E7ypj1p4FfMbOHgcuBPwEws+vMbLO7fxf4HPAoMAHsqz4nVUmML2/fc4ijM1ppSvpv0DO2k6I8\ngf7K+uI7Xbes3f0F4LcbPP+J0ONdBGPW0kfrSmsG+su9c8v65RWfkg4aWf5Sp0Hjq5Jn4cV3brn7\nsdjOH3Gek7QoSgqSaIlk/aqwX/pxUw9l3jenbHyRk8V1TlKwTkmcQUXBQ0S6pQurE/IwlKO1wUW6\n1M/u9ryoDb/oxhOSV0n1yMU1NKmWdQHk4aqwX/qVlKN7aJ9KPTwiyVHLOiVxBpTtew6xtFRhXWlN\nbPsUkWJTjku+qGWdIo0ZSZGohyc/1Atyqka9clk6R6tlLRLBIE+Ra0VDAyLJUMs6RXFNt5qanWd4\neEgBREQ6ol6QzhyZnuPIdDYWm1KwTkm7bqhOu1/UnSUiUShJsrXw+TQL51Z1g4tIbNS7I5IMtaxT\nElc3lLqzepOlBBKRftLa462Fz6dZOLcqWKeoVTdUN2Ml6s5KhgK5yODKWje4gnUGZe1DIiIijfXr\nol5j1jnT6D7K6sqKTgtDiEgjO7es5323Prj8OG0K1ilSkE1Xu9viKZBL3jS6mG/1vLR2+trsrAoZ\nOVib2aXAu9z93Q22XQFcCSwAN7v7fdGLKGEKIMmqndReWlqVyP1tRaRY+nVOjjRmbWZ3AB8Fhhps\neyXwR8AG4GLgY2a2qpdCDppmWYiaUx2v52aOASiTvk+UsJe8+oWWanWu+91H0y5jvp/n5KgJZgeB\nq2kQrIFfAw66+3F3nwaeBs6MeJyBdfraNZnqgimq4aHGS2RqSpzkTaPAcXQmmFWii/z8axmszWyr\nmT1R93OOu3+hxdtKwFTo9xlgLI7CDjoFkOSFWyCaEicirfRzLnbLMWt33wvs7XKf0wQBu6YEHG33\npnK51O4lA2VkJOi0qK+X06qt7Sj1pTo+2dTsMYaHh9n34YsB2L77wPI49W37Dzf9H7Sjem5sajYY\ndoijflTHp/rk9Ru5dMe9y48Bliontq8YGV5+vhOq484MD0c7T3QriWzwbwO3mNlqYA3wOuDJdm+a\nnJxJoCj5detVbwJOrZfFajDptr7K5ZLqOCToMlwC4NpdE+zcsp7jC4vL248vLDb9H7Siem6sUX1H\npTpubmw0SA+anJw5pct7bHRVx/WmOu5CJfo5uRu9LIpSqf4AYGbXmdlmd/8JsBt4GPgmcKO7H+vh\nOCJ9oWEGKZIfTM6mXYSBMDw8tNy6TlLklrW7PwQ8FPr9E6HHdwF39VY0aUTzJOOxc8t6rhyfWH5c\nozHqZDSrb0nOq8ujPPPD6eXflQ0ev3ZrNcRJi6LIwFJg7i/Vd/LCF/PhFbgAXvNKjUHHqTYlrl8U\nrEVC1HMhRTI8FCSZaWgnGWOjq5cDdqrZ4CJFpsDcX6rv9KhXI37hKZ79oLtuiYgU1LrSGvqQ+zRw\nwgvQZHq5URERyb7xbRtYV9JKiEl6zStLfek1Uje4iEiBafghfmnMblCwFhER6VK/8wDUDS4iIpJx\nCtYiIiIZp25wERGRLvU7F0AtaxERkYxTsBYREck4BWsREZGMU7AWERHJOAVrERGRjIucDW5mlwLv\ncvd3N9h2B3AeMANUgEvcfbr+dSIiItJepGBdDcabgMNNXnI2sMndn41aMBEREQlE7QY/CFwNnHI/\nFzMbBs4A7jSzR8zsvT2UT0REZOC1bFmb2Vbg2rqn3+PuXzCzNzd524uB3cCu6v4nzOwxd3+ixaGG\nyuVSh0WWqFTH/aF6Tp7qOHmq42xpGazdfS+wt8t9Pg/sdvc5ADN7EHgD0CpYi4iISBNJZIMb8IiZ\nDZvZSuB84DsJHEdERGQg9LI2eKX6A4CZXQc87e73mtnngEeB48A+d/9ub8UUEREZXEOVSqX9q0RE\nRCQ1WhRFREQk4xSsRUREMk7BWkREJON6STDrWXUBlT3AmcA8cLm7P5NmmYrEzN4I3OruG83sl4B9\nwBLwJPCH7q6EhYiqMx0+C7wGWA3cDHwX1XGszGwEuBP4ZYKE1j8gOFfsQ/UcKzN7OcHMnV8nqNt9\nqI5jY2b/G5iq/vp/gY/RRR2n3bK+BFjl7huAG4DbUy5PYZjZDoKT3OrqU7uAG939AoKV596RVtkK\n4t3AZLU+fxP4FMHnV3Ucr98Cltz9fOAm4KOonmNXvfj8M2CWoE51voiRma0BcPeN1Z+tdFnHaQfr\n84D7Adz9W8D6dItTKE8D7+TEkrBnu/uB6uO/An4jlVIVxxeBD1UfDxNMU1Qdx8zdvwpcVf31tcBR\n4BzVc+zGgU8DP6r+rs9yvN4AvNjMvm5m3zSzc+myjtMO1muB8N24Fqtd49Ijd/8KsBB6KryO+8+A\nsf6WqFjcfdbdf2ZmJYLAfRMnf59UxzFx90Uz2wfcAXwefZZjZWbvIegleqD61BCq47jNAuPufjHB\nUM7n67a3reO0A+M0EF6Adtjdl9IqTMGF67UEPJdWQYrCzH4eeBD4nLvvR3WcGHd/D8HqiHcBa0Kb\nVM+9ey/wFjObAH4V+AugHNquOu7dU1QDtLv/PXAEeEVoe9s6TjtYHwTeBlDtFng83eIU2mEzu7D6\n+K3AgVYvltbM7BXAA8AOd99XfVp1HDMz22Jmf1z99QVgEXhM9Rwfd7/Q3d/s7huBvwF+H7hfdRyr\n91LNyTKzVxEE5we6qeNUs8GBewiu6A5Wf9ftNONXyy58P8FtS1cBfwd8Kb0iFcKNBN1WHzKz2tj1\nfwJ2q45j9SVgn5k9BKwkqOPvoc9ykirofBG3vcCfm1ktIL+XoHXdcR1ruVEREZGMS7sbXERERNpQ\nsBYREck4BWsREZGMU7AWERHJOAVrERGRjFOwFhERyTgFa5EBZWZjZnZP2uUQkfYUrEUG1zqC5SVF\nJOMUrEUG127gVWb25bQLIiKtKViLDK4/Av7J3f9d2gURkdYUrEUG11D7l4hIFihYi4iIZJyCtcjg\nWiD9O++JSAcUrEUG14+B/2dm30y7ICLSmm6RKSIiknE9dYGZ2RuBW919Y93z1wFbgcnqU1e5+1O9\nHEtERGRQRQ7WZrYD+D3gZw02nw1scffDUfcvIiIigV7GrJ8G3knj6R/nADea2cNmdkMPxxARERl4\nkVvW7v4VM3ttk837gU8BM8A9ZvZ2d7+v2b4qlUplaEhTPkVEZGB0FfSSmrZxh7tPA5jZfcBZQNNg\nPTQ0xOTkTEJFEYByuaQ67gPVc/JUx8lTHSevXC519frYg7WZjQGPm9nrgeeBi4C9cR9HRERkUMQR\nrCsAZnYZ8BJ3v7M6Tj0BzAPfcPf7YziOiIjIQMrKPOuKulySpW6t/lA9J091nDzVcfLK5VJXY9Za\nwUxERCTjFKxFREQyTsFaREQk4xSsRUREMk7BWkREJOMUrEVERDJOwVpERCTjFKxFREQyTsE6h7bv\nOcT2PYfSLoaIiPSJgrWIiEjGKViLiIhkXFK3yJQO1bqzx7dt6Pg9U7PzSRVHREQySC3rnLnl7sdY\nWKywsFjhlrsfS7s4IiLSB2pZp6TWolYrWURE2umpZW1mbzSziQbPbzazb5vZITO7vJdjFNnU7HzX\nreSdW9azYmSIFSND7NyyvuNjbb35AWWQi4jkVORgbWY7gDuB1XXPrwR2AW8BLgSuNLOX91LIPGk0\nrSrqVKtm7xsbXc3Y6OoG7xARkSLqpWX9NPBOoP4G2q8Dnnb3KXc/DjwCXNDDcTKp17nOte7vdq3k\nozNzahGLiAy4yMHa3b8CLDTYtBaYCv0+A4xFPU4RTM3OnzQ2HU4SA/rSSn52el7j4yIiOZVEgtkU\nUAr9XgKOtntTuVxq95JMGRkJOhTqyz01e+yk57fvPrAclG/bf5jxay5g5YqR0DuGGBkZavj3T80e\nY6nCKdv3ffjipuXaevMDAOy9adPyc0EZlk4qQ6ca7U9ay9tnOY9Ux8lTHWdLEsH6e8AZZrYOmCXo\nAh9v96bJyZkEipKcZ6fngJPLHbSYg6B47a4Jdm5Zz/GFxeXtxxcWmZycYcdlZ3HleJCXNza6isXF\nyil/f3hf/3z0hY7rZ7F6YRB+faMydKrR/qS5crmkukqY6jh5quPkdXsxFMc86wqAmV1mZldUx6mv\nB74OHAL2uvuPYjhOZnQz17lZ9nYtSWx824aOFkTpZU51UIbhrjPI4dQufBER6b+eWtbu/o/Ahurj\n/aHnvwZ8raeS5dDOLeuXW8z1gTnKvt5364Ndv69ZYD1t7erlVnKnahcltcfdBnoREYmHVjCLoNVc\n506nVXXSoj597RqAjlvEWt1MRKSYtIJZRP2a5zw8FM+x9t60qeMxqPB65Y16CpIQZY10EZFBoWAd\ns0bBJmoAGt+2oas51s264XuhxVdERNKnYB1Rv1qA3R4njuCaRkKZkthERJrTmHVIr6uSFUEa497d\nHlP/JxEZNGpZF0zcLf5+jyF//8ea2ykiUk8t6xDNKY5+V684jgmddePr/yQig0bBukrTnk5I465e\nteO1C8L6P4nIIBq4YK3xzmyqBWkFYRGRUw1csG4mje5f6V64y/y5mWMpl0ZEpD8GKlhv33OIozNz\nTben0f2bRZ2uVx6nz2zf2PHF0tjoaobr76IuIlJghcgG72b1q6WK5vRmlS6UREQaG6iWtcZFsyGO\nvAFddInIIClEy7pTC23uOqV1qfPhpaVVHJmeW77oqo1d6/8nIkVViGDdSQurviWtJLLsUbAVEWks\nUje4mQ2b2X81s0NmNmFmv1i3/Toze7K6bcLMfjme4p4qyrzbWjaxpKPXRU3qM/e1SIqIFF3UMetL\ngFXuvgG4Abi9bvvZwBZ331j9eaqXQsah21WyJBm9LmpSG++uZe5rkRQRGQRRg/V5wP0A7v4toL5P\n+RzgRjN72Mxu6KF8bXUzP3psdDWnr12j7tYC06I3IlJEUYP1WmA69PuimYX3tR+4CrgION/M3h7x\nOB2ptbJ0os6+Xhefqe/yjmMxG31uRCTroiaYTQOl0O/D7r4U+v0Od58GMLP7gLOA+1rtsFwutdrc\n0kjdGHSzfe378MWRj1EEvdRxnE5buwbovjzbdx9Yzugvr3sR49dccMr+ap+FbvYd5T2tZKWei0x1\nnDzVcbZEDdYHgc3AF83sXODx2gYzGwMeN7PXA88TtK73ttvh5GT0WyMu1k3J6mVfRVUulzJTL7X/\nV7flOb6weNLj2vvD+3t2eq7jfdda07WWeu09tZXu1pW6HzLJUj0Xleo4earj5HV7MRS1G/weYM7M\nDhIkl11nZpeZ2RXuPkWQdDYBHACedPf7Ix6nI7XlMZUVnA9RlzNt1+UdJdlsanZeCWoiknmRWtbu\nXgGurnv6qdD2/QTj1rFrtrRo7URde6x51MXUKJO/9llQsBWRohqo5Ual2LpNNqv1wtS/Z2p2XsuZ\nikimFCZY6xaXAp3fOS3cZV57X/h50BryIpIduVtutFVrR4udFF8e5sh3cxe4fslimUSkc7lqWWu1\nKmmn0+S1cE9M+CIvvNJdv3tpNN9bRJrJXcu6FbUapBu1IF3/uRkbXc3RmblC9dRo/F0k33LVsta4\ntORB1CmESU09VI+USP7lrmVdpNaOZFMvPTThVdYaTSHU1EMRiSJ3wVokLoMybLJzy3quHJ9Yfiwi\n+ZPJbvBWiTZRV78S6Yfxay6INFST9BBPp1PaRCSb1LIWiVmroKiphyISRSZb1iJFFGeil6Z5iQyW\nTLasNc17HlVhAAAQyklEQVRE8iyLwzRZLJOIdC5zLWtNM5Giajcu3U0+hu4wJzJYMhesRYqs00Sv\nVt3ceb+gVRe+SPcyF6y18IlId8LBL4uBMItlEsmbSGPWZjYM7AHOBOaBy939mdD2zcAHgQXgs+5+\nVyf7rX2hlRUrRdVNN3cz9fOm0w6E3d4kRN33It2L2rK+BFjl7huAG4DbaxvMbCWwC3gLcCFwpZm9\nvNeCigyKVt3ctVZqs+70pMayO20dt3td3rvwReIQ5QI7ajb4ecD9AO7+LTML91e/Dnja3acAzOwR\n4ALgS53uXJmrIp2rBefwkqVXjk/wme0bEz1u7YTT7uJALWmR3kVtWa8FpkO/L1a7xmvbpkLbZoCx\nTnaqDFeR7vI2wi3VH0zO9qmEJ0zNzp/UUq7/Dte3pJWTIhJN1Jb1NFAK/T7s7kvVx1N120rA0XY7\nvG3/4eVWwW37DzN+zQURiybNlMul9i+SnsVRz6etXdNwX1OzxwC457bNQHDjkJpf+Lm1PP2DKRYW\nl4ChWP/fteOG9zk1e2z5OwvwTz+dPeU7vHLFyPL2lStGKJdLTf+2buiznDzVcXKmZo+x+f1fPXTv\n7e/ouBs5arA+CGwGvmhm5wKPh7Z9DzjDzNYBswRd4OPtdnh8YfGkx5OTMxGLJo2UyyXVaR/EVc+L\n1aAX3lfQSg2uia/dNcHOLevZcdlZy8lm4ccLi0vLr+nVleMTy0G4ts9wWQBWjAzxqpeN8swPgw63\n2ne4vnyTkzMN/7Zu6LOcPNVxckLfnTd1876o3eD3AHNmdpAguew6M7vMzK5w9+PA9cDXgUPAXnf/\nUbsdqntMJJos3KRjaanS9DuchfKJ9FMS0xUjtazdvQJcXff0U6HtXwO+1u1+9YUWCTRKsuzkVpdj\no6s5Mj2XyEVveJ/hsiwtnegK7+Q7HDWBtHby2/fhi7t+j5JWJStq352Fxcqj3bwvk2uDi0hj7YLh\n+LYNsV7RhzPMm5XlyPRc2zKJyMlriez78MVdfTEyFaz1pRbpXprfm/rM73pxtmw1U0TyIonPauaW\nGxWR5rq52UccwuPQnczbTqp84Slg4Qz4djQdVPqtk4V/onxHMtWyFpHexR0sm3W9147Tbhy9mfpW\nd9zjy+Es9tocb5E0bN9ziKMzcyxVgtyPKNSyFpGetMr2jqtlG27hh9dguPzjD3L5xx/sef8icWk2\nK6KWhxl1qV0FaxFJRLvuwG4DeauLglZTZdplxuuuYBK3JKYrqhtcRFpKcgy69njnlvWRWuBLjRPV\nW2axi/Rb+LMddVqlgrWI9KRZMO9kXnhNo+DdTri1PjU7H7klowQ0iVv4OxHXhaO6wUUkMc26A6Os\nWFifaR6+ccnCYuWkoFu/gloztSQ03bJTsk7BWkRSEQ7knQbv7XsOsfXmBwB4dXn0pG31AXd4KPjR\n+g0St27yHGqfbaDjKZCNKFiLSCZ0m5TTrkW+VGk+pl1P9ySQJI2NrmZ4qLcltTVmLSKJSbpVO1yd\nsjpcfVALuOEWdrMx8HZjiVpXXJpJI89BwVpEUlEfBDsJisFJ8tRFJXQTIOmXKMmQcVzwqRtcRHLh\nxLztpeWW810fuIi7PnDRKa/tZAw8PJaoYD848jqvXi1rEcm9Ri2XTgJw7TWN3q8pXdJIN1MS49R1\nsDazFwH/DSgDM8B/cPef1r3mDuC86vYKcIm7T/deXBEZVCdOkskng0Xp6uyWxsTTEcdFWBo9MVFa\n1lcDf+vuHzGz3wFuAq6te83ZwCZ3f7bXAoqI1IyNrmYk4o0QOlELoC8trUrsGJKeRhdhebloihKs\nzwM+Xn18P/DB8EYzGwbOAO40s1cAe939z3sqpYhIlzo5+cax+lozeQkCkg8tg7WZbeXUVvNPgFqX\n9gwwVrf9xcBuYFd1/xNm9pi7P9HqWOVyqdMyS0Sq4/5QPSdn34cvTnT/tVZ7uVzitLVrlh9HcfRn\n8y3fPzV7rKf9Jy2r5erFJ6/fyKU77l1+DNH+D0l/DhtpGazdfS+wN/ycmX0ZqP1VJeC5urc9D+x2\n97nq6x8E3gC0DNaTkzOdl1q6Vi6XVMd9oHpOXpJ1/Oz0HBCcjxar3aWRj1Vp/v6gO3YJgGt3TWRu\nQZYif47HRoMhjsnJmYb/h371iHR7MRRl6tZB4G3Vx28FDtRtN+ARMxs2s5XA+cB3IhxHRKRv2t3S\ns1vdrKAmycnrVK16UYL1p4FfMbOHgcuBPwEws+vMbLO7fxf4HPAoMAHsqz4nIpIb9TcO6Ub9Cmr1\notzIROKXp/9D1wlm7v4C8NsNnv9E6PEugjFrEZFciHP+bPiOYOHHYZ1M/1GSWvzq67L+/5DV+fVa\nwUxEpKrbm4k0E74jWP3dwaS/pmbnOw7AcQ+FxEnBWkQkZuGWeda7V4usUfDN6xi2lhsVEamKs7t5\nuM3aLZ3fuESSFP4/pLWUaCfUshYRScC60hrWldZEfn+Wu2TzolECWbtu8biGQuKmlrWISALqW85x\nJIsp4awz4XoKB95+rPmeFAVrEZEMynKXbJFl9UJIwVpEJKOy2B2bB+Fu7ryMSbejMWsRkT7oZgpR\np/vIa2ZzktqN9Wd1TLodBWsRkYRFTRYLr6IWR8KZgnt+KViLiORUHK31osnTEqLdULAWEUlYHAGk\nfh9RWtpHpuc4Ur2zWFjRWtyturp7WfM9TQrWIiJ9EMdYaS/7aHdzkSwp2sVDHJQNLiKSQ3FmNvej\nK70++OaxdZsmBWsRkT6IIzi1u2NUKzu3rOd9tz64/Lim2UIhaS7AMjU7z8Jihe17DkU6fhEvBCJ3\ng5vZpWb2+SbbrjCzvzazR83s7dGLJyIicUi7G7zTZLjwxYOS506IFKzN7A7go8ApS9Wb2SuBPwI2\nABcDHzOzVb0UUkRE4vP9H88st5zjyp5uNc4cToY7Mj2nIBxB1Jb1QeBqGgRr4NeAg+5+3N2ngaeB\nMyMeR0REmugmszkciOu7zxslrk3NznNkei6RRK+FxcryeHujcq4YaXPLsgHUMlib2VYze6Lu5xx3\n/0KLt5WAqdDvM8BYHIUVEZHoTl+7htPXtr8T2PbdB7ruim7Vzd0oANf234ruOHZCywQzd98L7O1y\nn9MEAbumBBxt96ZyudTuJdIj1XF/qJ6TpzqOZqQaMJ+dDoJqrR5rzzeq14XFStv6Dgf32/YfZvya\nC055zVIFhoeHGB4aYmFxafl4W29+AIC9N20KvXoICPa3csWI/t8kkw3+beAWM1sNrAFeBzzZ7k2T\nkzMJFEVqyuWS6rgPVM/JUx1Ht7hYWc60Brh21wQ7t6zn1qveBJw4D49fcwGb3//V5ffVXtfM8YXF\nkx43/P9Uqi3pUAP72l0TLFbLEn7P2OgqjkzPsWJkiB2XnVXI/3e3FyC9LIpSoXbpA5jZdWa22d1/\nAuwGHga+Cdzo7sd6OI6IiMRgfNsGXvPKzoJEN+PGnSSpLVWCn06Mb9vA6WvX5PKGG0mJ3LJ294eA\nh0K/fyL0+C7grt6KJiIicet0MZWx0dXLrdtOssRbBdbwuPPIyDALiyda4soM74yWGxURGTCdLFva\nTet2+55DHJ05dc3xdn4wOdt0ffO8ruGdFAVrERFJVLib/NPXX7j8/KvLoymWKl8UrEVEpGfrSmta\ntoTDrfnhoeCnqLezTILWBhcRGTCddi93+rpux53XlU7M9VYSWWcUrEVEJLJmNwJpRWPR3VOwFhGR\nxDUL0ArcndGYtYiIRKZx5/5Qy1pERHqicefkqWUtIiKScWpZi4hITzTunDy1rEVERDJOwVpERCTj\nFKxFREQyTsFaREQk4xSsRUREMi5yNriZXQq8y93f3WDbHcB5wAxQAS5x9+nIpRQRERlgkYJ1NRhv\nAg43ecnZwCZ3fzZqwURERCQQtRv8IHA1MFS/wcyGgTOAO83sETN7bw/lExERGXgtW9ZmthW4tu7p\n97j7F8zszU3e9mJgN7Cruv8JM3vM3Z/otbAiIiKDqGWwdve9wN4u9/k8sNvd5wDM7EHgDUCrYD1U\nLpe6PIx0S3XcH6rn5KmOk6c6zpYkssENeMTMhs1sJXA+8J0EjiMiIjIQelkbvFL9AcDMrgOedvd7\nzexzwKPAcWCfu3+3t2KKiIgMrqFKpdL+VSIiIpIaLYoiIiKScQrWIiIiGadgLSIiknEK1iIiIhnX\nSzZ4z6qrne0BzgTmgcvd/Zk0y1Q01elznwVeA6wGbnb3e9MtVTGZ2csJpin+urs/lXZ5isjM/hjY\nDKwE/tTd/yLlIhVK9Zx8F/DLwBJwhbt7uqUqDjN7I3Cru280s18C9hHU85PAH7p704zvtFvWlwCr\n3H0DcANwe8rlKaJ3A5PufgHwm8CfplyeQqpeFP0ZMJt2WYqqumrim6rnizcD/yrVAhXTJmDU3c8H\nPgLcknJ5CsPMdgB3EjSaIFjl88bquXkIeEer96cdrM8D7gdw928B69MtTiF9EfhQ9fEwsJBiWYps\nHPg08KO0C1Jgm4AnzOwvgXuB/5FyeYroBWDMzIaAMeBYyuUpkqeBd3Linhpnu/uB6uO/An6j1ZvT\nDtZrgfCtMxer3TASE3efdfefmVmJIHDvTLtMRWNm7yHovXig+tQpN7iRWJSBc4B3AX8AfD7d4hTS\nQWAN8D2CnqL/km5xisPdv8LJjaXweeJnBBdHTaUdGKeB8AK0w+6+lFZhisrMfh54EPicu//3tMtT\nQO8F3mJmE8CvAn9hZq9IuUxF9FPgAXdfqOYEzJnZy9IuVMHsAA66u3His7wq5TIVVTjWlYDnWr04\n7WB9EHgbgJmdCzyebnGKpxo0HgB2uPu+lItTSO5+obu/2d03An8D/L67/yTtchXQIwR5F5jZq4BR\n4EiqJSqeUU70dh4lSOQbSa84hXbYzC6sPn4rcKDVi1PNBgfuIWiRHKz+rntfx+9Ggu6VD5lZbez6\nrbW7oonkhbvfZ2YXmNm3CRoa21plz0ok48Cfm9nDBIH6j939hZTLVDS1z+z7gTurPRd/B3yp1Zu0\nNriIiEjGpd0NLiIiIm0oWIuIiGScgrWIiEjGKViLiIhknIK1iIhIxilYi4iIZJyCtUhBmdlrzewf\n0i6HiPROwVpERCTjtCiKSAFUbx95U/XXVwPfBm4mWMJwAvjXBMtHXuLuz5rZfwR+j2B5ySXgd9z9\ne2b2nwnu/rMIfNXdP2JmLwE+BfwKwdKTH9ca8yL9pZa1SHGcC1wFvI7gzkm/RXCnqtvd/d8APwH+\nffUObO8ALqw+/5fANjP7l8BvuvuvAhuAXzKz1QQXAY+5+3rgQmCnmf1Cn/82kYGmYC1SHN9w92eq\n62XfDVwE/NDdH6tu/z/Ay9x9Bvhd4HfN7GPAZoIW9g+BF8zsEeA64IPuPk/Q0v4DMzsMPAS8GHh9\nP/8wkUGnYC1SHOF75Y4Ax+ueqwBD1Vum/i+C+8nfB+wjuD3tIvBG4IPA6cCjZnYGwXni3e5+lruf\nBZxHcCc3EekTBWuR4thoZq80s2FgC/BXnHyD+5r1wN+7+x3AXxPcpnbEzM4kaDkfcPftBHcCMoJ7\noW8DMLOfAw4TjIuLSJ8oWIsUxw+BzxN0d/8A+AYnbsdXUyFoFQ+b2ZPA1wkC9Gvd/XHgUeBJM/sO\n8A/A/wT+BHiRmT0BfJPg3uiaEibSR8oGFymAajb4B9z9rWmXRUTip5a1SDFUOLUVLSIFoZa1iIhI\nxqllLSIiknEK1iIiIhmnYC0iIpJxCtYiIiIZp2AtIiKScf8fsKZTxTyGD8UAAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 4 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let's take a look at how supersmoother handles this. Using the standard algorithm to smooth each plot, we can do the following:" ] }, { "cell_type": "code", "collapsed": false, "input": [ "from supersmoother import SuperSmoother\n", "\n", "phase_fit = np.linspace(0, period, 1000)\n", "t_fit = np.linspace(0, 50, 1000)\n", "\n", "model = SuperSmoother()\n", "y_t = model.fit(t, y, dy).predict(t_fit)\n", "y_phase = model.fit(phase, y, dy).predict(phase_fit)\n", "\n", "ax = plot_data()\n", "ax[0].plot(t_fit, y_t, '-k', lw=1)\n", "ax[1].plot(phase_fit, y_phase, '-k', lw=1);" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAesAAAFkCAYAAAAAFROsAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XV4FMcbwPHvxbAIdsVdhmJFggXX4FKkUAgFUrx4kWIt\nFCqkOCS/QnEoFC9OcJdSaLF2gUIFCjRYEkIJkfv9EWkSope73CV5P8+T57nc7u28DJd9d2ZnZ3QG\ngwEhhBBCWC8bSwcghBBCiMRJshZCCCGsnCRrIYQQwspJshZCCCGsnCRrIYQQwspJshZCCCGsnF1q\nPqyUqgV8oWla4zjvjwI8Ab/ItwZqmnYjNWUJIYQQmZXRyVopNQ7oBTyPZ3M1wEPTtEvGHl8IIYQQ\nEVLTDX4LeBvQxbOtOjBRKXVCKTUhFWUIIYQQmZ7RyVrTtK1AaAKb1wMDgSZAPaVUG2PLEUIIITK7\nVN2zTsR8TdMCAJRSu4GqwO6EdjYYDAadLr4GuhBCCJEhpSjpmTxZK6VcgMtKqfLACyJa18sS+4xO\np8PPL9DUoYgY9HonqeM0IPVsflLH5id1bH56vVOK9jdFsjYAKKV6AI6api2NvE99BAgGDmqats8E\n5QghhBCZks5KVt0yyFWcecmVctqQejY/qWPzkzo2P73eKUXd4DIpihBCCGHlJFkLIYQQVk6StRBC\nCGHlJFkLIYQQVk6StRBCCGHlJFkLIYQQVk6StRBCCGHl0nWyHut9mrHepy0dhhBCCGFW6TpZCyGE\nEJmBuRbySBP+QcGWDkEIITK9qB5OryFuFo4k40q3LeuZay4QGmYgNMzAzDUXLB2OEEIIYTbpumUt\nRHomrRGRUUgvp/mlqmWtlKqllDoSz/vtlFLnlVKnlVLvp6aMmGIOKJvk4YqdrQ47Wx2TPFxNVYQQ\nFicDJ0V6Ir2cacPoZK2UGgcsBbLEed8emAM0BxoCA5RSb6QmyIS45MiCS44sSe8ohBXyDwrGPyhY\nkrMZSJ2KjCY1LetbwNtA3GW+3gRuaZrmr2laCHASaJCKcqL/8KJObkKkdzFbI3G/0/I9F+lJzF7O\nZ4Gv5CLJTIy+Z61p2lalVPF4NjkD/jF+DwRckjqeXu+U4Db/oFeEhoVH/z5r/SW8hjdg5cfu8e7v\nOcMXgGWTWyRVbKaSWB0L00lOPdvb2Ua/jkjYr9DrnRi74DihYRFrzEd9z8Xrkqpj/6BXCe4n54fk\nScn5IrdzVgCeBASn+LMiecwxwMwfiPk/5QQ8TepDCS10HtECCY/1XkhoWKILo4dFnuxk8fT/yGLy\naSO59TyuR1UGeB2JTsyhYeGMnBN7+EdS3/PMKqk6jnnOGDnnyGtjWuT8kLSUni/CInuIor7P8dW7\niC2lFzTmeHTrV6CMUiqXUsqBiC7wM6Y6eHIGlEk3okgPXHJkwc429l0kGThpfnJ+MD2vIW4Uyy+t\naXMyRbI2ACileiil+kfepx4N7AdOA8s0Tbtv7MFjnrzyOGdNckCZjEwU6UlUwo6ZnGXgZOokdsEj\n5wfzkQtN80pVN7imab8DbpGv18d4fxewK7nHCQ0NTXR7Sk5cfzwI5NW/AYS8DIJCztHvyzOtwtpE\nfRdlQI7pycWOZUi9m49VTIqi1+tp0aIVffp44upaM8H9Zg2uEyux//PPP/j67qVEiZK4utYkS5aI\nFsnPhzbz4OYplky8lBbhC2FSckFpPpM8XBngdST6tRDphVUk62vXrrF06Up69+7O1au30Ol0BAYG\n8ODBAx48uE/13PfYu3c35Rf2pHjxknh7L8XbeyHbt2+hceMmrF27Ek3TqFmzFv+EF+bBzTO45C3C\nRx+Nxdt7KSAz7AjrJcnZ9BKrU2n9mY98l83HKpJ1wYIFGTJkGBs2rKVGjco8fvwIGxtb8ufPT/78\nBciXLz9NmzZnxowvaNGiIa1aNeG99/px+vSP6PV6AJ49e8qoz1by/NJZilRsRtmaHbm4cyI7d27n\n8rPC0aMUZ665IFfUQghhBHPcTpRblMljFck6yoYNW/n33xfky5cfR8f4RxZu376XQoUKvbY9Z85c\nlK5UH13utwAoVciZIYu+pnfvHnQc6g3Ymzt8gfzhCesn302RHlnVqlsFCxaiVKkyCSZqAKXKJbg9\n7mhEV9ea9O7dhzsnl2Brk7zHvkTqyGMxQmRc5vj7lnNG8lhVsjaFuI+9jBkzAT8/Px7+ekjuVZmZ\nPBYjRMZljr9vOWckX4ZL1nHZ29uzePESLh1aRcDjvy0djlWRxQ6sw/Pngdy5c9vSYQghrFiGS9Ze\nQ9xeuyelVDnGjBqJ7f39kqDMSCZFMM6WLZuoU6caI0YM4d69u8yb9xXXrl21dFjphvxNpw1z/H3L\nOSP5rGqAmTl5ePSlbt3qNMjXDBd9EUuHYxXMcZ9IbjWk3LNnT/Hw6EvOnDlp0qQuoaFhLFniw+rV\n6xOddyCzkwSd9szx9y3njOTJcC3rhOTJk4exYz/i+LbZPHoWlOn/0OVekfXw9/enSJEiTJr0MUeP\nnmHjxm3Mn78YD493OHHiWKx9pRX5OhmgJDKDTNOyBrhvVxV0On7/aS/2NdpZOpwMSR6LSTl/f3+K\nFi0GQIECBSlQoCAAy5at4f33ezN37mLc3VsB8Pu1Ezjlyk/kLL+ZWsxVniD5cyjcu3cXvf4NHBwc\nzBlehmSOv285ZySPUS1rpZSNUup/SqnTSqkjSqlScbaPUkpdjdx2RClV1jThpo7OxobKzT/gxtnv\nCHjywNLhWJTcK7IeAQH+5MyZ87X33dzqsW7dJkaPHsaWLRsBuHp6G4fWT+PJk8dpHaZVieoZSszZ\ns6d5/vz5a++/+25X2rRpzq1bN80VnhAmZ2w3eEfAQdM0N2ACMDvO9mqAh6ZpjSN/bqQmSFOZ5OGK\nY+5ClHLtxOUDi5mx+gdLh2RRsrqTZfn67uX27d/w93+Gs7NLvPtUrVqdzZt3MH36VLoP+Zjnzx6Q\ns4CiZacehIeHx/uZzCjuRaefnx+dOrWhSpU3GTSoHwcO7CMkJISwsDDu3PmNLl260a5dC1auXIbB\nkHDSl9sOwloYm6zrAvsANE07B8RtmlUHJiqlTiilJqQiPpOzs9VR0rUDoa9ecvHod5YOx6LiGzkv\n0sbz54EMGzaI3bt3EhDgj4tL/Mka4M03y7N9+x4uHN5AyMtAqrQcwcsXASxevCANI7YuSS2d++ef\nv1OpUmXOnr1EzZp1mDdvNpUrl+WDDwaSJ09eBg4cyo4d+1m7dhW9e3fHz8/PQv+SjOXFixeEhYUZ\n9Vm5MEqcscnaGQiI8XuYUirmsdYDA4EmQD2lVBsjyzE5lxxZsLGxpWaHcWjntnPmzClLhyQyoRUr\nlhEcHMzZs6e4c+c2hQsn/oRCiRIlOXHkCBUavkeWLFnYs20TPj4LOXv2TPQ+me1kl1jP0F9//Unh\nwkXJmzcv/fr1Z/fuA+zbdwSlytGnjycAZcqUZc+eg5QtW44mTepy8OD+144jg9dS9r0aMWIIXbq0\n5+nTJ2aOKvMxdoBZABBzzk8bTdNi9snN1zQtAEAptRuoCuxO7IB6fcJTjJrSyo/d8ZzhC7mKMGr1\nKjw9Pfnxxx/Jnz9/mpRvSWlVx5ldUvX84sULlixZzMcff8z48eOZMWMGFSuWScZxFTUadwWgatXy\nrFy5gkGDPLl48SJ6vR5bW12yys8I9Hqn6H/vssktADAYDNy6dYvDhw+zatUqGjRoEKsu9PpKuLpW\neu1YCxbMoXPnDvTu3Zt27Y7i5eVFtmzZGLvgePR98VnrL+E1vEEa/MusR1TdJed79fLlS4KDg7l+\n/QpVq1alTZtm7Nq1C6VUssvzD3qVZDmZmbHJ+hTQDtiklKoNXI7aoJRyAS4rpcoDL4hoXS9L6oB+\nfoFGhpJyYZF/gNWr16VHDw8q1mxOnS7TWP5R8zSLIa3p9U5pWscxZabFPZJTz0uX+lC9ek1at36b\nCxcu0a/fkGT/30R9d/38AqlZswEdO3bhnXd6sH79Fp4EvIzelpFF1fEXA+sAEf/eadOmsH37FsLC\nwqhfvyE9evSmbdv2ya6L8uWrcejQScaPH02VKlXx8VlGSOh/3bkhoWEZvl5jivk9Ts73asGCOcyd\n+xWhoSEcPXqWTZs2UK9efXx8vqFhw8ZJlhcxYDCivTdyzpFMMeg1pRclusQGVyREKaUDvIHKkW/1\nJeI+taOmaUuVUj2AUUAwcFDTtGlJHNJgqT+EsLAwKri6UbRiM3zXzLRIDGlBknXaSKqeg4ODqVnz\nLVavXs9bb1VNdXkhISF06tQGu7wVyV2uLRCx4lxGPtnFreOQkBCKFn2DEyfOU6pUaXQ6XaqOv3nz\nd0yd+hFDhozgkKbDIZsT383qm9qw05WoOp655gK/3Yu441mqkDPDO5Zly5ZNFC5cmBIlSlKkSDEc\nHBzo378PFSpUpEiRonTu3A2A06dP0r9/Hz78cAJ9+76faHlxy8nI398oer1Tir6oRrWsNU0zAIPj\nvH0jxvb1RNy3ThFLnNRtbW0pWb0Dv/+0J83KzGwy+z2/mDZsWEeFChVNkqghYu77JUtW4Fbfjbcc\nS5CncAWTHDc9efLkCbly5aJ06aRvJSRHly7vULNmbYYOHcD58+dwyOaEb7P8tGjRyiTHT89OnDjO\nwoVzKVWqNHfu3Ob+/b/Jn78Ajx49YtSosZQv/9/3z82tHrt2+dKrVzdu3tSYPv1z7OziTzmTPFwZ\n4HUk+rV4nVXNYGaJwRyD5xwjT+GKPL1/g/5f7EvTsjMDmSntPyEhISxcOJdRo8aa9LgFCxZi+dJv\nuLhrFk/++jnTneweP35Enjx5TXrMokWLsX37Htr0n0/THp8wZswIvvnmfyYtIz2IOx/D77/foVWr\nNmzcuJ0ffrjMnTv32bhxOxs2bOHNN8u/9vkSJUqyZ89Bbt26Sc+eXQkI8E+wLHmUNHFWk6wteVK3\nz+pIQVWP8ztnJ/nYQWYbcWsOmbUOt2zZSLFixalRo5bJj92kSXPqdRzNxd1z+PPPP0x+fGtmjmQN\nEb1uy2e8x8rPPdm1y5eVK5cxceJYox9NSq9iJtE7d25TvHiJ6G329vaULFmK2rXdErz94OKSk2+/\n3UyJEiVp3boZv/9+J03izmisJllbis/ohgBUajqQkm84MG3aFAtHlLHITGkRwsLCmDfvK0aPHme2\nMgqXrUXFet0YOLAfISEhZivH2jx65GeWZB1TsWLF2b37AJqm8d57PeKdGS2ju3//b/bv30O1ain/\nO7azs+OLL2bTt29/2rRpztmzr1+sJ2feh8x6oQ9WlKwteVK30YGdnT3Ll6/h4MH9rFyZ8OB1ee4y\n5eJ2b2XGOvz++63kzavHza2e2crwGuLGrtVe5MqVi88//9Rs5VibvXt3Ua7cm2Yvx8UlJxs2bEGv\nf4MOHVpx//7fZi/TGngNcWN636q8914PPD0HGJWso3h6DmDRoq/p168XGzasM2GUGZ/VJGuw/D2L\nXLlys27dJry8Pufw4YOvbZf7r8Z79TIIyJx1GB4eHtmqHpvqkcpJsbGxYeHCr9m6dROHDvmatSxr\nsG3bZq5evcLQoSPSpDx7e3vmzFlIhw5v06pVU65cuZz0h9I5g8HAqFFDKVmyFMOHj0718Ro3bsr2\n7XuZM2cW06dPTdG0uZnxQj+KVSVrS/lmfBO+Gd8EiBgQsXz5Wj74YADXr1+zcGTpS3BwMHfv/vXa\n+18MrIXv0gEMHvw+wf9mvu7DvXt3kzVrVho3bpYm5eXJkwcfn28YMWIoDx7cT5MyLeHJk8dMmjSe\nRYu+Jnv27GlWrk6nY/jwUUyf/hndunXA13dvhu6e3b17J5qmMXfuYpNdbJYtq9i37zAXL16gT5+e\n3Lx5I8mxAJnxQj8mq0rW1jJXda1atZkx40t69erGw4f/rc4l919fFxQUxM6d3zNoUD8qVixD8+YN\noq+UDQYDd+/+xaZNG3B2dsHR0Yk9S4YT8EDDhrBMUYcGg4E5c2YxatQ4s7eqY6pTpy59+ngyePD7\nGXZA1KlTJ6latVqqumVTo337TqxZ8x2jRw9n45zeaFcyztTFnjN8oy8+bt7UaNKkGdmyZTNpGblz\n52Hjxu3kyZOHunVdGTzYM1ONtUgpq0rW1uTtt7vSq9d7dOzYmp9/vhT9vqW76q1BYGAAW7ZspG/f\nXlSurFi1ajm1a9fl5MkfyJ07D3379qJFi4aULFmIVq2asmnTBiZMmIyX11ymT/+Mc9tmsGNOZ9at\nW23pf4rZHTrkS2hoaPR61Glp1KiIbvc5c2aledlp4fLln0z2vLqxXF1r0vr9OeiLVeHinvkMnb7S\novGYw/37f1OgQAGzHNvBwYG5cxdx69ZfBAYGMnz4YEJDQ+PdN7M3loydbjRTGDVqLPny5cfDozu7\ndx+gSJGilg7JYp48ecz+/XvZtet7zpw5TZ06brRt24HZs+eTO3ee6P1mzpzF33/fo2xZRZkyZXFx\nib1Oc+vWbWk32If7ty/x2WfTKVCgAE2aNM+Qs5wZDAbmzv2KUaM+xMYm7a+LbW1t8fH5hqZN61O3\nbn2zDm6zhMOHD/LJJzMsHQYueQpQufkQCpSpw66VHzO0S41Yk4OkR08CgoGI2S0fPLhPw4ZNzFqe\ns7MLy5at4b33etC3b0/+97/l5MiR47X9MnNDyajpRs3AYtONJseSJd6sXbuKXbt8E1x32NoZM93o\nw4cP2bNnJ7t27eCnny7SoEEj2rZtT/Pm7iaph3PnztKnTw82btzO6hMRA9DSe7KOWc+nTp1gzJjh\nnDp1AVtbW4vFdPjwAUaNGsahQyfJm9e8jzilBb3eiePHz9G9+9tcvHjNonUbJWr2rValH/PJJ5PZ\nsWMfxYoVt2xQRoo79eeBZcOZM2dBmtxuePXqFaNHD+PmTY21azeh1+vNXqalpHS6UekGT8JY79Pc\nCHkLN7d6vP/+exn+nsrdu3/x9deLadfOnbp1XTl37jR9+nhy+bLGihVr6dy5m8kuWGrVqs2sWXPp\n1esd7t+/m+FGec6d+xXDh4+2eDJp0qQ5nTt3Y9iwgSkaeWvNNm3aQJcu71i8bqNE3R7r1KkLw4eP\nplu3jvzzzz8J7p9eBqQ9+PNX/P2fUanSW2lSnoODAwsX/o8mTZrTunVTfvvtZpqUmx5Isk4GnU7H\njBlfYmtry0cfjcVKeiNM5vbtWyxYMBd390Y0a1af69evMXz4KK5du8X//recdu06xNslZQrt2nWk\ndI0OnPxuCn5//ZphRnlevHiB27dv0aXLO5YOBYCPPprCs2fP8PFZZOlQUi0sLIwtWzbStWt3S4cS\nL0/PAXTu3I0ePToTGBhg6XBSLOLesA12tjqCftvPoEFDsbe3T7PydTod48dPYuTID2nfvhXnz59L\ncN/0ctFjCkbds1ZK2fDfqlvBwPuapv0WY3s7YAoQCizXNO0bE8RqEVGtPTs7O5YuXUnbtu74+Cxi\nyJBhQPpdUerp0yds2bKRb79dy8OHD2jduh0TJ36Mm1u9NP3DBKjWsCt+T55zae9cqlRZmeB+6amu\n582bzdChw3FwcLB0KMB/C364uzfih3uO6AuXSxf1GNdY79M8/P0n8ubVp8lEKMkVty7Hjv2Ix48f\n0bt3xPKlWbNmjbXd2nuRcjtn4cnDvzhz5iQLFvhYJIaePXtToEAB+vTpwaxZ82jbtr1F4rAWxras\nOwIOmqa5AROA2VEblFL2wBygOdAQGKCUeiO1gVpC3Of6HB2dWLduI19/vZjdu3daOrxkG+t9Gs8Z\nvoSHh3PixDEGDepHjRpv8cMP5/j440+5fFnDy2suDRs2TvNEDRFX8qVqvE2uAooft0/n2bOnsWKP\nStLpZUKE69evcfHiBd59t7elQ4mlSJGizJo1j+ObP+dVOn7e/dZPB+nWzTpb1VF0Oh2ffeaFXq9n\n0CDPWCOc08vzwtdOb6FPn/dxdHS0WAxNmjTnu++2MWnSOJYs8X5te0LnhIzY4jY2WdcF9gFomnYO\niDny4E3glqZp/pqmhQAngQapitKKFCpUmNWr1/Phh8O5cuXndJFADAYDV05tplatKkyePIHq1Wtw\n/vxPfP31Cho2bGwV9/1sbXRUaz2CKlWq0a6d+2uTq6SXExzAggWzGTBgiMmfSzWFtm3bk6twJc4d\nXGPpUIziHxTMXzd+oE0b629l2drasmjREoKCnjN27EgMBgNjvU/zxwPrHUwbpVOVV/zz2xk8PQda\nOhQqVXqLXbt8Wb16BVOmfBQ97iI9nRNMwdhHt5yBmDdjwpRSNpqmhUdui7kOWiCQ5Igkvd7JyFDM\nZ97oxnQatzP6dZRmzRrg7e1Np65dqf3OLByyOTNr/SW8hlvnNcnf9/7g8vGNnDh6AFdX1zSdnCO5\nwg0ANvj4LGTu3Lm0b+/O7t278Q96BYA+13+Jz97O1iq/LwC3bt3i+PGjrFixDGdn64tx7ILjKLd3\nObFuDAMne7P16/GWDinZxi44zr8vggh99ZJ1x/34asTrSzJao127dtC0aVPmzfsC2+yNyO2cNfLR\nqNjnFWvy1Vdf4ePjzZtvlkh65zSg11fg7NkzdOrUie7dO7Jhwwbs7f5rZMQ9J0SfN6z0PGEMY5N1\nABCzFqISNUQk6pjbnICnJMFaH91yyRFxzzFufI0bt6LMW9u4uPsrar79MSGhYVb5b5i55gJPH/6O\nS/4yrDvxnOLFra/7M+ZV8cg5R5jk8T6OjrmoXbchb7UcTd6ilQkJDcPONuIiY1yPqlZZ1wBffPEF\n773nSXCwzipjDAkNI6tjblzbTWDft59z7N2G6eaZ4JDQMF74PyC7Sz5Cw8Ktsn4TsmrVd7Rv706W\nwg8p7doh+nlha/w3+Pn5cenSJerWbWpl8dmxfv02vvrqC2rWrM369ZvxuvsMiH1OiGhxR6SjiPOJ\ndU6gktILCWO7wU8BrQGUUrWBmLPZ/wqUUUrlUko5ENEFfsbIcqzaznXeGAzh3Lt20Gq/EADPH/+J\nU55ilg4jRTp27Exrj6lc3P0VAX6/A9Y/e9y9e3fZunUr/fsPsnQoScpVUNGww1A8PN5J9BEjazLJ\nw5UXT/7CMXchq/57i0+ePHmo/850fruwk58PLeNp4AtLh5Sg/fv30LJlS6u8jWNnZ8eECZP58MPx\ndOjQiuDHt6z6nGBKxibrbcBLpdQpIgaXjVJK9VBK9Y+8Tz0a2A+cBpZpmpZuVxNIbL5yOzs7Krl1\n4v6vx2I9v2pNgxueBb7iRcA/OOXKZ7UnuISmEfSZ9j4Vmwzg3NZptK9iHaOqE+PtvQBPT89YM7pZ\nM1WtKd269aBPn3d5+fKlpcNJlqd3f6ZEOev8HifFKdcb1HlnJgH/3ObP82utdjT+rl3f8/bbb1s6\njER1796TefMWcXzLl/wb+CTWtow6LalR3eCaphmAwXHevhFj+y5gVyriSjcKl6nBtdNbmDlzGlOm\nTLN0OK/xDwrmZeBjsjhadwKJ7+p45poLFFT1QGdDu47tOeR7gFKlyry2nzU80uXn58emTRu4fv26\nxWJIjkkervT74nD06/Dwaty6dZORI4fi4/ONVY5niBIWFsa9Wz/i2ryPpUMxSkTdB+DafgK39n/C\nunWr6dnTup4Y8Pd/xvnz59i+fSvWfv3WokUrhg0ewNq147naphgVK1aK3pYRW9syKUoqzRneiOO+\nu9mzZycrVkQ8Tm4tI8SjRkv++/wRWXJYd7JOTMGybri18qRh44Y07tgvwYn+LenrrxfTsWNnsy14\nYEp5nLOSxzniuV8bGxsWLPDh999vM3eul4UjS9xPP12kRNGCfPeVh6VDMUrU2Az7rI64dZ7Ep59O\n5fHjxxaOKrYDB/ZTt249nJzSx8Cs0aPHMXXqdLp2bc/+/XstHY5ZSbI2gTx58vDtt5uZM2cWnuPn\nW83jBKEhr/D/5zb/BviRzdm659iN73ZDzO6sbxdNoVmvmTy+f4sJEz6MNYucpS+Onj17ypo1K/jg\ng5EWiyE1smXLxqpV61m7dhXff78VsK5bOVF8fffStGkLS4dhErnzFaNXrz707duT4GDLX9hH2bVr\nR7p4LC6mjh07s3btRsaOHYm390IMBoPVLLdsSpKsTaREiZKsWvUt+zfM4m/tpKXDAcD+wUFOrB2N\nvlgV7LPksPjFgzFiDipzyFWc6u3G8+OPP7BgwRzAOp61XLZsCS1atKJo0fQxiC++E1m+fPlZtWo9\nEyaM4dKlHy0UWcI2b/6OxV8v5XpQaUuHYrS491InTpxK3rx6RowYYhVTGAcFBXHixDFatGhp6VBS\nrHr1GuzZc5BNmzYwevQwXr16Fe9+1ngRmlySrE2oWjVXmveawdUjSwl/FcizwFcW/WJcvRoxSL90\nrS4Wi8FUopKyzi4btTtPYdWq5WzcuN7SYfH8+XOWLfua4cNHWzqUVKtUqTJz5iyiT5+esRZWsfQJ\nbulSH2bOnIZb109xylPEYnGYQsyLTxsbGxYvXsIff9zhyy9jL/VpiTo/cuQQVatWTzcDJOMqXLgI\nO3fu5/HjR3Tr1pEbNzRLh2RSkqxNzD5ncQqq+lz2XUToK8uN0Hj58iVnz56h+/hN5MxXGjtbm3Q5\nMjK+VqCjS17Wr9/CJ59M5s8blm0Frlmzkjp16lGmTFmLxmEqrVq1oVT1tpzdOpOX/76weG9MQIA/\nn38+g5079+OYp2jk5DkZR7Zs2Vi9+js2b97Ehg3rLBrL7t07aNOmnUVjSC1HR0dWrFhHgwaN6NCh\nJfv27Ym13dK3zFJDkrUJRbX+ytXrxauQEC4c+dZiX4zTp09SoUJFFoxxJ49zVnI7p+/RkXG7EJUq\nx7Jlq9m3biYBfncsEtPLly/x8VnIyJFjLFK+uVRr9A4ub5Tk0t55GMLDLXaCe/78OVOnTsTNrS4z\nN96Ofr/bxPT7oEl8F596vZ5vv93E9OlTOXHiGJD2SSU4OJiDB31p3bptmpVpLra2towePY7167fw\n4Ycjoi+CrOGWWWpIsjYDO/usVGzowT3ttMW+GIcO+dK0afM0L9ec4k6KUqdOXRbMncf5bTMIefEk\nzXsONmxAJbrzAAAgAElEQVRYR8WKldJsrd+0Mrl3Daq6Dyb0VRCn9nwTfYIb4HUkzWI4ffokjRu7\nERYWxqJFX6dZuZZStqxi6dKVDBzYj9Fem9M8qZw8eQylypEvX/40KS8tVKlSje3b9+Dl9Tne3gst\nHU6qSbI2oZitv4oVKhESHMSTv3+1SCwHD/pGj5z1GuLGsskZYxRtXB06vE2JCvW5dvjrBAeVmENI\nSAiLFs1j5MixaVZmWsrl7EizHlP59dJR/rp2KM3KDQoK4qOPPmTw4PeZOfNLFi78HzO//QXHrP+t\nBrfxs/Tf+otP3br1+fjjT9mxbCLBQc/StOzdu3emu1HgyVG6dBl27NjHt9+uxvDHLmxtSLeTpUiy\nNrGo1t/k92ryVrOB/LxvPiPfTtt1d2/fvsWLFy9iTRKQEST0OEZZ19bcvXGOWbM+S7NYtm7dRJEi\nRalZs1aalZnWsmZ3pnnPT/nl+Gqe3r+R9AdS6fTpkzRqVIeAgACOHTtDixatYm230UX8ZGTvvPMu\n/fu+x7lt03n5PG16i0JDQ9m3b3eG6AKPT6FChfn++32cOHGU64f/h1M2Y5fEsCxJ1mb0ZtVG5C9W\ngWnTJqdpuYcOHaBp0+ZWPRuVKS39uDtXrtzk229Xc+XKz2YvLzw8nAUL5jBiRMa6Vx1T9IVR9jeo\n3Hww57Z8wp+/nDJLt2xQUBATJ45l0CBPZsz4ksWLl5AzZ67o7VH3b3M5ZSWXU1aTl29tQgq0IG+R\nSvzsu4j+sw6bvbxz585QsGBhihUrbvayLCVPnjxs2bKT588ecmzTZ+lmet2YJFmbWNzWX82Wgzh4\n0JfDhw+aveyoxz1idoFnFvny5WPq1E8ZMWIoz54luchbquzevRNHR0caNrTO5Q1NJWpATv7Stane\nbjw3z2wweRlnzpyicWM3nj17xrFjZ3B3j92ajjkoKKeTQ4ab6CI+Op2OcvU8CA56yl/XzD9OICOM\nAk8OR0cnzh/z5a3Senr27Mrz59a0oljSUpyslVLZlFJblFLHlVK7lVJ549lnvlLqglLqiFLqsFLK\n2TThpi9eQ9yYP6YF8+YtZtSoD3j69EnSH0ql0FcvOX/+HA0bNjJ7WdbmnXfepXr1Gnh69o61sIop\nGQwG5s37ipEjx2aanguAPIUr8G/AP/Rq8Nqfu1GCgoKYNGkcAwf2Y/r0z/H2XkquXLlNcuz0bpKH\nKw4O9rzlPoxrx1by8OEDs5UVHh6eYe9XxydLliwsWbKCEiVK8fbbbXn06JGlQ0o2Y1rWg4GfNU1r\nAKwG4uvjrQa00DStsaZpTTRNC0hNkOldgwaNaNeuAxMmmLfb1D8omN9+/ZG33qqCs7OLWcuyRjqd\njs8/9+Lly5csXepjljIOHz5AaGjIay3AjCjmgEl9LkcqN+hB8+YNGDTIk59+uhhr35RM4nH27Gma\nNKnLkydPOHbsDC1btk5WDOlxUFBquLxRkqKVWtChe1+zzXB26dKPODk5UbasMsvxrZGtrS1eXnNp\n3Lgp7du7c/fuX2Ytz1QT3BiTrOsC+yJf7wOaxdyolLIBygBLlVInlVJ9UxdixjBp0idcvXqFbds2\nm+X4Ud2F93+7gC5XebOUkR58tOQ8hesMZN68r/j1119MemyDwcDcuV8xYsQYbGwyxx2kqAGTXkPc\n2Ld+Pj/8cJnKlavQr58H7dq5s2vXDsLCwpJ1rBcvXjBlygT69+/Dxx/PwMfnm2S1pq19HXNzKl2r\nG/6P77F9+xazHD+iVZ3xu8Dj0ul0fPTRVHr37ku7du7pYrazRM84SilPpdSVmD+ACxDVUg6M/D2m\n7MACoCfQEhiilMpYw5KNkC1bNhYvXsKkSeO5f/9vs0wnaDAY+OfOjxR/s6ZJj5veOOcuyOTJ0xgy\npL9JH+c6c+YU//zzkPbtO5nsmOmNi0tOhgwZxvnzP/P++wNZvHg+tWpV5eLxzTx6+hSDwZBgK9DD\nozsPHjzg2LEzGXbksaksGdsYO1sdWbI48O2qFUyePAE/Pz+TlmEwGCLvV2eOLvD4DBr0ARMmTKZT\npzZcvGieZ9pNNsFN1B9Xcn/Kli27pWzZsjUiX7uULVv2SpztNmXLlnWM8fuXZcuW7ZXEcTONTz75\nxODu7m7oO32fod+n+0167Kb9vA3ZnPIawsPDTXrc9KTj2B2GjmN3GMLDww3t27c3TJgwwWTHbtGi\nhWHp0qUmO15G8e4ob0OBsnUNdg7ZDXYOWQ2TJ0+OtT0kJMSwYMECQ4ECBQyvXr2yUJTpT79P90ef\nI8aNG2fo2rWrSY//888/G4oXL56pzxdRduzYYciaw8Xg/t5nJj3uh/OPGdqO3m5oO3q74cP5x+Ju\nTlHuNeaBs1NAa+AHoBVwPM52BaxXSlUDbIF6wMqkDurnl75G5hmrf/9hfP/9Dh4c306Jqq1M+u/2\nv3uJImVr8OjR89e26fVOGb6OI24FRAwsGzX3KJ9/PpfGjd1wc2tM7dp1UnXsS5d+5OrVayxf/i1+\nfoHRvSKvTx2Z8es5Ln3hclRvO5ZX/wZgeHKVH3/8KboOTp8+ycSJ48idOzfffbedZ89eAql7bCaz\n1HFYWEQPhZ9fIEOHjqFJk7osX76Wdu06mOT469Z9h7t7m0x7voipdu1GuLYbz7HNXnxRzh5Pz4Em\nOW5IaFis1zHrVK9P2Zrhxtx48wEqKKVOAO8D0wCUUqOUUu00TfuFiIFnZ4AjwMrI9wRgb29PxebD\nuX5yHc8e3TPZc6sd+k3mgu9SCpXJ3F3gMen1embPXsAHHwwkMDB1YxznzZvN0KHDcXBwMFF0GUfU\nILDsji6M9WzN3bt/ce/eXQYM6MMHHwxkzJhxbNmyk3Ll0nZyoPQu5mOg2bJlY948byZOHMuTJ49N\ncvw9e3bJ7YhIM9dcIGfB8rh1/4LPvWbz9deLTXJcUw6QTHHLWtO0f4Fu8bw/N8brOcCcVEWWgeXO\nV5Sytd/hp33zqVR+UaqP9+zZU87sWkDBsm58O3ekCSJMnyZ5uEbPXx31h+Hu3gpf371MnjyB+fO9\njTruL79c58KF8/j4fBP9XnpducdcogaALTvwgGvXr9O0aT369u3PvHneZM+e3cLRZQy1atWmQ4dO\nTJo0PtZ30Rh//PE7Dx78Tc2atU0UXcaQ3SUfnQfP4ZtvJhIaGsbQocNTfUxTDY7MHENarcwkD1dK\nu7Yh7FUQbSrZpvp4Yz5fTd5ib1Gt7Ti++PaSCSJMv+IbOTxt2mecOXOK3bt3GnXM+fNnM3DgkOik\nk95X7zGnLNldqNKoF/v3H2X8+EmSqE3so4+mcuHCefbv35uq4+zbt5sWLVpha5v6809GELMF/OWI\ntnz//V5Wr17OggXW0+aUZG0hOR2zka9Iea5cuZyq4wQGBvDTye3kKVzRRJFlPI6OjixevIRx40bx\n8OHDZH/u5s0bfPnlTI4ePUSfPp5mjDD9i+qyDXjxiqLVOmboqSstKUeOHMybt5hx40bh72/8Yh97\n9uzigaG0yZ9ISc9iXugXLFiI77/fy4YN65g/f3aqjpvQmgYpJcnagvIVq4i39wLOnjXuD0bTfsXd\nvTFN3KpQtmanTDlxRFwJ/WHUqFELD4/3GDVqaLInmJg5cxqzZ39Jnz6esSaZycwTdSRGehzSRt26\n9XF3b8XUqRON+vyjR4+4evUKOfKXl9s5icifvwDbtu1mzZpVrFq13NLhSLK2FK8hbpSq0pxiNd5l\n4MB+KZ6OdOfO7XTs2Irhw0dDia7obOwy7cQRyTVmTMSzqqtXr0hy3xs3NM6fP8OkSR8zcODQ17Zn\n5ok6hOVNnTqdkyePc+iQb4o/6+u7l3wlqmDQ2cuFVQzxXeh/teU21TtMZd68r5J13jAnSdYWpNPp\nKFa+HidOnCNr1qzUr1+LzZu/S7Tl9+zZU6ZM+Yhp06bw3Xfb6N69ZxpGnL7Z29uzePESvvjiU27f\nvpXovosWzcPTcyAjRoyJd5YtU3VtZSTS45B2HB2dmDt3EWPGjEhxd/jevbsoVbGemSLLeJzzFGTL\nlp3Mnetl0Ra2JGsLiprZxtnZhc8//4pVq75l0aL5dO3akdu3f3tt/8uXf6JZswY8fvwIX9+jVK5c\nJfo4NjY6SR7JULas4sMPJzB06ABCQ0Nxda382opo9+7dZd++3Xh6DrBQlOmX9DiknQYNGtG8ecsU\ndYc/f/6cU6dOsmjaILmwSoaoc3TJkqXYunUX8+Z9xYoVqRuJbyxJ1hYS3/296tVrcODAMerUcaNP\nn3dj7b9hwzreeacTH3/8Kd7eS8mdO0+CxxGJ69u3P05OzkybNpk///ydvXt3x9ru47OQHj08Yq2p\nLIQ1+vjj6Zw6dSLZ3eFHjhyienVXXFxyyoVVEuKeW0uUKMnWrbtYuHAuy5cvTfN4jJnBTJiRvb09\ngwZ9wFezvfhw8Sk+61+DKVMmcPz4UbZv34tS5SwdYrpnY2PD/PneNGwY8Yxp1PrXDx8+4OefL7Fx\n43qOHz9nyRDTLendSVuOjk7MmbOQ4cMHc+zYGVxccia6/549O2ndOvMt3GEqUQm7c+d2GAyGNO19\nk5a1hSR2fy9HjhzY2mfhzzu/0rFja+7f/5v9+4/Em6jlPqFxChQoyKJFX1OoTA3OXb5NWFgY/fv3\noVevd2jXriP58xewdIhCJEtyu8NfvXrFoUO+tGrVBpBxF0mJeT6N+bp48RJs3bqLxYvnc/To4TSL\nR5K1BSXUDTVzzQWyOuk5uWEi9nkrsHLlt4muTy3dWcZp0aIVVZu8R/CLALy9F6LT6Wje3J2hQ0cA\npluHVghzi+oOP3hwf4L7nDp1glKlSsuFaDLFvK0Y9xZjsWLFef/9QezZY9xES8aQbnArVbRSc7I7\nv0HNeo1jrZ2c0AISwjivyErAk7/x8VnA/v1HKVKkqKVDEiLFokaHDxs2iD17DlKwYKHX9tm1awdt\n2phmERAB7u4tadfOnQCX+tg7ZGPh2FZmLU+StQUllHAnebjy272A6NfGHkckbuaaC9g6OBEW8pJq\nzUe+lqhlwgiRntSv35D8b7agYXN3rvx4gaxZs0ZvG7PoBJu2bOX40biLJIqExLfWQEylSpWha9ce\n+MzvA8DCsalbLCgpRneDK6U6KaXWJbCtv1LqB6XUGaVUG+PDy5wS634x2ULmAgBb+yzUe9eLctWb\nA/91fcsoe5EeVazXDZe8RZg167NY7//zx1WyO+elePESFoosfUrqFuO4cRPJW7QSAD/9dNGssRjV\nslZKzQdaAK+tGqGUyg8MA6oD2YCTSqkDmqa9Sk2g4r9HCaJey4Cy1HkWGPGVzFuoLJN717BwNBmf\n3MIxv4AXr6jQeAAbN46jZcs2bLoQsZ7yjcvHyVdaVthKqaS+q/O2/kLtLp9y+8fvGTr2E04d2GG2\nWIxtWZ8CBgO6eLbVBE5pmhaiaVoAcAuobGQ5mZKM8E47NrrYS9hF9VzI/4FIb6Iu5m2zulCj1WCG\nDx+E35NnPPJ/wd83TpOvdB3pJTKTopVa8NfNi/z++x2zlZFoy1op5QnEXSC5j6ZpG5VSjRL4mBPg\nH+P3QCDhocyR9HqnpHbJVHI7R9xvilkv80Y3ptO4ndGvU0rqODb/oFeEG8DWVode78TYBcejey5m\nrb8U7/9Bckg9x88/KKInwxT1I3X8Onu7/5a7fLNaY3KG3uTU8TUULFcfh6xOOOYujL2dbbLrTuo4\nafNGN6bdmO+xc8jGiGGDWb16KQsXLjRLWYkma03TlgHLUnjMACISdhQn4GlSH/LzC0xhMRlbWGTS\niFsvLjkc4n0/KXq9k9RxDBGtkHAAnLLb4+cXSEhoWPT2kNCwBP8PEiP1HL+Y9T1yzpFU9VZIHcdv\nXI+q0QOixvWoypSAbjz4vh9P7l0nf5k62NnqGNejarLqTuo4+Wwi+5d79uxHgwa1+OCDMdEzTCYm\npRdD5hgNfh6YqZTKAmQF3gSumqGcDE3u66W9pEZ/CmHtYt7SeRRkw1vuwzi35ROqtBwhczGYSXjk\nukv58xegVau2rFy5jNGjx5m8nNRMimKI/AFAKTVKKdVO07SHwALgBHAImCiDy0xHZh0yjYTuScsE\nM+YhYwDSXmF9DvTFquA+ZC3O+uLyFIkZxH1yZ8iQ4SxbtoSXL1+avCyjW9aaph0DjsX4fW6M198A\nllmaRIhkkqSctqS+zS/mhfwkD1f6fXEY+6yOABTLL/egTWms9+nXLoCUKkeVKlXZuHE9vXv3NWl5\nMt2oEDFIz4XISKLup0qPhnm45MjyWo/R0KEj8PFZSHh4uEnLkhnMRKYlSTltSX2nvVxOWXka+FJ6\nNcwgqlUdt27r1KmLk5MT+/btoXXrtiYrT1rWQgiRQXkNcSOXU9akdxQpEnOGw7hd4TqdjqFDR7B4\n8XyTlinJWgghMjC5tWNexfI7vVa/bdq05+HDh5w/f85k5UiyFkIIIVIgqacb7OzsGDx4KN7eC0xW\npiRrIYQQIoWSesyze/denD9/ht9+u2mS8iRZCyGEECaWI0cOevfuh4/PYpMcT2cwGJLey/wMMrWd\necn0gWlD6tn8pI7NT+rYNPz8/HBzq87p0z+i1+tjbdPrneJbCCtB0rIWQgghzECv19O+fSeWLfs6\n1ceSZC2EEEKYyZAhH7Bq1TKCgoJSdRxJ1kIIIYSZlCpVhpo167Bhw9pUHcfoGcyUUp2ALpqm9Yxn\n23ygLhFrWRuAjpqmBRgdpRBCCJFODR06giFD+vPee57Y2RmXdo36VGQybgFcSmCXakALTdOeGBWV\nEEIIkUHUrFmLfPnysWfPTtq372TUMYztBj8FDAZeG82mlLIBygBLlVInlVKmXXpECCGESGeipiA1\n9gmsRJO1UspTKXUlzk91TdM2JvKx7ESsZ90TaAkMUUpVMio6IYQQIgNo2bI1AQEBnDlzyqjPJ9oN\nrmnaMmBZCo/5AligadpLAKXUYeAt4EpiH9LrZa1Vc5M6ThtSz+YndWx+UsemN27cWJYuXUyHDq1S\n/FlzLJGpgPVKqWqALVAPWJnUh+QBfPOSSQ7ShtSz+Ukdm5/UsXm0bNmRyZOncPLkD9SrVyNFn03N\no1uGyB8AlFKjlFLtNE37BVgNnAGOACsj3xNCCCEyrWzZstGvX3+jFviQ6UYzCblSThtSz+YndWx+\nUsfm8+TJY2rVqsqzZ09lulEhhBDCGuXOnYcuXbql+HOSrIUQQog0NGjQByn+jCRrIYQQIg0VK1Y8\nxZ+RZC2EEEJYOUnWQgghhJWTZC2EEEJYOUnWQgghhJWTZC2EEEJYOUnWQgghhJWTZC2EEEJYOUnW\nQgghhJWTZC2EEEJYuRQvkamUcgHWAk6AAzBa07SzcfbpDwwAQoEZmqbtNkGsQgghRKZkTMt6FHBA\n07RGQB9gccyNSqn8wDDADXAHPldKOaQuTCGEECLzSnHLGpgLBEe+tgf+jbO9JnBK07QQIEQpdQuo\nDFwwOkohhBAiE0s0WSulPIGRcd7uo2naj5Et6DXAiDjbnQD/GL8HAi5JxKHT652SEa5IDanjtCH1\nbH5Sx+YndWxdEk3WmqYtA5bFfV8pVQlYD4zRNO1EnM0BRCTsKE7A01TGKYQQQmRaxgwwKw9sArpq\nmnYlnl3OAzOVUlmArMCbwNVURSmEEEJkYsbcs/6MiFHgC5RSAM80TeuklBoF3NI0badSagFwgogB\nbBM1TXtlsoiFEEKITEZnMBgsHYMQQgghEiGTogghhBBWTpK1EEIIYeUkWQshhBBWzpgBZiajlLIB\nvImYNCUYeF/TtN8sGVNGopSqBXyhaVpjpVRpYCUQTsTo/KGapsmABSMppeyB5UAxIAswA/gFqWOT\nUkrZAkuBsoABGETEuWIlUs8mpZR6A/gRaEpE3a5E6thklFIX+W8OktvA56Sgji3dsu4IOGia5gZM\nAGZbOJ4MQyk1joiTXJbIt+YQMTK/AaADOlgqtgyiJ+AXWZ8tiZh2dzZSx6bWFgjXNK0eMJmIp1Gk\nnk0s8uLzayCIiDqV84UJKaWyAmia1jjyx5MU1rGlk3VdYB+ApmnnAFfLhpOh3ALeJuJLAFBN07Tj\nka/3As0sElXGsQmYGvnaBghB6tjkNE37HhgY+WtxIiZYqi71bHJegA9wP/J3+S6b1ltAdqXUfqXU\nIaVUbVJYx5ZO1s5EzHgWJSyya1ykkqZpW4lY9SyKLsbr5yQ9BaxIhKZpQZqmPVdKORGRuCcT++9J\n6thENE0LU0qtBOYD65DvskkppfoQ0UvkG/mWDqljUwsCvDRNcyfiVs66ONuTrGNLJ8a4U5PaaJoW\nbqlgMriY9eoEPLNUIBmFUqoIcBhYrWnaeqSOzUbTtD6AAr4hYmbEKFLPqdcXaK6UOgJUAVYB+hjb\npY5T7waRCVrTtJvAYyBfjO1J1rGlk/UpoDVAZLfAZcuGk6FdUko1jHzdCjie2M4icUqpfIAvME7T\ntJWRb0sdm5hSykMp9VHkr/8CYcAFqWfT0TStoaZpjTRNawz8BPQG9kkdm1RfIsdkKaUKEpGcfVNS\nxxYdDQ5sI+KK7lTk730tGUwGFTW6cAywNHJt8evAZsuFlCFMJKLbaqpSKure9QgipuGVOjadzcBK\npdQxIpbkHQH8inyXzcmAnC9MbRmwQikVlZD7EtG6TnYdy3SjQgghhJWzdDe4EEIIIZIgyVoIIYSw\ncpKshRBCCCsnyVoIIYSwcpKshRBCCCsnyVoIIYSwcpKshciklFIuSqltlo5DCJE0SdZCZF65iJhe\nUghh5SRZC5F5LQAKKqW2WDoQIUTiJFkLkXkNA/7WNK2zpQMRQiROkrUQmZcu6V2EENZAkrUQQghh\n5SRZC5F5hWL5lfeEEMkgyVqIzOsB8KdS6pClAxFCJE6WyBRCCCGsXKq6wJRStYAvNE1rHOf9UYAn\n4Bf51kBN026kpiwhhBAiszI6WSulxgG9gOfxbK4GeGiadsnY4wshhBAiQmruWd8C3ib+xz+qAxOV\nUieUUhNSUYYQQgiR6RndstY0batSqngCm9cDi4FAYJtSqo2mabsTOpbBYDDodPLIpxBCiEwjRUnP\nXI9tzNc0LQBAKbUbqAokmKx1Oh1+foFmCkUA6PVOUsdpQOrZ/KSOzU/q2Pz0eqcU7W/yZK2UcgEu\nK6XKAy+AJsAyU5cjhBBCZBamSNYGAKVUD8BR07SlkfepjwDBwEFN0/aZoBwhhBAiU7KW56wN0uVi\nXtKtlTakns1P6tj8pI7NT693StE9a5nBTAghhLBykqyFEEIIKyfJWgghhLBykqyFEEIIKyfJWggh\nhLBykqyFEEIIKyfJWgghhLBykqyFEEIIKyfJOh0a632asd6nLR2GEEKINCLJWgghhLBykqyFEEII\nK2euJTJFMkV1Z3sNcUv2Z/yDgs0VjhBCCCskLet0ZuaaC4SGGQgNMzBzzQVLhyOEECINSMvaQqJa\n1NJKFkIIkZRUtayVUrWUUkfieb+dUuq8Uuq0Uur91JSRkfkHBae4lTzJwxU7Wx12tjomebgmuyzP\nGb4yglwIIdIpo5O1UmocsBTIEud9e2AO0BxoCAxQSr2RmiDTk/geqzL2UauEPueSIwsuObLE8wkh\nhBAZUWpa1reAt4G4C2i/CdzSNM1f07QQ4CTQIBXlWKXUPusc1f2dVCv5aeBLaRELIUQmZ3Sy1jRt\nKxAazyZnwD/G74GAi7HlZAT+QcGx7k3PXHOB4OBgnvnd5Z/ff+Llo5vcvv0b//zzDy9evMBgMMT6\n/N0b51i0aD4//HCOV69eGRXDk4BguT8uhBDplDkGmPkDTjF+dwKeJvUhvd4pqV2siq1tRIdC3Lj9\ng17Fen/sguOEhhkIDw/jw682Ub9MOEc2f4/280lsHbKRzSkPNoTx7jFvAgMDCQwMJCQkBCcnJ4LD\n7LB1yEbw88eU7N2TSZO2cPPmTerXr88333wTb515zvAFYNnkFtHvRcQQDsCs9ZfwGp78jo74jicS\nl96+y+mR1LH5SR1bF3Mk61+BMkqpXEAQEV3gXkl9yM8v0AyhmM+TgJdA7LgjHquKSIoj5xxhkocr\nwa9eceXQEu5eP4xzrjf4t6Eb3ds14GA5d5zzFo2+9xzzOeuQkBCmLz/BzT/+IfTVv9jaOTBzZk8A\nAgL8Wbx4PhUqVCRfvny4udXDza0ederUI1++fISFGV6LKyQ0LNbrlNR1fMcTCdPrnaSuzEzq2Pyk\njs0vpRdDpkjWBgClVA/AUdO0pUqp0cB+IrrZl2madt8E5ViNqGedo17Hvd8cHhbKnV/OMnr0Kvbt\n24Mhqx73gd+wcmr76H20RCZDsbe3J2t2Z7I7xy5zkocrzs4ufPTRVMaNm8S1a1c4deokW7ZsZOzY\nUej1emxzKirU7RLreJM8XBngdRQwpGgEOcijZUIIYQ1Slaw1TfsdcIt8vT7G+7uAXamKLB0a3aUC\n7fpN4eYP2ylfphjV23di2LBReO99kOJjTfJwpd8XhxPcbmtrS+XKVahcuQqDB39AWFgY169fZdDE\neexaMgzXbR9Tq1ad6B9e/M2L5/68fFmHrFmzJiuGpC5KhBBCpA2ZFMUIES3VI9GvAwL8WbHiG5Ys\n8SFLrhI06jqBVZ/HfLz89WSdnOlF8zhn5XHAy2Q9U21ra8uOn0JQdT0o69aT3Pb+VMkXyLlzZ1i4\ncC737j8gp74Y5TZNp2rVatStW5969RpQrZorDg4OKfr3CyGESFuSrI3kkiMLL4OeMXPmNNasWUHj\nxs3YtOl7Vhz1T/rDKWCjI8XPVOt0NuTOVwwPD1c8PPoAkCtXNp4+/ZfAwADOnj3NyZMnmDLlI27d\nukn16jWoV68+devWp0qVakxc+gMQcUER86LEnIyZI10IITILSdZGCAsL49fzO/jp6Dq6d+3C/v1H\nKVasOABe5V/f39gE5DXELUXPWMdt8cdkZxfxX+3k5Ezz5i1p3rwlAP7+zzhz5jSnTh1n7NhR/Pnn\nH3OoNiIAACAASURBVOQs8CZVm/bBYKgjk68IIYQVkGSdQj/8cI4JEz7EyckJ3737efPNeLKzCaU0\n0ac0ubq45KRly9a0bNkagCdPHtNl8AwOrJ1C7f1fkTXfW1SMM2DNHGQgmxBCJExW3YohsVnJHj16\nxMiRQ+nXz4PBgz9g27bdZk/UluCz+w7Fq7ajaf9vqNNpAjqdju0L+9OlSwe+/noxt2//ZvIyU7qS\nWGpnjxNCiPRGWtZJCAsLY/XqFXh5fUbnzt04deoHnJ2td0I2U93z1el06AuV5uj3K3n+/DnHjx/l\n4MH9LFo0H0dHR5o1c6d5c3dq13Yz6QC1Px7Is51CCBGXJOsY4nbFXrx4gfHjx5AtWzY2b95J+fIV\nLBRZ2onvvrejoyOtW7eldeu2GAwGrlz5mQMH9vPZZ9O4desWDRo0onlzd5o2bcEbb6R8zZaoMkPD\nDMnqxpcucyFEZiPJOlLMZ4on+xzi31vf4+u7j6lTp9OlyzvodHHXK8m4EkuYOp0u+vnuMWPG4+fn\nx6FDvhw86MvUqRMpWbIkzZq506JFSypVegsbm+TdaXHJkYXHAS+TTMTy7LcQIjPKdMk6sUeEDOFh\n/Hn1IIfObsDj3e5W3+VtDfR6Pd2796R7956EhIRw7twZDhzYz+DB7xMQEECzZi1o3rwlDRs2wtEx\n4en1opJ01H1rScJCCPEfGWAGGAwGKjr/yZEVQ7j3y1F279jFjBlfSqJOIXt7e+rVa8C0aTM5ffpH\ndu7cT4UKFVm1ahmVK5eLNUjt+fPnr60ulhyTPFyxi1xE5VmgcSuQCSFEepOpWtZjvU/zNPAluZz+\nm27zzp3bTJgwhgcPHlC/04fkK1aRihUrWTBKyzPVILUSJUrSv/9g+vcf/NogtYcPH1C0aLHIx8ba\n4D2qPkPmngCSnoDFJUcWnga+NEmMQgiRHmSIlnVKHuUJN0R0uQYHBzN79pe0atWEBg0ac/DgcfIX\nr5Sp7k2npahBanPmLOTyZY2//vJj1ar15MqVm2nTplChQiku75+PISjl86gLIURGl6la1lH3RR/c\n+ZnK1T6gVvXKHDx4gsKFi1g4ssxlnM8ZIKIFX6FCRcaMGc/ff9+j19CpnN29iHv961OoUOFEjxF1\n0SWEEJlBpkrWQYFPuX5sBU/uXqdpl+Es/3JErO0yL7XlFCxYiCqNPDi5fTbNmtUna9Zs1KxZi5o1\na1OzZm3efLNC9JSpOZ0ceBzwMnowWtS9a/n/E0JkVBkiWSenhTV1yVHObv6EPEUq0vC9BSyf2iYN\nIhMpMWdEIxjRCIPBwG+/3eKHH85x/vxZli9fyv3796lWzZWaNWvxR0BuQuwLY58lu6VDFkKINGFU\nslZK2QDeQGUgGHhf07TfYmwfBXgCfpFvDdQ07UYqY41Xcp67vXbtKuvnDiJvMVfKN+yLvV2GuFWf\nbiV1caXT6ShdugylS5ehR49eADx+/JgLF85z/vxZ7v+8jR8vXqJCoz78b/QMPph/Ki3CFkIIizG2\nZd0RcNA0zU0pVQuYHflelGqAh6Zpl1IbYGpt2bKRyZPH8/mnX+B7541kz5IlzMPYSU3y5MmDu3sr\nDv7mQrlW7rxR9Tbndi+ijCqFvmQtnPIWY+LiED4b2sKc4QshhEUY28SsC+wD0DTtHBD3jFsdmKiU\nOqGUmpCK+JIU9dytna0u1ok/JCSESZPG8eWXM9m8eSdduryDS44s5HHOKvc2M4Dc+UvSynMO3Uf6\nkN0lP88e3GDNrD5Ur9+Cd4bPJiQkxNIhCiGEyRibrJ2BgBi/h0V2jUdZDwwEmgD1lFJmvUHskiML\nLjmyRD/C9fDhA95+uy2//34HX9+jVKhQ0ZzFixRI6OIqufyDgmN1o385oi3l6nSmRtvR/HJNo4iq\nw9VTm6hatTzTp0/l5s2k777IKl5CiP+3d9/hTZbrA8e/6aClpS2jkS1QxIcpq6XMMmRLWQJS2UNQ\nBA8oVRSRw1FUqKCAooclS/iBehwIKAhFZItshEdlqOA5UKCUUiilbX5/JCmhdKYjabg/1+VlkvfN\nk6cvSe7cz3R29jaDXwNs145001qn2tyfo7W+BqCUWg80BNZnVaDRmPlSlNlxd78zN/rCH8fp3HkE\no0aN4tVXX71rbeqlUzvZ/RquIC/XOD+V9jcvSpPb+kTO3Z7WhG4sVZyo58LuKq9atQrUDOlMzZDO\nRPZ5kCVLltC792M89NBDDB8+nH79+lGiRIl7yrW+f/Lr+jjLdXZlco0Lnlxj52JvsN4JhAOfKqWa\nAkesB5RSAcARpVRt4Abm7HpxdgXGxNi/NWKK5Qv890Ob+XnzYpYsXEj79p24fDnB7jJdjdHol6dr\nnJ+s/165rc/t5JS7blufb1velWvmlc3KlKlIZOQUxo+fxPffb2L16hU8//wLPPZYOE8+OZiQkCZp\n872tmbq1PNuV7nLbZeJM19lVyTUueHKNC15ufwwZ7FmfWSll4M5ocIBhmPupS2itFyqlIoAJmEeK\nf6+1npZNkaa8vjG2bv2ewcOG0OrJt1n91qA8leWKXOXDZ92+c0Fk27THrE3YJf2Kceq8uXemekX/\ne5rZL1y4wNq1q1m1ajlubm74V21F4EMt8Cxe+q7nSLB2bnKNC55c44JnNPrlarlMuzJrrbUJeCbd\nw7/aHF+Nud8632W0a9aGDd8wZtyzBHd/meIBFWXXJheW0Uh+63th+or9WT63bNmyjBs3nrFj/8G+\nfXuZMOUdjmwfh0/JClSsFUY5/y4FUmchhMirIr8oyoUL/+OVVyLpMug1bnlXc3R1hANNHhSclnln\n9WPNYDAQGtqUuh2epfajT3PpjyP89cs21s4ewentIcR518MY1IQ4N1knXgjhHIr06iCxsVfo168n\nQ4YM58N/jsjTKGPhGqwzA7Jjne+danLjgWoNaddvEocPn8SzfCh/ndjO9wtHsPermUSM/RdXrlwu\nhJoLIUTmilxmbR0MlJCQwJNP9iUsrC3jx08EMm4iFa6lIOfI+/r6UrPRo3iWDeHWjatcOL2fP07u\nIyRkAY0bB9O9ey+6du1G6dJlsiwno64aR3PGOgkhcq5IZdbWbOjWrSRad+pGjRoPM23adNnWUqSJ\nGtM8RwHJdr637Y886+NePiUJqt+BXd9/w5EjmoEDhxAdvYWQkPr07duDlSuXcfly/mbcMt9bCJGZ\nIpdZx/59kmPRCwk0lmP27Hl3zaOWrEHkhjVIp3/fBPh6ERufmHbc19eX7t170b17LxISEtiyZRNf\nf/0lU6dOplGjxpaMO5wyZbLOuB1JthMVomgrUpl19wbF+Onrt6hQoynRG75I2zJRiMJiDdyLFi3j\nyBHNoEFD+eGHaJo0MWfcCxcu5GJMjF3BMf3qbPnF2iJl3VJUCFH0FJlod/bsGQYM6Euzx8ZStU4r\nvL29HV0l4aJy2kKTPuPeunUzM+Yt4dSx3Xj5lqTT3lZMm/gUISFNcHd3BzLvO7Z3gxMhxP2hSATr\nixcv0q9fTyZMiOSXm7UdXR3hIvKz28TX15fw8J6cuF6FE2evEH/pLLcuHGTSpBeIiblI586P8dhj\n3UhJLoa7h2e+vW5O5HRKmxDCeTllsLbNPuLjrxER8Th9+jzBsGEjHVwzIbIW9VwYvV78mtLlglgw\nawQAZ86cZsOGb5g1ayaHjx6nYo0Q1lUeRtu27dPWKi/ogCozJYQo2py6z/rmzZsMGfIkjRoFExn5\nsqOrI0SOpJ/rXa1aEM8++xzr12+m3bB5+JdTLF/+MY88ohg8uD8rVizl77/P53iOuBDi/uO0wTol\n+TYjRw7GaDTy9tvvyPQsUeRNX7Efj+KlqFyvMw26T+HAgWOEh/dkx44faNu2Oev+PZZff97In3/+\nkW1ZMs1LiPuLUzaDx8bfYP83s3ikWgDvv78gbXCOEEVBTvvCS5YsRd++/enbtz8pKSk8OeE9zh7f\nTufO7ShRogRhYW0JC2tNy5Zh2S7Ekl91EkI4J6cL1m8s28eBjfNISrxO9TZT8PQs3ME4QhSUrPql\n3d3dWTP3BeAFTCYTJ0+eYPv2aNasWcWECeOoVi2IVq1aExbWhtDQZjJvWoj7jFMFa5PJRPQXc7kR\nd4HQ3lPx8Czm6CoJka9y0idtMBhYEn0VaMgnnzzL7du3OXDgZ7Zvj2b27JkcOHQIf2N1AqvUZ9wb\n53h3UkSRWnNAlj4VIvec6hP+1luv437jL5o9PgVPL2+ZZiIE4OnpSWhoU0JDmxIZ+TJTF/3I/n27\nufTnEbZ8OouqC16ibJV6jBnamz3nAggIrMw7z7ZwdLXTSHAWIu/sCtZKKTdgPvAIcAsYqbU+ZXM8\nHJgCJANLtNaLsitz/fp1LFi8lPCnP8DbN8Ceagnh9HIasLJq5p42shWjYpOpWCOEBZFtGTtzPf89\nc4hjx47y/TffkpqaQuKJjmnN5uXLV8iv6mcot8FYmvCFyD17R4P3BIpprZsDk4BZ1gNKKU9gNtAB\naA2MUko9kFVh27ZtIzLyH7SLmCqBWtz3sloe1DoK3HaaV/ESpQiq15Z3332f9k8toMUT0wkJCWXT\npm9p06YZTZs25B//GMOqVSs4deo3TCZTruuU09Hn2Z0nS58KgV0zOextBm8BfAugtd6rlLJtr64F\n/K61jgNQSu0AwoDPMiusX79+/PvfH/P1UfNgMmkuEyLnrJnq9BX7SUkFb/9y7LxoYPHi4aSmpnLy\n5An27t3N9u3beOedt0lMTCQ0tJmlab0Zdes+YtdATusXTnaZsmTSQuSdvZm1P3DN5n6KpWnceizO\n5lg8kGW6vGbNGlq1al1gGxkIUZTYbt+Z3bgN20z1XEzCPcfd3NyoXbsOw4aN5KOPFnPgwHE2bdpG\nt27dOX36FOPHj+Xhh6vw+OPhzJgxnW3btnL9enyO6xqXcOuuTDn9Zzh9Jp2bv00IcYe9mfU1wM/m\nvpvWOtVyOy7dMT8gNqvC2rZtS+Tc7WkbGcxcfZCo58LsrJrIjNHol/1JIs/y4zqX9vfOsKy4hCQA\nvpgZDkDk3O1px6qV9+f3c3Ekp6QChkzrYTTWpmHD2jz9tHk51NjYWHbv3s2OHTuYN28WBw4cQClF\ny5YtadmyJWFhYWmva1tmXEJS2mcW4O9LCfd8hj097qyR4OnhjtHol+nflhvyXi54co0LTlxCEuEv\nfLVr3aweOW5GtjdY7wTCgU+VUk2BIzbHTgI1lFKlgATMTeBR2RV4OznlrtsxMTn/dS+yZzT6yTUt\nBPl1nVMsQc+2LHOWav5NPH52NJMHBfNiRMO0udu2t5NTUtPOyZ4HISGtCAlpxYQJcOvWLQ4fPsSe\nPbtYtGgJAwYPx7dkeXwCytEt5iDzpz3N7M+Op9UFwMPdQIVAX06dNze4WT/D6esXExOf4d+WG/Je\nLnhyjQuOzee4WW6eZ2+w/gLooJTaabk/TCkVAZTQWi9USj0PfIe5mX2x1vq/2RUoOwMJYZ/8Xk/c\ny8uLJk1CadIkFJjAwClrib1wlsT4Sxze8SX16s2mXFAjSlQK4YFqjfDw9CY11ZTpZ1jWOxf3m4KY\nrmhXsNZam4Bn0j38q83xb4BvcluufKiFMMvoQ56TH7QBvl5cvpaYr33CPv5GivkG4uFuYMHqmVy6\ndIkNG9YxY97HHNn0PoGV61KpdhsSEkJz9Bm29wvM+gW4dGqnXD9HBq0KZ2H9HCenmHbn5nlOtSiK\nECJr2QXDqDHN83WDD+sAMVuBgYEMHjyMo9cViTeucfLwDk7t/4o6dd7H2788xsq1+NT4J40bh1Ct\nWpAESiEsrJ/NAF8vlk7tlKsPhlMFa/lQC5F7jvzc3DJ5UblOOyrXaUeVB7w5/dsvxPx1gu++28j0\n6dO4fPU6xko1iejRjvr1G1C/fkPKli1n12vJTBFRVBTEe9WpgrUQImuFHZhtm94XRLbN8lwPz2Is\ne2vkXY+NeesrLp3XJCXdZNGif3PkyCG8vLxp1649ISGhBASUJDS0GYGBgVmWbZvhR87dzosRDXNU\nfwnworDZvlet0xWtUlNSMLi52bUcsARrIVxMfgf0zJrera+TVT+6b4AR3wAjky3nmkwmzpw5xZYt\nm5m7+DNu3Yzn2oVnqFixItduF8c3oCyvv/gULVuGUaxY3jbysfQLAvd+aQpRmCLn7+Ly1Xh+XPUS\nSYnxcCacZcsW56oMexdFEUIIgLuWPk0v/SIpBoOBoKCHeOqpZ2j1+Eu0H/gGWp/lvfc+oHqDDviX\nqUhU1FvUq1eDsWNH8+23G0hMTLxrMRXbNRhGztjKyBlbC/xvFCKn0i/8Exd3lUPbVrJt+QR8S1Wk\nSa/XOHM194OpJVgLIQpEduuAWwO5h4cHjRoFE1SvLXVb9GHjxi1ER++iQYOGfPTR+9StW4NRo4YS\ne3YvPh6pGbxS1muSZzcyPqfrnguRU9YfsNHRW2jTpjnXLp+nVqvBNOz6PH5lKtOodd9clynN4EKI\nLBVEP3lG/Xq2GXiFChUZOfJpRo58mpiYGDZu/IbZHy5j19fvceNERzp06EqnTl1IzWRPkoxGsQtR\nWG7fusn+zYvYsugI7733AZ/sI+39aO+0SsmshRB5EjWmeabzwu1Z4zx9Fm40Ghk8eBgdBr3J4+OX\n0bNnT77++gtq11H89NVbXDi1j9j4G3bXX/YkEPlpz55d7Fw1gcY1SrFt2y52/OmXLz8cJVgLIQpM\nZv3Z9mzoETWmOXNf6MSQIUNYsWINncd8TNmgYH7/6T9899EIdm1YwG+//ZpWvm0dMmMdhCZbdoq8\nSkxMZOrUyTz11FBef/1t5syZj79//m35LMFaCOEQtoE8p8E7cv4uRryxCYBqlYw8WK8DLfq/TbO+\nr5OaaqJD54489lgHVq5cRkrSDdwMsn6DyH/pxzkcOnSA9u1bcf78ObZt203nzl3Tjlnf22BuAs9u\nCmRmpM9aCOEUcrvc8ORBwQx/2zwSvETpStQKG0KXfs8QWj6WVatWsHnrVspVD2V3mBdNmzbHYDBk\nWpZs2Snscfv2bWbPnsmyZUt4882Z9Oz5eIbnBfh6ERufmKcltSVYCyEKTEFntW6W+OtmufHa0KYA\n7L9oJLlKH86f+IFho54mwMediIiB9OsXQYUKFYHsB6HJuuIiM3EJt4i7eJZOnV6kfPnyREfvtHtl\nvpySYC2EcIj0QTAnQdE8EOzeDDmjjMXLpyRBjXvQPnwgXesaWLVqJW3aNKNRo2AiIgaSfLtgv1yF\na3p96V5O7FrLmQPrCOs+mpVzX82y1Qby5wef9FkLIYqEOyPGU9MGgy16qR2LXmp3z7m2feCvDg6h\nUaNg3nnnPQ4dOknv3n1Zvnwpq2ZEcGzLh8T+fRJ/n7ytliaKjrzMq9f6JGvmjeXyuWO0GjiLOk26\nZBuo84tk1kKIIi+jzCWjbNvHx4d+/SLo1y+Cc+f+YsTEKI5unsf/firOu0kR9O3bn0qVKgOyrri4\nIyUlhQ8/fJ8PPniPSZOmsONCVQyGwh3nkOtgrZQqDqwEjEA8MERrfSndOXOAFpbjJqCn1vpa3qsr\nhLhf3dlUJH++JCtVqswjYf2p1+oJ+od6snbtatq3b0X16jXwqRiCW+m6FPcrU2DrikufuGPk9kfY\nqVO/MW7cM3h5efHtt9FUqVKVYw5Y8c6ezPoZ4LDW+l9KqSeAV4Hx6c5pBHTUWl/JawWFEMIqwNcL\nd/f8bXY0GAwEBzchOLgJyZV68t/TB/n7152c3rQSkymVv4Pb07PRROrUqZuvrysKX0Yr52X2oyk1\nNZVFiz5i9uyZTJw4ieHDR+Hm5rieY3uCdQtghuX2t8AU24NKKTegBrBQKVUWWKy1/jhPtRRCiFzK\nScaa/hx3D08qPdyE1e+NZ1RUNAlX/0c9v1MMHNiPwEAjERED6d27DyVLlsq2bMmciybzznCnmTBh\nLCkpKWzYsIWgoOqOrlbWwVopNYJ7s+YLgLVJOx5Iv0SLDzAXmG0pP1optV9rfTSr1zIa/XJaZ2En\nucaFQ65zwVk6tVOBlm/N2o1GP0r7e1PavyqzXh3FzJlvsmXLFpYsWcJbb/2Lrl27Mnz4cNq1a5dp\nthV7/VZaWRmJS0jK8rijOWu98uK959vS68V1abfh7n+HHTt2MGbMGI4dO0ZUVBTjx4/H3d39nnIK\n+n2YkSyDtdZ6MXDXpptKqc8B67+iH3A13dNuAHO11omW87cC9YEsg3VMTHzOay1yzWj0k2tcCOQ6\nF7yCvMZXriUC5u+jFEtzqfW1GjZsxrx5zYiNvcJ//vMpzz8/kdjYK/TvP4D+/Qfw4INV7i7MdPfz\nbZmbY807iI2fHe10C7K48vs4wNc88j8mJj7t3+FWwlVqh3Ti6vljPNxiKAN7vc3gwa24csX+Neez\nk9sfQ/Y0wO8ErGupdQG2pzuugB1KKTellCfQEvjZjtcRQohCk92WnlalSpVmxIjRbNnyI8uWreLq\n1Vg6dmzN44935/PP13Lz5k0AUk1kuiuYKDzZTdWKu3iaH5b/g+IlSrJz509UrRuGm9u92bSj2ROs\nPwTqKKV+BEYC0wCUUhOUUuFa6xPAcmA3EA0stTwmhBBFRma7idmqV68+b74ZxaFDJxk8eChr166m\nQYOaPNpjAFf/9xsmU8aB356NTET+SkpKonjMVvZ+PpX6HUbz7dpFlCjhvE3/uR5gprW+CfTL4PF3\nbW7PxtxnLYQQRcKdqWHkOoB6e3vTo0dvevTozfnz5xjw3JscWP8O7h5e/FW/A6M6V8FoNN71nJys\nEy2D1PJf1Jjm/PzzTzz6aEuqVKlK+OgPKFHygbTjzjq/XlYwE0IIi8y29MyNihUr0aHncNoO/5A6\n7Z4iKe5PmjVrxMCB/Vi37kuSkpLyqbYiJ2z3K9+3by/Dhw9i8OAIJk6cxIoVa+4K1DntCnEECdZC\nCJHPJg8KxmBwI7ByPbZv/JxDh36hW7ceLF68gAYNavHaa69w9eIfjq6my7MG38v/PUVwy3Y888wI\nYlIr0m7oPHr06F1oS4XmBwnWQghhkZN+6pxyM9zZFaxECT/69x/Al19u4JtvNuHl5cWez6dybN0U\nPvlkOdevZzzy2jYrFPb572972Pv5NKrXa8Hu3Qeo2aQ73r53Zhzb/ps781gCCdZCCFEASvl5U8rP\n+57Hg4KqM3nyVA4e/IXx4yfy3XcbadiwDuPHP8tPP+3FZLqzwpazNskWBXFxV/l7zwJO/riM0F4v\ns+bD1ylWrFi2P4DyoyukIEiwFkKIApA+S08/hcjDw4NOnbqwfPlqduz4iaCghxg37mlatWrC/Pnz\nuBEfe0+Zedkx6n4y8p+raNC4MSVKlKDHmA8JergBULR/AEmwFkIIBytbtizPPTeB3bsP8M47czhx\n4jhrZg/n2JaPSE267nRNss7syy8/Z9OySTR8dAgzZszGs9i9rRtZyc+ukPwkW2QKIYSTMBgMNG3a\nnKZNm3PtWhyd+45g88LRjL24niFDhmMymYrUoKjCFBt7hUmTXuDYsaOE9nqVUuVrAHdPe8vL9DxH\nk8xaCCEKQW4Hi/n7B9Cix/P0fm4JtWrVYezY0Xw5/xl2fvMR+/aZ+7alWdxsy5ZNtGnTnAceKEvH\nEXPwe+ChTJu6nbVPOjuSWQshRAHLaGvGnLiTFXYhtkRT9u3dyeVzvzBw2DDKB/rjV7U1pcpWIzGx\nEd7e2Tf3utoiK9evX2fq1Mls27aFDz5YQMuWYZYA7Xoj6CWzFkKIIsDg5kbgg/VRzSMY/NIypk+f\nyfkzR9m1bh4NGtRkypRJ/PqrdnQ1C82ePbtp27Y5ycm32bZtFy1bhgHOPf0qLySzFkKIApYffaW2\nZbw6OITpK/bTsOtEAAK9r1P8xiF69+5GUFB1Bg0aSnh4z3uy7cuWXcXSK0oZd2JiIjNmTOezz9YQ\nFfUenTt3veecrJq5i8LfmBHJrIUQohDkR19pZmUElKnAK6+8xsGDvzBq1Bg+/fT/0rJtrU8C3NV/\n6+zTljLriz969DAdO7bm7NkzREfvyjBQuyrJrIUQogjKKFv39PSkW7fudOvWnT/+OMsnnyzn8cfD\nqVYtiJLV2+BWpgHunvcG+8JYJS198M1NhpucnMzcubNZtOgjpk17kz59nrjvRsUbrKvlOJjJVTc6\ndxauvJm8M5HrXPDkGt+Rk+br27dv8913G1mx4mN27N5HpVqtef2FQZQrV968xefKnzl1/hoA1Sv6\nM3lQMEajH0OnfZdt2fbU1SqrckdFRZOcYqKMvzejOxoZN240vr5+zJnzARUrVsqX+jia0eiXq18b\ndmfWSqleQB+t9YAMjj0FjAKSgTe01uvtfR0hhBD2s2bbh2MrYAj6jT+PbmZ85Mt4u9/Gx8eHwIda\nUbxCKMX9jdkXlgfW7D27rgDryPn4S39y9uAhwud9xcSJLzNs2Ejc3O7fnlu7grVSag7QETiYwbFy\nwDigMVAc2KGU2qy1ln3hhBDCgXwCylKz5UA83Afh71OMvsHufPrp/7F65fMEGKvQYswQTp8uSWBg\ng1yXnVWWbzt17fK1RDzcM04qr1y5zLefvIk+8D0Y3ChdQbFh/WaCgh7KdX1cjb2Z9U7gC2B0Bsea\nADu11reB20qp34FHAOce0SCEEEVMbpqoJw8KZvjbW4E72W1oaFNCQ5uS+mBPzv32E0eOHGbOnNn4\n+BTH5K+o1rgnkfPzfwR1coqJUVHRLIhsm/bYunVf8vLLkfTq1YfydbrgXbIi3t7FJVBbZBmslVIj\ngPHpHh6qtV6rlGqTydP8gDib+/FAQCbnCiGEKCRl/DNeOMXdoxhVarUgakwkJpOJF2d/zo4t6/jx\nk4k07DyWM2fK4e7ujru7Ox4eHri5uePu7mZz253Yawm4uWccUl4Z2JjBUz/jj2PbKO5vxCegHHEX\nTzFq1DKSkm7j5ubGr7+e5OOPVxISEprWZ21dhcyV5kvbK8tgrbVeDCzOZZnXMAdsKz/g3u1jdx6l\ncwAACFZJREFU0jEa/bI7ReSRXOPCIde54Mk1to+7pfn5yjVz/7H1Oloft943VqhOrbAhlKpQk1+2\nryDiyBpSUlJITk6+6//W27eSkklJSaZWq8FMXWCib3N/jh8/zrFjxzh+/DjHjx8nITGZB6oFE3P2\nZ24nJeJvrEr3AQP5Z9QiShof5PDhQzbzwg2Audnc08Nd/r3Jw2hwS2Y9Wmsdke7xssBmIATwBvYA\n9bPps5bR4AVMRtAWDrnOBU+usf0i5+8iLuFWWv+xdfR3ekajH+EvfJV2P7PzrKav2M/BQ0f4ceXz\nFPP2oX69utSsWZuaNWuiVC2UqsXLS45iMBhwczPc9fpX482hIf12ota+bdumcldSaKPBMf/sSYv0\nSqkJwO9a63VKqbnAj5gXXXlFBpcJIYTjRY1pzvQV+9OmamXFw/1OUM3O5EHBjPpfPF3HrWTJ5Mcy\nnANtwoDJlLOVuKLGNJcNStKxO1hrrX8AfrC5/67N7UXAorxVTQghRH7L6dKnAb5eadltTvqMA3y9\nwNcrw0Btu2Kau7sbySkpafcLY0EWV3D/TloTQoj7VE6WPo0a05wy/t45WiI1cv4uYuMzXnc8K+di\nEu4aSJb+9YvqOt4FQYK1EEKIAmW7E9aHz7dOe7yS0deBtSpaJFgLIYTIs1J+3llmwrbZvJvB/J+r\nbmdZEGQjDyGEuM/ktHk5p+fltt+5lN+d+d553YnsfiHBWgghhN1slxLN6QIm0hedexKshRBCFLjM\nArQE7pyRPmshhBB2k37nwiGZtRBCiDyRfueCJ5m1EEII4eQksxZCCJEn0u9c8CSzFkIIIZycBGsh\nhBDCyUmwFkIIIZycBGshhBDCyUmwFkIIIZyc3aPBlVK9gD5a6wEZHJsDtADiARPQU2ud/W7nQggh\nhLiHXcHaEow7AgczOaUR0FFrfcXeigkhhBDCzN5m8J3AM4Ah/QGllBtQA1iolNqhlBqWh/oJIYQQ\n970sM2ul1AhgfLqHh2qt1yql2mTyNB9gLjDbUn60Umq/1vpoXisrhBBC3I+yDNZa68XA4lyWeQOY\nq7VOBFBKbQXqA1kFa4PR6JfLlxG5Jde4cMh1LnhyjQueXGPnUhCjwRWwQynlppTyBFoCPxfA6wgh\nhBD3hbysDW6y/AeAUmoC8LvWep1SajmwG7gNLNVan8hbNYUQQoj7l8FkMmV/lhBCCCEcRhZFEUII\nIZycBGshhBDCyUmwFkIIIZycBGshhBDCyeVlNHieWVY7mw88AtwCRmqtTzmyTq7GMn1uCVAF8ALe\n0Fqvc2ytXJNS6gHM0xQf1Vr/6uj6uCKl1MtAOOAJvK+1XubgKrkUy3fyIuBhIBV4SmutHVsr16GU\nCgXe1lq3VUo9BCzFfJ2PAc9qrTMd8e3ozLonUExr3RyYBMxycH1c0QAgRmsdBnQG3ndwfVyS5UfR\nv4EER9fFVVlWTWxm+b5oAwQ5tEKuqSPgq7VuCfwLmO7g+rgMpdSLwELMSROYV/l8xfLdbAB6ZPV8\nRwfrFsC3AFrrvUCwY6vjkj4FXrPcdgOSHVgXVxYFfAj819EVcWEdgaNKqS+BdcDXDq6PK7oJBCil\nDEAAkOTg+riS34He3NlTo5HWervl9kagfVZPdnSw9gdst85MsTTDiHyitU7QWl9XSvlhDtyTHV0n\nV6OUGoq59WKT5aF7NrgR+cIINAb6AE8Dnzi2Oi5pJ+ANnMTcUjTPsdVxHVrr/3B3smT7PXEd84+j\nTDk6MF4DbBegddNapzqqMq5KKVUZ2Aos11r/n6Pr44KGAR2UUtFAA2CZUqqsg+vkii4Bm7TWyZYx\nAYlKqUBHV8rFvAjs1For7ryXizm4Tq7KNtb5AVezOtnRwXon0BVAKdUUOOLY6rgeS9DYBLyotV7q\n4Oq4JK11a611G611W+AQMFhrfcHR9XJBOzCPu0ApVQHwBS47tEaux5c7rZ2xmAfyuTuuOi7toFKq\nteV2F2B7Vic7dDQ48AXmjGSn5b7sfZ3/XsHcvPKaUsrad93FuiuaEEWF1nq9UipMKbUPc6IxJqvR\ns8IuUcDHSqkfMQfql7XWNx1cJ1djfc++ACy0tFz8AnyW1ZNkbXAhhBDCyTm6GVwIIYQQ2ZBgLYQQ\nQjg5CdZCCCGEk5NgLYQQQjg5CdZCCCGEk5NgLYQQQjg5CdZCuCilVFWl1BlH10MIkXcSrIUQQggn\nJ4uiCOECLNtHvmq5WwnYB7yBeQnDaKAu5uUje2qtryilxgIDMS8vmQo8obU+qZR6B/PuPynAV1rr\nfymlSgAfAHUwLz05Q9aYF6JwSWYthOtoCowGamHeOakb5p2qZmmt6wEXgP6WHdh6AK0tj38JjFFK\nPQh01lo3AJoDDymlvDD/CNivtQ4GWgOTlVLVCvlvE+K+JsFaCNfxvdb6lGW97BVAO+C81nq/5fhx\nIFBrHQ88CTyplHoLCMecYZ8HbiqldgATgCla61uYM+2nlVIHgR8AH6B2Yf5hQtzvJFgL4Tps98p1\nB26ne8wEGCxbpu7BvJ/8emAp5u1pU4BQYApQBtitlKqB+XtigNa6oda6IdAC805uQohCIsFaCNfR\nVilVTinlBgwCNnL3BvdWwcBvWus5wE+Yt6l1V0o9gjlz3q61jsS8E5DCvBf6GAClVHngIOZ+cSFE\nIZFgLYTrOA98grm5+xzwPXe247MyYc6K3ZRSx4DvMAfoqlrrI8Bu4JhS6mfgDLABmAYUV0odBbZg\n3htdpoQJUYhkNLgQLsAyGvwlrXUXR9dFCJH/JLMWwjWYuDeLFkK4CMmshRBCCCcnmbUQQgjh5CRY\nCyGEEE5OgrUQQgjh5CRYCyGEEE5OgrUQQgjh5P4fs3a0khKZ94sAAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 5 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Notice there are some issues here. In the upper plot, the data: the data is sparse in some locations, which makes for a poor model. Also, at the edges there is some extrapolation which does not look good.\n", "\n", "In the bottom, phased plot, things look a bit better. But there is a discontinuity between the left and right edges because the algorithm is not considering points across the boundary.\n", "\n", "We can improve on these by passing a ``period`` argument, which enforces periodicity of the smoothed model:" ] }, { "cell_type": "code", "collapsed": false, "input": [ "model = SuperSmoother(period=period)\n", "\n", "y_t = model.fit(t, y, dy).predict(t_fit)\n", "y_phase = model.fit(phase, y, dy).predict(phase_fit)\n", "\n", "ax = plot_data()\n", "ax[0].plot(t_fit, y_t, '-k', lw=1)\n", "ax[1].plot(phase_fit, y_phase, '-k', lw=1);" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAesAAAFkCAYAAAAAFROsAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd0FNXbwPHvpkAqoYXQO1y69NB7UUEU7EKUokgTKVKD\nKALSQUABERTBV0QFVAQD0juC9DZAAqFIL0ko6fv+kfJbQkjZbJlsns85npPs7M59uG7mmbnVYDQa\nEUIIIYR+Odk7ACGEEEKkTZK1EEIIoXOSrIUQQgidk2QthBBC6JwkayGEEELnJFkLIYQQOueSlQ8r\npfyByZqmtUzx+mCgF3Az8aX3NU07k5WyhBBCiJzK7GStlBoOdAPup3K4NhCgadohc88vhBBCiARZ\naQY/B3QBDKkcqwOMVkrtUEqNzEIZQgghRI5ndrLWNG0VEPuUw8uB94FWQBOlVAdzyxFCCCFyuiz1\nWadhtqZp4QBKqbVALWDt095sNBqNBkNqD+hCCCGEQ8pU0rN4slZK+QBHlVJVgIckPF0vTuszBoOB\nmzcjLB2KMOHr6y11bANSz9YndWx9UsfW5+vrnan3WyJZGwGUUm8CXpqmfZPYT70FiAI2apoWZIFy\nhBBCiBzJoJNdt4xyF2ddcqdsG1LP1id1bH1Sx9bn6+udqWZwWRRFCCGE0DlJ1kIIIYTOSbIWQggh\ndE6StRBCCKFzkqyFEEIInZNkLYQQQuicJGshhBBC5yRZCyGEEDonyVoIIYTQuWydrIfN282webvt\nHYYQQghhVdk6WQshhBA5gSRrIYQQWSKtnNZnrf2sbSLsQZS9QxBCCCGsLts+WU9cdoDYOCOxcUYm\nLjtg73CEEEIIq8nWT9ZCZGdJzYbT+jWycyRCZI20clpflp6slVL+Sqktqbz+glLqH6XUbqXUu1kp\nw5Rpv0hgQF1cnA24OBsIDKhrqSKEsDvp/xPZibRy2obZyVopNRz4Bsid4nVXYCbQFmgO9FZKFcpK\nkE/j45kbH8/c6b9RCB0KexBF2IMoSc5WIHUqHE1WnqzPAV0AQ4rXKwPnNE0L0zQtBtgJNMtCOfKH\nJxyO6dNIyibEpCQuLCs2NtbeITgk01bOexHRcq22ErP7rDVNW6WUKp3KoTxAmMnvEYCPueWYCnsQ\nRXTkfd599x3y5cvP1atXePHFLsDjfX7SF5g1RqORqVM/5/TpU1SuXIVHjx4xduxnGAwp78tEVtw4\n/y8h//6Ob6laxEZF8LB7TWatPElsnBFISOjSxWOepJudU6dO8uGHfalevSbLln3HX39tok6denaO\nzrFER0dzYtMCDAYDNy+folTlJiDXXouzxgCzMMDb5Hdv4G56H/L19X7qsbAH0cTGxQPwzx9TKVa8\nFC1bVuPhw3ACA4cTEPAG3t7/+7yzsyHdc+ZEGa2PBQsWsGnTej788EN27tzJ0qVLefbZNnTs2NHK\nETqGjNRzxxpxLPh4DmXqvMjD8BuEHvmLtwf6Ur/1W8nvcXVxlu/wU6RVL8PmbCc2zkh05H06denL\noAHv4eLiwttvv83nn39K+TajMTg5sXhMOxtGnP1k9LvXvXt3Ih/cpWi5WkTHxnN8zypiYiZTtGhR\nK0eYs1gjWZ8GKiil8gEPSGgCn5beh27ejEj19YTmwoREffc/jfu3L9Os/3S6dvWna9deDB06kLff\n7s7XX3+X/OQXl/hk8rRz5kS+vt4Zqo/79+8zceLnLFz4HXXr1qdDh5dp1qwNAwZ8QMWK1cmXL78N\nos2+MlrPH3/8CQ06DCBf6frExhkpW6cTe38ezdBeL3HucsL3ePibteQ7nIr06jgmNg6Ayyc341dC\n0b//0ITXY2Lo3LkDh7b+yDPNu0rdpiGj3+Pz50P4448/OHToFLNWniT4Sji5PPJSr+mz7N+xFVdX\nVxtEmz1l9kbcEvOsjQBKqTeVUu8l9lMPAdYDu4HFmqZdzXIhRiMnti6iSrMAPu7un/z6hAlT0DSN\n335bmfya9PmZ75tv5lOvXn3q1q2f/Npzz3WgffvnGDVqmB0jcxybN//NtWtXKaEa4OOZGxdnAz4F\nivLdosW8/35PPFziZOBkFgQG1MUY+5Bz/6xk3vRxya+7urqyePFSTu79jUsXztgxQsfxySeB9Os3\nEA8Pj+TXKjR4DScnZxYv/tqOkTmeLCVrTdMuaJrWKPHn5ZqmfZP485+aptXXNK2upmnzs1JG0uCF\n8GuniYt+QNW6bR877u7uzsyZcwgMHM6KFT/KNIIsiIyMZPHihQwePPyJYyNHfszBgwd44YX2PHjw\nwA7ROY55875kxIhAnF0SnjqSZjW0atWGtm3b8+fXH3DnWoido8zebpzZTtEyNalRo+Zjry/acJkq\nzXuxc8XHvD14ip2icwwhIcHs37+X997rC/zvWu3q4syvP3zLnDmzmDr1c+Lj4+0cqWPIFiuY+Xjm\nJvTQGqo0eAknJ+cnjtepU48vv/yaSZPGs33jb6meQ0aUp++XX36iRo1nqFSp8hPHvLy82LJlN0WL\nFmXIkAEYjUY7RJj9hYQEc/LkMTp2fJFp/Ro9MQhy8uQZVKz7PFtWjOfu3Tt2ijJ7i4+PR/tnDZXq\nv5Dq8eJVWvBMu/5s+/0r9u7dY+PoHMe33y6ka9d3HnuqTrrxLFu2PCtWrCIoaB0zZshNkSVki2T9\nIPwW1y4cZfmXY546wrtVq7asWLGak9u+J+x6sCyWYoalS79LvktOjaenJ198MY/g4GC+/XahDSNz\nHL/+uoIuXV4ld+7Um7ldXV1Z98MMur7WhX793pObIjNs27YFZ9fcFCpZ9YljSU9/xSr68+03i+jd\nuzvXr1+zQ5TZ28OHD/n55+X06PH0Na+qV3+G5ct/ZdmyJWzcuN6G0TmmbJGsQ0/soIRqgJeXV5rv\nU6oS9Z/tw5ENc8njLgMbMiMkJJgrVy7TtGnzNN/n7u7OggWLmTZtErdv37ZRdI7BaDTy228r6dz5\nlcdeT+0Je+zYz7h27Rp//vmHLUN0CKtW/cLQD3ozvX/jVI8nPf21bt2OV155nfHjP7FxhNnfpk0b\neOaZWhQrVjzN9/n5FWbBgsUMGzaY6OhoG0XnmLJFso6/dZAJI/tk6L35yzYEgxOXzux77HUZdJa2\n339fxQsvvIiLS/oTBMqXr0Dnzq8wY8ZkG0TmOE6cOE50dDS1a6ff4uPq6sonn4xnwoRPiImJsUF0\njiEqKor169fxwgsvZej9gwd/xJYtmzh+/JiVI3MMSd2Jv/22KnGNi8elduPZqFETKlZULF/+Q5rn\nFGnTfbK+dOki588H06xZy3TfO3HZAeLiQTV6i4N/f8P9+/eTX5dBZ2n77bdVTzzxpWXo0JH89ttK\nDh8++Njr8of3dL/9tpIXX+yS4cVlWrRoRZkyZZk/f66VI3Mc27ZtplKlKhQpkrE5vt7eeRgyZBhj\nxoyQgVAZFBP1iK1bN/P88xlfdyEw8BOmTJkoXQ5ZoPtk/fvvq3n++RcyNV/Pr2xdipapzvjxY60Y\nmeM4c0bj3r271K/fIMOfKViwIBMmTGHAgPeJjIy0YnSOIakJ/KWXXs7U56ZOncX8+XM5efKElSJz\nLL//vppOndJ+qk759PfOO72IioqScRgZdPnMPurVq0/+/AUy/JkaNWryzjs9GTLkAxmHYaZskKxT\nb25JjekatX8s/4YNG4LYunWz7NCVjt9+W0mnTi/h5JS5r0Pnzq+gVGUmT55gpcgcx6FD/5IrVy6q\nVaueqc+VLFmKMWPG8cEHfaTPLx2RkZFs2BBEx44vZupzLi4ufPnlAqZPn0xw8FkrRecYwh5Ecfbo\n1kzfdAIMHjyMa9euPdEcLl2UGaPrZJ006Klx46YZ/kzS4BEfn7zMnDmXwYMHEB4eJjt0PYXRaMzU\nDZEpg8HAlCkzWbny5+QpMPKHl7rVqzPXBG7qrbcC8PPzY9asdBcCzNG2bt1MlSpVKVy4SKY/W65c\nBT76aCQDBvQhLi7OCtFlfxOXHeDRwwfcDD3KibDMLyWaK1cu5s5dwPjxY7l06WLyOaWLMmN0nazX\nrPmNjh07ZWjQU2patmxN27btGTNmpIUjcwzD5u3m3U/+j8jISLM3NyhYsCDTpn3BwIF9+GTRDvnD\nS0V8fDx//LHarKcRSLgpmjlzLt9/v/iJMQLif37/fRWdOnU2+/M9e/bG3d2dr76aY8GoHMv1kP3k\nL14FNw/z1qyvUqUqffsOZNCg/jJGIJN0nax/+21Vpi9wKfujxo4dz549u2ha/DrT+jWSAVApXDi+\njU6dOmdpR61nn32eBg0asW7pOOLj5akkpX/+2Uu+fPlRqpLZ5yhcuAgTJkyhd+8eyU8lIsGwebsZ\nMnsrGzduyHQTuCknJydmz57H/Plz2LFjmwUjdAyBAXW5dnYXxSs1yVJ3Yv/+A3n48CGzZ8+QLspM\n0G2yDgk5x61bN/H3b5il83h5ebFw4XeMHPkR+/bttVB0juHe/UiCj23jpZcy3wSe0owZcyjt58n5\nA6vkDy+F1at/tUgdd+nyKq+99iYDB/aVp5IUrgT/S7Vq1fHz88vSeUqUKMnXX3/H++/3kJHLKYSH\nh3H70nFUjYx3S6bG2dmZr7/+luXLf2DNmt+lizKDdJust2zZTKtWbXB2fnJ50cyqVasOQ4cO54cf\nlkifaqKJyw5w52owAGuOxGb5fK6ursyYMYfgA7/j7WaNzdyyJ6PRyNq1a7LUPGtq8OBh3Lt3j61b\nN1vkfNmZaSvZJW1PpqYSpaVZsxa88UY3pk+XZTJN7dy5A9/ilcjl5pnlc5UsWYoRIwL55ZefLBBZ\nzqDbZL19+1aaNWthsfO1aNGKNeuCiHz0kNg4I72nbbHYubOrm6FHKFSmTpaawE2VLFkKT59CXLtw\nxCLncwSadho3N3fKli1nkfON/HofufzqEBS01iLncwT37kdyJfhwhtZiyKhu3d4mKGgtsbFZv5F1\nFNu3b+G9bl2euuRzZjVp0pw9e3YxpHNpi53TkekyWcfGxrJ7906aNm1hsXOWK1eBstUac3DtdOlX\nJaH/6c7loxQqXcOiTdZzpk/i6IY5hIScs9g5s7OdO7fRrFnaS7hmVrlnWhMUtI516/606Hmzm7AH\nUdwOjyT8zlXi42L4ZZ/l9qcuW7Y8VatWY/ToYTIvOJGlH6D8/Pzo338g3bq9nryAlXg6XSbrI0cO\nUbRoMQoVKmTR867/dQnG+Di0nT/k+D6S6Oho7vynUVbVtuh527d/juHDR/P6613kDxDYvn0bTZo0\ns9j5wh5EEeuSh6VLlzN06AccPXrYYufOTpKm/ADcuniUgiVrWKyFKMk33yxh3769rFjxo0XPmx1d\nuXKZu3fvULVqNYue98MPh1KjxjMMHjzAoud1RGYla6WUk1JqgVJqt1Jqi1KqXIrjg5VSxxOPbVFK\nVczM+XfssPzTCCT0q7Z8dSShx9bn+MEjBw/+S54Cxcjtbt4UjLQEBHSnQoWKrFu3xuLnzk7i4uLY\ns2cXjRtbJlmbzkldeyyejz4axezZMy1y7uzs1sWjFCr9jMUHNXp756F0w56MHjsux6/PvmPHNpo2\nbZ7phZPSYzAYmDRpOvv27eHYsaMWPbejMbfmXwJyaZrWCBgJzEhxvDYQoGlay8T/zmTm5JZubjEV\nbXCn9DPPc/jvBTl6LvDOndsoUuYZq53/3Xf7MHXqpBz9dH306GGKFCmS5RHKT/PKK6+xZ89Ojh3L\neWMEkqb8ODsZuXP5OOUqmbdOQHoKlayKd/6iOX4p0m3btlh0TIApNzc3+vcfyKRJn0mXQxrMTdaN\ngSAATdP2ASlvaesAo5VSO5RSmVqR5OHDhxw8+C8NG6a+vZ0lVPB/lXvXznIt9JTVytC7Xbt2EDjg\nDasN7GjVqg3+/g1y9KpbO3Zst2gTeMo5qT4+eRk7djxDhgwkPj4+x60h4OOZm/iI/8jt5oWnj6/V\nyvF/vh+zZk3j8uVLVitDz4xGY/KTtbX06PEeB45ovNZfRuA/jblzbPIA4Sa/xymlnDRNS5r8uRz4\nCogAViulOmialubwVV/fhObYDRv2ULt2LcqUyfxydhmxemonOg9fQ6XGb3Lj2K8ULNjX4n1depVU\nx5GRkRw+fJCOHdvj7W35ZvAkM2dOo0aNGowaNYwiRTK/BGR2lVTPBw/uo2fPnsm/W0L+PG6PlTFg\nwPssXbqYbdvW4+yc97FjjszX1xtnZwM3L5+gcJnqODsbrPLvDnsQDe6F6NOnD3PnTufbb7+1eBl6\nlVSfrw75nvuPYqlTp5pVr5UNX+jPjtUz8fH5hFy5clmtnOzK3GQdDpj+ZZgmaoDZmqaFAyil1gK1\ngDST9c2bCSM5N27cSp06/sm/W4OPZy6q13+O3T8G8euvf9CiRSurlaUXvr7eyXV6+PBBSpUqQ2Qk\nREZar55z5/bh9de7Ehg4lilTckbfqmk9HzjwL+PHT7XodzkucVCV6TmHDx/DqFEf0eit2Tg5O1v1\nb0cPkuo4Ls7IjYunGNjrFd56q6HF/90JYwQSLmt3fBqwZk1Pdu06QMWKyqLl6JHp9/h66Cl8i1fm\n1i3rdmnlLqBwz1OY2bPn0b17L6uWpQeZvbk0txl8F/A8gFKqAZA8MkAp5QMcU0p5KqUMQCsgw53D\nx44dpXr1GmaGlXFOzs6MGvUxEyeOo9fkTbw7JecsMmGrOgYYOHAIP/60gj4TVtqkPL24fv0aMTHR\nFC9ewuplNW/ekiJFinLsn3U5asGfaf0aERd+3ux17TMjt7sXffsOZNKk8VYvS28unz+Ol18Fq5aR\nNHhSNe7KZxMn8ujRI6uWlx2Zm6xXA5FKqV0kDC4brJR6Uyn1nqZpYSQMOtsCbAeOa5oWlNETnzhx\nLNPbCGZW0vrhHTu+iNFo5OrZPVYtT2+OHj1CjRrWG1xmqmDBglSq34kj2/7PJuXpxbFjR6hevabF\nmw1Trn0PCSNqi9Z+jdO7fiIqKirHDJy8efMmd+/epUKFTE02ybCUYwR69erNwYMHOHgwZ9QvJK50\n+J+Gj18Fm3yv8haugF8JxXffLbJ6WdmNWc3gmqYZgb4pXj5jcnw5Cf3WmTJw+nquXb9F6dJlzQkr\n05ycnAgM/ISefQZQqJy/TcrUg+PHj/Dyy6/ZrLwqDbuwek5PTp8+RaVKlW1Wrj0dPXrEZq0XAEVK\nVcHHrxyhR4KoWKqbzcq1p+PHj1KjxjMWn05kynQ9Bg8PD4YMGc7EiZ+xcuUfVitTT+Lj44i4fZE8\nvmWsWk5gQN3kVSW//WoGL7/8AgEB7+Dtnceq5WYnuloU5c71EPL6lbbqH19KKw464+5dgNAjf9F3\npuPvtBMXF8epUyepVs2yixuk5VGcC+Xqd+Gzzz7OMVMzjh07arPWC0i42FVt2o3g/Svp19E2N7v2\nZo+bv7feCuDy5Yts3rzRpuXay1tN8uPmmQ83dw+rb86TtKFH5cpVaNOmHVOnTrJqedmNrpL11Utn\n8SpQyqZlGgwGqrZ8j7N7fyby/t1035/dp8eEhl6gQIGCeHnZZsRwUl9UqWee5+Axjb/+yhlrWh87\nZruuhiQlSlekTLXmTJz4mU3LtZfTp09SqVIVq5aRstvB1dWV8eMnMXr0MKKiHH98wMmTJ8nvV9rm\nKz6OHTuelSt/zpFrCDyNbpL1xGUHCLt5Ec/8JWza5zZ/SHO8C5SgRLU2uF9z/ESiaaeztK+yuZyc\nXWnRZSBjxoxg0Bebs/UNT3rCw8O4desWZcpYZvOOzKjV8m3Wr1+XI/pVNe0UStm+W6Vdu+eoWLES\nX3012+Zl29qpUyfI52fdJvDUFCxYkI8/HsfQoQOJi/vfXg7Z/WEpK3STrAEibl/Eu6Btn6wBnAyg\nGr7G3r272bVrh83LtyVbPI2YMh2ks2Dce/j7N+ToNsdea/nMGY2KFSvatDsHEp4CZw9tx8cfj2P4\n8CGPXeQcTXx8PJqmUamS7W88ASZOnMLChfMIDb1gl/Jt5fTpU+QtZJtrcspWjDfe6Iq7uwffffeN\nTcrXO90k69Hd6hBx6yL5CpW0et9Ialxc3Rg/fjIjRgwhOjr6qe/L7vthJzyN2PYCZ7q5/KefTuTM\nwSAuhZ61aQy2lJCs7ZNEAF599Q08PT1ZssRxR9RevnyJPHny4OOT1y7llyhRkn79BjJy5FCHHodx\n6tQJJgzsbJctLA0GA9Onz2b69MlcvfqfzcvXG90k6//+u4JrLjd8C1pv2cD0PP98R0qVKs3cubNS\nPW66kUJ2nR5z+vRpu47IXrThEhUbvsHhDfOZ8P0/dovDmk6fPmXXZG0wGJg6dRbTp0/m2rWrdovD\nmhJaiOw7s6Bv3w+4cuUya9b8Ztc4rOXRo0dcuXKZcuXK2y2GChUq0qPHewQGjgCy/8NSVugmWZ8+\nfZK8hUrbpexFI1qxaESr5IvcokULOHnyhF1isabY2FhCQs5RoYJtV2BK2bxVqkZ74mJjOLl/vU3j\nsJUzZ+wzLsCUUpV4++0ejBmTqaX5s43Tp+3TX23K1dWVadNmM2bMSD6c+bfD9aWeO3eWUqVK233p\nzw8/HMrJk8fpNXJOtn9YygrdJOtTp06Rz0Z9I2kpVqw4o0aNZejQD57o80u5SEJ2c/58CIULF8HD\nw8NuMQQG1MXV1YVa7ftxePMSbt26ZbdYrCWhGdz+S1IOGjSMY8eO8PffGV6TKNsIDj5H+fLWXVUr\nI/z9G9C27bMc3PSdvUOxuAsXzlOmjP2nAbq5uTF16iy2rp5DbHTOXdlMN8la007x7mut7NI3klK3\nbu/g6por1W3xTPtfsxs9NB1CQh2WKluZV1553eH6/CIiIrhz5zYlS9r/xtPd3Z2pU2cxYsRQwsLu\n2TsciwoNvUDp0rYfpZyajz/+lNCTuwjWDto7FIvpNWED85ZvpVQpfdRxs2YteLZNK87tXZFtH5ay\nSjfJ+sKF85Qta/upLqlxcnJi5sy5zJgxhZCQYHuHYzEJi0jYbiR4ekaN+pizZzVWrPjRYaZknDp1\ninLlKuDs7GzvUICEdcPbt3+Ojz4a5FA3RRcunNdNsv5qTTDV2/Tl37WzGLtwq73DsZiIu//ppo4B\nxo37nMunthF1x3GuyZmhm2QdGnpBF08jScqXr8DQoSPo27fXY6PDU1ubObtI6Oezb1+qKXd3d+bP\nX8y4cWMIv+0Yoz1PnjypiyZwU2PHjufMmdOsWOEYU+YiIyO5ffsWxYoVt3coyfzK1adQmXps+mWm\nQ9wU3QmP4u7NK5Qpo59kXbBgQRq+MJAdK6cQHh5m73BsThfJOjIykjt3blOkiHX2sDbXu+/2oWBB\nXyZPnmDvUCxCL60Xpjc8VapUZciQ4exYNYX4uFg7R5Z1J0+e1NUNETx+U+QILUXnz5+nWLHiumm9\nSBrLUqNld5yjrmf7m6KkrUHv37vGusMP7B3OY1bMGcorL3Vk2DDHainKCF0k69DQUF398SUxGAzk\nr9mDb7//ga1bs/8WmpcuhVKihH5aL5K8+24fnHJ5su/v7+0dSpadOHHC7qOUU1OlStXklqKYmBh7\nh5MlwcHBlCpV2t5hPMbHMzf58+ZhwYJvE2+KzqX5fr13+8THxRD14A558vnZO5QnjBv3OadOncz2\nN0WZpYtkHRISoqsmcFNunnlp0nkoAwf25ebNm/YOx2zh4eFER0dToEABe4fyhM9/+Jdn2g0k9NhG\n+n262N7hZElwcLAuRtCmplev9ylQoCDTpmXvDRKCg4N11ZdqqkqVqnz00Uj69OmV5uJKehYYUJfo\nB7dx88rP2B4N7R3OE9zd3ZNvinp/9rOub3osyaxkrZRyUkotUErtVkptUUqVS3H8BaXUP4nH303v\nfOfPn9fNqMPUFClbi9dee5OBA/sQHx//2DG93yEnCQ0NpUSJkhbfX9lScnvmpUbb/qxfPol791Lf\nUEXvdW00Grl8+TLFihWzdyipMhgMzJ49n68XfUvAsPn2DsdsISEhNttG1xw9e/bG17cQU6ZMtHco\nZnOKvo1Pfn11S5pKuCkaxY6VU4iLzd4tRRll7pP1S0AuTdMaASOBGUkHlFKuwEygLdAc6K2UKpTW\nyfT8ZJ20Ys6IEYGEhYUxe/aM9D+kQxcuXKBEiZL2DiNVSdMw/MrW5c1XX6ZPn17Exma//uvw8DCc\nnJx0vQevr68vjV8cwo5V07hy5bK9w8m0YfN2szpor+6erE3HYSTdFP3yy09s2PBXqu/X+0pcEXeu\n4p2/iL3DSFPPnu/hkacABzY8uXa43m/szWFusm4MBAFomrYPMJ30Vhk4p2lamKZpMcBOoFlaJ0t4\nstZfsjZdXnTqT0f49ttlLFmyONtt8zhs3m4+/zpIt8kaEjZTcTIkrB0eHx/P6NHDMBqNj/3R6f0C\nd/XqVYoX188I5adx86tGqVod6dr1Ne7fj7B3OJkWfueq7pJ1SgULFmTx4qUMGtSf06dPJb8+bN5u\nek/bovuVuMLv/Id3Pv0+WUPCTVH1dh9wJfgQ33yTfVuKMsrcZJ0HCDf5PU4p5WRyzHRcfQTgk9bJ\nEpJ1aTNDsZ3ChYvw/fc/MmTIAI4ePQzoP4Ekibh7TZeDy1JycXFh0aLv2bt3NwsXzkt+PTusy375\n8kVKlChh7zDSlFR3Zeu8iFOeUvTu3SNbtWLcu/+IiHvXddsSZ6pePX8+/XQi3bq9zu3bt+0dTqaU\n942j35st7B1GmiYuO4CTqyf1XhrD51Omsn596q0YjsLFzM+FA94mvztpmpbUmRuW4pg3kHonZKKQ\nkBBq165GgQLeab3N5r4Y0pLOw9ck/wzQtm1zFi5cSPfub9HhvS+IjUuowqnLDzFtYJoNCHYT9iCa\nOzevUq3aG/j66quOkyRtJ+nr642vrzd//bWORo0aUbJBL4qUr49vPvfk97q6OOvy33Hr1lXKlSun\ny9iSuLokzLgwGAy0fW0Ix9ZNYvLkT5k9W/97Mw+bs537Ybdwze3JV3+c0e3fm6n+/Xtz6VIIffp0\nZ8OGDTg7G8ifx4074Qk3+EnXFb0JDg6mdu1q2eK77OHjx0u9JjB4cH82bNhArVq1CHuQMLhPz/Fn\nlrnJehd7tp/LAAAgAElEQVTwAvCLUqoBcNTk2GmgglIqH/CAhCbwaWmdzGg0Ehfnys2b+muS8/FM\nWMTeNLamTdvSvft7fL1wJHU6j8cllzsxsXG6jD9pzuTDsOus3HWLFi30GyPAoJlbCAyoi6dnAZq9\nOoZVC0fSLGAWMYW9cHFOGBw3/M1auqzrEydOU758WV3GlmT4m7XoOTlhGuKobvUIe+FbOnRoy/z5\ni3jlldftHF3aYmLjeHDvGh4+hXX795aaQYNG0r37W7z77vuE+72EwWBIXrJYj/8Go9FISEgIPj6F\ndBlfkuFv1qL3tC0AzBrWlTXVvOjY8QWe7TWT2LiE+k26nuhRZm8kzG0GXw1EKqV2kTC4bLBS6k2l\n1HuJ/dRDgPXAbmCxpmlp7tNXpkwZ3Y5SfpoPPhhE66b+nNz6TbZYq/Zh+E3y5C9s7zAypXCpypSv\n14V9Kz8lJjpS9+uyh4ZeoFw5+y86kxbTLoSJyw7g45OXOXPmM3bsKE6cOG7HyNIXGFCXqIgbeOUr\nrPu/N1POzs7Mn7+IdX9vI/jgOmLjjLruOrt58yYeHh54een/qdT0mvDCCy/yzjs92bB8kkMumGLW\nk7WmaUagb4qXz5gc/xP4M6Pn09OSdik9bWlRg8HA559PpXqtWjy8fgr4X3NW0oAoPSxLei8impio\nh8THRTOhT2t7h5OqwIC6yXfIphfhwIC6BF8JJ+xGCIbL68H3OXuFmCEJK8Tpd0rR09SuXZdRo8Yy\nbNgg/vxzQ3KXhB7FPrhJQT/9D+JLycvLmw7dx7H8i/74la1HtcrldHvDERqaPb/HAB98MJg///yD\nq9o2SlZtqds6Nocu/io7depk7xDM4uXljf/z/dj751wePdLv1m2Pwv7Dp0AxXbdepPXUXLVlL1as\n+JHb/51N9bgepmkYjUYuXgzV/UUu5c1Qkq5d38ZojGfZsiV2iCrjIu5cxTufvqcUPc30IS9RplZH\njm9eyOhudewdzlNdvBiq6wcoUyn3anBxcWH69C84ue17chv023rxtLUk0qKLZN2jRw97h2C2n2YP\noV2rJvTr954uR9WGPYgi4s5/uOXJfhe4pCbb3B55qdO+N8c3zGTYK+XtHFXqbty4joeHJ97e+m86\nLJDHjQJ53B57zcnJiS++mMeUKRPZuHG9nSJLX8Td/8ij8/m/aSlXrwtRD8P4+OORum2qDQ29kG2S\ndWpq1apDhVrt2PzjWN1OTTx3LvUHj7ToIllnd7NmfUlk5CMGDHifuLg43UznSpru9ODuf3jofM5k\neirVbk3e0k1o2rqdLvdmvnDhQraYfpgWpSrx/fc/8sEHfdi/f5+9w0lV3MObLJ7wlr3DMMvEZQdw\ndnHFv8tYVq/9m4kTx+kyYYeGXtB9C1F6Nq5cyHOtGtGt2+s8fPjQ3uE8QZK1neTOnZtvv/2BGzeu\n07rTm0RFRetqPvCDu1fxzKvvp5HUth5N2s0oaQBfjeZvUbRcLZ57rjW3bt1Kft/t8Ehuh0faOuTH\nXLgQovuFOpKktc1rvXr+zJ27gC6vvsI7I59cGcqeBs5Yz/2HUfj6+to7lCxxdfOic+9p/P13EDNm\nTLF3OE/ITs3gT2MwGJg6dSbFihUnIOCNx64XeiDJ2o7c3d1ZuvQnbt28zomti+wdDvC/Psn7d//D\nK19R3dw8ZIZpX7bBYKBuu/fo3PkVXnzxWSIjI58Y3WwvoaHZ/8k6SZs27WnaeRibf/qMkydP2Duc\nZFcuncfDp7Cux16kxfTmc0Lf1vzyyx+sWvULS5boZ/OauLg4Tp06QeXK+ts5LrOcnJyYPXsepUuX\nZuDAPk+siGhPkqztzMvLi9avB3I9+B/Cr2nci4i2+xfDyRDPg7tX8MzmzeDwv9Xihg0bRfnyFfny\nyy/sHVKy48ePUqlS9r/AJSlWoS61Wr7N8OGD+ejLnXb/HgPcvHiEAsWr2juMLDG9+SxUqBBLl/7E\n5MnjuXbt8dmt9koqR44colAhP91uRpNZLi4uTJo0ndDQCwQFrbN3OMnOnTuT/ptSkGRtYY/iXanS\nvCfHNs4nPs7+A86i74Tg4eOHh1febDmNIanJNuVyoxMnTmHRogW82Thf8nvt9e+Li4tjz57dNGzY\nxC7lW0PYgyj8KrUkJiaa4CMb7R0OADfOH8S3dG17h2FR5ctXICCgB59+GmjvUADYvHkjLVu2sXcY\nFpUrVy4mTZrOmDEjiI22b3cZwO3bt7l6Nc2lR1IlydqCkhJKUdUEV4+8HNz+q90Hml0+8w9+ZeuS\nP49+FxMxR/HiJfjggyG89vZ7yYN07NUMvn//PkqUKImfn59dyre0pO9xnNEJ1bw3/278jqiH9h1V\n2+vztYTdOE+B4lV5bXSGl3DQndTGCwwa9BH79//Djh3b7BTV/2zevJFWrRwrWQM0a9aCOnXqsn/z\n/9n9mrxx43qaNm2e6c9JsrYCg8FAzbbvc/afVTx6EGHXvtTY28eoUM3+i7NkVcrBZgC9e/cl4u51\nboUetmtsf/21lmeffd6uMViLX4mKFKnQgH0bl9q1KfzGhSPkL1YFZ1fHuukE8PT05LPPJhEYOJy4\nuDjAPoMm7969w+nTp/D3b2jTcm2lQI3XOX84iLDbV+16TV6//i+zrheSrC3INKFUrqQoVKYOIQfX\n2C2eixdDuXXrJgWKVbRbDJaUcuEUV1dXJk8Yx5k9P+LsZL9m8KCgtTz3XAe7lG0Npt9jgHL+r3H5\n5FauXrXP/tfD5u3mduhBCiU2gf/8eUe7xGFNzz/fEW/vPKxa9YvdBk1u376Vhg0b4ebmlv6bsyEv\nH19K13yes3tX2C2GqKgotm/fSps27TP9WUnWFpaUUAID6lK58RuEHllHv472mbP499/radWqLTMG\nNGXxmHZ2icHaXnyxC8a4GCKuHLFL+SEhwTx69Ihq1WrYpXxrSfoeh16LwM0zHyWrt0Xb87NdYjEa\njVw59y9+ZWvjlD0HgqfLYDAQGPgJU6d+TlxsjF1icNQm8CSBAXWpWP9Fblw4yKv+9lm8aNeuHShV\niYIFC2b6s5KsrahosVKUqtyEefPm2qX8v/8Ool27Z+1StjWk1t/n5OREzRbdOLxlqV0WmNiyZRMt\nWrTKttOJ0pPUklGubmeunt1DSEiwzWO4e/08Li658HKAGQ1padSoCWXKlKWk4WTya7ZqLTIajQ45\nuCylgvnzU63Ry0yZMtEu5a9fv4727c3rMpNkbWU1mr3J0qXfcv36NZuWe//+ffbt20uLFq1sWq49\nrJg7jOKFvPnjj9U2L3vr1k20bKnPDVKyIunGKGkwTi53b0rX6sg7fT+yeSzBJ/dSsHQtm5drD6NH\nj2XchAnExSTUe9IGN9Z28uQJPDw8KFtW37vGWUKl+p3455+9HD1q27EuRqORDRuCaN/evA2JJFlb\nWMqnP6+8hXjrrbdtcidnOjdz+/at1K5dhzx5fKxerr0ZDAY++WQ8EyZ8SlSU7UZ6RkdHs2vXTpo3\nb5n+m7OhpFHhScrW7sSlcwc5duxoGp+yfAzXQw5SsFRtyhTNw6IRjn3zWbNmbQqXqsL5w2ttWu6m\nTX875E1nStP6NeKLQa0YPHgY48aNtWlr3PHjR3FxcaFiRWXW5yVZW1FS4h40aChBQetsuhpUUJDj\njlBOTdOmzalYUfHtt7ZbInPv3t1UrFiR/PkL2KxMe3Jz92Bs4Bg+/XSMzS5ykQ8jCLsRTMES1W1S\nnh4s+/oLQg78RvSjcJuVuX79Otq10/cWtJYUENCdq1evsHnz3zYrMyhoHc8+28HsLrNMJ2ullLtS\naqVSartSaq1S6omecqXUbKXUAaXUFqXUZqVUHrOicxA+PnkZPPgjPvvsY5uUFxcXx99/B5ndN5Jd\nffLJBObMmcHdu3dsUt7atX/w3HOONzI5iemo8AJ53PDxzJ18kdu0aYNNYqhZ8CYFS1Qjt5tbtlzU\nxxzly1egWKUmnNn7s032GLhx4waadprGjZtatRw9cXV1ZezY8Xz66Rib7ZYYFLQuS7NGzHmy7gsc\n0TStGbAUGJPKe2oD7TRNa6lpWitN02x3i6hT77zTi/PnQ9iyZZPVykhajnP//n/w8ytCyZKlrFaW\nHlWsqOjY8SVmzpxq9bLi4+NZt+5POnbMnnuxZ1Rq0+XGjh3PuHEfP3aRs9bymEFB6yhbtfFT9zp3\nVJUavcGVU9u4f/c/q5e1YcNftGrVmty5c1Ydt2//HAUKFGT58h+sWs6webvp+/lvXLlyifr1G5h9\nHnOSdWMgKPHnIOCx4YNKKSegAvCNUmqnUir7blZtQbly5eLjjz9j3LiPkxc+sCTT5TjHzfyWZ5/N\nOU1aph75tmHJ0h+sPmp5//5/yJ+/AOXKVbBqOXphOhYj6SL344/LrFpmVFQUW7ZsonhFf6uWo0ff\nf/ISFeq/hLZzmdVbFP7660+efdZx1gnIKIPBwKefTmDq1M+5f/++Vcu6pO2lTZv2uLi4mH2ONJO1\nUqqXUuqY6X+AD5D0pByR+LspD2AO0BV4FuinlMo5HU5p6NDhBby8vPj55+VWexIxGuMJPr4rR/7x\nAbh75aNKwy5MnDjOquX8/POPDv9U/TSPX+QSliG1xh7uv/22kqpVq+HulS/9NzugWk1fIfxGMHv3\n7rFaGffv32fPnt20aeOY6zCkp2bN2jRp0oyvvppt1c1TLml7snxNNmR2oIhSaiUwWdO0/UopH2Cn\npmnVTY47AR6apt1P/H0KcEzTtLTaGvS3A7uV7Nu3j5dffplWPefhmsvNoouVdB6+hgtH1uNx/xA7\nduzAySnnjR/sPHwNsTFRHFk5hOXLl9O4cWOLl3Hu3DkaNmzI6dOnKVAgZwwuS023bt0oW7Ysjwq2\n4XToXQAqlcrHtIHNsnxuo9GIUopFixbRrFnWz5cd9ZqwgXOHNxJ1eTt79uyxylz+GTNmsHHjRv76\n6y+Lnzu7CA0NpXbt2rTtNRfPPAUtvoBU07cmc+CPady8FoqXl5fpoUz9DzXnmXwX8DywH3gO2J7i\nuAKWK6VqA85AE2BJeie9edO+GwXYStmyVahXz5/9m1dQqdHrFv13+3jm4uKx9SyYPZ3btx88dszX\n19vh6zihKyAenFyp3KQrH344iHXrNln8IrdgwSIKlW/K0C8PPLFIS06o5yRDhoyiTZumdPmgGpDQ\n3xkTG2eRf/+RI4eIi4unUqWaT5wvp9RxXJyR0lVbcujMBhYvXsqLL3ax6PmNRiMzZsxk+fKVObaO\nATw88tO16zv8tfF7Gr84xKL/7onLDnD2n9VUbPQmYxbuf6xLw9c3c6uomfPoNR+oqpTaAbwLjANQ\nSg1WSr2gadopEgae7QG2AEsSXxOJ8lZ5meB/13A//I7FRnoOm7eb6zduEHH7P5o0yZlPIqYq1W5D\ndHQMv/++yuLn3rVrB0XL54yRyWkpUaIk3bp15+HZNU9sspJVO3fuoHnzlg67MlxGTOvXiOkDmvDp\npxMYP97yawicO3cWg8FA5cpVLHre7OjDD4dwSdvHxfOaRc8bFxvDrYtHKVIh65ujZDpZa5r2SNO0\n1zRNa6ppWhtN024kvj5L07Q1iT/P1DTNX9O0Jpqm2W7iazbhU6AIJaq24syenyx63shbZ2jSuGGW\nBjFkZ6ZTjca8U59PPhnPpEnjLTo1IyzsHsePH8OtQDm7b7WnBx9+OIS//w7C+OCaRUdsJ4xQbmux\n82VnTZs2RynFkiWLLHredevW0L79czn6hijJ3N/PUt7/NY5vXWLRqXLPV3fGM29h3D29s3wjm/M6\nNXUgMKAulRq+yrWzu+nWLPMLuqcm7EEUZ49udeh5vxlhOtWoWbMWFCtWnJ9/Xm6x8//55x8UKlUd\nnN1sMgdW7/Lk8aFPnwEc2fYjYJkpXOfPh6Bpp3LEUrkZNWrUWObO/YIHDx6k/+YMiIuL44cfvufl\nl1+3yPkcQaka7bh/9wpXQo5Z7JybNv1N8fJ1LHIjK8naTnwLFqRSvReYNWtals81cdkBwu9c48aF\nIwRH2WeHL70aMWIMM2ZMsVgT4sKF83imSWeLnMtR9OzZm2sXjnD3+nmLnG/+/Lm8807PHDfvNy3V\nqlWnQYNGFluhb8OGIO7H5OLn/bZZEETvAgPqkitXLio1ep2rh3+1yDnDw8P48cdllK5qmW5JSdZ2\nVKVhZ4KC1hISci7L57p2bh9FKjTAzcM+W7/plb9/AypUqMj//d/SLJ/r7Nkz3Lt3j/nj3rN4H212\n5uXlxYiPPiL+UlCWp3DFxsby668/07Pn+xaM0DEMGzaK+fPnEBGR9TWm/vhjNRVrPytN4CZ8PHNT\nrd6z/PffFXbs2Jbl823ZsomqVaux5PMeTwxENYckazvK7e5N7979mDZtcpbOExhQl1sXD+NXpnaO\nTx6pbaM5cuQYvvhiOo8ePcrSubdt20yrVm1wcnJ6YmWvnK5Hj3fZvG0Hd66FZKl74NSpkxQpUgQ/\nPz8LR5j9KVWJFi1as3Dh/Cyfa//+fRQqVc0CUTkWJ2dnPvpoJJMnT8jy+vfbt2+1aFeOJGs7SUoq\nvXv3Zdu2LZw+bf6A+aFf7uDOldOUVTljG8HMqlmzNjVr1ub77xdn6Tx79+6hQYOs3yE7Ig8PD+q0\nfB1td9YGTf77737q1q1voagcz5Ahw1m8+Oss9V1fu3Y14enc3VcGSaaiS5dXCQu7l+VNPvbs2UXD\nhk0sFJUka7saNm8345YeY8CAQUyd+rnZ57l7LQQ3r/y4eea1YHSOZcSIQObO/cLsZQWNRiN79+6W\nZJ2GH+Z+Sti1M9y/GWx2C8+BA/9Qp049C0fmOMqXr0CDBo354YclZp9j//5/yFtEERdvkEGSJpIe\noJydnRkxIpCpUz/no692mTVg8saNG9y4cYMqVapaLD5J1jrQvXsv9uzZafZ61tcvnqBAcct9KRxR\n1arVqF+/AStW/J9Znz9/PhhnZ+fkzVFSa27P6dzd3anR7A3O7fvZ7HPs27dHnqzTMXDgYObP/5LI\nyEizPr9r13aKla1h4agcS4cOnQgLC+PGRfO2Nd63bw/16/vj7OxssZgkWdtR0mAcDw8PunXrztdf\nf2XWeS6d/RffktUkeaTj/ff7s3DhfLM2Utm2bSvNmrWQATnpqFj7Oe5cC+bQoX8z/dnz50N4+PCh\nLNKRjpo1a1OjRk2zR4Zv3bqZz4a8I4Mk0+Dk5ETv3v04uHU5t8IyP9Zl9+4dNGhg2aWOJVnbieku\nWROXHaB3736sWfN7ct91RuerfvbdHm5dOkH+Es9Ic1Y6/P0bULRoMbMWl9i6dbPM+80AZ9dcVG/6\nhlndOtu2baFFi1ZyQ5QBgYGf8OWXs7h3726mPhcaeoGIiAiqVKkqgyTTcclQlQf3rnHzwqFMX1u3\nbt1My5aWvV5IstYJX19fBgwYxOjRw4iPj8/w566FnsQzXxFyueexYnSOwWAwMH78ZGbOnMaNGzcy\n/LmYmBh27dpBs2YtrRidY5jWrxG/LvgETTvNrl07MvVZuSHKOKUq8dJLLzNy5EeZ+ty2bVto3rxl\njtzkJ7NcXHJRqUkAx7d8Q0xUxp+uL14MJSzsHlWrWnazSfk/ZiemS2MmNUP17t2X+/cj+PPP3zN8\nnpJu/+FXupY0Z2VQtWrVeeWV15k48VMgYy0YBw/+S6lSpfH19bVBhNlf7ty5mTx5Oh9+2I/o6OgM\nfSY2NpadO7fLDVEmjBkzjp07t3Ps2NEMf8b0hkjGXaQtMKAuRSo0II9vaQo8yni3zrZtW2jWzPI3\nRJKs7ShlM5SLiwvDh49m+vTJ3Ap7yO3w9AeQbN26mTKV6klzViYMGjSUoKC1XLiQ9opbSYl869ZN\n8sSXSe3aPUeZMmUzvNTrwYP/UqJESZlfnQkeHh4MHDiYadMy1uWQcEO0jebN5YYoI5KavlXDN5k+\na2by3u3psVYLkSRrnWnduh13Hxq5ejZhw/m0+kru3r3DmTNn8C1R2VbhOYR8+fLTo8e7zJ49I0Pv\nl+ZZ83z00Si++GIGMTEx6b5XbojM8/bbPTl69EiGBvQdPnyQokWL4+dX2AaROQ7vgiUpWaHOUxej\nMW2di4uLY8cOyy6GkkSStc4YDAYatHuHs3tWYDQ+2Xdt+sXYsWMbDRs2YubAFtKclUnvv9+fdevW\n8N+Vi2kuDBH1KAJNO039+g1sGJ1j8PdvQKlSpfnll/QXSpEbIvO4ubnx4YdDmTRpfLorbkkdZ45p\nt+KSeVNZuHAeYWH30vzMkSOHKFKkKIULF7F4PJKs7ehpfUalKtXHJZc7If/+kebn5Y/PfPny5adC\nnQ4c3byYmNi4VFswwh5EEXzqAP7+DWRTCTMNGzaSmTOnER4e9tT3hIXd49Spk/j7Z33P35yoa9e3\nuXHjBj//vPypYzCGzdvNkuW/yfUikwrkcaNAHjfKlatAp06dCQwckeb7t27dTPPm1qljs5O1Uqqz\nUirVFSaUUu8ppfYrpfYopTqYH17OZDAYqN3hI879s5Kw2/89dixpbrbRaExM1q3tFGX2V69NVyLv\n3+bSic3JryVd7JKm1l0/f4gYrwp2jDJ7a9CgES1btmbo0A+f+p4dO7ZTv74/bm5uNozMceTKlYvp\n079I2Ls9JvVWoqhHEdy5fl5uiLLgk08msGPHNvbv3/fY66bX5LVr19CypXWuyWYla6XUbOBz4IkJ\nkUqpwsAHQCOgPTBJKZUrK0HmNIEBdXHP40vZOp24tO/75OYt07nZbw+ZSlxcHOXLSyIx1/1IA9Vb\n90Hb+QPvP1f6ieN3r57h4rENlK7sb/vgHMi4cRM5dOhf3ho8J9WnvpUrf6Z9++ftEJnjqFu3PvXq\n+bNv/XePdesk3Xwe2rGSohUb4u7ubscosx/T1k9PT09Gjx7LyJEfJY/DML0mvztqLnFxcVYbwGfu\nk/UuoC+pJGugPrBL07QYTdPCgXOArG2XSQXyuFG/1ZtcvBjKzp3bnzj+z9/LGD9+kiwgkUX5i5Sn\ndNUmLFw477HXAwPqcu3MTlTD15k+RPavzgoPDw+mTp3FP0Hzn+hXPXbsCPv37+O11960U3SOw692\nAJdObePOtfPJ3Tp3IyK5cTuMkINrKVf/VVk4KYtee+1N3N3dWbduzRPHju7+gw8+GGS1OewuaR1U\nSvUCBqV4ubumaT8rpVo85WPegGkHVQTgY3aEOZiTswsDBw5hxIghdOjQiUaNmrB/71GuntuHc/xD\nOnToZO8QHUKZai1YsWImjx494ty9wuQpWIpVq37h5oWDNOk8zN7hOYSWLVsTFxvDnj++4GuXw3Tv\n/i4hIcEsWPAV3bv3wsvLy94hZnvuXj4Ur9KKY5u+5uGVmlxoWoCI25e4fHonnnmL4JnX8oOechqD\nwcDAgYMZPXo4hw4dpGHDxuzafoRboYe5GnKY1q3bWq3sNJO1pmmLgczuKxhOQsJO4g2kuyaer693\nem/JUZydE56Y+/fvTcGCPpw9e5b33nuHiPsPMBgMbNu6hcKFM7fLltTx48IeROPk5MS6xUNYvboM\nMxb9yf4tUzEYnNgYH01M1CMKlaiQ6XqTek5dvRdHc/vKSX74YQlff/0V169fJyYmhv3790sdW8AX\nQ1pyLvQ6V07vpErhCDp0aMOt23fJ5e5NpSbdcHF24oshGW+ilTpO3VtvvYq7uwtnz57lgw/e5969\nhBHiCxYsoEKFklYrN81kbaZ/gIlKqdyAG1AZOJ7eh27ezNiE85xi8vsJA0Fu3bpP69YdaN0aatas\nz//tDMfNIw9KPZOpOvP19ZY6NpHQ15QwNW7wrK0EBrRl+4V83L0fj2/pmlSpVJ52lQ00adJE6tkC\nJi47gFfB0ngVLI175Wo8WysfHTt2IjIyEnd3d6ljCymYPx8FG73AuO41ufKoILkK++PsmjCTwccz\nV4brTeo4bU2btqVp07aULVuJhX/fxCOPL126tMr09zgzspKsjYn/AaCUGgyc0zRtjVJqDrCDhD7x\n0ZqmZWzNQZGmBg0asfJg5vdWFRkTGFCX0GsJf2yfvdfcztE4rqJlqvHCCwlzWGXAk3V4eHiQt1xL\noqIzv8OcyLgWLVrx4/4tNinL7GStado2YJvJ77NMfl4EZH5rI5EuWfzEMgID6tJ72pbkn5PIsq3W\n8bT6FtZT3NeT4Cvhyb+ntfiPME/SaPCkn6353bZGM7gQ2YIkZtuS+rY+05v5wIC69Jz8vzUEShWW\nPmhLGjZvt01vgCRZC2FCWi6EI3EyQLwR2ZXPSnw8cycnbGvXryRrkWNJYrYtqW/7kVYNy0tK0raq\nW1kbXAghHFQ+bzecZN0kizNducxWTeGSrIUQwkFN69eIfN6y5ro1lSrsbZNWI2kGF0IIBybdD5Zn\nj9kNkqyFEEKITLL1OABpBhdCCCF0TpK1EEIIoXPSDC6EEEJkkq3HAsiTtRBCCKFzkqyFEEIInZNk\nLYQQQuicJGshhBBC5yRZCyGEEDpn9mhwpVRn4BVN07qmcmw20BiIAIzAS5qmhad8nxBCCCHSZ1ay\nTkzG7YBDT3lLbaCdpml3zA1MCCGEEAnMbQbfBfQFntjPRSnlBFQAvlFK7VRK9chCfEIIIUSOl+aT\ntVKqFzAoxcvdNU37WSnV4ikf8wDmADMTz79FKXVA07RjaRRl8PX1zmDIwlxSx7Yh9Wx9UsfWJ3Ws\nL2kma03TFgOLM3nOh8AcTdMiAZRSm4FngLSStRBCCCGewhqjwRWwUynlpJRyBZoA/1qhHCGEECJH\nyMra4MbE/wBQSg0GzmmatkYptRTYA8QASzRNO5W1MIUQQoicy2A0GtN/lxBCCCHsRhZFEUIIIXRO\nkrUQQgihc5KshRBCCJ3LygCzLEtcQGUeUAOIAt7VNC3YnjE5EqWUPzBZ07SWSqnywBIgHjgO9Nc0\nTQYsmClxpsO3QCkgNzABOIXUsUUppZyBb4CKJAxo7UPCtWIJUs8WpZQqRMLMndYk1O0SpI4tRil1\nEO6LxycAACAASURBVAhL/DUEmEQm6tjeT9YvAbk0TWsEjARm2Dkeh6GUGk7CRS534kszgdGapjUj\nYeW5F+0Vm4PoCtxMrM9nga9I+P5KHVtWRyBe07QmwBjgc6SeLS7x5vNr4AEJdSrXCwtSSrkBaJrW\nMvG/XmSyju2drBsDQQCapu0D6to3HIdyDujC/5aEra1p2vbEn/8C2tglKsfxCzA28WcnEqYpSh1b\nmKZpvwPvJ/5aGrgL1JF6trhpwHzgauLv8l22rGcAD6XUeqXUJqVUAzJZx/ZO1nkA09244hKbxkUW\naZq2Cog1ecl0Hff7gI9tI3IsmqY90DTtvlLKm4TEPYbH/56kji1E07Q4pdQSYDbwf8h32aKUUt1J\naCXakPiSAaljS3sATNM0rT0JXTn/l+J4unVs78QYDpguQOukaVq8vYJxcKb16g3cs1cgjkIpVQLY\nDCzVNG05UsdWo2ladxJWR1wEuJkcknrOuh5AW6XUFqAm8D3ga3Jc6jjrzpCYoDVNOwvcBvxMjqdb\nx/ZO1ruA5wESmwWO2jcch3ZIKdU88efngO1pvVmkTSnlB2wAhmuatiTxZaljC1NKBSilRiX++giI\nAw5IPVuOpmnNNU1roWlaS+Aw8DYQJHVsUT1IHJOllCpKQnLekJk6tutocGA1CXd0uxJ/l+00LS9p\ndOFQErYtzQWcBH61X0gOYTQJzVZjlVJJfdcfAnOkji3qV2CJUmob4EpCHZ9GvsvWZESuF5a2GPhO\nKZWUkHuQ8HSd4TqW5UaFEEIInbN3M7gQQggh0iHJWgghhNA5SdZCCCGEzkmyFkIIIXROkrUQQgih\nc5KshRBCCJ2TZC1EDqWU8lFKrbZ3HEKI9EmyFiLnykfC8pJCCJ2TZC1EzjUHKKqUWmnvQIQQaZNk\nLUTO9QHwn6ZpL9s7ECFE2iRZC5FzGdJ/ixBCDyRZCyGEEDonyVqInCsW+++8J4TIAEnWQuRc14CL\nSqlN9g5ECJE22SJTCCGE0LksNYEppfyByZqmtUzx+mCgF3Az8aX3NU07k5WyhBBCiJzK7GStlBoO\ndAPup3K4NhCgadohc88vhBBCiARZ6bM+B3Qh9ekfdYDRSqkdSqmRWShDCCGEyPHMfrLWNG2VUqr0\nUw4vB74CIoDVSqkOmqatfdq5jEaj0WCQKZ9CCCFyjEwlPWtN25itaVo4gFJqLVALeGqyNhgM3LwZ\nYaVQBICvr7fUsQ1IPVuf1LH1SR1bn6+vd6beb/FkrZTyAY4qpaoAD4FWwGJLlyOEEELkFJZI1kYA\npdSbgJemad8k9lNvAaKAjZqmBVmgHCGEECJH0ss8a6M0uViXNGvZhtSz9UkdW5/UsfX5+npnqs9a\nVjATQgghdE6StRBCCKFzkqyFEEIInZNkLYQQQuicJGshhBBC5yRZCyGEEDonyVoIIYTQOUnWQggh\nhM5Jss6Ghs3bzbB5u+0dhhBCCBuRZC2EEELonCRrIYQQQuestUWmyKCk5uxp/Rpl+DNhD6KsFY4Q\nQggdkifrbGbisgPExhmJjTMycdkBe4cjhBDCBuTJ2k6SnqjlKVkIIUR6svRkrZTyV0ptSeX1F5RS\n/yildiul3s1KGY4s7EFUpp+SAwPq4uJswMXZQGBA3QyX1WvCBhlBLoQQ2ZTZyVopNRz4Bsid4nVX\nYCbQFmgO9FZKFcpKkNlJatOqzJ1q9bTP+XjmxsczdyqfEEII4Yiy8mR9DugCpNxAuzJwTtO0ME3T\nYoCdQLMslKNLWZ3rnNT8nd5T8t2ISHkiFkKIHM7sPmtN01YppUqncigPEGbyewTgY245jiBlv/TE\nZQeIuHeL+3cuEfPoHr6+RQgJKYmXlzdeXl64u7tjMPzvHujymX18+eV+/P0b8MwztcyK4U54FGDM\nyj9DCCGEnVhjgFkY4G3yuzdwN70P+fp6p/cWXXF2TkimKeMOexD92OvD5mwnOjqaiNuXCPhwN8r3\nEavWbuP6lRB8CpUht4cPF49c560dXxMeHk5ERASxsbF4e3sTFeeCcy53ou7fpuzbXQkMXMnZs2dp\n2rQpixYtSrXOek3YAMDiMe2SXxs2ZzuxcfEATF1+iGkDM97Qkdr5RNqy23c5O5I6tj6pY32xRrI+\nDVRQSuUDHpDQBD4tvQ/dvBlhhVCs5054JPB43AnTqhKS4qCZWwgMqMut65fYvWIcsdGPKFa6Eg0q\ntmDW5M/4bvt9crl5Jfc9m86zjomJ4bNvd3A29AaxUQ9xds3NxIldAQgPD+Orr2ZTtWo1/Pz8aNSo\nCY0aNaFhwyb4+fkRF2d8Iq6Y2LjHfs5MXad2PvF0vr7eUldWJnVsfVLH1pfZmyFLJGsjgFLqTcBL\n07RvlFJDgPUk9Ikv1jTtqgXK0Y2kuc5JP6fsbw6/Fcr2f3eydv5hwsLCKFCqIdVa9OCb4a2S3/Pn\niacvhuLq6oqbRx488jxeZmBAXfLk8WHUqLEMHx7IiRPH2LVrJytX/sywYYPx9fXFOa+iauNXHjtf\nYEBdek/bChgzNYIcZGqZEELoQZaStaZpF4BGiT8vN3n9T+DPLEWWzRiNRlqVj+KnryZw79o5+vd5\nj2cD+1OjRk1GLNib6fMFBtSl5+TNTz3u7OxMjRo1qVGjJn37DiAuLo6TJ4/TZ/QX/LnwA+qu/gR/\n/4bJ//HwPx7eDyMysiFubm4ZiiG9mxIhhBC2IYuimCHhSTVhevnIt2rx559/8NVXX3Dnzv+3d9/h\nUVXbw8e/k0IglTbSQdpv0yGQECChShdEFBGkSJemNIMIIq9ewQsRFFAUKYKoqFyveLEQWjQCoRm6\nsEGkCCpSkxAJpMz7x2TiEFInZSaT9XkeHmbmzJyzcjKZNXufvde+Ts2GD1N74GxmTOqY5T5yUl60\nnG9JrsUm5GhOtaurK/87lIgKHsL/tRlEWfcYmlWIY+/eKJYufZNLf/xJaWMN6m14FX//5gQHtyUk\npB3NmwdQokSJnP/wQgghCp0kaxv5epbg7NEIQkKm4evry8SJU+jZsxczlu/N1+O4GMj1nGqDwYWy\nFWowZEgAQ4YMA6BMmVLcuHGbuLhY9uzZzc6dPzJ79ov88stpWrQIJCSkLcHBbWnWrDkzV+wHzF8o\nLF9KCrpVbUuNdCGEKC4kWdvg+PFjbP4glKTEBFa+8xbBwW3TplpllGxsTUBh49vkao61dYs/fXJ1\nczP/qn18fOnSpTtdunQHICbmJlFRu9m1K5LQ0ClcuHCe0pXq4//QMEym1lJ8RQghHIAk61yIjY1h\nwYJ5/Pe/GwgNncnQocNxdXUt0GPmNtHnNrn6+ZWme/eedO/eE4Dr16/Rb9xrbP1oNq3C36BkhaY0\nSjdgrSDIQDYhhMicrLplJbOqZCaTic8/X09wcCDx8fFERu5j+PBRBZ6o7eHdb87yoH9vHhq9ktZ9\nZ2AwGNi4dDT9+vVh+fJ3+PXXM/l+zNyuJJbX6nFCCFHUSMs6G8ePH2PGjGncvn2bDz74iICAlvYO\nKUv5dc3XYDBgrFKH779aw61bt4iM/J5t28J5++3FeHt707lzN7p06UarVm3ydYDa+T9lbqcQQqQn\nLWsrMfF30rpjY2NjmDVrOk888QiPP96f8PAIh0/U+SGjVb28vb3p2bMXixYt5cgRzfLlqyldujTz\n5r1Cgwa1GTFiCOvXf8Rff/2Vp2NCzrrxrX9PQghRHEiyTmXpik1MSmHwpHkEBwdy+/ZtfvxxP8OG\njXTKLu/MZLWql8FgoEmTZkyb9gKbN0cQFRVN167d2b59K8HBAXTr1oGwsNc5fPggKSkpuTomZH/t\nOrdd5kII4QyKXTd4ZlOETCYTVy8c4VTUetwMyaxd8zEtWgTaI8QixWg0MmDAIAYMGERiYiJ790ax\ndWs448aNIjY2ls6du9KlS3fat++At3fm5fUsSdqShKUAixBC/ENa1sCZM6c59L/XOLbjfao16MDh\n/XslUdvA3d2dkJB2vPLKXHbv/olNm8Jp2LARa9euokmTevcMUrt16xYmU+5XAbPuMr8Zdze/fwQh\nhHBIxaplHbpsNzfiEijjYy63efv2bRYvXsiaNSuZPPl5arR/HhdX12LV5Z2R/BqkVrNmLUaPHsfo\n0ePuG6R2+fKfVK9eI3Xa2MMsm9KW8W/+CGRfgMXPy4MbcQn5EqMQQhQFTpGsc1P9KsVk7nLdvn0L\nM2Y8T5MmzdixYxeVK1eR6UAFyDJIrWfPXphMJu7evcsvv5xm8+ZveOWV2Vy4cI4y1fxpHPKkvUMV\nQgiH4xTJOqdi4u9wO+4qx79fRcTq8yx/ZymdOnWxd1jFzvR3owDzl6uGDRsxbdoL/P77JQZPeJk9\n37zNpdFtqVKlapb7sHzpEkKI4qBYJetzx77neMRKHmzWk8dHvkynTsH3bJe61PZTuXIVmnUYws6N\nC+ncuS0lS5aiZcsgWrZsRcuWrahfv2FaydTSPiW4FpuQNhjNcu1afn9CCGflFMk6Jy2sEdPf4njE\nSlo/8S98jQ8yZ0Rwtq8RhWvRpA4wqQMmk4kzZ35h//697Nu3h9WrV/DHH3/QvHkALVsGcT62LInu\nVXH38LR3yEIIUShsStZKKRdgGdAEuAOM0lqfsdo+BRgJXEl96Bmt9ak8xpqh7NZcTk5OJixsHt9v\nXEvLvrPxNT6YNppY2Ed2X64MBgN16tSlTp26DBw4GIBr165x4MA+9u3bwx+Hv+Sn6IM07DCM96a+\nxsTFuwojbCGEsBtbp249CpTQWrcBZgAL021vDgzRWndM/VcgiTo7N2/eYPDg/uzZE8XeXVEYqyog\n94tdiPxja1GTcuXK0a1bD2LLdaNejzn0GLmIP05GUlfVJnrzO5w+8DUz39lSgJELIYT92Jqsg4HN\nAFrrvUD6uTYtgJlKqR+VUjPyEF+2MiqPCfDzz8fp2rUDderUZcOGrzAajfh5eVDOt6Rc23QCZSvW\nosfIRQyY/C6efhW5+ecp1i0YRou2XXnyuYUkJibaO0QhhMg3tiZrXyDW6n5yate4xXrgGaATEKKU\netjG4+SIpTymZTWmL7/8D48/3ovp02fyr3/9G3d394I8vMiFzL5c5VT6uuDzJ/WiXuvHCew1lRPH\nNdVUa47t2oC/fwNeffVlTp/OvlNHVvESQjg6WweYxQLWtSNdtNbWhaAXa61jAZRS3wD+wDdZ7dBo\nzLwUZXZcU69BpyQnc2DrKiJ/j2bbtm00a9bsnuetmdPN5mM4g7yc4/xU1tdclCa38YQuiUwbn2As\nU4qw59rds7+aNStTL7A79QK7E9qvOqtXr+axxx6mTp06jBgxgv79++Pt7X3ffi3vn/w6P45ynp2Z\nnOOCJ+fYsdiarHcBvYENSqlWwBHLBqWUH3BEKdUA+Btz63pVdju8csX2pRGTUz/Af/wijNvxMUR8\nt4OyZcvlaZ/Oxmj0cZjzYfl95TaexKTke25bXm+9v+ux5spm5cpVITR0NpMnz2Dbti2sX7+OqVOn\n8fDDvXnqqaEEBrZMm+9taalb9mdd6S63l0wc6Tw7KznHBU/OccHL7ZchW7vBvwQSlFK7MA8um6KU\nGqiUGq21jsE86CwCiASOaa0323icHFkwrjXl4yL466Im4JEZlC1briAPJ/IobHwbm8YNZNeFntHg\nNXd3d3r0eJgPP/yUnTv3U7t2XSZNGkdISCDHdm3gz8u/yypeQgiHZ1PLWmttAsale/iU1fb1mK9b\n57uMSou+9dYbrPt8I637z8PkUkJWbXJiGY3kt7wXsku2FSpU4NlnJzNx4iT27dvLlNlvcCTyWTxL\nV6ZK/XZU9O1RIDELIUReFflVt77/fgcrVrzLw0Pn4OFV2t7hCDvK6eA1g8FAUFArGnWZQI+J62jY\ndhBxf53h80UjGTDgMY7tC+fundtSzlQI4TCKdLKOjj7AuHEjWbVqHfMn9crTKGPhHCwzA7Jj6TJP\nMbnwQE1/OvWfweHDJ3GvFMRvJyLZtmIke79awMCJr3L9+rVCiFwIITJX5MqNWlo7p05phgwZwFtv\nLaN1a3PpUCl24vwKco68l5cX9Zo/hHuFQO78fZPLvx7g/Ml9BAa+T4sWATzySF969uyV7ZiI3KwC\nV1gcMSYhRM4VqZa1pTUUd/Mvuvd6mJdffpVu3eQ6o/hHTgevWXeZW3/Jszzu4VmaWk27sHvb1xw5\nohk8+GkiIrYTGNiUJ57ow0cfreXatfxtcct8byFEZopcyzopMYH9G+fSLORRnnzyqXu2SatB5IYl\nSad/3/h5eXAjLiFtu5eXF4880pdHHulLfHw827dv4X//28icObNo3rxFaou7N+XKOe4sBLn+LkTR\nVqRa1jMHt+Do1rfxMz7IF6vC7B2OKIYsiXvlyrUcOaIZMmQYP/wQQcuW5hb3ihUr+OvKFZuSY/rq\nbPnF1nrsQgjHUaRa1kuXvklC7GW6D38Dg0FWzhIFI6c9NOlb3Dt2bGX+0tWcORaFh1dpuu1tyyvP\njyYwsCWurq5A5teOs1s9TghRvBWZZL1tWzgrVrxHx6fCcHOXgWQi7/LzsomXlxe9ez/KiVs1OHHu\nOnFXz3Hn8kFmzJjGlSt/0b37wzz8cC+Sk0rg6la4tepnDQlgTFhE2m0hRNHjkMk6fevjzJnTPPfc\nOD744BOCglrZMzQhshT2XDv6Tv8fZSvW4v2FIwE4e/ZXvv32axYuXMDho8epUjeQTdWG07Fj57Ra\n5QWdUGWmhBBFm8Nfs46Li2Xo0IG8+OLLkqhFkZB+rnfNmrWYMOE5vvlmK52GL8W3ouLDDz+gSRPF\n0KEDWLduDb//finHc8SFEMWPQyfrlJQUxo8fTXBwW4YMGWbvcITIk7nrDuBWqgzVGnen2SOziY4+\nRu/ej7Jz5w907NiGTcsncuqn77hw4Xy2+5JpXkIULw7ZDW4ZEbtgwVxiYmJYtWqdnSMSIudyei28\ndOkyPPHEAJ54YgDJyck8NeUtzh2PpHv3Tnh7e9OuXUfatWtPSEi7PC9OI9MahSjaHC5ZW0bF/nFq\nN9t3rWXPzihKlChh77CEyLOsrku7urry2ZJpwDRMJhMnT54gMjKCzz77hClTnqVmzVq0bduedu06\nEBTUWuZNC1HMOFyyBoi9ep6j29/jsWcWYDQa7R2OEPkmJ9ekDQYDqyNuAv58/PEEEhMTiY7+icjI\nCBYtWkD0oUP4GmtTvkZTnn3tIm/OGIibm0P+KWdISp8KkXsO9xc+a0gAAR/MpW7LviyZ+VT2LxDC\nybm7uxMU1IqgoFaEhr7InJU/cmBfFFcvHGH7hoU8+P4LVKjRmPHDHmPPRT/8ylfjjQnB9g47jSRn\nIfLOpmStlHIBlgFNgDvAKK31GavtvYHZQBKwWmu9Mif7DV22m+SkRK5eOEzz9gNsCU0Ih5bThJVV\nN/cro9oy5kYSVeoG8n5oRyYu+IY/zh7i2LGjbPt6MykpySSc6JrWbV6pUuX8Cj9DuU3G0oUvRO7Z\nOhr8UaCE1roNMANYaNmglHIHFgFdgPbAGKXUAzndcfT2DyhbsTa+5avaGJoQRVtW5UEto8Ctp3mV\n8i5DrcYdefPNt+k8+n2Cn5xLYGAQW7ZspkOH1rRq5c+kSeP55JN1nDlzGpPJlOuYcjr6PLvnSelT\nIbBpJoet3eDBwGYArfVepZT1aJn6wC9a6xgApdROoB3wn+x2+pvey/mfdxK9bx9lypS1MTQhihdL\nS3XuugMkp0BJ34rs+svAqlUjSElJ4eTJE+zdG0Vk5Pe88ca/SUhIICiodWrXemsaNWqCu3vuq6pZ\nPnCyaylLS1qIvLM1WfsCsVb3k5VSLlrrlNRtMVbb4gC/7HaYmJhI1NdLafHwVEnUoljLTTUz65ri\nF6/E37fdxcWFBg0a0qBBQ4YPH2V+3sXf2Ls3ir17o1i//mMuXDhP8+YtaNnSnLwDAgLx9vbJUawx\n8XfuqWmePjFnVPNcSp8KkXu2JutYwPqv2ZKowZyorbf5ADey2+GgCbPxLluN0pUbsGD9QcKea2dj\naCIzRmPOPoBF3uTHeS7rWzLDfcXE3wXgywW9AQhdEpm2rWYlX365GENScgpgyDQOo7EB/v4NGDvW\nXA71xo0bREVFsXPnTpYuXUh0dDRKKUJCQggJCaFdu3Zpx7XeZ0z83bREDPD71fi0+5a/YXc317Tt\n7m6uGI0+mf5suSHv5YIn57jgxMTfpfe0r3ZvWtgnx6MubU3Wu4DewAalVCvgiNW2k0BdpVQZIB5z\nF3iW61meP3+ePVs+otWT/wYgMSmZK1fibAxNZMRo9JFzWgjy6zwnpyY9632ZW6nm78STF0Uwa0gA\n0wf6p7VUrW8nJaekPSd7bgQGtiUwsC1TpsCdO3c4fPgQe/bsZuXK1QwaOgKv0pXw9KtIrysHWfbK\nWBb953haLABurgYql/fizCVzh5vlbzh9fFeuxGX4s+WGvJcLnpzjgmP1d9w6N6+zNVl/CXRRSu1K\nvT9cKTUQ8NZar1BKTQXCMQ9gW6W1/iOrnU2YMIEpkyZxIqkKIN1jQuRGftcT9/DwoGXLIFq2DAKm\nMHj259y4fI6EuKsc3rmRxo0XUbFWc7yrBvJAzea4uZckJcWUaRe31DsXxU1BTFe0KVlrrU3AuHQP\nn7La/jXwdU73d/bsWZYvX8uslTI6VAjI+I88J9d7/bw8uBabgJurId++9Hr6GinhVR43VwPvr1/A\n1atX+fbbTcxf+gFHtrxN+WqNqNqgA/HxQTlKzLZ+gFk+ANfM6Zbr18gcb+EoZg0JoP/zq/AuXzMq\nN69ziKIo27Ztw81NSooKkZ3skmHY+Db5usCH9QAxi/LlyzN06HCO3lIk/B3LycM7OXPgKxo2fJuS\nvpUwVqvPBuMFWrQIpGbNWpIohUgVumw3KclJHPz2TW5eOZ+rPwyHSNaVKlXiypU4+aMWwgb2/Lu5\nY/KgWsNOVGvYiRoPlOTX0z9z5bcThId/x9y5r3Dt5i2MVesxsE8nmjZtRtOm/lSoUNGmY8kUMFFU\nZPVePbH3K7z8cl9G2yGStRAiZwo7MVt3vb8f2jHL57q5l2Dt66PueWz8619x9ZLm7t3brFy5nCNH\nDuHhUZJOnToTGBiEn19pgoJaU758+Sz3bd3CD10SyfSB/jmKXxK8KGwZTVe0+DvuGkd//Izvd0Tk\ner+SrIVwMvmd0DPrerccJ6vr6F5+Rrz8jMxKfa7JZOLs2TNs376VJav+w53bccReHkeVKlWITSyF\nl18F/jV9NCEh7fK82t6YsIhMPzSFKEyhy3ZzPSae/ZuWUr1xZ2rXrpvrfdhablQIIQDuKX2aXkz8\nnXtatwaDgVq16jB69DjaPv4CnQe/htbneOutd6jdrAu+5aoQFvY6jRvXZeLEZ9i8+VsSEhKYNSQA\nN1cDbq6Ge2owjJq/g1HzdxT4zyhETlm/Vy1fEE0pKRzasozEO/HUCXrSplK7kqyFEAUiuzrglkTu\n5uZG8+YB1GrckUbB/fjuu+1EROymWTN/3nvvbRo1qsuYMcO4cW4vnm4pGRwp65rk2Y2Mz2ndcyFy\nyvoL7Jkzpwlf+wLxN/8g8NFZuLrbNpVRkrUQIkth49vke9d6RoncuhVeuXIVRo0ay8aN3xIVFU1I\nSHtOR29mw6JB9O3blw0bPiU2NoYUE6RksC5JRqPYhShMJpOJjz5aS69eXXmgTita9XsFN/eSNk+r\nlGvWQog8ySyR21rjPP31ZaPRmDZV7M7tOLrXj2f9+s+YMnUyZao2pnqjh3CpE0AZH0+b4pdBaCK/\nvThAMXXqc6w8d5YeIxZwM7lcnvcpLWshRIHJ7Hp2Rtf1shM2vg1LpnXj6aefZt26z+g+/gMq1Arg\nl/3/Jfy9kez+9n1Onz6Vtn/rGDJjGYQmS3aK/HDnzh3mzn2FVq38qVatOuHhEZSrWDNf9i0tayGE\nXVgn0Zy2wkOX7cbV1cC/n2lNzapGUgxdqN64C7euX+S3Y9vp0r0rDevVZeDAwSTfLYe7h6fUbxD5\nLqPKeMePH2PChDFUr16D8PDvqVmzFvDPezsp2WSuApjNFMjMSMtaCOEQshpVnhHrpO5dtir12z3N\niNmf8txzU9m6NZytK0YR/d1ioqJ2YTJlff06P8uziuIlOTmZJUsW0a9fb8aOncDatZ+kJWoLPy8P\nXAx5q5MvLWshRIEp6FatiyH1/9QbLw9rBcCBv4wk1ejHpRM/MHzMWPw8XRk4cDD9+w+kcmXzgkHZ\nDUKTuuIiM5ZxDr/+eoZnnx2Lh4cHW7b8QLVq1QvsmJKshRB2kT4J5iQpmj8kDfc9nlGLxcOzNLVa\n9KFz78H0bGTgk08+okOH1jRvHsDAgYNJSrSt7Kko3uauO0BiUgoXjoTT4Z1PeGnmTEaNGouLS+Yd\n1fnxhU+6wYUQRcI/071S0gaDrXyhEytf6HTfc60HsL00NJDmzQN44423OHToJI899gQffriGT+YP\n5Nj2d7nx+0l8PWUhoeIir/Pqk5LucmjzW1w4tpUnJi5mzJjxWSbq/CItayFEkZdRyyWj1ranpyf9\n+w+kf/+BXLz4GyOfD+Po1qX8ub8Ub94dyBNPDKBq1WqATOkS9/vhhwi+Wz6VFFcj7Qa+zsLnuxfa\nsXOdrJVSpYCPACMQBzyttb6a7jmLgeDU7SbgUa11bN7DFUIUV/+MGM+fwWBVq1ajSbsBNG77JAOC\n3Pn88/V07tyW2rXr4lklEJeyjSjlU67A6orLNXH7sOVL2OXLl5kz50X279/HvHlhbDvjVwCRZc2W\ntvs44LDWuh3wIfBSBs9pDnTVWnfUWneSRC2EyA9+Xh6U9bV9RG1GDAYDAQEtWbDgTXqOX4tv3Yf5\n6+Jpdn86g23vj+SHjW9z/PixfD2msI+MKudl1S2enJzMqlXL6dChFVWrVicyci/duvUozJDT2NIN\nHgzMT729GZhtvVEp5QLUBVYopSoAq7TWH+QpSiGEyKWctFjTP8fVzZ2q/9eS9W9NZkxYBPE3iAjT\niAAAEZlJREFU/6SxzxkGD+5P+fJGBg4czGOP9aN06TLZ7ltazkXboUPRhIZOwdPTky+//JZ69erb\nNZ4sk7VSaiQwOd3DlwFLSzkOSN8f4AksARal7j9CKXVAa300q2MZjT45jVnYSM5x4ZDzXHDWzOlW\noPt3dTWPNDcafSjrW5Kyvg+y8KUxLFgwj+3bt7N69Wpef/1VevbsyYgRI+jUqVOmg4tu3LqTtq+M\nxMTfzXK7vTlqXHnx1tSO9J2+Ke023Pt7OHnyJK+//jqnTp3i7NmzzJ8/n6FDh2Iw3DsDoaDfhxnJ\nMllrrVcBq6wfU0p9AVh+iz7AzXQv+xtYorVOSH3+DqApkGWyvnIlLudRi1wzGn3kHBcCOc8FryDP\n8fXYBMD8eZScOgfbcix//9YsXdqaGzeu89//bmDq1Oe5ceM6AwYMYsCAQVSvXuPenZnufb01c3es\neQWxyYsiHK4gizO/j/28zCP/r1yJS/s9xF37DdX8Ia7/doTKDbrg4deUyMhPKVOmLFev3iqQOHL7\nZciWa9a7gJ6pt3sAkem2K2CnUspFKeUOhAA/2XAcIYQoNNkt6WlRpkxZRo58hu3bf2Tt2k+4efMG\nXbu25/HHH+GLLz7n9u3bAJmuCCYKV3ZTtWKvnGP3ZzMpX6kW+/Ydxr/T0zRo3ZcyZcoWYpTZsyVZ\nvws0VEr9CIwCXgFQSk1RSvXWWp/APPAsCogA1qQ+JoQQRUZOlgZt3Lgp8+aFcejQSYYOHcbnn6+n\nWbN6PNRnEDf/PI3JlHHit2UhE5G/kpOTqXg3mj1fzKFxxxFsXLMIb2/H7frP9QAzrfVtoH8Gj79p\ndXsR5mvWQghRJORmSc/0SpYsSZ8+j9Gnz2NcunSRQc/NI/qbN3B18+C3pl0Y070GRqPxntfkpE60\nDFLLf2Hj2xAdfYCuXTvg4+ND96fnU6bCg2nbHXV+vVQwE0KIVLldTCQjVapUpcujI+g44l0adhrN\n3ZgLtG7dnMGD+7Np00bu3r2bT9GKnIiJv0NM/B1u3Ypj9eoV9O7djcGD+zNu3ES+/PKbexJ1Ti+F\n2IMkayGEyGezhgRgMLhQvlpjIr/7gkOHfqZXrz6sWvU+zZrV5+WXZ3Lzr/P2DtPpWZKv3vsl9RvW\nY+fOSEpW70iHoW/Rr9+T943ydmRSblQIIVLlZ3ezi1Ue8Pb2SRs1/uuvZ1i//iP2fDqHqlWr8bHf\n0/Tp0zfD66WO2iVbVKSkJHNqz2dcOLKFp6YuZ/6kXvcNNrP+neflUkhBk5a1EEIUgDI+JSnjU/K+\nx2vVqs2sWXM4ePBnJk9+nvDw7/D3b8jkyRPYv39v2trbjtwlWxT88cfvRG98leu/HaX94DDmT+oF\n/NMtnpn8uBRSEKRlLYQQBSB9Kz39YDE3Nze6detBt249uHz5Mp999gnPPjsWNzc3nnpqKH9TD3DN\nch8iY6P/33rC185g0sTx1PJog4uL+Txar2FeUDXfC4q0rIUQws4qVKjAc89NISoqmjfeWMyJE8f5\nbNEIjm1/j5S7t4pUUrG3X345zQ8b5tGk3UCmTXshLVHnVE6m7NmDJGshhHAQBoOBVq3asHTpe0RH\nH8PNJYWtK55h4sRn7ukiF/czmUysX/8RvXt3pUaznlRq2BW4N/kW5fnt0g0uhBCFILeDxXx9/Qju\nM5UWnUdSz/M0Eyc+w7VbKTxQoxn7AlwJDGzJ9HejAOkWj4uLJTR0CsePH6X78AXEpJQjOSXjrm5H\nvB6dE5KshRCigNl6rfSfJNyDG96t2Ld3F9cu/szg4cOpVN4XnwfbU6ZCTRISmlOy5P2D2dJzxmve\nhw5FM2bMcNq27UB4+Pe8+cXPxFxyvlWZpRtcCCGKAIOLC+WrN0W1GcjQF9Yyd+4CLp09yu5NS2nW\nrB6zZ8/g1Clt7zALTUpKCsuWLeWpp/rx0kv/j4ULF+Pp6Vmku7qzIi1rIYQoYPkxf9d6Hy8NDWTu\nugP493wegPIlb1Hq70M89lgvatWqzZAhw+jd+9H7WtvXUlcVS6+otbivXLnCs88+Q2xsLJs3R9y3\n4llWXd1F5WdMT1rWQghRCPJj/m5m+/ArV5mZM1/m4MGfGTNmPBs2fJrW2tb6JMA9c7Udfd52Vitl\nRUZ+z0MPhdC4cVO++uq7+5cmdVLSshZCiCIoo9a6u7s7vXo9Qq9ej3D+/Dk+/vhDHn+8NzVr1qJ0\n7Q64lGuGq/v9yb4wKqVlVTksJxITEwkLe51PP/2YpUvfo337jvkZnsMzOMhUAJOzLnTuKJx5MXlH\nIue54Mk5/kdOuq8TExMJD/+Odes+YGfUPqrWb8+/pg2hYsVK5iU+P/qJM6kDsmpX8WXWkACMRh+G\nvRKe7b5tidUiq/2OCYsgKdlEOd+ShI1vw4UL5xk7diS+vr4sXbr8vhXMiiKj0SdXhcltblkrpfoC\n/bTWgzLYNhoYAyQBr2mtv7H1OEIIIWxnaW0fvlEZQ63TXDi6lcmhL1LSNRFPT0/K12lLqcpBlPIt\n2ARoab1ndynAeuR8TPwdNm3ayAsvTGXixCmMHTsBF5fiefXWpmStlFoMdAUOZrCtIvAs0AIoBexU\nSm3VWsu6cEIIYUeefhWoFzIYN9ch+HqW4IkAVzZs+JT1H03Fz1iD4PFP8+uvpSlfvlmu951VK986\nAV+LTcDNNfNG5Z2EeG78fpL4mMtcPRfNoQ0X+PjjDfj7t8h1TM7E1pb1LuBL4JkMtrUEdmmtE4FE\npdQvQBPAsUc0CCFEEZObLupZQwIY8e8dwD+t26CgVgQFtSKl+qNcPL2fI0cOs3jxIjw9S2HyVdRs\n8Sihy/J/BHVSsokxYRG8H3rvdecjRw4RvnIS5y/8hodnaR6o6c/2r3/Ex8c3X49fFGWZrJVSI4HJ\n6R4eprX+XCnVIZOX+QAxVvfjAD+bIxRCCJEvyvlmXDjF1a0ENeoHEzY+FJPJxPRFX7Bz+yZ+/Ph5\n/LtP5Ny5Sri7u+Pm5o67uxtubm6pt91xc3PDYDBkOUjNMhju5tVLePpVwMXFNa2lDeZSoatWLWfR\nogXMmxfGl4ddcXH3wuDiypKNp5xqvrStskzWWutVwKpc7jMWc8K28AFuZPcio/H+tVxF/pJzXDjk\nPBc8Oce2cU3tfr4ea06slvNoedxy31i5NvXbPU2ZyvX4OXId/Q+tJzExkcTERJKSktJuJyYmkpyc\nbF4sw8WVBu2GMW+dJxMeeZALFy5w/vz5tH+RW/cRe+UcJUp5U7FuG+7E36DTlpe4eDkWv/LVKO+Z\nwJ49e6hduzbfTt9EUnIKAO5urvL7Jg+jwVNb1s9orQeme7wCsBUIBEoCe4Cm2VyzltHgBUxG0BYO\nOc8FT86x7UKX7SYm/k5aq9Yy+js9o9GH3tO+Sruf2fPAXEnstbV7OXjoMLs/m4mLwUCVKpWpWrUa\nVapUpVq1alStWp3Pd12jTKX/4/rFY8Rc/Y0Snn40bNiI3eEf4+lTji3/WYGHh0danJZr2+m7yp1F\noY0GB0yp/wBQSk0BftFab1JKLQF+xFx0ZaYMLhNCCPsLG9+GuesOpE3Vyoqbq+GerurMuLi48PLw\n1owJS6DH+A9YNasXrq73L0u546L5enmlukEYa7UE4IEqvnToPwsgLVFb4sysKEpxZXOy1lr/APxg\ndf9Nq9srgZV5C00IIUR+y2npUz8vj7TWbU6uGft5eYCXMcNEbV0xzdXVhaTk5LT7hVGQxRkUzwlr\nQghRjOWk9GnY+DaU8y2ZoxKpoct2cyMu47rjWbl4JZ6kZBNJyab7SqBar0MtJFkLIYQoYNYrYb07\ntX3a41WNXnaMqmiRZC2EECLPyviUzLIlbN2adzGY/znrcpYFQRbyEEKIYian3cs5fV5urzuX8fln\nvndeVyIrLiRZCyGEsJl1KdG56w7kqIUs16JzT5K1EEKIApdZgpbEnTNyzVoIIYTN5Lpz4ZCWtRBC\niDyR684FT1rWQgghhIOTlrUQQog8kevOBU9a1kIIIYSDk2QthBBCODhJ1kIIIYSDk2QthBBCODhJ\n1kIIIYSDs3k0uFKqL9BPaz0og22LgWAgDjABj2qts1/tXAghhBD3sSlZpybjrsDBTJ7SHOiqtb5u\na2BCCCGEMLO1G3wXMA4wpN+glHIB6gIrlFI7lVLD8xCfEEIIUexl2bJWSo0EJqd7eJjW+nOlVIdM\nXuYJLAEWpe4/Qil1QGt9NK/BCiGEEMVRlslaa70KWJXLff4NLNFaJwAopXYATYGskrXBaPTJ5WFE\nbsk5LhxynguenOOCJ+fYsRTEaHAF7FRKuSil3IEQ4KcCOI4QQghRLOSlNrgp9R8ASqkpwC9a601K\nqQ+BKCARWKO1PpG3MIUQQojiy2AymbJ/lhBCCCHsRoqiCCGEEA5OkrUQQgjh4CRZCyGEEA5OkrUQ\nQgjh4PIyGjzPUqudLQOaAHeAUVrrM/aMydmkTp9bDdQAPIDXtNab7BuVc1JKPYB5muJDWutT9o7H\nGSmlXgR6A+7A21rrtXYOyamkfiavBP4PSAFGa621faNyHkqpIODfWuuOSqk6wBrM5/kYMEFrnemI\nb3u3rB8FSmit2wAzgIV2jscZDQKuaK3bAd2Bt+0cj1NK/VK0HIi3dyzOKrVqYuvUz4sOQC27BuSc\nugJeWusQ4FVgrp3jcRpKqenACsyNJjBX+ZyZ+tlsAPpk9Xp7J+tgYDOA1novEGDfcJzSBuDl1Nsu\nQJIdY3FmYcC7wB/2DsSJdQWOKqU2ApuA/9k5Hmd0G/BTShkAP+CuneNxJr8Aj/HPmhrNtdaRqbe/\nAzpn9WJ7J2tfwHrpzOTUbhiRT7TW8VrrW0opH8yJe5a9Y3I2SqlhmHsvtqQ+dN8CNyJfGIEWQD9g\nLPCxfcNxSruAksBJzD1FS+0bjvPQWv+XextL1p8TtzB/OcqUvRNjLGBdgNZFa51ir2CclVKqGrAD\n+FBr/am943FCw4EuSqkIoBmwVilVwc4xOaOrwBatdVLqmIAEpVR5ewflZKYDu7TWin/eyyXsHJOz\nss51PsDNrJ5s72S9C+gJoJRqBRyxbzjOJzVpbAGma63X2Dkcp6S1bq+17qC17ggcAoZqrS/bOy4n\ntBPzuAuUUpUBL+CaXSNyPl7809t5A/NAPlf7hePUDiql2qfe7gFEZvVku44GB77E3CLZlXpf1r7O\nfzMxd6+8rJSyXLvuYVkVTYiiQmv9jVKqnVJqH+aGxvisRs8Km4QBHyilfsScqF/UWt+2c0zOxvKe\nnQasSO25+Bn4T1YvktrgQgghhIOzdze4EEIIIbIhyVoIIYRwcJKshRBCCAcnyVoIIYRwcJKshRBC\nCAcnyVoIIYRwcJKshXBSSqkHlVJn7R2HECLvJFkLIYQQDk6KogjhBFKXj3wp9W5VYB/wGuYShhFA\nI8zlIx/VWl9XSk0EBmMuL5kCPKm1PqmUegPz6j/JwFda61eVUt7AO0BDzKUn50uNeSEKl7SshXAe\nrYBngPqYV07qhXmlqoVa68bAZWBA6gpsfYD2qY9vBMYrpaoD3bXWzYA2QB2llAfmLwEHtNYBQHtg\nllKqZiH/bEIUa5KshXAe27TWZ1LrZa8DOgGXtNYHUrcfB8prreOAp4CnlFKvA70xt7AvAbeVUjuB\nKcBsrfUdzC3tsUqpg8APgCfQoDB/MCGKO0nWQjgP67VyXYHEdI+ZAEPqkql7MK8n/w2wBvPytMlA\nEDAbKAdEKaXqYv6cGKS19tda+wPBmFdyE0IUEknWQjiPjkqpikopF2AI8B33LnBvEQCc1lovBvZj\nXqbWVSnVBHPLOVJrHYp5JSCFeS308QBKqUrAQczXxYUQhUSStRDO4xLwMebu7ovANv5Zjs/ChLlV\n7KKUOgaEY07QD2qtjwBRwDGl1E/AWeBb4BWglFLqKLAd89roMiVMiEIko8GFcAKpo8Ff0Fr3sHcs\nQoj8Jy1rIZyDiftb0UIIJyEtayGEEMLBSctaCCGEcHCSrIUQQggHJ8laCCGEcHCSrIUQQggHJ8la\nCCGEcHD/H5z6sN651EjuAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 6 }, { "cell_type": "markdown", "metadata": {}, "source": [ " " ] } ], "metadata": {} } ] }