Skip to content

Instantly share code, notes, and snippets.

@wakusei-meron-
Created April 14, 2019 13:39
Show Gist options
  • Save wakusei-meron-/d09888e4d7a78632e49cfdf3f1c93329 to your computer and use it in GitHub Desktop.
Save wakusei-meron-/d09888e4d7a78632e49cfdf3f1c93329 to your computer and use it in GitHub Desktop.
Desktop/PPP/python/audio/Untitled2.ipynb
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "# parameters\ndt = 0.01 # サンプリング間隔\nsr = 1. / dt # サンプリング周波数\nf1 = 4 # 周波数\na1 = 2 # 振幅\n\nt = np.arange(0, 2 * np.pi, dt) # フレームごとの時間\n\ny = a1 * np.cos(2 * np.pi * f1 * t)\nplt.plot(t, y)",
"execution_count": 121,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 121,
"data": {
"text/plain": "[<matplotlib.lines.Line2D at 0x13b51afd0>]"
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/plain": "<matplotlib.figure.Figure at 0x13b03e0f0>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAD8CAYAAAB+UHOxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzsvXmwbddZJ/ZbZ7r3DZIsWUISsoUYTBk3uBle2UB3OnRhKGFoGwwkuEiDu0OUonCgQ1JgoEIn/JF0F9VNisaBiKEZimAcwI06qNuYZjCmMfaT5wHbsjxIsmw9S5b0pnvPtPLHPmuftdf+xjO9d+/bX5VK7977nb3XWetb65t+37dCjBEdddRRRx1de9S70gPoqKOOOuroylCnADrqqKOOrlHqFEBHHXXU0TVKnQLoqKOOOrpGqVMAHXXUUUfXKHUKoKOOOuroGqVOAXTUUUcdXaPUKYCOOuqoo2uUOgXQUUcddXSN0uBKD0Cim2++Od51111XehgdddRRR0eGHnjggc/GGG+x8F7VCuCuu+7C2bNnr/QwOuqoo46ODIUQPmHl7UJAHXXUUUfXKHUKoKOOOuroGqVOAXTUUUcdXaPUKYCOOuqoo2uUOgXQUUcddXSN0toKIITw3BDCn4UQPhBCeH8I4UcInhBC+PkQwoMhhPeEEL563fd21FFHHXW0Hm0CBjoF8D/FGN8RQrgOwAMhhDfFGD+Q8XwLgOct/nsxgF9c/L+jjjrqqKMrRGt7ADHGx2KM71j8+zyADwK4o2B7OYDfjBW9FcCzQgi3r/tuiqazOV77Zw/izR8+Z/7MX3/0Cfztp58x87/nkafwwCc+Z+b/6LkLrvF85pkD/If3PmbmP38wwe898Ais13tOZnP8zts+ielsbuKPMeL1Zx/GxcOpeUx/9J7H8Pj5AzP/n3/ocXzssxfN/Gc//iTe9+jTZv4PfOoZ/M1DT5j5P/nEJfzp337GzP/EhUPc9+5Pmfkvj2f43bd/0rxms3nE6972SRxOZ+Z3vOGdj+DpSxMz/x+//9P41FOXzfx/9eBn8ZHPnDfzv+vhp/DOT9r3zYOPn8dbPvJZM/+nnrqMN77/02b+py9P8AfveMTMP57O8bq3fRKzuW3NYox4/dsfxuWxfc3+/bs/hc9eODTzr0sbzQGEEO4C8FUA/qb40x0AHs5+fgRtJZGecU8I4WwI4ey5c/ZDM1G/F/DLf/mQSxBe+ctvxd3/51+a+V/2C3+F7/zF/2zm/8Z/9Rf4vl97m3mzv+rfvh0/+NvvwDMHts37v973AfzP/++78a6HnzLx/9ZffwI/8Qfvxeve/rDODOBtH3sSP/Z778H/fv8HTfyfuzjGD/0/78B/95sPmPjn84hX/du345t/7i9M/ADwXb/01/i2f/MWM/9Lf/4v8V/f+1Yz/z/6hbfgn/76WbOS/OHXvRM//DvvxGNP2w7Qf/XHH8KP//578ecfssn4fe9+FK/5g/fi3r94yMT/kc+cx//4u+/GT7zhPSb+y+MZ7vmtB/DKX7bP0ff+yt/gm37uzWb+b3/tX+E7/i/7vnnJv34z/ptfLY8SeTz//W89gEtjm6Hyk294L3709e82G3+/+paP4TV/8F78u3c+auL/iw+fw4/9/nvws2/8kIn/M88c4H/4nXfih377HSb+TdDGFEAI4TSA3wfwz2KMdnO6oBjjvTHGMzHGM7fcYqpmLseBL7r5FD567oKJ/4LDql2Mr/631RJIdM6o2T/4WDV9D52zWcQf+2z1XT/1lM3i/swzFd/j523jeezpiv+TT14y8ae5/5BxYz22GM9kZvdgViWrNfb05Ur5PvI524H+oU9X3/ljxjV7+HPVXKa10OjxZ6q1eszI/4knqud/9HGrDF1sfE6jpy6NTXyJ8n0zd+6bJy/a3pW+g3XfJL6Hn7StcVLu1n2c1iyttUZp33zY4VWtSxtRACGEIarD/7djjH9AsDwK4LnZz89Z/G4r9MW3nMZHzUKwVBQHE/1weCITxkcNh0Mu7NbNuOS3KbH0ioeMSi99h3PGEM1Di401NR7QSZCNDk/je1q8pPyQshxEeejKE2YCYDYkgOjiH08rJfZx44GbQjPnD2wGS70G8K2ZlfL9ZVHIn84U16cNSix/plWul2Oz8SdZs/Ino8AaJvvYE9UcpbXWyHpmbZI2gQIKAH4VwAdjjP+aYbsPwPct0EBfC+DpGKM9yO2ku24+hXPnD00Heq79HzUs7MOZFWyxiPN43sMG/vwAtFrcT1w8dPGncVj5H3nSZ62mObUqgNxCslh7Ob/FesvX1fKdc6/QOkeHi01u5U9elUUmAODhhbFh56/4njLmAPI5tYS9Hsn4P/20Lhf5OlnmKJc1C3/ujVvnKB3o5n2wWAPvPrOGBdM+syqMTdAmUEB/D8A/BvDeEMK7Fr/7SQB3AkCM8ZcA3A/gpQAeBHAJwD/ZwHtZun6/+lrnD6bYH/ZF3vNZnP2CwbrKLbALh/rmeibjP28IN13OlJY1PJXG5OY3WpPpO1jGXz2/mpfxbI7JbI5hX7YzmnM6xbNP75n5zxvWoLHGhu/glYkYY+1leNfAO6fblgkAuHg4ww0n5TVryLVp36yzz3T+nMc+p759UK+BmX+1fXZxPMN8HtHrBdPn1qG1FUCM8S0AxJHGyqz9oXXfZaVTe9XXung4xS3XyYdJLjgWlEvOc+FQ9zAuOp/vHc8qh8/FsZN/wWdFAeXzcvFwimedHJmebx1Tc071NSjH43n+BUNC8fJkVofhLDJR8a02p1b+xHdpPMNsHtFXDpPyO99wcmjmv2iYowtO/vVkQuefz2M9DvuczszjyZ/r5QeAS5MZTu9tv1nzsawEThO3Dcsh57lgQOlccAryBafHcDid18lTq6BdcFqHF4rDROdfzovNOszn1DdHFi/sgtOaXGs8BpmIMdbjsFqH6Xt6rU/AeECv9Z2d+8wiE4erP9/Cf2kyq0OU2/LCLmQKwJLb8n6HTdA1rwC8loPXMvG6shfXsFbt1mSyfGzWqvc7r2Ptea1JrxfmXQPTGjs9ksPpvFakdutwYX2ObYdJPo9ez9OkVLdsoa/jObtlwhjm83ttFd88AgcTPa7vndNN0PFUAIscgNVySN6xx0IPwSdovWCMfS42n5X/gvP5k9kch9N5xX84NUHyzmdzZNuMs+WcOtfAYh021sxhTdrndLkGnsOqF4xe5EE2HofX1gtVYv2SAcqa+NO/Vf7GPrCE1VbbN15+65xeaMypJS+05LcYHQeTOebRuWa5XBs9Vc+cboKOpQKocwAma3KG267fX/zbwD+eYm/Qw40nRy5L6bbr942Hyazm91jbVv6LGT9QucKWz7jm6GDi48+eb7GgLx5OcfPpPfcBfatxDS5ka2BSSAfeNW7yaxb94XSG8WyezdHm5bqxxsbv7B3P6b0BrtsbuOXaP6d2w+y26/ddhtlt1+9jMotqRXbKzbnkeuyb003QsVQA1+0tUUAaXTiY4KbTIwx6wSz4p/cGOLXXd1not96wb7NkFoJ26w1GBXCwfL4l1ng+488/z9FsHnFpPKv5rRZ6/XzjZl/yG6y3wylO7w9w2nGY7A97uOHE0GhN+tbg/GFzDSzjSfzTeawhpBylw6NeA6McedbgvHMNLhxOccv1HpmY4PTeAKf3B65983lGw2adNfB4GNZ9k3JzVv7E49lnm6BjqQByFJBGFxeWyam9gdnSqA6fodlVHvYDbjo5clufpvGMl5aMJdaY86fxefitc+rjn+LzrttzhdVO7/kUwOm9IU5b13gRYrn1Or/16bVu858tz89/5mg6m+Ng4vMYml6eTa6v3x/g1KhvlolTe/1qnxnzQieGfVx/wrbP1lmD8XSuYu9zzzz/2TKe/H3amDoPYAN0ctRfHCa2ST81qg4Tm2U1w8nRAKf3+iZL6eLhtOLftx1WebjCNv5KED9vAXfVYo3p+YnfevjcahTklCy71WEdVofDAKdGNuuwmtPqMLHyp8PHKhPDfsCNp0buNbMdJsUaKN/hgpO/VmDONbjVEQK6uNg3njk9tTC0bB5DpTBO7/VNyKrGvjGOP/HnP/PjafLr+6y5L7U5mi5yc1b+TdGxVAAhBOwNeqprDQAH0xn2R33sDW38h9MZTgx72B/2bc+fzHBi2Mf+oI9DAxIgVS/fcGKIySyqsMvEn7D22juSh7Dkly2ZJX+FC9din2lOlvyOORr2TN0uDybzBX/fyJ8937gG+y5+3xy11kCZo9YaK/yHNb9vDU7t9THoBducTmc4Meq79sH+sI994748rPmtz5+j3ws4vTfA4XSmhkLba6aswdQ3p0v+tGaKTExLmbB3EF2HjqUCAID9Yd/UCuJwMsf+IB3QtsNkf9jH3qBvgnYdTObYH/awP+zVQqHxj/o9nBxVFczqgZspDAv/QcF/YDx8an7lO6cD8/TeACHoCia9wzenywPaugZ7C4VhXYN9F79vjtr8NoWxLf7pbI7pPFb7YLiCXFvWeLqcU5NMTBcyYR7PDPuDajzzqDcW9K5Buc/0NVhXJnbTDuLYKoC9gc16O5zOsDfsOTyAOfYGid+C+Z5hb5AE2crfw96gWhrtO6Qxmw/ogl/3AAoFo22UxZxUB3pPVTBpTL45XfAPbB5APqcmmZgs+S1eWILVptyTepiUa6BZ9NNSyW+WP63R3nAxR5Y5nSzkemA1tLI1MHkA84zfKBMLI6L6WbO4vXPqXbOS3yYTp/b66Bu9sE3QsVUALmtv4QGYLJncNfVYSgvBt7imyVoFoH6H1gG9YQ+gpWBUj6H6u9Xaq63PYVqD7XgAHov+MLNWq5/1OW3yXx0ewHX7RqU9WSptswcwnWNv4QFYDafl8+0eQCVDPg+g+tnmqV6/v5pF7zWcrB5ACpN1HsCadDV6ADFWDdJ0/ur5gM0DCGFZ/Gb2GE76BPn6WvBt1mey3qweyXY9gLnLAziY+L2wnN/sARjjybVF7+TfH9os7oN8DQwWd4wR4+l8Zx7AeDZXCxZX8QBGg55Zabs9gEm5xpoSzveBTa43QcdWAVisvdk8YjKLK3sA9vhzdpgYBKfyGOweQIrdWvkBvwdwctTHsB8Mzy88AON4UsjIk6TdM3oA6fBJMqEXXjU9ANMarOABWK3PwxU9AKvFnXsAlph7+n77Dg+gygH0VvIA8ney/JOm4WTxAJoeg9eit4WYlmtsSxrvdR7AZshymNTW6lY9gHntAQAW6zB5DA7rc7h6zsDqASRrb5segEVhxBgLD8Aeftgb9BCNCcJ1PAA9T1Il+k+Yrc/mYWLlt1rcB7n1afAAltaq1wPor+QB5N+J5Xd6AJXnn/MbLXqnB7A/7GHkmlM7wnATdGwVgMUDqC2lQW8lD8CUIJzM6hxA9bO+GRsegOGAXsUDsB8mTWtv2x6AWWHUME3bmu1l7r72HXbhAaT4efpZ5l8mCAe9sHkPILc+DdDXNB/7DsOp9ADUXFjhAVhi6N4cwL5nDaYz9HvBnOhfztEC+mqe036lMAxyvQna1JWQvxZCeDyE8D7m798QQng6hPCuxX8/vYn3SuTzAGx1AA3rc+iwTDIPwGSZ5B6A4YD2egDDfsD+sGeCaR5eZR5AI2ewggdQfSdjOGGLOYA0n/l34vnzNdDldD0PwLAGLQ9Afv50VnU/TeOfR2CqGk7b9QAOpjPfGpTjMSaZrTH90gOwoOc2QZvyAH4dwN0Kz1/GGL9y8d/PbOi9LO15PABjbHI8myPGZZwufwb/jqYHoPPPF8m7FT0AI38qlrPmAJK1d8U9gAZipUoQWorl0kbM38m+o/QADHO65/AAkle45/AAegELxa3L6aoegBVZlXsAluK9g2lzn1Vj9HlVXg/AckB7PYD9YR+Dfq/ywowwU+uclvDpI+UBxBjfDODJTTxrU+TyADLLSnJNm4gVq0W/mgewv6oHYOSvvocO09xlDsBkfRb8gHyHaoKZpjXOn8FRywOwzKnbA/A8v5IJa4X7rj0A/bDNZMgj1w4PYLwjD6D6HoazZfH3Ud87p0fTA7DQ14UQ3h1C+A8hhL+z7ZfZrM+mZaLBNA8K6zN/hvSZVXIAXg9gz+GRpGfbLJPc+jRYe5n1aXF9awVjtT7JNeDfcbiC9bmKB+DNAewPswPdYNGn73o1egAaTLPhARjkug617iAHYPbCFjJRfQ9DdGEBM+31gnNOj6AHYKB3APiCGOPfBfBvAPw7jjGEcE8I4WwI4ey5c+dWfqFNSzc9AEDevIeFpVTx8wvVsD7NCcjVcgDJNbXFSjMPwHBAN6zPHXgAZi/MYL01+bfjAYxX9ADSuKweAIAFosS2BpX1aeBf0wNI34nn93kAZW1I/p2kz1RoPq8HYIVnz9wewJLf7wEcKxRQjPGZGOOFxb/vBzAMIdzM8N4bYzwTYzxzyy23rPxOm5bO49u6ZVJaShW/Lsh5rNGMAjJbn/PMojdguBceQxqX5YBuWJ/uHIDPotf6uHg9gCb/4nBwW582D8CjYBrW5BY8gKX1uQUPoJ7TTK4Fucg9bYtne5jL0MDjAdhzc8kDMHthC5lI49JzZ/ka+zzbY+cBhBBuCyGExb9ftHjvE9t856hf9XGRXNMUOx71+9jr6wd02tijQQ+jvn6YjBv8C8E0xDJzft36XFomJrzx4vlm/knG39ctn/I7a4I8zqy90UA/TBrPNxy4bv4ZNR7lO8+qORr0e+gbYJqJP43L4mF4+ZM82zyG5hyp8fZZe06lAy7nN3nas1nr+ZYxefhTLUZ6j2mf5fzGWg9gsW8M4w8BGPQCbr5uDzeeGon8m6LBJh4SQvgdAN8A4OYQwiMA/jmAIQDEGH8JwHcB+MEQwhTAZQDfEy03W69BSRAm8zn2en2SJz8choPqMk4Jntbg7xv4Zzn/YjxKEVIS5MQ/NfIDVZxe5Z/m/D0Vjtd4/qCHyVw7DCu8dL8XMDLxt+dI+g6H5JzaDp+Rhb82CvLx6AdoevawH/Q5LQ4HaTzpO+RK+JJyocp4NnPJxGHxneexqpLvpwtqifHn/IB1Tpf7UBpTbRT0s30m8M8XFf2ufTZtzqm6BtM5To6q49K0b7Iw37Dfs+3LfuWR/Pjdz8eP3/18kX9TtBEFEGN8pfL3XwDwC5t4l5WS4ExmEXvMt0yLPuwHkyAnoRr1exguFldKGidBHvaXlon0/Nm8KiwbLizJfi+ogjmZxvq7Dg2CPJnN62KWYb+nXl4yns3rubFslMksH09QN+IkmyPTGmSHz6oHujSm9LdKJoLKn96fK1XTnGb8ljUbZgpGn9OY8dueD7QP9D5nOM0yuTbtm+U+K39H8y/WYLDcl9I+S0aGVYbS331zVMq14/mDntoDbDxbGgW7pGNbCVwLgjE8UAuaNZxgeX4WTrCMp96Ig3yz+6xDi6B5DvRxYd1Opnbrdtjv1UpNGg/QnFNRqTq9qvqwyrw2k8IY9E2HT/qMy6LPQjTDgX6gHxZz6pGJoUUmpnMMegG9XnAp1T3rGjQ8bbvhNOr3TYZTPh6r4dTwhA1rUHrOqpJfwcNI/Luk468ALIeJ25KxWpM5v25N5uNJn9M272TqtGSmsXmgO61b7+GTnsE/P7O4B/oc0XOqewx72ZqJ1qTTK5zNI+YRvjVwW5/NMJxNyVdzMxr4PYzqd5tbA9JjMBhCjTUQ+ZMM2eU6t7hN+yw3nExzWngMiuGUP3+XdGwVgMmaJDwASXCasVKHNWk8fPLxpM9pgnY4c1oyzvBDadF7rWFAnqMy/qzxu63JzAOwWJN5QtRiTZZrtjVr0uu1LdAzpsNw2lQwwIbXwOkxHDb2pcFwKtdAmdMYY8tzlhRMGlPTEPKtcecB7Ji2Yk3W4QdjbLIRfvB5GOn/kuUQY6ysQ094ID/QDdak3/LJFYxtDYb9gBCyORXDakuPwWNNNhSMkT+9xxRiyvIw+pwuY/S2vEqRAzBYk6NsPLO5goabNXMGgLwGuZz6LHrnPnN7bbY5nc0j4gpeW1Mm7KFWu0zQSfdt0vFVAH0dwZHH3F2WTxablNEMy6RxsiYl/knGDyRrUhdkjwKowgmhfo+OMoouREljowwMa+D0GJYhnWWMXkZuLcMDFuRWOw8jz2n6215mTWqooRZyS0GUTBpr0MNUQVaVYTsAIhprMps3xg9o6LY8Ua7zN/aZwxCyIsPGxJqJ+yyBOVxr4MurNOZ0YNlnS69tl3TsFYDFNR00kl8Gj2Fgiw8vBdOGHhgv8M9Do6Dl4Q1gWfsgUTP2aUuW5ZbMVLMmC/58nNx4huVhZUrq5igdyxr4wnBWa7Lk18ID6Tat3EK3hIAa1qSD3yLXFf9yPBW//J0TZDF9zo+G00M6VjTc0jBbfAfFcKLWzDunvn1jTEp3HsDmyCT4C8sqDz9Yk8aeBGTCP6sHeuEBaILW8hgs8DRv0njWxDMDsjWZewxWpZqPP/1OGn96tilcQeUYhBAKZU1K/F7kVrI0G7UVDo+hCsP5vDZAD9GU/FoIqOVhGMJ2q3gAaUw+fpvh1JxTPURjlYk0prLeRjKccv5d0rFVAB5LpuK3J3WHVmsys1YBgzXZ8hjkzX5IeQzGqlUzf2Yd2qzJWSMkVfHz7zicNvHSaYzS+ENAXWimjifb7B5rcjmmzVqTtIdhgdbaMehjYk617+DnX6KMKn5Lxb3NcDqk5nSTXlttyNn2WUoaW2Uijan0hEXDadqhgDZKfsvHkoD0WZPthKLRmuznCA5p/NWz9rLNK/FXSeMsAWmwJinrUIPk+azJ2Io/a7UVy/CD3Zq0FvHUc7ola9Jr3abPrGpNWg0Vt9fmgvr6wnDlGtgPdFsYrvYKzR6GTyaqMa3gCXcewObICk9b1ZKxWJMtS8ZsTdoqe5ceSe5heMIb27Img4N/1rKUtHYcq3ht5vDAijkA6+FD5ww2b022E+u+vI3ZYzDmAHoLr80UhqPmyMJvlOu8LiH93+TJZ56w3tph1pZTwbA57DyAzdK2LJlUMVl9xhabXM+a3LzHkB9WalLXaU3mc2pbAyJebXy+KwxnzatQOQAnosSalE7Plzyk+TxiOl8nr2KT672B40An19i2ZhbDqbUGiuFEybWpot+ca6OUvCNnYDCEOg9gw2SKP8/alozV8knv2D6iZPMegwciuBtrsolA0TZv6/kG69AahmvP6Qo5AIfHoCn5trWqW5Pj6bz2Ci1huBKzbnl+S8E4rFv7nNpqK1oFlEoYrkTPqTJB8FsMp5Zca0q18wA2R1ZESSnIWmwy19KbtyZLjyHYchJGa3JSHCbb8JJKjHs+To3fEobL18BqTTa8Ni0Mt+KcWqu3vR5DSyYsYbXZfOkVWqzPrD2IbQ2oymGZf2+QKwCnJ2zdZ84cQLMfkyPMl4pMGcMp9b9Ka2DaZ9ka7JKOsQIwWJPU4SZYAm1LxmpNGsMDK3oAZmuSSH4BvPU2nc2LPjcGa9JpHa7jMaQxqUgvz+FDzqls3aZxpOd749uSNUnxA4bcVpYXAnSkVLsWw+ox2HIA+b4xe8I924G+cg5g1dCsMqe1oTVoegDqnA66OoCNkdWSSYtjtSbzYg3dmlxWTFb/N1roxsKukt+e/LIlaVtNtoyxzBxllD+H43cldWe+MFyOegLs1mRTCTuQXkYPo1V4xViT5RpsIwfQDMPpB3qjBXnPJxNpTBoazpNr87bjoDxhr0wAvGHTRiUZ0HBdEnizZBX8UUMwN29NJsgiYIhNOpNNFL/FmiwLu7jv0IqtGnMArpARwa/NaR5O0DbvIbFmXi9pGygj65yyeRvmHclrax0+ygHaguIqubPUtqDXCxhotRVFgtMScmntM0NSdy8LuWwDuWVVqq01MxhOh7MjnAQOIfxaCOHxEML7mL+HEMLPhxAeDCG8J4Tw1Zt4r0Q2azLWi1N9RrccvNZkM1yhJfx8vYBYy4S1Jn05AKpgJn8O/Y68sZjBmiQTkDavLX1G4/co+fS3QW+71qQVIti2bmVrcnmZij0MR6+BbK2OPHI9ja01c+0zY21FHnKRk7opD9NsmMfdW0F52tVzNrPP0meOcg7g1wHcLfz9WwA8b/HfPQB+cUPvZcliTU5ncwyza+8sCbxBefgIizqdzQt+pU1AOnzMIaNmiEmPTcbF843W5Hxe8BsO9DlVOazNacVvCcNNM35AV5LVGjQPK3FO51V4o/baNJmYpzm1VZW25lSxDqfz5mGihu3S83v2+pbpLGLQK58vf4fEX41NBitM501DyDKnHsMp1Y0sv4Nxn2X8gCEUapzTVCOQGxHS81NTx3xOd0UbeWOM8c0AnhRYXg7gN2NFbwXwrBDC7Zt4N0cWa7I8TAZKq91SkAeaNVkI8kCxfKZrHiYDxZpMgp8EeaBYk6wgM/ylIA8U/vQdhj27kqwUTMbfU6zJeazHn76DqjCy8agyUc/pcg1ka5KeU+47tNagZzt80nPT56yeauKXDafYUqqq0m6tgWI45WtgrBtYfgdln82KfabkMaaF4TRQ+Nueuc1jyOd0V7QrlXMHgIeznx9Z/G5r1O8F9MIKySknvyex40c/KII/bXoAmqVRhgc0C52CRObPaT+/dMUt9wHEBvrBGx5Q14zEoCvjcRxu7fCA7XAw5wCY8IP98NE9YTIEpMX0HWvQTtzrhV0NmRjoMjHoBXuujQvDaXLtDcMZw2plWHCXdNUlgUMI94QQzoYQzp47d26tZ6mbd94+oDVEycARMprOS0tJ9zD6DfSDnNSdFOGBkbLZy/CAdqAvrU/n4dPzHT65tafOaWlNKjdwTeft+LNbySsyAWQhFzVJ6wsP1BZ9EaKxWpMjJWQ0X1xpuX4YTrboXYbTrOkVmmSiBQ3WvTZrGK6lVI387TVQ9tkx9gAeBfDc7OfnLH7XohjjvTHGMzHGM7fccstaL9UOdNo1VXIGLmuyDG/4XWWAT+q2DhPlFrQyPGA90AeFh6GGKxxJ4yllcXuTwJqC6dutyfIw0axJfk41a7IIDyhhu/owUazJ5XhsYbvSiEifVXMADiVZGkKWOfUYTpNZbO0zKQxHhVoBYQ0Kfs1rS/xmz7wI5e6SdvXG+wB83wIN9LUAno4xPrbtl2rdMVtJXYM16Y59lgnSwLk2AAAgAElEQVRIzfUtPJL0e5q/TNJq1mRxoCvWZCnImjVZCrLJmiwPEy2BR6yBmIAsDgc9AdlOcGqHYWpPXY1fS+oWh486p9xhwh0+RYJTUUilEZHeoYEPWuAJhX/QMoQ8SWbN024n+qv3akldWyh0iQxz7rNiDTTDKZ/TXdFgEw8JIfwOgG8AcHMI4REA/xzAEABijL8E4H4ALwXwIIBLAP7JJt6rkSXk4ov3FoKsWZOFIKuWT0uQs6TuHsVPHw5sbLJAS+jWJCP4mvXZK6w3ZrPX7anLw0f1kpoH+qXxlOcvk8ZqUrftkaQ7dXvEBqWsT8Bg0bcOH+6A5g50Lf7s9do8IZe2B6B7zqVMaIaQz3Aqx1/9fo79YfuaRQptl55Dj7+cUy2py3kM9jXYFW1EAcQYX6n8PQL4oU28y0PDfk++C5Vq7eAQZK0tbFuQlTuBCQVT/V4+TMocANdOOR3c5eHA8rdCOkHkpwRZmtPZ3H/4TGbLO43TmLT7a0+MmuOR7tSlwhtA5Rns9ejDZFDIRPUceU5ba6Ae6L4kc6m0OblbhoCaB64op0RuS1qDymtbPQ+j3albhVrbhhP/ncs8jKxU23BoxRAqobvKGpS1Hrukqy4JvElSWymUsUmN341maCsYMak7axfYAEKlrjeh2BJMWw7A2jCPEmTpQC8Pt/RvbQ3K8IAcAmrPqVxV6g/DleNPv+f4cz57UZEtoVgivbQwXDme6rOWlij2NSj7N2n7jE7qKoZZo9JYX4O81kPjL1u66GvgSxqX+3KXdKwVwFCJTU5bSVrdNV03BwBISd0548ry4YEc/qZVlbYtei35xVk+cvy5daAz/LT16UWgWBLrZRjOExZU5rRVtKQkdZkEpBrTd+Z52mE4G5AgvYNTSFTRkiW06QFDtGo9NMOJqPUApLBXG3kG6GG4Et2mGTalh2Gt9dglHW8FIFgyJfwNMFiThCDLVaVtmCkgWA5EgrPi5w+T8jAEZMEHsuSX0ZpMgqlZk+RhIqzBSglIYk5lL6+dV/HBTPUD1xO2K6GyamKdOUz0JHMhpw4PQPbamvHz9FmXRa9W6nrRcM59xhlmaozed6DXHkAqHOMMIWINdkXHWwFIgk9Zq054mmpNtgRZtiYpjwSQCk7aeOnq9z74m4pBN4a9yHCCsAYlHhswQgSN1m3i94bhqByAtNmpHICUUOwFNGo90ns5/pxPGw91QEsHeikT6R1WiGP6t7u2QqkD8HnCXsOpHcoF9PqZMtem5wyq52oN86g12BUdewVgRUsA28kBDAoFk35PERXe0PjLZnYSP3u7lFGQ02c12Kj18Ckhjun5EmpolcIulzVZhHQsXpsrB9BKMmshIzpxz1uTTYhj+oxmfbahr4rSLj1hMQREFVDK+6xs4AfwhtO4lRdS5pTzMDbUEoVdA04mijsfdknHWgFIgkYLsh+epid17dZkWTJvSUCSHsaGrMmyDW76t0eQpQOd5pc8DFphqJj1gW9OXQc6cUENIFiTxc1P1nbQ1tqKOszXaKWgz+mobMehdSd13oxXGiquWg+DIdTYZ2pSt40yqvh5maBrPZSQToFW4z2Mtlzvio61AljlMHEVLRmsSW94oBxPxc9/B+qwkjDoDUFWk7pMeIBN6vrCA2WCMz2f9zC4kJGS1HWE4dpFS7qFXlrD0vO9RUvekEtZ65HeoYbhjEldutaDlwmu1kMynMrEvRqGm9NQ3I3us16vBbbY3Jy2+XdFx1wB8JYMe/gw/JQgm2KT1OEgCTJ1+FgFX7NM5s2iJas12UxyrhAeUF1lH2zUmwPwheEKi96QpB161mDFWg/vge6F4ppzAFSthyATVK2HdqfuZNaGjeZjpb4DaZiJ+8axz2aF0la7gTLGpRK263IAGyYpqUshUCRrkhRk1ZpkwgOS9UlZJt5YplGQ05jcIRfHYWKxJltFQpxHwiSZ9aSuAyJYFC3piXXGK7QqGLUOgMtVOedUgZmaw3BM4Zg3LwTIYAV3Upfw8qQ5LXMS8vOb/HpSlzYuuzqAHZOUAyhbHQNauIKwlAx9X3z8/tgkhWaQDvQSa+wOuTgFWcwBUElmycOowyH2MJwXKVXx+8IDVB5GVjCENSnUVuS1HmlMHH8do99SHUDZIRawwUapdhms4US08AYEw2nqNZyYPI9oyDWPSlGuGU/VA8XdFR1rBWATfJs1WfYQT/z530qi7gMAhMNn6o1N0oIsCVp576itUtcWHqBhoII1SSSZpcNnmRC1WZOp1sOV1GWUqnVOLXcs5M/XrMnSKwRsYTgr7LJEhqXvoEN9rTJBy1D+rNZnSpimoQmhz3Aq9pnT0ErfR4JnVzw2D4BSGLuiY64A/AUtAG1NkiEj7UAnKhrzd1P8rtgk5yqzIZTYSjRZIILWamna8jGEB4w5AC4hCtDWG+lhWCCCmwjDGWs90mdEhUTwu0IulqSuMQzHHW5c+2Wq1sNSSOXZZ1StRz5W6ju4x1MoYW1OexnYYskv51U6D2DDJCeB6dgqQAsaLchavNcZHmCTzLbDQe37MidyAEqlbl60VH0HYU4JQbZVlRrj1WRSmp+jsoujxp9+7w3D0UlmW9FSGp/r8DF5bb4cgDcMV/JX76aUMLHPhDlKtR7evIo3DNeQCS2pW9R6pDFZDTmNn1qDXdExVwDSxmofJtIduZQgS/ySIEswzZx/oPK3D4dBTz5MWoLc6+HQKch6zqA5R7rCaPJz1iQZTlis3yHxDjLM11PmtEj4Le/InbH8jfErd+pOCpmovo+MViv55TlNct0ck47EyuSuJyG3CCWc5pRUwsQ+E9agNsx6tn0GtOfUtM+y8fQWhpPMX+wzRa7L3v6DXiBlFKDXYFd0zBUA30aWin2mf1OteUlsr2CZUIK8fL5UVdq2ZGR+T0y/LcjVHNkFWWppzVVAcvzl5SWJP42V42/GhxdzJORtyDmVMOik9cl7JS6ob5FkTp+R5rS8Kzblqmj+6r15+GE44Nuil73xgUW4QrmFziqnNMxU8BjSmjnyKm20Hb+P05jKoiuppXUJM03vmLAy1FbaVUtrJWcwCOTft0nHWgEMFhuluo6gSZIgUwtLCfJAOEwoON6S32aZqIVmheCnd0iHCWVNsgc0IcgDyZVl8iq+JDM/pxNqTnuC0iZQTAPhMIkxLvI2DqVdeGH9XkAv8Pyk19YPSofY4rDqye018lbHFb8gE0RIZyDwS3MqrQFpaElrTHgAotImPHmWv6i3AeQ5LTvEJn72QJ+3FcygxyuYI18IFkK4O4TwoRDCgyGE1xB/f1UI4VwI4V2L/35gE+/VKFla5AFNCbJwmCThGBH8lOCQ1mdPtibLyt60CfjDpy3IldLjBXnUOnx4S2ZMhYxES6m92SXBp/n1NSAre8m8DRUP5+c0tToeOpT2ZDZvWejVnNL8JQoojU9UMET4QVLa5fMHgpdHIbEkw2k8dRpO0/Ya1HMqeNqUTIgewMC5zwblHEmGEzWnAv+UXgMPUmpXtPaNYCGEPoDXAvgmAI8AeHsI4b4Y4wcK1t+NMb563fd5KC1yFftu/o1ErAjJKfLqvAF/mJBwvAEv+NV7m5ZGvxcQAs/PJwht4Y2KPwjPJywf6TBhkq6SNQxwtRW2NVhak4LCoCx6CunFtL7gxpM+Q62BZEHTStWWlE5jujimcxJcrQcbriBrPZaGU7n+okVPJoFprzCNtT0eAsWkeGHTosJ9KCiY6r1MGI5dA++ccjLBG2b9otZjV7QJD+BFAB6MMT4UYxwDeB2Al2/guWuTlJAjrcOehGYgEpyCZUIJsmbJlIITQsCw16tvJKLGRFkmnhzAQHR9qSSzjCgpBXm4SEpT1iTZktu7BoY8DN1egwo/JA+jzS/PqX0NSmRYGt9Y6MdEVm8LSWNKwXigu+KcrshvnVNqDQYC/3wBGCjDcBx/GlN7ToU1IJUwvwYUdNcrE7uiTSiAOwA8nP38yOJ3JX1nCOE9IYTfCyE8dwPvVUm0TAjrcCBYGpQgLwXTe5jQ8LfZvI3TH4gWOhNrdLqyclWs/TAhD7fF+2iMOH+YSGswIhUGFa7glTYdfuCVtm8N5GK5dT2AQU+wVqe00ta8MCr0KMq1cQ2ovFA9p6TXthjPIJcJIcREGBEhhEpOpRwAFdIxNnVM34HP5fkMpzERRtwV7eqt/x7AXTHGFwJ4E4Df4BhDCPeEEM6GEM6eO3durZdKySAy/lwLmoACsoYTyKSxQfAdySMWneCw6CXXl0tYSuOhFEx6d5u/HX4YOtdAVBgCdJdErJA9XPjDrXoHBZUVlDZ1mEgho1kkrUk+qUskOCWlPW92iK34paSuINfknPJrICoMAmYqKm3CqxJRQARMU87D2A0nMsksKRhin+2KNqEAHgWQW/TPWfyuphjjEzHGw8WPvwLga7iHxRjvjTGeiTGeueWWW9YamIQHFptUEfwSbFR0ZY2HCdcPZDSQr/OjYpmSu0+7sh5rlb9nucxhAMs5pb5D3YveOKdj6nAQwgNUHsY2nmYYjmvVUHeIdShhKkGo9Ymh4slSrUfZ7kMaD5eUrt5t9ADqsBoxp8SdDxKsk8sjcfxcFe1Q82yJ7yy3+6BybdyctpPMWjFeaUTsijbx1rcDeF4I4QtDCCMA3wPgvpwhhHB79uPLAHxwA+9VSYIU0gUqOqyTKvqhNwohyFLOgEhwpvFJSVralZVimRQKiIez0fhn3ppsj19PlJcQRJ6fQqDoh0kTIihZkzQcj5tTqtYj8XPhBwoiKCbWiTUQUUBzumhpHkH2uCKLlsQYPRXTlzwGYp+JXhu/L8l9RnQnTWOScwBtuRb3mcMLo2o95Nxcew12RWujgGKM0xDCqwG8EUAfwK/FGN8fQvgZAGdjjPcB+OEQwssATAE8CeBV677XQjKkkHJNpYQiZX36BLknYMQ5QZbdfSKkoyV1W/hnpyArrnJb8CWECIHEsqwBCdPk14AMARlrPQB+TjmlLWHEycNEwYhTcyoiVjiLfj7HXq8Jh6OMCJPhRCltY/X28kAXlDYVajUiw9LnpTAZuc/EUKvdcOKSzJ56m13R2goAAGKM9wO4v/jdT2f//gkAP7GJd3moDgEZQzQS3EwuWrIJcnqf3GyuCAFJcDMieTQcBBxMGNd0St0HIFcC0zmDCiNewtbICsj68KEPE6rVMWA/TCxhOGt/f64vOwdlZa3PgaRUOQSKlGRuz6nUqoGqck3P2it2POeRVPy20KkYomFucav4+RBTLtepuE5cYxKm2eZPtR5kU0QhT1LmYSSo72Q+x+lhc6I1JNaVqAEAjnklsFT0QyezhHiyUAFJ5hiYmH4VQ5cSooTrSDw/tTpuW5My3ridUPTBTMU8BpngXPAzeRgXP4WsEop+SIijYH3yCUV680oJSL4XEFMIxoaM2kpb9hhomCnA50koGQIUw2kDayCDJ9pySh3oXO6Mi7lz+0zKGXA9t7xQX664jspJ7IqOtQKQUTeEIEtwM0KQRVeZCG9UY6LRABRipfo87e5TmPj0eW8nShGdwIV0GJQOpWA4fg5mCjAeg6C0JYTIsLAmOX7vYUJ1iE3fQfQAyMpeuwcgoYbK7qRpPNXf6BAKJRPpbxR/2SHWuwbLUKuExCrDcPScUjIBpFClc59JYAtKIa2wBtxZcZRRQFctaYJPtTqu/sYfDjnCQsox1KihsuScwQPzhw99oHOwUbmi0duGgMKgy4nvsqGVhBGvSvhp61OcUwo1ZMwBhBDYsBrVShnglTZ1SVD6DhR/6hBL4frFw6RQGKN+1TuI7HE1bzePk5QwhUFXlTanYES4NWE4SQqjNAyYZmrcvuEOaAqVlMZkvUAm8XvbgwBceLnthe2KjrkCEOLDZDGIAaVD9h0XECut2CHtOi5jmbbYJI9YkeBmRNGSAkGkYKn5eEt+avwcP2UNa/wAvQbkhTDCgS6FmMi8jZgELnMAdFiNU9ojof8+3cCvhxi54jo6zwOADD1KHgBVGcvBXqvn2/Iw8hozoVCmIp5ChgHCPiNCuWlMUg6A5JdaijNrwHUl6HIAWyAJUjiZEm12DRWQuSAkjLgIKSTQA7KrbEMBLSsm299BbFvgQawwFZD5eMvvQI0/H2/5HTw5humsXbRkCfNZY+hU0VLil5LG7Tnlwg80f4oPU0Q38OPDCeMZdetbstDpMXFKmw5VEjkJkV9CAVHj5y10yQOwNkVcGmYUsqo9nrrWwwH1LbuTpvFXf7N5VbuiY60ANEghK8ie2CEDH2MPE0Zw2MOEQw2l8IaxU6QkyBJGnMI/A8KccoeJdw3IvEdVFZujhqTDhPOSvF4Vh5SasIeJrLSpHIALs64kdan4dnoWNSZK4Un8FJQ4/Y0aTz4GIDOcpAPaus8YL4xV2ix0l16D2Zzj13IGzBowZ8VRLgS7akmuaBS0tLFwDICAEef5XYcJl/ziBJk53FhBrmPonGXCJAiZzUhtRI6fRkvIhwkb3pDm1IjJZg8T5vDxIlZEZJjotdFrwHXTpDDx6Vltfjq8wfJTh5uYOxMMIY9X1aeRVSx0l1PaznobNsm8MJzIMJyAbuMQg1eqEOxYK4Clq2wTZClON+YOE+Zy6Do5VYZoBlxSl+FnDgfqvgGJf+nB2EMuXNI4H2/J365LkObUm2NoKwzpAhbucOAsbgqDnj5PbVxuDbg2AVxOQkINUUgp+cDl23GQY5rydQB2mZBlqKz1SO+QK41LMIEWOrWCJ3xKmEv0i/uAaPdR1xhxhlCXA9g8SfhkOr4twd9oQeYsdN6iVxKKxt4+XOEYG99mDiutOyaXA2CTui34npQz4FFGHFqCipXyGHEmScuggHjrk/MAGAXDygSf4KQw4nWtB+cBsK0daA+ADj/QHWIrfjoMx3phTA6AOtw4L4yD1vJhteRpt72Sze4zhycszakRkLIrOt4KQLiAhbI+pQtYJEHeBEach7P5DnS2apUTZKkqUygE4ypjufCDdU41JBalALiqTKrWI30HKVxh7SzJ13owMsHkhbjDhGp1nP/MeQBupU1Y2xU/o7QJoyM9ixoPddk5lwOQwnYe6C6ntNl9xhhOrELSDCe2GI8e03FvB31FSKuALAVTuoCFE2SuxJuLHWohGpJf6GbaEnzmAhZWISmVsV6YJtWdFOCrSqkS++pvnIJpK2EupEPVegC6F0YduNIFMtQta94kc/689njsSpKG7iohmlbYTpGJYjzJcLLyV9+JvoCFV3oMdJctiHTuMy6UyyltJrSZaj3al8L79tmu6FgrAG8FJJBcTbsgc3Azqmti4pcKZshe7samXOn5QDs5xaOS6BANJ8haL3e2AtK4BlpVKRkC4orrGNeaiw/zhwkvE+nvze+gJPrZA71U2nSYT0TdzHnrk/NiPGFBCrGSDCduTrk1EFtBkK0UHPuMuQSH3WcL1BBrODFrUL6D7RCrXILThYC2QBqkkJp01hVkBJmHmzGHCVs3wBwOCmqIg2mW34G3PlM8ubA+GUGWMOJUYzHvGtQoIG4NCEtJOqApdMUqSlvqNkpeFuLg5zDiPGJFyqusoLQ5mWBDQHbDiWt1zLWzSM/ot8Jqzn0meIXp7+V4AMJwYpFh9BpI3Unz5zU+Q3Qn3RUdawWwiiCPGFdwTBSOVe+QS9Tpvi/OHIDITyM+SteULRxjBJlDJXmTtGJx3bRd5CQV142ntNKWcP3l+Ct+DYnVVto0QqR6RnlZCA/15eaU9gC4EJN4AQtROCYqbcKr8raCqL4DH3On1mDQoyvQkxFRgi1YZBUTCtVzZzbDiUeG0YaThM7Ln5cTVW+zKzreCkBJfvHhBDrWSMcymcOEEWQW/sbEJrnbn1iII+uapoQojXwqx8QmjYWYPl0BKawBUTiWPsPxl4cbwFvcE+L+A/H5Qozey09hxFmII4PTl3IS+fsT1a2OndBdLinN1s8Qa8DdXMfFt0dMbx8KGQYIOH3hwJXag1D7DGjLtVRvUz2v5NcS95yS7DyAjZN0AcuYQEsAkqDRgszDzYTDTSoSMt5GpbumxYE+3ZEgszF9Xx6Gz9sQc8pZ3FPatVahu8R3kIvx6O/MHejcGnDWJ1upy/Bz8GYrHFo0nAhkWPUO3uvxeQzCPhPbfRBKW6qgNxpOXHdS3tBiZEIAW3AAk13QRt4aQrg7hPChEMKDIYTXEH/fCyH87uLvfxNCuGsT77WQJAgcpJBL0lL8PNyMOdzUyt52DmBGYMQ565C7gIVvH00fPlIFZPV3OqHowYhTMNP0GS78QOYAmAtYZOguLRN9T9GSgFhJ76f4rYnvJcy0ROnQBzQPe5WVdvl8yXDiDSHfGnCGE7fPeKgvlzRW9hnX5pzzwoyerboGhUxwtR67orUVQAihD+C1AL4FwAsAvDKE8IKC7b8F8LkY45cA+DkA/3Ld91qpcgW9gmw73NLzuT435OE24BEr6Xk5cbcnLROKjPU5tQkyZ61ygszBNOfziBkRcpEOH24NuFoJqpUyUB2QVphpGpO1NUU9HqkVBJcodx4m7bwNByRQUENsnQEXCmWUpCMHIBlOlCdcxfQdXiGTh5HCZFRxHZcD4AwnqUMsQIEt5DVo7TNm3++KNvHWFwF4MMb4UIxxDOB1AF5e8LwcwG8s/v17AL4xlGbWlmjAwM1419RnmXBtaqkSe4BPKCYlxbmm5WfGTF9zLoHHx0rp8MAyaWwLbyRBHhUJPwkjzicI+QOdy8Pw3UnpeDJX9k8l+iVrFSCSusqB3koyM2EyKSRV8dsOE+nmOt6z5W6uo+V6wOD0uSInFjXE5IX44ro5eqGNGlLlmkPpFN+B607K8UutJqjxcEbErmgTCuAOAA9nPz+y+B3JE2OcAngawLOph4UQ7gkhnA0hnD137tzag+Mw4hT+GZARIlwhmFeQOY+hbHWc+AEBb2xE3SyTzDQCpW2t0glOnp8XZA4jzkEE+Wpp3gOwti6u+R3WJ4cR52CjHEacw6BzGHE+D5MQKFxYkPMAmvxcrQfAG07cBebsGrAWPe8BlEZE+g4eDD1rCLFhOG4N5DlthVqdSpvzYHZFV10SOMZ4b4zxTIzxzC233LL289iLIaZcgYpTkIVwAmspMTkDTsEA7Ri65pq2LPTaY+D4bYLMWVZcgjONkUN8cJvX2sAPAHsBy5g5HLgLWNhaj8XvSlSPVozHYcS5w6GsjOVqPeo1Y6xPs9Jmaj3SZyR0W0kSGo4vxnPkeVhUD1OXwBlCitLmQ63WOZX5W54845nvijahAB4F8Nzs5+csfkfyhBAGAG4A8MQG3q0SW5XJCHLlyvoEmWzxyhwmHJ6ZanVcPZ9LNmkxfZsgsxBElp8ZD1O0lD7DQV89bQKmHKyTLa7jDge+EyXtkTBKkukFxOH0uWZzI8VatRbX8c3v0njsSltqZ+GF1nJzyoMnmLAdUxFPK3laTt2GE2sIyfwczJTdl0fYA3g7gOeFEL4whDAC8D0A7it47gPw/Yt/fxeAP43UhaZbIFFwnDBNUpCFEnUyydyjr/OjSvgrfg5upiSnjILsPUz4SmPaWk3v5NaADatxKCAvdNfFz7ctAGgvjGt1DPBzSt1GlcZbjj//e/18xVotwwlccZ2ktKWb6za1BjzMlNlnHkOO6UKrodusBzTvCcs5Ay7EdGRzAIuY/qsBvBHABwG8Psb4/hDCz4QQXrZg+1UAzw4hPAjgRwG0oKLbohHjOkpNqrwJRe4woWKZEuqGjn36LA022aRAFjlBbl+byY1HOEyIzTtbwN9YjDi3Bg7rk4OZjjiEC3uYcEqYy0kw/Nzh0KM9Bq6aXG8ex8X0uZCU3eKuKo3thhML9RUa8nEhIKq4TiocA1aByq5nOHG5OU7BcEnjXdFgEw+JMd4P4P7idz+d/fsAwHdv4l1eYgupGEHm4GY8AoUPJ5AVkAx8jIPLqUVCXAjIHB6gLaUJY332egF94oCWBJnC6XMQx/QMHjVkhyx6lTZ1kUc1fvpwGE8ZmWBx+kx4YCCvQWkYcMV13jmV+Edc6NRpOLEoI6Z2g0XbZQd6v9evf88ZctI+6/dCq0OsXoxnM5zGzJxyhhOnMHZFV10SeNMkxXtdhw9btMSjhrhYJtBO4EmuMtCG8E1nkRZkxvVle6BwrSCYEFP6DCfINIqmDevkeuOnMfIJQrvHwCUUuQtY2FqPxTupmD5n3VbvL3MG8hqwiBXOA2BDTDYviTMigBSiaaOGqFqPNEaPRe+FW3NIpjETRpSqsbn6n/S8kr8ar9MLI4rrKMMp7esjWwh2tRNVQSgJslT0w5Woc1WxXGwVIDYvE06Qqkq5qthqvEZBZpNl2mFigzgCtMXNwUxrfrKbJoNA4YrrmMOEK65jFQxjcWuHiXkNVCQWF66woZLSO/kcg+2A5sIhaUy8oeXYZwIwoBqz0XBiQ5XMPlMMJ65rrXcN2K67nQLYDlHJJkmQuTa1rCALfcQ5xEr1d4pfSma1BZ/DxANSVSkd02fbR3MJPyNeOo2JTYg6MeL8GtBKmNuI1d+pOeWVNvUdJKXNVZVSrY6r59kQK1xxHXfnQ/W79pxKCBSqFkPi55U2049JDJ1Sz+e9MM6QA2jwhORps1BZJgfgXoOp3RDaBV0DCqBdoq4JPtd9k1vU6u9t603ip5QSV5cAUK4pXxVb/Z1GP3AY9Bam3BmiEeeUwOlL/BROP8ZYheHY3vJc+IE/TFqbkTlMUgfSEprKd4hNz6e8PAE1xBwmtCFBHOjCYULh7jlUUhoTm8Nw4PS5NUgKiTKcuAr69LxyTO59JuTaqH0GtOeUhXUqa2BF2+2Kjr0CoDDikiCzXQ2ZBCEH4dOTWW1+qtUxd6CPlcPNWlXKXcAyYVpNAEuvpzGeqaAwqDVgWl9UY6JRQ1SrY0BArMzmrV79+Rjbvdwj25oCoDwAvj0Iyc8lmbmcAVO8B6SwGudV0XJBjV/ib1e5VvxcpS61BleLIeYAACAASURBVGMucb+QOxIO7dhnbHfSDe6z/O/l89trsJhTEgHIe2FdEnhLJCW/ONglGQJiY5P0YcIJMocR16pQqQNdEuTS4h4zSi+EAKqdhezutw9oSZCpHICMQW8rGAktwV4WoiGrCOuNu+2qel7b+uSax1H8bK0HFwIS5lRM6rKJew7pZVPaUq0HZzhxtR7cNZhjpoU3V1w3mSo5A6p2w7nPqFoP1nCSvCoiVHmlYaDHXgHQCUjBte63EzXVZ/jOkvkzE3El6hLs0n2YiMksmyADS1RMOZ7qefSYuD40bHiALYChD3TeWqXHTxbXsQlCASEiQQqJOeUOZ8DuMUhhRICDdRJrIIbVCKUtHD6UIaTBTEt+qdUxbzgp+6wFnlDmtGV4KIYcuc/a/KzhxOTagMor8EQjdkHHXgFQh48kyJT1KQmyCDeTYJ3Ow4SyHORklk2Q02faMFMpT9KO94p1AIK1al0DScGwRT/K4UDxS9ah9XBY5lVsCkm9EcxoccuhTd6r4uaUx6zT36FlnSseDMAgq6TEvRWJxRTXsYacUODIXdZCIaUkOHSFAqLluvMAtkQ0mkFOcJYXsIiCzBZ4MG0FhBJ1r+BLlhKl9DhBpixuUZAJxIcm+J41ICGIknXLWNBSVSnHL4aAjHPKrgGTZA6BxoinVsdlrQcAspBKruyVlDCjtLkGhEborkVpUzF0qRiPNoQ2oGAkqK9gOFmbzaUx8aHWzgPYClGCL2fq29aeDBtlkkdTvg96PoacfxMewBKBQhxuRM4jjamd/FrkSZgkpyeWKcerGbQEi4mXLO7lO1KrY2+bAB9GPJIJTvZ6wRnd6jh9phzPmDnc0pi4ZnOchc6FmLiKeE4muMT3bB4xn+f7RjecqAJBrqV4PuacX9xnrUS/ss+MYbv0GVeo0qkwdkHHXgGQiBVRkNuHg01hGAWZg5uxbXYZhAhzuCWMOCWYXJyRuoBFst6opKts+VCCL88pp7QlJFa+eeWiJSGpS1qfXIiGSTILnSi5NaBuyJrO6MMNYJR2XVXqm1OunQXfa4iQiYFzDQackmQuCRI8BilvQ1XQi/uMaMnNK4B252Cu1iO9gzWEukvht0MkYkW0JtvW3iqHCVu0xLqaShUq4Y5z1YMURlwTZK6ikRJksbKXmVPvGpStGrjeRMAyUZ2PSUQxCb3cpc6SdFWp3VrlFEw1Tqq4jj7c0pj4RmRcWK1tRHD8lOEkz+mKa5B9hxgj+JbcNOqGr/VI+6b9HSh+znDiaj3Sd6Ar+mmwRdXRlNk3nQewHfJWQI5qQZub+FnMt3KYUAUhUh8aUsEwQkMV/UiCTLqygiCTGHGpF5ATsVIf6PP2YSL3629bn5IXZp1TsapUbB1hUzDVZ6i8h7DGVHGd6FUJiBUGpsnWGTDFeNUzszUT12CRpM0s6NlcUPK1hU6FaIR9RnwHdh+QxXW0J5++g9WQA+jEumRc7oKuAQXQhnVuUpApjHhdtCQcJpRSkppaUb3c+eQUbXGzISAy/MALMo0RTwf0BhAr0hqIMf2cX84xAHSYzNMKQoWZGpPSAGNxi16eVAjGe1XleKq/MfFqll+Q60YISM4ZAM05lYEEzBq4YZ18TJ8urpP46RYn0r70GEK7oGOvACiMuFitR+QA6oQowT8irU8+rieVqPv46SRz9Zn25uWSXzW/MbyRxsTHMrn4M50D4KpK05jz8ed/K58PNA8HrSo2H0M+JklhtIvr6DVIGHEqqcuuAYMRlw+r9pxytR5DAuorJfrpPI++BpQXJqN6lt9BGg+bt/GCJ6Z+uZaSwO1rPPl9OSIVTOcBbJWoBJ4El6OuktPgdTlP/lkZBUQkvxyCzzWDS+Nst5umL5xJ/NR4WMSKE6VDhqSkw6Ge0/YacJeRAM0EnnxY0V4bd0ENX1wnHNCE1yOtAV2pK3gAfbrFNisTVM5AsdBZw8mYAxDXgARb6DJBwYNlFFAZqlQMIQLezK8B3RRRkglWqR7FHEAI4aYQwptCCB9Z/P9Ghm8WQnjX4r/yusitEoXTry9tkLpvEoIvhoAI61OMPxvhZjVGnICzSfhkCvLHo4DoEnU+x8AX13FhLKqPe3o3xQ8wayCFgIzhBwojrhWmlfzpM16MOB9P9lufLZlgLqhJ4+GUtrXAUbToJU/YaTj51kDLnTlCOoTSruaUV9pemaD4qXs9dkXrqp3XAPhPMcbnAfhP4K96vBxj/MrFfy9jeLZClODU1qQAN8vdZSm8UYeAprklI4dDSv40JsnSoJta8dae9zChLHQOmkbD3+gLatLzS4y4tAY1VHbaVtoS6qaJQBG8Nsr6nOseBtVim5tTLoQiH+h25BZdXCeFBdseg1TgSIU25TqDdkhHWgNvzoAynObz6l4Pij9dwELKNbfPBnS1tG9OeX7qKlJpPLugdRXAywH8xuLfvwHg29d83saJxojz1ueQtD5T0tjmataHlZF/Po+LugH7ga5VKFpDTDW/8aalip8ugJGsW2C9NZDDcO1EuXTTEmmtGvitrYjTc7xrYC32S2Oi4ts8PzEeqSMrYTiNhQ6xpMfgnFP5jggixKR00uQtdJ/hJKLnWmvg8wCkfbYLWvfNt8YYH1v8+9MAbmX49kMIZ0MIbw0hiEoihHDPgvfsuXPn1hxeBuskkk3W+LNY3Ue4mrUgk21w25ZSEmRfzN1n0etJXcfzmeI6tmhJig9L7ZeN4QTKQpfyMMuksS2+LRXXifFe8gIZCYLo8Ri4nIFPafcCV+tBKW3doicNIavCkNaYyM1pCVRKrrkK/TQmqrjOYzhpa0CuMbPvd0HqpfAhhD8BcBvxp5/Kf4gxxhBCJPgA4AtijI+GEL4IwJ+GEN4bY/woxRhjvBfAvQBw5swZ7nlmog4T0ZqkDmgh/kxhxLWeKdx4JJim9TaqxE/1cmctHxIjLsPZvPyJJx8/QHtJVJJWig+Ta6D0Mqr4s8NEuZqvxIinVhPynFKwTv4woapW94bCgU6EHzSlHWOsUUJc87vEn8ZQj0foWyPtM66dCEArAO4CnIrHts/SOEnYqFART/GzXh5TXOcxnCRo8C5IVQAxxpdwfwshfCaEcHuM8bEQwu0AHmee8eji/w+FEP4cwFcBIBXAponCiNsQIm1LQ0SgEElma/xZsnzSeym4mQ8jLlnobctEho36PIb6Ri2jVyVBCqXDpAEDFS+0SR5DNPGncdLjsYdoxjNJYfRwcTxr/G4ym+P0Pr1FKcTKhOmln54PNA+0yVSOb6cx1OOX5lTwwqRWDZRRQF3iQ4b5BC8yjcljcVPFdaJcEzfXaWvQSjILz98Frfvm+wB8/+Lf3w/gD0uGEMKNIYS9xb9vBvD3AHxgzfeaScbpS+EBKqkroYCoBGebv8aIUzkGKQTkgJuRSV0RzubrQ8MV10keSTWG5pwGJfxArgGVsBy0cfoTKQE5aB8m2s1MpcVt4aeKkHiFEYgqV+8ay+EHoG3RSzmMNIYlv2QIUTF6PcfQhIH6DKeppoQLi7tuNSFY6GtDfZU1IPfZFQwBrfvmfwHgm0IIHwHwksXPCCGcCSH8yoLnywCcDSG8G8CfAfgXMcadKQBv+IES/CVkkbcmqXg165oWgrOEUPIwzWkhyFKnSBKCKMDZKDyzCGdjMOK8q0wkFGdVO2u6Z4ovBETmDKTrFAmor3SlZfp9Q2lPU4jJMacipJCO0Uvhh3ZCUQ5JAe1EvOSRALRXxfV74vjFDq5EzoBrmd02nLQQUKm00753JmkdUF9tDSj49FUdApIoxvgEgG8kfn8WwA8s/v2fAXzFOu9ZhyiM+NLSsKETxCKkgYRmsCXwpJxE+j0NWRRcTbK3jyPJPI84xY0nK67r9/qLMUmuctsLE9scEHM6FuaIhPqKl5FQMFAtBFQo7ZS4F919R3x4QBfjWUNSwKLQTAsBFXKqhYAoi956c500p9Qta5LnXL23NJxkz7lUGFOVv4cL02njd1LehrxHQ1mDZDj1MyPnKIeArnqiY43VonEtW3Oe/N9S3xdrfDu9gyxMM/JLPVnSe+nDhH9++65Vudtoemb9HaZCp0smsS4V2JTPl5R2Heab53PKW/QhhFbxm5Tor34fGI+E/w5U3YAYrqD6ygjFeBQCRUqI5uOu+ZUQUCnXbKdLMhcmKW1+X8pNDolcnlRcR9SS+Irx+DWgi/3kRH8+7iX/lfMAjr0CIJOui2IN+n5cmh+Qk19W6xNoW9zT2pqUcgA2eF3FT/f3l63bdrxaE+TyO3CHGyn4kjUs5W2Ew4QsxjOugQop7NFK2Oox1A0CJQVD9ZXhrOF+IC5g4S/9oZS2KBNUcZ10uFGoIUFp04aT0xNW+Mt9k8YjteOg7yeQUEDtfSbVGaRn1vzCPtsFHX8FQB0mQrEGiREXYpkURnxcx4cFwWnEShWLvnX4KB4Gc5sTPx4qByD1rm9b3BrMFGjPqZaApIrxNtE6AkAL1qkhsUqFIVUm1/zGOoM0ThqyqCjJYg2k+DZQhtWES4IoqKwAWfQqbTokZZhTKjcnwjTbCkZscUKhhgSPxJszAIrvLBTv7YKOvQIg0QkruGnVs3jBsRYhAZUFQhemcTHx0vXVLB+qs6RUAFONJ7+AZSwc0FRxndZlMY0hH4/XA5Au2ijHI4WAgOS++9bYLRPGVhYAfSMY1+kSoC3uVeZUahCYjzv9W+On5JQ64Pq9gF5hOKlKsjSctH3GrpmQ6M/GLzUITOP0KQwqVNl5AFslUvAVSGTF34a/8RdJ0CEduYmUI1zBoIYka2xcCDJ3PwGwPNDLC1jYVhAU6kYKGVFFP0wfd4A+0GVrWMgZiDF0yjrkFEaPTCiKl4UYW1mk57QOEwGmScm1VIxHJmlFi54O6Wg5BvJAF1A642L8+bNKKg0nLRfWDvMZvDBj7q/mL26u03poVc8t9lmnALZHS/hbGX5QwhWNBKcuyE0IomZ9luEHH2qotnwE9AO5Edl4Mo0Q4flpmKY+p801kHIeFX9zDbj57Pfa/FKtRzUmOgTEzmmhtLUiJFZpS4dJK/4sX+OZePIxScV+QBGqFGP6bX5JJsgWJzO502VpOEn3AVRjCvQaS/uMqreR9hkRtuObLrYNp7GkJEm48tFuBnfVE1miLkEiCYx4dfWfIMiloClJ3VF5QAtdFtPzrZedpO9gvZ+g4qcx4hL8DSgLtXSvqrz9SQ1XFJuRO5zF4joxnuxN6tqgwYnfk2MYFRhxqdNlNX5uDZjxEBewSPcTjAZOJU8V1wkykcbk8oRXmFOPwigT91NFhii4slTYRa2B5Dnvgo69AqAx5VISeIFOmDYPXClRU1ZxLgtmJEFu83sFn3XfB3TbAo6/xulPbQpgibop51S26MfFhS1quGLa3IzSxdktnL4hBEQX+/EKg2xzIIRQqDWQ8jB5cd2kDiN6cP3CnJItSOS+OG1+PSxYwpulw62K6dtyBgBhOBkMIU/epuztY4H6AjRUlqKlx2DbZ7ug468AiEWaCq5ywoh7FqkFN1Mai5W3P0l909N3IAVZgPyR8XAhXJGPo/qM3OoYaHoMcusIKpxgaFtQ8jPhh+o7tIvrpPBDGXOX7idI34EOqwnxZwKCqOHu03M1hVEbKsbEOoe7d0F9vfFtpcipZTiZ2mvkXqQWAipzc5ZQazNPBcheHrBc27pBoBAqBpqG01G/D+Cqp0GxSIBu0VPJI1GQiwNXw+lz1qq1j/iypJ3fvPOIGiMu3U+Qv3diFGQKI14VIfGHbZtf6rJIhB/mUeybTqF0pBL7spBKs+hLjLjUTqTid8pEYUGrCoNcA5/SNqGGipwBJxPUzXVakVMLKivcTwC0lbzUnC79vgm3NuwzAv6tKe30mYliRHBrcKUuhAeuCQVAY8S1A71EJ4ixzNIycSd1LfFkKmSkIETm6TBZjMdovWmCzHVMlRKuQNsLc0EWhcQ9QClt3iNJ/C6valDKhJxYrxArjnhyYXHrxX4E8klqx1Fbn2UeRpGhIgwnzSkl11rYrhw/dz9B4id7+0j7bNo0ImT+5s11lvYgwFIR1eg/LtTKgC2kOd02HXsFMCIOHymWCaDVl0WyfADaY+A6XVb8dGzS2gqi7sliRIgsrVubK6taw9QBPeOLkKh49Vg4HGqMeAGXk/MwhAfg4pcteg6xwiM+/BDEnE9TGGWX29TpUsvDlBaui1+bUyLmLnW6pDxnyRrm0HAyTHN1w8nSIBDIPACDB5OPO/37SjaDO/YKgA0/KPHkMtkkCXIJN9OSXyVGXE9O0a6pCvlbPFe1JgsIn17kRM2pIZ48Lw8HKTzQQ9nbR57TMqFoWAPygJZgmr6kbo4R1xRGCVfWFQYTflBCRs0DV28FYUUBpTFZw4Ikv3A/QcVfKAwtBNRjegcpB/TUuA/4feY0nI5wO+irniiMuFTQArQRIlNFkEu42XQWWWs4PZ/sasgdDgMa1mkt+jGHH+bNBKQL1ilYn2RxnbIGo37xnYU6AIDCiMsKo0SgpM9K4QeKX0KsAMsDWut0yXphqtIuwg9KCMh6HwBVDyPlGNJ3aNUNOJS27uX5wBatfWbgB5ZzpM1pKddayKhcg/p+gs4D2B6FEKobtYrwg4pOMKIrABpuJml1Dm7GF9lUEMRkTVpimflzdYvey58QKEtBroqQNIXkOKCLNdAqJqkQkJq3abRq4BsEpu/gghQ657TNL1v0JU5fDT9QayDkVfpMuEJbgxLqK+2DFgrI8Py14NOGW9+ApSVvDQHVa+ZcAy2HsQta680hhO8OIbw/hDAPIZwR+O4OIXwohPBgCOE167xzFWoX/ejhh3bISDusSjSGXcFYCmCADCOuhYAKS0PvHcR5DDbXd6YIMnfJu7gGhZK0Ke3SutXCcGVSWh4PhQJSi+uKNeCRVXTexjqn2v0E1AUsUkdWqrhOKt5L72gW+/F5ofQdWqFWxdCy3k8ALGSiGE/1e2XfJPCE2qW3ya/XbjTXQLufYBe07pvfB+AVAN7MMYQQ+gBeC+BbALwAwCtDCC9Y870uasXQtaQugbtXQzrGroxAG242nsmCXMZ7pcu5q/E04WaWCkigjUBR4W/l4WZMSqcxyWG4ttKWIYXNNZDuJ6j5W0VUPq+weo4vpMO3d6aVtpakbSX6lZh+qYR1pdeMuYvw6TLxPbV4eY4kM+mZ0w0C0/jz4rp6nynw4zIHwM8pvQbe3NyRTQLHGD8YY/yQwvYiAA/GGB+KMY4BvA7Ay9d5r5dKSJ5moVMQPr2gpRkfll3fHsre9ZIgl0gm6aKN/PdlLNMaAtIsn3I8WisL6gIWfU6pkI7PA+Cs4cTfQIYp1m0bsmgNwxVKW00C25R22Y5DhR4Xt7KlBoH+OVVCQI05lQ/08kYti0yUxXUif3ZzXeJP7yXHMyjkWsuFFXOqXtRUdCXQ7ifYBe3izXcAeDj7+ZHF73ZGZIxes2TKBk+K9TkpXFmtCKlVDKJ4JEBmmShwthJupt9PQMcytfsAautWia2md5dzqhUJlW0CZAudgoEK/GU4QUsy94viOm0NuPiw1k55ulrOQEWgMNaqFyqr8bfCfJoHULYTUTznZs5A2WdFotyaWLejgJr7Um3rXnQl0MKCuyD1TuAQwp8AuI3400/FGP9w0wMKIdwD4B4AuPPOOzfyTKqSVj6s2gUqouVDYL7lpHG7aEkLbwDLA3fTFv0S1ZOSzLbDbRmS0pNZrfYUmsVNdNPUUUCeOS3j1br1CVRzs9fr1w0CpaRxem76XPW9lPBDGR825hgszexyPi1cAdCwTo+Sl9pTV/xNQ0hqTpe+Q2k4yfBsOlRpVdrWObXytwwnRWHsglQFEGN8yZrveBTAc7Ofn7P4Hfe+ewHcCwBnzpyJHJ+HWhhxJZY57AccTnK0gd/ysbiyMUaEEFT+5Z23ZXzYdkCrlk/LWpUVzLIIqWkNi3OUQVlni06XMj9x+IjFeG0P4PQeL94l1Ffz8vLah72BzRqu+G1zWoYTLPfXpvHkn+P4ywtYtAaB6R3W1hE1f54zMMh1qbRFD4AwtDQ4N0CE1dh7MZqGkHkNjLUbLaNA4d8F7eLNbwfwvBDCF4YQRgC+B8B9O3hvTSXcTOqbXvH74s/VxRZ265Oq1NUsK4AID2g3DyV+Bc5WHuh6l8VCYVisySwHYBH8AdWOQ0ANUbc/uaC+wv0EaTz52Mdaa4oeE09W7+xtWvTWdhza/QRAktNiDTSL2ynXZVLXKhPpO1gNp3o8gkxQ+0xsEFgaQkqMvuRXQ0Bl7kzZZ7ugtd4cQviOEMIjAL4OwB+FEN64+P3nhxDuB4AY4xTAqwG8EcAHAbw+xvj+9YbtoxJuNlWszxZO3wkp1OLPLbiZwVXO+TVBblWVGi6pr8bRTOpyYyqL67Qb09K7lwlRWSGlv7nbcZThBzGe3GuFBS0eQJ7U1WQijQPI5tR4B69m0bOtIySLOLOgNSMivTuNJ91PIOeq2v301RBQvgYGpV2NPTOclPGk59bjUWCp6bmJv/o9F9OnlbY1Z2DJw2yb1BCQRDHGNwB4A/H7TwF4afbz/QDuX+dd61BumWgXbST+Zg8RJZ5MXMBycsRPbRm/VT0SIt6rjafis4UfasE31hmUGHGtJB9oxtxNSeNSac+Udhy9dv8m0brtVRhxaxiupYQ1mSD4AT08YG0FwYcTNKVXeiS2A7quM9AKuxw5g7bhZJ/TEXp2z9m6z1o5A9lL4pS2VlzXggYfVQ/gqFCO09es4fQ3V1HRoAkplLosAtkFLJmlYUpAZq6mjd9mcZedIk1J3X6vxW8Nq5nWYFCG1eQk8GjQhplq42lcwGKd0+nyQLRYqwmBZel0WY27eUBrl7C3Dh8lpDOujQh9DUZZmExrJ5Ke5YFbD1s5AKdcG8KC6blLfvs+s3TpbfLLFn3qSlAXgh31SuCjQjncTEMC1Px5bFKJ95aIlbHQZRGg4Wba/QT52C2WVeID7K2Fp3ObIFffIbSsT42/nWOQIHzL8IN2P0H1fN99AEtrMosnW0JAWRLVclilOVVvx+IOEwU1NC7n1FhcpzWnS2NqyYTKb1fa1C1rXrnWxgMU+0yBmVZ8tsQ620Zd8UqWubkrHwK6JhRAjmawWD4lRlyz6EuM+HSu96IHss1uhCDmsUxP+EGHp/kFuWHRW+Y0Cz9Y1iCHIGrWcMmfPqP1ocnHrsWfS6Wtt44IBb/ikXCtIIR7kPMYvXY/AdDMk2jN6dJ3aCltLaxWJu4dhpO6z1prYNxnKyrt6Uxv617y58+hv0NmOHUewG4oh48t2xzY4tWAntSl4GYmSGEWczcdJhk8zVI4lh8mkiCXN3yZLHrigNZx9zZYavpbCd9zQxYVSCGQF/EoHgaBEDFBEGc2mWhjxHWLPp9Ti0Wf12JYLPpcaS+TxjaZSA0CtX02j80eVxp/PhZ1nxFy7VHa41l1Cx3faoLxGIxhMq2gcxd0TSiAfNJry0dZpJTYrFu2OtAJ3pi+nmOo/jbOvoNkDZc5Bm38vV51nV9uDefj5L6Dl388bR5u0hqMiDXTYZ0OREkrDyM3pyvjybqH0YxX+2VIPxzyOV3KtdVr882pbQ2Wz5/NbfusGovNE26tgQHMkfNbOsoC+T7TO8o2xuOcI62gcxd0zSkAm+WTucr1/bt21I3a5qDElM+UNgdUzsBh3WrWcPpMfROSKZ7cayFWrHNqiicTXps2p7MCI26Z0/w7WDyA/Dto1nDFv5Q7k1FQeD2c15Y+48rb9AkYqArTLHIGVus2QSgtqJu50UL3es4leEJpEEgV10n85c11tjkNDQ9D4982XRMKoBl+0GOZOUbc0rKVst48GPHpTG5c1oab6a0sKr7MklEaTg0b39kQT87iz8sYvQaVLSwl7Za1DHOfxsiPf7kGMUZMlGZw7jklMOJafDs9N73HE36YzOX7CdJncn5Ay6vkCBTbnJYdaLXCrnY7EUtexTan1CU4rn1mlYkaMSiPp/pMll+0oufmxT7rPIDtUn64WWKZFVxuXscxAT3+nD97osSTyfiwgT+v4hQhkQXCRYtvV+9oV+paL/S28o8da5BDELXmdOn5QLV5U6dLU97GOqfFBSxq64iiE6XmkfR7ASE0azc0dEgOV7bEk0cZAkVrELh8vmMNBr749jIMZ5vTtAZ57sm7xhb+cYbS0azzUb9Xz03dEkU8K3LDSc+FbZvWKgQ7KpS3/rVMejqYZvNoQriUcDNLp8tqLLkgG5LGRteUah2hCXJucVtc0zw+bEt++QQ/L1oyhYDSZp9GBFgOt3YCT6sOB5q4e491q3W6DCFg2GsqSa1PPKm0lTGVISO1GK8FnpChu7UHlkJASmVyGkv1GVtl7zKko93TkVBASUlqXXrbISlNCTcKHGdyg8D0HfLQLyAbQtuma8IDGGXJsvT/PUMSdTybO/lnldcwnYv8e0VSt+Lvm8az5Jdd8RDs/EB7jrTww2hAzan8HfLxVPzynHr5AeBwNjOPJ3+2eQ2mtjmtx5PzD/nnp8/Uz585+a1z6uQ/XGHfTGbRt2+m80WPHzs/kOZIShr3G/yHypyO+sTzBZlIn7HKRPoO7TmV37FNujYUACH4FhTNeDq38feXm92CWW8J5lRua0sJpsRf34OcH+iKYO4VB67Gnx/Qh8Y5cq3BwiOZz5dhOIl/L5ujw9nM8PzFGsxsc9pSwgo/qeQVV7+a02rsh0pzuiX/cjwhyOGHxuFjmdNBD+PprMFvNpwsa0AYWhb+w1xOJdQQobTluoEeesHOn95h3ZctfsN33jZdOwpgtrQC0u8kfmChAIwbxcPfFsyZz/IxCuahWzDTZp/p/P2m4AMO69MzR/nh0LdZ6LVlZTwc5vMqZCEeJpmST5+z8LsOB681uYrXVnhV2gHa5rdZ0FajIPF7DC2rYbNc41n9f4uFnhs2ln1z6DCcSkMrH+eVoGtCAewN+pjNI6azed3nXw7RVEJ+snw+xQAAF4FJREFUOM35ecFPz6r4Zw7+pbVnG0/Gr4QH9gb9ehMeTmeqm7mXKYzDiX74VM+f1fyAfJg0xmNag2xOF++R3P0m/9zBP1sqMIl/2FQA1Rrw/L1e1TCvuQbKnA57jefrh0nfLEMVf6+e++WcynJ6OK3i86Y1GC7ltF4DaY1Jfn48+601kOU6l4n0f5NcTzxz2s/mVA8ZlWswWNTgXCm6JhRA0zV1xJOnTld2BQ8gIY20LouJP/1f8wD2Bm3rUKK1XNlFkplrT13yW7ywhldltFZX5T808O8VYbtqTn3x4Y17APmcTucYGQ4flwfQrxrmTefRzJ+e7Q2d2viXazCdzTGPXk/bF9q8Gvi3TdeGAqAE0+vKbuPwmS2TXxJ/CKHhalpdU7dgrqowHPHqlCQHlocqxw/448PuePLUmej3HiZb5m8kmZ0KBlBgo5RcWw0h55xuw9BqgSGMYTjPnO71l6FTU96m7zO0tk3XhgJwCvLejg6fxkYRXGsgCZotZwAkQcsE08K/RUuGmlMRweH1qvpb5nd6bekzqyhJwOGFOfj3hv2GFzYa6DmDNHbPHDUsekHJ516ez9CaZc+XDSevJ9zgX8ET3rRMbJvWensI4btDCO8PIcxDCGcEvo+HEN4bQnhXCOHsOu9chch4siX+PDHGMgfO2CcRr1YFs4gPe+LJJkumFcv0xZMtsc/07LQGWs4gPduUMxhS/MLhk8WTLXmbfq/qvmld4/S8Zt5m83Na809sRkGet7GMP43dn7fx5GG2ww80LXrLGuRelT93ZsjzOGVi27RuIdj7ALwCwP9t4P2HMcbPrvm+lYi0TJwWvdWV9aAfmuOxxZPN1me/jA9vPmTkQT+UXs9AuNIyjT/xm5BbBUJE489j+hbrNv095/cgSkw5g0EfT1+eZPwGa9XpMcxjVZxWYdwda+ZGz+lKfh3klsUDqD5TgQ9mi5sAtTVYN3d2ek8+Uq82D2DdKyE/CEB0I68GIhOKzhDNJuFsOU7fomDS88ZTW84g5wdg3ux5eOCG0VDk31tBIQHLOXIdPo6cQVLA+TvF5888h0nPHOZLzytDLlb+VZLGFoUBLL+z5ftW/EuvR/baths6zWP6FkMrjcm6L9Pf10ka752y81sS99umXamfCOCPQwgPhBDu2dE7a1oeDpUgWwpmKn6jRe8MMaXPVOPRk19Nfrsg566prW5g5uIHFptxYj9M0pxaPYZGyMgbVrOs2cR+mIz6vSa/sRYjJb4th9WhJ2+zAj+w/M5WpX0wqdZArTNYKOg8dGryhI38yXDy7LO9FfcZYM+d5aFTC38qcDyc6vU22ybVAwgh/AmA24g//VSM8Q+N7/n7McZHQwifB+BNIYS/jTG+mXnfPQDuAYA777zT+HiZUhwvubJ7A1mQG/xTD64/t+htrqbdkum7PIYymaXXDTQLx/R4dTGnHn7jYZjzA8bEfeYBSHOUYvrjRusIZUzDpgdgiuk7QkZl+MEST875r9/Xww9AsrgN8e1h00K3zE/On74Tz98n+G37xso/avHrc5qH4awyseTf7Jxum1QFEGN8yboviTE+uvj/4yGENwB4EQBSAcQY7wVwLwCcOXMmUjxeKl1Ns3VrjX06ESXp76vEk805g0G/iD/7XF9L24LEa21bkPhNPVYKpaoVzIycCiB9xhUeKMJ2eky/h/MHU3eYL30PV0zfpFSdSrgI23nCfO5KYAMMtPp735zLA5Z5Eo+XN55W1eHTuXzxUsXfVMIWGCiA2jDQcgbbpq2rnxDCqRDCdenfAL4ZVfJ4Z1QKpiXhClRxW0t8OFV97iKe7OW3FMwAi5j+zB6uaClVM//M5wG4n++P99rDaiscPt74c9ayxPWdZwa5HjQPH2/i3ssPyPsgL3B0z+k2Da0Vng8YobjlnB5xGOh3hBAeAfB1AP4ohPDGxe8/P4Rw/4LtVgBvCSG8G8DbAPxRjPE/rvNeLzUgf8aSfGAZmwxBv7UnwSit8LQE4TNDCofO2OewmTPQ+fuIserkaIWBAst4shkGmsWTZf5mPNkCcczHkz9DGlNzTi38M/saLCB/dthoJRPLjrKOOTXCQIElVNYFA/XIxHSZY5BCrQmn34B1WmP6BlgqsMzDWPn3Wvy26ur5PJphoMAyH3mkYaAxxjcAeAPx+08BeOni3w8B+LvrvGddciNQEn8WMtKQTpX1ZitQWfI7Ypl9p7W6SJZ5PAZgVYvb1mQLqCxbKx7bM55Bv4f+IqZfPkN6xypILK9F7+Z3WKtA7gEYY/Qp/GDNGazoAWhhxPQZl5w6LfS9QR/PXJ66LXpLR9n879Y1aM3pUfYAjgrlKCBPLNOKlkifOZz4Dtw0HhP/oJlQtCFKfK4sUFncHljn4cQo+Dm/J1ZqfH76TMU/Q9/QZKtWks45tdQlNMdjzBn0q6aFB2OjEVHKqeZVFfUnduSWUSYa/DaES+kBWJFVZkOr3GeGvE0DzGGE1h5Oqlbw9jmyKdVt0zWhAE4s0AYHkzkuT2bYVxAxvV7lmh5MZziYzOrPi+8Y9XEwrZ4PQH3HiVG/Hk8+Ruk7HExmOBjb+PcXZf8XD6fm5wPAhcMpZvOo848Wczqd4fJ4bvq+QHU4HExm9c+253vWoBqPmX9SPR/Q12zfuQbl80+M5O2W/v65S+PqfY45ssxpet7BZIbLljXI941hDUb9qp9+xa/LBFDN6eGkGn8IuteT1qDeN4Y5auwzbQ3SGqd9rM3p4js+dXlcf157PpCtgWGOtknXhgJYLOLlifNAH9s2ClAJwuUF/7Af1IueT2T86fPaeC7ngm88cD93aWLiT4L+xMWxbTyLv18e2+Y057cIfgXVxXINLGs27OPy2Kbka/5ss1uVsHUNqsPKYRQMm2tgndPzBxNMZgalnR8+hgO9XrPFd9bGH0JozKll3yT+NB4t1FqvgcMQujyxK/kTwz6m84hnDnyGU71mRoVRz6lhjrZJ14QC2B80D5+Thkk/WQimyp9ZezYFM2hYk9qYToyaCkPjT3//nFEwTw6b/CdHcnw4Pe+ScU5PjpqCr40nP0ysSjhZ3AfGNT4xGtTPH/SC6o6fLJSw5TuPZ3NcWBwm+pxWf1+ugVHJX5yYxwPY16D2MFacU9M+y+bUy5+PUeI/GC+VvFWurWvQ2mfK3k/8Fw+rsNTJ4TGHgV4N1OsF7A97uDyZ4dLYpnX3R/36cLNaMpfG00oBmPh7uDSuxgPYLJPDLKRjCVcAwJNWa7LwACyuMgA8czCpQkZGy+fS2KEkh9UamL22hcK4NJ4a+Xu4vFgDC//+Yjy1NWkMoTx5yWfRm9es4Ld6GNY1GParPMql8dQ8RydGaU6n5hBQGo/Ha7s8nqEX9BzAiWEflybLfebdB+o+a+0bq9K27bNt0zWhAIBlyOXAcfj4LPo+Li/cfc3KACpLJFkyo0FPTVguN/ukfp+F3+uaWg+f9B2t/HuDFB+2z+l+faA7lfBkblLyebjCwp97eenz4vhHXiVcbUdzGK48fIzjOX8wxXQeVes2hFB5wouw2gmDXC9Dm848jFHJ72eesClkNKoS688c2PbNyWLNVC9sRaVtXbNt0zWjAE6OBpXF7XA1k4Vucn2HfVxeWEpWy+fyeIbL46l5PADw5MXDxs+b508KQ3GVS8E3hnS8a3C5tlYNh09SwuNpHdKS+QfuNZ7MqsPEkuc56Zyj9B3Na+bk93oMwOLAnUxdoc1LE59ce9bg5GLfVEaBTSEBq3tV1pCOeY2LfWn5Dtuka0YB7A97K1mfByskFE8YijtODKv48PlDW7hiaaFXlsy+UjdQWhrmBOQFqzXZc/EDlfA/dWmCGPXwSeJfJhRtc+pJ3Nf8Di8PqNbAJBOjYk6Nh4PZa1vMSS0TCv+w38OwH8zjAZZhuMvmNei5Evc5eMI6p5cXCsY0nmwN+ouKfW08iT//PPt8577pPIArRCdGfVwcT3E4tcHTcpSOOQRUJ0QtIaDqmU9dmjgPn8MqnKKEjNJhYE1OLVFDNv7RovCq5jce6NbnA8vDwR5/7uOSy1rtmREu6flANUdWjyHx5z+r/G5r9dDED1RzmsZjtdDPH1TQYItcL2P09jDcwcSOGtof9jGPwNOXjfsmW4OThpBR+o5Wud4v19jK3+UAdksnhwNzZj/xpPizDf2w5Le61kBlCVgEP3c1fSEjH6rHan2mkI7Hksn5PeEBc/x5tIhXG3MGJ0eDGvLnOdCfuDh2Ke0nLo5NeR53OGFYrrFtTj1rtp+tmUWuq1Dr1ByGc4eAsjmyygRQrYFVISX+/GdtPLVcK6HK/qLGaBli6kJAO6H9XPCtB3pyTY2xycks4vzBxBzLBCrrzRUCumQ8fApB1gpsWtakcY6s8eeKf1DzW70wa4FN4nElFJ1rkOdVfPFn3xpbD59Bv4dRv2dGGSUeT/z5pHON85COtw7AKhNA2ge2sCBQrYHVkEv8pnqeIqa/bwxLec6ibdI1owBODHu1B2BN0l48rFr5eqw9b0jncxcn5o2S+D2WzOcujrE/1ENGw34Pg16oMeXWwyTx2zZvb+XnWze7FzUE2Od0P+M3jSdfY5dMeOS051qD/cacetfAFlZ75mCKebQnmQHgKaths/Kc+vlN4x8s+QHPnHY5gJ3SydEAF43wvYq/7+JPwnLhcOo6fC4c2vDSKYRzwZg0TmO4aIyHp89cMLaOSDw1v/E7u/hHfv7ZPJrzPMkiNM/pGmvs4b84tkGDgUouPHN00rnG+xm/VU49zz859MnpOmvg8TCsMpFqjFz7xnm2bJOuGQVww4nlHbfPOinfd7sKf86Tf5Z9fsb/LAu/czynRoP6AHnWyZHKn79jNOhh32Dher9DPo5nndDHlH/nG5z87jndxvidMtHvBVy3uCDEMp7yuSa588r1lvdNQ4acz3fLhOH5+a1qFv78Haf3BhgYuns2voPxHduia0YB3Jht3hsNB6KX/6aM56ZTm+fPhdEynl4v4MbFZ240Clkax00nRypaIvEBQAi2w6cxp6d0/nxeTHPq5L/x1Ij8N88/JP/N0XV7g/ruact48nHY+atxXLc/UOPV+fPLf3PkntOTzjl18jfH45Qhw74Z9Hu1LFv2Wc5nkYl8HP1eUK/x3DZdMwogF5ar4zDx8e8P+zi1cHnNh8NJ52FSC7LvsLrhxNBk+XjX4MaTvs3uPkycSrgxHsPhEEKox2E+TJz8NzkVRhr3qN+r5Ukcj9cQytfYxO+b08a+MfCfHPXrHk9WuXbPaWY4WShfY4uhtU1aSwGEEH42hPC3IYT3hBDeEEJ4FsN3dwjhQyGEB0MIr1nnnavSTaf2sn/7BM2tMAyCkGt+q+C4LZMVBfnZZv7VxgPYYp83nfYePksey3e40amQcgvbepg8u14Dm3X4bOeauZX8gu/ESMfEl8+1hES8+yznscxpHpJ69mmdP4SwguE0NI8nf65XJqz7bJu0rgfwJgBfHmN8IYAPA/iJkiGE0AfwWgDfAuAFAF4ZQnjBmu91U77ZLcmgtUI0hs2ebz6r4PS84QSnILs9gAX/wJCsBJpz6gkxAcYQUyNEo3+H/EJu65x6+ZOs+dfAF7az5gwSv9XwzMdtCTHd5FwDr6ede5pWwyPJmpU/5XqshtmqHsCVbgUNrKkAYox/HGOcLn58K4DnEGwvAvBgjPGhGOMYwOsAvHyd965C3g2eC+/1+/rmyq90tArakt+2eRNqxZqcOrU44OzhhGU82ca/EGQjkuFG5+GTr5klxNQMVziVsHFOk9KwJtZTMt0eT14t/qxdKVrynzL2oLEeauXzq3fYUTeAXa6pd0l0am+hhI3PT2OyjieN43qrEk5rdoWvgwTWvBO4oH8K4HeJ398B4OHs50cAvHiD7zXRF99yGt/74jvxZbdfb+K/fn+AH/yGL8bJYV/F0Cf6yZc+H59++hC3Xb9v4v/Z73ohHvjE5/Dld9xg4v/xu5+PP3rvY/j6L77ZxP+9L/4CzOcR3/yCW038d3/5bfjI4xfw3V9D6fE2/f0vuRmv+Ko78I++8vNN/C98zg34r848B1/7Rc828d/xrBN41dffhefedNLEP+z38M9e8jxMZ7a2BQDwMy//O/jIZy7grmefMvH/H6/4CvzFh8/hq+8ko50tevU/fB5uv+ER/JdfeouJ/9u/6g6cu3CIb33h7Sb+l3zZrXjPI0/hlS+608T/4i+8Cd/51c/BNxll4ktvO41XvuhOvPA5Nhm96dQI9/yDL8INJ4YmLy+EgB+/+/l48uIhbjm9p/IDwL/8zq/Aux5+Cs+//ToT///yrS/Af3z/p/Fio9x9/9ffhWE/mOfoW194Oz7+xEW84qvvMPH/F196M77jq+4w82+TQoxRZgjhTwDcRvzpp2KMf7jg+SkAZwC8IhYPDCF8F4C7Y4w/sPj5HwN4cYzx1cz77gFwDwDceeedX/OJT3zC94066qijjq5hCiE8EGM8Y+FVzaQY40uUl70KwLcB+Mby8F/QowCem/38nMXvuPfdC+BeADhz5oysnTrqqKOOOlqZ1kUB3Q3gxwC8LMZ4iWF7O4DnhRC+MIQwAvA9AO5b570dddRRRx2tT+tmIX4BwHUA3hRCeFcI4ZcAIITw+SGE+wFgkSR+NYA3AvgggNfHGN+/5ns76qijjjpak9ZKAscYv4T5/acAvDT7+X4A96/zro466qijjjZLVx6H1FFHHXXU0RWhTgF01FFHHV2j1CmAjjrqqKNrlDoF0FFHHXV0jZJaCHYlKYRwDsCqlWA3A/jsBoezazrq4weO/nc46uMHjv53OOrjB3b/Hb4gxmgqPb+qFcA6FEI4a62GuxrpqI8fOPrf4aiPHzj63+Gojx+4ur9DFwLqqKOOOrpGqVMAHXXUUUfXKB1nBXDvlR7AmnTUxw8c/e9w1McPHP3vcNTHD1zF3+HY5gA66qijjjqS6Th7AB111FFHHQl07BTA1XD/8DoUQvi1EMLjIYT3XemxrEIhhOeGEP4shPCBEML7Qwg/cqXH5KUQwn4I4W0hhHcvvsP/dqXHtAqFEPohhHeGEP6/Kz2WVSiE8PEQwnsXjSbPXunxeCmE8KwQwu8t7k3/YAjh6670mEo6ViGgxf3DHwbwTahuHns7gFfGGD9wRQfmoBDCPwBwAcBvxhi//EqPx0shhNsB3B5jfEcI4ToADwD49iO2BgHAqRjjhRDCEMBbAPxIjPGtV3hoLgoh/Ciqi5qujzF+25Uej5dCCB8HcCbGeCTrAEIIvwHgL2OMv7JohX8yxvjUlR5XTsfNA7gq7h9eh2KMbwbw/7d3/zwyxWEUx7+HVQgJBRHJFqvS2i00KyIkQohXQKEWiZbGOxCdZpdILAprW7EJtWJFoqDSWCErUfjTCI7i/iQSirk74Td37vkkk7kz1Zli7nPvM0/m+VA7x3rZfmv7aTn+RPMX4PV337XgxufyclN5dOpKSdIkcAKYq52ljyRtAw4C8wC2v47ayR/GrwD8bf9wp04+40TSFDANPKmbpL3SPnkGrAHLtrv2Ga7SLGv6UTvIEAw8lLRSVsV2yR7gPXCjtOHmJA22ePo/GrcCECNC0lZgEbhg+2PtPG3Z/m57H80K0/2SOtOOk3QSWLO9UjvLkA7YngGOA+dKe7QrJoAZ4JrtaeALMHK/SY5bAWi1fzj+jdI3XwQWbN+vnWcY5bb9MXCsdpYWZoFTpYd+Fzgs6VbdSO3ZflOe14AlmhZvV6wCq7/dOd6jKQgjZdwKQPYPV1Z+QJ0HXti+UjvPekjaKWl7Od5MM1Twsm6qwdm+aHvS9hTNd+CR7dOVY7UiaUsZIqC0To4CnZmMs/0OeC1pb3nrCDBygxBDrYQcNba/Sfq1f3gjcL1r+4cl3QEOATskrQKXbc/XTdXKLHAGeF566ACXylrQrtgN3CxTZRto9lh3cpSyw3YBS831BBPAbdsP6kZq7TywUC5GXwFnK+f5w1iNgUZExODGrQUUEREDSgGIiOipFICIiJ5KAYiI6KkUgIiInkoBiIjoqRSAiIieSgGIiOipn6iTmm6gP7haAAAAAElFTkSuQmCC\n"
},
"metadata": {}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "N = len(t) # データ数\nyf = np.fft.fft(y) / (N / 2) # N等分して、正の部分だけ扱うので2倍する\nfreq = np.fft.fftfreq(N, d=dt) # 周波数\nyf_abs = np.abs(yf)\n\nplt.plot(freq[:int(N/2.)], yf_abs[:int(N/2.)])",
"execution_count": 122,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 122,
"data": {
"text/plain": "[<matplotlib.lines.Line2D at 0x13b634cf8>]"
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/plain": "<matplotlib.figure.Figure at 0x13b4d0da0>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAD8CAYAAABw1c+bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAG21JREFUeJzt3X2QXNV95vHvM2+SRhJohIYX6x2jNcYBRDyR7QWvwRtkmTjI2ThrESfBKXtVyUJi70u2wNmCLC7XehNXYrvMLlbZWuytIEJsE2u9skExJDghgEa2zKsxQsZoJmDJegW9jWb6t3/0ben2aEZ9e6ZHPffO86ma6r7n3r5zLgxPH8459x5FBGZmNnW0NLsCZmZ2Zjn4zcymGAe/mdkU4+A3M5tiHPxmZlOMg9/MbIpx8JuZTTE1g1/SQkkPS3pW0jOSPjbCMZL0eUnbJT0p6RdT+26U9ELyc2OjL8DMzOqjWjdwSboAuCAivi9pNrAVeH9EPJs65jrgD4DrgLcBn4uIt0maC/QCPUAkn31rROybkKsxM7Oa2modEBGvAK8k71+T9BwwH3g2ddhq4KtR/hZ5TNKc5AvjamBzROwFkLQZWAVsON3vnDdvXixZsqT+qzEzm6K2bt3684joznJszeBPk7QEuAJ4fNiu+cDO1HZfUjZa+WktWbKE3t7eeqpmZjalSfpp1mMzD+5KmgV8Hfh4RBwcS8VqnH+tpF5Jvbt372706c3MLJEp+CW1Uw79v4yIb4xwSD+wMLW9ICkbrfwUEbEuInoioqe7O9P/rZiZ2RhkmdUj4MvAcxHx56McthH4nWR2z9uBA8nYwAPASkldkrqAlUmZmZk1SZY+/iuB3waekrQtKfsEsAggIu4CNlGe0bMdOAz8brJvr6RPAluSz91RGeg1M7PmyDKr5x8A1TgmgJtG2bceWD+m2pmZWcP5zl0zsynGwW9mNsU4+OuwfddrPL5jT7OrYWY2Lg7+Otz58It84v6nml0NM7NxcfDXYWCoxGDJi9ObWb45+OsRUOOZdmZmk56Dvw6lCEpOfjPLOQd/HcItfjMrAAd/HYKg1voFZmaTnYO/DqUoryZjZpZnDv46ROA+fjPLPQd/HSICz+Y0s7xz8Nch8OCumeWfg78OpfDgrpnln4O/DuHBXTMrAAd/HXwDl5kVgYO/Ts59M8u7mitwSVoPvA/YFRG/MML+PwI+lDrfm4HuZNnFl4DXgCFgMCJ6GlXxZnCL38yKIEuL/25g1Wg7I+LPImJ5RCwHbgX+fti6utck+3Md+pC09p37ZpZzNYM/Ih4Bsi6QfgOwYVw1msR8A5eZFUHD+vgldVL+P4Ovp4oDeFDSVklrG/W7mqUU4Qa/meVezT7+Ovwq8I/Dunmuioh+SecCmyX9KPk/iFMkXwxrARYtWtTAajVO4Ba/meVfI2f1rGFYN09E9Cevu4D7gRWjfTgi1kVET0T0dHd3N7BajeNHNphZETQk+CWdDbwL+GaqbKak2ZX3wErg6Ub8vmbx4K6ZFUGW6ZwbgKuBeZL6gNuBdoCIuCs57NeAByPiUOqj5wH3S6r8nnsi4juNq/qZ5+mcZlYENYM/Im7IcMzdlKd9pst2AJePtWKTkRv8ZlYEvnO3DiVP5zSzAnDw1yPCj2wws9xz8NehMqPHj2Y2szxz8NehcvuWc9/M8szBX4dSKXl18ptZjjn461CJe9/EZWZ55uCvQ6Vv30/sMbM8c/DXIaL61cwsjxz8dfDgrpkVgYO/DpW+fQ/umlmeOfjrcLKP38wsvxz8dQi3+M2sABz8dajEvXPfzPLMwV+HSkvfj2wwszxz8NfB0znNrAgc/HWotPjdx29meebgr8PJwd3m1sPMbDxqBr+k9ZJ2SRpxvVxJV0s6IGlb8nNbat8qSc9L2i7plkZWvJn8yAYzy7MsLf67gVU1jvleRCxPfu4AkNQK3Am8F7gEuEHSJeOpbLOdHNxtckXMzMahZvBHxCPA3jGcewWwPSJ2RMQAcC+wegznmTQ8uGtmRdCoPv53SPqhpG9LektSNh/YmTqmLynLLQ/umlkRtDXgHN8HFkfE65KuA/4GWFbvSSStBdYCLFq0qAHVarwY9mpmlkfjbvFHxMGIeD15vwlolzQP6AcWpg5dkJSNdp51EdETET3d3d3jrdaEqNy4VfK0HjPLsXEHv6TzJSl5vyI55x5gC7BM0lJJHcAaYON4f18zuYfHzIqgZlePpA3A1cA8SX3A7UA7QETcBXwA+H1Jg8ARYE2Um8aDkm4GHgBagfUR8cyEXMUZ4j5+MyuCmsEfETfU2P8F4Auj7NsEbBpb1SYfr7lrZkXgO3frUOnb90PazCzPHPx1cIvfzIrAwV+Hkw19J7+Z5ZeDvw4npnM6980sxxz8dfAKXGZWBA7+Ong6p5kVgYO/Dn5Im5kVgYO/DicXYnHym1l+Ofjr4AVYzKwIHPx1KLnFb2YF4OCvg6dzmlkROPjrUDoxuOvkN7P8cvBnlA57t/jNLM8c/BlVN/Kd/GaWXw7+jNJR7xa/meWZgz+jdFePu/jNLM8c/BmlW/mezmlmeVYz+CWtl7RL0tOj7P+QpCclPSXpUUmXp/a9lJRvk9TbyIqfaembt5z7ZpZnWVr8dwOrTrP/J8C7IuJS4JPAumH7r4mI5RHRM7YqTg7psPd0TjPLsyxr7j4iaclp9j+a2nwMWDD+ak0+UdXV07x6mJmNV6P7+D8CfDu1HcCDkrZKWnu6D0paK6lXUu/u3bsbXK3xS/fr+5k9ZpZnNVv8WUm6hnLwX5Uqvioi+iWdC2yW9KOIeGSkz0fEOpJuop6enkmXrJ7OaWZF0ZAWv6TLgC8BqyNiT6U8IvqT113A/cCKRvy+Zqhq8buP38xybNzBL2kR8A3gtyPix6nymZJmV94DK4ERZwblQfXgbvPqYWY2XjW7eiRtAK4G5knqA24H2gEi4i7gNuAc4H9KAhhMZvCcB9yflLUB90TEdybgGs6IcB+/mRVEllk9N9TY/1HgoyOU7wAuP/UT+VQ1q6fUvHqYmY2X79zNqHpWj5lZfjn4M6qe1ePoN7P8cvBn5MFdMysKB39G4emcZlYQDv6MfAOXmRWFgz8jP7LBzIrCwZ+RH9JmZkXh4M/Ij2wws6Jw8GfkWT1mVhQO/oyqgt99/GaWYw7+jNJh70c2mFmeOfgzKlW1+M3M8svBn1F6QNePbDCzPHPwZ5SOes/qMbM8c/BnVP3IhiZWxMxsnBz8GfkGLjMrikzBL2m9pF2SRlw6UWWfl7Rd0pOSfjG170ZJLyQ/Nzaq4mdaydM5zawgsrb47wZWnWb/e4Flyc9a4H8BSJpLeanGt1FeaP12SV1jrWwzVU3ndO6bWY5lCv6IeATYe5pDVgNfjbLHgDmSLgDeA2yOiL0RsQ/YzOm/QCatqrn77uQ3sxxrVB//fGBnarsvKRutPHfc4jezopg0g7uS1krqldS7e/fuZlfnFNXP6nHym1l+NSr4+4GFqe0FSdlo5aeIiHUR0RMRPd3d3Q2qVuN4Vo+ZFUWjgn8j8DvJ7J63Awci4hXgAWClpK5kUHdlUpY71QuxmJnlV1uWgyRtAK4G5knqozxTpx0gIu4CNgHXAduBw8DvJvv2SvoksCU51R0RcbpB4knLd+6aWVFkCv6IuKHG/gBuGmXfemB9/VWbXEp+Vo+ZFcSkGdyd7LwQi5kVhYM/M0/nNLNicPBn5Ec2mFlROPgzclePmRWFgz+jqumcTn4zyzEHf0a+gcvMisLBn5EXYjGzonDwZ5TOes/jN7M8c/Bn5D5+MysKB39GVbN6mlcNM7Nxc/Bn5Ec2mFlROPgzqn5IW9OqYWY2bg7+rDyd08wKwsGfUfXz+J38ZpZfDv6M/MgGMysKB39Gns5pZkXh4M+o+gauplXDzGzcMgW/pFWSnpe0XdItI+z/C0nbkp8fS9qf2jeU2rexkZU/k/zIBjMrippLL0pqBe4ErgX6gC2SNkbEs5VjIuI/pI7/A+CK1CmORMTyxlW5Oaof0ubkN7P8ytLiXwFsj4gdETEA3AusPs3xNwAbGlG5yaRqIRYHv5nlWJbgnw/sTG33JWWnkLQYWAo8lCqeLqlX0mOS3j/mmjZZegqnY9/M8qxmV0+d1gBfi4ihVNniiOiXdCHwkKSnIuLF4R+UtBZYC7Bo0aIGV2v8Su7qMbOCyNLi7wcWprYXJGUjWcOwbp6I6E9edwB/R3X/f/q4dRHRExE93d3dGap1Znlw18yKIkvwbwGWSVoqqYNyuJ8yO0fSxUAX8E+psi5J05L384ArgWeHfzYPvAKXmRVFza6eiBiUdDPwANAKrI+IZyTdAfRGROVLYA1wb1SPfL4Z+KKkEuUvmU+nZwPlSfVjGpz8ZpZfmfr4I2ITsGlY2W3Dtv9khM89Clw6jvpNGlUt/lLz6mFmNl6+czejqumcbvGbWY45+DOKqoVYmlgRM7NxcvBnVMl9ydM5zSzfHPwZVbp3WiWP7ZpZrjn4M6p077S0yC1+M8s1B39GlaxvldzgN7Ncc/BnVGnlt7bIg7tmlmsO/owqWd8iP53TzPLNwZ9RpFr8zn0zyzMHf0Yn+vhb5Bu4zCzXHPwZVVr8LZIf2WBmuebgz6jkFr+ZFYSDP6OTg7ue1WNm+ebgz6h6cNfJb2b55eDPqGpw17lvZjnm4M+odGJw1w9pM7N8c/BnVIn68uCumVl+ZQp+SaskPS9pu6RbRtj/YUm7JW1Lfj6a2nejpBeSnxsbWfkz6eQjG1o8uGtmuVZz6UVJrcCdwLVAH7BF0sYR1s79q4i4edhn5wK3Az2UG81bk8/ua0jtz6CTffx+ZIOZ5VuWFv8KYHtE7IiIAeBeYHXG878H2BwRe5Ow3wysGltVm+vErB55cNfM8i1L8M8Hdqa2+5Ky4X5d0pOSviZpYZ2fnfQqYd/iG7jMLOcaNbj7f4ElEXEZ5Vb9V+o9gaS1knol9e7evbtB1WqcEwux+JENZpZzWYK/H1iY2l6QlJ0QEXsi4liy+SXgrVk/mzrHuojoiYie7u7uLHU/o9JLL3o6p5nlWZbg3wIsk7RUUgewBtiYPkDSBanN64HnkvcPACsldUnqAlYmZblzsqvHS+6aWb7VnNUTEYOSbqYc2K3A+oh4RtIdQG9EbAT+UNL1wCCwF/hw8tm9kj5J+csD4I6I2DsB1zHhIgIJhIhwX4+Z5VfN4AeIiE3ApmFlt6Xe3wrcOspn1wPrx1HHSSEAkbT4h5pdGzOzsfOduxmVImiREO7jN7N8c/BnFEG5q0fu4zezfHPwZ1QKkIT8PH4zyzkHf0ZBlPv4Bb5118zyzMGfUQRJHz9u8ZtZrjn4M6pM52yRH9lgZvnm4M+oVGnx+5ENZpZzDv6MIsrz+OUVuMws5xz8GQWVrh6P7ZpZvjn4M4pkOue0tlYGhtzXY2b55eDPKCJoEcyc1sqhY4PNro6Z2Zg5+DOq3MDV2dHm4DezXHPwZxQkLf6OVg4fH6LkyfxmllMO/ozKOS86p7URAUcH/YhOM8snB39G5Tt3Yea08pOsDx1z8JtZPjn4M6rcuTuzoxXA/fxmllsO/ozKN3CVB3cBDg04+M0snzIFv6RVkp6XtF3SLSPs/4+SnpX0pKTvSlqc2jckaVvys3H4Z/OilJrOCXB4wF09ZpZPNZdelNQK3AlcC/QBWyRtjIhnU4f9AOiJiMOSfh/4U+CDyb4jEbG8wfU+4wZLQVtrS6qP3y1+M8unLC3+FcD2iNgREQPAvcDq9AER8XBEHE42HwMWNLaazbfv8ABdne3M7PDgrpnlW5bgnw/sTG33JWWj+Qjw7dT2dEm9kh6T9P7RPiRpbXJc7+7duzNU68zad2iAOZ0ddFYGd93Hb2Y5VbOrpx6SfgvoAd6VKl4cEf2SLgQekvRURLw4/LMRsQ5YB9DT0zPp7o7ae3iAC7tnnejqOeyuHjPLqSwt/n5gYWp7QVJWRdIvA38MXB8RxyrlEdGfvO4A/g64Yhz1bZp9h47T1dlxYnD3kAd3zSynsgT/FmCZpKWSOoA1QNXsHElXAF+kHPq7UuVdkqYl7+cBVwLpQeFcGBgs8fqxQebObKejtYW2FnHYXT1mllM1u3oiYlDSzcADQCuwPiKekXQH0BsRG4E/A2YBfy0J4OWIuB54M/BFSSXKXzKfHjYbKBf2Hx4AYE5nR/KgtlYP7ppZbmXq44+ITcCmYWW3pd7/8iifexS4dDwVnAz2JsE/d2YHUH5sg6dzmlle+c7dDPYeKgd/V2c5+Ds7Wn0Dl5nlloM/g32HjgMnW/yzprV5OqeZ5ZaDP4NKV09XZzuAF2Mxs1xz8Gew/9DJwV2As2e0sycpMzPLGwd/BnsODTBrWhsdbeV/XBd2z+TlPYc57kXXzSyHHPwZbN/1OovP6TyxfdG5sxgsBT/dc6iJtTIzGxsHfw2lUvDDvv0sXzjnRNmyc2cD5S8EM7O8cfDX8JM9h3jt6CCXLzgZ/G88dybg4DezfHLw1/DDnfsBuDzV4u/saGP+nBm84OA3sxxy8Nfw6It7mDWtjYvOnVVVfvH5s/n+y/uImHQPEjUzOy0H/2kcOHKcbz35z/zq5RfQ2qKqfSvfch479x7h6f6DTaqdmdnYOPhP474tOzl6vMRvrlh8yr73vOV82lrEt5785ybUzMxs7Bz8o9i59zCf/dsf885l87h0wdmn7J/T2cG7Lz6Xe554mT2vHxvhDGZmk5ODfwSvHjjKjf/7CSTx6V+/bNTj/ug9b+LwwBCf+n/Pua/fzHLDwZ9ybHCIDU+8zHWf/x6vHjjK+g//EvPnzBj1+GXnzeamay7iGz/o5zMPPk+p5PA3s8mvoWvu5tGRgSG27dzPw8/v4pvb+vnZwWMsXziHz/zG5afM5BnJx//1MnYdPMqdD7/Ik30H+ONfeTMXn3/WGai5mdnYZAp+SauAz1FegetLEfHpYfunAV8F3grsAT4YES8l+24FPgIMAX8YEQ80rPbDvH5skOODJY6XShwfivL7oRIHjw5y8Mhx9h8ZYN+h4/TtO0LfvsO8tOcQ23e9TimgvVX8yzfO4zO/sZSrLppHspJYTS0t4r//m0t5yxvO4k+/8zyrPvs9fmlJF9cvn887L5rH4nM6M5/LzOxMUK2+aUmtwI+Ba4E+ymvw3pBeQlHSvwcui4jfk7QG+LWI+KCkS4ANwArgDcDfAv8iIk67iklPT0/09vbWfTFv+q/f5thg7QenzWhvZUHXDBbN7eSSN5zFZQvm8I43nsOsaeP7H6B9hwb46607uefxl3lpz2Gg/CjnN50/m4vPP4uFczu54OzpnH/2dM4/azpdnR1Mb2/xF4OZjZukrRHRk+XYLEm3AtgeETuSk98LrKZ60fTVwJ8k778GfEHlNFsN3BsRx4CfSNqenO+fslSuXre+92ICaG9tob1VtLe20NbawlnT2zh7RjtzOjuYM6OdOZ3tExK2XTM7WPuv3si/e+eF/OTnh/jHF/fwTP8BfvTqa9zXu3PEVbs6Wls4a0YbZ01v56wZ7cye3sb09tbyT1sL09pbmN6WbLe3MK2tlWntLbS2iLYW0drSkryWt9taq7dbW0Rb68njWiQkaJFoEUggJeVwYv/JY8rHkdo+cVwLJ96P9LnKP+HKP2p/wZlNDlmCfz6wM7XdB7xttGOSxdkPAOck5Y8N++z8Mde2hg9fuXSiTl0XSVzYPYsLu0+OEUQE+w8f55UDR3n14BFePXCM/UcGOHhkkINHj3PwyPETXVK7XzvGwGCJo8eHODpY4ljyOlSwweMTXwgntjVsu7K/+kCNsK/WuRhtf9U5R/79jPo7Rq5DLfV8/U3El2Xmetbxq5Xxquo7Z8bj6jhp5iOb9O9ybmcH9/3eO+o469hMmsFdSWuBtQCLFi1qcm0aTxJdMzvomtnBJW8Y2+Dv8aESx5IvgqEIhkrB4FAwWAqGSiUGh28PJceUyq/Hh0qUovwlFEAp4uR2lLdPeU0dR/I6/LigujwiqCxVUD4DVHoUT3x1JQVRvTnq8VX7T9mX4TOpbYZ/ro7fz/D9wz5XSz2zfrMeWt85Mx48IfXMftKJufas58z47zL7r8588OzpZyaSs/yWfmBhantBUjbSMX2S2oCzKQ/yZvksABGxDlgH5T7+LJWfaspdWC3jHosws6ktyzz+LcAySUsldQBrgI3DjtkI3Ji8/wDwUJS/NjcCayRNk7QUWAY80Ziqm5nZWNRsOiZ99jcDD1Cezrk+Ip6RdAfQGxEbgS8D/ycZvN1L+cuB5Lj7KA8EDwI31ZrRY2ZmE6vmdM5mGOt0TjOzqaqe6Zx+ZIOZ2RTj4Dczm2Ic/GZmU4yD38xsinHwm5lNMZNyVo+k3cBPx/jxecDPG1idPPA1Tw2+5qljLNe9OCK6sxw4KYN/PCT1Zp3SVBS+5qnB1zx1TPR1u6vHzGyKcfCbmU0xRQz+dc2uQBP4mqcGX/PUMaHXXbg+fjMzO70itvjNzOw0ChP8klZJel7Sdkm3NLs+E0XSekm7JD2dKpsrabOkF5LXrmbWsdEkLZT0sKRnJT0j6WNJeWGvW9J0SU9I+mFyzf8tKV8q6fHk7/yvkkelF4qkVkk/kPStZLvQ1yzpJUlPSdomqTcpm9C/7UIEf7Ig/J3Ae4FLgBuShd6L6G5g1bCyW4DvRsQy4LvJdpEMAv8pIi4B3g7clPz7LfJ1HwPeHRGXA8uBVZLeDvwP4C8i4iJgH/CRJtZxonwMeC61PRWu+ZqIWJ6awjmhf9uFCH5SC8JHxABQWRC+cCLiEcprHqStBr6SvP8K8P4zWqkJFhGvRMT3k/evUQ6F+RT4uqPs9WSzPfkJ4N3A15LyQl0zgKQFwK8AX0q2RcGveRQT+rddlOAfaUH4CVvUfRI6LyJeSd6/CpzXzMpMJElLgCuAxyn4dSddHtuAXcBm4EVgf0QMJocU8e/8s8B/AZJVmzmH4l9zAA9K2pqsPQ4T/LftxVsLJiJCUiGnakmaBXwd+HhEHCw3BsuKeN3JanXLJc0B7gcubnKVJpSk9wG7ImKrpKubXZ8z6KqI6Jd0LrBZ0o/SOyfib7soLf7Mi7oX1M8kXQCQvO5qcn0aTlI75dD/y4j4RlJc+OsGiIj9wMPAO4A5kioNtqL9nV8JXC/pJcrdte8GPkexr5mI6E9ed1H+gl/BBP9tFyX4sywIX2Tpxe5vBL7ZxLo0XNLP+2XguYj489Suwl63pO6kpY+kGcC1lMc2HgY+kBxWqGuOiFsjYkFELKH83/BDEfEhCnzNkmZKml15D6wEnmaC/7YLcwOXpOso9w9WFoT/VJOrNCEkbQCupvz0vp8BtwN/A9wHLKL8VNN/GxHDB4BzS9JVwPeApzjZ9/sJyv38hbxuSZdRHtRrpdxAuy8i7pB0IeXW8FzgB8BvRcSx5tV0YiRdPf85It5X5GtOru3+ZLMNuCciPiXpHCbwb7swwW9mZtkUpavHzMwycvCbmU0xDn4zsynGwW9mNsU4+M3MphgHv5nZFOPgNzObYhz8ZmZTzP8HCkPgBpAxD0MAAAAASUVORK5CYII=\n"
},
"metadata": {}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "plt.plot(freq[:50], yf_abs[:50])",
"execution_count": 123,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 123,
"data": {
"text/plain": "[<matplotlib.lines.Line2D at 0x13b6fc4a8>]"
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/plain": "<matplotlib.figure.Figure at 0x13b1ffb00>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAD8CAYAAABw1c+bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzt3X2UXHWd5/H3t566Ow8kgXRiyCMKahjlQdqgC6u4IBNZBd3xaNBR9OhmddXRmXN2js6cI7O4Z8Zd98yoRx3IYsaHVVBR3DgniqgojhpNB1EEDIQAJhHoDnlO6rm++8e9t7u60911q7o6VTf9eZ1Tp+o+1O1foOpTv/reX92fuTsiIjJ7pDrdABERObUU/CIis4yCX0RkllHwi4jMMgp+EZFZRsEvIjLLKPhFRGaZhsFvZivN7B4ze8jMHjSzD06wj5nZp81sl5n91sxeUrftBjN7NLzd0O5/gIiINMca/YDLzJYBy9z9PjObD+wAXu/uD9Xtcw3wAeAa4FLgU+5+qZmdCQwCA4CHz73E3Q/OyL9GREQayjTawd2fAp4KHx81s4eB5cBDdbtdB3zJg0+RbWa2MPzAuAK4290PAJjZ3cB64Lap/ubixYt9zZo1zf9rRERmqR07dux39/44+zYM/npmtga4GPjluE3LgT11y3vDdZOtn9KaNWsYHBxspmkiIrOamT0Zd9/YJ3fNbB7wTeBD7n6klYY1OP5GMxs0s8Hh4eF2H15EREKxgt/MsgSh/xV3/9YEu+wDVtYtrwjXTbb+JO6+yd0H3H2gvz/WtxUREWlBnFE9BnweeNjd/3GS3bYAbw9H97wMOByeG7gLuNrMFpnZIuDqcJ2IiHRInBr/ZcDbgAfM7P5w3d8AqwDc/WZgK8GInl3ACeCd4bYDZvYxYHv4vJuiE70iItIZcUb1/BtgDfZx4H2TbNsMbG6pdSIi0nb65a6IyCyj4BcRmWUU/CIt2Pn0UbY/odNVkkwKfpEWfOqHj/C3dz7Q6WaItETBL9KCY8Uqx4vVTjdDpCUKfpEWFMpV8mUFvySTgl+kBcVylXxJwS/JpOAXaUE+7PE3uqy5SDdS8Iu0oFCuAVCs1DrcEpHmKfhFWlAI6/sF1fklgRT8Ii2ITuzqBK8kkYJfpAXFsNSjE7ySRAp+kSZVa06pGga/evySQAp+kSbV1/VV45ckUvCLNKk+7PMljeqR5FHwizSpUDeEU6UeSSIFv0iTVOqRpGs4A5eZbQZeCwy5+4sm2P7fgLfWHW8t0B9Ou/gEcBSoAhV3H2hXw0U6pX4kj3r8kkRxevxfANZPttHdP+HuF7n7RcBHgJ+Mm1f3VeF2hb6cFooV9fgl2RoGv7vfC8SdceJ64LZptUiky0WXawCN45dkaluN38zmEHwz+Gbdage+b2Y7zGxju/6WSCep1CNJ17DG34TXAT8bV+a53N33mdkS4G4z+334DeIk4QfDRoBVq1a1sVki7VWoKPgl2do5qmcD48o87r4vvB8C7gTWTfZkd9/k7gPuPtDf39/GZom0V32pp6BSjyRQW4LfzBYArwT+X926uWY2P3oMXA38rh1/T6STol5+OmVjPgREkiLOcM7bgCuAxWa2F7gRyAK4+83hbm8Avu/ux+ueuhS408yiv/NVd/9e+5ou0hnFMPgX9mVV6pFEahj87n59jH2+QDDss37dbuDCVhsm0q2iIZyL5uYU/JJI+uWuSJPy5SrplDG/N6Nx/JJICn6RJhXKNXozKfqyaY3jl0RS8Is0qVCu0pdLB8GvHr8kkIJfpEmFco2eTJrebFqlHkkkBb9IkwrlKr3ZVBj8Gs4pyaPgF2lSEPxp+nIplXokkRT8Ik0qVKr0ZdM6uSuJpeAXaVK+FPb4w5O77t7pJok0RcEv0qRCuRbU+HNpAIoV1fklWRT8Ik0qVEZ7/KBr8kvyKPhFmlQISz29YfDXX6ZZJAkU/CJNKlSCUo96/JJUCn6RJhXKwaieqMevIZ2SNAp+kSa4O/mRcfxhqUfBLwmj4BdpQqlaw51xJ3c1qkeSRcEv0oToEg09mboav3r8kjAKfpEmRLNv9eWCSzaAgl+Sp2Hwm9lmMxsyswnnyzWzK8zssJndH94+WrdtvZntNLNdZvbhdjZcpBOikO/NpOnJqMYvyRSnx/8FYH2DfX7q7heFt5sAzCwNfBZ4DXA+cL2ZnT+dxop0WlTq0cldSbKGwe/u9wIHWjj2OmCXu+929xJwO3BdC8cR6RqFkVKPxvFLcrWrxv9yM/uNmX3XzP4kXLcc2FO3z95wnUhi1Zd6NI5fkirThmPcB6x292Nmdg3wbeC8Zg9iZhuBjQCrVq1qQ7NE2i/q8fdk06RTRi6ja/JL8ky7x+/uR9z9WPh4K5A1s8XAPmBl3a4rwnWTHWeTuw+4+0B/f/90myUyI6Iaf1Tm6cumKajUIwkz7eA3s+eYmYWP14XHfBbYDpxnZueYWQ7YAGyZ7t8T6aSox9+bDd46fZp+URKoYanHzG4DrgAWm9le4EYgC+DuNwNvBN5rZhUgD2zwYGaKipm9H7gLSAOb3f3BGflXiJwio8GfDu9V6pHkaRj87n59g+2fAT4zybatwNbWmibSfUZG9YwEf1rBL4mjX+6KNCFfN44fgl/wahy/JI2CX6QJI6N6MqM1fo3jl6RR8Is0oVCpksukSKUMYGTCdZEkUfCLNKFYro3U9wF6cwp+SR4Fv0gT8qXqyFBOCH7BW9RwTkkYBb9IEwqV6siJXQiu2aMevySNgl+kCdF8uxGd3JUkUvCLNCFfrtEzPvjLVYLfLIokg4JfpAmFcpXeTF2NP7wmf7GiOr8kh4JfpAnFcnVkAhZA1+SXRFLwizQhX67Sm5kg+HWCVxJEwS/ShEK5NnY4Z1bTL0ryKPhFmlAYV+rRLFySRAp+kSYUylV66ks9mnBdEkjBL9KEoNQz0cldjeqR5FDwi8RUrTml6tgav07uShIp+EViKlbGTsICwSUbQMEvydIw+M1ss5kNmdnvJtn+VjP7rZk9YGY/N7ML67Y9Ea6/38wG29lwkVMtGqtfX+rRqB5Jojg9/i8A66fY/jjwSnd/MfAxYNO47a9y94vcfaC1Jop0h0Ilmn1Lwzkl2eLMuXuvma2ZYvvP6xa3ASum3yyR7jN+onXQL3clmdpd438X8N26ZQe+b2Y7zGxjm/+WyCk1ValHNX5JkoY9/rjM7FUEwX953erL3X2fmS0B7jaz37v7vZM8fyOwEWDVqlXtapZI20Qnd+uDP50ychldk1+SpS09fjO7ALgVuM7dn43Wu/u+8H4IuBNYN9kx3H2Tuw+4+0B/f387miXSVoVwpq36UT3RckGlHkmQaQe/ma0CvgW8zd0fqVs/18zmR4+Bq4EJRwaJJMFoqWfs20YTrkvSNCz1mNltwBXAYjPbC9wIZAHc/Wbgo8BZwOfMDKASjuBZCtwZrssAX3X3783Av0HklChMUOqB4LINBc27KwkSZ1TP9Q22vxt49wTrdwMXnvwMkWSKwr3+sswAParxS8Lol7siMY0M58yNK/Xk0hrHL4mi4BeJaaJx/KAJ1yV5FPwiMY0Ef2aC4FePXxJEwS8SU6FcI50ysmkbs743p+CXZFHwi8QUzLebIhypNkLj+CVpFPwiMRXK1ZPq+xCM648u4CaSBAp+kZjGz74V0cldSRoFv0hMQY//5LdMdHLX3TvQKpHmKfhFYpq01BNOuF5UuUcSQsEvElOhUj3pAm2ga/JL8ij4RWLKlybu8WvCdUkaBb9ITMHJ3Qlq/DlNvyjJouAXialQqdIzQY+/J6MevySLgl8kpmK5NnGNXz1+SRgFv0hM+SmGcwLkSxrVI8mg4BeJqVCunnSBNtDJXUkeBb9IDO5OoVwdKevU6wuvz6/gl6SIFfxmttnMhsxswjlzLfBpM9tlZr81s5fUbbvBzB4Nbze0q+Eip1KpWqPmJ1+LH0bX6UJtkhRxe/xfANZPsf01wHnhbSPwzwBmdibBHL2XAuuAG81sUauNFemUaNrFnszkNf5oTl6Rbhcr+N39XuDAFLtcB3zJA9uAhWa2DPhT4G53P+DuB4G7mfoDRKQrFcMyzkSlnl79clcSpl01/uXAnrrlveG6ydaLJEp+ktm3oC74VeOXhOiak7tmttHMBs1scHh4uNPNERkjKvVMVONPp4xcJqXgl8RoV/DvA1bWLa8I1022/iTuvsndB9x9oL+/v03NEmmPwkipZ+K3jGbhkiRpV/BvAd4eju55GXDY3Z8C7gKuNrNF4Undq8N1IokyVakHNOG6JEsmzk5mdhtwBbDYzPYSjNTJArj7zcBW4BpgF3ACeGe47YCZfQzYHh7qJnef6iSxSFeKevwTXasHgpO++bJ+uSvJECv43f36BtsdeN8k2zYDm5tvmkj3GK3xT/wluTeb1rV6JDG65uSuSDcrhmP0J7pIG4QTriv4JSEU/CIxRGP0JxrVA5pwXZJFwS8SQ9SbnzL41eOXhFDwi8RQCCdSn7TUk1PwS3Io+EViiMo4E12rBzSOX5JFwS8SQ6FSJZdJkUrZhNv7sumRbwUi3U7BLxLDZNMuRvpyOrkryaHgF4khX5p42sVIb3itnuAnLSLdTcEvEkOhUp10RA8EJ3cBiir3SAIo+EViKJSrU5d6dE1+SRAFv0gMhXJt0uv0gCZcl2RR8IvEkC9X6Z1kKCeMzsyl4JckUPCLxFAsN6jxR/PuKvglART8IjEUGgznVPBLkij4RWLIl6cezjl6clejeqT7KfhFYig0KPXo5K4kiYJfJIaGwR/OxavglySIFfxmtt7MdprZLjP78ATb/8nM7g9vj5jZobpt1bptW9rZeJFTpVCuxTu5q3H8kgANp140szTwWeDVwF5gu5ltcfeHon3c/S/r9v8AcHHdIfLuflH7mixyalVrTqlai1fjV49fEiBOj38dsMvdd7t7CbgduG6K/a8HbmtH40S6QaNpF2F0HL9G9UgSxAn+5cCeuuW94bqTmNlq4BzgR3Wre81s0My2mdnrJ/sjZrYx3G9weHg4RrNETo1G0y4C9GbU45fkaPfJ3Q3AHe5e/+pf7e4DwFuAT5rZ8yZ6ortvcvcBdx/o7+9vc7NEWhddZ3+qUk8qZeTCK3SKdLs4wb8PWFm3vCJcN5ENjCvzuPu+8H438GPG1v9Ful6j+XYjmoVLkiJO8G8HzjOzc8wsRxDuJ43OMbMXAouAX9StW2RmPeHjxcBlwEPjnyvSzZoJfvX4JQkajupx94qZvR+4C0gDm939QTO7CRh09+hDYANwu4+diWItcIuZ1Qg+ZD5ePxpIJAliB38uTb6sX+5K92sY/ADuvhXYOm7dR8ct/90Ez/s58OJptE+k4wphmE91dU4IPhg0qkeSQL/cFWkgCvNoyOZk+rIpBb8kgoJfpIF8zFJPb1YTrksyKPhFGhgt9ejkrpweFPwiDYyc3M01qPHnFPySDAp+kQY0jl9ONwp+kQZGgl+lHjlNKPhFGiiUa6RTRjZtU+7Xl0uPnA8Q6WYKfpEG8uUqvZkUZlMHf2/Y4x/7G0aR7qPgF2mg0exbkegibsWKev3S3RT8Ig00mn0rMjrhuur80t0U/CINFCrVKS/JHNEsXJIUCn6RBgqleKWe6JIOCn7pdgp+kQaCHn+cGr+mX5RkUPCLNFAo16acbzfSp+CXhFDwizSQL8Wr8feOnNzVqB7pbgp+kQYKlSo9zYzqUY9fupyCX6SBYtxST3gRNwW/dLtYwW9m681sp5ntMrMPT7D9HWY2bGb3h7d31227wcweDW83tLPxIqdCvtxcqUcXapNu13DqRTNLA58FXg3sBbab2ZYJ5s79mru/f9xzzwRuBAYAB3aEzz3YltaLnAKFcrXhBdpApR5Jjjg9/nXALnff7e4l4HbgupjH/1Pgbnc/EIb93cD61poqcuq5O4VyteG0izA6jl+jeqTbxQn+5cCeuuW94brx/szMfmtmd5jZyiafi5ltNLNBMxscHh6O0SyRmVeq1qh542vxw+hlm9Xjl27XrpO73wHWuPsFBL36LzZ7AHff5O4D7j7Q39/fpmaJTE90meWeTOO3Sipl5DIpBb90vTjBvw9YWbe8Ilw3wt2fdfdiuHgrcEnc54p0s2LM2bcimoVLkiBO8G8HzjOzc8wsB2wAttTvYGbL6havBR4OH98FXG1mi8xsEXB1uE4kEaIef5zhnNF+6vFLt2s4qsfdK2b2foLATgOb3f1BM7sJGHT3LcBfmNm1QAU4ALwjfO4BM/sYwYcHwE3ufmAG/h0iMyLfbI8/lyavWbikyzUMfgB33wpsHbfuo3WPPwJ8ZJLnbgY2T6ONIh0zOtF6vNNhvdm0rscvXU+/3BWZQhT88Us9KYoVBb90NwW/yBSiUk+ca/VAWOpRj1+6nIJfZArRyd3YpZ6MTu5K91Pwi0whKtvELfX05hT80v0U/CJTiMo2GscvpxMFv8gUCi38gEs9ful2Cn6RKRQqzdX4+3LpkfMCIt1KwS8yhZEef4zLMkM4jr9cxd1nslki06LgF5lCvlwll0mRSlms/eeEl2Y+kq/MZLNEpkXBLzKFYrlGb4wrc0YuXrkQgJ89tn+mmiQybQp+kSnEnYQlcsnqRSyck+UHDz0zg60SmR4Fv8gUgvl24wd/Jp3iVS9Ywj07h6hUdZJXupOCX2QKcefbrXfV2qUcPFHmvj8cmqFWiUyPgl9kCoVyjd4mSj0Ar3j+YrJp44cPq9wj3UnBLzKFfLna1MldgPm9WV723LO4W8EvXUrBLzKFYpM1/siVL1zC7uHj7B4+NgOtEpmeWMFvZuvNbKeZ7TKzD0+w/a/M7CEz+62Z/dDMVtdtq5rZ/eFty/jninSzQrkW+wJt9a5cuxSAHz481O4miUxbw+A3szTwWeA1wPnA9WZ2/rjdfg0MuPsFwB3A/6rblnf3i8LbtW1qt8gpEYzqaf6L8coz5/DC58znByr3SBeK84peB+xy993uXgJuB66r38Hd73H3E+HiNmBFe5spcurtGjrGnoMnWHXW3Jaef9XapQw+eZBDJ0ptbpnI9MQJ/uXAnrrlveG6ybwL+G7dcq+ZDZrZNjN7fQttFOmIz92zi95Mmre/fHXjnSdw5dolVGvOj3cOt7llItPT1pO7ZvbnwADwibrVq919AHgL8Ekze94kz90YfkAMDg/rjSKd9cT+43z7/n289dJVLJ7X09IxLlyxkMXzejS6R7pOnODfB6ysW14RrhvDzK4C/ha41t2L0Xp33xfe7wZ+DFw80R9x903uPuDuA/39/bH/ASIz4bP37CKbTrHxFc9t+RiplHHlC5dw785hShX9ile6R5zg3w6cZ2bnmFkO2ACMGZ1jZhcDtxCE/lDd+kVm1hM+XgxcBjzUrsaLzIQ9B05w56/3cf26VSw5o3dax7rq/KUcLVb41eMH2tQ6kelrGPzuXgHeD9wFPAx83d0fNLObzCwapfMJYB7wjXHDNtcCg2b2G+Ae4OPuruCXrva5Hz9Gyoz3vHLCqmRTLj93MT2ZlEb3SFfJxNnJ3bcCW8et+2jd46smed7PgRdPp4Eip9K+Q3nu2LGHN790Jc9ZML3ePgQzcl1+7mJ+8PAz3Pi68zGLd11/kZmkX+6K1Ln5x48B8N4rzm3bMa9cu5S9B/M88ox+xSvdQcEvEnr6cIGvbd/DGy9ZwfKFfW077pVrlwCo3CNdQ8EvErrl3seouvNf29jbB1h6Ri8XrFig4JeuoeAXAYaOFvjqL//AGy5ezsoz57T9+FetXcr9ew7x1OF8248t0iwFvwhw608fp1yt8b5Xtbe3H3ntBcvIpVO8Y/N29h8rNn6CyAxS8Musd8/OIb78iye57qLlnLO4tevyNPLc/nlsfsdLefLAcTZs2sbQkcKM/B2ROBT8MmuVqzX+4bsP885/2c7qs+bw1+tfMKN/77JzF/PFd67jj4fyvHnTNpV9pGMU/DIr7TuUZ8Ombdzyk9285dJVfPt9l7FsQftG8kzm0ueexZfftY79R4u86ZZfsOfAicZPEmkzBb/MOj946Bmu+dRP2fn0UT59/cX8/Rte3NIsW626ZPWZfPndl3L4RJkNm7bx5LPHT9nfFgEFv8wiRwtl/se/PsS7vzTIikV9fOcDl3PthWd3pC0XrVzIV//zyzheqvDmW7ax/Qldy0dOHXP3TrfhJAMDAz44ONjpZshpwN351eMH+NrgHrY+8BSFco23v3w1f3PN2lPay5/Mw08d4e2bf8Xw0SIvXbOI917xPF71giW6tIM0zcx2hJfAb7yvgl9OR0NHCtxx316+MbiXx/cfZ15PhtddeDYbXrqSC1cu7HTzxjhRqvC17Xv4P/fu5o+HC7xg6Xzec8Vzee0FZ5NN60u5xKPgl1nn4PESO548yPYnD7D98QP8Zu9hqjVn3Tln8qaBlVzz4ucwJxfrmoQdU67W+M5v/sjNP3mMR545xvKFfbxpYCWveP5iLlixkHRK3wJkcgp+Oa0dOlFi19Axdg0d47f7DrP98QM8OhRcAC2bNi5YsZCXP/cs/uySFTM2Ln8m1WrOPTuHuOXe3Wx/4gDusKAvy2XnnsW/P6+fy89dPCO/LpZkayb4u7sLJLNSteYMHS3wx0MFnjqc56lDBXbvP85jw8d4bOgYzx4fnbx8fk+GS9Ys4vUXL2dg9SIuXLmwK2r305FKGVeuXcqVa5dy4HiJn+3az08fHeanj+5n6wNPA7B8YR9rl53B2mXzeeFzzuCFy+az5qy5+lYgsajHLzOuVnOOFiscyZc5nC9z6ESZZ48XefZYiQPHSyOP9x8r8vThAs8cLVKtjX1dLpyT5dz+eTyvfx7nLpnH85bM5dz++Sxf1Ddrws7deWz4GD95ZD/3/eEgO58+yu7hY0T/qXqzKc5dMo9VZ85hxaI5rFzUx4pFc1gR3vflkv2BKFNTqUea5u5Uak6xUqNYrlKo1CiUqxTLNQqVKoVycDtRCm75UpV8tFyscCy8HS9WOF6scrRY4VixzOETZY4WK0z2MkunjEVzciyel+PMuTmes6CXsxf0sWxhL8sW9LJsQR9nL+jjjL6MRrpMoFCusmvoGL9/+ii/f+oIjw4dY8/BE+w9mD9pnt/5vRn65/fQP68nuJ/fw+J5PZw1N8fCOVkWzgnuF83JsaAvm/hvTrNN20s9ZrYe+BSQBm5194+P294DfAm4BHgWeLO7PxFu+wjwLqAK/IW73xXz39G0qX4FGQWP4+OWg9DzkXWOOyPLHi7X3EefEy0TrY/2CZZrPvY5VXdqNafmQRnD3am6U60F+9Zqo/uM3NeCIK6FgVytjj6nXHWqtRqVmlOpengfLJertXBdjXJ0X3FK1Rrl8FaqBNuKlRrFSpVSpUapWps0nBuZk0sztyfD/J4Mc3syzO1Js3xhH/N65rGgL8uCvixn1N0v7Mty1rwcZ83tYUFfltQs6bHPhN5smhctX8CLli8Ys75Wc/YfL7LnQJ694QfB8NEiw8eKDB8t8tAfjzB8tMjRYmXSY/dkUszvzTK/NzNym9eTYV5Plnk9afpyGebm0szpGb3vy6aDWy5FTyZNXy5Y7smk6A3vMxqp1HENg9/M0sBngVcDe4HtZrZl3Ny57wIOuvu5ZrYB+J/Am83sfILJ2f8EOBv4gZk9392r7f6HAFz9T/eSL8/IobtGyiCTSpFOGZm0kUkZ6VSKbDpYzqZS4fpoXYpcOsUZuSy5tJFNp0ZuuUyKnvCWywT75cI3aG929I3aE97PyWWYk4ve2Gnm5NL0ZtIK7i6UShlL5veyZH4vl6xeNOl++VKVgydKHDxR4vCJMgdPlDmUL3HoRJkj+TJHCsE3uaOFMkcLFfYfPRF8sytVOFGsUqrWJj32ZDIpG/O6il6L0esvel0Gr1Mbsz1T9xrOpGxkn0zdcjplwbq690I6Fb5XRt4zRtqC90zKRvcZc7Ng/7QZqVTwvktZ8N82bcE+qZF7EvWNNE6Pfx2wy913A5jZ7cB1QH3wXwf8Xfj4DuAzFvxXuA643d2LwONmtis83i/a0/yx/v4/vYipXofR/5bo/8/IPVa3zrBwW7Q+yDUL10HKovXhuvA56dTo+tF9gu3pcF0qxcgLpv5FEz2OtmVSRipV9yJNRS/A5Ly4pPv15dL05fo4u8UZx0qVGvlSNfggKFXIl4LSYFQKLJSDx9G3y6h0GN2XKrWRb5ylSo1iuHysWAm+oVaCb7Glum+rlapTrtVGvu12CzPGvs9HHo++x6PHaTMs3K8+LxbP7eHr73n5jLc1TvAvB/bULe8FLp1sH3evmNlh4Kxw/bZxz13ecmsbeMPFK2bq0CIygaiHvmBOtiN/391HSpqjpc/wQyH8gKiFZdNqWAqNyqW1sHRardsePK6F2xh9rjvVam2knFsdU5YdLfGOPK+ulOtev0xY9g0eu48+rrozv+fUDLTsmuGcZrYR2AiwatWqDrdGRJLAzMhljJwuO9aUOP+19gEr65ZXhOsm3MfMMsACgpO8cZ4LgLtvcvcBdx/o7++P13oREWlanODfDpxnZueYWY7gZO2WcftsAW4IH78R+JEH40S3ABvMrMfMzgHOA37VnqaLiEgrGpZ6wpr9+4G7CIZzbnb3B83sJmDQ3bcAnwe+HJ68PUDw4UC439cJTgRXgPfN1IgeERGJRz/gEhE5DTTzAy6dERERmWUU/CIis4yCX0RkllHwi4jMMl15ctfMhoEnW3z6YmB/G5vTTmpba9S21qhtrUlq21a7e6wfQXVl8E+HmQ3GPbN9qqltrVHbWqO2tWY2tE2lHhGRWUbBLyIyy5yOwb+p0w2YgtrWGrWtNWpba077tp12NX4REZna6djjFxGRKZw2wW9m681sp5ntMrMPd7o99cxss5kNmdnvOt2Wema20szuMbOHzOxBM/tgp9tUz8x6zexXZvabsH3/vdNtqmdmaTP7tZn9a6fbMp6ZPWFmD5jZ/WbWVRe+MrOFZnaHmf3ezB42s5mfcioGM3tB+N8ruh0xsw91ul0RM/vL8H3wOzO7zcx6Wz7W6VDqCecFfoS6eYGB68fNC9wxZvYK4BjwJXd/UafbEzGzZcAyd7/PzOYDO4DXd9H/73phAAAC9ElEQVR/NwPmuvsxM8sC/wZ80N23NXjqKWFmfwUMAGe4+2s73Z56ZvYEMODuXTce3cy+CPzU3W8NL/U+x90Pdbpd9cJM2Qdc6u6t/qaone1ZTvD6P9/d8+FVj7e6+xdaOd7p0uMfmRfY3UtANC9wV3D3ewkuV91V3P0pd78vfHwUeJgZnBqzWR44Fi5mw1tX9FTMbAXwH4FbO92WJDGzBcArCC7ljruXui30Q1cCj3VD6NfJAH3hZFdzgD+2eqDTJfgnmhe4awIsCcxsDXAx8MvOtmSssJxyPzAE3O3u3dK+TwJ/DdQ63ZBJOPB9M9sRTmvaLc4BhoF/Cctkt5rZ3E43agIbgNs63YiIu+8D/jfwB+Ap4LC7f7/V450uwS/TYGbzgG8CH3L3I51uTz13r7r7RQTTdq4zs46XyszstcCQu+/odFumcLm7vwR4DfC+sNzYDTLAS4B/dveLgeNAt52TywHXAt/odFsiZraIoIpxDnA2MNfM/rzV450uwR97bl8ZK6ydfxP4irt/q9PtmUxYDrgHWN/ptgCXAdeGdfTbgf9gZv+3s00aK+wh4u5DwJ0E5dBusBfYW/fN7Q6CD4Ju8hrgPnd/ptMNqXMV8Li7D7t7GfgW8O9aPdjpEvxx5gWWccKTp58HHnb3f+x0e8Yzs34zWxg+7iM4ef/7zrYK3P0j7r7C3dcQvNZ+5O4t977azczmhifrCcsoVwNdMaLM3Z8G9pjZC8JVVxJMzdpNrqeLyjyhPwAvM7M54fv2SoJzci1pOOduEkw2L3CHmzXCzG4DrgAWm9le4EZ3/3xnWwUEPde3AQ+EdXSAv3H3rR1sU71lwBfDERYp4Ovu3nVDJ7vQUuDOIB/IAF919+91tkljfAD4SthJ2w28s8PtGRF+UL4a+C+dbks9d/+lmd0B3Ecwf/mvmcaveE+L4ZwiIhLf6VLqERGRmBT8IiKzjIJfRGSWUfCLiMwyCn4RkVlGwS8iMsso+EVEZhkFv4jILPP/AW4wVtnzRzKTAAAAAElFTkSuQmCC\n"
},
"metadata": {}
}
]
}
],
"metadata": {
"kernelspec": {
"name": "python3",
"display_name": "Python 3",
"language": "python"
},
"language_info": {
"name": "python",
"version": "3.6.3",
"mimetype": "text/x-python",
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"pygments_lexer": "ipython3",
"nbconvert_exporter": "python",
"file_extension": ".py"
},
"gist": {
"id": "",
"data": {
"description": "Desktop/PPP/python/audio/Untitled2.ipynb",
"public": true
}
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment