{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## ThinkDSP\n",
"\n",
"This notebook contains\n",
"\n",
"Copyright 2015 Allen Downey\n",
"\n",
"License: [Creative Commons Attribution 4.0 International](http://creativecommons.org/licenses/by/4.0/)"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false,
"jupyter": {
"outputs_hidden": false
}
},
"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": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"\n",
"from thinkdsp import decorate"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Do we hear phase?\n",
"\n",
"This notebook investigates what effect, if any, changes in phase have on our perception of sound.\n",
"\n",
"I'll start with a simple waveform, a sawtooth, and move on to more natural sounds."
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false,
"jupyter": {
"outputs_hidden": false
}
},
"outputs": [
{
"data": {
"text/html": [
"\n",
" \n",
" "
],
"text/plain": [
""
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"from thinkdsp import SawtoothSignal\n",
"\n",
"signal = SawtoothSignal(freq=500, offset=0)\n",
"wave = signal.make_wave(duration=0.5, framerate=40000)\n",
"wave.make_audio()"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false,
"jupyter": {
"outputs_hidden": false
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO29eZQc2XWf+d3aF9SC2rDvO9BNNkg0emOvWJrS+LA5Y8kixwtli8PxjCQvOvKIGo0lmZbOUGMfU+MjjS2ORIu2ZZE0LY/aFuUKAN1oNNnobqDZGyILhSrsW0TWvq9Zb/7ISCBRyKwtM/b3nZOnMiMjM1/eing37+/duFeUUmg0Go1GEzRK/B6ARqPRaDS50A5Ko9FoNIFEOyiNRqPRBBLtoDQajUYTSLSD0mg0Gk0gKfN7ACuhpaVFbd261e9haDQajaYIvPfee71Kqdb520PpoLZu3cr58+f9HoZGo9FoioCIXM+1XUt8Go1Gowkk2kFpNBqNJpBoB6XRaDSaQKIdlEaj0WgCiXZQGo1Gowkk2kFpNBqNJpAUxUGJyLdEJCkiF/I8LyLyL0WkW0Q+EpFPZT33JRHpcm5fKsZ4NBqNRhN+ihVB/THw2QWe/wlgl3P7CvCvAESkCfgN4AngMPAbIrK6SGPSaDQaTYgpioNSSp0B+hfY5RXg36o0bwONIrIOeBk4oZTqV0oNACdY2NFpisx/eOcG568t9K+LF3/8o6t8fGvI72EEhj944zKd1ojfwwgESil+77UurvaO+T2UQDCTmuPL3z7PD7t6XfsMr9agNgA3sx7fcrbl2/4QIvIVETkvIud7enpcG2icGBqf4R//+QW+e+7m4jvHgLtDE/zmf0nwZ+/f8nsogaA7Ocr/+ZcX+a8f3fF7KIHgo1tD/HPjEoZp+T2UQHDuaj8nO2xGp2Zd+wyvHJTk2KYW2P7wRqW+qZQ6pJQ61Nr6UMkmzQp4rdMmNadyGzyGnEzYAOgm02lOaHs8gJFIOyZtjjRGwqayrITndre49hleOahbwKasxxuBOwts13iAYeoJKBvDmZA1ae5PyPoAAX2+ZKOUwjAtnt3VSk2FeyVdvXJQrwJ/y8nmexIYUkrdBdqB4yKy2kmOOO5s07jM5EyK051pqVRPQGm58+zlPiB98sUde3iS928MAnpCBrjSM0pXchTQ5wvAhdvD3Bma5OUDa1z9nKK4PhH5U+AFoEVEbpHOzCsHUEr9a+AHwE8C3cA48Led5/pF5J8C55y3+ppSSq/Ye8APu3qZmEmlH+jzjdc7k8zOpQ2hzXFf3gNtD3gwutYOOx1dlwgc2RcCB6WU+uIizyvg5/M89y3gW8UYh2bpGAmLuqoy6irL9ARE2h5tdZXMpOb0BER6Qt7WUsvtwQltD8AwLfaureOizmgEoN20eHxrE021Fa5+jq4kEUNSc4qTHUle2ttGeVlJ7CWtjNx5bP8aSkSYi7k9hidnOHu5l+P71yBoSSs5PMn7Nwc5fmAtAHNz8bbH1d4xLtmj9+zhJtpBxZDz1/rpH5vm+P61zgQUb37U3cv4dIrjB9Yiou3x+sUkMynF8QNrEH2AcKLDRinurbfE3ByccJJnju93V94D7aBiiZGwqSgr4fk9rYgIMf9BiGHa1FWW8dT2ZkBiL2kZCZuWVZUc3LQaQUeUhmmzpbmGvWvrAb0GZZg2+9fVs6mpxvXP0g4qZiilMBIWn9nZwqrKsnQEFeMzLi132ry4t42KspJ0xBDj38hTsylOX0ym5c4SSUeU8TUHI5MzvHW5l5cPrL130WacJc+ekSneuzHAcZez9zJoBxUzLloj3OyfuB+ex1zS+vGNAfrGpu+dcGmH7e+Y/OSty32MTacetIe/Q/KV0509ablzvyN3Eu/j4+Q9udP99SfQDip2tJsWkpUeGvOAgfYLFhWlJTy/O12dJO4Rg2Fa1FaU8vSOZgBE4i15tpsWLasqOLh5NeJ4qBibA8O02NRUzd61dZ58nnZQMcMwbT69eTWtdZWAMwHF9JRLy502T+9spq6qHAAhvvZIzSlOJGxe2NtGZVkpQKyz+KZm09mdR/etobQkqypbTD326NQsP+ruSydXSa4qdcVHO6gYcbN/nMTd4QfC8zhLWp32CDf6xx+0R4wjqA9uDtA7Ov2gfBNje5y93Mfo1OzDx4ePY/KT051JplNznmTvZdAOKkZkroY/lnWAxXlCbr9gO3Jn271tcV5zaTdtykuFF/bcL8bsze/kYNJu2tRWlPKUI3dCvH/QGaZNU20Fh7Y2efaZ2kHFCMO02LOmjq0ttfe2xVnSMhIWn9q8mra6qnvb4rrmopSi3bR4akcL9Y7cCRl7xM8gcxm5c08bVeWl97bHVRKfnp3j9YtJju5re1DudBntoGJC/9g05671P5QeGtcI6tbAOOad4ZxyRRwnoK7kKNf7xh+yR1wlrfdvDtI7OvXw+UI8z5e3r/QxMjXL8f3eZO9l0A4qJpzqsJnLkx4aw/PtXjHU+eVa4lo5IdOE7yEHRTwnZCNhUV4qvLi37YHtcXXY7aZFTUUpn9nlXu+nXGgHFRPaTZv1DVUcWF//wPa4SlrtpsXuNavYliV3QpwnIJuDmxtpq696YHscJa10ryObJ7c3PyB3giOJx8sc9+TO53e3PiB3eoF2UDFgfHqWN7t6nFpzD+rHcbyyY2Bsmnev9ueUK9ITULzscWdwgo9vD+WxR/wiqO7kKFd7x3IXQ5X4ScAf3hokOfKw3OkF2kHFgDOXepmanct5gMVxDerUxSRzivz28GFMfpKRO3M1n4ujPe5lu+bodRTHC9uNhE1ZifDSHu2gNC5gJCwaqss5nCM9NJYTkGmxrqGKRzc0PPRcHCMGI2Gxs20V21tX5Xg2fpKWYVo8tqmRtQ1VDz0Xx/Ol3bR4cnszDTXli+9cZLSDijgzqTlOdSQ5sq+NstKH/91xk7QmplOc6epxaqs9nC6bXnOJD4Pj07x9pT/vxZdxK557d2iCD28N5ZWz4na+dCdHudIz5ou8B9pBRZ5zV/sZmpjJmx4at1+EZ7p6mJyZy9tsLW7V3V+7mCQ1pxaxh7dj8pN72Z0LnS8xsofh9H466nJr93wUxUGJyGdFpFNEukXkqzme/4aIfODcLonIYNZzqaznXi3GeDT3MRI2VeX3i6HOJ24TkGHaablzW56r4WPmsA3TZm19FZ/IIXdCDCdk02ZHay0723LJnfGrNNJu2nxiYwPrG6t9+fyyQt9AREqB3weOAbeAcyLyqlIqkdlHKfUPs/b/ReBg1ltMKKUeK3QcmodJp8taPLurleqKPOmhMZK0ZlNznLpoc2RvG+U55E6I1yL45EyKNy718FOf3khJnuoAcao0MjQ+w9tX+vifntued584XZZhDU3y4c1B/tHLe3wbQzEiqMNAt1LqilJqGvgO8MoC+38R+NMifK5mES7cHubO0OSCxR3jJGm9e62fwfGZBfX0OF3382ZXLxMzqUXsEZ8I6rVOm9k5tfj5EpPj40RHRu70R96D4jioDcDNrMe3nG0PISJbgG3Aa1mbq0TkvIi8LSKfz/chIvIVZ7/zPT09RRh29DESFiVZvZ9y4VHV/EBgmDaVZSU8l0fuhHhJnoZpUVdVxhPbmvPuEydJyzBt2uoq+eTGxvw7xchhG6bFtpb8cqcXFMNB5Zri8v0LvwB8XymVytq2WSl1CPgfgd8VkR25XqiU+qZS6pBS6lBra/4JRnMfw7Q5vK2JptqKvPvEZUJWKn01/LO7WqmpyK9sxyVimE3NcbIjLXdWlOWfBuIiaWXkzuMH1uSVOyE+1d2HJmY4e7kvb7arVxTDQd0CNmU93gjcybPvF5gn7yml7jh/rwCneXB9SrNCrvWO0WmPLFrcMS6SlnlnmNuDE4umy8ZlzeX89QEGxmfyZu9lEwd7/LCrl/Hp1NLOlxh47NOdybTc6VFr93wUw0GdA3aJyDYRqSDthB7KxhORPcBq4GzWttUiUuncbwGeARLzX6tZPpn00GOL6MdxiaAM05E75xX/nE9cIijDtKlYRO6E+BTPNRIWdZVlPLk9v9wJ8bkswzBtWlZVcnDTAnKnBxScxaeUmhWRXwDagVLgW0opU0S+BpxXSmWc1ReB76gHf37sA/5AROZIO8uvZ2f/aVaOYdrsX1fPpqaaBfeLzYScsDm0tYnmVZWL7ht1c6Rb3Vt8ZmcLqyoXngLiMCGn5hQnO5K8uIjcCfH4QTc5k+J0Z5LPPbZhQbnTCwp2UABKqR8AP5i37dfnPf7NHK97C3i0GGPQ3KdnZIr3bgzwD47sXnTfOEha1/vGuGiN8I//yv5F943DmkvH3RFuDUzwiy/tXHTfOFROeO/6AP1j0zlb0cwnDpL4W5d7GZteOLvTK3QliQhyssNG5SmG+hAxiKAMc+npsnGo7t5uWk6r+yXYIwYRVLtpUVFawvN7Fk++ikMEZZg2qyrLeHrHwnKnF2gHFUEM02JTUzV719Ytum8clhiMhMW+JcidEA/J00jYHNqympYlyJ1Rn5AzcuczO5sXlTsh+g47LXfavLCnlcoyb3s/5UI7qIgxOjXLj7r7OL7/4d5PuYj6Injv6BTnrw8s+WLDqE9AN/vH6bg7vCQ5C6JfPPeiNcLN/ollZKtFWwJ+/8YAvaPTvmfvZdAOKmK80dnDdGpu6RNQxNegTjly57LsEeEZ6F6vo6U6bKJdacQwbUSWXgw16tXd202n1f0S5E4v0A4qYrSbFk21FXx6y+ol7R91SavdtNm4upp96xaXOyH6EVS7abF3bR1bmmsX3xkiXzy33bT49ObVtNYtLndCtCXPtNxp8/SOFuqqvO/9lAvtoCLE9Owcr19McnRfG6VLTA+N8oQ8OjXLD7t7lyx3QrQnoL7RKc5fy9/7KRdRDhhu9o+TuDu8rGy1KP+gu2SPcr1vPBDZexm0g4oQb1/pY2RqdslyFkRb0jpzqYfpPK3u8xLhNZf7re6XcXxEOK16sd5PuYiyJJ7J7szV6t4vtIOKEO2mRU1FKc/sbFnya6IcQWXkzkNLlDsh2msuhmmxobGaA+vrl/yaKEeU7abFnjV1bG1ZotxJtCMoI2FxcFMjbfUPt7r3C+2gIsLcXLoY6vO7W6kqX156aBRPuOnZOV67mOTI3tyt7vMR1eruY1OznOnq5dgyi39GdULuH5vm3LX+ZctZUU16vT04wYXbw4HJ3sugHVRE+PDWIMmRqeWfcBGVtN652sfI5OyyT7ioRgxvdq1A7iS6ktapDjstdy5D3oPoVho5YaZrd/rZ+ykX2kFFBCNhU1YivLRn+b8Io3jGGaZNdXkpz+5autwJ0V1zMUybxppyDm/N0+o+D1GNoIyEzfqGKh7ZsHS5M0MUj49202Zn2yq2t/rX+ykX2kFFhHbT4sntzTTULC89NIprUHNz6eoAK5E7oxhBzdzr/bRmWXJnhoiZg/HpWc5c6uH4gaVnd2aI4oXtA2PTvHutn5cDlL2XQTuoCNCdHOVKz9iK0kOjOCF/dHsIe3j5cidEM2J492o/w5OzK7RH9CStM5d6mZqdW5GcFcUfdK9dTJKaU8uWO71AO6gIsNTeT7mIoqRlmBalJcJLi/R+ykUU11wM06KqvITndi2/OkAUi+caCYuG6nIe37Y8uROieVmGkbBYW1/Foxsa/B7KQ2gHFQEM0+aTGxtY11C97NdGMYIyEjZPbm+isSZ/q/u8RCyCylQHeG5XK9UVyy/+GbWIcjY1x6mOJEf2tVG+ArkzahHUxPTSWt37hXZQIccamuSDm4MrTg+N2gR0uWeU7uToiuWKqC0xfHx7iLtDk4UdH0Uek5+8e7WfoYmZwo6PCBnkza4eJmfmAinvgXZQoedEx9J7HeUmWmnmmd5PK5E7IXqL4IZpU1oii7a6z0fUJC0jYVNZVsJzu5eX3ZkhapdlGAmbuqoynti+fLnTC4rioETksyLSKSLdIvLVHM//rIj0iMgHzu3LWc99SUS6nNuXijGeOGGYFttaatnZtrL00HQEFZ1TzkhYPLqhgfWNy5c7IXprUEbC4vDWJlbXrkDuJFoRlFIKw7R4dlcrNRUrayYepUojabnT5sjelcmdXlDwqESkFPh94CeA/cAXRSRXb+3vKqUec25/6Ly2CfgN4AngMPAbIrL0ujQxZ2hihrOX+zh+YHnVAbIJnuq8cpLDk7x/Y7CgdNkoSZ5Xe8e4ZI8WVPwzSpKWeWeYO0OThaVTR8hhn7s2wMD4zLJqd3pNMdzmYaBbKXVFKTUNfAd4ZYmvfRk4oZTqV0oNACeAzxZhTLHgdGeS2QLTQ6M0IWd6HRVSriVKEYNhrjy78x4RkrTaTYuSJba6z0eUqrsbCYuKshKe2x2M3k+5KIaD2gDczHp8y9k2n78qIh+JyPdFZNMyX6vJgWHatNZVcnBT44rfI0qSlpGw2dpcw64Vyp0QrTUXI2FzYH09G1cv3uo+H1GStAzT5vGtTTStUO6E6FyWkZY7bZ7d2ULtElrd+0UxHFQulWj+f/C/AFuVUp8ATgLfXsZr0zuKfEVEzovI+Z6enhUPNipMzqQ43Znk2P7C0kOjEkENT85w9nLviqoDZBOVCCo5MsmPbwwULN9EpXjutd4xOu2RgouhRkXyNO8Mc3twIlC9n3JRDAd1C9iU9XgjcCd7B6VUn1Jqynn4/wKfXuprs97jm0qpQ0qpQ62twQ1JveKty72MTacKLu4YlQn59YtJZlKqKOVaojABnUwkUYqCJ6CoTMiZi9mLcr5Ewh42Jctode8XxXBQ54BdIrJNRCqALwCvZu8gIuuyHn4O6HDutwPHRWS1kxxx3NmmWQTDtFlVWcZTO5oLep+oSFpGwqZlVSWPbSosx6YkImsuRsJic1MNe9YsrdV9PqIiaRmmzf519WxqWrncCZnjIwr2sDi0pYnmVUtrde8XBTsopdQs8AukHUsH8D2llCkiXxORzzm7/T0RMUXkQ+DvAT/rvLYf+Kekndw54GvONs0CpOYUJztsXtjTSmXZ8qsDPEAEIqip2RSnLyY5tn/pre7zIREIGUYmZ3iru4/jy+z9lIsImIOekSneuzFQNDkr7Pa40TfORWsk8PIeQFFWx5RSPwB+MG/br2fd/1XgV/O89lvAt4oxjrjw/o0Bekeni5IeGoWspLcu9zlyZ3HsEXJz8MalHqZTc7z8SBHsEQFJ61SHnZY7i3F8RCDCvi93Bje9PEMwr87SLEi7aVFeKrywp/C1uEiccKZFbUUpT+8sTO6EaFTvbjdtmmsr+NTmwi8pjEKWZ7tpsampmn3rCpM7IRoRpWHa7F1bx+bmwuROL9AOKmRkin8+vaOFuqrl9X7KRdjTiFNOq/sX9rYVLneStsdciO0xNZvi9YtJju5bU7DcCYS+eO7o1Cw/6u7j+P7CsjszhL3ySu/oFOev9weutXs+tIMKGZfsUa73jRdNPw57Ft8HN9NyZ7FaVYdd0nr7Sj+jU7O8/EiR7EG4j483OtNyZ1GPj6K8kz+81pF0Wt0Hf/0JtIMKHYZpIVJgdYAswi5ZGKZNeanw4gqLoT5MuCVPw7SoqSjl6R0rK4Y6n7AXzzUSFk21FRxaZqv7fIQ967XdtNjQWM2B9ctvde8H2kGFjPaExcFNjbTVVRXl/cKcRqyUot20eGpHC/VFkDsh3BLOXEbu3LP8Vvf5CPMa1PTsHK9dTHJ0X+HZnRnCHEGNTc3yZndvQbU7vUY7qBBxe3CCC7eHi6ofhzmC6kqOcq1vvKhyRThO29x8cGuQ5MhUUbOzwix5vn2lj5HJ2eLag/Da48ylHqZng9v7KRfaQYWIE2ZxroZ/gBBPQEUphjqPME/IhmlTVlJMuTPcEYORSMudn9lVHLkTCHXxXCNhs7qmnMe3hqdhhHZQIaLdtNnVtortrSsvhjofCXHMYCRsDm5uZE19ceROCLekZSQsntrRTEN1ceROCO+aS0bufH538eROCG/W60ym99O+NZQFtPdTLsIz0pgzMDbNu9f6i371d1jXXO4MTvDRraGiyxVhjaC6kyNc6RkrenZWWCOoD28NYg9PuXK+hJF3rvQzPDkbmuy9DNpBhYTXLiZJFdj7KRdhTdI6ca/3k56QIR1dAxxzYX0hjA7bSKRb3b+0p8jHB2G1h0VVeQnP7gpXoW3toEKCkbBYW1/FJzY2FPV9wxoxGAmLHa217Cii3AnhlbSMhM0nNzWytqF4cieEt9KIYVo8ub2JhpriyZ0QzqzXubl076fndrVSXVE8udMLtIMKARPTKd641ONKemgY11wGx6d5+0q/O62qQxhBWUOTfHhz0BX5Jl2rMVwW6U6OcrlnzJXjI4wR1Me3h7CGJwPd2j0f2kGFgDe7epiccSc9NIwR1D2506UJKGwe6oRT/LMYvbDmE0bJM1MM1Y1eR2E8X4yElZY7i5jd6RXaQYUAI2FTX1XGE9uLczV8NqGcgEybNfWVfGJDceVOCKekZSRstrcUX+6EcEYMhmnziY0NrG+sLvp7h1FxMEybw1ubWF1Aq3u/0A4q4MxmpYeWu5IeGq7q3ZMzabmz0Fb3+QhbGvHQxAxnL/cV3Oo+H2Fbc7GHJ/nAJbkTCN11g1d6RulKjroSXXuBdlAB59y1AQbGZ1w74UpCpmm92dXLxEzKNT09bBHl6xeTzM4p15rPlYRsQjac7E7Xjg/CdXxk7HEshOtPoB1U4DESFhVlJTy325300LBp6oZpUVdVxhPbCu/9lIuwSVpGwqKtrpLHNja69AnhirAN02JbSy0724ovd0L4iucapsUjG+rZ4ILc6QXaQQUYpdLpoc/ubKG2sijNjx9CQlS9ezY1x8kOm5f2tlFR5s6hGyZJa3ImxelO9+ROCFdEeU/uLEKr+3yEaQ0qOTzJ+zcHQ1V7bz5FOctF5LMi0iki3SLy1RzP/5KIJETkIxE5JSJbsp5LicgHzu3VYownKph3hrk9OOFqemiYKkm8dz0jd7poD8ITQb11uZfx6ZSrzefCtCZ3utNduRPCpTiccFrdhzG9PEPBP8tFpBT4feAYcAs4JyKvKqUSWbu9DxxSSo2LyP8C/F/AzzjPTSilHit0HFHESNiUCBzZ5156aJgUi3bTpqKshOeL0Oo+LyGagNov2NRVlvHUdnfkTghXaR/DtGlZVcnBTe4VQw1TRGmYNluaa9i9xh250wuKEUEdBrqVUleUUtPAd4BXsndQSr2ulBp3Hr4NbCzC50Yew7Q4tKWJ5lWVrn2GSDjWGNKt7i0+s7OFVS7JnRCe4rmpOcXJjnSre7fkTshU1nDt7YtGWu5Muip3QngqjYxMzvDW5V5X5U4vKMaRvQG4mfX4lrMtHz8H/GXW4yoROS8ib4vI5/O9SES+4ux3vqenp7ARh4AbfeNctEZclSsyhOGE67g7wq2BCdeLXYZF8vzxjQH6xqZdTx9ORwzBt8fZy32MTadcP1/CEkGd7uxhJuXOxexeUoyforncc87/oYj8DeAQ8HzW5s1KqTsish14TUQ+VkpdfugNlfom8E2AQ4cOheEYKYjM1fBuL3CG5YQzEulW90fddlCExB6mRUVpCc+7lN2ZISxrLkbCYlVlGU/vcE/uzBAGe7SbFi2rKvjU5vD0fspFMSKoW8CmrMcbgTvzdxKRo8CvAZ9TSk1ltiul7jh/rwCngYNFGFPoMUybvWvr2Nxc4+rnhEXCaTdtDm1ZTYuLcieEY0JOt7q3eXpnM3VFanWfD0GYC7hBUlmt7ivL3C2GGoZKI1Oz6ezOo/vWFK3VvV8Uw0GdA3aJyDYRqQC+ADyQjSciB4E/IO2cklnbV4tIpXO/BXgGyE6uiCW9o1Ocv97vSXgeBknrZv84HXeHPUmXDUMacac9wo3+cW/Sh0MQYb9/Y4De0WlvzhcI/C+Ys5f7GJ2a9WR5wG0KlviUUrMi8gtAO1AKfEspZYrI14DzSqlXgX8GrAL+o7Ngd0Mp9TlgH/AHIjJH2ll+fV72Xyx5rSPJnHKn+Od8wiBpGS71fspFGCIow7QdudP94p9hKDRiJGzKS4UX3MzudAiDJN5u2tRWlPL0jiK2uveJoqRDKaV+APxg3rZfz7p/NM/r3gIeLcYYokS7abGhsZr96+pd/6wwTMjtpsXetXVsaa51/bPCMQFZfGrzatrqitv7KRdBl7TScqfF0ztaqHdZ7oTgXyc3d0/ubCtqq3u/0JUkAsbY1Cxvdve60vspF0GvnNA3OsX5a/0etqoO9prcrYFxzDvDntkj6BfqXrJHud437pmcFfTz5f2bg/SOFr/VvV9oBxUwzlzqYXrWnd5PuQj6L8JTF9Nyp1fpshJwTet+q3vv7BFca6SzGQGOudD7KRdBP1+MhEVZifDCnvD1fsqFdlABw0jYrK4p5/GtHqWHBn4CstnQWM2B9e7LnRCCCci02b1mFdta3Jc7IQT2SNgc3NxIW737cicEWxLP1O58akczDdXuy51eoB1UgJjJ6v1U5krvp4eRAHuo8elZ3uxKF0P16mr4IEcMA2PTvHut39Pin0GWtG4PTvDx7SGPa80Fd02uOznK1d6x0F+cm412UAHinSv9DE/OerjeEuxKAWcu9TA1O+epnh7kUjan7rW699IewY0YTpiZi9k9Pl8CapB7vZ88kju9QDuoAGEkLKrKS3h2l/vpshmCPAEZpk1DdTmHtxa/1X0+ghxBGabFuoYqHnWh1X1eAixpGQmbnW2r2O5Cq/t8BPmyV8O0+OSmRtY2eCN3eoF2UAFhbi6tHz+/u5XqCu/SQ4M6Ic+k5jh1McmRfW2eyZ0QXIc9MZ3iTFeP58U/g1o8d3B8mneuepndmSaoa1B3hyb48NZQaFu750M7qIDw8e0hrOFJz5uLBVXSevdqP0MTM573sklXdw+ePc509TA5M+f5+kJQJa1THWm50/PjI6CVRu5ld4a4OWEutIMKCEbCorREeGmvt+mhQY2gDDMtdz7nodyZIZj2cOTObd7JnRDcSiNGwmJtvcdyJ8GNoAzTZnure63u/UI7qIBgmDaHtzaxurbC088NoqSV7v1k8+wub+VOcK6DCpg9ZlNznLpoc2RvG+Ueyp0QzAl5YjrFG5fcbXWfjyD+oBsan+ZA6+kAACAASURBVOHtK32Ri55AO6hAcKVnlK7kqD/6cQCbmX18e4i7Q5Oery9ARsIJFu9e62dwfMaX6gBBlLTevCd3+mSPgHns1zptZudU5NafQDuoQHAvPdSH6xcy7ilIJ51hplvdH/UhXTaIay6GaVNZVsJzLvd+ykUQIygjYVNXVcaTLra6z0sAIyjDtGmrq+STGxv9HkrR0Q4qABimxSMb6tnQWO35Z2cCqCBNQkbC4vA27+VOCN6ai1Lp4p/P7mqhpsK9Vvf5CJqkNZu5mN0HuROCV919csY/udMLtIPymeTwJO/fHORln/TjTBpxUM65q71jXLJHfdPTgxYxmHeGuT044WN1gGAVzz1/fYCB8Rnf7BG06u4/7OplfDrleTajV2gH5TMnOmyUh8VQ53M/ggrGaXev+KcP608QvNI+hmlRInDE4+zODEErnttuWlSUud/qPh9Bq+5uJCzqKn2SOz1AOyifMUybLc017F7jT3rovTUoXz79YYyEzYH19WxqcrfVfT6CltVoJGwObW2i2eVW9/kIkj0yxVCf3dlCbaX3cicES/JMzSlOdiR5cW8bFWXRnMqj+a1CwsjkDG9d7vW8OkA2QVqDSo5M8uMbA/6mywZoArreN8ZFa8SXbMYMQZqQE3czcqeP9iAY5wrAe9cH6B+bjkzvp1xoB+Ujr3f2MJPy/mr4bDKOMQiy1qmOpCN3+jkBBWdGzlQH8PX4CFBadSa784iPxVCDJAG3mxYVpSWR6f2Ui6I4KBH5rIh0iki3iHw1x/OVIvJd5/l3RGRr1nO/6mzvFJGXizGesGCYFi2rKji42aPeTwsQhDmo3bTY3FTD3rV1vo0hSNXd202Lfev8kzshWBFUu2lxaEsTLT7JnRCcCCp9MbvFMzubWeWT3OkFBTsoESkFfh/4CWA/8EUR2T9vt58DBpRSO4FvAL/jvHY/8AXgAPBZ4P9x3i/yTM2mON3Zw9F9ayj1MT00KNfpjkzO8FZ3n69yJwRnAuodneL89QFf5T0Ijj1u9I2n5U6/5ayAZHletEa42e9ndqc3FCOCOgx0K6WuKKWmge8Ar8zb5xXg28797wNHJD0LvQJ8Ryk1pZS6CnQ77xd5zl7uY3Rq1vcT7l6auc8n3RuXephOeV8MdT5BiRhOOdmdfqcPB6V4rpHwN7szQ1Cquxumjfh0MbuXFMNBbQBuZj2+5WzLuY9SahYYApqX+FoAROQrInJeRM739PQUYdj+0m7a1FaU8vSOFl/HcS9Jwudp2TBtmmsr+PQWf+XOoKy5GKbNxtXV7Fvnn9yZwX9rpLMZ966tY0uzN63u8xGUSiPtpsWnN6+mtc4/udMLiuGgcv2kmP8fzLfPUl6b3qjUN5VSh5RSh1pb/bkGoljMzaWrA7ywp42qcn8Vzfuljvwbw/TsHK9fTPoud0IwIqjRqVne7O7l+P61vsqdEIziuX2jU5y/1u97dA3BqDRys3+cxN1h39UXLyiGg7oFbMp6vBG4k28fESkDGoD+Jb42crx/c5De0alAHGD3Iyj/OHulj5EAyJ0QjDWXM5d6mPa41X0+glA891RHkjnlbWv3fASh0kgmu/NYBKuXz6cYDuocsEtEtolIBemkh1fn7fMq8CXn/k8Br6l0nPwq8AUny28bsAt4twhjCjRGwqKsRAKRHnp/Dcq/s84wLWoqSnlmp79yJ+B7xAJpezTVVnDIZ7kTgiFpGQmLDY3VHFhf7+s4AEoCkGZuJCx2r1nFthZ/5U4vKNhBOWtKvwC0Ax3A95RSpoh8TUQ+5+z2R0CziHQDvwR81XmtCXwPSAD/Dfh5pVSq0DEFmczV8E/taKahutzv4fgeQWXkzud3t/oud4L/pZ+mZ51W93u9bXWfD78lrbGpWc509XLM5+zODH5HUP1j07x7td/35BmvKEoCvVLqB8AP5m379az7k8BP53ntbwO/XYxxhIHu5ChXe8f4O5/Z5vdQHsCvk+6DW4MkR6YCc8JlZzX6MR++c7WPkcnZQKy3gP8TckbuDMrxgc+S56kO25E7g2IPd/H/J1rMuNf7KSDpoeJzCGWYNmUlwosBkDvB/4jSMG2qy0t5dpf/cif4XznBSNg01pTz+Fb/5U7w32EbCZv1DVU8ssF/udMLtIPyGMO0eGxTI2sbqvweCpBdLNafs85IWDy5vZmGGv/lTvC3gWPQ5E7wN2lk5l7vpzWBkDshc3z4Y5CJ6RRvdvVw/ID/2Z1eEYz/eky4OzTBh7eGApGdlcHPYrHdyRGu9IwF0h5zPtjjo9tDWMOTgbKHn6UJ37nSz/BkMLI7M4j4c2xA+mL2yZm5QGQzeoV2UB6SSQ8Nkn7sZ7uNdjOTLhucE87P4rmGaVFaIrzkU++nXPhZPNdIWFSVl/DcruBc9+jnhdxGwqKhupzHtzX58vl+oB2UhximzfbWWna2+dP7KRf3JmQfTjojYfPJjQ2sa/C+1f1i+DEHGQmbJ7c30Vjjfav7fPhVPDeT7frcrlaqK4Ihd4J/F3KnW90nfWt17xfx+aY+MzQ+w9tX+gKUjZTGr6QAa2iSD28OBiZbLYNf0v7lnlG6k/61us+HX2tQH9+TO7U9AN692s/QhH+t7v1COyiPeK3TZnZOBU4/9qvU0Qmn+Gfw7OFP8VwjgHIn+BcxtDtyp1+t7vPhV/FcI2FTWVbCc7uDkd3pFdpBeYRh2rTVVfLJjY1+D+VBfFpzMRI221uCJXeCf8VzjYTFoxsaWN8YLLnTrzUXw7Q5vLWJ1bXBkTszeG2NtNxp8eyuVmoqotv7KRfaQXnA5EyKNy71cGz/Gkp8LoY6n3uj8fCsG5qY4ezlPo4dCEZ1gGz8iCiTw5O8f2OQlwOUrZbBjwjqSs8oXcnRQGXvZfCjeK55Z5g7QwHL7vQI7aA84IddvYxPpwK3/gT+rEGd7kw6cqe2B8CJDie7M4jHB37Iv8GUO8Gf4rntpkVJDHo/5UI7KA8wEhZ1lWU8ub3Z76E8hB9rLu2mRWtdJQc3BUzuxJ/iue2mzdbmGnYFTO4EfMkaaTctHtlQz8bV/rW6z4cfxXMN0+bxrU00BVDudBvtoFwmNac42ZHkxb1tVJQFz9xer7lMzqRb3QdR7gTvI6jhyRnOXu4NbHUArytrJIcnef/mYCCja/C+eO613jE67ZFARtdeELwZM2K8d32A/rHpwOrHXq+5vHU5LXcGLXtvPl7Z43RnDzMpFcj1J/C+0sjJjiRKEdzzxeNafPcv7g+mPdxGOyiXaTctKkpLAtH7KRdeRwztF2xWVZbx1I7gyZ3gffHcdtOiZVUlj20KRjHU+YjH1efaTYvNTTXsWeN/q/tceF08t9202L+unk1NwZM7vUA7KBdRSmEkLJ7Z2cyqymCmh3q55pKWO21e3NtGZVlwqgNk42Xx3KnZFKcvJjm2v833Vvf58LI/1sjkDG9d7uXlAGZ3ZvAyaaRnZIr3bgwENpr0Au2gXOSiNcLN/olg68ceSjg/vjFA39h0oOUKLyWtty73MTadCux6C3hbqzEjdwb9fPEqfjrVYaflzgAfH26jHZSLtJsWEvD0UC9/pxqmRXmp8MKe4BT/nI+XE7Jh2tRWlPL0zmDKneCtwzYSNs21FXxqczDlTvC2eK6RsNnUVM2+dcGUO71AOygXMUybT29eTWtdpd9Dycv9YrHufo5SinbT5ukdLdRVBaP3Uy68Kp6bcno/vRBguRO8q+4+NZvi9YtJju5bE1i5E7wrnjs6NcsPu3o5vj+Y2Z1eUZCDEpEmETkhIl3O34d++ojIYyJyVkRMEflIRH4m67k/FpGrIvKBc3uskPEEiZv94yTuDgdeP/ZqzaXTHuFG/3ggL1bOxqsciQ9uDtA7OhVouTMbt3/AnL3cx+jULC8/Emx7eLUG9UZnD9OpePV+ykWhEdRXgVNKqV3AKefxfMaBv6WUOgB8FvhdEcm+QvMfKaUec24fFDiewBDE3k+58ErCMUw7LXfuD2Y2Ywav0u4N06a8VHgxYMVQ5+PVj3cjYVNTUcrTO4JdDNWr0k9GwqKptoJDW+PT+ykXhTqoV4BvO/e/DXx+/g5KqUtKqS7n/h0gCQR3EaJItJsWe9bUsbWl1u+hLIhXEYORsDi4qZG2umC0us+LB5JWWu60eGpHC/UBljvBm0ojmVb3L+wJTqv7fHhRPHd6do7XLiY5ui+42Z1eUaiDWqOUugvg/F3w56CIHAYqgMtZm3/bkf6+ISJ5F2tE5Csicl5Ezvf09BQ4bHfpH5vm3LX+wMt74E2a+a2BcS7cHg52dpaDF8Vzu5KjXOsbD4V840WlkfdvDtIzMhV4tQG8iaDevtLHyORsKOzhNos6KBE5KSIXctxeWc4Hicg64N8Bf1spNeds/lVgL/A40AT8Sr7XK6W+qZQ6pJQ61Noa7ADsVIfNXEjSQ72IoDJyZ9DXn8AbexhmuhdWEIuhzscLydNIWJSVBF/uBG/WoIyERXV5KZ/ZFWy50wsWvXpUKXU033MiYovIOqXUXccBJfPsVw/8BfB/KKXeznrvu87dKRH5N8AvL2v0AcVI2KxvqOKRDfV+D2XJuDoBmTa72laxLeByJ3gjaRkJm4ObG1lTH3C5E/cddqa1+1M7mmmoDrbcCbi+KJeRO5/fHXy50wsKlfheBb7k3P8S8OfzdxCRCuA/A/9WKfUf5z23zvkrpNevLhQ4Ht8Zn57lzKWewBb/nM/9MbozBQ2MTfNuSOROcF/SujM4wUe3hkIRXYP7EvDlnlGu9o6FQu4E94vnfnhrEHt4KvDZjF5RqIP6OnBMRLqAY85jROSQiPyhs89fA54DfjZHOvmfiMjHwMdAC/BbBY7Hd85c6mVqNjzpoW5LOKcuJkkFtPdTLty2x73sztA5bHdov9fqPiTHh8tZr0bCprREeGlPOI4PtymoQJxSqg84kmP7eeDLzv1/D/z7PK9/qZDPDyJGwqKhupzHt4UjPdTtCcgwLdbWV/GJjQ0ufUJxcd0eCYsdrbXsaA1g76cFcG1CNi0+uamRtQ3BlzvB/eK5hmnx5PYmGmpCIHd6gK4kUURmU3Oc6khyZG8b5aXhMK2bay4T0ynOdPVwPMDFP+fjpqQ1ND7D21f6Q5EsksHN6u53hyb48NZQaNQGcLd4bndylMs9Y6E6PtwmHLNoSHj3aj9DEzOhSKfO4Oaay5muHiZn5kIj7wGuFs89ddFOy51hOj6cv24cH/ezO0PkoJy/bkRQRiKd3Rnk2p1eox1UETESNpVlJTy3OzzpoW6uuRimTX1VGU9sD4fcCe4WzzVMmzX1lXxiQzjkTnB3zcUwbba31rKzLTzFUN22xyc2NrC+sbr4bx5StIMqEul0WYtnd7VSUxHM3k+5cOuEm03NceqizZF9a0Ijd4J7xXMnZ1K8cSm4re7z4VbEkJY7+8IVXeNe8Vx7eJIPbg6GSu70gvDMHAHnwu1h7gxNhiY76z7unHDnrg0wOD4TuhPOLUnrh129TMykQre+4FZ199c7k8zOqRCeL2mK/QMmTBeze4l2UEXCSFiUBLz3Uy7ciqDaTYuKshKe2x3sqh/zcdMedVVlPLEtuL2fcuFWjkS7adFWV8ljGxsX3zlAuJXr025abGupZWdbuLI73UY7qCJhmDaPb22iqbbC76EsCzfON6XSV8M/t6uF2oC2us+HGxPybGqOkx02L+1to6IsXKecG2uUYZU7wZ2s16GJGc5e7uP4/vBku3pFuM6WgHKtd4xOeyRU2VkZSlxYczHvDHN7cCJ06wvgTpr5e9cHGBifCaV848aay4+6exmfToXyfHEj6/V0yOVON9EOqghk0kPDtt4C7pxwRsKmRODIvuAX/5yPGxGUkbBDKXdClqRVRIMYpk1dZRlPbQ+X3AnuRJSGadOyqpKDm4Lb6t4vtIMqAoZps39dPZuaavweyrJxY83FMC0ObWmieVVwW90vRrHsken99JmdLawKmdwJxa+ckJpTnOxIt7oPm9wJxf8BMzmT4nRnMpRypxeE7wgJGD0jU7x3YyC04XmxJ6DrfWNctEbCa48ihwwdd0e4NTARyugaiv8D5r3rA/SNTYfq4txsii0Bn73cx9h0KrTni9toB1UgpzpsVEh6P+WkyKVbwtLqPh/FlnCMhOW0ug/nBFTstHvDtKgoLeH5EMqdUPwIykhYrKos4+kd4ZM7vUA7qAJpNy02NVWzb114robPpthLDIZps3dtHZubwyd3ggsTkGlzaMtqWkIqdxYzglJKYSRsnt7ZTF3AW90vRjHskcpqdV9Zpns/5UI7qAIYnZrlR93pq+HDmh5azMoJvaNTnLveH8rsrAzFTCO+2T9O4u5waKNJKK4EfNEa4Ub/eLjtUcRfMO/fGKB3dDrU54vbaAdVAG909jCdCk/vp1zcd6uFn3EZuTOs6wtQ3KxGI2S9n3JSRAnYMG1H7gxfdmeGYkqeRsKmvFR4YU845U4v0A6qANpNi6baCg5tDU8x1PlkJuS5IvwiNEybDY3V7F8Xnlb38ynmGpRhWuxdW8eW5uC3us9HUe2RsPjU5tW01YWj91MuiiV5ZrI7n97RQn3I5U430Q5qhUzPzvH6xSRH97VRGuL00GJJWqNTs7zZ3Ruq3k+5KNYE1Dc6xblr/aGOrqF4EvDN/nHMO8Pht4fzt1B/fcke5XrfeLijaw8oyEGJSJOInBCRLudvzivNRCSV1e791azt20TkHef13xWR0NQJevtKHyNTs6HW06F4DdjOXOphejZkvZ9yUpzKCacuJplThH59oViS1v1W9yG3R5GK5xpm+uL+YyGr3ek1hUZQXwVOKaV2Aaecx7mYUEo95tw+l7X9d4BvOK8fAH6uwPF4hpGwqC4v5TO7wtP7KRfF+kVomBara8p5fGu4r4YvVgSVkTsPrA+v3AlFtEfCYveaVWxrCa/cCcXL8jQSNgc3N9JWH1650wsKdVCvAN927n8b+PxSXyjpnyIvAd9fyev9ZG5OYZg2z+9upao85OmhRZiAZlJznLqY5Mi+NZSFqPdTLoohTo5Pz/JmV7oYapjlTijOhDwwNs27V/sjEF0XZ03uzuAEH98eCmVtRq8pdDZZo5S6C+D8zZeeUyUi50XkbRHJOKFmYFApNes8vgVsKHA8nvDhrUGSI1O8/Ej4w3MpgqT19pU+RiZnQ7++AMVZczlzqYep2blIrC8Uo3LCyQ7bkTvDbw+KUDw3I+9F4Xxxm0WLg4nISSCXq/+1ZXzOZqXUHRHZDrwmIh8Dwzn2y/tfF5GvAF8B2Lx58zI+uvgYCZvSEuGlPeE/wIpR2ccwbarLS0NZDHU+xVhzMUybxppyDoc4uzNDMSIoI2GzrqGKR0PU6j4f9+LhQs6XhM3OtlVsb9W9nxZjUQellDqa7zkRsUVknVLqroisA5J53uOO8/eKiJwGDgL/CWgUkTInitoI3FlgHN8Evglw6NChYvdPWxaGafHk9iYaasKfHlqof5pzroZ/bndL+OVOCl9zycidRyMgd2azUntMTKd4s6uHnzm0KfRyJxTusAfHp3nnaj//83PbizamKFPoGfQq8CXn/peAP5+/g4isFpFK534L8AyQUGnN4HXgpxZ6fdDoTo5yuWcsMvpxoZLWx7eHsIYnI7G+AIVPQOeu9jM0MRMNOYvCi+ee6ephcmYu9Nl7GQq9LONUR5LUnIrM/OE2hTqorwPHRKQLOOY8RkQOicgfOvvsA86LyIekHdLXlVIJ57lfAX5JRLpJr0n9UYHjcZ1M76ewtXbPR6GVE9pNi9ISCWXvp1wUuubSblpUlZfw3K7wy51QeFJAu2nRUF3O4W3hlzuh8PPFSFisrY+G3OkFBTWoUUr1AUdybD8PfNm5/xbwaJ7XXwEOFzIGrzFMm09sbGB9Y7XfQykKhU5ARsLmiW1NNNaE5hK2BSkkgsoUQ312VyvVFeGXO6Ewe8ym5jjVkeTI3jbKIyJ3lhQgAU9Mp1vd//SnN+neT0skGkeNR9jDk3xwczBS2TeFTECXe0bpTo5GzB4rl3Au3B7m7tBkpOSbQiStd69FS+6Eworn/rC7l8mZuUgdH26jHdQyyBT/jNYBtnJJyzDT9jgWIXsUUjy33bTSre73RkPuhMIkLcO0qQxpq/u8FFB5pd20qKsq44nt0ZA7vUA7qGVgmBbbWmrZ2Rad9NBCIigjYfHIhno2RETuhMKy+IyExeFtTayujYbcCSuXgJVSGKbFs7taqKkIX6v7fKzUHmm5046U3OkF2lJLZGhihrOX+zgegeoA2az0uo7k8CTv3xjk5Yhk72VYqYRztXeMS/ZoZLIZM6zUYZt3hrkzNBmZ7L0MKz33z18fYGB8JnL2cBvtoJbI6c4ks3MqUno6ZK25LHNKPtERjeKf87nXfmSZ/UdOONmdUTs+Vlo814ig3An3f9DNLdNjG6ZNRVl4W937hXZQS8QwbVpWVXJwU7iLoc5npZJFu2mzpbmG3WuiI3fCyi9cbjdtDqyvZ+PqcLa6z8dKI6h20+bQ1iaaQ9rqPh8rsUem99OzO1uorYyO3OkF2kEtgcmZFKc7kxzbvyZy6aErOeGGJ2c4e7k3cnInsKLiucmRSX58YyBy8h6srHjutd4xOu2RSGV3ZljJmm3i7jC3ByciGF27j3ZQS+Ds5T7GplORPMBWsuZyurOHmVQ0r4ZfSfHcUx1JVFSKoc5jJWn3JyKZ7ZpmJRdyG6adljsjcnG/l2gHtQTaTYtVlWU8vaPZ76EUnZU0LDRMi5ZVFRzcHC25E1ZWPNcwLTY31bB3bZ0rY/KTlRTPNRIW+9bVs6kpWnInrCyCMhI2h7Y00RIxudMLtINahNSc4mSHzQt7Wqksi0Z1gFws9YSbmk1xurOHo/vWhLrVfT6W659GJmf4UXf0sjszLFcC7hmZ4vz1gUjKe9ks1R43+8fpuDscyejaC7SDWoT3bwzQOzoduWy1DMudgN663Mfo1GxkT7jlSlpvXOphOhWdYqjzWW7EcKrDjqzcCcsvntt+r/dTNI8Pt9EOahGMhE15qfDCnmimh8oyYwbDtKmtKOXpHeFudZ+P5VZOMEyb5toKPr0lenInLH/NxUjYbFxdzf514W51n4/lZr0aCZu9a+vY3Bw9udMLtINagEx66NM7WqivCn/vp1wsJ4LK9H56YU9bJHo/5WI5E9D07ByvO72foih3AvezGpew6+jULD/s7uX4/rWRlDtheRFl3+gU56/1Rza69gLtoBbgkj3K9b7xyMoVsLwT7v2bA/SOTml7OJy90sdIhOVOWJ7DfqOzh+mItLrPx3KK557qSKZb3Ud8Pc5NtINaAMPRj49FOD10OSecYablzhcjVh3gQZYuaRmmRU1FKc/sjKbcCctbczESFk21FRyKqNwJy5OAjYTFhsZqDqyPptzpBdpBLUB7wuLg5kba6qv8HoprLPWEy8idT25vjqzcCUuPoO7Lna2RlTth6RHU9Owcr11M936KUqv7+SzVHmNTs5zp6uX4gWhmd3pFdI+kArk9OMGF28ORvNgwm6WecF3JUa71jUdeT19q8dwPbg2SHJmKfHbWUh3221f6GJmcjf7xscQ12zOXHLkz4seH22gHlYcT99JDoyvvwdInoDjInbD04rmGaVNWIry4J8py59IlYCNhUV1eyrO7oit3plni8ZGwaawp5/Gt0ZU7vaAgByUiTSJyQkS6nL8P/TdE5EUR+SDrNikin3ee+2MRuZr13GOFjKeYGAmbnW2r2N4arWKoD7O0NRcjYfPYpkbWNkRX7oSlR5RGIi13NtREV+6EpVUaycidz++OttwJS4ugZu71floTabnTCwq13leBU0qpXcAp5/EDKKVeV0o9ppR6DHgJGAeMrF3+UeZ5pdQHBY6nKAyMTfPO1f7IR0+QVdpnAe4MTvDRraFIZ2dlWMoE1J0c5UrPGC/HwR7O34X89Ue3h7CHo53dmWEpq0nvXu1neHI2FseH2xTqoF4Bvu3c/zbw+UX2/yngL5VS4wV+rqu8djFJai6axVDns5SIIVP8Mw56+lKK52aqAxyNwQ+YpVR3bzctSkuElyKd3ZlmKZVG2k2LqvISnt0VzYv7vaRQB7VGKXUXwPm72BH6BeBP5237bRH5SES+ISJ5qymKyFdE5LyInO/p6Sls1ItgJCzW1lfx6IYGVz8nCCxlzcVIWGxvjVar+3wsRdIyEjaf3NjAuobotLrPx1KquxumxRPbmmisiU6r+3wsVjw33ere5rldrVRXRFvu9IJFHZSInBSRCzluryzng0RkHfAo0J61+VeBvcDjQBPwK/ler5T6plLqkFLqUGure79MJqZTvHGpJ5K9n3KxWAQ1ND7D21f6YxFNZpNvOraGJvnw5mDks9UyLHYZVHdylMs9Y7E5PhaTgD++PYQ1HL1W936xaHtHpdTRfM+JiC0i65RSdx0HlFzgrf4a8J+VUjNZ733XuTslIv8G+OUljts13uzqYXIm2lfDZ7PYCfdap01qTsViPQ4Wt0em1X1c1hcWW4PKyL/H4nZ85HneMG1KSyRyre79olCJ71XgS879LwF/vsC+X2SevOc4NSStM30euFDgeArGSNjUVZXx5Pbo9X7KxWJrLu0XbNrqKvnkxkbvBuUjixXPNUyL7S217Ih8dmeaxdZc2k2LRzc0sL4x+nInLF48t920OLy1idW10Zc7vaBQB/V14JiIdAHHnMeIyCER+cPMTiKyFdgEvDHv9X8iIh8DHwMtwG8VOJ6CmL2XHtpGeUzSQxdac5mciZfcCQtHUEMTM5y93MexGFUHWKjSiD08yQc3B2MTXQMLFs+90jNKV3I0NuqLFywq8S2EUqoPOJJj+3ngy1mPrwEbcuz3UiGfX2zOXx9gYHwmlvpxrhPuh129TMykYrO+AAtLOKc7k8zGJLszw0JrlPdauz+i7QFZ2a4xOj7cJh5hwhJpNy0qykp4fnd80kMXIFFmLQAADblJREFUWgQ3EhZ1lfGRO2HhygmGadNaV8ljMZE7YWGHbSRstjbXsCsG2Z0ZFiqe225aPLKhng0xkTu9QDsoh0x66LM7W6itLCiwDBX50sxnU3Oc7Ejy4t42Ksric5jkk7QmZ1Kc7kzGSu5Mk3vNZXhyhrOXezl+ILq9n3KRL4JKDk/y/s3BWFwr6CXxmXkWIXF3mNuDE7HTj/OdcO9dH6B/bFrbw+Gty72MTafitd5C/gjq9YtJZlLxye7MkM8eJzuSkW517xfaQTm0mzYlAkciXgx1PvlOOCNhU1FawgsRL4Y6n7z2MG3qKssi2+o+H/mquxsJm5ZVlRzcHK9iqPkkYCNhsaW5hj1r6nwYVXTRDsrBMC0ObWmiZVXeYhaRJNcJl+n99MzOZlbFSO5M87Cklcr0foqZ3Am5JeDJmRSnLyY5tr8tuq3u85Ar63Vkcoa3uvs4vj8+2Z1eEa+zLQ83+sa5aI3EMjzPtebScXeEWwMTscxGyjW//PjGAH1j07GTsyC35Hn2cp8jd8bw+HD+ZgdQpzt7mE7NxfJ8cRvtoEiH5xCfq+GzyTUBGQkLETgaM7kT8tjDtBy5Mz7ZnRlyXRdmJCxqK0p5emd8sjvvkdMeNs21FXwqZnKnF2gHRXp9Ye/aOrY01/o9FO/JseZimDaf3rya1rp4yZ3wsKSllMJI2Dy9s5m6CLe6z8f8SiPZcmdlWfyKoc4vnjs1m+L1i+nszrjJnV4QewfVOzrF+ev9sQ3PZd5Pwpv94yTuDsdS7oSHI6hOe4TrfeOxlLPg4TWX928M0DsaT7kTHr5u8OzlPkanZmN7vrhN7B3Uax1J5lT0W7vnY37WmhGj3k+5mC9pGaadljv3xyubcT7Zx0d5qfBiTIuhzl+DMhJ2Wu6MWXanV8TeQRkJiw2N1RxYX+/3UHxhfsRgmBZ71tSxtSWGcicPS1pGwuJTm1fTVhftVvf5yHbYmezOp3a0UB9DuRMeLJ6baXX/wp62yLe694tYO6ixqVnOdPVyLMbpofdPOEX/2DTnrvXHWq7IlrRuD05w4fZwbKNreLC6e1dy1JE7Y2yPrKzXD24N0jMSj1b3fhFrB3XmUg/Ts3OxKv45n2zJ4mSH7cid8bVHBkU6moR4F//MjqDaL8Q32zVDtuLQblqUlUjsLmb3klg7KCNh01hTzuNb45semj0BGabN+oYqHtkQT7kTHlwEN0ybXW2r2BZTuRMeXKM0EjaPbWpkTX085U6YZw/T5qkdzTRUx1Pu9ILYOqiZe72f1lAWk95PuchIOOPTs7zZ1RO74p/zyXz3gfFp3r0Wv1b388kcH3cGJ/j49lDs7ZGJobrsEa72jsU6uvaC2M7M71zpZ3hyNjatu/Pi+KIzl3qZmp2L9foC3JdwTnUk063uY358ZH6rGGam15G2B9y3x7EYXszuJbF1UEbCoqq8hGd3xa86QDaZE+7c9X4aqst5fFuTvwPymWx7rK2v4tENDf4OyGcyDvvc9X52tMan1X0+su3xyU2NrG2Ir9zpBbF0UJneT8/taqW6It7podmLvnFqdZ+P7OK5x2PU2j0f2WuUWs56MM087mqDFxQ0G4nIT4uIKSJzInJogf0+KyKdItItIl/N2r5NRN4RkS4R+a6IVBQynqXy0a0hrOFJfcLBAxOwtseDxWL1egtkNdzQ9iDbGtoeXlDoz+ULwP8AnMm3g4iUAr8P/ASwH/iiiOx3nv4d4BtKqV3AAPBzBY5nSRgJi9IS4UhMr4bPJnPCVZaV8NxufTV8xh71VWUcjrncCfcd9pr6Sj4Rc7kT7ttje2stO2PU6t4vCnJQSqkOpVTnIrsdBrqVUleUUtPAd4BXJP3T/SXg+85+3wY+X8h4loph2hze2sTqWk8CtkCTOeGe3dVKTUXcej/lwLHHkX1rYi93wn2HHb9W97nJSMD6WkFv8OIM3ADczHp8y9nWDAwqpWbnbc+JiHxFRM6LyPmenp4VD2ZkcobK8pLYZyNlqCor5bMH1vJ3PrPV76EEgtU1Fby0t42/9dQWv4cSCNY3VvPsrhb++hPaHgBbWmr4zM4WfubxTX4PJRaImt+7eP4OIieBXD8Xfk0p9efOPqeBX1ZKnc/x+p8GXlZKfdl5/DdJR1VfA84qpXY62zcBP1BKPbrYoA8dOqTOn3/oo5aFUir2C+AajUYTBETkPaXUQ3kMi2o6SqmjBX72LSD758ZG4A7QCzSKSJkTRWW2e4J2ThqNRhNsvJD4zgG7nIy9CuALwKsqHbq9DvyUs9+XgD/3YDwajUajCQGFppn/9yJyC3gK+AsRaXe2rxeRHwA40dEvAO1AB/A9pZTpvMWvAL8kIt2k16T+qJDxaDQajSY6LLoGFUSKsQal0Wg0mmCQbw1K59FqNBqNJpBoB6XRaDSaQKIdlEaj0WgCiXZQGo1Gowkk2kFpNBqNJpCEMotPRHqA6wW+TQvpi4U1abQ9HkTb40G0PR5E2+NBCrXHFqXUQ835QumgioGInM+V1hhXtD0eRNvjQbQ9HkTb40HcsoeW+DQajUYTSLSD0mg0Gk0gibOD+qbfAwgY2h4Pou3xINoeD6Lt8SCu2CO2a1AajUajCTZxjqA0Go1GE2C0g9JoNBpNIImEgxKRz4pIp4h0i8hXczxfKSLfdZ5/R0S2Zj33q872ThF5eanvGWSKbQ8R2SQir4tIh4iYIvL3vfs2hePG8eE8Vyoi74vIf3X/WxQPl86XRhH5vohcdI6Tp7z5NoXjkj3+oXOuXBCRPxWRKm++TeGs1B4i0uzME6Mi8nvzXvNpEfnYec2/lKV2jFVKhfoGlAKXge1ABfAhsH/ePv8r8K+d+18Avuvc3+/sXwlsc96ndCnvGdSbS/ZYB3zK2acOuBRne2S97peA/wD8V7+/p9/2AL4NfNm5XwE0+v1d/bIHsAG4ClQ7+30P+Fm/v6sH9qgFPgP8XeD35r3mXdJ9AwX4S+AnljKeKERQh4FupdQVpdQ08B3glXn7vEL6BAL4PnDE8eCvAN9RSk0ppa4C3c77LeU9g0rR7aGUuquU+jGAUmqEdOPJDR58l2LgxvGBiGwE/jvgDz34DsWk6PYQkXrgOZyGo0qpaaXUoAffpRi4cnwAZUC1iJQBNcAdl79HsVixPZRSY0qpHwKT2TuLyDqgXil1VqW91b8FPr+UwUTBQW0AbmY9vsXDk+e9fVS6w+8Q6Q6++V67lPcMKm7Y4x5OOH8QeKeIY3YTt+zxu8D/BswVf8iu4oY9tgM9wL9xJM8/FJFad4ZfdIpuD6XUbeCfAzeAu8CQUspwZfTFpxB7LPSetxZ5z5xEwUHl0jLn587n22e528OAG/ZIv0hkFfCfgH+glBpe8Qi9pej2EJG/AiSVUu8VOjgfcOP4KAM+BfwrpdRBYAwIy7qtG8fHatJRxjZgPVArIn+joFF6RyH2KOQ9cxIFB3UL2JT1eCMPh9P39nFC7gagf4HXLuU9g4ob9kBEykk7pz9RSv2ZKyN3Bzfs8QzwORG5RloCeUlE/r0bg3cBt86XW0qpTFT9fdIOKwy4YY+jwFWlVI9Sagb4M+BpV0ZffAqxx0LvuXGR98yN34tyRVjUKwOukP61klnUOzBvn5/nwUW97zn3D/DgIucV0ouEi75nUG8u2UNI68a/6/f3C4I95r32BcKVJOGKPYA3gT3O/d8E/pnf39UvewBPACbptSchvV7zi35/V7ftkfX8z/JwksQ54EnuJ0n85JLG47dBimTUnySdWXYZ+DVn29eAzzn3q4D/SHoR811ge9Zrf815XSdZmSW53jMst2Lbg3RmjgI+Aj5wbks6wIJwc+P4yHr+BULkoNyyB/AYcN45Rv4/YLXf39Nne/wT4CJwAfh3QKXf39Mje1wjHU2Nko6c9jvbDzm2uAz8Hk4Vo8VuutSRRqPRaAJJFNagNBqNRhNBtIPSaDQaTSDRDkqj0Wg0gUQ7KI1Go9EEEu2gNBqNRhNItIPSaIqEU835A+dmicjtrMdvufSZB0Ukbz1AEWkVkf/mxmdrNG5T5vcANJqooJTqI309ECLym8CoUuqfu/yx/zvwWwuMqUdE7orIM0qpH7k8Fo2mqOgISqPxABEZdf6+ICJviMj3ROSSiHxdRP66iLzr9MvZ4ezXKiL/SUTOObdncrxnHfAJpdSHzuPnsyK2953nIX3h7F/36KtqNEVDOyiNxns+Cfx94FHgbwK7lVKHSbfu+EVnn/8b+IZS6nHgr5K7rUfm6vwMvwz8vFLqMeBZYMLZft55rNGECi3xaTTec04pdRdARC4DmVYMHwMvOvePAvuzGo/Wi0idSvfjyrCOdJuLDD8C/oWI/AnwZ0qpTIuDJOmq2hpNqNAOSqPxnqms+3NZj+e4f06WAE8ppSbIzwTpumgAKKW+LiJ/QbqW2tsiclQpddHZZ6H30WgCiZb4NJpgYgC/kHkgIo/l2KcD2Jm1zw6l1MdKqd8hLevtdZ7azYNSoEYTCrSD0miCyd8DDonIRyKSAP7u/B2c6KghKxniH4jIBRH5kHTE9JfO9heBv/Bi0BpNMdHVzDWaECMi/xAYUUotdC3UGeAVpdSAdyPTaApHR1AaTbj5Vzy4pvUAItIK/AvtnDRhREdQGo1GowkkOoLSaDQaTSDRDkqj0Wg0gUQ7KI1Go9EEEu2gNBqNRhNItIPSaDQaTSD5/wHIWfc1tWKDvgAAAABJRU5ErkJggg==\n",
"text/plain": [
"