{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Waves and Stability\n", "\n", "Copyright (C) 2010-2020 Luke Olson
\n", "Copyright (C) 2020 Andreas Kloeckner\n", "\n", "
\n", "MIT License\n", "Permission is hereby granted, free of charge, to any person obtaining a copy\n", "of this software and associated documentation files (the \"Software\"), to deal\n", "in the Software without restriction, including without limitation the rights\n", "to use, copy, modify, merge, publish, distribute, sublicense, and/or sell\n", "copies of the Software, and to permit persons to whom the Software is\n", "furnished to do so, subject to the following conditions:\n", "\n", "The above copyright notice and this permission notice shall be included in\n", "all copies or substantial portions of the Software.\n", "\n", "THE SOFTWARE IS PROVIDED \"AS IS\", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR\n", "IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,\n", "FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE\n", "AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER\n", "LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,\n", "OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN\n", "THE SOFTWARE.\n", "
" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Wavelength and Wavenumber\n", "\n", "Let $j$ be the wave number, then\n", "$$\n", "j = \\frac{2\\pi}{\\tilde{\\lambda}}\n", "$$\n", "where $\\tilde{\\lambda}$ is the wavelength.\n", "\n", "A sinusoidal wave is expressed as\n", "$$\n", "A \\cos \\left( \\frac{2 \\pi x}{\\tilde{\\lambda}} \\right).\n", "$$" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false, "jupyter": { "outputs_hidden": false } }, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAD4CAYAAADhNOGaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOy9d3hjZ5n3/3kkWZJtSW6S6/SSKamQIaTMpJCQRgn1R+BlSVhY3l1gWVoS2N2XzpLQl7ZLqKEGyBISICGEJENm0sikJ9MzJRlXSWNbxep6fn8cHVljS7Kkc44k75zPdfka++jo6BlbOvdzt+8tpJSYmJiYmBy/WBq9ABMTExOTxmIaAhMTE5PjHNMQmJiYmBznmIbAxMTE5DjHNAQmJiYmxzm2Ri+gFrxer1yxYkWjl2FiYmKyqHjssccCUkrf3OOL0hCsWLGCHTt2NHoZJiYmJosKIcThYsfN0JCJiYnJcY5pCExMTEyOc0xDYGJiYnKcYxoCExMTk+Mc0xCYmJiYHOfoYgiEED8UQkwIIZ4t8bgQQnxDCLFfCPG0EOKlBY9dJYTYl/u6So/1mJiYmJhUjl4ewY+BS8s8fhmwNvf1HuC/AIQQ3cAngZcDZwCfFEJ06bQmExMTE5MK0MUQSCnvB46WOeUK4CdS4WGgUwgxAFwC3C2lPCqlnATuprxB0cSO2/+Lh3/zFYgGjXoJExMTE/3JZnjh8T/z2Hf/kYmpsO6Xr1dD2RDwYsHPR3LHSh2fhxDiPSjeBMuWLatpEfbdt7Jp5hHY+XlYdR6c+HpY/2po667peiYmJiaGkc3Ci4/Ac7fCzttYFhnDJ+1Mjl0NnWfq+lKLprNYSnkjcCPApk2bapqm86s1X+Y/nn2Ym88ehed+C7f/M/zhQ7DqgpxRuBxazciUiYlJg8hm4cijuZv/7yA8CjYnrL2YuziLDz7Rx2OrNun+svUyBMPA0oKfl+SODQPnzzm+1ahF+DxOfj4zROqCd9Ny4Sdg9CnFIDx3K9z2Xvh9C6x+xaxRcHYYtRQTExMTBSlh+DHlPvTc7yB0BKwOWPtK5V50wiXgcPPQ7c9hdRyhza7/bbtehuB24P1CiJtREsPTUspRIcRdwH8UJIgvBj5u1CJ63U4AApEEAx2tMHia8nXRp2Hk8dk/xL67wGqHNRfl/hCXgtNj1LJMTEyON6QsuOfcBtMvKPec1RfChZ+AdZfNu+f4wwl63Q5DlqOLIRBC/BJlZ+8VQhxBqQRqAZBS/jdwB3A5sB+YAd6Ze+yoEOKzwKO5S31GSlku6awJX+6X6A/nDMHsfwCGTle+XvlZOLIj9we6Ffbccax1XncZ2NuNWqKJgUgpyWQlNuvibJ+RUpLOSloW6fpNUKIQz+aiEFOHwWJTohAX/Ktyb2ntLPlUfziBt5kNgZTyrQs8LoH3lXjsh8AP9VjHQqjWdCKUKH2SELD0ZcrXxZ/Lxet+q3gKu/8Ancvgqt9D14p6LNlER268/wA3P/oi9330/EYvpSZ+8bcX+MY9+3joYxdisYhGL8ekGqSEP/87PPQt5ea/6nw471pY/6qK85IT4TgnDRkTrl40yWI9UD2CiXAZQ1CIxQLLXq58XfIFOHAv3PIu+NHlijHoWW3gak30ZvdYmIOBKOF4CrezpdHLqZqnX5xmPJRgciZJj8uYnaGJAWSzcOe18Oj3YNPfwyv+X02VihPhRD68rTfHlY/pdc2GhqrGYlFyBlf9HtJxxRj49+i8QhMjCUSUv/vwVKzBK6kNdd3BaLLBKzGpmGwW/vAvihE46/3wqq/WZASiiTQzyUx+M6s3x5UhsNssdLfbmQjHa7/IwClw1R9AZuHHr4Lxnfot0MRQAhHlBjqySA2Buu5ALRsZk/qTzcBt74PHfwJbPqKEmkVtIT01imFUsvi4MgQAPpejNo+gkL6NcPUfQVgVYzD6lD6LMzGUvEcwufgMgZQy7xH4I6YhaHoyafjte+CpX8D5/6qEg2o0AjAbxTA9Ap3o9TgqzxGUw3cCvPMOaGmDm16j1AGbNC3ZrORoLqQyPKXBI2wQwWiSRDqrfB8xQ0NNTToJt7wTnr0FLvwknH+dJiMA5KMYvR7TEOiCLh6BSs9qxRg4O+Anr4MXHtHnuia6MxVLkckqDemLMUdQ6MUETI+geUkn4NfvgF23wyX/AVs+rMtl8x6BQUUCx58h8CiGQKlo1YGu5fDOO6HdCz99PRx6QJ/rmuhKsODmuRhzBIVrNj2CJiUVg5vfBnvvhMu/DGcVrZiviYlwAptF0NVm1+2ahRx/hsDlIJnJMh1L6XfRjiVw9R3QMQQ/eyMc2KrftU10QY2rD3W2LsocgerFDHW2mh5BM5KMwi/eAvvvgdf8J5zxD7pefiKUwOd2GNY/ctwZgl6PUoerW3hIxTOgJJC7VypviH1/0ff6JppQd9GnLu1gPBwnmYu3LxaOTMZot1tZ5WsnYJaPNheJMPz8zXBoG7zuO3D61bq/hD+SMCxRDMejIai2qawaXL1Kaal3Ldz8Vthzp/6vYVIT6i76lCWdSAnjocWVMB6ZijHU1YrP5TDLR5uJ+DT89A3wwsPwhu/BaW8z5GUmQnHDSkfhODQEhXpDhtDeA++4HfpOgl+9HXbeZszrmFRFIJLAahGcOKgIeR1ZZOGh4akYg52teN0OAhEdc1wmtRObVIpERh6HN/8ITn6TYS8VMD0CfZn1CAzcEbZ1wzt+p4jY/ead8Mwtxr2WSUUEI0m62+0s6WoDFl/CeGQqxlBnKz3tdhLpLNFkptFLOr6JBpWy8fFn4S0/g41XGPZS6UyWYDSJzyB5CTgODYHLYcPZYjHOI1BxdsDb/weWnQm//Qd48hfGvp5JWQKRBF6Xg4EO5cO0mEpIZ5JpJmdSikeQKx80w0MNJDIBN70aAvvgyl8qqqEGEowmkdK4ZjI4Dg2BEIJet9OYHMFcHG74P7+BFVvgd++Fx24y/jVNihKIJPG67DhbrHhdjkVVOaR6L0u6WulxKeWDwahpCBpCaFRREzh6EN72K1h7keEvqaolmzkCnfG5dWwqWwh7u/KGWXMh/P4DsKN2xe1doyG27wvouLjjB9UjABjqamVkevEYAjWfMVTgEfjDZuVQtUzNJPnNjhdrz6+Ex+DHl0NoRPH2V52v5/JK4o/kuopNQ6AvvW6dZCYqpaUVrvwFrL0Y7rgGRp+u6TLX37mbj/22tuce7wQjSXrald30UKdzUXkEahirMDRkegTV87snhrnmlqd53h+p/slSKl59eAze/ltYcY7+CyyB6hE0fWhICHGpEGKPEGK/EOJjRR7/mhDiydzXXiHEVMFjmYLHbtdjPQvR63YwUe/yQZsDXv9daO1W3lCZ6hvado6GmNCzK/o4IZpIE0tl8tOdhjpbGZ6KLZrf48hUDKtF0Odx5kNDAdMjqBp18/fcSKj6Jz/xM3j+HmWs7bKX67yy8kwYLDgHOhgCIYQV+DZwGbAReKsQYmPhOVLKD0kpT5NSngZ8E/htwcMx9TEp5Wu1rqcSfG4HoXiaeKrOlRdt3fDqr8H4M7Dtq1U9NRBJ4A8nSKazhBNpgxb4vxO1h0DdTQ92tpJIZxeNrv/wZIx+jxOrRdBitdDZ1mJ2F9eA+jvbNRqu7onTw3DXv8Lyc+Bl7zZgZeXxhxN0tLbgsFkNew09PIIzgP1SygNSyiRwM1CuluqtwC91eN2aUaf81C1PUMiGV8NJb4L7vwhjz1T8tF2js7sYs2KkOtQ5BOpueqhTmVe9WEpIR6biDHXNztjuabeboaEaUN8HhZ+lBZESfv8vigd/xbeUAVV1ZiJsbDMZ6GMIhoAXC34+kjs2DyHEcmAlcG/BYacQYocQ4mEhxOtKvYgQ4j2583b4/X5NC843lTVqV3XZF5U5pVWEiArfvA0xYIsYdSfoK0gWw+KZSzA8FWNJ56wh8LocZmioBtTPTVWG4MlfwP674aJPQfcqQ9a1EP6wsc1kUP9k8ZXALVLKwpjMcinlJuBtwNeFEEUHAUspb5RSbpJSbvL5fJoW4atkiL2RtPcoI+vGnobtX6/oKTtHQnlJ84CpPlkVwRIewWLoJUhnsoyF4gzONQSmR1A1gUgCIZSYe0WhtdAI/OnjsOxsOOM9xi+wBMqs4uY3BMPA0oKfl+SOFeNK5oSFpJTDuX8PAFuBl+iwprL0NtojANj4WjjpjfDXG2D8uQVP3zUa5pShDsDUo68W9ffV06783TtaW2i3WxeFIRgPJ8hk5TGhIa/LboYHq0RKSTCSzH+GFvQKpITffxAyyYaFhJRlSPzhRF4s0yj0+N89CqwVQqwUQthRbvbzqn+EEOuBLuChgmNdQghH7nsvcA5g+BDgHpcDiwB/o4XHLvuS0oH8u38qGyJKpDM8749w1movFmEagmoJRhJ4nDbsNuXtLoRgcJHIUatrLPQIelxKscNiU1BtJKFYmmQmy5a1SjRhQUPw1M2w7y646JPKAKoGEYqnSaSzhg2kUdFsCKSUaeD9wF3ALuDXUsrnhBCfEUIUVgFdCdwsj63Z2wDsEEI8BdwHXC+lNNwQWC2CHledewmK0d4Dr/6qMvP4gdIhon3jEdJZyUlDHrrbHaYhqJJAJJkvHVVZLE1lIwVzCFTMXoLqUb3/tX0u+j3O8pVDoVH403Ww7Cw44//WaYXF8Rs8olLFpsdFpJR3AHfMOfaJOT9/qsjzHgRO1mMN1aLryEotbLwCTnw9bL0B1r0K+jbOO2VnbveyYcCD12VvjnUvIvwFXcUqQ52tPPXiVIlnNA+zzWSzoQFvQS/BQEdr0eeZHIv6mfG6HGwYcLOzVC+BlPCHDyojJ6/4dsNCQioTBo+oVDkuO4tBxyH2enD5l8HpyYWI5vcI7BoN4WyxsKKnXZHHMJPFVRGMJPI3T5XBzlYmZ1LMJJu7J2N4KkZ3u502++yerUcVnjM9gorJV465HWwY8PC8P0IiXaSP6Olfwd4/wYWfaGhISEU1YEZ7BMetIWgajwCUecev+gqMPgkP/ue8h3eNhljf78FqEeZgkhpQBOeO/SAt6VocvQTDk7FjwkIwuzs03weVU9hUuHHQQzor2Tc+R2oiPAZ3XgtLz4SX/2MDVjmf2aH1zZ8sXpT0epRYezbbJDIDJ74eNr4Otl4PE7vyh6WU7BwJsWFAGahiDiapjmRamU+tVgypqMnXZh9QMzIVOyYsBBQokJqeYaWog4k6W1vyn6WdhQljtUooHxIyrou3GibCCew2C55WXaL4JTluDYHP5SCdlUzONNGH6fIvK9LVBSGikek4oXiajQNuQIkPJ0yZiYo5mrtZet3HhoZmu4ubd2SllJLhqRhDnW3HHG932GhtsZoeQRUEworooMUiWNHTjrPFcmzl0DO/gb13wiv+H3jXNG6hc/DnegiEMGZovcpxawjUutymyRMAuHyKMRh5Ah78BgC7ckmtjbkRi2oznHkTqIy5PQQqvW4HVotgeGqmEcuqiKmZFDPJzDyPABTDZlaPVU7h8HerRbC+3zNrCMLjiirwkjPgzH9q4CrnMxGOG95VDMexIfAZOcReCye+Hja8FrZ+ASZ259+s6/pzoSE1PmwmjCtiNkl4rEdgs1ro9zib2iMYLhhIM5eedocZGqqCwJzKsQ0DHnaNhpHZLPzhQ5CKweu+0zQhIZWJkPFdxXAcG4Jeo4fY14oQSuLY7oLb3svukUmW97ThcigxwllD0GTrblJUgzk3WQxKL0EzN5UVziGYi7eZih0WAYHwsYZg44Cb6ViKyb/dDHv+CK/4d/CubeAKi+M3eGi9ynFrCHz1GGJfK65eeNWXYfgxTn7hp2zIeQNgGoJqCaqhoWKGIDeXoFkp1kym4nXZTY+gQqSUuabCWa9ww4AHL9O47v04LHkZnPW+Bq6wOIl0hqmZVF4t2UiOW0PQZrfhctiad1d14htIr3s170z+knM6g/nD3e12RR6jWdfdZAQiCZwtFtrt813+oc5WxkJx0pnmlGoYnozhbLHQ3W6f95jX5eBoNNk8VW9NTCiuyEsUNmWt73fzuZYfYknNwBXNFxKCWW/W9AgMpu4jK6tBCJ57ySeYwcFrDn4Oskrzi9UiTJmJKlBGVBavuhjsbCWTlYw36XtgeCrGYGdr0bX3uOxkspKpWPWT7o438rX4BTdU177buNT6KLd3Xw2+Exq0svKoUxTNHIHBeOs5xL4Gnp5y8snU1XQefQoe+lb+uCIzYYYFKsEfSczTGVIZavKmspGp+c1kKmaIsHLmTqgjMgF3XMMBxwa+FbukgSsrT76r2AwNGUtvkxuCXaMhttq3INe/Cu79PPj3AsrOxrwBVEYgksRbJLQCBXMJmjRhPFyJIWji92+zcIwhkBL++GFIRnngpM9wcDJBtEl7cuoxq1jluDYEvkYMsa+CXaMhNgx0IF71NbC3wW3vhWxGGUxiGoKKCBYRnFNR6/ObMWEcT2UIRJJlDEFOeM5MGC9IIC84Z4fnboVdv4cLPk7/6lOREnaPVTnDuE5MhJVBOj2u4hsZPTmuDUGv20k0mWnKHUEmK9k9Glba4d19yuyCI4/CQ9/OK5DWW2ZieCrGY4cn6/qaWshmJcFocl5XsUqb3UZ3u70pDcFImdJRWJweQSieYuueibq/rj8nL9Elp+GOj8LQ6XDWP7Mh162/s5rRlXXEH07Q3WanxWr8bfo4NwRN2ksAHA5GiaUybMzponDym2Dd5bD1epY5IiTSWSJ1NmBfu3svb/3ew00bU5/LdCxFJivndRUXMtjpbMrQkNroNlSkmQyUKWtWi1hUMwk+ddtzXP2jR/OyH/UiLy+x9T8gHlKqhKw2hjpb8Tht1c0wriP+OnUVw3FuCBo+xL4M6uAMVVoCIeCVn4V0nLOHfwTUv7t4eDJGMp3l63/ZW9fXrZV8bLjMh2mos7UpDZsqfVEqNGSxCHra7YtmiP2u0RC3PqlMsK337zsQSXBqawAe/wlseif0rgeUSXVKh3GzGoL6NJOBToZACHGpEGKPEGK/EOJjRR6/WgjhF0I8mft6d8FjVwkh9uW+rtJjPZWianw3bIh9GXaOTmO1CNb0umYPetfAS9/BysO/YakYr3ueYDwcRwi45bEj7BtvzrhqIfmu4hLJYlBCL8NTsaZTcx2ejGER0N9RumLE63IsGo/gS3ftyX9f7ybOQCTBu9M/B5sTzr3mmMc2DHjYPRom04T9GMrQeuMrhkAHQyCEsALfBi4DNgJvFULMH7MFv5JSnpb7+n7uud3AJ4GXA2cAnxRCdGldU6WoDSb+Juwu3jUaZrWvHWfLnEaX864DYeUjtt/UPaQ1Ph3ndacN0Wa3HfPBblYq9Qhmkhmmm6wef3gqTp/HWTY+3OOyL4ohRX87eJR7d0/wd2cuB2Bsur7vW29oJy+fuR/Ofr/StV/AxgEPsVSGw8FoXde0ENmsJBBJGD6QRkUPj+AMYL+U8oCUMgncDFxR4XMvAe6WUh6VUk4CdwOX6rCmiuhqs2OziKZsKlMqhjzzH/AMEDv9PbzO+iCZkafqtp5wPEU0mWHDgJv3nLuKP+8c5/EXmjtxPKs8WtojUAXdmm0uwfDUTMlEscpiGFIkpeSGP+2mz+PgmkvWIQSM1bFST0rJ38d/woytE856/7zH1c9Y2RnGDWAqliKVkYaPqFTRwxAMAS8W/Hwkd2wubxRCPC2EuEUIsbTK5yKEeI8QYocQYoff79dh2UqctRnFuyajSUan47OJ4jk4z/8wU7Kdk3bNn2ZmFOO5D2+fx8m7Nq/E67Jzw527my6kUkgwklSqRdrKh4ag+UpIR6biJfMDKj0uO8Focw8p+suuCR47PMm/XHgCbmcLXpeD8en6GYLo7ns4x/IMT618tzIOdg5r+1xYLaLp8gT1GlGpUq9k8e+BFVLKU1B2/TdVewEp5Y1Syk1Syk0+n0+3hTXV7OIcuwqG1RfD2tbFTdY3sHLqQTi0vS5rUt35Po+TdoeNf37FWh45eJSte/UxykYQiCQUbSZL6aEeswNqmscQZLOS0elYyYohFa/LQTyVJZosMnu3CchkJV+6azcrve28edMSAPo8jvp5BFJiu/fTHJFe/OveXvQUZ4uV1b72pishVfMoi8kjGAaWFvy8JHcsj5QyKKVU77bfB06v9LlG04x6QzsXMAQA97ivYNLqhb98SumWNBj1w9ufG+jz1jOWsay7jS/+aU/TCp8FIsmyYSFQRPycLZamKiH1RxKkMnLB0JCqqBpswqo3gFufGGbveISPXrwun+vo9zjz3qXh7PwdTv/TfDX1Jro73CVPa8bKIbWARR2gZTR6GIJHgbVCiJVCCDtwJXB74QlCiIGCH18LqEN57wIuFkJ05ZLEF+eO1Q1fE8pM7BoN43M7ypaOedweft76NqXJbPcfDV+T+uFVq1jsNgsfufgEdo2G+P3TI4a/fi0EKtByF0Iw2NnKyHTzGAI1X7FkAUOQ7y5uQkMQT2X42t17OWVJB5ef3J8/3udx1scjyKTgns8S8qzld9nNJZsKQUkYj07HmWqisbX+SP3kJUAHQyClTAPvR7mB7wJ+LaV8TgjxGSHEa3OnfUAI8ZwQ4ingA8DVueceBT6LYkweBT6TO1Y3fG4nwWiiqaSISyaKC/C67PwmfS54T4B7PpOfcWwUY9NxOlpbjqlies0pg2wY8PCVP+8lmW6e359KMJpY0COA3FyCJvIIyg2kKaSZp9X97OHDDE/FuO7S9ceop/Z7nEzNpIinDA5nPfFTOPo8f1v1frJYSsqMAMWH2TeYiVCCNrs1P5DKaHTJEUgp75BSniClXC2l/Hzu2CeklLfnvv+4lPJEKeWpUsoLpJS7C577QynlmtzXj/RYTzX43A6kpO7djqVIprPsmwjn299L4XM7GI+mka/4dwjsgad+aei6xkPxfFhIxWIRXHvpOl44OsPNj75g6OvXQiCcLHsDUGm2ATWz8hLlwwLNqkAajqf49n372bzGyzlrvMc81pfzKA3t3UnOwNYbYOmZPOk8c8GCgbwhGGkeQ+CP1GdEpcpx3VkMszITzZIneN4fIZWRJSuGVPKJwlWXK9opW78AKeNc7vFQPP8hLuT8E3y8fGU337hnX1NpNkUTaWKpTNHJZHMZ6mwlEEkav0utkOHJGB6nDbezpex56sCaZusu/t79B5icSXHdpevnPaZuJgwNDz3y3xAZg4s+RSCapLvdjrVMwYDP7cDrcjRVCelEqH7yEmAaglmZiSYxBGrSqhJDALmwwEWfgtAwPPo9w9Y1ForTX6SUTQjBdZetJxBJ8oPtBw17/WoJ5mcVLxwaGmyyyqGRqRhDXW0Lnme3WehobWmq7mJ/OMH3tx/kVacMcPKSjnmPqzkmwwzBzFHY/nU44VJYfta8ofWl2DjYXAljxSOoT6IYTENQ4BE0R3fxzpEQdpuFld72sud5C3WSVp4Lqy+EbV+B+LTua0pnsvjDiXmhIZWXLuvi4o193Hj/gaapYPFX0FWsMjugpjneA+XmEMzF67I3VWjom/fuI5HO8pFXFp/61Zd7DxnWS/DA1yERggs/AYA/kqxoM7BhwM2+iXDT5Lr8ofrpDIFpCGaH2DeJ3tCusRDr+tzYFpCezVeMqJ7MRZ+E2CQ88A3d1xSIJMnK8qVs1166jplkmu9sfV73168F1SB5yyiPquQH1OSE3hqNYggq2w32uBxNkyx+ITjDLx55gbe8bCmrfK6i53icNpwtFmM8gulheOS7cOqV0HcioHw+KqnF3zjgIZWRPO+P6L+uKoklM4QTadMQ1BOHzUpHa0tTKJBKKdk1Gl4wLASzBiy/Gxw4FU56Izz8HQiP6bquuT0ExVjT6+ZNpy/hpw8d5shk42+oecG5MmWDKv0dTixC0fdpNKF4inA8vWAzmYqviYYUfeXuPdisgn+5cG3Jc4QQ9BtVQvrX60Fm4fyPA8rnyV9BCTHMhmKbITxUbMay0Rz3hgByTWVN4BFMhBMcjSYXrBgC6G6zIwTHio5d8G+QScJfv6jrusamj+0hKMUHLzoBBHzt7n26vn4tqDfH7grKR1usFvo8zTGXQF3DQqWjKl6XvSn0hp4bmea2J0d45zkr8+GfUvR5nPqHhvx74YmfwaZ3QZcibhdOpEmmsxXlCFZ627HbLE1hCNQwtVk1VGd8bkdTeARq+dpCPQQANquF7rY58eGe1XD61fD4TRDUL0RTqDNUjsHOVq46azm/feIIexo8/i8YSeBx2nDYrAufjCpH3XhPRk1YV5oj6HE5CMXTDY9tf+muPXS0tvCP561e8Nz+DgM8gns/Cy1tcO5H84fyIyor8AptVgvr+txN0UtQz6H1KqYhQJWZaHxYQH0Trq/AEADFBfPOvRasdrjv87qtazwUp8UqKmrOeu/5a3A1gUx1IJKsKFGsogyoafx7YLhKQ6DudhtZOfTwgSBb9/h57/mr6WgtX/IKSohxIqSjWN6Rx2DX7XD2P0P7bN9CPjxYoV7PhgE3u0bDDRfxq+fQehXTEDArM9HoN8Cu0RBLulor+jCBstOZFx9298GZ74Vn/wdGntRlXWOhOL1uZ1nxNpWudjv/97xV/GXXODsO1bVJ/BgCkURFiWKVwc5WRqdjDddNGp6KYbeW74QtRB1sHmxQwlhKyfV37qbf4+Sqs1dU9Jw+j5NkJsvkjA4zIKSEv3wS2rxw1vuOecifH1pf2e9y44CHo9Fkw3uK/GFlxnIlYU29MA0BigsWT2UJN7ghamcF0hKFlEwUnvMBaO2Cez6ty7rGQ3H6qpDD/fvNK/G5Hdzwp8bJVAciiYpCAipDXa2kMrLhIcLhyRiDnZUZXZi9yTVq3X/eOc6TL07xwYvWzh+iVIJ8L4EeeYLn74VD2+C8a8FxbG4tUKVeT7N0GE+E4/Qs0ASnN6YhoDlGVsaSGQ4FolUZAq/LQSCcnH+zdXbAlo8oH5IDf9W8trHp+IL5gULa7DY+cOFaHj00yX17JjS/fi0Eo8myQ+vnogq8NXpAzchUrOJEMczKFDfCI0hnsnzprj2s9rXzptOXVPy8fC+B1jxBNquo73YuU3JjcwhEElgEZeUlClnfJJpDE+H6TSZTMQ0BhSMrG2cI9oyHybPK94EAACAASURBVMqFO4oL8bodxFKZ4nr0L/sH8CxRvAKNu/LxUKIqQwBw5cuWsrxHkamu9zzYVCbL1Eyq4pAANM+AmmqayWA2NNSIEtLfPj7M/okI11yybsG+l0JU71JzwnjnrTD2NFzw72Cb/7dW5lE4Kt5Zd7S2sKSrteGVQ/4Kex/0xDQEFHgEDUwYVyotUUheZqKYAWtxwgUfh+FcIq1GIok0kUR6wdLReS9vtfCRi9exeyzMbU/WdcREfnfcU0FHqYoq8NZImYlkOstEOFGVR9DusNHaYq17CWk8leFrf9nLqUs7ueTE/oWfUIBaDaMpNJRJwb2fg76T4OQ3Fz3FH66sq7iQZphNUM+h9SqmIQB8LuWX3kiPYOdICJfDlp+hWwnzmsrmcsqV4F0H93y2ZpnqfA9BDQMyXn3yACcOevjq3XtJpOsn6JYfWl/FrsrtbMHjtDW0l2BsOo6UVNxMpqKMrKxvaOinDx1mdDrOdZeuO0ZmuhLsNgtel11baOjxm+DoAUVKwlL8NlZpM1khGwY8HAxEiTVo6lsmKwnWcWi9imkIAE+rDbvN0lBDsGs0xPp+d8VJQqhgMInVpnxQgvvgyZ/XtK5KewiKYbEIrrt0PUcmY/zikfrJVM8mCavbDQ51tTXUIziS62OoJjQEuVxRHUNDoXiKb2/dz7kn+Dh7tXfhJxRB04CaZFRpmlx2Nqy9uORplcpLFLJxwE1WKqHaRhCMJsjK+paOgmkIAKXt3dfAIfbZrGT3WJiNg5WHhaAgt1EuUbj+VbDkDNh6PaSqv8nNnUxWLVvWejlrVQ/func/kTpVZeVDQ1Uki6HxcwnUPobqDYG9rnpD3/3r80zNpLj2knU1X0MZWVnj5+3h/4LIOLzy01DCG5FS5irHqjUEimJqo8JDs81ki9AQCCEuFULsEULsF0J8rMjjHxZC7BRCPC2EuEcIsbzgsYwQ4sncV+3BbI00coj9kckYkUS6qoohUOQThFggpCWEIlMdHlEEuaqkEp2hcqgy1cFoku9vO1DTNaolUIXyaCFDnc6GGgI1LDVQoeCcSj09golQnB9uP8RrTh3kpKH5MtOV0tdR4+zimaPwwH/CulfB0jNKnhZJpEmks1XnCJZ0teJy2BpmCBrRTAY6GAIhhBX4NnAZsBF4qxBi45zTngA2SSlPAW4BCsVwYlLK03Jfr6VBNNIj2DmqSEdXawhsVgtdc2UmirHiHFjzStj+VUWhtArGp+N4nDZa7ZXViBfjtKWdXHZSP9+7/0BdbljBaBKHzUJ7lWse6molHE8TiuvQ6FQDw1Mz+NyOimUxVLwuB0ejybo0w33j3n2kMqVlpiul3+PkaDRZfe5o+1chGYEL/1/Z06rtKlaxWATr+90N6yVohLwE6OMRnAHsl1IekFImgZuBKwpPkFLeJ6VUhVweBiovOq4TikfQmKqhnaNhLALW9S0sNjcXn8tRWcXIRZ+EeEjZTVXBWKi6HoJSfOTidcRSGb51737N11qIQFgZRlJtEjNfQtqghPHIVLzqsBAoyeJMVjIVM9aAHQpEuflvL3LlGUtZscC8jIVQPcyqenemj8AjN8Kpb4XeDWVP1aLguWHAw+6xcEO6zBuhPAr6GIIh4MWCn4/kjpXiXcCdBT87hRA7hBAPCyFeV+pJQoj35M7b4ff7ta24CL1uJ5MzqYaId+0aDbHS217TrruozEQx+k+Gk9+khIcilf/+xkKJmvMDhazpdfHaUwf5n8ePGN5t7I8kqg4JwGxsvlEJ42p7CFTqNbv49qdGSGclH3hFaZnpSumrZVLZtq8AMi8zXY5aKsdUNg56iCTSDWkunAjFcTttFXdp60Vdk8VCiLcDm4AvFRxeLqXcBLwN+LoQoqh8oZTyRinlJinlJp/Pp/vaFizFNJCdI9VJSxTirWYwybnXKgnjh75V8fXHq+wqLsfJSzoJx9OGlzoGI5UNrZ/LUAObyqSUiiGosnQU6tdUdjAQZbDDWXZAUaXkm8oq7SWYHlZkpl/yduhcuuDpWgxBXmpiVP9pfwtR76H1KnoYgmGg8C+zJHfsGIQQFwH/BrxWSpl/x0oph3P/HgC2Ai/RYU1V06gh9tOxFMNTMU2GoOLchu8EOOkN8Oj3laTbAmSyivZOrYniuaz0KnN4DwejulyvFJXOqZ2L1+XAbrU0xBAEIkmS6WxNHoEv7xEYa2APBaOaQ0Iq/dXKTDzwn8rQmc0fquj0QFiRl6hFuG1dnxuLUEK29WaiziMqVfQwBI8Ca4UQK4UQduBK4JjqHyHES4DvohiBiYLjXUIIR+57L3AOsFOHNVVNo4bY71Y7iqssHVXxunIyE5WWZp57jZJse/g7C54aiCTIZGXejdfKih7lJnIwYJzufzYrORpNVtVVrGKxCAY7GzOgRg1HVdNVrJKXojbYIzgU0M8QdLS24LBZKvMIwmPw2I+V3EDnsoqu748k6K5RuK3VbmWlt70hlUP1HlqvotkQSCnTwPuBu4BdwK+llM8JIT4jhFCrgL4EuIDfzCkT3QDsEEI8BdwHXC+lbIghUH/59U4Y1yItUUjVIa3eDbDhtUquIDZV9tRxjaWjc1nS1YZFGOsRTMdSpLOyJo8A1AE19TcE1c4hKKSjtQWrRRgaGpqeSTE5k2JFT5su1xNC0N/hZLySjdeD34RsGrZ8uOLrK/ISte+sGyE1IaVsmEdg0+MiUso7gDvmHPtEwfcXlXjeg8DJeqxBKz2uCmryDWDnaIjudnvNccHC7uLlPRXu1s69RtEfeuS7cP51JU/TIi9RDLvNwpKuNg4GjDME6oCWWjwCUG7E9+/TvxhhIVQvpBZDYLEoQ4OMVCA9lDPeKyp9j1VARSMrI3549Adwyluge1XF1w7UIC9RyIYBD394epTpWKri+SBaiSYzxFKZRZsj+F9BS270Y71zBLtGw2wYcFdd6qjirUU5deAUWHc5PPxtpaS0BHl5iQ793pgrvO35m4oR+MPKzbBW9cbBzlYmwom6V48NT8VwOWx4Wmvbm/UY3FSm/s1W6hQaAiobYv/QNyGTUGTVq6DWPJGK6qHvrqNXMJH7XdRbZwhMQ3AMvjoPsU9nsuwZD9ccFoKC3Ea1u8Fzr4H4NDz6vZKnjIXiWC2iaqmGcqzoaeNwYMawElL1ZthT401gqKsVKXUamlIFw1PKQJraNwT26t8DVXAoMIMQsLRbn9AQzM4uLvleiAbhb9+Hk94I3jUVXzcvL1GjVwizObt6hofyXcWuRZgj+N9EvYfYHwhESaazNVcMwWxVRNUyxEMvVQS7HvwWJCJFTxmbVkrZ9JyUtKKnnXDCuBLSYL5ssLabQH5ATZ0H2Y/U2EOg4nU5DE0WHwpGGexo1bW+vc/jJJlWZkcU5eHvQGoGtny0+OMliCTSxFNZTR5Br9tBd7u9rkNq8l3FpkfQWHrdTvxah2VUgbrb0GIIWqwWutsrbCqby7nXQuwo7Phh0YfHdeoqLmRFroT0kEF5gkAkWdVUqrkM5pvK6u8R1NJDoKIIzxk3d/tgIMpynRLFKmUH1MQmlRzWxiugd31V11XLaLXkCIQQ+WH29WLWIzANQUNRPYJ6zdndORrCbrWw2ufSdB31JlA1S18Gqy6AB78Byfk74LFQXLdEsYqabDwUNGbHHYwqU6mqkfMuRBV8q2cJaTSRZmomVVPpqIrX5SCeyjJjkI7+YR17CFTU91ZRQ/DIdyEZVkKYVaKlmayQjQMe9oyHSWfqky/yhxO0WAWdbfVJThdiGoICet0OUhlZ2lXVmV2jYdb0urDbtP0Zqmoqm8t510LUr9Rpz2F8Oq6LvEQhS7rasFqEYR5BLVOpCnHYrPjcDobrGBoa0VA6qtJjoMyEWjq6UseKISiYXTw3HxMPKWGh9a+G/pOqvq4aJtVqCDYMeEims4ZWuRUyEY7jq0EjSw9MQ1DAbOK1PnkCLdIShVQlMzGX5WfDii1K52Zq9gMZTaQJJ9K6h4bsNgtDna2GVQ5prRYB5YZcz9DQER0MwYJDijRwMPe30j80pHYXz1nz325UChlq8AagUIa89g0BFEpN1CdP4A8n8On8easU0xAUkJeZqEPlkD+cIBBJsGGgesXRuWjWoz/vWoiMwRM/zR+aHUijf7zSyBLSYFRbtQjUf0BN3iPQlCMwTmbisAGlo6BsCnra7ceGhhJhRQtr7SUweFpN1/WHEwgB3TXmiVRW+1y0WEVdDUEjegjANATHoIpp1aO7eJdGaYlCfG4HM8kMM8kaJ4Ct2AJLz4TtX4O0YlDGNIyoXIiVPW0cMqiENBBO1lw6qjLUpRiCeuWKhidj2CxCk7SAkQqkBwNR3UtHVfo8cwbUPPoDJVF83rU1X9MfSdLTbsdm1XZ7s9ssrO2tX8J4ItyYrmIwDcEx1FNvaKdGaYlC8mGBcI27QSGUD15oGJ78BaBtVvFCLO9pJ2JACelMMk0sldElNJRMZ+s2/nFkKkZ/h1NTma5aRmxEd/GhgP6loyr9Hc7Zno3kjCInsfoVsGRTzdfUIzyosmHAU5chNalMlqPRpOkRNAMuh402u7Uu3cW7RkMMdDjp1Oi+wuxIRn9Egyez+hUwdLoyASqTYmxa+R3oXTUEsyEGvRPGqiHUGhoarLMcda1zCAqx2yx0tLYY4hEcCs7ky3715hiP4LEfwUwAziste1IJ+hoCN4FIwvDNofp3Mz2CJsHnrs/Iyl2j+iSKoWCIfa0eAeS8gutg6gV4+leMh+K4HTbaHbrIUR2DmnTUuxojENWnWqTeA2pqnUw2F6/LGL2hQ8Fo5TpWVdLncRCMJknEIkrBwspzYdmZmq6ptau4ENVjN7rDuFEjKlVMQzCHXrfxIyvjqQzP+6O6hIVAx/jw2oth4FS4/8tMTEV1k5+ey9JupYT0sM69BHqVDQ7VcWRlOpNlLBTXlChW6XHp3xk/NZNkyoDSURXV45x5+McQGVeaHDUgpcQf1jc0BMYbArVAxfQImoR6eAT7xiNkslI3j0C3CVVCKB/EyYOs9d9lSFgIlG7oJV2t+bJEvVBj+rUqj6p4Wm24HLa6hIbGQnEyWampmUzFZ4DwnNr4p3czmUpfhxM7Kdof/SYsOxtWbNZ0vWgyQzyV1e2G2tVuZ6DDaXjlkGrAzRxBk9DrdhqeI5iVltBeOgrKjbWrTaf48LrLoe8k3hD5Jf1u4zocl/e06z6XIJgXnNNmCIQQdSshVfsV9AgN9RgQGlLzOHrNIZhLv8fJm61/xT4zBuddo2xGNKCXV1hIPWYTqB6BnuuuBtMQzMHndhCOp4mnjGnVB6ViqM1u1TXuqqm7uBCLheyWa1guhzk3/YD265XAiBLSQCSB22nDYdNe3TLY6axLjkDtYNbDI/C6HEzHUrpKaB8KGlc6CtDfbuGfbLcz0XGKIneikdlmMv1uqBsHPDzvjxp6T5gIx+lqa9GsMlAruryqEOJSIcQeIcR+IcTHijzuEEL8Kvf4I0KIFQWPfTx3fI8Q4hI91qOFepSQ7hoNsa7frauqp6bu4jkEll3M3uwQ547+CLLG6Kys8ColpHqWaAaiSd0Eu9ReAqPR0yNQd5NHdSzLNbJ0FKBz3/+wRAS4r+9qzd4AFOoM6ZMsBsUjyGQl+yeKq/TqgdJM1phEMehgCIQQVuDbwGXARuCtQoiNc057FzAppVwDfA24IffcjSgzjk8ELgW+k7tew5gdYm9MwlhKyU4dK4ZUvG794sPjoRTfSr+ezsjzsPv3ulxzLqr4nJ7hoYCOScLBzlamZlKVz4KukSOTMXra7bTatb/tdcsVFXDQwNJRMmnE9q+y27KaB3iJLpf0G6DgqYZwjewnaGQzGejjEZwB7JdSHpBSJoGbgSvmnHMFcFPu+1uAC4WirHQFcLOUMiGlPAjsz12vYRjtEQxPxQjH07obAp/LUf1MghKMheL8IXsm8Y5V8NcvgQEdtmryUc8S0mCNQ+uLUa8SUmUgjXZvAIzpLj4cjOo6nvIYnvkNTB7iNs/bGdPpveuPJBV5iXb9PILlPe20tlgNTRg3Ul4C9DEEQ8CLBT8fyR0rek5u2P000FPhcwEQQrxHCLFDCLHD7zdupuzsEHtjDIHarq5X6aiK120nqkVmooCxUJwsFhJnfQjGn4E9d+qwwmNZ0tWqqJDq6RHo2Eg0lB9QY6wh0DqQppBZ4Tl9QkP50lEjKoayGdj2Zeg7mRd95x0rM6GBQCRBd5t2eYlCrBbB+gG3YQljteR1sXsEdUFKeaOUcpOUcpPP5zPsdbrb7VgMHGK/cySEELC+X5+KIZX8blBLU1mO8WllRKVr01uhawXc/0XdvQK1hFSvuQSpjDLpSjePoMt4j0BKyfBk83oE6t/GkGay526F4H447xr6O1oZmy4zsrIK9AwPFrJhwMPO0ZAh+lOhWJpkRr+S11rQwxAMA0sLfl6SO1b0HCGEDegAghU+t65YLQKvy7jZxbtGQyzvbtO9YzffXazDTWAspOiiW20tytDwkSdg/180X3cuK3radZOZUBOket0Eet1ObBZhaFPZ1EyKWCqjSzMZQJvdirPFotvISvVvs1LvHEE2C3/9Ivg2wPrX0N/hJJHOMh3TPgckEElolp8uxoYBD+F42pACAjUf2dsgCWrQxxA8CqwVQqwUQthRkr+3zznnduCq3PdvAu6Vimm9HbgyV1W0ElgL/E2HNWmi12Pc7OJdYyFdFEfnou4m9NgNjofis13Fp1wJHcvgrzfo7hWs6GnjcFCfElJ/WN9qEatF0N/hNLRyaFiHOQSFCCF0rR5TVUeXdOlsCHbdBoE9cO5HwWIpPZegBvyRhCGjHmelJvRXIjUiwV0tmg1BLub/fuAuYBfwaynlc0KIzwghXps77QdAjxBiP/Bh4GO55z4H/BrYCfwJeJ+U0rhi3QrxuYyRmYgk0hwOzrChX39DoGdYYDwUp18doG2zw+YPwpFH4cBWzdcuRM8S0qDOHgGoA2oWjyEAHWZTFHDYgIH1ZLNw/5ehZy2c+HqA/BS8oiMrqyQQThoSGlrf70YIY6Qm1HxkI4bWq+iSI5BS3iGlPEFKuVpK+fncsU9IKW/PfR+XUr5ZSrlGSnmGlPJAwXM/n3veOiml/lnJGuh1Ow0JDe3WYVh9KdTYuB65jbHpObOKX/J2cA8q7ryOqJVDeiSMjegoHepsNTQ0pF5br9AQqPOrdfIIgjP6J4r33AHjz+a8AcXA9JcaWVkl0UROhtyAWHu7w8by7jZDSkjVTedizxH8r8PnVhQRM1l9QyF5aQkDQkMtVgudOshMxJIZQvH0sfFKmwM2fwheeBAObde40lnyg+x1yBMEo/rISxQy1NXKWChu2PDy4akYzhZFHkQvvC6HrjkCXcdTSqkUHnSthJPelD+s7oS1egR6Da0vxYYBD7vG9DcE/nACZ4sFtwFKv5ViGoIi9HocZLJS1w5NgJ2jYTpaWxg0SNXT63JorhpSP4zzBOde+nfg6oOt12u6fiF6lpAGIkkcNgsuHT9MQ52tZKU+IYtiqKWjeg4r73HZCUaTZDVuYqZmkkzHdC4d3fdnGH1KKUCwzv6dHDYr3XNHVtaA3nmiuWwc8HA4OEM4rj2pXYjaTNaIofUqpiEowqy+v77hoUcOBjltaadhf3A91CfVaVH9c41VSyuc80E4tA0O6aNB1GK1sLSrlUMB7SWkatmgnr/bwXxTmTGGQM9mMhWvS9nETGmswDmYF5vTyRBICVu/AJ3L4NQr5z3c53FqDg0ZPdzltGWdAOw4NKnrdRstLwGmISiK6qrqmTAemYpxwB9ly1qvbtecix4yE2VHVG56Z84r+IKm1yhkeY8+g+wD0aTuO0E1dq8Kw+nNyFSMJTrmB4D8vGat4SH1b6KbvMTePyllyOdeA9b5obA+j0O7R5DLjRhVffOyFd3YbRa27Qvoet2JBncVg2kIiuJzKTdBPT2C7bk3z2YjDYHLrnnN+dBQsfBVSyts/rDiFRzcpul1VFZ6lV4CrSWkRjQSDXYYN6AmnsoQiCR1rRiC2bCI1vLnQ4EZ/VRHVW+gawWc+taip/TPHWJfA4FwQnd5iUKcLVbOWNHN9v36Khs0uqsYTENQlFmPQD9DcP8+Pz63g3V9+nYUF+J1OYgmM8SStVfgjofiuBy20rH2068G94Dywdah/n95TxvRZEZzpUswmtA1UQzQarfS025n2IDQkFqWqndoyJf3CLT9Pg/lSkf1kPRmzx1KbuDca4t6A6B4oIFIkpSGxHwgkqBLZ3mJuWxZ62XveCQfQtVKPJVhOpYyPYJmxNlixe206eYRZLOSB58PsmWN19CEkE+HXoLxUJy+cvXMLU7FKzj8ABy8v+bXUdGjhDSblQQjxtSPGyVHbUQPAcyGhrSGCA/pVTqazcJ9X4DuVXDKW0qepnqgWjZfyohKY7wBFdWj375fn/BQvpnMNATNiZ4jK3eOhjgaTRoaFoIC5VQNN4Gx6Xjx/EAhL32H0ldw339o9grUWbhaVEhD8RTprMzfBPVksMOYpjI13KS3R9DZ2oLVIrR7BIGoPvmB3X9QhAvPvfaYSqG5qFVqWnbagYjxIZYN/R562u1s36dPeGh2RKWZLG5K9BxiryaXNq8x1hDMCs9p8QgSC88qbnHClg/Diw/Dgftqfi1QdtzKIPvaDYERw0hUhrqUpjK9xcZGpmJYRIlcjAYsFkF3u12TRzAZVUpHNVcMZbNKuXH3ajj5zWVPnZWZ0GIIjPEKC7FYBJvXetm+P6i5RBcaP7RexTQEJfC5nbp5BNv2+Vnf7zZcVEoV26rVI8hm5bE6Q+V46TvAs0Rx+zXcJPUoIfWH9ZeXUBnsbCWWyjA5o2/t+JGpGH0eJy0GxLO1ykzkK4a0GoJdt8PEc3DedWW9ASiQmdDoEdRj5u/mNV4CkQS7x7TrDjV6aL2KaQhKoHgE2g1BLJlhx6FJw70BgJ52bVLUwWiSdFYu7BGA0m285cNw5G/w/D01vZ7KCm+7ptCQ2lVsSI7AoAE1es4hmItWmYnZ0lENhiCbVYQKe9bCyW9a8PSuthbsVkvNHkE0kWYmmamLIdiyVpHB16N6yB+KYxEYEtasBtMQlMDndjCTzGgeVfi3Q0dJZrKG5wcA7DYLHa21y0yU7SEoxkv+DjqWavYKVvS0czhYewlpwMCO0vyAGp1LSIenYrpqDBWi1SM4GJjBImBpt4b17fwdTOyE8z+W1xQqhxCCXg29BEaGB+fS3+FkTa9Ll34CfyRBd7tD1/nltWAaghLMzi7W5hVs2+vHbrXw8pU9eixrQXwamspKdhWXwmZXxMOGd2iaV7AiV0Jaa0grGE1iEdDZZkyOAPT1CDJZydh0XPdEsYrXZdeULD4cjDLYqaF0NJtRvAHvurzCaCX0e5w1h4aM7iqey5a1Xv528CjxlDax5IlQ45vJwDQEJcmPrNTY5LJ9f4BNK7p0GU5eCVqaysbDJXSGynHa/1FkAzRUEC33qoPsa8sTBAzcVXW1tdDaYtW1hNQfTpDKSMNCQz0uB7FU7d7soYDGOcXP3Qr+3XD+dRV5Ayp9Hc6aN15G5omKsWWtl0Q6q1luotFD61VMQ1ACPUoxJ0Jxdo+F6xIWUtESFhifVuKVVbnX1hZFNmDkcdh7V02vq7WE1B/WX15CRQjBYKdTV49AlawwLkdQey+BlJKDWkpHVW/AtwE2Vu4NwKxHUEuIsN4ewctX9tBiFWzTmCdo9NB6FU2GQAjRLYS4WwixL/dvV5FzThNCPCSEeE4I8bQQ4i0Fj/1YCHFQCPFk7us0LevRk3xoSMNcArXp5Ny1xs1YnouWCVVjoThel6P6zsxTc7ONa+w2XtLVis0iapajDkaNrRYZ6mrT1SNQO5WNyhH0aBhiPzWTIhRP1+4RPPs/ENib8waqex/1e5zEUooMerWohsAoeYm5tDtsvHRZF9v21p4nyGZlXXofKkGrR/Ax4B4p5VrgntzPc5kB3iGlPBG4FPi6EKKz4PFrpJSn5b6e1Lge3ehsa6HFKjR5BNv3Behut+fH3NUDn9tBJJGuKXY5FkrUVteuegWjT8Ke6mcL2XKD7LWEhvSWlyhkqNOpq96QUc1kKlo6zA9qKR3NpBVvoPdE2HBF1U9Xy5ZrqRzyhxN0tbUYUo5bii1rvewcDdXsgU/OKFV6i94jAK4Absp9fxPwurknSCn3Sin35b4fASaA+m2Ra0QIoYysrNEjkFKyfX+As1f3YKljRYAWCe3xSrqKS3HKlcrAkRq9Ai0lpEbJS6gMdbYSjCY1JwZVRqZidLS26Do7oRCvBr2hw1pKR5+9BYL7c5VC1d9atHQXN2JnvTnn6T9Qo9zE7IjKxnYVg3ZD0CelHM19Pwb0lTtZCHEGYAeeLzj8+VzI6GtCiMabxgJ87tqH2O8djzARThgqO10MLU1lY6F4dYniQqw2pXFo7GnY/ceqn74iJ0ddbXx4Jml8/fisHLU+XsGwgT0EMBseqckjqLV0VPUG+k6G9a+u+nWhwBDU4BHUo6t4LicPddDR2pJXFq6WZtEZggoMgRDiL0KIZ4t8HeP7SeUTXPJTLIQYAH4KvFNKqUoMfhxYD7wM6AauK/P89wghdgghdvj9+srAlsLndtZcNbQtp0WyuY75AahdZkJVQdQkeXDymxU5ga3XKw1FVbCip42ZGkpI1V2vkaEhVY5ar4TxiAEDaQpR+0lqmUlwKFBj6egzv4ajB2r2BmBW9beWATX16iouxGoRnLOmh237AjUluPMewWIwBFLKi6SUJxX5ug0Yz93g1Rv9RLFrCCE8wB+Bf5NSPlxw7VGpkAB+BJxRZh03Sik3SSk3+Xz1ubn2emoXntu2L8AqX7uhO79izFaMVBcWUN3xmkNDMOsVjD+jiI1VQV6FtEqpCdVwGDWMBAo8Ap3yBMOT+g+kmUtPjd3Fh4PR6lVHM2n46xeh/xRY/6qqX1PFyxmmGwAAIABJREFU2WKls62lNo/AgHkUlbB5jY+xUJzn/ZGqn9sMQ+tVtIaGbgeuyn1/FXDb3BOEEHbgVuAnUspb5jymGhGBkl94VuN6dMXncnB0pnqN9EQ6wyMHFdnpejNbMVKdASs5q7haTn6TIitQpVew0lvbIHvV8zHSI+jzOLEIfUJD07EU4USawU6Ddadc1Yc186Wj1SaKn74ZJg/C+R8HjTLrtQyomUmmiSYz+bBoPVFDv7V0GfvDCVwOG232xg2tV9FqCK4HXimE2AdclPsZIcQmIcT3c+f8f8C5wNVFykR/LoR4BngG8AKf07geXen1OJCy+qTbY4cmiaeyeU2SeuKwWWuSmRjPTybTuDuxWBWvYOI52DVvX1CSoc5cCWmVKqTBqPGNRC1WC/0epy6GYCQ/h0CnEZAl8LkcVYeGJnOlo8t7qlhbJqV4AwOnwbrLqlzlfPo8TsarLNBQtbWM9ApLsbS7jRU9bTUZgmZpJgPQZIqklEHgwiLHdwDvzn3/M+BnJZ7/Ci2vbzSFFTjVxM637Q9gswjOXF0fWYm51NJdrBoCXSoYTnoD3P8l2HqDUkZYQczYZrWwtLutakOgegRG148PdrbqEhpSr2FUD4FKLaEh9XdfVWjoyV/A1GG47IuavQFQPIKdo6GqnqN6Pt4G3VQ3r/Xy28eHSaaz2G2V762bYUSlitlZXAb1pljtXILt+wK8ZFmnYeWBC1FLd/HYdII2uxW3Hmu2WJWGIv8u2HlrxU9b3tNWdY4gGE3idtpwthgr4THU1crItA4ewbTaQ2B8aGg6liKZrjw8p4blllcaGkon4f4vw+BL4YRLalnmPPo6nAQiiarCsYE65InKsXmNj5lkhideqE5uolm6isE0BGXJy0xUsbs+Gk3y7Mh0Q8JCKl539d3F47nSUd1GaW58vSIzsPUGRXagAmopIfVHEnW5AQx1tjI6FSejcRjJ8GQMu82Ct93YNas5k6PRyt8HhwJRLAKWVTqw/smfw/QLcMG/6uINgOIRSFndZ25WebQxN9WzVvdgtYiqx1eaHsEiQb3BVCOE9cD+AFJSV32hufhcjqrLR8dCGprJimGxKF5BYI8iQlYBK73tSglpFWsPGtxVrDLY2Uo6KzUPKzoyFWOww2l4k2EtekOHgjMMdbVWFt5IJ2HbV2BoE6y5qNZlzkPNUVVTOeSvQ8FAOTpaWzh1SQf3V5EnmEmmiSTSDR9RqWIagjLYbRa62lqq+vBv3xfA7bRxylCHgSsrj8/tIFylzMTYdFz3sYlsuEKRG/hrZV6BmqQ8VIXURL0aiWabymqfpAa5gTQG5wdgVjiwOkNQRcXQEz+F6RfhAu2VQoWoN8ZqegkCkfrLS8xl81ofzxyZYrrCSXbN1EwGpiFYEF8Vs4sLZSWqFm7TEfUmUKkBy2YlE2GdPQIo8Ar2KmJkC1BLCanROkMqeg2oGZ6M5RvUjKTafpKqSkfTCcUbWHIGrJ5XK6KJ/MjKKjyCQLj+XcVz2bLWS1bCg89X5hU0UzMZmIZgQXrdlWukHwhEGZ6KNTQ/ANWHBZReCUm/x4A35frXQN9JSl9BpryqpFpCerDCyqFUJsvUTKo+HkF+ZGXt8ykS6QwT4USdPAJVb6iy98DkTIpwPF2ZxtDjP4HQsO7eAEB3m50Wq6jOEDSgq3gupy1VikMqDQ81y9B6FdMQLIDPXXl3sao5Um99oblUuxuc7SEwIF5psSiyA0efh2d+U/ZUtYT0cIWGQE2E1mPea7vDRmdbi6bQkNq9XY9u8za7FWeLpeLNgCr4t2KhHoJUXPEGlp4Jqy7Qusx5WCxC2XxV0UsQiCQaVjqq0mK1cOaqnornGPtzUQbTI1gkqEPsK6lk2bYvwNLu1srL7wxC3WVUehPQtYegGOtfrcgP3P/FBb2CFT1tHKywhHS2bLA+ScLBjlZNHsFwvpnMeEMghMDrclTcDFmx6ujjN0F4VNdKobn0d1Q3stIfrk/l2EJsWevlxaOxijYyE+EENougy4DxqrVgGoIF8LkdJNPZBYdlpDJZHj4QbHhYCGarJyr1ZMamlfM0y0uUQghFfuDoAXj6V2VPXV7FIHvV46lXWGCoq5UjkzOkM9mavo4crU8zmUpPFTITauno0q4yHkEqBtu+CsvPgZXn6rTK+VQjMxFLZhomLzGXauQm/DltpHpK1Jej8SIXTc5sL0GcjtaWkuc9+eIUkUS6IfpCc3HYrHictoo9grFQHCEMjleuu0yRIfjrDXDSG6GluNEpLCFdyENR49/1CA2BMknt7p3jrPm36ofvqFiEQSG4Ivhc9oo9mIOVlI4++n2IjMEbv2+YNwCKzMR9eyaQUi7Y19LoHoJCVnoVkclt+/y8/czlZc+dCCfyaqvNgGkIFiA/xD6cYE2vu+R52/YFsAg4e3XjDQGoTWUVhoamlRGVhpbfCQGv/DT85Ap4+Nuw5SNFT1NDEwcD0QUNwexNoD67wXdtXonX5SCroalslc9VvcRzjfS0O3j6yHRF5x5eqHQ0MqFoCq25CFZu0WmFxenvcDCTzBBOpPE4S2++oD7qs5UihGDzGi93PDtKOpMtWznoDycYqNOGoBJMQ7AAlXYXb9/n55QlnXS0lX/j1guvy5EX41oITQNpqmHV+Uq+4P6vKHOOPYPzTlGTlYeDM7x8VXmtpkAkid1mqZuUx5KuNt53wZq6vJYeeN12gtEk2awsG4JQS0dfd9pQ6Yvd8xlIzcCl1xuw0mNRy5jHp+MLG4Jw83gEAFtO8PKrHS/y9PA0L102b4R7nolwglOXNq7XaC5mjmABVPetXBXDdCzFky9ONbxaqBBfNR6B3l3F5bjk85BNw92fLPpwNSWkgZy8hG6yGP/L8LocZLKS6Vj5Jqej0WT50tHhx+CJn8GZ/wTetQas9Fj6qphUli8YaJLqm3NWexGCskPt05kswWhzJLhVTEOwAG6HDYfNUjbp9tDzQbISNjdBfkDF56q87HUsFNcuP10pXSvgnA8oE61eeHjewzarhWXdbRU1lQUiyYbJCiwGeirsJ1E7uVd6iySKs1m48zpo98G51+q+xmJUM7tY9Xqb5X3Q1W7npMGOsmWkR6NJpARfE8wqVjENwQIIIej1OMqOrNy+30+b3cpLyriC9cbrslckMxFPZZiaSdUnNKSy+UPgGYI7ry0qPbG8p60imYlgEzQSNTOzMhPlQ4RlVUef+TUceRQu+hQ4PTqvsDhqMr2SRs5AJEFng+Ul5rJ5rZcnXlCKR4rRbF3FYBqCivAtUIa3fV+As1b1VKVFbjSVdherIS/DegiKYW+HV34GRp9SNGvmsMJbWQmp0lHaHDvBZqTS98DhYInS0UQY7v4EDJ2u5HTqhLNFGa5UkUfQhJuBLWu9pLOSh58PFn282XSGQKMhEEJ0CyHuFkLsy/1bdEsshMgUTCe7veD4SiHEI0KI/UKIX+XGWjYd5TodXzw6w6HgTEPVRotRaXexbiMqq+WkN8Kys5UkZGzqmIdW9CysQprNSoKRZN1KRxcjlRqCkqWj938ZIuPK0JkaB9LXSr/HWVGOQKnHb67bxunLu2htsbJtX/Hw0ESTdRWDdo/gY8A9Usq1wD25n4sRk1Kelvt6bcHxG4CvSSnXAJPAuzSuxxB87tIewbYmkZWYS767eAH3esxIeYlyCAGX3QCxSUWHqIDCEtJShOIp0lnZdLvBZqKztQWrRSzYXXyomNhc8Hl46Ntw6ttgySYDV1mcvo7KmsoCkQS+JpFyVnHYrJyxspttJeYTNFulE2g3BFcAN+W+vwllAH1F5AbWvwJQB9pX9fx60ut2MDWTIpGeH8/evt/PQIeT1T5XA1ZWGlV7ZaHOUlXut25VQ4UMnAKnXw1/uxEmducPr8zdlMqNrax3D8FixGIRdLfby3oEUkoOBaPzx1Pe9a9gc8JFxau7jKbf46gwNJRsyvfAlrVeDvijRedcT4QTdLS2GD5Vrxq0GoI+KeVo7vsxoK/EeU4hxA4hxMNCCPVm3wNMSSnVjMoRoEwhc+NQS0jnhioyWckD+4NsXuNtuhLGntwM30o8gtYWpRO5IVzw7+BwwZ+ug1xOYLDTmRtkXzphXG95icWKMra0tEeglo4ekyjedzfs/ROcdw24++uwyvn0e5SRlekyIyvjqQyRRLop3wOq1Mz2IuGhiVDzjKhUWdAQCCH+IoR4tsjXFYXnSSWzVyq7t1xKuQl4G/B1IcTqahcqhHhPzpjs8PsrU/jTi1JNZc8MTzMdSzVdfgCUhJu7ApkJpXRUxxGV1dLeoxiDA1th9x+BykpIA3l5iebbDTYTXld5j2B2YH0uUZxOwp8+Bj1r4OX/VI8lFqXX4yQry3u0+aRrExqCE/pc9LodRXWH/JHmGVGpsqAhkFJeJKU8qcjXbcC4EGIAIPfvRIlrDOf+PQBsBV4CBIFOIYS6FV0CDJdZx41Syk1Syk0+X32F3QplJgpRrf05TdQ/UIhvgd0gKKGhvkZrnmz6e+jdqIQjUko4YIW3vWyOIGh6BBXhdTkIRssYgpzSaz5H8Mh/Q3A/XPIFsDXOyFbSS6AaiWYQnJuLKjfx4PPBeZIkE+H44vMIFuB24Krc91cBt809QQjRJYRw5L73AucAO3MexH3Am8o9vxko5RFs2xfgxEFP096MvBXMUhg3YjJZtVhtinTB1GF46JuA0ktwODhTsoQ0EElgETSNjG+z0tNuLys1cihXOrqkqw3C44qe0NpL4ISL67jK+ajFC+NlOvoDeY+guZLFKltO8HI0mmTnaCh/TErZVEPrVbQaguuBVwoh9gEX5X5GCLFJCPH93DkbgB1CiKdQbvzXSyl35h67DviwEGI/Ss7gBxrXYwg97XaEONYjiCbSPP7CZFOGhVQUj6B8onA8lKh/6WgxVp0HG69QZI6nh1npbSeWypRsKgpEEnS327E2iYxvs+J1O4ilMkRLNDcdDERZ0tWmlI7e82lIx+HSL9R5lfPJ6w2VqRzK54ma0COA2UjB/QV5AqXJM9s0Q+tVNBkCKWVQSnmhlHJtLoR0NHd8h5Ty3bnvH5RSniylPDX37w8Knn9ASnmGlHKNlPLNUsrKxxLVEZvVQk+7/Zjd9SMHg6Qyki1rGj9/oBRel71sjHVyJkUynW28R6By8edAZuHu/7+9cw+Pqj7z+OdNwkzuCcmEhEswFLEWFalGLIK0Aq6uWsELeClqW1fdSvtsa/Vpu3atvde1tV23rYItlbrW69pKn1bxUhSxgsWtRS4VUi5KS0ImgWRIyP23f/zOGSaTmcyEuZ2T+X2eh4dk5szM+8vMmfe8v/fyvSu4VREtT5Au0Xq3c0yyMnJUsK+lU5fr7n8L3n4UZt8KlSNO4SWdyqLYkpXBPFGRMz8H40ryObmmJKhcCM5sJgPTWRw3VSX5QXk50NtC3rwc6uucM1YiHF+xl0BX9DET9v5r2nsIolE+GeZ8HrY+zQe73gGil5CmS7Te7QRFiiJcECil2OvvYEpFPjx3BxRXw7w70m1iRGzJyqZhcgT+I7oM00kd/eGcO83H5r2HONqjz8FgJ79xBO5kXNh++2u7/MyaUuGoWuBw7F6Clo7IV4N22O2YiABgzr9B6STGvX4X3lwVVbayxUQEcVEVjAiGOoLWjh4C3X0s6FmnJ4wu/Dp4o2tupJvqUu+wEYETu4rDmTutip7+ATbt0eMmbIfsJFEaMI4gbqos7WKAA21HaTh4xHHdxOHYXwLREsYZ6yoeDk8hXPAtpOkdbineEFX/1YkzZpxI5TCD5/a2dFBMJ2fvvh8mnQUzrkq3ecNSUzb8mAm/A8sww5lVV4EnNye4PWQPr3Ragts4gjixI4KBARV8U52gTzwcvhhjJhrbtESl08JUpi+GE+ZyU++jNB9sGnJ3Z08fnT39ZmsoDuz980gRwR5/J5/N+w3eLr8e95HmeUKxiL015PyosMCTy1lTxrLBGjfRHOjGk5dDaYGzNMGc9c47mKoSL30DisNHe9nQ4MdX7OXkGueE0ZE4NoY4siNoau+isijFEpXHgzWHqGjgCIsOrx5SQmp6COLHk5cTVb+6ff92Pp37HP2nL9MTRh1GTVk+HT39BLoiC+v4A+6ICueeWMVfGwMcbO/SpaMOFFNy2DeAc7HLvZrau9iwy8/cEysd92aGE2v6pFYmc+iJVHMqO2uXcI28QMuetwfd1WzmDI0IrV89dGvo7He/T494yT0/M/OEYlFTGr2XoKtXaxo7fWsIjg2k3NDgd5xovY1xBHFif+DW72ympaOHuQ7fFgJrzIQ3L2p3caNTegii0DrrdgIUMuaFrwTnEIGJCEaKL1I/yc61nNK5iTXl10HxuMwYFoPhegmOTfB0/sXA9PGlVBR52LDLH4wInIZxBHFi76M/8396CobTE8U2VcN0Fze1d1HtpERxGLUTJ/GDviWUNb4BO4IyFiFzhpx3QjmRIfOG+rpRz3+ZPWoCDSekT3BmpNhFDJHGTDhNq3g4cnKEOSf6eK3Br8dLmIjAvdgfuHebApxUXeysksth8EVRV+vu66e1o8fREcH4snyeZgEHC0+EtV+FXj3StyXYSOT8q0EnoOcNhUSFGx9AWndzd+91TKoqz5xhMagZRsTebdNnzz3RR3Ogm0OdvY7rKgbjCOKmyJtHkUf3DMx1cDdxOL6SyNMn7cYWJzuCvNwcJlSU8Ej5cmh7D16/H9BfAiX5eY7u4XASlUVaT6O3fwACjbD+Xg7XLuDVgdOH6hA4iAKPHo8eaWvomB6FOxxB6CgaJ0YxxhGMAFvX99yT3LEtBNb+cIStIfsqy8lbQ6CnYr7YeSKcchls+CEcft/0EIwQexZPa0cPvHQ39PfwxrQvAnq4n5OpKcuPvDUUcNcY8gnlBUyt0k7XceXaGEcwIqqKvXhyczh7SkWmTYkbX7GX9ghjJoLjJRwcEYB2BPtaOlHnf0Pf8MKd+ANdrkgSOgXbaR5peB3+8hjMXs62oz5yc4TaCmc7gurSyJKVzUe6Kc3Pw5vnnqjQ7jsyEYHLWTh9HNfMqqXQ46xmkOGoijJmoilTovUjZIqvkKO9/TTJOJh3O2x/lo8cWuPYQWNOxFfsoYJ2Jr60HMpq4dzb2dvSwaSxBc7rIQmjOoqIvRu6isO5/IyJfHhyueNkbQHc843mAG6el/mpjCMl2EsQ6GZieUHw9qb2LrwO7HAM54QQ/eKauV+A9zexfNdKVuXMAJzXBOVEfIW5/HjM/YzpaoVlL4C3mL0tHYPlKR1KTWk+zYFu+gfUoJHj/oDzu4rDmTGpnF/fOifTZkTE2ZcDhoSJ1l3c2N6dWYnKOLGTmXv9HZCTS+/ilewf8LFs71eh/UCMRxsAxv/pe5yTu51XP3gnTJhpTR3tZIrD8wOgc1gDaujn13+kOzhCxZA4xhGMcqJ1F2uJSmdvC4FOso3JPSZkf6i/gFt6b8Mz0AlPXq81dg3R2fIUnjd/yiMDF7KxRKuOtXT0cKS7T+sQOJxokpXNR5zZmOVWEnIEIlIhIi+KyC7r/yHD+UXkPBF5O+Rfl4gstu57WET2hNw3MxF7DEOJJrPZ2N7l+PwAEExo2gI1zUe62alq2XrWd2H/m/D8lzJsoYM5sAXWfA4mn8NDBTcGK23sv2WdS7aGYHAvQVdvP4Eud4yXcAuJRgRfBl5WSk0DXrZ+H4RSap1SaqZSaiYwH+gEXgg55A77fqXU2+GPNyRGpDETSintCBxeOmozpbIoKFBjr6Pv5EVaxGbzKnhrdSbNcyadrfDEJ6BgLCxdzdjSomBj4R7bEbggIqgu01/2oZVDfjNrKukk6ggWAfZZuBpYHOP4K4HnlFKR1UYMKcFXMri7+LDTJCpjUOfTjkApFewq9hV7YcFdMHU+/P522L85w1Y6iIF+ePrTunnsqv+B4nFUFXuCM5r2tXSSmyNMGlsQ44kyj6/IS16ODNoacltXsRtI1BFUK6XsjF0jUB3j+KuBx8Ju+7aIbBGRH4pI1HdWRG4Wkc0isrm5uTnaYYYI+Io9g5rKGl1SOmpTV1lIV+8ATe3dIXOGPJCTC1f8HErGwxPXQWCodkFW8vI3YPc6uPg+mKQrqyqLjg2e2+OS0lGwJSsHK5X5A+7qKnYDMT8JIvKSiGyN8G9R6HFKD41XUZ4GERkPnAasDbn5K8DJwFlABRB1w1cptVIpVa+Uqq+qcs+IBycQPn3ymDKZO04kewtjj7+DliM9ePJyKPFaZa+FFXD1o3D0EDz1SeiPPLs+a9j6DLz+I6j/NJxxXfBmX4mH1o4eBgYU+1o6XJEfsBkX1lQWHENucgRJI6YjUEotVEqdGuHfs0CT9QVvf9EfHOaplgK/VkoFz1Sl1AGl6QZ+AcxKbDmGSIRPILXl8pw4/CoS9pfWvpYOmo904yvyDC57rTkNLv1veO+PsPbODFnpAJq2wbPLofZsuPCeQXdVFmlhpbajvez1d1LngtJRm5rS/EGaBH4XjaB2C4nGhmuAG6yfbwCeHebYawjbFgpxIoLOL2xN0B5DBOwxE919esxEY5s+kdySI5hQXoAnN4c9LToiiHglOGMJzP4svLkC3v5V+o3MNEcPweOfAG8pLP0l5A3+krT/Zn9tDLimdNSmpmywZKXfheMlnE6ijuB7wPkisgtYaP2OiNSLyM/sg0SkDqgFXg17/KMi8g7wDuADvpWgPYYI2HupdrKwsb2LyiIPnjzn7xGDXUJawD5/5/AD5xZ+HabMg99+Hv7x5/QamUkG+uF/b4K2/doJlNQMOcS+en5rXyvgjoohm+rSfALdfXR09wGWVrHZFkoqCX0TKKValFILlFLTrC2kVuv2zUqpfwk5bq9SaqJSaiDs8fOVUqdZW03LlFJHErHHEJnw7mItUemOaMCmzioh9R/pjq5DkJsHV/5CK249vgw6/Ok1MlO88l1oeBEu+k+YfHbEQ2znuXnfIcAdPQQ2di7Lzm01m+mzSccdl4SGhLCvnmxH0Njmnh4CG7uEtLUjxtVgkU+XTHb6reRxX9pszAg7fgvr74UzroczPxX1MPuL8619h1xTOmoTlKy0tof8DpV7dDPGEWQB9kljJ4xdGRH4iujqHaC3X8VWJpswEy75Eex9DV68Kz0GZoLmd+HX/woTz4SLvg/DzI0qLxhDbo4Q6Oqj1iWlozbh3cXNLpw86nTc82kwHDdVwYigh+6+flocLlEZidAql7i+BGZeA7NugY0/gS1PpdCyDNHVBo9fC2MKYOkjkDf83yQnR6iwHKgbpo6GUh3iCOzxEqZiKLkYR5AF5I/JpdibR3Og+5hEpUt6CGxC97Tj1iK44Nsw+Rw9b+fAlhRZlgEGBuCZW+DQXliyGsomxvUwO5JysjxlJIq8eZR482hq6wrqapgcQXIxjiBL8BVr7eKDAauHwGURgV1CCsekF2OSOwaWrtbzdp74hJ6/MxpYfy/sfA4u+C7UxT/f3o6knC5PGYnqMt1LYLqKU4NxBFmCr1g3ldk9BG7bGrJLSGGEXwLF43TyONAIT3/K/cnjd5+HV74Dp18Ls24a0UPtv5ubSkdtaiylMjvPZcpHk4txBFlCVYkeM+G2OUOhTPEVkSMwtnCE+8OTzoSLfwC7X4EX/0NvrbiRpm3wzE0w/nS45L5hk8ORCG4NuSxHAMe0i+3KN5MsTi7GEWQJet5QD03tXXjycigvHJNpk0bM7Kk+zjxh7CDJwrg543o46ybY+FNYfQm07k6+galioB82/AhWngd5+TrCGTPy8s+Zk8uZWlXERBeVjtrUlHk5GOjmoBURxKwcM4wIZwvWGpKGr9hL29Fe3m/tpKbU+RKVkbhx7hRunDvl+J/gont1aenz/w4PzIEFX4NZN0OOg6+HDv4Vnr0V/v4WfOjjeqJo8bjjeqpLZkzgkhkTkmxgeqgpzad/QPFuY4CS/Dzyx5jxEsnEwWeAIZnYCdZt/2h35bZQUhCBDy+DW9+Aurla3ezhi6Hlb5m2bCj9ffDafbDiXGjdA1eu0mWix+kE3I5dQvrO39tMM1kKMI4gS7AThe+1dlLtsq7ipFM2Ea59EhY/oPfdH5gDb/zUObmDgzvg5+fDy1+Hky6E5W/CqVeMOCcwmrA74d9r7TQVQynAOIIsITS5VlNqTiREYOa1sHyjHlS39ivwi38Gf0PmbOrvg/XfhxXz4PA+WPIwXPUIFBv9jdBOeJMoTj7GEWQJoeG028ZLpJTSCXDtE3DZCmjeAQ/OgT/+WCdo00nTNvjZAvjDN+Hki3UUcMpl6bXBwfiKvcEiAdNVnHyMI8gSfMYRREcETr8abt0EHzgPXrgTVl0I/l2pf+3+Xnj1XljxUT1GeslqHQkU+VL/2i4iN0eCFzNmayj5GEeQJRR4ciny6EoLt00eTRul4+Gax+Dyh8C/Ex6cC6/fn7rooHErPDQf1n0Lpl9qRQGLU/NaowA7t2WayZKPcQRZhH0CZW3VUDyIwIylsHwTTJ2vG9BWXQDNO5P3Gv298Mo9sPJjEDigq4GuXAVFlcl7jVGIndsyEUHyScgRiMgSEdkmIgMiUj/McReKyLsi0iAiXw65fYqIbLJuf0JEzOZfCrFD63EmWRybkhq4+ldw+c+gpcGKDv4r8ejgwBZ46Dw9JuKUxXo7avqlybF5lGNfwJhkcfJJtKFsK3A5sCLaASKSC/wEOB/YD/xJRNYopbYD9wA/VEo9LiIPAjcCDyRokyEKvmIvFUUeo/UaLyJaC3nKPPjdbVrbYNMK8JYc/3O2NEBBBVz1KHzokuTZmgUEt4ZMsjjpJOQIlFI7gFhdqrOABqXUbuvYx4FFIrIDmA9cax23Grgb4whSxg3n1HHeyaYUccSUVOuxDtuege1rAHX8zzV1Psy7AworkmZetvDxGRPo6h1gYrn7RmQ4nXSMmJgIvB/y+37gbKASOKyU6gu5PepgdRG5GbgZYPLkyamxdJQeBaSfAAAEw0lEQVQze2ols6eafejjQkQ3dZ16RaYtyVpqKwq57fyTMm3GqCSmIxCRl4CaCHfdqZR6NvkmRUYptRJYCVBfX5/AJZnBYDAYQonpCJRSCxN8jb8DtSG/T7JuawHKRSTPigrs2w0Gg8GQRtJRPvonYJpVIeQBrgbWKKUUsA640jruBiBtEYbBYDAYNImWj14mIvuB2cDvRGStdfsEEfk9gHW1/1lgLbADeFIptc16ii8Bt4lIAzpn8PNE7DEYDAbDyBF9Ye4u6uvr1ebNmzNthsFgMLgKEXlLKTWk58t0FhsMBkOWYxyBwWAwZDnGERgMBkOW48ocgYg0A/uO8+E+wJ9Ec9yAWXN2YNY8+kl0vScopYaMF3ClI0gEEdkcKVkymjFrzg7Mmkc/qVqv2RoyGAyGLMc4AoPBYMhystERrMy0ARnArDk7MGse/aRkvVmXIzAYDAbDYLIxIjAYDAZDCMYRGAwGQ5Yzah1BNJ3kkPu9lk5yg6WbXJd+K5NLHGu+TUS2i8gWEXlZRE7IhJ3JJNaaQ467QkTUcNrabiCe9YrIUut93iYiv0q3jckmjs/1ZBFZJyJ/tj7bF2XCzmQiIqtE5KCIbI1yv4jI/dbfZIuInJHQCyqlRt0/IBf4G/ABwAP8BZgedsytwIPWz1cDT2Ta7jSs+Tyg0Pr5M9mwZuu4EmA9sBGoz7TdKX6PpwF/BsZav4/LtN1pWPNK4DPWz9OBvZm2OwnrngecAWyNcv9FwHOAAB8BNiXyeqM1IgjqJCuleoDHgUVhxyxC6yQDPA0skBjiyw4n5pqVUuuUUp3WrxvRYkBuJp73GeCbwD1AVzqNSwHxrPcm4CdKqUMASqmDabYx2cSzZgWUWj+XAf9Io30pQSm1Hmgd5pBFwC+VZiNa5Gv88b7eaHUEkXSSw/WQg8corZnQhtZEcCvxrDmUG9FXFG4m5pqtkLlWKfW7dBqWIuJ5j08CThKR10Vko4hcmDbrUkM8a74bWGZpo/we+Fx6TMsoIz3fhyUd4vUGhyEiy4B64KOZtiWViEgOcB/wyQybkk7y0NtDH0NHfOtF5DSl1OGMWpVargEeVkr9QERmA4+IyKlKqYFMG+YWRmtEEE0nOeIxIpKHDilb0mJdaohnzYjIQuBO4FKlVHeabEsVsdZcApwKvCIie9F7qWtcnDCO5z3ej5aC7VVK7QF2oh2DW4lnzTcCTwIopd4A8tHD2UYzcZ3v8TJaHUFEneSwY9agdZJB6yb/QVlZGJcSc80i8mFgBdoJuH3vGGKsWSnVppTyKaXqlFJ16LzIpUopt8rbxfO5/g06GkBEfOitot3pNDLJxLPm94AFACLyIbQjaE6rlelnDXC9VT30EaBNKXXgeJ9sVG4NKaX6RMTWSc4FVimltonIN4DNSqk1aH3kRyy95Fb0B8y1xLnme4Fi4CkrL/6eUurSjBmdIHGuedQQ53rXAv8kItuBfuAOpZRrI9041/xF4CER+QI6cfxJl1/UISKPoR26z8p9fA0YA6CUehCdC7kIaAA6gU8l9Hou/3sZDAaDIUFG69aQwWAwGOLEOAKDwWDIcowjMBgMhizHOAKDwWDIcowjMBgMhizHOAKDwWDIcowjMBgMhizn/wEFeF4tUi6OlAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "kdx = np.linspace(0, np.pi, 20) / np.pi\n", "plt.plot(kdx, np.cos(2*np.pi*kdx / 0.001))\n", "plt.plot(kdx, np.cos(2*np.pi*kdx / 1.0))" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "def ampFTBS(jhx, lmbda):\n", " return np.sqrt(1 - 2 * lmbda * (1 - lmbda) * (1 - np.cos(np.pi*jhx)))\n", "\n", "def ampBTCS(jhx, lmbda):\n", " return 1 / np.sqrt(1 + lmbda**2 * np.sin(np.pi * jhx)**2)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false, "jupyter": { "outputs_hidden": false } }, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, '$|s(j)|$')" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEYCAYAAACUdWs9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOydd5xcVd3/32f6zu5s7303yaaHCARCiYAIAkoogjSxQrAgKD4oNkBRAX3Q50Ee5CdFVJQmKKE3CQHSQ+puNptsS7b3Ov3e8/vjzia7yc6W7GzNeb9e93Xbufeem+zMZ875NiGlRKFQKBSKSGCa7A4oFAqFYuagREWhUCgUEUOJikKhUCgihhIVhUKhUEQMJSoKhUKhiBiWye7AZJKcnCzz8/MnuxsKhUIxrdi6dWuLlDJlsHPHtajk5+ezZcuWye6GQqFQTCuEENXhzqnpL4VCoVBEDCUqCoVCoYgYSlQUCoVCETGUqCgUCoUiYihRUSgUCkXEmBaiIoR4QgjRJITYHeb8dUKInUKIXUKIdUKIEya6jwqFQqGYPi7FTwIPAX8Nc74SOEtK2S6EuBD4E3DqBPVNoVAoRo2UEjSJ7tOQ/r5FRwY09IBubAeNfRnQIagjAzpSk8igDppEajoyKEE3ttElUie0liAlyNC+BAjth7CmOUm4bE5E32taiIqUcq0QIn+I8+v67W4Asse7TwqF4vhG6hLpDaJ7guju0NoTQPdo6N4g0hNE9wbRvZrRzhs67tMMIfFpoE9y6ZFxeP60EJVR8nXg9XAnhRCrgFUAubm5E9UnhUIxxZGaRO8NoPX40XsD6D0BtJ6Asd0bQAutdXffEhzwq39aMg79n1GiIoQ4B0NUzgzXRkr5J4zpMU4++eTp/iehUCiGQQY0tC7/4aXbWOvdoe1uP3qPIRTTXiSmADNGVIQQS4DHgAullK2T3R+FQjH+6L4gWocPrdOP1ukj2OFD6/QZ4hFaS09wsrsZHpNA2MyY7CaEzWwsVtPhtdWEsJgOb5tD+xYBfdtmYdzHEloLY41JIEwCBCCMtRB9DzY2hMMc8VeaEaIihMgFXgSul1KWTXZ/FArF2JFSovcECLZ7DeFo9xqi0e5D6zAERHonVzCEzYzJaTGWKAsmpxWTw4KIMhv7DmMRDrOxtvdb280Iy7RwwB0V00JUhBBPA2cDyUKIGuAuwAogpXwEuBNIAh4WhhQHpZQnT05vFQrFSNG9QYKtXoJtXkM0+tbtXrR2n+H1NIGYnBZMMTbMMVZMMVZM0VbM0Ye3Tc7QvtOKyWmZkaIwVqaFqEgprxnm/A3ADRPUHYVCMUKklOjdfkM4Wj2HBCTY6kFr8xrG7vHGBGaXDVOsHXOsDbMrtMTaMLkMATG7bJhirAizEomxMi1ERaFQTF36hCPQ7DGEoyUkGiERGdfRhllgjrNjjrNhibNjjrdjjrUfOmaOs2OKthq2BcWEoERFoVCMCN0bJNjsIdjiIdDsJtjiObRI//gIh7CbsSTYMSc4MMfbsYTW5ng7lniHMbpQgjGlUKKiUCgOIXWJ1uEj0OQ2BKTZbYxAWtzo3YHIP9AssCQ6MCc4sCQ6Bm4n2BFRFoRQojGdUKKiUByHyIBOoMVDsMkdEhA3wSYPgRYPBCM76hB2M5bkKCxJDiyJxtqc6DDWsXY10phhKFFRKGYwujdoiEaTm0BTSESa3Wht3ogG+g0QjuQoLElRh/ZN0VY12jiOUKKiUMwAdE9IPBqNkUegsZdgkxut0x+5h5jEYdFIcWJNCYlHSpRh21DCoUCJikIxrdC9QQKNIdE4JCBu9K7IiYcp2oIlxYklOQprqjMkIlFYEqOM6G2FYgiUqCgUUxDdZ4hHsNEQjUCTm2BDL1oExcOcYDdEI8WJJTXq0LY52hqxZyiOP5SoKBSTiO7XCDZ7CDT0HhKOQKMbrcMXmQeYCE1RObGmOQ3hSHUaU1a2yOd9UiiUqCgUE4AM6oZrbmNvaPrKmMKKmMHcJIzpqjRDPCypzkNTVyqViGIiUaKiUEQQqekEW70EGnsJNIS8rhp7CbZ4IBKeuoOJR1pIPFSKEcUUQImKQnEMSF0aOaxC01WG51UvgWYPaBEYeggOG8rTnFjTog+Lhxp5KKYwSlQUiiGQukRr84aEo/ew8bzZDcHIiIc5wREaeYSEI82JNcWJsCrxUEw/lKgoFBilZINtoeDAQyMPI0VJpCLMzXF2QzDS+wlIqlMZzBUzCiUqiuMKGdAJtnoOC0dzyO4RqWkrwOSyHhINa1p0aPrKicmhPm6KmY/6K1fMSHR3wPC2CgnGeKQnMUVbD09X9Zu+MjlVnIfi+EWJimLaIjXdMJaH0rEHm0Mp2Zs96L2Ry6hrirZgSY0+yuPKHGOL2DMUipmCEhXFlEbqEq3bP6B2R5+ABNu8oEcuK6IpxjogxqNvW4mHQjFylKgoJh2pS7RO38CSs62GeGhtka8caI63Hw4ODKUnsaaqaSuFIhIoUVFMCLo3SLDNi9YeqlHe5iXY6kVr8xJs90bMSH4Iszgc55EShTVFpSdRKCYCJSqKMSOlRHqCBDt8aB0+QzjafWgdoXW7F90dHJdnm1xWLMmhNOwph1OymxMcqviTQjEJKFFRDImUEunV0Lp8aJ1+tE5faPET7AyJSIcP6dfGrQ/CYcGSbNTxsCaHij+FFuWmq1BMLdQn8jhF6hLdHUDvCaB1+9F6AuhdfmO724/W5Ufv9qN1+ZD+yNo0BsMUbcGSGIU5yXG4amCiISQmp6pTrlBMF5SozACklEifhu4OonuChli4A8Z+bwCtN4Deb9F6jHUky8kOi1lgSXBgTrAbYpHowJwYEo4khxpxKBQzhGnxSRZCPAF8DmiSUi4a5LwA/he4CHADX5FSfjyxvRwZUpeg6ciAjgxKZFBHBjRj36+jBzSkX0P6daRfQ/dpyNCie4PoXg3pDaKH9qXHEJIJFYhBEFYT5ni74VkVb4iHOcGBJd6OOdGB2WVTNg6F4jhgWogK8CTwEPDXMOcvBOaEllOBP4bWEeetb/0PmY68fkcE4tBaIETf2oQJEyZh6rdtxiSmYZJAiwlznA1zrM3IXxVnx9y3hIRETVEpFAqYJqIipVwrhMgfosklwF+llBLYIISIF0JkSCnrI90Xh4gm3pYc6dtOCsJhxuyyYYqxYXZZjW2XDbMrJCCxxraIUoKhUChGxrQQlRGQBRzst18TOnaUqAghVgGrAHJzc4/hUZM8zxSGgO4noHvxa158uge/7sWvefDpbnyaB5/mxqcba6/mNvaFCb81Cs3mRDhiMDtjsEa7iHLF4kqIJy4hkaSURNJSU8jKSCXW5Zzs11QoFFOcmSIqI0ZK+SfgTwAnn3zyqBUiEpKi6UE0aSy61AjKAJoMoMkgQT1IUPrR9ABB6Seg+wnKAEHdT0D3GYv0E9B8+HVv6JgX/RjKClqlhtUfAH8X9Bw+HgTaQ0tVv/YBYcVni0Z3uDBFx2JzxeGMSyA2KYnk1BTS09PIzU4jJTkBk2kaTvMpFIoxM1NEpRbI6befHToWcexnJVNcvRMpdaSUICW6pqNpGlIPomlBtEDAWIIaejBIMBA02gSNJIhSFwhpRkozEhMIM7owoZkEuhCAhpQBIAgyGNoOgAwY29Jv7E8wVhnA6usAXwd0GseCQFtoKQu1CwozPpsLLSoWc0w8jvhEXEnJJKWmkpGVTn5eNlnpyUp4FIoZyEwRldXAzUKIZzAM9J3jYU8BWP6FyyN+z6AWpL2jnra6Ctrrquk6WEdvQyu+Zj96RxDhtmLxx2CWcQRssfhtcXhtsfhsTqBPaHxI6Q+tfUesvaCH1tJ7aD1eU3kWqWHpE5+OA1BjuOS5MeYoN2EIj9fmQnfGY4lNxJmYTEJaGhlZmeTlZVOYn4XDrhI5KhTTjWkhKkKIp4GzgWQhRA1wF2AFkFI+AryG4U68H+O766uT09Njw2K2kJKUQ0pSDiwO367T005d1W6ayovprPgIX/UB9AY39nYLTk8suiURr8NYPI4svI4kdPPgX8zGKMuHlB6QHqTu6bftNrZ1N1K6kbobpBuOYYot7DtLjZg+4WmvgurDU24lgI7Aa40h6IzHEpeEMymVpIwMMrMzmTUrj/zsDMwWlcNLoZhqCMNh6vjk5JNPllu2bJnsbkSEDm8HVRUfU1+8ma6yPegV1TgPtpLS5gRzCp6oFNzOVDxRqbidabijkpGmkWflNUTIGxKZXtB7QXbhoB4rLaB3EdT9uINmgnL8v+w1THgc8UhXIvbEVOJSM0jPzqKgMI+5c/JwRjnGvQ8KxfGKEGKrlPLkQc8pUZkZojIYUkpaPC2U7fmIhh0bcO8pwVZeQ3qth+QuE15HEr3OdHqj0+l1ZtIbnY47OgN9FGJzJDa6SLHsId5cSgzVWGjGEzTRoUXTHIynPWDH7weTPn65wiTgsboIxiRhS0wlLi2T9JwsZs0uYO6cfCU4CsUYUaIShpkuKuFo97ZTUrGBg5vfp3fXDpxlteTXBIh3gy5MeKJS6InOoifGWLpjcvDb44/pWQKNZEslmbYSMmx7yLDuwWnuxKeZaTenU2+bQ42eTq03hpZeE4HuTqyeTqKCvRF+awMdgccWi+ZKxpGcRkJGFll5uRQVFTI7P1tNqSkUI0CJShiOV1E5El3qVHZUUrzrPzRt+gB2lZJd0U1OM/T5Z/mtMXS7cumOyaHblUOXKw+fI/GYnpdgPkimrZgs226ybLtxmkOuZAhIXQA5p+BOPZEKPYeKVklDbR0dTQ142prRu1qxu9ux676IvHt/gsKMJyoREZ9KTGoGqTm5FM4qZP78WSTGuyL+PIViuqJEJQxKVMLT2NvIx/vXcmDd22gf7yRnfyeFDWDq9+fit7rois2jy5VHV2w+XbH5BC2jD5BMsBwg27aTbNsusmy7sZvch09Gp0Lucsg9DfJOg7TFYLbQ1NLB/vIqqqsO0lRbS3dTA4H2Ziy9bTgD3RH4FxiI2xpD0JWCIzmdxKxscgsKmD9/DrnZaco1WnHcoUQlDEpURk59Tz2bytdQ/cEbmLbsYna5h7zmgW0kgl5nOl1xBXTGFtIZV4jbmTaq5wg0Uq37ybHtIMe+nTRrGWbRz/5ic0HuqZB3OuSdCZmfAMtAD7eeXjdl+6qpqKii4WAtnY11+NuasPS0RlxwfCY7vugkLEnpJGTmkF1QwLz5s5lTkKOm0hQzFiUqYVCicmzoUqe0rZRNu9+k6f13iN9ZxZIKnTj30W391hg642bRETeLjvjZdMfkwCiSalqFm2zbLnLt28i1fUys5Qglszoh5xTIXwEFnzRExhze0aCr203pvgoq91fRcPAgXY11BNubsPe2Yte8I+7XcASEBW90MubEdGIzsskqKGDe3NnMLcrDapkWnvwKRViUqIRBiUpk6PJ3se7gh+xet5rAuk0s2Othdt1he0x/gmZHSGDm0B5fRLdrdCKTYK4h176VPPtWMm17MIsjyhTbYoxRTMFZhsikLYIRTk/VNbSwp3Q/1RWVNNfU4GmuR3Q14/R2YIpQoGhQmPE4kzAlphOXkUNWQT7z5hUpsVFMK5SohEGJSuQJ6kG2N23nw92v0vaft5m1u40TKiT2MCXqA5YoOuJm054wl7ak+bij0kf8LKtwk2vbTl5IZA4b/PvhTIbCs6DwHJh1DsRlj/qdPF4fe0or2L+vgrrqaroaatDbm3D0tmKVkUmXY4hNMqbEdOKzcskuKGD+giKKCtU0mmLqoUQlDEpUxhcpJWXtZby773UOvPsKudvrOLlM4hpilslni6MtYS5t6SfQkbQQnz7SmBmddGsZ+fbNFDg2kWCuYdBs/UlzYNanjCX/TLDHHMurAaAFNcoP1FK6Zz8HKyrpqKvB31KPvac5Yt5pAWHFE52MJSmDhKwccmYVsmjRXApyM5WDgGLSUKISBiUqE0t5Rzlvl79B2Xv/JndrLafulcQPYofpQyLoicmiffYZdGYvo8UdhRxhppg4cx2Fjo0U2jeSZi1DiEH+zk1Ww7Ns9qdhznmGO3ME6sbouk5NXTPFJWVUl1fQVnsQX3Mdtu5mHJpnzPcH8JtseGNSsSVnkJidS97sWSxcOIecTOWNphh/lKiEQYnK5CClZF/HPl4vf5V9773EnI+bWF4qiR3m+zZgjcK9/BI68k+jriMKT/fIpp6cpjYK7RspdGwgy7YbkwijTK5MmH0uzDkfCs8GR+yo3msk1NY3sWt3GdXl5bTW9IlNE44IOQl4zQ78rlTsKZkkZedSMGcWixfNIyMtKSL3VyhAiUpYlKhMPlJKdjTv4NWy1Rx47xU+saOHU8okUf5hrouNRVxwDe0Fp1NTL2g+MDJXYYfoosCxiVmO9WTbdh5t6O/DZDFiY4o+A3M+A8lzIjKKCcfBuiZ2795LdXkFrTUH8DfXYe9uitg0msfiJOBKxZGaRUpOHgVzClm8aB6pyceWKUFxfKNEJQxKVKYWAS3A2tq1vFbyL9xr1nLG7iBLyyXmYf5E7YsWYl95Na3pJ1Jd2kXN3nZ0bfi/a7vopsCxidmOdWTbdgyMhzmShAIousAQmbwzjoqNGQ90XedATSO7i/dSvb+c9tqDBFrqcPQ2Y9Mj4yDgtsYQjE0lKiXTEJuiWSxZVERyohIbRXiUqIRBicrUpc3bxmsVr/Hm9ufIWFfOWbt0ChuHvkY4o4i7eCXRl11JYyCZiu3NVO9uJeAbPnml3dTNLPsGZjs+HHqKDIwAzNnnwtwLjaky57GlqzlWdF1nf1UNJcVlHCyvoKP2AFpbA1G9LVhkmJHXKOkvNsk5uRTOmcWiBUWkpiRE5P6K6Y0SlTAoUZn6SCkpaS3hxX0vsnP9ak79uJcVxcPbXxwnnEDitdcQde751FX0Ur6ticodLfjcw3/pRlm6mW1bS5HjA9Kse4ee9RImY5ps7kUw7yJILBzdC0aQQDDI3rJqSkvLqK2sorPuAHpbA053K+YI1cLxWKIJuFJwpGSQmJVDbmEBCxbMVg4CxxlKVMKgRGV64Q64ebPqTV4seRbH+l18aodkSaUcNMiyD3NiAvFXXEnC1VdhSkundm875R83U7GtGW/v8FNILnsnRdZ3KXK8R6KlZvhOpsyHeZ81BCbzxHG1w4wUvz9AcWkFZaX7qKuqoqv+ILK9EaenLWJi4zPZ8cUkY01MJz4zm6yCPIqKZlE0WwV1zkSUqIRBicr0paS1hOfLnmfjxy9z+lY3n9ohSRgqW77JhOvTnybxS9cTddJJSF1SW9bB/q1NlG9rwtc7/Agm2dVBke1t5pjeIMbcNnwnXZkhgfmsERMzRPqYycDr81Oyp5y9e/dTV1lJT0MNsqMRp6c9YmKjYcIdlYiISyE6NZOU7BzyZuWxYN4c5SQwjVGiEgYlKtOfbn83q8tX81zx0yR/XMn5H0uWVA39N21fMJ/E679E7GcvwmSzoWk6NaXt7N/cSMX2Zvze4W0w2antzLW/S6HvBWymEbgDO+ING8z8i43AS2vUSF9xwvF4fZSUllO2t5z6qiq6G2oiPrIB8Fqi8EcnYU1IIzY9k/ScHAoK85k/N5+Y6NFnu1ZMHEpUwqBEZeagS531dev5R+k/2Lfjfc7fpnP2Tkn0EB655pRkEq+9lvirr8aSYBiggwGNA7vbKNvcSNWuFrTA0F+iFptgVk4Hc21vkd3+HIIRGMqtTiPYcv5Kw9A/DvEw40HfNNq+snLqq6vprD+I3taIw92KNUIOAjBY5c4M0rKzyS/IY+6cfGJdSnAmGyUqYVCiMjOp7qrmH3v+waslL3LyDjcXbtXJbQ7fXjgcxF9+GYlf+Qq23NxDx/2eIBXbmynb1EBNaTvDfVRi4q3MLexgnvV14uv+BYEh0gX0YbYZI5f5K42RzAR7kkUCLaixv6qG0tL91FRV015bQ6CtAVtPM45gZDII9MdjiSYQnYAlPoXopFSSMjLIysmisCCXvOx0lSttAlCiEgYlKjObLn8XL5S9wN9LniJ5TwMXbZGctG8Iw74QuM47j6QbbyBq8eIBp3o7fezb3EjZpsYRBVpmFLqYN6uN2fI1bOWrwTMCG4zJYmRXXnAJzPscRE//KPjGplaKS8uprqii+eBB3M11iM4WoryRs9v0JyjMeO2x6NGJ2OKTiElOJSktjczsTPLzssnJTlOOAxFAiUoYlKgcHwT0AG9WvcmTu5+ko6KUC7fonLNz6Kh956mnknTDDUSfeQbiCA+u1roeyjY2sHdjI70dQ0e8W+xmZn8imfmzWsjoegVR+gp01w3faWGGghUhgbkYYlJG8qrTBr8/QNn+A+wvr6Su+iAdDbX42powd7fg9HcxXj5zGia8thi0qDjMrgQc8UnEpSSTlJpKWkYauTkZZKenqtHOMChRCYMSleMLKSXr69fz591/Zmfles7dLrloi07SEAMPx4IFJK26Edd55yHMA79odF1SW9pO6YZ6KrY1ExzG/hKf5mT+aenMzW8i+sDLsGc1tFcN33FhMrzHFlxqTJPNMIE5kp5eD2X7qqioPEDDgYN0NtXjb2/G3N1KlL8rYrVtwqFhwmd1EnS4EM44rK44ouMTcSUkkpCcREpaMulpKWRlpOCMcoxrX6YqSlTCoETl+KW4tZjHdz3OexVvc9oenYs36uQ3hW9vKygg6aZVxH32swjr0a7Bfk+Q/R83Ubq+nvr9g9R16YcwCfIXJ7HgjAxyk+owla42BKalbPiOH2cCcyRen5/95QepqKymvqaWzsYGPK1NyO42bJ6OiOVKGyl+kw2/LRrdFo1wxmB1urDFuHC6YomOiyM2Lo74xHgSE+NJSYonKTEeh338U/yMNzNCVIQQFwD/C5iBx6SU9x1xPhf4CxAfanOHlPK1oe6pREVR2VnJE7uf4JX9L7OwIsAlGySLqsN/JqxZWSStWkX8ZZcibIN/OXQ0uildX0/phoZhp8diEuzMOz2DBWdk4gpWQcm/oeQlaCoZvvPCZJRRXniZITAzwAYzVppaOiivrKbmQB1NdXV0tTTha2+Fnnas3s6IZYMeCwFhJWCxo1kcSKsDYXMg7FGY7VFY7HasjijsUU7sTgd2RxRRTieOKAdRUQ6iHFFERdlD+3YcdjtRDhsOhx2b1TJhWQ2mvagIIcxAGXAeUANsBq6RUpb0a/MnYJuU8o9CiAXAa1LK/KHuq0RF0UddTx1P7H6CF/e9SG6Nn0vX65xaFv6zYcnIIOnGG4i/4gpMYcRF1yUHS9rYs66Oyh0tQye5FJC7IImFKzLJX5yEqW2/IS4lL0HjruFfQJiNCpcLLzeCLaehF9lE0NnVQ0VVDTU19TQ3NNHR3IS7vZVAdzuitxObr3vCRzuRRMOELsxIIdCFCSlMgEAK0W/dh0BPzODu//3vUT9nJojKacDdUsrPhPZ/BCClvLdfm/8HVEgp7w+1f0BKefpQ91WiojiSJncTf979Z54ve57kRi+Xrtc5szh8pmRLWhpJN60aUlwAPN1+9m5soOSjetrrhwr9h+g4G/PPyGTBmZm4Eh3Qsh9K/gXFIxQYk8VwU154mZGTLEpFro+Grm43B2vrqa1toKWphfbWVnrb2/B1dRDs7UJ4urH4e7EHPeNu3xlvulwZ/PyxR0d93UwQlSuAC6SUN4T2rwdOlVLe3K9NBvAWkABEA5+WUm4d5F6rgFUAubm5J1VXV0/AGyimG83uZp7Y/QTPlz1PbKuXSzcYwZTWMMH2lvR0km9aRfznPx92WgwMZ4HGyi5KPqxj35ZGgv7wxn0hIHdREotWZJG7KAmTSYxeYMw2o7Llwsth7gVgdw1/jWJEaEGNhuZWGhtbaWpupb2tna6ODtxdXXi7Own09qB5e8HrxhzwYAl6sWq+KSVEXa5Mfv7Yn0Z93fEiKrdhvM8DoZHK48AiKcMXoFUjFcVwNLubeXz34zy/93liOnxcskHn3O0SWxhxsWZmkvztbxF3ySWIYeIh/J4gZZsbKfmwbtjYl5hEOwvPzGL+GRlEx9mNgy37oPhfxjISG4zFYUTwL7rcKDxmU5HpE40W1Ojo7qGltYO29g66u3ro6e7F3duLu6cXn8eN3+Mh4PWg+X1ofi+634cMBiAYQGihRQ9i0oOYdA2z1DBJ7ZjE6ngWlZFMfxVjCM/B0H4FsFxKGdanR4mKYqQ09Dbw2K7HeGHfC8R0Bbhkg8552yS2MNlJrHm5pHznFmIvuhAxAuNpU3UXxR/WUbapkeAQ9V9MJkHB0hQWnZVFVlH84Riapj2HBWYkXmRWpxHBv/ByYyRjPT5dY2cSWlDDFwjg9QUIBjWCwSB+fwBd19F0HU3TDmWFkFIipY7DYWdWfvaonzUTRMWCYag/F6jFMNRfK6Us7tfmdeBZKeWTQoj5wLtAlhziBZWoKEZLTXcNf9zxR16peIXYbo1L1ofEJYwO2IuKSLn1FmI+9amjgigHw+8NUrapkeIPamk52DNk24R0J4vOymLu8gzsUaFRkZTQWAzFLxoC01Yx/EvZYw3by6LLofCcCalqqZjeTHtRARBCXAT8D4a78BNSyl8JIX4BbJFSrg55fD0KxGDkpPuBlPKtoe6pREVxrFR0VPCHbX/gnQPvkNAtuXS9zqe3h7e5RJ1wAinfv43oU04Z0f2llDRVdbP7g1r2b24cMrDSYjNRdEo6i8/OIjnb1f8mUL/DEJjd/4LOA8M/2BFnRPAvusxIGTPF0vUrpgYzQlTGAyUqirGyq3kX//vx/7KxYSNJnZLL1xkpYCxhNCB6xQpSb/sejvnzR/wMnztA6YYGitfW0t4wdJLKjFlxLDori1knpmK29Jt2kxJqt8Lu0AhmJKliohJhwUrDiyx/BZhU6hKFgRKVMChRUUQCKSXr6tbx+62/Z2/7XtLaJF/4UOeM4jDJK4Ug9uLPkXLLrdiys0b1nLqyDna9X0vl9mZ0PfxnN8plZcGZmSz6ZBYxCUfYS3QdDm4MTZH9G3qHSCXQR3SKkYds4WVG+WQlMMc1SlTCoERFEUl0qfNqxav8YdsfqO+tJ6dJcvVanWX7Bv+MCauVhOuuI/kbN2GOH10sSW+nj5IP6yj+oG7IqH1hEhSekMyis7MHGvYPdVqD6o+MEdZQygQAACAASURBVEzJSyPLphyTdlhgcpaDqk1/3KFEJQxKVBTjgU/z8fc9f+exnY/RHeimqEZy7RqNBQcHb2+KjSX5pptI+OJ1mOz2UT1L13Sqdray6/0aakrbh2ybkBHNkrOzKDo1HZtjEHdnLQCVa40RzJ6XwTt0DjMAXBlGipiFl0HOqUpgjhOUqIRBiYpiPOnwdvDIzkd4tvRZgnqAT5RLrlsTvmCYNSuLlNu+R+xFF43IU+xI2ht62fV+LXvX1w9ZEtnqMDPvtAwWn5VFQnr04I2CfqhYY9hfSl8BX9fwHTgkMJeqEcwMZ0yiIoQYvmD32JBSykmpmqNERTERVHZW8rstv2NNzRqELvnkbslVH+gkh/medpywhLQf3oHzxE8c0/P63JJ3ramhrW7olDA5CxJZcnb24Yj9wQj6oPw/xhTZ3tfAP7SrMwAx6YaRf8ElygYzAxmrqES+PNsRSCkn5SeNEhXFRLKhfgP3b7qf/R37sQYkF22RXLZexxnGJOK68AJSv/9fozLm90dKSd2+DnatqaViezNyCMN+bLKDRWdlM//0DBzRQ7gRBzyw/11jBFP2xsgEJjrFqGS5YKXhRabclKc9avorDEpUFBNNUA/y4r4XeWjbQ7T72nG5JVd8aARQDuaGLGw2Er/yFZJWrcIcE2aqagT0tPso/qCW4g9q8XQHwrazWE0UnZrOknOyScqKGfqmAQ/sf8cQmL1vQGDoUREAUQlGoOX8i41ASxXJPy1RohIGJSqKyaLT18kjOx7h6dKn0aRGRqthbzklTLp9c0oyqd+7jbhLLxlR2pdwaAGd/R83sWtNDY2VQ9tJMufEs+ScbApOSMZkHuaZhwTm3yMfwdhiYM55xihmzvngiB3FmygmEyUqYVCiophsyjvKuW/TfWyo3wDA/AOSL7+jUdg4eHvHokWk/eTHOD9xbPaW/jRVd7HzvRr2bWlED4b/HohJsLPwk1ksPDOTKNcIUrgEvIYNpuQlwwYzEiO/2WZE8M//nDGSiUkdxZsoJholKmFQoqKYCkgpeefAO/x282+p761HSMknd0mueV8nMcwP/tiVF5P6/f/Cmjb2L19Pt5/iD+soXltLT3v4mBezxcSck1NZfE42qXkjHFUEfVDxPux5CUpfBc/Qbs8GArKXwbyLYO5nIaVoZM9STBjjKipCiGggDUgGooBWoEVK2TCmG08ASlQUUwlP0MNjux7jz7v/TEAPYPcbhvzPbRw8YaXJ6ST5W98k8UtfGrKGy0jRNZ3KHS3sfK+Gun0dQ7ZNK4hlyTnZR6eDGQotCNUfGjEwe16BnhF+RSTNhqILjBFMzqlgnhRnUUU/IioqodK+nwPOB1YAC4DBfBE7gXXAWuA5KWXVqB40AShRUUxFqruquXfTvXxU+xEAKR2SL76nc1rp4J9VW34+aT/5MTErVkSsDy01PexaU0PZxoYhk1k6Y20sWJHJohVZRMePInBT16F2S0hgXob2ypFd54g37DBFF8Dscw3Dv2LCiYioCCHygJuB64EUBheSwZCh5T3gUSnlcyO8btxRoqKYqkgpeffAu9y36T4a3YaBZWG1zlfe1skLEzwZ8+lzSbvjR8fsgjwY3t4Aez6qZ/faGrpavGHbmUyCwk+ksPjsbDJmx40ueFNKox5M6atGoGX99pFdJ8yQc4ph5J9zPqQtNMplKsadscappAA/wyjB2zfG3gF8BGwJbbcA7YAXo5xvApAPLAst5wAuDHEpBn4ipXx5LC8VCZSoKKY67oCbR3Y8wt9K/kZQBjHpkvO2Sa5aqxMzyHe8sNtJumkVSV//+qhTvgyFrksO7G5l55oaDpYMnR8sKSuGxWdnUXRKOlb7MQQ9dtbA3tcNI3/lB6CHd4EegCvDGL3MPg8Kz1KjmHFkrKLSjVHzvRL4C/C0lHLfKDvgAD4LXAusBEzA7VLK343mPpFGiYpiulDWXsY96+9he7PxK97lNpJVnrtt8EzI1rxc0n/604hOifXRlw6mdH09gSHSwdiiLMw/LYNFZ2URn3aMpYu9XVD+riEy+94aoaEfECbIOglmnQuzzjG2VdBlxBirqOwB7gX+LqUcc8oWIcRs4A5gv5TyvrHebywoUVFMJ3Sp8+K+F/n91t/T5TfcdAvqJV9/S6MoTHkU13nnkfbjH2HNyIh4f/yeIHs3NrBrTc2wdV5yFiSy+Kws8hYnh08HMxxa0LDDlL0BZW9CU8nIr7W5oGCF4bZceBakzFNTZWNgrKIihirJO4ZOjct9R4MSFcV0pMXTwm83/5bXKl8DQEjJ2Tsl167RiRvku104naR8+9skful6hDXyv9allNSUtrNrTQ1VO1sY6lMdk2hn0SezmH96Js7YMXqsdRw0Ri/73obK9yEwtLANIDrVEJn8MyH/k5A0S4nMKFBxKmFQoqKYznxU+xH3bLiH2p5aAKI9xpTYeR8PPiVmnzOH9J/fjfPEE8etT91tXnavraXkwzq8PeFtISazYNaJqSz6ZNboDfuDEfTBgfWGwJT/Z3SjGDASYOadDvlnGAkwU+arLMtDMNaRym+AncBWKeWecejfpKFERTHdcQfc/HHHH/lryV/RpeH6W1AvufFNjdn1g18Tf+UVpH7/+6MuDDYaRpMOJjEzmoUrspi7PB17VIRiULrqDHEpf89I4e9uGd31jjgjJibnFGOdeSLYh8mFdhwRiSzFEiiVUi4MHfsNsB3D82uPlHLcMxmPB0pUFDOF4tZi7l53N6VtpQAIXfLp7UZU/mBeYuaEBNLu+CGxK1eOfZQwDE3VXex+v5ayzY1oQ8S8WGwmipalsfCTWSOP2B8Jug6Nu40pssq1UPXRyJJf9keYIXUBZJ8EWScbhv/kouM2EHOsovIosAQw992kn9AA+IASDIHpE5odUsoRlI2bXJSoKGYSAT3AX4v/ysPbH8av+wGI7ZVc/x+ds3YP/jl3nracjLvuwpafP+798/YGKF1fz+61tXQ2eYZsm5LrYuGKTOYsSxu8SuVY0AJQty0kMB/AwU2js8f0YXVC+hLIXGqsM5YYDgDHgZdZxG0qQojHMIRmEdA/d3X/mx3gaKGpGPXDxhElKoqZSFVnFXetu4uPmz4+dGxBtc6Nb+hkDRJiImw2kr/5DZK+/vWIpHsZDqlLava2s3ttLZU7Woas82K1m5mzLI0FZ2aSmucan1FV0A/1O4wUMgc2GIt36DQ1YTFZDWFJX2SMbFLnG/tx2TPKEWDcDPVCCBMwFzghtCwNrdP7Nev/gE5gPfBv4C9SSv8xPzwCKFFRzFR0qfPs3mf5/dbf4wkaowJLULJyo+Tyj/RBc4nZZs8i4xe/GFdD/pH0tPso+aiOkg/r6O0In8wSICk7hoVnGqOXIQuJjRVdh+ZSOLgRajYb69b9Y7unLQaS50DyXEieDYmzDI+zxEKwuyLT7wlkwr2/QlH4fQLTJzZzgb5xrAQ+Bs6RUo6g8ML4MJioBAIBampq8HrDp6SY7jgcDrKzs7GOg3upYmpR21PLz9f9nPX16w8dS2+T3PCmzpKqwT/78ddcTeptt2F2TdyXna7pVO9uZffaOg6UtA78KXoEZquJwqUpLDgjg6yiBMSxxr2MBncb1H5sxMnUbDGmz0Zr/A+HMxkS8iA+D+JzIC7HGNm4MiA20zg/xTzRpoRLsRDChjFldgHwLYzMxr+WUv5shNdfAPwvYAYeGyxwUgjxBeBujD/JHVLKa4e652CiUllZicvlIikpadwNmJOBlJLW1la6u7spKCiY7O4oJgApJS/ue5HfbvktvX0GailZUSz50ruDx7ZYUlNJv/NnuD796YntLNDV6mHPR/WUfFSHu3PoyYzYZAfzTstg7vJ0YpOiJqiHGPnKumoNcanfaUyfNeyE7jAud2PBZDHiamJCizMptCQaCTYdcRAVbwR42qINLzWrEywOsEYZ10f4u2xKiMqAhwqxFGOksqfPo2yY9magDDgPqAE2A9dIKUv6tZkDPAd8SkrZLoRIlVI2DXXfwURlz549zJs3b0YKSh9SSkpLS5k/f/5kd0UxgdT31PPz9T/no7qPDh2LcUuuf0/nnJ2Dfw+4zj+f9J/9FEtKykR18xB9o5fiD+s4sLt1yKBKBGTPTWDe8nQKP5F6bDnHIkFvKzQVQ2OxESvTVGoky/R3T05/ABCG84DZZnixCRFaTJBxAlz/r9HfcYzeX4WRNrCHRKINsEgphy28LYQ4DbhbSvmZ0P6PAKSU9/Zr8xugTEr52Ej7EU5Ujocv2+PlPRUDGXTUAiys0ln1hk7GIKm1TLGxpP3gduI+//lJ+7HV0+5lz7p69qyrp7t16Klpq93MrBNTmLs8g6w58RMzPTYUUkJ3A7SUGUtrObSVG+uOAyNPmDkeZJ0EN/5n1JeNVVT8wNMYU1V7R/30gfeyAl8FfojhNfa6lPKGEVx3BXBBX1shxPXAqVLKm/u1+TfGaOYMjCmyu6WUbwxyr1UYGZfJzc09qbq6esD54+XL9nh5T8Xg1PXUcee6O9lYv/HQMWtAcsVHOis3SsyDhJM4T1tOxi9+gS0nZwJ7OpA+z7E9H9VRsb0FLTh0iFxMop2iZekUnZpGUuYUDF7UNSNQs73SSDvTWQOdB4xjXXXQVQ++cYzOyF4GN7wz6svGKiofAqcDOvAB8AzwTyll6wgfLoCzgauBy4FEoBe4Xkr50gjvMRJReQUIAF8AsjGKgy2WUob1DZzKI5U33niDW2+9FU3TuOGGG7jjjjsGnF+7di3f/e532blzJ8888wxXXHHFqO4/Vd5TMXkM5iEGkNsk+cZrg0fkC4eDlFtvNfKImSdpiimEtzfA/i2N7FlXT1P18NNLSdkxFC1LY86yNFyJjmHbTxn8buhtgp5mY+1uA3ersXg7Q0sH+HrA3wv+Hgh4jCXogaFi03NOha+/NeoujdmmIoRYCfwao8pjX9GtfcBWjBQuffVU/EA8Rj2VAuBk4BMYqfMFxpf+/wPukVKGKTU06PNHMv31CLBRSvnn0P67wB1Sys3h7juUqOTf8epIu3dMVN332bDnNE2jqKiIt99+m+zsbJYtW8bTTz/NggULDl9fVUVXVxf//d//zcqVK5WoKI6Z6q5qfvrhTw+l1QcjIv+iLZJr1kpsgaO/IxwnLCHzl7/EPmfORHY1LK11Pexd38DeTQ3DGvcB0gvjmLMslVknphIdF7m6M1MSXTMCPjW/ITBSN6bkkIZdxZk46lsOJSojClWVUq4WQryM4bl1A0Y54bmh5Zqhnh1aVwBPAH+WUh6Le8RmYI4QogCoxRj1HOnZ9e9QX/4shEgGikLPnXZs2rSJ2bNnU1hYCMDVV1/NSy+9NEBU8kMR0KYp5mqomH7kxebx5AVP8mTxkzy0/SGCehBpErx6imDzHMmqN452P/bu2EnF5Z8n+ZvfIPnGG8cl+/FoSMqM4fTPz2b5pYXUlLazd2MDFduaw5ZCbqjopKGikw+e20fm7Hhmn5RK4dKU0ZVEni6YzMZinZjR2YjzH4TS1L8OvC6ESMSo5ngGcAqQASQDdgwDfAuwF6M65IdSyjFFGEopg0KIm4E3MewlT0gpi4UQvwC2SClXh86dL4QoATSMImAjmqKbatTW1pLTb946OzubjRs3DnGFQjE2zCYzX1/8dc7MOpMff/hjytrLAGhKEPzyahNn75J89T8Q5en3JR0I0PLgH+h+620yf/0rHP1+9EwWJrOJ3IVJ5C5Mwu8NUrm9mb2bGqnZ0za495iEun0d1O3rYO2zZWQUxlH4iRQKl6YQmzyBLsoziGNKqiOlbANeCC0TgpTyNeC1I47d2W9bAreFFoVCcQzMTZzL0599moe2P8STu59EIkEI1iwRbC+U3PiWiWV7B/7695WWUnnlF0i68QaSv/UtTBOQ6mUk2BwW5i7PYO7yDHo7fZR/3MS+zY00VITJmiyhvryT+vJOPvrnfpKyYyg8IZn8Jcmk5I5TipgZyPGZYnMEDGXzGG+ysrI4ePDgof2amhqysrImrT+K4wub2cZtJ93GiqwV/OTDn1Dfa8xYd8QIfnsZLC818Y13zDj710vRNFof+X90v/MOmffeS9TixZPU+8GJjrOz5JwclpyTQ1eLh/1bm9i/tYnmA+EN/K01PbTW9LD51Sqi4+3kLU4if1ESWXMTIp/kcgZxTBPyQogbhBB7hBB3Dt9aMVqWLVvGvn37qKysxO/388wzz7By5crJ7pbiOGNZ+jJeWPkCFxdefPigEGyYb+Lmr+tsWHz0iMS/v5yqq66m6YEH0H1D5/KaLGKTozjxM3l84cfLuO4Xyzntslmk5g2dkqa3w0fJB3W89sddPP5fH/DS/2xj21sHaK3t4XgudDgYx5ql+J/AZcB3pJQPH3HuMxjpVLKBYuB3UspnI9DXiDOVXYpfe+01vvvd76JpGl/72tf4yU9+wp133snJJ5/MypUr2bx5M5dddhnt7e04HA7S09MpLi4e8f2nynsqpgdvVL7BLzb8gu4jIsNP2qfznXfsODuOTmVvKywk895fE3XCCRPVzTHR1eKhckcLFdubqd/fMXQEfz+ccTay5yWQPTeR7HkJ08td+RgZj9T3ZcAsIFdKWdvveBpQDkRx2PNLAj/t7/47VZjKojLeHC/vqYgcDb0N/OTDn7CpYdOA406v5NtrnSzbOshUkslE4le/Qsott2CyTx/PKk+3n6pdrVTtbOFASStB/8jrEMYmO8gsSiBzdjwZs+OIS4macfaY8RCVTsAupXQccfzHwC+BKuDbwGLgHoxptmVSyu1MIZSozPz3VEQWTdf4S8lf+MO2PxDUgwPOnVgp+N7bDuytR4uLraDAGLUsXTpRXY0YwYBGbVkH1btbqd7VQlfL6DKYR8XayJgVR1p+LOmFsaTkxk5ebrIIMR6i4gU8UsqEIx+Ekeb+Uinly6FjtwP3A49LKW8c9cPGESUqM/89FeNDSWsJP1z7Q6q6qgYcj/JKbt+QwqL1DUdfZDKR9LWvkvyd70yrUUt/pJR0Nnk4UNLKgZI2ave2j2oUA0Yux4SMaFJyXaTkukjOjiE5Owa7c/qUohgPUanGsJm4pJTu0LFMjAzCzUB6yMUXIUQC0AqUSymnRvhtCCUqM/89FeOHO+DmN5t/wwv7jo4sOLMmhpvfEJiaj85QaZs1i8z7pp6H2LGgBXUaK7s4WGoITGNlF7p2bIb7mEQ7iRkxJGY4SciIJj7NSXyqkyiXdcpNn405on4QNmKIymXA30PHvhhavyv7KVUoDX0XoHxiFYoZhNPq5O7T7+aMrDO4a91dA4z4H2b3sPWLkl9um0vOmoF5aP3l5VRdfQ1JN9xA8renTlzLsWC2mMicE0/mnHi4GAJ+jcaKTmr3ddBQbkTtj3Qk09Pmo6fNx4HigTHbNocZV3IUsUkOYpOjiEmwEx1vJybBgTPWhjPWNqWm0451pPJp4C2MEcg3gSDwJOACviClfOGI9j2ALqWMHWuHI4kaqcz891RMDPU99dzxwR183PTxUecubcnnupc6kU1Hp/uzFxWRed+9UyIafzzQNJ3Wmh4aK7toqOyksaKLzuajPeXGisVuxhFtwRFtxe60YnOYsTrM2OwWLDYTZqsJs8WEySwQIrSYDM+1omXpwz/gCMalSJcQ4mHgGxwu/CkwUrMskVIG+rXLBg6gpr+mFMfLeyomjqAe5NGdj/LIzkfQj8iMm6JFc//2ecS8NUi6IYuF5G98g+SbVk16DrGJwOcO0Hygm6YD3bTW9NBS00N7gxupT3y8S2qeiyt/tGzU143H9BdSym8JIbYDXwNygd3Azf0FJcRnQuvyY33W8chwqe9/97vf8dhjj2GxWEhJSeGJJ54gLy9vknqrUIDFZOGbS7/JsvRl/PCDH9LkPlx4tdncy9dO2soti87mrKeK0fqPWoJBWh56iJ7//IeM++7FUVQ0Cb2fOOxOK9nzEsmedzg7sBbQ6Why01bfS1tdLx2Nbjqa3HQ0ukftCDAaxqOA2biXExZCrANOBe6SUv5yXB82SoYcqdwdN74Pvzt84Z2RpL5/7733OPXUU3E6nfzxj39kzZo1PPvsyGNM1UhFMZ50eDv42Uc/Y03NmqPOLbHP4s6NWeivH11xUFitJN/yHZK+9rVJr9cyFZBS4u7y093qpavFQ1erl94OHz3tPno7fHi6/bi7/MfsHJBeGMfnf3DSqK8bl5HKCB9swyjqtQcjNb1iBIwk9f0555xzaHv58uU89dRTE95PhSIc8Y54HvzUg/yj9B88sOUBAv1K5u70lXP9yXXcd9KXyX7kVbSWlkPnZCBA8wO/o+edd8m4717sBQWT0f0pgxCC6Dg70XF20gsH/6ErpcTvCeLtDeJzB/D2Bgh4NQI+Db9XQwvoaEENLSjRNYnUJbqUoIMrKfLR/+MqKlJKP/DgeD5jJjLa1PePP/44F1544UR0TaEYMUIIrpt/HUtTl3L7+7dzsPtwklRP0MOt/J2rfnoR17/uwf3m2wOu9ezYQeVll5N62/dI+OIXEapuUFiEENid1lCcy+Sn61f/U9Ocp556ii1btnD77bdPdlcUikFZmLSQ5z73HBfmH/3D59nG1/jupw5i+dUdmOPjB5yTXi+Nv76XA1/5Kv6amonqrmKMqPzN4RjC5jHejDT1/TvvvMOvfvUr3n//fezTNEJZcXwQY4vh/k/ezykZp3DfpvvwaYczGO/v2M91lv/j5w99l0VPfETPfwbaWtybNlG58hJS7/gh8VdeOeUCARUDGXakIoTQxnkJDteH442RpL7ftm0bN910E6tXryY1NXWSeqpQjBwhBFcUXcHfL/o7+bH5A855gh5+UHwvT3wxhZRf34PJNTAVve5203DnXRy86SYCjY0T2GvFaBnJ9JcY50VNwR2BxWLhoYce4jOf+Qzz58/nC1/4AgsXLuTOO+9k9erVANx+++309PRw5ZVXsnTpUlVvRTFtmJs4l2c+9wwXFVx01Ll/7nuBb1qfxf6Ph4k+/fSjzveu/YCKi1fS+fLLqo7JFGXcXYqnMir4cea/p2LqIqXkn/v+yX0b78Ov+weci7HGcM/pv+CkDa00/ua3SLf7qOtd559P+t13YUlMPOqcYnwZyqVYjRIUCsWkIITgyqIr+ftn/06uK3fAuZ5AD997/zYenXWA3BefJ+rko2Mput96i4rPXUz3O+9MVJcVI0CJikKhmFTmJc7j2c89y3l55x117m8lf+OmkruwP3w/qT/8IeKI5JNaWxs1N3+H2h/8AK1z8pxrFIdRoqJQKCadGFsMD5z1AHeccgcW00Cn1O3N27n69Wsp+8xcCv71Io5Fi466vmv1y1SsvISeDz6cqC4rwqBERaFQTAn6giX/csFfSI8emDm3zdvGTW/fxJM975D79N9JufUWsAwUn2BjIwdvvJH6O+9C6+mdyK4r+qFERaFQTCmWpCzh+c89zxlZZww4LpE8tP0hbln7Paxfu5aC557FPkjyyY7nnqPykkvo3bhporqs6Me0ERUhxAVCiL1CiP1CiDuGaPd5IYQUQgzqmaBQKKY+8Y54Hj73Yb699NsIBgY7rq1Zy1WvXEVlmiD/n8+TtGoVHJHGJVBby4Evf5mGX/4K3RP5+iWK8EwLURFCmIH/Ay4EFgDXCCGOquojhHABt2JUppzWvPHGG8ydO5fZs2dz3333HXX+ySefJCUlhaVLl7J06VIee+yxSeilQjF+mISJb5zwDR457xHi7QNTuNT21PLF177IS9Wvknrb98h/+h/Y8vOPukf7U09RcemluD/eNkG9VkyXNC2nAPullBUAQohngEuAkiPa3QPcD4w5Edbiv4xv/exdX94V9pymaXz7298ekPp+5cqVA7IUA1x11VU89NBD49pPhWKyOT3zdJ6/+Hm+v+b77GzZeei4X/dz57o72dmykx+d8iMK/v0vmn//P7T99a/QL/4uUH2A6uuuI/GrXyXl1lswqZRG48q0GKlg1Lc/2G+/hiNq3gshTgRypJSvDnUjIcQqIcQWIcSW5uajy5tOBfqnvrfZbIdS3ysUxyvp0en8+YI/c9Xcq44698+yf/Ll179MY7CdtB/dQd5f/4I1O3tgIylpe+IJKi+7HM+OHRPU6+OT6SIqQyKEMAG/A74/XFsp5Z+klCdLKU9OSUkZ/84dA4Olvq+trT2q3QsvvMCSJUu44oorBiSgVChmIjazjZ8u/ym/PvPXOMwD64Dsbt3NVa9cxcb6jTiXLaPwpX+TcO01R93DX1FB1TXX0vTAA+g+31HnFWNnuohKLZDTbz87dKwPF7AIWCOEqAKWA6tnsrH+4osvpqqqip07d3Leeefx5S9/ebK7pFBMCBfPupinLnqKHFfOgOPtvnZWvb2KJ3c/iXA6Sb/zTnKfeBxLZsbAG+g6rY8+RuXln8ezcyeKyDJdbCqbgTlCiAIMMbkauLbvpJSyE0ju2xdCrAH+S0q5hWNkKJvHeDOS1PdJSUmHtm+44QZ+8IMfTFj/FIrJpi8p5Y8/+DHv17x/6LgudR7Y+gC7WnZxzxn3EH366RSuXk3T/ffT8fw/B9zDX15O1dXXkPT1r5F8883K1hIhpsVIRUoZBG4G3sQoTfyclLJYCPELIcSMS887ktT39fX1h7ZXr16tkkMqjjtibbE8+KkHB3U7fqv6La599Vqqu6oxx8SQcc895Dz6KJa0tIE36Ru1XHY57m3KQywSTAtRAZBSvialLJJSzpJS/ip07E4p5epB2p49llHKZDOS1PcPPvggCxcu5IQTTuDBBx/kySefnNxOKxSTQJ/b8f+d+3/E2mIHnCvvLOfqV65mzcE1AMSsOJPCl1cTd/nlR93HX1FB9bXX0Xjf/SquZYyo1Pcq9b1CMSM42H2Q7733Pfa27z3q3LdO+BY3nXATJmH8ju754APqf3YnwYaGo9pa83LJ/OUvcS5bNu592YP/4QAAG35JREFUnq6o1PcKhWLGk+PK4W8X/W3Q4l8P73iYW9+7lW5/NwAxK1ZQ+PJq4q+84qi2geoDVF//Jep//nO0np5x7/dMQ4mKQqGYMURZorhvxX38YNkPMAvzgHNrDq7h2levpaKzAgCzy2XYWh5/7GgPMaDj6WeMei1r1kxE12cMSlQUCsWMQgjB9Quu59HzHyXRMbAqZFVXFde+ei3vHXjv0LGYM86gcPXLg8a1BBsaqPnGN6m97fsEW1vHve8zASUqCoViRrIsfRnPfPYZFiQNTG/UG+jllvdu4eHtD6NLHQBzTDTpd95J3t/+ijUv96h7db322v9v787Do6rOB45/3ywkQIAQIAoEJESDLAooCoilIAKyFLBaZJO1oogKpShQHq24L1Q0NIIs2oqSiKKWFgLVEoQWoeCPYtnBhECoCAkQ1uzn98dMQiYzk0ySWTLk/TwPT2buOXPve4HkzbnnnveSMnAQ51Z/Tk2eh3aFJhWl1DWraVhTPhzwIUNi7FceLNq9iOnJ07mYe3XexLIa/y9ETJpoV/m4ICuLH+fO5dj4CeSkpno8dn+lSUUpdU0LCQzhpR4vMfvO2XbzLMnHkxm9bjRHs44WbwsIDeW6p5+m1aerCHFwh+Tl7dtJHTqM0/HxFObmejp8v6NJpZoqr/Q9wKpVq2jXrh3t27dn1KhRDvsopa4+VdLRPEtKVgqj1o7inydsH0Vcu317old9QpPfzkBKrbY3ublkLPwjqUOGcmnbNo/H7090nYqTdSr7b/bsGo62B/Y7bSsoKCA2Ntam9H1CQoJN6fvDhw8zfPhwNm7cSMOGDTl16hSRkZEuH1/Xqaia6seLPzIteRr7z9h+DwrC9NunM6H9BERsV+jnHjvGyefncWnrVof7rD9oEJGzniG4At+D/kzXqfgZV0rfL126lKlTp9KwYUOACiUUpWqyonmWwa0H22w3GBZ8t4BZm2dxJd92VX2tli1psXwZzd54ncAI25EOwPm1a0kZMJDMP/0Jk5fn0firO00q1ZArpe8PHTrEoUOH6NGjB926dWP9+vXeDlMpvxUaFMord7/CzC4zi1fZF0k6msS4pHH8ePFHm+0iQoMhQ4hJWkf4Q/bPdSm8dIlTr71O6i9/WaMviWlS8VP5+fkcPnyYTZs2kZCQwCOPPMK5c+d8HZZSfkNEGNd+HIvuXWRXN2z/mf2MWDuC7376zu5zgQ0a0HTe89yQsNLhRH7O4SMcGz+B9KemkZue7rH4qyt/KX3vdWXNeXiaK6Xvo6Ki6Nq1K8HBwURHRxMbG8vhw4e5Q+sVKVUhdzW7i8RBiTyV/BRHzh0p3n4m+wy//vuvmdt1Lg/G2pdzqdO5M9GfruJs4iecfucdCi9csGm/8Pe/c3HTJiLGjaXRo48SGBbm8XOpDnSkUg25Uvp+2LBhbLKWj8jIyODQoUO0bt3aB9Eq5f9a1G/BRwM/oneL3jbb8wvzmfftPF7e9jJ5hfZzJRIURMSY0cQkraPB/ffbtZvcXDKXLuOH/vdxZuXKGjHfokmlGnKl9H3//v1p1KgR7dq1o3fv3rz55ps2D+5SSlVM3eC6vN37bR7r+JhdW+LBRB776jHOZp91+Nmgxo1p9uortPokkdBbbrFrL8jM5KcXXiRl8C84v37DNb0qX28p1tL3SqlSvkr7irn/nGt3F1jzsOYsvGchNzW8yelnTWEhWX9Zw+kFC8g/dcphn9D27WkyfRp1777b7vZlf6C3FCulVAX0vaEvKwasoFndZjbbT1w8wZh1Y9h4bKPTz0pAAOH3DyNmfRKNH38cqV3brk/23r0cf2QyaQ8/zKVt266pkYsmFaWUcqBNRBsSBidw+3W322y/nH+ZacnTWPr90jKTQUCdOjR56kliNqwnfPhwu1piAFd2fsex8RNIGzWai5s3XxPJRZOKUko5EREawdK+S/lV7K/s2uJ2xTFryyyy87PL3EdwZCRNX5hH67+uoV7fvg77XNm1i+OTHyX1/l+StWaNX0/oa1JRSqkyBAcG82y3Z5nbda5dQcqk1CTGrx/PT5d+Knc/ITExRC2Mo9Vnn1G3588c9sk5cID/PTOLI337kfHeEr98hosmFaWUKoeIMOLmEbzX9z27hZJ7M/cycu1I9mTscWlftTu0p+WSJbRKTCDs5z932Cf/5ElOL1jAkV69OfH0M1za/m9MYWGVz8MbNKkopZSLujbtSuKgRFo3sF0TdvrKacavH8+6lHUu76t2p060eG8xrVZ/Rr1+/cDBXWAmL4/zf/0rx8aN44e+/Tgdt5CclOr9LBdNKtVUeaXvf/Ob39CpUyc6depEbGws4eHhxW2BgYHFbaUXTSqlqqZF/RZ8PPBjekb1tNmeU5DDrC2ziPu/uOInSrqidvv2RMW9Q8yG9TQcNQoJDXXYL+/ECTLefZeUgQNJ+cUvOL3wj1zZs7fajWB0nYqTdSrxjzm/ZdAdpi6+x2mbK6XvS1q4cCG7du3i/fffByAsLIyLFy867FtE16koVTUFhQW883/v8MHeD+za7m15Ly/f/TJ1gutUfL/nznHuiy85m5hAXtqxcvsHhodT967u1LnjDmp37EhIbCwS5NkKXGWtU/Gb2l8ich/wDhAILDPGvFaqfQbwayAfOA1MNMakeT1QNyhZ+h4oLn3vLKkkJCQwb948b4aoVI0XGBDIjC4ziAmPYd6382zKuHx97GuOJx1n4T0LaRrWtGL7DQ+n0YTxRIwby6Wt35L1+edc+Mc/MDk5DvsXnDvH+XVJnF+XBICEhhJy442ExLSmVusYgps1JahJE4KaNCEgLAypVYuAkBAkJAQJDHS4z6rwi6QiIoFAPNAXSAd2iMgaY8y+Et12AV2MMZdFZArwBmBfn9oPOCp9v337dod909LSSE1N5Z57ro58srOz6dKlC0FBQcyePZthw4Z5PGalaqqhNw7lhvo3MC15GmeyzxRvP3j2ICPXjuTt3m/TKbJThfcrAQGE3d2DsLt7UHD+POfXr+fC+vVc2v5vKChw+jmTnU32nj1k7yn/xoGQm2+m9ZdfVDi2svjLnMqdwBFjTIoxJhdIBIaW7GCMSTbGXLa+3QZEeTlGn0hMTOTBBx8ksMRvHGlpaezcuZOVK1cyffp0fvjhBx9GqNS1r1NkJxIGJdCmYRub7ZnZmUzaMIm/pfytSvsPrF+fhsOH0/L997npn1to+vJL1OvXj4B69aq03xo7UgGaA8dLvE8HupbRfxKQ5KhBRCYDkwFatmzpdAdlzXl4miul74skJiYSHx9v93mA1q1b06tXL3bt2kVMTIznAlZK0SysGR8O+JA5W+aw8fjVOdncwlzmbJlDyrkUnuj8hN1DwSoqqGFDwh94gPAHHsDk55O9Zw+X/r2DK9/v5sp/dlOQkeHyvpzdFFCl+Ny+Rx8TkTFAF8DhDeDGmCXAErBM1HsxNJeVLH3fvHlzEhMTWblypV2/AwcOcPbsWbp371687ezZs9SpU4eQkBAyMjL417/+xTPPPOPN8JWqseoE12FB7wX8cdcfWfrfpTZtS/+7lNSs1EpP4DsiQUHU7tSJ2p0sl9eMMeSfOk1uyg/kpKSQezSN/FOnyM/IID/jNOZKNoU5OZicHEx2NgEhtdwSR0n+klROAC1KvI+ybrMhIvcCc4GfG2Mcz2r5gZKl7wsKCpg4cWJx6fsuXboU3yacmJjIiBEjbKqc7t+/n0cffZSAgAAKCwuZPXu20wl+pZT7BUgAT932FNENonl+6/PkFuYWt3197GtOrD9B3D1xXF/3ercfW0QIvi6S4OsiqVvil01HjDGQn+/+GPzhlmIRCQIOAX2wJJMdwChjzN4SfToDnwH3GWMOu7JfLX1/7Z+nUr60+/Rupm2cRma2bbmVxrUbE9c7jlua2D97xR/4fel7Y0w+8ASwAdgPrDLG7BWRF0SkaHXfm0AY8KmI/EdE1vgoXKWUAqBjk44OJ/AzrmQwYcME1qeu91FknuMvl78wxqwD1pXa9lyJ1/d6PSillCpH07CmfDjgQ2ZvmU3y8eTi7TkFOTy9+WlSslKY0nGKXz6syxG/GKkopZQ/qxNch7d7v82EDhPs2hbtXsSszeWX0PcXmlSUUsoLAiSAGbfP4MUeLxIUYHuRKOloEpM2TCLjiuu3A1dXmlSUUsqLht04jGX9lhEeEm6z/fuM7xm5diQHzxz0UWTuoUlFKaW87PbrbmfloJV2JfRPXjrJw0kPs+n4Jt8E5gaaVKqp8krfp6Wl0adPH2699VZ69epFenq6D6JUSlVWi3ot+GjgR9zV7C6b7Vfyr/DUxqf4894/++Uz6/3m7i9v+8NDgz26/99+4rwWUEFBAVOnTrUpfT9kyBCbRYwzZ85k7NixjBs3jo0bNzJnzhxWrFjh0ZiVUu5Vr1Y94vvE88aON0g4kFC83WCYv3M+qVmpzO02l+CAYB9GWTE6UqmGSpa+r1WrVnHp+5L27dtXXJm4d+/edu1KKf8QFBDE77r+jt91/R2BYlvgcfXh1Uz5agpZOVk+iq7iNKlUQ45K3584YVuVpmPHjnz++ecAfPHFF1y4cIHMTNtVu0op/zHy5pHE94knLDjMZvv2k9sZs24Maef94/FQmlT81Pz58/nmm2/o3Lkz33zzDc2bN7cpf6+U8j89mvfgo4Ef0TzMtir50fNHGb1uNDtO7vBRZK7TORUnyprz8DRXSt83a9aseKRy8eJFVq9ebfOceqWUf4oJj2HloJVMT57OrlO7irdn5WQx+avJPNftOe6/6X4fRlg2HalUQyVL3+fm5pKYmFhcmbhIRkYGhYWFALz66qtMnDjRF6EqpTwgIjSCZf2WMbi17Q1D+YX5PLf1Od767i0KTaGPoiubJpVqqGTp+7Zt2zJ8+PDi0vdr1ljqZG7atIk2bdoQGxvLTz/9xNy5c30ctVLKnWoF1uKVu1/hyc5P2rV9sOcDZmyaweW8yw4+6Vt+UfreU7T0/bV/nkpdCzYc3cDcf84lp8D2MVHtGrVj4T0LiawT6dV4/L70vVJK1WT9W/Xng/4f0Ci0kc32fZn7GLl2JAfOHPBRZPY0qSillB+4pcktJAxKILZhrM32U5dPMTZpLMnHkp180rs0qSillJ8oejZLz6ieNtuv5F9hWvK0alHaRZOKUkr5kbrBdYnrHceYtmNstheVdnlx24vkFeb5KDpNKkop5XcCAwKZdecsnu32rF1pl08PfcrUr6dyPve8T2LTpKKUUn5qeJvhvNvnXbvSLt/++C0Pr3uY4xeOO/mk52hSqYYmTpxIZGQkHTp0cNh+4MABunfvTkhICPPnz/dydEqp6uSu5nexYsAKu9IuKVkpjF472mZVvjdomRYn0mdv8ej+o177mdO28ePH88QTTzB27FiH7REREcTFxfHll196KjyllB+5seGNfDzwY6YlT2P36d3F28/mnGXShkm80OMFu9X5nqIjlWqoZ8+eREREOG2PjIzkjjvuIDjYf56xoJTyrEa1G7G8/3IGRA+w2Z5XmMecLXN49z/veuXOME0qSil1jQgJDOH1n73OlI5T7NoW7V7ErC2z7Fblu5smFaWUuoaICI93epxXf/aq3RMjk1KTmLRhEplXPPfsJb+ZUxGR+4B3gEBgmTHmtVLtIcCHwO1AJvCQMeZoZY9X1pyHUkpVd4NbD6ZZ3WZMT57O2Zyzxdt3n97N6HWjie8TT0x4jNuP6xcjFREJBOKBAUA7YKSItCvVbRJw1hhzI7AAeN27USqlVPVy23W38fHAj4luEG2z/cTFE4xZN4at/9vq9mP6RVIB7gSOGGNSjDG5QCIwtFSfocCfra8/A/qIiHgxRrcZOXIk3bt35+DBg0RFRbF8+XIWL17M4sWLATh58iRRUVG89dZbvPTSS0RFRXH+vG8WOimlqrcW9VuwYsAKujbtarP9Yt5FHv/6cTYd3+TW4/nL5a/mQMlVPOlAV2d9jDH5IpIFNAIySnYSkcnAZICWLVt6Kt4qSUhIKLP9+uuvJz093UvRKKX8XYOQBiy6dxEvb3uZ1YdXF29vUa8FnSM7u/VY/jJScRtjzBJjTBdjTJcmTZr4OhyllPKK4IBgft/998zsMhNBCA8JJ75PPA1CGrj1OP4yUjkBtCjxPsq6zVGfdBEJAhpgmbBXSimF5c6wce3HEVUvivCQcFrWd//VGn9JKjuAm0QkGkvyGAGMKtVnDTAO+BZ4ENhoKrnSxxiDn07HuMTXpbGVUr7Vp2Ufj+3bLy5/GWPygSeADcB+YJUxZq+IvCAiQ6zdlgONROQIMAOYXZljhYaGkpmZec3+4DXGkJmZSWhoqK9DUUpdg/QZ9aWeUZ+Xl0d6ejrZ2dk+isrzQkNDiYqK0jIvSqlKKesZ9f5y+ctrgoODiY6OLr+jUkopO35x+UsppZR/0KSilFLKbTSpKKWUcpsaPVEvIqeBtEp+vDGlVuvXAHrONYOec81QlXO+wRjjcPV4jU4qVSEiO53d/XCt0nOuGfScawZPnbNe/lJKKeU2mlSUUkq5jSaVylvi6wB8QM+5ZtBzrhk8cs46p6KUUsptdKSilFLKbTSpKKWUchtNKuUQkftE5KCIHBERu8rHIhIiIp9Y27eLSCvvR+leLpzzDBHZJyLfi8g/ROQGX8TpTuWdc4l+D4iIERG/v/3UlXMWkeHWf+u9IrLS2zG6mwv/t1uKSLKI7LL+/x7oizjdRUTeF5FTIrLHSbuISJz17+N7Ebmtygc1xugfJ3+AQOAHoDVQC9gNtCvV53FgsfX1COATX8fthXPuDdSxvp5SE87Z2q8esBnYBnTxddxe+He+CdgFNLS+j/R13F445yXAFOvrdsBRX8ddxXPuCdwG7HHSPhBIAgToBmyv6jF1pFK2O4EjxpgUY0wukAgMLdVnKPBn6+vPgD7i30/4KvecjTHJxpjL1rfbsDyJ05+58u8M8CLwOnAtPBfBlXN+BIg3xpwFMMac8nKM7ubKORugvvV1A+B/XozP7Ywxm4EzZXQZCnxoLLYB4SLStCrH1KRStubA8RLv063bHPYxloeJZQGNvBKdZ7hyziVNwvKbjj8r95ytlwVaGGPWejMwD3Ll3zkWiBWRf4nINhG5z2vReYYr5/w8MEZE0oF1wJPeCc1nKvr9Xi59noqqNBEZA3QBfu7rWDxJRAKAt4DxPg7F24KwXALrhWU0ullEbjHGnPNpVJ41EviTMeYPItIdWCEiHYwxhb4OzF/oSKVsJ4AWJd5HWbc57CMiQViGzJleic4zXDlnROReYC4wxBiT46XYPKW8c64HdAA2ichRLNee1/j5ZL0r/87pwBpjTJ4xJhU4hCXJ+CtXznkSsArAGPMtEIql8OK1yqXv94rQpFK2HcBNIhItIrWwTMSvKdVnDTDO+vpBYKOxzoD5qXLPWUQ6A+9hSSj+fp0dyjlnY0yWMaaxMaaVMaYVlnmkIcaYnY535xdc+b/9JZZRCiLSGMvlsBRvBulmrpzzMaAPgIi0xZJUTns1Su9aA4y13gXWDcgyxvxYlR3q5a8yGGPyReQJYAOWO0feN8bsFZEXgJ3GmDXAcixD5CNYJsRG+C7iqnPxnN8EwoBPrfckHDPGDPFZ0FXk4jlfU1w85w1APxHZBxQATxtj/HYU7uI5/xZYKiK/wTJpP96ff0kUkQQsvxg0ts4T/R4IBjDGLMYybzQQOAJcBiZU+Zh+/PellFKqmtHLX0oppdxGk4pSSim30aSilFLKbTSpKKWUchtNKkoppdxGk4pSSim30aSilFLKbTSpKFUFIrLU+nyV+Q7a/mZtm+mL2JTyBU0qSlVNR+vX7x20dbJ+/Y+ngxCRCBHJF5FMEQn09PGUckaTilJVEwgcBL4rudFaK6uohLjHkwowwBrLemNMgReOp5RDWvtLqSowxtzupKlolHLCGJPhhVAGW7/+zQvHUsopHako5RlFl8W8cekrCOgP5APrPX08pcqiSUWpShKRZc4m6SkxnyIiwSLypIhsF5HzInJSRD4XkdZuCqUH0BDYWvToX2t8p63xufonwU3xqBpML38pVXlljUaKksopYCdwK5bS4gJcB9wP3CEibY0xF6sYh92lLxG5AUu5+stAHSAby6OuSwoD6gI5wDkguYpxKKUjFaUqw3qHVQfr292l2kKAm61vfw9cwjKaKPoh/issl6qigIfcEI5dUjHGpBljrufqQ6jijTHXl/wDfGZtW2bdtsQNsagaTpOKUpXTBstTAXOA/aXaOnD1KsAhoJcxZqu56jOuJoDYqgRhvYR2M5BijCkdB8At1q//ddDWsYw2pSpFk4pSlVP0A3mfMSbfSVseMNIYk+vg80VzH1V9St4vrF/t7vqyPjK3jfXtf0u1BQFtHbUpVRWaVJSqnKLEsdtBW9F8ytfGmKNOPl80Se+s3VVl3UrcFsuIqQDYV6rtZiDE+npPFWNQqpgmFaUqp6zV8kVtfy3j80WXpSr9A11E6gE9gQvAN2Uc47AxJrtUW1FSTDPGnK9sDEqVpklFqcpxOFIREcFypxeUWmVfok9LIALLpS9HIx1X9QdqAX93comtKA6dT1Feo0lFqQoSkUjgeuvb0kkhGmgAFOJ8FFI0kkk1xlwosd/mInJFRFaWOl536/blpfZT3ir6opGKo7pkZSUcpSpNk4pSFVeUFI6VXGxoVTQCOGSMuezk8w7XtxhjTgBLgOEichOAiEQDfwG2Ao8V9RWRACz1vgqBdU6Oo3d+Ka/TpKJUxbkySV9WeZaiPo4+/yqQC8wRkQbAWiATeMAYk1ei351AJLDDGHOq9E5EJIKrBS1L3/nViKsjLZ2kV26lK+qVqjhXVtK7klTs+hhjTorIIuBJLLcDNwa6GmPOlerq6qWvi0BqqbboEq9/KiNOpSpMRypKVVxZI40yk4qI1OfqD3Vnk/TzsZSx7wIMNcaUTgrgelLZY4wpvRamsMTrd0RknJN9KFVhmlSUqgBrCZaiBYWl7/wKB1pa3zobqdyKpf7XOWNMmpM+c7F8bwYCZxzEEIVltJRujHF2nLLmU77HcrtzHpbkdIOTfShVYWL/S4xSyhkRuQP4N5bijA0djAKquv/pwAJgGvAisNYYM6pUn8eARcB7xpjH7PeilO/oSEWpiulv/fqtBxLKEOAPwMvGmDhgIfCQiLQt1VUfyKWqLU0qSrlIRDoBU61v/+Lmfd8GrMRSOfhZ6+YFWErXP1+q+2ZgHvAPd8aglDvo5S+lyiEik7FMntezbtoF3OWg9Ell9x8FbAeOAb1L7ldE3gBmArcaY/T2X1XtaVJRqhzWleyjsBR//AJ4TetlKeWYJhWllFJuo3MqSiml3EaTilJKKbfRpKKUUsptNKkopZRyG00qSiml3EaTilJKKbfRpKKUUspt/h/Kat8IZ4MFcQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "jhx = np.linspace(0, np.pi, 100) / np.pi\n", "\n", "#lmbdas = [0.1, 0.2, 0.4, 0.8, 0.9, 1.0, 1.2, 1.4, 1.8, 2.0]\n", "lmbdas = [0.1, 0.2, 0.5, 0.6, 0.75, 0.9, 1.1]\n", "\n", "for lmbda in lmbdas:\n", " plt.plot(jhx, ampFTBS(jhx, lmbda), label='%g'%lmbda, lw=4)\n", " \n", "plt.legend()\n", "plt.xlabel('$j h_x / \\pi$', fontsize=24)\n", "plt.ylabel('$|s(j)|$', fontsize=24)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Dispersion\n" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false, "jupyter": { "outputs_hidden": false } }, "outputs": [], "source": [ "def alphaoverj(jhx, lmbda):\n", " return np.arctan((lmbda * np.sin(np.pi*jhx)) / (1 - lmbda + lmbda*np.cos(np.pi*jhx))) / (lmbda * np.pi * jhx)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false, "jupyter": { "outputs_hidden": false } }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/andreas/src/env-3.7/lib/python3.7/site-packages/ipykernel_launcher.py:2: RuntimeWarning: invalid value encountered in true_divide\n", " \n" ] }, { "data": { "text/plain": [ "Text(0.5, 0, '$j h_x / \\\\pi$')" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEYCAYAAABfgk2GAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dd3hUVfrA8e+ZmfRGQhIICSGhh15CEZUuCqugqyK47mJvWNeGbS2ru64ruCpIsYEF9GcFFVGKgCCCQTpIMbSEFlpILzPn98edhBAChMzcTDLzfp7nPnNn7p257w3hzZn3nnOu0lojhBDC+1k8HYAQQojaIQlfCCF8hCR8IYTwEZLwhRDCR0jCF0IIH2HzdABnEh0drZOSkjwdhhBC1CurV68+rLWOqWpbnU34SUlJpKWleToMIYSoV5RSu8+0TUo6QgjhIyThCyGEj5CEL4QQPqLO1vCFEKKmSkpKyMjIoLCw0NOhmCYwMJCEhAT8/Pyq/R5J+EIIr5ORkUFYWBhJSUkopTwdjttprTly5AgZGRkkJydX+31S0hFCeJ3CwkIaNmzolckeQClFw4YNz/sbjCR8IYRX8tZkX6Ym5yclHSHOg9aaYkcx+SX55JXkkVeSR0FpAXkleeSXnnwtvySfFg1aMDBxoKdDFqKcJHzhlRzaQWFpIfml+RSUFlBYWkhBacEpS35J/qnPnfsWlBjr+aX55Jec+lhQUkCpLq1WDFGBUZLwfdy8efO4//77sdvt3HrrrYwbN+6U7RMmTODtt9/GZrMRExPDu+++S7NmzQCwWq107NgRgMTERObMmeNyPG5J+Eqpd4HLgUNa6w5VbFfAa8AwIB+4UWv9mzuOLeqPUkcpxfZiCu2FxmNpIUX2oiqfF5U6H+1FFJYWUmgvNB7PsF5QWkChvbA8uRfZi84rNquyEmwLJsgWRLCf8RhkCyIqMIoEvwSCbEGE+IUQbAsm2C/4lMcQv5BT1kP8Qpi0dhJz0+ea9JMU9YHdbmfs2LHMnz+fhIQEevTowfDhw2nXrl35Pl27diUtLY3g4GAmT57Mo48+yieffAJAUFAQa9eudWtM7mrhTwcmAu+fYftQoJVz6QVMdj4KN9NaY9d2Sh2llDpKKXGUnPZY4iihxF5S/lqxo7j8eYmjhGJ78Snrxfbi8n2KHcWnvFZ5vcheVP5a2XrZY3VbxlUJsAYQYA0g0BZIkC3o5Lo1iLDgsPIEHWgNPLluq3q9LLGXL35B+Fv83VrzDbAGuHS+ov5btWoVLVu2pHnz5gCMGjWK2bNnn5LwBwwYUL7eu3dvPvzwQ1NjckvC11ovVUolnWWXEcD72rif4i9KqQZKqTit9X53HL+iE3nH+OiH/6AAlEKjUMq4wKGVMv5TK4VSFpRFGYuyoCwWsFqwONeVRWGxWlFWCxaLFe18v8O5aBQOZcEBaIwk69CO05ay10sdpac8t2u7se5wUKpLsTuc+znXy5J25ceq1isndTMoFP5Wf/wt/vhZ/crX/a3+5esB1gBC/ULxt/qXJ+iy7YHWwPLnZcm67PXy7baA8v0CbYHlrwVYA7Co+tW/wGqxYnfYPR2GAJ77ehOb951w62e2axLOM1e0P+s+mZmZNG3atPx5QkICK1euPOP+77zzDkOHDi1/XlhYSGpqKjabjXHjxnHllVe6HHdt1fDjgb0Vnmc4Xzsl4SulbgduB6NmVRNZRzN588S3NYvSTaxaY0GhMLpBWVFYUCcflcKqLFiVBZuyYrVYsSorNosVi7Jis/hhtdiwWmwEWP2xWv2xWfyw2QKxlj+34Wfxc77PVr74WfxOW6/46Gf1w08Zj+Wvlb1u8StP6GXr/lbjuU3ZvL7XgzvZlA27loQvqufDDz8kLS2NJUuWlL+2e/du4uPjSU9PZ+DAgXTs2JEWLVq4dJw6ddFWaz0NmAaQmppao7urx8c2540O/wSHBg1aO4zPdjjQ2oHWGu3QaIcd7QCHoxSH3YF22HHYHeCwo+2lOBx2dKkD7ShFl9rBXorDXgKlpWAvRdtLobQUZS9B2Yuw2EuxOIqx2ktQjmKs9mL8KMWfEgJUCQGUGOvOJVAVE0AxgRTjr84vMWhlAf8Q8A9DBYQa6wFh4B9mPAaEQUAoBIQb64EREBAM/hEQGG48D4wA/1CQJG4Kq8WKXdvRWssfSg87V0vcLPHx8ezde7Kdm5GRQXx8/Gn7LViwgBdffJElS5YQEBBwyvsBmjdvTv/+/VmzZk29SfiZQNMKzxOcr7ldYEAw/bu7/tXHHYpLHeQXl5JXbCevqJTcolKOFBqPuYWl5BSVklNYQk5+IUX5eRQV5FJUkEdxYS6lRXnYC/OgJJ9gighWhQRTRAiFBKtCQkoKCS0oItJWRISliHDLEUJVJsG6gEBHHgH2PCznqiErq/MPQAMIagBBkVUsURAcdfIxuKHxx0KS2FlZlRUAu7ZjU3WqXSVqSY8ePdi+fTs7d+4kPj6ejz/+mJkzZ56yz5o1a7jjjjuYN28esbGx5a8fO3aM4OBgAgICOHz4MMuXL+fRRx91Oaba+k2cA9yjlPoY42Jtthn1+7rG32bB3+ZPg+Caf0ap3cGJwlKO5RdzPL+YY3klzvUS0vOLOZpbzJG8Yo7mFRmPucXkFJUCmgBKCCefMJVPKAU0CSwiPqiEJgFFxPoVEW0rINKST7jKJ1TnEJR3DNux3aiCY1B4HJzfjk5jsRmJv2wJiYaQGOcSDSGxxnporLH4h9T8B1BP2SzGfy27tmOrW1+kRS2x2WxMnDiRSy+9FLvdzs0330z79u35xz/+QWpqKsOHD+eRRx4hNzeXa6+9FjjZ/XLLli3ccccdWCwWHA4H48aNO+Vib41jcvkTAKXULKA/EK2UygCeAfwAtNZTgLkYXTJ3YHTLvMkdx/UFNquFqBB/okL8q/2ewhI7h3OLOJxbzOGcIrJyi8jKKeJQTiF7TxSRllNEVnYhh3KKKHWcWjmzWRSNwgNpEuNP83AHySGFJAYW0sS/gFhbHlEqh8Di45B/xFjyDsOBDZCXBYXZVQfkH+pM/o2MJayxc4mr8BhnfNvwEmWtervDDlYPByM8ZtiwYQwbNuyU155//vny9QULFlT5vj59+rBhwwa3x+OuXjqjz7FdA2PdcSxxboF+VhIig0mIPPtXC4dDcySvmIMnCtmfXciB7ALnYyGZxwtYsa+EL7JLKbFbgVDn0ogGwX4kRAaR0CCYxLhgmkYF0zQyiGYN/Ejwz8Ov4LDxhyD3oLHkZTnXD8HBjbBjIRTnnB6QfyiEN3Eu8c6lCUQ0hYgEYwkINeNH5nZWi5HlpWumqEvku6YPs1gUMWEBxIQF0CE+osp9HA7N4dwiMo4XkHmsgMzjBWQcy2fv0QK2Hcph0dZDFJeeLP1YFMRHBpHUMIxmDRuR1DCEpMQQmseE0DQqGD+rs3tlUS7kHIDcA3BiP+Tsq/C4D/740dhWuawUFOlM/onQIBEimxmPDZoZ6wFhZv24zkt5DV+6Zoo6RBK+OCuLRREbHkhseCDdEiNP2+5waA7lFLHnaD67j+Q5H431OWv3caLwZAvXalEkRgXTIiaEFjGhNI8JoWVsO1q27EVEcBVzettLIWc/nMiE7Aw4vsd4zN4LR9MhfTGU5J36nuCGEJkEkckQlQxRzU8uITG1drG5Yg1fiLpCEr5wicWiaBwRSOOIQHomR522/VheMTuP5LEzK4+dh/NIP5xLelYeS7cfPuWbQXRoAK1iQ2ndKJSWjcJo41wiGjSFBk1P+1wAtIb8o3B8FxzbDcd2nVwyVsGmL079huAfBg1bQMOWxhLdylgatnT7heWyFr5ZA+GEqAlJ+MJUkSH+RIb4n/btwO7QZBzL54+sXHYcymX7wVy2H8rl898yyS06mSQbhwfSpnEYbePCSGkcTtu4MFrEhBqlIaUgpKGxxHc//eClxca3gmM74cgfcPQP4zHjV9j4OVDhgnVEU4huDTFtnEuK8RjUoEbnXVbDlxa+qEsk4QuPsFoUzRqG0KxhCAPbNip/XWvNvuxCth3MYduBHH53Liv+OEKx3Wit+1sttGoUSru4cNo3Cad9fATt4sIJCaj062zzh+iWxtLqklO3lRQafwAOb3cuWyFrK+z+GUoLTu4X1gRiU6BRO4htD43aG38IbAGcjdTwRV0kCV/UKUop4hsEEd8giAFtTg5EKbE7SM/KY8v+E2zZf4LN+0+w8PdDfLo6w/k+aB4dQsf4CDomNKBTQgTtm4QT7H+GX3G/QCN5N6o0CtPhgOw9cOh3yNpiPB7aDCunQdkMnBab8W2gcUdo3AniOhnrQSe/xZTV8KWXjm871/TIU6ZMYdKkSVitVkJDQ5k2bRrt2rVj1apV3H777YDRCHr22We56qqrXI5HGT0m657U1FSdlpbm6TBEHaa15sCJQjZlnmDjvmw2ZmazITObgyeMxGxR0LpRGF2aNqBz0wZ0adqA1o3CsFpqcOHWXmpcKD64AQ5shIObjPEHOftO7tOgGTTpAnFd+MEfHvr9Pb4Y/gWtIlu56YxFdW3ZsoWUlBSPxmC322nduvUp0yPPmjXrlAFUJ06cIDzcGH8yZ84c3nzzTebNm0d+fj7+/v7YbDb2799P586d2bdvHzbbqQ2Yqs5TKbVaa51aVUzSwhf1llKKuIgg4iKCGNzuZFno0IlC1mdksz7jOOsyspm36QAf/2rMaRLsb6VzQgO6Jjage7NIuiVGElmdQW1WG8S0NpYOV598Pe8w7F8HB9bDvrXG+ubZWIODoFEM9o9GQlyqcY2haQ9o1NEoNQmvV53pkcuSPUBeXl75vEvBwSfH0BQWFrptPiZJ+MLrxIYHMrhdYPkfAa01e47ms2bPcdbsOcZve44zdWk6duco4xYxIfRIiiI1KYoeSZEkRgVX/z9YSDS0HGQsZfKPYtv0EWyehj0qCXYtgw3/Z2yzBhjfApr2hKa9oWkvCI1x49mL03w3zvg25k6NO8LQl866S3WnR540aRITJkyguLiYRYsWlb++cuVKbr75Znbv3s0HH3xwWuu+JiThC6+n1MkLxFd2NWYgLCi2sy7jOKt3H2P17mPM3bC//FtAo/AAeiY3pGdyFL2To2gZG3p+LazgKKxNusJmKL3kWYjpDNmZkJlm9BDauwpWToWf3zD2b9gSEntDYh9IutAoDcnkdD5j7NixjB07lpkzZ/LCCy8wY8YMAHr16sWmTZvYsmULY8aMYejQoQQGBrp0LEn4wicF+Vvp3bwhvZs3BIwBZNsP5bJq11F+3XmUlTuP8PU6oz4fHepPr+YNuaB5Qy5sGU1Sw3N/Azitl05EvLG0G2E8Ly0ySkB7VsCeX2DLN7DGebej8AQj8Sf3NZYGNbs3hHA6R0vcLNWdHrnMqFGjuOuuu057PSUlhdDQUDZu3EhqapWl+WqThC8ExgCyNo3DaNM4jL/2blZeBlqZfpQV6UdY8ccRvl1vTPAaFxFInxbRXNTK+AMQG3Z6q+ucI21tAZDYy1jA6B2U9TvsXm6UgHYshPXGvU2JTILkftBigPEYfPoAN1H3VGd65O3bt9OqlXFR/9tvvy1f37lzJ02bNsVms7F7925+//13kpKSXI5JEr4QVahYBhrZoylaa3YdyWf5jsP8/MdhFv5+kM9/M7qEtm0cxsWtounbOoYeSVEE+lnPf6StxWL09W/UDnreZowiPrQFdi41lk1fwm8zAAVNujqvGwyG+FTjgrKoc6ozPfLEiRNZsGABfn5+REZGlpdzli1bxksvvYSfnx8Wi4U333yT6Ohol2OSbplC1IDDodm07wQ/7chi2fbDpO06RrHdQaCfhT4tommdeISP9j7C5MGTuSj+ItcPaC+Ffb8Zk8r9sciYOkI7jJvRtBgErYYYg8tCXE8K3qAudMusDdItU4haYLEoOiZE0DEhgrv7tyS/uJRf0o+wZGsWi7dlsXjXbkKS4fEv1jK0RUMGtImlR1IkNmsNb8ZutTl79vSE/o9BwTFIXwLb58P2H4x5g1DG9jZDoc0wY0SwEBVIwhfCDYL9bQxs26h8moiFfzTggWXQIMTGe8t3Mm1pOhFBfgxoE8Pgdo3o1zqGsMAqZgitrqBIaH+lsTgccGAdbPsets6FBc8aS8NWkHI5pFwBTbpJzx8hCV8IMyRGGTdquXdgc/r8ZSA/bctiwZZDLPr9IF+t3Ye/1UKflg0Z0q4xl7RrREzY2efmOSuLxajrN+kK/ccZXUC3zoXfv4Hlr8OyV437B7QbDu2uhIRUSf4+ShK+ECaoOFtmaICNoR3jGNoxDrtDs3r3MeZvPsD3mw7yxJcbePKrDfRoFsVlHRoztGNj4iKCXDt4RLxx4bfnbcb00dvmwebZsGoarJhodPNs/2foeA006iDJ34dIwhfCBGX3tK3cS8dqUfRMjqJnchRPDEth68EcvttwgO83HeD5bzbz/Deb6d4skmEd4/hTxzgaR7g20IbgKOhyvbEUZsPvc42poX9+A5b/z5gGutNI6Hjtme87ILyGJHwhTFCd+fCVUrRtHE7bxuE8eElr0rNy+W7jAb5Zv59/frOZF77dTI+kKK7oFMewjnE0DHWh7ANGj54uo40l7whs/hLWfwoLnzOWpIuhy1+M0o+bbwgj6oYadhkQQpxNTebDbx4TytgBLfnu/otZ9FA/HhzcmqN5xTw9exM9/7WQG99bxZdrMsgrcsOUyyENocetcMv3cN9aGPCkcfvIr+6EV1rD7Htg76/GeABRY/PmzaNNmza0bNmSl146fcTvnj17GDBgAF27dqVTp07MnTsXgJKSEsaMGUPHjh1JSUnh3//+t1vikRa+ECZw9Z62zWNCuW9QK+4d2JLfD+QwZ90+5qzdx4OfrCPIbyOXdWjMVV3jubBldM2me64oKhn6PQp9HzGmeVj7IWz8AtZ8ALHtoNsY6Dyqxnf/8lV2u52xY8eeMj3y8OHDT5kt84UXXmDkyJHcddddbN68mWHDhrFr1y4+/fRTioqK2LBhA/n5+bRr147Ro0e7PNpWEr4QJnDXPW2VUqTEhZMSF84jQ9qQtvsYX67J5Nv1+/hyTSaNwgO4qmsC13RPoGVsqGtBKwXNLjCWS/9t9O1fPQPmPWaUfDpcDT1uMXoDiXOqzvTISilOnDgBQHZ2Nk2aNCl/PS8vj9LSUgoKCvD39z9lKuWakoQvhAnMuKetpcIF32euaMei3w/x+eoM3vopnSlL/qBrYgNGpjblis5NCK18u8fzFRgO3W80ln1rIO1d2PCZ0epv2gt63m5MBGd1YSxBLfnPqv/w+9Hf3fqZbaPa8ljPx866T3WmR3722WcZMmQIb7zxBnl5eSxYsACAa665htmzZxMXF0d+fj6vvvoqUVGuz6EkNXwhTFDWS8ese9oG+lkZ1jGOd27swYrHB/LEsLbkFpby+Bcb6PHCAh75dB2/7TmGW6ZOadIVhr8Bf99itPxzD8Hnt8D/OsJPE4yun6JGZs2axY033khGRgZz587lr3/9Kw6Hg1WrVmG1Wtm3bx87d+5k/PjxpKenu3w8aeELYYKyFn5t3NM2NiyQ2/u24LaLm7Nm73E+WbWXr9fv49PVGbRtHMbonolc1S2ecFdG9oJRw7/gbuh1J+yYD7+8aZR6lv7X6N3T5x5jZs865lwtcbNUZ3rkd955h3nz5gFwwQUXUFhYyOHDh5k5cyaXXXYZfn5+xMbGcuGFF5KWllZeHqopaeELYYKa9NJxlVKKbomR/OeaTqx6cjD/uqojflYLz8zZRK8XFzLu8/VszMx2/UAWC7S+FP42G+5cbozeXT0dXu8Gn93i/rtL1VMVp0cuLi7m448/Zvjw4afsk5iYyMKFCwFjIrTCwkJiYmJITEwsv/tVXl4ev/zyC23btnU5Jkn4QpigPOG7sYZ/PkIDbFzfK5Gv772Ir++5iOGdm/DV2kwuf2MZf35zObPXZlJc6nD9QI07wFWT4YH1Rut/2zyYchHMvA4yfHu224rTI6ekpDBy5Mjy6ZHnzJkDwPjx43nrrbfo3Lkzo0ePZvr06SilGDt2LLm5ubRv354ePXpw00030alTJ5djkumRhTBJl/e7cHOHm7mv232eDgWA7IISPl+dwQe/7Gbn4TyiQ/25vlczbuidWOVNXGqk4Disegt+mWTM6Nl8APR//OSNXmqJTI9c9fTI0sIXwiRWZfVYC78qEUF+3HxRMgv/3o/3b+5J54QGvLFoOxe+tIgHP1nLhgw3lHuCGkC/R+CBjXDJ80Z5590h8MGffb7FXxfIRVshTGK1WF3uh28Gi0XRt3UMfVvHsOtwHtN/3sWnaXv5ck0mvZKjuOWiZAanNMLiyoCugFC48H5jNO+vb8Oy/8Hbg4x5+gc+bdzZS9Q6aeELYRKbstWpFn5VkqJDeHZ4e1Y8MYin/pRCxrECbv9gNYMnLGHmyj0UlrgYv3+Ikfgf2AADnzLu1zu5D3xxOxzb7Z6TOIO6Wq52l5qcnyR8IUxSV1v4VQkP9OPWi5uz5JH+vDG6KyEBNp74cgMX/WcRExdtJ7ugxLUDBIQaUzfcvw4uvM+YrnliD/jhaaPu72aBgYEcOXLEa5O+1pojR44QGHh+117koq0QJun/SX8GJA7gmQue8XQo501rzYr0I0xZks7SbVmEBtj4S69EbrkomdhwN1zgzc6ARS/CullG3b//45B6s9tG7paUlJCRkUFhYaFbPq8uCgwMJCEhAT+/U39mck9bITzAarHWaj98d1JK0adFNH1aRLNpXzaTF//BWz+l897PuxjVoyl39GtBfAMXbtQSkWB05+x9F/zwFHz3qDF9w2X/hhYDXY7fz8+P5ORklz/H20hJRwiT1IcafnW0bxLBxOu7seih/vy5azyzVu2h38s/8uhn69hzJN+1D4/rZAzgGjUTSgvhg6tg1vWm1/d9lSR8IUxSn2r41ZEUHcJLV3diySMDuKF3M75au48B4xfz8Kfr2Hk4r+YfrBS0/ROMXQWDnoH0H2FST1jyMpR4b0nGEyThC2GSutYP312aNAji2eHtWfboAMZckMTX6/YxeMISHv7UxRa/LQAu/jvc8yu0vgx+fNHo0ZO+2G2x+zpJ+EKYxGax1dsafnXEhgfyjyva8dNjRuKfs24fA8cvZtzn68k8XlDzD45IgJEz4K9fgnbA+yPgizsg77D7gvdRbkn4SqnLlFJblVI7lFLjqth+o1IqSym11rnc6o7jClGXWZW1VmbL9LTYMGfif3QA1/dK5IvfMhnw38U8O2cTWTlFNf/gFgPh7hVw8cPGjdcn9jDuwVtHexbWBy4nfKWUFZgEDAXaAaOVUlUNo/tEa93Fubzt6nGFqOvqcy+dmmgUHsjzIzqw6OF+XNU1ng9+2U3fl3/kle+3cqKwhv34/YJg0NNw508Q1Ry+uBVmjYLsTPcG7yPc0cLvCezQWqdrrYuBj4ERbvhcIeo1b+mlc74SIoP5zzWdmP9gXwalxDLxxx30fflH3lqaXvORu7EpcMsPcOm/IH0JvNkb1s6U1v55ckfCjwf2Vnie4XytsquVUuuVUp8ppZpWsR2l1O1KqTSlVFpWVpYbQhPCc7y9hn8uzWNCmXh9N7659yI6xkfw4twtDBq/hC/XZOBw1CBRW6xwwVi4+2do1AG+ugtmjYacg+4P3kvV1kXbr4EkrXUnYD4wo6qdtNbTtNapWuvUmJiYWgpNCHNYLb5Rwz+XDvERfHBLLz66tReRIX48+Mk6Ln9jGcu21/AibFRzuPFbZ2v/R6O1v+Vr9wbtpdyR8DOBii32BOdr5bTWR7TWZVdv3ga6u+G4QtRpVuVbNfxzubBlNHPGXsRro7pworCEG95ZyU3vrWLbwZzz/zCLxWjt3/ETNEiET26A2WOhKNf9gXsRdyT8X4FWSqlkpZQ/MAqYU3EHpVRchafDgS1uOK4QdZrV4p398F1hsShGdIlnwd/78fjQtqTtPsZl/1vKU19t4EhuDXr0xLSGW+bDxQ/Bmo+Mu21lrnZ/4F7C5YSvtS4F7gG+x0jk/6e13qSUel4pVXYDx/uUUpuUUuuA+4AbXT2uEHWdTdm8aqStOwX6WbmjX4vyUbuzVu2l/yuLefun9PO/9aLNHwb9A26aC45SeGcILH8dHG64haOXkdkyhTDJ/YvuJyM3g8+Hf+7pUOq87QdzeHHuFhZvzSI5OoSnL09hYNtG5/9BBcdgzr1GTb/FILhqKoT61vVAucWhEB7ga/3wXdGqURjTb+rJezf1QCm4eXoaN763ij+yzrMmHxQJIz+AP00wbrYy9WLYvcKcoOshSfhCmMRX++G7YkCbWObd35en/pTC6l1Gff/f320hr+g8SmNKQY9b4NYFxsCt6X+C5a9Jn30k4QthGm+bLbO2+Nss3HpxcxY93J8ru8QzdUk6A8cvZvbazPO7g1VcJ7h9MaRcDvP/YfTkKTxhVtj1giR8IUzirbNl1paYsAD+e21nvri7D7Fhgdz/8Vquf2sl28+nG2dgBFw7Ay79N2z9Dt4aAId8t5OgJHwhTOLrI23dpVtiJF+NvZAXruzA5v0nGPraT/z7uy3kF1fz25NScMHdcOM3Rgv/rUGw6Stzg66jJOELYRJfmS2zNlgtiht6N2PRQ8bEbFOXpDN4/BK+33Sg+mWeZn3gjqXQuAN8OgYW/tPnum5KwhfCJDLwyv0ahhplns/uvICwQD/u+GA1t85II+NYNW+8Eh4HY76Gbn+Dn16Bj0dDYba5QdchkvCFMIlMrWCe1KQovrnvIp4Y1paf/zjCJROW8tbSdErt1Wix2wLgitdh2CuwYwG8fQkcTTc/6DpAEr4QJrFZpFummfysFm7v24L5f+/LhS0b8uLcLVwxcTnr9h4/95uVgp63GXfVyjtk1PV3LTc/aA+ThC+ESaxKumXWhoTIYN76WypTbujG0bwirnpzOc99vYnc6vTdT+4Lty6EkGjjVoq/fWB+wB4kCV8Ik0gNv/YopbisQxzz/96PG3o3Y/rPuxgyYQkLt1RjrvyGLYwJ2JIugjn3wILnvPZiriR8IUxiUzYc2oFDe2fyqIvCA/14fkQHPruzD6GBNm6Zkca9s9Zw+FwzcQY1gL98Cp82D5AAAB4wSURBVN1vgmUT4LOboMSFG7HXUZLwhTCJ1WIFkFa+B3RvFsk3917M3y9pzfcbDzB4whI+X51x9i6cVj+4/FUY8gJsng0zroC8I7UXdC2QhC+ESazKmfClp45H+Nss3DeoFXPvv4gWMaE89Ok6xrz369m7cCoFfe6Fke/DgQ3wjnf14JGEL4RJbBYbIC18T2sZG8and1zAc8Pbs3rXUYa8upQZP+86+3112w2Hv80xplt++xLI8I6p2iXhC2GSsha+9NTxPItFMaZPEt8/2JfUpCiembOJ66atIP1s0y8n9jIu5gaEGuWdbd/XXsAmkYQvhEmkhl/3JEQGM+OmHrxybWe2Hsjhstd+YsqSP848YCu6pZH0o1vDrNGw5sPaDdjNJOELYRKp4ddNSimu6Z7Agr/3o3/rGF767neunvzzmW+mHhprTLzWvJ9xo/Slr9TbufUl4QthEqnh122x4YFM/Wt33hjdlb3HCrj89WVMXLSdkqpa+wFhMPoT6DQKFv0T5o2rl331JeELYRKp4dd9Simu6NyE+Q/25ZL2jXjlh21c9eZytuyv4kYpNn+4cjJccA+snAJf3AalxbUftAsk4QthEqnh1x8NQwOYdH03Jv+lGweyCxk+cRmvLaiitW+xGP30Bz8HGz+DWaOgOM8zQdeAJHwhTGJTzpKO1PDrjaEd4/jhwX5c1iGOVxdsY8TE5WzeV6m1rxRc9AAMnwjpP8L7V0L+Uc8EfJ4k4QthkrIWvtwEpX6JCvHnjdFdmXJDdw7lnKW13+2vxgCt/WuNG6XnHPBMwOdBEr4QJpFeOvXbZR0aM//BfgzreJbWfsoVxhw8x/fAO0Pg6E7PBFtNkvCFMIn00qn/IkP8eX10V6b+9Syt/eb9YcwcKDoB714GBzd7KtxzkoQvhEmkl473uLT9qa39KydV6skT3x1u+s5Ynz4MMlZ7JtBzkIQvhEnKa/iS8L1CWWt/yg3dOXjCaO2/vrBCaz82BW6eBwHh8P5w2PmTZwOugiR8IUxSXsOXko5XKavtD+0Qx4T5lVr7Uclw8/cQkQAfXg1b53k22Eok4QthkvIavly09TpVtfbLa/vhcXDjXGjUDj75C2z4zNPhlpOEL4RJymv40i3Ta1Vs7Zf15Nm0LxtCGhrTKzftBZ/fCmnveTpUQBK+EKYpH2krLXyvVrG1fyiniBETlzNh/jaKbaFww+fQ6hL45gH4+Q1PhyoJXwizlI+0lRq+TzBa+325onMTXl+4neETl7HhYDFc9xG0uxJ+eAoWvejRmTYl4QthEinp+J7IEH9eva4Lb/8tlWP5xVz55nJeXpBO4Yi3oOsNsPRlj860afPIUYXwAVLS8V2D2zWiR3IUL367mTcX/8H3mw7w8tXP0T0gHH55E4py4IrXwVq7KVha+EKYREo6vi0iyI+Xr+nM+zf3pLDEwTVTV/Js0V8ovugxWPsRfHYjlBbVakyS8IUwibTwBUDf1jF8/2Bf/ta7GdNX7GZAWm/+6PYkbPkaZl4HRWe5r66bScIXwiRSwxdlQgNsPDeiA5/eeQEBfhYG/dyeWU3GoXcugQ9qb3plSfhCmERa+KKyHklRzL3vYu4d2JKnd3XmEfUQ9n1r0bU0vbJbEr5S6jKl1Fal1A6l1LgqtgcopT5xbl+plEpyx3GFqMv8LH6A1PDFqQL9rDw0pA3f3HcROxoO4K+Fj1CUlU7pW0PgaLqpx3Y54SulrMAkYCjQDhitlGpXabdbgGNa65bAq8B/XD2uEHWdzJYpzqZt43A+v6sPQ/40kjH2p8nJPkr+lMGUZK4z7ZjuaOH3BHZordO11sXAx8CISvuMAGY41z8DBimllBuOLUSdJfe0FeditShuvDCZV/9+C682fZ2jRfDb+5fz0YLxphzPHQk/Hthb4XmG87Uq99FalwLZQMPKH6SUul0plaaUSsvKynJDaEJ4jtzTVlRXgxBNp74ObmuXwq1NGvDRzunYS0rcfpw6NfBKaz0NmAaQmprqufHHQriB3NNWnMvO7J18svUTZu+YTW5JLq0iW/FIq/vpFZqK1c/P7cdzR8LPBJpWeJ7gfK2qfTKUUjYgAjjihmMLUWdZlAWFkha+OIXdYeenzJ+YuWUmK/avwGaxMaTZEEa1HUWXmC6YWe12R8L/FWillErGSOyjgOsr7TMHGAOsAK4BFmntwRmEhKglVotVavgCgJziHL7Y/gUf//4xGbkZxAbHck+Xe7i69dVEB0XXSgwuJ3ytdalS6h7ge8AKvKu13qSUeh5I01rPAd4BPlBK7QCOYvxREMLr2ZRNWvg+bveJ3Xy05SNm75hNfmk+3WK78UD3BxiYOLC8625tcUsNX2s9F5hb6bV/VFgvBK51x7GEqE+sFqvU8H2Q1pq0g2m8v+l9Fmcsxs/ix9DkodyQcgMpDVM8FledumgrhLexKqu08H1IiaOEH3b9wIxNM9hydAuRAZHc0ekORrUdVWtlm7ORhC+EiWwWm9TwfUBeSR6fb/ucD7d8yP68/SRHJPPMBc9wefPLCbQFejq8cpLwhTCRVVllpK0XO1xwmJlbZvLx1o/JKc4htVEqT/Z6kosTLsai6t5UZZLwhTCR9NLxTntz9jJ943S+2vEVJY4SBjcbzE3tb6JjTEdPh3ZWkvCFMJHU8L3LjmM7eHvj28zbOQ+LsjC8xXBu6nATzcKbeTq0apGEL4SJbBab9NLxApsOb2La+mks2ruIIFsQN6TcwN/a/43Y4FhPh3ZeJOELYSJp4ddvaw+tZcr6KSzPXE6Yfxh3db6Lv6T8hYiACE+HViOS8IUwkdTw66ffDv7G5HWT+WX/L0QGRHJ/t/sZ1WYUof6hng7NJZLwhTCRjLStX9YeWsuktZP4Zf8vRAVG8XDqw1zb+lqC/YI9HZpbSMIXwkRWJSNt64MNWRuYuHYiP+/7uTzRj2wzkiBbkKdDcytJ+EKYyGqRGn5d9vvR35m0ZhKLMxYTGRDJQ90fYmSbkV7Toq9MEr4QJrIqqeHXRbuydzFx7US+3/U9Yf5h3Nf1Pq5PuZ4QvxBPh2YqSfhCmMhmsclI2zrkQN4BJq+bzOwds/G3+nN7p9sZ034M4f7hng6tVkjCF8JEVmWlSBd5Ogyfl12UzTsb3uGjLR+h0YxuO5pbO95Kw6DT7rTq1SThC2EiqeF7VpG9iJlbZvLWhrfILc7lihZXMLbLWJqENvF0aB4hCV8IE9mUzJbpCQ7t4Nv0b3l9zescyDvAxfEX80D3B2gd2drToXmUJHwhTCQ3QKl9q/av4pW0V9hydAvtGrbjxQtfpGdcT0+HVSdIwhfCRDK1Qu3Zlb2L8avHs3jvYuJC4njp4pcYmjy0Tk5T7CmS8IUwkUytYL7somymrp/KrC2zCLAFcH+3+7kh5YY6deORukISvhAmsinplmkWu8POZ9s+Y+LaiWQXZXN166u5p8s9Ptfz5nxIwhfCRNLCN8evB37lpVUvse3YNlIbpTKu5zjaRLXxdFh1niR8IUwkNXz3OpB3gPFp45m3ax5xIXGM7zeeS5pdglLK06HVC5LwhTCR3MTcPUrsJXyw5QOmrJuCQzu4q/Nd3NThJq+b3MxskvCFMJHcxNx1q/av4oWVL7Azeyf9m/bnsR6PkRCW4Omw6iVJ+EKYyGqRhF9ThwsO80raK3yb/i0JoQlMGjSJvgl9PR1WvSYJXwgTyUjb8+fQDj7b9hn/W/0/Cu2F3NHpDm7teKt0s3QDSfhCmEjm0jk/245t47kVz7E+az29Gvfiqd5PkRSR5OmwvIYkfCFMVHbHK6219CQ5i8LSQqaun8r0jdMJ8w/jXxf9i8ubXy4/MzeThC+EiawWK2CUKazK6uFo6qZfD/zKcyueY/eJ3YxoMYKHUx+mQWADT4fllSThC2EimzL+i9m1HSuS8CvKKc5hwuoJfLbtMxJCE5h2yTQuaHKBp8PyapLwhTBRWQu/1FGKv9Xfw9HUHUszlvL8iufJKshiTLsxjO06VvrU1wJJ+EKYqKyMIz11DCeKT/DyqpeZ/cdsWkS0YEL/CXSK6eTpsHyGJHwhTGSzOEs60lOHnzJ+4tkVz3Kk4Ai3dbyNOzvfKd96apkkfCFMVNbC9+WboOSV5PHfX//L59s/p2WDlrw+4HXaR7f3dFg+SRK+ECYqq+H7agv/1wO/8vTyp9mft5+bO9zM2C5jpVXvQZLwhTBRxV46vqTYXswba95gxqYZJIQlMOOyGXSJ7eLpsHyeJHwhTOSLNfxtx7bx+E+Ps+3YNq5tfS0Ppz5MsF+wp8MSSMIXwlS+VMPXWvPRlo94dfWrhPqHMnHgRPo17efpsEQFkvCFMJGv1PAPFxzm6eVPsyxzGf0T+vNsn2flVoN1kEsJXykVBXwCJAG7gJFa62NV7GcHNjif7tFaD3fluELUF75Qw1+euZwnlj1BXkkeT/Z6kuvaXCdz4NRRFhffPw5YqLVuBSx0Pq9Kgda6i3ORZC98RvlIWy8s6ZQ4SpiwegJ3LriTqMAoZv1pFqPajpJkX4e5WtIZAfR3rs8AFgOPufiZQniN8pG2XlbSyczN5NGlj7I+az3Xtr6WR3s8KvPV1wOuJvxGWuv9zvUDQKMz7BeolEoDSoGXtNZfVbWTUup24HaAxMREF0MTwvPKa/heVNL5cc+PPLn8SbTWvNLvFS5NutTTIYlqOmfCV0otABpXsenJik+01loppc/wMc201plKqebAIqXUBq31H5V30lpPA6YBpKamnumzhKg3ymr43nCbwxJHCa+tfo0Zm2eQEpXC+P7jaRrW1NNhifNwzoSvtR58pm1KqYNKqTit9X6lVBxw6Ayfkel8TFdKLQa6AqclfCG8jbe08A/mHeThJQ+zNmst17W5jkd6PEKANcDTYYnz5OpF2znAGOf6GGB25R2UUpFKqQDnejRwIbDZxeMKUS94Qw1/1f5VjPxmJFuPbeXlvi/zVO+nJNnXU64m/JeAS5RS24HBzucopVKVUm8790kB0pRS64AfMWr4kvCFTygfaVsPW/haa97d+C63zb+NiIAIZv1pFkOTh3o6LOECly7aaq2PAIOqeD0NuNW5/jPQ0ZXjCFFflY+0rWc1/PySfJ5a/hTzd89nSLMhPH/h84T4hXg6LOEiGWkrhInqYw1/94ndPPDjA6Rnp/NQ94cY036M9K33EpLwhTBR+UjbelLDX5a5jEeXPIrVYmXK4Clyj1kvIwlfCBPVl5G2Wmve3/w+E1ZPoFWDVrw28DXiQ+M9HZZwM0n4QpioPvTSKbIX8fyK55nzxxwuaXYJL1z4gkxn7KUk4QthorreS+dIwREe+PEB1mat5e4ud3NnpzulXu/FJOELYaK63Etn27Ft3LvwXo4WHmV8v/EMSRri6ZCEySThC2GiutpLZ1nmMh5a/BChfqFMv2y63FTcR7g68EoIcRZ1sYb/f1v/j3sW3kNieCIz/zRTkr0PkRa+ECaqSzV8h3bw6upXmb5pOn0T+vLfvv+Vi7M+RhK+ECaqKzX8InsRT/z0BD/s/oFRbUbxWM/Hyv8YCd8h/+JCmKgu1PCzi7K5/8f7WX1wtYyc9XGS8IUwkadr+AfyDnDn/DvZnbObl/u+LJOf+ThJ+EKYyKIsWJTFIyNt04+nc/v828kryWPq4Kn0jOtZ6zGIukUSvhAmsyprrbfw12WtY+zCsdiUjemXTadNVJtaPb6om6RbphAms1lstVrDX565nNt+uI1w/3A+GPaBJHtRThK+ECazKmut9dL5YdcP3LPoHpqFN+P9oe/LPWfFKSThC2Eyq8VaKy38L7d/ySNLH6FjdEfeufQdooOiTT+mqF+khi+EyWqjhv/h5g/5z6//oU+TPrza/1UZUCWqJAlfCJPZlLk1/Lc3vM1rv73GoMRBvNz3Zfyt/qYdS9RvkvCFMJnVYk4NX2vN5HWTmbxuMkOTh/Kvi/4lo2fFWclvhxAmsyqr2/vha6157bfXeGfjO1zZ8kqeveDZ8lG9QpyJJHwhTGaz2Nxaw9da87/f/se7G9/l2tbX8lTvp7Ao6X8hzk0SvhAmsyr39dLRWvPqb6/y3sb3uK7NdTzR6wlJ9qLaJOELYTJ31fArJ/snez0pk6CJ8yJNAyFM5q4W/sS1EyXZC5dIwhfCZO6o4U9dN5Vp66dxdaureaLXE5LsRY1IwhfCZK720nlv43tMXDuR4S2G848L/iE1e1Fj8psjhMmslpqPtP3494+ZsHoCQ5OG8nyf5yXZC5fIb48QJqvpSNuv//iaF1e+SP+m/Xnx4heln71wmSR8IUxWkxb+wj0LeXr50/Rq3ItX+r2Cn8XPpOiEL5GEL4TJzreGv2r/Kh5Z8gjtG7bn9YGvE2ANMDE64Usk4QthsvNp4W8+spn7fryPZuHNeHPwmzLrpXArSfhCmKy6NfzdJ3Zz14K7CPcPZ8rgKUQERNRCdMKXSMIXwmQ2i+2cI20PFxzmjvl3oLVm6iVTaRTSqJaiE75EplYQwmTnuuNVXkkedy+4m6OFR3n30ndJjkiuxeiEL5GEL4TJznbHqxJHCQ8teYhtx7bx+sDX6RDdoZajE75ESjpCmMxmsVXZS0drzT9X/JPlmct5uvfT9E3o64HohC+RhC+Eyc7Uwp+2fhpf7viSOzvfydWtr/ZAZMLXuJTwlVLXKqU2KaUcSqnUs+x3mVJqq1Jqh1JqnCvHFKK+qWq2zLnpc5m4diKXN7+cuzvf7aHIhK9xtYW/EfgzsPRMOyilrMAkYCjQDhitlGrn4nGFqDcqz5a59tBanl7+NN0bdee5Ps/JzJei1rh00VZrvQU41y9sT2CH1jrdue/HwAhgsyvHFqK+qDjSdm/OXu5bdB+NQxrzv/7/w9/q7+HohC+pjRp+PLC3wvMM52unUUrdrpRKU0qlZWVl1UJoQpivbKRtbnEu9y26D7u2M2nQJBoENvB0aMLHnLOFr5RaADSuYtOTWuvZ7gxGaz0NmAaQmpqq3fnZQnhKWQv/8Z8eZ2f2TqZeMpWkiCRPhyV80DkTvtZ6sIvHyASaVnie4HxNCJ9gs9hwaAeLMxbzZK8n6RXXy9MhCR9VGyWdX4FWSqlkpZQ/MAqYUwvHFaJOsCpjHvuRrUcyqu0oD0cjfJlLF22VUlcBbwAxwLdKqbVa60uVUk2At7XWw7TWpUqpe4DvASvwrtZ6k8uRC1FPDEwcSKG9kLu7SPdL4VlK67pZKk9NTdVpaWmeDkMIIeoVpdRqrXWV46JkpK0QQvgISfhCCOEjJOELIYSPkIQvhBA+QhK+EEL4CEn4QgjhIyThCyGEj5CEL4QQPqLODrxSSmUBu8/zbdHAYRPCqct88ZzBN8/bF88ZfPO8XTnnZlrrmKo21NmEXxNKqbQzjTDzVr54zuCb5+2L5wy+ed5mnbOUdIQQwkdIwhdCCB/hbQl/mqcD8ABfPGfwzfP2xXMG3zxvU87Zq2r4QgghzszbWvhCCCHOQBK+EEL4iHqX8JVSlymltiqldiilxlWxPUAp9Ylz+0qlVFLtR+l+1TjvvyulNiul1iulFiqlmnkiTnc61zlX2O9qpZRWSnlF173qnLdSaqTz33uTUmpmbcfobtX4/U5USv2olFrj/B0f5ok43Ukp9a5S6pBSauMZtiul1OvOn8l6pVQ3lw+qta43C8YtEv8AmgP+wDqgXaV97gamONdHAZ94Ou5aOu8BQLBz/a76ft7VOWfnfmHAUuAXINXTcdfSv3UrYA0Q6Xwe6+m4a+GcpwF3OdfbAbs8Hbcbzrsv0A3YeIbtw4DvAAX0Bla6esz61sLvCezQWqdrrYuBj4ERlfYZAcxwrn8GDFJKqVqM0QznPG+t9Y9a63zn01+AhFqO0d2q828N8E/gP0BhbQZnouqc923AJK31MQCt9aFajtHdqnPOGgh3rkcA+2oxPlNorZcCR8+yywjgfW34BWiglIpz5Zj1LeHHA3srPM9wvlblPlrrUiAbaFgr0ZmnOudd0S0YLYP67Jzn7PyK21Rr/W1tBmay6vxbtwZaK6WWK6V+UUpdVmvRmaM65/wscINSKgOYC9xbO6F51Pn+vz8nm0vhiDpHKXUDkAr083QsZlJKWYAJwI0eDsUTbBhlnf4Y3+SWKqU6aq2PezQqc40GpmutxyulLgA+UEp10Fo7PB1YfVLfWviZQNMKzxOcr1W5j1LKhvH170itRGee6pw3SqnBwJPAcK11US3FZpZznXMY0AFYrJTahVHjnOMFF26r82+dAczRWpdorXcC2zD+ANRX1TnnW4D/A9BarwACMSYY82bV+n9/Pupbwv8VaKWUSlZK+WNclJ1TaZ85wBjn+jXAIu28AlKPnfO8lVJdgakYyb6+13ThHOestc7WWkdrrZO01kkY1y2Ga63TPBOu21Tnd/wrjNY9SqlojBJPem0G6WbVOec9wCAApVQKRsLPqtUoa98c4G/O3jq9gWyt9X5XPrBelXS01qVKqXuA7zGu7L+rtd6klHoeSNNazwHewfi6twPjgsgoz0XsHtU87/8CocCnzmvUe7TWwz0WtIuqec5ep5rn/T0wRCm1GbADj2it6+232Gqe80PAW0qpBzEu4N5Y3xtySqlZGH+4o53XJp4B/AC01lMwrlUMA3YA+cBNLh+znv/MhBBCVFN9K+kIIYSoIUn4QgjhIyThCyGEj5CEL4QQPkISvhBC+AhJ+EII4SMk4QshhI+QhC+8llLqLec8+a9Use0b57aHPRGbEJ4gCV94s87Ox/VVbOvifFxrdhBKqSilVKlS6ohSymr28YQ4E0n4wptZga3A6oovOuefKZtm1vSEDwx1xjJPa22vheMJUaV6NZeOEOdDa939DJvKWveZWuvDtRDK5c7Hb2rhWEKckbTwhS8qK/XURjnHBlwKlALzzD6eEGcjCV94JaXU22e6YEuF+r1Syk8pda/zhvcnlFIHlFJfKKWauymUC4FI4OeyWxI648tyxlfdZZab4hE+TEo6wludrRVflvAPAWlAJ4zpZxXQCLgK6KGUStFa57oYx2nlHKVUM4xpjfOBYIz78WZXel8oEAIUAceBH12MQwhp4Qvv4+wJ08H5dF2lbQFAW+fTZ4A8jFZ4WYK9FqP8kgBc54ZwTkv4WuvdWuvGnLzJxyStdeOKC/CZc9vbztemuSEW4eMk4Qtv1AbjjkhFwJZK2zpw8pvtNqC/1vpnfdJnnEzOrV0JwlkWaguka60rxwHQ0fm4oYptnc+yTYgakYQvvFFZstystS49w7YSYLTWuriK95fV2l29O9AVzsfTeuc4b+XXxvl0Q6VtNiClqm1CuEISvvBGZUl9XRXbyur3C7TWu87w/rILtmfaXl1n646ZgvFNww5srrStLRDgXN/oYgxClJOEL7zR2UbRlm37+izvLyu11DjZKqXCgL5ADrDkLMfYrrUurLSt7A/Wbq31iZrGIERlkvCFN6qyha+Mu7t3cj49ZfRthX0SgSiMck5V3xCq61LAH/jhDGWjsjikfi9qjSR84VWUUrFAY+fTygk7GYgAHJy59V72DWCn1jqnwufGK6UKlFIzKx3vAufr71T6nHONri1r4Vc1z8/Z/hgIUWOS8IW3KUvYeyoOdHIqazlv01rnn+H9Vfbf11pnAtOAkUqpVgBKqWRgNvAzcGfZvkopC8b8OQ5g7hmOIz10RK2ThC+8TXUu2J5tSoWyfap6/7+BYuBxpVQE8C1wBLhaa11SYb+eQCzwq9b6UOUPUUpFcXLytso9dBpy8huKXLAVbiUjbYW3qc4I2+ok/NP20VofUEpNBu7F6FIZDfTSWh+vtGt1yzm5wM5K25IrrB88S5xCnDdp4Qtvc7YW+lkTvlIqnJMJ90wXbF/BmOo4FRihta6csKH6CX+j1rpyX39HhfXXlFJjzvAZQpw3SfjCazinTSgbzFS5h04DINH59Ewt/E4Y8+kc11rvPsM+T2L8v7ECR6uIIQHjW0aG1vpMxzlb/X49RpfREow/HM3O8BlCnDd1egNDiPpJKdUDWIUxEVlkFa1nVz//AeBV4H7gn8C3WuvrK+1zJzAZmKq1vvP0TxHCc6SFL7zJpc7HFSYk++HAeOBFrfXrwBvAdUqplEq7ys1ORJ0lCV94BaVUF2Cs8+lsN392N2AmxgyWTztffhVjeuNnK+2+FHgOWOjOGIRwBynpiHpNKXU7xoXUMOdLa4A+VUxXUNPPTwBWAnuAARU/Vyn1MvAw0ElrLV0oRZ0nCV/Ua84RrtdjTHT2JfCSzD8jRNUk4QshhI+QGr4QQvgISfhCCOEjJOELIYSPkIQvhBA+QhK+EEL4CEn4QgjhIyThCyGEj/h/qtwQGHOLsQMAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "jhx = np.linspace(0, np.pi, 100) / np.pi\n", "\n", "#lmbdas = [0.1, 0.2, 0.4, 0.8, 0.9, 1.0, 1.2, 1.4, 1.8, 2.0]\n", "lmbdas = [0.25, .33, 0.88]\n", "\n", "for lmbda in lmbdas:\n", " plt.plot(jhx, alphaoverj(jhx, lmbda), label='%g' % lmbda)\n", "plt.legend()\n", "plt.xlabel('$j h_x / \\pi$', fontsize=24)" ] } ], "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.8.3" } }, "nbformat": 4, "nbformat_minor": 4 }