Skip to content

Instantly share code, notes, and snippets.

@duartenina
Created October 16, 2019 21:44
Show Gist options
  • Save duartenina/6d31fd6062fedfbe359c292dc15e1a91 to your computer and use it in GitHub Desktop.
Save duartenina/6d31fd6062fedfbe359c292dc15e1a91 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [],
"source": [
"def func(t):\n",
" return np.sin(10*t) * np.sin(5*t) * np.exp(-t/5)"
]
},
{
"cell_type": "code",
"execution_count": 23,
"metadata": {},
"outputs": [],
"source": [
"ts = np.linspace(0, 10, 1000)\n",
"dt = ts[1] - ts[0]\n",
"\n",
"fs = func(ts)"
]
},
{
"cell_type": "code",
"execution_count": 24,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7f07fa133748>]"
]
},
"execution_count": 24,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD7CAYAAABpJS8eAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOy9eZgk513n+X0zjryzso6s6jr6PiS1bqktW/iSbQS2Acs7eHZtGGCYWbzaQSzs8Dzg8ezDXPDALMdjBhu8gp0ZFsaP8WB7ECCwsbEtfEiWZN1q9X3W1VWVlZV3RkbGu39EvJGRkW9EZne90S1Vvp/n0dNV1aGMqKyub/zi+7sIpRQSiUQi2fnEbvQFSCQSieT6IAVfIpFIRgQp+BKJRDIiSMGXSCSSEUEKvkQikYwIUvAlEolkRBAi+ISQ9xJCThBCThNCPsb5+zFCyF8SQl4ghLxCCPlpEeeVSCQSyfCQ7dbhE0IUACcBPAjgMoCnAXyEUvqq55iPAxijlP4yIaQA4ASAXZRSY1snl0gkEsnQqAJe4z4ApymlZwGAEPJZAA8BeNVzDAWQJYQQABkARQDmoBeempqi+/btE3CJEolEMho8++yz65TSAu/vRAj+PIBLns8vA3iz75hPAngMwBKALID/hVJqDXrhffv24ZlnnhFwiRKJRDIaEEIuBP2dCA+fcL7m94l+EMDzAOYA3AXgk4SQHPfFCPkoIeQZQsgza2trAi5PIpFIJIAYwb8MYLfn8wXYkbyXnwbwBWpzGsA5ADfzXoxS+iil9Bil9FihwH0qkUgkEsk1IELwnwZwmBCynxCiA/gwbPvGy0UA7wEAQsgMgJsAnBVwbolEIpEMybY9fEqpSQh5BMCXACgA/jOl9BVCyMPO338awH8A8F8JIS/BtoB+mVK6vt1zSyQSiWR4RCRtQSl9HMDjvq992vPxEoAfEHEuiUQikVwbstNWIpFIRgQp+BKJRDIijIzgmx0Ln3nqIqqtgf1eEolEsiMZGcF/7IUlfPyLL+F3v3LyRl+KRCKR3BBGRvBfvLwFADi7VrvBVyKRSCQ3hpER/PMbttBf3mzc4CuRSCSSG8PICP5yqQkAWCxJwZdIJKPJyAj+Rq0FAKi2TLTMzg2+GolEIrn+jITgdyyKYs3ARFoHAJTq7Rt8RRKJRHL9GQnBL9UNWBQ4NJ0BABRrcu+KRCIZPUZC8MtNu/Z+70QKALBZl4IvkUhGj5EQ/ErTtnBmxxIAgHpLevgSiWT0GAnBrzoR/q6xJACgZshuW4lEMnqMhOCXXcGPA4AcryCRSEaSkRB8ZunM5KSlI5FIRpcREXwnwncEX0b4EolkFBkpwc8lNSQ1BTUp+BKJZAQZCcGvttpIago0JYZ0XEXNkJaORCIZPUZC8CtNE5mEvc0xHY8+wm+2O6jLSiCJRPI6Q4jgE0LeSwg5QQg5TQj5WMAxDxBCnieEvEII+YaI8w5LpWkiywRfVyMX/I9/8SXc+m++hE3Z0SuRSF5HbFvwCSEKgE8BeB+AowA+Qgg56jsmD+D3AXyAUnorgH+83fNeDTXDRCZuC34mrkZah08pxRe+twhKgWcubEZ2HolEIrlaRET49wE4TSk9Syk1AHwWwEO+Y34MwBcopRcBgFJ6RcB5h6ZudJDUFABAKq6gFmFZ5oYnqj+xUo7sPBKJRHK1iBD8eQCXPJ9fdr7m5QiAcULI1wkhzxJCflLAeYem2e4gqduCn45Ha+lcKtbdj5e2mpGdRyKRSK4WVcBrEM7XKOc89wJ4D4AkgO8QQp6klPYtmCWEfBTARwFgz549Ai4PaBgdJPO24Cc1Bc12dBH+WqXlfrwkl61IJJLXESIi/MsAdns+XwCwxDnmbymlNUrpOoAnANzJezFK6aOU0mOU0mOFQkHA5QENT4Sf0GJompaQ1+VRathdvYemM9ioyqStRCJ5/SBC8J8GcJgQsp8QogP4MIDHfMf8BYC3E0JUQkgKwJsBHBdw7qFoeDz8hBpthL/lLFfZN5keagxzs93BRx59Ep/+xpnIrkkikUgAAYJPKTUBPALgS7BF/HOU0lcIIQ8TQh52jjkO4G8BvAjguwD+iFL68nbPPSyNdlfw41oMzXYHlPpdJzGUGgaUGMHCeHKozVpPny/iO2c38Bt/81pk1ySRSCSAGA8flNLHATzu+9qnfZ//JoDfFHG+q4FS2mvpqAosCrQ7FLrKSz9sj1K9jXxSw3hKR7Vlot2xoCnB99VXlrqVPMWagclMXPg1SSQSCTACnbYt0wKl8Hj4ivP1aGydUqONsaSG8bRmfz4gyr+w0a3qWZRJXolEEiE7XvAbztwc18PX7G+52Y4mcVtutDGW0pBPsYXp4T7+YqkBNWY/aSxuSsGXSCTRsfMFv90r+HHnz6gSt11Lx4nwG+ERfrHWwpGZLABgcwjPXyKRSK6V0RH862bpGBhLasgn7Qh/0DydYtXAgULaPlYuV5dIJBGy8wXfZ+nE1WgtnVqrg0xCRd6J8LdCInxKKTZqBubySSS02ED7RyKRSLbDzhf8gAg/KkunbphI6yrSzrC2sDEOdaODlmlhIq1jPKVLS0cikUTKzhd8f9LWifBbEXTbdiyKZttCUleQjtvnC1u2UnTsnom0jnxKH6puXyKRSK6VHS/4deP6RfjsaSKtq4irCjSFhO7PZZM1J9M6xlOatHQkEkmk7HjBb7b9ZZlM8MVH+HVH3IedzFms2YPW7Ahfk0lbiUQSKTte8Ps9fJa0FR/hs6cJZuekdTU0wvcuV78elk7Honh5cUuOcJBIRpSdL/h9VTpOhB9BWSbbpJXUPNu1hhD8bEJFNqGiEvHqxd/96in88O99E3/14nKk55FIJK9Pdr7gB0b44i2dhj/CH7Bdi0X/2biGjK7CMC20O9GNbv726XUAwNdOXNeFYxKJ5HXCzhd8o4MYAXRngFmUjVfM0kl5PPxwS6cNJUaQ0GJIOWWc9ZCqnu1ycrUCAHh1Sa5elEhGkZ0v+M5oZELseTVM+I0IyjLrjqWT0oezdKpNe7k6IQRp5yYR1frFrUYbZcdCkpu4JJLRZDQEX+9OgY7FCNQYiUjw+yP8UA+/ZQs+AE+EH43gM5G/bT6HctMc6jwnVip49sJmJNcjkUiuPztf8I0Oknrvt6mrsUgEv+ar+c8MsHSqTRPZhC303Qg/GkuHNXndNjcGAFgZsGDd7Fj4wU88gR/9g28PnAckkUjeGIyG4Du+PUNTYjAiSI42nKg57TxRpOMKakbwdq1qyyP4bBRDRBE+E/yjczkAwEo5XPBPOH4/ABnlSyQ7hJ0v+O1+wY8qwq/7SkDTcRUdiwaOcag0u5YOu0lEFeGzpq5bZh3BHxDhv7bcFXyv+EskkjcuQgSfEPJeQsgJQshpQsjHQo57EyGkQwj5kIjzDoNt6fgEP6II3zAtaApBzFlowsQ8yNaptkxkEvZUzZRTyhmVh88ifDZ7/0qlFXr8hY0alBhBLqHiUrEeeqxEInljsG3BJ4QoAD4F4H0AjgL4CCHkaMBx/xH2svPrBi/Cj0cU4bdMy23sArxRO1/EKz0efsQRfs1ALqEil1ChK7GBYxwubTawK5fA/kIGl+UmLolkRyAiwr8PwGlK6VlKqQHgswAe4hz3cwA+D+C6dv14F5gzorJ0DNOCrnbf0vTACL+NbLzr9wPRlWUW621MpHUQQpBPadgaMMZhrdLCTC6OQkZ3h7xJJJI3NiIEfx7AJc/nl52vuRBC5gH8TwA+LeB8V4WdtFV7vqar0Vg6LbPjLlgBupYOL2pvdyw021a3LFOPNmm7WTMwnra3cNmz98NFfK3SQiEbx0Rad4e8SSSSNzYiBJ9wvuYvS/kEgF+mlA70KwghHyWEPEMIeWZtbW3bF9dod9xxCgxduV4RfnDUXnWaoDKOpcM6bqOK8DfrBvJJO18wltIGDmpbq9qCP5mJo1gz5MA1iWQHIELwLwPY7fl8AcCS75hjAD5LCDkP4EMAfp8Q8kHei1FKH6WUHqOUHisUCtu+uGa74zZCMaKydGwPv9/S4Y1LYDYPi/ABe+xDVKsX7RJQW/DHBwi+YVoo1gwUMglMpnW0O9Tt0pVIJG9cRAj+0wAOE0L2E0J0AB8G8Jj3AErpfkrpPkrpPgB/DuBfUEr/h4Bzh0IpDS7LjKhKxxvhsxsNz6bpTsrU3K8lNSWy1Yu1luk+TeST4ZbOhmPhMEsH6Fb5REHL7OA/ffUUzq5VIzuHRCIRIPiUUhPAI7Crb44D+Byl9BVCyMOEkIe3+/rboWVaoBRI8Moyr2OVTp1n6QRF+BFcF9Bb859Payg12oE2zZpTsjmV0T2CH52P/xfPLeF3/u4kHvnMc5GdQyKRAOrgQwZDKX0cwOO+r3ETtJTSfyrinMPAouWE6uu0jbJKR+m3dHh7bVnUz3x+wC4XjSLCb3cstMxugjif1GGYFhrtjpss9sLsnom07t7A1qvRRfjfObsBwG7warY77kRTiUQilh3daeufhc+IK7FIlpi3Or2Wjq7GoCmE20xVb7HZ+X4PP4LFLC12c7HPNeYkb8sNvi9fbtqCn0tq7rGVCD3848v2uOaORd0RzhKJRDw7W/B9ow4YuhqLZNFIq91blgnY5Za8skwmwt6EclQefnfRik/wm/zE7Vaj7R7HGsMqAceKYKnUwH37JgAAF4fo6v3aa1fw0//lu5HmFSSSncjOFnxm6VyvpK0vwgfsKZi8Uktm6fR6+LFIqnSqvgg/l7T/ZMLuh0X+uYTmJnqjivBrLRPlpolj+8YBAItDdPX+0udfxNdOrOHPnr408FiJRNJlRwt+M8DSiSxp2+5N2gL2nHteWWY3wr9+lg4T71yCWTrBEb6m2H0BmhJDQotFFuEvb9kCf2Qmi1xCHTjGYbXcdJPKz5wvRnJNEslOZUcLfsOwRf16TcsMjPA5Hn7N6EBXYj3HJzTFfSoRCYvOM06CODfA0ik32xhLau6WsGxCC53rvx2WSvbUztmxBBbGU1gcsI3rtRXb45/JxXHqiizjlEiuhh0t+G6EzxF806KwLLHdo0Eefj3Aw0/Fe68rKkuH5RAycafTdkDSdqvRdp8CACCbUCNrvGJjmmfHkpgfTw60dJjH/9ZDU1jeaqAj+GcokexkdrTgd6t0+jdeARDu4xsdq0/w7SUonAi/1XHr9BkJTUErkqRt270WAG4iNtjDb7tPAfbxWmQePhvMNpXVMZWJu01fQVwq1hFXY7h7zzjaHYorlfC5/hKJpMtICL7fV2e18iJLMymlfaMVACfCD/Dw030RvoKmKV7wu7t2baHXlBhSuhLo4fsFP5dQI/PwS3UDuhpDUlMwldFRrBmhUfuFjRp2T6SwezwJYLgkr0QisdnRgh+YtGURvkDBNy0KStHv4ceDq3S8NfiA3SDW7lCYgp88GhxrK5fQQjx807V9AGc3b0QR/mbdwHjKzhdMpnVY1L4JBHGx2MDu8SQWHMGPelZ/pdl2y3slkjc6O1rwA+vwFfGWDnta8At+aITfZ+nY/6/o8QosL+B9+hhLagM8/O61ZRNqZJZOqd5GPmmPb5jIxAGEz+1ZqzSxayyBmVwCACK1dLbqbbzrt76OD37qWzJXINkR7GzBD6nDB4C2QGFlTwt++4hV6fjn1tSNTp+lw55ERJdmsmQyW70I2LX4PA+fUopyo+2L8KOr0inV28in7HNNOXN7gsY4dCzqTPGMIxO3N3cVa9E1hP3d8VWsVw2cWK3IRe6SHcGOF3xdjUGJ9Y7sjyJpawRF+HEVlKKv+qZmcCJ8NRrB5239CrJ06kYHpkV7PPx0XEGdc9MSgW3p2EI/6UT4QYnbjVoLFgWmsnEQQjCe1iId6vbcxa7Iv3i5NPB4SinOrFWFV39JJKLY0YLfNPpHIwMeS0dghN8yWYK4vw4f6B+RXGt1+soy48zSEVya2Wx3+gbI5ZJ8wWdf80b4SV2BRcUmuRmlRhvjaftck5nwUczrFfvrBefGMJGORzpe4eRqBfftm0AhG8erzryfMP7s6Ut4z29/A5/46qnIrkki2Q47W/DbFl/wVfFVOoERvjsiuTdqt6t0+ssygSgifKsvwg/y8L1jFRjumGfByUtKKUp1A2PJ7upFQoItnfWqM7Y5awv+ZFqPVPAvbzaweyKFW2ZzeG158FC3//H8IgDgs9+9KDeESV6X7GjB51kZQDRVOq0gDz/eH+Gbzrhiv6WTjEjwm5yGsFxCRbnZ7rMfvIPT3OtiTymCffxm20K7Q91zKTGC8ZSOjSrfpnEF34nwxyMUfMO0sFJuYmE8iSPTGZxZq4aKeK1l4tkLm0jpCq5UWpFXD0kk18KOF3zebPV4BB5+WJUOgJ4RyWw+fnCEL97S6fPwkxooBao+q4nV5rMBa0A3whc99qE7QK57bWFRO5uhU/BE+BsRCf7yVgOUAgvjdgdwy7RCz3VitYJ2h+In7t8LADgtt3dJXofsaMFvchaYA4Cu2AJzXTx8d5F5Vyzd+fR6/2gFIJoIv8/DDxigxovwUxFF+Mzm8g6Qm0jr2AixdBJazH3fxlM6Kk0zklHXLEJfGE9hPj+4yeuMM9fnPTfPAAAurNeEX5NEsl12tOA3gpK2EVg6Az18TyRdN3rHFTNYhC86kuZW6QTM03GXnyT6BV90A1J3bHP32ibSOooBjVfrVQNTmbg71I09hUTRI9Cd8WMPdQMQOtjt7HoNmkJw95480rqC8xuD5/pvh1rLxIUNeVORXB1CBJ8Q8l5CyAlCyGlCyMc4f//jhJAXnf++TQi5U8R5B8FbYA4AmmILhtERJ2CupaP4q3ScNYeeCL/qbrvy1eFH5uFbfU86QTPx2edZT+MVu2nxVjVuB3bj80b442k9sNN2vdpy/Xuge1OKYuwDKw2dysYx73b1Bov42bUq9kykoCkx7J1MD7XI5VrpWBQf/NS38MBvfR0vXBpcLiqRMLYt+IQQBcCnALwPwFEAHyGEHPUddg7AOymldwD4DwAe3e55h6HR7vQtMAeijfD9wspKL3sifNfS6Y3w4xF12jaMEEvHJ5blhr3sXPXcuHjfgwh4uYyJlI7Nen8yGQA2qgamnNJNAJ5tXOIj/PWqgbhq20djSQ1pXcHyVnBX76ViA3sm7CeB3RPJ0JvDdnnu4iZOXamCUuCzcgmM5CoQEeHfB+A0pfQspdQA8FkAD3kPoJR+m1LKulieBLAg4LwDCazDd5O24krnuhG+v9O2PzoelLQVPTGzZfbf+LojkvsjfK9/D3QtHdFlmXWOpTOe1tGxKFfEizUDE2mv4IcvctkO7GmC2UfTuQSuVIKbvFbLTewas58EprPhx26X5y7aUf3de/J4eoglMIulBh74za/hl/78hciuSfLGQITgzwPwhhmXna8F8c8B/I2A8w4kyNJhtkskoxV8EX5Ci4GQrrgB/H22QHSdts22xW28AtA3577cbPfYOQCQ0qKpw69ynnTGnTELfh+fUnuswrhH8JktFcWs/o2q4TaCAXZl0FqZL+Its4ONmoFdznyfQjaOUr3tJvJFc3y5jJlcHA8cmcaZterAsRd/8p0LOL9Rx+eeuYwzsnpopBEh+ITzNW7oTAh5F2zB/+XAFyPko4SQZwghz6ytrW3rwpqchiPAM0tH6GgF+5fb7+ETQpDWVV+E37/PFrBzC4SI72i1k7a915WNqyCE7+HnfBE+ew/roqt0OE86TND9pZk1owOjY2HSK/gRe/jec01n44GD2q44N4JdY3H3WCC4gWy7vLpcxi2zOdyxewyUAi9d3go9/usnrriVRk+e3YjkmiRvDEQI/mUAuz2fLwBY8h9ECLkDwB8BeIhSGvivjlL6KKX0GKX0WKFQuOaLopQG1uFrkYxW4Ef4gB3J99ThswjfJ/iEECRURajgtzsWOhbti/BjMYJMXO2zQ/yD0wD7BqkpBPWI6vC9TzoTzlydTZ/gFx3xnEh3k7bsSSSqCN+bIJ7OJtw+AD+rZftGwCZ4Tufs/+9KWfwkT8uiOLtew5GZLG7elQUAnF0Pjto3qi28tlLBP3nLXoynNLx4KfzmINnZiBD8pwEcJoTsJ4ToAD4M4DHvAYSQPQC+AOAnKKUnBZxzIK2AJCoAqM4wNbERPr9KB7Aj2N46fKf+nNcUpsWEWjqNgJ0AAH+AWtm33pBhr2oUnLRtmVBipKd3gXn0mz5Lh1k8E+neOf1ANBG+3z4qZOOoGR1uL8KKI+y7xhxLJ8NGN4v38ddrLRimhYXxJGayCcTVGC6ElICyHcC3z4/h9oU8XhhiCNx2qLZMfOv0eiS9EZLts23Bp5SaAB4B8CUAxwF8jlL6CiHkYULIw85hvwJgEsDvE0KeJ4Q8s93zDiJoFj5gR9K6GhOetFVipKe6heFfglJrmUjpSs+4YkZcjaElsNOW3TzinPeBN0/Hv/yEYT+lCI7wWx2kdcVNjAJdS6dP8J0ySW+Eryp2FU3QXP9rpWV20DKtnp0AzKbhibhbs5+zbRPWCRz0RLAdWPPX3FgSsRjB3skUzoc0eR13hr7dPJvFTTMZnFuvRTrN8xc++xx+/I+ewn/8m9ciO4fk2hFSh08pfZxSeoRSepBS+mvO1z5NKf208/H/Sikdp5Te5fx3TMR5w+BtefKiKzGxZZkdixvdA3Z0XPONVvBX6DDiqiI02dc07O+R9z7kkr2WjtmxUG2ZPWMVGFEIfp2z9SutK9w596z7llk+DHvfrtgIn1UIeXMZYTbNylYTCS3mvm9s+qfflhLBUsk+P+sN2DuZHhjhT2XsfcF7JlJomVZkFUTLWw185fgVAMB/f/byUJvb/urFJXzmKTls7nqxYzttmWfOszIA25cW+djZcmbv80j7xNLedsW/roQWE+rhsx25PGvLb+kwoeNH+Kr4OvxWp69Sic2594sli/gnMn7BF7+Ni71etifCD7ZpVspN7Mol3CeVuKogpSvYrIu3mhZLtrgzwd83mcKFYnDUfm69hoOFDABgz2QaACJrCvvuObtE9P949yFsNdp4+nz40phXl8p45DPP4eNffMm9UUiiZccKPutm9VfCMDSFCI/w/XN0GKm42mPp8CJbRlxVxHr4IdZWLqn1RPgl5+NxXxQN2BG+6E5b3l5fdn5/WeZGzV527r9RBs313w7uADlPLmM6xKZZLTfdhC1jPBXcMbwdlreaSOuKe217J9NotoOj9gsbdeydtBvCWGNYVCMZXri0hbgawz97234Q0r0BBPEnT15A3Flg//lnL0dyTZJedqzguwPKAgVfdIRvDR3hVzn7bBlxVXCEH7DmEXA8fE90zARqLMX38EXP0qm3Otz3gSeWxaqBiZTe4/cD0UT47jwhz5NOPqVBU0hwhD/WK/j5lNaXhxDBZs1wN4MBcMX8PEfEay0T69UW9jqR/Xw+CUKASxGNbn5psYRb53LIp3QcLGTw0mJ4gviJk2t4103TeN/tu/D0+aK0da4DO1bwWTNKUIRvJ20FCn5YhK/3RviVJt8nB2xhFin43b2+fEun2jJdr7XkWBB5nqUTV/u2dm0XO8LvvxFNcEYkb9Z7u2wZUXr4XkuHEIJCpr8Wn1KK1XKrT/DHnRERoinW2z3VQ7udwW6XODYN8/bZTUFXY5hMxyMpFwWAk6tV3DybAwDcMT+GF0P6Ay4V61gsNXD/wUnct28CGzUDZyOcMPqt0+t40699Bb/6V69Gdo43AjtW8AdF+KKTtnaEz/fl7Z2wHTeCsbtZ+0UVsCN8kZYOm63Pi/D90yZLDVtk8xxLJ6mJj/DtaiVOhJ/W+sRyo8YXfHuRi+AIn2PpAEAh11+Lv1lvwzAtt8uWETYEbjsUay1MeJ7A5pyonbdw5WLRFtB9ToQPsAYy8UnbUt3AVqON/c65js7lcKXSCtxt8Jwz9O3YvnHctScPAHh5MboegV/96+NYq7TwR98851YuhfGdMxv4rS+diORneCMZAcG/PknbUA9fV2Fa1H2iqDTNnpI/L3HRSdsQS8c/QC00wteVCBag8KuVJhxLp+NJRG4GCH4moaJ6HZK2gCOWvvEKrCSzT/BT/TctEWzW2j2lqboaw2wugUucYW0swt/jRPgAMJOLu41iImHjoPdN2YJ/cNpOFJ++wm8KO7FShhIjODSdwf6pNGKku1NANJeKdRxfLuPn3n0Iaoy4qyiDWNlq4qf/63fxya+dxr997JVIrulGsWMFf3DSVqylY5jhVTqA7VlTSlFuhEX4gssyQ8pTx3wz8Zng+0crsP8/iuFpvGql8bQOi/YORQuO8DUYHUvoe1ZuthEj/dNMp7NxrPnWL66U7ch6ps/D11FutntuWiKwB8j1/nwWxlO4XOyP8C8U6xhPaT1PKjMDhsBdKywRvM+5uRwqDBL8Kg5MpRFXFcRVBXsmUpFtCfv6CbsC6B/ds4D7D07iK6+uhh7/589eQrNt4d03T+Pxl1eEFwXcSHas4NcNE4QE1+FrCkHbFNt4FRThZ9yZLyZqRgcWRYiHH03jFd/Ssa+LzdPZarSRS6hQOA1hSV2BYVrCBMyyKGpGp2+8BNCtEmKVOoZpodI0+RG+8/+LjPLZDdnfGDedTaBYM3qswJUtZ44OJ8KntH9W0XZoGB002p0eDx8AFiaS3Aj/4kbdLcVkTGfjWK+2hqqRvxpYDmG3Uwk0n08iqSnBgr9axk3OaAgAODSdCTx2uzx3qYRCNo79U2ncf3ASZ9ZqobuQ//KFZdy3bwKPvPsQDNPC116LrmS00mzjZ//b9/DxL74k1GIOYscKPquE8Vd1MHRVERzhBwv+mEdYWYIxLMIXO1ohvPEK8Fo6Bte/BzxbrwRdG3udDMdyY4LG/NOSO1YhRPAFjn2oNM0+OwfoDkfzJm5Xyk0Q0u2uZbCblshKHXYDnPS9D7vHU1gpN/ueci4W69g7ker52nQuAUohfBfwarmFfEpzA4tYjOBAIc2N2qstE5eKDXcWEGDnGS4W65FU6ryyWMZtc3Yy+U37JgAAz17g9whsVFs4sVrBu26exp0LeWTiKp4Z0E+wHX77yyfx1y8t4zNPXcSfPnkhsvMwdqzg11r8ChCGLrgOv2UGl2UywS81DNc+4c2rAaIry+TdjPx7bUuNNnRjjBgAACAASURBVPKckkyge8MQlbitcbZdMVg3Leu2ZeLkFzrA9vABsUtQyk3+PCE2737FswhlZauBQibuDuRjsPdRZNKPNaP5+yQWxpOgFFguda/L7FhYLHWXsjBYP4FoH3+13MRMtvcp59B0huvLn1y15/scmekK/sJ4Es12+KL4a6HZ7uD0WhW3zY8BsGcK6UoMzwTsEWDNYvftH4cSI7hrdz7w5iDi2v782cv40XsW8JYDE3j0ibORjr0AdrTgB48vAMTX4dsRPv8Gw375eyP8oKStXZYpKtJptjuIqzHu3B7Xw/ckbXldtgCQdIRZmOAHrHkE+kcTsMdvv5UB2GOeAbERfjkowndsG+/mq6VSE3PO6GEvboRfE2fpsPfB/6TDbBSvrbNUaqJj0T7BZw1i/uTzdlmttNzxE4xDhQwWS42+gXMnnYFuN+/KuV+bZ3uDBfcIXCrW0bEoDjlJ5ISm4Lb5XKCIP32+iLgac28Qd+3O48RqJZLdBt88tY5qy8RDd83hI/ftwUq5iWciurkwdqzgV1tmYMIWiGC0gtkJnKXjRvj1NrepxwuLxEVF+c2AEdGAbdMoMdLj4QdZOknBC9a7S2A4EX6618NnFsq0zzYBuhG+aA+f9/NhtfYrPYLfcGfNe2GCXxLo4Q8UfE/i9qLPU2cwUV4NmO1/raxxuo0Pz9gi61+68tpKBSldwcJ4931bcPcGixV83vtw1+5xvLy0xf39f+Z8EXfuzrvB25FdWXQsinMR9Ah85+wG4moMbzkwie+/ZQa6EsOXX1kRfh4vO1bweZubvGiih6eZFncWPtDr4fOWhHsRLfhBW78Au5kol1A9VToGtyQT8K45FCOsrOKHd1NOagriasyN8Fk0Ou0TFO//L9rD51k6uYSKlGe3LaUUi6UG5vL919W9yQv08AMEf1cuATVGeiJ8JnTekkwA7oz/9Yq467IsiiuVVt8N+VBAaeaJlQqOzGR7njqHWRR/LbBksvdJ587dY2i2LddaYtQNEy8vlXGf4/MDwGHnezi5Kj6h/MyFTdy5kLdHhsRVHNs3jm+eXhd+Hi87VvDtuu3+iJARxXjkoAg/4QhYudF2pz5OBVybu9dW0CNks21xu2wZbBZNu2Oh1GhzbRPvdQmL8DnLTxiEEExlug1Cq+UWMnGVe3NwPXyhlg4/WCCEYNdYwi3FLNYMtEyLa+lkE/ZGMZH7djfrBmKkP/+jxAjm8sme6PhCsQZdifVVD2lKDPmUhvWqOEunWDdgWrQvwt87mYYaIzjlE/yTqxXc5PHvAft7GktqEUT4DSQ1pSf/c9duu9HrBd8ymOcvldCxKO7dN+5+7UDB7hE47bs5bBfDtPDK4hbu2ds911sPTeG1lUokY7UZO1bwizWjpyPRj91pK86XC4vwATviK9XbWK8a0JVYYFmmG+ELKs0M2vrlva5yo431aguU2o05PNwqHWEefngn9Fw+geUt+5d/tdLs84cZ2bj9MxZl6VgWdUZE8//tzI4l3AifjSrmCX4sRuyfuUDB36gZGE/p3HzM7olkz3iFS8U6FsaT3BLbqUxcqOB3N371/ow0JYaDhQxOrHTF8kqliY2a0VOSyZjPJ7FYEuzhb9axZyLVU623ZyKFfErDC5d6Z/2w5fD37O6KcFxVsG8yLTzCP7deg2lR3DLbfR/edmgKAPDtM9FF+TtS8NsdC+WmGRitAk4dvqAIv2NRmBaFrgQLaz6lodQwsF5tYTLTPwSMERce4YcLfi6hYavRxmqZX0/OSAouy6y7SVu+4O8aS7rCeqXc5Pr3gN23oMQIqi0xwlo1TFCKwE7oXbmk6+G7o4o5gg/YHcslgd22Qd3GgF2a6bVDzq/X+/x7xlRGFyr47EmMZ7ndMpvFa55RBq8s2R/fOpfrO3ZhPBmJpbN7ovfnQwjBnZztX9+7sImDhXTf8MDDMxmcvCI2wj/lvB6zvQDgtvkxjCU1fPOUFPyrgv2SBf1yAGKTtiwXEBbhT6bjWK8a2HAEPwgW4TcFRfi24AdfVyEbx2q51beX1Q/LA4jqtmWWTtBeABZJU0qxwkkIMgixd/OKivCD5ugw5seTWHVq3hdDInwAGEvpwpO2QUHM7okU1qsGai0THYvizFrV9Z/9TGXiQm0DNoyNd1O+eTaHpa0mtpzfyVeceTlHOYI/l0/2lJZuF0qpI/j9N747d+dxcrXiPmlSSvHcpRLu2TPed+zBQgYXN+pCm9VOrlYRI3B3FQC2Nfe2Q1N44tRaZJNDd6Tgs2YX3lx3hqbEYFpUSN1r2D5bxmw+gaVSA+u+5dh+uh6+KMG3ApO2gB1VrZSbWHIepYOsExbhi2oKC6vSAWzBN0wLq+UWFjcbfQ1EXjJxVZiH7/ZJBFhuBwtpWNSeU3NuvYpsXMV4gHU4ltSwJbIOv270bfxiMHE/sVrBxWIdLdPqqXP3Yls64q6LPR36m88AuM1Vr63Ykf1Li1vYP5XmNh7OjiVQaZnCpp8WawZqRsedKOrlrt1jsGh3YNv5jTqKNQN3cwR//1QapkWF5hdOX6lg72S67+n7nTcVsFpu4fiy2CcKxo4U/KBqBi+sSUpEty2zX4IarwD7sX+13MRquRkq+F0P//pYOgvjSXQsihculaDECCYDksndKh1REb5dxhr0nrGo+dtn1mFRYH8hzT0OsBOkoiL8QZ3QLCI7c6WKk6tVHNmVDbTn8klN6GiFYs3o2/jFuMUZS/zacqXb2MTxyQFbmKstU9jNe7XcxERa5/ahsOtiVs7Li2WunQN0f+bePoftwOb+8yL8OxbsxO3zjo//tLOs5d69/YJ/wPm3J7I08+RqtcfOYTxwpAAA+NqJaMY5CBF8Qsh7CSEnCCGnCSEf4/w9IYT8J+fvXySE3CPivEGwcr6grlGgG42LEXzH0gkR/Ll8Eha1/U5eJMS4nmWZADCft38ZnjpXxK5cgpvkA4CEKjZpW2+ZSIV0QrMqjsdfsuuS90/x7QnAjvBFlWWyUctBlg775T99pYqTq5XAKBpgeRsxgm9ZFJv1dmCEP59PIhNXcXy5jBculaDGCI7M8N+zQkbsknVeSSZjJpfAwngST53bwMpWE4ulBu5YGOMey8pbRSVueSWZjKlMHIemM/gHxy//xsk1zOTi3PeM/dvz9xNcK4Zp4fx6jXuu6VwCR2dz+MaJNSHn8rNtwSeEKAA+BeB9AI4C+Agh5KjvsPcBOOz891EAf7Dd84bB6pGDkmkA3Fb4tgBhZeIcFuHPeqYp7p8MjlajKMuMD4jwATuqOhjg+QJ21UlCiwlL2lYDtl0x9kykkNYVfOW4Pdkw7D3LJMQJ/qBO6JSu4sBUGn/90jJK9TZuChBVoBvhi7ANK03bmw/y8GMxgrv35PHtM+t48uwG7lgYC7TLprL2a4hK3F4pN7kJW8b9Bybx1Lkivvqa/bN8hxPF+pl1RleI8vFZL4K3wcvLe26edm9ET5xcwzuPFLhPaxNpHfmUJizCP79hV+gcnuYHC++6uYCm2RE+4A4QE+HfB+A0pfQspdQA8FkAD/mOeQjA/0dtngSQJ4TMCjg3l1OrVeeHNNjSGVSp07HowAYtN2kbMFoBgLtmDgAOTgeL17BJW7NjDZV0bg6I8L2Pu94SMR4pXR0qwh/Gg60bJrcGnxGLEXd70v6p/soJL8Mmbbca7YHJMDdpG1CWCdiC9ZpTavim/ROBx42ldFA6eM5PpdkeaK9s1Gxx5s0TYrz75mmcWavhexdLuP/gZOBxbvPVAB/f7FhDjQZeLbcwE/LU+uDRGZTqbfzrL76M+XyyrwafMZ2NQ4kRN58UBKUUG0PcrC5v1jGV0QMrwT5w1xzaHYof/YNvo9Iy8Y+P7Q58rf1T6aEE/0qlOXCi7CmnxPNwQLDwiw/ehMceeRvUkJzgtSLiFecBXPJ8ftn52tUeI4Rqy8Tfvrzi1rQGwSL8MDFvdyy867e+jkc+873Q12LReJils28yhUI2jmxcdf1DHsOWZf7Lz72A+3/9q271QxCDqnSUGMH7btsFAPiBozOhrzXMTPz//M1zuPPffRlPnAx/JA1afuLlh++wY4Ifuj08NsgmBidtv3NmA3f/+y/jN/72tdDjygHLT7z8kHNdN+/K4ugs348GeofmBbGy1cTb/++v4QOf/GboDdwtRAgR/PffPou0rkBXYvjQvcHiNTWEpdPuWPjAJ7+FN//aV/HKUvAmqo5FsVZtBVZRAcB7bplx7Yufefv+wJyHqsQwk41jaStc8H/+s8/j3l/9Ch594kzocReLdSxwEraMW+fG8IO3zmCx1MBbDkzgGMe/Zwwj+I8+cQb3/dpX8ZE/fDL0Z3lytdJXoeOF12chivDfuOHgXZ3/FjfMMfaBhHwUtu2DPXv2XPXFpHUFn/zxe3osFB6aYl9SmIf/2rJd8XCxWEclZC2hMYSHTwjBl37hHTBMq2+yopdhPHzDtPDYC0sAgCdOreFH7pzjHtfuWDAtGhrhA8Dv/M934ZF3V3HrHN9bZST1waOb//TJC7Ao8PnvXQ58dAec5SchHj4A/NPv24c7FvKBni9jmAj/C9+7DIsC/+Wb5/GLD94UaL9Vmm0kNSX0Z/SmfRP4i599a19Dj5+8b98Aj89/7zJK9TZK9TaeOLmG99zCv+myyaFBHj5g++V/8/PvACH8RCWDlQWHWTpPnFzDq079/B8+cRaf+PDdAddlbyYLatgD7KDizz56P86uV7llj15mB5RmnlytuP/2P/n3p/ETb9nnVpD5uVRsuF21Qfzuh+/G0+eLOLZ3IvRneWAqjS98b9F5Mu2XzdVyE7/95ZOYycXx3XNFPP7SMh66ix/Tnr5SxZ6JVGgxRVSIiPAvA/CGEwsAlq7hGAAApfRRSukxSumxQiFYMIIghOCdRwqhyTSgK6xhd2JWSgbYjSxBDOPhA7YX6F927cf18EMsHe/MlFdD9nOGLT/xktSVgWIPsAg/WFhbZgfnnc1Hz/u6GP1UA/bZeiGE4N6946HiCwCZuIZGO9zzZNMRjY4VGrGWG8EL5r3cuTsfGm0D3hHJwYL/1Lki9k2moCsxfPccf2Qv4BmNnA62mgB7dk6Y2AO29ZhLqKGC/9cvLmM8peFH7pzDE6fWA/MQrH+jkA3/dz2e1nHvAFEF7OKGsAj/i88tQlMIfv/H70G5aeLLr/KHjbHx0P6mKz8JTcHbDxcCbxoMlrgN0oHHnl9Cy7TwmZ95Cw4W0vjjb58PfK2TqxUcCvDvo0aE4D8N4DAhZD8hRAfwYQCP+Y55DMBPOtU6bwGwRSldFnDua2YYS8e7Co4JGY9hPPxh6Xr4wZH0ec+j5dmQygGWYE0M+Mc8LMkBe20vbtRhUdYx2Qi9mdaNTug006vBnZgZYOts1gycXa/hw2+yY46wzUqVVvCT3NXiCn5AhN+xKJ49X8TbDxdwx8IYvhswox3w7gQIjqSvhkI2fLzCsxc38eb9k3jXTQUUa0ZgYMEmmYZF+FfDnKfhjsdzFzdxdG4M7711F8aSGr4VMGxseYs/HvpaYdVZZ9f5/3Yef3kZt83ncLCQwQfvmsdzl0puQ5qXdsfCuYAKnevBtgWfUmoCeATAlwAcB/A5SukrhJCHCSEPO4c9DuAsgNMA/hDAv9juebeLPkSE7/U4w345ho3wh0GNEcRIuKXDvMTb58dCa5bZU0JCwHUBdoQflrQ961zXu2+eRseiocm3QUnbq4HNxA9KjrLreo8zgjZsd2q5Ebxg/moZS9pPAEHNVyvlJmpGB7fM5nDn7jyOL5cDI+nNuoGEFhsYiQ7LVCYeODGzWDNwYaOOO3fn3Q1R/jEEDDbJNMzDvxpYwx1vEYplUbyyWMbt8znEYgT3H5jEt05vcF/HXbkY4uFfDfucootza/2BX7Vl4sXLW3jgyDQA4MFbZ0Ap8Pec1YjnnRk6QQnbqBGiBJTSxymlRyilBymlv+Z87dOU0k87H1NK6c86f387pfQZEefdDsNE+GvVFvZNpkAIQndgGp3BjVfDQghBQgtfZH6xWEc2oeLITDY08dYY0tIZlpQenrQ96/wyPHCTbcWd3wi2waotc2DSdlgGRfjsl3//VAr7plI4cyX4aa3c5M/Cvxa8exB4sBvi/HgSBwsZNNtWoJ1hDwMMt5CuhqmQCJ91n965MIaF8SSyCRWvLvEj/LAu22uBNV/xgoWLxToqLRO3Ofbjm/ZPYLHU6Fk5yWC25yB7a1iSuoK5sQQ3cfvM+SI6FsVbDtiVUTfNZDGTi+MfOE8fbHJoUElm1OzITtth0IZovFqr2NUHEyk9dPXaMEnbqyGuxkLLMterdqPLdM7+pQ2KCpktNChpOyxJLdzSWdlqYCypufmAiwE2WMeiaLat0Dr8q2HQTPxuPXYKB6YyOBfwWA6wfbZiBF9XY0jpSmDSlm13ms8ncdCxDM5wIkggvMv2WiiEzNNhNuGhmQwIIbhlNofjAZbOaqWJybQ+MM8yLF3B7xfxl53ci3ddIWDvrPVzsViHEiMDizeuhv2FtPu06OXJs0WoMYJ79toJYkII3naogG+dXu8r0Xxtxa7Q4XXZXg9GVvDjQ9ThrztdsRNpHcWQmmWRlo59beER/kbVwGQ6jkImjnaHBgoKu2mIivCTerils1a136+pTBwx0psD8cISv4OqdIZl0NarS8U6ZnJxJDTFHtAV4hGXG21hlg7gTMwMEnwnip3LJ3DAM7KBR9EZjSyKqYyOSsB4hbPrNWTjqtuRe8uuLF5bqXDfs0FNV1cLE+hlzpPOS4tb0JWYW5BxdC4HQuyv+7lUtDeRiaxl3z+Vxtm1at/78NS5Ddy5O99ThPD2w1Mo1dt9BQLHl8vYP9U/Q+d6MbKCP5Sl4xX86xnha+GLzDfrBsbTmvsYHSSsLBpP6gI9/JAIf63SQiFjN89MhsxcZ/tsB1XpDIvr4QdE+Jc3u2sIZ8cSqBsdt97eC6VUaIQPOBMzAyydxVIDE2kdKV3FVEZHJq66TyN+NuvBo5GvhW7zVf/P6OxaDQcKabei5vBMFnWjgyVOvuhKpSUsYQvAmckT41o6ryyWcdOurBtYZeIq9k+lXQvKy0XOWOTtcmAqg3LT7NGCumHipctbeLOvAe+tTh/QP/hGHR9fLuPoEBVxUTGygj8oadvuWKi0TIyndExmdLfTkYfoCD+hhte7F51tXmx+SdCjOXsNEdVDgO3hN9qdwOh4zTMnKGwEb+06R/jexqDZfHAE2TItGB1rqLLMYbHHK/CDhcXN7npEQgjm8onAOTLFkFn410JYt+259Rr2T3W7wdkkzlOcrU+rIbsKrgX7fUj23VwopXhpcQu3zfc2ut02N8YV/PMbtZ7udhHs5wxRe+5iCaZFcZ9P8AvZOG6ZzfXMtt9qtHF5szGwoz1KRlbw3cargEi66um4HBTht4YYj3w1hEX4bIjWZFrHVDY4SgM8Hr6gyo6EroDS4Aoir+AXssGC7y4/Ee7h8yPp9WrLFbiuZdAfrQ6ahX8tjIUsQfEvQJ/PJ7mRrWFaqDRNoUlb9nNa9/2MGkYHi6WGazEBdoQP9JeztjsWrlRa7gwcUcyOJbDsex8ubzaw1Wi7/j3j9vkxLG01e0YtbNYMlOptHJgSK/js9bw+/nfPFREj/Cmbbz88hWcvbLoWJlsEc0tId3bUjKzgD5qWyUr8MnEV+aSOctMMjGxbZge6GhvYVDIscTUW2HhVbrbdIVqsCiRo3smwjVfDkgpZglJrmagZHU+ErwfOanH32QqK8NmNgxfhG6aFUr3tXlfYgK5hxipcLUETM9kCdDatFLATlrwIvzTEWIWrJShYYNGrN8KfSOuYTOt9gn+lYq/FFJkYBZzmK9/Ph0Xxt/nsEHYDeNET5TNBPhAyUvtamM8noSmkJ8J/+nwRN+/KcW3Atx2agtGx8JTTUMeaEf3fw/VkdAV/gKVT9sxFzyRUdCwa6F8bpiXMvwcQWpbZbcDRXWEKWpTNEqzCqnRC1hyy5exsuBeL8Hk3SXefraAIPxYjgUtQmBXHIvxCNg5Cuh2iXtjPXFRZJgCMpTTu0LZSvY260XEtHcAuzyzV2+77wyjWB+93uFrYz8n/FHYuQCwPTWf6lpGzKHxQ9/jVMjeWwJVKs6dz+uWlLagx0rcL97Z5J3F7uSv43ZuW2EoY1dnRy2b7tzsWnrtY6rNzGPftn4Cuxlxb56lzRRwopIWVsF4LIyv4g5K2rMQvl1C7lkGAR9wSLPhhZZnMHhhLaYirChJajJuABICmyap0BCVtHYFucMYrsAFhbELpVDoOw9kt7KdmhO+zvRaClqAwQWO/ZJoSw2Ra5ya6KwNm4V8L+aQOw7T6fp4skveO7p0PqEFnFWIiBT+hKchyxiuwktX9Pjvk8EwGp1Z7K3WYLRa04vFaYbsjVj0/o5cWyzg8k+17Ws0mNByYSuNFT2PY2bUq1BgJHIu8He7ZO47nLmzCsihevLyFRrvjNqf5SWgK7ts3gW+eWkez3cF3zxXx5v3BU0yvByMv+EERvmvpJFQ3kg6qAjFMS5h/D4SXZXpvRPafWuDETBbhJ0QlbTW2BKX/PWOloWycAPuT9/RRb4lN2gLBS1CYoE15atgL2QS37b3r4Yu1dID+iZlM8L2WDhN8v62zXut9ehJFgbPq8Nx6Hbtyib4KqsPTWZSbZs8TAUt8C4/wnfeBNcxRSvHK4hZuC9iUdedCHi/6Ivw9kylhvQFe7t0zjkrLxPGVMr5yfNXdQxvE2w9P4cRqBZ/+xhlUW6Y7AfZGMcKCz6Zl8n1576q77IAqkJYZvmTkaglL2naTybaQjCW1YA/fyS2IGreadNccciJ89uSRZIJvi9MmZ6xAdcA+22shaAmKP8IH7LkvvAh/a4hZ+FdLULcta7ryWzoAR/ArvbaUKKYycaz5Iny7uqW/O9Wt1PHYOstbTaR1xS2LFQVbxMM2TC1vNbFRM3B7wNTUOxbGcKXSworzxPHqchlHIupkfceRApQYwWMvLOGvX1zGm/dPhO5q+NC9C8jGVXziK6ewdzKF+w/ICP+GQAiBrsQGWjrZhIpMXOv5mh/D7AiP8IMsHXYjYjZTLkzwjfDlJ1dLmIfPxJIJ3HjIpEiW9E0Lqh4CnEXmnBsyi2C9YjmdjXM9fP/3IIJ8gOAvlRpIaLEem2Y6m4DKWQCyXm1BjRGh1wXYm6/8ls55X0km49BMf2nmcqmJ2XxSWLECY24sgbSuuItCWPQeNNH1dme/xAuXSyjV7TlAd+yOJjFayMbx9sNT+H++cRYXi3X85P37Qo+fzMTxez92N374jll86sfuiXTW/TCIvTW/wdDV2GBLJ9718IOGcxmmhbggnxxwqnQGWDoZ19JR+6I0RrNtCfPvgW7yl9dt6xdLZmXwIvyaYSKuxoR2QWYTKrfUcq3SQjah9ni/M7kE1qstdCzas8O33GhDV2NCuyBZ9Ofvhl4sNTDnE0slRrBrLOFG/4z1aguTGV24WBQycXyz0q0TLzfb2KgZ2McR/EImjrGk1hvhl5vCK3QAOxg7NJPFqSv2zeW5S5vQFBK4/PzWuRwSWgzfObPh/uzuDFkytF3+1ftuweJmA7fPjw1cGgQAD9w0jQdumo7seq6GkRZ8TSGBEX652YauxNzkFhC8uq8l2sMPsXTYTYdVuOSSWuD8labZESpeqQERftwjlu6kSK6HP3jb1dUStARlrdpyRwQwprNxWBTYqLZ6xgJsNdrCo2hmbfmbr/w1+Ix5TmnmetUQbucA9lNPuWk6W9EUXHBmve/jNCwRQnDYV6lzcaOG9w3YRnat3LIri8dfWoZlUTx3oYRb58YC/y0nNAVvOzSFrxxfRceiSGgx3L0nOsG/aVcWf/cv3xnZ60fJyFo6gJ24DYrwq03TjaIHDecSHeEnVAWGaXFLGitNE5m46kamYR5+Q7SlE1KHv1Vvu1E9EL78ozbEtqurJRPXAj38KV8ZHBN5v48fieAHefilBreKZGE8hcucCD8SwXfeF1bqe84Zdrdvij9h8vBMxq3FL9UNbNbbocvlt8N9+ydQbpp49uImnr9cwpv2hW/K+qE7ZnF5s4E/efIC3nXTtND80E5ipAVfV2OBjVcNo+NGtINa96OI8Nnr+qm22j2LQ3IJDeWA5dxNwclk5uHzxj6UGkaPWGpKDJm4GmjpiKrBZ7CkrX9y6HpAhA/01+JHIfgpXYEaIz3NV812B+tVA3OcDtWF8SRWys2eJ8/1SjSCz6phmIXEFuvsneCL+KHpLIo1AxvVFrdBSyRvdpKb/+4vX4FhWviBW3eFHv+BO+dxz548kpqCh995MJJr2gmM9G0wLGlb90THmhJDQouFlmWKmlcDdGfftNpW32NstdV98gCAXFKFRe2v+7v97KSteA+fG+FzxDKf4peM1j03U1GwSpGa0fs+rFdamDrUW844ExLhi1rkwSCE2N22nvfBOwffz/x4EpTaJY97J9OglNqWTlZsSSYA7HOqcc5v1HDf/gmcX69hdiwROIrDW6nDvgee3y+C+XwS9+4dx7MXNrF7Iol7B+zCVWIE//3h70PHosJmWu1ERvqdCUvaNtq9opTW1b4OSAYbrSCKhBvh9wurPc2xN8IHwG1wEu3hq0oMuhIL8PBN17dn5FNaYFmmcA+fswSlZdpTMf2djSxa5kX4ecERPuDYbp4Iv1uDz4/wAbi2TrlpwuhYfU8pIpjPJ6HGCC44Vs6pK1UcLAR3p7ItTaeuVPHyYhkJLcYt4RTFb37oDnzkvj1DV7coMSLFfgAj/e5oIRF+w+gVy7BZ8IZpCf2HxiJ8Xmkm8/AZ7jwdTnK0YXSENV0xElqMX6VTN/oj/KTOnSNT9d20RMDrhmYD7yZ9Yqmr/G7brYa4bVde8im958bXrcHvF3y2ku+ys7Gp2zgmXvBVJYaF8STOb9RhdiycWK2ETnLclUsgE1dxcqWC713cxB3zBn6A2QAAGdlJREFU+UiamxgHChn8+j+6HXdEWHEzamzrp0UImSCE/B0h5JTzZ99zFyFkNyHka4SQ44SQVwghP7+dc4pEU0jgAhR/hB+23i+K0Qr26/afz7ZuvJYOv+wP6P8eRJDS1cCyzGEtnWorGg8f6O2G9s/38VLIxnu6bTuWPQtftIfPzu+td18qNRAj/A7VXWMJxEg3wl91Sk1FjiD2crCQwYmVCs6u12CYFo4GlD4Ctj31pn3jeOyFJby8uIW790ohfqOxXZX6GICvUkoPA/iq87kfE8AvUkpvAfAWAD9LCDm6zfMKISxpay/Z7opSUldRDxmeJjbCD07aVpq9SduwLuBmu4OEYMFP6krf+9DuWKgZnZ4qHSDE0mn25iFEkOVE+Exk/RE+YPv43gi/HEHTlfdcbPcrAFwuNTCTS3CjY02JYXYs6Qr+ZXfmTjTWyV278zh9pYqvn7AXbg+a5PjBu+ex1WjDtOhQNeiS1xfb/a17CMADzsd/DODrAH7ZewCldBnAsvNxhRByHMA8gFe3ee5toymxwGaqpi9hmtIU7tAwgEX44oSVnZdXDVP1bWRym8I4s+BFl2UCztYrX4Qf1KE6ntKx1WjDsqjrwVJKUTVM4e34PA9/w+2y7Y/wp7NxvLbS3YUaRZctYyYXx1aj7da7B9XgM+bHk66lc3mzARLwNCCCe5w57r/95ZOYG0sM3LX6/ttn8dzFEibTOu4ZkEiVvP7Yblg64wg6E/bQdjJCyD4AdwN4apvnFYIeUodvR/iDLR1KKYzO9YnwOxZFzej0RPhBJaOU0ogsHQWNdu+5/IPTGGNJDRbt7VCuGx1QCuERPs/DZ6ORgyL89arhLpmOUvBZ3T+b63Op2MDuieCIfWE8iUtFJ8LftIeZRZWMvHtPHmldQcu08P1HZwaOSdCUGP7tB27Fz73nsPCRCpLoGfiviBDyFULIy5z/HrqaExFCMgA+D+AXKKX9a+a7x32UEPIMIeSZtbW1qznFVTOoSsdbnpaK871rZgkJ9fCdqNzv4Xvn+zCCqnRapgWLilt+wuAlr1nJoT/hyZZueydFuqMh4mKFNRvvXwazUTUQV2PcmT3TuTg6FnVvCq7ghwzCulZYqeeqU1+/tBUu+Aem0lgpN1FptrG4Gf40sF1Suopf/9E78ODRGTzyrkORnUfy+mBgmEUp/f6gvyOErBJCZimly4SQWQBXAo7TYIv9f6OUfmHA+R4F8CgAHDt2jJ9RFURQlY5lUTTbVo8dktIUd1OTFxaFi12A4kT4viodnuDH1RjUGOnrMnXXG0Zg6fiXZgT5395u273OkEDv2GmRsM5d7/uw5nSo8iLR6axTi19uYTqbcCt6RM6cZ7Al3yvlJiZLDVAK7AkR/Jt32YnTEysVnL5SxbtvjnYOywfunMMH7pyL9ByS1wfbVanHAPyU8/FPAfgL/wHE/m37fwEcp5T+zjbPJxQ7wud1qPbvgk0GWDpGBILvlmX6I/xmf3RMCOEu/2gI3mfLSOr9C9ZdS6dP8G3xLNZ5Eb7Y61KVGJKa0mvpVA1Mcvx7wI7wAeBKxa6Cccsf0+KrYVjC9cJGHRedGe9h9eu3OJUyXztxBRs1o2+Pq0RyrWxXpX4DwIOEkFMAHnQ+ByFkjhDyuHPMWwH8BIB3E0Ked/57/zbPKwRN4Q8pY8Lu9/B5lo67wDwKD98X4bujkX3RcSah9g12430PIuDlMtjO1f6krdbz9wD/piUK/0z8jVorcGmI223rVM+sVw1oCkEuKb75PBNXMZ2N49x6DWedGe97QyL8ubEExlMa/vTJiwDsNX4SiQi29a+bUroB4D2cry8BeL/z8TcBvC6zO7pCuB6+uylK6xV806J9JZhGBIKfcD18n+C70bFP8DmDw3jfg6hr66/ScbZwBXj4m7Xuzaga8D2IIOubib9RNXDLLr5YFtxu25ZzbAuTab79I4IDhbQr9lMZPXSvKSEEbztcwF++sIS4GsPRWRnhS8Qw0p22QUlbZoekeiwdts+1V+xYYlXsLB37x+K3Tlh07F/Bl02ofUnbqDz8lFOH7x3WttWwewP8deW5pIYY6Z2Jz8tDiCKf0twEMaXUsXT4wqqr9vIRr6UTxbwaxoFCBmfWanh5cQu3zo0NvLH8b+84gFxCxcPvPCjclpOMLiMt+EFJW2ZZJH0RPgDUfSWJboQvdOMVvyzTv/yEkeXMgo/Kw8/ENXScpDbDPymToThbmnoE37exSySTmbhbe19p2TNoeDX4DHvzlRPh1wxMRuDfM26bG8NWo43XVipDWTS3zY/h+V/5AfyfDx6J7Joko8dIC76uxmBatG+kLovikzpH8PsifPGWjqrEoMRIf1kmW34S74/w/Y1XvJuWCHjLYMohY4XH0zrX0hE9PA2wrRK20pDtgQ2rupnOJbDGIvxKKzDBK4J33lRwP/7+W4brUL3R6/AkO4+RHo/MLAijYyER6wojayzyimXQej+WWBXulaux/qQti/D1/qStP8JvRhThM8EvN01MO4FqqR4i+L7BYZWWCV2NRdJINJmOo1hrwbKoG7nvChl3PJ2N4+RKBe2OhZVyEwsR1rvP55P4vY/cjc26gbtlh6rkBjHSgs9smHand4xCw7CF1jtLh33sH5HMonCRu2MBu/nKX5ZZa9ndv/7IL5uwk7aUUtcbbkQU4bNGL2+Ev9VoB47VHU/p7pgA9j2IHqvAmMrosChQarSxvGV3qoaNJJjJxbFWbeHyZgMWBRZCKmdE8COy1l1ygxl5SwdAXy1+3eBE+K6Hz7d0RCZt7dfrj/BrLZPrfWfiKtod2uP5R1WWyZtZs9Vo941VYIz7ln9EMTiNwRK0G9WWu9B8lrNVijGdTaBjUTx/aRNAeDOURLITGGnBdy0dX3KUZ4ewTk6/pcOOFdl4BdgWEa8skyf4XV+9K8IsaSvaavKfi1KKUoiHP5HWUawbblXPVqMdSYUOANeDX6u2sLLVRD6lhVpabB79P5xaByAFX7LzGXHBt+0Pf2kmt/FKU3v+jsFEWbSwxtVYX1lmLWBTFC+R2mx3QIj4G1HWZ+nUjQ4M08J4QHI0n9JhmJZ7Ayo12m59vmhYNL+42cDyViPUvweAOxfs+vYvfG8RcTUmfL2hRPJ6Y6QFXw8of+RFxyxS9I9IbkUU4cfV/i7gYEvHFmGvzcJGI4tuJPJH+CwhOxEg4hNp+9rYrJqwBO92mc8nESPAxWIdl4oNd11gENO5hDuY7Pb5MSiyKkaywxltwfckbb00jA7iaqxHAILKMpvMwxedtFWVvrLMSjM8wvdW6tQjGI0MwN1UxRq9WMllkIfP5ukwH3+zbkQW4etqDPPjSZxZq+LMWhWHZ4LX9THYYLL33z4byTVJJK8nRrtKRw0QfN9oZKCbwO2zdNoRJW21/uUsNcPkDh1jUb+327ZpiF1gzlBiBBlPo5cb4QdYOuzrxZo9e36r0XZn7ETB3ok0vn5iDaZFcdMQgv/x99+CB4/O4G2HpiK7Jonk9cJIR/hBSdu60UHKJ5axGEFCi7kVPIyW2YGmEOF2gB3h+y2dDrfCJcupnGm0xW+78p6PefhM8PMBUTsT9826gXKjDUqDjxXBgULavSkfngnf3gTYVt07jhRkk5NkJJCCD/TttW0Y/F2wKV3tt3TaYtcbMuJajNtpy7d0+mvj60Y0lo59vu6QskFz5PPuADXDvTmMp6OL8L/v4KT78ZEhInyJZJSQlg766/CDVgPy9rm2zI7wpisASKhKTx2+YVowOlZfly3AX+/XaHfczVmiySY0d5TDZr0NQoJXA46ndMSIPX54s87m5kcX4T94dBf+2Vv34/6Dk9wl4RLJKDPagh9o6ZhcO4Q3C170AnNGQou51UJAt8OXZ+noagxxNeaOXmDHR1VmmImrbrS+WbMHpwVZWkqMYCaXwNJWw30aCCrhFIESI/iVHzka2etLJG9kRjoECk7aWu44ZC9sNLCXZrsjvCTTPZcnXzBo6JjXZgGCSzhF4D3XMFU3u8YSWNlquuMO5vKy3l0iuRGMtOCzxit/hN8wTCQ5Nk1KV/vr8E0rEuskpatoti10nEme7hz5QMHvXYJSDWjSEsFYUnPXGtqCH+7Jz40lsbzVxFKpCU0hkawRlEgkgxlxwQ9I2rY7PYPTGCldQa3Fs3TEv43uKAfniaI2IMLPxHvXHFZbZmQjDKYycWzWDbQ7FjZrgztnZ8cSWCo1sFhqYHYsKStiJJIbxLaUihAyQQj5O0LIKefPwLmvhBCFEPIcIeSvtnNOkcSDLJ2AGvZUXO3x1QG70zYKwWeWErN1Blk63tp4s2Oh2bbcJinRFLJxUGpX6KxXW6Ez5wFgNp9Ey7TwyuKWtHMkkhvIdpXqYwC+Sik9DOCrzudB/DyA49s8n1CC6vAbASWNaV3pG4/cNK1IGpzSrLPXeaIYtBrQ66uzp5CoplKyfayLpQbWqi13CFkQBwtpAMDZ9drAYyUSSXRsV/AfAvDHzsd/DOCDvIMIIQsAfgjAH23zfELROBE+pTRwLAGvDj+qCN8/ymGgpZNQ3ZsCK5nkdeWKgAn+y4tboHRwEvboXHel39HZwev9JBJJNGxXqWYopcsA4Pw5HXDcJwD8EoD+BbI3EF5ZZsu0QCl/+mU6rqBmmD0LvI0Ik7aA19JxovYAmyaX0FwPn0X4USVtC87c+RcvbwHAwKh9OptwRw+/80gh9FiJRBIdAxWBEPIVALs4f/WvhzkBIeSHAVyhlD5LCHlgiOM/CuCjALBnz55hTnHN8Kp0GiGLQ1K6CkrtmwK7ITTbHSQiTNqyCL+7z5Z/c8nEVXfrVbUV3aJwoBvhv3CpBGCw4APA5//378PFYn2ogWYSiSQaBioCpfT7g/6OELJKCJmllC4TQmYBXOEc9lYAHyCEvB9AAkCOEPKnlNJ/EnC+RwE8CgDHjh2jvGNEQQiBrsbQ8lg6LCkb1HgF2PYKE3y7LDOCpK3WG+HXDBMJLQY1oHs0k1BhUfsG4T4NRCT4CU1BIRvHqStVAHYVziAK2bh7o5BIJDeG7SrVYwB+yvn4pwD8hf8ASum/opQuUEr3AfgwgL8PEvsbQUKNoenx5VlEzduUxBuRHFWnbV+EP6CRio022Gq03aeBqJK2AHDI2WG7dzLFLWGVSCSvP7Yr+L8B4EFCyCkADzqfgxAyRwh5fLsXdz1IaAqanpk1zZAIn3niNU/zVbMdzSwddsOpeSydMMFntfDFmuF6+VFF+ADwpv0TAIB79wRW4kokktcZ21IESukGgPdwvr4E4P2cr38dwNe3c07RJDQFTbM/wg9qvAK6SVGzY8G0aDQRvnN+1tm7FbI3FuhOq9ysd4eUDaqP3w7//K37UW+Z+In790Z2DolEIpaRfxZPaL27Y5lnntT7o3YW4bPELuvQjaTxynnCYH78VqONXKjgd1cJbtYNxNVYZPPwAWAspeH/+mE5pEwieSMx0qMVgDBLJyTCd24KbHxxFI1XsRhB1jMuoTxQ8O2EaLFmoFgzMJHWhe+zlUgkb2yk4KuKL8IPTtqmfbXxzAqKIsIHgJxnSNkgS2csqYEQZ9FILbq9sRKJ5I3LyAt+XIu5i8iBblkmvw6/18N399lGkLQFbMEvN+za+kGCr8QIxlM6inUDxboRqX8vkUjemIy84Cc0BS1PhM/8+aDhaUA3wnefBjj2jwjGkirKjTbqRgemRUMFH7D3x25UnQhfCr5EIvEhBV/rtXTcsQQBKw69xzTatvBHtTs2l9BQbrZdW2eQ4M/lk1gqNbBZb2NiwIx6iUQyeoy84Ce1WE/Stm6YiKv8jlYlRuy9tu6M+mD7RwRs0ciwgr9nIoUTqxVsNdqYysiuVolE0svICz6vDj9MwNPx7ojksASvCFjStlQfTvD3Tqbcm9duZ1iZRCKRMKTg+y0dwwwdFeAdkcwsnagWjYwlNdSNjrsLdnrALJq9k2n34z2TUvAlEkkvUvBV29JhI4/rrU7gREqArTnsjfCj8/DtGwkbUjadCx9SdpNnEqWcOy+RSPyMfKdt3DP1MqEpQ0T4iiv0bBtVVJYOq7Q5vlxGQou5N4Ag9k2l8dNv3YeF8VQkzWASieSNzcgLvjvmuG0LftB6Q0Y6rnIi/Gjexhknon95cQszucRQnbP/5kdujeRaJBLJGx9p6ThNUyxxWzM6w0f4bbuiR4lFM8KACf561XA/lkgkkmtFCr7a3VwF2GWZYR5+WlfdWTqDnga2y+xYAiyol4IvkUi2ixR8d1WhXc5Ya4WLeCquuN69fWx0rlhCU9xSzJtmMpGdRyKRjAZS8Jml40T4jQFJ254Iv21GlrBlvP2wvfT7rYemIj2PRCLZ+cikrWcZuWVR1Nsd7lgFRkpX0WxbMDsW6kb4sSL4zQ/dgZ95+37csZCP9DwSiWTnIyN8J8JvtDtomh1Q2h2SxmMsaf9duWmi3upEHuEnNEWKvUQiEcLIC35c7Xr4w8zGyTtz5kt1Y+DIYolEInk98f+3d7+hdd11HMffn+QmTZq/LW23NJ3bCnW1FGdHkM3BUNvJ/IPZAxWFjSKCT3SuIkj1ic9kiIh7IIPSbRYsE4mDFRnTWQc+K9Z1aLdMOuq6ZsuW2Nk27ZIlWb8+uCfpbUiTuHvuPbnnfF5Q7jknJ/f3PU345Hd/58+vqsCXtF7S85JOJa+LzmgtqVfSkKRXJQ1LuquadtM0fx3+7Afzj0Zeagx/LuDPJw81c+CbWaOotoe/HzgaEduAo8n6Yh4FnouI7cDtwHCV7aam8qTt3MnYpcble5LHDl+YnOH85PR8j9/MbLWrNvAHgUPJ8iHg/oU7SOoG7gEeB4iI6Yg4X2W7qWmvuCzz6gTmSwzpJD36sYtTTM1ccQ/fzBpGtYF/Q0SMAiSvmxbZZyswDjwp6YSkg5I6FtkvEx3JCdpL78/O30HbseRJ23LAnzn33jXrZmar3bKBL+nPkk4u8m9whW2UgDuAxyJiF3CZ6w/9IOnbko5LOj4+Pr7CJj68tpZmWpubmJianT9p277Eg8fmA/9dB76ZNZZlr8OPiD3X+5qkdyT1RcSopD5gbJHdRoCRiDiWrA+xROBHxAHgAMDAwEAsV18aOttKTEzNcHFq+YlGSs1NdK4pcebcZQB6PZWgmTWIaod0jgB7k+W9wDMLd4iIt4Gzkm5LNu0GXqmy3VR1tZWYmJrlYjKVYPcyvfb1Ha2cHi8Hvnv4ZtYoqg38R4B7JZ0C7k3WkbRZ0rMV+z0EHJb0D+ATwE+rbDdVXW0lLr0/y4XJGZoEXUuM4QPc2NM2P96/zlfpmFmDqOrRChFxjnKPfeH2t4AvVKy/BAxU01Ytda1pYWKqPHdsd3sLTcs87rivp/zkyib5KZZm1jgKf6ctzI3hl3v4vSsYork5mSC8d20rrSX/F5pZY3BacXUMf6V3zu7s7wFg64ZVc3WpmdmyCv+0TIDutmRIZ3KGnhWMyX9m+yb27dnG53bcWIfqzMzS4cAHOtckJ23fm+YjyXDNUlqam9i356N1qMzMLD0e0qF8Lf2VgLP/nZx//LGZWd448IGNXWsA+OBKsLHTV92YWT458Ln20sq+Xge+meWTAx/o721fdNnMLE8c+MCWdVdDfufmngwrMTOrHZ+hBCTxs698nHOXpucnODEzyxsHfuJrAzdlXYKZWU15SMfMrCAc+GZmBeHANzMrCAe+mVlBOPDNzArCgW9mVhAOfDOzgnDgm5kVhCIi6xquS9I4cOZDfvsG4D8pltMIfMz5V7TjBR/z/+vmiNi42BdWdeBXQ9LxiFi1E6fXgo85/4p2vOBjTpOHdMzMCsKBb2ZWEHkO/ANZF5ABH3P+Fe14wcecmtyO4ZuZ2bXy3MM3M7MKuQt8SfdJ+pek1yTtz7qeWpN0k6QXJA1LelnSw1nXVC+SmiWdkPSHrGupB0m9koYkvZr8vO/KuqZak/T95Pf6pKSnJOVu0mlJT0gak3SyYtt6Sc9LOpW8rkujrVwFvqRm4FfA54EdwDck7ci2qpqbBX4QER8D7gS+U4BjnvMwMJx1EXX0KPBcRGwHbifnxy6pH/geMBARO4Fm4OvZVlUTvwbuW7BtP3A0IrYBR5P1quUq8IFPAq9FxOmImAZ+CwxmXFNNRcRoRLyYLE9QDoH+bKuqPUlbgC8CB7OupR4kdQP3AI8DRMR0RJzPtqq6KAHtkkrAWuCtjOtJXUT8FXh3weZB4FCyfAi4P4228hb4/cDZivURChB+cyTdAuwCjmVbSV38EvghcCXrQupkKzAOPJkMYx2U1JF1UbUUEW8CPwfeAEaBCxHxp2yrqpsbImIUyp06YFMab5q3wNci2wpxGZKkTuD3wL6IuJh1PbUk6UvAWET8Peta6qgE3AE8FhG7gMuk9DF/tUrGrQeBW4HNQIekB7KtqrHlLfBHgMrZyLeQw4+AC0lqoRz2hyPi6azrqYO7gS9Lep3ysN1nJf0m25JqbgQYiYi5T29DlP8A5Nke4N8RMR4RM8DTwKcyrqle3pHUB5C8jqXxpnkL/L8B2yTdKqmV8gmeIxnXVFOSRHlcdzgifpF1PfUQET+KiC0RcQvln/FfIiLXPb+IeBs4K+m2ZNNu4JUMS6qHN4A7Ja1Nfs93k/MT1RWOAHuT5b3AM2m8aSmNN1ktImJW0neBP1I+o/9ERLyccVm1djfwIPBPSS8l234cEc9mWJPVxkPA4aQzcxr4Zsb11FREHJM0BLxI+Wq0E+TwrltJTwGfBjZIGgF+AjwC/E7Styj/4ftqKm35Tlszs2LI25COmZldhwPfzKwgHPhmZgXhwDczKwgHvplZQTjwzcwKwoFvZlYQDnwzs4L4H4XYpv+8pvsVAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.plot(ts, fs)"
]
},
{
"cell_type": "code",
"execution_count": 27,
"metadata": {},
"outputs": [],
"source": [
"fft = np.fft.fft(fs)\n",
"fftfreq = np.fft.fftfreq(len(ts), dt)"
]
},
{
"cell_type": "code",
"execution_count": 34,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"(-20, 20)"
]
},
"execution_count": 34,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAD4CAYAAADo30HgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deZhkdX3v8fe3qqv3hVl6FmaG6QEGZBEIDktUvERQiEtQlIRoInk04uM1ZLmP8cKD97mQxBvikuR6r5BgxJBgXJ5H0NFEEPACbogDgs4wM6GHaZhmVmbrvbvqnN/945zqrurppU5NdZ/Tcz6v5xmr+tQ5XT/K6vM5v/WYcw4REZGiTNwFEBGRZFEwiIhIGQWDiIiUUTCIiEgZBYOIiJSpi7sAlVq6dKnr6uqKuxgiIgvK008//apzrjPKMQsmGLq6uti0aVPcxRARWVDM7KWox6gpSUREyigYRESkjIJBRETKKBhERKSMgkFERMooGEREpIyCQUREyigYRGYxUhhh446NRF2i/oneJ9gzsGeOSiUydxQMIrP46uYHufVHt7Ll1R2Rjrvpkf/GX/3wS3NUKpG5o2CQVCn4BfJ+PtIxLx/uA2Dv0YFIx3nk2dcf8RjfI+9FK59IrSkYJFVu/+ntfOLxT0Q6xvM9APJ+IeK7OXznRzriS5u/xO/+++9GfB+R2lowayWJ1MKewT0MjEW8ig9P7p5f+UneOYeZw3de5PLtGVS/hMRLNQZJFd/5ka/i/TAQCn7lJ3kvDITI71VF+URqrSbBYGb3mNl+M9tcsm2xmT1sZi+Ej4tKXrvFzLrNbLuZXVWLMohUwvM9Ci5ak1DxJO9FuPoveNUFQ8EvRHofkblQqxrDPwNXT9p2M/Coc2498Gj4M2Z2NnA9cE54zJ1mlq1ROURm5Dt/vAZQqeKJuhDhuGJ/hE+0k7xqDJIENQkG59wTwKFJm68B7g2f3wu8q2T715xzo865nUA3cHEtyiEyG9/5ka/Ix2sMEZqSiiHiIp7kPd9TjUFiN5d9DMudc3sAwsdl4fZVwK6S/XrDbccwsxvNbJOZbTpw4MAcFlXSwnPRT7zFK/gofQyFYo0hajA4TzUGiV0cnc82xbYpp5Q65+52zm1wzm3o7Ix0ZzqRKfnOj3TlDxOjkqKcsIs1hqhNScVgiDrLWqSW5jIY9pnZSoDwcX+4vRdYU7LfamD3HJZDZFw1NYZikBSidD6Hxziq689QrUHiNJfBsBG4IXx+A/Dtku3Xm1mDma0D1gNPzWE5RMZV08dQnIvgeZUfl69yVFIxhBQMEqeaTHAzs68ClwNLzawX+J/AHcA3zOxDwMvAdQDOuS1m9g3geaAAfMw59bbJ/KimDb+4vxfhuGL4RK0xTLyXR45cpGNFaqUmweCcm24O/xXT7P8p4FO1eG+RKHznj3cMV6qqeQzjV/7RrnmKcyxUY5A4aeazpEpVM5/H5zFEb0qqtsagYJA4KRgkVarqY6B4so6wJEax87nKPgbNZZA4KRgkVTznRR6uOrFWUpSZz2FTkkYlyQKkYJBU8f3jmPk8H8NVVWOQBFAwSKoU5zFEmUA20e5fRTBUMfO59D1F4qBgkFSppnO32BwUpSnJ0wQ3WcAUDJIq1Zx4fRd96Gkh/P3HM49BJC4KBkmV8QXxItyTYXzmc4QaQ3GuhIu4VtL44nsRlwYXqSUFg6RKVTUGol/FF5fPiNrHoBqDJIGCQVJlYgntymsMroolMaptSlIfgySBgkFSpbo+hrDDOso9n8c7n6Mtn13N0FiRWlMwSKpU01RTzTH545zHoBqDxEnBIKnhnJs4yUe4+i/ebCdSLWO887i6piTVGCROCgZJjdITe3U1hih9DMUaQ3WzrFVjkDgpGCQ1qg0GV8UieoUq+xiKNQ3VGCROCgZJjdKTbaSmpPHmneiT4qI2Jel+DJIECgZJjdKTbZQTr3NV1BiO9w5uEVeAFaklBYOkRmmNIcrMZ0f0PoaJ+0NrVJIsPAoGSY1qawzFmc9RZjFXW2PQqCRJAgWDpEa1fQzuOOY+EKHz2Xf+eGe1agwSJwWDpMbxj0qKsoheWGOwCM1PpcGlGoPESMEgqVFaS4h29R+93X/ivaIvoxH1vURqTcEgqVFWY4jSlFTFPIZqm5KKVGOQOCkYJDWqbao5nqYkIjQllY6UUo1B4qRgkNQ47j6GCCOMqhkBVXpzHtUYJE4KBkmN0pNtNRPcogxX9aoIobLy6Q5uEiMFg6TG8fYxVNNhPfl9Z6JRSZIUCgZJjfnsY/CqCKFqJ+CJ1JqCQVLjePsYosxiriaESm83qmCQOCkYJDXmtcbgV9H5rOGqkhAKBkmNslE/Vc1jqGbZ7Qg1Bg1XlYRQMEhqVDsqiWqakqqYxazhqpIUCgZJjdITdGl7/myOb+Zz5bWT6oNLpLYUDJIaVc9jqKrzOXofQ7Wrv4rUmoJBUqP6zt3oS2H7RO9j0CJ6khQKBkmN6kclRb/pzvHOfFYfg8RJwSCpUe3MZ+ZpVJL6GCQp6ub6DcysB+gnWJi+4JzbYGaLga8DXUAP8NvOucNzXRZJt+onuLnwMUIwlIwwKniVvZfmMUhSzFeN4Teccxc45zaEP98MPOqcWw88Gv4sMqeqvVEPVdzzuXQl1kKFtRPNfJakiKsp6Rrg3vD5vcC7YiqHpEjVaxFZcdnt6hbRqzQYVGOQpJiPYHDA983saTO7Mdy23Dm3ByB8XDYP5ZCU88pO1pXNYyg9WbsId2MrXRKj0mBQH4MkxZz3MQBvcM7tNrNlwMNmtq3SA8MguRHglFNOmavySUpUU2MoPVm7KBPcqmhKqr6pS6S25rzG4JzbHT7uBx4ALgb2mdlKgPBx/zTH3u2c2+Cc29DZ2TnXRZUTXDU1htKTdaTO55L3yldTY9CNeiRGcxoMZtZiZm3F58Bbgc3ARuCGcLcbgG/PZTlEoMrbbZY1JVV3a8+KQ0jzGCQh5ropaTnwgJkV3+vfnHMPmtnPgW+Y2YeAl4Hr5rgcIlWdeMuakiL0MVS1VpJmPktCzGkwOOdeBM6fYvtB4Iq5fG+RyaoZ9VNeY4jSx1DSlFThPAbVGCQpNPNZUqOaRerKO58rrzG4shCK3tGtGoPEScEgqVHN/Q5q0cegez7LQqNgkNQohkHWslWteBopGEr2rXRUUrGTOkr5ROaCgkFSo3gVnsvkIl/FO2eRh6s6Z0D02kkuk1ONQWKlYJDUKJ6gm5zP6KHuio4Zvw+zy0KEYHDOD4+BglfZcNVRbzQoX12TagwSKwWDpEbxKrx1dJChXT+LdAwuG70pqRgMFZ7kh/JDALTkWjTBTWKlYJDUKF6Ft/k+w5nKvvrFY5yri9z57FwwGtyr8CQ/VBiiMdtILptTjUFipWCQ1Che/bf5jsGMVXaMP1FjiNSUVGWNoTnXTNay6mOQWCkYJDWKV+Etvs+QVfbVn+hjqIs08znoYwhqDJX2MQwWBmmua8bMVGOQWCkYJDXG+xh8n6EKawwDQwcByHo5otQYfHyc1wTAYH6womNUY5CkUDBIani+h2G0+o7BCmsM/d0PA7DOHY1cY3B+IwAD+f6KjhkqDNGSayFjGdUYJFYKBkkN3/lkLUuz8xnOWEVLXAwUgqv9Bq+6PgbnNTCYH6jomKH8EM11qjFI/BQMkhq+88lYhhbfUTAjXxid9Zj+sBmo0asD8yMs1+0FweA3MlSIEAy5ZtUYJHYKBkmNoMaQoSk8uQ8OTXl/qDJ9hWBuQVOhAZiYazCbxsJeznD7cV4jg5UGQ2GIpromspaNtGCfSK0pGCQ18n6eukwdzX5w0h0aenXWY/rzwzT4PvVePQADFTYLeQyxgZ6gxlBpU1LYx5DL5Mj7+YqOEZkLCgZJjYH8AG31rbSEcxOGhg9VdozvaPCDUUyV1hgGzWjxHfiNDBUqG5U0mA+Gq7ZmcvSPHKnoGJG5MNd3cBNJjL6xPtpyrTSHzTSDwwdnPaa/MESr79MaLrpXSY0h7+UZzWRoccGQ1aHCgYqOKfgFmnPNtG1/kP6mxlmPEZkrCgZJjb7RMBiKTUkjh2c9ZsAboc33aQ/7JSoJhsHRo0AwX8JRWefzUGFinaQ236evwiU7ROaCvn2SGv35ftpzLbSEJ/mh8AQ+kwF/LAiGsM2/kqakwcF9ADT7Duc1MuwNztqZXJwE1+x5tIdrOamfQeKiYJDU6B/rp62ueaLGUEEw9LsCrb5PR3gTnUpqDAODwWinVt8HvxHPFRjxRmY8phg4TcNHaAv7QAZGK5sYJ1JrCgZJjSAYWmgudj5XcOIdwKPNOU4Kg6GS5S0Gw9FOLWGNAWBgbOZAGW9KGjxIe1i+/rDmITLfFAySCgW/wGB+kPa6JlrDpqTDoxX0MQBtns8iF1zxVxQMYad2i/NxfrBeUv8sy2IcGQ1GIbX176fNC2o0/f27Z30vkbmgYJBUKJ7Q2+qaaHCwvFBg1+DMJ96hkaOMZIwO36eNPPi5yjqfw07tVt/Hec0AHJxlBNTLfS8DcIpv401JfaoxSEwUDJIKfWN9ALRlgyv4rnyBl0Znniuw46XHAVg7VqDJxsA1MDg2e41hIPy9Lb7DH+0EYOfRnTMe09PXQ1t9G4vcRDCoKUniomCQVOgfC5py2uqCNv+1LstOxmYcLdS9J7j959jYSpoZBb+BwQomqw2OBiHU4vu4QgcNmWa6j8x8j+mX+l6iq70Lc/5EH0MF8yxE5oKCQVJhPBgywZpHaxs76c8YR45MfyXffXAbDb5j/+g6mhgFv7GiGsNg+F71vgHG8qa1FQXD2va14BfGy9g/MvvMbJG5oGCQVCg2JbWHNYauji4AXur96bTH7BjcTZfLMEArjTaK8xsqG66aH6DZ98lYsIzG8sa17DiyY9r9hwvD7BncMx4Mzc1LyDhHXwXDaUXmgoJBUqE46qc9GwTD2WsuA+CpnQ9OuX++MMJmr5/T6xYz7BrI4jC/vqJRSUOFYVp8nxwe4FjReCqHRg7xUt9LU+6/ae8mAM5Zcg74Bay+hXbfcUTBIDFRMEgq9BztoSHbwLJcGwBLl5zBeX6OHxz81ZT7P/nsPRzNGBevfDPDBCurZgtN7BuavUN471gfS7ygnyCHx5ntvw7AQz0PTbn/D3b9gOa6Zi5ZeQk4HzJ1nOIy9Iyqj0HioWCQVOg+0s2pHaeSLd5oJ1PHm5eez5aMR/eOh4/Z/4FtX6fN9zn7rA8yTNDmnxvr4NDIIY7OciXf4w2wLh8sZ5GjQEeukwuXXch3dnyHQjhRrmgwP8ijLz3KZasvoz5bD34BMlnWNyzmhcKA7ssgsVAwSCp0H+5m/aL1wYkXwLJc+4ZP0uY77vjRrRRK7ub2k2fu5mHvEL/b/hqsroMhFwRD42hQ2+jp65n2fUZG+9ltjlXhMkc5Cni+4wPnfICevh7u3XJv2f53Pnsnh0cP84GzPxBs8AuQqeP09nUcycDBvl21+QBEItDqqnLCOzp6lP3D+zn9pNMhXD6bTB2LFp/Gn626kr/Y8yg3/dvl/M4Zv83OQ1u5a99PWOeMD1/9j2w75JEP/0waxsJgONrD+Z3nT/leL/f+FGdGs78IOEgOj+G8x7Vr3sxb1r6Fv3/m7zk4cpBLV17KE71P8PXtX+e9Z7yX8zrPC35BMRiWnQ+Hn6b7pcdZet7vz/VHJFJGwSAnvC2vbgEIagwD4aiiTBaA69769+S/91/533uf4Efb7wHgEmvkjrffS2PLUna/uKckGJqpy9TNWGPoCTuSvfwqoJumrMfuI8OYGXdcdgdt9W3c9/x9/Ovz/0rWsrzvNe/jzy/684lf4HtgWdavvRy238PmV37EpQoGmWcKBjnhPfzywzTVNbFh+QboeyTYmJn46r/vN+/kt/peYUfPYyxuX8OatZdBONS09/DQeDDU41jRvIath7ZO+17b9v+SjHO8MrgWMrCqLUvv4eHg+Gw9t7/+dm76tZvo7e9lTdsaljQtKf8FYY1hyYoLOLcAj+zfxB/W8LMQqYT6GOSENuaN8ehLj3L5mstprGuc6GPIlF8Ttbav4vzz3s+arjeNhwJA7+FhcvXBqKQcBU5tuYBNezcxXBie8v1+eGQb57kcRwrtAKxur6P3cPk9HJY2LeWCZRccGwoQ1BgydWDGVUsuYAtjvLR7U7X/+SJVUTDICe2r277K4dHDXLv+2mDDNMEwnVcOD9PR2gJAzgqsqLuQUW+UJ3c/ecy+e/c+x7aMx+tazh6vZZzcVscrR6YOkSmFNQaAt130xzT4jjsfv6Xy40VqQMEgJ6zth7bzhWe/wBtXvZFLV14abBzvfM5W9Dt6Dw+zuK0VgKaMj42dRkdDB/d333/Mvg/87DMArFjyLvIEv39la5Z9faOMFrzKCl0SDMtWXcQNrafzH2N7efCnn67seJEaUDDICemxXY/xwYc+SHt9O7e//vaJF8ZrDLMHg+87eg8PsbQjqDEsb8mw90ie97/m/Ty26zG2Hdo2vm9//x7uO/Qsl1srR+ouHK8xrGgN3qfYzzD7m3plZfvI2+/hQi/LLdv/hXsf/Bj5/Mx3ghOphdiCwcyuNrPtZtZtZjfHVQ45cRwYOsADLzzA+//9/dz0g5tY3rKcL1/9ZZY1L5vYKUJT0sbndjM45nHmqsUALGvJ0HNwkPed9T4WNSzitp/cNn5f5s989wYGDD560cfpeXWQ+vpg7sMZnQ2Ywb/97OXK/iPCCW5F9c2L+b/Xfos30sxn9z3Bb913MV/+j4/wUu+TmvwmcyaWUUlmlgW+ALwF6AV+bmYbnXPPx1EeSS7nHCPeCCOF4N9AfoBDI4c4OHyQgyMH2T+0n51Hd9J9pJtXBl4BYF3HOm65+BauO+M6ctncpF84MY+hlOc7jg7nOTI0xv7+UZ588SBffOJFXruqgze9Zjk8CueuaObvnunju88e4ZaLb+UTP/w4d/zsDs4cOMQDY3v4UOuZ9Ddczv3PPMVHT10Ku4JRSb+zYQ1f/vFORgsev3HmMrqWtrCsrYHWhjqspKMbKGtKKmo7qYvPf+BJfvyzz3HX1q/wtwd+wt8++hNO8h1nZVvoalzG8ublrGhfTWfbKbS2LKWtuZOWluW0tawkV99U0/9P5MRncVx1mNmvA7c5564Kf74FwDn319Md07SuyZ1+2+nzVEJZiDLOOC3fwSmujVw2w+RzLsDygW3UuTH2NZ/BrtkXSmVpSz2tDLF0+EXGss38Mr96/LVnF5evm3TBoU6KlfBTG/tZPLaHwdwSDjWsYl/f7E1Ana0NrB0K1m56pe28Y153DgqeT787xNZcP31ZtQTL7Db/weannXMbohwT1zyGVUDpXP9e4JLJO5nZjcCNAI1djfNTMllQ6soubBw7c4fZyQz3cm5vBBqBvbAowhvVt4dP9o5vyjrwStJn86L9E88BmtqBPNBT+Xt1FN+nZ+rXx1uZrOy/vTBVCopUKa5gmOpbfEzVxTl3N3A3wOted6F78vofH3PIMRWecIOj/IQx6WWOPXDKIoCb4j2Y4j1Kn7opNhZ/muJ9XXFhtwrLc+yvmaIcbvJ7uyl/n/On/xymq0sWa5l2zHuE/+v8oBPVeZjzcV4heO77QVOOH/wcPA/3873x18wbw/LDWH4YCsGjFYaCx7EBGDxAZij4Vzf0KuaCfgM/k+PI4vPZsfxqfnHSW9hxNMPTLx+me/8A56/u4F8+dAkdv/gH+P4nOWfkS7zrkjP51LtfO81/ZYnBV+Ezp8HbPgsXfzjY1r8X/vG/4LI5+L37sc4zyo/ZvxXuvBTe+2U499pZ3+JrT73Mzff/iu2L/oyGs67Cf8fn+dOvP8vG53azsqORDV2LWb+slaWtDSxuqae5PktDXYaGXPBYX5cha0bGjEyG4NGMjIGFj8VtFr5ucEytyib9ac6UN1GPnfyrJjejHfv6zPtLZewPon9ucQVDL7Cm5OfVwIx3ZjfLkGtQW6lM4hXgcA8c2Eam9ykWv/AIi7f8FRc1fB7e/EnctR/moef38dGvPMM/Pr6DT7QGIeKR4fWnLa3sPYpt/t5Y8OgcbPxjGDmKffhRmBwKANlgUhyTVlOdzq+fFkx2c14eMnU8uGUvG5/bzZ9csZ4/vmI92YxOijJ/4mqk/Dmw3szWmVk9cD2wMaayyEKWrYOlp8NZ74C3/AV89Mfwhz+A1RfB9z6Bfecmrj5nBe8472Tu/UkPY2PByd0jy8qTKmyeLJ7ki8HQ/Qi88BBc8T9g+TlTHzM5TGaxoiMoiwtnPt/12A5O62xRKEgsYgkG51wB+CPgIWAr8A3n3JY4yiInGDNY/Tr4vW/CZR+HX9wHP/o73v7alQyOeRweDDqBC2Q4uaPCGuh4MIRX/098FtpXw0UfruCYyoKhoS7L0tZ6zC/gk+H5PX1cfe4KhYLEIrZhDc65/3DOneGcO80596m4yiEnKDN48yfhrHfC459mXS7okO4fGsFhZDNZOtsaKvtdxXkF3hgcehF2PQmXfATq6qc/ZnKYVODkk5ow5zFYCIbPrlnUXPGxIrWk8W5y4jKDq/4X+HnW7rgPgIGhETyyLG9vrPxq3Cw40XtjsD28R/RZ75z5mGy0piSAlR2NmPPoGw268VctUp+axEPBICe2k06B06+kYdu3aK03BodH8chycqX9C0XF227+5/eg8yxYvG72/SFiMDSRcR59Y0EwrFaNQWKiYJAT37nvxfp6ubKtl6HRETyMFZX2LxRl6oKT/J7nYO3rZ99/ulFJe56D+z8C/3QlfOtjcGD7+EsndzRQZz4HBoLZ2Ss7NHdH4qFgkBNf1xsBuKihh6GRMQpkac5VtrrquGw99O2GkaPQeebs+2eygJXXGJ7fCF+8Iqh11DXC89+Guy+HnT8EoDlsfXp12GNZWwONUcsoUiMKBjnxta2A5iWcycuMjI7hkSUT9ZufrYd9m4PnS9dXfkxhNHh+9BX41kfh5AvgT56DP/gu3PQ0dKyBb/4hjA6QIZjoeGTYZ+VJ6l+Q+CgY5MRnBsvPYeVIN+YKeGTIRJ1Fm60LJtIBLJ1iQttUco0TwfDYXwfNSu/5J2gK18doWw7XfAEG9sKTd5GzIBhGfKMppz9NiY++fZIOy89l2ciL5FyeAtkqgiHsM8i1QPuqyo7JNUN+CPIjQbPRue+BRV3l+6y5CLoug83fJEvQt5D3M5q/ILFSMEg6dJ5Jzh/lZA7gUcWJtxgMi9bOvIBQqWIwdD8Mo33w2vdOvd9r3g4HttIxuBOAMVdFjUakhhQMkg7NwbpIS+wonpt6Se4ZFZe4aDyp8mNyzZAfhpefDDqbu9409X5nXA3AygM/AiDvTMEgsVIwSDqE7fpL6KNAlmy1TUmN7TPvVyrXFNQY9m0JRjJlp1mzclEX1LdyUn83AGN+BrUkSZwUDJIOzcHtORfZAAUyZKptSmqIEAz1YY1h//OwbJrF9iBomlpyGm39OwDIO/UxSLwUDJIOTRN3yvHIRm9KKl7tN7RVfkyuGY72wsA+WH72zPsuOZ22gRcBGPNN9x6QWCkYJB1KgiHvMtGbkjLhvaOjNiX1BfehZtlZM++7ZGJuxKhfRflEakjBIOlQ10A+E0wa86oZrlq8y16UpqRcyVpHsw1xXTJxP/OqJuCJ1JC+fpIaI7kOAF5xS6L3MeSHg8dINYaSYJhtNFPHRHAUqpmAJ1JDCgZJjUI2WJRuhzs5+qifQhgMDR2VH5MrWdaiaZZgqG8df1pVjUakhhQMkhqN+T4Advirorfh50fCX1JFjSFTVx4SU2mYCIZCNRPwRGpIwSCp0ZQ/BIQ1hshNSUPBY9ThqgDZCu4UVz8x2smnigl4IjWkYJDUedGtiH7iLYQ1hkjDVcNaQjY3+75lNYYqJuCJ1NA0UzFFTjzf3/BFXnjyOwzTOL9NSdkZ7g1dVDdxUx7PqY9B4qVgkNQ40HkJnym0AEQ/8Y53PlcRDHUVNCWVlKeqmdkiNaSmJEmN0jCIfOIt3s4z6sxnqKwpqUQwKinSISI1pRqDpEbpyTbyifd37guWt8hEuN3meB9DBTWGEprHIHFTjUFSo/RkG3k4aEPb7MtaTFYMhroK+hhK+BquKjFTMEhqlAbDvCxS5xeCx0o6n0sUqlnkT6SGFAySGqVX4fNyQd66PHhcf1Wkwzy0iJ7ES30Mkhql59p5OfEuXgd/urnie0Q7y2DOD+5JraYkiZFqDJIa5TWGeTrxnrSGSpdK9cO1nLRWksRNwSCpcVzDVeeBCye5eU639pR4KRgkNcqCIYEn3mIwZM3XqCSJlYJBUqP0XJvEE+9g11sBGHb1urWnxEqdz5IapWGQxBPvoctu551Pn89BOjQqSWKlGoOkRtKbkjLZHLtcMMQ1ieWT9FAwSGqUdjgn8Yq8bNSUkkFipGCQ1Cg91yaxKam0SBquKnFSMEhqZI9nraR5UFqmrP4yJUb6+klqWNL7GCyGCXgiU1AwSGokvQ1fwSBJMWfBYGa3mdkrZvZs+O9tJa/dYmbdZrbdzKKtMCZSpUzC2/CP634RIjU01/MY/s4599nSDWZ2NnA9cA5wMvCImZ3hnPPmuCyScpn5Xl01ouO6X4RIDcXRlHQN8DXn3KhzbifQDVwcQzkkZcpOvEmsMSR8Ap6kx1wHwx+Z2S/N7B4zWxRuWwXsKtmnN9x2DDO70cw2mdmmAwcOzHFR5USXne8b9USU9CU7JD2OKxjM7BEz2zzFv2uAu4DTgAuAPcDniodN8avcVL/fOXe3c26Dc25DZ2fn8RRVpPx+DAk88c77jYREpnFcfQzOuSsr2c/Mvgh8N/yxF1hT8vJqYPfxlEOkEkk/8WpUkiTFXI5KWlny47uBzeHzjcD1ZtZgZuuA9cBTc1UOkaKk349BwSBJMZejkj5tZhcQNBP1AB8BcM5tMbNvAM8DBeBjGpEk86F0NnEST7zqY5CkmLNgcM79/gyvfQr41Fy9t8hUFtLM5wTmlqSIZj5LamQT3lRTtvprEpNLUkPBIKmxENrwi3mQ1PJJOigYJDWSPlwVJsqV0OJJSigYJDWSPlwVJvpBVGOQOCkYJDWSPlwVJvpBFAwSJwWDpEYm4cNVYaImk9SmLkkHBYOkRibhw1VhoowJzRwhYkQAAAgBSURBVC1JCQWDpEbSh6vCRBOXagwSJwWDpMZC6GPQcFVJAgWDpEZpH0MS78cApcNVk1k+SQcFg6TGQuhjmBiuGnNBJNUUDJIaZfMYEnrmLdZk1McgcVIwSGqUts4ktammmAdJvMOcpIeCQVIju4CaklRjkDgpGCQ1FsKoJK2VJEmgYJDUyGQWwDwGDVeVBFAwSKqMLzmR0BNvRsNVJQEUDJIqxaaapJ53M+pjkARQMEiqJL1zN6t5DJIACgZJlaQva10sVlI7xyUdFAySKhOdu/GWYzqZhAeXpIOCQVIlkzHMkjuBrNjEldTOcUkHBYOkSsYs0VfjEzOf4y2HpJuCQVIlm7FEX43rfgySBAoGSZWMJftqXH0MkgQKBkmVjFmir8bHR03pL1NipK+fpErS+xjGh6smuIxy4lMwSKpkw1FJSTU+8znJhZQTnoJBUsUs2R27urWnJIGCQVJloTQlmf4yJUb6+kmqZDPJDgZNcJMkUDBIqpgldzkM0HBVSQYFg6RKNuHDVTMarioJoK+fpErS+xh0BzdJAgWDpEpGw1VFZqVgkFTJLJDhqsoFiZOCQVIl6aOSip3jSV0WXNLhuILBzK4zsy1m5pvZhkmv3WJm3Wa23cyuKtn+OjP7Vfja501/ATKPzCzRo5KSHlySDsdbY9gMXAs8UbrRzM4GrgfOAa4G7jSzbPjyXcCNwPrw39XHWQaRimUt2R27GTPd1lNid1zB4Jzb6pzbPsVL1wBfc86NOud2At3AxWa2Emh3zv3UOeeAfwHedTxlEIki6aurZhJeo5F0mKs+hlXArpKfe8Ntq8Lnk7dPycxuNLNNZrbpwIEDc1JQSZdgVFJyz7wZ04gkiV/dbDuY2SPAiileutU59+3pDptim5th+5Scc3cDdwNs2LBh2v1EKpVZADOfk9zUJekwazA4566s4vf2AmtKfl4N7A63r55iu8i8yGYS3pSUUR+DxG+umpI2AtebWYOZrSPoZH7KObcH6DezS8PRSB8Apqt1iNRcxpLdlJTNJLtGI+kwa41hJmb2buD/AJ3Av5vZs865q5xzW8zsG8DzQAH4mHPOCw/7KPDPQBPwvfCfyLz44BvWMZL3Zt8xJu+5cDVnrWyPuxiSchYMDkq+DRs2uE2bNsVdDBGRBcXMnnbObZh9zwma+SwiImUUDCIiUkbBICIiZRQMIiJSRsEgIiJlFAwiIlJGwSAiImUUDCIiUmbBTHAzs35gqiW+k2Yp8GrchZjFQigjqJy1pnLW1kIp55nOubYoBxzXkhjzbHvU2XtxMLNNSS/nQigjqJy1pnLW1kIqZ9Rj1JQkIiJlFAwiIlJmIQXD3XEXoEILoZwLoYygctaayllbJ2w5F0zns4iIzI+FVGMQEZF5oGAQEZEyiQ4GM/uMmW0zs1+a2QNmdlLJa7eYWbeZbTezq2Iu53VmtsXMfDPbULK9y8yGzezZ8N8/JLGc4WuJ+TxLmdltZvZKyWf4trjLVMrMrg4/s24zuznu8kzHzHrM7FfhZ5iYO16Z2T1mtt/MNpdsW2xmD5vZC+HjogSWMXHfSzNbY2b/z8y2hn/nfxJuj/55OucS+w94K1AXPv8b4G/C52cDzwENwDpgB5CNsZxnAWcCjwEbSrZ3AZvj/hwrKGeiPs9JZb4N+Hjc5ZimbNnwszoVqA8/w7PjLtc0Ze0BlsZdjinK9SbgwtK/E+DTwM3h85uLf/cJK2PivpfASuDC8Hkb8J/h33bkzzPRNQbn3Pedc4XwxyeB1eHza4CvOedGnXM7gW7g4jjKCOCc2+qcS/ys7BnKmajPcwG5GOh2zr3onBsDvkbwWUqFnHNPAIcmbb4GuDd8fi/wrnkt1CTTlDFxnHN7nHPPhM/7ga3AKqr4PBMdDJN8EPhe+HwVsKvktd5wWxKtM7NfmNnjZnZZ3IWZRtI/zz8KmxPvibtZYZKkf26lHPB9M3vazG6MuzCzWO6c2wPByQ5YFnN5ppPU7yVm1gX8GvAzqvg8Y18Sw8weAVZM8dKtzrlvh/vcChSArxQPm2L/OR13W0k5p7AHOMU5d9DMXgd8y8zOcc71Jayc8/55lr35DGUG7gL+MizPXwKfI7hISIJYP7eI3uCc221my4CHzWxbeCUs1Uns99LMWoFvAn/qnOszm+prOrPYg8E5d+VMr5vZDcA7gCtc2EhGcGW2pmS31cDuuSlhYLZyTnPMKDAaPn/azHYAZwBz1vlXTTmJ4fMsVWmZzeyLwHfnuDhRxPq5ReGc2x0+7jezBwiawZIaDPvMbKVzbo+ZrQT2x12gyZxz+4rPk/S9NLMcQSh8xTl3f7g58ueZ6KYkM7sa+O/Abznnhkpe2ghcb2YNZrYOWA88FUcZZ2JmnWaWDZ+fSlDOF+Mt1ZQS+3mGX+SidwObp9s3Bj8H1pvZOjOrB64n+CwTxcxazKyt+JxgUEeSPsfJNgI3hM9vAKar6cYmid9LC6oGXwK2Ouf+tuSl6J9n3D3ps/SydxO04T4b/vuHktduJRgRsh34zZjL+W6Cq8dRYB/wULj9PcAWgtEqzwDvTGI5k/Z5TirzvwK/An4ZfsFXxl2mSeV7G8Hojx0EzXWxl2mKMp4afgefC7+PiSkn8FWCJtd8+N38ELAEeBR4IXxcnMAyJu57CbyRoGnrlyXnzLdV83lqSQwRESmT6KYkERGZfwoGEREpo2AQEZEyCgYRESmjYBARkTIKBhERKaNgEBGRMv8fGXR0643O7qUAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.plot(fftfreq, fft.real, fftfreq, fft.imag)\n",
"plt.plot(fftfreq, np.abs(fft))\n",
"plt.xlim(-20, 20)"
]
},
{
"cell_type": "code",
"execution_count": 35,
"metadata": {},
"outputs": [],
"source": [
"ifft = np.fft.ifft(fft)\n",
"ifftt = np.fft.fftfreq(len(fft), fftfreq[1] - fftfreq[0])"
]
},
{
"cell_type": "code",
"execution_count": 38,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"(0, 10)"
]
},
"execution_count": 38,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAD7CAYAAABnoJM0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO29e5QkZ3nm+bx5v1VmXburr2o1tCSEsLi0kWw8gI3lFcJr4T3MWNjYGtZerQ4Ge4zPGWvsszs+ts8Ma+9hYY+xtRrMwIwZYxaLQTAahJHFYmMDat0QUqtR05K6q7u66155j8jIfPePiC8yMjMiMqsrIzMr8/2do9OVWVHxfZ3q+p5478TMEARBECaX0LA3IAiCIAwXEQJBEIQJR4RAEARhwhEhEARBmHBECARBECYcEQJBEIQJpy9CQES3E9EZIjpLRPe5fD9HRF8momeI6Dkien8/1hUEQRB2D+22joCIwgB+AOA2AEsAHgfwXmZ+3nHN7wLIMfPvENECgDMAFplZ39XigiAIwq6J9OEebwZwlpnPAQARfQ7AnQCed1zDAKaIiABkAGwAMLrdeH5+no8dO9aHLQqCIEwGTzzxxBozL+zkZ/ohBIcAXHC8XgJwS9s1fwrgIQCXAEwB+AVmbnS78bFjx3Dq1Kk+bFEQBGEyIKJXdvoz/YgRkMt77f6m/wHA0wAOAng9gD8loqzrzYjuIaJTRHRqdXW1D9sTBEEQ/OiHECwBOOJ4fRjmk7+T9wN4kE3OAngJwA1uN2PmB5j5JDOfXFjYkXUjCIIgXAX9EILHAZwgomuJKAbgLphuICfnAbwDAIhoP4DrAZzrw9qCIAjCLtl1jICZDSL6IIBHAIQBfIqZnyOie63v3w/gDwF8moiehelK+h1mXtvt2oIgCMLu6UewGMz8MICH29673/H1JQA/04+1BEEQhP4ilcWCIAgTjgiBIAjChLPnhOBbZ9fwzIWtYW9DEARhbOhLjGBQ5Ks1/NInvwMAeOnf3wGzUFkQBEHYDXvKInj+Ut7+erWoDXEngiAI48OeEoIfrhbtr19ZLw9xJ4IgCOPDnhKCpc2K/fV5EQJBEIS+sKeEYCWvYSpuhjXWS+IaEgRB6Ad7SwgKVRxfSCMWCWG9KKMMBEEQ+sGeEoK1oo6FqTjm0zGsl0QIBEEQ+sGeEoJ8pYZcMoaZdAwbIgSCIAh9Yc8JQTYZwVQigmK164AzQRAEoQf2jBDUG4yCZiCbiCITj6KgiRAIgiD0gz0jBIVqDQCQS0ZNi0CrDXlHgiAI48GeEYJ8xbQAsskoMnFxDQmCIPSLvSMElkWQTUSQSURQ0upD3pEgCMJ4sHeEoNJ0DWXiEej1BjRDxEAQBGG37Bkh2LaEIGvFCACIe0gQBKEP7BkhsF1DlkUAAMUBZQ5Va3WUJEtJEIQxpS9CQES3E9EZIjpLRPd5XPN2InqaiJ4jov9vp2vYweJEBGlLCAoDsgh+94vP4qbffwQr+epA1hMEQRgkuxYCIgoD+ASAdwK4EcB7iejGtmumAfwZgJ9j5tcC+Oc7XUfVDaRjEbvx3CAsgkaD8eCTF8EMfOuHa4GvJwiCMGj6YRG8GcBZZj7HzDqAzwG4s+2aXwTwIDOfBwBmXtnpIhXdQDIaRihEyAwwRrDssAKeu5j3uVIQBGFv0g8hOATgguP1kvWek+sAzBDRN4joCSL6lZ0uUtbrSMXCAGDHCAoDKCp7Za1kf31hU2YgCIIwfvRjZrHb4GB2WedNAN4BIAngn4jo28z8g46bEd0D4B4AOHr0qP1+Ra8jaQlBKhax3mvsfvdduLRtWgTH59O4sFHpcrUgCMLeox8WwRKAI47XhwFccrnmq8xcYuY1AN8EcLPbzZj5AWY+ycwnFxYW7PedFkEyav5ZrQVfR7BmzUZ+7aEcVgoyDEcQhPGjH0LwOIATRHQtEcUA3AXgobZrvgTgnxFRhIhSAG4BcHoni5RrdSQtSyAeNbddGYQQFDQkoiEcmUlis6yDud3Ycee3/vppfPjzTwe8O0EQhN2zayFgZgPABwE8AvNw/zwzP0dE9xLRvdY1pwF8FcD3AHwXwCeZ+fs7WaeiG0hZlkA8EgLRYCyC9ZKO+Uwcs+kY6g2201j9uLBRxhefuogHn7yIS1viThIEYbTpR4wAzPwwgIfb3ru/7fWfAPiTq12jrNexmI0CAIgIyWgYFT14Idiu1DCdimIuEwMAbJR15FJR35957lIzu+h7S9s4OJ0MdI+CIAi7Yc9UFjuDxYAZJ6gOoNdQvlLDVDyKmZQlBKXucYKLDivg/EbJ50pBEIThs2eEwBksBoBENDyQrKFC1UA2GcFcOg4A2Ch1T1m9uFlBKhZGLhnF+Q1JORUEYbTpi2toEJR1w04bBYBENDSQGEGhWsNUIoqZtOkO6sUiWNos49B0Eg1mma0sCMLIs2eEoFJrcw3FwgPJGipUDUwlmhbBeg8H+8WtCg7NJFHW6lgvihAIgjDa7AnXUK3eQK3OdtYQYMUIAhYC55zkZCyMeCSErXJ319DydhUHp5OYTcewWRYhEARhtNkTFkHZyg5KtsUIgm46p+6v5h9MJaJd16w3GJtlHfPpGJiBjVdECARBGG32hBCoNNHWGEEYqwFX+hbs8ZhmfGAqEena+jpfqYEZmE7F0GBgs1xDo8EIhdw6cQiCIAyfPSEEZd08fFPt6aMBu4bsGQhJ82PKxCMoVv1dQ8oVNJOOgmFaCPlqDdNW+qkgCMKosUeEoNM1ZApBsOmjyiKYclgE3VxDm1YMYToVA1n9+DZKugiBIAgjy54IFqvsoNY6glDgWUPKDaRiBJl4d9fQdsWyCFIxzKRVEdpg4gTVWh2nl2VmgiAIO2NPCEFZdxGCAaSPliyXlBqNmekhRrBpFZzNpKLIJU1LIt/FndQvPvDZJ/HOj/89nr6wNZD1BEEYD/aEEFSsAzkZbXqyktEwdKOBeqO3bqBXQ7sATcV7cQ2ZT//TqRgycfPnilrw9Q5GvYG/e8Ec/Pa15y4Hvp4gCOPDnhACV4tgADMJ7HWjremjfq2ot8o1hEOEbCKCTNy0CAYxUtPZyuLZi9uBrycIwviwZ4UgHjG3rhvBBYxtS0SNyExEUG+wr0tqq6Ijl4yCiJC2LIJSwPUOAPDKuikEB3MJXJD+RoIg7IA9IQQVl6yhmBKCeoBCUKsjEiJ7LTUr2e8Jv1A1kLWCy+mYmq08OIvg1lfN4dJ2FY0eXWYvXM7j1MsbQW5NEIQRZ08IQdmloCweMUVBCzCFtKzX7bGYQDN7KN9FCDLWdaEQIR0LD8QiWC1oCBHwI4dy0I1GTz2RdKOB2z/293jP/f+EK/lq4HsUBGE02RtCUDMQi4QQdlTnNi2C4GIE7TMQlBD4BYyLVcO2HADTnTSIGMF6ScdsOo7DMykArTMRvPjeUjO76B9/uBbY3gRBGG32hBBU2mYRAM0YQZBFZe0zEHoJ/uatttWKdDyCoj4AIShqmM/E7GlovYzIdAaVn78k9QeCMKn0RQiI6HYiOkNEZ4noPp/rfpSI6kT0np3cv6zXWzqPAo5gcYAxgrJeR9LhjlIWQcGnLqCoGZhyWART8UFaBDEcyCUAmB1Qu/HD1SKmEhFcv38KL63JJDVBmFR2LQREFAbwCQDvBHAjgPcS0Y0e1/0fMIfc74h2Fw3QdA0FGSOo1Iw2i6B78LeoGbZgAJZFMIAYwXpRw1wmjlwyinCIsNlDjOCV9TKunU/j6FxKJqkJwgTTD4vgzQDOMvM5ZtYBfA7AnS7XfQjA3wBY2ekC7dPJAEewOMC5xV7BYq8nfGZuCRYDpngMIli8XtQxl44hFCLMpKI9BYuv5KtYzCawbyoeeCdXQRBGl34IwSEAFxyvl6z3bIjoEICfB3D/1SxQdrEIBlNH0LquajXh1WaiWjMrnVUsAeitP9Fu0Yw6CpqB+YzZ22gmFevJIriS17A/m8DCVByb5RpqAbrZBEEYXfohBG6N9tuT2D8G4HeYuevjOxHdQ0SniOjU6uoqADOf3ytYrAUpBG3rRsMhJKIhFDX3GEFBU91KW11DpYCDxXZ/I6vJ3Ww61rXRXbVWx3alhsVcAvMZawynjNUUhImkH0KwBOCI4/VhAJfarjkJ4HNE9DKA9wD4MyJ6t9vNmPkBZj7JzCcXFhYAdGbvAI700QCFwG3dTDzi2TuovVspYFZDq4K4oMi3DdCZTcew0WVEpqobUBYBAKwVxT0kCJNIP4TgcQAniOhaIooBuAvAQ84LmPlaZj7GzMcAfAHAB5j5v/a6QEWvtzScA5wxgoBdQ23r+vn8VezAWUeQjIWhGY2eK32vhubcBHPdXiyCy9tKCOK2RbA6ICE4c7mATzx2NvDBQoIg9MauB9Mws0FEH4SZDRQG8Clmfo6I7rW+f1VxASdmsNgrRhDMYcLMKOsGkrFWrfTLAmrOOG7GCFSwuVKr2zGGftO0RJoWwVZZR73BLUV4Tq5YweH92YS9x0EFjH/vi8/i1CubiIQI/+vbXjWQNQVB8KYvJxMzPwzg4bb3XAWAmf/lTu/v5xoKyiIwGowGoyVrCPAXAvVk3m4RAIMRgqzDImgwsF2pYTbtPhlNBZPn0jE7I2sQrqGybuApa17CN86sihAIwggw8pXF9QZDMxqeWUNBCYG6r3JBKfxcQ24xAtUuO8g4QbtFoAbibFe8C9+2rJGauWQUyVgYsUjI9/p+cXo5j3qDcWg6ie8tbfm29BYEYTCMvBC4jakEgEg4hBAFFyzWrHXj0daPaKdCkBzA3ATbErHWtSej+QlBRcdUPIJIOGT/jN/1/eLcqlnB/DOv3Y+SXu85LvHpb72ED/6XJwNNDhCESWXkhaBszwTodKvEI+HACsqaFsHOYwROF1DK4RoKikLVQIiAtLVWtgeLYLtcQy7VjGVkExHkK4OZmxAOEd7yqnkAwPn17hXN60UNv//l5/GV7y3jEZm+Jgh9Z+SFoGJPCQt3fC8WCQVnEXi6hsK+QpCMhhENNz/W5EBcQzVk4hEQmYHhXmYlb1VqmHYKQTI6kNnKl/NV7JuK4/hCGkBzoI4fz7R0SV0PbG+CMKmMvBC4TSdTxCOhAGMEdXsNJ+l4BNVaA4ZLFW6hanQEhBMDsgicmUqqnsDPItgs65hONgPJg3INrRQ07JuK49BMEkTAKz30OPr+xTyIgJsP5/DCZemSKgj9Zs8IQXuwGDD998HFCBr2Gk5URlDJ5Qm/qDWnkykGYRHkq62N7poxAm9Xz3a5zSJIRAcSLF7JV7EwFUc8EsbBXLKnsZrPXtzG8fk0XnMg25MrSRCEnTHyQlBxmU6miIWDtAi8s4YA9+E0xWqtpeEc0FpHEBRFrWZbAQCQiIYQDdMOXUMR38lr/WKtqGFhymyVvS/bW7O7sytFXL84hSOzKayX9IF0cxWESWLkhUAFi91dQ0EGi71dQ4D7QPqiZthzihXJAbiGSlodqXjz8yEi5JLeT/iNBmOrzTWUTZiuoSDTOY26OUJzn9XSYiHTXQgaDcbFrQoOz6RwzZw5fU2sAkHoLyMvBOoAdXMNxYKMEdR2bhG0t6AGHEIQoGvIrSmfOtjdKOoGGowWiyCXjMJocKCCtVbUwWxaAgCwMBXvmj66VtKgGw0cmk7imlkzwHx+YzBDdJgZKzLLWZgARl4Ihhcs9ogRJLwtgpLeOp0MGEwdQUWv24VriqyPRbDtKCZzXg/4B5h3i3r6X8g0hWCjpPu2v764aY7cPDSdxMFp06V0uYfpa/3g//r6i3jzv3sUX3xqaSDrCcKw2DtCEHWpI4iGA0wf9XANxXxcQy5ZQ9FwCJEQBfqk7WoRJKOePn9VVTydarqGmmM4g/O/b1XMthaq7YVqdufXIO+iNXv50EwS06kYQmRaFkFTbzD+0z+9DAD49LdeDnw9QRgmIy8EFbugzMU1NMRgcfuBycwoap2uIcC0Cip60F1SW/fplw6qDmSnayjt4/LqF00BMtdV7a/94gSXLCE4OJ1EOESYTcexXgq+J9KLKwVslWs4NJ3E9y/lpVOqMNaMvBCU9TqiYbKbzDkx00cDChbXvILF5oHbbhFoRgO1Orc0nFMkYmFUasEcsA3Lr98uBGalsIcQqAM52TpJDQDKHrMW+sGWtR/lhupFCFbyGpLRsJ2WO5+JYbUQvEVwetmsV/jFW46i3mA8vyz1C8L4sieEwC11FADig7AIvGIEbcHfktbZZ0gR5HAatc/2Fhwqa8gtC0gdyE7XkHItBWkRbFvDclRsYqGHOQirRQ0LU3G7ano+Ex9Il9SX18ogAt71ugMAgO9f3A58TUEYFiMvBCXNsHvotBOPBicEVStrKBZu/YjikTCiYeo4MO0+Qy6ilYyGA4sR2FlVbYKV9ckCaj+QAYdFEOBYze1KDalY2Ha39WIRrBU1exYzYFoEg3ANnd8o42AuiWvmUsglozhzuRD4moIwLEZeCMp6HSmPPv6xcJC9huqIhMjuzukkHY/Y08gUKmbgFiNIRMOo1ILZZ7POotMiANyzgLbKNaSt1tOKlE8QvF9slWst7qhENIypRMRXCFYLmi0YADCXiWNtAK6hl9dLODqbAhHhmrkULljZS4Iwjoy8EJR0P4sg2O6j7fEBRTrW2Yrank7mIlrJaBjVgFxDKoiZcKkjANzbTGyWay1uIcBZHxFsjCDrEALAKirzcw21CcF8Jo5KrR6o5QKYRWvH5s0CtiOzKZxfH0ztgiAMg5EXgrJWd80YAsxArm40AqmG1Yw64i4dTwEzDtDuGiq5tKBWJGPBuYbKHt1Z/SyC7Yre4hYCzLYUIQrYNdTW3wgA5jIxbHikg9bqDWyWa3aaKdDMOFIB7yCo1upYL+k4mEsCAI7OprC0WUE9wLnTgjBM+iIERHQ7EZ0horNEdJ/L93+JiL5n/fePRHRzr/c2LQJv11CDzbGS/Uar+VgE8QhKurtF4Jk+GlSMwKMpXzZp7sMtc2jL5UAmIqRj3rMW+sFWpbWtBQDM+aSDrlsC4bQIemmxvVuUq0pVQF8zm4LRYCxvi3tIGE92LQREFAbwCQDvBHAjgPcS0Y1tl70E4G3M/CMA/hDAA73e3zdGEOC4Sl/XUDzS4UKxp5N5WQQBuYaUwLRXFvvGCCqdQgCYf68g00fzFcMWKMVcJmYf+O2oA9lpEdgttgO0CJSrap/VHO/wjOkiWpI4gTCm9MMieDOAs8x8jpl1AJ8DcKfzAmb+R2betF5+G8DhXm/ulzWkhCCIgLFm1DuKyRSZeBjFtifSYVsEbr2GAPcn561yDblk51D7VDyMYoCuoaJmIBNvdw3FsVHWXd0uq0WzlYSbRTCQVhjWuos5888r0ndIGFP6IQSHAFxwvF6y3vPiVwH8915vXvGpIwhWCBodNQQKc25xZx1BiNBR2AWY/vegKlOb6aOt66p6hvYDk5mxXdEx42IR+M1j3i3MjJJuIBNv3edcOgZmc1BOOyo7aCEzWCFY6RACM1YwqB5HgjBo3E/YnUEu77k67YnoJ2EKwU943ozoHgD3AMDRo0cR1g27mrcdleMfiBB0ixFonemjace4SCcJyyJgZtfv7wavpnyRcAiZeOcc4rJeR63O7q6hWHCuobJeB3NnMH3OqhHYKOktLiCg6aIZhkVAZIoUYApkJh7BsgiBMKb0wyJYAnDE8fowgEvtFxHRjwD4JIA7mdlz8CwzP8DMJ5n55Pz8AhrsPpQGcFgE9f4fXv6uoQiKutGSrVTUOjuPKhLRMJgB3afL5tXilT4KwHUmgV1V7OIaSvvMY94tXllVc2nzkHerFl4taMjEIy3xj6lEBETuQfB+sVrQMJeOtdSQ7M/GxTUkjC39EILHAZwgomuJKAbgLgAPOS8goqMAHgTwy8z8g15vXLcOWi+LID7EYDFz82kccO882r7PagBFZXbWkItLaioR6YgRbKmqYq9gcUAxgqJHCw5VNewWMG6vKgaAUIgwFY8EahGsFzW7Q6riQC6JyyIEwpiyayFgZgPABwE8AuA0gM8z83NEdC8R3Wtd9r8DmAPwZ0T0NBGd6uXeDSuAOIoxAqC1Cter8yjQTO3UAogTlGtmU76oSwW0m0Ww7dJwTpGKdWZD9QuvFhxzljto3cUiWC92uosAU8SCFILtSq3DYtqfTUiMQBhb+hEjADM/DODhtvfud3z9awB+baf3bSiLwLOgzHx/8FlDzZbN+6z3iprh2nAOABIRNZwmGIugPXVUMZ2K4uW11rGObg3nFJl4OHCLoN1qmk5GESJg3WUmwVpRw/GFdMf7fmM4+8F2pWanjCoWc3GsFDTUG4xwqL9xHkEYNiNdWawyCrvVEQThe+8WLAbQkjnkKwRqSlkA7TDMrCoPIUjGOrJx1GvPOgK9blti/UR9Vu1tukMhwmw65jpsxnQNdVoE2USwQpCv1DoqrxdzSdQbPJDOp4IwaEZcCJRraAhZQ74xAnM/Ba15GBWr3hXQiaiKEQQgBC6zCBTT6Si2yq2tqDetJ+8ZF4vAnr4WgFXQDBZ37nUuHe9wDbm1l1AMwiLoEILsYMdkCsIgGW0haHQRgqALyrx6DVlFUe0WgVeMQFkEQVQXV2r1jlkEiplUDHq90RLU3ijVkIlHXAf9KEunHMA+Cz4Fd3OZWMe4SiVY81NeQhCMC6tWb6Ck1zuE4EDOFAJJIRXGkdEWAjtGMFjXEDP3ZBGop1yj3kBRMzoOD4VtEQQgWOaYSvd9qqIxp3tos6xjJu2+T/X3CiKFVH1WbhPc5jLxjhiBqiGYT3daLrlkNLBeQ9t2DKX1M9pvWQSSQiqMIyMtBHU7RuDvGup3+mitzmDuHFOpcAaLAdhD4r2EIG4Hi4OyCNw/H9VGwtmpc6OkY9bFLQQ0s7OCsFz8Kq/n0rEO37uKGbhZBNlkFLrRCOTzVELQ/v9yLh1DNExiEQhjyUgLgXINeVkE8YBcQ2rGgVfWULotfdTr8FDYweIghECvIxn1cg11tmw2LQJ3IVDZWUG0mShqZgzFrbJ6PhNDoWq0zJZYc2k4p1AzDYIoKvP6fxkKkZVCKo3nhPFjtIXAcg15BUODihF4zStWpGJhEDUtgu5CYFkuQaSP+lgE6sB3uob8LAJ1nyBiBCWfGMqsVV3sjBOo1tRzmc69qkH2+Wr/BUv9v2wfoAOYcQKxCIRxZKSFoN5gTCUiCHnkbQcVI7CFwMM1RGRWt6on0i2XGcBOkkGnj/rUETj3B5hBWLcaAiDYYHFR8668nnOpLl4r6ohFQq5tO1SabiGAOEHeR9QP5JIiBMJYMvJC4HW4AsGlj6oKYC/XEGA+bW9aLpdhuobKuuFpEUy3xQiqtTpKeh2zHsFilZ0VRPpoUat7CoFqI+GME6wVNCxk4q6upCmrxXYhQIvAXQjM6uIgJuIJwjDZ00IQCZvjFQNzDXlYBIBZmatcLvZTpEuRFuAUgiCazjU8K4tjkRDSsbAtWEoQvGIEQQeL21tQK+ZcXEOrLn2GFE2LIAAhKPsLgV5vuFZBC8JeZrSFgP2FADAPu8BcQx4xAsAMxG71aBGEQ4RomPo+nMaoN6DXG551FoApWMo1pITLO2soOIugpHkX3Hm5huZcAsWA0yLov2toq1JDKhZ2rbNQcwmWt8Q9JIwXoy0EDfe++U5i4dBQXEOzDotgu1JDIhryvT4RCffdNeQ1lMbJTDpq79OuKvawCOIRa4B9AI3nClXvYLEqcFtzzC6+kq9i0SriaidQi8ClqlhxcFoVlUnmkDBejLwQdLcIwn2vI+jVNbTlcLl022c8Gu67a8gWAh+LYCbVrNrdKHu3lwCaA+wDyRrSDddiMrXuXDpmTySr1urYKOk4kHUXgkzMnEkQhEXgJwQHlEUgAWNhzBh5IXBL43MSjwRgEdhC4HfARlHUDOhGA6tFrWWKlhuJaKjvbaj9ZhEo9mcTuJI3n7TVLF63lExFKqAOpCWfrCEA2JdN2FW7qp/Pgemk67WhECETiwSWPur1b24uHUMsHMIlsQiEMWOkhYDh7XdXBBMjsKZ++cQIptMqI0fHSl7Dvin3p1dFIhrue/poLxbBYjaB1aIGo97A8nYVsUjIHsHoRioAi0AzzPGYXhYBAByeTuLilnnAqifugx6uIcB0DwXhGnLrPKoIhQiHZ5I4v152/b4g7FVGWgiAHoQgHILe5wNWFX51swgAYLNcw0pBw74uFkEyCNeQ3l0I9ucSVvtkHcvbVRzIJXznJqdi/bcIvFpQOzk8YwpBo8G2D94rRgCYAeOgXENuQ3sUx+bTeGmt1Pd1BWGY7H0hCNI15GMRKAvg0lYF66XuQpCIhvofLO7BNWS3T85XsbxVsbtoepGKhVu6qvaDYtV9KI2TQzNJ6EYDayUNlyzLQPnk3XAbw9kP/GIEAHBsLo1X1stSSyCMFSMvBN1cLkG6hvyCxYdnzEPq6QtbYAYWPAKbikQ03Pf0UXU/v/RRZx990yLwPlzNe0VQ7vM+i3bnUe99HrLiAUubFZxbLeFALuFr6QThGtINs2W3rxDMp1Cp1bFSkAE1wvgw8kKwP+v/pB1I+mgPweL92QSiYcKT5zcBoKtFEI/03zVU7sUiyDUtlyv5aleLIB0Po9znpnOqLsHPIjgya46GfGW9hBdXinj1vozvPU3XUH/3ud2lMBAAjs+b+zq7Uuzr2oIwTPoys5iIbgfwcQBhAJ9k5o+0fZ+s798BoAzgXzLzk73c288iYGZEIyEUSwaMegMMgBlgmG2k0faarZ9R18Hjexc3TdfEdqWGfLXmfg82s5r+/sU1AGbK40trpdb7O37mwkYZZ64UcHo5v7M9trzf+r3HXlgBADx1fgtX8prHPcxF/uArzwMAzlwu4GvPXfa8/2MvrKJSq+OhZy7ZP9uyv7b7d36Gra/BjMfOrAIA/uaJJZxdKbbcQ92/1jBF8rf++hn7/+8nHjtr/39uv/9Dz1wCAHz0a2dc11Wv0fb/zP3fiPl3+eGqebj/ySNn8MLlQuvnY/2Msm6ev5THW1497/lvUxD2ErRbXycRhQH8AMBtAJYAPA7gvcz8vOOaOwB8CKYQ3ALg48x8S7d7xw+c4KP/88c9DkhBMCECCGY9AtmvzTedr9uvg/M1taEQnJIAACAASURBVLbrNttbOH+meY8QEf63n30Nbr/pwOD/soLQBSJ6gplP7uRn+mERvBnAWWY+Z23icwDuBPC845o7AfwnNlXn20Q0TUQHmHm5281/7Z8dd/1lhvXL/PFHXwQAfPi265rXWVkx7T+jG2ZfHpU003pwNA+EP/pvp6EbDfzRu29yOUSah8zvP/Sc7Z756L+4ueXa9vv/2y89h/WSjj//pTda3+88mMy9m1lLiVjY9XBr/v2B//DNc3jszCo+/f4fRcoqslLXOu//+19+Hs9c2AIA/PU9t9oumpZ7Wl//8VdfwKMvrODrH36r9x6t90u62TYiFPLeI4Hw+VMX8CePnMF/+42fwMFcsvWgdtz/D778HD5/agkA8OUP/gRO7M94HuL3f/OH+OOvnsHpP7jdN5YAmIVnqVgEYY8utorHXljB+z/9OB78wI/jjUdnfK8VhHGiH0JwCMAFx+slmE/93a45BMBXCA7kEvid22/wXfzCRhnffXkDv/GOE77XPXNhCz//Z9/CH737dfjFW476XvuDKwX892cv4323XuN73f5sAnd/6rv45Vuvwf/0xsO+1z59fgv/5bvn8c7X+T9Fnl7O484//RbufdtxfPhnrve99jvnNvDYmVW89cSCZ6tuAPjdd96AX3jg23jbdQu45fic7z1vPjKNR19YwbG5NCJh7xDSi1cK+PlP/CPeedMiPvoLr/e9pwq6H55J+QZif+FHj+Lzp5bwmgNZ3HQo65vm6uw35CcEL62V8K7/++/xxqMz+M+/+mbfe9pjKrtkqgnCuNGPYLHbb1a786aXa8wLie4holNEdKpRKXRdvNf00YefXUaDgf/87Ve6XqvVvOcVO3nriXl8/cNvw++96zVdr01EzV5D3VxxX3zqIvR6A//PN891vbZSqyMRDfmKAADccnwOX/nQT+DP3/fGrvtUGUjdMoe+8OQSKrU6HnzqYtc0TpWOmu7y5P6ma2bwtd96K/6yy4EN9D6c5q++ex5lvY5/OLuG71/M+17bba6EIIwr/RCCJQBHHK8PA7h0FdcAAJj5AWY+ycwn9y10D8bFIqGeeg09e3EbgGlBdDtgNaOBuE8mjoKI8Op9Gc820E6SsTAabM5D9uN7S1v2Hi53GZRe0et26+hu3HQo19O1vbaifvKVTfvrp89v+V5b1GpIRsO+Fobiuv1Tnl1HnfQ6nOY759Zx3X4z0+cfzq75XrtdMUWlW1sTQRg3+iEEjwM4QUTXElEMwF0AHmq75iEAv0ImtwLY7iU+0Au9po9e2DTbAhQ1w+6544Vm1HuyCHaCul+3WoKX1ko4Mmvm1L94xT9FsVKr+6aOXg3pePe5xcyM08sF3PG6RQDAuVX/ffoNpblaehlOU28wzlwp4K0nFnAwl8CZy/4WwXalhnQsjGgPgiUI48Su/8UzswHggwAeAXAawOeZ+TkiupeI7rUuexjAOQBnAfwHAB/Y7bqKXgrK6g3G8lYVJ6zc9G7FQJrRm2toJyirwa/xXK3ewEpBw49ZfvylTf/mZhW97tsP6WpQwuLXb2izXENRM3DymllMxSM416Xlgt9Qmqull1bUL62VUK01cMOBLK5bnMILl/1djd2qigVhXOnLYxozPwzzsHe+d7/jawbw6/1Yq51YJIR6g1FvsGdWyHpJg9FgvO5QDi+uFLtbBLWGbzHZ1dDLlLIr+SqYgdcdnsb/+8RSd9dQrXfXUK/0Mrf4woZpXR2ZTeH4QtrOv/eiW+fRq6GX4TQvWBbADYtTeHFlCt86u4ZaveH5xO/XeVQQxpk9bwPbA+x93EObJfOwuG5xCoA5BtEPzaj79hm6GtSTu18HUtV188hMEvOZOFa6CEFZN/ruGuplSpmyVA7PJHF8IYOXVv0tgqLmPYvgaunFInjF6hJ6fCGN6/ZNoVZnnN/w7hyar9S6DkIShHFk7wtBDwPs1WAWFTRc6yoEAbiGIt0H2Nvtl6eT2J+N92ARNLrm0O+UXoLFS1a85dBMEgenE7hS0NBoeAfBgxCCXobTLG1WMJeOIRWL2C0sLvq428Q1JEwqe14I1IGt1f182qYQHJxOIh0L9xAsHo5raHmr2X55MZuwB7R4UdX7Hyy2LQKfYPHSZgXZRATZRBT7s2aba7+B7kG4hnoZTrO0WcYhqzmg+lPNPHBjq6KLEAgTyZ4Xgl5cQ8oimE3FsDAVx1rR+9ACgskasl1DPhbBakFDIhpCNhHFvmyia1C7XDN8O49eDXYdgY9FsFbUsN/qaqp6QV3xsV6CyBoCuncgvbhVsbua7p+KIxIi25pxQywCYVKZCCFQQ9unUzHMZ+JYLfg/aWu1RgAxgu6uoc1yDbPWPOH9UwlslHTUfDKiKrrZiqKf9BIs3ijpmLGmnKnusCs+n2kQWUOA/3AaZsbFzYrdLjwSDmExl/B0DWlGHdVaQ4RAmEj2vhCEzQPGL4V0o6xjKh5BLBJqGTrvRZCuIb86gq2yjmlLCGbSUes9771WAggWxyMhhAi+U8o2y7o9oU1ZBmoucjtGvYFKrY5MvP8HrJ9FsFHSoRkNHHTMPT40nfRMybVbUIsQCBPI3heCHiyCrXIN09bBmk12H2gSpGtI84kRbJZ1WwCUIKi2B+0ws5U+2l8hIKKuc4s3yzXMWhbBgjWHwSueUbLukw7AIsgmoyho7kKpMsOcbcwPTSftgHw7eXsWgfc8Z0EYV/a8EMR7jBHMWL/g2UTUtzcOMwdaUOaXPrpZrtn7dM5EdkOvN9Bg9NTeYqf4zS1mZmw6Ps9oOITZdMwzJbdkTycbbIxgxbJQ9jkGGy3mEriSr7pmOIlFIEwye14IerEIipqBrFWAlE1EUNQMz3THWt0cRNJLr6Gd0FuMQHcIQcx+zw2V3tlvi0Dd02tucUEzYDTYtggAYDYds+Mw7SghGHSwWGWGLTj6Fh3IJWA0GGulTtESIRAmmbERAs0nRlDSmtk12WQUzEDR44m3l3nFV0MiorKG3PdZbzC2KzXbElCFTV6uIRVr6HeMAICva8gZeFfMpmJ2ZlY7xUAtgijylZprE0FloSxMOS0CM17g5sZSsRgRAmES2ftC0ENBWVlvpi+qilTlE26nOa+4vx9NJBxCJESeFoF5oDUP2KZF4L5Pe15xABZBOu7tGlL7mU03D8zZdMzTclFCoD73fjKViMBosKu4rhY0pGLhFktEzWt2ixOov9eMVBYLE8ieF4JeYgRl3WERdOla2cvg+qvFnEngvk91kCqXSyoWRiwc6uoaCsIiSPZgEcw4LIKZtI9FYH3OmUCEwLvf0EpBw76p1nbWi5YQuFsEOkLU/PchCJPEnheCXmIEJc1pEZi/6J4WgfXE3u86AsDMHPJKH233URMRcqkotj0sAmVZBGIR+ASLN1yEYDYdxWa55hp3KQToGvIbTrNaqLa4hQDThRULhzwsArOquNuQH0EYR8ZHCDxiBPVGa5plNunfrCwo1xBgWgRebajVfpwulJlU1NMiKAcYLE76BIvVfmYcweKZVAz1Brt+psoimAqojgBwtwhWC1qHEIRChP25OC5vd9YSODO2BGHS2PtC0CVGoJ7A07E2i8AjhTRw15BH+qgdVHUIQS4ZRb7iLljq7xVE+mg6FvG0XDbLOsIhsp/GgaY7a8NFtIp21lAwlcWAu6ivFrSWjCHFgax7LYFZzCduIWEy2ftC0MU1VLYOolRcxQi6BItrwWQNAaZryCtGoJ5qpxw+6mwiaruM2gkyRpCKhz2bzm2UzMwm50xhWwhc0jKLmtHzmMqd4tWKulqrI181sC+b6PiZxVzCtavrZkksAmFyGR8h8HAN2ZWtlkWgnrhLHsFQ2yIIIkYQCXtmDbm5hnJJ7+I39cTe78E0AJCKRqAZDdRdfP7OYjJFUwg691qoGoEEigHvYLFbDYHiQC6B5e1qR8rpdqWGnFgEwoSy94XAetL0GmCvnmyVLz1mpXF6PfEGnzXkkT5qCUHacbBnk1FPyyVIi0C5cdwCxmYbjFYhsFNdXTKHipqBqQACxYC3ReBWQ6BYzCWgG42OtFxnMZ8gTBq7EgIimiWivyWiF60/Z1yuOUJEjxHRaSJ6joh+czdrutzfd4B92e51E7GvN1souB/I6qCODdg1VKyaw1uc4zbNXjruVdCVALOGkj6tqDfLut0hVeEbI6jWArMIvIbT2BaBixAcsIrKlh0BY82oo6zXpYZAmFh2e9rdB+BRZj4B4FHrdTsGgN9m5tcAuBXArxPRjbtct4VYxFsI1MhF54GZjkfsIGY7yiJIBGARxH2CxYVqraPoKpuIgLmZgumkpBmIhCgQwVJWiZsQbJRqdmM8RSoWRjwS8rQIgkgdBbyH06g5Du11BECzqMxZS6CqiqfFIhAmlN2eIncC+Iz19WcAvLv9AmZeZuYnra8LAE4DOLTLdVuIRULQPSaUlbXWGAFgCoFXnrxqMZEIKkbgYYkUqp0Hphqk7uYeKuv97zyq8JpSxszYcnGhEBFm0zHXKWVuf69+4tZvaLWggQgt/ZAUbtXFdkqsCIEwoez2tNvPzMuAeeAD2Od3MREdA/AGAN/Z5bot+LmGlEXgPDTTPnnyynUTRIwgGQuh6rHPomZ0WASquMwtc6ii1wMJFAOOucVt8Qy3hnOKmZR747kgg8WA5T5zcQ3NpWOumUpzGXNSmdM1tFmS9hLCZNP1N5SIvg5g0eVbv7eThYgoA+BvAPwrZs77XHcPgHsA4OjRoz3d2881VHbpfmk2VevSdG7gWUO1DtdE1qfmoaT3f0ylQqXatlsEbg3nFLPpmGcdQVDBYsA9xdYsJutMHQWAcIiwP5tosQhUYz/JGhImla6/ocz8017fI6IrRHSAmZeJ6ACAFY/rojBF4LPM/GCX9R4A8AAAnDx50r1XdBuma8g/fbTFIoiHcWnLq3VDsJXF1VodzNyShw+YT86HZ1Mt76kqaDfXUEWvBxIoBrznFtuzn9OdB+ZsOtYxD5iZzRhBgBbBTDqKl9ZKLe+tuLSXcLKYS7TECNaK3ummgjAJ7Pa0ewjA3dbXdwP4UvsFZJ54fwHgNDN/dJfrueKfNWQgHKKWgz0Vi9guo3Y0o45YJNRxUPeDRDSEBpszD9opaEZLtS7QdA25VReX9XpL3KOfeAWLt+wOne4WQXuMoFozaxGCGFPpXLe9fuFKvorFbO9CsFLQECLTbSQIk8huheAjAG4johcB3Ga9BhEdJKKHrWveAuCXAfwUET1t/XfHLtdtIRYJ+dQRmEFV58Gejkc8YwRarWHPDug3flPKzKyh1gPTDha7uIbKujEAi6BVgJoWgXuMoFA1UHNYZmqMZJAWgWqBrVJsjXoDqwUNiy5VxYoD2daistWChtl0vCV1VxAmiV39hjLzOoB3uLx/CcAd1tf/ACDQ3zDfGIFudDw5+3XX1Ix636eTKeKOKWXOdse1egPVWqMju0blybsFi8t6HQengxICd4tAZde4xggyzaIy1dqh2XAuSCGIo95g5K0Yy3pJR4Ph2l5CsZhLoFKrI18xkEtFXVtWC8IksecriwHTn+8VIyjrdTv4qUjFzX77boVa1VojkNRRoDmlrH2AvVt7CcDMk88m3KuLywFmDSWiIYSoM1i8UepsOKeYs6wEp3soyKE0Xusql89+P4tAFZXlzcyhlUK1ZbaxIEwaYyEE3SqL3SwCoDM9ErAsggBSR4Gma6h9XfvJ2WUoSjbZWTAFtA7b6TdEhHS8Mz9ftWp2i5+4tZlQsY0g6whUuwu17hWroZyfa2ixrZZgJS8WgTDZjIcQ+FUWa50HZso6mNz6DQVqEXgMsFcxALcDM5d070AaZEEZYKZltscmzIZz7oHfuUynRbBV6Zxd0G/aLQIlBPt9nvCd1cX1BmO9pPtmGQnCuDM+QuDjGkq3+95VnrxLlW+QFkHSFgJ315Cby8XNNVRvMDSjEViwGHCv2N1waTinaM5YdgiBat0Q4EB4FbheLyoh0BAOkW8G0MJUHCEyLYKLmxXUG4wjMynP6wVh3BkPIQiHOvzuCrfCK+VbH7xFELLWaHMNaT6uIZeCKRXoDip9VK3bXrG75dJwTqEsBXUgA80gdzZAIVCHupoxcDlfxULGPwMoGg5hYcqcVPbyulmDcGw+HdgeBWHUGQ8h8LMINLcYgXdTtUHECNqFQB24bmmWbjMJ7BbUg7YISjVPiyASDmG6bbTmVllHMhoOZIqaIhoOYd9UApe2zMDvlXwV+3Pe8QHFYs6cVKaE4FoRAmGCGR8h8Ok11Jk1FLa/1442CIvA6C1rCLCCxW0FZUHOK1a0CwEzWz37vZ/uZ1OtRWVb5ZpdFBckB6ebQvDKehlHZpJdf+bITBLnVkt4aa2EVCwswWJhohlrIWBmj6whH9dQgBaBuq+XReAqBIkoKrV6y9+v2UgvyK6erZZIvmqg7tFwTjGbbm08t12pDWQO8MFp8+leNxpY2izjeA9P92+6ZgYXtyr4xplVHF9IB1JJLgh7hbEQgnjYdA21jx9U4xbbLQJ7ApdLdbFWawTSZwjwcQ1pBmLhkKsAqUZozkO5MkCLQH2mWz20ajbbPTizhgZlESRxcauC8xtlNLg3f/8t184BAF5aK+HkNbNBb1EQRpqxEAKvucVlvXMWgfO1m2uoWqsH5tP2ChYXqp0tqBV2B1JHwHgQrqFsMop6g+2aB3XAtw+lcTKbjmHNGSwuD8YiODSdhG408N2XNgD0JgQ3LE7hyKzpQvofbz4Y6P4EYdQJzrcwQGwhMBotT9Xt84oVKXsmr1uweBAWQWeMwFMIrA6k2y1C0Dl1rd845wGnYhHf9hKKfVNxbJQ0GPUGIuEQNss6Xp+cDmyPihP7MwCALz9zCQDwqvlM158JhQh/9b/cirWijtcfCX6PgjDKjIdFEG4KgZPmk3Ok43q3AfbMZn5+UL2GouEQwiHqTB/1metrdyB1BG69LJ1+MtVmiagnfb9WzfuyCTTYLO5ShVrzU8FP/brxQBZEwD+dW8e18+me5wocnkmJCAgCxkUILCug3TVkB1XbYgRqgH27EKgOpkFZBIBZVOZqEXi0alauIadF4DZjod+o7KDNshICs2e/qiB2Q7V1uLxdxYYlBvs8BsT0k+lUDDcfNg/0n7zed0ieIAgujIkQeFgELvOKFZl4pKOyWBWlBZn3noiGOtpQ+7mGci5ziysDcA2p7KCNkikA60UdqVjYN1NJNXq7kq9ipdC91UM/+cM7b8K/OHkYH/jJVw1kPUEYJ8YuRuDEbV6xIuUywN4eUxmgRZCIhu2sH0XBxzXkNpOgWDVAFKxraC5tHuBq6MtaUcN8l8Et6tC/UtAQtT5Dr5GR/eZ1h3P44/fcPJC1BGHcGAshUAd3+3AauxWDSzM3twH21QFYBJl4pMMlla8anmmW8UgIsXCoxTWUrxrIxCIIBThIRWUHOS0CP7cQYE74CocIV7ariIXNvUmhliCMPmMhBF7poyXbNeRiEbgMsB+ERZCKhVvSVusNc65v1qXPEGDGM7LJaEt1sZ8rqV/EI2Fk4hG7UnitqOHIrH9jtnCIsJhNYGmzbH+G0tVTEEafsYgRxD2zhlSw2MUiiIdRbLMIlEURpEXQPiZTzSLwa8xmziRwuIa0zrGWQeAsEDNdQ90zgF61L4Ozq0Us56uYTkUD/SwFQegPuxICIpolor8lohetP2d8rg0T0VNE9JXdrOmGZ4zAOnCTLodR2iVGoNI6g7QI2l1D6oB3a0GtaG9FXagagc4BVighqNUb2CjpPfn7X72QwQ9XSji3WsSxOWnkJgh7gd2eePcBeJSZTwB41HrtxW8COL3L9VzxzBrSDSSjYdeWxGkXX/0g0kfb1+2lVXMu2SkEQbuGAHPoy3pRx6WtChoMHO6hmdur92VQqdXx7XMbOL4gQiAIe4Hdnnh3AviM9fVnALzb7SIiOgzgXQA+ucv1XPGMEeh1u69QO+7BYvN1oK6hWLglbbVpEfi5hlpnEhSqg3ENLeYSWN6uYGnT7OzZixDccGDK/voNUqwlCHuC3QrBfmZeBgDrT69qno8B+NcA3HtF7xLPymLN8Mx7T8cjqNTqqDsG2NsWQUBtqNW6Ja3ZzE0FgVUrCTeyida5xUVtMBbBsbk0Nss1PHdpGwB6muJ18+Fpu6fS26W4SxD2BF1PEyL6OoBFl2/9Xi8LENHPAlhh5ieI6O09XH8PgHsA4OjRo70s4VNH4D3XtzmcxrCfrlV+fyKgNtSAKQSGNWoyEQ33ZBEo1xAzg4iQH5Br6Jo58+D/h7PrCIfInvXrRzhEePS33461QvcsI0EQRoOupwkz/7TX94joChEdYOZlIjoAYMXlsrcA+DkiugNAAkCWiP6Smd/nsd4DAB4AgJMnT7LbNe3Yff7bKnbLuuFaQwA0awtKWr0pBLXgWzeoVNaybnY5zfcQI8gmozAa5myFSJigGw1f4egXqovnN3+wisMzSUTCvVlKh6aTODTd3Y0kCMJosFsfyEMA7ra+vhvAl9ovYOZ/w8yHmfkYgLsA/J2XCFwtqtVCezO3kuZjEbhMKRvECMimAJnrKpdPxkOwAGfjuZo9Nczv+n5xdDYFNa/lpoO5wNcTBGE47FYIPgLgNiJ6EcBt1msQ0UEieni3m+uVRET1+W91DVVcppMp3KaUqa6ebumm/UId4GpgfaFaw1Q84jtsvTmTwHBYEMELQSIaxqsXzJbOP/7qucDXEwRhOOzqNGHmdQDvcHn/EoA7XN7/BoBv7GZNNyJWW+kOi8BlXrHCnlvsyByq1Opmi+oeXSBXgypuUzUM+Yrh6xYCWmcSqCd01QsoaD521+vxX5+6iPe86fBA1hMEYfCMRYsJwHyKr9TaYwTeFkEm3mkRVHQjULeQua55f+XiyVdrXQO/zg6kRsO0err1/ekXrz2Yw2vFLSQIY83YCEHcpc9/SfOxCFzGVZZ9soz6RfugmXyl1tUimE5aLaHLup0Z1a0TqCAIQq+MjRAkoiFoDovAqDegGQ2kot0sglbXUJDxAaA56nHbGv2YrxpdM2z2We2dV/JVqJo5vyHygiAIO2Esms4B1uQvR/po2RIFr8ri5tzi1qyhoF1DyiJQk7+2y3rXwG8iGsZsOobl7SrWSxpyyahdOyEIgrBbxsgiaB34oqaTeVYWx1qzd4DBuIai4RAy8Qi2ymaB2FpJ950DrFjMJnB5u4pELDyw+IAgCJPBGAlBqCVGULKH0rgf7OEQIREN2SmjgOkaGkTFbi4ZxVZFR1EzoBuNng72A7kElreryCYjmB9QxpAgCJPB2PgXEu2uoS4WAWDGCYpaq2soaIsAAKZTUWyXa1gvmnGCXlJB9+cSuJyvYnm7ioUBzQEWBGEyGC8hcLMIfA72dDyCslMIanVf4egXM6kYtio1rFtjIHuxCA7mEtgo6XhlvYxrpc+/IAh9ZMyEwGER+EwnU6RikZYpZeUBBIsBIJeKYrOsY82yCHpJBb1+MWt/rXoACYIg9IPxEYJIqEUI/OYVKzLxcFvWkBF4+igATCdN19Ba0bQIehGCm480i7p+9JjnIDhBEIQdMzbB4mTs6iyCLSufn5kt11DwQjCXiWOjrONK3hSC2XR319C+qQR+46deDaPBOCrtnQVB6CNjIwTtMYJiDxZBNhnFhY0yALNhXYP9g8v94kAuAWbgmQtbmE3Heq4J+PDPXB/wzgRBmETGyjVUqdXtyV+qh5DXPAJATf4yC7sK1p+DSB9VA16++9KGDG8RBGHojI8QWE/+atxkUTMQj4QQ9ekkqmYBM7Pd+2cQQqBcO5VaXdw8giAMnbERAhXkVXGCXub6ZhNR1OqMaq1hWwSDmPx1jSP986aDWZ8rBUEQgmdshKC9ZUSx6j2mUqF6/Dgnfw3CIgiHCD923Bz08tbrFgJfTxAEwY+xCRan27qJljSj6zjH5uQvpxAEbxEAwAO/8ia8sl7Gaw6IRSAIwnAZH4vA6ilkj4DsRQhaZgEPLlhsrhPFTYdk4IsgCMNnjISgdeJYbxaB5RqqGAN1DQmCIIwSuxICIpolor8lohetP11LXolomoi+QEQvENFpIvqx3azrhooRqEKyomYg0+MIyO2KaREQwXO0pSAIwriyW4vgPgCPMvMJAI9ar934OICvMvMNAG4GcHqX63agnv6LjhhB92Bx0zWUr5oWRChE/d6aIAjCSLNbIbgTwGesrz8D4N3tFxBRFsBbAfwFADCzzsxbu1y3AxUjUK6hQtXAVI/B4s2SGSweROqoIAjCqLFbIdjPzMsAYP25z+Wa4wBWAfxHInqKiD5JRH1vn5mON9NHa9a84m4WQSwSwkwqipVCFYVqTeIDgiBMJF2FgIi+TkTfd/nvzh7XiAB4I4A/Z+Y3ACjB24UEIrqHiE4R0anV1dUelwDikRDCIUJJM3pqL6HYn01gpaBhvaT31PxNEARh3Oh6UjLzT3t9j4iuENEBZl4mogMAVlwuWwKwxMzfsV5/AT5CwMwPAHgAAE6ePMnd9ufYC9KxMEqage2KmQqqgsF+LEzFsZKvYqtSw48cnu51OUEQhLFht66hhwDcbX19N4AvtV/AzJcBXCAi1TrzHQCe3+W6rpijJ+vYLJtCMJPqLgS2RVDUMS9D4QVBmEB2KwQfAXAbEb0I4DbrNYjoIBE97LjuQwA+S0TfA/B6AP9ul+u6kopHUNYNbFozBmZ66vMfx/J2FUXN6GlAjCAIwrixq+goM6/DfMJvf/8SgDscr58GcHI3a/WCGkavhs3MpLoLwf5swv76mjnpBCoIwuQxNpXFgDnpa6OkY6PUu2vIefhfMyuzgAVBmDzGSgjm0jGsFTVslXWEqLeW0jc7AsQn9meC3J4gCMJIMlaJ8/NTcawXdWyUdOSS0Z6qhGfSMXz8rtejwYzEAAbXC4IgjBpjJQRz6RiMBuP8Rrmn+IDiztcfCnBXgiAIo81YuYYWpsyssKoIiAAABGBJREFUnx9cKWC6h/iAIAiCMGZCMJc2heBKXsNiLtHlakEQBAEYMyGYn2q6gw7PSCqoIAhCL4yVEBzIJu2vj8wkfa4UBEEQFGMlBDlHXOANR11n5AiCIAhtjFXWEADc/7434anzm7hRhsILgiD0xNgJwe03LeL2mxaHvQ1BEIQ9w1i5hgRBEISdI0IgCIIw4YgQCIIgTDgiBIIgCBOOCIEgCMKEI0IgCIIw4YgQCIIgTDgiBIIgCBMOMfOw9+AJERUAnBn2PkaEeQBrw97ECCCfQxP5LJrIZ9Hkemae2skPjHpl8RlmDnzo/V6AiE7JZyGfgxP5LJrIZ9GEiE7t9GfENSQIgjDhiBAIgiBMOKMuBA8MewMjhHwWJvI5NJHPool8Fk12/FmMdLBYEARBCJ5RtwgEQRCEgBlJISCi24noDBGdJaL7hr2fYUFER4joMSI6TUTPEdFvDntPw4aIwkT0FBF9Zdh7GSZENE1EXyCiF6x/Hz827D0NCyL6Lev34/tE9FdElBj2ngYFEX2KiFaI6PuO92aJ6G+J6EXrz67jGkdOCIgoDOATAN4J4EYA7yWiG4e7q6FhAPhtZn4NgFsB/PoEfxaK3wRwetibGAE+DuCrzHwDgJsxoZ8JER0C8BsATjLzTQDCAO4a7q4GyqcB3N723n0AHmXmEwAetV77MnJCAODNAM4y8zlm1gF8DsCdQ97TUGDmZWZ+0vq6APOX/dBwdzU8iOgwgHcB+OSw9zJMiCgL4K0A/gIAmFln5q3h7mqoRAAkiSgCIAXg0pD3MzCY+ZsANtrevhPAZ6yvPwPg3d3uM4pCcAjABcfrJUzw4acgomMA3gDgO8PdyVD5GIB/DaAx7I0MmeMAVgH8R8tN9kkiSg97U8OAmS8C+D8BnAewDGCbmb823F0Nnf3MvAyYD5MA9nX7gVEUAnJ5b6JTm4goA+BvAPwrZs4Pez/DgIh+FsAKMz8x7L2MABEAbwTw58z8BgAl9GD+jyOW//tOANcCOAggTUTvG+6u9h6jKARLAI44Xh/GBJl67RBRFKYIfJaZHxz2fobIWwD8HBG9DNNd+FNE9JfD3dLQWAKwxMzKOvwCTGGYRH4awEvMvMrMNQAPAvjxIe9p2FwhogMAYP250u0HRlEIHgdwgoiuJaIYzMDPQ0Pe01AgIoLpBz7NzB8d9n6GCTP/G2Y+zMzHYP6b+DtmnsgnP2a+DOACEV1vvfUOAM8PcUvD5DyAW4koZf2+vAMTGjh38BCAu62v7wbwpW4/MHJN55jZIKIPAngEZgbAp5j5uSFva1i8BcAvA3iWiJ623vtdZn54iHsSRoMPAfis9bB0DsD7h7yfocDM3yGiLwB4EmaW3VOYoCpjIvorAG8HME9ESwD+LYCPAPg8Ef0qTKH8513vI5XFgiAIk80ouoYEQRCEASJCIAiCMOGIEAiCIEw4IgSCIAgTjgiBIAjChCNCIAiCMOGIEAiCIEw4IgSCIAgTzv8PP9Mhm3SpduwAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.plot(ifftt, ifft.real)\n",
"plt.xlim(0, 10)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.3"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment