{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## ThinkDSP\n", "\n", "This notebook contains code examples from Chapter 1: Sounds and Signals\n", "\n", "Copyright 2015 Allen Downey\n", "\n", "License: [Creative Commons Attribution 4.0 International](http://creativecommons.org/licenses/by/4.0/)\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Think DSP module\n", "\n", "`thinkdsp` is a module that accompanies _Think DSP_ and provides classes and functions for working with signals.\n", "\n", "[Documentation of the thinkdsp module is here](http://greenteapress.com/thinkdsp.html). " ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "# Get thinkdsp.py\n", "\n", "import os\n", "\n", "if not os.path.exists('thinkdsp.py'):\n", " !wget https://github.com/AllenDowney/ThinkDSP/raw/master/code/thinkdsp.py" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Signals\n", "\n", "Instantiate cosine and sine signals." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "from thinkdsp import CosSignal, SinSignal\n", "\n", "cos_sig = CosSignal(freq=440, amp=1.0, offset=0)\n", "sin_sig = SinSignal(freq=880, amp=0.5, offset=0)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Plot the sine and cosine signals. By default, `plot` plots three periods. " ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO29d5xkV3Xv+11V1bmqcw4zPaG7J2tG6hlJKAMCCYMEBtvC4Cvb8LCv4Tld+1q+fh+TbD+4fr5gbGMsY2wwYAkTBQgTFFBAM9KgydPTYXpC5xyqc6j9/qhTPa2ezpXOqVrfz6c+XXXqhHV619lr77V/e20xxqAoiqIodsMVbwMURVEUZTnUQSmKoii2RB2UoiiKYkvUQSmKoii2RB2UoiiKYks88TZgMxQWFprq6up4m6EoiqIsw89//vN+Y0xRuOdxpIOqrq7m+PHj8TZDURRFWQYRuRKJ82iIT1EURbEl6qAURVEUW6IOSlEURbEl6qAURVEUW6IOSlEURbEl6qAURVEUWxIRByUiXxCRXhE5u8L3IiKfEZEWETktIjcu+u5hEWm2Xg9Hwh5FURTF+USqB/VvwH2rfH8/UGO9PgD8I4CI5AMfBm4GjgAfFpG8CNmkKIqiOJiIOChjzHPA4Cq7PAh8yQQ5CuSKSBnwZuDHxphBY8wQ8GNWd3Rh81xTH//ja6eYmJmL5mWUCGOM4fFXrvKFFy4xNx+ItznKBpgPGB597iL/dbYLXX9O2QixyiRRAbQt+txubVtp+3WIyAcI9r7YsmXLpg1p6vHzzRPtnG4f5rPvuZGaEt+mz6XEhtn5AB9+4hxfPXYVgO+c6uRvfukAO4u17OzOzFyAP/jaSb5/uguA1+0o4MNv20tdqZadsjaxEknIMtvMKtuv32jMo8aYemNMfVHR5lM8vf+O7Xz5fTczNDHDA3//It860b7pcynRZ2Rill//15f56rGr/M7dO/jMuw9xdWCct3zmBf7ppxeZD2iL3K5MzMzx/i8d5/unu3jk/l18/O37ON81yls+8zwfeeIcIxOz8TZRsTmxclDtQNWiz5VA5yrbo8ptOwv5/u/ewf7KHP7g8VM88o3TTM3OR/uyyga51D/OOz77Ii9fGuRvfukG/ud9u3jghnJ+9Ad3cU9dEf/vDy7wS5/7Gf4prejsxsjELO/9/DFeaO7jk+/cz2/ftYNfu2Urz/yPu/nVI1v40kuXectnntdQu7IqsXJQTwD/zVLz3QKMGGO6gB8CbxKRPEsc8SZrW9QpyU7nq++/mQ/es4PHXmnjcz+9GIvLKuukZ3SKd3z2RYYmZvjK+2/hnTdVLnxX5Evjc++9ib/5pRt49eow/340InkplQgxOD7Drzz6Emc7Rvnse27kVw5fC8nnZaXy8bfv44u/eYSO4cmFsK1iDwbHZxibtk+jIVIy8/8AXgLqRKRdRN4nIr8tIr9t7fIk0Aq0AP8M/A6AMWYQ+DjwivX6mLUtJnjcLv74zbt4w65ivvTSFe1F2YgvvXSZ0clZvvZbt3JkW/5134sI77ypkjtri/jCC5e17GzEP/30Is29Y3zh1w9z376yZfe5o6aI1+0o4J+ea9WysxF/93Qzd3zyaduUSaRUfO82xpQZY1KMMZXGmH8xxnzOGPM563tjjPmgMWaHMWa/Meb4omO/YIzZab3+NRL2bJT337GdwfEZvvlqRzwuryxhanaerx67yht3l6wpYvntu7bTPzbNN17VsUQ7MDU7z+PH23jz3hJurylcdd8PvX4nff5p/vN426r7KbFhcmaeb/y8ndtrikhPccfbHEAzSQBwy/Z89lVk8/kXWgnooHvc+faJDoYmZvmN27atue+t2wu4oSqXR59rVcGEDXjiVCfDE7P82i3Va+576/YC6rfm8Y/PXmRmTqcOxJvvne5kdGqO99y8eZV0pFEHRTBc9H/dsZ3WvnGeaeyNtzlJjTGGf33xMrtKfdyy/frQ3lJEhP9+13auDEzwg7NdMbBQWQljDF966TK1Jd51l92HXr+TzpEpVdPagK8cu8qOoixuXiakHi/UQVm8ZX8Z5Tnp/PPzrfE2Jal56eIAjT1+fvP2bYgsNwvhet60p5TtRVn847MXdSJoHDnRNszZjlF+7dbqdZfdXbVF7K/I4bPPXtQJ2HHkbMcIJ9uGec/NW9dddrFAHZRFitvFb9y2jaOtg5xpH4m3OUnLF168TEFWKg/cUL7uY1wu4bfu3M65zlFeaOmPonXKavz7S1fwpXn4xUPLzrVfllAv6srABN87rT3gePHVl6+SnuLinTdWrr1zDFEHtYhfOVKFN82jvag4cWVgnKcu9PCrN2/Z8CDt2w9VUJKdptMF4kT/2DTfP93FO2+qJCttYwlq7t1dQl2Jj79/pkXHgOPA2PQc3znRwdsOlJOTmRJvc16DOqhFZKen8NDhKr5/pouO4cl4m5N0/NvPLuMW4b23bN3wsWkeN++7fRsvtgxwun04CtYpq/H4K23MzAc2VXYul/DB1++kpXeMH53viYJ1ymp8+0QH4zPzvGcTZRdt1EEt4TduDyrH/u3FS3G2JLnwT83yn8fb+YUDZZRkp2/qHO8+sgVfuod/+9nlyBqnrMrcfIAvH73C7TsL2Vns3dQ5fmF/GflZqfyXCl1iijGGrxy7yt7ybG6ozIm3OdehDmoJFbkZvHlvCd98tUPDDTHkGz9vZ2x6bl3S8pXwpafwxt0lPNvYp2UXQ37S0EvXyBS/duvmW+Bul3B3bRE/berT6QIx5ETbMA1do7YTR4RQB7UM9+4pYWB8hnOdo/E2JWl48mw3u0p9HKzKDes8d9cVMTg+w+kOFbrEiq8cu0JFbgZv2FUc1nnu2VXM0MQsJ9uGImSZshZfOXoVb5qHBw6uX5QUS9RBLcMdNcFs6c/qnKiY4J+a5dUrQ9xdF14FB3BnTREi8MwFLbtYMD49x9HWAd56oAyPO7zq5M7aItwu4Wktu5gwOTPP90538uDBcrwbFLbECnVQy1DoTeNAZQ4/beqLtylJwdHWQeYChjtrV0+Nsx7yslI5WJXLs1p2MeHYpQFm581Coy4ccjJSuGlrHk9f0LKLBS9fHmR6LsCb95bG25QVUQe1AnfVFvHq1SFdsyYGPNfUR2aqm/qtkZnBfk9dMafbh+kfm47I+ZSVeb65nzSPi/rqvIic7/W7imnoGqVrRFW00ebFln5S3S4OV9snc8RS1EGtwF21RQQMvHhRJ35Gm+ea+7h1ewGpnsj8HO+pK8aYoONTosvzzf0c2ZYfseSir7fGsZ7RXlTUeaG5n/rqPDJS7ZEYdjnUQa3AwapcfOkeftqoD0o0uTIwzpWBCe5YI/P1Rthbnk2hN5VnteyiStfIJC29YxEtu5piLxW5GToOFWUGxqY53zXKbTsjV3bRQB3UCnjcLu6oKeSnTX2a3y2KPNcc7KHeWRv+GEYIl0u4q7aY55pVshxNXrDK7vadkSs7EeENu4t5saXfNmsSJSI/uzgAwO3qoJzLXbVFdI9O0dQzFm9TEpbnmvqozMtgW2FWRM97z64ihidmOdmmWSWixfPN/RR609hVuvqaXRvlnl3FTM7Oc7R1IKLnVa7xQnM/2eke9lXYb3LuYtRBrUKoVf/TJg03RIPZ+QAvXRzgztqiiE8SvGNnES7RqQLRIhAwvNjSz+07C3C5Ilt2t24vID3FpVMFooQxhhda+nndjkLcES67SKMOahXKcjKoK/HpWEaUePXKEGPTc9wZAYnyUnIyg5JlLbvocL5rlIHxmYjIy5eSnuLmth2FPN3Yq+H1KHBlYIKO4Ului+DYYbSIiIMSkftEpFFEWkTkkWW+/5SInLReTSIyvOi7+UXfPREJeyLJXXVFvHJ5kPHpuXibknA839yP2yW8bmdBVM5/d10xZzpG6PVPReX8yUxoWZO1lnXfLPfsKqZtcJKLfRpejzQLZWfz8SeIgIMSETfwD8D9wB7g3SKyZ/E+xpg/MMYcNMYcBP4O+OairydD3xljHgjXnkhzd20Rs/OGly5qPDzSPNfcx6GqXLLTo5Pi/+46K0SrvaiI83xzH3Ulvk0n9l2Leyy5uar5Is+LLf1U5GZQXZAZb1PWJBI9qCNAizGm1RgzAzwGPLjK/u8G/iMC140JN1XnkZnq1qwSEWZwfIYzHSMRVe8tZU9ZNsW+NM0qEWGmZud55fJQ1HpPEEzavKvUpw4qwswHDD+7OMDtOwttmRx2KZFwUBVA26LP7da26xCRrcA24OlFm9NF5LiIHBWRt690ERH5gLXf8b6+2FU4aR43r9tRwLNNGg+PJM8392FMZOXlSxER7q4r4rmmPl1OPIK8fGmQmblAVB0UBENQr14dZmZOyy5SnO0YYWRy1hHjTxAZB7WcG16pJn8I+LoxZvEEhy3GmHrgV4FPi8iO5Q40xjxqjKk3xtQXFUWvUluOu2qLaBuc5FL/eEyvm8g839xPbmYK+6Msc71tZyH+qTkae/xRvU4y8XxzH6luFzdvi26KnBu35jEzF6ChS1cViBSh8afX7YjOuG+kiYSDageqFn2uBDpX2PchloT3jDGd1t9W4FngUARsiiivswYTX7k8GGdLEgNjDM8393HbzujLXG/cEswRd+KqzoeKFM8393PT1jwyU6ObAfvQluDSKyeu6vIbkeLFln52l2VT6E2LtynrIhIO6hWgRkS2iUgqQSd0nRpPROqAPOClRdvyRCTNel8I3Aacj4BNEWVbQRbZ6R6d9Bkhrg5O0DM6HZNWXGVeBgVZqVp2EaLXP8WFbj93RCDz/FqU5WRQmp3Oq9q4iAiTM/McvzzE7VFSzUaDsJtAxpg5EfkQ8EPADXzBGHNORD4GHDfGhJzVu4HHzGsHcnYD/yQiAYLO8hPGGNs5KJdLuKEqV1vhESLkLA5VRSYD9mqICIe25GorPEIcvxz8P966PTaV3I1bczmhCxhGhONXBpmZD9g+/95iItJHN8Y8CTy5ZNufL/n8kWWO+xmwPxI2RJtDVbn8/TMtjE/PkWXTxb2cwsm2YdJTXNSWeGNyvUNb8vhJQy8jE7PkZEZH0p4snGofJsUt7CnPjsn1DlXl8eSZbvr80xT5nBGWsisvXxrE7RKORHnsMJJoJol1cnBLLgEDZ3Qp8bA51TbM/oqcsFdgXS+hZeRPtWsPOFxOtQ2zpyybNE9slmi4cauOQ0WKU+0j1Jb4oj52GEnUQa2TGyqDD4qOZYTH7HyAs52jC04jFhyozEFEhRLhMh8wnGkf4YYYlt3e8hxS3MIJfe7CwhjD6fZhbqi0d3LYpaiDWicF3jS25GdySh+UsLjQ5WdmLhDTSs6XnkJtsU/HMsLkYt8Y4zPzC421WJCe4mZPWTavXtGyC4ergxMMT8xyIIZlFwnUQW2Ag1W52oMKk5NWmC2WlRxcKzudbL15Qo2zG6pi2wo/tCWP0+0jOtk6DE61B4cmDmgPKnE5WJVL18gUPaOafHSznGobptCbSmVeRkyve2hLLsMTs1wemIjpdROJU+3DeNM8bC+MjbglxKEtuUzOzutk6zA43TZMmsdFXYTX7oo26qA2wMGFiYPai9osJ9uGuaEyN+Z5wA7qpM+wOdU2woHKnIiv/7QWocnWOh9q85xuH2FveTYpMRImRQpnWRtn9pRlk+IWDfNtktGpWS72jcV0/ClETbGPrFS3lt0mmZqdp6FrNC5lV5mXQaE3TRsXm2RuPsCZjhHHjT+BOqgNERqwPamD7ZvibPsIxhBTBV8It062DouGrlHmAiYuKrDQZOuTWnaboqVvjMnZ+ZiPHUYCdVAb5GBVLmfaR5gP6GD7RglJheM1UHuwKpeGrlEmZ+bX3ll5DdcEEvFphd+4JY/W/nGGxmficn0nc7otJJDQHlTCc3BLLuMz8zT36oDtRjnVNsy2wixyM1Pjcv1DW/KYCxjOdupk641yqn2EYl8apVFaoHAtQoljNUS7cU61D+NL87CtICvepmwYdVAb5KCVP07DDRvnVPtwXMJ7IULX1rLbOKfahrmhKvbilhAHKnNwu4RXdRxqw5xuH2F/HMQtkUAd1AapLsgkNzNFW3IbpGtkkp7R6bjOZC/ypVGZl6ETdjfIyOQsrf3jcW1cZKZ62FXq0zHEDRIStzgxvAfqoDaMiHBDpU7Y3SjxHsMIcWhLnlZyG+SMTSZ5HtoSfO50/Hf9xFPcEgnUQW2Cg1W5NPX4GZ+ei7cpjuFEW2yzYK/EIWuydfeITrZeL6Ekuwcq4tu4uHFLHmPTc7T0jsXVDidx2mpcxLthuFnUQW2CUGbzUOEraxPrLNgrcXBhsF3DfOvlZNsw2wuz4r5USaiS1RUF1s+p9mEKvWmU5cRH3BIu6qA2wUHNbL4h4pEFeyWCs+mFk21aya2X0+3Dtii76oIs0lNcnO8cjbcpjuF0+wg3VObETdwSLuqgNkFeVirVBZrZfL209AazYMdzkD1EmsfNjiIvDV1aya2H7pEpekan4z7+BMHJ1rtKs7Xs1onfytziVIEEqIPaNHvLczivD8q6sItAIsSeMq3k1stJm5Xd7rJszneNalb6dXCmI5i55YADM0iEiIiDEpH7RKRRRFpE5JFlvv91EekTkZPW6/2LvntYRJqt18ORsCcW7C7zcXVwAv/UbLxNsT3nOkfISnXbZqLg7rJsev3TDIxNx9sU23O6fRiPS9hTFl9xS4g95dmMTM7SpSKXNVkQSCRzD0pE3MA/APcDe4B3i8ieZXZ93Bhz0Hp93jo2H/gwcDNwBPiwiOSFa1Ms2G09sI3dmlFiLRq6/dSV+mwzUTBUdg1dWnZrcaZjhLpSH+kp8RW3hNhTFlwuQseh1uZ0+zCVeRnkZ8Unc0skiEQP6gjQYoxpNcbMAI8BD67z2DcDPzbGDBpjhoAfA/dFwKaoc62S0wdlNYwxNHSNLvy/7MBuq5LTslubhi6/rcqurlSfu/USFEg4t/cEkXFQFUDbos/t1ralvFNETovI10WkaoPHIiIfEJHjInK8r68vAmaHR1lOOjkZKTRoD2pVOkem8E/N2aqSK/CmUZKdppXcGvT5p+kfm2aXjRa586Z5qC7I1PHfNRidmqV9aDLu8w7DJRIOarm4zdIRzO8C1caYA8BPgC9u4NjgRmMeNcbUG2Pqi4qKNm1spBARdpf5tJJbgwYrFBPqtdiF0GC7sjKh8LVdxp9C7FaRy5o0WWVnp8bFZoiEg2oHqhZ9rgQ6F+9gjBkwxoRGpP8ZuGm9x9qZXaXZNHb7CWjqlRW50B2sSEKhGbuwuyyblt4xpud06Y2VuFZ29qrk9pRlc3lggjHN5LIiocjOLps1LjZKJBzUK0CNiGwTkVTgIeCJxTuISNmijw8ADdb7HwJvEpE8SxzxJmubI9hTls3EzDxXBifibYptaejysyU/E2+aJ96mvIbdZdnMBYymzVmFhi4/xb40Crxp8TblNVwTKGkvaiUudI3iS/dQ7tAMEiHCdlDGmDngQwQdSwPwNWPMORH5mIg8YO32uyJyTkROAb8L/Lp17CDwcYJO7hXgY9Y2R6BCibVp6B61XXgPrqnBVMm3Mg1do7ZsgYfGVc5r2a3IhW4/u0uzHZtBIkREmrXGmCeBJ5ds+/NF7/8U+NMVjv0C8IVI2BFrakq8uF1CQ9cob9lftvYBScbkzDyX+8d524HyeJtyHdUFWaR5XNq4WIHZ+QAtvWPcUVMYb1OuIyRQUqn58gQChsZuP79447J6M0ehmSTCID3FzfbCLG2Fr0BTj5+AwVYKvhAet4tdpSpyWYlL/ePMzAdsWXYiotlAVqFjeJKx6Tl22WzcdzOogwoTVRStTOj/YscQH1wrO02bcz2hsttl47K70D2qa0Mtg93LbiOogwqTXWU+OoYnGZnUlEdLudDtJyvVTVVeZrxNWZbdZdkMTczSM6opj5ZyodtPilvYXuiNtynLsrvMx9RsgMsD4/E2xXZcsBR8dSXqoJKeUAjkgvairuN816itUhwtRUUuK9PQNcqOIi+pHntWEQtCCR2Huo4L3aNsLcgky2bK2c1gz1+fg9ijldyyGGO4YLMUR0sJhUB0wu71XOjy226C7mJ2FnvxWAIl5bVc6PI7foJuCHVQYVLsSyM/K1WFEkvoHJlidGrOljLlENnpKVTmZaiDWsLQ+Azdo1O2HsNI87jZWezVslvC5Mw8lwbGE0IgAeqgwmYh5ZFOGnwNoZDnHhtXcqAil+W4sJAmx96VnCr5rqepx48x9hUmbRR1UBFgt5XySBVF1whVHHZLcbSUPWXZXO4fZ3JGUx6FCKU4snMPCoLjUD2juq7XYuyaWmyzqIOKALvKspmeC3CpXxVFIRq67ZniaCm7y7IJGGjs0RBtiIauUQqyUimyWYqjpei6XtdzodtPRoqbLfn2VM5uFHVQEUDXF7qehq5RRwzUqsjlei50B9eAsnuanJCDOt81EmdL7MOFLj+1pT7cNlXObhR1UBFAFUWvJZTiyM4KvhCVeRl40zxadhbzVpocJzQu8rNSKc1O1x6UhTGGC92j7HZA2a0XdVARIKQo0kouyLUUR/Z/UFwu0ZRHi7g8MM70XMDW6svF7CrzLYg6kp1e/zRDE7OOaFysF3VQESKoBtMHBa4N1DqhBwVW2pwuv6Y8IhgiAucsdFdb4uNi3xhz84F4mxJ3rqU4csZztx7UQUWI3WU+ukenGBqfibcpcaehy94pjpZSW+rDPz1H18hUvE2JOw1do7hdQk2JPVMcLaW2xMfMXEDXZGPx9ABnNC7WgzqoCBGSdaoaLFjJ2TnF0VJqi4OVcZOWHRe6R9lRlEWaxx1vU9ZFreVIm7XsuNA1SllOOrmZqfE2JWKog4oQ+qAECQ7U+h01D6PWSqqpDirY+7X7BN3F7Cz2IgJNPboy8gWHiFs2gjqoCFGanY4vzZP0D0qff5qRyVnqHBIiAsjLSqXIl5b0ZTc6NUvH8CR1DqrkMlM9VOVlJn3kYmYuuMBkIo0/QYQclIjcJyKNItIiIo8s8/0fish5ETktIk+JyNZF382LyEnr9UQk7IkHIsG4fbK3wkOVfK3DUv3XlniTvvcbun+ntcK17AgKRQLGcWW3FmE7KBFxA/8A3A/sAd4tInuW7HYCqDfGHAC+DvzvRd9NGmMOWq8HwrUnntSW+GjuTe5WeMhB1zjMQdUUB8sukMTpqpzbuPDR2jfOzFzyKvkaHZI/caNEogd1BGgxxrQaY2aAx4AHF+9gjHnGGBOS2RwFKiNwXdtRU+JjcHyG/iTODdbc6ycvM4VCr7MGamtLfEzMzNMxPBlvU+JGU08wTU5Fbka8TdkQtSU+5gImqRcvbOrx43EJ2wqz4m1KRImEg6oA2hZ9bre2rcT7gB8s+pwuIsdF5KiIvD0C9sSNkFAimcN8TT1j1JT4bJ8mZyladtDcM0ZNidcx6ssQNVp2NPeOUV2YZdsFJjdLJO5muV/zsnESEXkvUA/89aLNW4wx9cCvAp8WkR0rHPsBy5Ed7+vrC9fmqBAKjTQn6WC7MYamHv9CZe8kahaUfMlZdhCs4GuKnRXeA9hR5MUl0JTEGSWaHfrcrUUkHFQ7ULXocyXQuXQnEXkj8GfAA8aYhRiYMabT+tsKPAscWu4ixphHjTH1xpj6oqKiCJgdeYp9aWSne5K2JdczOo1/as5xYxgAORkplGanJ+1g+8jELL3+aUdWcukpbqoLspK2cTE1O8+VwQl2OrBxsRaRcFCvADUisk1EUoGHgNeo8UTkEPBPBJ1T76LteSKSZr0vBG4DzkfAprggItSW+JLWQS0IJBz6oNSUeGnqTc6ya+4NiVuc56AgGL1I1rJr6R3DGBzZuFiLsB2UMWYO+BDwQ6AB+Jox5pyIfExEQqq8vwa8wH8ukZPvBo6LyCngGeATxhjHOigIhoqaesaSMq9byEE59UGpLfHR3DOWlAtPhnofTm1c1JZ4udw/ztRs8i082dLrTPXleojIanLGmCeBJ5ds+/NF79+4wnE/A/ZHwga7UFfi5T9enqXPP01xdnq8zYkpzT1jFGSlUmDzhe5Woq7Ex/RcgLbBCaoTTA21Fk09fjJTnafgC1FT4iNgoLVvnD3liSW1XouQgq+6IPF+s4kl+bABtUk82N7U63dsiAiSWw3W3Ounpth5Cr4QoewXzUkY5mvqGWNbAir4QB1UxKlJ0rxuxhhaesYcHWYIlV0yTrYOTQ9wKtUFWXhcsjBhNZlodnjDcDXUQUWYQm8qeZkpSdeS6xqZwj895+hKzpvmoSI3I+kqueGJGfocquALkepxsb0o+ZR8U7PzXB2ccOzY4Vqog4owwZx8vqR7UBYEEsXOreSApMynuCCQcHDjAoL2J1vD8JqCz9lltxLqoKJArVXJJZOSr9mhedyWUmfldUumFVqvTQ9wduOittjH1cEJJmeSR8kXcshO7v2uhjqoKFBb4sM/NUfPaPLk5Gvq8VPoTSMvy1k5+JZSU+JjZj65Vmht7gmugOxUBV+IulIvxlyTXScDTT1jQQVfgqpO1UFFgVA8OJnWqGnqHUuIVlwyLjzZ3DvGTgfmT1xKMgqUmi0FX4o7MavyxLyrOJNslVxQwed3fIgIgiu0AjR2J1cr3OljhwBb8zNJdbuSy0H1+h0fVl8NdVBRoMCbRkFWatI8KB3Dk4zPzDt+kB2CK7Ruyc9MmrQ5Q9byMIlQyXncLnYUJ4/IZXLGUvAlQORiJdRBRYmgGiw5WuGJIpAIkUwrtF5bYDIxKrnaJHruLvYltoIP1EFFjboSnyUBTXwln9Nz8C2lJolWaG1KsDxutSU+OoYnGZuei7cpUSfRFXygDipq1JT4GJueo3NkKt6mRJ2mnjGKfGnkZjpbwReitsSbNCu0Nvf48aZ5KMtJjLyR19ZkS/wecFPPGCluYWsC5uALoQ4qStQmkaIoOFCbOK24kAozGcquqcfPzmKv4xV8IUJCnWRIV9Xc409oBR+og4oayaLkCwRMcKnwBEq1srM4uEJrMqyM3NyTGNMDQlTlZ5LmcSXFXCin509cD+qgokRuZipFvrSEH7DtGJ5kcnY+YcYwILhC65b8zISv5AbGphkYn0mosnO7hB1Fia/km5yZp21ogtoEagpUDeUAACAASURBVBguhzqoKFJT7E34UEOiCSRC7CxO/JWRQ7/NRGuF15R4E773e03Bl1jP3VLUQUWR2hIfLQmeky+RK7lL/ePMJnBOvuYEbVzUFHuDc/MSWMl3bXpAYj13S1EHFUV2FnsZn5lPaCVfc88YJdlp5GSkxNuUiBJS8l1JYCVfU88YvjQPpQm28nOo0k7kEO01BV9mvE2JKhFxUCJyn4g0ikiLiDyyzPdpIvK49f0xEale9N2fWtsbReTNkbDHLiSD5DW4EmviteKuKfkSt5Jr7vWzsyRxFHwhkkHJ19LrZ3uhN6EVfBABByUibuAfgPuBPcC7RWTPkt3eBwwZY3YCnwI+aR27B3gI2AvcB3zWOl9CsPCgJGglFwgYWnrHFvLXJRI7irxIgiv5mnvGEnKQfYuVky+R14YKKvgS77lbSiTc7xGgxRjTaoyZAR4DHlyyz4PAF633XwfeIMFm24PAY8aYaWPMJaDFOl9CkJeVSqE3NWEflI7hSSZmEkvBFyIj1U1VXmbCll1IwZeIlZzHHVxdN1EbFyEFXyJGLpYSCQdVAbQt+txubVt2H2PMHDACFKzzWABE5AMiclxEjvf19UXA7NhQU+xL2FBDy4JAIvEqObBUmAlaySWquCVEIq+umywKPoiMg1ougL1UtrbSPus5NrjRmEeNMfXGmPqioqINmhg/akq8tPQkZk6+RFmJdSVqSny09o8l5Oq6Cw4qUcuu2Evb4CQTM4mn5Eu0BL+rEQkH1Q5ULfpcCXSutI+IeIAcYHCdxzqammIv/uk5ukcTT8nX3JtYOfiWUlPsZXbeJOTquomWg28pod7Fxd7EU2E29yZ+Dr4QkXBQrwA1IrJNRFIJih6eWLLPE8DD1vt3AU+bYJfiCeAhS+W3DagBXo6ATbbh2iqfiRcqau5JrBx8S6lJ4HRVzT1jCZWDbyk7rfGZRAzzNfckh4IPIuCgrDGlDwE/BBqArxljzonIx0TkAWu3fwEKRKQF+EPgEevYc8DXgPPAfwEfNMbMh2uTnbim5EusB8UYQ3NvYuXgW8rOBFZhNvcmVg6+pWwtyCTFLYnZMOwdY2cCl91iPJE4iTHmSeDJJdv+fNH7KeCXVjj2L4G/jIQddqTAm0Z+VmrCTRoMKfgSOQ6emeqhMi8j4UQuoVV0E7lxkeJ2sb3QS0uC9aBCq+j+4qHKeJsSExK/j2gDahJwGeprg+yJW8lBgpddAjcuAHaWJF4uzJCCL9HLLoQ6qBhQYz0oiaTka+lJbBVYiNoSH6394wml5EuWPG41xV6uDk4wOZM4owbJsIruYtRBxYDaEh/+qTl6/dPxNiViNPX4KfSmkZeVmAq+EDuLvczMBWgbmoy3KRGjpXeMrFQ35Qmq4AtRU+zDmGCvI1FIhlV0F6MOKgaEBtsTKVQUFEgkfisuEVdGburxs7PEl7AKvhChXkYijf8294wl/Cq6i0mOu4wzoXGaRFGDGRPMwZcMYYYdxQlYyfWOUZsEjYutBVl4XJJQjYtETc68EuqgYkChN5W8zJSEmZPRNTLF2PQcOxN8DAPAm+ahIjcjYaYJDE/M0OefTopB9lSPi+rCrIQRSkzNBhV8yVB2IdRBxQARCebkS5Ae1MIquknQCoegyCVR5tMkeg6+pdSWeBOm99vSayn4tAelRJqdCaTka0mySq6m2MvFvjHmA84vu0TPn7iUncU+rgyMMzXrfCVf6LlLhtB6CHVQMaK22MvI5Cx9CaDkCyr4UslPcAVfiJpiH9NzAdqHnJ+Tr7knqOCryM2ItykxoabYS8BAa5/zc/I19fjxuJJHwQfqoGJGqLeRCPHw5gRdpHAlQjH/RAjzhRaYTHQFX4iFVa0TYPy3uTeo4Ev1JE+1nTx3GmcSJSefMYaWnsTOwbeUhZx8CVDJNfX4kyY0C1BdmInbJQkx/tvc408qgQSog4oZRb40cjJSaHJ4D6p7dAr/9FxSxcF96SmU5aQvZM9wKiMTs/T6p5Nm/AkgzeOmuiDT8VLzqdl5rgwmxyq6i1EHFSOCSj6v4yu5UEt0Z5I9KDUlPhodXsldS5OTXGVXV+pzvINKthx8IdRBxZCaEh9NvX5HK/kWJOZJ9qDUWXJlJyv5QuOfyTR+CEGHfMXhOflCDcNka1yog4ohtSVehidm6RtzrpKvpXeM/KxUCrxp8TYlptSWBJV8VwacqwZr6vGTmUQKvhB1JcGcfE6eD9XcG1TwVSeRgg/UQcWUOqv109jt3HBDU48/6VrgEAwTgbNz8oUUfC5Xcij4QtRaZefkEG1TzxjVSabgA3VQMWXhQXGogzLG0NQzxq7S5AozQHAulAg0dju3Fd7Uk1x53EJszc8k1eNyfOMimcQtIdRBxZBCbxqF3lTHPigdw5OMTc8t9CaSiYxUN1vznasGG5mYpWd0OunGDgE8bhc7i7yObRhOzc5zZWA8qaYHhAjLQYlIvoj8WESarb95y+xzUEReEpFzInJaRH5l0Xf/JiKXROSk9ToYjj1OoLbE59gHJWR3XRI+KBAsuwvdo/E2Y1OE7E7GxgU4W8l3sW+MgEme9FSLCbcH9QjwlDGmBnjK+ryUCeC/GWP2AvcBnxaR3EXf/7Ex5qD1OhmmPbYn+KCMEXCgGuyC5aBqk7iSuzww4ci8bqHxl2R1ULUlPrpGphiZnI23KRsm1DBMxtB6uA7qQeCL1vsvAm9fuoMxpskY02y97wR6gaIwr+tY6kp8TM7O0+bAvG5NPX7Kc9LJTk+Jtylxoa7Ux3zAODKvW2O3n+x0D6XZib2K7krUlTo3k0tjt59Ud3DpkGQjXAdVYozpArD+Fq+2s4gcAVKBi4s2/6UV+vuUiKyoXRaRD4jIcRE53tfXF6bZ8aPOwUKJxm5/0rbA4Vpo04mhosZuP7tKs5MmB99SQvOHnKjku9DtZ0exN2lW0V3MmncsIj8RkbPLvB7cyIVEpAz4d+A3jDEBa/OfAruAw0A+8CcrHW+MedQYU2+MqS8qcm4HrMahUvPZ+QAX+8aoK82Otylxo7owixS3LIQ6nYIxhsae5G5cVORm4E3zOO65g1DjIjnLzrPWDsaYN670nYj0iEiZMabLckC9K+yXDXwf+H+MMUcXnbvLejstIv8K/NGGrHcg3jQPlXkZjmvJXeofZ3beLIRKkpEUt4sdRV7H9aA6R6bwTyWn+jKEiFBb4jwl3/DEDN2jU0lbduH2GZ8AHrbePwx8Z+kOIpIKfAv4kjHmP5d8V2b9FYLjV2fDtMcR7HKgoujCgoIveXtQ4EwVZqOl4EvWVniIkJLPSanGFp67JC27cB3UJ4B7RaQZuNf6jIjUi8jnrX1+GbgT+PVl5ORfEZEzwBmgEPiLMO1xBLUlPlr7xpmZC6y9s01o6vbjdgk7ipNvoHYxdaU+OoYn8U85Rw2W7OrLELUlPoYclmosmRV8sI4Q32oYYwaANyyz/Tjwfuv9l4Evr3D868O5vlOpK/UxFzC09o+xyyFjOhe6/WwrzCLN4463KXHlmlBijJu2Xjftz5Y0die3+jLEQtl1j1Hsc4aa8UKSqy+TTxZiA5yo5GvsGU3aMMNinJiTL9nVlyGcmJOvsXs0qdWX6qDiwPZCLx6XOMZBjU/P0TY4ya4kzSCxmIrcDDJT3Y4pO1VfXqPQm0ZBVipNDim7UO7LZG5cqIOKA6keF9uLshxTyS2sAZXED0oIl0uoLXGOyCWkvkzWMYyl1Dpo4cn2oWDuy11lyVt26qDihJMelGQfqF1KnYOUfMmuAltKXamP5h6/I1KN6XOnDipu7Cr1LbSQ7E5jj5+MFDdVeZnxNsUW1Jb6GBifod8BarDG7lE8LmFHUfLOX1tMbYmP8Zl5OoYn423KmjQurF6tDkqJMbUOSpvT2O2ntiT5FrpbiWtqMGeU3bYkXOhuJUITzZ3w3F3o9lORm4EvidWX+quNEyF5uVMqOQ0RXaPWquScEKK9oGX3GmoclJMvqOBL7rJTBxUnKvMyyEhx2z6vW//YNAPjM0kdZlhKkTeN/KxU249DjU3P0T40mfSV3GKy01Moz0m3fcNwZi5Aa9940jcu1EHFiaAazP553a4N1KpMOcRCXjebl13TwhpQWnaLqS310dgzFm8zVuVi3xhzAaMOKt4GJDNOWOWzUVVgy1JX4qOp29553VQFtjx1JT4u9o4xO2/fVGMXFvInJnfjQh1UHKkt8dE/Zm81WGO3n/ysVAq9qfE2xVbUlgbVYO1D9lWDNXb7yUp1U5GbEW9TbMXusmxmrAnMduVCt58Ut7C9KLlzX6qDiiNOEEpc6PFTV+JL2lQrKxEqu4au0ThbsjIXukepLfWp+nIJe8uDZXeuw75l19jtZ0dRci5SuJjkvvs4E1KDNdjUQQUChuYkX+huJXaX+RCBc532rOSMMUH1pYpbrmN7kZf0FJdtyw5UORtCHVQcKfalU+RL41znSLxNWZb2oUkmZub1QVmGzFQPO4q8ti27Pv80QxOzWnbL4HYJu0qzbVt2IxOzdI0k7yKFi1EHFWf2lWfbNtQQeoB3lyX3QO1K7C3Ptm0rvLFHxS2rsbc8m/Ndo7YUuYTKbneSCyRAHVTc2VeRQ3Ovn8mZ+Xibch1nO0es1qZWcsuxrzyHrpEpBmwocgmNjSW7Cmwl9lXk4J8KZum3G6EVkLVxoQ4q7uwtzyFgrslK7cSZjlFqir2kpyT3IoUrsTDYbsNe1JmOUSpyM8jPUvXlcoTK7qwNw3wN3X586R7KcpJzkcLFqIOKM/sqQg+KvSo5YwznOkbYX5ETb1Nsy97y4P/GjpXc2Y6Rhd+Wcj21JT7cLrHlONS5jhH2leeocpYwHZSI5IvIj0Wk2fq77BrYIjIvIiet1xOLtm8TkWPW8Y+LSNI19ypyM8jNTOFch70elO7RKQbGZ9inDmpFcjJTqMzLsF0PanRqlkv949q4WIX0FDc1xV7bld3MXICGLj/7K7XsIPwe1CPAU8aYGuAp6/NyTBpjDlqvBxZt/yTwKev4IeB9YdrjOESEfeU5tmuFn2kP2qMOanX2lefYrnEREt1o2a3OHhuKXJp6/MzMB7TsLMJ1UA8CX7TefxF4+3oPlGD/9fXA1zdzfCKxtyKbxm4/M3P2Sb1ytnMUlwTn+ygrs7c8m8sDE/inZuNtygJnLYepPajV2VueQ59/ml7/VLxNWeCMVXYHtOyA8B1UiTGmC8D6W7zCfukiclxEjopIyAkVAMPGmNCKfe1AxUoXEpEPWOc43tfXF6bZ9mJfeQ6z88ZWefnOdoywo8hLZqon3qbYmlBL97yNWuJnOkYoz0mnwJsWb1NsjR1FLmc6RvCle9haoIuDwjoclIj8RETOLvN6cAPX2WKMqQd+Ffi0iOwAlhsBXHFSgjHmUWNMvTGmvqioaAOXtj+hSs5OA7ZnVSCxLuxYyQUFElp2a7HHKjs7NS7OqkDiNazpoIwxbzTG7Fvm9R2gR0TKAKy/vSuco9P62wo8CxwC+oFcEQk10SuBzrDvyIFszc/Em+bhrE0m7PaOTtHrn2avVnJrUpwdzAZilzFE/9QsrSqQWBfZ6SlsLci0TcNwZi7AhS4/B1QgsUC4Ib4ngIet9w8D31m6g4jkiUia9b4QuA04b4JTuJ8B3rXa8cmAyyXWgK09HpRQZauV3PrYW55tm1Z4qCenPaj1YadsICqQuJ5wHdQngHtFpBm41/qMiNSLyOetfXYDx0XkFEGH9AljzHnruz8B/lBEWgiOSf1LmPY4ln3lOZzvGmU+EP/UK2c7RhG5FgJRVmdfeQ7NvWNMzcY/G0hIIKGV3PrYW57DlYEJRm0gclkQSGgPaoGwRsCNMQPAG5bZfhx4v/X+Z8D+FY5vBY6EY0OisK8im6nZAK19Y9TEOQP1mY4RthVm4U1TgcR62FuezXwgmD38hqrcuNpypmOEUivsqKxNqBHW0DnKzdsL4mrL6fYRstM9bMlXgUQIzSRhE0ItXjuMZYRmsivrw05ld0YFEhvCTiKXkLhFBRLXUAdlE7YXZpGe4oq7UGJgbJrOkSlNk7MBKvMyyE73xL2SG5ue0wwSG+TakjfxLbvpuXkudI9qBoklqIOyCR63i91l2QtjCPHirA6ybxgRS+QS57I73zmKMbC/UhsXG2GvDQRKTd1jzM4bbVwsQR2UjdhXnsP5zlECcRRKhBzkXg3xbYh95Tk0dPuZnY9fNpAzKpDYFHvLs+Mucjmj2T+WRR2UjdhXkY1/eo6rgxNxs+Fc5whbCzLJyUiJmw1OZG9FNjNzAS72jcXNhrMdI5Rkp1Hs02UaNsK+8hzmA/HN5HKmY1gFEsugDspG2GH5hjMqkNgUof9ZPFdHPqPZPzbFwnMX77KrVIHEUtRB2YjaEh8pbonbgzIyMUvb4CR7VSCxYbYXeYMilzg1Lsan57jYN6bhvU1QlR9c2PHE1aG4XH96bp7Gbj/7K+I7RcGOqIOyEakeF3WlvrgN2GoGic3jdgn7K3I4cXU4Ltc/3xUUSGjvd+OICDdtzeP4lfg4qMZuvwokVkAdlM3YX5HLqbbhuAglFrIQaCW3Keqr8znbMcLkTOwH20Prd6lMeXMcrs7jUv84ff7pmF9bBRIrow7KZhyuzmN0ao7GOAzYnmofpiI3g7yspFvYOCIcrs5jLmA40Rb7lvjZjhGKfGmUZKtAYjPctDUfgJ/HoRd1tmOEnIwUqvIzYn5tu6MOymYc2RZ8UF6+NBjT6xpjePnS4ML1lY1z09Z8ROCVS7Gv5E61D2sLPAz2VWST5nFx/HJsnzsIpjjarxkklkUdlM2ozMukPCedl2P8oLT2j9M/NsPN6qA2TU5GCnUlPo5fiW3Z9Y9Nc7FvnMPVWnabJc3j5oaqXF6JcQ9qanaeph6/iltWQB2UDTmyLZ+XLw0SXJEkNhxrHVy4trJ5Dlfn8+qVIeZiOGE31Nu+ebuWXTjUb83jXIzHEF+9OsTsvOFwdV7Mrukk1EHZkMPb8unzT3NlIHYTdl++NEChN41thVkxu2YicnhbPuMz8zR0xW4M8WjrAJmpbg3xhcnh6nzmAoaTbbFTYh5tHcQlwd+Ncj3qoGzIzTEehzLGcOzSIDdvz9c4eJiEWsKxDNEeax3kpq15pLj1cQ6HG7fkIUJMx6GOXhxgX0UO2emauWU59BdtQ3YUecnPSo1ZJdc+NEnXyJSOP0WAspwMKvMyYlbJDY7P0Njj55Y4r2WUCORkplBb7IvZfKjJmXlOtg1r2a2COigbIiIcrs6LWQ/q2CUdf4okR6rzeeVybMYQX740AKCNiwhRX53Hq1eGYrKy9YmrQ8zMB7hVHdSKqIOyKYer87k6OEH3yFTUr3WsdYBcq/WohE99dT79YzNc6h+P+rWOtg6SnuLiQKWmyYkEh6vz8U/P0dgd/THEl1oHcLuEehVIrEhYDkpE8kXkxyLSbP297j8tIveIyMlFrykRebv13b+JyKVF3x0Mx55E4uZtwVZVLMJ8L18e5HB1Pi6Xjj9FgiPbgo/B8cvRDxUdbR3gpq15pHq0rRkJbtpqlV0MpgocbQ2OP/l0/GlFwv1VPwI8ZYypAZ6yPr8GY8wzxpiDxpiDwOuBCeBHi3b549D3xpiTYdqTMOwu85GV6uaVKIf5ukemuDIwoSGiCLKjyEteZkrUGxfDE9b40zYNEUWKyrwMSrPTo964uDb+pM/daoTroB4Evmi9/yLw9jX2fxfwA2NM/BY8cgget4ubqvOjPg51bGEMQyu5SCEi1FfnR10oEZwrBzfrGEbECJZdXtTL7udXgvOfdPxpdcJ1UCXGmC4A62/xGvs/BPzHkm1/KSKnReRTIpK20oEi8gEROS4ix/v6+sKz2iEcqc6jscfP8MRM1K7x8qVBvGkedpfp+FMkOVydx+WBCXr90RtDPNo6SJrHxQ1VOv8pktRvzaNzZIqO4cmoXeOl1n5r/El7UKuxpoMSkZ+IyNllXg9u5EIiUgbsB364aPOfAruAw0A+8CcrHW+MedQYU2+MqS8qKtrIpR3LEatX80oUww0vXxqkvjoPj86hiSihtEPRDBUduzTAjVvySPO4o3aNZKR+oeyi14s62jrI/oocvGmeqF0jEVizVjLGvNEYs2+Z13eAHsvxhBxQ7yqn+mXgW8aY2UXn7jJBpoF/BY6EdzuJxYHKHFLdLl6J0oMyMDZNc++YysujwL6KHNJTXFEL0Y5MzHK+a1TTG0WBXaU+vGmeqDUuJmbmONU2zK07NLy3FuE2m58AHrbePwx8Z5V9382S8N4i5yYEx6/OhmlPQpGe4uaGqpyFeUqRJuT4VCAReVLcLg5V5UVNDRacZ6Vjh9HA43ZxaEtu1BoXxy8PMRcwOkF3HYTroD4B3CsizcC91mdEpF5EPh/aSUSqgSrgp0uO/4qInAHOAIXAX4RpT8JxZFs+5zpGGJ+ei/i5j10KzqHRpaajw+Ft+ZzvHGVkcnbtnTfIsUsDpHqCFakSee6oKaSxx0/bYOT1XEdbB/C4hPqtOv9pLcJyUMaYAWPMG4wxNdbfQWv7cWPM+xftd9kYU2GMCSw5/vXGmP1WyPC9xpixcOxJRI5sK2AuYKIS5jvWOsiNW3QOTbS4u66IgIGnGnoifu6jrYMcrMolPUXHn6LB/fvKAPivs90RP/dLrQMcqMwhS8ef1kRrJptz87Z8slLdEX9QBsamaege1RBRFDlUlUt5TjpPnumK6HlHp2Y51znCLRqajRpV+ZnsLc/mB2cjW3bj03Ocbh/R8ad1og7K5qSnuLl3Twn/da6b2QiuMfTkmS6MgTfvK4nYOZXXIiLct6+M55r6GZ2KXJjvZy39BAzcopVcVLlvbymvXh2mZzRyUwWOW3n+dPxpfaiDcgBvPVDO8MQsL7T0R+ycT5zqpLbEy67S7IidU7meXzhQysx8gKcbVhO4boxvn+ik0JvKEZ1DE1Xu318KwA/PRS568VRDD6ke10JKJWV11EE5gDtqC/Gle/jeqciEGzqGJ3nl8hAP3FAekfMpK3OoKo/S7HS+H6Ew38jkLE9f6OVtN5Tr3LUos7PYx85iLz84ExkHNTMX4LunOnnTnhIyU3X8aT3oL9wBpHncvGlPKT863830XPjLUX/vVCcAD9xQEfa5lNVxuYT79pXy06Y+/BEI8/3gTBcz8wHecUjLLhbct7eUY5cGGBibDvtczzT2MjQxyztvrIyAZcmBOiiH8NYbyvBPzfFcU/hhvu+c7ORgVS5bCjIjYJmyFr9woIyZuQBPXwg/zPetEx1sL8zS5d1jxH37SgkY+EkElJjffLWdQm8ad9QURsCy5EAdlEO4fWchuZkpfP90Z1jnaekd43zXqIb3YshNW/IoyU4LW83XMTzJsUuDvP1QBcG57Uq02VueTVV+Bj8IU0U7ND7D0xd6eftBDc1uBP1POYQUt4v79pby4/M9TM1uPsz3xKlOXAJvPVAWQeuU1XC5hPv3lfFsY19YE66fOBlsnDx4UBsXsUIkWHYvtvSHNeH6e6c7mZ03/KKG9zaEOigH8dYD5YzPzPNs4+ZCRcYYvnuqk1u2F1CcnR5h65TVuH9fKdNzAZ4KI8z37RMd3Lgll60FWRG0TFmLN+8tZXbe8EwYZff1VzvYVepjT7mqZjeCOigHccv2fAqyUvnu6c2Fis52jHKpf1zDe3GgvjqfIl8aT26y7Bq6Rmns8fN2FUfEnENVuZRkp2160m5L7xin2oZ5103ae9oo6qAchMft4r59pTzd0MvEzMZDRU+c6iDFLQtpXJTY4XYJ9+8r5ZnG3k2F+b59ogOPS/iF/Vp2scblEu7bG1Ribua5+9aJdlwCD2hodsOog3IYbz1QzuTsPE9tcOJnIGD47qku7qotIiczJUrWKatx/74ypjeh5gsEDN852cmdtUUUeFdc01OJIvftK2NqNsC3T2xMpBQIGL71agd31hZR7NOw+kZRB+UwjmzLp9iXxuOvtGGMWfdxL18epHt0irdpeC9uHNmWT2VeBp95qpmZufWnrTp6aYDu0SkN78WRW7bnU781j7/5UeOGxBJHWwfoHJlSccQmUQflMNwu4bfu2sELLf08cWp9rbn5gOFvftSIL93DvXs09168cLuEj7xtL829Y/zLC5fWfdxXj10lK9XNvbu17OKFiPCRB/YyODHD3/6ked3HfePVDnxpHt6kz92mUAflQH79ddUcrMrlo989v64Z7p9/vpVXLg/x0Qf2aoqVOPPGPSXcu6eEv32qaV1rDX33VCffO93Fb96+jYxUXVojnuyryOGhw1v44kuXae7xr7n/q1eHeOJUB287WK7LomwSdVAOxO0S/ve7DuCfmuVj3zu/6r6N3X7+5kdNvHlviabHsQkfeWAvgvDR755bdb+2wQn+1zfPcGhLLr/7hpoYWaesxh+9qZasVDcf/e75VUPsA2PTfPArr1KSnc6fvHlXDC1MLNRBOZTaEh8fvGcn3znZueKCeDNzAf7g8ZNkZ3j4q3fs1+wDNqEiN4Pff2MNP2no5UcrZMqenQ/wu4+dAOAzDx0iRbMP2IICbxp/eG8tL7T086Pzyz938wHD7z12koHxGT733ptUlBQGYf3qReSXROSciAREpH6V/e4TkUYRaRGRRxZt3yYix0SkWUQeF5HUcOxJNn7n7p3Ulfj4s2+dXTYR6d893cz5rlH+6h37Vf1lM37z9m3Ulfj4yBPnlpWd/+1PmjlxdZi/+sX9VOVrzkQ78d5btlJb4uXj3zu/bFaXT/24iRda+vn4g3vZpzkTwyLcZtlZ4BeB51baQUTcwD8A9wN7gHeLyB7r608CnzLG1ABDwPvCtCepSPW4+OS7DtDrn+IvvtdAx/AkQ+MzTM/Nc+LqEJ999iLvuqmSN+0tjbepyhJS3C7+4h376ByZ4qPfPcfJtmEGxqYxxvCzi/38w7Mt/HJ9paoubYjH7eLDb9tL+9Akf/z10/zkfA/9DD+Z8AAACMNJREFU1ljwUw09/P0zwbL7lcNb4myp8wlrxNwY0wCsFTo6ArQYY1qtfR8DHhSRBuD1wK9a+30R+Ajwj+HYlGwcrMrlN2/bxudfuMTjx9te811FbgZ//rY9KxypxJvD1fm895YtfPnoVb52vB2AzFQ3xsC2wiw+8sDeOFuorMRtOwt5+NatfPnYVb5rqWkr8zIYnphlb3k2H3twX5wtTAxiIemqABbXnO3AzUABMGyMmVu0XUfxN8Ej9+/iyLZ8hiZmGJ+eZ3J2nunZeR44WE52usa/7czHH9zHr91STdvgBG1DE1wdnGBkYpb/fvcOVVzanI8+uI9H7t/NmY4RTrYNcbJtmO6RKT79K4dUtRch1nwCROQnwHIxoj8zxnxnHddYrntlVtm+kh0fAD4AsGWLdp0X43G7NIznUESEulIfdaW+eJuibIKMVDdHtuVzZFt+vE1JSNZ0UMaYN4Z5jXagatHnSqAT6AdyRcRj9aJC21ey41HgUYD6+vr1p1BQFEVRHEkstKuvADWWYi8VeAh4wgQnETwDvMva72FgPT0yRVEUJQkIV2b+DhFpB24Fvi8iP7S2l4vIkwBW7+hDwA+BBuBrxpjQDMU/Af5QRFoIjkn9Szj2KIqiKImDbCThqF2or683x48fj7cZiqIoyjKIyM+NMSvOjV0vOj1dURRFsSXqoBRFURRbog5KURRFsSXqoBRFURRbog5KURRFsSWOVPGJSB9wJYxTFBKcKOxk9B7sgd6DPXD6PTjdfnjtPWw1xhSFe0JHOqhwEZHjkZBAxhO9B3ug92APnH4PTrcfonMPGuJTFEVRbIk6KEVRFMWWJKuDejTeBkQAvQd7oPdgD5x+D063H6JwD0k5BqUoiqLYn2TtQSmKoig2Rx2UoiiKYksSwkGJyH0i0igiLSLyyDLfp4nI49b3x0SketF3f2ptbxSRN6/3nA65hy+ISK+InHWa/SJSJSLPiEiDiJwTkd9z4D2ki8jLInLKuoePOu0eFn3nFpETIvI9J96DiFwWkTMiclJEor4UQpTuIVdEvi4iF6zn4lYn3YOI1Fn//9BrVER+f1UjjDGOfgFu4CKwHUgFTgF7luzzO8DnrPcPAY9b7/dY+6cB26zzuNdzTrvfg/XdncCNwFkHlkEZcKO1jw9ocloZAAJ4rX1SgGPALU66h0XH/SHwVeB7TvstWd9dBgqjaXsM7uGLwPut96lArtPuYcn5uwlO6F3RjkToQR0BWowxrcaYGeAx4MEl+zxIsHABvg68QUTE2v6YMWbaGHMJaLHOt55z2v0eMMY8BwxG0e6o2W+M6TLGvGrdh5/gYpcVDrsHY4wZs/ZPsV7RVCVF5XckIpXALwCfj6LtUb2HGBPxexCRbIINzn8BMMbMGGOGnXQPS459A3DRGLNqRqBEcFAVQNuiz+1cX5Et7GOCK/yOEFzBd6Vj13POSBKNe4glUbXfCh0cItgDiRZRuQcrNHYS6AV+bIxx3D0Anwb+JxCIvMnXEa17MMCPROTnIvKBKNi9rH3L2HHdPuu8h+1AH/CvVqj18yKSFR3zX2vfEjuW3WcTddJDwH+sZUQiOChZZtvSVupK+2x0e7SIxj3EkqjZLyJe4BvA7xtjRjdt4dpE5R6MMfPGmINAJcGW8L6wrFydiN+DiLwV6DXG/Dxc49ZJtH5LtxljbgTuBz4oIndu3sQ1icY9eAiG6//RGHMIGAeiOTYezWc6FXgA+M+1jEgEB9UOVC36XAl0rrSPiHiAHIKhr5WOXc85I0k07iGWRMV+EUkh6Jy+Yoz5ZlQsX8a+pXYst89Gy8AKxzwL3BdJo1eybyU72Pg93AY8ICKXCYZ5Xi8iX46G8UvtW2LHsvustxyMMaG/vcC3iG7oL1p1UvuiHvjXCTqsaBHN5+F+4FVjTM+aVkRrkC1WL4Iti1aCg3Ghwby9S/b5IK8dzPua9X4vrx3MayU4eLfmOe1+D4uOqyb6IololIEAXwI+7eDfURHWQDaQATwPvNVJ97Dk2LuJvkgiGuWQBfisfbKAnwH3OekerO+eB+qs9x8B/tpp92B9/xjwG+uyI5o/tli9gLcQVHldBP7M2vYx4AHrfTrB7mQL8DKwfdGxf2Yd1wjcv9o5HXgP/wF0AbMEWzXvc4r9wO0EwwKngZPW6y1OKgPgAHDCuoezwJ878Xe06Pu7ibKDilI5bCdYYZ4Czjn4eT4IHLd+T98G8hx4D5nAAJCzHhs01ZGiKIpiSxJhDEpRFEVJQNRBKYqiKLZEHZSiKIpiS9RBKYqiKLZEHZSiKIpiS9RBKcoGEZGCRRmZu0WkY9Hnn0XpmodEZMVceCJSJCL/FY1rK0q88MTbAEVxGsaYAYJzUhCRjwBjxpj/L8qX/V/AX6xiU5+IdInIbcaYF6Nsi6LEBO1BKUoEEZEx6+/dIvJTEfmaiDSJyCdE5D0SXB/qjIjssPYrEpFviMgr1uu2Zc7pAw4YY05Zn+9a1GM7YX0Pwcmb74nRrSpK1FEHpSjR4wbg94D9wK8BtcaYIwSXrfi/rX3+FviUMeYw8E6WX9KinmAmihB/BHzQBJPQ3gFMWtuPW58VJSHQEJ+iRI9XjDFdACJyEfiRtf0McI/1/o3AnuAyOgBki4jPBNfAClFGcKmFEC8C/0dEvgJ80xjTbm3vBcojfxuKEh/UQSlK9Jhe9D6w6HOAa8+eC7jVGDPJykwSzHsGgDHmEyLyfYK50o6KyBuNMResfVY7j6I4Cg3xKUp8+RHwodAHETm4zD4NwM5F++wwxpwxxnySYFhvl/VVLa8NBSqKo1EHpSjx5XeBehE5LSLngd9euoPVO8pZJIb4fRE5KyKnCPaYfmBtvwf4fiyMVpRYoNnMFcUBiMgfAH5jzGpzoZ4DHjTGDMXOMkWJHtqDUhRn8I+8dkzrNYhIEfB/1DkpiYT2oBRFURRboj0oRVEUxZaog1IURVFsiTooRVEUxZaog1IURVFsiTooRVEUxZb8/4BYYA2g6P7VAAAAAElFTkSuQmCC\n", "text/plain": [ "