{ "cells": [ { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "![alt text](https://github.com/callysto/callysto-sample-notebooks/blob/master/notebooks/images/Callysto_Notebook-Banner_Top_06.06.18.jpg?raw=true)\n", "\n", "
Sound Demo 2
\n", "

\n", "

Michael Lamoureux
" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "slideshow": { "slide_type": "skip" } }, "outputs": [], "source": [ "%matplotlib inline\n", "from matplotlib.pyplot import *\n", "from numpy import * ## for numerical functions\n", "from IPython.display import Audio ## to output audio" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## What is sound?\n", "\n", "Rapid changes, or **vibrations** in air pressure are heard in our ears as sound. \n", "\n", "A computer creates sound by sending a list of numbers to a speaker, which then **vibrates**." ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Random numbers\n", "\n", "A list of random numbers, will create a sound. \n", "\n", "What do you think is will sound like?\n", "\n" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "slideshow": { "slide_type": "skip" } }, "outputs": [], "source": [ "Fs = 44100\n", "random_sd = random.randn(Fs)\n" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "slideshow": { "slide_type": "slide" } }, "outputs": [ { "data": { "text/html": [ "\n", " \n", " " ], "text/plain": [ "" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "## Random sound\n", "\n", "Audio(data=random_sd, rate=Fs)" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "We can plot a few of those numbers, using a plotting command. " ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "slideshow": { "slide_type": "fragment" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAD8CAYAAABjAo9vAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzsvWewZdlVJvgde++796U35UspqeQNEqoSCNNIGloSTAximAGhjoaehghNdE8PmqZpAdMD0dMGZoKJoQMXIESjpgGBACHRLQkGI4sMqpKqSqZKUvnMMmle5svnrjlmz4991t5r77OPue/dZ27mWREV+eq64/Ze+9vf+tZanhACnXXWWWedXTvm7/cJdNZZZ511Nl/rHHtnnXXW2TVmnWPvrLPOOrvGrHPsnXXWWWfXmHWOvbPOOuvsGrPOsXfWWWedXWPWOfbOOuuss2vMOsfeWWeddXaNWefYO+uss86uMQv346AnT54UZ86c2Y9Dd9ZZZ50trN1zzz2XhBCnmj63L479zJkzuPvuu/fj0J111llnC2ue5z3e5nMdFdNZZ511do1Z59g766yzzq4x6xx7Z5111tk1Zp1j76yzzjq7xqxz7J111lln15h1jr2zzjrr7BqzzrF31llnnV1jdt049mma4713n0Wed60AO+uss2vbrhvH/ulHVvCOP74fX35qbb9PpbPOOutsV23Hjt3zvNs8z/uI53lf8Tzvy57nvX0eJzZvm6a5/DfL9vlMOuuss8521+ZRUiAF8C+EEJ/3PO8QgHs8z/tLIcRX5vDbc7NcSAomy/f5RDrrrLPOdtl2jNiFEE8LIT5f/L0O4AEAt+z0d+dtonDsad559s466+zatrly7J7nnQHwSgCfdbz3Ns/z7vY87+6LFy/O87CtjJB659c72w9Luq1iZ3toc3PsnuctA/gTAP+bEKIUoRRCvFMIcacQ4s5TpxqrTs7d8g6xd7ZP9sTKFl70M3+Orz6zvt+n0tl1YnNx7J7nRZBO/feEEO+bx2/O28ix07+ddbZX9uTqCGku8NTqaL9PpbMZbDTN8IZf/Bg+99jl/T6VmW0eqhgPwG8BeEAI8f/u/JR2x7rgaWf7ZdNi0GXXYQ7F3zx4Hv/09+7Z79PYll1YH+Nr5zfw4ALutOaB2L8VwA8BeL3nefcW/333HH53rkYMTNZRMZ3tsZHUNrsOd4t/9+gVfOiLzyjxwiLZ1lRKo7MFRIM7ljsKIT4JwJvDueyqdYi9s/0ycuzXY9YzBY2nWY5eGOzz2cxmo0Q69nQBn9t1k3mqHPsCIofOFtvIuV2PYy8trj3JFu/aR4TYO8e+d/blp67iXZ94pPXn6dlci1TMJM3wC3/xILam6X6fSmcOU1TMAjqIndq0cOh0DxbJiIrpEPse2gfufQo//+EHW3/+WqZi7j93Fb/6kYfxuceu7PepdOawSTHorkdFFiH2xXTsEigt4oK8sI59kmTIctE6KEP85qIi9vvOruLclS3ne2mBiq5HDncRTCP2fT6RfTBCu4uYoDXuOPa9t3FCKKjd5zUVs0sntMv29j/4An7lbx5yvqd3I4s3AK8Hu56DpyT1nCwkYieOffHOfXEdezpbYGPRg6ejJFNRetvoHizqtV3rdj3LHXXwdPGcY8ex74NNktkCUsr5LeAAA+T5Vw0wurbrEREugiXXcYJSusDBU6JisgVU9CysYyfE3rb2C4GlBXxGACRqSCsWpQ6xH2ybXsfB0ynTse+3pVk+U3mADrHvgxFib0t/aR56/wfYdizLRCXiS/OOYz/IRmg1XVRUsQOja04OAGL/yFcv4vt//dN47NJmq89vdTr2vbdZEXu24HLHtIaK6QqcHWyjwOH1+HyIhpocgIm3Pk4AAFe2pq0+36li9sHGM3LsNKcWdXJluahEfIqK2f+505nDrucEpSQ/OBw7zR9C4k2mdez7f+6z2sI69gmpYmbUsS/qdjjN88rdSRc8Pdg27UoKHAhVDD2HzUm7DO2OY98HI469raPOFljuKIRALqoRXxc8PdiWXMc69uQAZZ7SItMWsStVzAI+t8V17DPr2It/F/Ah0TU2yR0XcQBeD6brse/P8cdJhrf8xqdx39nVPT+2Cp4eAMRO82ezZU2lDrHvgymOvSVK1c2sF+8hpQ00UtYFTw+07XeC0qWNCT776GV8+pGVPT/29AAhdjqXrUk7xK6qOy4gfbvAjn02xK546F2cXP/kd+/Bz3/ogbn/bofYF9v2u6QAAYJL65Md/9blzSkub7ZTlfBjH4SSAnQu1wNi33Gjjf2wNMsbUaxt9Gx2M3j6paeu7goy0Tr1hgSlBRyA14NN9jl4SuPn0sbOHfs7/vh+eB7wmz98Z8tjH5x67MmMHPsoWdxaMQvp2Pnq317uuPuIfX2cKnnXPC1romL2YDfS2fZt3xF74ZguzsGxX9mawp+hXxpd+0GgYmhxaauKGS0wYl9IKmbMimG1RUG7jWqFENKx7wpiLxRADQlKByA+dc3bxiSdeQztd60YTcW0p1AqfyuvzoCu+jxwQIKnMyD2NMsXugn5Qjp2E7G3GzCKitmlh7Q1zYpCXfMfwE2LUtoh9j0xIQRe9/98FL/z6cdK713amFQiwf0Ons6TikmzvHWpbMDsebrflsygY+eVVDvEvkfGEXt7jn13k3jWx3Kw7AaX2CQZq3P8X31mvbJ42HbtUw9fwv/35Wfm+puLYGkucHF9gsdXyg1Pfvi3/g6/8BdfdX5v36mY4vlf3prueCxkuWgNIIQQaj4cCComb595OmKf6RD7HpmB2Ntmnu5ygtJaUYdiNxB7UyONKsd+eXOK7/6lT+CDX3x6rufzrk88iv/wV1+f628ugtG4WxslpfdWNieViHi/M0/JuQqBmRQtLpuFiuGfOxCIvXh+bVQx3Pl3iH2PzODYZ220sWuIvXDsu4HYW8odbSS1WfDBO53MtiVZfiA40722STHuaBHnlmaiEpXud60YftydBlDTLG99HXz3ehAQO82fURvEbviY/T/3WW1BHbu+0W1XU90ab5cQ+4iomN3j2Ku20VWLFr3O79c8LM2qK01ey6YRexnxJSzYZtt+O/aEOaZLGztH7G03Hvy4BwEIKI59BsTeC/2FrC+1kI6dygkA7XnLvIKumJcRittNjr3KmSqduzXj6FrHFS31tn0++fWJ2MlBuxB7lrsRuxBi30sKcMd0cYdJSlkuWlNK6QFD7MkMmaeE6g/1o45j3yvbHmLfbSpGooB5ByqBZlUMLVr2IqcQezpfx57UNP24lq2OY08qHDtf6PdLtZQZiH1njj3JRGswxRf/g+DYZ8k8pZK9h/vhQo71hXTsHLHPWgRst4Onu5GgxHXswnH+OjPVfJ3+fzJvKibPnTuTr51fxwfvn2+g9iAZjbu1cdkxpBVUDH9t36gY9qx2WlYgy/PWc8hw7DWA596zq3j5v/6LuceCbKNzGCfNcQLi2A/1w4WkHRfTsSezT5bdpmL2ArHbf9uv2Yhw16iYzK3X/93PPI6fet/9cz3WQTJC7BuT1HjOeS7LKrtQ6XQbCq55G40Dz5tD8HQGuWPb4OnjK5tYG6c7pomajFNDWw2onaiY5Q6x752Nt4XYd5ljH+2+Ksb+m6wpeDqas2NPstx5nUkmsDFJF7I0chvjgGKDJblQkLDJse/XfSHkfHK5t2MqJs1E6z7DfPGri8nQWNptB8rBSJOWnd4/1It2RcK827aQjn2yA1XMbvGcKkFphkEwTXNnIM62vMGxzzt4KoSodQBZLpwTNctzCAFstKyet2jGKUCujKmrYDjdRl2jeRuNj5uO9HdcViCbQcduIPYax75XtY6m7Hyask8JDB1e6hD7nhl3VG1REDm93ZIurW9DFfNrH30I/8Ovfarxc9yZu2pDVwVPsxq540/+8f341MOXnMf75EOX8M0/99e4sDZ2vp9UyB1p7q47OOhrwbiT5gsyjakmjn2/gqf0rG443MfK5k6pmNk59jjwa6mYpjyNeRnfQTQh9tE0Q+B7WIqCA1GZclZbTMeezp4VRkHH3QueUuNbd4DTZefXJrjQglfMDMReniBVtWLyCsQ+TjL84d1n8cmvux37uSsjpLnAqkP9Qefguk46/nqLXcgiGkfkXBmTtqRi9q8ImDyHQ/2wVXJOlVEsoS2Yovsy6NU7R1Lt7Pb9SbIch/qyoG0TYt+aZliKAoSB3yH2vTIzeNqyCFhOn99dxA60R+1Z3i6Lr5Fjr1TFkNzRfIMQdRWKokFftbvRtWvc1M+1itgNKoYj9uK6nY6dq2L2yT/Q8xrG4Y5S+2ctNjdN2XFbIPZd59gzgaODCEALxJ5kWIoDhL7Xcex7ZeM0UzWhW9et2CMdO9C+Xkyat6sGyRcvJwVSUWueXp9YiF3TRlWOnepQu9+n79nva8e+eIh9kmb46ffdX0k/yc9wxM6CpzUVDA9C8JTGzCAOZtKh2zZr6WuF2OOgtoPSXjWKSfIcR5akY2/Ssk+SDP3IR+B7HWLfK5skOYax3FLNSsXsFs+5NkrUYtMesbcLRDVx7FWIh/yuTcUoxF7h2EkK1lSb5lpC7F8/v4H3/N3Z2r6gfKfo4thdz/MgUTGDYs5sF7WTMKDtZdB9GfTCelXMHgVPk1Qox96UfZrkApHvF4j9OnXsnuf9R8/zLnie96V5/F6TjdMcg14AYPYEpV2RI2Y5NqcZjg9jAO3rYkjE3nw+/BpdqhtFxVQgdjt4So63CkURmqm6t5SEZWv2M7G4jn2jgX4CbFVMmWMHys9+msnvLEXB3B37M1fHWKlRL13dSmTpXIbYXefY1rIZZYm0gAzjoJaKaerpOy9L8xxHl+QcbULsWZ4jDDwEvg8xQ1zhoNi8EPu7AbxpTr/VaOMkU4h91g5Ku4EKyCkcG8hB03bxyDLRatDw33NNKoV4Sohd/r+tY2+qa6OomEqO3R3sWmTEvqHq6Vc7oEmaw/NkmjnPPuX30V4siWcexMHcA/c/9p4v4N/816843zu/Nsar/t1f4jOPXJZOyvcQh35xTttz7LMia4XY4wbEnrnH77wtyQQOL7Xj2NNMIPB9hIHchi8aap+LYxdCfBzA5Xn8VhubcMTOJtXq1tRAVdzmkaD01OoIf/i5J0qvkyM7NjNi16UC6sxQxdTIHascbZmKkY59WnGvKHjqulekjADK5RMWmWMnBFf37KZpjl7o48ggMhF7TYYlodb+LiD21dG0Mt/g4voEaS7wzNqocFLMsW8TsaeKipmNYx/2ApmxWknttZsHO7UkyzEsAqJNmadpLhD6HoKCX617dm9952fwnz/9WOn1X//Yw3jv3Wd3csrbtoXk2Dli54Phzb/6t/j1jz7i/A6NxZ1Mrg/c+xR+8k++WBoUhIBPFI697QBtqwbgSM8V0GwMnqa5IU1sUsVsTauDp5wKsqmYfIGpGB13qKNicvTCAIf7kaWKqa6JQve4H/lz3y1muahEnrRzSFLZxSgKfMTBDhH7rFRM2o7b30tVTBT6GMSB2pVWfjaXi2HoE2Kvvmf3nVvFV55eL73+/i88iQ/PuclNW9szx+553ts8z7vb87y7L168uKPfkhHrAL5nDoYLaxOcvVJuWwawyPsOJhchC1sLTAoJQuxt68VobrH+801yxypVAUdILg12FRWzUYPYOTq1v0/vLSRiVxx7HRWToRf60rHzzFN2n0qIvfj/pXj+iD3LRaUunXaukyxXfPFOqRhNZ6JVrgbdl2HB7Vc59r3IPKXyyZHvYdgLGxF7lueIgmbELoTAOMmcTEHOSjbvte2ZYxdCvFMIcacQ4s5Tp07t6LcmxZY49H3DUSdZ7iypCuhBs5PgKT0jm7NetxB7W1VMWwSUNVQI1IuW/P+vnV9Hnpt1szkds9aI2KtVMeYic+0FT2s59iRHL/JxeCk0EHtdeVqiu3YjeJoJUVkHaKoQe46koBWiArHXSQ/rjD/vNpdC92XQC41zKv/u7gdP6d5HQYHYW3HsHLFXxJsKatJ1T7NczL2yaltbWCqmHwWGxlQIqTCpcirk43aCCqpUJpSheWJGjr1JDfCBe5/E186vG++7fpsHT89e3sIbfvHj+MRDlwxHws9ZqWIadOxuxM6pGPcOYb1FF/iDZuTYm6iYOCDEPjvHPncqJmumYqZZjjTLEfo+ejvm2PX5t7kWAjjDBjVOxsbvbhmdexj4OLIU4XJDJynJsfsIfN84R9voPtu5IvSdhUbsnue9B8CnAbzA87xznuf96Dx+t8omaa6SB+wsyKqiWvNIUKpK0V8pBsmNR/oA2icoEV9ddU7/x/u/hN//7BONZXt58JSu/+ooMSYfP2eVoNQgd0yLrf47P/6wclh1i8wiI/bNNog9zSTHvhQZqhiTYzfHBo3LXUPsFY6dI/Y0FwYVU/Xcm6xJnZXnAv/uv34FZy9vFZ9vi9jNedC2JMcsRg42CjzcemyAc6tuylafUzvETvPKidjFgiN2IcRbhRA3CSEiIcStQojfmsfvVtk4kRMs8D3lvGhCVjmVeahidA9R27FPsBTJoJo8l3bHaMq4SzOBSTEx1WsuKobVwdG/mRslBsaOdPjKBCWF2HN85tEV/NyHHsTfPHi+uLbqLNhFVsVstOLYCyqmHxk12evkjvT/vShondjT1mTwNHU6QhU8LUosh3NRxdQj9mfWxnjXJx/FR796QR0bYBx7hWPn8+Br59fx4p/9Czx6aXNb51h57pmmYm47voSnVse1z7rEsVfMaY3YXfkl27/XO7UFpWIkYud1HGgQVSF2GpM7Cp5WUDErm1OcWI4RFhxmWx6/iWOn8rhNcke+lSUnk1qp4y4qpqpoFQ3GNNPt3j718Erp+CXEvkA69nf88X34adYUZKNYzOoWZYrtDAup7VZS1vu7gqdx6CPwdgYqXJbV8Lt0HpMsR5rnCOegiuGlLVzXonqKTs17OWgInqp5IATOXdnCKMnwuUfnq56mcwsDD7cdGyDLBZ6pKR+hOPagXhWjEXt555Tm+b61BFw4xy6EUFtin3HsU4bYXQhG0RU7CJ5WUjGbU5xY7qlB0LYmexPHnhQDwygp4JI7MsUPVxjwRYxv2ddrknG4WiDLhZp0yrHXLDJ07I2J+xkcJLv37CoefEZL1DYadjGAljsuFY5q7JCFuhx7L/CNsTovSyvGozxu4VxTUUbsc6BiXJdSduw5At9DL2yH2FM23vizmYcliorxceuxAQDg7OVR5ec1x16viiGk7g6euh3+XtjCOfYkkyiFELsaFAz9ugJKOaMrtmv0bG0lwsrGBCeHsUJErRG74hZd2ziZlZoUcjX9nXrEzpOejOBpWubYXRONqwV4kbKHLmzgwtrYcGJVCVF1+uqDYisbU2Oxo4BxHf88SaTccSmSjmrkQuyOkgJR6CPwvF3RsQPuLEodPM1KHPtcqBjHOKQsW3VfcoGIc/tVmcxs/NI1PfjM2rbOsco0FePhtuNLAIBzFdJoQN7bNhy7kpU6xk0uxLYVSDu1hXPsdCPjsAie5ibHDrjpGEXF7AA1VWVyrmwQFeOVzqXN79XJCqVjZ6/XUDEGYs/N3pQUtRdC1ModeZ1qjtgB4NOPrNRTMY4kqDr72vl1vPlXPrnnnHyWC1zemhqTTnHsNeNjmuboRYGiFlyJXE4qJtidKoF1jl0HTyWdF/qaitm+3JHtHB2LlEbsekcY+b6SWbbh2OkYDzy9Ntddn6JifB83HVmC5wFnr5iI/a++ch5/cs85ACjoK69RFUMUp2vXlOWio2LaGt1g2iZlDsfucirzaGatgqfsYQkhsLI5kVSMT8ikpSqG8eG2aYpJNCN2FRjWv1VC7MUAHCUZslzA99zIjTt2jth9D/jbhy7VBk/zAuUAwMak2Vl/+amruO/cVTy5Wr0l3g27sjWFECY9peWOTVSMj76F2OuaNieZUCBktxC7SxnD5Y6ZhZznwbG7EDuNFU7FRKHPdgruXRxXxdA1XdlK5trcOlGIXZ7PTYf7OHd5C5M0U2P6tz/1KH7zE4/Icyroq50g9iwXpazvvbLFdexFxFo5wFTfPFeS0jyoGIXY2URaG6dIMoETwxgRBVpqqJgP3v80fumvv278nruwlylX0+dQzbHnDPHYJWQJUdCid3wYY5qVBx1Hf1muj3378QEeu7TVKHc8WhRZWmuB2FVLuT1GNSRPHbFdjEpQqqNi0gwxo2IUx84TlBwlBezd5byMxrIri5LOY1qoYgLfU4h9u9UdkwaOXVExU01RycQoz3i/dB3F6WTCrHb6wDZ49vd9/hweXykrajTHLs/l1uMDnLsywo+8+3P4sfd8AYAcFzxZKjA4dvc9a9Kxy2N3jr3R6Gb5nmdy7HkDYi8mwU5KcNLX+LaLyqaeXO5pVUxN8PTPv/yMKgxUl3HHg8JG2d5GKkajH5eOnRa948MYQpSPvWEj9uJ4y/1QKizY8UvB00zgSNGhpg0V46LR9sKocNaIaZBdOz/bJolE7BQ85Vwymb1IUVKT782XiqEYDKDVOea56phBkucKqbrOsa0Z+RS1VIw89jSTx21KjMoMxK4/89UZeXYhBH7ij+7DH919rvQezUmao7ceW8K951bxtw+t4OsXNgAAlzamRuVSA7FXOOcmHTuwP5LHhXPsKoPM9+B7s3PswPZRu1LFsEDkyqZEfyeWY0Q+cezVv5/lOQv06kFkG/1GkknUXNcxijezTgwqRn+G6CNC0ieGPfX73AxVTCZ0WnjR3iw1aKFqxN6GN6dnZ08KIQT++oHzu5aJSI59Wjh0ezGrsklWqGJqqJhS2d5MI/Z5Xo9RLsLFsTPETk5qXmV7ATc4orHCEXsUeIiDelUM32XS33Hg40FHYa06I/mna3FOWPAUAG47NlDnc3F9gjwXuLw5YfRorlgB+m2X0fOWIMgt/3Wh+d22hXPsdLNIY0oDzKBiahA7/42Zjy1oq6kfICH2E8OeChI11Z62y/U6Kzaq68qRZUJJxuoCrXXBU03FFOUPlmP1+9w2WNU7ztMP4wCTNKstApbluiZ9K8ROzsc6h/vPXcWP/qe78ZlHq7sZ7cRWWDr5JM1ULXbXuZAJIVTZXuLYCZlmtcHTTDn2edZj52PYqYopYipJliOh2uK+B68ittLGuOOqAyNbCQWipX4+CutFBZySpL9vPzHAuRljL3U7YC53BCRiB6S67uoowYX1ibEoUNnepnrsfCG376stxd5LW1jHThFrJ2Kv4dj5b8xqOniqJ9KlwkmcZKqYOo6dR8rrEpR44lWaC/QjklK6ZVX0O3xw80WQgqfkcE8uS8RuO6ItWxVT/AYhdiN4WhrIuWpk0NQFnl+3fQ5Ue2e3Ep1WNnVQbjTNDMRe5Xx0BqmvVDFjC7GHvufk2HuhpGLm2ROZOxonFcMWzTSTWZSeJ3n2eSD2VlRMWuwUGlQxBmJnJRjaVkm1j++ucUSIXZ7Lq599HC+5+TB+9NueDUDLK8sce0OtGHbveRIg39HshzJm4Rw73Xjf84xsvkaOnd3bHVMxBscuHfuxYawRe80M5u3w6uSOPPCS5XkrxJ4zjjJn6H0QB5pjtypR2o5oc5rB8+RWOMlzRsUEpWQpV0kBQrNtUAo5RPuztJV3ScjmYRyxjxLt2HuhX0mjKcfOEpRGTO5IjSymaY7/+88fxLv/9lH1vV7oI/B3Fri3jTuakSt4qlQxchwQtxyH/rbljvyYLqVHiYrJJQ0VhfU7WZdctxf6MwOwlNGXttEYI878WSeG+OCPfTte9axjAHRCFF8c2qlicva3udt1fWavbOEcO6HT0C/K9jqpmHrEvu0u7cXXDMe+OcHRQYQo8BsDLYCJSpTMq6ZB9TTNkQmJFOn7rt+U52dy7LkQ8DyJfmxVzAkLsT90YR1/dt9T2JykGEQBosBDVjgFz5P1xCcWYredYC6gdhZJhQLCPG83FTMqtvK7VUDpEnPs4yRTu4tjg7gGscv71wt99MNyghJx2EmW40NffBof+9rF4ntyUQ5aBE8/88hKa01/IxVjlO3N1djshf6cqJjy+zZiV/elQT/PgQn93Y8CY3ytjRP8d7/8SXztfDXvXldUz0bsZKeWZeG+B59eU58TghZDTcVUqmKYL+DjlfubDrG3MHpAge/NpmM3uhDtkIphD3BlY6rQb6Acex1il7VYeIs5F5Lj15UVqgbPq6Bi1A7AlD5muUDgeegzx74xTuF7UN3ap1mO/+vDD+Lv/+LH8WPv+QLufvwKBr1QyfMSNjkp2KiuxTqXtDhP33PHDUbTzPiOQuy2Yy9iGLuVjm1SMbnuWTuscewJIXZZHqAX+kbwlDoUTdMcq1uJ1pETFVOMjSpQcWVzirf+5mfwnr8rt150mYHYa0oKkNyRHHsc+Nuv7sipGMd1UMnjkaWKaaJieGs8+rsf+cYYOnd5hC8+eRUPPF2tlOGgpnzuxLF7xuunDkmAQ4g9ZRJfG7H//mefwFMW7z82ELs7BrEfZQUWzrFz3jhgRcDo3yNLkZNj5/N123JHBxVzaWOi0C9xmHYvUNf5TzP3ICAjB8MnZlihheZNhg2OXQj4vod+5KvFyJa+JanAf/rUY3jV7XJLet/ZVSz3QoSBX+wuilreka/OxT6uvj8U1Hajwv/2lz+BX//Yw6XrrlLm2MXW5mWXNiYqxsCpmGODqJmKKaimpTgwKAdK2R8nGdbG2rGT9j3wCuRn7Rx/7kMP4NyVLTy2sgkhgKdWqwtTcTOpmBrEnlHZXvm8ox0gdn5MV7JVaoxZ4vb1QljVtYjHmmhM9cLAojOai7TR8V3ghxYVG7GTiOChQvKYZIKBR82xX9lK8L//6Rfx/nufNL5vIHbmwPlz7qiYFpYJvZpKpYF8nbb+J4axczsrhFCr9Xa5TpXpZ1AxGrEDMqhbh4hchZvqOXaJku0dCjcePDV07ByxF4Muy8zGxuM0wyjJ8K13nMRzTg4BSD6dELuqMxIERSYdP28reCrk8eLAd9JRZy9v4eGLOnmEts72wN8Ljp3qhYwTrYqpo2LIMRD6XIq0Y5e7GrlYrmzKrFaekcgRO39+T66O8M6PP4IP3PsUnihqmF9Yb+nY2RiuKylAwVOO2LdLDZgJStWBV3vvAAAgAElEQVRgBJAB3SQTiso4PoxxZauiV4KlivE9iaxdVTPrdsP1iN1NxUSBj+PD2JmLwBE7+RRbWsrBh4HYa7KR98IWz7EXzoAK9CietnjgJ5bjSrkjpfxvtz2ei4q5OkpwtEjKAVCJqvX5lxcHF3+n1T5CSa8ipgJyfZbr2OUkkfepH2oqJhNmFuLVYrINewG+5Y4Txd+hurdpwc/SQmAUCWP3kXhJv5CI2RNQKmoELm9qfpsGfyl4WpP0sVPbmqbYmmaqwt+o4Ng9Dzi8VIfYC469iCEsxYE6T6rdHQc+LqxJmkdV/UtkfRmi6bhDpO8/dGEDT6xIx35+rV0aPXccjYidOVgK8G7Hmsv26tfGU5mqT4702CDGlU131yIeI0oViDGDp+p6auaW3bDDPLfCWVtUDACcKnZv6tyLZ8117ETv2rQXBzociHSIfUazOXY7En68ArFTvQxg++3xaLzw7df6OMGhvnbsUeA36tgBczLWFvbKpTwy8D0EDocphM5A5KoCCp76nnRGtBgR+o8LbTFJCwdxiG957kkAUrNOiJ0UFZQ9uGVIAzmCk/8GnofQ90st5uh6V7f05OYBYm5bu4jYSRFDOuZxkmF9kmI5DgtVTIPcMdSIfcyCp0FBNxDi5rVaepyKcVAoD13YwOMzIna+W3KWFFBUkKRiIqaK2a1GG/w5bk0zpGzOHR/GuLzldux2dcegKEPAn8WkBWK3u6lxUwlKftnlnT5sOnZ6Li7EXnbs7Bw7ueP2jVZCm5qgh3p82DM6yJMJAYU6txs8tdF2kuUYJzmWi9ZfAAq0av5+kumC+zQhOeqvi+LT8ULfd+4GjDRvxlHyScKDp4T+KRuQHO2wF+A1zzkBz+OIXe4AIgdij0MzuMXzC2LHAkRJK3w7nlaoYuhcx7sQdKJM4VuOSsc+mkrEPuyFJWfCjcsdgYKKoXFQOM449NX1yWQuSaNRdUfAlN3S9x++uKHqm1xYm7QqGpU3UDGak86VHBOQVMy2qztmfKwBv/bRh/Cv/+zL6jWDirER+7AZsVMJi9D35TxyOMc6ZZGWKjp2wJSgFLZH7LxWDNF1PDkRoG5uZdVPXZmJvbDFc+wsYs2z+RQVM4wxSrLSBOVUzPZrxZj8OMnkDMfu+yUd+89+4Ev4J797DwDGsVdoXsnsAl5VHLuNogyFgSg7dqrASI76SuHYB3GIY8MYP3jXbfh7zzulOfYsR8AqAxI6XLLkaHRvfE8GT+1ronrnVzhir1DFkKPaDbkjOZebj0qZm6RiMiz3ZcC4WhWj5Y6ApGK0rM+kqwA5yXlSk1JMsbFByHBrmuG+c1fV91oVULPGh22qbG9B5VG5i7hmV8LtNz/+CP76gfOVx8xygbsfu4JPP7zifH+UpMpJA8DxQWTQcObvyvOhPAxJs/oGOGgVPGULhG28bK9tpIwhepJAV1icB6CpGBtsTFKdlGcETx000l7awjl2laDkm0XAOBUDlCWPmRBqtW4Knv7qRx7C9//6p0qv2xw7HWO5rx17HJYDh0+ujlVpWlUhMnEPAjK+OIySTGpq/XICDZ1TrFQsQr2e5wK+56EfaiomLQKqtEVeJY49ltfw89/3cvzAXbfJHAHlFDQVQw66H5kTj55L4MvAV1XS0fpY9wpVVExJFUOIff4TglAy1coZFSqW5V5Y0Giitoco6fT7dvCUyfoA+Xw5ylfBUwfHDkhH/LzTywCACzUt28i4OqwueJrlkqojVUyvJcf+23/7KP7svqesY+rvCSGKcgX6NZuKmRYZr4BE7LISavnYHLFnwq0Aa0PFJOozNVSMi2MvHDtRMoqKCST9CTDHPnU49mL+cwdu6tg7uWOjEdoOfQ++xbETOgXKKDAXOiJeFzwVQuCP7j6LrzpKhnIqhpd6PdzniL28nU8zrY2lYzeqYnggKskKJOyVtpk62u8VtS7YJGFUDKGJPBcGAifHPij6eJJpxC4Dbz2G2H2PqBgTwcnvycYK9gTcZDww8fpVckdFxewCx06TlnT840Tqzo8NIsQ1dUGmFhXDs3lTCp6WELumrYhjN6gYy0nceeY4AOBCUYf83rOr+PkPP+BcaOjeLffCygQln/mwwJ8teDplBeDIOKigeA4HIDYVkzIqhgDXqkMZY9c6UvEkB51RFzyl96pKdMhaOdWO/cbDchdHz5V3UFqr4NjHSabG0qRiTne1YlpYypBKyJoXUCW5quazQggVOKkLnj58cROPrWw5Jzf/yUmqE1uWe0wVE5RRdZrpoKZG7JyfrlbFAHKSEPVU1o4XiN1K287yXOrYiwXBqIHhlakYQuz6OrQqJvC17n1jkiL0fUS+uaWn8wg8d5yBB/iI16fv21tV+uxubGFpKz2IAxUAvbI1xbEBKwnhmIiu4Gk581QvjkJoqo5KCgCWTJE5EAC464zMJbiwLnd4P/ruz+E3PvaIk5qhcXSoH5acDcVahowiJKQaVeQY2CbjQs3xHJ5hzO+bTEbTahwqDnfFEUAlEJMrjr1QgDmomPrgqRnHMq9HOBUxgObYbzhiOnbe87RaFcOpGHfwdLcyqOts4Ry7maDkG9v5urT+LGdUTM2qT7yic8sozIe14aBiosArDSxec4XO11DFtOTY5WB3Tza7/ZjSsfuekc6uVDGBhdjjCsReKBsIqW5NM5Vq7UpWkoqGsvPgqPLyponYyyUFKMC8e4i9FwUqyejqVoIjg0jRFa5yCBx9A5KK4a3xbCoGAK6OyLEH8BViZ8+1+D5RMHcViP382gT/9HfvUYFel5xRO/aopIqh+3nIiv0A7XXsqUWzyOu04zkmqk8yoe4BUTGxhdhdPLtduC4oZIa50Pdrkpjzx2V1OnYeyLXtzjPH8c+/8/l4/QtOA+DBU0/ttJQqpkTFZDjc17s/Mu4rOsTewszgqbmdjwOfNbuwuWjt/Ood+4Xi98pcq61Bpu2ZGTwto1XOfRM6N3Tszlox5oQhTW2VKkY5dlbEKBMobWtlcSOfUTGkirEQexG/sDvcb05SFVRyyd+CAs2XEDsrB3xFIfYKx06IfRccO6GqpShAP/SxPk6wPklxdClWVIyriFtJFRObckc7eAro7XsvdHe7pzHwytuP4fgwxq3HljCIA3zkwQu479xVvPrZ0tG75IzKsfdCjJMcX2Lp9nQ/OeCYVceeOKgYszWenGNTw7Fr9Lo1TVVyG8AQu8Ox8zIYaTE+I4sWo+PU0agqC93JsVc79jj08fbvfB4OFfeLlC9RIOle39MNaGywMU5yDHsBQt+rDJ52qpgWxrlcXgQsSYuiPQ71ATloomKqHPv6OMHdj19WjRSqaA+gyFgsHvYhYwKVVQekTOC/2TbzlCwoJGA2bUPIQEuuMvV9GTyVTjpnx/cZYr/SgNjtRsijaSZ3RtbORD8XlN4D3FSMnVymP7t7CUqjaaZ00v04UAlBVMgNqKBiCjQWMyqGnF+Sm1p/Mtq+x9yxWzLFKPDwk296Af7wbd8Mz/Nw+lAPn330MnwP+L5X3qI+ZxunYgDgR979Ofyb//IVea7FGLDVWnQuTY5GCOmwy7Eim2O3i8LlOLIUqnMmEAEwxO6gYnh/AtpR2t3IdEOLOiqmhmNn9XKqjKpQji2KLPR9ladRomKSDL0wKOJY5fnAz30vbeEcu9rye1YRsKIGSuCgYhSqbaBiVrcS5EInr5SQN29cwRo08AkUOxx7munOQ3TsqkGgv2O+VlUrxq5ap/XyLBBVoGshRCH7lPfO8+RApeQablGhsqEqd5RxuTFJ1QLquse+R1SMeZ48Y5UWk6oEpdFuBk+TDP3Qh+d5WIoCPFMoUAzH7qBixql0wjS+VN/TpAgSuhD7SCN2JxWTZFiKAhwdxHjeDYcAAKeLAN6dZ47jlmIcuhw73TtC5RfWJ7g6Ig09IXYe+2GIvYEa4FnPrtfpOtJMGJ9JMoGlWFYGpWune0LZ2TZi58XwZK2j3Fkul8s3x0mGN//KJ3HP41eM36J559pxUaG2OiPgpzJPi/8P2IJg1y8apzl6kVzUO7njDkwF6QKzCBjxe6oejEETyH8VFVMRPKUBPygctT1AqMgVoBs0+J6Jdu3ECkAnXtDf9H0y1/nYv6G1vW56KFaIPVevq+CpSmfXiJ0KlgHy/G21gELsuSnlk4Fcv7Qz4bGPyJGgNJpmUk0T+IqKqdKx61oxuxA8TTJVT30pCvD0VSlDPTrQjVJcjm9zkhoLOO97SkFCukfHBtTQu3DsrKSAIXec6nMhO10oNN7w4hvUuHJSMcXv8HPatILOy0zppMr2FvGPuiQo5SBLVIyN2PV/9PkokM2+LxadxXosJjGMAxVfsa+Dfl8hdgug8eDppY0J7jt3FV94wnbsJmL/h+/6LH73M4+zc6tH7PT8adxpxK6/xxG77qoVSMdeUba3q+7YwnhJTdlHUr6epITYzS0coG9yE8dODmZQoDG7mFcuBIaqe06O9bGc7NwpurTmXO/rrhXjomLMY6uErEqOnTrBa55RBU8ZPZWx2ty0GNiKGDqe7M9qJt8QbxpZCxjPCHYtQJvTFMM4xNFBhNVNQuxlB0LVCIHdmRCjYusMEE8uj310KVKOOclyvOfvnjDKH2wU2alkqu/pNJP1zlnFzBsK1L020qoY31FSYCvJMLDuPX33DS++Ub3nDJ5mJmIHtApHcew9kyIE5DN3NTHnRjuWMqWonWMuynLVaTEHB3GI+86uAgCeWwSGgSL71KJi+P1QRed8D4FNxbDgKc0vO1fF5tjvPbuKr1Cd9eIZ1VmkHHtm/H/AFgReuporpXoWFeMqYLaXtnCOnRweIVFetjcMPJVhx28sT+KRv+Ee1DRAhz03x57lWkI2TiVi53ViADjRKiERKpQFWBy7M3jqQOyBV9pF0O/ZiF0GT030wxERvx+2hp2ORzLN0NeqGEBn5PHrVDsp30PkyG4cTTMMeoEsBkWI3UFLEe3A2/nN0yZJrlAy5TwAJhXz+MoWfvp9X8R/uf9p9f7GuB6xcyqG6BRSUnCO3dax83MAgLe++jb82ze/BLefGDDE7nDsxZimyqLPOjFQMR/NsbMaRtZiPk6ySqprqhB7efzTmOFNMfi/ceBjKQ7wWFHU7MU3HVbfPz6MS6oYm95RtWJsxJ5xxy7/thvqaFWMjt0QOJumzVQMUS9ljt1E+vQ+LTb9SCL2qiJgnWNvYVwV43tSEiWDPcLk2B1UTOigabjRgFkqkJKrLIHqdznNSpNdHqOcTk/BUzMNnNMYDmmlk4pxIHZrN6J2BlRt0fOMe5LlQkm46hG7r+RsvH47HUvWVTHRFiBjH5FfXoA2pxKdHh1EZSrG0j8Dks4Yp1mruimz2KjgtQGNugFJxVDwjBpx8Lr+m9NqxE6F0ugZEJ1C+vMqHbvk2M0peMfpQ/ih15yRx6ijYor7/W13nMJv/aM78b2vuAXjRO6wXKoYGgN0jr/5iUfx2l/4aG2LO9shpblQ40AqpooFgFWSjAJP3Zvjw1jdCwDGoq6uwwrI6loxZjKhQuxs91tC7NaumH+WksjqLLQQu4tjB3j1UV1mome1HLRzXvbaFs6xKwcSmEg0SUnu6HLspvOrSlCiB0B0i6v1m43Y+eQBJDIq63/NyD59374m13WSqSBohY7d7lKT5TJQyqkYanXHA2lAWRFDx6MkFDu4SuUN3KqYssYdkFUhB3Fg1OV2BU9p0hwdxBCivjbIdkyiZOJ95b++J2WDhBJJ288R4cYkMxw7IW2qS8QzT28oUtN18DRwUzHTtETFcKP36lQxcejjv3nRDUods1noxwFTx86rOwLAF564gmfWxs5dURXHnrLeu1Ril39umkq6g8bTC288ZNCUbsRujqESx67mji7foKmYxPotvXPIiqCsQvFtgqdWrRiuiuE2slRb0rEHcPUqCHyvQ+xtLGPIMGBOPMlyRKEu2mNqbk3HXsUv0iCgCeWiVAYGx56UELurbK9GHax287SeY7ePrVUxM1AxBToPbcRuoTdbw07HS/K8FBgEtPTSRXdRglK5K1KGQSwVIKsKsZeRIaHT44Xued4VHseppj8IWR5ZiuAXFBKg5Zi8SujmJDUc5ZIaB5le/AKLY1fBU18XoDNyIfISFcONzq/OsdOzpGe4NU0VunXq2ItzPFuUCXaVuK507JlG7ESbAJy6yRUVAwAvYjQM4K7JXs5mzY1eo/bin+bViF2hc6bB5+WTm+SO9P7I5tgrqBj6tx8F6EUWYi+e8yAKugSlNmaU7WUoKMmFkQLsqhUeN1ExFDytQOxZLhRtMZoWdbz7NhXjLtsLmEWtZkfs5VKm/FrcmaeA70MFlKl1nqJiggbEnsmJFvmyA5BKSyfn7UDsvnLsFmIv0OmxQYTVrcSIN/CBP1aIPTL+f17Gee1+cd2UPEP38IoLsY9TFXsBYPDfpPUnSejpQxQ8LTj2wAeBPj42uELHZbRTcvU05SICQDv2zUmq7icHHbxWDACcuyLVQOuTMs0zVcFTB8dOjp0FYOlzSVHWY1Dh2I8PI2xOs8qs63pVDHPaqZtj503i7cWJ162pMo3Y3Rw7B3X8nBRiN+hV+a9sAt+pYhqNODlSiQDS2ZMqxqV8odUzdLzHjQYBBRNtxCKEUFJI0rEftqkYB1p1BUxHDYjdfo0Qe5UqhiacmXlqb2ut4Gkdxx7oImCBhfbqdOx0ni7EPiyCp2kusD5J1cJgVwUEtLOdd52NSZqXOPYjxSKiK15KVMkR4WaNKobqFL32+afxE294Pl5+6xEAFsdOOnYjQSlVCqwqG/ZCg2P/X9/zBbzrE4+oXShVjSRp48Yk08HTfpmKIUqNHKqr8Tvvt2u8zoKnvKmL7Txpx/uimw4Z33/uKamQefAZ3ZC6jNipuqM7QSnJcnVeVYg9Y8oZ7eDFjjl2GpM2x64RO5/TFK8LOiqmjdlFwADp7JMsRxyagUIym4qplDuSKoaoGEeWJwW7xoksAlYKnlpJRIJxkYZjZzzePDh2F2K3g6fEsbdVxdBEsxUVOvO07NgDT1Ia5SJgGZaiUDmbzUmqFmm7eBQAHB2Wa1zPw0bTcvD06BI5dguxF4hbCIHNqfms+zxBqZDSHRlE+Gevf55CdmujBHGRDOUsKeDQsdu2xGrSAMDnH7+CLz55tYzYY4bYnbVizGdI1kTF8OBqlucGFWM3SplmAlEoOfbQ93AHkzoCwMtvOwoA+OKTV9VrtmRWZ56aO+8py6iukjvynYMdAE6yZrmjVsUUc9OiYqjxtU6g44jdtwQR8t+laH8ce3Xk5oAatXvzLLUHbYf1Fo7zXfJfomKqgqf0AJYqqRhJa/QjH5sT2TuTS8qAQhXjUIsAVuNb1rihShXje0zRQ4k/FSUF6NpUIgcLnvKKl60Re5GtmrIJIQNniVPuSOfh+1oVI4RQwbOtqaQyeJCXyqwmhebe9z01aQgdzVvyOEp08JSe89ESFVNw7KxUKw+c8+9KKkYvfoCuJzNJc+VcXfXYxw0cOyC3/7zODtEM6n57JhWzMUm1CMBZBMw8nguxE8ARAkawPc0EBrEGRzZiT4oM3O+/8zZ8w61HDYksANx8pI+TyzHuO3sVeA2K39EAJ8sFUq9QxVjlP3jPU82xJ8YY43ODdnq8+J5dpM02QvSKYycqhrX3A3jwVBeUk8HT8nwYxIHKCN5Lmwti9zzvTZ7nfdXzvIc8z/upefxmlVGRIMDkjilluE4VU1UgjEwHT6upGN8DTi738LXzsl57SRUTaKcGmJyqETxljt2pY8+EIcejTNsmKsZsZi1U6QX1mhAl9NaE2O3Ph0E5eGqqYnzlFMg2C3SqKCOW6AFoZ7LF5I7A9hD7e+8+i688teZ8b5xkilvvh+TYbSpGTkRyeq7SEX2rVSBHgxwVE+8eWCUF0oJScMU3uA3iQJX3BTRizYr7Rc9muVdG7HbLRvvc5DW6ELtw/p3mQlE5uSjr2Eka+9JbjuAH7rqt9Lue5+Hltx7F/edWjd8ECOCUETuBB95og9MrfOHn50rOecp2hVVle8nCSo5dvq4C+iUde7mkAD3nhaViPM8LAPwqgO8C8GIAb/U878U7/d0qy3KhAlGhhdhjpopx1ophyRUuo+2eomIctWIC38OrnnUMn33kMgBzu0vH4E6NBxj5ij5Sjj1wq2LyHFGodyCau66nYvjrWYGC+X1KM6HQIzmyqszTNM9VpinAHHtN8JRUMXQ8AEpbPYxDQ71jdN8p/h5Nd4bYhRD4mfd/CX/wuSdK7+W5wCTNlUMnB390SR5LlzImVYx0ehuuFohFmQVyjNxpUFkFQKN3m4qh57/UgNhlaWGNqim3gIYB0QU8eErj7JBVTpr/S+bk2NPyc6Fjx+zZ0uaDdy6KwnqX8vJbj+ChixvqntIciwvH3lQrxq5PwxcmvoPkfYnpHJsTlMySAuRL6HWF2K2evL2wrIqh816KQvX61a0EP/aeLxgxht2yeSD2VwN4SAjxiBBiCuAPALx5Dr/rtMxA7CbHzhG72cZL/muXArWtCbETZ33nmeNaeeBQxfBj2CoIMuXYo3JCE30/ZE6Seok2IXb1utCInWuocwOx625AtpmZp2bgLQz8UvBUyR1Z2z26f4Q4BzGjYjKJ2IkWoYmrqJihiY6abGuaIsuFoiJcQdeJRbUpjr1A7ITY6BZTFyRqB2jLQge9AFcLdG9L6WjxoHvmW8FTus5+A2IfxmaHJAIxisJQVAwPnpapGNrdlhF7dfDU/jtlHPvUcv5UEbLJeX7DrUchBPDlgmen8asQe1ZV3bFw1FZFSd6EhM8jAgia/89Lhe5ss1UxocWxH1+2qBjGsfdDicxpp07gccAQ+4e/9DT+7L6n9gTBz8Ox3wLgLPv/c8Vru2KcI+Y1UBLKWHNmns4YPKUiYHaCUuHYqdMNgLKO3TczQDnvxyWOXCrlVMVkZtPpJh27E7GXOHadcARohFqlY7cDdCp46uuG1TblxIPadP+IIx4wxD5N84KzDdX/A7pY2CFHH8kqy3OB7/iFj+J3Pv0YLm1Mi++VFwQbJduO3aWaWB+nCl0OLcrq+CBWbezsJBaiYOh6NWKHuk4AjaoYagbCrzVJNQ3Cq036nqZi4mLxpfwgVQSsOJ8TwxjLvbCWYwcsx55pKoY7J16WOvLL95DbywrF0P1F8+6UAROid3j5beqLYCJ2fWyO2PnrYwuxT5K8tKjZFviy9nqJiiGO3cqtsHXsAMsjEYTYA7Xwvf/eJ/Gck0O87JYjtecxD9szVYzneW/zPO9uz/Puvnjx4rZ/Jy26mAP6xudCJyg5y/a2zDzVJQWoVozpVHIh640///Qh5XgOVSH2rIzYOYqkU+hHVVQMNRwoHENRLqGyuqPDsVNAMmCBKDN4aupzuQXMUakCUkzuGFkLqJGgRJI6xZun6jiGY881x0wTd6soPUCItw1iXxsnuLg+wYNPr+NSUVXQlRSiUHIxCU8ekhP1lqOyPK4Lba6NEicVA8it+fmi7K+9KPQsxG6XFFCLTBuO3UbseW7U5gEkfz2Mw2LHkqFXqHHomanM0yJ4euORPg71w2aOnZUw5iUFbOdPc6eJijm53MOxQYTHL2/K+5Hr8WvWY9e7TCkEgPp/vqgYiN3BsfNEJTuY6zJeEiRUfkZe09FBDM/TyYXkxONCxw6UOz3R833i8hY+++hlfM8rbnb2XZ23zcOxPwmAR0puLV4zTAjxTiHEnUKIO0+dOrXtg0llipk4oDh2pmO3JYdAMxUjq9PpyVDis4V0lL7v4c5nSdReduzFdx3dXFxOqhKxFxl4pHYJHSiZX0uJimGBUrsIGG+TBlQgduaoFF9coEsZPDVjGTTPqYWfvAcFYp+WqZhJmiMXZdprlKRYigMmJ2xG7ITSz6+PsUKO3YH0OcICgBfeeBh/9eN/D68qniV37AQ818apqppY59htKR05dJrwdj12QuFNHPsgDo1G4ARieG0esmEv1Ig9NJ+xrYS66Ui/ErG7Yh+AHD+2rJY+TwtAExUDyGvWDpB2rgEL1mtVTJLlxq5NLiJujt0lm52kmiZqQuyAueNwJSj1Q9brNtd+hZ43ofnccux/cs85CAF87yt2jcwwbB6O/XMAnud53rM9z4sB/CCAP5vD7zqNl52lIGCSCtX6TiP2styxKXhKPH1VJ50811mbdxVty6gVGJlqr0bReIOKKTubXhiUdgYAFGUSsa08r6tOVkXFpHkuM089u2yvUE6GdwOyLXAMcIXYWXCQro/L75QGmaFwQDoeOiaheDr2hFExS0W1PPl6M2InZ35+bYKLiopxIPap6dgBWXSLEBRtxQGdPbo+rkbsJ5ZjhRhtjp2uU6liqoKnTTp2i4pJCyomzwU8T88BQFJFm9O0QKcmFRRZAXCF2CfVOnbX3xFbmMmkdLVArw3KEwBGiVseI3JlntoInVd3BMwYQZK7EbvdiLzOQjbG+bgA5LNYioPSbiDyfQXw6Hzouohq+5PPP4lvuPUIzpwcNp7DPGzHOnYhROp53j8D8BcAAgD/UQjx5R2fWYWlzDGpKDZ1PAnK0XSg7PyqSjeUJJO2Yxf62D/8mjN4zsll5QDItCqnjNjtgJ7nyQlXlXlKFSzpWnW5hByBHxjXZqORPKd4BMqInXWsB6o5dvW3lbVo5gsQYtfyO71rkq8R4uRyR8UxE8fO5I6DGRE7NX0+v6YRu8ux0yJRh5KjQKobbjzSxzNrY6yNNGK37xOpJOh73Gh3Q4uhCmCLGRF7FCAtnFvoe0VhtNyQoZIt90JsTmQjFJ5MBuhnSFz8rccGOHt5VKq2CJTROBmNSbuwFU8IaoPYeYnb1HLsuVfsMhmlqTM8fYP2AapVMZpj131Z2zh2V30Yus9LUYClKFA9UeU8lDt4AnikWad5SQv3k6sj/NBrntV4/HnZXBKUhBAfAvChefxWk+XMMfFuRoCcRJS45C4pUNvSjgEAACAASURBVFbMcKOovg6AmnROLngKd4g3vfTG0m/YGXNVwVOAnLVbFZNkQtaHKf6fS+jsNGygzPGmea6Cp3bZXhuxDytUMfw8AT0xooCpFlQad/k86R7Q8xnGodGJCUCJY5cJRBqxt+HYiVe/vDnF06uSGnEjdjOG4rK4cOw3H+3j3rMasXteORZxYqhL0toa6Z6F2Ol9RcUk5vVXGZWwoHr2gNwl0bPlRlRMP2KOPTRB0FIc4Hd+5Jvw8tuO4EtPXlXFwLiZOnYTLYeBLI8wNV7XVExTdidgIfZMq2JyIQChm7XQtaqiZsXCRd3S0jw3ELvBsU91pipx4q0QuyVxBPRckICDLUqZXlyPFI6dJLIqeMqe7xtfUvYXu2ULl3nKVR3KsatqbHrbyykQoqSJlqhsjVdsYWkyJA46J2gIfNhotY5jl0XLqmrFSIRGC4n8bDl+QIuWPWhzoVU8Kosv0zphQDpaz2uB2C1+1rUz4rVLQksZRIh3EGuHrQKqPVMVQ4idGm63UcUQxw5AdcxxcewqeFoTRItCH5gANx2RAdW1wrEvx2Ep6EUp5kCZirE59qACsbfJPAWArSRVi0SSFrJAr+zYL29KR03ZtJw+I/u2550EABzqR0bwkczg2HnwNKMxaWndU12/pakeC+BG7L3QDFry4ClXq62N0iJnxUff842a+UkFiKJCZ204dhs00vkALipGxxyOVCB2en7Pv2EZz94jGgZYQMdO2mxA3/yx5dgj3zML+JPT8aprswC6UUBY40Ab1FwlisJIc7acjUrNr1LFMAdqVrwrJ0K4tsDTYqvIOfZc6Pv2P77qVjz39LAieKp/z67l7WqEYDTaYFmw/+pPv6ikbYNY9/7csuR+5IjXRokqFtWL/FaInegXAPjqM+vF75W/N1a8dvUEp3t8+lAPvgdFxbjuEadiSnJHW8e+TY6dV5E83C/iNgUVYyP25Z4MtD6zNsZ3v+wmACbYse1wpSqmTMVQ0+mgoAdNVYymYprS9gG5mK1ZDrAXBgqABb65aybEPohDpWOXVSTDSsTOlUQbqhhbsypG+RB2HQSoBlGIpShgi5LOZj3cN/vcquBpMb7ftIdoHVhUx85QLKC5a77aOp2yL0v9NgdPy4idl6WtM3JqrvZiNCB8r5BOVpQJoONxvtpomMGvzeHY40A2LJ6muRE8pfRqVa1uGOP1L7zBeR2u4Knm2KuDp1IVI9+bJBl+77NPAJB0z3I/VJNPOfaCXuAV+w4vyWHZj9qVPF3ZmKr7SL/jpGJYtm+V0X1c7oc4vBRhfZxgc5KVNOyA5dhtKsbWsVuNNkbTdlSMqsk+yXQWZparOkDcBnGAZ66OkWQCd9DiGJKevTxul3syK5KraAA3FZMxZVng2Rx7rp5rW45dlQigoKvVyCVgtKlu9Sed/zjJEAUyYMl3HKSMm2a5AQg2Z0HsFhvAXyO1FgXTKXcGYIjdaiLz3FPLOD6M8T17pIYhm4cqZk+NB/9oDNGE1Tprs7sP+UHqk1oVPKUB7nnl0rNcp11ntv7aFTyl7XdV0hGgKafY4Uz5ouWSO/bY4sIXBJqMTddAxyOzFRUy89RE7FxXTWie0MtPvOH5+Mi/fC16YVBSxdhyx7VxovrI2p3fq2xlc6IcGZmLipm0QMl0fsu9UDmOdUcVT8Dk2EvBU0XF1CP2OloI4F2UUnWPNWI3j7ncC5VTpsqKXAxgm1Zy2L1Dy4hdJ6DJ2vy2Y1dUTAvnKZ+r5sABE+nTPJW/rVUxvKNUFPhq4eXn0WfVV8mIiplVFUNGACwKPPQjrVKSC4meG0uRLviVC6laet4Nh/D5n/n7pUqXu20L59h5owga2DYVY9ct5zSK73s1wVPNmVUWuWrg2HlnHXm+5Ug9oTByupw24sczEpR8nWhUVwcH0GhRCJPCIeTT1EmGjkdmVwaMWIcbXj0PMEsKUAeiY8NYqYcoG3LTUsVMUtmvU1IO0rH3o6BVB6WVjSnuOL1sOIc6xF6vitHxh8P9CGujpJKKOTbUUtcyxx4Y//JkOkAi9l7oN+4AaUezlWjEnhbVMG1wzM+RO/bId09zWkA3rGYbpo7dpBSVKsbSumvpX/PY6rPgqeLYWe9XO/dCl0jQNfCjwMNhC7GnuXD2iZ2Niil8C8/j8GXzEM/zTComM0sBH1mKDI69yVfspi04FWMHT/XrRgcl5nTqgqfUNxWQ2tRpVkb9Tc9qyZLpGcFTVphfnqe7K5L8Xo6gFxrBL4Vi8vJOgqMRG/3Q/ZoksyB2xjGqBKXq3YO6x4Gmj2iQc7qBsiFtud801QoHQpJSgdAmeDrBqUM9nD7cw7krIxzqh+7g6dTcMdVd97AnHft6kaB0YjgofbYXBkX2ZlqZoFSmYopzSbJGGgbQ9240zdSzllpuAZvXJ8c+jAPcdKSvjh80IvZqx05BUtVIJZBafwOxp0zu2Baxp9WI3YgLZZqKGVqI/VA/wvp4XX0vzQTrRavPb2MbVAxfDP/BNz0L31gksS1FLHjK2APA4dhbzLPdsoVD7DxzkmR7NPm5ZtdE7PJfLYV0/zaVJQBQahZhp3BXGW0FyXGZiL3Q0zIJXBXHrouAac6PnIUtQePXLn+flfv1Nf9LC1WbAcc/EynE7qJiLI6dJSgRFbMUmfghDn2mDS87dtIE2w2CXTZJM6yNU5wYxqrX6C1Hl5wlBcZphpglsbmMHNOhPlExibOhCtmJIdVyN3+THIymYuTr5Jy3WMOPOhtE2pkZtf3TzBE8lb/33NPLilOP2LOy7ZAV8CObptop8RK5AFRuBX8uRkmBlsHTcVKD2AMZE5AlsDliL6SfBcfOHSmdhwZWLHi6DSqG39sX3HgIby44cp4wlma5sQB0jn0HZpTttVpZcSrGxY/7nnQ8dTp2jpCNDuotOXbVMo1JoshUogXbnlepYmhgcImhQiNWUSh57fq8+AAOto3YGRVjcewmFUMJShRcRi1ip/PTqhCdoEQOZhbETl3vTyz3cMNhyXnffHRJpqdbzn00zQwH4jLiTIc9GTxd3ZKO3UXFADqAat9Tm2O3g6erW0kpa9llSwqxpwZlN0nKjp3OkccbeqFfWnTI6hC7HfvQZZn9EhVDRfiA9nJHhdiVmoaDEQ1mXFTMOMkQhT5OLsdYH6eGdLLnmCOzIPaqBtZkxqKUmYj9MHfsYn+pmIVz7FKHba6qY4uKqUpQUpx2TfCUy524hpd+o6mAT89CDC7ETpM19OsROy9vELBmydzZ1QVPAVNXPs3MqnV1xrfvdoJS4POaPNZWnUnVSNJmByvjwNdlBliCEjl2xbGHQaPccWWDHLvm8W8+2i+u13zQk7QZJets3AAvvukwnlkbY3UrqXHsPeN7ZPSseiyeAuj7RPRRk5GD3Zxmpe5LVY79uSxQN+wFlcHiNo69xLEHhdzR0LHPJnfshQGSoiS0a/xyOiQp1F38+ramGeLAU/ePEtSSLFcNUHgDcM2xt6FiNJByn7tfWab4yFJklBRoip/spi2cY8+YDjvwTMceMqograVi3J5dUjEalRrKGsqsbAqe2o7dIXfsR9pBuhpUAxqx812Ia5tJNUP4dpsHiVyIvU3w1FVSwEDsigO1ELsHRsWY9WDI4lBz7L0wUJwtBVtJ7hiHfmPtaprUJ5djPP+GQzjUC3HrsYFxvWRteozStS73QvxP33IGb3zJDcX/u79HVExV8JQcnefJoDEBhIvrE5xabnbsSu44zYxxO0qy0liksrIvvFE3kX77dz4fv/SDr3T+NlExLlUM8dm2KobAiC2JpM+1yTztqxK3mVGPnYxEAkFgIXaLY6eFnEonp5koiReAWRE7zTf3Z3l1S2piTmZTMW3m2W7Z4jl2R9leQrAGjcKpGCtB6fJWgn/+h/eqbTzZNMvR44uDo/Rv07iNCt7crgAHsEYPURvELrP8eD3vvsOx05aPj8PYQuzKsRfH91tsEetKCvDqe7xsr++h4EbrqZg49FX9GKKbkixXDoYQexT4Tq6cm0Lswx7ectdt+Ng7XqeQ6DTL8e8/+BX85888DqAoV9CgjIiL4OBSJLNf/8NbXol/+M23V+r9qflCpdyRUT+SBpTVOS+2ROx+saCPpqmx0xxNy1TMy289gnf/47vwuhecVq/dcnQJL62o/12F2KeZ0OULLPUKFUqrru7YjooB5MJbh9hD35cJSqmp+R9NU0SBr+7fxXWN2F3gZ30bqpgqp8zLZXD2AJCOfWOSIi3yDPYTsS+cKibNmComsFQxoX7dlaBEVMwnv34RuQC+5xU3G5MgSYWBkKfOxaH+YXmeh37oMx6uLHfU5W9rdOwZIXa9iNHCxbeZpHfn6M1EPxpd02Rs6v0ImIil1Bov4EWa9MSn72i5Izn2cvtAomKoTPIkzRXCJ8ceh76Ruu6ylU05qU8sxwh8D8eHsVEb+/33PoUrm1O84tajsnl0A2KPAh/Dni4fsBQH+Hff+7LKzyvEXlUrhjkTv1BkrY1lad02jh2Aartm1x2yWxp6nofXsvHcZFHgYxAHRlo+IJ05OUibY4+Ceh17KyqGnC9D7LE1ZuWxPNVWEbCDp57h2EXRpGPJybEnpWNUmSt4yk1Vfk1lXKEfcY5dnt/aOO3kjrMab+1GN86uFVPWsct/fa/IPC3+397mc1VMHPpW6d92jh1AqZ4EGVEDHLH7vqfrurDBRCUFKKhkInZzweHyMKCMEnXm6XYRu+mkooBnnmq5ow5qFxz7uIJjZxSLROyyyww5GGo3GIfNiP382gRLUWCoVvR2OcOoUJP8+HvvReB7KkOwyqLAr1TAuOxNL70RlzamuMGq8tmzVDGAfN55LhR91NaxU4Cdr/+TJMfh/s4dx9GlCKu2Y89y9CK/yOLUSVEAFIig5yI7BOlAdTtVjAOxB2XEzoOnvE9CLuRxTgxl44sL6xOmrjHFC8BsqpjIAo226abxcqG1OXZA7lRdJR/20hbOsafMAWruuHDsLKjqSuKx61fb9dZ5MIS3huO/0eZh9cLAGTylydBnNbpVIoYQ8KF/W3HsrDqfrbiRn0NBxXDEbqJEWgB1glJ7PS+gETgN3OVeVJY7MoRCn6+kYqw6NL3C0a+PpayQ7nGbImCPr2zh9uMDI6itK0Pm2JymuPNZx3D/k1cxTXO89gX1TV7ectdt+ObnnKj9DLdbjw3wU9/1wtLrto4dICpGUwcnW3DsgKYW7WS3eTiOI4MYq1tlx77cD6Xc0ELsoQUiluJAoVeg3W6QxqdE7HkJmGjE7ku5YyK7H5nNX2S9ohPDGBfXJ2q+09zi7Cb1rG2zmwgVFeP+rGo0UtBPNscOyHGf77NjX0COXRgcHMB07KEeEGmeYzTN8L7PnzOoGO6wbMTOa2aEgRm44/VmmmwpDmBXr+MPmXheXrHR5tllJT3f6ICjVTHcsecy5Zk5NiPZo1jMOC86q46d/n7JzYfxB2/7Ztx15lipBWDGOMWIPZeIJSyp87O23VHgKVUM70jVaxE8fXxlE7efMJOH6PfXxgmEAF7/otP487d/O974khvw+hfWUxXfesdJ/INvur32M23s1mNLCH0PNxdVIgGd9UyOvTViL6hFPkZGc3LsR5ciXB3ZsaZCkVXEPgA9jsPANxbRpSgwSwpsA7Hb85LXfMryHNMsK2SbJiAA5OJ4cX2ikvZ4oJ6c7vo4RRw0Z/nS9QHVHLtqxl40F7EzTwHp2LnIYz9sIR27Cp7aHDtzgmku8JcPnMePv/c+PHpJ9lf0LWRrI3YqIgToHoxkms5pPseliCcxFEgiJHWE2eDY1RiErpMXAQt9v+hj6Q6eVlIxbBGczODYbXQkz93DNz/nhBEgTZjckXe2okO45IV8SxwFMng6STOsjxPFrwOaihEVmcJ5LvDE5S2csRw7/T4h0WEc4jmnlvEbP3Qnfvg1ZxqvfR72kpuP4Ev/5xuNRYeynpVjnwmxm459nGRzUV0cHUROxE5tJrUqhpcU0J+1m1+0KwKmu2ZlRU1z3wAS2rlSdcc4NBPLqF/v6cN9XFwfq3nGxxvFdjYmSSt+HdBUTCPHngmpimGf4zXZuchjP2wxHTvVirHkjtwJppnAVsGtUVTc98ytJEeDWa7b6wFUK8Z8nx+zznhiDf1Gj/HqnO5R+marXkyaCwSBhzvPHMN3vug0jg/jIjAbWIhdFEoFtiOITCoGkAOVHPvMckfH511yRyNbtbhGO3AK2Ihdytaeviq7FVEAin5DCHe9ekD2OJ2kOW4/Yda5JsdO3YGaJI67ZXbpAp+omI0JoqCZ7ycjrtmOG80FsQ/cHDtx2pTLwalIPgcGcaicnE2pVBlvouIK/isde+AXJQXyUqIVja9TBWJPFc2p7/kw1jGpNvw6wCXTLTh2R60YoEDs+xw8XTjHTkFFQA9sCtJxuSPXv+pyuR5efNNhVRt5amlxAa7V9o3325btBeTgMupJ+GZddSWp4uoSSxlDJQVecvMRvOsf3aX7k7LALJ2XPaFMKkZz9Eru2IqK4cGs8jBRTbuL+8bbBgLcsZedqh0ou+P0Mh6+sGFUdgRMzbDLHl+RTSVsxE4BZ47YD4IFvtxlXFqf4MSw11oOR8jZrnE0F459KcbVrcTYFSWpLl9tUzFRYKJromJoMWhj5Hwnaa7b7bHvcjkztQWMQ7M0gnLsh3q4uDFRY5urngaOgHqTNXPsWoiQ5MLYoRw2HHu7ebZbtnCOnQclVM/TJMd3vfRGhcyCQHZQomAhd+z/9ntfil/4/pcDMBE7DQydeeptWxXTL1WAMymVkP1tZyQCsg2fdNjlx8OllPJ7BcXETsvOPKV7QkHm2ROUyp+3i4DxNmH8Oy60bNfefu7pZWxOMzx8cQOHGcfO+UyXPVE49mcdtxB7RFSMROxtim3thQVeQcW01LCr7xVAxabr5oEIjw4iTLPcAAvTTCAKTSqGzw9jdxhrjr2qiqRtvFE5lR92I3ZJQY0SmTHMxyH9ffpQD0kmVN/bPhtbvOVjW8QeMRrIZTQmZcDYXMyopaOmYlodclds4Rw7L9vr+x6GcYBXP/s4fvEtr1CfiQixFw5wxBw7wB6Oo+50zFZsF8fe5mGZXVYEIubAeTZpFcfOFQi29WObislV0ST6uC13pH+3WwTM5dhdRcA4QqlF7IYE0Fe1TcZJ7kbsFY79sZVNGaA8akoN6fle2XKrcvbLfF/KHS+uz+bYKZPabhAzr+ApAINnd3HsVAZ32AsNJ9wPfdVBqU1lRwCGbJdaQLp6jNKCRmWT7WYygA5AP7U6AmAWwOO0zKyIvZJjtzJPbWRP2acyeLp/7vVg7FFnMBvJfujt344bDveNhxgUHDuhDApk0tcCX6Z3u7qxm7Viyhx7Kx07o2KoUJDe4jEqxtc11jlidylpyEocu2DO2/eQZ8IsKcAmyXYbbTipmOJ9krnZ8i4KKi05aJDIQcWQcY49Vsiugoq5vCXVJxV1WjRiPxjDnAdPX3Lz4dbfoyS2EmKfE8cOSMd+81Gp4CEkGoWajiTJ4DAOjSznpThwSv/qrIzYTXpHJ7r52EhTjKYZbh0MLAmu27FHPiX9CcQh0UmiVdYp/90qjj1ioNBFPy33QmxNMzkf9o+JWTzEzjsoAcCzTgxLQSqq86KoGCuVnmqCm017bY7dqhXTsrojYAVPC0kU19iHfjvE7pooUkrJFxwdfVc7EgcVE87o2Okznuf+PClf6FzthAxytgOHKsamYk4ux4qC4Yi918CxP7GyVQqcAg7EXlHnZa8t8CS1sLI5nQ2xO1QxwPw4dgBYZZJHahEZB57K/NXNxwPjuJxjb5MfATAde5Kr+eyS19J1b00zLPcCJ8d+uriPT18dq9dVpnSg5cKtETtTkblM04PCiPep90PZp7cr2zujZbloRM08Yw0AxlOTigHKyS82Yrc7KLUtKQBIuoRnnnL1C+/AHgacYy9XbHRy7JHvDJ7SdQPlsr1AwbGTY29xDeSY63jTMPCV3DG3ypTSouSiQXoGYpe6aELtXO7I0ZFtQgg8trJZCpwCejt+5YBx7L7v4dLGBFkuWksdAbeOHWgXK2kyQuzUq1MI2cQjsqgYA7FbCqwkzbE6StRvNVkvKiN2F/UXFHQoUTGmBFf+bSN2WaZDz2GiTmZVxVQ5ZZJZUgGzclXPQAWFO8c+g1FxrDrj+leAcezsaqnwFAB85KsXWO0SrophyHiW4GkYYFo8XB081Q6Wq3e0Y2fXyDTDtnH+nr7nMx4dsKo7+nRPdA3sWaiYus9Gvtcod2wTPAV0KzeDiqkJnm5MUqyPU9x6bKn0Hn1vdctdq2a/LPA8hSxPH+43fFobxXtsxz4P1YWiYkYJPnDvk/jyU2sAUNAYJsfej/yStLYXSY59ZWOCEy0XK14EjHTsHBTYtWKo0Ql37Lw37SAOcO5KQcUwsMTLXrcOnjYWAZO/M7L8Bb+2SSpLLHclBVqaEKKVfpc6KLlUMWQ0aB+5uIF//Nufw//8Hc8BoAeA1NByxC7/bfOwyJkRIoksJQwP0NDvGU09ajj2XmTLHXP1e74DsfuMf5+lCJjaYdR8Ngx0PR2+wPDvuRKU7DR7QDv2tsFTKhjm0oLT5Dxoqhjf93DuilTy3H68vNOoMiopsCuIvaBiLq5P8Mt/83W8sZACU8Ywcewbk1TJRvnus1fs2i5tTPHcU+0aNnue7AY25qoYg2PXY3aSSh37sBcau0eeNHfj4b66rzyGtT0qpkDsDRw7AcFyuWYfG5MUQrQDgbtlC4XY69Qi3KiDkgqeOhw7FaIiB3H2shwYZq0Yl9yx+Tx5hTkVPGWDNWSSKt60lyytuc6lKDDqjGdsADlVMWwSKmVPGyqGoZ4qo9Zl8vxza6vcrGP3Pb0YveyWowCAm49oJFvn2NdVt6WyYyfHkQsodcdBMP4MzpwsxwYqvxeUE5SA+XDsS7GU6N3z+BUkmcCTqxr5xqGW/G5NMxWr4GCBksgurk9woihh3MZ6oS8Re7EDr6oVQ/WGbCqGz43Th3sqnsLpzTDQZa/bB0+LYzfUiqGYQ7nXbVBcV1ePvbWRw2vagtJEsBOUbA5Y0jXyPR180YMqyYRK3MhmqBVDtTDGRalVjtKNPqYBR+ztCo6VOfa8RJsYzawdlMoswdO6z8oMX5I7Wjsin6gYV+YplS3W5/ma557AJ97xOjzvhkPsc8WW3cGx242vy8coFpYDEjgF9LM4udybqYJk6O9eghIg6Zi7H7sMAHiSURomx84RO52XZ8j/2lIxAHHRmaqpYiJ2zXNT8uEwDky5I0PgNzJaK2IcOyVZ2Z+vMy5Fdhn9Di87za0f+Yp/7zj2lkaouWklpKALOXRSkXgGxx6obR4APEOOnaliANZIYqaSAhqxJ0WXFT5Y+a4grJE7upBmmWPX+nFyrEYza6/soNuoF6jbVFTn2FmA2ZY70kJWp2O3n+Ntx+0M0u0hdkAjNJcqZ7+MfMCzT7anYQANVHYjQQmQdMxm4aioG5HWsctjbk0z9SxpnIXWbohq07exftQGsXugtWzIqn7K9/Rxb2C7PA6itsOxcwDmsjiwHXsZsY+TvJSJvde2UI69Tt/NjR4K3XwnFVOkS5OTPL82Ll6niLqpyJilbC/vdJQWckRdRoCjd71tNMsMV8sSqVwB7SRG04z1Ii1z7C7k3TZvIvC92lZnVFYVgLpO/h5QHzxtWqB5XQ7bCLFXId+eQuwHJ4xE9+dZDolm0/eSTKjGz/r1+UzfIw41SxQSx14g9qlu6M0ltByxnpxBwtkrgJWqMePIPOUAZLkXOnXsAIxa+HHoGQsDjaH5cezkW1LjXNV1qaYoHWJvbVQoqzl4at58VwEv4tgJsRMY4rViAJaAU7zfZhHm7bkkx854PxuxB2WOXdW2rnDsQkBVPXzk0iaeXfC1CrE7gqdNCUcuszMCXe8TFZPnZaoLqOLY9SJXZ3WIXXdbanDsByRwCuhn8ewZ+HVA1z6y6sTNLWWdsk9tx0ngBwC2JhqxBxbHTnZyOItj13rvko7dkf1JXa14uQGyG23E7lTFzIdjp+TGzQZVTFePfQbLWlIx9D5pb8m4Uyb+0M5q5LViAC09nCVBiVDqqKheFwUa+fKkJJNjb6eKUbuBaY6LGxOsj1OlRlAc+1wRex0Vo7fqtryLFo9ZqJiqz5FjT7Ic7/jj+/DIxY1GKiY+gI6d7s+ZGRE7lRSgnRyN43khdpI8fuPtx9RrimNPGWJ3qGIMKmaG4Gm/0HvXqWK48x72dDtJwIwj3XBYLyimjn0bqpgGjp2SG0nuWG6JWOjY91nuuFiOvSZxhxs9WGqYTFZSxTAqhr/OfyNhOm2gJcfOMuuoYUbkmAyNqhhX5inrF/nwBVln3nbsVTp2slkQe929jgJPORt760mxiqXIETwNdAnjOrOrOz58cQPvvfscPvLVi9gYpwh9TwWqq757UDTsAKdiZuTYC4UWMTG8teI87OhAOuRvf95J9VocUKMNzbErKoYlvXHne3wGjl3pvXMzgQ/gYMSkYgBWpMtw7Dx4yjNP9fm1LwJGu8nqexsHPqNiyoh9ShRTx7G3M82x139OUTEWYrc5YE7F6NdNtJBYiN1r8bCW4iKJIdHB08ARPOVIxVTF5Op928iRjYpqiADw3NNExcjPcLmj7wietkUSgW/WwC6/7xnBZVetmFrE3kTFWIj96VUZB7m8KXcqh/ph5fM4yFTMLFJHQDqPLNOInRz7vMrCUtDzO1jbQELsRPltTFImd9TnRYj4UC8slfaoM1kBNXdnnlpN0QG9QIfW/ASA04c4FeMZVIyWO84HsQMStBBip0xUfl3A/DpcbdcODpxpYXlLxE431K4xwu9zJWJnmacAU8XMVCuG6BJH8JTRG1wtQ4j9w198urYhhoHYL25gEAdK7kUTxFUpz/VakzVx7JHvG8FlV4LSjqiYwETspLG+vDnFaJpV0jD8WXv/5gAAHEdJREFUGAfJsQe+h1OHZpM6AvJeJgyx9+eM2N9y1214/o2H8MIbdWEy2Uhd/v44yTFN8xIVw0HKLDQMoBE7gFLmKV84yOieBb45PwH5rE8MY6xsTqU4gaH67api6kppRIGnOHYXYgekPLRz7C2tLnGHW9VD8QxVTD3HHlqInSbVLHJHyqzjSpjANwuC0fiklPH/5fc/jxcXlf+qgqeAnGwPX9zEc04N1XX5hUTR4NOdcse2iN2rpW3CQGez5kIYO6laVYx1j6vMDp5SPZCVjSmyXFRq2AEmdzxAVMxrX3AKL7rpUPMHLVPB0zwvtVachx0dxHjdC2Qv2OPDGJc3p0bg0W5KbqhiQnLs7QOnAAVPcxU45b/pKTkljWu9U63Spd9wuI+VzSliJlSIWfmO9q3xWiD2wFfd2Uocu6qDky9uow3P877f87wve56Xe55357xOqspoK9p0w1wPxX4tVlSMm2O3C1CpImAt7pideRr5Zrbp0WGEYRzglqNL0GV7c2xNU+RCdwZyOT6ukX/4woaRxk1o3VV3g/9WW30tLzfsfp/LHc3a1G1a4zXtvPziemzHfnlzqqiYKovV8Q8OYv/h15zBv3zjC2f+HtU+ol4EKsFtFxwHVUvkGbtU+ZE4dlMVI/+eRcMOUPBUc+yu2kT0GiliAE7F2I69p97nRfZmzTzV363h2EMfW0m1jt0+//2wnXLsXwLwfQA+PodzaTRCzY2qGMdDsb8ShWahMPVdKyKvEnBm7KAEsOBpoCVYQeDhcD/CF372DfiO559Sr1N5UkBrtOs49tWtKZ5cHbkdu4N2MWiSWRB7nY7d53JHYdXUrqZidD2e5vMgWSoAPHWVOPZpqY1e6RjFfRoeIB37do3nVBgxml0IzlFxsijUaHzValiiqRKNiLeD2DfGKZ65OsaJ5Z5TvUV/c+oqcvDvgJY82hx7NCNiv/nIEp59cojnna7eWcngaXWtGLL9TFDa0agXQjwAtAsozsPSmsQdbi76wD7HOAgUYl/uhdiYpIgDX33ORuzZDBx7UAz4UZIVfRG1g6SBYG+ncyHLk5rX4eDYi8n1ladlFT7u2H1PHoe6KeWirGP3vPZBt0YdO8s8zYTZWKAfBQh9z8ltzkIlUCwEYFTM5hTLvbAVYncVIVs0o/s0SXNjsd1NxE46dqDcO5aXqaBzOTkjx96PAsVTf+sdJ5yUIc1Bvji75I4A8Mrbj+ELT6wWGdPa+c/KsR8ZRPjIT7y29jO8CY9Lx07Wcewtra3k0EnFWN+JQtkqbpLkODaMSrUdbLnjLPXYAWq2kSEtUFZVGVzeaINQQN11kJTyS09eBaAVMfR5HSyVDtGWjs2C8l5442Hc4iiLq8490PXYZQahHtRv/abb8fJbjzgXfRWgbsFrEWWW5QLPXB0j8D1cHSVIs9yo3V76XkhOYfEdO6HTSZIXoKGZLtiuuaiYqwUVQ6oY1SCd1WKZlYrpsQD6q599gmWblgPwvHcpD4xy+4E7b8MP3HkbALPeUxS6OfmdGN8t2I7d7OR2gB2753l/BeBGx1v/SgjxgbYH8jzvbQDeBgC333576xPkphx7w4B2bfHte9wrHMY4zdALA5xYjg3ErGvFmJmpbZ8VNbQm7rmqBgU9/Mzh2OsQ+/3nriL0PTznJKNiPGsBydgkrFhY6uyX3vrK2vd5PXY7eHrL0SXcctS9KEQzIHaS3F1cnyDNBV544yE8+Mw6NqfZwgVPt2u0YE7SzKBidmOrT+3x+lGgnNSFNVk/RitT9Jg6NojhebOXSaB6Rt94+zGDauHggHPs6jWH3NE2XVLAV01d2nLsbcxo7VhKUFoQxC6E+M55HEgI8U4A7wSAO++8UzR83Glty/baPJ2r6xI9nI1Jhn7kY9gLjdT1EhUzQ60YQDrgUVFSgG8Jbd5cN4UWpYQqF79NiH1lc4o7Ti9bbfD0QFPcevG2K017p2bWY2+faddWFQPohI+nrkoa5qW3HMGDz6wDqK4TAxxMueN2TSH2NC8luM3bvu8bb8EtR5dw6lAPtx2TiVQPnpf3mxZJj9EmNx9dwkf+xWtnTroiB/itd+ikKJv6o7lhOHbmtKuMx8lmpWLaGJ9z9q7TVYBvP2x+V7sHlrakYvhDX7YKF5HRw1kfJ+iFAU4OY+dKnNjB05aTqR8GUhVT9Dx16ckBqK3iNMtLCVVOuWOsz/H5N5iNDbg80Q5GbQexNxmvx57PkELd1KWGGwVPiV9/2S1H1Hu1wdMDmHm6XeMcu89UMbshpxvEIV73Qil9vL1w1g8UXZWI1tLZzPL4Z04OZ46z0W7g21i2q++7W+QtG4i9rGO3jX+GdofzpGLiloh9P+WOOxr1nuf99wB+GcApAB/0PO9eIcQb53JmDstbomb+/uGlEFdHSVkVE5BjT3H6UA9vfMmNeM6pYel9WxXTdhWmvqeJkju6HetSJJsDb4zTEmJ3XacM8AJCoBS5Jx07/65dtneeKI92QwCUFK+NeZ6suteUeQro4Ck59pfeopNoaoOn1xBiD5Vjz3YdsXM7shTh+DDGoyuydAUtknpMbd9ZvuElN2BrmuKVtx1VrwWeKa/ViF0/wygo52rYZmSe7gJi54tKXfB0P+WOO1XF/CmAP53TuTRaXQ0VbvyGHupFAEYlKoYj9tuPD/ADd91mvF+mYuTrbXnNo0uRqm0tU/Pdk9HzPCz3QqyPE5X08P+3d66xcVzXHf+fmdkHuXyLFEVZ1oMSrdhym9RR/GgdJ5Ec1zZcq23SwG1RO4gBoUUrJA3a1KnQwm0/JUb6qUANpwkapG6cFK0RfwlgOwha9IOTOobsOI3lR+ogdiQ6kR8iKZHikrcf5t7ZO7Mzs7PcufPYPT9A0HI4mj26986ZM+eeR1X6lcMWBhG5bwPrGzi4w6/YwzZprYBCT9UVE8w87eLaNS0ENA61efqzt1cxWnN8vtx4H3v/bJ4qBbe2vinrs5iLigmye2oYb674WwwG19RW2D5ax/Gb9vuOBUN1vc1TnysmvswF4M/dMGGxV3RXTFujjWJsnpbKFbORMDJFtwTVzd/mitEs9rCnuVq0wVoxSY2UXZND+Im0dPSypGFW6mjdwdJq00t6UGVdoxaG2kANumJmRmpe13bPqiK/gk/bFeNrtNHF63hXFnvT3TzdPlbzNuuADq4YeYOFdXAqG47mY9et2iwUx17pjlHNrQFtbaUclWNZQR+7VOzaHOp+8yi82HXbwq/u34aPXLPL12WpV/yumGiLvbRx7FnT2jyNn1h9cYzJWtPBe0BvcRWm2IO1Yrqp7ggAuyaHfe2z4vzKo/UKzq82cWFtAxYB8zMNnF5civx/1h3XaglGIvzNsUNtm7xB6yppZcckuOVkteqOXdzoVSeZxV6xLaysNfH2xUuYHK7CtggTQxW8dWE9shY7AHz4ylm8uXzJ10O1rDiBqJhqQMGaRK0xPeQwDYs9jGA1UaU0G76omc6KXY+K2T8zgi987N2pyqkX/mpvtFGMzNNSKfbEHZQ0BeNZ7BFRMQBCq9JFVXdM+hS+fKoV6ucW/49+fXYt9nWv5rUKE4xSlPWqjf0zI20LXN8oDFp16oZJUa+jYpHXF3ZTdGex/9EH9/tCNaOoOrIjzYV1zEklPdWo4q0L8ZmnO8br+OTNC4nlKTL65mlNt5yzsNhlGz99bYUlE6VBlMU+ornT9Nj5KDwjxkCcP5Dcx17azdOsSRpy6LPY68piD/exA+EbK8pa8IpcJWykrVChYoCs4hhjsY/VHfzs7VW3S03Nxq8tTOPUT9+OzJp87+5JXw3qMIKuGFMWOwBZoKq7xgJ337A30Xlq83RtddPbU9jWqOGVn6/E+tj7CT1BabjqtBWqM4my2P2x5u7faVukNgXL97b72IcrdsdIp6js1LRo5RGEJxxa5Oa9ZDA9kZTqzkia1q+/zo15Pnb/OfpTvxaiQEdrDoartlcqttuOKHpjZlu32ENme7RewdLakmexf+jgdq/aXhgP/E7nV8u2OHb1c4qLTQ8J3RRmLEiVSPbOxXVMDLnZjVMN1yXTDxEvSWhZ7BuwrdbazcKHq7o9DWtWs97MOk2CJaeVQp8abmW1njiygF+srHW4TueQyF5oPVjbr09EqMnghkInKBUJVd0xaWs8QPexR1vs9RCL3bIIC7OjeFEmZ7h1V5LLOjnsVnBcubThrzgX6Ypp4uKljdAyt1shGO5owmJXyRnqrcaEz7fquAWXlteaXgu3uYk6tjWqmdUoyhulSNxaMXqElfnYh8nhCkbrjm8DM6yPbhoEwxgP7RzDv9x7Ha6f3+Yd271t2Iuvj8K0K8arABvx/69XLKnY84tNKVVUTLPLZtZAy8ceVDr6a1qYxQ4AB2dHcPqs26VoMyR7NQ4iwi7pjnEbavh93TpKsS+vNX03UC+0xbObCHe0W5YkYManWHUsvHXBDbdTiv3EkQV89d7rUv+uoqLmbHV9AzZproAM7l4iwuE9k74G3CbWkrqefu8SEW5cmO56XUWV9k2LTs3YVfkCQ1+fiFJZ7Em7GOmWjNpgC+rkTj52ALhidhTfePo1nFte69qHDLgbqCq6JVjdUWe07hYhO7dyCbtiim51g2253dQpYF2lXVIAcMsT69+RJhXbgpx2jMu3r6lGtav+mmVHuV42hYzjdtJ/+4rjS/e8z3f/mIqKsa34HrtJOTAzgvnphjFXXbDhfRBVMprDHROStIOSXiluyKshHR0VE2mxy82604tLsixtdxOlW+yVGMWq3ioW31nFwdnuO+yEYVvhDTdSjWPXfL9pX1uhP4AnhgdHmevoys62SOvyk833By3mNDJPw9CL2PXCLYd24JZDYXUL06Gl2KMs9uw2t6MolWLfSBiZojYoa06rultbB6UEFrtSsi+eXYIQ7VZ/J9QGqt60OspiB4CltWZqVoZjkW+c7Jhwy62ih+EBZiwU3WU2MRQd3tjPBPvVtjoo5fOubyoq5vKpIa+6ZJGpdlDcyhXDFntCklZ3VBZNrWJ7k9DeaCM+jh0AZkZrmBiu4PTiclu6cxKUW8Wx45sj6GF7aXX8schvsZspAtZKnEn72gq/xT6git0OKPaEzcBNoVd3TJOH/sB4d81U6NRPQBmKXFIgIUnj2NXva46lPT395ySx2IkIV8yO4vTZ812HOwLAtXuncPOVs7h653gr8zTkKa9nUKYVFaOXMQDMFAFztPhq/TvSRJ8bFe44aLRZ7Fa+PtxgpFVaWIG3zKLi7XFEJRBKQ7HMPU8zJXHmqabYo9qw+XzsMQWCDs6O4qXFZQghug6vm2xU8U/3HMbMqN7TMSwqpmWJNlJS7JZsj6cwUV9Ebd6tZmCxE8UX/epnnKCPvYNiMY23lvPMwMmRTuGmSp/k6YoplWJParFbMvur5sS4YpzOrhgA2DZSxdJaE+sb3W+e6uwcH8KOsTrmp9s7zegKK6364VFNrdMuAga4zUqAVqPtNFE30fhQpRTWnAl880iEd+0YwxWzI14D56wxFRVTFrzN0wiDUEXFcIJSQpIWAVPn1CqW9vT0/96XeRpjsavNzJW1Zk8TNdmo4qm/PBr6O5/FnlKZ2WCyh+djT9GKUPsG5y+6zY6HKukvJ+XPHNSNU8D/dmnbhCvnxvD4n34gN3nUssozASdPOiUoee5fdsUko+mV7e18rmOT3xUTl6AU0w9R1WtZXmsaSwhpVG3v/5SWxW5b5HsVNBEVoxb2O1Kxm4gbVvM3PqChjkC7xZ43SoZOxbj6lU6tHWs5b24DJVPsmzJJKImv27bI54oJ+rtIazEW50JQ9bxX1prGfGZE5FntaSnHKIs9Tb+sZ7Gvmlfsg22xt89jnpio7V8mOsWxK9cu9zxNSLOL7M+KbcmomNbmWxCvbVaMj10pq+W1ptGJUn729Cx2y/cgUjdjmg8n9ZA4f9Ht/JRWRI+Op9gHNNQR8CvQIuwzmIqKKQudevZ67l+22JOxsbmZWLnaFqFWafU8DHsgqM2POB+754pZbRqdKGWxp+Vj3zlRx86J1uaano2bFirsTlnsadW50WEfu98yLIIybVUMLZX6SI3Emae8eZqMjc3kg+VYJBs/u3+HWapKacRFxQxpFvuIwXC7tC32P7/lIOSWBID2jkppoK55nn3sRgmLbsoTU9Udy4JS3JGKvZL/5mnJFPtm4tjZE0cWMD/jhhbWHCvUFaMmJlFUzKUNo3GpKkkpLYs9WHnOiMWuXDGrBl0xbLG3JSjljanqjmWhU6OTWkTARpaUSrE3u2iY/HvX7fY+Vx0rdBHWHAsWxSs75YrZSnXHbmhtnpqZkrgEqa2iHh5LnsVuwBWjLPYBVuwkOwuZXoNJGXgfu+dq4ZICqfB3x66OjAWPo+pEuGIcC/WKHRtlo1uhZhW7csWYKTXqGKgIqG7s86vrkQ/PXjmwfQT33rgPHzw4k/q1y4SJPIStwlEx8eGeyhXDCUoJsSyChe4HS1nmQVTkTBx631GTHXvmpxvYOV431hygVSsmvesrWd+5uO7riZkmFdvCX91xlZFrlwnHIqyhGGn8rdZ4+cuSBzXb1QlR96rSGXnG+ZdKsW+VesUOnQRlscehuxdMztPdN+zFXdfu7nziFjFSK0breWrKhcS4BFsd5olqcrKtUctZknzoVKvnyLu242+PHcL+mZEsxfIxEHfj/XceCi0gVZHZqXHoyUsmX60si1C3zDVnVlZWupmnrbExsXHKtGgVnspfsR/YPoL//osPeY1kBo1O4Y6NmoO7b9iboUTtDIRi15vh6lQdO7acAOC6X4YqbtfxMjdPNtMar3UtU3sDjIuJcNVeGFSlDsgmNlTskgoDodijeN+eSZxbudTxvOGqq9iL8Bq8VUy8yuuKfaiDS4vpjSJZ7IMOEeHzH303Du+ZzFuUSAZasZ84upDovHoBdrl7xYSPXXfFpNX5iQmnaBb7oPPR9+7KW4RYShXumBfKzVDmm8pEByVV9x5gH7tpnA71SRhGhxV7AjzFXuJ7SoU5pv1wUtcdZleMUVp7JHzLMp3hVZKAIpTh7BUTFjvQsiR589QsXoJZeZcgkyGs2BPQD66YRtXG9fNT+KXLxlO9rnpQDHEcu1G8PRJDCWxMf8F3YwKG+sEVY1t45PgNqV9XRWuwxW6WIpUUYIoPP/4ToHp5ljkqxhTsiskG5YrhzVMmCazYE9DaPOWbKoi3ecquGKOoB2iZ3YFMdvSk2InoASJ6gYieI6JHiWgiLcGKxBAr9kgqbLFngqnNb6Y/6dVifwLA1UKIXwbwIoDP9i5S8RjqgwQlU6ia7BzHbhaHE5SYLuhJsQshHhdCNOWPTwEodjrWFmGLPRqlcNhiN4vDJQWYLkjTx/4JAN9K8XqFQSktjjRrpxUVwz52k5go4sb0Lx3vRiJ6EsCOkF+dFEJ8U55zEkATwMMx1zkO4DgA7N5tru64CZQrhi32djgqJhuUxc7hjkwSOip2IcTNcb8noo8DuAPAUSGEiLnOQwAeAoDDhw9HnldEhvogQckU7IrJBs9i59RTJgE9vT8T0a0APgPgA0KIC+mIVDw8VwxbS22ocEfePDULJygx3dCr1/gfAIwCeIKIThHRgynIVDhUghIb7O20XDHsYzcJV3dkuqGnu1EIcSAtQYoMu2KiUZun3GjDLKaqczL9Ccd5JIBdMdE4FqFesThawzCcoMR0Ayv2BHCCUjQV22I3TAZUuKQA0wWs2BOgXDFlbmZtinrFxmidFbtpbC4CxnQB35EJ4ASlaP7kyAGcW17LW4y+x7PY2bhgEsCKPQF1x0bVYZdDGPumG9g33chbjL7H87FzHDuTANZUCbAswiPHr8c8KzAmJypcK4bpAlbsCblm92TeIjADjLLY2RXDJIG9xgxTAhwv3JFvWaYzvEoYpgS06rHnLAhTCniZMEwJaNVj51uW6Qz72BmmBBy9cjveWFrA7Fgtb1GYEsCKnWFKwNz4ED794SvyFoMpCfxexzAM02ewYmcYhukzWLEzDMP0GazYGYZh+gxW7AzDMH0GK3aGYZg+gxU7wzBMn8GKnWEYps8gIUT2X0r0cwA/2eI/nwbwixTFSYuiygUUVzaWqzuKKhdQXNn6Ta49QoiZTifloth7gYieFkIczluOIEWVCyiubCxXdxRVLqC4sg2qXOyKYRiG6TNYsTMMw/QZZVTsD+UtQARFlQsormwsV3cUVS6guLINpFyl87EzDMMw8ZTRYmcYhmFiKJViJ6Jbieg0Eb1MRPflKMflRPQdIvpfIvohEX1SHr+fiF4nolPyz+05yPYqEf1Afv/T8tgUET1BRC/JvzPtzE1EB7UxOUVE54noU3mNFxF9mYjeIKLntWORY0REn5Vr7jQR/XrGcj1ARC8Q0XNE9CgRTcjje4noojZ2D2YsV+Tc5TxeX9dkepWITsnjWY5XlH7Ibo0JIUrxB4AN4BUA8wCqAJ4FcFVOsswBuEZ+HgXwIoCrANwP4M9yHqdXAUwHjn0ewH3y830APpfzPJ4FsCev8QJwE4BrADzfaYzkvD4LoAZgn1yDdoZy3QLAkZ8/p8m1Vz8vh/EKnbu8xyvw+y8A+OscxitKP2S2xspksV8L4GUhxI+FEJcAPALgWB6CCCHOCCGekZ+XAPwIwGV5yJKQYwC+Ij9/BcBv5ijLUQCvCCG2mqDWM0KI/wLwZuBw1BgdA/CIEGJNCPF/AF6GuxYzkUsI8bgQoil/fArALhPf3a1cMeQ6XgoiIgAfA/A1E98dR4x+yGyNlUmxXwbgp9rPr6EAypSI9gL4FQDflYdOyNfmL2ft8pAIAE8S0feJ6Lg8NiuEOCM/nwUwm4Ncirvgv9nyHi9F1BgVad19AsC3tJ/3SbfCfxLR+3OQJ2zuijJe7wewKIR4STuW+XgF9ENma6xMir1wENEIgH8H8CkhxHkA/wjXVfQeAGfgvgpmzY1CiPcAuA3AHxPRTfovhfvul0soFBFVAdwJ4N/koSKMVxt5jlEURHQSQBPAw/LQGQC75Vx/GsC/EtFYhiIVcu40fhd+AyLz8QrRDx6m11iZFPvrAC7Xft4lj+UCEVXgTtrDQoj/AAAhxKIQYkMIsQngizD0ChqHEOJ1+fcbAB6VMiwS0ZyUew7AG1nLJbkNwDNCiEUpY+7jpRE1RrmvOyL6OIA7APy+VAiQr+3n5Ofvw/XLZtbtOmbuijBeDoDfBvB1dSzr8QrTD8hwjZVJsf8PgAUi2ictv7sAPJaHINJ/9yUAPxJC/L12fE477bcAPB/8t4blahDRqPoMd+PtebjjdI887R4A38xSLg2fFZX3eAWIGqPHANxFRDUi2gdgAcD3shKKiG4F8BkAdwohLmjHZ4jIlp/npVw/zlCuqLnLdbwkNwN4QQjxmjqQ5XhF6Qdkucay2CVOcbf5drg7zK8AOJmjHDfCfY16DsAp+ed2AF8F8AN5/DEAcxnLNQ93d/1ZAD9UYwRgG4BvA3gJwJMApnIYswaAcwDGtWO5jBfch8sZAOtw/Zn3xo0RgJNyzZ0GcFvGcr0M1/+q1tmD8tyPyDk+BeAZAL+RsVyRc5fneMnj/wzgDwPnZjleUfohszXGmacMwzB9RplcMQzDMEwCWLEzDMP0GazYGYZh+gxW7AzDMH0GK3aGYZg+gxU7wzBMn8GKnWEYps9gxc4wDNNn/D/zRIZIz3PGEAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plot(random_sd[0:200]);" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Periodic vibrations\n", "\n", "A vibration that has a periodic vibration will sound like a tone to us. \n", "\n", "The sine function has a periodic cycle, and can be used to represent such a tone. \n" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "slideshow": { "slide_type": "skip" } }, "outputs": [], "source": [ "Fs, Len, f1 = 44100, 3, 440 ## sample rate, length in second, frequency\n", "t = linspace(0,Len,Fs*Len) ## time variable\n" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "slideshow": { "slide_type": "slide" } }, "outputs": [ { "data": { "text/html": [ "\n", " \n", " " ], "text/plain": [ "" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "sine_sd = sin(2*pi*440*t) \n", "\n", "Audio(data=sine_sd, rate=Fs)" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "We can plot this sine sound, using a plotting command." ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "slideshow": { "slide_type": "fragment" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAD8CAYAAABzTgP2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJztvXt0I9d54Pn7AJAAHwDfBNlkP9XdklpqqWV3ZFmWZTu2HMmTRPass2tnJrbP2tFqY2eT2d2ZKCc7WU8yPpPNniS72Tj22mNvlMSxxxvba40t22Mpfsnyq2VJ3a1HP9QPkWwS4BMAHwAI4O4fVQWi2WTzgXrcAut3Dg6BQlXhXlbV97rfd68opQgICAgICLAIed2AgICAgAC9CBRDQEBAQMBVBIohICAgIOAqAsUQEBAQEHAVgWIICAgICLiKQDEEBAQEBFxFoBgCAgICAq4iUAwBAQEBAVcRKIaAgICAgKuIeN2A7dDb26v27dvndTMCAgICfMUzzzwzpZTq22g/XyqGffv2ceLECa+bERAQEOArROTyZvYLQkkBAQEBAVcRKIaAgICAgKsIFENAQEBAwFUEiiEgICAg4CoCxRAQEBAQcBW2KAYR+ayIpEXk9Drfi4j8pYicF5GTIvKamu/uF5Ez5neP2NGegICAgIDtY5fH8DfA/df5/gHgkPl6CPgEgIiEgY+b3x8B3isiR2xqU0BAQEDANrBFMSilvg/MXGeXB4G/VQY/BjpFZBC4EzivlLqglCoCXzD31Y7TYxk+/f0LnE/Pe90U17g8vcCnv3+B50fmvG6Ka0zmCnzmqYs8dW7K66a4Ri6/zN/96BLfODVOpbIzlvotlMp88cQI//jMKMVSxevmaIdbBW5DwEjN51Fz21rbX7fWCUTkIQxvgz179jjTynX4zstpfvNvT1CqKP7822f5h998HXfs6XK1DW5zeizDf/N//4iFYpmQwF//i9dy/60DXjfLUSYyeX7lr55iMlcA4N/+8hE+eM9+j1vlLAuFEu/+xI84k8oB8IG79/HRX73F41Y5S7mi+NCjJ/iBqfy/cWqcT7/vOKGQeNwyffDN4LNS6lNKqeNKqeN9fRtWdNvGfKHEv/nSSQ72t/ON33kj3W3N/Ot/PEmp3LhWRqWi+L0vnaQ9FuHr/8M9HB3q4JEvnySbX/a6aY7yvz52moVCiX98+PW8/UiS//D4S1yeXvC6WY7yF98+y9l0jk+/7zgfuHsff/P0JX58YdrrZjnK535ymR+cm+KPHryFP3jHzTz5cpp/fGbU62ZphVuKYQzYXfN52Ny23nZt+E8/G2EyV+Bj7zrKzYMJ/pd/djPn0/N84/SE101zjCdfTvPClSyPPHATt+zq4GPvOsrc4jJ/96NNVdP7kpfGs3zrhRQPv+kGju/r5t+/81YiYeGvv/OK101zjOn5An//k8v88zuGue9IkkceuInBjhh/+eQ5r5vmGMvlCp/47ivcub+b37hrLx96435uH+7g4989T3mHhNE2g1uK4THgfWZ20l1ARik1DvwMOCQi+0WkGXiPua8WKKX43I8v89q9Xbx2rxE6evstA+ztaeUffvKqx61zjs/95DIDiRi/ctsuAG4d6uANB3v4/E9fbdgY9Od/+irNkRDve/1eAPoTMd55bIjHnr/SsJ7SV54dI79c4b970wEAYk1h/uVde3n6lWkuTDbmWNo/vZxmPJPnoTceQEQQEX7z3gNcnl7kB+cmvW6eNtiVrvp54EfAjSIyKiIfFJGHReRhc5fHgQvAeeDTwG8BKKVKwEeAbwEvAV9USr1gR5vs4PRYlgtTC7z7tcPVbeGQ8K47hvjxxWnGM0sets4Zsvllfnh+il89totIeOX2ePdrhxmdXeLnr8562DpnqFQUj58a574jSTpbm6vbf+34bpaWyzz5UsrD1jnHV5+7wu3DHRxOxqvb3v3aYUTgayfHPWyZc3z1uTF625t5840r4ej7jiSJxyI89vwVD1umF3ZlJb1XKTWolGpSSg0rpT6jlPqkUuqT5vdKKfVhpdQNSqmjSqkTNcc+rpQ6bH73MTvaYxdPvJRCBO6/5epB13ccHUQp+N6ZxrMwvn92kuWy4r4jyau2/+KNSUIC3zvbeH0+NZZhar7I227uv2r7Hbs7SSaifLMBw4aTuQKnxjK8fdW9nUzEOL63qyFDpcvlCt8/O8V9RwauMnqikTD3HUnyTy+nG9Yj3iq+GXz2gh+en+K2oQ662pqv2n6ov51kIlrNamgknnwpTVdrE69ZlXXV0drEHXu6GlIxfOdMGhF40+GrFUMoJPziTUmePj/dcMkGVtjk3kPXJnK86XAfL41nmVkout0sR3l+ZI75Qol7D/Ve8929h/qYW1zmxfGsBy3Tj0AxrEMuv8yzI3O84eC1N5GIcM/BPn74ylRDDVhVKorvnEnzlpv6Ca+Ruvemw32cGsswPV/woHXO8Z0zk9yxu5PuVQYAwF0HuskVSg0nML5/dpKetmZu2ZW45ru7zXu+0bKTfnBuipDA3Tdc+0zffbAHgKfON56xtx0CxbAOP704Q7miuGcNxQBw7+Fe5haXeeFKxuWWOce59Dxzi8u8YY0HB+DNN/ahVGM9PJmlZU6Ozl3jLVjcdcAQGD+5cL36TX9RqSh+cG6Kew/3rZm7f9tQB+3RCE+/0jjXGYz79uhwJx2tTdd81x+PcWMyzg8b6N6uh0AxrMMPz08TjYR4zd61C9ksq+PpVxrHqrIGll+7Tp9v2WUIjGcuN84A9MnROZRav8/JRIwDvW0NZT1fml5geqHIXQe61/w+Eg5x5/7uhrq3F4slnhuZ4x7TM1iLuw/28NOLMxRKZRdbpieBYliH50fnuG24g1hTeM3v++JRdne3cHK0caaLeObyLN1tzeztaV3z+3BIODrU0VBTZDz76hwicNvujnX3ed2BnqoH2Qg8b96zx3avX71/14FuLkwuNEzY8IUrWcoVdc3YWS3H93ZTKFU4O9GYqbpbIVAMa1AqV3jhSoajQ53X3e+24U6eH2mcUNLPX53lNXs6EVl/aoDbd3fy4ni2YayqZ1+d5VB/O4nYteEFi9fu7SJXKHFxqjEExvMjGVqbwxzsb193H+vePzXWGPe3ZcwcHV7fALjN/O75BjL2tkugGNbg/OQ8+eUKt1/HigQjFjs2t9QQVtXcYpELkwvrhs4sju3uYLmseGk851LLnEMpxXMjcxzbfX0D4NYhY4D2hSuNMQD93MgcR4c61kwwsLjF7PPpBlEMJ0czDHbE6I/H1t1nuKuFrtYmTo02Rp/rIVAMa3DSvDGODm2gGIYNgXKyAR4eq88bCcnbze+fa4BCt8vTi8wuLm84IeINfe00R0INISSLpQovXslueJ0TsSb297Y1jMdw0gwNXw8R4ehwZ0M8z/USKIY1ODk6RzwaYV9P23X3OzrcgQicbIBw0ssThjV888C16Yu1DCRi9MejVUXiZywPYCMDoCkc4uaBOKfH/O8xvDyRpViuVI2a63HLrkRD9DmzuMyl6cVN9fm2oQ7OpnIsFRsjVLpdAsWwBqfGstwylNhwGt72aIQDvW2cboCU1ZfHcwwkYtcU861GRLh5MMHLE/4PJZ2ZyBISrhtrtziyq4MXrmRQyt8D0NZ1W6t+YTVHzVCp3wvdrJTyjQwAMIy9ckU1XN3KVgkUwyoqFcXZiRw3D2784ADcOBDnXMr/QvKliRw3DcY33hGjz+cn531fDfzyRI59vW3rZp7VcutQgmy+xOisv+fHOjuRIxoJsbt77cyzWm41Banfa3WstSZuGtj4/rY85rMN8EzXQ6AYVjE2t8TScvmqicWux6H+OJdnFn3tei6XK5xP57hpgzCSxeFknGKpwuWZRYdb5ixnU7lNCQuAI6ah4HdL8kwqx6Fk+3UHni2sZ8DvqxaeTc3T0dJEXzy64b7DXS20NIU5l/J3n+slUAyrsCyFw8mNwwtgWM9KwSs+nqb4wuQCy2W1aSF5oykwzvo4nLRYLHF5ZpEbk5tThjeY4Sb/C8ncpo2e3vZmOlubOOfzPp9L5TicbL9uGrZFKCQc7G/nXNq/97YdBIphFWdNS+Fg/+YeHkuB+Nn1tAaeNxtKOtjfjgi+Hmc4l5pHKbhxYHMGQCLWxEAixis+FpJzi0VS2UJVsW+EiHCov53zPraelVKcTeU4tMk+gzFJpp+fZzsIFMMqzqVyJBNROlrWL3iqZW9PG01hqSoUP3IuNU84JBzo3ZyQbGkOs7e71dcPzxlTqd24yfAZGArxvI89Q+sePbxJzxAMA+lsOufbQfd0rkA2X+LwJhIMLA4l46SyBTJLjblA02awa6Ge+0XkjIicF5FH1vj+X4vIc+brtIiURaTb/O6SiJwyvztx7dnd5Vx6ftOuNhipjAd62309AH1hap7dXS00RzZ/OxxOxquDen7klcl5miMh9mxiENbiYH8759Pzvp2z/0w1TLr5+/twsp25xWWmfZqZdHabfQZ8/UzXS92KQUTCwMeBB4AjwHtF5EjtPkqp/10pdUwpdQz4feB7Sqna6SrfYn5/vN721EOlojifnufQJsNIFoeS7Zz1cUzywuQCB/o2b1GBEXN/dXrRt5lJF6YW2NfTuqlBWIuD/e0sFsuMZ/MOtsw5LkzO09ocZlfH+tW/q7GeBb8Oxlrt3kooyVIifo4C1IsdHsOdwHml1AWlVBH4AvDgdfZ/L/B5G37Xdq5kjIykG/qvX9i2mgN97YzNLvly/qBKRXFxaoEDvVvr8/7eNkoV5dv0zYtTC+zfYp8P+nwA+tLUAvt62jY1CGtxyLKefWr4XJiaJx6L0Nt+/fqcWoY6W4hGQg277vVmsEMxDAEjNZ9HzW3XICKtwP3Al2o2K+AJEXlGRB6yoT3b5vK0kX65f4OK59Xs722lomDEh+mbVzJLFEqVLXsMliK5OLXgRLMcpVSucHl6gf2bHFOx8LtiuDi1wP6+rd3b/fEobc1hX15ngEtTi+zv3ZoyDIWEvT2tXJr23/NsF24PPv8K8MNVYaR7zBDTA8CHReTetQ4UkYdE5ISInJicdGZ5Sevm37dFS9KaOuPilP9upAuTRp+3aj3v87FiGJtbYrmstuwl9bQ1k4hFuOTDPi+XK4zMLm3Z6BER9va0+bLPYNyfG01tsxb7etq4NO3PPtuBHYphDNhd83nY3LYW72FVGEkpNWb+TQNfwQhNXYNS6lNKqeNKqeN9fdeuU2sHl6YWiEZCDCQ2H4OFFaHqx4fHcpdv2KIl2dPWTDwW8aViuGC2+cAW+2wJST8W9o3MLFKuqC0bAAD7elur3rSfyC+XuZJZ2rKhB8Yz/er0YsOswbFV7FAMPwMOich+EWnGEP6Prd5JRDqANwFfrdnWJiJx6z3wduC0DW3aFpemF9nX07bhHEmr6Ww1CoEu+tDCuDC1QHs0sqmq0FpEhP29bb5UDBe36SUB7O1p5bIPr7Nl/W5HSO7taWNk1n+JBiMziyhlhHq3yr7eNorlClfm/DmGVi91KwalVAn4CPAt4CXgi0qpF0TkYRF5uGbXdwH/RSlV+1QlgadE5Hngp8DXlVLfrLdN2+XS9MK6q5dtxD6futuXphfZ19u6pRishW8Vw9QCiViE7g0mDFyLfT1tjM4usewzIWmFDLcaPgNjzG25rBjP+Csbqxoa3kYoyZIDOzWcZMsYg1LqcaXUYaXUDUqpj5nbPqmU+mTNPn+jlHrPquMuKKVuN1+3WMd6QbmieHV6cVtWJBhC0o+KYXRmkd1d21OG+3vbuJJZIr/sr2ysS9ML7NvigKTF3p5WyhXlO0vy4tQCna1NG86euxaWkPSbEWAJ9e0809XwsA9DaHYQVD6bjGeWKJYr23K1wbBKrmTyvhKSFTPddDMzba7F3p5WlMJ3Kaujs0vbVoZ7e/wpMF6dWWTvNq+z9Uz4LYR2cWqRztYmOlu3rgyT8RixppAvjT07CBSDiTW4tt1QknXc6Kx/BEY6V6BYrmxbMVjCdcRHfa5UFGOzSwx3t2zr+H3mdfabkBydXWJ4m9e5Px41hKTPlOFIHcowFBL29bT57jrbRaAYTCyBvpUpEmoZ7jIEzYiPrGdLoO/u2p6QHO6ylKF/+lxVhtv0GPriUVqawr7K0qkqw21eZxF/CsnR2cVtK0Mw7m8/3dt2EigGk9HZJcIh2XKqqoUfhaRVkLddj6E/HqU5HGLUR+mbVWW4zT4bKav+ykyqVxmCYfj46d6uVBRjc9tXhrDSZ79OIFgPgWIwGZ1dYrAjRiS8vX9JVUj6KKwyMmM86EOd23t4QiFhyGcCw1KGdggMv2Ddk/X1uZUxHwnJdK7AclnVrQznC6UdOctqoBhMRmcXty0gwZ9C8tWZRZKJ6KaWtlwPQ0juHGUIppD0UVZSvV4SGP+vXKFEdqlkV7McxR5l2GKeyz/X2i4CxWAyOrtUDQdtF79ZkiOzi9seU7EY7mr13bhKvcpwqLOFXN4/luSoLcrQGkPzhxFgPYf1PNN+DA/bRaAYgGKpwkQ2X5d1AcbDM+aTBwfqq2Gw2N3dwsxCkYWCfyzJevs8ZN4nYz4RGCOzi2ZmUT2eofE/84unZK/H4J9n2i4CxYBRw6BUfTcRGA/P1HyRpaL+tQzFUoXxbL6urA3wn1U1MrP9ug0Ly/L2j5CsbxAWVpShX67z6OwSve31KcOOlibaoxHf9NlOAsWAPW6ncbwlMPS3MK7MGcpwu6mqFn6yqpbLFcYz9QtJP/UZ7AmTdrU20doc9o2XZIcyFBHfhYftIlAM2ON21h5vDXDqjB0DkrBS5OaHh2d8Lk9FUXcoqbutmVhTyBdC0pq+Y/c2C/osRIShTv8kGozOLtb9PIP/kivsIlAMGLHicEgY3MKSh2uxElbR/0aylFe9iqG33RCSflikyFKG2616trCEpB9CSRPZPKWKqttjAP8kV5SrNQx29Nlfabp2ESgGDGt3ILH9GgaLvnarlkH/h2d0dpFIHQV9Foa77Y8K0WpBnw0CY8gnKav29tkfyjCdy7NcVrZ5DH5K07WLQDFgTzwS/FXLMJHJk0zECG9x7Ym1GO5qYdQH4ypjc0uEhLo9Q/CP9bwyfmaHkGwls7RMLq93mq69ffZXmq5dBIoBKx5Zv0UF/olJXsks2SIgweizH8ZVrswZyrBezxCMzKSZhSKLRb0tydHZRURgsNMeZQj6Z2OtjBnaE0oC/ftsNzteMVg1DEM2WBfgH0tyPJNnsI6Cp1oGO1rILC1rLyTHbVaGgPbrMozP5elrjxKNbD9t08JK0x3V3Aiwo6DPYqdWP9uiGETkfhE5IyLnReSRNb5/s4hkROQ58/WHmz3WaVJZI1Nl2EYhOb1Q1HpdBqWM1bjsEpK7TGtU9xW+7FSGVSGpucCw1zP0R3LFlUye7rZmWprrV4YdLf5K07WLuhWDiISBjwMPAEeA94rIkTV2/YFS6pj5+qMtHusYljCzw9WGlfh1KquvkJxZKFIsVWwTGIMdhpCc0FgxKGWkbe6yXUjqLTAmMvnq9amX3vZmmiMh7cMqEzYqQxEjW3Eiq3ef7cYOj+FO4Ly5TGcR+ALwoAvH2sJ4xrjgdgtJna3nqjK0SWBY/zudwyqzi8sUShXb+twfj9IUFu2F5Hgmz4DtQrJgy/mcwk5vGIznROfn2QnsUAxDwEjN51Fz22ruFpGTIvINEblli8c6hnXBB+wSktWwir4CY0Ux2PPwJBP6h5IspWVXn0MhoT8eI6Vxn7P5ZeYLpWqozw4GO2JMaHxvg6UY7HmeAQY6YozP6XudncCtweefA3uUUrcB/xfw/231BCLykIicEJETk5OTtjVsIpMnHovQHo3Ycj5L8OgsJKtekk0CI9YUpqetWes+T1RDhvYJjMGOmD/6bKOQHOxo4YrGQnKxaMx6a5eXBLCrI0Y6l6dUrth2Tt2xQzGMAbtrPg+b26oopbJKqXnz/eNAk4j0bubYmnN8Sil1XCl1vK+vz4ZmG9iZqQLQ2hyho6VJawtjPJOnKSz0tkVtO+dgZ0xzL8lom11jDADJjpjWY0l2e0lgWM+pbJ5KRc9KYLu9YTCiCRUFk/N6h9DsxA7F8DPgkIjsF5Fm4D3AY7U7iMiAiIj5/k7zd6c3c6zTGDFY+ywq0N+SHJ9bIpmIEbKhuM1isKNFa2V4xVKG7TYqw4RxnXWdLmHcIS+pVFFMLegpJJ3xkqwxNH3vb7upWzEopUrAR4BvAS8BX1RKvSAiD4vIw+Zu7wZOi8jzwF8C71EGax5bb5u2wngmb6sVCWifxWD34BxYylDjPjugDAc6Yiwtl8nm9azfGM/kETEGyu3CmkJF1ww0JzwGK+Sqa5+dwJbAuhkeenzVtk/WvP8r4K82e6xbFEsVpuYLtsYjwXA9T45mbD2nnYxn8hzb3WnrOQc7WsjmSywUSrTZNF5jJ1cyeXbZ7Bla981EJk9HS5Ot57aD8bklM3vKvqHE2qy724ZtO61tjJvhMzuf6cGE1Wd9DR+72dGVz6lsHqXstS7AiGPrWuRWqSgjt93GTBWoHXTX8+EZzyzZbwAk9O7zRNaBMKnm1vN4Nm9Oi15/cZtFoiVCS1NY6/Cw3exoxTCRtT8eCSvWSlrDfO+ZxSLFcoXBOmdVXY3O2ViViiKVKdiuDK3rrOsAtJ0FfRbdrc00h0NaXmcwFFa9MwavRkQY7IxpqwydYEcrBifikcb5zHl0NLQkrQFiOwckAXaZ59NxAHp6wVCGdoeS+uP6KsOVaU/s7XMoJCQ7otp6SVfmlmyt27AY7Ihp+Tw7xc5WDA7EI0Fvd9vuSm+L/kTUPP/O6XNzJERve1RLjyGbL7FYLNveZzBi7jpeZ7DCZ/b3eSDRouXz7BQ7WzFk8rRHI8Rj9g4cVtPbNLQw7J4OwyIaCdPbrqclaaUZ2t1n45x6pibbXcRYy0CHnmGVpWKZucVlR67zrs4Y6VxhxxS57WjFMOFA2iasFLnp+PBcySzRFBZ62pptP/euzhhXNOxzOme0KdlhX9qmRTKhp5Cshgyd8BhMxaBb/YZTniEYyrBcUTumyG1HKwYnMlUsdLUkJ8xJ1ezM57dIJmKkNQyrpLMFwiGhx8ZKbwujZkW/PlvhraTNA7Fg9LlYrjCzULT93PVgKWgnnmmdkyucYIcrBmc8BtDX3R6fs39A0iKZ0DPensoai9XYsYzpagY6YswtLmuXmpwyM+KsAXI7GdB0BmFLQdudlWScU9/kCifYsYqhVDaL2xy4iQCS8Vg1hKET41l754aqJRmPmdNbayYkcwWSCfu9BdC3EjiVy9PTZqyfYDe6Ws+WMnTCS6qmoGv4TDvBjlUM0wtFKgr6nVIMiSiTuQJljSYbU8rI53dMGSb0rN9IZ/P0OWA5w4rA0E1IGn12SBlqWr+RyuaJRyOOVN53tTbRFJaq8ml0dqxicDIGC4bCqSiY1miwam5xmWK54mCfDUGkm1WVdtJj0FRIGn125jr3tDUTErQbT0rn8vQ5dJ1FjPU3dOuzU+xgxWC5nc7cSNZDqZOFkTIFdv8O6nOhVGZmoeiYkBzQdJGiVDbv2L0dCVv1G/pcZzA81aRDniGYY2iaGT1OsYMVgykkHbqRrIdSJ0sy7WAMtva8OsXbJ3POGgBt0QjxWESrVc3KFcWkgx4DGNdaNyGZyjmnDMHss2bK0Cl2rGJI5wqIGAucO0HVetbo4amGzxxShl2tTTSHQ5r12bnsHAvdUlan5wuOjp+BMZW3TkJSKUUq64Iy1Og6O8nOVQzZPL3tUSI2TklcixWH1enhSZvWs1OhJBGhPxHVavB50uHwGehX5LaiDJ3rc38iVv3f6kBmaZliqeLYgDsY91AuX2KxqOf6G3ayYxWDkzFYWInD6jRYlc7mScQitk5JvBrdrConUxgtdPMYqpXejlrPUabmiyxrMkVEOuf8dbY8bZ0MH6ewRTGIyP0ickZEzovII2t8/y9E5KSInBKRp0Xk9prvLpnbnxORE3a0ZzOkHB6oAj2FpJMPDuhX5JbK5omEhO5WZ0KGYAxAT2o0j47TiRXGuY37yBrD8Rqnswxrz63T/e0UdSsGEQkDHwceAI4A7xWRI6t2uwi8SSl1FPhj4FOrvn+LUuqYUup4ve3ZLOlc3tHwAlhCUo8HB4zxDqf7bKT0adTnbIG+eNSRKUAsrNTkqXk9pohIZfPm+JmTikGv5Ap3lKHZZ02UoZPY4THcCZxXSl1QShWBLwAP1u6glHpaKTVrfvwx4OmigMvlCtMLRUcHJMEQGDrl9DudzgeGVZUrGEt86oBhADjb5wHNLMl0Lk9Pm71Leq7GenZ0MXyczjKElcF8ncLDTmHHnTMEjNR8HjW3rccHgW/UfFbAEyLyjIg8tN5BIvKQiJwQkROTk5N1NXhqvoBSzrqdYMQkdYnDKqVcEZLJapGbHgLDUIZOe4Z6KYZUtuDowDPoV8w4mSsQj0VoaXZu/MwYnwtpc52dxNXBZxF5C4Zi+L2azfcopY5hhKI+LCL3rnWsUupTSqnjSqnjfX19dbXDDbez9vw6xGFnF5dZLivHBYZ2QjKXd2VcBfTpc9rhfH6AnjZjUkJdwoZGMomz19mqftbFS3ISOxTDGLC75vOwue0qROQ24D8CDyqlpq3tSqkx828a+ApGaMpR3BioghWrSgeB4UaminF+ffqcXzYWbnFaGfaYM7fqIjDcSDIIh4Q+jVavczrL0CKZiGrjJTmJHYrhZ8AhEdkvIs3Ae4DHancQkT3Al4HfUEqdrdneJiJx6z3wduC0DW26LulqPNL5gVjQIw7rlpe0Eof1vs+TLqQwgl5C0po12OmQIVhTRHh/ncEKnznf5/6EXskVTlH3NIRKqZKIfAT4FhAGPquUekFEHja//yTwh0AP8NciAlAyM5CSwFfMbRHgH5RS36y3TRuRzhUIiWHpOUl1tlENLAw3BucA4tEILU1hLYRk2oXiNotkIqpFLcPUfNEcP3O+z/2JGCMzi47/zkYoZUwB4sp1jsf4bjbt+O94jS3z0yqlHgceX7XtkzXvPwR8aI3jLgC3r97uNClzSmInFm6ppaet2QwxeC8wJh2uerYQEW0sSTeK2yySiRiXp70Xkm4ZAMbMd2KsAAAgAElEQVRvRDlxacbx39mI6qzBLvQ5mYiyUCwzXyjR7sD03rqwIyuf3YjBAoRCos2cMqlsno6WJkerni36NSnsS7kUMgR9JpVLOzxpYC3JhB4LM6VcGj+r/Q0d7m8n2aGKIe+KRQV6CUk3BCTos/ZzOlegKSx0OVj1bJFMRLVY4tOtxArjN/TIuqvODeVK+Eyf5Aon2ZGKIe1SPBIgGddjUjknF25ZTdL0kpTydvU6ywBwsurZQpfV69LZvDF+1ua8MuxP6JFcYRkh7oSS9LjOTrPjFEOxVDEWbnHJY9AmxJB1URkmYiwtl8nmva1+drvP4P0066lswdFZg2tZmVTO2z47PWtwLUEoqUGZnHcvBmv9jtchhmrVs1vKsEMXgeFu+Ay8X6TIjYI+C11qVtwcP2uPRmhrDnvuJTnNjlMMbsZgYcXd9jIOa1U9u6YM45bA8D727NZ11mW+JDemw7Doam0mEhLPM9DcHD8DfaIATrLjFEO1uM3tEIOHAsNtZahDn/PLZTJLy671OdESIRoJeT5H1KQL82FZWFl3Xsfb3Rw/A8zFqALF0FC4sdRjLSvutncPj5tpm1CTueGhVWUJKydX9KrFqN/wdiW35XKFqfmia54h6DGDsJtjSbAz1n7ecYohncsTDokrWRuwMkDnpfXsxupWtbQ2R4jHIp5akm7NDVXLgMepyW5NAVKL1wszWeNn7vbZuM5eZ905yY5TDFYM1o0URoDO1iaawyGPrWfjt92ynsH71evcmhuqlv5E1NNQ0krIcOdYz27NGlxLfzxKoVQhu6THmiNOsAMVg3sxWDCn6k14G4dNZQuuZW1YeD13UMrF3HYLK5TklSXpdpgUjD5nlrzLunN7/Kz2txp5AHrHKQY3Fm5ZjdfWsxvz868m6fESn6lcnuZwiM7WJtd+c8Cs38h5tHrdpIuTBlpYXqhX19orL6n2txuRHacY3Fj3eDXGfEnehlXctKhgZVDSK+t50lzr2Zy51xWqq5p5dK1T2YI5fua+kPRqADrtiZfkfUKJ0+woxVAoGQu3uBleAHPuIA9jz2lzNlk36Y9HWS4rZheXXf1di5QXXlK1yM0767mv3flZg2vxWkhaCsnN+9tSQl5nYznJjlIM6az7WRtgWJK5fInFovshhkpFMTnvvsfgtSXphZfkdZFbKlfwJGQIeDaelMoW6Gx1d/yspTnsedad09iiGETkfhE5IyLnReSRNb4XEflL8/uTIvKazR5rJ24u3FLLypwy7t9Is4tFo+rZ9XEVby1Jt6thwfv6DcMzdFcZWll33hkAedcjAOD9uKHT1K0YRCQMfBx4ADgCvFdEjqza7QHgkPl6CPjEFo61DTcXbqmlGnv2IJy0MiWxNx6DFw/PUrFMLl9yvc9W/UbKoyK3tAceg4jQ52H1c8rFmZJr6Y97m5rsNHZ4DHcC55VSF5RSReALwIOr9nkQ+Ftl8GOgU0QGN3msbbhdAWzhpZBcKfRyt88r2Spe9tl9S3LAo7z+QqlszBrsQZ+TiahnHsOki2ur1BJ4DBszBIzUfB41t21mn80caxtuLtxSi5fVz15kbQDEmsJ0tjZ5IiS9KG6z8GqCtUkXV25bjVdFbpWK8sRLAqq1SY1a/eybwWcReUhETojIicnJyW2do7UpzO3Dna5VPVskWiI0R0KezLCa8qDq2cJwt90Xkm6ue7ya/kTUk1CSVyFD8M56nlksUqooT7yk/niMYrnCnEdZd05jh2IYA3bXfB42t21mn80cC4BS6lNKqeNKqeN9fX3bauhvv/UQ//jf372tY+vBmGDNm1qGdM79rA0LryxJN9c9Xs2AmZpcqbhrSbq5itlqvMq686K4zSLp4bihG9ihGH4GHBKR/SLSDLwHeGzVPo8B7zOzk+4CMkqp8U0e2xAk494ISa+yNsCwqjwZY8jmaY6E6Ghxr+rZIpmIUaooZhaLrv6um6uYrabfo6w7q89uZ2JB41c/160YlFIl4CPAt4CXgC8qpV4QkYdF5GFzt8eBC8B54NPAb13v2HrbpCP9iagnsWevsjbAGpR033pOZY3iNjerni0sS9Lt6bdT2TyRkNDt8vgZeGc9p730GDSYNdlJInacRCn1OIbwr932yZr3CvjwZo9tRPrjMX5wdsr1353M5jnY1+v678LV1nNvu3sPrzGDrjde0tWFfR2u/W7KnALE7fEz8M56tjxwT8bPglBSgB0kEzFyhRILLk6w5mXWBqykBbstMLyYNNBiRUi6HVZxd02CWryynlPZPF2tTUQj7o+fxZrCJGKRhl3JLVAMLmEJSTctDCtrw+26DYv+qvXsdojBO4/BmLjPm1CSV8rQyrpz/Tq7vKTnarxei8JJAsXgEtUQg4sWRtrDFEaoiT272OeFQolcoeSZwGgKh+hpcz9NN53zThlaWXduW89pl9dWWY1X44ZuECgGl6jOHeSiVZXysAIYVmK/blpV1ewcj7wksJa7dK/P+WVz1mCPPAbwJusu5cHaKrV4veaIkwSKwSVWUvrc9Bi8y9oAiEbCdLc1uxp7XumztyEGN0NJkzlvPUNwv+K77NGswbV4veaIkwSKwSUSLRGiLsdhvczasDAWKXLTS/KuuM0iaQoMt/BiecvVuL187fRCgXJFeXydvV1zxEkCxeASRhzW3akDUtk83W3NnmRtWPS7LCQtj8Fb6znK1HyR5XLFld/zcm4oi/54jHkXs+68Hj+DlShAI9YyBIrBRdyeFsPI5/dOWAAkXZ6SOZXNE42ESMRsKdHZFkmXs7FSHk6HYeF2kZsOXpLbfU7n8vz5t89yYXLe8d8KFIOL9MfdXeLTy9x2i2QixuS84fa7gZXC6EXVs8XKIkXuGAHpXIHmcIjOVvenALFwu8hNBy/J7T5fmFzgL588x7gL41eBYnARt+OwXua2WyQTUcoVxfSCe5ak9312N9HAWtN7JynDVDaPCK5W1K/G7TVH3PSSAsXgIsmEe3HYckUx6XEBENQUubmkEL0sbrOw/uduZSalPKz0tuhzeSK9dC5PT1uUprB3IizWFKajpcmD8Jnz1zpQDC7iZvXz9HyBivJ2cA68CDF4Hz7rbm2mKSyu1aykst4bAIlYhFiTe2s/pzUYPwN3xw1T2QKtzWHao86PnwWKwUXcFJLVGKzHD4+bynC+UGKhWPbceg6FhP64exloOijDlaw7l5ShBl4SuDsthnWd3QgZBorBRdyMw+qQtQG11c/O99kK3Qx0eNtnMKdLcKHPS8UyuXzJ01oVi6SrytB7LwmM+9u9sST3vKRAMbiIm3FYr6fDsGgKh+htb3bFqqrWMHg8xgDuTRGR1uQ6A/SZ6284TalcYWq+4HmYFFay7txYcyTlYpZhoBhcxM04bCpbMLM23F+4ZTVureS2ogy9t54HOmKurP2sQ9qmRdKl6zw1X0QpXfpsVT87u2KfUopUNu+aN1yXYhCRbhH5toicM/92rbHPbhH5joi8KCIviMjv1Hz3UREZE5HnzNc76mmP7rgZh01n8/S2R4l4mLVhkXRpFsqJjCUkvbck+xNRV9bf0CVkaLQhykKxzLxbfdbAM+yvjhs6+0xnl0rklyu+CSU9AjyplDoEPGl+Xk0J+J+UUkeAu4APi8iRmu//Qil1zHztgJXc3Ik965DPb9HvUlgllc0Tj0ZocyFrYyMsoeV0aEUnIelWcoVuyhBwPArgdmi4XsXwIPCo+f5R4J2rd1BKjSulfm6+z2Gs7TxU5+/6lv5ErDobppMYUxJ7/+CA8fBMzxcoOTx3UDqX92x969VYLr/TtQzpXMGYAqTFe2XY71JyhQ4TJVr0uzRu6LYyrFcxJJVS4+b7CSB5vZ1FZB9wB/CTms2/LSInReSza4WiGg23MjcMIamHYuhPxKgomF5wNg6rS6YKuGhJZg1l6GXVs4VbQjKdzRMS6PGw6tnCNWXo8ljShopBRJ4QkdNrvB6s3U8Zk5KvOzQvIu3Al4DfVUplzc2fAA4Ax4Bx4M+uc/xDInJCRE5MTk5u3DNN6XchDrtcrjA1X9TCogL3QgwTmTwDmiiGfpf6nNbMMwR3lGFfPEo45L0yjEbCdLY2OT6G5rbHsKH/qZR623rfiUhKRAaVUuMiMgik19mvCUMpfE4p9eWac6dq9vk08LXrtONTwKcAjh8/7tuVMWqXu2zva3fkNyZz+gzCQm39hnOWpFJKKy8pHo3Q2hyuDog7RSqX5+aBhKO/sVnazT47PZ6U0mDak1rcWMktnc3T0dJErMmdKfTrDSU9BrzffP9+4KurdxDDx/0M8JJS6s9XfTdY8/FdwOk626M9ybjzWQxuzqmyGdzwGGYXl1kue7twSy3VDDSHLcl0tqDNuIpba44YM+jq0Wew1n52Xhm62ed6FcOfAPeJyDngbeZnRGSXiFgZRm8AfgP4xTXSUv9URE6JyEngLcC/qrM92lOdVM5BgWEpHV2sqp62ZkScnYXSEka6hJLAyEBzss/zhRLzhZI2niFYlcDOW8+6eIbgTp3OhMvTntSVyqCUmgbeusb2K8A7zPdPAWsGA5VSv1HP7/uR/mooybmHR6dqWIBIOERvu7NVsRMarNy2moGOGD9/ddax869UeutjPScTMU6Nzjl2/mKpwvRCUZtxFTA888mcUf0ccmjcI53Nc0NfryPnXgvvq592GPFohJamsKPu9kQmTzgk9LR5X/Vs4fQslGnNwmewMsGaU4vFpzUbSwKjEtjJPk/O65OqapFMxChVFDMOVT9XKsr18FmgGFzGiMM6G5O0lvR0ynrZDk7PHaRb+AwMS75YqpBZcmaxeN3GksAQkkvLZXIOZd3pVNxm4fTkmDOLRUoV5WqfA8XgAU5PyaxTdo5FfyLm8LhKnp62Zpoj+tzS1SI3h661FY7U6Vr312TdOUE1fKaRMnR6ckwvlKE+T9EOot+MSTpFKptnQKMHBwzreWq+yLJD1c8pzQYkoTYbyzmB0dIUJq7BFCAW/Q5n3a0UeulzrZ2u30h7MFFioBg8wErpcyoOq1MFsIXVHqcUotvpfJthJTXZGYGRyhW0qXq2cFpIprJ5IiGhu1Wf8bOVNUcCjyGgDvrjURYdqn7OL5fJLC1rqBicXcnN8JL06nN1ugSH5kuayCxp2GfnPYY+zcbPopEwXa1NjhkAVijSzcWYAsXgAVbs2YkbqRp31iiFEZwtctNp4ZZaYk3OTpcwnskzqMFqdbW0RyO0RyPOeUkurkmwFZws7EtlC/S2N9Pk4hT6gWLwgMGOFgCuzNl/I1lCSDch6eSg5NR8kYomC7esxqlsrErFWrilxfZz14uTU8tfySxppwwBBjtijDvkGaayedez7QLF4AGDDk7JbN2cuj08PW3GpGdOCEkdq54tkh3OWJLTC0WWy0q76wyGR+zEva2UMidK1E8ZDnS0ODbF+ngmz67OQDE0PJb17ISFMT63BOinGMIhMdd+dk4Z6jauAlbBl/19toSQjmGVwY4WR+7tbL7EYrGs3b0NxvM2vVAkv1y2/dzjmaVqlMEtAsXgAdFImN72KBPZJdvPPZ7J0x6NEI812X7uehlIxBzJ6R/P6KkMwVwsPlegbPNi8Tr3ebAjRjpn/8JMljIcdNl63gzWdbC7lmGpWGZucdl1AyBQDB7hVExyXNMYLDhnSY5n8kQjIbo1mgLEItlhLFI0NW+vwLAUrJYeQ2eMckVVp6+wC72VoTlumLHX2LP6HISSdggDHTHGHRh8nsjombUBhsAYn1uyvX7Dys7RKZ/fIhl3ZrqE8YyRz9/bpt+AuyW47TYCxqvhMx3HGJwZN6yGDF0eVwkUg0cYHoP9oaQrmTy7NHxwwOjzQrFMNm9v/cb43JK2ytCp6ueJjDENs075/BaW9Wy34TOeySOiXyo2OKcMr3iUTBIoBo8Y7Gghmy+xYGORW7Fk5PPrKiSrAsN2d1tfZejUfEl6hwwtIWnvdZ7ILNEfj7qaz79Z2qIRErEIEw70GdwPGer3H94hDDogMIxpNtyPR24Wq112WpLlimIim9dyQBKMRYpCDixSpHPIsKOliZamsCOhJB3DSBaDHS1VC98urmTydLc1u7akp0VdikFEukXk2yJyzvzbtc5+l8yV2p4TkRNbPb4RcSImuTIgqefD48QAnZXx43Y632axFimyc4xBKaVl1bOFiDDoQC3DRCbPoIYpyRaDnc702Yv6nHo9hkeAJ5VSh4Anzc/r8Ral1DGl1PFtHt9QWA/1lTn7hKR1rl2aCoz+eJSQ2Osx6JypYjHQEWPCxjGGucVlCqWKtgYAGELS/gwdfb0kcCbT8MrckicRgHoVw4PAo+b7R4F3uny8b7EGJW31GKp53noKjEg4RDJhr8BYqfTWs89g/5rAula31zKQsLcSOJdfZr5Q0jZMCkafp+YLFEr2FblNeDQ3VL2KIamUGjffTwDJdfZTwBMi8oyIPLSN4xGRh0TkhIicmJycrLPZ3hNrCtPT1sy4zQIjbk5ipiuDNqfpVr0kjQWG3cuaWoWROlvPuzqNqUDsKnKb0DhV1cIa57KryM0qbvPC6NlQMYjIEyJyeo3Xg7X7KSM5fb0E9XuUUseAB4APi8i9q3fY4HiUUp9SSh1XSh3v6+vbqNm+wO45ZcYzS9oOwloMdrbYOuA+nskTawrR0aJfpbfFQCLG7OKybdMl+MFjGOxooaKwrcjNH322N2XVyzDphopBKfU2pdSta7y+CqREZBDA/Jte5xxj5t808BXgTvOrTR3fqNgdk9Q9awOM8Y8rNha5jWeW2NXRomVxm4XdixRdmVsiHBL62vXL57dYGUOzV0jqOFGihd1puisFfRoqhg14DHi/+f79wFdX7yAibSISt94DbwdOb/b4Rmawo8XWwWcjn1/fBweMPhdKFWYXl20533hG31RVi36bF4sfmzUW6IlomM9vYV0TuzziK3NGcZuOEyVarNTp2KUMjfN4UaNT7531J8B9InIOeJv5GRHZJSKPm/skgadE5Hngp8DXlVLfvN7xO4WhrhYyS8vk8vULSd2L2yyssQC7FOL4nJ7TMNdid5Hb2NwSQ11693kwYW8x49jcEsl4jOaIvsqwLRqho6WJsVm77m3vxpLqGqVUSk0Db11j+xXgHeb7C8DtWzl+pzBsPtxjc0vcNFBfjHwiYxa3aR5KqrWqbh3qqOtcpXKFdM79ueq3ysraz/aEkkZnl3j9DT22nMspEi0R2prDjNlkAIzOLlafF50Z7mqxrc9jc0v0tkddL26DoPLZU4a7WgEYnan/RhqdXTTO2a33w2OFGOywJCeyeSpK71RVgM7WJqKRUNUCrIdiqUIqm2dY05RkCxFhuKuVUZus59HZJd8oButZrJcRD5VhoBg8pNZjqBfrAdxtKhtd6W2L0hwO2dtnzZWhISRbGLFBYExkDGU4rPl1BuO6jMzU3+dSucJEJu+LPlvK0I7kCi+VYaAYPKSnrZlYU8gWC2NkdpFwSLRO5wMIhYShrhZbvCRL6OiuDAF2d9tjPY/OGX3WfYwBDCE5ZoOQTOUKlCrKJ31uYbFYrju5olxRXJlb8kwZBorBQ0SEoc4WewSGDzJVLOyynkdmlxDRc0Wv1ezuarXFerYGNv0SVskVSmSX6ptBeNT8v/mhz0NmiK9eYy+dy7NcVoHHsFOxKw7rl8E5MKxnO4Tk6OwiA4kY0Yj7g3NbZbjLmGY9s1SfJTlqKUPNx1VgRZDXawSMVpWh/p5hddywzmd61GMDIFAMHmNXFsPo7BK7u/V/cMCwnmcXjblv6mF0ZskXYSSgem3qtSTH5ow1CXRO27RYEZL19xn0nvbEwgp31dvnajJJEEramQx1tTCzUKxrwZ5CqcxENu8jj8G0JOv0GrzM2tgqVeu5zrEVwzP0iTK0zXpepD8e9YVn2NHSRCIWqb/PM4HHsKOxHvJ6vIbxOaOGwW8Cox7FUCxVDGXoIy8J7LGehzRPVbVItESIR20Qkj5JVbWwBt3rYXR2ib64NzUMECgGz6mmrNZxI3kdj9wqVlhlpI4+G/Mt+afPna1NtNcpJMsVxficfzxDETMDzYYxBr8YPWDVMtSpGOa89YYDxeAxwzbEJK3BPb+MMXS1NtHWHK7LY/BL3YaFVctQz3VOZfO+Sdu0qDdNt1xRjGf85zGMzi7WlaY7MuOtMgwUg8f0thkDifU8PKOzi0RCQjKu72ybtYiIKTDsUIb+Ehj1jDFY4Ua/Wc8jM9sXkitpm/7q84K5lsJ2WKlhCDyGHUsoJAzXWcswOmusw+CHGgaLeoXkyIxR0KfzNMyrseo3tiskX532Tz6/xXBXa11C0nou/OQlrWQmbe/+tjzDQDHscIa7W3m1zrCKX0IqFru76xOSo7PGWrh+Uoa7u1vrqoq9PL1ASPwTPoPaUOn2hKQ/lWF99Rs61G3456lqYPb3tHJpemHbQnJkZtE3mSoWu7sMITmzUNzW8SOzi74SkAC7u+pL0700vchQV4svahgsrGu0XcPHj8pwb08bAJemF7Z1/KvVqV4Cj2FHs7enjVy+tC0huVAokc4V2Nfb5kDLnKPezCRjcM5nyrC7fiG5r8df13lvj9Hn7QrJi9NG3YaflGF7NEJfPMqlqe31+dLUAuGQeJpM4p//dgOzv9eyMLYuMKwHbr/PFMOeOoRkLr/M1Lz/lKElJC9vU0heml6snsMvtEUj9NchJC9PL/iuzwD7elq5NLU9A+Di9ALDXS00eRgmreuXRaRbRL4tIufMv11r7HOjiDxX88qKyO+a331URMZqvntHPe3xK1WrahsPz0XzGD9akiJwcXLrfbYeuAM+UwytzREGEjEubOM6zy0WySwt++46A+zrbavep1tBKcXFqQXfGT1gPI/b9ZIuTXnvGdarkh4BnlRKHQKeND9fhVLqjFLqmFLqGPBaYBH4Ss0uf2F9r5R6fPXxO4HhrlbCIdmWJWkpk329/rKqYk1hdnW0cHFqfsvHXjCP2d/bbnezHGf/NoWk5U3u9aFiONC7PSE5s1Akly/5ss/7ettI5wpbnupGKcUlDZRhvYrhQeBR8/2jwDs32P+twCtKqct1/m5D0RwJMdTZwsVthJIuThkzjLY217VKqycc6GvblvV8cWoBEXwZYtjftz3FYBkN+3zY5329bUzNF8lucW1zSxnu95nRAyse/FYV4uR8gYVi2fPrXK9iSCqlxs33E0Byg/3fA3x+1bbfFpGTIvLZtUJRFiLykIicEJETk5OTdTRZT/b1tm0zlDTvO2/BYn9vGxcnt56NdWlqgV0dLZ7NI1MPB3rbmFtcZnaLiQZW+Mwv1e21VIXkFu9va39/egxWeHhrxp61//4+b73hDRWDiDwhIqfXeD1Yu58ynu51n3ARaQZ+Ffh/azZ/AjgAHAPGgT9b73il1KeUUseVUsf7+vo2arbvMAartiEkpxd9GVIBQ0jmCiUm5wtbOu7i1AIH+vwnLGAlSeDiFi3JVybnGe7yqTI0r9VWPaVLZqqq37LPYPspq5Yy3K/7GINS6m1KqVvXeH0VSInIIID5N32dUz0A/Fwplao5d0opVVZKVYBPA3fW1x3/crC/nZyZerpZMovLzCwUfelqw4pVtJUBaKUUFzSIwW6XqmLY4qD7ufQ8B/v9aQDs6TYTDbaoGM6n59nb0+aL6bZX0x6NkExEeSW9tTG0C1MLNIXF87Un6g0lPQa833z/fuCr19n3vawKI1lKxeRdwOk62+NbrIf+bCq36WPOTxr7HvCxxwBsaZxh2hyQ9DprY7vs7jYSDS5sYdC9XFG8MjnPIZ8qhlhTmOGuFs5vUUj6WRkCHE7GObfVPqdyHOht97yiv95f/xPgPhE5B7zN/IyI7BKRaoaRiLQB9wFfXnX8n4rIKRE5CbwF+Fd1tse3HE7GATiX2vyNdNbc98aBuCNtcppdnS1EI6EtWVWW4jyU9KfAaAqH2NfTuqXrPDq7SLFU4VC/P68zwI3J+Jb6XCxVuDS1wGGfXmeAQ/1xzqfnqVQ2Hx4+m85pcW/XlcqilJrGyDRavf0K8I6azwtAzxr7/UY9v99I9LQ109XaxLn05j2GMxM5WpvDvpsOwyIcEg72t3NmC17SmQlj3xuT/hWSNw0kOH0ls+n9LYF6g4+t50PJON87O8lyubKpwq1L0wuUKsrXyvBwsp2l5TKjs0vs2USW0WKxxMjMEr/22t0utO76BJXPmiAiHOrfmlV1Lp3jUDJOKCQOtsxZbhyIV4X9ZjibytHV2kSfT6YYX4vDyTivziyyWNxcjvv5SeOe8HNY5cZknOWy2nRmkvUc6GA9bxer7Zs19qw+6+AlBYpBIw4l2zmbym06M+nMxDyHfSwsAG4aiJPOFTY9T9TLEzkOJ+OI+FsZKrX5sOG51Dz98SgdLU0Ot8w5LCG5We/wbCpHSOAGj9M26+Gg6e2c3eR1tsKkhzXwhgPFoBGH+tvJ5jeXmTSzUGRqvuDb8QWLmwYSALw8kd1w30pFcXYix00+77N1zTbrKZ1N5bQQFvVwQ187Idm8kDyfnmdPd6sv03MtOlqaGEjEOLdJZXguPU9zJKRF3UagGDTiRlNIvji+sZDUybqoh5u2ICTH5pZYKJY57HPFYAi80Kas52KpwpmJHLcMJVxomXPEmsLs62njzCYMADCeAb8bPQCHB+K8tEkD4MxEjhv62glrEBoOFING3Go+/KdHNx6YfOGK8YD53Xrui0fpbmvelGKw9vF7n8MhYzxpM17S2VSOYrnCrbs6XGiZs9w8mNiU0ZPNL3NxaoHbhjtdaJWzHB1KcC6VI79cvu5+SilOj2U4MqiHARAoBo2Ix5o40NvGybGNFcPJ0TkGO2L0+2hpy7UQEW4aiG9KYJwcnSMkK56Vn7l1KMGp0cyGqYwvmNlLtw75XzHcNtzByMzShuNJp8cap89HhzopVRQvbXB/X8nkmV4ocvtuPfocKAbNODrcUd954/cAAAvCSURBVH0wrsfzI3PcNqzHTVQvtw138uKV7IZW1XOjGQ4n47RH/Tdh4GqO7e4kmy9tOGXC6bEs8WiEvT6cI2k1lgfw/Ojcdfez7v+jDaAYrGf01AbP9MmROXN/PbykQDFoxtGhDsYzeSavMwA9t1jk0vSiNjdRvbxmj2FVXU8hKqV4fmSOY7sbo8/HdhvzRT43cn0heWosw5FdCV+nJFscHe5ABE6ObCAkRzMMdbbQ3dbsUsucY7AjRm97Myc3CA8/P5qhKSzcPKhHmDRQDJphWUnXE5LWTdYoQvKOPYaQfPbV9YXkxakFMkvLDdPng/3ttDWHr6sYCqUyL41nGyKkAsb8QQf72jf0GE6NZRrGGxYRjg51cGoDxXBydI6bBhLazAsVKAbNuHWog0hI+NmlmXX3ed4UJo0iMPriUXZ3t/DsyOy6+1gC9NiexlAM4ZBwdLjjuorh+ZEMhVKFO/d3u9gyZ7ltuJPnR+bWrdVJZ/Ncnl5sGAMA4PbdnZxL58gsrb0eRbmiODWqlzIMFINmtEUj3L67kx++Mr3uPj+6MM2NybivC55Wc8fuLn5+eX2B8czlWdqaw76eImE1d+zp4sUr2XVX+frRK9OIwF37r5lNxrf8wr4upheK604u97R53999Q6+bzXKU1x/ooaLgxxfWfqZPjWXIFUpaGQCBYtCQN9zQw6nRuTUtjMViiROXZrn3cOM8OAB37u9mIpvnlTWmo1ZK8b2zk7z+hl4tcrzt4p6DvZQqqioMV/PjC9McGUzQ0do4BsAbDxtrqXz/7NqLbf3w/BQdLU0c2eX/zDOLO/Z00dIU5ofnp9b8/qlzxv/inoP6PNOBYtCQuw/2UlHw04vXhpN+cmGGYrnCvYcba7Git9zUD8A/vZy65rtXJucZnV3izTc2Vp+P7+uitTnMd89cu4xJfrnMM6/OcteBxvEWAIY6W7ihr40fnLtWSCplKMnXH+hpKAOgORLidQe6eWodxfD9c1PcsitBT7s+838FikFD7tjTSawpxA/OXWtVff/cJNFIiF/Yp4/baQdDnS3cNBDnyZeuFZLfPWP8HxpNMUQjYe6+oZfvnpm8JoT2k4szFEsV7r6hsRQDwBsP9fGTi9PXpCdfml5kbG6JNxxsvD7fc7CXC5MLjM0tXbV9vlDi55dneeMhve7tQDFoSDQS5t5DfXzz9ATlmgIopRRPvJTirgM9vp5DZj3eenM/Jy7Pklm8OoT2Ty+nOdTfznCX/3P5V/OWm/oYm1u6Zg6hb54ep605zBs0Ci/YxZtu7CO/XOGpVV7D46eM5eMt77GRePONRp++eXriqu3feTlNqaJ4k2YRgLoUg4j8moi8ICIVETl+nf3uF5EzInJeRB6p2d4tIt8WkXPm36562tNI/Mrtu0jnCjz9ysrD87NLs4zMLPHgsV0etsw53n5kgHJF8Z9PXqluG5tb4kcXpnng6OB1jvQvbz8yQCQkfPnno9Vti8USXzs5zn1Hkg1pANxzsJfutma+/OxKn5VSfOXZMV6zp7MhDYCD/e0cGUzwlWdHr/IOv/LsGAOJmFYDz1C/x3Aa+OfA99fbQUTCwMcx1nw+ArxXRI6YXz8CPKmUOgQ8aX4OAO47kqSnrZnPPnWxuu0zT10gEYvwS7cMeNgy57htuIMjgwkeffpS1VN69OlLCPBrrx32tnEO0ReP8tab+/lPJ0bI5Q1P6Us/HyOXL/HeO/d43DpnaAqH+K9eM8S3XkgxMrMIGHH28+l5fv11ez1unXO893V7OD2WrY4dvjI5z3fOpPmvjw9rN6ZSl2JQSr2klDqzwW53AueVUheUUkXgC8CD5ncPAo+a7x8F3llPexqJWFOY//ae/XznzCTffjHF985O8q0XUnzgDftpa4ApIdZCRPjwWw5yLj3Po09f4sxEjkefvsQ7jw2xuwGmhFiPD7/lIHOLy/zZfzlLOpvn/3ziHMf3dmlnRdrJB+85QDgk/Lv//AK5/DIf+/qLDHW28Cu3N6ZnCPDu1wzTH4/yx19/kflCiY8+9gItTWHed/c+r5t2DW5ImCFgpObzKPA6831SKTVuvp8Aki60xzd88J79fP3kOA/93QnAWK/ht958g8etcpZ3HB3grTf180dfe5GQQHdbM488cJPXzXKU24Y7+cDd+/ibpy/xdz++THM4xB+/81ZfL0a0EQMdMX7v/pv446+9yLE/+jZKKT77gV/QpvLXCVqaw/zxO2/l4b9/hts++i0qCj72rlvp1SgbyWJDxSAiTwBrxS7+QCn1VbsaopRSIrLuVJMi8hDwEMCePY3pYq8m1hTmcx96HZ/6wQXKFcXDb7qhIWPOtYgIf/Xrr+EzT11gMlfgQ2884PsZZDfDv/3lI+zvbeNMKsev37mHmzWZftlJPnjPfnramvnRK9P86rFdDTnQvppfumWA/+cDv8A3Tk3wxsO9/PJteo4XymaXkbzuSUS+C/zPSqkTa3z3euCjSqlfMj//PoBS6j+IyBngzUqpcREZBL6rlLpxo987fvy4OnHimp8KCAgICLgOIvKMUmrdRCELN9JVfwYcEpH9ItIMvAd4zPzuMeD95vv3A7Z5IAEBAQEB26PedNV3icgo8Hrg6yLyLXP7LhF5HEApVQI+AnwLeAn4olLqBfMUfwLcJyLngLeZnwMCAgICPMSWUJLbBKGkgICAgK2jUygpICAgIMBHBIohICAgIOAqAsUQEBAQEHAVgWIICAgICLiKQDEEBAQEBFyFL7OSRGQSuLzNw3uBtVfMaFyCPu8Mgj7vDOrp816l1IZzfPtSMdSDiJzYTLpWIxH0eWcQ9Hln4Eafg1BSQEBAQMBVBIohICAgIOAqdqJi+JTXDfCAoM87g6DPOwPH+7zjxhgCAgICAq7PTvQYAgICAgKuw45SDCJyv4icEZHzItIw60uLyGdFJC0ip2u2dYvIt0XknPm3q+a73zf/B2dE5Je8afX2EZHdIvIdEXlRRF4Qkd8xtzdyn2Mi8lMRed7s878ztzdsny1EJCwiz4rI18zPDd1nEbkkIqdE5DkROWFuc7fPSqkd8QLCwCvAAaAZeB444nW7bOrbvcBrgNM12/4UeMR8/wjwv5nvj5h9jwL7zf9J2Os+bLG/g8BrzPdx4KzZr0buswDt5vsm4CfAXY3c55q+/4/APwBfMz83dJ+BS0Dvqm2u9nkneQx3AueVUheUUkXgC8CDHrfJFpRS3wdmVm1+EHjUfP8o8M6a7V9QShWUUheB8xj/G9+glBpXSv3cfJ/DWOdjiMbus1JKzZsfm8yXooH7DCAiw8A/A/5jzeaG7vM6uNrnnaQYhoCRms+j5rZGJamUGjffTwBJ831D/R9EZB9wB4YF3dB9NkMqzwFp4NtKqYbvM/B/AP8GqNRsa/Q+K+AJEXnGXOseXO5zpN4TBOiPUkqJSMOln4lIO/Al4HeVUlkRqX7XiH1WSpWBYyLSCXxFRG5d9X1D9VlEfhlIK6WeEZE3r7VPo/XZ5B6l1JiI9APfFpGXa790o887yWMYA3bXfB42tzUqKREZBDD/ps3tDfF/EJEmDKXwOaXUl83NDd1nC6XUHPAd4H4au89vAH5VRC5hhH5/UUT+nsbuM0qpMfNvGvgKRmjI1T7vJMXwM+CQiOwXkWbgPcBjHrfJSR4D3m++fz/w1Zrt7xGRqIjsBw4BP/WgfdtGDNfgM8BLSqk/r/mqkfvcZ3oKiEgLcB/wMg3cZ6XU7yulhpVS+zCe139SSv1LGrjPItImInHrPfB24DRu99nrEXg3X8A7MDJYXgH+wOv22NivzwPjwDJGjPGDQA/wJHAOeALortn/D8z/wRngAa/bv43+3oMRhz0JPGe+3tHgfb4NeNbs82ngD83tDdvnVf1/MytZSQ3bZ4ysyefN1wuWnHK7z0Hlc0BAQEDAVeykUFJAQEBAwCYIFENAQEBAwFUEiiEgICAg4CoCxRAQEBAQcBWBYggICAgIuIpAMQQEBAQEXEWgGAICAgICriJQDAEBAQEBV/H/A1iDwm4iFR1fAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plot(sine_sd[0:500]);" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "By clipping the tops of the sine wave, we get a square wave.\n", "\n" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "slideshow": { "slide_type": "fragment" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAD8CAYAAABkbJM/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJztnX2wXOV93z/f3XsvYLANBEWIFxWRyM6ITEocDfY0ieMGiAXJWLgzdWGmtuJ6RqFj8tpOIpe+uDNNS93YyXhKoCRhKk8cExpC0GRUE6HmZZqGGIExRtgyAkOQLJCCY3CMkO7d/fWP85y9y2Wv7kV79/mu9jyfmTt79uw5+zvP3d3ze57fqyKCQqFQKBSWouW+gEKhUCicGhSFUSgUCoVlURRGoVAoFJZFURiFQqFQWBZFYRQKhUJhWRSFUSgUCoVlURRGoVAoFJZFURiFQqFQWBZFYRQKhUJhWUy5L2AlOe+88+KSSy5xX0ahUCicUjz88MN/GxGrljpuohTGJZdcwp49e9yXUSgUCqcUkp5dznHFJFUoFAqFZVEURqFQKBSWRVEYhUKhUFgWRWEUCoVCYVkUhVEoFAqFZbEiCkPSJkn7JO2XtG3A698n6a8kHZP0r5dzrqRzJe2S9GR6PGclrrVQKBQKJ8fQCkNSG7gVuAbYANwgacOCw74J/Bzwa2/g3G3A7ohYD+xOzwuFQqFgYiXyMK4A9kfE0wCS7gI2A0/UB0TEYeCwpJ98A+duBt6TjtsO/BnwKytwvYty36MHeerw349SxGuReP8PXsi6887MJ3MBTRzzN79znN998FnmOt08AiU2X34B37PqrDzyBvB3acyzucYMnP2mGX76H11Cq6VsMvv55neO89mGjXnUrITCuBB4ru/5AeCdK3Du6og4lLafB1YPegNJW4GtAGvXrl2m2NcTEfzS3V+i0w2U6bOOgJePzvLx912WR+Dr5Ae/+PuP0g0aM2aAP3j4OT6162tZx7zv+Zf5Hx/cmEfgAD6/93k+uetrQJ7POqJ6fPfbVvG93+1RlPc3cMyj5pTI9I6IkBSLvHYHcAfAxo0bBx6zHLoBnW7wS1e/jZ+7cv3Jvs0bYuN/eoBjc/lmPwuJqMb9i1e9jZ+/Ks+Y3/Wfd/PK8bksshbj8YMvc+HZZ/CX2348i7x/+0df5p6HD/LqbIfTp9tZZC7kO8eq//mX/sNP8NYzpkcu7/OPH+LG330k6+x+IcdmOwB88d9dzTlnzoxc3jiMedSshNP7IHBx3/OL0r5hz31B0hqA9Hh4yOs8IXPd6kNuZ1xKTreVzywygE6aEuVcPZ8x0+borPcH9fg3XuKyC96STd4V676Lo7MdnnnxO9lkLqSemJw+nScwspWm9J3uSc/hhmYuyZ5q5/mCt1vV/3au4xvzqFmJb89DwHpJ6yTNANcDO1bg3B3AlrS9BbhvBa51Ueov9lTGu+dUW70vtYN6zDntradPtzl6vJNN3kIigq//7XdYvzqfyWAm3bC6Rj356mwHCWbaeRRGPfHqhu/7PZtu3NOZxlzfO+acH/SIGdokFRFzkm4C7gfawJ0RsVfSjen12yWdD+wB3gJ0Jf0CsCEiXh50bnrrW4C7JX0EeBb4wLDXeiLqL1fWFUarZV2+1r/lnGM+Y7rFq7M+hdGNatynTeUzDdWzbefN89XZDqdPtVEmx009CbGuMNJvK9cksF7JOMc8albEhxERO4GdC/bd3rf9PJW5aVnnpv0vAleuxPUtB8cKY7rdsi5fXSYpr8LIP+bxUBjdbOYogPYYjHm2m3cSWMuZLSapyadeRk5lWr5WsmRdYfRMUrnChYDTp9ocNSoMhxmuPQaz7aOzHc7I6HCvx+ycEM11uky3lW1VNZV8GJO8wigKI+HxYbR6syAHEfnNcKfPeBVGPeHNqSRrUc77SO4IrZ7T27jCmOtG7yaeg9okNck+jKIwEnMWH4Y5Ssqwwjhjus2rRqd3fQNrZxzzODiAX53tclpGhTE1Bo7+2U43W4QU9Dm9i0lq8skdglfLGgsfRlant3eFUd+0M+qLeR+GcYlxbK6T1YcxFiuMTmSLkIJ5k5Qz8nHUFIWR6NQ+jIxL2Ol2i1njFKwXJZVzhTHT5lVjHkY3syMUxuPmWUdJ5aK3qrLmYXSzh8nXcieVojASc6YoqfFweueTefpUi6OznZ7/JDf1/SuvwqgejfrCFiXldAAfn8u7whiH4IZRUxRGwuHDmGqZTVKOxL2ZapbrKolSjzlX5AyMx43kaG6nd7qzeJ3eeX0Y0yXTuzk4fBjuFYbFJJVuWq5s7zA4vTUGOQmvmsJqvYl7kdVi0C4mqebg8GHYS4P0nN75ZNY3rVfnPArDkazYxCipcTBJzXa6mZ3etcIoK4yJp15GZvdhGKvVWsJqZ7wrjPq3nDVxrxcllU3k6zg2mzlKagyU5Fw3LGG1xYfRADqG6JnpthqXuFcXv3P5MLoGJVmLskZJzeX1YYzLCiOrxSDJKqVBGsCsIw+j1RqT8uYZI4bMM89uT0nmk1krZFtkWDeY7QSnTTUrYqjKw8gfVtspPozJp9Prh5HZhzEOUVI5FYbZPOMcs2tu0NTs9ioPI7+SLD6MBuDwYcyMS+Je1kJ81aN7hZFTYYzNmC0FF7OJfB2zHY8Po4TVNoCOqTSI097pSNyTOeu55/RuUFhtPSexrKrMKwxH4l5ZYSyBpE2S9knaL2nbgNcl6dPp9cckvSPtf7ukR/v+Xk7NlZD0cUkH+167diWudTEcmd5TrRadbths245aUm1zXaX54IZ8Mt29ITpGv03HuMTInYchKSXjTq4PY+gGSpLawK3A1cAB4CFJOyLiib7DrgHWp793ArcB74yIfcDlfe9zELi377xfj4hfG/Yal8OcwYdRO+RmO8HMVMZpfsKRxDbfTCibyNcwX3ywOT4Mixmut8LIJvJ15M7DgEpRlrDaE3MFsD8ino6I48BdwOYFx2wGPhMVDwJnS1qz4Jgrgaci4tkVuKY3jMOHUTdrcmWGdhymirpkhOlHVf+rsyrJ2odhG7MjGu61sh3kzsOA1EWzKIwTciHwXN/zA2nfGz3meuBzC/b9bDJh3SnpnBW41kVx+DDq2Y/LjzFfSyqfzPpG7Ss+mH/M7hatjoKLPZOUubx5zigpqMY9ySapsXB6S5oB3gf8r77dtwGXUpmsDgGfXOTcrZL2SNpz5MiRk76GOVPiHmCrJ9V1mKTMNxJH7on75ukIbhiHPIzZ1KI1J1Mtb7mfUbMSCuMgcHHf84vSvjdyzDXAIxHxQr0jIl6IiE5EdIHfojJ9vY6IuCMiNkbExlWrVp30IOpZgSMz1BWG5wi3dPswwqAw3C1aLWG1Y9A0ymGSmmoXH8ZSPASsl7QurRSuB3YsOGYH8KEULfUu4KWIONT3+g0sMEct8HG8H3h8Ba51URwrjCnzCsOTxFY9+qKkqses5hnzzdOTe+I3SeUuDQLVJHCSS4MMHSUVEXOSbgLuB9rAnRGxV9KN6fXbgZ3AtcB+4BXgw/X5ks6kirD6mQVv/QlJlwMBPDPg9RWlvnnmXMJO98ohu8tkNKc3hLVFq9kklbuku9Ss0iBQrzAm14cxtMIAiIidVEqhf9/tfdsBfHSRc78DfNeA/R9ciWtbLh4fRu30dhXiqx5z2rbdN8+u4ebZMivJ+l+d0yQF1f/YqjC63V4kYi7axYfRDObDag31811RUo7ig26FYbh5ulu0OpzeUP2PXSapiKrg4nTmQbu7aI6aojAS9TKykbNtg0nKNQlrYpSUw/QI1QrD57dJ12DwYZQVRgOY61b2TkcGsHu2nfNGUoty+zAaNTEwZLdDnfWcVWQPRwkYqLtoTq4PoyiMRKcb+WdgYzPbzifT3g/DsKpqmaOkOobsdqj+x24lmd1vU0qDNIO5bv6sULln285+GG4fhiOU2JyHkXu27bx5OkKJAaZbreLDaAJzna5theEuk+HISXCbKiw+DHOF3twmqZZ8EUM9c6thVVVMUg2g9mHkZL6KqbtkhKEonTF6BvIqyTonwTUxCNvN01/GPvOQS6Z3U+h0I/vy1V0mw1oaxHUjMfhtKnm+ENP5vid55baNY3ZMDKCaHEywRaoojJpu5Hd6z9u2m7Nstzv6u/mVJKQQU/eYc0+IWr6wWteY28aVZA6Kwkh0uoYvlzliyNOiNck2m2dyf9aSb1Xlmm23jYl7jgRNSCvJYpKafCIi+5JdZh+Gs4qprS2toTQIeENMfbPtcYiSyiu31fKtJHNQFEaiE/l9GPNRUlnF9rDUVRoTJZnbGVrNPPPKrHGEEsN45GHknhi0jCvJHBSFkeiG58sFxigpR07CmFSrdfir3DfP3LPttrGukm1VZVSSOSgKI9HthmXWCcYaQ44Wre5VlWm23RoDk1R+Jekbs6tCr4yRYTkoCiPhiZJqXuJeb1XltucbQkztK4wGlclwVehtS7bJUA6Kwkg48jDmM4Cziu3RxPLmjhatkGaeNh+GL6zWlZPgND2WKKklkLRJ0j5J+yVtG/C6JH06vf6YpHf0vfaMpC9LelTSnr7950raJenJ9HjOSlzrYnTDYe+sZTewlpR55tmkrOeuq/ig0QHsqtDrND3mYGiFIakN3ApcA2wAbpC0YcFh1wDr099W4LYFr//jiLg8Ijb27dsG7I6I9cDu9HxkdI1htfbEvQb1w7D5MIwmqY4pMsxbfDBdg+NzLiuME3IFsD8ino6I48BdwOYFx2wGPhMVDwJnS1qzxPtuBran7e3AdStwrYvSjTDMwJqXuOeODJu35+eV6ywN4ijpDuZyKEYfxgTrixVRGBcCz/U9P5D2LfeYAB6Q9LCkrX3HrI6IQ2n7eWD1ClzronS6YanmWcnOKrZHN6rIsJzjrgvx2R3A2U0VzYsMcxbic5WAabV8AR05mHJfAPAjEXFQ0ncDuyR9NSL+ov+AiAhJAz+FpGS2Aqxdu/akLyLCMANz+zAMqyrwRgz1mgkZakn58m3qFUZeuWMRVmuYBJZaUifmIHBx3/OL0r5lHRMR9eNh4F4qExfAC7XZKj0eHiQ8Iu6IiI0RsXHVqlUnPYgqSuqkTz8p/A7g/D8ocGc9+zK9mxYZ5rTnO5VkiZI6MQ8B6yWtkzQDXA/sWHDMDuBDKVrqXcBLEXFI0pmS3gwg6UzgJ4DH+87Zkra3APetwLUuStdYGsRZ3jy3LR9q84zZnt+g6Blv1nNWkT2cfcwnWF8Mb5KKiDlJNwH3A23gzojYK+nG9PrtwE7gWmA/8Arw4XT6auDe9KFOAb8XEZ9Pr90C3C3pI8CzwAeGvdYT0Q1ji1bjzdNhknLOwnxRUvPhrblxRMOBNyfBNTFwViXOwYr4MCJiJ5VS6N93e992AB8dcN7TwD9c5D1fBK5cietbDp1ucNqUJ0rKVrnVsKoCb2Od+WZCzYkYqm9gTTLDuZSk0z+Xg5LpneiG5wcF3llY7hsnkNqVZhcLpDHnH7LVGTpvz29Q7olLSRp7gOSgKIyEpZaUu/ucYczgTuhq5pjB1QMkq8geYRpz5ejPKjIrRWEkHE7vXotWoz3fMdu2NhMymeGq8ubZxQL9Jqnm2PN9pkdfmHwOisJINLFFa9dQcBGqG5czPt+iMKzNhKpHR4tW95gdv+likmoAEWGJ2QZvqW+Leca4bLeO2R5Wm1fuODj6HWOO8AWyjJqiMBKO8ubz/TCyiu3hqNALKdzSmN1uGLI5lNgXGeYLJfY5+iv5WcVmoyiMRJXE5vFhOG8knsQ9XwZw17TCaLW8CZpg6m9tX1U1q2XBqCkKI+Hph+ENq+0YE/ectm1XORRnCZj6GnLSRB9GHVgwqeVBisJIVIX48spUzyTlXGGYQkyNocSWZEXrzdNT0n0sugzmbsVr7lk/aorCSDh8GOCNqnCEEoPXVBHhSdyTjO1KjeYZe796U6j8pEZKFYWRiMjvFARvwxWrScpphrN8zt4SMNU1NKccSseUe+LuWT9qisJIOMqbgze5qetSkmbbtq+ku9me76ifZRpz/fWyRUkVH8Zk4ywZ4U3cyy/XatvuGiPDGlY/q92SzZZva9FqLvczaorCSFTx+c1qJuSrJeXtMmjz2zSsfpa1vLmtaVT1WKKkJpxu5LfxgtcB7DLPWLOeTZ+zM7ih45oMjUVkmKegaMn0nnBcPgzrj8pqkmpWprezfpazUZY7D8Ph6IcSJXVCJG2StE/SfknbBrwuSZ9Orz8m6R1p/8WS/lTSE5L2Svr5vnM+LumgpEfT37Urca2LYctJMNcYcvltnP0wfPWzGlaV2BwBCAYfxoSXBhm6456kNnArcDVwAHhI0o6IeKLvsGuA9envncBt6XEO+FcR8Ujq7f2wpF195/56RPzasNe4HJyVW53JTbZaUs5yKA0rb94xNcpyfs5hMknVX60SJbU4VwD7I+LpiDgO3AVsXnDMZuAzUfEgcLakNRFxKCIeAYiIbwNfAS5cgWt6w3Qjfwge+JObbGUybKsqX3lz583TUz/LZ8/31ZIqeRhLcSHwXN/zA7z+pr/kMZIuAX4Q+Ou+3T+bTFh3SjpnkHBJWyXtkbTnyJEjJzcCasfgSZ9+0jhj1V0mKa9t2xRWa27R6poYgGeVYfdhlBXG6JB0FnAP8AsR8XLafRtwKXA5cAj45KBzI+KOiNgYERtXrVp10tcQ4XQMZhcLeBP3rC1aTZFhvqxnUzScMSehnpAod4+bkoexJAeBi/ueX5T2LesYSdNUyuKzEfGH9QER8UJEdCKiC/wWlelrZLhqSbXsOQn55TqT2DpdV4ip7ybiq59VPTq+3/WExFVLqpikFuchYL2kdZJmgOuBHQuO2QF8KEVLvQt4KSIOqfrl/g7wlYj4VP8Jktb0PX0/8PgKXOui2Gbb7iiphuWehMlX1cz6WT57vq1F64TXkho6Sioi5iTdBNwPtIE7I2KvpBvT67cDO4Frgf3AK8CH0+k/DHwQ+LKkR9O+fxMRO4FPSLocCOAZ4GeGvdbFcLVzrGR6awwVJZkHd05C83wYvpLuMLk+jKEVBkC6we9csO/2vu0APjrgvP8LDPwmR8QHV+LaloOrmic0OXEvv1zwJe65/Tau+lmV/PyyuyaTVOmH0QBcZQSgrjGUXSzgrSXlDCV2maR8vdt9pkfw5CTUk8BSS2plKQqD+Ru2rex1A/MwrGY4V7Ki0QznjZIy+jBMtaQm1YdRFAZ93bkM/42qTIbTJGVKYnPePJtWxt5U9qZnz3coDFdtuAl3eheFgW/5Wsu0Je65TFJG84yzRavL9Njtmir09poJZRft62/Tc3pnF52FojCAcJqkrI11jOYZp5K0mGd8s05fFYPq0ZKHYetvUz2WFcYEM7/CyC/b2w+jiSXd8/d5Bq+vyl1LyjE5CFd/m9qHUZzek8u8D6OBOQkNS2KrTBX55dZRUq5CfK5GWeAJMXX1tyktWhtAffPytWhtljPU2SfBGRlWyc8umo4pQbPO/bA4vU3f7V5YbTFJTS69ypamH5W1+KDFJGUOMTXlnoDH/FgV1swu1hoxZIsALFFSk4/Xh2GuMeQqk2HrDeFx9DtLRtgKa8pnz3f1t3GOOQdFYdBfS6ph8fkNzEnomGbbzpIRNtOj0Z7fMQV0FB9GA+iVBnG1aG1cyQiz38YYbulJYvNFAIIrSsrVcrl6LKVBJhinD6MtX12ljs0xaKyr1HXVDDOapIz1s8DXD8NZDsX1mx41RWEw/yNuWotWa+Keta5Sfrkt+W4k/siw5vkwSpTUBFP/iF3JTc4QU1f9LGeyoqdCr7E3hL34YHbRdLuu7Pbiw1gSSZsk7ZO0X9K2Aa9L0qfT649JesdS50o6V9IuSU+mx3NW4loH4a0l5Yuo6NjsvMa6SvaSEdlF22bbTnu+r4x9kj+hGmNohSGpDdwKXANsAG6QtGHBYdcA69PfVuC2ZZy7DdgdEeuB3en5SHCWN3fNtiPCFmLaNuZhdN0lI2z2/Oxirfb8ju27XfIwluIKYH9EPB0Rx4G7gM0LjtkMfCYqHgTOTj27T3TuZmB72t4OXLcC1zqQ+SipUUlYHJlqDHkd/d5yKE0re+32YbhWGM7PuURJLc6FwHN9zw+kfcs55kTnro6IQ2n7eWD1ClzrQOy1pEzJXOBTkq66Ss5yKOC8eTarHIqz1wuUFq1WUk/wgR+BpK2S9kjac+TIkZN6/5mpFm9bfRZnnbYiLc7fEFW12uxirW1p3c5QZ3y+qxCf1Z5vmhg4xzypUVIrcYc8CFzc9/yitG85x0yf4NwXJK2JiEPJfHV4kPCIuAO4A2Djxo0n9Sl93/lv4U9+8cdO5tShabU8YbW9VZUzic1wI3M5gJ1RUmEqPui053dMZezbRtNjDlZihfEQsF7SOkkzwPXAjgXH7AA+lKKl3gW8lMxNJzp3B7AlbW8B7luBax07qu5zTpNUwxzAdvOMa8zZxVrrZ7lCxusxT2qU1NArjIiYk3QTcD/QBu6MiL2Sbkyv3w7sBK4F9gOvAB8+0bnprW8B7pb0EeBZ4APDXus44mqs04sMa1gGsKtFq1NJukrAOFcYLr/NpNeSWhGjfUTspFIK/ftu79sO4KPLPTftfxG4ciWub5xxJe7Nm6Tyy3Y6gF32/LbVAezq9TIvPze+Cr3z8ieRU8LpPcm4Evc6Rqd3y+n0Du/N06cks4u1lskIV2mQkodRGCWuxD1nSXdXNmwtr2mJe+5yKL62tNnFlgZKhdHiKj5oTdwz3TydCZrzjXXyy/aVQ6lNj9lF+3wYpZZUYZS4Sn07uwzKZKpwmuGcLVpd5VC8Yy79MEZBURhmXKW+rV0GTbPt+t9sbdHaIPOMjOYZd77NpIbVFoVhxubDcJZDMc0861mfpaS7sx+GqxWvuYFSKW++8hSFYcZV6tuZuOdK6HK24nXb812teMEz5iilQUZCURhm2i2fjRdM9vzebDuvXGcZ+5bRnm9rxeses2ky1FJp0VoYEbZM7zpKynjztDm9rVFSns/aWg7F0mXQM2bwtl0eNUVhmKmjpHLPSJzlzV2x6l6/jbtCb3651jGbSsCAt+3yqCkKw4zLSdZTGM5aUqbEPWumt8k8Y23R2qBkRahbFkymxigKw4wrYshZ3rxX6jv7mF8rPycuJWltxWuMDHPVkgJfU7QcFIVhxhcxVD22jDWGckeHjYUPw6Qkm9ai1dUDBHx+yRwUhWHGVSbD2g+jriWV++bpTFY0NVCa99tkFQt4i0y6eoBANe4J1RdFYbiZv3nmlTseDuDm5GHI9Dl3xsBv44kM8+SeQKreUExShVHgWrZbK7eazXBNVJLWtrSmBmEOJQm+6g05GEphSDpX0i5JT6bHcxY5bpOkfZL2S9rWt/+/SfqqpMck3Svp7LT/EklHJT2a/m4f9L6TgKsEdP0jtsw8TaaK+dl2Xrng92E4Jwa+ku7ZxQKpekNRGAPZBuyOiPXA7vT8NUhqA7cC1wAbgBskbUgv7wK+PyJ+APga8LG+U5+KiMvT341DXufY4lthVI/O8gm5f1RhnG27PuexUJKmplHeKCmL6JEzrMLYDGxP29uB6wYccwWwPyKejojjwF3pPCLiTyJiLh33IHDRkNdzyuGabTudoa4WrR2jD6PWUbknnk4l6U3cc0ZJlVpSi7E6Ig6l7eeB1QOOuRB4ru/5gbRvIf8C+N99z9clc9SfS/rRIa9zbHHNtsfDJOVZVTUpSmocouE8DcLcmd6TqTCmljpA0gPA+QNeurn/SUSEpJP6L0m6GZgDPpt2HQLWRsSLkn4I+CNJl0XEywPO3QpsBVi7du3JiLfiKgE9Dk7v3Mt2pwPY7cNwzLYlIVMhPm+U1OQm7i2pMCLiqsVek/SCpDURcUjSGuDwgMMOAhf3Pb8o7avf46eBnwKujPTNiohjwLG0/bCkp4C3AXsGXN8dwB0AGzduPOU+pWZGDNXX4AqrzSq2kmkPJc4qtocria3qh+GMkrKIHjnDmqR2AFvS9hbgvgHHPASsl7RO0gxwfToPSZuAXwbeFxGv1CdIWpWc5Ui6FFgPPD3ktY4lrV6UVF65TmeorUWrsX6Wq9dzx7iSrOW6Ci66akmp+DAW5RbgaklPAlel50i6QNJOgOTUvgm4H/gKcHdE7E3n/3fgzcCuBeGz7wYek/Qo8AfAjRHxzSGvdSxx2Xmt8fmmGkPeMhnVo61plPHm6Svpnl0sUH2/J7UfxpImqRMRES8CVw7Y/w3g2r7nO4GdA4773kXe9x7gnmGu7VTBntDVoE5s1jGb8m2cjn7wJbG5mkZB6YdRGCEyOUOtJSPqBkq2iKGsYpPM5tWSgvrmmV9umDruQemHURghLtv2OJSMyG+SaqAPw5h7Usk1tWjtemtJTWqUVFEYZmy27TrT22mSapB5Rq7IMGMeBvhyEpyJe6WWVGFkuMIt5xP3sooFfF0GxyG73VZLynXzNNRV6hpNj1CZeTuTqS+KwnBjS2Lr+kxSrrLX1ux2k6Pf6beBdPNsUHBDJbeYpAojop7tNqldqSsyLKxRUtVjk3qAQPX9dlVidkZJFZNUYSTYoqTGwCSVP0rqtfJz4qrc6lYYjhDT+qfk9NuUsNrCSGi7biTOWlLuMhlWH0ZeuR2j6RHq2XZemW4zXEv5KzfkoigMM64S0OOQ6Z19zMaIoVqky/RoM8+0mtVlsJZbSoMURkLvRmJrrNOkMhnVo+NGIskSn+8uPtg2mKTqABJX8cHiwyiMDF9dJWOUlLktrbNya5PK2INpzL3ghqxie0xyefOiMMzUN09blJQ1xNQTJWV1hjaoURZ4EvfcUVKlvHlhZLiS2Jzlzd0OYGeZjNyWCmfvdvD0t/ZHhnm6DOagKAwzriQ2Z+KerUzGGGQ9N634oKM3hLtCb/FhFEaGK4mtEz7bdttkkuoaV1XguZG4TVLtVv7eEG4lWRRGYWS4W7RaTFINDCWGZM83+W2cTu8mRQBC8WEsiqRzJe2S9GR6PGeR4zZJ2idpv6Rtffs/Lulg6rb3qKRr+177WDp+n6T3DnOd44ytEF83aMnzo6pF2pyh1lLfeWU6s9vB0xui/lq5lKSry2AOhl1hbAN2R8R6YHd6/hpSb+5bgWuADcCldygEAAAP90lEQVQNkjb0HfLrEXF5+tuZztlA1fv7MmAT8Jt1j+9Jw1VjqBO+nse27HZzyQhHQpczux08/TA65jGX8uaLsxnYnra3A9cNOOYKYH9EPB0Rx4G70nlLve9dEXEsIr4O7E/vM3H4yl6HNbEJHM5Qf5mM/C1azcUHjXkY1vpZRWEMZHVEHErbzwOrBxxzIfBc3/MDaV/Nz0p6TNKdfSatpc7pIWmrpD2S9hw5cuSkBuFERgewza7d8qww/DWGDPZ8t9/GkuntVxi5Q4lzsaTCkPSApMcH/L1mlRDV1OmNfjNuAy4FLgcOAZ98g+cTEXdExMaI2Lhq1ao3erqd+XaleeV2ur6bCHgcg84WreAac/Xo82E4ghuqR5+S9LSlzcHUUgdExFWLvSbpBUlrIuKQpDXA4QGHHQQu7nt+UdpHRLzQ916/BfzxUudMGr66SmGbaUNKbmqYqcLhDHV3n2tJzGXuoOReSRYfxuLsALak7S3AfQOOeQhYL2mdpBkqZ/YOgKRkat4PPN73vtdLOk3SOmA98IUhr3UsaRl9GK6ZNrhqDFWPtk5shhuJO5TYOWbfxCB/l8FcLLnCWIJbgLslfQR4FvgAgKQLgN+OiGsjYk7STcD9QBu4MyL2pvM/IelyKlPWM8DPAETEXkl3A08Ac8BHI6Iz5LWOJbae3kYfBngKtDnLoUDtDM0r010OxdHf2q0wHF0GczGUwoiIF4ErB+z/BnBt3/OdwM4Bx33wBO/9q8CvDnN9pwLzWc955XbDZ8sHjz0/7A5gZ9MoX3/r7FWJxyAarkRJFUZCy9XTu+v1YchQoM2exGZYVbnNcJYoKWMVAyjlzQsjZCppjOzmmfCapLy27axiezjG7HYAO/pbu/02jra0uSgKw0x9055zREk5TVKmhC6ZyqGAxxkaZpOUo6T7XMft6J/csNqiMMy026YyGd2wmWbAc/Ps2ldVPnt+E8uhTJlqgzjMcLkoCsOMa4XRCXfinuPm6btxgscZ2hkDH0bu2fZcz+mdVWyPVkvZV1W5KArDjKsfhj9xL/8sLCJsBenAY9ueN0nllVvjCZ+ulq5t2wojfxBLLorCMFMrjLnMwepuk5Tj5tmxj9lZP6s55VBqU+eUzYfRotONiczFKArDTK80iCNxz5npbXAMdsNnmoFkz3eF1Zo+a0/4dKUxXEqyVlST6McoCsOMpHQjyesB7obXnu/qb20cMjKW+naNu20o6V77MKbavlVV/3VMEkVhjAGVwsgrs+u255vyMKyOfkdP79oB3KDeEO5M7+l2WWEURsiUYYUxFrWkDDcSqw/DUurbnMRm8WHUYbU+HwaUFUZhRLSlZibuNax+lrOZkK+7Yn5H/5zZ0V98GIWR0m47agy5E/eaVz/L0qLVnm9jMD2OjQ9j8mqcF4UxBjhWGG6TVLuV/+Y5DvWz8ifu+ZVk7pn2nNmHUVYYhZHiKcTnS+YCVxXTsJlmoDbP5JU5Dvk2+dsPex39rtyqHBSFMQa0W2pe4l7L0FjHnXtiCqv1jtmTYwS+WlJTJUpqMJLOlbRL0pPp8ZxFjtskaZ+k/ZK29e3/fUmPpr9nJD2a9l8i6Wjfa7cPc53jjstU4Q0xzV9Lqon2fHf9LGdJ97bNhzG5UVLDtmjdBuyOiFuSItgG/Er/AZLawK3A1cAB4CFJOyLiiYj4Z33HfRJ4qe/UpyLi8iGv75TAlQHsNlVkzwA2J+65zHDeRln5o+HmzCap4sNYnM3A9rS9HbhuwDFXAPsj4umIOA7clc7rocqw/AHgc0NezymJRWG4I4YMM88wO70dVUzt4dOWEjDufhiV3Nnc2bgZGFZhrI6IQ2n7eWD1gGMuBJ7re34g7evnR4EXIuLJvn3rkjnqzyX96GIXIGmrpD2S9hw5cuQkhuBnyqAw7LWkDA5ge+KeI5TYrSQNmd61P9CVuDfJK4wlTVKSHgDOH/DSzf1PIiIknex/6AZeu7o4BKyNiBcl/RDwR5Iui4iXF54YEXcAdwBs3LjxlPyEfKYKr217LvMMzJ245ykN4kvaA5hut4jIO0HpFR90KYx2g30YEXHVYq9JekHSmog4JGkNcHjAYQeBi/ueX5T21e8xBfwT4If6ZB4DjqXthyU9BbwN2LPU9Z6KTLU9sequxCbwzDzdZjiHPb+KDMsrs5/pJPz4XJczZtpZZHYibKsLmOwVxrBfpR3AlrS9BbhvwDEPAeslrZM0A1yfzqu5CvhqRByod0halZzlSLoUWA88PeS1ji2OxL3jc11mjHcSR2Mdd4ipy57vNEnNTCWFkXE1OWc2t5ZM78W5Bbha0pNUN/5bACRdIGknQETMATcB9wNfAe6OiL1973E9r3d2vxt4LIXZ/gFwY0R8c8hrHVscoYfH57q9H7MDS2Od8JpnfJFhRoWRVrHH5/LdPDsdr8KY5BXGUGG1EfEicOWA/d8Aru17vhPYuch7/PSAffcA9wxzbacSjsS94x2vwmgZGutUUVJZRb4GR+XWMOee1N+xnBFD9hyj0g+jMEociXuVSSqPTXkQjSxvrvwmqY7Zb9Pvw8hFp+v2YVRj7pTSIIVR4MjDcJukXGUyGhclZR5zrTByrjDGx4dRFEZhBNRN43MREZVJymifcfgwul3sUVIeM5zfJHUs4wrDXTOs1JIqjJR2Znv+bFoqO1cYMjTW8UdJeSq3Os1wM6YVhqvwIJQoqcKIyb3CqEMc3VFSnt4Q5kxvR82wsXB65xt3p+vtV1/7T0p588JIyF0apHZAOvMwHPb88Sjp3qxkRZ/T27/CKCapwkjIPdvuKYwpX5SUJes5vD4MR4tWd4ipJazW7MOYnuDSIEVhjAG5o6TmFYZzFuYJMfX2AGleGfvp5ADO6fSe63btrXhhvqbVJFEUxhiQXWF0OoBXYTS2RWvkbRzlNkmdZllheJMVp0pYbWGU5FYYx3o+DK89P3tYrTnE9E2nVYUVjs52ssl0R4Z5fBhda2HN4sMojJSq+GC+H9Q4hNU6nN5HZzucPu0b85tPrxTGy0fnssnsdM21pAwrjDlzcMPUBLdoLQpjDGi3Rc7WEPNRUs7SIPlnYN96ZZaz3zSTVWY/bz59GoBvvzqbTaZ7VdVbYWT8gnfN5c3LCqMwUipnaL4f1Dg4vXO3aO10g2+/OsdbzpjOJnMhb6lXGK/mW2F0x6T4YE6T1NyYVKsteRiFkdBUp3fOTO96Vn+2UWHUK4yXM64wKpNUNnGvY8awwnBHw7VaQipRUoUR0cjEvcxO72+9Ut2k32pUGG89o1phfDvjCiPGxOk9O5dxNWkeM1S/6dliknotks6VtEvSk+nxnEWOu1PSYUmPL/d8SR+TtF/SPknvHeY6x53ciXvHxsAkJZF1zC8dTSuMN43BCuNoxhXGGPRub7eUPXHP6cMATwXqHAx7x9gG7I6I9cDu9HwQ/xPYtNzzJW2g6sR3WTrvN+uWrZNI7i9XL0rKXBokZz7Ct476Vxh1lFTOFUZVodd785xpt/K2aDX7MACmW63iwxjAZmB72t4OXDfooIj4C2BQi9XFzt8M3BURxyLi68B+4Iohr3VsaWamd94xj8MK44zpNlMt5Y+SMhuep9vK6vR2555AHfk4eT6MoVq0Aqsj4lDafh5YvULnXwg82HfcgbRvIqnt+Vd/6s+zyKtn216TVN4x1wrDGSUliTefPsVn//pv2PXEC1lkPvviK6w778wsshZjZqrFHz5ygL/c/7dZ5D374itcuso75qmWuPeLB/l/T72YTeZ73r6Km39yw0hlLKkwJD0AnD/gpZv7n0RESDrpKePJni9pK7AVYO3atScr3sp7Lzuf/Yf/PmuY6QVvPYNzjLPtTZedz9NH8o951VmnZZM3iJt+fD0PPztosT0a1q8+i3+68eJs8gbxL9/zvdnH/AHzmG/8se/hkb/5u6wyV7/l9JHL0DB2ZEn7gPdExCFJa4A/i4i3L3LsJcAfR8T3L3W+pI8BRMR/ScfdD3w8Iv7qRNezcePG2LNnz0mPp1AoFJqIpIcjYuNSxw1rk9gBbEnbW4D7Vuj8HcD1kk6TtA5YD3xhyGstFAqFwhAMqzBuAa6W9CRwVXqOpAsk7awPkvQ54K+At0s6IOkjJzo/IvYCdwNPAJ8HPhoR+Sq2FQqFQuF1DGWSGjeKSapQKBTeOLlMUoVCoVBoCEVhFAqFQmFZFIVRKBQKhWVRFEahUCgUlkVRGIVCoVBYFhMVJSXpCPDsEG9xHpCnfsH4UMbcDMqYm8HJjvkfRMSqpQ6aKIUxLJL2LCe0bJIoY24GZczNYNRjLiapQqFQKCyLojAKhUKhsCyKwngtd7gvwEAZczMoY24GIx1z8WEUCoVCYVmUFUahUCgUlkVRGICkTZL2SdovabG+5Kccku6UdFjS4337zpW0S9KT6fGcvtc+lv4H+yS913PVwyHpYkl/KukJSXsl/XzaP7HjlnS6pC9I+lIa839M+yd2zACS2pK+KOmP0/OJHi+ApGckfVnSo5L2pH35xh0Rjf4D2sBTwKXADPAlYIP7ulZobO8G3gE83rfvE8C2tL0N+K9pe0Ma+2nAuvQ/abvHcBJjXgO8I22/GfhaGtvEjhsQcFbangb+GnjXJI85jeOXgN+jasw28d/tNJZngPMW7Ms27rLCgCuA/RHxdEQcB+4CNpuvaUWIiL8AFvbG3AxsT9vbgev69t8VEcci4uvAfqr/zSlFRByKiEfS9reBr1D1g5/YcUfF36en0+kvmOAxS7oI+Engt/t2T+x4lyDbuIvCqG4mz/U9P5D2TSqrI+JQ2n4eWJ22J+7/kNoC/yDVjHuix53MM48Ch4FdETHpY/4N4JeBbt++SR5vTQAPSHpY0ta0L9u4p4Y5uXBqExEhaSLD5CSdBdwD/EJEvCyp99okjjuqjpSXSzobuFfS9y94fWLGLOmngMMR8bCk9ww6ZpLGu4AfiYiDkr4b2CXpq/0vjnrcZYUBB4GL+55flPZNKi9IWgOQHg+n/RPzf5A0TaUsPhsRf5h2T/y4ASLiW8CfApuY3DH/MPA+Sc9QmZB/XNLvMrnj7RERB9PjYeBeKhNTtnEXhQEPAeslrZM0A1wP7DBf0yjZAWxJ21uA+/r2Xy/pNEnrgPXAFwzXNxSqlhK/A3wlIj7V99LEjlvSqrSyQNIZwNXAV5nQMUfExyLiooi4hOr3+n8i4p8zoeOtkXSmpDfX28BPAI+Tc9xur/84/AHXUkXTPAXc7L6eFRzX54BDwCyV/fIjwHcBu4EngQeAc/uOvzn9D/YB17iv/yTH/CNUdt7HgEfT37WTPG7gB4AvpjE/Dvz7tH9ix9w3jvcwHyU10eOliuT8UvrbW9+rco67ZHoXCoVCYVkUk1ShUCgUlkVRGIVCoVBYFkVhFAqFQmFZFIVRKBQKhWVRFEahUCgUlkVRGIVCoVBYFkVhFAqFQmFZFIVRKBQKhWXx/wFTPBSKwjjS7QAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "square_sd = minimum(.1, maximum(-.1,sine_sd))\n", "plot(square_sd[0:500]);" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "slideshow": { "slide_type": "slide" } }, "outputs": [ { "data": { "text/html": [ "\n", " \n", " " ], "text/plain": [ "" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "## What does a square wave sound like?\n", "\n", "Audio(data=square_sd, rate=Fs)" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "The Sine and Square waves have the same **pitch** (frequency) but different **timbres**." ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Beats\n", "\n", "Playing two tones at similar frequencies results in **beats**. \n", "\n", "This is useful for tuning instruments!" ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "slideshow": { "slide_type": "slide" } }, "outputs": [ { "data": { "text/html": [ "\n", " \n", " " ], "text/plain": [ "" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "sine2_sd = sin(2*pi*440*t) - sin(2*pi*442*t)\n", "\n", "Audio(data=sine2_sd, rate=Fs)" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "We can plot to see the beats, or variations in amplitude. \n", "\n", "Note we plot the whole 3 seconds of the sound. Two beats per second. " ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "slideshow": { "slide_type": "fragment" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAD8CAYAAAB+UHOxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJztnXu81VP6xz+PLpSSS0dK5SQVFZKjm1uIbkPuU+4GScwYDHK/0zDjN5OQyGBco1FRRKlJLumU7hdOCV3ouCXSJD2/P84+M9tuX777u+7ru96v1361z95rf9f69t37+6y1nuf5PMTMCAQCgUDy2M70AAKBQCBghmAAAoFAIKEEAxAIBAIJJRiAQCAQSCjBAAQCgUBCCQYgEAgEEkowAIFAIJBQggEIBAKBhBIMQCAQCCSUmqYHkI+GDRtyaWmp6WEEAoGAM8yePfsrZi6J0tZqA1BaWory8nLTwwgEAgFnIKJPo7YNW0CBQCCQUIIBCAQCgYQSDEAgEAgklGAAAoFAIKEEAxAIBAIJRdgAEFEzIppKRIuJaBERXZ6lDRHRMCKqIKL5RNRRtN9AIBAIiCEjDHQLgKuYeQ4R1Qcwm4jeZObFaW16A2iVenQG8HDq30AgEAgYQngFwMxrmXlO6vkGAEsA7JnRrB+Ap7iK9wHsTESNRfs2xfLKH1A6ZAJKh0zAP9+PHHJrJePmrv7vuSyv/MH0cGKzdSvj1vGLUDpkAnr9bTr+s+UX00OKzcbNW9DtnikoHTIBf31jmenhCLFw9fr/fr+mLl1nejhCPDSt4r/nsvq7n0wPRwoksyYwEZUCmA6gPTN/n/b6qwCGMvOM1N9TAFzLzNtkeRHRQAADAaB58+YHf/qpXTfY0iETsr6+9I5e2KFWDc2jic+WX7Zinxtey/reyqF9NY9GjE+//hFH3jdtm9dv/k1b/O6wFvoHJMBfJi3D8KkV27z+3nVHo3GDOgZGFA9mRovrJmZ9b/ndfVBjO9I8ovj8+J8taHfLpG1er1u7Bhbf3svAiPJDRLOZuSxKW2lOYCKqB2AMgD+m3/yLhZlHMnMZM5eVlETKZtZGrps/AOx70+vYulWeMVUJM+e8+QP5z9M2vvx+U9abPwDc/upivDDrM70DEuC+SUuz3vwBoOs9b+H7TT9rHlF8ct38AaDl9bnfs40tv2zNevMHgI2bf0HrPL8jF5BiAIioFqpu/s8w87+yNFkNoFna301TrzlDlK2evR35Yuf7cVZz7+tLNYxEnM53T8n7/rVjFuC7jZs1jSY+X6zfhAenLs/b5oBb39A0GjEufXZOwTauTDLyTZQAYPMvWzFh/lpNo5GPjCggAjAKwBJmvj9Hs/EAzklFA3UBsJ6Znflf27qVcdPYhZHa2v5leH/F15HaPTRtOX7+Zavi0YjR9ubXI7XrcPubikciTpd78huyak55+F3FIxFj08+/RP4N2O5zeuztFZHaXfrsHMjcSteJjBXAoQDOBnA0Ec1NPfoQ0SAiGpRqMxHACgAVAB4FMFhCv9ooZmYfZfZjkv4j34/ctpXFy9uNm7dg4+bojt4pS75UOBox/vHOJ5Hbzv70W/xi8VbjvjdFM8oAcMxf/61wJOLcOWFJ5LYH3zlZ4UjUISMKaAYzEzMfwMwdUo+JzDyCmUek2jAzX8rMLZl5/2zOX1vZ9HPx0SQP5tjHNc3TMSKWbI2maXtz9n3ZXFzwpL1fudteWVy4URq27qGv27Cp6M/M+PgrBSMRZ/Azs4tq/82Pm53xAaYTMoELUMyMppr7JtkZundjxG2sdNrcWPz5qyaOUQbsvNk89d7KWJ+zccuh013RtrHSOWvUTAUjEWfigi+K/szh905VMBK1BAOgiHcr7LrZLF4TOzDLuptNHKMM2HmzuXncolifO+khu3wBm7fE9xet+najxJGIEze3x8XcgGAA8iCylXPGY3bdbPoMezv2Zwc/Y7dfoxhsigiq3PCf2J+d+/l3EkcizskPvxP7s4f92a6Zc9SAj2w8/4E7YcdAMAB5Ed3KibtVIRvRGfxrC4tfDqtC1JlrU0TQIXeJOQ5tMgILV8dfYQL2rDJ/+M8Woc8P+dcCSSPRQzAAOfj2R/GZ4m+LiLhRyRUvzBU+RhwHnwpsdubq5sQH48+6ZfLJVz8KH+PpmXbMnNvnSPoqhi2Wh0+nEwxADroOLd6hlck8S2ZoY+euET5GHAefbGTNEkX8IbKYvNjesNRiOeov04SPIbLtYhsDHrVj4heFYABysOlnOVZ842axJaUoPwouaW1iqKTsZBF/iCwufErOSubV+eLG3RZMh1HKckbPWvmtlOPoIBiALHwjYfunmr7DZkg7VhwuknSjAap0d0zyyL+jZWYmicue/dBo/x9+Ju9md9fE6IlXKpDpjBb1JegiGIAsHP+AvJu2jP1REd5dHk36IQqFdHdUInt2uGjNeqnHK4Zpy9yWRU5HZjjqqBnRM6Jt5+J/uuGrCgYgC7LjeU05hWyWDCiWF2d/LvV4Jldm5/1jltTjvbvcrpwTEUxFA339Q/yQ3Gy8UyFv4qWSYAAyULEPGScDVwZPvLtS+jE3GJIkvnaMW+F1OjnjUTM5JyJ5DLmIk4Erg7NHfSD9mC5MwIIByGDYWx9LP+bzs+TOXqNyx6vFacxE4U8vzpN+TFOs36jfmNmUiCbK2Qoyq02JKS5eKz8ybNxc+xXvgwHI4G+T5RsAn5i0SH/44jpFzuerX9JvzE4d8Z6S45ooFrP0iw3a+1SBqpn6laPtnywFA6CJBav0Oh1XGnY+y+SE4WoSnt4wEItfsU6NBv41L85XclwT6Fagdb1WsQjBAKSxQmGBiuOH63U6nviQuizRZZpnfl8oDD+1RYJAlNcX6d07VxnJpHvmLCsnIxumfGZRkVUS8nEiWkdEWb2dRNSdiNanFYy5WUa/simmWIrtfKdwf7vn36YrO7ZuVn2rT8Hx3x9VautLNbIjmdKxvapeMVzzkt0rM1krgCcA9CrQ5u20gjG3S+pXKusURDWkYzrT0UVUhzjqlFU+93H5kSbpmM45cRHVv0mbhBSzIcUAMPN0AN/IOJYpdGwFzNBUI2COxOxM06gOcfxKcvy3Sc7UpEGjYyKjy4f1DwWh0i6h0wfQjYjmE9FrRNROY7+RWF6p/gt3juIZYDUna5jVzoxYXD5QhY6b5pr1eqQ6xsxZpbyPi/9ZXEnGuKgIlc5kjcWFYnQZgDkAmjPzAQAeADA2V0MiGkhE5URUXlmpb8/09EfUhOf5ig6pa12JNDoitKZ6JP9wtYZ97WVf+hFiCtjtB9BiAJj5e2b+IfV8IoBaRNQwR9uRzFzGzGUlJSU6hgdArgCcSXyJagH0FTw54zH1xkxXHQNVYaY+oqtgk66t3zhoMQBEtAcRUep5p1S/idxDUF0rWGdUi2pOeViPg3bDJjeUG6NwmeJM2p816lqpVp8d+6H9mbqqkRUG+hyA9wC0IaJVRHQBEQ0iokGpJqcCWEhE8wAMA9CfLZqq6pQ5Vl0reKCmvVPArKKmS+j8qqvOzn1fo+/nlnGLlB5fZ/nGnzbbUR42E1lRQAOYuTEz12Lmpsw8iplHMPOI1PvDmbkdMx/IzF2YWV/sXQR8qka0RIGmSS4uedqfYvFr16tbOa38Wk6hERtQIZqWC93JbSp5eFqF6SFkJWQCw4wkgA989o26G5sq/Z9c3DZeXTTIJU/rW5UB9s42bUJ3Ts6wt4IBCKSQXW+gGhPqlqq4TUN4XjoqZ5u6RdMefTtUTivEJ1+HpDkgGABtkQDp3DVBzc1t5NvLlRw3H6pqDvskB6Cb+9/8SMlxvzUQKffKPDU1jy9/3mwpTVtIvAGYuED/jUZV0YsHp+o3AE+//6n2Pl3CxARDFfdOWqa9z98/p+ZGvXC1Pl9ZNTod6FFJvAFwQbPbZu55ban0Y5qqpPT2x/ITD5+d+Zn0Y5riuQ/8ORcTDNLsC4pC4g2AL1gUVSuMqaLpKhQub9fsy6gmhOjmxtQEQ6VCb1yCATCE7CzXjz3KANWVNZuJCzVco3K95Bh3kxMM2UETL802U6LVRhJtAEyU0avm98/JjaE3qTfyucJw0EA85knWN5q10pzC7Ihpcn1b147RlwCWicl7TjYSbQAelvzFKobPv5E7q9Glm5ONv76h3zmoCpkVnFSF+5rAZNTMPz0KNHjdsvoAwQAEhBk7V16onqqw0qiMk3guD0z5WNqx4rBFom7PWk1S075jmzJoog2AaUKFsG0xHf9/o0RZkOdnmd1r/mCl0zWalOBTASAZBANgEFnZpwtX+xPxcc0Yu2ZILjP4GTl+pv9sMZ/LIGs77aXZ6ovZuERiDcC6DeaXtLJ+oDJnrXFZ+oX+xJpAfmSFHc7XUDCnEPe9LiffZKiCvJVi2bxFn6R2IRJrAJ54Z6XpIUjDpAO4mntf98cRLEPywLQvQyayQ0rjINPPZBqb6hAk1gA8FBzAUnlrqXjyli0+kZcl/EBV6fEUi4wCLj7lmNiATduciTUAthCke/+HjBuvDGRk746a8YmEkYizaE3YmqvGp0Q/WciqCPY4Ea0joqyb0VTFMCKqIKL5RNRRRr8+MG+V2PbNPAu2f6oRzRa96sWgyySbOwxJUahg/U9iPo0xc4IDOBNZK4AnAPTK835vAK1Sj4EAHpbUbyx0loAsxFWCYnQ2OLWq+dZCrZOkM/tTsQze+YITFJk8Ol2szoFNMfi2bHfKKgk5HUC+oON+AJ7iKt4HsDMRNZbRdxxs2WoAxMPb3rNIYvbvk+3Y9zZN5QZ/Ys3vmWjPBGP4VDurasXBlnKXunwAewJIz4pZlXrNCDbNmn3iyffip+zLzFqVwZsCZUJVFTExgU0TDJ8QXfnLwjonMBENJKJyIiqvrJSvzx6wkw8+sStr9arRc2N/1pQEdC4+/lJvSUobsU0u/SdLCgXpMgCrATRL+7tp6rVtYOaRzFzGzGUlJSVaBmeauGqaa9f7IzY2xIJY83S+3+RPHL9PNYLj1r3+90dhMpkNXQZgPIBzUtFAXQCsZ2Yjoi82pLVnMjKmc8vGZLa4ceefBUlpZYwujxf9svIr+wqnT/soXr7JreMXSR6JODasSmSFgT4H4D0AbYhoFRFdQESDiGhQqslEACsAVAB4FMBgGf3Godygrnku4srdPiIYFaGCl+fY42A3gU8ZwM/MtE+G+fLn423NrfzavgnGkrXmt+ZqyjgIMw8o8D4DuFRGX6LcPXGJ6SF4zTVj5uP0Q5oVbpiGDTOhbCxcvR7t92xQ1GcWeCTM9+jbdiSz+cr9by7DY+ceYnQM1jmBVRMyI+1DNMFHFQ+8Vbye/7UWpfmnE2SQ7WPyEjO1r9NJnAGwlbjOLR+wcSsLACYtKj4U9FMLtxoAYNoyf5ygxfqZQsnS3AQDYAlvLSvuZjPj468UjUScYrd0QmU29Vz9UnFx55ssCVPMxtrvisvk/4eFwRK2kCgDsMbiGq1XvFDcD/Qui30Zi9eGbTbbKNbNMucz+4IlqvlzkbUBHn/HXl+GaUmIRBkAmyQgRFli8U32QY9S9ovxT9g8ay6Wv75hr6zHhAVmy4bKZOJCs+eSKANw3yR/ipbYzMQFduicyODtj6PvndsiAZ2LYmaboiJygWhcN8ZsAmSiDIDt2BoOqZJ3Kuz1ZQDF6UbZPsGwNdoqyWwwnDcSDIBFrPo2mo/CZl9GsdwrqdarKqJeExcYFiOs1VaiTpZskrO2kcQYANvUJrMRVRJi8pL4SpW6iFr4ep4FBceTQtRomI2b7c9mjhqia2uIsS0kxgC4kAAWVRLi5nH26ZpkYpu6ZyA6o2d9XriRYaJWOpsw336H8ScGNZcSYwD+8obd+7O+4dN2w/ebCu+df2hx2GSx3PqKXXLW2RAtpGQTJutHJMYAvG1x4pSPRFkBuCJP8Mz7nxVs40oyWxIDDWzn/jfNhdwmxgD4wi+W1BKVwSxHtomiJB69IVBBTCfvLU9Oha+fNvuTl6GKYAAsY97n+aMWfIoAuupFO8riJYm/TfFna+6L9fklIaYsdcMomyQRBiBukRITPFWgrm7c4jEmKDQD2xhmaNoptDVXse4HTSMRp1CS3o1jF2oaiTimohQTYQBmrnBjqwEAxszJX70pbvEYE7w635/i6Pn4waMiMM9/UNjfYQtXv5Rfevs7hxR2V35tJhJIVkWwXkS0jIgqiGhIlve7E9F6Ipqbetwso9+oPDLdDQedb1yXp86va87IxXnCiJc7NGsuxGOWy1n4yvC3zOhnCRsAIqoB4EEAvQG0BTCAiNpmafo2M3dIPW4X7bcYQgSQGbbkcVhv+tmdbTkgf3nEW1+xPy8jnSAJYR9j55pZLctYAXQCUMHMK5h5M4DnAfSTcNzE4kLWsigvzrY/2SidZ2bm3hr58DO35AaWWqwkWyy5VpLrvi+uZkBSkWEA9gSQ/mtelXotk25ENJ+IXiOidrkORkQDiaiciMorK/2pYlQM7+YI1fv4S/NFpGXhQjazr1z/cvatORcnHt/8uDnr628tNV9u0QV0OYHnAGjOzAcAeADA2FwNmXkkM5cxc1lJSYlwxy5+qXOpSr7gQIp+Jj6FrfrC8srsDsdvNma/mdrM8zl+E7Yrs9qCDAOwGkCztL+bpl77L8z8PTP/kHo+EUAtImoooe+CxKnrapoFq7MLpLnooJvigHBdVKIK3LmKKUekCLlu9F/nWBnYTMU6/St8GQZgFoBWRNSCiGoD6A9gfHoDItqDiCj1vFOqXy0piXdOsF/XxGfumGBv6cpiWbdh233l6R/5s01ZKAcloBYT///CBoCZtwC4DMAkAEsAjGbmRUQ0iIgGpZqdCmAhEc0DMAxAf9YUB7i2QLZgQC3ZZs2VG9zQAMrk75O3zaL1qfxlwCxOGgCgaluHmVszc0tmviv12ghmHpF6PpyZ2zHzgczchZnfldGvz2zIUKB0LW4+H5MWuVky8sXZ2ybpzXREzyiTacv8dZL6pJelmkRkArvI0i9+vR84x7FQw3QyjZdLKfq+8pAj6qVRyIwEerHcvWAJU3htAHKFiLnAXzKcW6NmuKMBlEmYkNlHpibQ0i/czQ0Y++GvYk5w07gwwYiK1wYgVzSNC2RuLUxc4Oa2CeDvjOw/W/wRszNZlESU2zOqg/38i7szjvKVercUvTYAN+RIeAnoxaekr/RkvFwJey7y4FR/toRc5vF39IZ6e20AVn0bkpBsYHNaMt6qbzcaHIk449I0W4ZOLFwoxmZcTJL0Hd0rfa8NgC/4tNXwpiOVs3IxPC3sc5nj0hxTPJRLcDXE2BTBAFjMuxVVKqbrvvfnS32XR4lhrlNd6tKHEOPPvq5aWb67PCj/FoO3BsCHL/Xj76wEAAzzoIzft6mIrHwS0QG9rEhpArlUBSwXb6XKP1452v0yozpX/N4agBkV7s8EJqd0dLIlILnGVM8Sj370qApYPqlrV7j1lapIIB+SwHT6Lr01AKMcFE7zmWvHzPdiVQZUVQdbXun+rLmaJ95daXoIgTR0ivJ5awCmLfNHpMsHfv6FnY7PTueFWZ/h3tf9kBvOlBwJmOfljMQ2ldTU1lMgFls9WNJWM3auvi+2Sp70SDVzRY7aAC7icua/KbxdAfjCvzTOBlQz9DW34+Z9xLV6xvmY7HiIsQm8NAA+zZr/9KL7UQ3VhBmafbhWzzgf14yZb3oIzuGlAfAhAigQCCSXLzTVMZFiAIioFxEtI6IKIhqS5X0iomGp9+cTUUcZ/ebigbfcj5sPBALJ5dmZevxMwgaAiGoAeBBAbwBtAQwgorYZzXoDaJV6DATwsGi/+Zi18luVhw8EAgGlDNMUCipjBdAJQAUzr2DmzQCeB9Avo00/AE9xFe8D2JmIGkvoOxAIBAIxkWEA9gSQLvi+KvVasW2ksOlnf4TTAoFAQCXWOYGJaCARlRNReWVl8clc29e07pQCgUCgaHTIWsi4W64G0Czt76ap14ptAwBg5pHMXMbMZSUlJUUPhoiK/kwgEAjYRo3t1N/LZBiAWQBaEVELIqoNoD+A8RltxgM4JxUN1AXAemZeK6HvQCAQCMRE2AAw8xYAlwGYBGAJgNHMvIiIBhHRoFSziQBWAKgA8CiAwaL95uPgvXZRefhAIBBQSq92e2jpR4oWEDNPRNVNPv21EWnPGcClMvqKwiVHtsSFT5Xr6i4QCASkMviollr68dJj2nnvXU0PIRAIBGLTrkkDLf14aQDq71DL9BCkcfNvMnPqAgF5NN+1rukhSOOS7npmzTrQ4QAGPDUAPnFG5+amhyCNm4Ixs45rerUxPQRpDDrCHwOgi2AALGeHWjVMD0Eax+y7u+khSOG8bqVotmsd08OQQp/2/iTkN6jrz8pfF94agB1r+3Pj9IXShjuaHoIUzu66F67rvZ/pYUhhO01bDYHo9GzXSFtf3hqAv55+oOkhBNLwaauhZUk9HNdW349UNRce1sL0EAJpnNdN3/Xw1gAcvJf7kUDVDrqj2hSfEW0b/Q/xx5cBADVr+PPTOa2sWeFGlnPVsa1ND0EaXVvupq0vf77FGZTU3970EIS5+6T9AQAXHb634ZGIs0vYn7WOffeoDwBok/rXZc47tBQAcH2ffc0OxDG8NQA+0LbJTgCAbvs0NDwScao1mq7tFX6gtuBT2GR16PcZnfcyPBK3CAbAYnbdsbbpIUjn9LKmpocgxJGt/7cd17bxTgZHIs5RnkRlpVPXo6g5HXhtAGoHaWjr2K2e21tz6c7sy47ex+BIxNnJo4TJalyPauqxn16j7PUd8o5+7UwPIQC/HHR7N6z33+c99vMnEsinqCaXuaS73kmF1wbg5I7ubjdkOrFPOkhJATUt9N7fn2SjOmn5JT6tMF0O0z27iz/7/rq3Ff35BmehlsOhepmyCed1KzUzEAnss3u9wo0CWsnUmmm6i7uaQOenIoCqudLhFWcdzQms7t4hPad3+1/rgbdr4rbDMZ2LDg+JR6Z5YMBBv/rbZcmRFhkZ5scf2MTQSNwjGABLyVy9+JR4dP6hbhqA7lkS8hrUcdORWuZR0aTMMrCZBiGQG6G7ChHtSkRvEtHHqX+zfquIaCURLSCiuUSktVJLQ8ejTnxkj512MD2EWAw6ctu4+ftOPcDASMTZ3dFr4DMdmu2svU/RaeUQAFOYuRWAKam/c3EUM3dg5jLBPovC5f1AH8iW+OVqqF6XvbdN0e/o0Uxadwhi4NcM6a0/SVLUAPQD8GTq+ZMAThQ8nnR6tddTW1Mm+++ZvRpQ5xbu6Rv1P8R9nZl8+LTCvDjLCsd2/uB4LkY6Lq4AGjHz2tTzLwDkCiZmAJOJaDYRDRTssyhczKb9fY4v9bUGZgii7OLg/7/v5PJbuOgXyOXwdbGSnglHfEEDQESTiWhhlke/9Hapwu+c4zCHMXMHAL0BXEpER+TpbyARlRNReWVlZTHn4g1H5lD/dF16IB2fksNc49YTst8cM52pLtCqUXYhu34dQiRQFAoaAGbuwcztszzGAfiSiBoDQOrfdTmOsTr17zoALwPolKe/kcxcxsxlJSXuyyDHYfua2WcCLofqZdL3ALeSwwZ0yi1n3aSBWw7VEw50N6kwKq5LjuhCdAtoPIBzU8/PBTAuswER7UhE9aufAzgOwELBfouijkc3Tl/Yu8St5LB8uQs3OrbdoKvgeCA6pmRFRA3AUADHEtHHAHqk/gYRNSGiiak2jQDMIKJ5AD4AMIGZXxfstyjuD9XBjOCT9HPpbrljy3u1cy/QIBcXhOpgRhh4hJmaHzVFPszMXwM4JsvrawD0ST1fAcDoHVhnhR1RWhWQTTjhwCYYP2+NptGIcerB7moxZZIvdNXVsNZsnNu1FKNmfGJ6GJEopF+0HQFbc3klLaOToQg/f9JL87BzXXciUe48sX3e9zN1T2zGh6psvnFggVDD5ru5own0uwIZ5bedENSAC5EIA+AS7XLkAFRzUHP3QvVycevxbu2d+4BPiZHbF1Bj7eewgq4uggGwjHrbC+3KOcWZjsj45osAquZoR6prdWyuP9lIFYXCVn0seCObxBgAn774LhBlW9wVue4LDist2CZX8p5t1A83Res4zaCvzI1foASuPNbdghcuMixDbthlWkYIWW3rkVz3xUeaiUgpBp9Cu68wuC2XGANwWKuGpodQkKg65tlUKW3j8H38SeKLkiGbK3nPRVyosHV7xHKvuzsQiGAyWCIxBsAFLowYgz3IgRlag7rRthpKHYo6cZ0jW0czyi5UB+veJprP5V4H5LpNboUGA2ARhUL0qvHJuXVziATSxh+OaWV6CNKIOmsuK3VPQVcnwQA4iE+JR7b/QAvlZaTz2zK7pa+TGAiRpKi6OCTKAJzS0Z/MVJvp2S66rontq5li6steaHmt42LUPms7EqHlOrcYXgEn6ir/4Rg3QvWiYPMiYHB3f/6fi6n5m0ua2EXuO83evXOfVjI9DetIJcoA2OzcuqHPfkW1H35GR0UjESeqLyNgL6ZvTPm4umdxIoPndStVMxAJNNm5jtH+E2UAbJbBHdC5cLZpOl2z1Kd1lVC8Qz3FTjBsrj1RrLjjWQ6EtZoiUQbAZop1VvlUatGUFG4hivFlVGPr3vnZXZN7E9yngMJukrHz2xpIFO2a5BfAM8XlxxSfofnnU/dXMBJxbJ7RJ5XOhiSg00mcAfh7/w6mh+A1V/f0R3KjdaPiZ479PCq3aPPeuQ8UqmegAyEDQESnEdEiItpKRGV52vUiomVEVEFEQ0T6FMVU6bV8HNg03gz4rC7F+Q10cFqZP6G2NWNs5/iUo3HpUfZFc90cs/zmTjvYlw9w8F7urwAWAjgZwPRcDYioBoAHAfQG0BbAACIyFvy6o4WJITfF/FKfX6Aghgl2rx+vQPouEaUjAsUTt9qUjQV9zow56bnnZHvDWk0iZACYeQkzLyvQrBOACmZewcybATwPoJ9Iv77RMWaRlygqla7wl9NC3WZVFJPNbDtxRfeObOOPOKFMdPgA9gTwedrfq1KvZYWIBhJRORGVV1ZWKh+cDfi0bRCXYyzbmrv7pPjOXNvKduYraJ8UgiREdgoaACKaTEQLszyUzOKZeSQzlzFzWUmJGqt9SXfYpDJ4AAAPr0lEQVT75ZRdxEb/SlyOixECWs1llu2d1y5QOjEfNm4D+YDIBEMmBb8ZzNyDmdtneYyL2MdqAOkqWU1TrxnDth+oCHEdyCq44lh/1CYb1ot/49vVoxyN4RYV9und3t7s5GLpu39j00MAoGcLaBaAVkTUgohqA+gPYLyGfnNikyNYNCw1rgNZBbbG8+umGNE11dQU3F48oKk9sh6ilbOu6GGu8lYmUetlqEY0DPQkIloFoCuACUQ0KfV6EyKaCADMvAXAZQAmAVgCYDQzLxIbtj8cf4CYDML+Fq0ARLm2V3EaL4HCiGpG1altTwJZa0GxvdMP8SdEWRaiUUAvM3NTZt6emRsxc8/U62uYuU9au4nM3JqZWzLzXaKD9glRB7BPpQhPNVgcO50/HSc+UzwiYvUt1fTYL1rlrCTQuIFZ4TUbSVwmcDXHtfXHYWkDnSQUdrHF4Tigk3iC3bUWZHkC8ZLZAmqxaSsqsd+O64tUR7QZG1RObznBHl+EKLsJOICr8ckfYoN8ig26ObLo38meynGJNQB7WVCMXFbB6sfOyanCoQ2fbniBX/MbQT+VDGQFO9gQAt5op3jZ8ipIrAGwIVJDVijYofs0lHIcGwiOYHmMOEtO0SAbVpjt95QzwbDBANhEYg2ADcgKRxVJ9LGN/oeYXR4PlniDMO1n6tXejlhzm7C9BrVu/LlzxOA3B4QfiAz2lzQ7A8wXupE5Q7zOIz/TjhaFg7qMbSvcRBuAq46zI1JDBk0amNtXlOXLsIH6EmeILRr6o8HzqEE/ky1ZszK48HC7FHwTbQBM/kBHnHWw1OP93WDK/r57iCXoBOQjImWRjS4Ga1Df1q+d1OOZLMRSy7KwXLtGkyCO2lduotAhEuLw4yLboW5D2KHryHIAV2NSsXY3yduCZ3ZObn3kTIIBMIRPGbyyOeFAM2GHfz5FvkLjxYYK3tuk4SOK7AlGgzrBEVxN4g1A3BJzgSrO6Sp/NmUqRLevgnj3wYaUZ1VEhtmkPOsij5wtd9tXBok3ACYKX8uQTcjGyR31FyT3SVpbRdEQn2ab9xmo2nZ1T38CNY60RB8qncQbABN7m0MVbDUAwI199a9mdleU1dhqd3/KXepG1URAVI0zDmd2FtdlysaTv+uk5Lj52KGWfdu+iTcAJthbUS1frwqRCMoYF0vH5v7smfs0a965rprv9BGt/MmeFyEYANihC+Qi9RUW1mmjObRUZVH6Zy7srOzY2Qiyx4XR7Wc6sYN5PaVsBAMAYOTZ5sXUXOTZi7qYHoI0VK3KAKBbS3Mx9LJRaSgzOcijVdm1ve3KAK5GtCLYaUS0iIi2ElHOuygRrSSiBUQ0l4jKRfpUgc7Z5sNnqt3aeFbjbNOnamQq0TnbbLST2poKp2gMNLj/dLX5IDplGWxdlYmuABYCOBnA9Ahtj2LmDsyc6On20YorNHX1aLap2lj6yD8vUDsB0GnMVGfqDzSUo2EToiUhlzDzMlmDSQKqE8BskLmWRa/2e2jp59mL1K+a7j1Fj16SiUgdV9Elc926kb0Rbbp8AAxgMhHNJqKB+RoS0UAiKiei8srKSk3DA55WPHPyjWEatId0GTNVeRnpnFZmR71jGeiIMqpjYchkXEade4jpIeSkoAEgoslEtDDLo18R/RzGzB0A9AZwKREdkashM49k5jJmLisp0Zc4cZiGsLDbTpArapWLUeeq32Xr1U7P7FwHOurm+rQyu+hw9VsnL1/aTXkfgLo8g3Sa7WpvlGHBbz4z92Dm9lke46J2wsyrU/+uA/AyAP1ZGBZwhoYvGwAcs5/6QiS6itDIVk31mX8N1nPT1HHt22jayrrGMn1+3Si/kkS0IxHVr34O4DhUOY8Th21SsC6guqqWjv3/alSHUB7UzJ+wSV0rJtVSHboMWVxEw0BPIqJVALoCmEBEk1KvNyGiialmjQDMIKJ5AD4AMIGZXxfpVxVvXJFzZyqQxj/O17enqVqqo1tLfRmhqlVOdW4zqdzO3HNnO0Mm4/C4xt9KHESjgF5m5qbMvD0zN2LmnqnX1zBzn9TzFcx8YOrRjpnvkjFwFajUn3lC8xfhlcsOU3bso9qoDWX1FZVbJ7orwg3opG47U+eqDABuUqgIbLsxC3sSaaicQXXXfNP0KUnr1d+rMWb7eCQ4N3pQV639qTRme+2mt1Lf+QYUgW0hGIAMjt43zG7zYUI0rb3EovPpjL5Y700TACb8QY0xa7qLvZEmtqNqm/HSo1oqOa5MggHI4P8UpJ+b0oJRsef8tGZhM5WYUE9t18SfldnkK4+UfsxbjjdToKlWDflG4Mpj7VdlDQYggwZ15UcFPGwolPHPCrJP69ZWpwCaj98f7U/hGdkMP0N9Ul42VGyhmSjQBADvDjlG+jF1ZRqLEAyABkxVhapT259syit6tJZ6PN1O+XTuPVWuYe67f2OpxzOJqYS5kvpyRfR0O+XjEgxAFqZcJX9pa4qddpA3Y1cZWVQI2fu0JsvznV7WTOrxTGYZj5RY59an7PIpV3U3PYRIBAOQhZYSteGn/am7tGPF4b3r5C1tTUcWnXawPD0dX6QZLjasaHmsxES9v/VXK/9cCJnRZq6svoMBUEypYknbQuyosGqXbu4+WU4t5cfOMa9IftWxcra0TEsZyDSkpmvmyoo2q+fQby4YgBzMuqGH6SFIo12TnYSP8drlh0sYiRiypDR6KJaXiMJlkpzaNjgaZZS87KFBv0oX0685yvQQIhMMQA5kOIWW3tFLwkjEeVFCktB+jcWNiAxciK2OgoyZsy1CeYfuIy6n8YhEX4IIb0u4eZsIL45LMAB56NxCTCfe9JK2GtHQzaa72JPO/qfjxGKr3xlytKSRiDPu0kOFPq+rYI4ObFjJAOLSzWd10aP4K4tgAPIgUvT8hj77SRyJOCIzrKmGHdnpiM6cbdJmOdAj9c4lt8df7Y4VNISyuaR7/FXmHf3aSxyJeoIByIPIrOQiy+qN9hQIsbNNxnrezcfF+pxtRhkAzu26V6zPfXRnb8kjEUMk6qWDZYZQpFi8a9Fldv2yLWTOTccW/ZlDSndRMBJx4ixP37vOni2TauJma9tmlAHgtpgzRl0FeYrh9T8WHyhwY1/7jDIA1Iwx+bPNKEfBvm+RZcRx6JgQGYvCnScWH0bZuIE9WybpzLy+uPyGCw9roWgk4hxTpABh3BWQavbdo/hAgQs1lJeMQ5ybuY1GuRCiBWHuI6KlRDSfiF4moqxrOSLqRUTLiKiCiIaI9GmCxbf3jNx2cPeWVi8Di9E+t/VGAwCNdiou1f5GhZrvoow6rzhZChV6VbIoJorm4TM7KhyJGNttRzisiOgmF2f/gPgK4E0A7Zn5AAAfAbguswER1QDwIKoKwrcFMICI7P01ZqFu7ZqRI4JMJ+YU4oKIM+HWjepZfaMBgOV394nUToVqpWyiFkH55J5o52yKYqJoeluuYRRV+fasLs2dnP0D4hXB3mDmLak/3weQLVe/E4CKVGWwzQCeB9BPpF8TvBBhW6eYlYJJosxW3rjC/ptmje0IzxWI1DqidYkThV+6tWxYMEJp8pVHWL26rGbl0L4F20Q13qaJsgqOs7VqCzLN1u8AvJbl9T0BfJ7296rUa86Rb/b1zpCjjUklF0vtmtth9o25M51tn2Wm07Xlbvi/32Yvtt5817p46nedNI8oPvlyFB4/rwz77G53gfF0Ku7KPclYcOtx1sT9F6JB3VqY9MfctcKjGDubKWgAiGgyES3M8uiX1uYGAFsAPCM6ICIaSETlRFReWVkpejipEBFWDu2Lq3v+Lxmpbu0aqLirt1Xx5VHYrd72WHF3H+xd8j+toit6tMbKoX2dmGWmc9JBTbHwtl+vvib98QinUvKrWTm07zYJYotv74mj93VLKqFmje2wcmjfX+n7d2y+Mz65pw/q72D31mImbfaov82q+Y4T2zt/8wcAYmaxAxCdB+BiAMcw88Ys73cFcGt1wXgiug4AmPmeQscuKyvj8vJyofEFAoFAkiCi2cwcSfFQNAqoF4BrAJyQ7eafYhaAVkTUgohqA+gPYLxIv4FAIBAQR9QHMBxAfQBvEtFcIhoBAETUhIgmAkDKSXwZgEkAlgAYzcyLBPsNBAKBgCBCXktmzqppy8xrAPRJ+3sigIkifQUCgUBALm4GrwYCgUBAmGAAAoFAIKEEAxAIBAIJJRiAQCAQSCjBAAQCgUBCEU4EUwkRVQL4NObHGwL4SuJwTOLLufhyHkA4Fxvx5TwAsXPZi5lLojS02gCIQETlUbPhbMeXc/HlPIBwLjbiy3kA+s4lbAEFAoFAQgkGIBAIBBKKzwZgpOkBSMSXc/HlPIBwLjbiy3kAms7FWx9AIBAIBPLj8wogEAgEAnlw2gAUKjZPVQxLvT+fiKytQh3hXLoT0fqU6upcIrrZxDgLQUSPE9E6IlqY432Xrkmhc3HlmjQjoqlEtJiIFhHR5VnaOHFdIp6LK9dlByL6gIjmpc7ltixt1F4XZnbyAaAGgOUA9gZQG8A8AG0z2vRBVZlKAtAFwEzT4xY4l+4AXjU91gjncgSAjgAW5njfiWsS8VxcuSaNAXRMPa8P4COHfytRzsWV60IA6qWe1wIwE0AXndfF5RVAlGLz/QA8xVW8D2BnImqse6ARiHIuTsDM0wF8k6eJK9ckyrk4ATOvZeY5qecbUFWXI7MutxPXJeK5OEHq//qH1J+1Uo9Mp6zS6+KyAYhSbN6VgvRRx9kttQx8jYja6RmadFy5JlFx6poQUSmAg1A120zHueuS51wAR64LEdUgorkA1gF4k5m1XhehgjABrcwB0JyZfyCiPgDGAmhleExJx6lrQkT1AIwB8Edm/t70eEQocC7OXBdm/gVAByLaGcDLRNSembP6nFTg8gpgNYBmaX83Tb1WbBsbKDhOZv6+ernIVRXWahFRQ31DlIYr16QgLl0TIqqFqhvmM8z8ryxNnLkuhc7FpetSDTN/B2AqgF4Zbym9Li4bgCjF5scDOCflSe8CYD0zr9U90AgUPBci2oOIKPW8E6qu3dfaRyqOK9ekIK5ck9QYRwFYwsz352jmxHWJci4OXZeS1MwfRFQHwLEAlmY0U3pdnN0CYuYtRFRdbL4GgMeZeRERDUq9PwJVdYj7AKgAsBHA+abGm4+I53IqgEuIaAuAnwD051SYgE0Q0XOoisJoSESrANyCKueWU9cEiHQuTlwTAIcCOBvAgtR+MwBcD6A54Nx1iXIurlyXxgCeJKIaqDJSo5n5VZ33sJAJHAgEAgnF5S2gQCAQCAgQDEAgEAgklGAAAoFAIKEEAxAIBAIJJRiAQCAQSCjBAAQCgUBCCQYgEAgEEkowAIFAIJBQ/h9JsW2bxP4mbgAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plot(t,sine2_sd);" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## A chirp\n", "\n", "Bats use a chirp, or tones that rapidly increase in frequency, to use as sonar.\n", "\n", "We can use a chirp to explore the range of frequencies that we can hear.\n", "\n", "Humans (with very good ears) can hear from 20 Hz to 20,000 Hz. " ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "slideshow": { "slide_type": "slide" } }, "outputs": [ { "data": { "text/html": [ "\n", " \n", " " ], "text/plain": [ "" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" } ], "source": [ "f = 20000/(2*Len)\n", "chirp_sd = sin(2*pi*f*t**2) ## t**2 creates the increasing freq's.\n", "\n", "Audio(data=chirp_sd, rate=Fs)" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Count the pulses\n", "\n", "We test how well you hear, by testing various frequencies.\n", "\n" ] }, { "cell_type": "code", "execution_count": 73, "metadata": { "slideshow": { "slide_type": "fragment" } }, "outputs": [ { "data": { "text/html": [ "\n", " \n", " " ], "text/plain": [ "" ] }, "execution_count": 73, "metadata": {}, "output_type": "execute_result" } ], "source": [ "freq = 10000\n", "\n", "n = random.randint(1,5)\n", "test_sd = sin(2*pi*freq*t)*sin(pi*n*t/Len)\n", "Audio(data=test_sd,rate=Fs,autoplay=True)" ] }, { "cell_type": "code", "execution_count": 74, "metadata": { "slideshow": { "slide_type": "fragment" } }, "outputs": [ { "data": { "text/plain": [ "4" ] }, "execution_count": 74, "metadata": {}, "output_type": "execute_result" } ], "source": [ "n" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## More on chirps.\n", "\n", "We can vary the parameters, to vary the sounds" ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "slideshow": { "slide_type": "fragment" } }, "outputs": [ { "data": { "text/html": [ "\n", " \n", " " ], "text/plain": [ "" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Fs = 44100\n", "Len = 3\n", "t = linspace(0,Len,Fs*Len)\n", "fmax = 15000\n", "chirp2_sd = sin(pi*(fmax/Len)*t**2) \n", "\n", "Audio(data=chirp2_sd, rate=Fs,autoplay=True)" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "
Thank you for your attention.
\n", "\n", "![alt text](https://github.com/callysto/callysto-sample-notebooks/blob/master/notebooks/images/Callysto_Notebook-Banners_Bottom_06.06.18.jpg?raw=true)" ] } ], "metadata": { "celltoolbar": "Slideshow", "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.5" } }, "nbformat": 4, "nbformat_minor": 2 }