Skip to content

Instantly share code, notes, and snippets.

@drvinceknight
Created March 10, 2018 21:53
Show Gist options
  • Save drvinceknight/f3163de47741f5de5a3a452e243e0591 to your computer and use it in GitHub Desktop.
Save drvinceknight/f3163de47741f5de5a3a452e243e0591 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"%matplotlib inline\n",
"\n",
"time = np.array([0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0])\n",
"height = np.array([1.1, 2.2, 3.1, 4.0, 4.8, 5.5, 6.0, 6.4, 6.9, 7.1])"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"array([ 6.67272727, 1.04 ])"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"line = np.polyfit(x=time, y=height, deg=1)\n",
"line"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7fb030039978>]"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAD8CAYAAABXe05zAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAGbtJREFUeJzt3Xt0VfWd9/H3lxgkIBAUtCYQwKpx\nFMTAEW+tbbWK4/QRvNT7OLbOMNOZUds6ceR5qnas1kuso+20nTK19qZWqwhMQYJtneVoFQ0iEi5R\nqtwCCgLhGiAk3+ePc6IhnHB2yLnsfc7ntRbLZGcn5+te4eP2c357b3N3REQkOnrlegAREekeBbeI\nSMQouEVEIkbBLSISMQpuEZGIUXCLiESMgltEJGICBbeZfcPMFptZvZk9aWZ9Mj2YiIgklzK4zawc\nuAmIufsooAi4MtODiYhIcod0Y78SM2sB+gJrD7Tz4MGDfcSIET0cTUSkcMyfP/8jdx8SZN+Uwe3u\njWb2ILAKaAbmuvvcA33PiBEjqKurCzSsiIiAma0Mum+QqmQQMBEYCZQB/czs2iT7TTazOjOr27Bh\nQ3fmFRGRbgjy5uQXgffdfYO7twDTgDM77+TuU9095u6xIUMCne2LiMhBCBLcq4DTzayvmRlwLrA0\ns2OJiEhXUga3u88DngHeBBYlvmdqhucSEZEuBFpV4u53AndmeBYREQlAV06KiESMgltEJGIU3CIi\nB6luxSZ++OLyrL9u0CsnRUQkYcO23dz3/DKefXMN5aUlXHfGcPr3Kc7a6yu4RUQC2tvaxuPzVvHg\n3AZ2tbTytc9/mhvPOZa+vbMbpQpuEZEA6lZs4vYZi1m6biufPW4w377oJD495LCczKLgFhE5gI+2\nx2uRZ+av4eiBffjRNWP5y1GfIn49Ym4ouEVEkkhWi/zzF46l36G5j83cTyAiEjLzV27i9umLWbJu\nK585Nl6LHHtkbmqRZBTcIiIJYaxFklFwi0jB61yL/MPn4qtFwlCLJBPOqUREsiTstUgyCm4RKUgd\na5FPDejDD68ey4Wjw1eLJKPgFpGC0trmPD5vJQ/WNrBzT/hrkWSiM6mISA91rEXOOvYI/u2ikzj2\nyP65HqvbFNwikvc+2r6b+59fxm8jWIsko+AWkbzVuRb5+88dw03nHBepWiSZaE8vItKF+Ss3c8eM\nehav3cqZnz6CuyZGsxZJRsEtInmlcy3yH1dX8Vejj45sLZKMgltE8kK2a5HpCxqpqW1gbVMzZaUl\nVE+oZFJVeUZeqzMFt4hEXrZrkekLGpkybRHNLa0ANDY1M2XaIoCshHfK4DazSuCpDpuOAe5w94cz\nNpWISAAda5GjBhyatVqkprbh49Bu19zSSk1tQziC290bgFMAzKwIaASey/BcIiJdam1znpi3kpr2\nWuTsY7jx3OM4LEurRdY2NXdre7p199/yXODP7r4yE8OIiKTSsRY545h4LXLcUdldLVJWWkJjkpAu\nKy3Jyut39ynvVwJPZmIQEZED2bh9N7c+s5BLf/wn1mxuZlDfYl59byPXP/YG0xc0ZnWW6gmVlBQX\n7bOtpLiI6gmVWXn9wGfcZtYbuAiY0sXXJwOTASoqKtIynIhI51rknBOO5E/LP2LX3jYg+28Mdnyd\nXK0qMXcPtqPZROCf3P38VPvGYjGvq6vr6WwiUuDeXBWvReobP6lFrn/sjaQ1RXlpCa/cdk4OpkwP\nM5vv7rEg+3an474K1SQikgUbt+/m/jnLeLouvlrkB1dV8aWT46tFcv3GYBgECm4z6wecB/x9ZscR\nkULW2uY88foqauYs63K1SK7fGAyDQMHt7juAIzI8i4gUsGS1SLLVItUTKve5+AWy+8ZgGOjKSRHJ\nqY3bd/PAnAaeqlu9Xy2STK7fGAwDBbeI5ER7LfJgbQM7du/t1kU0k6rKCyqoO1Nwi0jWLVi1mdsD\n1CKSnIJbRLKmu7WIJKfgFpGMa21znnx9FTWJWmTy2cdwUxbvLZJvdNREJKNUi6SfgltEMmLTjj08\nMGcZv3kjXot8/6oq/o9qkbRQcItIWqkWyTwdSRFJm7dWN3H79HoWNW5RLZJBCm4R6ZZkz1o8+/gh\nPDBnGU/VrWbIYapFMk3BLSKBJXvWYvVvF1J8SC/27G3j7z6rWiQbdHRFJLBkz1psaXOszZl982c5\nXrVIVnT3CTgiUsC6unXqnr1tCu0sUnCLSCCtbc7AkuKkXysvoFuqhoGCW0RSemt1Exf/6BWamlvo\n1en9xkK7pWoYqOMWkS5t2rGHmtr4RTRDDjuUR648hbY258G57xTsLVXDQMEtIvtpbXN+80b8Ippt\nu/byt58ZyU3nHkf/PvGq5OKxQ3M8YWFTcIvIPhaubuL2GfW8vWYLpx9zOHdNHKU3HkNGwS0iQPJa\n5KIxZbqIJoQU3CIFrrXNeeqN1TxQu4xtu/Zyw1kjufmLn9QiEj4KbpEC1rEWOW1kvBap/JRqkbAL\nFNxmVgr8FBgFOPBVd381k4OJyP6S3SfkYFZ0bN6xhwdqG/jNG6tUi0RQ0DPuR4A57n6ZmfUG+mZw\nJhFJItl9QqZMWwQQOLxVi+SHlMFtZgOBs4HrAdx9D7Ans2OJSGfJ7hPS3NJKTW1DoOBeuLqJO2bU\ns1C1SOQFOeMeCWwAHjOzMcB84GZ339FxJzObDEwGqKioSPecIgWvq/uEdLW9XcdaZLBqkbwQ5JL3\nQ4CxwI/dvQrYAdzWeSd3n+ruMXePDRkyJM1jikhZF/cD6Wp7W+JJNF/43v/wdN1qbjhrJH+85XNM\nPKVcoR1xQYJ7DbDG3eclPn+GeJCLSBZVT6ikpLhon21d3SdkYeLeIlOmLeL4o/oz+6bP8q0vnagu\nO0+krErc/QMzW21mle7eAJwLLMn8aCLSUXuPfaBVJZt37KFmbgNPvq5aJJ8FXVVyI/B4YkXJe8BX\nMjeSiHRlUlV50jci29qcp+pWc/+c+GqRr541kq9rtUjeChTc7v4WEMvwLCJyEN5e08TtMxazcHUT\n40cezl0TT+KETw3I9ViSQbpyUiSiOtciD19xChNPUS1SCBTcIhHTXos8MGcZW1WLFCQFt0iE7FOL\njDicuyapFilECm6RCGjauYea2gaeeH0VR/Q7lIcuH8PFVVqPXagU3CIh1tbmPJ1YLbJ1116+cuZI\nvn7ecQxQLVLQFNwiIbVozRa+NaNetYjsR8EtElC6bqmaSuda5N+vGMMkXaYuHSi4RQJIxy1VU1Et\nIkEpuEUC6OktVVNZtGYLt8+o563VTZw6YhB3TRzFXxytWkSSU3CLBHCwt1RNRatF5GAouEUCKCst\noTFJSHd1S9VU2tqc385fzX3PL2NLcwvXnzmCb5x3vGoRCUTBLRJA9YTKfTpu6PqWqqmoFpGeUnCL\nBBDklqqpNO3cw4NzG3h8nmoR6RkFt0hAXd1SNZX2WuT+OQ007dyjWkR6TMEtkkH1jVv41vRPapF/\nu+g0TixTLSI9o+AWyYB9a5HeqkUkrRTcImnU1uY8M38N981ZRtPOPfzNGfFaZGCJahFJHwW3SJp0\nrEViwwdx10TVIpIZCm6RHtqys4UH5zbw63krOaJfb7735TFcMla1iGSOglvkIKkWkVxRcIschPrG\n+EU0C1apFpHsCxTcZrYC2Aa0AnvdXU98l4LUXos8Pm8lh6sWkRzpzhn3F9z9o4xNIhJibW3OM2+u\n4f7nl7F55x6uUy0iOaSqRCSFjrXIuOGD+OXE8ZxUNjDXY0kBCxrcDsw1Mwd+4u5TO+9gZpOByQAV\nFRXpm1AkR7bsbOF7LzTw69dWMqhvbx788hguqSqnVy/VIpJbQYP7M+7eaGZHAi+Y2TJ3f6njDokw\nnwoQi8U8zXOKZI1qEQm7QMHt7o2Jf643s+eA8cBLB/4ukeipb9zCHTPqeVO1iIRYyuA2s35AL3ff\nlvj4fOCujE8mkpCNh/SqFpEoCXLGfRTwXGK50yHAE+4+J6NTiSRk+iG9nWuRvz59ON88v1K1iIRa\nyuB29/eAMVmYRWQ/mXxIr2oRiSotB5RQy8RDejvXIjWXncylY4eqFpHIUHBLqKXzIb1tbc6zb67h\nPtUiEnEKbgm1dD2kt2MtMrailF/eoFpEokvBLaHW04f0bmlu4aG5DfxKtYjkEQW3hN7BPKQ3aS1y\nXiUD+6oWkehTcEveWbx2C3fMWMz8lZsZW1HKL746nlHlqkUkfyi4JW+oFpFCoeCWyOtci1x7+nBu\nUS0ieUzBLZGmWkQKkYJbIqlzLfLAZSdzmWoRKRAKbomUtjZn2oJG7nt+KZt2qBaRwqTglsjoWItU\nVZTy86+oFpHCpOCW0NvS3MK/v/AOv3x1BaWqRUQU3BJeHWuRjTv2cM1pFVSff4JqESl4Cm4JpSVr\nt3LHjHrqVIuI7EfBLaGyXy1y6clcNk61iEhHCm4JBXdn2puN3JuoRa49bTi3nH88pX1753o0kdBR\ncEvOqRYR6R4Ft+SMahGRg6PglqxTLSLSM4GD28yKgDqg0d2/lLmRJJ8tWbuVO2fW88aKzZwyrJTH\nrh/P6KGqRUS6oztn3DcDS4EBGZpFQmj6gsaDfvpMR6pFRNInUHCb2VDgr4B7gG9mdCIJjekLGvd5\n3mNjUzNTpi0CCBzenWuRa06r4F/Or1QtItIDQc+4HwZuBfpncBYJmZrahn0e0gvQ3NJKTW1DoOBW\nLSKSGSmD28y+BKx39/lm9vkD7DcZmAxQUVGRtgEld9Y2NXdre7utu1p4aO4ntcj9l47my+OGqRYR\nSZMgZ9xnAReZ2YVAH2CAmf3a3a/tuJO7TwWmAsRiMU/7pJJ1ZaUlNCYJ6bLSkqT7uzvPLWjku7OX\nsXHHbtUiIhnSK9UO7j7F3Ye6+wjgSuCPnUNb8lP1hEpKiov22VZSXET1hMr99l26biuX/+RVvvn0\nQoYOKmHmP32GuyeNVmiLZIDWcUuX2nvsA60q2bqrfbXISgb0OUS1iEgWmHv6W41YLOZ1dXVp/7kS\nHqpFRNLLzOa7eyzIvjrjlm5bui5+b5FPVoucqtUiIlmk4JbAVIuIhIOCW1JSLSISLgpuOaCl67Zy\n54zFvL5iE2OGlfKz62OcPLQ012OJFDQFtyTVuRa575LRXB5TLSISBgpu2Ye7M/2tRu6ZFa9Frh5f\nQfUE1SIiYaLglo8t+2Ard0xXLSISdgpuYeuuFh5+4V1+8eoK1SIiEaDgLmDttch3Zy/jo+3xWuRf\nzq9kUD/VIiJhpuAuUPvUIkMH8ujfqBYRiQoFd4HZtquFh3//Lj//k2oRkahScBcId2fGW2u5Z/ZS\nPtq+m6vGV1CtWkQkkhTcBaDhg23cPqOe199XLSKSDxTceaxzLXLvJaO5QrWISOQpuPOQahGR/Kbg\nzjOqRUTyn4I7T3SsRfqrFhHJawruiHN3Zi5cy92z4rXIladWcOsE1SIi+UzBHWENH2zjjhn1zHt/\nEycPHchPr4sxZphqEZF8p+AOsekLGpM+qHfbrhYe+f27PJaoRb578WiuOHUYRapFRAqCgjukpi9o\nZMq0RTS3tALQ2NTMbc++zfyVm6ld/AEbVIuIFKyUwW1mfYCXgEMT+z/j7ndmerBCV1Pb8HFot9u1\nt41fvbaSk4cO5L9Ui4gUrCBn3LuBc9x9u5kVAy+b2fPu/lqGZytoa5uau/zac/94lmoRkQLWK9UO\nHrc98Wlx4o9ndCrh6IF9km4vLy1RaIsUuJTBDWBmRWb2FrAeeMHd52V2rML2zofbKOldtN/2kuIi\nqidU5mAiEQmTQMHt7q3ufgowFBhvZqM672Nmk82szszqNmzYkO45C8L23Xu5Z9YSLnzkf9m4Yw+X\nx4ZSNrAPRvxM+95LRjOpqjzXY4pIjnVrVYm7N5nZi8AFQH2nr00FpgLEYjFVKd3QfhHNPbOWJlaL\nDKN6wgkcrtUiIpJEkFUlQ4CWRGiXAOcB92d8sgLxzofxi2heey9+Ec3U62KcotUiInIAQc64jwZ+\nYWZFxKuVp939d5kdK/9t372XR37/Do+9soLD+hzCPReP4spTK/TGo4iklDK43f1toCoLsxSE9lrk\nu7OXsn6bahER6T5dOZlFHWuR0eUD+c9rx1FVMSjXY4lIxCi4s6BjLdLvUNUiItIzCu4Mcnf+++11\n3DNrCR9ujdcit16gWkREekbBnSHvfLiNO2cs5tX3NqoWEZG0UnCnWeda5O5Jo7hqvGoREUkfBXea\nqBYRkWxRcKdBx1pkVPkAfnztOMaqFhGRDFFw98D23Xv5/h/e5Wcvv69aRESyRsF9EDrXIlfEhnHr\nBZUccdihuR5NRAqAgrub3v1wG3eoFhGRHFJwB9SxFunbu4jvTBrF1apFRCQHFNwpqBYRkbBRcB9A\nx1rkpDLVIiISDgruJFSLiEiYKbg76FyLXB4byr9ecIJqEREJFQV3Quda5EfXjGPccNUiIhI+BR/c\n23fv5Qd/eJdH22uRiSdx9WnD+e+Fa7npyQWsbWqmrLSE6gmVelCviIRCwQa3u/O7t9dxd5JaZPqC\nRqZMW0RzSysAjU3NTJm2CEDhLSI5V5DBvXx9vBb505+T1yI1tQ0fh3a75pZWamobFNwiknMFFdxd\n1SKdV4usbWpO+v1dbRcRyaaCCO72WuSeWUv5YOsuLo8N5dYLTmBwF6tFykpLaEwS0mWlJZkeVUQk\npV6pdjCzYWb2opktMbPFZnZzNgZLl+Xrt3HNT+dx45MLOLxfb5792pk8cNmYLkMboHpCJSXFRfts\nKykuonpCZabHFRFJKcgZ917gFnd/08z6A/PN7AV3X5Lh2Xpkx+69fP+P7/Lo/x64FkmmvceuqW3Q\nqhIRCZ2Uwe3u64B1iY+3mdlSoBwIZXC7O7MWrePu38VrkS+PG8q//mXXtUhXJlWVK6hFJJS61XGb\n2QigCpiXiWF6avn6bdw5czGvLN/IiUcP4IfXjNVFNCKSdwIHt5kdBjwLfN3dtyb5+mRgMkBFRUXa\nBgyicy1y18STuCZgLSIiEjWBgtvMiomH9uPuPi3ZPu4+FZgKEIvFPG0THoC7M3vRB3znd0t6VIuI\niERJyuA2MwMeBZa6+0OZHymY5eu3c+fMetUiIlJwgpxxnwX8NbDIzN5KbPu/7j47c2N1bcfuvfzg\nj8t59OX36FOsWkRECk+QVSUvAzlPxfZa5O5ZS1i3RbWIiBSuSFw5uXz9dr49czEvL/+IE48ewH9c\nXcW44YfneiwRkZwIdXCrFhER2V8og7tzLXLZuKHcplpERAQIaXBv2rGHW59ZyPAj+qkWERHpJJTB\nfcRhh/L0P5xB5VH9OaQo5X2wREQKSiiDG+CksoG5HkFEJJR0OisiEjEKbhGRiFFwi4hEjIJbRCRi\nFNwiIhGj4BYRiRgFt4hIxCi4RUQiRsEtIhIxCm4RkYhRcIuIRIyCW0QkYhTcIiIRE8q7A05f0EhN\nbQNrm5opKy2hekIlk6rKcz2WiEgohC64py9oZMq0RTS3tALQ2NTMlGmLABTeIiIEqErM7Gdmtt7M\n6rMxUE1tw8eh3a65pZWa2oZsvLyISOgF6bh/DlyQ4Tk+trapuVvbRUQKTcrgdveXgE1ZmAWAstKS\nbm0XESk0aVtVYmaTzazOzOo2bNhw0D+nekIlJcVF+2wrKS6iekJlT0cUEckLaQtud5/q7jF3jw0Z\nMuSgf86kqnLuvWQ05aUlGFBeWsK9l4zWG5MiIgmhW1UC8fBWUIuIJKcLcEREIibIcsAngVeBSjNb\nY2Y3ZH4sERHpSsqqxN2vysYgIiISjKoSEZGIUXCLiESMuXv6f6jZBmBl2n9wbgwGPsr1ECGg4/AJ\nHYs4HYe4dB2H4e4eaC11RoI7n5hZnbvHcj1Hruk4fELHIk7HIS4Xx0FViYhIxCi4RUQiRsGd2tRc\nDxASOg6f0LGI03GIy/pxUMctIhIxOuMWEYkYBXeCmV1gZg1mttzMbkvy9W+a2RIze9vM/mBmw3Mx\nZ6alOg4d9rvUzNzM8nJVQZDjYGaXJ34nFpvZE9meMVsC/N2oMLMXzWxB4u/HhbmYM5NSPQnM4r6f\nOEZvm9nYjA7k7gX/BygC/gwcA/QGFgIndtrnC0DfxMdfA57K9dy5OA6J/foDLwGvAbFcz52j34fj\ngAXAoMTnR+Z67hwei6nA1xIfnwisyPXcGTgOZwNjgfouvn4h8DxgwOnAvEzOozPuuPHAcnd/z933\nAL8BJnbcwd1fdPediU9fA4ZmecZsSHkcEr4D3A/syuZwWRTkOPwd8EN33wzg7uuzPGO2BDkWDgxI\nfDwQWJvF+bLCUz8JbCLwS497DSg1s6MzNY+CO64cWN3h8zWJbV25gfh/XfNNyuOQ+F/AYe4+K5uD\nZVmQ34fjgePN7BUze83MsvZc1iwLciy+DVxrZmuA2cCN2RktVLqbIT0SygcphJmZXQvEgM/lepZs\nM7NewEPA9TkeJQwOIV6XfJ74/329ZGaj3b0pp1PlxlXAz939e2Z2BvArMxvl7m25Hixf6Yw7rhEY\n1uHzoYlt+zCzLwL/D7jI3XdnabZsSnUc+gOjgP8xsxXEu7yZefgGZZDfhzXATHdvcff3gXeIB3m+\nCXIsbgCeBnD3V4E+xO/fUUgCZUi6KLjj3gCOM7ORZtYbuBKY2XEHM6sCfkI8tPO1zzzgcXD3Le4+\n2N1HuPsI4l3/Re5el5txMybl7wMwnfjZNmY2mHh18l42h8ySIMdiFXAugJn9BfHgPvgnhkfTTOC6\nxOqS04Et7r4uUy+mqgRw971m9s9ALfF30X/m7ovN7C6gzt1nAjXAYcBvzQxglbtflLOhMyDgcch7\nAY9DLXC+mS0BWoFqd9+Yu6kzI+CxuAX4LzP7BvE3Kq/3xFKLfJF4EtjngcGJLv9OoBjA3f+TeLd/\nIbAc2Al8JaPz5NnxFRHJe6pKREQiRsEtIhIxCm4RkYhRcIuIRIyCW0QkYhTcIiIRo+AWEYkYBbeI\nSMT8f+NSmQ2mxlQHAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7fb030070160>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.figure()\n",
"plt.scatter(time, height)\n",
"plt.plot(time, np.poly1d(line)(time))"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/home/vince/anaconda3/lib/python3.5/site-packages/ipykernel/__main__.py:1: RankWarning: Polyfit may be poorly conditioned\n",
" if __name__ == '__main__':\n"
]
},
{
"data": {
"text/plain": [
"array([ 5.35903504e+03, -2.20342165e+04, 3.22968022e+04,\n",
" -1.24946605e+04, -1.97891743e+04, 3.03477926e+04,\n",
" -1.92241427e+04, 6.74687107e+03, -1.34570302e+03,\n",
" 1.50051460e+02, -5.55531337e+00])"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"high_deg = np.polyfit(x=time, y=height, deg=len(time))\n",
"high_deg"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7fb0300111d0>]"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAD8CAYAAABXe05zAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAHsBJREFUeJzt3Xl4VdWh/vHvygAJEBKGMIVAkCEI\nAQmESWsFrGK1FmqtRUVAqTjU4dZeruW2v6u3tbe94rWTU6kgKAKiAtrWilZRqlRIQpglMiecAAlg\n5vmc9fsjkSIEciI5Z5/h/TxPHpKTzcnLfpI3i33WWttYaxERkeAR4XQAERFpGRW3iEiQUXGLiAQZ\nFbeISJBRcYuIBBkVt4hIkFFxi4gEGRW3iEiQUXGLiASZKF88adeuXW1KSoovnlpEJCRlZ2cft9Ym\nenOsT4o7JSWFrKwsXzy1iEhIMsYc8vZYXSoREQkyKm4RkSCj4hYRCTIqbhGRIKPiFhEJMipuEZEg\n45PpgCIioW5Njov5a3MpKK6iV0IscyenMjU9yS9fW8UtItJCa3JczFu1nao6NwCu4irmrdoO4Jfy\nbvZSiTEm1Riz5bS3UmPMv/k8mYhIAPF4LAeOV/DnrQX8bM2OU6X9hao6N/PX5volS7MjbmttLjAC\nwBgTCbiA1T7OJSLiGI/Hsv94BTsLSth+uITtrhJ2FZRSVlN/3r9XUFzll3wtvVRyJbDPWuv10kwR\nkUDm9lj2FZWzw9VQ0DtdpewsKKGitmFE3TYqgot7dmRqehJpSR1JS4rnziVZFJRUn/VcvRJi/ZK5\npcU9DVjuiyAiIr5W7/awt6ic7YdL2OEqYUdBKbsKSk9d9oiNjmRIr47cOKo3aUnxDOsdz4DEDkRF\nfvmq8n9cM/hL17i/+LtzJ6f65d/hdXEbY9oA3wbmnePzc4A5AH369GmVcCIiTfFmRkdtvYc9hWUN\nBe0qZburhE+PlFJT7wGgfZtIhvaKZ9qYZIYlxZOWFE//xA5ERphmv/4XX8upWSXGWuvdgcZMAX5o\nrb26uWMzMjKsdgcUEV84c0YHQExUBD+cOIAuHdqyo6BhNL37SBm17oaS7tA2iqG9Op4q6LSkePp1\nbe9VSfuLMSbbWpvhzbEtuVRyM7pMIiIOm78296wZHdX1Hv7v3c8A6BgTRVpSPLMuS2m43JEUT9/O\n7YgIoJK+UF4VtzGmPXAVcJdv44iInJvbY3GdZ+bG+rkTSe4cizGhU9JN8aq4rbUVQBcfZxERaZLb\nY/nLtgL+8P7ecx6TlBBLny7t/JjKOdqrREQCVr3bw+qcw1z1mw95cMUWIgzMGN+XmKgvV5c/Z3QE\nAi15F5GAU+/2sGZLAU+v28uB4xUM7hHHM7eO5JqhPYiIMIzs08mxGR2BQMUtIgGjzu1hdY6Lp9ft\n5dCJSi7u2ZHnpo/k6iE9vvTi4tT0pLAq6jOpuEXEcbX1HlZtPszTH+wl/2QVaUkdWXDbKK4a0j3k\nX2j8KlTcIuKYmno3r2Uf5pl1+3AVVzG8dzyPXj+USYO7qbDPQ8UtIn5XU+9mZWY+z36wj4KSakYk\nJ/DYd9KYMChRhe0FFbeI+E11nZtXGgv7aGk1o/p24tffHc7lA7uqsFtAxS0iPldd52bZxjye+3Af\nhWU1jE7pxBPfu4TLBnRRYX8FKm4R8ZnK2vrGwt7P8fIaxvbrzG+njWD8RSrsC6HiFpFWV1FTz9JP\nDvGnf+zneHktl/bvwlO3pDPuIi3Abg0qbhFpNeU19bz4z4M8/48DnKyo5fKBXXngyoGMTunsdLSQ\nouIWkQtWVl3Hkg0Hef6jAxRX1nHFoEQeuHIgo/p2cjpaSFJxi0iLnH4Tgx4dYxjRJ4EN+05QUlXH\nxNSGwk7vo8L2JRW3iHjtzJsYHCmt5siOowzt1ZGXZo9heO8EhxOGB+0OKCJesdby2F93nXUTA4Di\nyjqVth9pxC0i51VcWcvqHBcrNuVzvLy2yWMKznNzA2l9Km4ROYu1lk0HTrIiM5+/bj9Cbb2H4b3j\nSYiNpriq7qzjeyXEOpAyfKm4ReSUE+U1vL75MCsy89lfVEFc2yi+n5HMtDHJDO0V3+SNesPtJgaB\nQMUtEuY8HsuGfSdYnpnHOzuPUue2jOrbifk39ue64T1p1+ZfNfHFHtjhfBODQKDiFglThaXVvJp9\nmFcy88k7WUl8bDTTx/Xl5jF9GNQ97px/L9xvYhAIVNwiYcTtsazfU8TyjXm8t7sQt8cytl9nfnz1\nICYP7UFMdKTTEcULKm6RMFBQXMXKrHxezTqMq7iKLu3b8IOv9eP7o5O5KLGD0/GkhVTcIiGq3u3h\n/d2FrMjM54PcQjwWLh/Ylf+89mKuGtKdNlFaxhGsVNwiISb/ZCWvZOazMiufwrIausW15d4JA/j+\n6GSSO7dzOp60Aq+K2xiTADwPpAEWuMNa+09fBhMR79XWe/j7p8dYvimPj/YexwATUrsxbXQykwZ3\nIypSo+tQ4u2I+3fA29baG40xbQD92hZxwOkbPPVKiGXmpX05UV7La9mHOVFRS6/4GB68ciA3ZSRr\nUUwIa7a4jTHxwNeBWQDW2lqg6XWvIuIzZy5+cRVX8T9v7SbCwDcu7s7NY/vw9YGJREbozjKhzpsR\ndz+gCHjBGHMJkA08aK2t8GkyEfmS/317d5MbPHWLi2HBjAwHEolTvLnwFQWMBJ611qYDFcBPzjzI\nGDPHGJNljMkqKipq5Zgi4ausuo6n3t/DkZLqJj9/rLTpxyV0eTPiPgwcttZubPz4NZoobmvtAmAB\nQEZGhm21hCJhqqSqjsUfH2ThR/spra4nJiqC6nrPWcfpWnb4aba4rbVHjTH5xphUa20ucCWwy/fR\nRMJTSWUdCz8+wAsfH6Csup6rhnTngUkD2VdUrg2eBPB+Vsn9wMuNM0r2A7f7LpJIePq8opaFHx1g\n8YaDlNfUM3lodx64ciBDe8UDMKx3w5/a4Em8Km5r7RZAr36I+MDJilqe/8d+lmw4SEWtm2uH9eD+\nSQO5uGfHs47VBk8CWjkp4pjj5TX86R/7eemfh6iqc3PdsJ7cP2kgqT3OvTOfCKi4RfyuqKyGBev3\nsfSTPGrq3Vx/SS/umziAgefZSlXkdCpuET8pLK3muQ/3s2zTIWrrPUwdkcQPJw2gv3bnkxZScYv4\n2NGSap77cB/LNuXh9li+k57EDycOoF/X9k5HkyCl4hbxkYLiKp77cB8rNuXjsZYbRjYUdt8uKmy5\nMCpukVbmKq7imXV7eTXrMB5r+V5Gb+6dMEBbqkqrUXGLtJL8k5U888FeXss+DMBNGcncM6E/vTup\nsKV1qbhFLlDeiUqeWreHVZtdRBjDtNF9uGdCfy1FF59RcYt46cy9sGddmkLusTJW57iIjDBMH9eX\nu6/oT4/4GKejSohTcYt4oam9sH/51qdERxpmjk/h7isuoltHFbb4h4pbxAvz1+Y2uRd2l/Zt+a/r\nhziQSMKZbkQn4gVXcVWTj2svbHGCilvkPEqr63jolS3n/LxegBQnqLhFzmHD3uNc85v1vLG1gMlD\nuxMT9eUfF+2FLU7RNW6RM1TXuXn87VwWfXyAi7q25/V7LmVEcsJZs0q0F7Y4RcUtcpodrhJ+9MoW\n9hSWM2N8X+Z982Ji20QC2gtbAoeKWwSod3t47sN9/Pbve+jcvg1L7hjDFYMSnY4l0iQVt4S9g8cr\neGjlFjbnFfOt4T15bGoaCe3aOB1L5JxU3BK2rLUs25THY39pWEjzu2kjmDJCl0Ik8Km4JSwVllbz\n8OvbWJdbxOUDu/L4jcPpGa+pfRIcVNwSdt7afoSfrt5OZa2bR68fwozxKUREGKdjiXhNxS1ho7S6\njkff2MmqHBfDe8fz5E0jGNBNtw2T4KPilrCwYe9x/v3VrRwrq+GBKwdy/6QBREdq/ZkEJxW3hLRz\nLaYRCWZeFbcx5iBQBriBemtthi9DibSG8y2mEQlmLRlxT7TWHvdZEpFWosU0Eup0qURCysHjFfxo\n5RZytJhGQpi3xW2Bd4wxFvijtXaBDzOJtJgW00g48ba4v2atdRljugHvGmN2W2vXn36AMWYOMAeg\nT58+rRxT5Ny0mEbCjVfFba11Nf5ZaIxZDYwB1p9xzAJgAUBGRoZt5ZwiTdJiGglHzRa3MaY9EGGt\nLWt8/2rg5z5PJnIeJVV1PPrmTlZrMY2EIW9G3N2B1caYL45fZq1926epRE5z5g0MpozoxZoclxbT\nSNhqtrittfuBS/yQReQsa3JczFu1/dQd1l3FVTzzwT4S49pqMY2ELQ1TJKDNX5t7qrRPFxVhVNoS\ntlTcEtAKiquafPxoSbWfk4gEDhW3BKx6t4f2bZu+mtcrQdP9JHypuCUgnSiv4baFmyivqSfyjOl9\nsdGRzJ2c6lAyEedpybsEnK35xdyzNJsTFbU88b1LiIowX5pVMndyqu62LmFNxS0BZWVmPj97YweJ\nHRpmjaQlxQOoqEVOo+KWgFBT7+a//7yLZRvzuGxAF/5w80g6t9fmUCJNUXGL446WVHPPy9nk5BVz\n1xUXMffqVKK0oEbknFTc4qhNB05y78ubqayt5+lbRnLd8J5ORxIJeCpucYS1liUbDvLYXz8luXM7\nlt05lkHd45yOJRIUVNzid1W1bn66ejurclx84+JuPPn9EXSMiXY6lkjQUHGLX+WfrOSul7L59Ggp\nD101iPsmDtA2rCItpOIWv1n/WREPrMjB7bEsmjmaiYO7OR1JJCipuMXnrLU8++E+nliby8Bucfzx\ntlGkdG3vdCyRoKXiFp8qr6nn31du5e2dR/nW8J48fuNw2rXRt53IhdBPkPjMvqJy5ryYxcETlfzs\nuouZ/bV+NN6QQ0QugIpbfGLtzqP8eOVW2kZF8NLsMVzav6vTkURChopbWpXbY/nNu5/x1Lq9DO8d\nz3PTR2kLVpFWpuKWVlNcWcuDK7bw4WdF3JTRm59PSSMmOtLpWCIhR8UtrWJXQSl3L83mSEkVv/xO\nGreM6aPr2SI+ouKWC/bGFhcPv76N+NhoXrlrPCP7dHI6kkhIU3HLV1bn9vCrt3az6OMDjEnpzFO3\nptMtLsbpWCIhT8UtX0lRWQ33LdvMxgMnmXVpCj+97mKitRWriF+ouKXFtuQXc/dL2RRX1fKb71/C\nd9J7Ox1JJKx4XdzGmEggC3BZa7/lu0gSSNbkuL50v8evDejK6hwX3To23FpsaK94pyOKhJ2WjLgf\nBD4FOvooiwSYNTku5q3aTlWdGwBXcRWvZOWT2iOOFXeOo5NuLSbiCK8uShpjegPXAc/7No4Ekvlr\nc0+V9unKqupU2iIO8vbVpN8C/wF4znWAMWaOMSbLGJNVVFTUKuHEWQXFVU0+fqSk2s9JROR0zRa3\nMeZbQKG1Nvt8x1lrF1hrM6y1GYmJia0WUJzTpUPTo2otYRdxljcj7suAbxtjDgIrgEnGmKU+TSWO\n++u2I3xeWceZax9joyOZOznVkUwi0qDZ4rbWzrPW9rbWpgDTgPettdN9nkwcs+ijA9y3fDPpyQk8\nNjWNpIRYDJCUEMuvbhjG1PQkpyOKhDXN45ZTPB7Lr9/ezYL1+5k8tDu/m5ZOTHQkt47r63Q0ETlN\ni4rbWvsB8IFPkoijaurdzH11G29uLWDG+L48cv1QInUTX5GApBG3UFpdx90vZbNh3wkevmYwd19x\nkXb2EwlgKu4wd6y0mpmLNrG3sJwnb7qEG0Zq+bpIoFNxh7E9x8qY9UImxZW1vHD7aC4fqGmcIsFA\nxR2mMg+eZPbiTNpGR/LKXeNJS9KeIyLBQsUdht7ecYQHVmyhd6dYltw+huTO7ZyOJCItoOIOM0s2\nHOTRP+8kPTmBhTNHa88RkSCk4g4THo/l8bW5PPfhPq4a0p3fT0snto1u5CsSjFTcYaC23sPDr29j\ndY6LW8f24edT0jRHWySIqbhDXFl1Hfcs3cxHe48zd3Iq907orznaIkFOxR3CCkurmfVCJrnHyph/\n43C+l5HsdCQRaQUq7hC1t7CcmYs28XllLQtnZjAhtZvTkUSklai4Q1D2oZPMXpJFVIThlTnjGdZb\nc7RFQomKO8Ss3XmUB5bn0CuhYY52ny6aoy0SalTcIeSlTw7xyBs7GN47gYUzM+jSoa3TkUTEB1Tc\nIcBayxPv5PL0un1cObgbT90yUnO0RUKYijvI1bkb5miv2uzi5jHJ/GJKGlGR3t4DWkSCkYo7iJXX\n1HPvy5tZ/1kRD101iPsnDdAcbZEwoOIOUoVl1dyxOJNPj5Tx+HeHc9NozdEWCRcq7iC0v6icGYs2\ncaK8ludnZDBxsOZoi4QTFXeQ2Zz3ObMXZxJhDCvmjOOS5ASnI4mIn6m4g8i7u45x//LNdO8Yw5Lb\nx5DStb3TkUTEASruALYmx8X8tbkUFFcRHxtNSXUdw5PiWThrNF01R1skbKm4A9SaHBfzVm2nqs4N\nQHFVHREGbhnbR6UtEuaanfBrjIkxxmwyxmw1xuw0xvy3P4KFu/lrc0+V9hc8Fn7/3l6HEolIoPBm\nxF0DTLLWlhtjooGPjDF/s9Z+4uNsYc1VXNXk4wXneFxEwkezI27boLzxw+jGN+vTVGHueHkN0ZFN\nL6TplRDr5zQiEmi8WhttjIk0xmwBCoF3rbUbfRsrfB08XsF3n90AQJszlq7HRkcyd3KqE7FEJIB4\nVdzWWre1dgTQGxhjjEk78xhjzBxjTJYxJquoqKi1c4aFLfnF3PDsBsqq61l513gev3E4SQmxGCAp\nIZZf3TCMqelJTscUEYcZa1t21cMY819ApbX2iXMdk5GRYbOysi40W1h5f/cxfvhyDolxbVlyxxj6\naY62SFgxxmRbazO8OdabWSWJxpiExvdjgauA3RcWUU63YlMed76YzYBuHXj9nktV2iJyXt7MKukJ\nLDHGRNJQ9CuttX/xbazwYK3lt3/fw+/e28MVgxJ55taRtG+rqfUicn7NtoS1dhuQ7ocsYaXe7eGn\nq3fwSlY+N47qza9uGEa09tEWES9oeOeAytp67luWw/u7C7l/0gAeumqQ9tEWEa+puP3sRHkNdyzO\nZLurhMempjF9XF+nI4lIkFFx+9GhExXMXLSJIyXVPDd9FFcP7eF0JBEJQipuP9maX8wdizPxWMuy\nO8cxqm8npyOJSJBScfvButxC7l26mS4d2rDkjjH0T+zgdCQRCWIqbh9bmZXPvFXbGdwjjhduH023\nuBinI4lIkFNx+4i1lj+8v5cn3/2Mywd25dnpo+igOdoi0grUJD5Q7/bw/97YyfJNedyQnsSvvzuc\nNlGaoy0irUPF3cqqat3cv3wzf/+0kHsn9Gfu5FTN0RaRVqXibkUnK2q5Y3EmWw8X8/MpQ5kxPsXp\nSCISglTcrSTvRCUzX9hEQXEVz946imvSNEdbRHxDxd0KdrhKmPVCJnVuDy//YCwZKZ2djiQiIUzF\nfYE+/KyIe5dmk9CuDSvmjGVAtzinI4lIiFNxX4DXsg/zk9e3MbB7HItvH033jpqjLSK+p+L+Cqy1\nPPPBPuavzeWyAV14dvooOsZEOx1LRMKEiruF3B7LI2/uYOkneUwZ0Yv5N16iOdoi4lcq7haornPz\nwPIc3tl1jLu+fhEPXzOYiAjN0RYR/1Jxe+nzilp+8GIWm/M+55Hrh3D7Zf2cjiQiYUrFfQ5rclzM\nX5tLQXEV3eLaYoHiqjqevmUk1w7r6XQ8EQljKu4mrMlxMW/Vdqrq3AAcK6sB4P6JA1TaIuI4varW\nhPlrc0+V9ulW5bgcSCMi8mUq7iYUFFe16HEREX9ScTchMa5tk4/3Soj1cxIRkbOpuM9w6EQFVXVu\nzpzkFxsdydzJqY5kEhE5XbPFbYxJNsasM8bsMsbsNMY86I9gTigsrea2hZuIijA8/M3BJCXEYoCk\nhFh+dcMwpqYnOR1RRMSrWSX1wI+ttZuNMXFAtjHmXWvtLh9n86uSqjpmLNrE8fIalt05jhHJCdx9\nRX+nY4mInKXZEbe19oi1dnPj+2XAp0BIDT2rat38YEkm+4rK+eNtoxiRnOB0JBGRc2rRPG5jTAqQ\nDmz0RRgn1Lk93LdsM1mHPuepm0dy+cBEpyOJiJyX1y9OGmM6AK8D/2atLW3i83OMMVnGmKyioqLW\nzOgzHo/l4de28d7uQn4xJY3rhmtxjYgEPq+K2xgTTUNpv2ytXdXUMdbaBdbaDGttRmJi4I9arbX8\n8q1PWZXj4qGrBjF9XF+nI4mIeMWbWSUGWAh8aq190veR/OOZD/ax8KMDzLo0hfsnDXA6joiI17wZ\ncV8G3AZMMsZsaXy71se5fGr5pjzmr81l6ohe/Ne3htDwu0lEJDg0++KktfYjOGs9StD62/Yj/HT1\ndiakJjL/e5doP20RCTphtXJyw97jPLhiC+l9OvHsraOIjgyrf76IhIiwaa5th4u588Us+nVtz6KZ\no4ltE+l0JBGRryQsintvYTmzXsikU/s2vDh7DPHtdGNfEQleIV/cR0qqmLFwIxEGXpo9lu4dY5yO\nJCJyQUK6uD+vqOW2hZsoq65n8e1j6Ne1vdORREQuWMjeuqyipp5ZizPJO1nJi3eMIS0p3ulIIiKt\nIiRH3LX1Hu5ems32w8U8dXM64y7q4nQkEZFWE3IjbrfH8tDKLfxjz3Eev3E4Vw/t4XQkEZFWFVIj\nbmstj7y5g79sO8K8bw7mpoxkpyOJiLS6kCru3/x9D0s/yeOuKy7iLt0EQURCVMgU9+KPD/D79/Zw\nU0ZvfnLNYKfjiIj4TEgU9xtbXDz6511cPaQ7//OdYdo0SkRCWtAX97rcQn68citj+3Xm9zenE6X9\nR0QkxAV1y2UfOsk9S7NJ7RHH8zMziInW/iMiEvqCtrhzj5Zx+wuZ9OgYw+LbxxAXo/1HRCQ8BGVx\n55+sZMaijcRER/LS7LEkxrV1OpKIiN8EXXEXldVw28KNVNW6eWn2WJI7t3M6koiIXwXVysnS6jpm\nvbCJo6XVvPyDsaT2iHM6koiI3wXNiLu6zs2dS7LIPVrGs9NHMapvZ6cjiYg4IihG3PVuDw8sz2Hj\ngZP8btoIJqZ2czqSiIhjAn7Eba3lP1dv551dx3j0+iFMGZHkdCQREUcF5Ih7TY6L+WtzKSiuon3b\nKMpr6nlg0gBmXdbP6WgiIo4LuOJek+Ni3qrtVNW5ASivqScywujuNSIijQLuUsn8tbmnSvsLbo/l\niXc+cyiRiEhgaba4jTGLjDGFxpgd/ghUUFzVosdFRMKNNyPuxcA1Ps5xSq+E2BY9LiISbpotbmvt\neuCkH7IAMHdyKrFnbBYVGx3J3Mmp/oogIhLQAu7FyanpDdP9vphV0ishlrmTU089LiIS7lqtuI0x\nc4A5AH369Lmg55qanqSiFhE5h1abVWKtXWCtzbDWZiQmJrbW04qIyBkCbjqgiIicnzfTAZcD/wRS\njTGHjTGzfR9LRETOpdlr3Nbam/0RREREvKNLJSIiQUbFLSISZIy1tvWf1Jgi4FCrP7EzugLHnQ4R\nAHQe/kXnooHOQ4PWOg99rbVeTcnzSXGHEmNMlrU2w+kcTtN5+BediwY6Dw2cOA+6VCIiEmRU3CIi\nQUbF3bwFTgcIEDoP/6Jz0UDnoYHfz4OucYuIBBmNuEVEgoyKu5Ex5hpjTK4xZq8x5idNfP4hY8wu\nY8w2Y8x7xpi+TuT0tebOw2nHfdcYY40xITmrwJvzYIy5qfF7YqcxZpm/M/qLFz8bfYwx64wxOY0/\nH9c6kdOXmrsTmGnw+8ZztM0YM9Kngay1Yf8GRAL7gIuANsBWYMgZx0wE2jW+fw/witO5nTgPjcfF\nAeuBT4AMp3M79P0wEMgBOjV+3M3p3A6eiwXAPY3vDwEOOp3bB+fh68BIYMc5Pn8t8DfAAOOAjb7M\noxF3gzHAXmvtfmttLbACmHL6AdbaddbaysYPPwF6+zmjPzR7Hhr9AvhfoNqf4fzIm/NwJ/C0tfZz\nAGttoZ8z+os358ICHRvfjwcK/JjPL2zzdwKbArxoG3wCJBhjevoqj4q7QRKQf9rHhxsfO5fZNPx2\nDTXNnofG/wImW2v/6s9gfubN98MgYJAx5mNjzCfGGL/dl9XPvDkXjwLTjTGHgbeA+/0TLaC0tEMu\nSMDduizQGWOmAxnAFU5n8TdjTATwJDDL4SiBIIqGyyUTaPjf13pjzDBrbbGjqZxxM7DYWvt/xpjx\nwEvGmDRrrcfpYKFKI+4GLiD5tI97Nz72JcaYbwA/Bb5tra3xUzZ/au48xAFpwAfGmIM0XMt7MwRf\noPTm++Ew8Ka1ts5aewD4jIYiDzXenIvZwEoAa+0/gRga9u8IJ151SGtRcTfIBAYaY/oZY9oA04A3\nTz/AGJMO/JGG0g7V65nnPQ/W2hJrbVdrbYq1NoWGa/3fttZmORPXZ5r9fgDW0DDaxhjTlYZLJ/v9\nGdJPvDkXecCVAMaYi2ko7iK/pnTem8CMxtkl44ASa+0RX30xXSoBrLX1xpj7gLU0vIq+yFq70xjz\ncyDLWvsmMB/oALxqjAHIs9Z+27HQPuDleQh5Xp6HtcDVxphdgBuYa6094Vxq3/DyXPwY+JMx5kc0\nvFA5yzZOtQgVjXcCmwB0bbyW/wgQDWCtfY6Ga/vXAnuBSuB2n+YJsfMrIhLydKlERCTIqLhFRIKM\niltEJMiouEVEgoyKW0QkyKi4RUSCjIpbRCTIqLhFRILM/wfKLHA0NwYSaAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7fb02ffe99b0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.figure()\n",
"plt.scatter(time, height)\n",
"plt.plot(time, np.poly1d(high_deg)(time))"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"anaconda-cloud": {},
"kernelspec": {
"display_name": "Python [default]",
"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.5.2"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment