Skip to content

Instantly share code, notes, and snippets.

@alexlouden
Created March 21, 2018 14:35
Show Gist options
  • Save alexlouden/e42f1d96982f7f005e62ebb737dcd987 to your computer and use it in GitHub Desktop.
Save alexlouden/e42f1d96982f7f005e62ebb737dcd987 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"# Some random data\n",
"N = 1000\n",
"t = np.linspace(0, 10, N)\n",
"np.random.seed(0)\n",
"x = np.random.randn(N).cumsum()"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD8CAYAAAB0IB+mAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzt3XeYlNXZwOHfmZ3tvVK2sJRFpAouiKIiQRQr9tg1jXzGxCSaGEtMsUSTGI1GscRo1BiNsccuNkRAXEB6W3YpC8s2ts/2Od8f78zsDDvbmLozz31dXszbzxDyvGdOeY7SWiOEECL0mQJdACGEEP4hAV8IIcKEBHwhhAgTEvCFECJMSMAXQogwIQFfCCHChAR8IYQIExLwhRAiTEjAF0KIMGEOdAGcZWRk6Pz8/EAXQwghhpQ1a9ZUa60z+zsvqAJ+fn4+RUVFgS6GEEIMKUqpPQM5T5p0hBAiTEjAF0KIMCEBXwghwoQEfCGECBMS8IUQIkxIwBdCiDDhccBXSuUqpT5VSm1RSm1WSv3Utj9NKfWRUmqn7c9Uz4srhBDiSHmjht8J3KS1ngjMBq5XSk0EbgE+1loXAB/btoXNvkMWPt5aQWVDa6CLIoQIEx5PvNJalwPlts+NSqmtQDawCDjFdtqzwGfArzx9Xig4UNfCSX/61LH91o/nMDUnJYAlEkKEA6+24Sul8oHpwFfAMNvLAOAgMMybzxrK1u2tc9n+quRQj3Pe2VBOcWWjv4okhAgDXgv4SqkE4FXgZ1rrBudjWmsN6F6uW6yUKlJKFVVVVXmrOEGtprnNZXvp1gp2VnQH984uK9f/ey2nPrDM30UTQoQwr+TSUUpFYgT7F7TWr9l2VyilRmity5VSI4BKd9dqrZ8EngQoLCx0+1IINTVN7S7bX5UeYsGDy5gzLh2zycRvzpnoONba0UVMZIS/iyiECEHeGKWjgH8AW7XWDzgdegu4xvb5GuBNT581FFmtmnvf28qKXdXUWdr5x/JS6ls6iIuKIDHa9X37ZXENn++o4sLHVjj23fzKBn8XWQgRopTR2uLBDZQ6EfgC2AhYbbtvw2jHfxnIA/YAl2itezZWOyksLNShli1z7d5aLliywmVfVmI0JqX47JensKqkhmuf+brPe+y+7yxfFlEIMcQppdZorQv7O88bo3SWA6qXw/M9vf9Q93Vpz3dcZWMb47ISiImM4OSCTEalx7GnxuJyzrisBIorm8hIiPJXUYUQIU5m2vrY9oPuR9okxhjvWpNJccdZE3scP3faSG5aMJ7qpnYaWzt8WkYhRHiQgO9D1z6zmtfW7WdaTjIXzsgBIMJk/BhKj492nJefEef4nBoXCcCxo1IpzE8D4NdvbBrwM7/efYiNZfUel10IEXok4PvI08tL+Wy7Mcw0OS6Key+Ywurb59NlNfpMpuUkO84dm5ng+BwZYfxPMjw5hpn5qYxMjuHDzRWO6/pz8eMrOeeR5d76GkKIECIB3wdaO7q48+0tju39tRaizCayEmMc+y6Zmev4rJRiRp4x0/aaE/IBGJEcgznCxHWnjKWlo6vH2H2Akqom8m95h69KavhkWwUfbD7oOCbNQEKIwwXVmrahYldVk8t2Qkyk4/MfL5zC/rpWhiXFuJzz4uLZtLZbSYo188OTx2C21fQzbS+JyoY2lxcGwLsbjYnMH26p4B/LS12OPf75Ln55+gTvfCEhREiQgO8Du6u7R9w8cMk0Thib4dj+9sw8t9dEmyOINhsTrMwR3YOespKMtv6KhlYmZye7XFNW2wJAfFTPiVmPfrqLS2fmkZsW1+OYECI8SZOOl2mtuf7fawFY9st5XDAjh+HJMf1c1Tt7+/7W8oYex3ZWGr8kHv6k2O21Nc3tbvcLIcJTWAX8prZOOrqs/Z/ogUNOQTYv3fPadXJsJNkpseyqanbZf+97W1mzp9btNTctGA/AeY9+KSN2hBAOYRPwu6yayb/9gJ//5xufPqek2gjMz1w702v3TI6N7NEJ+8TnJT3OO6kgg+yUWM49ZqRj3+vr9nutHEKIoS1sAr498L29oZxfv7GRTft9U/P9b9E+wJgp6y2JMWYaWjpd9iXY8vD87bLpZCRE8dgVM3juu7P44uZ55KXFcdd5k4mPiuDDLQfd3VIIEYbCptPWOcD/a9Ve/rVqr09y1LxcVAbAyJRYr92zprmd4somKhtbyUqMQWtNa0cX150ylnOmjeScad01emXr771q9igONbXz4NIdtHdaiTKHzbtdCNGLkI0ClvZOl3bu6qae49gt7Z099nnCapscFRsZ4ZhR6w3Fts7ZpVuMDNNNbZ10WjXp8X3n2bGP8Hl2xW6O+vV7tHZ0ea1MQoihJ2QD/pJPd/HE5yVc+NgKfvfWZj7Z1jMd/4E6764n+8f3twFwo63T1FsuKTTSMthr7w2txovKno+nN5kJRsC/592ttHVaqWzo+dITQoSPkA34zkMS/7liN5b2nrXbCi8vIP7EMqMjtcXLNem7z5sC4AjYDS1GB26S04Qud8YPS3TZ/sV/1+NpOmwhxNAVsgG/qa1nc01GQrTLdnm9dwO+nT0BmrdEmU2kxUfxj+UlvLX+QHfAj+37Oblprv0Iq3cf4qCXX3JCiKEjZAN+hZtgPmt0qst2aXVTj3M8kZ0SS0ZCFJcfN8qr9wVj0ZSG1k5ueHEddQOs4SulmDU6zWXfl8U1Xi+bEGJoCMmAX9vcztq9tcQethbscaPTuXPRJN68fg5gpB/4aEuF955raef86dle7bC1S3PqoLV3RA9kYtffLpvusn3X21vIv+Ud3pDx+UKEnZAM+E8sK6HTqrn8uO68NT+eN45vz8zl6uPzmZab4tj/92U9JzAdidaOLiztXaTE+WaFKuc8Ok8uK2FYUjTJ/TTpgPHL4MfzxvHsd2fx9LWF1Nt+Hdg7mIUQ4SMkx+Ev21HFrPw0jnEK7IvnjiEmsmeSsQ3767zyzFqL0Umc1s9QySP1g5PG8KTTy+nwXy+9UUrxi9OPArqHjYLxgiqtbuZQcxvHjkrr7XIhRAgJuRq+1priqiam5SY7xtkvPnlMr+3drR1WrzTr1DYbNWdvd9jaZSZG8/FNcx3bc8Zl9HG2eyanpiZzhIl593/GhY+tHPTIncqGVjYfkBw9Qgw1IVfDr2/poL3TyojkWBYdk82eGgvXzxvX5zWLny+i9F7PZt3aa/ipPmrSASNz5t+vLuT1dWXcfd5kj+5V1djm8jkraWAZPf/z9V5+9epGwHi5rfvNaR6VQwjhPyFXw7cHsszEaGIiI7h54QTio3u+1/79/eOwV3gzDxuuORgvF+1j5a4aR5bMVB816dgtmDiMJVcci1Le6xguq2vp9xytNeX1LY5gD1Br6ZCVtYQYQkIq4O+saOTT7caM2sPH3B/uhHEZ7LznTAAmjkw64mfe/MoGLvv7Kur8UMP3hjevn8Ndiya57Ktu7H8G7jNf7ub4ez/psd/bs5WFEL4TMk06ZbUWFjy4DDAmKk3PS+nnCogwKRZMHMa6vXVHlGDMue37jjc3A5DiozZ8b5mWm8K03BSuOj6fTfvrOftvy6lu6n+hlFUl7sfvl1Q1kRBjJtuLyeKEEL4RMjX8Okt308Ifzp/idkSOO2dOGU51UxvPfFna/8mHae3ouZhKZMTQ+Su1p15wl1jOWXunlQ+dOrZPKsjg/+aOBeC6F9Yy575PJGWDEEPA0IlO/Whs7U6lcNGxOQO+riDLCHr3vjf4cen2jtqhKspsQil4ftWePs97cfVex+crZ+fx/PeO46fzC1zOqbVIW74QwS5kAr59QtHhM0v7c9Tw7gRjLxftY09Ncx9nu9pR0TioZwUjrY2O7vL63jtu2zu7f8nYE7nFHrZwemn1wP/ehBCBETIB355QzHmy1UBERph48NvTAKMD9oYX1w342qe+cG0GeueGEwf17GDSVyK5upb+f8nsloAvRNALnYBvGx6YfASdpiOSuzsc99e1Drg9uqzWwolOE6AmjUzu4+zg9NgVMwA40MfQzN01FgBeWjzbZf9z353Fiz8w9j308U7aOmWBFSGCWcgE/PqWDpSChKjBDzyakp3M6ZOGMSs/jeqmNsdC5P2ptXQwJjOezb8/nbV3LBj0c4PBKUdlERmh2LS/we3x4som3tlQzsJJw5k9Jt3l2MnjMzl+bDonFWSw95CFm1/ZQJdVOm+FCFYhFfCTYiJd0gcMVHy0mSeuKuTO84zx6ev39Z9fp7PLSn1LB6lxUcRHm32WQ8fXYqMimJqTwuc7qtweP+/RLwFjIltvHr50OqcencWb3xzg7Q0HfFJOIYTnQirgDyR7ZF8KshLJSIhi6db+c+v8zxbY0hOGZqB3Nmt0GlvLG1hdeshlf5dVOxaS6eqjmSs1PoonriokKsLExjIjx87cP3/KhY+t8F2hhRCDJgHfSYRJcfzYDNbv6z8xmP2cC2YMfAhosLpgejYA++ssLvudRyzNyHNdPOZwESbFmMx4nlpeSmeXlT01FtbsqfXqegNCCM+ETMA/UNfC8OSBJQDrS0FWAvvrWmhxswaundaaZTurOH5MOglu8vQMNem2NBS/fn0T+be8w1JbkP7P1/sAuGvRJC6ckd3vfezLJ55mm/EMuF08XggRGCER8K1WzZ4aC/kDWAGqP8NtWSP7mn1aZ+mgpKqZeRMyPX5eMEiKMV5azbaX3PefK2LFrmrHouwzRqUOKFnbI5cZI36cO70HmrdfCOF7IRHw1+ytpa3T6rIq1JGyZ7vsaxbtbltTR356vMfPCwZmN+kgLv/7V47P/SWiszuxIIOcVNecOnVDfDayEKHE5wFfKbVQKbVdKVWslLrFV885qSCD+UcP8/g+afFGP0Bfbc8fbalAKTh6xJFn2Qw2Ny4Yz+NXzuixf/aYNIYNMFc+wFG2/DzT81KYMDzR0ekrhAg8nwZ8pVQE8ChwBjARuEwpNdHbz5mZn8bz3zvOK+3pWYlGcPvbJ8VUNbbx92UljolYVqvmoaU7efObA0wYnkRumudNSMHihvkFLJw8osf+Ry7v+RLoy3hbqooJwxOJjzbT3N7J5gP1bNovK2QJEWi+ruHPAoq11iVa63bgJWCRj5/pkdy0OMfs2Rtf/oZ73t3KlnJjUlJxVRMPLt3B/roWtpa7n6gUagY78mmCLeC3dViJjzbT1NbFWQ8v5+y/LfdF8YQQg+DrgJ8N7HPaLrPtC2oLJhpNQ1vLjeRoL3xlZIt0XigkL4Rq930ZbLrn0yYO54Lp2Vx2XB4t7Z0uk9j6St8ghPC9gHfaKqUWK6WKlFJFVVXuZ3v6m33WrH2kzr9tAb/KaeTOaz86wf8F84O7zpvMHy+ccsTXx0ZF8MC3j2Fmfhpf7651Obb5QHj8KhIiWPl6EPl+INdpO8e2z0Fr/STwJEBhYWFQJGJxV6tt77RSY1sZ6srZeQMeuTLUXDV7FGAkgusrg+ZAzB2f6ZKyYe8hSx9nCyF8zdc1/K+BAqXUaKVUFHAp8JaPn+mxU4/Ocny+9oR8AHZWNlLd1EaESXHnuZMDVDL/mZyd7GjaOlIPfvsYl+3yuhZaO7pkdSwhAsSnAV9r3Qn8GPgA2Aq8rLXe7MtneoM5wsR7Pz2J+y6YwlXHGzXezQcaqG5qIz0+6ogStIWjwzt8N5TVM+GO91ny2a4AlUiI8ObzvABa63eBd339HG87ekQSR49IwmrVxEZG8JcPt1PR0MbYzNCYbOUPESZFalwktZYO5k/I4mNbmoX7P9zOdXPHyotTCD8LeKdtsDOZFFOyk6loMDps7X+KgVl7xwJW3votbjxtvGOf1vDHDwa/hrAQwjMS8AfgV2dMcHyOMstf2WAopRiRHMukkclccVyeY//zK/teOF0I4X0SvQZgYgilUAike86fwtQcI9+RpY9spEII35CAPwCxURGcMXk4AGdOGR7g0gxt03K6F5mvbZbEakL409BP5u4nj115LAfrW0NihatA+vXZRzM8OYY/f7Cd9WV1KKV4aOkObl44oceauUII75Ia/iAMT44ZdKoB4SraHMG500YCcLC+lZtfWc/avXXc8cYmGZ8vhI9JDV/4nX3Ngc0HGqhoaGNydhKb9jew+UADq0pqiI82c9msvH7uIoQYLKmuCr+Lj4rAbFI8v8oYqXP17HwAtpQ3cPc7W7n1tY0BLJ0QoUsCvvA7pRQpcd19IfNtqSxufmWDY9+Hmw/6vVxChDoJ+CIg7DOWrz5+lGMRdWeLn1/j7yIJEfIk4IuAGJ1hBPzsFGMN3Mtm5fZ1ep/aOrtY/FwR6/bW9n+yEGFMAr4IiJ/ML+D4MemcP91YD+fnC8ZzwYxslv9qnuOcgY7TX7O7lg+3VHDPO1t9UlYhQoUEfBEQ2SmxvLh4Nlm2BdKzEmN44JJjyEntXkls2c6BLYhTWtMMGGkvinYf8n5hhQgREvBF0PnbZdMBMKmBZdOstCW0W7GrhoseX+mzcgkx1EnAF0FnZn4aAI2tnQM6v7LRNYNpZ5fV62USIhRIwBdBJzHGmA/Y2NoxoPOrGl2XYmyWxGxCuCUBXwSduKgIIkxqwDX8qsY2MhOjGZZkDO9sbhvYdUKEGwn4IugopUiLj6KycWCLqFc2tjF3fCa3nzUR6Dvg/+qVDSx69MsBl2VnRSPFlY0DPl+IYCYBXwSlnNRYXi4q6zehmtWqHTX8hOgIoLtJp6mtk8oG15fGf4r2sX5f3YB/BVz99GpOfWCZ9AuIkCABXwSluCgjeBdXNvV5Xq2lnU6rJisxmsQYY9H0l4v2AXD2w19w2l+XsaGsjt//bzNWa/fLo7ppYEtVltcbL4yt5Y18LUM+xRAnAV8EpevmjgOM9nm73765iTv/t8XlvCpb4M5KjOFo28pk//5qL5WNreyusVBn6eDcR77kmS93U1bb4riuoWVw7fy3vr6Bix9fyQeS40cMYRLwRVDKsnXA1jjNtn125R6e/rLUUYMHWF1q1LqHJ0eTEG3mVFsitln3fNzjnruqu38tNAxgBJBzc9Km/Q0AvLqmbDBfQ4igIgFfBKU0W878Tfvr+d/6Ay5t7je/soHPd1Shteal1fs4ekQSx+SmAnD9vHG93vODTd218y+Lq1mxq5rX1/UewNs6e7bbf737kCzUIoYsWQBFBKXUuCgiIxRPLCsB4E8XTnU5fs3Tq7nn/MlsPdjAT+cXEGEyZuVmJnZn3owymzht4jDqLB2sL6vjpa+7fxks+WwXSz7bBcD503PclqHpsI5dpaDW0kGdpcOxiIsQQ4kEfBGUIkyKjIRoR6fpMyt29zjnxdV70RpH2z24Bvz/LJ7N9Dyj5v+TF9fxv/UHmD8hiyizifecavv1LR0kxxodvj94rojUuEj+dNE0GlqMZp87F01iRHIslvZOfvrSN9z33jY6rZr7L57K0q2VzB2fSZRZfiyL4Cf/SkXQsjo1nWwtN9rQR6V3J1ezt6tnOQX5aHMEX902n533nOEI9gDTcpIBY13iP5w/xeU5JVXdbfsfbang5aIy/rG8lIO2IZ3jMhNYMHGY42Xyn6J9vLq2jHX76vjBc0Xc9rqs0CWGBgn4ImhFuEme9viVx3LUsESOGpbo2Jd2WPPKsKSei82PtOXdnzgyiZS4SBKju3/c7qoysm06d+Te9fYWLv/7VwCMsF3r/GIBKK8zXgif7xhYVk8hAk0CvghaJlu7/LisBMDIsTNheCIf/PxkUuIiHecNpD39jMnDeeH7x3H5rDyUUkRHdv/Tr7DV5Ettgf9wI5KNFM4Zh63Mtb3CmIFbbxlYzh8hAk0Cvgha9iaZGXkpgNFco2y1fnNEd+3fubbeG6UUc8ZlOK5fcsWxjmP2nD27qtxP8oqJNCaBJcdGEun0XPv5HVaZhSuGBgn4Imj98cIpvPzD47lqdj4ANc3dk7Duu2AqN8wvYOc9ZziC+GDMGp3GxzfNJdpsorG1A60172wo7/PloZRyqeV/uq0SAK3hd29t5uGPdw66HEL4k4zSEUErLsrMrNFptHYYuXHmjM1wHMtNi+PGBeM9uv/YzARGJMfQ0NrJm98c4ONtlcQ4NfX86cKpjlTNdpmJ3SOHLE5pmP9pG0V0w/wCj8okhC9JwBdBLyYygrd/ciJ5TiN0vGV3jYXdNRZHU01rR3fzzMWFOT1+PRzejj85O8kxWkiIYCcBXwwJk7OTfXJfkwKr7m7HP2vKCH78rXFsOdDgtqko87CAf1JBpkvAb+vsItoc4ZOyCuEpacMXYe2Ry2cAxlh8k4K/XDKNo0ckceGx7mffXlyYw80Lj3JsXz4rz+X4nf/bQlunrLglgpPU8EVYs8+w3VXVzJxx6Y4ROb0pzE+jMD+NSSOTaWnvJDfNtZnpha/2khoXxS9OP6qXOwgROFLDF2HNHvABRiTHDvi6ueMzWTh5BAAvfP84l2MfbpEUyiI4ScAXYc2ehhlc8/AMxpxxGZTee6Zj23QEw0SF8AePAr5S6s9KqW1KqQ1KqdeVUilOx25VShUrpbYrpU73vKhCeF9WYowj0KfFHXkGTOcO3prmdkmhLIKSpzX8j4DJWuupwA7gVgCl1ETgUmASsBBYopSSoQsiKJ1ckAngSLF8pOw5faoa2xxpF5btqOqxrq4QgeJRwNdaf6i1ticNXwXYhzYsAl7SWrdprUuBYmCWJ88Swlem21I3eJriePVt81l923wSY8ws/OsX5N/yDlc/vZqrn16N1ap5f1O5y7q6QvibN9vwvwu8Z/ucDexzOlZm29eDUmqxUqpIKVVUVSVZB4X/XT4rjwe/PY1LZ+Z6dB9zhImspBjuOHuiy/6S6mYK71nK//1rLf9ds6+Xq4XwvX4DvlJqqVJqk5v/FjmdczvQCbww2AJorZ/UWhdqrQszMzMHe7kQHjOZFOdPz8Ec4Z36zyWFuUwY3p2+ub3TyiHb2ryHmiWzpgicfsfha61P7eu4Uupa4Gxgvu7uqdoPOFeXcmz7hAgLozPi2Xawscf++GjpyhKB4+konYXAzcC5WmuL06G3gEuVUtFKqdFAAbDak2cJMZR878TRbve3u1kYXQh/8fQ37CNAIvCRUuobpdTjAFrrzcDLwBbgfeB6rbXMNxdhozA/zWVs/kc/PxmAQ83t7DtkobNLAr/wP49SK2itx/Vx7B7gHk/uL8RQppTivgumsK/WQoFtScYln+1iyWe7+PNFUynMT6OptZMpOb5JDCfE4SSXjhA+dKlTcrUZeSms3VsHQFltC7985TMAnrq6kFMnDgtE8USYkdQKQvjJwsnDHZ/317U4Pn//uSK2HZSc+sL3JOAL4SfDkmIcn19ZU+ZybNWuGn8XR4QhCfhC+MmYjIRejz362S5KellEXQhvkYAvhJ+My3IN+IWjUh2fqxrbOPPhL/xdJBFmJOAL4SexURFs/v3pPHzZdCaNTOKpawr54uZ5zBqdBriupyuEL8goHSH8KD7azLnTRnLutJEApMRFkZMay+rSABdMhAWp4QsRYL8/dxLQs8lHCG+TgC9EgCXGRHLB9GwsbZ39nyyEByTgCxEECoYlcqC+lTpLe6CLIkKYBHwhgoA9nfKNL68PcElEKJOAL0QQyE6NBeCTbZV0SGI14SMS8IUIAtkpsY7PBbe/R2uHJJcV3icBX4ggEB9t5q/fPsaxvbr0UABLI0KVBHwhgsR507O567zJANRK563wAQn4QgSRU4/OAsDSLk06wvsk4AsRROKijMnvzTImX/iABHwhgkhclLHI+StrymS0jvA6CfhCBJHICOP/ktsONvKXD3egtQ5wiUQokYAvRJB6/PNd/P2LkkAXQ4QQCfhCBLG/Lt0ptXzhNRLwhQgyX97yLS6zLX5uae/iYENrgEskQoUEfCGCTHZKLPdeMIV/fe84AEqrmgNcIhEqJOALEaQmjEjEbFK8s7E80EURIUICvhBBKiMhmqk5yeyukRq+8A4J+EIEsdS4KL4srqG9U8bkC89JwBciiH28rRKAV9eWBbgkIhRIwBciiN1tS6Z262sbA1wSEQok4AsRxK6cPSrQRRAhRAK+EEHuhm+NA8DSLgnVhGck4AsR5PLS4wGobGgLcEnEUCcBX4ggNzojDoCt5Q0DOn9jWT3f7KvzZZHEECUBX4ggN2lkMnFREfx79d4BnX/OI8s579Ev2XfIwg0vrqNJcusLGwn4QgS5mMgI5h89bNATsH74/BreWn+ADWVS2xcGCfhCDAFZidEcaup/ndsua3dmzdJq4wVxy6sb+dELayTrppCAL8RQkBwbSXN7F3tqmmnt6CL/lnd4fuXuHudVNXZ37LZ0GOvi7j1k4d2NB6lploXRw51XAr5S6iallFZKZdi2lVLqYaVUsVJqg1JqhjeeI0S4GpNpjNR55JNiDtYb6ZL/8O62HuftrGzs9R4VkmY57Hkc8JVSucBpgHOP0hlAge2/xcBjnj5HiHB21pQRTByRxJ4aiyNwt3R0UXlYEH9o6c5e72F/UYjw5Y0a/oPAzYBzA+Ei4DltWAWkKKVGeOFZQoQlpRSjM+Opbm6jpLq78/aJZSX8+o2NfLKtAoDmdqMZ57fnTOTxK2cwJjOe604ZS7TZxPLi6oCUXQQPsycXK6UWAfu11uuVUs6HsoF9Tttltn2S2FuII5QRH8W+QxbuenuLY98/lpcC8K9Ve7nn/Mm0tHdy9tQRfGfOaAAWTjbqWSt31bCjovfmHhEe+g34SqmlwHA3h24HbsNozjliSqnFGM0+5OXleXIrIUJaRkI0HV2ajq4u/nTRVLISo7n2ma8dx29/fRMA848e1uPa3LQ41stkrLDXb8DXWp/qbr9SagowGrDX7nOAtUqpWcB+INfp9BzbPnf3fxJ4EqCwsFDGjQnRi/SEaMfn847JprrJfaoFewevs9zUWN7bWE6XVRNhUm6uEuHgiNvwtdYbtdZZWut8rXU+RrPNDK31QeAt4GrbaJ3ZQL3WWppzhPBARkIUAAVZCUSZTWQmdr8Adtx9huPz0SOSelybmxZHp1VTXt/i+4KKoOVRG34f3gXOBIoBC/AdHz1HiLARHRkBwFHDEwGIjDAxLTfF8QJ47ruz+Gx7FdNzU3pcm5MaCxhj8nNS4/xXaBFUVDDNvissLNS5zMxvAAASRklEQVRFRUWBLoYQQamjy8oDH+1g8UljSI03avtaaw4bMOHW7upmTrn/M8f2s9+dxdzxmb4qqvAzpdQarXVhf+fJTFshhojICBO/WjjBEeyBAQV7gJEpsS7bT31R4tWyiaFBAr4QYSDKbOK2Myc4tpNiIgNYGhEoEvCFCBOLTx7L4pPHAHCol7w6ZbUWl3w8zpZuqeD9TTL2wlsqGlp5aOlO5tz3Cfd/sN0vz/RVp60QIgjddubRPLmshJUlNVjaO4mLcg0BJ/7xU8wmRfEfznTsW7PnEBc+ttKxfcH0bB749jF+K3Mo+HhrBSt31fDrsycCsLfGwsl//tRx/JFPi7nptPEDbqI7UlLDFyJM2VfFqmho5aLHVlBWawGg0+o6kOOX/93gsv3auv3sr5PhnYNxw4vreGp5qWO287tufin5I5upBHwhwszSG+cC3amUX1lTRtGeWm59baPb8xNiejYEfLa90ncFDEETbHMj3t5QjtWq6ei0AnDC2HR+dmoBYIyk8jUJ+EKEmdQ4o8O2ztIBQJstb/4XO3smV7O0d7KhrJ4p2clMz+se33/765t6besXPUXYmmp2HGxk6u8/5C8f7SAxxsy/fzCba0/IJyrCxNsbfN8/IgFfiDCTHOsa8EtrLD3OabW9BFYU1wBw8vgMXv/RHHbfdxaLjhkJwEdbKvxR3JBQ32L8Xb+/+aBjjeE2Wy0/JS6KJVfM4Pp543xeDgn4QoQZc4SJxGgzdS1Gm/EeN2vlFlc28c6Gcsc6ut8/cYzj2IOXHINJwYG6Fr7ZVydLJw5AXUs744cluOyLMXeH31MnDnNJleErMkpHiDCUHBdJvaWD19aWsaGsvsfxs/+2HICZ+anER0W4TPYymRRp8dG89PU+Hvm0mB/OHcO500YyaWSy38rvzopd1Xyzr44fneL7mvJgWK2a6qZ2Ljo2h2tOyOevS3cyNjOe+y+e5veySMAXIgylxEVS19LBjS+v7/O8vYcspMRF9difHh/FdtuIkyc+L+GJz0vYfd9ZPinrQFitmsv//hUAVxw3ytFsFQwOWdrpsmqyEmO44rhRXHHcqICVRZp0hAhDKbFR1FnaHRk4e1PR0Eaim1E6WUm+b34YjM92dI8a+ibI8v5vP2i8GIcFwd+ZBHwhwlCyrYbf2mHlO3PyHfvfvH4O37WtlmWX5Ka2nJ/eM+d+IK3d0x3k7313awBL4qq1o4srnjJ+eZxYEPhkddKkI0QYSo+PoqTK6JDNSIjm/Z+dREVDG9NyU5iWm8JvzpnIpN+8T3N7l9u8O6Mzegb8A3UtPZK0+UtDa4fRTGXp6DVtRCCs3VPr+JwQHfhwKzV8IcJQtlNgzk6JZcLwpB7pku0zbieOSOxx/YjkGMB4cdjd9rr7iVv+0NDSQXJsJFccl9djpnAg7apqAoKjOQck4AsRliZnd4+oyeplOKB9nPgJ4zJ6HLPPHP3eSd3NP/FRgavBNrR2khQTSUK02THOPRhsPdhITKSJ5b/6VqCLAkjAFyIsnTA2nYcuPYasxGgmZbsfThlrW2Er200zzeiMeL75zQKumzuWn3zLGAZZa+nZlPL8qj2U+iFlQENLB0mxZuKjzbR3Wunosvr8mQNRWtXMxBFJREYER6gNjlIIIfxKKcWiY7JZffupvQ5hfO57szhzyvBe2+VT4qJQSnHTaUdxSWEOK3bVuNSua5rauOONTVz6pJFpc2NZPfPu/8zRzOENz63czZo9hyjaU0t6fLSjnbw5CGr59ZYODja0+mVC1UBJwBdCuDUzP40lVxxLhKn/lL32ESgHnLJobrMNR6xoaGNnRSMPf7KT0upmvv+sd5Yxbevs4jdvbnakbu6yakeit8bWwAX8ji4rB+tbmXbnh5RWN5OeEDwBP/DdxkKIIc/e7FNWa+GDTQeJNJtckqsteHCZ43NpdTPVTW1keBAI2zutPV4cFx2bA7Z3U2VjK7lpgVms/a63t/Dcyj2O7Yz4vuc6+JMEfCGEx3JSjYD/3X92B+G+1vL469Id3H3elCN+3uYD9S7ZPZdcMYN5E7Ic+eb317VybIAmtC51Sio3OTuJC2bkBKYgbkiTjhDCY5luautaw9Qc9x3C/1q194if1dll5fwlKwB44/o53H3eZE6fNBwwfmmYTYr1AZxtG+c03v6t608k382chUCRGr4QwmOmXtr5T5803JGcbe74TI4dlUpNUxvPrtxDa0cXMbaRQINxoK7V8XnC8ESOye3O0x8fbWbW6DSKdh8a9H29ocuqKa7s7pTu7e8lUCTgCyG8ymxSjslPOamxnD89m8gIxZ8uMrJD/vsro3Zfa2lneFIMmw80sGl/PZfOyhvQ/auajIA/Mz/V7QsjNzWOTwK0IpfzCKSfzi8ISBn6IgFfCOEV9kD/+o/mcM4jRnrl4UkxPHjYguf2Fbd2VTZz/L2fOPZfUpjba414y4EGnli2i/svnkZlg9EZ/PtzJ7s9d3hyDNVNbbR3Woky+7fV2rl2//MF4/367IGQNnwhhFd8ctMp/PM7M5mSk8yPbas3Tc9L7XGePd3yQx/vcNlfZ1sVyp1bX9/Im98cYNP+eg7ZJnil95Lpc2RKDFobi7P7U0lVEz96YS0A//zOTL8+e6Ak4AshvCIvPY5TjsoC4KbTxrPj7jPc1rCn5CSTkxrL17trXfb31dFqP3b+khWssSUkc5e2GWBEsjFiyHlOgK8VVzay6NEvjXJFmx1/D8FGAr4QwuuUUr02pyREm1loG1Xj7L73trGqpKbH/sNnzb67sRyzSTlSPxwu2zZEdLebpRt95c63tzome71+/Ry/PXewJOALIfwuzdYckxRj5v2fncSs0Wlsr2jk0idXOc5ZuqWC51budkl3nJkYTWuHlcQYM6qXgf6j0+PJTIxmVYn/RurssM0qBhibGTzDMA8nAV8I4Xf2tMpTcpJ7pGbusmqsVs33nyviN29uZk+NBYDr541l0kgjS2eimxz9diaTIi8tztGGX93UhtXHKZMt7UbtPj0+qtcXUTCQgC+E8LvRGQkATMk2xtAvnNzdxHOgroVNB7oXVr/yH8aKUfOPHsZRw43c/LlpfS+0khYfxaHmdl5bW0bh3Ut5ftUe7v9gOy3tXV79HmDk9Glo7eTSmbl8fvM8r9/fm2RYphDC72aNTuO//3e8I4DnOeW9eeSTYkqqe2bUzE6JZaItD//50/tOV5AeH8VHWyoci7T/9q3NgNHR+8O5Y73yHexW7jL6HY4dlRoUq1r1RWr4QoiAmJmf5lg+MTLCxKbfnw7Af4r29RjBAzAsKYZzpo7k1etO4MIZ2X3eu7eUzu9uLPdqrvyOLivXPvM14PorJVhJwBdCBIWEaDPnT+8O5M4rcb3/s5MAo33+2FGp/baTz3FapesmpwlQ68vq+eN727xVZMfInHOnjeyzXyFYSMAXQgQNe9NOYrSZVbfO56QCI3CPy0wY1H2OyU0hNS6SgqwErrdNArN745sD3iks3UNG7eUMdsHd4CSECCsnFmTw0Mc7aWzrxGRSPHlVIbWWdsyDXCIwwqT47JfziIxQPdI12MfpH6l6SwcVja2MH5ZIs210TnyQt93beVzDV0r9RCm1TSm1WSn1J6f9tyqlipVS25VSp3v6HCFE6Dva1ilrFxsV0Wt7fH+SYyOJsy2svvLWb7Hhd6fxw7ljWL+vzmVs/2BYrZppd37IaQ8uo6KhleY2Y9RPWAR8pdQ8YBEwTWs9Cbjftn8icCkwCVgILFFKDT4PqhAirCREm4mKMHHVbO+uXjIiOZakmEhmj0kH4J2N5Ud0H+d8P2W1LY4mnfiooRHePH0tXQfcp7VuA9Ba23OSLgJesu0vVUoVA7OAlR4+TwgR4nbcc4bP7m2fuHXHG5u4fFbegNbrdVbT1L1sY1VjK/YBP2FRwwfGAycppb5SSn2ulLKniMsG9jmdV2bbJ4QQAZOVGOP4/PTy0kFfX+UU8HdVNbNuby1RZhP56cGbTsFZv68lpdRSwN0A09tt16cBs4GZwMtKqTGDKYBSajGwGCAvb2ALIAghhKdeXVvGD04eVLiipslo+4+JNPHnD7YDxq+G2CHSpNNvDV9rfarWerKb/97EqLm/pg2rASuQAewHcp1uk2Pb5+7+T2qtC7XWhZmZme5OEUIIr9t2sJHKQebMr7bV8E8q6I5Vw5Niejs96HjapPMGMA9AKTUeiAKqgbeAS5VS0Uqp0UABsNrDZwkhhMeuO6U7tcIra8sA2HawgTMe+oIfPl/U63XFlY088OEOTAqmZncvzp4UG/wTruw87Wl4GnhaKbUJaAeu0VprYLNS6mVgC9AJXK+19n7WIiGEGKRfLZzAj04Zy5Tffcif3t9O0e5aPtlmjDfZWt5AS3sXX+8+xOLni/jq1lNJti3JeNfbW2m0jcpJje9ebWv8sET/f4kj5FHA11q3A1f2cuwe4B5P7i+EEL7gnAbBHuzt3tlYzi/+ayRdW7uvlnm21asyEoxUD1cfP4qpOUYN/6SCDBYPsh8gkIbGWCIhhPATe7AH+KrkkCPgt3dZyU6J5bfnTCLCpPji5nnkpMYGdf77w0kuHSFEWFp923xudEqs5s7jn++irdNojW5s7SA9Icoxdj83LW5IBXuQgC+ECFNZSTEu2Tl7U2+bXVte19rrwulDhQR8IUTYyk2Lo/TeM7n/4mm9ntPQ0sHavbVsr2h0mbg1FA3t15UQQnhIKcVFx+aQnRJLSlwklzyxksbWTk4Ym86KXTXUt3Q4VrW69cwJAS6tZyTgCyEEcPxYI7Haxt8ZyX3X7a3l/CUrONTcwf66FlLiIqWGL4QQoSg93hiG+YPnjMlY/S2cPhRIG74QQriRlRTtsl1v6ejlzKFDAr4QQrgREzk0EqINhgR8IYToxV8unsaSK2YAkBIX1c/ZwU/a8IUQohcXHpuD1pqbFoznzKkjAl0cj0nAF0KIPiil+Mn8gkAXwyukSUcIIcKEBHwhhAgTEvCFECJMSMAXQogwIQFfCCHChAR8IYQIExLwhRAiTEjAF0KIMKG01oEug4NSqgrYc4SXZwDVXizOUCDfOTzIdw4PnnznUVrrzP5OCqqA7wmlVJHWujDQ5fAn+c7hQb5zePDHd5YmHSGECBMS8IUQIkyEUsB/MtAFCAD5zuFBvnN48Pl3Dpk2fCGEEH0LpRq+EEKIPoREwFdKLVRKbVdKFSulbgl0eXxNKZWrlPpUKbVFKbVZKfXTQJfJH5RSEUqpdUqptwNdFn9RSqUopV5RSm1TSm1VSh0f6DL5klLq57Z/05uUUi8qpWICXSZfUEo9rZSqVEptctqXppT6SCm10/ZnqrefO+QDvlIqAngUOAOYCFymlJoY2FL5XCdwk9Z6IjAbuD4MvjPAT4GtgS6Enz0EvK+1ngBMI4S/v1IqG7gBKNRaTwYigEsDWyqf+Sew8LB9twAfa60LgI9t21415AM+MAso1lqXaK3bgZeARQEuk09prcu11mttnxsxgkB2YEvlW0qpHOAs4KlAl8VflFLJwMnAPwC01u1a67rAlsrnzECsUsoMxAEHAlwen9BaLwMOHbZ7EfCs7fOzwHnefm4oBPxsYJ/TdhkhHvycKaXygenAV4Etic/9FbgZsAa6IH40GqgCnrE1ZT2llIoPdKF8RWu9H7gf2AuUA/Va6w8DWyq/Gqa1Lrd9PggM8/YDQiHghy2lVALwKvAzrXVDoMvjK0qps4FKrfWaQJfFz8zADOAxrfV0oBkf/MwPFrY260UYL7qRQLxS6srAliowtDF80utDKEMh4O8Hcp22c2z7QppSKhIj2L+gtX4t0OXxsTnAuUqp3RhNdt9SSv0rsEXyizKgTGtt//X2CsYLIFSdCpRqrau01h3Aa8AJAS6TP1UopUYA2P6s9PYDQiHgfw0UKKVGK6WiMDp53gpwmXxKKaUw2nW3aq0fCHR5fE1rfavWOkdrnY/xv+8nWuuQr/lprQ8C+5RSR9l2zQe2BLBIvrYXmK2UirP9G59PCHdSu/EWcI3t8zXAm95+gNnbN/Q3rXWnUurHwAcYvfpPa603B7hYvjYHuArYqJT6xrbvNq31uwEsk/CNnwAv2CozJcB3Alwen9Faf6WUegVYizESbR0hOuNWKfUicAqQoZQqA34L3Ae8rJT6HkbW4Eu8/lyZaSuEEOEhFJp0hBBCDIAEfCGECBMS8IUQIkxIwBdCiDAhAV8IIcKEBHwhhAgTEvCFECJMSMAXQogw8f9bXsiJP9y83wAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x10a102e80>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"_ = plt.plot(t, x)"
]
},
{
"cell_type": "code",
"execution_count": 33,
"metadata": {},
"outputs": [],
"source": [
"# The moving window parameters\n",
"window_length = 100\n",
"window_step = 50"
]
},
{
"cell_type": "code",
"execution_count": 34,
"metadata": {},
"outputs": [],
"source": [
"def np_rolling(x, window_length, window_step, func):\n",
" vert_idx_list = np.arange(0, x.size - window_length, window_step)\n",
" hori_idx_list = np.arange(window_length)\n",
" A, B = np.meshgrid(hori_idx_list, vert_idx_list)\n",
" idx_array = A + B\n",
" x_array = x[idx_array]\n",
" return t[vert_idx_list + int(window_length/2.)], func(x_array)"
]
},
{
"cell_type": "code",
"execution_count": 42,
"metadata": {},
"outputs": [],
"source": [
"# Our aggregation function\n",
"# x is shape (window_size, Y). Must return array with shape (Y, ) \n",
"def agg_f(x):\n",
" return x.max(axis=1)"
]
},
{
"cell_type": "code",
"execution_count": 43,
"metadata": {
"scrolled": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"(18, 100)\n",
"(18,)\n"
]
}
],
"source": [
"x2, y2 = np_rolling(x, window_length, window_step, agg_f)"
]
},
{
"cell_type": "code",
"execution_count": 40,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x10a7da048>"
]
},
"execution_count": 40,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlsAAAFpCAYAAACrn+1KAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzs3Xd4VGXax/HvmZKZ9N5IAiEU6YQqIFWaFGVdFVFpuhZs2HZdy+qK/bWs69pdXQEVu1JEECxIFaT3TiCNJKT3ZDLn/eNMJhOSQJKZSQK5P9eViylnznMiGH4853nuW1FVFSGEEEII4R665r4AIYQQQoiLmYQtIYQQQgg3krAlhBBCCOFGEraEEEIIIdxIwpYQQgghhBtJ2BJCCCGEcCMJW0IIIYQQbiRhSwghhBDCjSRsCSGEEEK4kYQtIYQQQgg3MjT3BTgKCQlRY2Njm/syhBBCCCHOa9u2bWdUVQ0933EtKmzFxsaydevW5r4MIYQQQojzUhTlZH2Ok9uIQgghhBBuJGFLCCGEEMKNJGwJIYQQQrhRi1qzJYQQQlxoysvLSUpKoqSkpLkvRbiJ2WwmOjoao9HYqM9L2BJCCCGckJSUhK+vL7GxsSiK0tyXI1xMVVUyMzNJSkqiffv2jTqH3EYUQgghnFBSUkJwcLAErYuUoigEBwc7NXMpYUsIIYRwkgSti5uzv78StoQQQogLnF6vJz4+nu7du9O7d29effVVrFbrOT+TkJDAokWLmugKWzcJW0IIIcQFztPTk507d7Jv3z5Wr17NihUrmDdv3jk/I2Gr6UjYEkIIIS4iYWFhvP/++7z55puoqkpCQgLDhg2jb9++9O3bl40bNwLwyCOPsG7dOuLj43nttdfqPE44T3YjXqgSt0DCOogdBjEDL/xxhBDiIjBv2T72p+S59Jzd2vjxzyu7N+gzcXFxVFRUkJ6eTlhYGKtXr8ZsNnPkyBFuuOEGtm7dyosvvsgrr7zC999/D0BRUVGtxwnnOR22FEWJARYC4YAKvK+q6uuKogQBXwCxQAIwVVXVbGfHE2gBaMGVYCkDvQFG/xNCL3HpECWWCoqT9xO46QWosGjjXPEiRPQEnUH70hurHtufG0Gnd3jPCLrzTKBKoBNCCLcpLy/nnnvuYefOnej1eg4fPuzUcaLhXDGzZQEeUlV1u6IovsA2RVFWA7OBn1VVfVFRlEeAR4C/u2A8cWQ1WGxbUCvKYNXjLh/CbPuyqyiD5Q828mzKWeHLIYxZK6AgDVBBb4LZ30vgEkJcsBo6A+Uux48fR6/XExYWxrx58wgPD2fXrl1YrVbMZnOtn3nttdfqdZxoOKfDlqqqqUCq7XG+oigHgChgCjDSdtgCYA0StpyXlwp7vrI90WmhZeLLEO66/8GvfXcjlgqVS5RTPG2cj54K9AYjyrhnIKgjWMvBaoGKci0sVXtuqfn4XO+d3g0Fp7WBK0phyV1wzf8gspfLvh8hhGhNMjIymDNnDvfccw+KopCbm0t0dDQ6nY4FCxZQUVEBgK+vL/n5+fbP1XWccJ5L12wpihIL9AE2A+G2IAZwGu02Y22fuR24HaBt27auvJyLz5mj8PHVUJyl3dIrL3LLrbetljQAdqodOVIWzSDdAUaO/jMDL73inJ87llGAr8lAmF8D/jWUuAUWXKXNnCkK5KXAe8Og+9Uw8jEI7ezMtyKEEK1CcXEx8fHxlJeXYzAYmDFjBg8+qN2NuOuuu7jmmmtYuHAhV1xxBd7e3gD06tULvV5P7969mT17dp3HCecpqqq65kSK4gP8Bjynquq3iqLkqKoa4PB+tqqqgec6R//+/VVZjFeHpG2w6DpAgelfQ5s+bhlGVVXaP/pDjdcfGNOZ+8Z0OudnYx9Zjsmg49CzExo2qOOarZDOsOlN2PQ2WIqh1zQY8TAENa5FghBCuNuBAwfo2rVrc1+GcLPafp8VRdmmqmr/833WJaUfFEUxAt8An6qq+q3t5TRFUSJt70cC6a4Yq1U6+pO2IN7DB/6yym1BCyCnqLzW1/el5PLqqkNsO5lV6/t5JdrnSi1WGhzgYwbCsIe0Xz0D4PJ/wP27YdBdsPcbeLM/fP+ANuslhBBCXGCcDluKVsP+Q+CAqqr/cnhrKTDL9ngWsMTZsVql3V/CoushKA7+shqCO7h1uNxiLTS9el1vIhxuB649ksEbvxzlmnc2ATD3sx0MfO4nCkstABw+XXXfPyGzyPkL8Q6B8c/BfTuh7yzYvhBej4eVj0FBhvPnF0IIIZqIK2a2LgNmAJcrirLT9jUReBEYqyjKEWCM7bloiE1vwbe3QdvBcPNy8K112VuDWa1qnbNPlTNU/p5GJvWKBCA60JOS8qq2DyXlFSzdlUJ6filT39tEeYWVnw5UTVyOemUNGfmlLrlW/NrA5H/Bvdug57Ww+R14vTf8/AwUSyURIYQQLZ/TYUtV1fWqqiqqqvZSVTXe9vWDqqqZqqqOVlW1k6qqY1RVrf3+k6hJVWH1k/DjY9BtCtz0NZj9XXb6m+f/wYiX15CaWwxoC9tv/mgL+1JyySvWZqr8PI38ddwlvHRtL769a0i1zydkFtof70vJ40BqHot3JDMgtmpJ3q+HXHzXODAW/vQ23LUZOo+Dda9ooWvty1Caf96PCyGEEM1F2vW0NBXlsPgu2PA69P8LXPsRGF1X6yS3uJzfDmdwKquIwS/8QmJWEaNf/Y1fD2Xw4foT9pktP08Dnh56pvaPIczXTLifyX6OQ6erh5v9KXmczivh8i5VM28nzhTiFqGd4br5MGc9tB0Cvzyrha6Nb0J5sXvGFEIIIZwgYaslKSuEz2+EXYtg1OMw6VWtAKgLbT6eWe35kp3J9sdGnY4825otP7Ox2nFf3jGYnlHa7NrmE9UnKVfv10pFdArz4bZh2q7Bvcm5Lr3uGiJ6wo2fw19+0h6vehz+0wf++ECrrC+EEEK0EBK2WoqiLFg4Rdt5OPnfWrkDRXH5MNtP5VR7/sqqwwR4GenbNoCTWYUOM1vVw1a7YG8en6Rted2VmIOHQce/pvYG4OeD2i3DSyJ8eXxSN24YGMPupNyG70psjJgBMHMJzPoeAtrC8ofgzX6w41OtzZAQQrQCiqIwffp0+3OLxUJoaCiTJ092+tyZmZmMGjUKHx8f7rnnnmrvbdu2jZ49e9KxY0fmzp1r/7mflZXF2LFj6dSpE2PHjiU7u3FrbJ9//nmnr7/S/PnzCQ0NJT4+nvj4eD744AP7ewsWLKBTp0506tSJBQsWuGzMShK2WoKcRPjfeEjdDdctgP43u22og6fz6Bbpx/HnJ9I+RCtY1ynMh9hgb05lFpFXbEGngLdHzRm1tkFegLZOK8jLgz/3jbbPdg1sH0SM7f0eUf7kFpdzKssFuxLrq/0wuOVHbX2bZ6BWif7tQbD3W7Baz/95IYS4gHl7e7N3716Ki7XlFKtXryYqKsol5zabzTzzzDO88sorNd678847+e9//8uRI0c4cuQIK1euBODFF19k9OjRHDlyhNGjR/Pii43bI+fKsAVw/fXXs3PnTnbu3Mmtt94KaMFw3rx5bN68mS1btjBv3rxGh8O6SNhqbukH4MNxkH8aZnwL3a5yyzAHT+exNzmXzIIywvxM6HQK867SWvzEBnsT4W8mPb+UtLwSAr08UGqZVYvwM2MyaH9kAr09AAiy/Xpp+yD7cYPjggGYt2x/08xuVVIU6DQWbv8Npn6s3YL9+mZ4bzgcWgGnNsO6V7UiqkII0ZwSt7j859HEiRNZvnw5AJ999hk33HCD/b0tW7YwePBg+vTpw5AhQzh06BCg9UO85ZZbANizZw89evSgqKj6P5S9vb0ZOnRojV6Jqamp5OXlMWjQIBRFYebMmSxevBiAJUuWMGuWVv1p1qxZ9tfrkpqayvDhw4mPj6dHjx6sW7eORx55xF4Z/6abbgLgk08+YeDAgcTHx3PHHXfYWwr5+PjwwAMP0L17d0aPHk1GRv1LBP3444+MHTuWoKAgAgMDGTt2rD00uopL2/WIBjr1OyyaCgYz3LwCInq4bagr/r0OgKgATzqF+QAwtGMIT07uxjV9o/luRxIWq8pPB9Lo07b2Qv86ncKQDsH8eiiDQC/tNmOFVQtT0YGe9uPiQn3oGObDLwfTWbn3NBN6Rjp9/ftT8vhk80luHdqeuFCfcx+sKFpo7TIJ9nwNa16Az6bZbssqWsPrWUul4bUQwvVWPAKn95z7mNI8SNsLqhUUHYT3AJNf3cdH9IQJ558ZmjZtGk8//TSTJ09m9+7d3HLLLaxbp/3s79KlC+vWrcNgMPDTTz/x2GOP8c0333DfffcxcuRIvvvuO5577jnee+89vLy86vWtJicnEx0dbX8eHR1NcrK2DjgtLY3ISO1nf0REBGlpaec816JFixg/fjyPP/44FRUVFBUVMWzYMN5880127twJaBXcv/jiCzZs2IDRaOSuu+7i008/ZebMmRQWFtK/f39ee+01nn76aebNm8ebb75ZY5xvvvmGtWvX0rlzZ1577TViYmJITk4mJiam1u/DVWRmq7kcWqGt0fIK0arCuylo5RSV8dTSffbnyTnF9tkonU7hlqHt8fcyEm4rYJpdVM6A2KBazwUwvnsEAKfzSgAos2i36CL8Pasd982cIXh76Pl2h2v+wH6/O4VFm0+xcNPJ+n9Ip4fe18M9f0CXyVpJDdWq9WFMWOeS6xJCiAYrydV+FoH2a4lrNhT16tWLhIQEPvvsMyZOnFjtvdzcXK677jp69OjBAw88wL592t8LOp2O+fPnM2PGDEaMGMFll13mkmtxpChKrXdLHA0YMICPPvqIp556ij179uDr61vjmJ9//plt27YxYMAA4uPj+fnnnzl+/Lj9+7j++usBmD59OuvXr6/x+SuvvJKEhAR2797N2LFj7TNvTUFmtprD9o9h2X0Q2UtbY+Qd4rah5i3bz3dnBZ4gH48axzk2j75xYN0Nwcd0C4dv92DQaf/jdI30ZUtCFm38q08v+3sZGdwhhMQGrNvam5zLiTOFXNm7TY33KivbH8soqPf57PRGuOw+OLIaKkq1H27B5+7zKIQQjVKPGSgSt8CCq7R/+Ok94JoPXDbTftVVV/HXv/6VNWvWkJlZtfv8iSeeYNSoUXz33XckJCQwcuRI+3tHjhzBx8eHlJSGtUSLiooiKSnJ/jwpKcm+Tiw8PJzU1FQiIyNJTU0lLCzsnOcaPnw4a9euZfny5cyePZsHH3yQmTNnVjtGVVVmzZrFCy+8cN5rqy3cBQcH2x/feuutPPzww/bvY82aNdW+D8f/Pq4gM1tNSVW1e/RL74G4EdoOOjcGLYDtp2ou8vM9q6wDQIQtLF3WMRh/r5rvVwrxMfHu9L68P0Pru/noxK58fvsgOoXX/FdIpL/ZPgNWH5PfWM+9n+2wP/96WxJT3lyP1araw9bxjEbW74oZCLO/hwG3gYeXVioi51TjziWEEM6IGagtZbj8cZcvabjlllv45z//Sc+ePau9npubaw9C8+fPr/b63LlzWbt2LZmZmXz99df1HisyMhI/Pz9+//13VFVl4cKFTJkyBdBCX+WuvgULFthf37JlS40QBXDy5EnCw8O57bbbuPXWW9m+fTsARqOR8nLt5//o0aP5+uuvSU/XdsBnZWVx8qR2t8NqtdqvfdGiRQwdOrTGGKmpqfbHS5cutTeVHj9+PKtWrSI7O5vs7GxWrVrF+PHj6/3foT5kZqupWK3w46Ow+V3oeR1MeRsMNWeYXDukWq1tzpiuYfx0IL3WnYZRAZ68N6MfQzueP/xd0aNqDZbZqGdQXHCtx0X4m8kpKqekvAKz8dz1wqzWqoX0lcf/9atdAOxIzLGHreScYg6ezqNjqA8GfQP/rRAzUPuKvxE+/hPMnwSzl2slI4QQoilV/jxysejoaObOnVvj9YcffphZs2bx7LPPMmnSJPvrDzzwAHfffTedO3fmww8/ZNSoUQwfPrzGTFRsbCx5eXmUlZWxePFiVq1aRbdu3Xj77beZPXs2xcXFTJgwgQkTJgDwyCOPMHXqVD788EPatWvHl19+CcCpU6fw9Ky+7ARgzZo1vPzyyxiNRnx8fFi4cCEAt99+O7169aJv3758+umnPPvss4wbNw6r1YrRaOStt96iXbt2eHt7s2XLFp599lnCwsL44osvaozxn//8h6VLl2IwGAgKCrKHzqCgIJ544gkGDBgAwJNPPklQUN3LaRpDadLdYufRv39/devWrc19Ga5nKYPFc2DvNzDoLhj3HOjcP6mYlF3E0P/71f581z/HsSsxh8s6hqDXub6G19m+2ZbEQ1/tYs1fRxJrKzNRl6Pp+Yz511oANjxyOVEBnvR86kfySyzcP6YTvx5Mp6isgiPp2m3E6YPa8uyfep7rlOeWvB0W/gk8A7QZLwlcQohGOnDggH2WRJzb3/72N2bMmEGvXr1cel4fHx8KChqxzKQBavt9VhRlm6qq/c/3WbmN6G6l+bDoOi1ojZkH459vkqAFcDBVa6vzwJjOLLxlIP6eRoZ3Dm2SoAXQNljb0VKfdVZbE6pud2bkl3Iso4D8Eq0o6dH0AnKLy+kSWbVbZ/GOhq0tqCGqL8z8DopzYP5krdaZEEIIt3r55ZddHrQuBBK23KkgQ/uL/MQ67bbh0PvdUhW+NmcKSrl1oTZLeMeIOIZ3Dm2ScR31jPLHqFfYevL8xeF2JlZVtk/PK+HLPxLx0OvoFe3PsYxCcovL8fc08M5NfQHoHH6e8g/1EdXPIXBNksAlhBAXKHfPajlLwpY7JG6BH/+hFdPMOATTFkGfm5r0En6xtdAJ9vY473opdzEb9XSJ8Dtnn0RVVfnyj0R2J+USE6Tdx3/4m93sT82jc4QPg+KCOZCaR3ZROf6eRib0jOSK7hHklVhYcyidU5lOVqmP6gczvoPibFgwGXKTzv8ZIYQQogEkbDVAeYWV4+e7JZa4RZvN2vQG5KfAFS/CJVc0zQU6OGpb27R87rAmH9tRiI8H2UV1N4Y+lJZvD1cD2gUxtls4OUXlrDtyho6hPozuUrVIc3Cctng/zM/E0fQCZn/0B/d+vqOuU9dfdD+YsVjrTzl/kgQuIUSDtaT1z8L1nP39lbDVAA9+uYvLX/2N3KLyug86sVar5QSg6KE4s+5j3ejQ6Xy6RfrZSzo0lwAvD/tOwtqk5pZUO/aWy9rbn3eJ9GNAbBBT+0czJb4NQzpoux7DfE32Yw6fznfNDzkJXEKIRjKbzWRmZkrgukipqkpmZmaNdkUNIaUf6iktr4Rlu7RF2YfT8+nfLpCisgq8TWf9J2w/XGsHY7Voxepim2dm6XBafrV+hc3F39N4znCamlMVtkZcEkqvaH8i/c0MiA1i1uBYdDqFl67tXe0z1/aLQa/ToVPghRUHWbH3NBNd0BJIC1zfwcdXa7OTs5eDv2sauQohLl7R0dEkJSU1qB+fuLCYzeZqrYkaSsJWPX29rWqm43BaPolZRTz45S5+fmgEHRx79VUWz0xYpwWtZui/t3BTAqm5JXSOqFlotKn5eRrJL7VQYVVr3QWZmluMToGDz0zAw9bketOjo895zgh/M3eO7ECFVeXD9SdYvT/NNWELILo/TP/WFrgmSeASQpyX0Wikffv25z9QtFpyG7GeUnOLCfAy4mnUczyjkI3HtNuDH9fWqy9mIAx7qNkaHT+5ROt51fF8DZubgL+nEVWFl1YerPX91NwSQn1N9qDVEHqdQkyQF6dz61+lvl5iBmgzXIVnbLcUXduQVAghROsiYaue0vJKifAzE+LrQVZhGZYKrYnoN9uTKC6raOarq93gDrVXdm9KfmZt8vS9tcdrfT+7sIwQH1Ot79VHhJ+ZNFtLoJLyCpbsTHbNugnHwLVgMuQ5WddLCCFEqyVh6ywfrDvOsJd+Yd2R6vfe0/NLCfU1Eeilha3knGIA8kss9p1/LUFlCJx7ecdaeyA2tcJSyznfzynWSjo0Vrif1n9RVVWe+X4/932+k231qOtVLzEDYMa3tnppkyRwCSGEaBQJWw7S80t4dvkBErOKeWfNMQDeXnOUg6fzSM8rIczXTKCXBzlFZSRlF9O9jVbR/ERmI5sju0HXJ1cCEOjt3r6L9TUlvmq9U4W15oxTTlEZAedofH0+4X4misoqyCux2CvVF5e7cKYxZiBM/6aqQK0ELiGEEA0kYctBel5V02aLVaWkvIKXVh5iypsbyMgvJdzPRKCXkfT8UtLySuylCBLOtIywVWFVKa/QAk1QCwlbgd4ezLuqO0Ct9bZyiy34ezb+Witb+Gw/lU2ZRZvVy8gvpcSVgavtpbbAlSaBSwghRINJ2HJQufanW6QfR9LyybPVhyq1WLFYVcJ8TQR4eZCaW4JVhY5hPoT7mTjpbBVzF0nKrroOn7NLUjSjUFtdrIz80mqvq6pKbrFzM1uXtg/CbNRxy/w/2H5Ka/nz4Je7uPrtjY2/4Nq0vVTbpWgPXKmuPb8QQoiLloQtB2m2ma2eUf5kF5WTddZMTLifudqMUXSgFxF+ZtLzXbwbrpEc1445E2BcrbIIaVJ2Me/9dsw+61RcXkF5herUmi2zUc8l4b6cvSb+QGqe6wsMOs5wLZDAJYQQon4kbDmoDE3tQ70BSLEtgq8U5meiXbCX/XlUgCehvqYaMzbNJd12HW/f1Jd+7Zq/oGmltkHaf7M3fjnCCysO8vavRwHIsRU7DXAibAF11thyy8aFtoO0wJV/WgKXEEKIepGw5WDLiSxigjwJ8tJmrxZsrF5DK8zXzPjuEUT4menTNoC2QV6E+ppJzy+tEcyaQ1ahNhN3uUM/wZYg1NeE2ajjYGo+ABkF2nVW7ugM93OupdBtw+LoFe1f4/WPf6+lBportB0EN32tBa0Fk7XgJYQQQtRBwpbNiTNaodIbBrbF11Yb6rfD1cs/aKFBz08PjeCbOUPQ6RRCfU1kFZYx5MVf2H7KRSUHGimrsAxvDz1mo75Zr+NsiqIQ7memzFaWorIcxP6UPAC62XZ1NpZOp9DG39P+fOQloUztH83CTSdZf+QMU9/dxLT3Nzk1Rg3tBmszXHmp2houCVxCCCHq0HJWUTezD9ZpRTev7hPFsfTquws3Pzaa0nKrPcQ4Lj53bIq8bFcKfdsGNsHV1i6rsKzFlHw429iu4Xyw/gQAS3eloFPAZNAT7O1R7b9hYz06sQveJgPPXd0Dk0FHqcXK2sNneP3nw/yR4KYQ3G4wTP8aPrkWFlwJs5aBb4R7xhJCCHHBkpktIDGriE83n+LqPlFE+nvaZ7YqBXl70NZhrZajUIegsMnWwqe5ZBWWtZiSD2d7fFLXas9/PZTB/tQ8urXxQ1Fq9kxsqHbB3rw6tTdmox5FUTAb9dw+PK5a0Kq8zepS7YZogSs3WQtc+WmuH0MIIcQFTcIWsDspF4C/DNUaifo4hK3Xp8Vj1Nf9n8lxVubg6XxSc5tv7VZLDluKotiD6bQBMeQWl7MnOZcubmyW3adtQLXnR9K0NWMFpRZybWU9XKLdELjpKy1wfTAavn8I9i+FkjxqbJMUQgjR6shtRGBvSi5GvUKncK1xs+PMlmMF9NqEnnUL7ONNJ3n4ii6uv8h6OFNQav8eWqJ1D48CYOXe03z+RyIAPaJqLmx3lbN7Lp621VEb96/fSMktIeHFSa4bLPYyGPs0/PAQbP1A+wJQ9GD2B88A7VdzgO1xgMPrAXW/r6tl/V3iFkhYB7HDmq3ZuRBCiPqTsAWcyiwiJtALk0H7i62y7tONl7Y972fPDltvrznG5V3C6B/btKUXrFaVjHytWXZLVbnmLdhHm33rEuHLpDrKNrjC2b836XmlqKpKSm5V42qXbiYozUVVdCiqFRUFtcNodG16Q3EOlORASa72ODex6jXruXtHYvKzBTBbEFOtcOp37VeDGWYtlcAlhBAtnIQttBkPx/IDJoOenU+Oxa8ejZwrA5ojl96iqqesojIsVtXpMgpN4bIOITx3dQ/+3Ccawzlu0Trr7CD13A8HCPOrCmD7U/NctqGhsNTCtvIuDMKITi2nHAMv5k1i3oxb6v6QqkJ5kS145Wrh6+xgZn/N9n7mUVBtrYgspdoMl4QtIYRo0SRsobXpGXDWTFSAV/3XPk3oEUGv6AD+b+VBADLdsRD7PCpbDYX7Ob+zz910OoWbLm3XLGPf9/lO++NTmUUuC1sv/3iI+Ruhr/Iog3QH+N3ale2J4TxSVoGnRx2zZ4oCHt7al/+5b1fbJW6BBVeBpRiwQnBHl1y/EEII92n1C+RVVSU9r7TajEdDvTO9H3eO7MBfx3UGtLVT7rZsVwqxjyy3z6JVNtEO9W35M1tNqXe0P0M6BBPpX/O/S2VAdcbR9HyyC8vYdlLb9bhd7czbFVPYrmp/Fvan5jk9RjUxA7Vbh0PuA4MXbH4frFbXjiGEEMKlWv3M1v7UPMoqrIS7IKTcc3kn3vvteJO073lv7TEAjmcU0KdtIBm2gOeKmlUXkyX3DAW0HYiJWUX8/ZvddArzZeXeVHsvTGeM+dda++P2Id6cOKPVaBvROZTfDmfYdqe6uPZazEDtK6QTLL0HtrwHg+507RhCCCFcplWGrZyiMob+369UWFWKyyvwMxsY2y3cJecO9TXZ/8J1J7NtrVhl7ahs268ttahpc/MxGega6cdSW/jacSqbNCcbiFut1cs6LJ87lG5P/gjAf6b1offTqzid68Ym5X2mw4Fl8NNT0HGMFr6EEEK0OK3yNuKupFwKSi0Ul2sLjW8a1I6YoNqLljbUhJ4RrDmUwWFbTSd3MRm137qU3BL+9NYGXlhxEA+9Du+61geJakJ9TaQ7eRuxsKxqJ2HbIC+8PKr+7eLnacDLQ09KTgkVVpXlu1NrhDOnKQpc9R9tV+LiO8Fa4drzCyGEcIlWGbZUh0KTc0Z0YM7wDi47959sdbleWnnQ9X+5OtDrtN+61JxidibmAOBt0rukGntrEO6nNRBvrDKLlV2JufbnJoP2+/H+jH4sufsyFEUhwt/M6bxiPt6UwN2LtvPdjmRnL7sm3wiY9Cok/QEb/+P68wshhHBaqwxbOUVVpRkeHn8J/l7nL/FQX7Eh3gD8dCCdlfsITogcAAAgAElEQVTc15w4p0i7bbhyb9UY2UVNX3LiQhXuZyItr6Ra8G6Il1YeZPqHm+3PL+sYAsC47hH0jtEq17fx9yQlp4RTWVpXAWdvW9apxzXQbQr8+jyk7XPPGEIIIRqtVYatbFtQGd89HJ3OtTNBjq19Nhw949JzO6pcq3W8CdaHXYzC/cyUlFv5z89HG/X5rSerei6+Pi2ev9fSNSDC38zOxBz+t0FrwJ1XfJ4Cpo2lKDDpX1oB1O/mQIWEbiGEaElaZdiqnNl668a+bjn/Lw+NIDbYi30peaiqyu6knEbPoNQmOaeY5JyaPRg9XVkN/SIXZiv++tpPhxv1+WCHjQgdQn1qraV1drmJxOyiRo1VL94hcOXrcHo3rH3FfeMIIYRosFYatsrwMxvcVr08LtSHQXHBJGUX88G6E1z15gZ+OpDusvM/uXgvCtX/Mn90QhdWPzjcZWNc7Bw3EpRaGr6wPMOhllpdnQbKKrT6Vx62P2eHT7t30wRdJ0OvabD2ZUjZ4d6xhBBC1FvrDFvF5W4vkRAd6MmZglKe++EAgL3opbOsVpXfDmdwXb8YxjmUqxjQPojoQNfsqGwNBrav6hiQktPwtVQpOcWYDDoCvYx1FsQddUkYAEvuuYz7x3TiaEYBhaVuupVYacKL4BMG390J5W4sOyGEEKLeWmXYyi4qJ8DTdYvia+MYfML9TC5bv5VXUo7FqtIp3IdHJ3Zlav9oAGIkaDWIr9nIl3cMBiAxq2G399LySjhTUMbfxl/CjifH1dnMelBcMCdemEjXSD96RwegqvDLwXSX3lKuwTMQrnoTMg7AmufdN44QQoh6a5VhK6eorEG9DxtjQPsgurfx00pLjOjAnuRcDrigdUtl38UQHxNmo56Xru3NiRcmEiqV4xusY5gPALtspTPqa8cp7fg+9eirWFmKY3CHYADu/WwHcz7Z1qDxGqzTGOg7Cza+Aac2n/94IYQQbtVKw1Y5gS4s91CbqABPls8dxiMTujCxZyQAG49lOn3eyl2IQQ63QaW2VuMEeXvQI8qPdQ2YdVy2K8Uelrq38av358xGPa9Pi6dP2wB+3JfGvZ/taNRasXob/xz4R8PiOVAmO1aFEKI5tcqwld0EM1uOwv3MRAV4st0F67YybQuzg32kLY8rDO0Yyo5T2fVeS3XvZ1ULz+u6fViXKfFRfHH7YGYObseyXSn8uC+tQZ9vEJMvTHkLso7DT/PcN44QQojzanVhy1JhJb/EQoCbZ7bONqRDMGsOpVPg5ALpyrpQwd5y29AVhnUKobxC5Z9Lz18MtNy2u9AZHgYdT0zuhodex97k3GrvbTuZxdfbkpwew679cLh0jtao+sTa8x8vhBDCLVpd2Mot1mpsBTbhzBbA1AExFJZVsPZwRqPPkVtUzv7UPCL9zbJGy0UGxQUTF+rNb/X4fVm6M8UlYxr1OrpE+rIvpXrYuuadTfz1q132puIuMfqfENQBFt8NJc6vGRRCCNFwrS5sVVaPb+qZrZ5R/uh1CvtTGv8X3q4kbWH2q9f1Ru/iyvetlV6nML57BNmFZefdJbg/Nc9eOHbm4HZOjdu9jT/bTtZ++7LPM6udOnc1Hl5w9buQlwSr/uG68wohhKi3Vhe2Ttm2+UcFeDbpuGajnrgQb6d2JK49nIFep9Aj2t+FVyaCvDywWFUKSi2oqsoXf5yyr41zlJ5fSrifiaPPTWDeVd2dGrN7Gz9Kyq30nrcKgO2nqq/nKyl34eL5mIEwZC5sXwBHXBjkhBBC1EurC1vHM7SdWXGhPk0+dlyotz3sNZSqqny7I5lx3cLrrFguGqeywO3tC7fR5YmV/P2bPbz+8xGs1qqZro82nGDZrhSCvD0w6HVO7wCNtzWrtlhVjqYX8Oe3N1Z7P6WWdkxOGfUYhHaFpfdCsWsK7AohhKifVhe2jmUUEuBlrFY6oamE+JjsdbIaKj2/lKzCMnu9JuE6lWVANh3PpNSiLYJfuOkk8U+vst9afN7WCaCyr6azekT5M2dEBwB21lLnKzHbxWHLYIKr34HCDFjxd9eeWwghxDm1qrBltaqsO5JBr+iAZhk/2MdEdlEZlkbsajuWUQBAXEjTz8hd7BxbN80eEmt/nFdiYdvJbCwVVvw9tWPO1HJ7sbEm9owAqt9CXP/3UYAbZrYA2vSB4X+D3V/AgWWuP78QQohataqwtfVkNknZxfy5T1SzjB/i44Gqau2CGupgqtbEuEOYt6svq9UL99MaevdvF8hTZ63FuvbdTXR8fIU9ZL0zvZ/Lxm0XpP1eVpaAmHt5R/su2bxi18yg1TDsIYjsDcvuh0LXtJASQghxbq0qbPVvF8iiWy9lfPeIZhm/sjZWQ5tSq6rKO78do3sbPyJswUC4TlSAJ4vvvoxFtw0653H/ndmfyzqGuGxcfy8jvmaDPWxNH9wOLw89OgWn67HVSW+EP70LpXnw/f3gzj6NQgghgFYWtnQ6hSEdQ/D0aFjlb1epXCc255NtVFjr/5dcZmEZGfmlXNcvWlrzuEl8TAAehur/O1zVu439sbeHnrHdwl0+bkygF5V/FPw9jSiKgo/JQH6Jm8IWQHg3bcH8gWWw52v3jSOEEAJogrClKMoViqIcUhTlqKIoj7h7vJasV7Q/Xragt/mE1ifx3d+OccfHW2s9fumuFK58Yz0/7EkFICrQq2kuVABw7+UdeX1aPADumv+JDdF+Tw06BZNB+7Phazbaw9aWE1m8/ONB1w88ZC5ED4Af/gp5qa4/vxBCCDu3hi1FUfTAW8AEoBtwg6Io3dw5ZkvmbTKw+sERAJzK1EpAvLjiID/uS6u1oObiHcnsSc7lySVaK5mmrg3W2gV4ebi908A1faMBrQREJV+zgYJSbc3WnZ9s461fj3E0vcC1A+v02u1ESyksmyu3E4UQwo3cPbM1EDiqqupxVVXLgM+BKW4es0ULsTWQzsivvquttpIQZxe2jAqUsNWUAryMbg9bl3cJ4y9D23PDwBj7az4mAyfOFFJhVfExGwD46YAbmlaHdIQxT8GRVbDjE9efXwghBOD+sBUFJDo8T7K9Zqcoyu2KomxVFGVrRkbj+wZeKEwGPX5mA2cKSqut2+r/7E81jk1yqLV0Ze82+HtKMdOmULkJwajXub2tk6IoPDG5Gy/8uZf9tbyScg6nFfDe2mP2wL29gZsq6m3g7RA7DFY+Cjmn3DOGEEK0cs2+QF5V1fdVVe2vqmr/0NDQ5r6cJlFYVsGCTSd5bvmBaq87znZVWNVqtZb+eWWrvfva5JbdO5Tv7x0KVNXgCm/CXaCH07RbhmsPZ5CWp/2Z2FFL4VOX0OlgyluACkvuBmvDa8AJIYQ4N3eHrWQgxuF5tO21Vq1yRut/G04Q7O2Br+1W0Q6H4pY5RWVYrCpXdI/gun7RBDdDxfvWKtTXRI8orf+kj8nAC3/uycJbBjbZ+D1tYxt02v+eHcN8yMgvJb/ETbW3AtvB+OfgxFrY+qF7xhBCiFbM3WHrD6CToijtFUXxAKYBS908Zot346Vt7Y/HdY9gy2NjUBQ4YCtcCpBdpK3hmtgrkpev6y0lH5rRDQPbEhPUdDtBF9iCXeWi+N62jgfJ7qgqX6nvLOg4BlY/CZnH3DeOEEK0Qm4NW6qqWoB7gB+BA8CXqqruc+eYF4Lnr+7JgNhAAPq2DcDTQ0+7IC8Op1WFrcwCLWwFuXmBtmh5grw96NcukNN5JQDEx2gzXcmu7pfoSFHgqje0oqeL7wJrxfk/I4QQol7cvmZLVdUfVFXtrKpqB1VVn3P3eBeKJyZ3Y95V3fmzbet/53BfDp7Os79fObMV6C2L4lujg6lVfxbiY7RgnphVhKqqvLTyINtOZrl+UL82MOFlSPwdNr3l+vMLIUQr1ewL5FurXtEBzBoSi16n2J77cyyjkCxbCYisQm19TpCs1WqVHNsCdY7wIdzPxKbjmWw9mc3ba47xwBe73DNwr6nQZTL88iyku6GYqhBCtEIStlqIwR20v1zXHz3DnqRcXvvpMIDb6zyJlunf0+KZ0CMCk0GHyaBnYs9Ifj2YwQfrjgNgdVcRUkWBya+ByQe+nAG/vQyJW9wzlhBCtBJKbZXLm0v//v3VrVtrb11zsbNUWBn+0q+k5JbYX/M06jnwzBXNeFWipUjMKmLUK2uqVZof1imE/80egFHvhn8zrX0FfnlGe2wwwazvIabpdmQKIcSFQFGUbaqq9j/fcTKz1UIY9DrmTelx1muyA1FoYoK8GNFZq0PXr522hmvdkTPsT8k718caT1EA258/Syn8+rz2qxBCiAaTsNWCjO0WXu25RC3h6P+u7cW70/ty7+Ud7a9tP+WmyvKxw8BgBkUPig6O/wrvDoWEDe4ZTwghLmIStlowqa0lHIX4mLiiRyQR/lXV7LefclNl+ZiBMGspXP443PIj3PgVWEpg/kRYfDcUuWE3pBBCXKQMzX0BorpALyPZRW6qFC4uChEOrYPc1jMRtMDluE4rdjP89n+w6U04vALGPQe9p9luOQohhKiLzGy1MEvuHsq70/sBcOfIDs18NaIlcmxInpxT7N7K8o48vGDsPLhjLQR1gMVzYMGVcOZI04wvhBAXKNmNKMQFKCm7iMyCMqa8tYHnr+5ZrQVUwplCfM0Ggn1M7rsAqxW2z4fVT4GlGIY9BEMf0HYuCiFEKyG7EYW4iEUHetEr2p8gbw/2JFet23r8uz2MfGUNV76xnuIyN7bc0emg/y1wzx/Q9SpY8wK8M0RrZi2EEKIaCVtCXKAURSHUx2TvOrA/JY9PN58CICW3hOV7Ut1/Eb7hcO2HMP1bsFq024rfzYHCM+4fWwghLhAStoS4gAV6G+1ha8NRLeD88fgY4kK8+WDdcSociqCm55fUeg6X6Dga7vpdu5245yt4sz/s+ARa0DIFIYRoLhK2hLiABXubyLSFraPpBQR7exDqa2LOyA4cPJ3PvpRcAFbsSWXgcz/zR4IbSzYYPWH0kzBnPYR2gSV3w/xJkHHIfWMKIcQFQMKWEBewIG8PsgrLyCkq44utiXQI8wG0xuYAp7KKAFi8MxmAExmF7r+osK4w+we46g1I2wfvXAa/PAflbpxZE0KIFkzClhAXsEBvD3KLy+1rtQbHBQMQE+gFwLPfH6DCqpJwRgtdOcVllFrcuHC+kk4HfWfCPVuhxzWw9iV4ZzAc+9X9YwshRAsjYUuIC1iYrwlVhWW7UmgX7MUDYzsD4G3S6hWfziuhw2M/cCgtH4DnfzjIiJfWNN0F+oTCn9+DmUu05x//Cb65DQoymu4ahBCimUnYEuICdlnHEAAOns63N6iuNPKS0Fo/czqvGW7nxY2EOzfBiL/Dvu+0BfTbFsDJ32Hdq5C4pemvSQghmoiELSEuYO1DvOkcrq3TOjtsfThrgP3xyvuHER3o6ZZrqLCqHE0vOP+BRjOMegzu3AjhPWDZXJg/AX55FhZcJYFLCHHRkrAlxAVuXLcIoGbY0uuqehZ2ifBjfPcI+/Myi9Vl4z/+3R7G/Os3ewmK8wrtDLO/14qhqlbty1IMm96WRfRCiIuShC0hLnC3DYvj1et6c0m4b4331j08is2PjQbg8Yld+VN8GwA+WH/cZeN//kciAFmFpfX/kKLAkHvBYAYU7Wv/d/CvLrDyMcg47LLrE0KI5iZhS4gLnL+XkWv6RaMoSo33YoK8CPczA6DTKYzuGg7ASytdU/vKsbdqdlF5wz4cMxBmLYPRT8DNK2HmUmg/Ara8B28NgI8mwZ6vwdKAECeEEC2QobkvQAjRdDrXMvvljPT8qiCUU1TO51tO0T7Em0ttJSjOK2ag9lUpbgQUpGvV57fNh2/+Al7BEH8T9JsNwR1cev1CCNEUZGZLiFbkkghfbhjYFl+Ta/6d5bgwPruojEe+3cP17//u3El9wmDYgzB3p9Zzse1g2PQWvNFXW0i/bzFUNHAWTQghmpHMbAnRyoT5msgvtVBhVastom+M4xlVYevQ6XxnL606nU7rudhxNOSlarNd2xfAV7PAOwz6TId+syAw1rXjCiGEi8nMlhCtjL+nEYC8Yudnh9LySqnMax+uP+H0+erkFwkj/gb37YIbv4KofrDh3/B6PHz8ZzjwPVRY3De+EEI4QWa2hGhlKsNWbnE5gd4eTp3rTEEpIT4mMgpKqVwr38bf7Owl1k2nh87jtK/cJNj+MWxfCF/cBL6R0GeG1iYoIMZ91yCEEA0kM1tCtDKOYctZGfla2Pph7jD+N7s/A2OD0OuduzVZb/7RMOpRuH8PTFukFUpd+zK83gsWXQ+HVoK1CfpACiHEecjMlhCtTICXFrZyXBG2CkoJ9TXRNdKPrpF+/HowgyPpjVu7papqreUrzktvgC6TtK/sk9pM146P4fBK8IuGDpeDdzBcMrH6zkchhGgiEraEaGXaBGhtexKzipw+15n8UjqFVZWT8DEbKCi1NDg4ZRaU0u/Zn3jjhj5c2btN4y8osJ1Wt2vkI5TvX45x0+uwY6H23vrXIDIeogdASGcI6aT96tdGK7IqhBBuImFLiFYmws+MyaAj4UyhU+dRVZUzBWWE+prsr/mYDJRXqJRarJiN+nqfKyFTC36PfbvHubBls/lkHtd/auSzroMZrOzUWgKhQGE67P4CSvOqDvbwgeCOEHpJVQAL6QxBcWAw1TmGEELUl4QtIVoZnU6hbZAXH6w/wd2jOjZ6kXxesYWyCishPlWf9zVrP1IKSi3VwlZ+STlbT2YzsnNorTNeeSXaLc38UtfsKDxuC5L/ORbBYA8TVJSB3gOuW6DNbBWkwZnDtq8j2q8nN2pBrJKi08pKhJwVwkI6gVeQS65TCNE6SNgSohXq3saPI+kFrD96ptEzSRkFWtNox5mtQC8teK3Yk8qMwbH2129dsJXNJ7J468a+TOoVyW+HMyizWBnbTWsflFNUzybW9VRQooW2TWUdyLrxK4IytlDRdij6yjVbvhHaV/vh1T9YWgCZR6sCWGUYO/YLVDi0DfIKqQpeoZdoj8uL4MxRaD9M1oYJIaqRsCVEK/Tkld1ZvDOFMwW19x08nVtCck4x/doF1nmOjHwtIDmGrQGx2ozPE0v20TM6gPiYAKxWlV1JOQAs2ZlM53AfZv1vCwAJL04CILuwarG+pcKKQe/cRmnHGbJ1JXF0iOvF5DfWM//mdEZeElb3B00+0CZe+3JkrYCckzVD2MHvtUKrjgyeMGupBC4hhJ2ELSFaoQBPI3qdQkZ+7WHruvc2kphVzOFnJ+BhqD34ZNiCWqhPVdiK8DcTHehJUnYxK/eeJj4mgNUH0igptwJwNKOA6R9uth9fXFaBp4e+2s7IvBILQU7W/yooseDlocegU9h0LJPsQi0YPvTlLrb+Y0zDdz3q9NoarqA46Dy++nuFmfDLM1ovR1SwlEDCOglbQgg7qbMlRCuk0ykEe3vUObOVmFUMwC8H09iakFXrMSv2pKIoEOZbvYjpb38bRfc2fizblcLNH23hH4v3AjCmaxhJWcXVAt4xW7sfx9uIrrilWFBajr+nkR5R/hxKyyc1V7vlmVlYxp7kXKfPX413MMTfCAYzoAAqlBe7dgwhxAVNwpYQrVSIj4kzBTWDTYVVtT+e88l2rn13E6qqVjvmcFo+K/ae5tah7fG31e2qpNcpRPp7kpxTzK+HMsjIL2XW4HaM6BxKWYUVqwq3XNYeqOqneDyjamekK4qt5pdY8DEZiA3x5tDpfN5be9xezHXdkTNOn7+GmIHarcNRj2vlJdb/GxLWu34cIcQFScKWEK1UqK+J9PwSCkst3PXpNvanaOUQUnNrzsr8sOc0H204YQ9d3+9ORafAnSM71nrucL/qJRMGtg+mfYiP/fnVfaLw0Os4nJbP3uRc1h89w+gu2lqqx77bS1GZhZLyCuIeXc4325Ia/L0VlFrwMRuIC/GmqEyrIj+kQzCdw334o46ZOqfFDNT6N85cAkHt4YvpkHnMPWMJIS4oEraEaKVigjw5mVnE2sMZ/LDnNFPeWo+qqqTlldQ49u5F25m3bD/fbk8GYMepbC6J8KtzbZW1+kQYfdsFMKRDMA+N7cxzV/egZ7Q/caHebDyWyX/XHcfXZGDu6E4AHEjN46MNCWTkl2JV4Yklexv8vVXObPW1LfCPDvTk39PiiQ32JjWn5vfnUp4BcOMXWumIT6+DIjeFOyHEBUPClhCtVFyID/klFuZvTACgvEJl68ls0vJqX8cFsOGYdgtud1Iu8TEBdR7naauxZTbquK5fNJH+nuh0CveO7sRNl7YDYFy3cPYk57JkZwqTe0fSK9qfp6d0p42/mV8PplNg21FYOTNVm18PpddaCT+/pBxfs4G+bQNZcd8wvpozGJNBT6ivyb6wPyWnmFOZzlfRr1VQnNavMTcRvpwJFteWthBCXFgkbAnRSsWFegOw+UQWk3pFAtqsUuXM1j8mdSXEYadht0g/DqTmU1RmIbe4nLZBXnWe+/6xnXhsYhf2zbuCl6/rXesxtw2Psz++vEs4iqIwc3AsE3tGsic5l5yic6/dslRYufmjP5jy1oZqr2uzc6X2hftdI/2I9NdaFIX6msgqLKO8wsqQF39h+Mu/2j+38egZPt188pxjNkjbQTDlLW1n4vcPwFnr3oQQrYeELSFaqXbB3vbH0y9th07Reh2m5ZVi1Cv8ZWh7Ft89xH7M8M6hHE3Pt+8mDPapuzyDn9nI7cM7oNfVXWLB11y1sL5nlL/9cYcwH0otVg6ermqpk1nLrsnkHG1tWVZhGVPe2kCZRSsvkVtcTkGphehAzxqfqawJllVYc6bpgS938vh3e9l2MrvOa26wXlNhxN9h5yew4d+uO68Q4oIiYUuIVirMoRhpuJ+JYB/tFltqbjHhfmYURbGHkztHdiAq0JPyCtW+gzDkHGGroRwX1MeFaCFwV2KO/bXDaQU1PnPCobfjrsQc0vJKWHMond8OZwDUGrYqZ+ocy0+UlGu3Kf1s4e/345mN/j5qNfJR6HEt/PQU7F/i2nMLIS4IUtRUiFbK21T1v3+Ir4kQHxMZ+aVkFpYRE6jdIjQZ9PbCpqv2nQbg9o+3ARDs7XyT5pev7UVGQWm1IqOdwn0BWLwzxf7a0fR8BncIrvbZE2c10p7x4WZ7Q2uA6MCatzkrw+Nxh8+ezi0hNsSbsgpb4dX0msHOKYqi3U7MOQXf3gH+0RDVz7VjCCFaNJnZEkLgazJoi8fzS0nMKq62Hquygny4X/Xipee6jVhf1/WP4a6zykcEeXvQxr/6WLXVAzuQmlftecJZi91rvY1om9nak1Q1a5ZiK3WRaRvjSHp+fS+//oxmbcG8Tyh8dgPkJLp+DCFEiyVhSwiBoihEBXiyJzmXMwWltA2uOSsUcVYAclw872rvz+xvf+xrMpBXUnOx/L6UPIZ1CuG2Ye3trw209Wb09tDbi5g6qpzZ+iOhal1WZkEZJeUV9t2PKe4qDeETCjd+pVWX/2walLoh1AkhWiQJW0K0Yi9d24uHxnYG4KZL29rrY7WrJWwF22pqjekazpK7L8NsK+/gDj2i/Fl026U8MqELfp7GGlXl80vKOXQ6n55R/lw/IMb++thu4QDoFKXW/odmox5fs4GdDuvBsovK7G2Lwv1MZBeVYbHdUnS5sC4wdQGkH4Cv/6I1uBZCXPQkbAnRik3tH8O9tmKi3dv42V/vEuFb41iDXscfj4/hnel96X2OGluuMqRDCHNGdMDf00jeWWFr3ZEzWKwqIy8Jw8dUNYM1rrsWtvJts1S1qZzdqtwgkF1Ybi/WenmXcFRV2624fHcqTy3d59LvCYAOl8OkV+DIj/Dj464/vxCixZGwJYQAqDYTFOtQFsJRqK8Jo75pf2z4eRqqzWxtPHaGuz7djkGn0LdtAD7mqoX+7YK9eWxiFz6aPaDO81Wu2+oU7oOf2UB2URkbj50hPiaAYZ1CAG2N2N2LttsLvrpc/1tg0N2w+R3Y8l/3jCGEaDFkN6IQwu79Gf04kJqPoYkD1bn4exo54lD64Y6F2m5IRdFm27zPquV1+/AO5zxfVKAnnIA2/p4kZReTXVTGycwiBncItq9DO+NQ16u8wuqegDnuGcg6Div+DoHtodMY148hhGgRWs5PVCFEsxvXPYL7xnRq7suoJq/YwvEzhXy1VdvB5+mhrRUrr9AWmNW2NutcIm0L/UN8TQR4eZCaU0Jqbgntg73ttcPSHepw5ZfUfUvSKTo9XPMBhHeDr2ZD2n73jCOEaHYStoQQLVrncB8A1h/V+jJW3jb0dagT9tI1vVh6z2X1Ol/lLJWqaov+tyRojaJjQ7yJDvTCoFM4llE1k3b2ejGXMvnADV+Ahzcsuh4K0t03lhCi2UjYEkK0aI9O7AqAl21Gq8A20/T5HYPsx0wdEEOv6Pot2u/RRmsN1L2NH90iqzYFDIoLxsOgw9/TyDtrjtlff/nHQ859A+fjHwU3fg5FZ7QaXOXF7h1PCNHkJGwJIVo0s1FPpzAfsgvLKbVUkFFQyn2jO9G9jf/5P1yLMd3CWfXAcK7s3Ya+7bSAFhXgad+l6O9VvT7X8j2pzn0D9dGmD/z5v5C8DRbfBVY3lZ4QQjQLCVtCiBYv0MvDvpBdVSEutPbdkvXV2dYS6LKOIdx7eUe+ubOq4fb/ZtXcyfj0sv389atdTo15Xl0nw9h5sO9bWPO8e8cSQjQpCVtCiBYvwMtITlE5xzO0noZxIT4uOa/JoOehcZdUq44fG1IzyP1vwwm+3pZEcZmbi5AOmQt9Z8Lal2HnZ+4dSwjRZCRsCSFavMqZrcrm07EhNSvcu9JXcwbzzJTuBHlX7/9Y2UfRbRQFJv0L2g+HpfdCwgb3jieEaBIStoQQLV6IrwdZhWWcyirC12TA11yz76ErDYgNYsbgWP41tXe111Pd1TfRkd4IUxdCYCx8cRNkHjvvR4QQLZuELSFEi9c/NgiLVWXF3rvPcGcAACAASURBVFSCfTzO/wEXObuZdUpOE+0U9AyEm74EFFg0FYqzz/sRIUTLJWFLCNHiDWofjLeHnpyicoJtVd6bgmPY8jDoOJSW32RjExQH0z6FnFPwxQywlDXd2EIIl3IqbCmK8rKiKAcVRdmtKMp3iqIEOLz3qKIoRxVFOaQoynjnL1UI0Vp5euiZMTgWAI8mbCUU7K0Fu94xAfRrG8jGY5lNNjYA7YbAVW9CwjpY/qBWiVUIccFx9qfWaqCHqqq9gMPAowCKonQDpgHdgSuAtxVF0Ts5lhCiFRvSIRiAVHcvUnfg72Xkk79cyqJbL2VUl1AOpOax+Xj1wPWfn4/UeM2lel8Pwx+GHR/DhtfdN44Qwm2cakStquoqh6e/A9faHk8BPldVtRQ4oSjKUWAgsMmZ8YQQrVdvW4X46ED37kQ829BOIQBMH9SOj38/yf1f7GR89wiiAz1ZtT+NLSe0dj8JL04CQFVVcorKCfR24dqyUY9B1jH46Z+QthcG3g4xA113fiGEW7lyPv4WYIXtcRSQ6PBeku01IYRoFH8vI4tuvZT/3NCnWcb38jDwzk39SM0tYf7GBJ5dfsAetBz9b0MCfZ5ZTVJ2kesGVxTod7P2656vYMGVkLjFdecXQrjVecOWoig/KYqyt5avKQ7HPA5YgE8begGKotyuKMpWRVG2ZmRkNPTjQohWZEjHkBq1r5pSjyh/Ft9ds+F15UL6mf/bwjPf7wcgPb/UtYMnbQEU7bGlVFvHJYS4IJz3NqKqqmPO9b6iKLOBycBoVbWv3kwGYhwOi7a9Vtv53wfeB+jfv7+s/hRCtGjxMQGM7hLGzwfT7a8Fe3tQUl7B2sNV/2AsKXdxtfnYYaA3gaUEUCG4o2vPL4RwG2d3I14BPAxcpaqq45z5UmCaoigmRVHaA50AmfMWQlwUogM9ARgUF4SHQcfxM4Xc/NEf1Y7JK7a4dtCYgTBrKQy9H/QecPAH155fCOE2zq7ZehPwBVYrirJTUZR3AVRV3Qd8CewHVgJ3q6rq5qZiQgjRNHpE+dsfj+0WDsCms3Yk5pWUu37gmIEw5ikYdBfs/gLS9rl+DCGEyzkVtlRV7aiqaoyqqvG2rzkO7z2nqmoHVVUvUVV1xbnOI4QQF5Kr+0QxbUAM917eibziqlB1afsg7hgRB1DtdZcbej+Y/ODnZ9w3hhDCZaSCvBBCNJBBr+PFa3pxWccQcm2halinEBb+ZSB/H98FRYG8EhffRnTkGQhD74PDK+DUZveNI4T4//buO77q6v7j+OvkJrmZZC8gYQ/Zsh2gAioqilp3qVZttda6a6u2Vmut7a91d7i3VuusuAUXONl7g4yETLL3+v7++N7cJCSQdW/ujXk/Hw8e3Pu9937PwduSN2d8jkcobImIdEFKVAgA95w1Fmegg4AAQ0igg037iwDYnFnMvnwPloFoMO0XEJEEi+9UZXkRP6ewJSLSBX/70Xge/8kkUmMbi62GBAWweHMOv3hhJac8tJQrX1jpfq28upaKag8sYQ0Oh5k3w96vYcfirt9PRLxGYUtEpAuiwoI4aXRys2vXzxkOwIcbswAoKG88RHrcnR8z96Elnml84iUQMxA++SPU13vmniLicQpbIiIedtzwhGbPAx12MdK80ipq6y32HCjn7TUZ/Oe7vV1rKDAYTvgdZK2HjW927V4i4jUKWyIiHpbsWsfVIL/UHtn6tEkh1OteWcNtb62nqLyLuxbHnANJY+DTu6HOizsgRaTTFLZERDwsJMjhftwvOpSy6jp25JTwm9fXtXjvku1dPKYsIABm/wEKvofVL3TtXiLiFQpbIiJeEBJk//U6NDECgJeX7QNgSEJ4s/dtzizuemPDToLU6fD5/0G1F3Y+ikiXKGyJiHhB32j7SJ/hSXbY+t/qDNJiw7hm1jD3ewYnhLMps5g1+wqpr+9C+QZj7MrypVmw7LEu9FpEvEFhS0TEC6YPjgNg0oAYAA6UVXP8iATiI5wARDoDGdM3is+35nLmv77iuW92d63BAUfBsJPhywegoqBr9xIRj1LYEhHxgj/MG8U/LzqSk5uUhbj4qIFMHRTLbaeO5OUrppPWpDbXrtyyrjc6+3aoLIKvHu76vUTEYwJ93QERkR+ikCAH88b1BeCB88cT7HC4129dMXMIAJv2N67XCnd64K/j5LEw9lz49hGYdiVEJrf9GRHxOo1siYh42VlH9ue0cSktrjetOl9T56GipCfcBvU1sOTvnrmfiHSZwpaIiI9MGhDD9MGxABSUVbfx7naKHWxXll/5LOTv8sw9RaRLFLZERHwkODCAV644irH9opod6dNlx/0GAoLgs3s8d08R6TSFLRERH4sOCyLfUyNbYK/Vmn4VrH/dPspHRHxKYUtExMdSY8PYlVvWtVpbBzvmWgjpA5/8yXP3FJFOUdgSEfGxSWkxlFTVsj2n1HM3DY2BY2+A7R/Bnm88d18R6TCFLRERHxufGg3A++szsSwPjm5NvRIikmHxneDJ+4pIhyhsiYj4WGqsfbTPQ59s57ElHtxBGBxmL5bf9y1s/9hz9xWRDlHYEhHxMWegw/34tRX7PHvziRdDzCD45C6o91AtLxHpEIUtERE/kl5QwdNffu+53YmOIJj1e8jeABve8Mw9RaRDFLZERPzAb+aOwBkYQFVtPXe9u4lb31znuZuPPhuSxsJnd0OtB0tMiEi7KGyJiPiBXx4/lDV/OMn9fPXeQs/dPCAA5twBBbth1XOeu6+ItIvCloiInwgNdvDxDTMZmhhBZU2dZ28+dA6kHW2fmVhd5tl7i8hhKWyJiPiR4UmRnDw6ibLqOs+WgTDGHt0qzYbvHvXcfUWkTQpbIiJ+Jiw4kLp6i6paD+8eTJsOw0+BLx+C8nzP3ltEDklhS0TEz0Q4AwEoq6r1/M1n3w5VxfDVQ56/t4i0SmFLRMTPhLvCVnm1h9dtASSNhnHn2VOJxfs9f38RaUFhS0TEz4QH20VO9xwo9+y6rQbH3wr1tfDF3zx/bxFpQWFLRMTPNIxsLXjqO+5YuNHzDcQOgkmXwqrn4cBOz99fRJpR2BIR8TPhzsbje57/Zg/pBeXU1nl4sfzMmyHQCZ/92bP3FZEWFLZERPxMw8hWg2P/7zPuW7TNs41EJsH0q+wjfDLXevbeItKMwpaIiJ8JCwpsce3V5fs8v37r6GshJBo++ZNn7ysizShsiYj4mf4xofzqhKG8cPlU97UDZdV8uCHLsw2FRsOMG2HHItj9pWfvLSJuClsiIn4mIMDw65NHcPSQeC6Yksq71xxLhDOQb3cd8HxjU6+AyBRY/Efwxs5HEVHYEhHxV44Aw19/NI4x/aIYkhDOjtxSzzcSFArH/RbSl8G2Dz1/fxFR2BIR6QmGJEawNauUak8f4QNw5AKIHQKf3AX1XiikKtLLKWyJiPQAJ49OJq+0iv8u3+v5mzuCYNbvIGcTrH/N8/cX6eUUtkREeoCTRycTHuxgV16ZdxoYdRYkj7PrbtVWe6cNkV5KYUtEpIdI7BNCbkmVd24eEABz7oDCvbD4Tlh6H+xb5p22RHqZlsVcRETELyVEOlm5p4DKmjpCghxtf6Cjhsy2R7e+/RcYBziC4ZKFkDq17c+KyCFpZEtEpIcIMJBZVMnd723yTgPGQP/J9mOrDmortYZLxAMUtkREeojiiloAln2f771Gxl8IDidgAAuWPQ6PnwArn4MqL5SeEOkFFLZERHqI+84bD8C27FI+3ZLtnUZSp8JP34XZt8OCN+CUv9kjXO9cC/eNgHeuh/1rvNO2yA+U8fhZW10wefJka8WKFb7uhoiI37rhv2t4a3UGxw1P4LnLumktlWVB+nJY+SxseBNqKyBlAkz6KYw9B5yR3dMPET9jjFlpWdbktt6nkS0RkR7ktlOPACAkqBv/+jbGHvE6899w0xY49V6or4V3r4d7R8DCayBjpY77ETkE7UYUEelBEiKdzBgWT1ZRpW86EBoNU38OU34GGatg5TOw/nVY9Twkj3WNdp0LIVG+6Z+IH9LIlohID5PUJ4SsYh+FrQbGQP9JMP+fcNNWOO1++/p7N8F9I+F/V8O+5RrtEkEjWyIiPU7fKLu4qdfqbXVUSB+YcjlMvgz2r7bXdq1/Hda8CImj7dGucefZo2IivZBGtkREephx/aOpt2DNvkKP3G9HTike2SxlDPSbCGc8DL/eCqc/BIHB8MHN9mjXW1fB3u802iW9jsKWiEgPM2VQLI4Awz8+3d7lkLR8dz5z7v+CV5bvc18rLK/uevhyRtojWld8Dld8ARMuhM3vwNMnwb+nwwe3wKd/1pFA0isobImI9DBRoUH8eFoaX+04QGF5TZfutS27BGgslLpidz4T7lrE++uzutxPt74TYN4D9k7GM/5hj2x99wgs+Rs8N0+BS37wFLZERHqgSQNiAMgr7drB1Hkl1QDU1dsjWY8v2QXYU4se54yAiRfD+PPBuH781FbBmpc835aIH1HYEhHpgRIinADkdjFs7T5QBkBGYQUARRX2SJmFxf7CCj7dku2Z9VxNDZxhHwlkHICB1S/BpoWebUPEjyhsiYj0QPGRdtjKK61u8VpHwlFGgR2y9uWXk1Nc6R4pK6qo4XdvreeyZ1ewcO1+D/S4idSpcMlCmPU7uOg1e5rx1Yth2ROebUfETyhsiYj0QPGuka3XV6Y3u/7GynQG3fo+OSXtq8OVWWyHrZySKqbe8wk7c+2Rrl25ZSzZngfA5swST3W7UepUmHETDD8RLl4Iw0+G938Nn9yl3Yryg6OwJSLSA0WHBgGwZFsumzOL3ddf+HYPALvzytu8R329RVZRJQmuUbKmvtiWi8MYADKLKjzR5UMLDoPzX7LXcy29D97+FdR1beG/iD9R2BIR6YECAgxPXGyff/vuusZpvtr6egBW7y1o8x6r9xVQU2cx2bXY/mCzj0hk6sBYMgu7oVq9IxBOfxiOu8UuhvryhVBd5v12RbqBR8KWMeYmY4xljIl3PTfGmIeNMTuMMeuMMRM90Y6IiDQ6cVQS/aJDm4Wh2jp7Cu4vH2xh5Z7DB66bX1sHwMmjk93Xbjt1pPvxrJGJpESHuBfPe50xcMKtMO9B2PkJPDsPyvK6p20RL+py2DLGpAInAXubXD4FGOb6dQXwSFfbERGRluIjgskra1wkX1Vb7378omtK8VCqauuJcAZyxvi+3HPWWG6YM5wrZg5hrit8nT6+L2P7RZFRWOGux9UtJl9qTyvmbIKnToT877uvbREv8MTI1gPAb4CmKxrnA89btm+BaGNMigfaEhGRJuIinOSX2TsIN2cW831e49TbW6szWJdeyNX/WcXYOz/i8SU73a9ZlkV+WTUXTEklIMBw0bQ0rpszDIAHzp/AqttPJCTIwVlH9gNg0absbvxTASNPtRfOVxTYgWv/6u5tX8SDuhS2jDHzgQzLstYe9FI/YF+T5+muayIi4kFx4cEccJV/uHPhRqJCg3joggm8/oujiHAG8tSX37NkWy4llbX867Od3P6/DRRV1FBcUUtFTR3JUSEt7hka7CA2PNi+f4STlKgQdnqjyGlb0qbBZR9DYIg9pbjjk+7vg4gHBLb1BmPMYiC5lZd+B9yGPYXYacaYK7CnGklLS+vKrUREep24CCcHSqs5/R9fsj6jiN+fdgTzJ9j/tj1vcirPf7ObWld1+KKKGl74dg+OAMNJo5MASIkKbbONIQkR7Mj1QdgCSBgOly+Cl86B/5wH8/9tV6AX6UHaHNmyLGuOZVljDv4F7AIGAWuNMbuB/sAqY0wykAGkNrlNf9e11u7/uGVZky3LmpyQkNDVP4+ISK8SHxFMdV096zOKANzTfmDvJmwIWmdPbLz+7Ne7ueiJ7wAYEBfWZhtDEyPYmVPq+Ury7dUnBS59H9KOgreugC8fVC0u6VE6PY1oWdZ6y7ISLcsaaFnWQOypwomWZWUBC4GLXbsSpwNFlmVleqbLIiLSoOk04KMLJhEX0Vgzq19046jVzGEJfHj9jGahyxFgGJYU0WYbQxIjKKuuI6u4G0pAHEpIFCx4A0afDYvvgA9vgfo63/VHpAO8VWfrfeyRrx3AE8AvvdSOiEiv1rdJoBp+UHBKiW4MYvERTkYm9+Gu+WOIdNorSKJCg3AGOtpsY2iCfV+vHE7dEYFO+NFTMP1q+O5ReP1SqPFhABRpJ4+FLdcIV57rsWVZ1tWWZQ2xLGusZVkrPNWOiIg06ttkzdWAuPBmrzUNUmmx9nRhhDOQRxZMAuBHE9u3b2lIon1fnyySP1hAAMy9B066Gza9DS/+CCoKfd0rkcNqc4G8iIj4r6ZH7TgCTIvXbz55BBHOQNKarM06ZmgcL/1sGtMHx7WvjQgnIUEB3VfctD2OvgYikuB/v4RnTrGnGPv09XWvRFqlsCUi0oM5Agy3nDKSI1OjW3396hOGtrhmjOGYofHtbsMYQ9/oUF76bi8pUaHMG5fC2vQiThyV1Ol+e8S48yA8Af67AJ480Q5ciSPb/pxINzM+213SismTJ1srVmjGUUTE35z0wBdsy7anEY2xNwNuu/sUggP94IjdzLXw0rlQWwUXvgIDjvJ1j6SXMMastCxrclvv84P/l4iIiL/Lb3IkUMO/0QvLqw/x7m6WMh4u/xjC4+GFM2HzO77ukUgzClsiItKmf100kauOH8LghMZF+AXlNT7s0UFiBtrV5pPGwKsXw/Infd0jETet2RIRkTZNGxzHtMFxHCitYleuff5igb+MbDUIj4NLFsLrl8F7N0H6SogbCoNmQOpUX/dOejGNbImISLtdduwgRvftA3hvGnFffjlfbs/rXMX64HA4/yUYPhfW/gc+vQueOx32LfN8R8VvbNpfTHFlDb//33que2U1+/1p5ywKWyIi0gEjk/vwxMX2euBfvLiKmrr6Nj/z7a4DDLzlPdILyg/5HsuyWLo9lw0ZRdyxcCMLnvqOhWv3k1nUiR+ajkDoPxVwlcKorYR3rofdX+qYnx6oqrbukMG+sqaON1amc+rDS7n4qWW8+O1e3l6zn/sXbevmXh6ephFFRKRDYsKC3Y9X7y1k6qDYw77/yaXfA7BmXyH9Y1o/i/HxJbv4ywdbml277pU1AOy859RWa4gd1qAZEBgCddX29snC3fDsaRA/AiZfBuMvgNDWy2WIf/nt6+v435r9bLrrZMKCm8eW+xdt4/EluwD7f18NNLIlIiI9WmiwgzlH2DW2vtl5oMXrW7NKuPSZZZRX1wJQXGkvpK+rP/So0hfbcg/52sb9RR3vZOpUe/3WrN/BpR/AzTth/r/BGQkf/hbuGwlvXw0ZKzt+b+lWS7fnAfDGqowWr23aX9zs+bWzh3HG+L6kFyhsiYhID/fkJZOJjwgmu6Tl2YR/encTn23NZck2+4dkwxRQbknVIe+XWXToMw4/2pjVuU6mToUZN9m/B4XCkT+Gn38CVy6B8efDhrfgiVnw2HGw6nmoLutcO+JVQ1xnc765Kr3Fa8bYZ3w+umAiRw2O4xfHDSY1NpSMwgpq2zHF3V0UtkREpFNiwoIpKGu5lsa4Zvw27S9id16ZuxhqbmnrYcuyLDIKK5gxLJ4Zw+zK9g2L8AHeWpXRucXyh5IyHk5/CG7aDKfea081LrwG7jsC3v8N5Gxp+x7SbRp2va5PL6Kqtg6Aoooa6ustsooqmT44lrljUnj5iumEBQdyREof6uotvt2V78tuN6OwJSIinRITHtys2GmDPQfshfC78sqaraPJK2l9kfPXOw9QXVvPnCOSePqnU/j9aUfwzKVT+OnRA7lm1lD2F1Uy+e7F1B9mGrJTQqJg6s/hqq/hso9g+Mmw8hn49zR45lRY/7pdlV58qqC8hkhnILX1FrtyyyivrmX8Hz/m/z7aQlZxJcl9Qpq9/8RRScSFB/NGKyNhvqKwJSIinRIbFtyi1lZJZQ178+2wdaC0mjzXaFa/6NBDjmzd+/FWAEYkRxLkCOBnMwaTGBnCnWeM5uKjBtr3Kqvmm10t14d5hDGQNh1+9ATcuAVOvAuK98Mbl8P9o2DxnVCw2ztty2FZlkVheTXTBtubMD7emO1eJ/jYF7soqawlNbb5pgtnoIPnL5/KX84e2+39PRTtRhQRkU6JCQ8mf0/zKvJbs0rcj7/ZdcAdkAYnhJN3iDVb6QUVzDkikemD41q8lhDp5OtbZnH0Xz9l1Z6CDh2g3SnhcXDMdXDUNbDrM1jxNHz1MHz5IAydY+9kHH4yBDi82w8BoLSqltp6i8kDYymqqOGBxS1LOoztF9Xi2ui+La/5kka2RESkU2LDgygor262nmprth22RqU0rrkKCQogJSrEPcoFkOHaml9cWUNuSRWTBx66fETf6FD6hAQ2+7zXBQTA0NlwwUtw/Xo47reQvQFeuRAeHAdf/B1KOrlwX9qtYZo6LjyY3582qtX3jGklbPkbjWyJiEinxIQFU1dvUVxZS1RoEAVl1dz/8TaCAwMYmhjBpkx7W36fkCASIp3klFRxzcurmTksnptfX8czl04hPtwJwOD48MM1RUKkk7zSaqpr6znu759x7exhXDg1zet/RgCi+sEJt8LMm2HbB7D8KfjsbvjirzDyNBhwNJTmQuo06DfRo01bloVlQUDmati/Goac0KuOHspy7VJNjgphfGo0jy6YxNasEt5cnc7PZgxmSEI44U7/jzL+30MREfFLDcVNC8qqiQoN4tpXVnPANRJhmtQg7RMaRHyEHareWbufd9buB+zdZeNT7cKicRHBHE58hJPckirWpheSWVTJrW+uZ2JaDO+t288NJw7HmA4WPe0MRyAccbr968BOezH9imdh09tea9LgroNvW/I3+PFrMGSW19r0J1nFdthKibIXwc8dk8zcMclcN2eYL7vVYQpbIiLSKbHhdkDKL69mIOHu4pMAtXWNU4thwQ532DpYQw2uqNDDh62ESCcb9xfzXZNF8uc++jXFlbVcPmMwUaFBnf5zdErcEDjpbgiOhM//CtQDxj6TcehsjzRRV1/Pne9sYmbAOuY4VmOwoL4WXr7QntacdqV9FuQPWMPIVtJBOw57GoUtERHplJjwxpGtpmckXjt7GFuzGit719ZZjGpSN6tBdnEl0WF2SGr4/VASI0P4pCiH9RmN1eSLK10V6itquj9sNRhyAnz5gF2ryxEMM2702DTfb19by+t1/dlYP5ATgjYSaNVCQCAkj4NP/gjf/huOvcFetB8U6pE2/U1mUSURzkAiQ3z0/XqIFsiLiEinxLqmEfPLqslx7TT869ljufHE4c1GtuZP6NvqmqzMokoKy+3djG2FpcEJ4VTU1PHplhxGJkc2e62oouYQn+oGTY8FumShx4JWeXUtr6+060Stsobz6IAH7DZ++i78bBFcvggSR8FHt8HDR8KyJ6C29TpmPVVpVS1vrkpnQmrPP8NSI1siItIpMeF2QCoor26c7nGtrWmYNnzliunukg7vXnMs1XX1LNqUzfLv89lfWMHAuHAinIEEOQ7/b/9hifaRLTV1FnPHJJNRWEFJk5Etn0qd6vFF6w1n+/3jwiN5bWU6926G42bPY2z/qMY2L1kI3y+Fz/4M7//aLlFx3G9g/IX2+rIerLy6luP+9hnFlbVcduxAX3enyzSyJSIinRLhDCTYEUBeaTXZroXMSZF22Lr99FH85eyxTBvUWNJhTL8oJqbF8Nu5IzkipY89slVR3a4pwGFJjaNZA+PCmTkswf3cpyNbXtIw4hcTFsz88X0B+GZXXss3DpphH7S94A27RtjCX8G/psC6V6G+rju77FFLt+e5N1sMS4xs493+T2FLREQ6xRhDWlwYjy/ZxS9fWgU07hqLcAZy4dS0Q+4STIkOoaiihk37i90jZIfTsBgfoF9MKJMHxrifX/XSqmZrxn4IGjYORIcFcfbEfhiDeySvBWPsgqs//wwueBmCwuDNn8MjR9s7Jet73n+bHTml7sd9o3v+ejSFLRER6bQhCc3XYrW10L1B3yj7B+iWrBLOnNCvQ232iw5l4EFrwDILKzt0D39XWNG4ls0YQ4Qz8NBhq4ExMPJUuHIpnPMMWPXw6sXw+EzY+iF48jBvL1uXbp+pOW9cCo6Abijr4WUKWyIi0mkHl3Rob72rpqMVx49IbNdnEiOd7t8HxTUPW/uLKtp1j56iyDWN2BBe+4QEtR22GgQEwJiz4ZffwlmPQVUpvHw+PDkHdn7aI0LXhoxi5o1L4Z8XebZIrK8obImISKddM2sYf5o/usOfG54U4X6cGtu+aaL/XX0Mz146hUBHAP1imn9mf+EPK2wVVlTjCLBHtAAiQwIpqezg2rQAB4y/AH61HE5/2D5e6IWz4NnTYPdXXui1Z+SXVZNRWNHqmYc9lcKWiIh0WnJUCD85aiB3nzmGt68+pt2fiw5rXIPlDGzfoc59o0Pdo2BBjgCW3HwCq28/EYB16UWH+yi1dfV8vDGr2TmO/qywvIZo1xQiNIStdo5sHcwRBJMugWtXwan32tXvnz0Vnj8T0ld4sNeesedAGQBDEyPaeGfPobAlIiJdtmD6APfRO+115oS+nDImudNtpsWFERMezPjUaP7z3V6KDzPy8/jSXVzxwkoWbcpudr2+3uKdtfuprPGvnXvZxVXNjjCKcAZSUtXFXZeBTpj6c7huDZz0Z8haB0/Ohv+cD5lru9hjz2mo2ZYY2bOrxjelsCUiIj7x4AVH8siCSV2+z1XHDaa6rp49eeWHfM/2bHt324o9Bc2ur95XyDUvr+bnz/t2hGdffjkDb3mPlXsK2JJVzOLN2QyObxzZiezImq22BIXC0b+C69bCrNth7zfw2Ex7Mf2612DpfbBvmWfa6qANGUXusJUQ2foRTz1Rz656JiIivd4A12L53QfKGot+HiTDVST08SW7uGBKKoMTIsgvqyav1P7BvnR7HpU1dYQEtW9K09O+2JYLwNNf886UWgAAFJBJREFUfs976zMBSOrTGDYiQwI9X7zVGQkzfw1TfmYf/fPVP1yHahsIDPFoRfz2WL47n3Mf/YYIZyDGtH04eU+ikS0REenRBsSFAbA3v3FkK72g3L2jr77eYuP+xjVd3+eVsWRbLhP/tIh312W6r//oka+7qcfNVdbUkeFa4N8QtKB5IdeUqBAKymuoqPbCdGdoNJxwGxx9NWAAyz7rcfdSz7fVigcXb+Oe9zezNasEsI/piQkLbvNUgZ7kh/MnERGRXiksOJCESKd7YbVlWRz7f59xzqN2eHp3fSZl1XXceOJwwF4TtHizvXbr8y05AIQHO9iRU0p9ffcuoN9zoIwJd33MI5/vdF8bnxrN4htncuHUNPe11Fg7UKYXHHqqtMuGnWSPaBmHfaj2wBnea6uJBxdv5/Elu9iWXeK+NiKp51eNb0phS0REerwBsWG8uiKdqX9ezAcbsgDY7qpCfu9HWwHci/F35JTy2gr7kOeSqlrCgx3cdtoRVNXW8+SXu7q132vTi6isqWfGsHgumJIKwMxh8QxNjGxWzLN/jB229nkzbHnpUO32WrQpmxnD4lnx+zk8d1n3tu1tWrMlIiI93oC4cFbsKSCnpMp9dBDAja+uIauoknnjUhiWFElMWBCvLt9HRZPdhzHhwe7z9+55fwvnTkolJtz764WKymtYvbcAY+CJiycT5AjgyLRoTnedhdhUmmtka0tWCbNGJnmvU144VLu9Mosq+eMZo1sUyv0h0MiWiIj0eA3rtg725qoMquvq3dNSCZFOSqpqCQ1ycIYr1CT1CWHSgBjOOtI+NujLHa0c+OwFJ9z3Oc98tZt+0aGEBDlwBBjOn5JGWHDLcZCESCeTB8Tw+sr0bumbr4zr37HyIT2FwpaIiPR4DSM/DSYcVPMr3lVG4IiUPq7fI/nbOeN46IIJ/P2ccTgCDH+YNwrAvUPR2/LL7MOm+8e0r4L+zOEJ7Moto6rWv2qCdcV7TTYoQPMdmD8kClsiItLjpTYJW3+YN4pZI5uftxjnmhacP8Eezbr0mEGEBDmYP6EfgxPselZRoUE4AgwHSqu93t+mlexvO/WIdn0mOcou8plT3D1hsDtc/Z/GKd/ESGe7z9bsaRS2RESkx2s6snXZsYMIdzafimsY2Zo1Momvb5nV6rqogABDTFgQB8oaw1ZZVS2WZfHk0l385KnvPNbfAldZij/MG9XuqbMUV9j6oZwDWV1b7358x+mjWHzTcT7sjXcpbImISI8X7yqA2SfEDlk/npbG8012tA1oEsb6Rh962i4u3El+mT1ylF5Qzug7PuLlZfu4+73NLN2e1ywgtFd2cSW3vrmu2ZFAmUV2YGoIUO3R8N6s4soO98Efbc9pLPUQ6AigT0iQD3vjXdqNKCIiPZ4xhhcvn+ZeKB8S5GDm8AT363Ht3OEWGx7MngN2eYWvXAvln/t6t/v17OLKZlOW7fHwJ9t5edk+jkyL4bzJdnmH3E4cSZMSZYfE9IKeP7K1IaOIW95cB8BF09I4d1J/H/fIuxS2RETkB+HYYfEtrn14/YwOjZj0iwnlm10H+GB9Jt99nw/A1ibFNjOLOh62Go4AymkyIlXkOnonOqz9fQt3BhIf0Vi8tSeb948vAbhwair3nDXWx73xPk0jiojID9bI5D6HnTY82B9OH0VabBiPLdnlHtmCxh2D69ILO9R+fb3FzlxXcdWPt3HCvZ+TU1zpPuewT2jHps4GxYex+zAHbvu71XsLWLuv8b/h9MFxPuxN91HYEhERcekTEsTxIxJYs6+Q7Ca7/h5dMAmAez/e2qEjfV78bg+fb7UPmT5lTDLf55WxfHeBe2QrqoNha2BcOLvySpvtZuwp6ustzvr318z/11cAjEyO5PRxLTcq/BApbImIiDSR0GR912ljU5g1MpEx/aK4+eQRVNbU8/w3u7nhv2vaFXh2uo4MAnjg/AkEGNiaVUxRRQ2hQQ6cgY4O9W3igBjySqvZuL+4Q5/zBxuaHAYOcNkxgwgI+GGWejiYwpaIiEgT8U0Wrf9q1lCe/ukUAI4fYS+4v/OdTby1OoPcVoqfrtidzzc7D7ifhzUpQRES5GBQfDibMospLK/p8KgWwMmj7fMdl2zP7fBnfW1rVkmz58kd2InZ0ylsiYiINNF0ZCupT2MgGJ4U2axUw7as0mafK6uq5ZxHv+HCJ751V6HPdxVI/fiGmQBMSI1h1d5CCis6F7Ziw4OJcAb2yMKm2QeVrGio5t8bKGyJiIg00XRkK6bJbsEgRwC3nDLS/XxbdvORmk+25LgfNwSLA2VVjErpw3DX2YxTB8WQX1bN6r2FRHVgJ2JTCZHOZkcKVdbUsakHTCtmFVcSHRbEwxceyaMLJnWo7EVPp7AlIiLSxKC4cMA+3Prg42OGJUa6Hx8ctj5Y33jOX15pNac+tJTFm3OIcxVcBZgyMNb1epV7WrKj4iOC3WGrqLyGUx5ayqkPL6Wksoaiiho+3JDZxh18I7u4iqTIEM4Y35e5Y5J93Z1upTpbIiIiTUSFBbHi93Ooa2XX4eCEcPfjjfuLsSwLYwyWZbFkWy5HD4nj650HyC6qZFOmPdrUdOpxUHw48RHBVFTX8eNpAzrVv/gIJ9uyS9ieXcLZj3xNSWUtALvzynn6q+95a3UGH14/g5HJ/jNNZ1kW+/LLSepF67Sa0siWiIjIQeIjnM3WazUICXKw7e5T+PG0NNZnFPHCt3vYkVPKL15cSVl1nbtq/Zom9biOGtJYS8oYw40njuCOM0Z3as1WQ99yS6p4fWU6ZVW1/PGM0QA8umQn3+6yF+ev2duxemDeVFdv8c9Pd7Alq4QTOjma19MpbImIiHRAcGAAPz16IACr9hTw3+V7+WhjNgATUqMJDgxg8Sb7+YikSE4Zk9Ls8xdNS3Mf29MZR6T0obiylvc3ZDIwPpzzp9j3em9dJplF9lqx11am+00trsWbs7lv0Tb6hAR2ejSvp1PYEhER6aBhSZEcNTiOvfnlzWpeDUuMYMrAGHJcZx/+98rp7uN6POW0sSkEOQz78itIjQkjJMjBpccMbPaelXsKmHT3Yiqq61q/STdq+G/x2E8mExzYO2NH7/xTi4iIdNGAuDD25pezObOYs47sxxtXHUVchJM7Th/tfk90WPBh7tA5UWFBTBoQA0BqrH2M0B2nj2baIHvx/ZwjEgHIL6tuUUjUF0pda8ompEb7uCe+o7AlIiLSCcOTIskrraagvIbhSZFMGmCHnWGJEQQGGM6c4L2jaH40sT9x4cEcNzzRfa2s2g41p41rnLbcX1jhtT60V2lVDY4AQ0hQ740c2o0oIiLSCXPHJHPXu5sASI5qrBlljGHTXXNxePEomnMnp3LuQeu+Ql3TlWP6RjEqpQ+bMov94tDqsqo6IpyBLcpo9Ca9N2aKiIh0Qd/oUPfjg3cuBgcGeDVstea+cydw04nDGZoYwfvXzWBAXBiLN2e3WsLCG/YeKOezrTlkFjUfTSuprCXC2bvHdhS2REREuii5lTIR3S0tLoxrZg9zjyBdOXMI6zOK2JLl/erylmUx8++fcekzy7nhv2uavVZaVUNkiMKWiIiIdMJc18HQ/nio8pSB9iL6gyvde0NpVW2rjxuea2Sri4wx1xhjthhjNhpj/tbk+q3GmB3GmK3GmJO72o6IiIi/efCCCSy6YSZhwf4XJgbG29Xub/jvWiprvFsCoqG8A8DB5b1KK2sJV9jqPGPMCcB8YLxlWaOBe13XRwEXAKOBucC/jTGeLTQiIiLiYyFBDoYlRbb9Rh8IcgSQFhsGwPLd+V5tK6fYDluD4sPJchVWbVBUUUOEphG75Crgr5ZlVQFYltVw5Pl84BXLsqosy/oe2AFM7WJbIiIi0gH3nTcegJ88tYytWd6bTswpsQPW2H5RHCirdhdTzSyqYPeBckb39Z9zGn2hq2FrODDDGPOdMeYLY8wU1/V+wL4m70t3XRMREZFuEhfeWFR1fYb3CpzmuqYRTxljr2F7e00GlTV1fL3DPqtx9sgkr7XdE7Q5rmeMWQwkt/LS71yfjwWmA1OAV40xgzvSAWPMFcAVAGlpaR35qIiIiBxGXHhj/a/deWVeayenpIrgwADmjkmmf0wot7y5nlveXM9F0+yf6wPjw7zWdk/Q5siWZVlzLMsa08qvt7FHrN60bMuAeiAeyACaVlvr77rW2v0ftyxrsmVZkxMSeudp4CIiIt7QJ7RxTOWfn+3g6515Xmknp7iSxEgnxhhGNFnD9uGGLBIinTgDe/ey7a5OI/4POAHAGDMcCAbygIXABcYYpzFmEDAMWNbFtkRERKQDDq7a/tgXu7zSTk5JFYmR9ija0KQI9/X8smr6+mFZjO7W1e0BTwNPG2M2ANXAJZZlWcBGY8yrwCagFrjasizfHz0uIiLSy1w7exgDYsP4eucBlmzPbfF6bkkV2cWVXPz0Mh758USmDY5r970PlFaxcX8xOSVVDE2wQ1bD7w1SokJb+2iv0qWwZVlWNbDgEK/9GfhzV+4vIiIiXXPjicMBKCiv5o1V6Qy+9T0+vmEmQxPt6b6p9yx218a68dW1fHXLLPdnv88rI7OogqOHxLd673s/3srLy+z9cEcPsUPa0MTmYWuyq7hqb6YK8iIiIr3AESl2+YV6C+bcv4TPtuQw677PmxUhzShsfq7hSQ98wUVPfId1cKVSl6bnLk4aYIeqIU3C1rRBsS0OzO6NeneVMRERkV5iZHLz4qv3LdrKrtzWdyhuzy4ht6SKmjo7TKUXVJAa23JHYa3r9emDY5nrKvvQJyQIgHMm9efec8d7rP89mcKWiIhILxAX4Wz2fENG6wdUl1bVctrDX1JdV+++tjmzuNWwVVxZw+i+fXjliqOaXd91z6kEBJgW7++tNI0oIiLSS/z5rDH866KJHD/i0KWWsosrmwUtgOe+2d3igGmwj+JpGMlqSkGrOYUtERGRXuLH0wZw2rgU9/qt1uw9UM7BWemrHQe4c+HGFu8trqhtVstLWqewJSIi0suMcoWtgXEtpwYvfXY59a2sh99/0OJ5sKcRo0JbjmxJc4qjIiIivcy8cSn0iwmlT0ggc+5f0q7PRIc1D1UllTVkFlW2Oo0ozWlkS0REpJcxxjAxLYahiZHs/utp7utTDlMTKzCgeWT45UurABieFNna26UJjWyJiIj0ct/eOpuaunriI5xs2F/EuY9+w8zhCYzvH8W8cX05+cElFFXUuN9fVFHD0u15/PL4IZw3RXW02qKwJSIi0sslNzm/cMrAWJb/bg7xEcHusxWPG55AYXm1+z0N67dG9T30QntppGlEERERaSYh0tnsEOvosCAKyhtHtjKL7LDVN1rnHraHRrZERETksOLCnezNL+eBRdvIKamktKoOgL46ZLpdFLZERETksPpG29OMD32y3X3NEWBIiHQe6iPShKYRRURE5LD6x7QcwXIGBuBQpfh2UdgSERGRw+oXbRc/Terj5NEFEwFQzGo/hS0RERE5rIaRrfMnp5KidVodpjVbIiIiclgx4cGsuv1EokODyCurAiDcqQjRXvovJSIiIm2KDQ8GICHCyZUzBzN/Qj8f96jnUNgSERGRdjPGcOupR/i6Gz2K1myJiIiIeJHCloiIiIgXKWyJiIiIeJHCloiIiIgXKWyJiIiIeJHCloiIiIgXKWyJiIiIeJHCloiIiIgXKWyJiIiIeJHCloiIiIgXKWyJiIiIeJHCloiIiIgXKWyJiIiIeJGxLMvXfXAzxuQCe7qhqXggrxvakfbTd+J/9J34J30v/kffiX/qju9lgGVZCW29ya/CVncxxqywLGuyr/shjfSd+B99J/5J34v/0Xfin/zpe9E0ooiIiIgXKWyJiIiIeFFvDVuP+7oD0oK+E/+j78Q/6XvxP/pO/JPffC+9cs2WiIiISHfprSNbIiIiIt2iV4UtY8xcY8xWY8wOY8wtvu6PgDEm1RjzmTFmkzFmozHmOl/3SWzGGIcxZrUx5l1f90XAGBNtjHndGLPFGLPZGHOUr/skYIy5wfV31wZjzMvGmBBf96k3MsY8bYzJMcZsaHIt1hizyBiz3fV7jK/612vCljHGAfwLOAUYBVxojBnl214JUAvcZFnWKGA6cLW+F79xHbDZ150Qt4eADy3LGgmMR9+Nzxlj+gHXApMtyxoDOIALfNurXutZYO5B124BPrEsaxjwieu5T/SasAVMBXZYlrXLsqxq4BVgvo/71OtZlpVpWdYq1+MS7B8g/XzbKzHG9AdOA570dV8EjDFRwEzgKQDLsqotyyr0ba/EJRAINcYEAmHAfh/3p1eyLGsJkH/Q5fnAc67HzwFndmunmuhNYasfsK/J83T0Q92vGGMGAkcC3/m2JwI8CPwGqPd1RwSAQUAu8IxravdJY0y4rzvV21mWlQHcC+wFMoEiy7I+9m2vpIkky7IyXY+zgCRfdaQ3hS3xY8aYCOAN4HrLsop93Z/ezBgzD8ixLGulr/siboHAROARy7KOBMrw4ZSI2FxrgOZjh+G+QLgxZoFveyWtsezSCz4rv9CbwlYGkNrkeX/XNfExY0wQdtB6ybKsN33dH+EY4AxjzG7s6fZZxpgXfdulXi8dSLcsq2HU93Xs8CW+NQf43rKsXMuyaoA3gaN93CdplG2MSQFw/Z7jq470prC1HBhmjBlkjAnGXsS40Md96vWMMQZ7Hcpmy7Lu93V/BCzLutWyrP6WZQ3E/v/Jp5Zl6V/rPmRZVhawzxgzwnVpNrDJh10S215gujEmzPV32Wy0ccGfLAQucT2+BHjbVx0J9FXD3c2yrFpjzK+Aj7B3jDxtWdZGH3dL7FGUnwDrjTFrXNdusyzrfR/2ScQfXQO85PrH4i7gUh/3p9ezLOs7Y8zrwCrsndWr8aOq5b2JMeZl4Hgg3hiTDtwB/BV41RhzObAHOM9n/VMFeRERERHv6U3TiCIiIiLdTmFLRERExIsUtkRERES8SGFLRERExIsUtkRERES8SGFLRERExIsUtkRERES8SGFLRERExIv+H5gtYGX6yCzKAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x10b4a2ba8>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.figure(figsize=(10, 6))\n",
"plt.plot(t, x, label='Data')\n",
"plt.plot(x2, y2, marker='.', label='Max 100, step 50')\n",
"plt.legend()"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.6.4"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment