Created
February 21, 2020 09:37
-
-
Save prl900/18092298ec5f772be87bde010394dca8 to your computer and use it in GitHub Desktop.
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| { | |
| "cells": [ | |
| { | |
| "cell_type": "code", | |
| "execution_count": 15, | |
| "metadata": {}, | |
| "outputs": [ | |
| { | |
| "data": { | |
| "text/plain": [ | |
| "[<matplotlib.lines.Line2D at 0x7f07fc431be0>]" | |
| ] | |
| }, | |
| "execution_count": 15, | |
| "metadata": {}, | |
| "output_type": "execute_result" | |
| }, | |
| { | |
| "data": { | |
| "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD4CAYAAAD8Zh1EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAZN0lEQVR4nO3dfXRc9X3n8fd3RpbkBz1gW5aN5SfAELsOxkShpEkpC3kwtAvtpmmhS9vdzdb7UNq0ZdtDTlrSZU+7Tcvp0zk0uyZNk7AsLmHTrkvdZRPCQ7YJ1CI2DrYxEbKRZVsP2NZItjTSSPPdP+bKTGQ9XI1mNHNnPq9z5lj33p/ufK+v9PH1b+79/czdERGR6IsVuwAREckPBbqISJlQoIuIlAkFuohImVCgi4iUiapivfHKlSt948aNxXp7EZFIevXVV99x96apthUt0Ddu3EhbW1ux3l5EJJLM7O3ptqnLRUSkTCjQRUTKhAJdRKRMKNBFRMqEAl1EpEzMGuhm9kUz6zWz16fZbmb252bWbmaHzOzG/JcpIiKzCXOF/iVg5wzb7wA2B69dwOfnX5aIiMzVrPehu/tLZrZxhiZ3A1/xzDi8L5tZo5mtcfczeapRJNKSqXH6BkfoGUjyzoURLo6MkxwbJ5lKMzaeZtyddNpJO7hDOhjS+rKBrTXUddm4fUsz29c15n2/+XiwaC1wMmu5K1h3WaCb2S4yV/GsX78+D28tUnrSaaft7fO89GYf3/p+H4dOJfKWxWb52Y8U16r62pIN9NDcfTewG6C1tVWXG1J2+gZH+E9ffY0X3+wjHjNuWNfI/f/sGtYtX8Kquhqa6mpYVlPF4kVxaqriVMWNeMyImREziJlhBqbklhzkI9BPAeuylluCdSIV5cU3+3jgqYMMJsf47D/fysff10J97aJilyUVJB+Bvhe438z2AD8MJNR/LpVm72un+dUnD3Bt8zKe+Lc3c93qumKXJBVo1kA3syeBW4GVZtYFfBZYBODu/w3YB9wJtANDwL8uVLEipSiZGucP9h3l+pYGnvp3H6B2UbzYJUmFCnOXy72zbHfgl/NWkUjEPPFKJ6cTSR75xHaFuRSVnhQVmYcLI2M8+nw7H7pmJT9yzcpilyMVToEuMg9f+FYH5y6O8psfu67YpYgo0EVydfbCCI+91MEd21YX5J5ikblSoIvk6IlXOhlKjfPAR68tdikigAJdJGf7T5xjy+p6rlmlWxSlNCjQRXKQTjsHT/Zzw3p1tUjpUKCL5KDjnYsMJsfYob5zKSEKdJEcHOg8D8AOXaFLCVGgi+Tg4Ml+6mqruGrlsmKXInKJAl0kBwc6+7lhXSOxmEZFlNKhQBeZo6HRMY71DHKD+s+lxCjQReboe10JxtOu/nMpOQp0kTk6eLIfgO0tCnQpLQp0kTk60NnPhhVLWLGsptiliPwABbrIHB082a/+cylJCnSROTiTGKZ7IKkHiqQkKdBF5uBgZ6b//Ib1VxS5EpHLKdBF5uBgVz/V8Rhb1mhALik9CnSROTh5boiW5YupqdJUc1J6FOgic9CdSLK6vrbYZYhMSYEuMgc9AyMKdClZCnSRkNJpp3cwySoFupQoBbpISOeGRkmNO6vr9UCRlCYFukhIPQNJAFY36ApdSpMCXSSkiUBXl4uUKgW6SEg9AyMA+lBUSpYCXSSk7kQSM2iqUx+6lCYFukhIvYNJViytYVFcvzZSmvSTKRJSdyLJ6gZdnUvpUqCLhNQ9MEJznfrPpXQp0EVC6h1I0qxbFqWEhQp0M9tpZsfMrN3MHpxi+3oze97MDpjZITO7M/+lihTPyNg4Zy+O6gpdStqsgW5mceBR4A5gK3CvmW2d1Oy3gafcfQdwD/AX+S5UpJj6BoNbFtWHLiUszBX6TUC7u3e4+yiwB7h7UhsH6oOvG4DT+StRpPgmHipq1j3oUsLCBPpa4GTWclewLtvvAveZWRewD/iVqXZkZrvMrM3M2vr6+nIoV6Q4Jh4qUqBLKcvXh6L3Al9y9xbgTuBxM7ts3+6+291b3b21qakpT28tUnjdiWAcFwW6lLAwgX4KWJe13BKsy/ZJ4CkAd/8OUAuszEeBIqWgZyBJdVWMxiWLil2KyLTCBPp+YLOZbTKzajIfeu6d1KYTuB3AzLaQCXT1qUjZ6BlI0lxfg5kVuxSRac0a6O4+BtwPPAscJXM3y2Eze9jM7gqaPQD8kpm9BjwJ/Ct390IVLbLQugc09ZyUvqowjdx9H5kPO7PXPZT19RHgg/ktTaR09A6MsOXK+tkbihSRnhQVmYW76wpdIkGBLjKLCyNjDI2O06yp56TEKdBFZqGHiiQqFOgis+hOaKYiiQYFusgsdIUuUaFAF5lFtwJdIkKBLjKL3oEk9bVVLK6OF7sUkRkp0EVm0Ts4wipdnUsEKNBFZtE/lOIKjeEiEaBAF5lFYjhFw2IFupQ+BbrILBLDKeoV6BIBCnSRWQwMp2hcXF3sMkRmpUAXmcHYeJrBkTF1uUgkKNBFZjCQHAOgYXGogUlFikqBLjKDxHAKgAbd5SIRoEAXmcGlQFeXi0SAAl1kBv1DowA06ENRiQAFusgMdIUuUaJAF5nBgAJdIkSBLjIDXaFLlCjQRWbQP5Ri8aI41VX6VZHSp59SkRkkhlM06pZFiQgFusgMNDCXRIkCXWQGGphLokSBLjIDXaFLlCjQRWagQJcoUaCLzCAxnKJRgS4RoUAXmUZqPM3Q6Liu0CUyFOgi09BIixI1CnSRafQP6SlRiZZQgW5mO83smJm1m9mD07T5GTM7YmaHzex/5rdMkYWnx/4lamadhsXM4sCjwEeALmC/me119yNZbTYDnwY+6O7nzWxVoQoWWSgamEuiJswV+k1Au7t3uPsosAe4e1KbXwIedffzAO7em98yRRaertAlasIE+lrgZNZyV7Au27XAtWb2j2b2spntnGpHZrbLzNrMrK2vry+3ikUWiAJdoiZfH4pWAZuBW4F7gcfMrHFyI3ff7e6t7t7a1NSUp7cWKYyJD0X16L9ERZhAPwWsy1puCdZl6wL2unvK3Y8Db5IJeJHISgynWFZTxaK4bgaTaAjzk7of2Gxmm8ysGrgH2Dupzd+SuTrHzFaS6YLpyGOdIgtOj/1L1Mwa6O4+BtwPPAscBZ5y98Nm9rCZ3RU0exY4a2ZHgOeB33T3s4UqWmQhaKRFiZpZb1sEcPd9wL5J6x7K+tqB3wheImUhMTxKw+JQvyIiJUGdgyLTUJeLRI0CXWQamZEWq4tdhkhoCnSRaSSGUxqYSyJFgS4yhWRqnGQqrS4XiRQFusgUJsZx0V0uEiUKdJEpTDz2r9mKJEoU6CJT0DguEkUKdJEpKNAlihToIlPQbEUSRQp0kSnoCl2iSIEuMoWE7nKRCFKgi0whMZyirraKeMyKXYpIaAp0kSloHBeJIgW6yBQU6BJFCnSRKSjQJYoU6CJTSAynaNTAXBIxCnSRKegKXaJIgS4yibuTGNL0cxI9CnSRSZKpNKPjaU1uIZGjQBeZRE+JSlQp0EUmUaBLVCnQRSbpHxoFFOgSPQp0kUkuTW6h2xYlYhToIpOoy0WiSoEuMolGWpSoUqCLTJIYTmEGdTVVxS5FZE4U6CKTJIZT1NcuIqahcyViFOgik2gcF4kqBbrIJBrHRaJKgS4ySf+QAl2iSYEuMsnAsAbmkmgKFehmttPMjplZu5k9OEO7j5uZm1lr/koUWViJ4RSNCnSJoFkD3cziwKPAHcBW4F4z2zpFuzrgU8Ar+S5SZKG4u/rQJbLCXKHfBLS7e4e7jwJ7gLunaPdfgM8ByTzWJ7KghkbHGUu7Al0iKUygrwVOZi13BesuMbMbgXXu/vcz7cjMdplZm5m19fX1zblYkULr12P/EmHz/lDUzGLAHwMPzNbW3Xe7e6u7tzY1Nc33rUXyLjGkgbkkusIE+ilgXdZyS7BuQh2wDXjBzE4ANwN79cGoRJHGcZEoCxPo+4HNZrbJzKqBe4C9ExvdPeHuK919o7tvBF4G7nL3toJULFJAGmlRomzWQHf3MeB+4FngKPCUux82s4fN7K5CFyiykBLDmtxCoivUcHLuvg/YN2ndQ9O0vXX+ZYkUx7uTW2iCaIkePSkqkiUxnCIeM5ZWx4tdisicKdBFskw8VGSmoXMlehToIlk0MJdEmQJdJIse+5coU6CLZBlQoEuEKdBFsugKXaJMgS6SpV+BLhGmQBcJpNPOgOYTlQhToIsELoyOkXY9JSrRpUAXCUyMtKiBuSSqFOgiAQ3MJVGnQBcJXBrHRYEuEaVAFwlcukLXh6ISUQp0kYC6XCTqFOgiAQW6RJ0CXSTQP5SiOh5j8SINnSvRpEAXCSSGU9Rr6FyJMAW6SCAzMFeoSbxESpICXSSggbkk6hToIoH+4VEFukSaAl0k0Dc4QlNdTbHLEMmZAl0EGBtP0zc4wur62mKXIpIzBboI8M6FUdIOzQ0KdIkuBboI0DOQBKC5ToEu0aVAFwG6g0BfrSt0iTAFugjQGwT6qnp9KCrRpUAXIXOFHo8ZK5cGgd75Crz9neIWJTJHeixOBOhOjLCqroZYzCA9Dk//G0iPwa8fhrh+TSQadIUuAvQOJmmeuGWx4wUY6IIL3fDWc0WtS2QuFOgiQHci+e496Af+Byy+ApY2wYHHi1uYyByECnQz22lmx8ys3cwenGL7b5jZETM7ZGbPmdmG/JcqUjg9A0ma62tg6By88Qxc/7OZ17F/gIvvFLs8kVBmDXQziwOPAncAW4F7zWzrpGYHgFZ3vx54GvjDfBcqUijDo+MMJMcyDxV972kYH4Ud92Ve6TE49NfFLlEklDBX6DcB7e7e4e6jwB7g7uwG7v68uw8Fiy8DLfktU6RwurMfKjrwOKzZDqvfC6u2wNpW+O7j4F7kKkVmFybQ1wIns5a7gnXT+STwD1NtMLNdZtZmZm19fX3hqxQpoImnRK8e74DuQ7Dj59/duOM+6DsKp79bpOpEwsvrh6Jmdh/QCvzRVNvdfbe7t7p7a1NTUz7fWiRnE4G+8eTfQLwGtn383Y3b/gVULYYDTxSpOpHwwgT6KWBd1nJLsO4HmNmHgc8Ad7n7SH7KEym8iUCv7/wGbP4ILFn+7sbaBrhuJxzbV6TqRMILE+j7gc1mtsnMqoF7gL3ZDcxsB/DfyYR5b/7LFCmc7sQIV1YPE0t0Qsv7L2/Q8n4YPAODPQtfnMgczBro7j4G3A88CxwFnnL3w2b2sJndFTT7I2AZ8FUzO2hme6fZnUjJ6RlI8oElwX8612y/vMHEuu5DC1eUSA5CPdPs7vuAfZPWPZT19YfzXJfIgukZSPJT1Z2QZOpAX/3ezJ9nDma6ZERKlJ4UlYrXPZBki3dAw/of7D+fUNsAy6+CM68tfHEic6BAl4rm7vQOjLAh1Q5rrp++4ZrtCnQpeQp0qWjnh1JUj19gRbIT1twwfcM126G/MzM0gEiJUqBLRetOJNlinZmFqfrPJ+iDUYkABbpUtJ7BJNtixzMLMwX66mCbul2khCnQpaL1JJJsi51gfGkz1DVP33DpCmhYp0CXkqZAl4rWMzDCD9kJbKar8wn6YFRKnAJdKtrZ/n42x04Ru3KGD0QnrNkOZ9shOVD4wkRyoECXilZ99ihx0jP3n0+YaNPzemGLEsmRAl0q2hWJo5kv5hLo6naREqVAl4qVTjvNQ29wMd4ADSHmZKlbDcuaFehSshToUrE63rnIdekOBq/YCmbhvmnNdjh9sLCFieRIgS4V6/WOLrZYJ4s23BT+m9a2Qt8bemJUSpICXSrWwLEXqLI0V2y9Pfw3bboFcHj7HwtWl0iuFOhSsepOf5tRqomt/+Hw37T2fbBoCXS8WLjCRHKkQJeKNDw6znXDBznTsB0W1Yb/xqpq2PAjcPylwhUnkiMFulSko28dZ2vsbVLrPjT3b950C7xzDAa781+YyDwo0KUinX39GwCsvD6HGYg23ZL58/i38liRyPwp0KUi1Zz8f1xkMY1Xz6H/fMLq6zOzGB1XP7qUFgW6VKQNA69yfOl2iIeaVvcHxeKw8UfVjy4lR4EuFaen6y02cJqLV34w951sugX634bzJ/JWl8h8KdCl4nQf/L8ANGy9LfedXOpH11W6lA4FulQcO/ES572OTdty6D+f0PQeWLpKgS4lRYEulSU9ztpz/8SRmu3ULFqU+37MMlfpHS/C+Fj+6hOZBwW6VJSB/U+yIv0Op1vumP/O3vvTcLEXDu2Z/75E8kCBLpVjPMXYN3+fI+kNvG/nL8x/f9fuhCtvhBc+B2Oj89+fyDwp0KVinPv2l1k+cor9V/0HrlpVP/8dmsFtn4FEJxz4yvz3JzJPCnSpDGMj8NIf8ppfzUd/8hfzt9+rb4d1N8NLj0BqOH/7FcmBAl0qQs+Lj7E81cPh636FNY1L8rdjM7jtt2HwDLT9Vf72K5IDBbqUPe97k9pvP0Kbb2HnXT+X/zfY9KOw6cfgW49A79H8718kJAW6lLXBQ88w/Be3khob563W32H5sprCvNHOPwCLwxc+DG/8fWHeQ2QWoQLdzHaa2TEzazezB6fYXmNmfx1sf8XMNua7UJG58KHzdP7NZ1n6tfvoSK/i6x/awyd+/M7CvWHzVtj1AqzcDHt+Dr75e5qmThacufvMDcziwJvAR4AuYD9wr7sfyWrzH4Hr3f3fm9k9wE+5+8/OtN/W1lZva2ubb/1S6dJpRob6OX/mBIO9J0ieOUbt8a+z6eIBqkjzXNUtrP75x/ihDasXpp7UMPzdr2XuTbd4ZjKMa3fCymuhoQXqr4SauswAXyI5MLNX3b11ym0hAv0DwO+6+8eC5U8DuPt/zWrzbNDmO2ZWBXQDTT7DznMN9P1f+zOaXn9szt8nxWehWzpG5kfHfOJrJ4YT83HijFPFOLWeZLFdfv/3W6zlzcZb4D0/wS23foyltfN4IjQX7nD6QKbr5Y1nMpNKT1ZVm5nKLl4NsSqIxTL/AJiBxbj0t2VG6L85C/83LEX2Y78F2z6e07fOFOhhxg5dC5zMWu4CJg+CcamNu4+ZWQJYAbwzqZBdwC6A9evXhyr+soKXreDckk05fa+UgnCh49ilgHJiYJaJ+dgi3OJ4LE66agletQSvWUasYS01KzZQ17yJq666hqvjRfx4yAzW3ph53f47mZmN+k9C4iQMnIbRC8FrCNJjkB7P/Onpd1/BkTPLBde7wraTklDbWJDd5jAYdO7cfTewGzJX6LnsY8dH74OP3pfXukQKqm515rXu/cWuRMpcmMuYU8C6rOWWYN2UbYIulwbgbD4KFBGRcMIE+n5gs5ltMrNq4B5g76Q2e4GJx+9+GvjmTP3nIiKSf7N2uQR94vcDzwJx4IvuftjMHgba3H0v8JfA42bWDpwjE/oiIrKAQvWhu/s+YN+kdQ9lfZ0EPpHf0kREZC70pKiISJlQoIuIlAkFuohImVCgi4iUiVkf/S/YG5v1AW/n+O0rmfQUaoWoxOOuxGOGyjzuSjxmmPtxb3D3pqk2FC3Q58PM2qYby6CcVeJxV+IxQ2UedyUeM+T3uNXlIiJSJhToIiJlIqqBvrvYBRRJJR53JR4zVOZxV+IxQx6PO5J96CIicrmoXqGLiMgkCnQRkTIRuUCfbcLqcmBm68zseTM7YmaHzexTwfrlZvZ1M/t+8OcVxa4138wsbmYHzOyZYHlTMPF4ezAReXWxa8w3M2s0s6fN7A0zO2pmH6iQc/3rwc/362b2pJnVltv5NrMvmlmvmb2etW7Kc2sZfx4c+yEzu3Gu7xepQA8mrH4UuAPYCtxrZluLW1VBjAEPuPtW4Gbgl4PjfBB4zt03A88Fy+XmU8DRrOXPAX/i7tcA54FPFqWqwvoz4P+4+3uA7WSOv6zPtZmtBX4VaHX3bWSG5r6H8jvfXwJ2Tlo33bm9A9gcvHYBn5/rm0Uq0IGbgHZ373D3UWAPcHeRa8o7dz/j7t8Nvh4k8wu+lsyxfjlo9mXgJ4tTYWGYWQvw48AXgmUDbgOeDpqU4zE3ALeQmVMAdx91937K/FwHqoDFwSxnS4AzlNn5dveXyMwRkW26c3s38BXPeBloNLM1c3m/qAX6VBNWry1SLQvCzDYCO4BXgGZ3PxNs6gaai1RWofwp8FvAxCzJK4B+dx8LlsvxfG8C+oC/CrqavmBmSynzc+3up4BHgE4yQZ4AXqX8zzdMf27nnW9RC/SKYmbLgP8F/Jq7D2RvC6b4K5t7Ts3sJ4Bed3+12LUssCrgRuDz7r4DuMik7pVyO9cAQb/x3WT+QbsSWMrlXRNlL9/nNmqBHmbC6rJgZovIhPkT7v61YHXPxH/Bgj97i1VfAXwQuMvMTpDpSruNTN9yY/BfcijP890FdLn7K8Hy02QCvpzPNcCHgePu3ufuKeBrZH4Gyv18w/Tndt75FrVADzNhdeQFfcd/CRx19z/O2pQ9GfcvAv97oWsrFHf/tLu3uPtGMuf1m+7+L4HnyUw8DmV2zADu3g2cNLPrglW3A0co43Md6ARuNrMlwc/7xHGX9fkOTHdu9wK/ENztcjOQyOqaCcfdI/UC7gTeBN4CPlPsegp0jB8i89+wQ8DB4HUnmT7l54DvA98Alhe71gId/63AM8HXVwH/BLQDXwVqil1fAY73BqAtON9/C1xRCeca+M/AG8DrwONATbmdb+BJMp8RpMj8b+yT051bwMjcxfcW8D0ydwDN6f306L+ISJmIWpeLiIhMQ4EuIlImFOgiImVCgS4iUiYU6CIiZUKBLiJSJhToIiJl4v8DW2mOjDOKiD4AAAAASUVORK5CYII=\n", | |
| "text/plain": [ | |
| "<Figure size 432x288 with 1 Axes>" | |
| ] | |
| }, | |
| "metadata": { | |
| "needs_background": "light" | |
| }, | |
| "output_type": "display_data" | |
| } | |
| ], | |
| "source": [ | |
| "%matplotlib inline\n", | |
| "\n", | |
| "import numpy as np\n", | |
| "from matplotlib import pyplot as plt\n", | |
| "\n", | |
| "# Sigmoid centered in 50\n", | |
| "def sigmoid(x):\n", | |
| " return (1 / (1 + np.exp(50-x)))\n", | |
| "\n", | |
| "x = np.arange(100)\n", | |
| "y = sigmoid(x)\n", | |
| "dy = np.gradient(y)\n", | |
| "\n", | |
| "\n", | |
| "# Plotting curve and its gradient:\n", | |
| "plt.plot(x,y)\n", | |
| "plt.plot(x,np.gradient(y))" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 16, | |
| "metadata": {}, | |
| "outputs": [ | |
| { | |
| "name": "stdout", | |
| "output_type": "stream", | |
| "text": [ | |
| "[47 48 49 50 51 52 53]\n" | |
| ] | |
| }, | |
| { | |
| "data": { | |
| "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD4CAYAAAD8Zh1EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAWC0lEQVR4nO3dfZBd9X3f8fd3792V0LNBK8mWhIVBuCgeG9EdguNMQ2PHAdKiZvJQmHpqd9yo44babWg6ZJIhKf2nTjJJ6il1gx2XxEOh2HUTTaKGqR0Se9pAJSoZGzAgiI0ESFoe9kog9LC73/5xz4rLsg93pbt795z7fs3s6J4H3fs9c5YPX/3uOb8TmYkkqfz6ul2AJKkzDHRJqggDXZIqwkCXpIow0CWpIurd+uC1a9fmli1buvXx6nUvvACnT5//+wwMwLvedf7vI7XpkUceeSkzB6fa1rVA37JlC3v37u3Wx6vXfeIT0ImG4vvfh7vvPv/3kdoUET+YbptDLpJUEQa6JFWEgS5JFWGgS1JFGOiSVBGzBnpEfCkijkbEd6fZHhHxuYg4EBGPRsRVnS9TkjSbdjr0u4HrZth+PbC1+NkJfP78y5IkzdWs16Fn5jcjYssMu+wA/iib8/A+FBFrIuKdmflih2qUSu3kmTGGj5/iyLGTvPTaKV4/NcbJ0TFOnhlndGycsUzGx5PxhEwYL6a0ftvE1k51XRkfvmI9H9i8puPv24kbizYCB1uWDxXr3hboEbGTZhfPxRdf3IGPlhaf8fFk7w9e5ZtPDfOtp4d59PlGx7I4ojPvo+5at2rpog30tmXmXcBdAENDQ7Ybqpzh46f411/5Nn/11DC1vuDKzWu45e9exuYLl7Fu5RIGVy5hxZI6F/TXWFKvUa8Ftb6gL4K+gL4IIiBMbp2DTgT688DmluVNxTqpp/zVU8Pcev9+jp8c5df//jZ+5m9vYtXS/m6XpR7SiUDfBdwSEfcBPww0HD9Xr9n17Rf49L37uHz9Cu75p9fw3g0ru12SetCsgR4R9wLXAmsj4hDw60A/QGb+Z2A3cANwADgB/JP5KlZajE6eGePf736C929azf3/7IMs7a91uyT1qHaucrl5lu0J/GLHKpJK5p6Hn+OFxkl+++c+YJirq7xTVDoPr0WdOx88wI9etpYfuWxtt8tRjzPQpfPwxWWX88rrp/nln3xvt0uRDHTpXL2cdb6w/HKuf9+GebmmWJorA106R/eMDXIi6tz60cu7XYoEGOjSOdszvpIrRke4bJ2XKGpxMNClczCesH98OVeeeaXbpUhnGejSOXg2l3KcOtvPvNztUqSzDHTpHOwbXw7A9tN26Fo8DHTpHOzPFaxklPeMHe92KdJZBrp0DvaNL+fKvtf9D0iLir+P0hydyD6ezGVcGa93uxTpLQx0aY6+k8sYI9je91q3S5HewkCX5mj/+AoAPtBnh67FxUCX5mjf+HLeHSe5KEa7XYr0Fga6NEf7x1c4fq5FyUCX5uDF7OcwA46fa1Ey0KU5mBg/v9Lxcy1CBro0B/vHlzPAOFfEiW6XIr2NgS7NwcFcwqY4xZLIbpcivY2BLs3B4RxgQ5zpdhnSlAx0aQ6OZD8bON3tMqQpGehSm8YTjtLPOjt0LVIGutSmV6hzhj42hB26FicDXWrTkewHcAxdi5aBLrXpSA4AsM4OXYuUgS61yQ5di52BLrXpcA4QJIMY6FqcDHSpTUfp5yJG6femIi1SBrrUpuZNRY6fa/Ey0KU2Hc4B1jt+rkXMQJfadDT7WW+HrkWsrUCPiOsi4smIOBARt02x/eKIeDAi9kXEoxFxQ+dLlbrnVAYv02+HrkVt1kCPiBpwJ3A9sA24OSK2Tdrt14D7M3M7cBPwnzpdqNRNwxSXLDqPixaxdjr0q4EDmflsZp4G7gN2TNongVXF69XAC50rUeq+iZuKHHLRYtZOoG8EDrYsHyrWtfoN4GMRcQjYDfyLqd4oInZGxN6I2Ds8PHwO5UrdMXFTkUMuWsw69aXozcDdmbkJuAH4ckS87b0z867MHMrMocHBwQ59tDT/DhcdupctajFrJ9CfBza3LG8q1rX6JHA/QGb+NbAUWNuJAqXF4Ej2M8A4axjrdinStNoJ9D3A1oi4JCIGaH7puWvSPs8BHwaIiCtoBrpjKqqMIznA+jhNRLcrkaY3a6Bn5ihwC/AA8ATNq1kei4g7IuLGYrdbgV+IiG8D9wKfyEzvj1ZlHKafDc7hokWu3s5Ombmb5pedretub3n9OPChzpYmLR5Hs58r4o1ulyHNyDtFpVlkOo+LysFAl2bxGn2coOYli1r0DHRpFt5UpLIw0KVZeA26ysJAl2ZxpJjHZb1XuWiRM9ClWRw+O+RioGtxM9ClWRzNflYxygUx3u1SpBkZ6NIsjmY/6+zOVQIGujSLEeq8g9FulyHNykCXZtHIGqvDQNfiZ6BLs2hknVXOsqgSMNClWRyjxho7dJWAgS7NYDThOHVWhx26Fj8DXZrBsWJC0tV+KaoSMNClGTSyBmCHrlIw0KUZNOzQVSIGujSDkbMduoGuxc9Al2bwZofukIsWPwNdmsExO3SViIEuzcAOXWVioEszGMk6FzDGQGS3S5FmZaBLM2hQY43duUrCQJdm0Mi64+cqDQNdmkEja6zyGnSVhIEuzaDhPC4qEQNdmoFDLioTA12agV+KqkwMdGkaZzI4gU8rUnkY6NI0GhR3idqhqyQMdGkaI1ncJWqHrpJoK9Aj4rqIeDIiDkTEbdPs8/MR8XhEPBYR/7WzZUoL780O3UBXOdRn2yEiasCdwE8Ah4A9EbErMx9v2Wcr8CvAhzLz1YhYN18FSwvl2NkO3SEXlUM7HfrVwIHMfDYzTwP3ATsm7fMLwJ2Z+SpAZh7tbJnSwrNDV9m0E+gbgYMty4eKda0uBy6PiP8dEQ9FxHVTvVFE7IyIvRGxd3h4+NwqlhZIww5dJdOpL0XrwFbgWuBm4AsRsWbyTpl5V2YOZebQ4OBghz5amh8jxYjkKq9yUUm0E+jPA5tbljcV61odAnZl5pnM/BvgKZoBL5VWI2usYIx+p85VSbQT6HuArRFxSUQMADcBuybt88c0u3MiYi3NIZhnO1intOAaWXf8XKUya6Bn5ihwC/AA8ARwf2Y+FhF3RMSNxW4PAC9HxOPAg8AvZ+bL81W0tBAa1Fjl+LlKZNbLFgEyczewe9K621teJ/BLxY9UCU7MpbLxTlFpGg1q3vavUjHQpWk0ss4aO3SViIEuTaPZoRvoKg8DXZrCyQxOUvOmIpWKgS5N4djZm4rs0FUeBro0hUY253FZY4euEjHQpSk0ig7dMXSViYEuTWGiQ3cMXWVioEtTGLFDVwkZ6NIUGj5+TiVkoEtTmBhycepclYmBLk2hQZ2VjFKLblcitc9Al6bQSG8qUvkY6NIUGjgXusrHQJemYIeuMjLQpSk0qLPGDl0lY6BLU2h26Aa6ysVAlybJbHboXrKosjHQpUlO0sdp+ny4hUrHQJcmaVDM42KHrpIx0KVJvO1fZWWgS5OM2KGrpAx0aZKJDt0xdJWNgS5NcnYudK9DV8kY6NIkE08rWuWdoioZA12apJF1gmSlY+gqGQNdmqRBjVWM0efUuSoZA12apJF1vxBVKRno0iQNal6yqFIy0KVJRrLuTUUqJQNdmuRYMYYulU1bgR4R10XEkxFxICJum2G/n4mIjIihzpUoLSzH0FVWswZ6RNSAO4HrgW3AzRGxbYr9VgKfAR7udJHSQmlOnVvzpiKVUjsd+tXAgcx8NjNPA/cBO6bY798BnwVOdrA+aUGdoI9R+nz8nEqpnUDfCBxsWT5UrDsrIq4CNmfmn830RhGxMyL2RsTe4eHhORcrzbeR4i5RO3SV0Xl/KRoRfcDvALfOtm9m3pWZQ5k5NDg4eL4fLXXcxDwua+zQVULtBPrzwOaW5U3FugkrgfcBfxkR3weuAXb5xajKaGKmxVV26CqhdgJ9D7A1Ii6JiAHgJmDXxMbMbGTm2szckplbgIeAGzNz77xULM2js08r8ioXldCsgZ6Zo8AtwAPAE8D9mflYRNwRETfOd4HSQnrzaUUOuah86u3slJm7gd2T1t0+zb7Xnn9ZUndMTJ27xiEXlZB3ikotGlmjRrKc8W6XIs2ZgS61aFBnNaOEU+eqhAx0qcVI1hw/V2kZ6FKLiQ5dKiMDXWpxLGtesqjSMtClFs0O3SEXlZOBLrXw4RYqMwNdKoxn8+EWa+zQVVIGulR4jRrjhB26SstAlwoTMy36+DmVlYEuFSZu+7dDV1kZ6FLhzbnQDXSVk4EuFc526A65qKQMdKkw0aE75KKyMtClgh26ys5AlwojWWeAcS5w6lyVlIEuFRrUWMWYU+eqtAx0qXDM2/5Vcga6VGhQc/xcpWagSwUn5lLZGehSYTj7GYwz3S5DOmcGugSMJgzTzwZOd7sU6ZwZ6BLwEv2ME6y3Q1eJGegScCQHAFgfdugqLwNdAg5nPwAb7NBVYga6BBwtAn2dHbpKzECXgMM5QI1kLV62qPIy0CWagb6OM/R5279KzECXgKP0+4WoSs9Al2h26BsMdJVcW4EeEddFxJMRcSAibpti+y9FxOMR8WhEfCMi3t35UqX5cyT7vQZdpTdroEdEDbgTuB7YBtwcEdsm7bYPGMrM9wNfBX6z04VK8+WN7OMYdYdcVHrtdOhXAwcy89nMPA3cB+xo3SEzH8zME8XiQ8CmzpYpzZ+Ja9Dt0FV27QT6RuBgy/KhYt10Pgn8z6k2RMTOiNgbEXuHh4fbr1KaR0do3iXqPC4qu45+KRoRHwOGgN+aantm3pWZQ5k5NDg42MmPls7ZETt0VUS9jX2eBza3LG8q1r1FRHwE+FXgxzLzVGfKk+af87ioKtrp0PcAWyPikogYAG4CdrXuEBHbgd8HbszMo50vU5o/h7Of5YyxMnw4tMpt1kDPzFHgFuAB4Ang/sx8LCLuiIgbi91+C1gBfCUi9kfErmneTlp0juSA3bkqoZ0hFzJzN7B70rrbW15/pMN1SQvGa9BVFd4pqp53mAHWY6Cr/Ax09bTM5tS5DrmoCgx09bRXqXOaPodcVAkGunram08qskNX+Rno6mleg64qMdDV07xLVFVioKunTczjss6rXFQBBrp62uHs5yLOMBDZ7VKk82agq6d5l6iqxEBXTzuS/Wxw/FwVYaCrZ40nPJdLeJcduirCQFfPejaXcpw674/Xu12K1BEGunrW/lwOwPa+17pcidQZBrp61r7xFaxklEvjZLdLkTrCQFfP2j++nPf3vU5fdLsSqTMMdPWkN6jxvVzGdsfPVSEGunrSd/rfwRjBlY6fq0IMdPWkff0XAnBlnx26qsNAV0/aP3ARm+Mka2O026VIHWOgqyft67/Q8XNVjoGunvNi4w0O15Y53KLKMdDVc/Y/NwJ4Q5Gqx0BXz9l/cISBHGNbnOh2KVJHGejqOfsOjrDtzAhLnANdFWOgq6cMHz/Ftw+OcNWZl7tditRxBrp6yp0PHmB0PPnYiWe6XYrUcQa6esbBV05wz8M/4OeHNvGeMb8QVfUY6OoZv/f1p4kIPv3hrd0uRZoXBrp6wlNHjvO1fYf4+AffzTtXX9DtcqR5YaCr8jKT3/zzJ1k+UOdT117W7XKkeWOgq9Jeff00O7/8CF9/4gifuvZSLlw+0O2SpHnTVqBHxHUR8WREHIiI26bYviQi/lux/eGI2NLpQqW5yEz+zzMvccPnvsVfPnmUX/upK/jUj13a7bKkeVWfbYeIqAF3Aj8BHAL2RMSuzHy8ZbdPAq9m5mURcRPwWeAfzkfB0mSnRscYPn6KI8dOcejVE/z1My/zradf4vmRN7hk7XL+xz//EO/buLrbZUrzbtZAB64GDmTmswARcR+wA2gN9B3AbxSvvwr8x4iIzOz4rXj37znIF771bKffVgug3V+GzHxz32z+vfFMMmFsPBkbT0bHk1NnxnjjzBij429955VL6/zIpRfxqWsv5ae3b2T5knZ+zaXya+c3fSNwsGX5EPDD0+2TmaMR0QAuAl5q3SkidgI7AS6++OJzKnjNsn62rl9xTn9X3Re0+QDP4OyefRFEsVzr66PW1/zzgv4aS/v7WDZQY3DlEtatXMr6VUu5fP0K6jW/HlLvWdDWJTPvAu4CGBoaOqfu/aM/tIGP/tCGjtYlSVXQThvzPLC5ZXlTsW7KfSKiDqwGnCxDkhZQO4G+B9gaEZdExABwE7Br0j67gI8Xr38W+Iv5GD+XJE1v1iGXYkz8FuABoAZ8KTMfi4g7gL2ZuQv4A+DLEXEAeIVm6EuSFlBbY+iZuRvYPWnd7S2vTwI/19nSJElz4aUAklQRBrokVYSBLkkVYaBLUkVEt64ujIhh4Afn+NfXMuku1B7Ri8fdi8cMvXncvXjMMPfjfndmDk61oWuBfj4iYm9mDnW7joXWi8fdi8cMvXncvXjM0NnjdshFkirCQJekiihroN/V7QK6pBePuxePGXrzuHvxmKGDx13KMXRJ0tuVtUOXJE1ioEtSRZQu0Gd7YHUVRMTmiHgwIh6PiMci4jPF+gsj4n9FxNPFn+/odq2dFhG1iNgXEX9aLF9SPHj8QPEg8oFu19hpEbEmIr4aEd+LiCci4oM9cq7/VfH7/d2IuDcillbtfEfElyLiaER8t2XdlOc2mj5XHPujEXHVXD+vVIHe8sDq64FtwM0Rsa27Vc2LUeDWzNwGXAP8YnGctwHfyMytwDeK5ar5DPBEy/Jngd/NzMuAV2k+kLxq/gPw55n5t4AP0Dz+Sp/riNgIfBoYysz30Zyae+IB81U633cD101aN925vR7YWvzsBD4/1w8rVaDT8sDqzDwNTDywulIy88XM/H/F6+M0/wPfSPNY/7DY7Q+Bf9CdCudHRGwCfgr4YrEcwI/TfPA4VPOYVwN/h+YzBcjM05k5QsXPdaEOXFA85WwZ8CIVO9+Z+U2az4hoNd253QH8UTY9BKyJiHfO5fPKFuhTPbB6Y5dqWRARsQXYDjwMrM/MF4tNh4H1XSprvvwe8G+A8WL5ImAkM0eL5Sqe70uAYeC/FENNX4yI5VT8XGfm88BvA8/RDPIG8AjVP98w/bk973wrW6D3lIhYAfx34F9m5rHWbcUj/ipzzWlE/D3gaGY+0u1aFlgduAr4fGZuB15n0vBK1c41QDFuvIPm/9DeBSzn7UMTldfpc1u2QG/ngdWVEBH9NMP8nsz8WrH6yMQ/wYo/j3arvnnwIeDGiPg+zaG0H6c5trym+Cc5VPN8HwIOZebDxfJXaQZ8lc81wEeAv8nM4cw8A3yN5u9A1c83TH9uzzvfyhbo7TywuvSKseM/AJ7IzN9p2dT6MO6PA3+y0LXNl8z8lczclJlbaJ7Xv8jMfwQ8SPPB41CxYwbIzMPAwYh4b7Hqw8DjVPhcF54DromIZcXv+8RxV/p8F6Y7t7uAf1xc7XIN0GgZmmlPZpbqB7gBeAp4BvjVbtczT8f4ozT/GfYosL/4uYHmmPI3gKeBrwMXdrvWeTr+a4E/LV6/B/i/wAHgK8CSbtc3D8d7JbC3ON9/DLyjF8418G+B7wHfBb4MLKna+QbupfkdwRma/xr75HTnFgiaV/E9A3yH5hVAc/o8b/2XpIoo25CLJGkaBrokVYSBLkkVYaBLUkUY6JJUEQa6JFWEgS5JFfH/AU/+RCNGrQGcAAAAAElFTkSuQmCC\n", | |
| "text/plain": [ | |
| "<Figure size 432x288 with 1 Axes>" | |
| ] | |
| }, | |
| "metadata": { | |
| "needs_background": "light" | |
| }, | |
| "output_type": "display_data" | |
| } | |
| ], | |
| "source": [ | |
| "non_zero_idx = np.nonzero((dy > 0.05))[0]\n", | |
| "\n", | |
| "\n", | |
| "fig, ax = plt.subplots()\n", | |
| "ax.plot(x,y)\n", | |
| "\n", | |
| "\n", | |
| "\n", | |
| "print(non_zero_idx)\n", | |
| "ax.axvspan(non_zero_idx[0], non_zero_idx[-1], alpha=0.5, color='red')\n", | |
| "\n", | |
| "plt.show()" | |
| ] | |
| }, | |
| { | |
| "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