{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Tutorial 7 - Fourier Transforms\n", "\n", "Let's play ..." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "\n", "import matplotlib\n", "import matplotlib.pyplot as plt\n", "from matplotlib import rc\n", "import numpy as np\n", "\n", "rc('text.latex',preamble='\\\\usepackage{libertine}\\n\\\\usepackage[libertine]{newtxmath}')\n", "rc('font',**{'family':'serif','serif':['Linux Libertine O']}, size=18)\n", "rc('text', usetex=True)\n", "rc('contour', negative_linestyle = 'solid')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# do Fourier series first" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# then fourier integral" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsIAAAFoCAYAAABdfZ3MAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAACKn0lEQVR4nO29eXwV5d3+f53sCWRlCwoKJxZFXGrAFcUtcbdWJe57C7S1trXtQ2rr02qrDw/RLvb32Ap172ZJ6l4qm61YdwgoggvmsMkqkARIQrYzvz/yvYf7TGa5Z87MWcj1fr14cc6c2c45OTPXXHPdn09I0zQNhBBCCCGEDDAykr0DhBBCCCGEJAMKYUIIIYQQMiChECaEEEIIIQMSCmFCCCGEEDIgoRAmhBBCCCEDEgphQgghhBAyIKEQJoQQQgghA5KsZO9AuhGNRrFlyxYUFhYiFAole3cIIYQQQogBTdOwd+9eHHLIIcjIsPZ9KYRdsmXLFowePTrZu0EIIYQQQhzYtGkTRo0aZfk6hbBLCgsLAfR9sEVFRUneG0IIIYQQYmTPnj0YPXq0rtusoBB2iYhDFBUVUQgTQgghhKQwTjFWDpYjhBBCCCEDEgphQgghhBAyIKEQJoQQQgghAxIKYUIIIYQQMiChECaEEEIIIQMSCmFCCCGEEDIgoRAmhBBCCCEDEgphQgghhBAyIKEQJoQQQgghAxIKYUIIIYSQgNA0DZqmJXs3iAUUwoQQQgghARCNRnHaaafh3HPPpRhOUbKSvQOEEEIIIQcj69atw9tvvw0A6OjoQEFBQZL3iBihI0wIIYQQEgC9vb36456eniTuCbGCQpgQQgghJABkIdzd3Z3EPSFWUAgTQgghhARAZ2en/rirqyuJe0KsoBAmhBBCCAmAjo4O/bFRCHd3d+Of//wntm7dmujdIhIUwoQQQgghAWAnhM8//3xcdNFFuPXWWxO9W0SCQpgQQgghJADshPC//vUvAMCCBQsSuk8kFgphQgghhJAAsBPCJDWgECaEEEIICYD9+/frjymEUxMKYUIIIYSQAKAjnPpQCBNCCCGEBACFcOpDIUwIIYQQEgAUwqkPhTAhhBBCSABQCKc+FMKEEEIIIQFAIZz6UAgTQgghhAQAhXDqQyFMCCGEEBIAVuXTNE2zXe7DDz/EpZdeihUrVgS2b6SPrGTvACGEEELIwUhnZ6f+WBbCvb29tstVV1dj27ZtWLp0KVpbWwPbP0JHmBBCCCEkEGTBKwvhnp4e2+W2bdsGANizZ08wO0Z0KIQJIYQQQgLAqxAmiYNCmBBCCCEkAGQh3N3drT82CuFoNJqwfSKxUAgTQgghhASAqiPslBkmwUEhTAghhBASAKpC2CoqkZXFmgZBQyFMCCGEEBIAcuTBrmqElRDOzc0NZseIDoUwIYQQQkgAxOsI5+XlBbNjRIdCmBBCCCEkALwIYbnZBh3h4KEQJoQQQghRoK2tDU888QS++OILpfllISw/thPCbW1t+mMK4eChECaEEEIIUeCHP/whbrvtNpx33nlK88viV84LqwrhzMxMr7tKFKEQJoQQQghR4NlnnwUArFy5Uml+WfyqCmH5sdyimQQDhTAhhBBCiAKFhYWu5vcSjZAf79+/3+0uEpdQCBNCCCGEKDB48GBX81tFI+zKp9ERTiwUwoQQQgghClgJ4Wg0igceeADvvPNOzHQvjrA8H4Vw8LBlCSGEEEKIAlbRiKeffhozZ84EEFv+zMtgOaMjrGkaQqFQfDtOLEmKIxyJRFBdXY3S0lJMnDgRdXV1CVlPQ0MDqqurUVFRgdLSUlRXV6OhocHTtgkhhBAysJAdYVnkWg2es4pDqAphILb+MPGfhAvhSCSCiRMnIhwOo7m5GfX19ZgzZw5qamoCXc+MGTMwa9Ys1NfXo6mpCc3NzQCAmpoaz0KcEEIIIQOHQYMG6Y9bW1v1x1ZiVdURtopQAIxHBE3ChbAQqnPmzAEAhMNhzJ49Gw0NDWhsbAxkPYsXL8bcuXNRVVWFkpISfbpYtra21vP7IYQQQsjBRWdnJ1asWBETcwAQE1EQhhoAdHd3m64n3qoRACtHBE1ChXBjYyMaGxtx1VVXxUyfOnUqAHVB6mU9JSUlqKioiJkWDof1x5FIRGnbhBBCCDm4ufXWW1FZWYnHH388ZroseOXGFypC2EtGGKAjHDQJFcLCga2uru73WjgcxuLFi9HS0uL7eqqqqtDc3Izp06fHzCvPI4tiQgghhAxc/vrXvwIAfvKTn8RMlwWvLHLdOsJ25dOMr1mtm/hDQoXwsmXLAJiLThFZUHFm/VrPvHnzAKCfQJbp7OzEnj17Yv4RQggh5OBn+/btMc+tBKuVWI23s5xxO8R/EiqEhTiVc7qCsrKymHkSsZ7a2lpUVlbqDrMZs2bNQnFxsf5v9OjRjuslhBBCyMGHX46wGyFsfE78JaFCWCX2oCJg/VhPTU0NysrKsGTJEtv57rrrLrS2tur/Nm3a5LhtQgghhKQv5eXlptNlwSsLVHl6vHWEjQ4wHeFgSbmGGmYur9/rmTFjBhobG7F8+XLH7eXm5iI3N9eXfSKEEEJI6mN13ldxhHt6epCdnd1vHq+OMIVwsCTUERaZXjtHV0QbglrPjBkzsHjxYiURTAghhJCBh7FsmkAlIyzXFPajagSjEcGSUCFcWVkJwDy2IKZVVVUFtp66ujrMmzePIpgQQgghlshCWBawKtEIWQh76SzHaERiSagQvvrqqwFYC9jKykolgeplPQ0NDaitrcWSJUv6vVZTU+OqmQchhBBCDl5kIWwleGWBKotfFUfYKG5V3WLiPwkVwlOnTkU4HO5XpaGhoQFA38A0mZaWFtTV1WHx4sVxr2fatGmYPn267ibLrzU0NLCOMCGEEEL6YeXkytPb29v1x1ZC2EoUG58zI5xYEj5Yrr6+HhMnTsTcuXMxffp0RCIR1NbWYvr06XpnOMHcuXNRW1uLkpKSmFaGXtbT0tKCZcuWYcaMGTGvLV68GCUlJYxKEEIIIQSAe0dYnm41j6ojzGhEYkm4EK6srMTy5csxbdo0zJ49GyUlJbqANVJVVYWSkpJ+rZTdrue9994DcKA1s9l2CCGEEEIA6xJoVhlhq5iEH44woxHBkpTyaULEqsxndIK9rKe+vt7V/hFCCCFk4GLlCFuJX1nIWg2Ws3psfM5oRGJJaEaYEEIIISTVsRKmVrEHPx1ho/ClIxwsFMKEEEIIIRKqjTMEVo6wakbYrrQaHeFgoRAmSWXfvn1sW00IISSlsIpDqDjCYh5N05TaLRufUwgnFgphklSmTp2KI444Ah999FGyd4UQQggBYC143WSE7SpDMBqROlAIk7j4+OOPceedd2L37t2ul92xYwcWLFiArq4uPProowHsHSGEEOKeeBxhIYSNYpeOcGqSlKoR5OBA0zSMHz8eAJCTk4PZs2e7Wv7VV1/VH69atcrXfSOEEEK8oGmaUitlp/JpbmoFs3xa8qAjTDzz4Ycf6o/ff/9918uvX79ef/zZZ5/5sUuEEEJIXNg5ufFEI+wcYSt32ew58RcKYeKZjRs36o937tzpevnt27frj7dt2xYzqIAQQghJBrLrCxwQv9Fo1NLVNRss5zUjzGhEYqEQTgPee+89TJkyxZPrGiRbt27VH3/88ceuhey2bdv0xx0dHdi3b59v+0YIIYR4wUqIWglkIH5HmNGI5EEhnAZceeWVeP311zF58uRk70oMshBua2vD3r17XS0vO8JmzwkhhJBEYyV47ZxaM5HrpjKEao1h4j8UwmmAqLPb1taWUvEBWQgDsG2HbYbsCJs9J4QQQhKNlSOrWgVCTKcjnB5QCKcBQ4cO1R+nkmtqFMJuS6iJ91JcXBzznBBCCEkWRkdYCFijeHUaOOcmI8zOcsmDQjjF2b9/P3bt2qU/91KvNyh27NgR89yNIxyNRvX5jzrqKADeBtwRQgghfmLlyNo5vGaP/XKEKYSDhUI4xdmyZUtMHKKlpSV5O2OgtbU15rkbkS7HPA477DDT9RFCCCGJRtURFtONkUUrIew1I8xoRLBQCKc4RnHoNocbJHv27AEAjBw5EoC7fRPLZmVlYcSIEQAohAkhhCQftxlhK1Fr5/oeDI7w0qVL8Yc//CHZuxE37CyX4qSyEBb7dvjhh2Pr1q2uHGEhhIuKivSMMIUwIYSQZGNXR1hl+kAZLHfmmWcC6Is3nnHGGUneG+/QEU5xUlUIR6NRvVzamDFjAHhzhGUhLKYRQgghycLKkXUreL021EiH8mltbW364w8++CCJexI/FMIpjlEcpkpGeN++fXou6vDDDwfgTggLgU9HmBBCSCqh6ghbCWSrqhGqbZTTIRqxfv16/bGxglS6QSGc4qSqIyxnfEV5N/kKUXV5CmFCCCGphJUQtRrEFo8jnJ2dbflaVlaW6f6kArIQ/uijj5K3Iz5AIZziGB3hVBHCQrQWFxdj8ODBALwJ4eLiYhQVFcWskxBCCEkWbh1hKyFsFLSapul3Uo2vmUUjcnJyTNefCshCON2bYVEIpzhCHIofhNs2xkEhO7pCCO/bt8/T8swIE0IISRWsnF+30Qij2LV7zcwRFuf9VHSE5br/qaJLvEIhnOIIISxKlHV0dCRzd3RkR3jQoEEA3DnCzAgTQghJRVQFr6ojLOIPZsvYOcJiuVR0hGXji0KYBIpwScvLywEA7e3tydwdHbFfhYWFnhxhMe/gwYNRWFjoenlCCCEkCFQdYavpRtdXFsJWr5mVVkvlaIR8vk73czeFcIojrrRE04lUEcJiPwYNGuTJERbOdkFBAQoKCgAAnZ2dKfmDJ4QQMnCIt0yalesrT1MZLJfK0Qg6wiRhCME5bNiwmOfJRojeQYMGeXKExfvIz8/XhTCQOtEPQgghA5N4M8J20QiVjHA6DJaTxW9nZ2e/AYbpBIVwiiME45AhQ2KeJxszIezVEc7Ly+u3XkIIISQZqApeq/JpdtEIlYxwukUjgPR2hSmEU5xUFcJW0QjjgcJp+fz8fGRkZOiucKq8P0IIIQMTt46wm8FyViLZbrBcqkcjzJ6nExTCKY5wTkXTilQRisK5LSgo0B1hQD3aIDvC8v+p8v4IIYQMTFQHxalGI+wywnSEkw+FcIojhKEQwqmSoZWjEfn5+fp01atC2REGKIQJIYSkBl4HxQmMrm9mZiZCoZDpa0II21WNSAdHmEKYBIZRCHd3d6dEKF2ORmRkZOjxCFUhbOUIMyNMCCEkmag6wqpOcUZGBjIyMmJeM8Yf0nWwXGZmJgBGI0iACMEoMsLytGQiRyOAA87u/v37lZY3OsJCSNMRJoQQkkycMr/G6U6OcEZGhi4YVTLC6eQIi2ZfFMIkEKLRqC4sS0tL9empIBblaAQAvfKDqhBmRpgQQkgq4ndG2MwRVimfZiaSU4Hu7m50dXUBOHC3WvXcn4pQCKcwsvM7aNCglBKLRiEcryOcSu+NEELIwEXVERbPnea3c4TtBsulaotlWZsIk45CmASCLArlxhOpIBbljDBwwBGOt2oEM8KEEEKSiargdZruJiMsr8PutVRAFr3FxcX9pqUbFMIpjBCLubm5KVdr15gRdhON0DRNf2/MCBNCCEklrIRwPNEI4Qi7abGc6kI4NzfX9d3gVIRCOIURolCIzdzcXAB97QyTjZUjrPJj6O7u1n/YzAgTQghJJYKIRghHWNO0mNdUohGplhEW5/m8vDzX44NSEQrhFCaVhbDR0XXzY5DjE8wIE0IISSVUHWHV6bIQVskIp0s0gkKYBI5xQJkopaIihP/5z39iwoQJeP755wPZt3iEsHhfGRkZ+nuiECaEEJIKBOkIG5dJx85y4jyfn5+vn/tTwaDzCoVwCmMcUCYcYVG2xI7/7//7/7BmzRpcfvnlgYhL+YoQcFc1QhbRotuOiFhwsBwhhJBk4lXwWk23c4RVBssxGhEsFMIpzFlnnYUdO3ZgwYIFANxFIzZv3qw/Xrduna/71dvbq3e3Ez8CL46w3JqZjjAhhJBkYDxvOVWBME53UzXCj/Jpr7zyCm6++eaktTWmECYJIzs7G8OGDUN5eTkAdSGsaRqampr05xs2bPB1v+TtG4WwSvk0o9MtP6YQJoQQkijWrFmDwsJCfP/739enqVaHiCcaEU/ViAsvvBBPP/00fvSjH6m/UQueffZZvP76666WkYWw0CUUwiQhiLyQUzRi69atMREDv4WwLHbpCBNCCElXfvrTn6Knpwe//vWv9WmqjrCTQPaaEVYdLLd06VLnN2jD559/jiuvvBJTpkzB7t27lZejI+wDkUgE1dXVKC0txcSJE1FXV5fQ9UQiEdTU1KCxsdHTdpOFqiO8ffv2mOd+C2HxB5+VlaX/iL1UjZAdYWaECSGEJBqzc06qOcJWGeH169dbvi8VPv/8c/3xM888o7wchXCcRCIRTJw4EeFwGM3Nzaivr8ecOXNQU1MT+HoikQhmzJiBiooKNDQ0oKSkJM53k1hUhbDxyk7OC/uBcaCc/JiOMCGEkHTB7Jzj1EHOStQalw+6fNq+ffsc3p09u3bt0h+vXbtWeTkK4TgRQnXOnDkAgHA4jNmzZ6OhocGVQ+t2PXV1da7FdqqhWjXCKITd3PJQwVg6DYjfEaYQJoQQkmjMzjlOUQchXo3tkq3mVxHC8joSVT5t586d+uPm5mbl5cQ5nELYA42NjWhsbMRVV10VM33q1KkAgNra2sDWM3PmTCxfvhxz5sxBZWWll91POqp1hI1/0H4LYTNH2E35NDrChBBCUgEVR9ipLbKVIyw7yFYZYbP4Q6IaashCmBnhBCHc2+rq6n6vhcNhLF68GC0tLYGvp6ysTG2H0Sc69+zZE/MvWbiNRoTD4ZjnfhFvNIIZYUIIIamASkbYKdebyIxwT0+P0vvq6elxrD/s1RE2E8JsqKHIsmXLABwQaDIirxuJRBK2HhVmzZqF4uJi/d/o0aN9Wa8XVIWw+IM+4ogjYp77hV00wk35NDrChBBCkoldNMIYW7DK9TpFIzIzM32rI2wU7mbCeP369SguLsY3vvEN8zf9//BTCNMRVkSIU7NBasKlVRGwfq1Hhbvuugutra36v02bNvmyXi+olk8TDnBFRYX+XNM03/bDr8FyZhnhzs7OlGsnSQgh5OBEFsLiPGkVWzCWPPOjs5zbwXJGIWx2zp03bx7a29vxhz/8AV988YXle/cjGiEMOhUTLFVJqBBWiT2oCFi/1qNCbm4uioqKYv4lC7fRCCGEe3t7fe1A41c0wswRll8nhBBCgkQ2icS51SkL7Ec0wk7s2jnCRgfb7Jwra6TFixebvOs+5KoTXhzh3NxcXZeIbrPpSMo11PCrpFm6lUZTQVUIixxzeXm5LlD9zAnHWzXCzBGW1+U2HtHd3Y0HHngA77//vqvlCCGEDGyEQAUOnL+sRKqVi+s03W1G2G4gnYojLJdCM/YVkJFNp/379yubUOKudF5envKd6lQmoUJYZHrjHcjm13rSDdXyaeIqb/DgwSgtLQXgb044CEc4FAp5vsUyd+5czJw5EyeccAK2bt3qallCCCEDl1AopD8W5y9VR9hYHcJKOLuNRrjJCJudcz/77DP9sRx/MGI817a2tlrOKyPMuJycHH0fKYQVEWXLzGILYlpVVVXC1pNuqJZPk4WwiHIEHY2It3ya23XIPPHEEwD6bnG99NJLrpYlhBAycJFv6as6wqqRiXiFsFkdYRUhvGPHDv2xGyGsWrVJiN6cnBw6wm65+uqrAVgL2MrKSqVIQ7zr8bucWKJQjUbIQnjw4MEx0/wgXkdY7L+8vPzcjSPc1dWFlStX6s/fe+895WUJIYQMbMyEsFMW2KpqhBtH2GpdZuuTX1PJCMvOrt1gOeO5VlUnCNGbm5urC2FmhBWZOnUqwuGwXgdY0NDQAKCvQoNMS0sL6urq+oW93a7HCpVBd6mEajRCXNUNGjQIhYWFAPx1hOMtnyaEsHg/Ai+O8Lp162KulkVpPUIIIcQOTdNizh9GR1g1CyzmNzq4XjLCmqbpA/jMohHG87/xnNvT0xPj7CbKEVapW5yqJHywXH19PSKRCObOnQugz8Gtra3F9OnT9c5wgrlz56K2tta0NbKb9chEIhG9BbPdaMpUJEhHeOPGjZgwYQJuvfVWx/2I1xGWR5zKeHGEP/nkEwAHGnKsXbvW11JxhBBC4kPTNDz33HMxg7hSAaOLaeUIW0UjVCMTblosy2JSRQgbz7nGpl8qQlici1WFsFlGGLB3hZ999lncfffdWLp0qdI2EknChXBlZaXe6riiogI1NTWora3t5+4CfTnfkpKSfq2U3a4H6BPAFRUVekkxoK8Vc2lpadoIYpWMcG9vr/7HPXjwYGVH+Jvf/CbWrFmDJ5980vHH4CSEnYSon46wOLCKTHhbW1vaRl9Sld/85jcYNWoUvvGNbyh3NSKEEMHChQtxxRVXYNy4ccnelRishLBVNMEpMuFlsJxdu2YzIWy1zwLjgDcVITx06FAA8TnC8nQzXn75Zdx///148803lbaRSLKSsVEhYlXms6t2oLoeoK/SRFNTk/I+piIqjrD8h+zGEZav1FesWIHTTz/dcl67aATQ92MwilwZq4ywWJ8bR3jz5s0AgHHjxmHEiBHYvn07NmzYgCFDhiivg1jz+uuv48477wTQ19r82GOPxe23357kvSKEpCrf/va3sXTpUrz99tt6icz//Oc/pvNGo1GsXbsW48aNi6nekCi8OsJWkQk/MsKy6DXLCDsJYWPk08oE6+3t1dc1dOhQbN68Oa6MsDzdbhl5/lQh5eoIE2tUMsJCCGdkZCA3N1fJEe7p6cG6dev0504XF3aOsPy6FVaOsJdWjaJc2siRIzFmzBgAwIYNG5SXJ/Y89NBDMc8feOABRk8IIZY8/PDDWLVqFV5++WV9mpXI/f73v4+jjjoK//d//5eo3YvBKCqFCeO2s5yTg5wMR1iYQR0dHabdWmXDadiwYQC8RSPk9tF20Qir834qQCGcRqhEI+R8cCgUUnKEN27cGHPL+/PPP7fdDzMhnJOTox/svAphL46wEMLl5eU4/PDDAfS9HxI/+/bt08vRvf322xg0aBA2bNjAAYmEEFPk84zcqEIWwvL5S1xof//730/A3vXHKNzEeVDV+Y3HERb/Z2Zmmk4HDpzzo9GobkCoCuFDDz1Un2bWpEo+zwrR7CUaIf9vZ9JRCBNfUIlGyEIYgJIjbIyMbNu2zXY/rBpiqDq6QTnC5eXlAJz3n6ixdOlSdHV1YezYsTj55JNx0UUXAQD+8Y9/JHnPCCGpiNzQSHYh5btIZnHHZI09cBLCTtEIUeHBD0fY6C7Lr8nzqwrh8vJy/QLETOCK83h2djaKi4sBuC+fJgSwSlMNRiOIL6hEI8QfsqiioOIIG7uxOQlJM0dYfu7k6PrpCIt9HTlyJEaMGAHAvqXkwcTatWuxadOmwNb/6quvAgDOPfdcAAcGJKbiqF9CSGJpb2/HCy+8EOM2yucSOasqVzLws8tpvFgJYSdhK4SwmOa2fJpZiTS7aIQ83SrOIRCfdVFRka4DzM7/sqEl5nPrCItzOB1hkjDcDJYTf9jCEbYTwmJUqRg56iQknYSwk6PrVD5N1RHu6OjQf/QjRowYUEL40Ucfxbhx4xAOh/Haa68Fso13330XAPSBk1OmTAEAvPXWW2ldPJ0QEj/f+MY38NWvfhW1tbX6tC1btuiPZSEsP05lIey2TJqY5jYaITvkKoPl5OnGfTY+FxcmgwYN0o0wO0fYixCWM8Ly/3bnBTrCxBdUMsLiRyBG64ofgl00QnSeOeaYYwB4i0YA6kLWL0dYHFAzMjJQXFw8YIRwe3s7fvCDHwDoczDuuOMO3wew9fb26vW2TzzxRAB9lTkKCwuxf/9+fPrpp75ujxCSXvzxj38EgJiBbvK5Ix2FsJMjbJwuXnNbPk0Wu1Y5ZON2xDJG19X4XNYAdgLXD0eYGWGScORohJXwMYpUFUdYCOFjjz0WQJ9DbJfZitcR9isjLA6uJSUlCIVCA0YIv/jii9izZw9KSkqQnZ2NVatW4f333/d1G2vXrkVbWxsKCgpw5JFHAug7mIu/kQ8++MDX7RFC0pORI0fqj+XzjIoQNkYMEo1RuFllhJ2iEcb5xfnZSgibiV3jNkKhkD6QDnDvCMtC2Oz8L86z+fn5ymVWBcZohEpGmEKY+IL8B2R1C8IohN04wkcddRSAvh+xXfvpeISw3IYx3jrC4oBaWloKADFC+GAu8SUGq82YMUMfwDZ//nxft7FmzRoAwNFHHx1zMD7++OMBwHfhTQhJH2QxO2rUKP2xfOyXzyFykwdZbMm3yeWs8d/+9reElFSzikY4ObyyEO7t7XXtCMtC2KoUW2ZmZsyx1yojbOcI20UjZGGaCEeY0QjiC7IQtopHxOMIl5eX6z8cFSFsjEaoZJjl1/x0hIEDGefu7m7lK9t05PXXXwfQN4hNDGT797//7es2Pv74YwAHLo4Exx13HAA6woQMZGQhLN89lE0M+RwiH9NlsSSLOvnxNddcgzvuuCPwC27VqhFWlR7ENLdVI8yEsN0ygP+OsFguJyfHlRCORqP658RoBEk48pWUqhBWcYTFYLlhw4bpotJOCIttGB3deIWwV0dY7HN+fr6+zoO1zfLmzZuxYcMGZGZm4tRTT8WZZ54JoG8Am5+3Fj/66CMAwPjx42OmJ0oINzU14fzzz8dxxx2H559/PtBtEULcIbu3siiWBa98zpEFkngcjUZjRJw4N8jCevXq1T7udX9UM8Kqg+WsBr65iUbIy9hFI8T6vDrCYrns7GxX0Qh5exwsRxKOfKvE6srLSgh3dXVZLiNHDFSEsFU0QohQO0dXvJaRkRFze0len1tHWEQjQqEQysrKABy8QnjlypUA+pzawYMH4+ijj0ZeXh727duHzz77zLftCEfYKIRFRnjz5s3YtWuXb9uT6erqwmWXXYaFCxdi1apVmDp1KlasWBHItgghzuzfvz8mbiabFfKxVsX5FdONTR7EdHndmzdvjnfXbfHqCFtFI1TLp6k4wpmZmQiFQnotYKMQFi6uVTk1J0dYFqZuHGH5e2X5NJIUnFxXcXAxCmHA/Mcg54FLSkriEsLiuYojbPZjcOsIG6MRAA56ISxuFYqsblZWlu7S+iUWNU2zjEYUFhbisMMOAwB88sknvmzPyJ/+9CesXr0aQ4YMwcSJE9Hb24sf//jHgWyLEGLP1q1bMWbMGFx66aX6NFnE7tmzRxdpshCWzwNmjrBRdJkJZKcup/HilBG2cnidyqeJOsGyqFUVwrIjLP9vzAiLylBWjrBTNQixnuzsbFdCWP5exfvlYDmSUJxKqMlXg0DfH6j4wzMTwu3t7fpVsCyE5cENMj09Pfr88WSEzX4Mbh1h42A54OAXwiKSIMQvAHz5y1+OeS1eNm/ejH379iErKwtHHHFEv9fFND8daJlHHnkEAFBbW4t58+YhFArhlVdeQSQSCWR7hBBr7r//fmzfvj2mo6TRzRWmhGxiOAlh4/nITAgH2TAI8N5ZzjiIzSicAcQIYRVH2LiM2Ib4X9URlqMRduZUvI5wVlaW/p44WI4kFKfucmY1fu3aLIsDWGZmJgYNGqS3WrRyhGWRGk9GOChHWIjiVKpV6Sdr164FEOvUivJm4rV4EU5vOByOcT4EX/rSlwAEI4QjkQjee+89ZGZm4uabb0Y4HNY72jU0NPi+PUKIPfKdJiGUjMdocW5RcYTFdKtohDxd1RTxitfOclbC1sopVqkjLLYrLwP0F8LicxJml50QFudUs89RHiznJSMsC1qnjLCcB6cjTOLGSWyaCWG7P3Lh/IpavE7RCL+EsHFZeRodYXM0TUNTUxMAxDi148aNA+CfEBbOq5kbLE8PQggvXLgQQF83u+HDhwMAvvKVrwAA/vnPf/q+PUKIPfLxXFQYMopYMY9VRtjssVVjCHndducSP/DaWU6OOpiVTxPTnYSzsVawLJKNQtgYjRAurt1gObtzqjxYTnaXnbqGGmsIA86OsFmuOJWgEE4zvAhhFUdYCGDxv5WjKtafk5MTU9pFZd/k15gRds+uXbv0C5dwOKxPFw7t2rVrfamfvH79egDAmDFjTF8PUgj/61//AgCcc845+rQLL7wQAPCf//xHb6lNCEkM8nljx44dAPoLYSG0zBxhTdNMB8sZBZeY3ypeEQReB8vJFR3MyqeJ6U5CWBbUZq+JZYH4ohFm51SzaATgHI8wtleW37eKEGY0gsSNakbYzBG2E8IiElFUVGQ5L2A9UA5QqxoRREZ4oAhhITxHjRoV8/2OHTsWGRkZaGtrw9atW+PeTrKEsKZpej3ks88+W59eUVGBL33pS+jp6cGSJUt83SYhxB754lM4wkZhZSeEjV1KrYSwmSNslzn1wt///nf8/ve/1587RSPclkOz6jhnNb+xVrCdI2wlhL06wnI0IicnR993p3iEXTTC6vuS9QqFMIkbLxlhu1qCRlfVKSukIoTtruLF8n46wgMlGiGEpzGykJOTg7FjxwLwJx4hhLBYpxHhRjc3N/v6OX/00UfYsWMH8vPzcdJJJ8W8dt555wEAXnvtNd+2Rwhxxo0jbObmWkUgVISwfC7RNA3/+te/9Lr3Xpg6dSq+9a1v6fWJvUYjZAFrVj5NTPcihFUHy6lkhFWjEaFQSHnAXDzRCHmAXSqRentEbHFbPg2wF7dGR9hJCJsJbUG85dPoCNsjhHBFRUW/10Q84tNPP417O06OcEFBAQ499FAA/uWSgb7oAwCceuqp/f4+Tj31VADAO++849v2CCH29Pb2xggjq4ywmSPc3d0NTdN8E8IvvfQSzjnnHFRWVnp6L/J+rFu3znQf3AyWk6MRbgfL2XWPsxosp1I+TdO0GLNK1REG4FoIuxksl8ql0wAK4bQjnsFyZn/g4raXUQhb/RjidYT9yghHo1F931PVEf7000/xs5/9zLeyX2YD5QRyTjgeOjs7sWXLFgDWQhg44BYL0ewHYnT6iSee2O+1k08+WZ/H79ulhJA+urq6YkpnGg0R8VwlGiHWZ/y9OjnFVkJ4/vz5ANRLqn322Wcx8S3Z2RbnDq/l06yqQHiJRjgNlnOTEZY/UychLDvC8vqcohFmGWHVaASFMPEFpysv8Qev6giLA4LIBgcdjVBxhHt7ex1Hru7Zs0cfGJaKjvDOnTtxyimn4Oc//zkmT57syyAvIYTNHGEhjsU8Xtm4cSOAvoPikCFDLOc7/PDDAQAbNmyIa3syQgifcMIJ/V6rqKjAkCFD0NnZqTcVIYT4h6ZpuPzyyzF8+HD9Atc4VkScF1QcYaDveK/qCDsNlhs2bJj+2On80Nraii996Us48sgj9Xnlc5pwto35ZWNDDStH2Fg1wswRTlQ0wqw0HdB3jnXjCDuZYAIzR1h1sFwq5oMBCuG0w2mwnJlQtbvSEwc6UVlCVQibRSP8coTl7VghYhH5+fkx6xJCONl1hB999FF9H7Zt24Y//vGPca9TdFkSnd1khDCNtwC9HIsQrT3N8FsI9/b26g1BRIMQmVAopOeGGY8gxH/eeOMNzJ8/H11dXfqgVaMQFiJJJSMM9B3vrSIQbqMR4hwFHLhgt2LRokUAYu8cyu9FDCo2CmGnqhFWg+XcOsLyeuTjrMpgOfE5mTnCxkFpdnWEjeI0yIwwHWHiK17+4FQcYVUhLA508VaNsFte3o4VRidbIIRwe3t74MXY7XjuuecAAEcffTSA+JtBRKNR/eB9yCGH9HtdiGOnE4QTQggLoWuFX9sTfPrpp+jo6MCgQYP0mIcREY947733fNkmIeQAjY2N+uNt27YBQL87WVYNNfbv3w9N0/qZIG4cYSchLE8XGV/jfglef/31fsvJ57TNmzcDsBbCRkfYLBphVj7NqpWynSMcCoV0MWx8TaxT3rZd1QhxzsvOzkZGRoanaISX8mlOd6rNxHMqQSGcZjgJYTNH2O6Wh7hKFoLSKScUZDQiIyNDSUzL+yfem6CoqEg/cCQrHtHa2oply5YBAB566CEAwFtvvRWXMN+5c6d+kC4vL+/3uhCm27dvj2s7dq6zjN+OsIhFHH/88ZajioVTvGrVKl+2SQg5wK5du/THVtEIO0fYqo2vWyFs1ZRDPifJx/Z77rkHRUVFePPNN/VpckMosc/yexE5aKtohEodYbeC16qOsHjd6TWVqhFGo8lLNCLI8mmMRhBfSJQj3NbWFtMPXWAnhFWqRtiVT5PX4eQIi4ObUQiHQiHHNtFB8/bbbyMajaKiogLnnnsuysvL0dnZiXfffdfzOsUAtuHDh5u2PS4rK9MPjPHEI4RTIqpCWOG3ELaLRQiOPfZYAMDq1av7ncAIIfHhRgiL47NsXMjnJGGoeHGE5end3d36eUg2cuRzzL333otoNIrvfve7pq+L5eTzn7HGsTivWjnCgLWwlcunZWZmxkQZVASy/L+8LqtohNhHM0fYeP5Xaajh1hH2khFmNIL4ip0Qlm9Nuc0IGwfLAf2v+gH78mnxOsLyelUdYbkjjsCpTbQdr7/+Op577rm4OrStXLkSADBp0iSEQiG9CoKY7gUhhM1iEUDfBYAfcQVVISy2tWfPHl8uOD766CMAB6IkZowZMwaDBw9GV1eXL2XiCFFB0zR8+OGHeOyxx/DAAw/gz3/+sz7Y6mBCFsLiYtoqGiHODaJij9ERFucRMyEs5rMaLGclkOXzl9n5QRatZpEKWdSL14W4FOdLq4ywmCa7tWbl07xEI8RyZuuS/zcKYXGuNHOEjUI4iPJpXuoI0xEmvmD3B2fVz9vOETYOlsvPz9fzSmY/iCCjEfJ6vTrCwAEhLJcBUmHBggWYMmUKrrjiCjz88MOulpURgldUPxAuZzzVDpyEMHDApY1HCItohJMQlqtK+JETFkJ4/PjxlvNkZGTorrBwkINE0zT87ne/w1lnnYXrr7/e15rJJD148cUXUVlZiWOPPRZf//rXMXPmTNxwww0oLy/Hbbfdpv9eDgZkISwG+hqNEqMQFmMyOjo6YpomyHcH43GExTrkbQPmwk4eqGbmCJsJYSEqxfnIWDXCOPjNabCcl8iE/L/Ta/I+mwlh4x1XlYyw22iEl4wwHWHiK3ZCWP7xq2aEjYPOMjIybB3koIWwW0fYTgi7dSoffPBB/fH//M//mEZDVPjwww8BAMcddxyAvtwrELwQTqQjDPgXj+js7NRrLR911FG284rPNBFC+Ic//CFuv/12vPbaa/jLX/6Ck08+GZ988kng2yXJp62tDddddx0uu+wyrFy5Ejk5OaiqqsL111+PE044AdFoFE888QSOPPJIPPXUU8neXV+wE8JiXIKxjrAQwrIjnJubG9MB1W3VCCshbBZtkO/cye6tfH5UiUaoOMJWUQd5ulVZNTdCWDUaIc6VZtEIY0a4s7Oz311OP6MRrBpBEopdi2VZPMp/pG4cYaf5VaIRKlUj4nWE/Y5GtLe3x7Tv3bp1qy5o3RCNRvVavuPGjQNwQNx99tlnniMXiRDCHR0d+gkwkUL4s88+Q29vL4qKijBy5EjbeRMlhF999VX86le/AgD89Kc/xaRJk9Dc3IwbbrhBPyERe6LRKJ588klMnjwZxcXFGD58OKZOnRpXVj4R7Nq1C2effTb++te/IjMzE7W1tdi6dSsWLVqEP/3pT2hsbMRbb72FyZMno729HbfccgtmzJiR9n8XshDev38/9u/frxsl4ndpF42QBZJsiojp4tjupyMs3/WzcoTtohFuhLBVIwyjeDWLOVgJZztH2CkaIQ+WE+cVq2iE8TMRywGMRgAUwmmHiiMsSqcIrIStpmn9bn3ZzQ/45wibLQ+oO8Iq0Qg3Qvidd95Bd3c3Dj30UFxwwQUAgCVLligvL9i6dSs6OjqQlZWlC8VwOAyg76Atn2zc4EYIexWmwg0eNGiQPuDQDr9KqMmxCLvaxQAwYcKEmGWC4p577gEAfPOb38S9996LF198EcXFxVi2bBnq6+sD3fbBwO7du1FdXY1bb70Vb775Jvbs2YMvvvgCf//733HyySfjhz/8YUoOeGxvb8cll1yC9957D0OGDMFrr72G//3f/9WdT8Epp5yC1157Db/4xS+QkZGBuXPnYvr06XGNLUg0+/bti9lf47Gpubm5nyNsFY2QhbDsCMtC2Di4yw9HWM5qWzXiEPssn1OsohF2g+XsHF635dPsusc5lU8Ty4pzpaZp+jQ7IWzW9Q/wp8UyB8uRhGInhK1EqlXUoa2tTT8Qyo6w12iEStWIRDjCXqpGiNq0p512GiZPngwgtq6mKqKl5+GHH647FPn5+Rg1alTM625REcKjR48GAM+5RTkW4SRIAf8cYSFqnWIR8jzr1q0LrE7066+/jtdffx05OTm4++67AfQ5Yt///vcBALNnz04rwZNompubUVVVhVdffRUFBQX4n//5H6xZswZvvfUWbr75ZgDAL3/5S1x77bUpJYZ7e3tx7bXX4u2330ZpaSlee+01/VhgRmZmJu6++24888wzyMzMxOOPP45f//rXCdxj7zzzzDMoLi7Gb37zGwB95xPj8V4WwrIjrGmabTTCyhGWB9AB7gfLmTnCO3fujNlf4zLycmYRAqMjbFU+TUyLZ1Cc6mA5pxbLxmiE/JkZM8JWcRF5GTGPHxlhOsIkIag4wkaRKf7A29vbY3Kv4rZXRkaGfptFnt9rNCKe8mnJygivWbMGQN+td5Hp9VLlQQhd0fJYIJ4HKYSF2P788889CTXVgXICv4SwyN2qCOHhw4ejpKQEmqYFVjlizpw5AIBbbrkl5vO+/fbbUVBQgJUrV+LVV18NZNvpTjQaxXXXXYcVK1Zg+PDheOedd3DXXXdh/PjxOOWUU/Dkk09i3rx5yMnJQUNDA775zW8me5d17r//frz44ovIy8vDSy+9pN99cKKmpkavFz5z5kz861//CnI34yYajeLaa69FNBrFzJkzARxwgzMyMvTfdUtLS79oRDQaxf79+/XzgFk0Ijc3N2bwlFUnNKPgtZouzidmjTZkcSwf780cYbPqClbRCDeD5ezKp6kMopP/d4pGaJrWLxoBHNADxjuuoVDIshutn44wB8uRhOLFEbYqiSbng2UHMBHRCCch7EfVCC9C+Oijj9arPHz88ce278WMIIRwT08Ptm/fDsBeCAsB29bW5rpiBuBuoBwQfxRDID4Tq45yMqFQSK8s8fHHH8e1XTPa29vx/PPPAwBuu+22mNeGDBmCm266CQDw9NNP+75tt/T09Fg6MMnigQcewCuvvIK8vDwsWLAAxxxzTL95ampqMG/ePGRkZODRRx/FY489loQ9jWXJkiV6HGbu3Lm2TrAZ3/rWt3DTTTfprrL4vaYiO3bs0B9Ho1H09PToQri0tFSvBiM7wiNGjNCXkSMUVo6wfJ4ydkLr6emBpmn6dCE2hcCzEsJmHdTk47OVEBbnPDtH2BiNUBks57Z8ml2JNBW3WN4+YO4Im51frTRDvIPl2GKZJA0vjnBeXp7+o5LFrVWbYrsqEypCuKenx3LgiGo0wo86wqpiUNO0mJzqqFGjMGjQIPT09LgWeVZCWIg8L0J4x44diEajyMjIwPDhwy3nKygo0E9MXuIRboWwcI62bdvm+oJBRgwuNH5mVgjnOAgh/PLLL6OtrQ1jx47FSSed1O/1G264AQDw7LPPmtbZTgTbt2/Htddei0GDBiE/Px9VVVX6hVwyiUQiuph8+OGHbZujXHbZZfjFL34BoM9pj6fGdrzs2bMHN998MzRNw9e//nXceOONrtcRCoXw+9//Hscccwy2b9+Om266yXPVmaCRjw3RaBSffPKJLm6HDBmiu7yyEC4rK9PPPXIu12qwnOxCGh1hoO8cYdUhTSUyIc4P8nnQKSPsJhohO7xW7Y/dlE+zE8jy/3ZVI8RFiyAnJ0d/zegIy+dXqwH2XjvLmUUjnDLCjEYQX/HiCIdCIdPcr1nFCEDNEbaLRgDWrnAqOsJffPEF9u3bh1AohHA4rP8PQC/rpYoQuhUVFTHT43GERSyivLxcP/BZEU9OWAhhEbFwYujQofrfmljWLa2trXrOT3zmTgghHMSAuYaGBgDANddcY5qTPu200zBmzBjs27cPL730ku/bd2L9+vWYPHkynnnmGXR1dSEajWLJkiU46aSTsHDhwoTvj0DTNNx+++3Yv38/zjnnHNx6662Oy/zoRz/CJZdcgs7OTtx6662Wt1WD5u6778bmzZtxxBFH4Le//a3n9RQUFOCZZ55Bfn4+Fi5cGFOOMZUwHhs2btxoKYTlzqPiHCIc5ezsbP34ayyfJp+nzIRwd3e3shAWz82ErHyekaMI8rxWjrAcM7CqGqFSDs1YPi3eznJ2LZZlIZyVlaULUKuMMICERiPoCJOE4MURBszFrbG9st28AiFQ7RxheV/c7KO83kRmhNetWwegzwkV+yVEmXArVRHC2SiExfN4hLBdLEIgRKyXNstuHeFQKKQLb6+VI8TnO3z48H5/h1YEFY3o6enBokWLAABf+cpXTOcJhUK45pprAPQ1W0gkXV1duOKKK9DU1ISxY8fivffeQ1NTE84++2y0tbVh6tSpWL16dUL3STB//ny88soryMnJwe9//3ulwZYiGlFWVoaVK1eirq4uAXsay7vvvov/+7//AwA88sgjphf4bpgwYYIupn/yk5/g7bffjnsf/cZ40fr55587OsKFhYX6sVZcuBYUFMQcr80cYT+FsJMjDJgLZDMhLaY7VY2wizRYDXDzu6GGLKyNQtioB8yqMjk5wn5GI5gRJgnBiyMMmMcdzEqnyfO6zQhnZ2frJ8CgHWE/6wgLITxmzBh9mhdHeM+ePfrFhRCIAhEj2LVrl+tb6l6EsBdH2O1gOSD+EmpWURI7hCP8ySef+Hr7+b333kNLSwtKS0v1tthmXHzxxQCAV155JaG1Y3/+859jxYoVKCsrw9KlSzFp0iSEw2G88sormDJlCvbu3Ytrrrkm4bnhaDSK//7v/wYAfPe739XrZ6swYsQIXTj+/Oc/T2j3vp6eHr3k2Y033ohzzz3Xl/V+7Wtfw9VXX42enh5ce+21jreaE43x2KAihIuKivRjrYhGGIWwkyMsD+4yE8JC5BkriXR3d8eITcBc8FpNF9s3E81OVSOMkQanznJey6fZVY2wE8JGR9guI+zGEbYbbO3FEWY0gviKn46wVTTCa/m0UCjkOGDOqY6wqiNsF41wWz5t/fr1AICxY8fq07wIYeG0FBcX9/tMi4uL9c/VbYwgEUK4t7cXW7duBeBNCHtxoIEDjrDRQbdj7NixyMnJQUdHhy/tnQWvvPIKAKCqqso2gnLKKaegtLQUu3fvTpjjt2nTJv1W+9y5c2PiK6ICw9ChQ/Hhhx/igQceSMg+CZ599lmsWLEChYWFqK2tdb38ddddhwsuuABdXV343ve+5/8OWvDQQw/h/fffR1lZGX75y1/6tt5QKIQ5c+ZgzJgxWL9+Pf73f//Xt3X7gTg2iGPU5s2bdSE8dOhQXQjv3r075hxhFML5+fmWjrDsQorp+fn5ulEiV5NQcYSN0+wcYU3TYqa7EcJmjrCbaITKdK9VI4wZ4YyMDEtH2M1gOWNGWFQGsULuV+C0DeMydISJL6gIYTORaSZunQbLuS2fBjhXjlAtn6bqCNtFIzo7O5VqzQoRJ0QdcECYuRHC4gRjlrGVYwRuRWMihPCOHTvQ29uLjIyMmBHiTvgVjXAjhLOysvTBh37GIxYvXgwAOP/88x23L5quzJ8/37ft23HPPfegs7MTZ555Jq644op+rw8bNkyvYTtr1qyYygBBommaPkDuzjvv1CsOuCEUCuGhhx5CdnY25s+fj5dfftnnvezP+vXr8dOf/hRAX2v1YcOG+br+4uJi/ft48MEHXY81CBJxIX7yyScDsHaE5TKMshB2ikYYHWF5UJbZdCshbCaaBXaOsErVCTHdqWqEMQLh1FnObn6v0QizjLCotKGSEbY6J1tFIwD7eIRZNIKD5UhCUYlGxOsIe60aIW87URlhs2iEXA5OxRUWpY7k9r6yI6xak9dOCMvT3YrURAhhcXIsLy+PqZ3pRDKiEYD/lSO6urqwfPlyAMCUKVMc57/ooosAAP/4xz982b4dW7Zs0cu1zZo1yzJ/e/311+PEE09EW1sbZs2aFfh+AcCCBQuwevVqFBYW4s477/S8nnHjxunLf+973wusWQpwYGBfe3s7pkyZgltuuSWQ7Vx22WWoqqpCZ2en3owlFRDHBichLH7Tos68ajTCKiOck5MTI9yEEBPmh7F8miyQVR3hrq4uywiAm2iEG0c43kYb8v+q0QhxjDZmc83MMFVHODMzUz8v2wlhN66z3TKpBIVwmuHVETYTt06OsNybXaAqhK1OZH5khKPRqJ6zNXOEMzIy9HiESgm1bdu2AYitlXn44YcjFAph3759Md2L7HCquuB1IFsihbBqxQhBMqIRgP+VI95//310dnairKxMSZRfcMEFCIVCeP/99z1XzFDld7/7HXp6enDGGWfg1FNPtZwvFArhvvvu05fx+p24QbieX/va1/Q7MV65++67MXLkSDQ1NQXaoa2hoQHz589HTk4O5syZozSwzwvC6c7KysILL7yQsLsHMvv27cPkyZNx5ZVXQtM0aJrWTwjL0QhZCIu/H2EsuIlGuBHCVo6wPN14vjOrIyyeexHCsiMsGx9Wg+Xclk+zE87idcC8aoRZgw6jI6xSPk3+XOSKGXLEQaWEmlNDDTPjiELYhEgkgurqapSWlmLixImeRwu7XY9f200mVsF3wL0jLESiEI0Cu9Gj8UYj/HCE5cFmZkIYcDdgTgjh8vLymP0QwlMMpnPCabCZ19JmXoSwPHBPBS8D5YDYphpuu9l1dHTo23XrCPtdOUJkfU855RQlYTR06FBMmjQJAALtJNbV1aV3ulPJz1ZXV2PKlCno6uoK/Pi2evVqLFy4EBkZGfjOd74T9/oKCwv1fPN9993nuVW4Ha2trfq+3nXXXUrdDOPh6KOP1r+3O+64wzHy5TcPPvgg3nzzTTz77LN4/fXX0dLSou+DqJPd0tKi10sfMmRITLQMOGCUiGOtmSMsx9CsBstlZ2d7FsJu4g5WzqeqIywPwLUaLGdVVs1r+TQrsS3/L2eEnRxhJ7dW/jxlQatSOcIuGgH0H+woL8NoxP8jEolg4sSJCIfDaG5uRn19PebMmYOamppA1+PXdpONnxlhIYRVG2pomuboCMsHRiOapvniCIv9CoVCloLcjRAW0QhZCANwnekNIhrR1dWln3hUhPDgwYP19+5mO25LpwnEZ7Rv3z7X3ezEBUZRUZHrbKnfjvA777wD4IBLpsLZZ58NIFghPH/+fOzcuRMjR460LOkmEwqF9AoOjz/+OJqbmwPbt9/85jcAgK9+9asxA03j4brrrsPpp5+O9vZ2/OAHP/BlnTI//vGPsW3bNowbNw4/+tGPfF+/GT/72c9w6KGHIhKJ4Oc//3lCtilYsmSJ/ri+vl4/JgwZMgQjRozQY3HydOEIC8Q8dtEI4MAdRKuGGk7RiFRxhOXqFHbRiKDKp7mJRnh1hOXPwq0QNmuoIT92O5A/FUi4EBbCU7gc4XAYs2fPRkNDAxobGwNbj1/bTTZ+ZoSdHGHj7ZGenh79B+wlIyzfdorHEZbzwVbunaoQ3rdvn74+4yAxtwPBVIWwm1vWwq3Ozs5WFoteBLdXIVxQUKDvl9tb8XLzEbe3p4888kgAfSfl3bt3u1rWDNkRVuWss84CEKwQFtngG264QTm7fe655+LYY49Fe3s75s6dG8h+ffHFF/jjH/8IAHFlg42EQiE8/PDDyMzMxLx583xtEvL222/j97//PYC+84DVMcxvBg8erNcqrqurw5tvvpmQ7QKxF4orVqzod4wy/t5VhLD4vcnRCODA+cTKEXaKRrjJCFvVBo5XCJs5wm6iEX5khO2iEV4Gyzk5wrKb64cjbFZLmI6wRGNjIxobG3HVVVfFTJ86dSoAKJfecbsev7abCgSRETYKYStHWHZpvUQj5GlWJyEVR9huoJxAtYSacIMLCgr6xSz8doS9RCNESbORI0fqB0snEimEAe8D5ty2VpYZPHiw/j4/+eQT18vL7Ny5U98Xs7bKVpx++unIzMzEunXrXLfiVmHXrl16BYWbbrpJeblQKKQPzvrtb38bSF3hRx55BJ2dnZg0aRImT57s67qPO+44fPvb3wYAfPvb346rfbegs7MTX/va16BpGm6++Wb9IiZRfPWrX8XNN9+MaDSKG2+8MSG1hb/44gs9+wsAK1eu1I9l4nduPFaVlZX1y3qLO4bG421BQUFM7XhxPrGqDuFnRtjogArcRCPkhhpyNELVEXZyiu2iEUbXN97BcioNNawcYfkCW1z02MXq7FosA+ZCmI6whHBjq6ur+70WDoexePFipVvZbtcTz3Y7Ozv1vKXb3GUQJNMRll1aqz9oOyGssryKI2xXQ1ig6gjL+WCjK+lmIFhHR4d+0nFyhN001XCTDzZux4sQdjtYDohfCLsdKCfwKycsYhFHHXWUqwFfhYWFeuONIFzhv/3tb+ju7sYJJ5yAY445xtWy1157LcrLy7FlyxbU19f7ul+dnZ14+OGHAfS5wUEMNrv33ntRXl6OtWvX+tKq+L777sOaNWswfPhwX2sGu+Ghhx7CYYcdhkgkkpAqEmvWrAHQdwGen5+PtrY2/e9U/M7l3/uQIUOQl5eHzMzMmLic0REWFBQUIBQK6cdscT4JYrCcUVzZubxOjrDY366uLtNohKojbFU+LR5H2Cie5f/NMsJuGmqYOcLyhQygdt40c3dDoZC+T4xGOLBs2TIAB0pTyYgvQKXeotv1xLPdWbNmobi4WP9n7BiWaLw21LDLCFsJ4f3798dcHctC2+rkp+IIywcXI24cYT+EsFU+GHDnCAshWVBQYCmmSkpKXDfVEELYjVPrRQh7HSwHeK8lLEcjvOBXTthLLEIQZE74r3/9KwDgxhtvdL1sbm6u7qr+6le/cj2Q0Y5nnnkG27dvx6GHHhrYGIvi4mJdsN53332uW53LrFy5Um9q8bvf/c5TrWM/KC4uxlNPPYVQKIQ//OEPejvvoBC/i2OPPRbHHXccAOh3GMwcYbmzphyPEI+NQlgcq41COC8vz3KwnDy4y0tGWBaxgLeMsHgfZkLYzhFWEbxehbBVG2f5fzeOsGpG2BhVkLsKWmGlM+zaLDMaISHEpplQKCsri5nHz/XEs9277roLra2t+r9ElCSyw65MiUqLZSEiNU2zFMKywJTjEUKc2mXr7MqnqVwVus0IWyG+a6cBXGal0wRuhLAci7C6SAiFQq5zwl4cYbcRjD179uifaTzRCLe/jXiiEYB/tYS9DJQTyELYT7G5bds2vPHGGwAORLjcMmPGDOTl5aGxsVFfV7xomqaXNvv2t78dc0vUb6699lqcc8452L9/P26++WZP7azFsj09Pbjyyitx5ZVXBrCn6px11lm44447APRVATETDX4hhPD48eNxwgknADhw7DTLCIs28ECsEBb11c0cYeDAMVuYDnl5eaad5WRH2KyhhnBjhRA0c4Tl7me9vb22jrAQmcYKE2IdVkLYqmqEXTTCa91hvwbLuc0IxyOErZa1a6pBR1hCJfagIoTdriee7ebm5qKoqCjmXzKR/5CsSsqoRCM6Ojr0H5XxPeXl5eliThbC4na+nQBVcYTtfgxuqkb4HY0wIgTl1q1bHU9YTvlggVu3NhHRCLk1tN13a4WXaERPT4/e2jpeRzgeIRyNRnUh7MURnjx5MrKzs7Fp0yZfu4e98MIL0DQNJ510kue7UEOHDtXdZFHhIV7+/e9/4/3330dBQQGmT5/uyzqtCIVCeOyxx1BYWIg33njDU0TizjvvxAcffIBhw4bpA9aSzT333IOhQ4dizZo1eOSRR3xbbzQa1e9wAeZCWCB+OyqOsDg2Go+3bh1hlWiEfJw1c4Tl41N3d7dtRlgWvPL/ZkLYLiPsplOc285yxjrCqi2WvWSEraIRMsIYtBLCdpWf7BxhCmGXxFuU3et6/Npu0NiVKXEzWE5knUOhUL8DnDzNzBEWBygz7Mqn+e0IBx2NGD58OLKzs6Fpmi5IrVAVwm4H4MUjhFW3Ec9AOcBbNGLjxo3o6elBbm6u5+2Kk3kkEvE8oOqTTz7Bnj17UFBQ4DqHC/T9FoST7Gc84tlnnwUA03bKbvjud78LAHjuuef0C494EG7wzTffrJ80g2TMmDF46KGHAPQ13HDzGf/xj3/EI488glAohD/96U+mv/FkUFpail/84hcA+kqryQPavNLb24sLL7wQhxxyCP7whz8AiBXClZWVMfNPmDAh5n8AOP7442P2USA+NzeOsJvBcmK9dkLY6AiLdds5wiLb3NXVFePOinXIg+WsDKZklE9TabEstu0mI+wmGmFViUeuEezGEWY0QkJkdO3EicrB1e16/NpuKmAnhN0MlpNrCJtVIzDLFAtH2KpihLzteB3hnp4e08LcwAFx7kfVCLtoREZGhrKodOsIJ0IIt7S0KI1OF/vi1XkUjvDmzZuVb1+LfHA4HFauhmFk5MiRKCwsRG9vr+cMqcgHT5o0yVVraRm/c8LNzc149dVXAQCXX355XOuaMGECqqurEY1G43ZEP/74Yz1jKgR2Irjllltw/fXX6/GGtWvXOi6zYMEC3HbbbQD6agefd955Qe+mK6ZNm4bjjjsOzc3NmDZtWtyVPV566SUsXLgQ0WgU3/ve97B9+3b9mDR+/PiYi7xQKKQfH8eOHYvGxka8/PLLuO666/R5vAhhcVx2Wz5Nzgi7cYTldYvfrpUQls9HdtEIsT3AudlFUOXT/KgaYZYRNnOE3UYj5HUYz+NGUS5DR1hCXJWa3UIU06qqqnxfj1/bTQXkH5obR9gobK3ywQIzR1gIYTtHWKVqhIojDFjHI/x0hO2iEYC6g6vqqgqR6nawnBshXFRUpJ8EVLajKuKtGDlyJDIzM9HT06N/nk7EWzEC6DuhxxuPiCcfLDjnnHMAAK+++qovOeF//OMf6OnpwTHHHINx48bFvT7R2ezRRx81bZuuyoMPPghN03DZZZfpdZwTQSgUwqOPPoqTTz4Zzc3NOPfcc/Hpp59azr9gwQJcccUV6OnpwbXXXpvwJhYqZGZm4pFHHkFOTg6ee+45nHPOOcq/HTPkgXft7e16h74RI0agtLQUeXl5OP300wGg3+dxwgkn4OKLL465EJSdYishLASs0RiRHWGVhhpyRlgWUbJAlgfLCWEoC1wR75PFsTgGapoWcy5xEsJiebPYglNnOdXyaSpVI9y0WFbJCJs5wsZohJMQltdhFNF2A/npCEtcffXVAKwFaWVlpVJEwe16/NpuqmD1B+fFEbYSwnaOsFchbCfUBfJrVvGIREUjAHUhHERGeP/+/fotKjdC2O12xDxeHeHMzEx9e6rxCOEIex0oJ4hXCMdTMUJwyimnIDc3F9u2bYu7pjHgXyxCcMEFF2DcuHFobW3FU0895WkdW7Zs0RtoJKPuel5eHp5//nkcddRR2LRpE0466ST85S9/iRnY1NHRgXvvvRcXX3wx2tvbcf755+PJJ5/0fMchaE499VQ8//zzKC4uxhtvvIGJEyfaCnw7xB2Eo48+GgD0ihviOQC8+OKLeP/993H33Xc7ru8b3/gGJk2ahGOOOUa/4+PkCAvsqkbYCeFoNKqfI7KyskznlStPmIleM0cYiDV07KpGiOmAc2mzeMqnWbm+ThUlVFssy9+Jl8FyVtEI+bMx3kGzcoRlAU9HGH2jn8PhsF7XV9DQ0ACgr0KDTEtLC+rq6rB48eK41uN2/lTHSgirZIQ7OzvR3d1t2V5ZYNZhxo0Q9lo1IiMjQ39/Vo6wSjRCRQhrmmYbjQD8d4TdVHQQzTTy8/MtL1iscCOExXvz6ggD7rPP4vb2l770Jc/bBOITwm1tbVi1ahWA+IRwXl4eTjvtNADxxyPa29vxyiuvAIg/FiHIyMjQowwPPfRQjHhU5aGHHkJXVxdOP/10nHrqqb7sl1vKy8vx2muv4ZRTTkFrayuuv/56HHnkkbj11ltxzTXXYPTo0bjnnnvQ29uLm266CS+++GLKOlCCCy+8EO+++y7Gjx+PLVu24MYbb7SMhFmxZcsWfPzxx3pZNhnxdwn0iRxRRs2J/Px8vP322/jggw90gWM83oo4oaoQdnKEgQPHfKNoltchr9voCFsJYWGeZGVlxUQF/IpGWJVPU8kOy//bNdvws8Wy1WC54cOHAwB27NhhenfLrJmGwEmXGPcrlUj4pXJ9fT0ikYje+jMSiaC2thbTp0/vVyZo7ty5qK2tNa1V6WY9XuZPZbw4wrLgbW1ttewqJzBrwOGXI+z0YxC3xPxwhDs6Oizzd62trfo+WQlhlYoI3d3duqBWdYR37txpOyAQiI1FuG1a4EZwxxuNANxXjkgFIbxs2TJEo1GMGjXKteNuRI5HxMOiRYvQ0dGBMWPGxAxeipebbroJJSUl+Oyzz3QDQJXW1la9usHMmTN92ycvDB8+HEuXLsXPf/5zFBYW4rPPPsOTTz6Jv/3tb9i1axcOP/xw/PWvf8WTTz6Z8iJYMG7cOCxcuBDFxcV499138atf/cp2/l27duHBBx/EypUrARy4+KqsrMRpp50WUyFC5Ne9kJmZGXPcMR5vRT1mN0LYro4wcOAcY+Ueu3GEZeEuzhlGIS1Ep/werKIRfnaW86vFsvx5ys1AvJZPE3dGu7q6TE0ku4iDlSMsbzdVf5MJF8KVlZVYvnw55syZg4qKCtTU1KC2trafWwv05XZLSkr6tUZ2ux4v86cyXhzhrKws/WDR0tKiHI1w6wjHWzVCXkc8GWGj8DdDxCKKi4stBwCqOJ3btm2DpmnIzs7GsGHDLOcD+lwZsS2n/K6XfLDAzaC8RAvh3t5ePabkZzTCbT7Xj1iEQAiOf//7354cV8Hzzz8PoK8lr58d2wYPHow777wTQF/XNjc1eR988EHs2bMHRx99NC6++GLf9skr2dnZ+O///m98/vnnaGhowP3334/Zs2dj0aJF+Oyzz3DNNdcE0u0uSEaNGqUL4J/97GdoamrC9u3bUV9fj3Xr1unzdXd344wzzsB//dd/4eSTT8by5cv1iy9xMfad73wHQF/N4ilTpvi2j0ZHeOjQoQD6n2+Mg+XMqkZ0dnb2qxcMHDjHGKMRbhxhMa+8HyIbLy8vn1/MBqE7OcJiuuzge2m0If9vF41wcoStOrfaOcJGYZqXl6ebSOKOpIzdOVzFEQ6y7ng8eBsmHSdClKrMZ1fYWXU9XudPVbw4wkCfS7p37140Nzd7GiwnDhxBVo2Q12/lmKpEI0Sb0D179qClpcVUoDrFIgA1ISwE7ciRIx3ziKKpxtq1a/H555/bDhbzQwg7OcJtbW36byyeroluSqht2rQJXV1dyMnJ0QW0VyoqKpCZmYm9e/di69atrj4rPwbKCU488UQUFBRg586dWL16NY499ljX6+jp6cFLL70EoE8I+813v/td/PrXv8aaNWswb948XHvttY7LbN26VRdo9913X0rlbYuKipLeIMNPbr31Vvz5z3/Gq6++ihNPPBHt7e3o7OxEbm4uXnzxRZx33nl48skn9bJoXV1d+Na3vqUfx4QQvuWWW1BVVYVDDjnE1+/L2A1UCCYzR1gWaGIfZCEsn1fk84lVNMKLIyyEcFdXl6kjLLe5z8nJQUZGBqLRqO1gOTOH16rcmtcWy3b5YbuGGlYRBDtH2EyYjhw5Ei0tLdi6dWtMxlxezosjnJOTk7IXqKlzVCPKeHGEgdjcbDIHyyXCEQacS6g5VYwADgi8nTt3Wu6PEMKqjqqqSE2EEBavFxYWxtUsRgjaDRs2OM4rl06zarWtSm5url4ecfXq1crLaZrmqyOck5ODM844A4D3eMQbb7yBXbt2YciQIZg8eXLc+2SkuLgYP/zhDwEA//3f/61Ue/mee+5Be3s7Tj311EDEOTmAqI4xatQoNDc3o7OzE/n5+ejs7MR1112HDRs24L777gPQF1EpLCzEu+++i40bNyI/P1+vCgH0/f6DvmgRv10zISzncMXfmSxujSJUrMsqGuElIyw7wrIQNhPjWVlZurg0ZoTNRKrs8MrCz2tDDbsWy3YZYbPPIjMzM2Ygm5sWy8CBToJ+O8Kpmg8GKITTEq+OsBgR2tLSomeErcRPssqnAf44woDzgDmnihFA32cm3q+VqBTTVRtDpKIQjicWARyo1R2JRBwjCiIfHG8sQiDcVzHwTYXPP/8cW7duRVZWVr9mA16Jt56wiEVceumlnmsaO/Gd73wHI0eORFNTk2OntrffflsffFVXV5eybs7BxNixY7Fq1SrMmzcP7777Lpqbm/HlL38Zu3btwpgxY7Bx40aMHDkS99xzT0wZtBtvvNHRGPADs2O3akY4NzfXVAibCeR4HGEz0WzmCMvntuzs7JiybIC9SFVxhOPNCLttsezU8U2lsxxwQAiblfPz4ghTCJNAiNcRVolGxOsI21WNsCufJr8eryPsJIRVHOFQKOQYj3DbmU1VpMqRC7eIbezatcu2XbUfFSOAvhM40Nex0KlTll8D5QRiUNn777+vvIxwg4877jjbv2c3CCH82muvucrgAn0OtRDCl112mS/7Y0ZhYaFeWusXv/iF5WfW1taGW2+9FZqm4aabbopxG0mwlJSUoKamBieeeCJyc3Pxl7/8JeaYeffddyM/Px/f/va3cccdd+Bb3/qWXjc4aMyqTtgJYVm8WQlh2Y316gjbRSMAZyEsO7luy6cZHeFkNdSwMprcOsLifGjmCMsXNUasWiyneg1hgEI4LTETwlYjRmVkYbhz504ABwY8GIm3fFqyq0YA6kLYLiMMOOeEg4pGuF2vTElJif492Q3Ki7eGsCA/P1+/EHDq8iaiEX4L4Q8++EB5GZEP9iMWIaisrERRURFaWlrQ2NjoatkPPvgA69evR35+fuBd0K655hpccskl6OzsxFVXXaUfCwS9vb245ZZb8PHHH6O8vFxvq0ySw/jx4/Hcc8/h7LPPxk9/+lN885vfBNAnhn7729/i4YcfjivW5IYnnngCxcXFuP/++/VpdoPljNONQjgrKwuhUEifrpIRlgfcielWQlica+wywmIfhLh06wirDJZzqiNsF5uQX7NrsWxlNHl1hO2iEXaOMKMRJCGYdYqRRaNKRnjHjh0AYFnlwGv5tERUjXAbjXCqGmHnCAPqQtitI2wnUDVNcx25kJGdbDvB7Vc0AoiNR9jhdzRCuFRr1qwxbe9phnCE/RgoJ8jKysK5554LAJg/f76rZYUbfN555/nmUFsRCoXw5JNPYvTo0fj0009xzjnn6OXndu3ahauvvhoNDQ3Izs5GQ0ND2rSfP5i54IIL8Oqrr+Lee+9NakRlwoQJ2L17N3784x/r0+wcYRk5nyvOJWI+JyFs5gjLJo1ZRtgpGiG7zwAcHWHVaIRfneXcOsJW51c/M8JeHGEKYRIIZld4KkWrxQlt165d+OKLLwBYC+FUdYSj0ai+T4mIRgDOQjiIjHBra6v+eXsRwvJ27Cpe+CmERQUMO0dYLp3mlyM8ZswYFBYWoqurS6mzW2dnJ5YtWwYAvjeHuPTSSwEAL7/8sqvl/v73vwMINhYhM2TIECxatAjl5eVYtWoVjjnmGEyYMAGHHXYY/v73vyMnJwd//vOfAxm0R9Ib4yA8WQiLgVpmTqOZEBbPvWSEZZNGOMLGTLJKNMJKCJsNVpMHyyUjGmHXYtmvjLA4H5plhL04woxGkEAw+8MWojErK8tyJL6IAGzZskXPcYpOMkZSoXyamSMsZ8sSFY0QFRHMBKWmaZ6jEdu2bbN0MIVALSsrs/28VbZjJ7j9yggDBxxhOyG8ceNGvXRavHEMQSgU0l1hlZxwY2MjOjs7MWzYMN9cacFFF12EUCiEZcuW6YMdnVi1ahVWrVqFnJychFZmOPLII/HOO+/gkksuQW9vL9asWYP29nYcf/zxePXVV00bGRFiRBbC4rGxBa8onWXlxnrJCIvawIC7jLCVGDcKWzMhKrByhEOhUIxwNpvfqY2yXUUJvzLCVnWEATVH2K6zHB1hkhDsHGG7PzZxpbdmzRp9ZL/oDmQkFcqnmTnCQpiHQiFHgWhXPi0ajerxkHgc4ebmZn0/Vas7DB06FDk5OdA0zfRgA8SXDxa4EcJ+iFLhCNtFI9asWQOgr5tWvKXTZNzkhN944w0AwOTJk32/zTxixAicdNJJANTjEX/9618B9IloUdklURx22GF46aWXsGHDBsyfPx/vv/8+Ghsb6QQTZWQhLB+TZbEkjvlOjrCbOsJy22SxXSshbNZQQ5xLhKg0ZoSNbq1VdQgx3ShcVesLy/87xSZUqkaoZISd6ggDfXcljWaU3TmcjjBJKHaOsIoQFkKlrKzMskxTvOXTzESs2/JpZo6wOPgNGjTIUcTYOcK7d+/WDypWrrjATggLkTl06FDlK96MjAw97uBXSTYznIRwc3Oz3kxjzJgxnrcjEEJYDIYzQ9T6nTBhQtzbk3HjCP/nP/8BgMDE3iWXXAIAenMMOzRNw1/+8hcAUGpwERSHHXYYLrzwQhx33HEp1TSDpD6y8JLz5LLwMWaBjbEEp2iEnSOck5MTY8DYlU+ThbRqNMKsg5ydEDYbROdURzjeqhF20Qi5prPAzhGWO60ajRo6wiRlsHOE7UqTGZ1POwGY6o6wSs1MOyEsYhFDhgxxbPsohHBra6tef1ngdqCcwEmk+pHdddqGaN06fPhwX2qQjhs3DkBf9MZqgGJQQlg4witXrrStY6xpGt58800AwQlhkRNeuHBhv78XI2+99RY2bNiAwYMH6wKakHTCSgibtfn1khGW2zSbOcJGIaxaPk01GmElbN3O79YRNgpxpxbLThnhzs5O/dho5wiHQiH97qYx3uXFEaYQJoHg1RE2ZmFF7VczvA6WU6ka4VRHWNURdsJOCKtWjAD6RLdYl9EVXr9+PQD3jqqTSPUqsN1sQ9wZENneeCktLdVvq4k2sEaCdISzsrKwfft228GBa9euxRdffIHc3FzfGmmY7ctRRx2F/fv364PgrHjssccAAJdffnng1SIICQK5dJuVI+wkhOPJCMtl0twIYaOotCqfZtVBzih4VR1kP1ssyy6sU0ZY3ic7RxiApRCOt8VyqkIhnIaY3epQEZn5+fkxGUS7Ufty+TRxFZkKjrBqDWHAvnyaasUIgVU8QohJu4sKMxLhCIt93rFjh+n34bcQBg4IXLN2x9FoVBfIfgvhgoICfPnLXwYA3fE1Q+SDRbOCIAiFQrjxxhsBAE8//bTlfLt379ZjEd/4xjcC2RdCgkY+RsnnlyAywlYl0cyEsFP5NIFT+TSnzK+VI2yXEbYTu14ywrIjbJURFp+P/B6t7oYKQ8NKCLPFMkk6Xh1h4MAtZMBeCAvHVdM07N+/H5qm6QcpFSEsX70K/Kga4SUasXv37n6vqVaMEFgJYREvcCsmE+EIl5WV6QdFs5rFQQjho48+GsCBQXEy69evR3t7O3JycvQ8sZ+IUmhvvfWW5Tz//ve/AQQXixBcf/31+vasSro9/vjj2L9/P44//njfy7gRkijkY5Q8bsMuIyzODW4zwnI0wsoRlsWaWUMNo/hTLZ9mVS84noywX9EIlfJp8ntTdYSNGWGV8ml0hElC8JoRBhBzO3j8+PGW88nRg3379sW4syrl0+R9Mj73wxFWiUaIGsnt7e0xZdcAd9EIILGOsKZpeuQinmoOoVDItpZwoh1hUdFh/PjxloM04+G0004DYO0Ia5qGhQsXAgCqq6t9377M4Ycfjq985SsAoLc0lmlra9Pb4t5xxx1JbZJASDzI5xz5OGvnCAvc1hF24wg7RSOM+yCOSVbl05wGy/mVEfbaYtnq/JqVlaWvS9URtqolTEeYpAxey6cB0E/O559/Ps466yzL+TIzM/UDXFtbW8wBTlUIG4WsH46wm2hEYWGh/lmJBiICv6IRwhH2KoQ3btzY77Xm5mZ9kJXb9RoRzuunn37a77VEO8LvvfceAGDSpEm+bU9GCOGVK1fGZNsFH3zwAbZt24aCggKcfvrpgeyDzMyZMwEATz31FDZs2BDz2m9/+1vs2LED4XAYN910U+D7QkgicBLCRhEqpqtmhJ0c4d7eXn0dshA2i0sIxLadHGHVqhF+lU+zqyNs1mLZ7q6wUTPYNdQADphIxvOmF0eYQpgEgl00wskRPvPMM7Fx40bMnz/fsY6rXEJNbolp5+aZXX0a99GPqhEqjnAoFNIrYwQhhJubm/WBeG4FqxCfmzdv7if4hUAtLy+PewCVcP2Ng9d6enp0cRaEI7xp0ya9NJvgnXfeAQC9zq7fjB49GmPHjkVPTw+WLFnS7/UFCxYAAM4666yEHJQnT56Ms846C11dXZgxY4aetf/www9x7733AgDuuecex6olhKQ6Imp022236dNUHGGVFssinifmt3OEAfOawfL23GaEExmNMGuxbJcRNnOEzTSAceyOU1xBnDdFrX2BSvk01hEmCSEeRxjoEwwqtULlEmoqA+UEVpUjEu0IAweubI0/aBGNUM0Im3WXE27wiBEjXAvWIUOG6A0/xHoEfjq1Rx11FID+Qvizzz5DT08PCgoK4sohGyktLdWz53JWNxqN6o5wUEI4FArh4osvBgD84x//6Pf6Cy+8AAC44IILAtm+GY888ghyc3OxYMEC/OAHP8C///1vXH755ejs7MRFF12EG264IWH7QkhQ/POf/8Q777yDqVOn6tPsMsICu0F08rzydKMQlh1h4IBZIjvC8n44RSOs6gi7HSznpo6w7Pq6iUaoZITF+5bfW7yOsF35NDrCJCHE4wi7QS6h5kYIW1WOUM0x+1VHGICjI+xlsJxw9uIRrKFQSG/va2xA4XUAnhnCEf74449jpq9atQpAn4PrdwMFs6zup59+ij179iA/P9/3ihEyshCW6wlHIhG8+eabCIVCuPLKKwPbvpEjjzwSDz30EADg17/+Nc4++2x89tlnGDNmDB5//HFmg8lBQWFhIU466STLwXKqjrCVEBbHfStH2OwupZUQVh0sZ+UIWwleu/JpcuvleAbLCZHsJiMsT1N1hK2EsBdHmEKYBEK8jrAqcgk1P4WwX53lVDBzhLu6uvTnqm6oyPR2dHRg165dAA4IWK+C1aoTm9cBeGYIIbxhw4aY/J4Qwscee2zc2zAibpOKUmUA8PbbbwPoG6wZZBTgrLPOQkFBATZv3qw70AD0MmXnnHOOcitsv5gxYwbmzZuHY489FkOHDsX111+PN998U/kijJB0REUIq2SEjdPFeoVgzMnJQSgUMq2WoOIIO2WEjcJWVSBbTffaYtmsjnDQGeGOjo6Y8RZeHGFGI0ggJNMRthsoJ7ASwqr76FcdYcDcERa1EXNzczFkyBCl9eTm5uqiWQw8E5URxAAxtwhHuKmpKWa6n9GIoUOHoqysDJqmxZTxSoQQfuedd/SD4qJFiwAAU6ZM8X17Mnl5ebjssssAAH/4wx8A9J3AHn/8cQAHypolmpqaGnzwwQf44osv8Kc//Umv00nIwYpdZzmB0REWjmd2dnaMy2vmCAvE+cKskYSXaIST82tVX9iq3JpKZYh4q0bIn5EfGeHBgwfry8jnTrZYJilDMhxhUcVA7iJkRbxC2K86woD5LR65Rq+bW9NCNIoyYB9++CEA4JhjjlFeh4xVNMJPIRwKhXRXWK7kEKQQPuqoo1BWVoaOjg689dZb6Onp0cuWJSKfK5pTPP3001i3bh2eeOIJrFu3DkOHDsVVV10V+PYJId4ywoLs7GyEQiFHpxjwTwh7jUY4lVszCuQg6ggDsZlpI24d4VAopJtEch1+u/JpVi2W6QiTQIinoYYbhNjcu3evJyEsO7q9vb36j8/JVfarjjBgHo0QtXvdDhITonHVqlXo7u6Ou0uaiEbIjnB3d7deUs2PaARwoIlKY2MjgL6LCSG2jzvuOF+2IZORkYFLL70UADBv3jwsXLgQO3fuxLBhwxLSOOKMM87A2Wefja6uLkyZMgXf/e53AQA//vGPlf9uCCHxUVhYqD9WdYQF4rmKkBXrNotGGKfl5eV5jkY4Ob9uHGE3VSNUHWHRQdUPRxg40JBKVEaSl6cjTJKOnRBWiS6oIrcoFkJYVDqww8wRlh/H4wj7EY0QjrDb9sWihe+7776LFStWoLOzE6WlpZ4Fq3CE169fr38+a9euRU9PDwYPHhxXe2UZUaVBlC9bsWIFNE1DeXm5fqHgN9dddx0A4NFHH8W0adP0aYkoFRYKhfDYY49h2LBh+Pzzz9HR0YELL7wQd9xxR+DbJoT0IbdbFo/NhC2AfoPdjI025OkqQjgjIwNZWVn95s3Pz3d0hN06v07l01SEsF3VCPHZODnC4hxtNo7HrSMMHPjOZCFsJ6CtHGEKYRIIZkJYiMYghHBLS4t+taniCJuVT5NFrWpGuKenp1+bZq/RCD8c4TPOOANAn5AUNWlPO+00z1UXRo4ciZKSEvT29urussgdT5gwwbeKAqeccgoAYPny5ejo6MBrr70GAIE2laiursakSZPQ2dmJLVu2IC8vDz/4wQ8C256RsWPHYuXKlairq8OTTz6JF154IZBudoQQc2QhLG6zq0YjzIQeoJ4RtmrgkZ+fb+k+O5VPcztYzo0QtotGiP2Sz4dimhwfEedoMw0QjyMs14NXKdHGwXIkIRi75QDBDJaThXC8GWGxf5mZmY6CRP4hG+MRbqMRfjrCo0ePxpgxY9Db24uf/vSnAOIb/BUKhfTYgl+5YzPGjRuHQw89FPv378fSpUv1gWtnnnmmb9swEgqF8Le//Q0nn3wywuEwnnnmmbjaRXvhkEMOwX/913/h5ptvZtMKQhKMihB2ikZ4dYSthHBeXl6MeJS3IYSnOG+JeZyErVNnuXirRtgJYXn/7YSwF0fYLBpBR5ikDIlyhOVbI+JH5jUa4Uaoy/MY4xFeHeH29nZ9Wa+OMABcc801Mc+vvvpq1+uQERndFStWAACWLVsG4ECu1w9CoRAuuugiAMAvf/lLLF26FABwySWX+LYNM8LhMN5++200NTXplRwIIQMDWQiXlZUBUHeE7TLFKo6wmMcYExDnR3kdxoywOG+pZn6tBLKKI2w2kM4qGiELYbkrrHgv4hxtFo0Qn41RCPuZEXZqqEFHmPhKosqnBeEIqwj1jIwM/T1aOcKqQnjw4MH6vMIJFkLYSwb39ttvx9ChQwEAX//613H44Ye7XofMySefDKCv5q6maXq9XRFn8AvR+nTRokXQNA1nnHEGxowZ4+s2CCFE4MYRNt4ltKsm4cYRNopCcX6U1+G2s5xxuiCejLDZYDmVaIT8nkQDITtH2BiNcOsI253HnVos0xEmvpLojHBzc7OnwXKyiHUr1MV8siMcjUZ1V1c1GhEKhfT2yBs3bkQ0GtXrCHtxhEeNGoUPP/wQ8+fPxyOPPOJ6eSMid9zY2Ih33nkHu3fvRm5urq+OMNAnrEX709zcXDzwwAO+rp8QQmTE+QNwnxG2m+4mI6ziCFtFI1QdYYFV+TQ3tYLtBtJZCWGr9yhjFY2wc2nFhYycEbbTGWyxTBJKMhxhN4PlzBxh8QNS3T/xQ5PFtHCDgdjSPE4I13bDhg3YsmULenp6kJmZifLycuV1yIwYMQIXXnhhP0fAC4cddhiOOOII9Pb26q5tVVVVILeR/vKXv+DZZ5/F6tWrdSeaEEKCQL5r5zYj7EYImznC8UQjjI6wU5k0gRvh7Kahhlw1wthi2e49ylgNlvPqCJudx50cYUYjiK+kY9UIPxzhvXv3Aug74LgR/LIjvHbtWgB9+dVUGUB17bXXAoBeOeKKK64IZDvZ2dm4/PLL9frFhBASFHIrc3EuURXCZtOzsrKQmZnpyhE23jkU88rrNcYzVKtGeI1G2L3mtmoE0F8IO5VP0zTNc0bYiyMchEnnNxTCaUiiHGExwGHv3r16vlal7my8GWF5PtkRFkK4sLDQVWkxIYTXrVunt0ceN26c8vJBc/vtt+vVLcaPH48bbrghyXtECCHxUVhYiI8//hiRSCRGPMoiVJwrVDLCZvleeV5ZCIt5VRxhIZZVoxFWjrBqNMLuNS9VI9w6wr29vXqe2I0j3NPTozvSdo5wT0+Pvn6xTatlUgUW1kxDjCNAgWAc4bKyMmRnZ6O7u1t3hEeMGKG8f0E5wiqutMyRRx4JoM9xFUI+lYTwiBEjsHTpUrzxxhu4/PLLU/oWEiGEqCKOvTKDBw/Wc6duohFWJdHMpovonMpgOaMQFoLXGJlwygi7iUZYCWE5/qAqhI3nfKeMsOzYuhHC8rnYzhEG+t67ccA7hTDxFeNtjlAoFMgfWygUQnl5OTZt2gSg74cpKibY4YcQdnKE3SBq8q5evVo/6B199NGu1hE0Rx55pOlJgxBCDiZkIexmsJyVI2wWjRD5ZBVHWMxjdKWdhK0xGqFaPs3utXijEUbHXSCfk2UDzc1gOflcbNdQA+jTJukkhBmNSEPEH5imafrVYxCOMNDX/UwwbNgwpQFiZlUj3A6Ws3OE3QrhiooK5ObmxnRVE22HCSGEJA55EJ1V+TQzp9jJETYTwsbzjVifnSMs8DpYzk00wriM3WC5aDTaz60GYoVwfn6+aWwwHke4vb0dXV1d+rk4NzfXtJuq0REW71O8H1aNIL4i/4iF6xrUVZcshFViEUDqOcJZWVmYNGmS/nzQoEGYMGGCq3UQQgiJH3kAm5O4lc8X4rHVADgzIWw1lkQWbU5CWDUa4cYRdtN1Tha9xq53QH8hbIaZI5yZmWkqaAVyBLG1tdXxHC7vk9iGfP6mI0x8Rf7Bi/C7+MML0hFWFcJ2VSNU989PRxgALr/88pjHfpQ+I4QQ4g4zR9h4XhDnOHm6LI7NnGIzIWxFPNEIP8qnuYlGyOcqcR61EsJmFSPk9ys7wk5Vk7KysvRzbUtLi+Nd51Ao1K9yhFOcIlWgEE5DRBkZoO8PTRacfl91jR8/Xn+s6qL66Qj7JYSnTZuG6upqHH744fjJT37ienlCCCHxI4tU4TpaCWFZ2IlzRygUinErzRxhp3OEvD2vjnBQ5dOsohGAdyFs5girDMqWm2qpnMONFa3EMtnZ2SltPnGwXJqSl5eHtrY27N+/P0Ys+i2EL7jgAv3x2WefrbSMH0JYHJza29v1afEI4aKiIixcuND1coQQQvxDFqziLqNRwJk1xJCXKyoqwq5du2Kmu3GE5TKgYhtGoSbEZiLLp8ktls2EsChLJu+r/BmJxiVGvDjCQN+AuU2bNqG5udnUpTfbTltbWz8hnMqxCIBCOG2RhbBci9F4eydevvSlL+HSSy/F559/jqqqKqVl/BDC4kAmd5OLRwgTQghJPrK5IYSbSjRCPnfI4k9UN5DFr5kQ/vrXv64/FnXbAfeD5VTLpwnRKp7LeWWr5hx2VSNkrBxhq6pO8TrCLS0tKC4uBmB/Djee+9OhhjBAIZy2iD+s/fv3W+as/CAUCuHFF190tYwfVSPMhPCePXsAUAgTQki60tbWpj8WQtDoCJtFI4xjYwRCTB9xxBH6NFkIL1q0CPPnz8f999/fbxl5G1YZYaPDa5bdld+LlVMcCoUQCoWgaZqraERGRoa+nMCqjnAQjjDQJ4TFOrxEIyiESSCYCeFU+WPzo7McHWFCCDn4kIWwwHheECLNLCMMxLrKQpDKTZLk9VVVVfW7mykLZVVH2C7qIM9vJZDF497eXleD5YA+4SuXPZNfk8+HVkJYbsLlNSNslec2246xmlUqD5QDkjBYLhKJoLq6GqWlpZg4cSLq6uoSup5IJIKamho0NjZ62m6qIAvhoGoIe8WuagSjEYQQMnAxE8Ky4M3JydFjBGZVI4BYISwYNWqU/lg0gbBCLsEmtuE0WM443al8mnG6/NjoLtuVTwOs6ywDffFFgVU0Qszf2dnpyRFubm5WuqtrJYRTxaSzIqFCOBKJYOLEiQiHw2hubkZ9fT3mzJmDmpqawNcTiUQwY8YMVFRUoKGhQb/SSVfMhHCq/LExI0wIIcSMH/3oRwCAa6+9Vp9mZeJYOcLyAHFBRkYGLrjgAgwePBgXXnih7T7I538hup3Kp8nbkV93ml+ezylvbNZi2WzfZBGr0iXVqyMsRyNUDLd0jUYkVAgLoTpnzhwAQDgcxuzZs9HQ0ODKoXW7nrq6OtdiO9WRXVe3sYOgoRAmhBBixk033YTVq1fjqaee0qfJglc4osbpsiMsBJqRf/zjH9i6datjzfvzzz8fp512GqZPn65Pc6oaYZzPKVMssHOEVQbLma1TFrFyzMPK5fXDERbnYruKHHSEHWhsbERjYyOuuuqqmOlTp04FANTW1ga2npkzZ2L58uWYM2cOKisrvex+ypFujrAfg+WEEJZrSBJCCEkfQqEQjj766BghJj8WjihgXTXi2WefxdFHH40FCxbErDsjI8OxdBrQJwzfeOMN3UwD1KMRZmXN5OluohEqneXMtmXcp7lz5+K8887DbbfdZvJuzQfLuc0Ii3OxnRGVro5wwgbLiT+46urqfq+Fw2EsXrwYLS0tjpGFeNdTVlbmar87OztjBJ2oXJBsZCEsfjCp5gjLVSM4WI4QQogTcnUEK0f41FNPxerVq33drmo0wmt2WH7spmqEcd/kDLVg2rRpmDZtmuV7Myuf5tYRFudfN45wupRPS5gjvGzZMgB9YtWIEK2RSCRh61Fl1qxZKC4u1v+NHj3at3XHQ7o5wm6vDIXYpRAmhJCBg2xiycbJYYcdFuh2g3KEZdHqxhE2Vo0QqAhYI14dYTkjrHL+ZTTCASFOzZxa4dKqCFi/1qPKXXfdhdbWVv3fpk2bfFt3PMiua6plhOX8sri6jzcjLIf8KYQJIeTgRK7+IJ/TTjjhhEC363bwm2pG2GywXDzRCBUBa8SrI2wWjbBzhBmNcKClpcVxHhUB69d6VMnNzU3JGniyI2yclmzkz6u7uxs5OTmehXB7ezt6e3v1q1H5NUIIIQcXhx56qOn0CRMmBLpdt1UjVB1heT6rrnNGgWxcTyo4wioNrQ76OsJ1dXV6ZxTVf4sXL3a1M36VNEv30mgqyEI41Rxhsw5Abmsdyz+2vXv3orW1FUBfZszvNtKEEEJSg9NOO01/XFlZia985Su46667Aj+/ua0jHE/5NONzMY9V0wy/HOGenh5dL7jJCEejUWzduhXAwRmNUFYUVVVVMSMsVZg0aZL+OBwOIxKJ2Dq6KgPZ/FpPuiMLYXF1mSp/bEYhXFhY6PoHkZubi/z8fHR0dKClpUUvkG5VNocQQkj6smDBArz00ksxlZ8yMzPxwgsvJGT7quXTrBxh8VwlGmF8LVGOMHCgoYnKevLz85Gbm4vOzk49FjqgoxGVlZVxlR6rrKxEJBJBJBLptx4RZTC2QQxyPemOmRBOFUc4IyMD2dnZ6O7u1p1gLz+I4uJiCmFCCBkAnHfeeTjvvPOStn23mV8rIezFETY22jDOK68jHkcYODDuRnU9JSUl2L59O3bv3g3g4HSEEzZY7uqrrwZgnt8VolYl0hDvesSXme6IP6yOjo6UqxoBHGhh2d7eDk3TPJVREd9ja2srhTAhhJDAcFsOzW0DDrt1mTnCflaNkJcRQlh1PcZzroojTCFswdSpUxEOh/vFKxoaGgD0VWeQaWlpQV1dXb+csdv1WKEy6C6VEfUV29vbUy4jDJjvH+BNCNMRJoQQEiTxVo0QwtJNNELVEY43Iyzu0gLuHWHjOVfFERbRCNYRNqG+vh6RSARz584F0Ofg1tbWYvr06XpnOMHcuXNRW1tr2hrZzXpkIpGI3oLZ7UC+VEM4rm1tbSntCLe1taG9vV2fLhdId6K4uBhAnxAWTj6FMCGEEL+Jt45wPI6w2WA5PzPCwAHhK6o/eBXCdudgRiMUqKys1FsdV1RUoKamBrW1taaD8KqqqlBSUtKvlbLb9QB9AriiogIVFRX6tNraWpSWlqatIBa3J2QhnKqOsAjn5+bmuqr4QEeYEEJIInDbWS6IjHBQ0QjggEgVFZhU9cLw4cP1x3l5ebbR04N+sJxfCBGrMp8QP/GsB+irNNHU1KS8j+mA7LgKUqnRhBDCbW1t+q0Ysc+qmAnhgVARhBBCSGKxyvxaOcJBZ4T9jEbIy7kVwocffrj++JBDDunX3lkmXesIsyBrmiJE5b59+/RONKkkhOXBckKsxyOEGY0ghBASFG4zwsbpVhlhs4YaxudmGWFZcPrpCItohKpLKwvhkSNHKm0j3aIRFMJpiuwI9/b2AkgtIWwWjXArhIX7u3PnTmzfvh1A7G0aQgghxA/cdpYLhULIysrSXVwvjrBxG0IIZ2RkWArhZDvCdlAIk4QiZ4TFjyeVhLAs1L0K4fLycgDA9u3bsW3bNgDOV6SEEEKIW9x2kAMQtxA2OsLC1DLOJ0cL4h0sJ4SwqjgdO3as/njUqFG288r9DeT/KYRJIMhCUwTTU0kI++EIjxgxAkCsEBbimBBCCPELt3WEjcuolE9TzRsbtym7t14dYSGmRRUnVUd47Nix+NnPfoaPP/4Y3/rWt2znFeuMp5FWMqAQTlPkjLC4ErUrdJ1o5MFy8TrCTU1N+joohAkhhPiN285yxmX8cITNlgFiRWu8GWGBG3F6zz33KM1ndIRZR5gEihC9mqaldDSivb1drxrhVqgLR1iI4IKCgpQS+4QQQg4OvDjCQQnhIBxhoxEVRLnVdHWEKYTTFLPGFKkkEv1whIcOHRozYKC8vNy2dAshhBDiBS/l0Mzc4Xg6y1ntix+OcDKFcKqXT6MQTlMyMzNjrrLy8/NdNasIGj8Gy2VlZeHQQw/Vnx955JH+7SAhhBDy/1CNRjh1fFMVzvK8To6wbHz55QgH4dKm62A5CuE0RnaAUykWARzYt3379nkWwgDw5S9/WX98zDHH+LJvhBBCiIzbOsKA+2iEH4Pl0sURjkajzAiT4Bk6dKj+ONWEcHFxMYC+Ui3xCOHKykr98bHHHuvPzhFCCCESquXTnDLCdg013LZxFgSREQ7SEe7o6NDjEYB5lDOVoBBOY+TmEmJgWapQVFQEoE8Ii7qFQhy74Stf+QoyMzMxduxYVFdX+7qPhBBCCOBNpJqVTwsiGiELYa+iMpGOcFdXF/bu3atPT3UhnDqhUuIaWQinWlkx2RFuaWkBcKBlshsmTpyIPXv2IDc3t99BhBBCCPED1cFy8VSN8EMIDxkyxOZdWGMcTB+EIyzv5+7du/VpxveTalAIpzHDhg3TH6daxzW/hDCQ+leThBBC0hurOIKfdYRVhbBdNEKORLohEY6wLK537txput1UJLVlOrElXRzh5uZmAN6FMCGEEBIkxkFoQtS5dYRVhbO8LjeD5bw6wokQwllZWfp7pBAmCUEWwqnqCPf29mLr1q0AKIQJIYSkJkZhqOIIJyMj7IcjnJeXF1i5VbGvFMIkIcjlxMaOHZvEPenP4MGD9R+zCM2XlpYmc5cIIYQQU4zVGEQTiGRkhI3zyZEDP4SwHKv0G7Gvu3bt6rfdVIUZ4TRmypQp+M9//oOmpiaceeaZyd6dGEKhEIqKivR8MEBHmBBCSGpi7FrqtY6wl2iEkyPc1dWlP/YjGhGkEBaO8BdffNFvu6kKhXCaM3nyZEyePDnZu2FKcXGxLoQzMjJSqgU0IYQQ4oRdftcsGqEapZDX5SSE5fKoXtsVy+VLEyGE0ykaQSFMAqO8vBwbNmwA0FdXONVLqBBCCCEybh3hUCiEzMxM9Pb29pvfS/c6ADjyyCPxpz/9CYcccojXt4EjjjjCdL/9RgjfHTt2xDxPZSiESWCMGjUK77zzjv6YEEIISSfcdpYD+txhIYTl6V5bLAPA9ddf73bXYxBNroAD+d0gENsRg+TTQQjToiOBIYvfVBvMRwghhDjh1hEGYkux2WWERS7ZKRrhN0Guv7CwEACFMCEAYoXwmDFjkrcjhBBCiAfsRKr8WBa/VkJYfmzlLBvn85PHHnsMI0aMwG9+85tA1g8ccITTabAchTAJjNGjR+uP6QgTQghJN+wGv8moOMJWQlglGuEHt912G7Zu3YoTTzwxkPUDsREMID3KplIIk8A466yzEA6HAQAnn3xykveGEEIIcYedI6xpmv44HiGckZFh6TT7jbFMnN+koxDmYDkSGCNGjMBHH32EHTt2cLAcIYSQtMPOEZbLmqlEI6wyxeK5qBccVDQiERiFcFlZWZL2RB06wiRQcnJyKIIJIYSkJUYHVXZrjzvuOP2xW0dYnse4fDqXGqUQJoQQQgg5SAiFQjFiWBazctbWauCclRA2c4TN1pVuMBpBCCGEEHIQkZmZiZ6eHv2x4JRTTsH3vvc9lJSUxIjleIUwoxGJhUKYEEIIIQOem2++GU899RSqqqpiplu5vaFQCL/+9a/7rUeeR6XWMGA9kC7dkFs5AxTChBBCCCFpwcMPP4zzzjsPF110Ucx0WaTm5OQ4rideRzjIFshBIypFCQYPHpykPVEnfT9tQgghhBCfGDRoEK677rp+060cYSviHSynso1URe4fkJOTE3i5Nj9IX/+dEEIIISRgZCEcjyPsVD7NbPl0Q/6shg8fnsQ9UYdCmBBCCCHEAjsn14x4HeF0jkYA0KMlP/nJT5K8J2qk96dNCCGEEBIgQUQjDlZHGAD++Mc/YsWKFTjnnHOSvStKUAgTQgghhFiQ6MFy6S6Ey8rKcO655yZ7N5RhNIIQQgghxAK/HGE7sXswRSPSDQphQgghhBAL4hksl5eXpz8eKI5wukEhTAghhBBiQTyD5eQGEwNlsFy6QSFMCCGEEGKBW0c4Go3qj62EMB3h1IFCmBBCCCHEArcZ4Y6ODv2x3FlNtcUyhXBiSbgQjkQiqK6uRmlpKSZOnIi6urqErKehoQHV1dWoqKhAaWkpqqur0dDQ4GnbhBBCCBkYyN3RVERqe3u7/nigtVhORxIqhCORCCZOnIhwOIzm5mbU19djzpw5qKmpCXQ9M2bMwKxZs1BfX4+mpiY0NzcDAGpqajwLcUIIIYQc/MjiVyUaIQthGUYjUpOECmEhVOfMmQMACIfDmD17NhoaGtDY2BjIehYvXoy5c+eiqqoKJSUl+nSxbG1tref3QwghhJCDG1n8uo1GyHCwXGqSMCHc2NiIxsZGXHXVVTHTp06dCkBdkHpZT0lJCSoqKmKmhcNh/XEkElHaNiGEEEIGFrJodRuNkLETu3SEk0fCLjuEA1tdXd3vtXA4jMWLF6OlpSXGtfVjPVVVVXoUQqalpSVmOUIIIYQQI7IwlQfOWRGvI0whnFgS5ggvW7YMgLnoFKJVxZn1az3z5s0DAEyfPt12vs7OTuzZsyfmHyGEEEIGBm6FaW5urul0DpZLTRImhIU4NXN8y8rKYuZJxHpqa2tRWVmpO8xWzJo1C8XFxfq/0aNHO66bEEIIIQcHboXw008/jaOOOqpfZSrVFst0hBNLwoSwHEWwQkXA+rGempoalJWVYcmSJY7ruuuuu9Da2qr/27Rpk+MyhBBCCDk4UKkUIXP88cfjo48+wpVXXhkznY5waqIshOvq6hAKhVz9W7x4saudccoH+7GeGTNmoLGxEcuXL1faXm5uLoqKimL+EUIIIWRg4JdDy/JpqYnyZUdVVZVjjMDIpEmT9MfhcBiRSMTW0RXRBjviWc+MGTOwePFiZRFMCCGEkIFNEEKY0YjUQVkIV1ZWorKy0vOGKisrEYlEEIlE+q1HRBmqqqoCW09dXR3mzZuHdevWUQQTQgghRAm/hGlhYaH+2OgI27nFJFgSlhG++uqrAZjnd4WoVRGoXtbT0NCA2tpaLFmypN9rNTU1rpp5EEIIIWTg4DYjbMXw4cP1x4xGpA4JE8JTp05FOBzuF68QoyrvuuuumOktLS2oq6vrlzP2sp5p06Zh+vTp/RzklpYWNDQ0sI4wIYQQQkzxS5jKQnjfvn0xr3GwXPJI6KddX1+PiRMnYu7cuZg+fToikQhqa2sxffp0vTOcYO7cuaitrUVJSUm/hhhu19PS0oJly5ZhxowZMa8tXrwYJSUljEoQQgghxJSxY8f6sp68vDz98e7du2NeoyOcPBIqhCsrK7F8+XJMmzYNs2fPRklJiS5gjVRVVaGkpKRfK2W363nvvfcAHGjNbLYdQgghhBAzfvjDH2LdunW4/PLLfVvnrl27Yp5TCCePhPvvQsSqzGfWGtnteurr613tHyGEEEKIoKCgAE888YSv67RzhBmNSCwJywgTQgghhAxkJkyYAAC44oorYqbTEU4evOwghBBCCEkAS5cuxVtvvYXzzz8/ZnpBQYH+mEI4sVAIE0IIIYQkgLKyMlx88cX9pn/5y1/WHzMakVgYjSCEEEIISSJyJ96Ojo4k7snAg0KYEEIIISSJFBUVoaamBocddhhOOOGEZO/OgIL+OyGEEEJIkpk3bx6i0SgyMuhRJhJ+2oQQQgghKQBFcOLhJ04IIYQQQgYkFMKEEEIIIWRAQiFMCCGEEEIGJBTChBBCCCFkQEIhTAghhBBCBiQUwoQQQgghZEBCIUwIIYQQQgYkFMKEEEIIIWRAQiFMCCGEEEIGJBTChBBCCCFkQJKV7B1INzRNAwDs2bMnyXtCCCGEEELMEDpN6DYrKIRdsnfvXgDA6NGjk7wnhBBCCCHEjr1796K4uNjy9ZDmJJVJDNFoFFu2bEFhYSFCoVDg29uzZw9Gjx6NTZs2oaioKPDtEf/hd5je8PtLf/gdpj/8DtOfRH+HmqZh7969OOSQQ5CRYZ0EpiPskoyMDIwaNSrh2y0qKuKPP83hd5je8PtLf/gdpj/8DtOfRH6Hdk6wgIPlCCGEEELIgIRCmBBCCCGEDEgohFOc3Nxc/OxnP0Nubm6yd4V4hN9hesPvL/3hd5j+8DtMf1L1O+RgOUIIIYQQMiChI0wIIYQQQgYkFMKEEEIIIWRAQiFMCCGEEEIGJBTChBBCCCFkQEIhTAghhBBCBiQUwoQo0NDQgOrqalRUVKC0tBTV1dVoaGiwnD8SiaC6uhqlpaWYOHEi6urqbNfvdn5CCCGExA+FcIpCYZQ6zJgxA7NmzUJ9fT2amprQ3NwMAKipqTH9XiKRCCZOnIhwOIzm5mbU19djzpw5qKmpMV2/2/mJf1RXV9t+zrygSV1aWlpQW1uLiRMnWs7D7y/1aGhowMSJE2NMhcWLF1vOz+8wuUQiEdTU1KCxsdFxvqC+p8C/U42kHE1NTVpJSYk2ffp0/Xk4HNamTp2a5D0beCxatEgDoM2cOTNmelNTkwZAM/sJVVZWaiUlJTHT6uvrNQDa8uXL456f+MPs2bM1APrvzIjb3yF/t4lj9uzZWklJiTZz5kytubnZdB5+f6nHzJkzNQDaokWL9GnTp0/vN03A7zB5NDU16d8NAK2pqcl23qC+p0R8pxTCKQiFUeqwaNEiraSkRJszZ06/18wOEMuXL7cUVwC0qqqqmGlu5yf+sHz5cq2kpMRWCPOCJjURJ2enz5TfX+phd6wzEzb8DpPD7NmztcrKSmUhHOT3lIjvlEI4xaAwSg+am5tNHWFx4Kivr++3TDgc1gDEOFhu5yfx09zcrIXDYW3OnDmWvzVe0KQmU6dOVToB8vtLPcRnbBS84lhqnM7vMDWorKy0FcJBfk+J+k6ZEU4x5syZA6Avu2gkHA5j8eLFaGlpSfBeESPz5s0DAEyfPj1m+rJlywD0fVdGSkpKAPTlnbzOT+Jn2rRpqK2tNf3MBW5/h/zdBs/cuXPR0NCAmTNnorKy0nZefn+ph/jOGhoaYvKms2bNAgDcddddMfPzO0wNysrKbF8P8ntK1HdKIZxiUBilB7W1taisrNR/qALx3YjvSkYcUOTvz+38JD4aGhrQ0tLS7wLGCC9oUo/a2loAwOzZsx3n5feXmojj5cSJEzF37lzU1taioaEBy5cv73dxw+8wPQjye0rUd0ohnGJQGKU+NTU1KCsrw5IlS/q9pnJ1Kn9/bucn3hFVBowXL2bwgia1mDt3LlpaWhAOh2PKGApBZYTfX2oyffp0LFq0CCUlJZgxYwbq6upQX19v6vDzO0wPgvyeEvWdUginGBRGqc2MGTPQ2NiI5cuXm/44VXC7nNftkFimTZuGGTNm2EYiBLygSS3q6+sBAFOnTsWiRYvQ1NSE5cuXA+j7TRpL4PH7S23Kysr0uzITJ07U3X4ZfofpQZDfU6K+UwrhNITCKDnMmDEDixcvthXBQmTZ/YDlzJXb+Yk3GhoaEIlEMHPmTN/WyQuaxCFOdldffbU+LRwO63dljLlTFfj9JZ7a2lrU1NRg0aJFmDNnDpYvX45wOIy6ujpPddP5HaYHQX5PfnynFMIpBoVRalJXV4d58+Y5OsHiFp/ZVaqYVlVV5Xl+4p6WlhZMmzYNd911FyKRSMw/ANi9ezcikUjMb44XNKnF7t27AfQ/6ZWUlOi/D5EnBPj9pSINDQ16FEJ83pWVlfoxVeT3BfwO04Mgv6dEfacUwikGhVHq0dDQgNraWixZsqTfidjYcUc4VlbfX2VlZcw63M5PvFNTU4OKigr934wZMwD0fb8VFRUxHcp4QZNa2J0Qxe9DPiHy+0s95syZE3PhIigpKdFjEnKHOX6H6UGQ31OivlMK4RSDwii1EG7i9OnT+w3oaGlpQUNDQ0zmdOrUqQiHw/0GZDU0NADoXyLI7fzEPSUlJWhubobWVzdd/7do0SIAfQN4NE1DU1OTvgwvaFKLSZMmATD/fIU4ln+H/P5Sj927d1u6d0OGDAGAuL4TfofBIO7GWBHk95Sw79SXasTEV8LhsBYOh2OmiU4qZo0XSHCINryiy478LxwO9+t4o2kHioCLbnSiJaRVBzO38xN/EO2zrT5nt79D/m6DQ/xGzArol5SUmE7n95daiOZBZg2Cpk6danos5XeYfERDDbsmNkF+T4n4TimEUxAKo9RBdLKy+mfV2Wb58uVaZWWlFg6HtcrKStMWzfHMT+JHdJaz+w55QZM6CCE1e/bsmGklJSWmXa/4/aUWzc3NWklJiVZZWRnzfYnf4aJFi/otw+8wuTQ1NennOvl3ZyTI7ykR3ymFcIpCYURIMCxfvlwrKSmJuaApKSkxdRd4QZNa1NfXa5WVlVpJSYkWDoe1qVOn2rYg5/eXWjQ1Nel308RnPH36dMv2vZrG7zAZCLFpNH5KSkpML1g0LdjvKejvNKRpmhZ/wIIQQgghhJD0goPlCCGEEELIgIRCmBBCCCGEDEgohAkhhBBCyICEQpgQQgghhAxIKIQJIYQQQsiAhEKYEEIIIYQMSCiECSGEEELIgIRCmBBCCCGEDEgohAkhhBBCyICEQpgQQgghhAxIKIQJIYQQQsiAhEKYEEIIIYQMSCiECSGEEELIgIRCmBBCCCGEDEj+f2ErG3pQX0/EAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig = plt.figure(figsize=(8, 4))\n", "\n", "def dispersion(k):\n", " return np.sqrt(k)\n", "\n", "def transform(k):\n", " a = 0\n", " n = 20\n", " rs = np.random.RandomState(1)\n", " for i in range(n):\n", " a0 = rs.random_sample()\n", " k0 = 20 * rs.random_sample() + 1j * rs.random_sample()\n", " a += a0 * np.exp(-(k - k0)**2)\n", " return a\n", "\n", "def wavepacket(t, N = 1000, x1 = -10, x2 = 10):\n", " \n", " \n", " kmin = -100\n", " kmax = 100\n", " dk = (kmax - kmin) / N\n", " k = np.arange(kmin, kmax, dk)\n", " a = transform(k)\n", " w = dispersion(np.abs(k))\n", "\n", " a = a * np.exp(-1j*w*t)\n", " \n", " eta = np.fft.ifft(np.fft.ifftshift(a))\n", " eta = np.fft.fftshift(eta)\n", " \n", " return np.real(eta)\n", "\n", " \n", "\n", "\n", "eta = wavepacket(100)\n", "\n", "ax = fig.add_subplot(1, 1, 1)\n", "\n", "ax.plot(eta, color=\"black\")\n", "#ax.set_xlim(1200, 3800)\n", "#ax.set_ylim(-0.02, 0.02)\n", "#ax.axis('off')\n", "\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.6" } }, "nbformat": 4, "nbformat_minor": 4 }