Skip to content

Instantly share code, notes, and snippets.

@ulikoehler
Created February 7, 2016 22:25
Show Gist options
  • Save ulikoehler/6198f9f154a36e06f2c1 to your computer and use it in GitHub Desktop.
Save ulikoehler/6198f9f154a36e06f2c1 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import numpy as np\n",
"from scipy.interpolate import splrep, splev\n",
"import matplotlib.pyplot as plt\n",
"plt.style.use(\"ggplot\")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Create a simple waveform to work with"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7fad0ed08198>]"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAEECAYAAADK0VhyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXlwXdd93z/nAQRJEAAXkNg3AuBO7TJlSZZNW5ZjO07k\ntOmJE09rd0nUcdx0maZNHHek6WQ6SZc0dZM0ntbx2K1b5zqOd1uxZYe2FkvWRklcxA3E9rByBwmS\nAPFO/zjvQSCI/d3t3Pv7zGgEPN537g/3vne+97ec31HGGARBEIR0k4naAEEQBCF6RAwEQRAEEQNB\nEARBxEAQBEFAxEAQBEFAxEAQBEEASv0YRGv9OeBDwLDnebfPc8xngA8AV4CPe5530I9zC4IgCMXj\nl2fweeDn5vtHrfUHgA7P87YBjwF/vtSBtdb7i7ZOAORa+o1cT3+R6+kfK7mWvoiB53nPAOcXOORR\n4Iv5Y18A1muta5c4/P7irBNmsD9qAxLG/qgNSBj7ozYgQexf7hvCyhk0An0zfs/mXxMEQRBigCSQ\nBUEQBH8SyEsgCzTP+L0p/9ot5GNd+wu/e573OPB4gLalBs/zQK6lb8j19Be5nv7heR5a65kvHfA8\n78BC7/FTDFT+v7n4JvCbwF9qrd8OXPA8b3iuA/MGH5jx0uMDAwM+mglmaorcf/wd1M47UB/+KEpZ\ns83kBLn/8mnUrjvJPPprvp4zDlRWVjI2Nha1GUsm96U/x1w6T+bXfxtVaj+qxhjM//kzzOUxMv/0\n307fuyhw7XoGQe4nT2Ke+haZ3/1PqLXlb73+t9/B/O13yfzeH6FWr17SWHI9/aOhoQHP855Yznt8\nCRNprf8v8BywXWvdq7X+h1rrx7TWvwHged53gdNa65PAZ4FP+HHelWL+5q+hbDXq0V+7aTJRq8rI\nfOJ3Mc/9EHP0tQgtFMwrz2EOv0LmY781LQQASinUR34Dzo5gfvitCC0UTG8X5utfIvOJT90kBABq\n/wdRzVsxX/18RNYJy0U50MLa+OkZmOEBcn/wb8h8+r+iqrfMfczLz5L77lfsU00mOWkVV568zOQE\nuU89Zp/8O3bOfczwALk/+G0yv//nqHWVIVtoceV6BsXUHz+OuuM+Mu/+4Jz/bsYvk3vit8g89m/m\nvY8zSfv19JOGhgaYP1IzJ8mZ6ZaI+d5fod798/MKAQB3PwAlpZgXnw7PMGEa8/T3obVjwQlE1Tag\n7rof8zdfC9EyoYA59gYMD6AeemTeY1R5BeqDv0zuO16IlgkrJVViYM6OYl59HvXwhxY8TilF5pc/\njvnGlzC5qZCsE8B6BeZ7XyXzCx9Z9Fj187+C+fGTmEsXQrBMmEnua/8b9ehHUaWrFjxOPfhe6DuN\n6TkZkmXCSkmXGPzg66h3PLKksILavhfKK+CwdM0IE/PTH0FLO6q1c9FjVfUW1L0PYn78ZAiWCQVM\nz0m4cA6176FFj1WrylA/90vkvvuVECwTiiE1YmCuX8f89Eeo9/7ikt+j3vV+cj/+XoBWCbMxT/+A\nzP65Y9Bzod75fsyzT2FyueCMEm7C/ORvUA+9D5UpWdLx6h2PwJuvYy4u1KRAiJr0iMGrP4X2naiN\n1Ut+j9r3Tjh5FHNuNEDLhAIm2wMXzsGeO5f8HtXaAesq4M3XA7RMKGCuXcW89AzqwYeX/B61Zi3q\nzrdjXjgQnGFC0aRHDH76t6j7372s96jVa1D7HsI881RAVgkzMc8+hXrgPUt+4iyg3vEI5pkfBGSV\nMBPz4tOwfS9qw9IfqgDUAw9jnvsRDlQvppZUiIE5fxa6T6DuvG/Z71X37ce89EwAVgkzMTduYJ4/\nsKwnzgJq37swh17GjF8OwDJhJuaFH5N5YPn3iG274dpV6O3y3yjBF9IhBj/7Cequt6PKlrYS8ia2\nboer45jBvsWPFVbOicNQXYOqaVj2W9W6Cti2B/PGywEYJhQwYxeh9xTsuWvZ71WZDOr+d0uoKMak\nQwwOPo+698EVvVdlMqi778e8/JzPVgkzMQdfWJHnVkDd9XZ49XkfLRJmYw6+gNp918oeqsDmDQ6+\nIKGimJJ4MTBjlyDbAztuW/EY6p4HMK+IGASFMcYK9l1vX/EY6o59mCMHMZMTPlomzMS8+jzcff/K\nB2hphxs3YKjfP6OEmzDnRldcWZd8MXjjJdh5O2pV2coH6dwFF89jRof8M0x4i74uKF0F9c2LHzsP\nqnI9NLfBEekpFQRm/AqcOIy67d4Vj6GUsqJ98Gc+WibMJPeff88+/K6A5IvBaz9D3bGvqDFUpgR1\n2z2YQxKTDgLzqg0RFduBVN31dltCLPjPkVehc9ctDemWi7rzPsxrL/hklDATMzIIE9ehqW1F70+0\nGJjJSTj6WlFPM9Psvgtz+NXixxFuwbzxEur24gQbQO2914aKJCbtO+bIQdQKEse3sH0vDPZJC5EA\nMIdfRe2+c8UPVYkWA04egbpGVNWGoodSu+6E44cwNyZ9MEwoYC5fguEsdOwofrDaBtuncWjOfZOE\nFWKMsWKwq3gxUKtWwfa90iI+AMzhV2DP3St+f6LFwLz5BmrXHb6MpSqroKYBuo75Mp6Q59gb0Ll7\n0YZnS0Ephdp1J+ao9JPylZFBmLoBDSvP6cxE7bxDVoz7jLkxCccPoXYvffX+bJItBsdeR+283bfx\n1O47MdK4zlfM0dd8E2wAdt2BOSL3yE+sV7Dy8MNs1M7bMSIG/nLqGNQ02EKKFZJYMTDXxqG/G5aw\nqcZSUXvutq6Y4Bt+i4HadQccP4yZktbjfmGOHIQinjhvoaEZJiekOs9HzJGDRXkFkGAx4MQRaNu2\n4gUyc9K+A4b6MVfH/RszxZizI3B1HBpbfRtTVW2A6ho4fdy3MdOMyU3B8TdQu30UbKVQO8Q78BNz\n4hCqiLVUkGAxMG++jtpZ3MWZjVq1Clo74dRRX8dNK/Ye3e771qJq522Y44d8HTO19PdA1QZU1UZ/\nx915G0gS2RfM5ITt+VRkEUayxWCHf/mCAmrbbsyJI76Pm0pOHIZte3wfVm3bjTkpgu0H5sQRVBD3\naOftmOOHpAzYD06fgPpm1Jri1oAkUgzMtXEYHoC2bb6PrbbtwRw/7Pu4acScfBPVucv/gTt3wamj\nsuGNH5w8Ap27/R93cy2g4Myw/2OnDHPiMGp78YKdSDGg6zg0t9uwjt907IC+LumBUyRm7CJcugCN\nLb6Prao2QsV6GOj1few0YYzJewb+i4FSCjp2YiTkWjTm+GFfvLdEioE59SbKxyqimag15VDXZF0z\nYeWcPArt25e9kc1SUdt2YU5KOK8oRofsIr7NtYEMrzp32c+BsGLM1BScPma94SJJqBgcRXUGIwaQ\nDxWdkFBRMZiTR1FBhB8KdOyCEzLRFIM5eQTVudu39QWzUZ27MKfeDGTs1NDfDRs3oyqqih4qcWJg\nclM2TNQeoBh07sLISuSisIIdQL4gj00ii2dQFCePBpMvKNDcDqPDtiOqsCJM1zFUuw+tXEigGDDQ\nB5XrfelHNC9bt8Hp41IJsULM5AT0nba7yAVFbSNMXLNbngorwnQdQ/nRM2oeVGkptHZIi5diOH3M\nt+9R4sQgyHzBNBs3Q6ZEKiFWSm+XbSC4ek1gp1BKQdt26Jbczkow167anMEK2yEvFdWxU0JFRWBO\nH0e1ixjMTdcxu1I4QJRSsHU7RiaaFWG6T6DaAvQK8qit2zCyEnll9J6CxlZfGgguhJLv0YoxVy7D\n+XPQ4E9FXuLEwPScRG31f33BbNTWbTY3ISyf7hPQ1hn4aWSiWTnm9AlUkGG8Am3boPuEhFxXwunj\n0NrhW0VeosTAXL8Go4O+9rqZD7V1O+a0xDpXguk+EYpg24nmpCw+WwHGx1j0QqiN1VBSCmdHAj9X\n0rAhIv+iIIkSA3q7oCF41xawE03facyNG8GfK0GY8ctw/izU+7/YbDaqcj2Ur4ORgcDPlTjC8gzA\neoniwS0b03XM13uUKDGwsejgww+A3Qu2umbFm0+nlp5TdnV4STCLzWZjPTiZaJaDuXgerl2FmvpQ\nzqfaOiWct0yMMdBz0teWO4kSA7r9vTiLodo6MT0nQztfErCCHd49KpQBC8ug+yS0dQa22Gw2qm0b\nplu+R8uiUDK9sdq3IRMlBqFPNC0d0NcV3vkSgDl9PJTkcQHVuk0Ee5mY3lOo1o7wTti6DXpPSW5n\nOfSetMljHwU7MWJgxi/DxfO2b1BIqJYOTM+p0M6XCHq7CFWwm7dCtseuTBeWhOk9hWoJTwxUZRWs\nqyQ32B/aOV3H9HT5fo8SIwb0dkFTa2ixaABa8hONbLG4JMyVMbg8BlvqQjunKl8HVRtsS3NhafSe\ngpb2cM/Z0sGU5A2WTBDeW2LEwPSfRjWH+wFWa8rtauQheaJZEr1d0Nzm+85mi9LSLh7cEjFjl+xW\npJvDE2wA1byVKckbLJ3eUzZM7SOJEQP6uwNfOj8XSiaaJWP6ukIXbMCeU3I7S6PvFDRvDV2wVUu7\niMESMRfPw8SErWb0kcSIgenrRkUgBrR2WJUWFqe3K/zwA/ncTq+IwVIwvf7HopdEcztTPSdlJfJS\n6D3le/IYEiIGZmoKhvpCWXk8GzvRiBgsBTvRhC8GtLZDb5dMNEuht8u2lg6bDZvs/y+cC//cjmF6\nTgXyPUqEGDCchQ3VqDVrwz9381a7ElkmmgUxE9fh7LBvTbWWg6raCKvK4Nxo6Od2DdPbFW5ZaR6l\nFCVtnVaMhIXp74amrb4PmwgxMH2nA7k4S0FVVMGatdJbZTGyPVDbGE6rkLloabern4V5MdeuwvnR\nUMuzZ1LS2omR3M6imGwPqsn/KEgixIBsRPmCAo2t0pZiEUxvNMnjAqqpDSP3aGEGeqGuOdzy7BmU\nbN0mYrAIZnLCPngGINiJEAPT141qbovs/KqxDdPfHdn5naC/GyK8RzS2YrLd0Z3fAUx/tA9VJa2d\ndgc8YX4G+2BLXSAediLEIKgY2pJpbLVPVcK8mGw3qrEtsvOrxlbIyj1akIjKswtk6hrh4jnbil6Y\nE9PfYz/LAeC8GJjLl+D6Vd9rbpeDamqVEMQCGGOgvyfSiYa6Rjg7Yt1sYU5MxOFWVVJi964e7IvM\nhtiT7QmsatJ5MaDvNDS2hdZhcU7qmmBkEHNjMjob4sz5M1BWZvcXiAhVusq2wZD+N3NiBbs7WsHG\nenDyYDU/QQp2qR+DaK3fD/wxVlw+53neH87693cB3wAK2aG/9jzv9/04d9RPMwCqbDVUb4GhbORf\npljS3w0RhogKqMZWzEBPNGsd4s75M7AqWsEGoEGKMRYkG5yHXbQYaK0zwJ8ADwMDwIta6294nvfm\nrEN/4nneLxZ7vlvo64YO/7Z+WzH5J5qohSmO2MRk+AsCb6Gx1YarhFuJjWC3kPvh61GbEUtsSPwa\nbNoSyPh+hIn2ASc8z+vxPG8S+DLw6BzHBRLHsRNNhMnjPKqxDaRaZW6yPbGZaIwk+uck6kqiaRpb\nYUAEe06yPdDQElhI3A8xaARmZnz686/N5n6t9UGt9Xe01rt9OG+kbShmoxpbMfLUOSfx8QzaRLDn\no787Ft8jNm2Ba1dtu3PhJmwlUVtg4/uSM1gCLwMtnueNa60/AHwdmHMnZ631fmB/4XfP86isrJxz\n0Km+01zZVEPV5mDcpuUwtWMPl7/yF/PaGgfKyspCt8/cmOTi6BCV23ajyspCPfcttqxbx8Url6ko\nUajyiqLHi+J6BsWl4Szlv/RRSiP8e8rKyqiqqmKseStrz5+htK4hMlviyPhIlpKO7axe4j3SWj8x\n49cDnucdWOh4P8QgC8xsONOUf20az/Muz/j5e1rrP9Nab/I875auVHmDD8x46fGxsbmfEnLHDmMa\nWpjv38PErK3AjF3k0sgwam151ObMSWVlZejXyvR3Q3UNl69fh+vXQz33nDS0MHbsMKqzeOc0iusZ\nBGZqitxgP+NVm1AR/j2F65mra+LKyaNkYhD+jRNT3SfJ3PMgE0u4R5WVlXie98RyxvcjTPQi0Km1\nbtValwEfAb458wCtde2Mn/cBai4hWDbZaFcez0RlMlDfLJUQszCDfZE0p5sP1dCCkcVnNzM6BOs3\nolavidoSS0OLLBCchcnl7MLWAMNERYuB53lTwCeB7wOHgS97nndUa/2Y1vo38of9stb6kNb6VWwJ\n6q8Ue14AM9CHitNEIzXStzLQh2pojtqKt5A+Urcy0GsfZGJCoQRYmMHZEVi7DrWu+PDmfPiSM/A8\n70lgx6zXPjvj5z8F/tSPc93EQC/Ux0cMaGqVBOUszGAv6u4HojZjGtXYSu7gC1GbESvMYNwE23oG\nxphoF5PGiQBXHhdwdgWymbhuN8IIcXP1xVAN4hncQhw9g4Ee2X9iJgPxCuVRuQGUgovno7YkNphs\ncD2JCjgrBgwP5Lv3hVUQtQSa7KImmWgs5sYNODNs+83EhaoNgEw0MzEDvagYedhKKVlvMJtsj51f\nAsRZMbAf4Bg9cZLfUaukRLbuKzA6CBs3o1ZFW1I6E5lobsbkpmAkC/XRbGgzH5LovxnTH3zXX2fF\ngMG+WCW9ppEE5VsM9EGcQkR5ZIHgDM4MQ+WGaLaMXQj5Hk1jpqbyHnaw6y6cFQNbshjDiaa+CTMk\nLXghnzyOpWC3iGdQIGaVRAWkdcgMzg5D1QbbEDNAnBUDBvvjOdHUN0ub5AJx9QwaWjEDItgQv/Ls\naRpaYKDX1tenncFwwnhOioG5MRm/xGQeVddkvRbBlizGKDE5TX0TDPVLoh+sZxBHwS6vgPIKW1+f\ncsxwPyqEuc5JMWCkkJj0fx/QoqlvtvsapBwzNQUjA4Fs3F0sal0llK6Ci5Lot2sMYijYkPey5cGK\noazdqS9g3BSDmOYLAFi/EW7csL3H08yZYajaiFodbJxzxUg4z1YSDfXHrpKogKpvwqT8HgGYoX5U\nCA9VToqBGeiLZ76AfOlifVPqJxoG45mYLCDhPODMCFRUodbEs7EidTacl3rEM1iAOHsGyEQDhcRk\nfO9RIW+QamLWRHA2tjIv3ffIXLkMkxOwflPg53JSDGxiMuYTTeo9g7549Y2ahRXsdN+jOHvYwPT3\nKNWJ/qF+qG0MpUeTc2JgpqZsK4oYJiYLqPpmeaKJvWfQLJ7BQG+sPQMqN4DJQYrzb2YoiwohRAQO\nioFNTG6IT+/1uahvSnUVhMnlYp2YBGBjNVwdx4xfidqSyIi7h23zbylP9A/3h/bg654YxDzOCcDm\nWrh0AROHnb2i4OwIVFTGNzFJfjOi2sbUegdvCXZ8xQBA1TWmekW/eAYLEPenGQCVKYGaeqvqaWQg\npn2jZpHqBOX5s7C2HFW+LmpLFqauya7ATStDWfEM5mWgL97hhzxpTlCaofgLNpCfaNJ5jwqJybiT\n5vzbdIO6mvpQzuecGLjgGQDpTlAOZWWiiTlmKIty4KHKCnZKw0QhNagr4JQYGGPin5gsUJ/etQZh\nxjmLQjyDqK1YnDTn30JqUFfAKTHg/FlYvcY2sYo5Ks1VEMPhrJgsmtp6ODtiGx+mDDOcDaXFQbGo\nkhK7te3IQNSmhE5YDeoKuCUGrkwyYDeiGB2ycb8UEeaKyWJRpaugusY2Pkwbg/3ufJfSuqI/pDYU\nBZwSAzOUDVUpi0GVrbZN60aHojYlXEJcMekLdY2p8+DMtXEYH4NNW6I2ZUmolPYoCqtBXQGnxMAp\nzwBSmUQ2w+4INuSTyGl76hwegJoGu9bCBeqb0tkWXjyD+TFD/aja+Mc5C6SyBW/IH+CiSeFTpxkM\n94mzWFQKizGiCLc6JQZOTjRp+xA75r2lUrAdu0fUNsLIgN1/IS1EEG51RgzMxHW4eN6WmjmCqm+y\nk2OacCivA1jBHs6ma6/dwfD63fiBWrMWKqrg7GjUpoRGFOXZzogBI4OwudaWmrlCbbr22jW5KZsw\nr22I2pQlo8rXwZpyW7acElwpK72JtIXzQmxQV8AdMXDNtQVUZRWoDIxdjNqUcDg7anfOinNH2bmo\na0xNHymTy9mafYcEGwqJ/nTcIxDPYEFcKiu9iboUdcZ0ULABVG0jJi3VKudGYV2VDb24RNo6zIbY\noK6AM2LgXPI4j23Bm46Jxm3BTsc9YsihxWYzUHWNqcm/hd2groAzYmDjnO59iAsJylTgqmeQponG\nlb5Rs0lTH6kz4TaoK+CEGNgGdVmbkHUMVdeYmlins55Bbco8Awe/R2zYBBMTmPHLUVsSPEPhNqgr\n4IQYMHYBMhmbkHWN2hR5Bo6G8thcA2MXU9EZ07audu8eKaVs0jsFoh12g7oCboiBq5MM2I6L585g\nJpPdGdNcu+pUv5uZqEyJXb+Shs6YjnrYkKL8W0TznRNiYIYHUDVulcIVUKWltjPmaMI7Yw4PwJZ6\nd/rdzKa+KfETjbk6DlevwMbqqE1ZGbWNqfCyw25QV8CNb67LngGkorzUOFqlUkCloXRxKAu1DjWo\nm414BoHixKfC2UqiPKlwb4cHnGoieAspKC81w241qJuNSkFlXpT7gTghBi7HOYH8Uvpkf4gZzjq3\nqnUmqjYF5aWue9i1DTAymOyGdRHuB+KGGJwdCX0Bhp9YzyDZIQhn69cLFBrWJbmP1FDWjX2P50Gt\nXgOVyW5YF+X3yA0x2FiNWrUqaitWTm2yJxpjjE0gOywGal0FrCqDi+eiNiUwXA+3AslfExJBg7oC\nboiBw08zkG9Yl8nY9RJJ5MI5KCtDlVdEbUlxJHiicbVB3WzsavHketlmUDyDBXH+aQYSPdG42u9m\nNolO9J8/C2srUGvKo7akOGoTnn8bDr9BXQEnxMB1zwBsJURS8wZO9sefi7oE17EPi2DHnaga1BVw\nQgwS4RkkuXTR8cRkgSS3sna2b9RsklxeGlGDugJOiIE80cSbRCQmIeGegdsJ/mk2VsP4FbuaOmlE\n1KCugBtiEMECDN9J8rZ9CfEM2FwH589iJieitsR3rGfgdvIYsKunaxoSKdpmKJoGdQWcEIMoFmD4\nzvREk6yGdWZy0lYTba6N2pSiUaWltoPpSAL7SDm618RcJNbLjvgelfoxiNb6/cAfY8Xlc57n/eEc\nx3wG+ABwBfi453kH/Ti3K9w00TS2RG2Of4wMQnWN/fuSQKHqq7E1akt8w0xch4vnodp9wQYSG84z\nQ/1k9r0zsvMX7RlorTPAnwA/B+wBflVrvXPWMR8AOjzP2wY8Bvx5sed1ktoEbryekCqVAolcLT4y\nCFvqUCUlUVviD0kt0464XYgfYaJ9wAnP83o8z5sEvgw8OuuYR4EvAnie9wKwXmudkMeUpZPEXc8S\nU6VSIIltkh3vGzWbJIaJzJWxyBrUFfBDDBqBvhm/9+dfW+iY7BzHJJ8klsW53vxsFnY9SLLuUeIE\nu64RRrJ2VXVSyBdhRJkfjV2gV2u9H9hf+N3zPCorKyOzx09utG/j6jM/iOzvKSsr8/3cY2eGWPu+\nRylNyD3Kde5gbGSAioqKRb+YQVzPILhyboTSXXewOua2Lvl6VlZycV0VFRNXyWypC96wELh+4Sw3\nmttY5+M90lo/MePXA57nHVjoeD/EIAvMzIg25V+bfUzzIscAkDf4wIyXHh8bGyvayDhgqjaRG+jj\n0qVLkTwBVFZW4ue1NMaQy/YyXrURlZR7RAZjYGygH1W1YcFj/b6eQTHV183U/e9hIua2Lud6mpp6\nLp86jlqzLmCrwiHXcxKqa337PFVWVuJ53hPLeY8fYaIXgU6tdavWugz4CPDNWcd8E/gHAFrrtwMX\nPM8b9uHcTqEqqqAkQQ3rLl8CDFSuj9oS31BKJWq1uO0om5B1IDNIWt4gygZ1BYoWA8/zpoBPAt8H\nDgNf9jzvqNb6Ma31b+SP+S5wWmt9Evgs8Iliz+ssdU0wmJAPcQzinEGQqI1uxi4CCiqqorbEX+oS\nVpkXYYO6Ar7kDDzPexLYMeu1z876/ZN+nMt17ETTj9qxN2pTiiYxbShmkyDPoNCGInmC3UTu9Zei\nNsMXom5QV8CJFciJoq4xcZ5B0rA985Nxj6JucRAYSVp4FnGDugIiBiGj6pqSM9Ek1jNogqSsB0nY\nGoNpqrfA2CXM9WtRW1I8ETeoKyBiEDZ1jclpWJdQz4At9XBuFHPD/T5Sdk/d6Ccav1GZEhtWScCD\nVVy8NxGDsElIw7q4xDmDQK1aZVsljyag4C0pravnIin7T8SkiaCIQcgkpjPmmWFYvzHyOGdg1DU5\nX60yLdgJWZg1G5WQRL8Z6o+F9yZiEAW1CQgVJawNxWwSsetZ0gU7KX2kYvJdEjGIgCTsh2yG4/E0\nExhJeOqMySQTFElYeBaHBnUFRAyiICkTTQySXkGRhPJSW+2VcMEeztpV1q4So4WbIgYRkITy0sSW\nlRZIQs/8pJaV5lHlFbB6DZw/G7UpK8ZWe8XjeyRiEAX58lKnn2iGBxLtGbB+I9yYtG68oySudfVc\nuL74bLg/8jYUBUQMIsA2rCuFS242rDNXx+HaVdgQfZwzKJRS7nsHMSlZDBLXE/1xaFBXQMQgKlzO\nGwxloaYelUn2x8flBKW5Og5Xx2FDddSmBIvznkH0DeoKJPvbHGPsE42bFUWJryQq4PJq8eEs1DQk\nX7Br3a3MM1NTMDoUm4Wbyf6kxJn6Jrc9g6THonF7C8w4JSYDxWUPO2brQEQMIsJlzyDp9evTuByC\nSHIbiplsroWL5zET16O2ZPnEpEFdARGDqKhrcnaiSXxZaYGaBhgdsu68ayRwd7O5UCUltt3GyEDU\npiybuDSoKyBiEBWba/MN6yaitmRZmFzOfvFi9CEOClW22paYnnWvYV3cJppAcbXqK2bVXiIGEeFs\nw7qzI7CuCrVmbdSWhIODE43J5dITJsLdqi8zGK9CDBGDKKlzMImclnxBHicnmgvnYG05am151JaE\ng6u5HfEMhAIuJpHj0m43NFycaFKSLyjg4sKzODWoKyBiECUulpcO9ceqAiJonJxohrKoBPckugUX\nG9bFqEFdARGDCBHPwAFc9QxiFH4IGtvepcSp9i5xXAciYhAlLj7RDKbLM2BDNVwdx4xfidqSJWOG\ns6jaFN06y/BSAAAWjUlEQVQjyOffHHqwGuqPnWCLGETIdMO6i+ejNmVJxDHOGTQqk7FtoF3yDoaS\n3bp6LlwL55mhflR9c9Rm3ISIQdS4FIYYsk214hTnDAOX2lKYyQlbTbS5NmpTwsW1thSD/VAnYiDM\nwKmJJm35ggIuTTTDWdhSZ9expAiXdqYzNybtep3aeDSoKyBiEDUudcYcjF+cMxRqGzHDbtwjk7ac\nToFah3IGw4NQXYMqXRW1JTchYhAxznkGKZxolEuLAwf6YheLDoUthfYuk1FbsjhDfbEUbBGDqKl1\nyDMYis9GHKFS2wCjg7bNQ9wZ7IMUioEqXQWbtjjR3sUM9sXyoUrEIGo218KFc7FvWDcd59wSrzhn\nGKg1a6G8Es6NRm3KoqTVewPs07YL4bwYJo9BxCBybMO6WttYLM6MDsGmLahV8YpzhoYDSWQzNWWf\njNO2xiCPK+Wl1jMQMRDmwoXy0rQmJvM4Ua1yZhiqNqBWx2PnrNBxQbBj3FFWxCAGuJBENoN9sVs+\nHyoutLIe7IWGlqitiAxV64Bgnx2B8opYdpQVMYgDLpSXpjV5nMcFz8AMpjhfANOeQazbu8S40aOI\nQQxwIdaZ2gVnBWobbagszgz2pVqwqVwPGLh8KWpL5iWu+QIQMYgHec8grk80xphYNtYKleotcGUM\nc+1q1JbMi/UM4jnRhIFSKv6iHePcm4hBDFAVVVC6Kr4N6y6eg1Vl1s6UojIldr1BTMN504KdYjGA\nfP4txuE88QyExYlzC96BeK6YDBvV0IIZ6I3ajLk5fwZWr0Gtq4jakmiJcUWRMUY8A2Fx4rzXrhno\nQTW0Rm1G9DS0QFzFYLA/3fmCPKqhGTPYF7UZczOW33ynckO0dsyDiEFcqGuK71qDgT5oTG/JYgHV\n2IIZiOdEE+fwQ6g0tEC2J2or5ibvFcS1BbyIQUywnkE8w0Qm24NKcf36NA5MNKlncy1cvoS5Oh61\nJbcQ9wS/iEFciOmiJmOMDY2IGOQnmouYazGcaIbEM4B8or++OZ7hvBivMQARg/hQaFg3cT1qS27m\nXD4xmeJKogIqU2IbjMUxVDSQzm6lc6EammOZ6I97KE/EICZMN6yLWwte8QpuQjW0YGIWKjJjl2Bq\nCtZvjNqUeNDYGk/PIOZJfhGDOFHXZHeqihFmoFfyBTNpaImfZzDYF+vEZNiohtb4CfaVyzB+Bapr\nojZlXkQMYoRqbIGBeH2IyfbYJy0BKFQUxeupM+7hh9CJYwlwtgcamlGZ+E658bUsjTS2YfrjJQbi\nGcwijhNNfzc0iWBPs2kzXL+GiVGPIpPtQTW1RW3GgpQW82at9UbgL4FWoBvQnuddnOO4buAikAMm\nPc/bV8x5k4pqasN87YtRmzGNyeVSu43ivGzaAlfHMeOXobIyamsAMNluMnffH7UZsUEp9VbeYPve\nqM2xZLtj72EX6xn8DvCU53k7gB8BvzvPcTlgv+d5d4kQLEBNPVw8F59maGdHoKISVb4uaktig8pk\nbHlgTLwDYwz090DMnzrDxib643GPIO8ZJFwMHgW+kP/5C8CH5zlO+XCuxKNKSmy1QVySX9keqSSa\ng1jlDc6fgbIyVOX6qC2JFzEK5xlj8rm3tqhNWZBiJ+gaz/OGATzPGwLmS5Ub4Ada6xe11r9e5DkT\njWpsi00lhM0XxPtpJhLiVFHU3x37SSYKVGMrJi7FGOdGoWw1qjLea3UWzRlorX8A1M54SWEn90/P\ncfh8Dfkf9DxvUGu9BSsKRz3Pe2ae8+0H9hd+9zyPypjEZsPgWscOcqMDlAfwN5eVlS3rWl4ZGaD0\ntntZnaLrvxQmO3Zy/ciXl309g+Da6CBm6zbWJuAe+Xk9c9t3MzbQR0VFReQlt5PHD3G9pYOKkO+R\n1vqJGb8e8DzvwELHLyoGnuc9ssDJhrXWtZ7nDWut64CRecYYzP9/VGv9NWAfMKcY5A2eafTjY2Nj\ni5mZGMyWOnI/e5qpAP7myspKlnMtp3q6mNr/QSZSdP2Xgtm4hVzfaSYmJpZ1PYMg13Uc9tzNjQTc\no+V+PhfCZEoxGMayfaiIF+PlTh6F+qZQPyuVlZV4nvfEct5TbJjom8DH8z9/DPjG7AO01uVa64r8\nz+uA9wGHijxvcmlsg2xP5LuemakpGMlKJdFcbKyGyUlyl24pnAsd09+NkrLSW1BKQUNrPPJvjqzV\nKVYM/hB4RGt9DHgY+AMArXW91vrb+WNqgWe01q8CzwPf8jzv+0WeN7lUbQCl7O5iUTI6COs3oVav\nidaOGGInmmZy/acjtcNMTsLokAj2PMRlMyJbSdQWtRmLUtQ6A8/zzgHvneP1QeBD+Z9PA3cWc540\noZSyZYL93bChOjpDpCfRgqiGFqb6u6G5Izojhvphcy1qVVl0NsSZxhbo7YrUBHNj0vYbi3G30gJS\n7hlDVGP0vVVMthfVIE+c89LYxlTUE01/d+xXtUZJLHoUDWWheguqbHW0diwBEYM4UvAMIsT0n4am\nrZHaEGdUcxtTPaeiNcKBVa2R0mjXGphcLjITjCP5AhAxiCUqDj2K+k6jmtujtSHONFnPwOSmIjNB\nPIOFUesqoXwdnBmOzoj+bifyBSBiEE8aWmA4i7lxI5LTm6vjcOkC1NZHcn4XUOUVZKo2wMhQdEZk\npQ3FojS3Q190iX4XGtQVEDGIIWr1ati4GYYHojGgvxsaWuzOXsK8ZFo7IKKKIjN2Ca5fs43zhHlR\nzVsxfRHmdhwK5YkYxJWmVky2O5JTm74uVLPkCxajpLUDE9VTZ36SiXp1bdyxYhCRYI9fhiuX7Q6G\nDiBiEFNUY4QLZvq7QcRgUUpaO6ObaBzoghkLogwTZXvzHrYb06wbVqYQ1dQW3UTT24WSSqJFKWnt\niK7qq79b8gVLoboGro1HstGN6e1yqghDxCCutHREsmDGTE3ZDW2kxcGiZLbUwdUrNn4fMqbvtDOJ\nyShRmYwVzSgerHpOQquIgVAs1TVwYxJzIeS2FENZ2LAJtaY83PM6iMpkbDitL9z1BubGJAz22hCI\nsCiquT0SL9v0nkK1RLhCfZmIGMQUpRS0tENvyBNNz0lUa2eo53QZ1boNE/bis2wvVNdK36il0rw1\n/O/RxHXb38uRNQYgYhBrVGtn+BNN7ykQMVg6rR2YnpOhnlIEe3motgi+R9keqG1ErVoV7nmLQMQg\nzrR0YMJ+ouk+gWqTiWapqNZOiESw3Qk/RE59C5w/YxdThoTpcStEBCIGsUa1doQ60ZipqXxZqcSi\nl0xtA1y+FGq1iuk5JZ7BMlAlJTaJHGZBhoOCLWIQZ7bUwbWrmEsXwjnfUD9sqEaVrwvnfAlAZTKh\n5nbMjRu2vbisA1kWqm0bpvtEaOcTz0DwlbCTyKZbYtErQbWEGJMe6IXqGtSateGcLym0dtpSzxAw\nkxP2wcqxtToiBjFHtW3DnA7piabnBEi+YPm0dmC6Q5poek7a8KGwLFRbZ3ieQW8X1DbYHmMOIWIQ\nc1T7Dszp46Gcy8aiZaJZLqptG4Q10Zw+Dlu3h3OuJFHbAGMXMVcuB34qc/o4auuOwM/jNyIGcad9\nO3QdwxgT6GnM5KRNHjsW54wFtQ1w/RrmwtnAT2W6jqHa3ZtookZlSmzItScE0e46Zr+3jiFiEHPU\nhmpYvdruoxokvaesayux6GWjlIL2HXYSCBBzddxu1CJtKFaE2roDcyrYewR5z8BBwRYxcIGt2zGn\nA55ouo6hOnYGeo4ko9q3Bz/RdJ+A5nZUqTsLmeKE6tyJ6Xoz0HOYsYu2bXVtY6DnCQIRAwdQITx1\ncupNaBcxWCmqfScmaM9AQkTF0b7ThlyD3BO56zi0dTrTtnom7lmcQlT7DkxXsElkc+pN8QyKYet2\n6OuyTeQCQsSgOFTVBlhXCYP9gZ3DnD7mZPIYRAzcoKUDBvsw168HMrw5NwpTN+wiN2FFqLXldker\nvu5AxjfG5BOTbk40cUF17Ao0VGROHkV1uHmPRAwcQJWttknDgPIG5pSdZGQLxeKwHlxAE81wFlaV\noTZWBzN+WujYAaeOBjK0uTEJ3Sehc3cg4weNiIEjqG17MMcPBzP4qaMSIvKDbXswJ4K5R+b4IdSO\nvYGMnSZUxy7MqYAEu+cU1NY7285FxMAR1I69mOOHAhnbHHsDtV0mmmJR2/fC8cPBrAk5dhjkHhVP\nYwtcuoC5dN73oc3xw6hte3wfNyxEDFyhYxd0n/Q9QWkuX7K169KTqGhU9RZYs9ZuG+ojxhjrGYgY\nFI3KlFgP7pj/D1bmhIiBEAKqfJ1d6ep324Pjh6FzF6q01N9xU4raHsBEMzoExkBNvb/jphS18zZ4\n8w1fxzS5KZuL2OZmvgBEDJxCbfc/b2COvYHacZuvY6aa7bfBMZ8nmuOHUNv3SILfJ9SO2zFvvu7v\noP09ULXRlq86ioiBQ6gdezF+TzQiBr5SyO34mjc4fgjkHvlHYyuMj2HO+9dLyhx9zXocDiNi4BLb\nb7MrKCf8WW9gxi7BuVFpTucjqroGVq+x+w74gDEmP9Hc7st4Qn5Dou23YY755x2Yw6+g9tzl23hR\nIGLgEKp8ne286FNVkTnyKmzfa7cFFHxD7b0bc+gVfwbr77brC2ob/BlPALDietQfMTDXr9s2FDvc\nFmwRA8dQe3ycaN54CXXbvf6MJUyj9t6LeeMlX8Yyh15G7b3Hl7GEt7CC/bI/fYqOH4LWdrsK3WFE\nDBxD7b0bc7h4MTC5Keva3iYTje/svA16TtqW00ViDr0s9ygA1JY626fIh60wzeFXULvdDhGBiIF7\nNG2F8SuY0aHixuk6DhuqUZu2+GOXMI1avQY6dsLRg0WNY8avQE+XzRUJvqNufxvm9ReLHsccfhW1\n924fLIoWEQPHUJmMDRUVGYYwEiIKFLX3HswbLxc3yNHXoHOnc3vpuoIfYmCG+uHaODS3+2RVdIgY\nOIi65wHMy88WNYZ5/UXU7SIGQVGYaMzU1IrHMK88h7rzPh+tEm6iYyecHS2qxNS89Czq7gec3L9g\nNu7/BWlk913Q37PiPXfNYB9cviTtkANE1dTDxs0rXoBmrl/HvPEy6p4HfbZMKKBKSlC33Yt55bkV\nj2FefhZ1zwM+WhUdIgYOolatQt2xD/PyT1f0fvP8j1H73mn7tAiBofY9hPnZT1b25jdehK3bUZXr\n/TVKuAl137swzx9Y0XvNUBbGLkLnLn+NiggRA0dRb3sH5qWnl/0+YwzmZz9G7XtXAFYJM1H3PoR5\n9XnM5PKbC+ZefBr1tncEYJVwE7vugPNnbOx/mZiXn0XddX9iHqpEDFxl1x0wlF1+VVHXMShdZRev\nCYGiNm22mxIdWl4i2YxfgaOvoe66PxjDhGlUSQnqbe9ctndgcjnMcz9E3ZechyoRA0dRpatQDz6M\nOfDdZb3PPPsU6r53SdOzkFD3v5vc099f1nvMs0+h9t6DWlcRkFXCTNT9+zHPH1jeArSjr0FZvoQ4\nIYgYOIx61wcwz/1wyXsj5y6et67tO98XsGVCAbXvndB9ArPETdhNbgrzo2+jHv6FgC0Tpmluh8r1\ncPCFJb8ld+C7qHd/IFEPVSIGDqO21EHHLszPfryk469//xuoe9+BqtoYsGVCAVW22or2U99c2hte\nfxEq18s2pCGilCLz/r9L7smvLqnbrDk3CieOoO7bH7xxISJi4DiZh38B8+RfL7oDmpm4zsRT30Q9\n8mhIlgkF1Ls/gHnpaczYxQWPM8aQ+5uvo97zoZAsE6a56z64ctlu9rQI5tt/iXrHe+1K8wRR1PZW\nWutfBp4AdgFv8zxvzqY5Wuv3A3+MFZ/PeZ73h8WcV3gLtesOqKnH/Og7qPd9eN7jzJNfpXTn7eTq\nmkK0TgBQVRtR978H89UvoD7+W/MeZ158Gq5dRb3toRCtE8Buh6ne/3fIff3/kPnt/zDvIjLTfxpz\n8AUyv/8/QrYweIr1DN4AfgmYN06htc4AfwL8HLAH+FWttfjAPpLR/xjzvb/CXLow57+bbC/mb7/L\n2o/9ZsiWCQXUox/FHH4Vc+LInP9uro1jvvJ5Mh99TFqKR4R68GEwuXmLMowx5Ly/QH3oV1DlyUvu\nFyUGnucd8zzvBLBQFmUfcMLzvB7P8yaBLwMSq/ARVd+EeugRcp/9j7dsfGOuXSX3hc+gHv01MtKU\nLjLU2nIyH/kn5L7w3zGXzt/0byY3hfnin6J234nqdHcPXddRmRIyH/stzLf+35wJf/P1L1nP7Z3v\nj8C64AkjZ9AI9M34vT//muAj6sN/H7V+oxWE/GRjLl0g918+jWpsTewH2CnufgC17yFy/+n3pteH\nmPErmM9/BnP5Euqj/zRiAwVV34T6e/+Y3H/+1PR+42Zygtx3PMyLPyHzz/4dqrSo6HpsWfSv0lr/\nAKid8ZICDPB7nud9KyjDhOWhMhn4R/8CvM+R+3efgE1bYHQI9fAvoj780USVwLmKUgr1i79Gbk05\nud//V7CxGs6Nom57G5nf/DSqTLqTxoHMA+/BbNhI7n/8B1hTDpMT0NpJ5l/++0S3B1lUDDzPe6TI\nc2SBlhm/N+VfmxOt9X5g/4zz09AgW/4tmX/97xf858rKypAMSQcrup4f/4T9T7iF2Hw+GxrgvT8f\ntRVFobV+YsavBzzPO7DQ8X76O/M9er4IdGqtW4FB4CPAr843SN7gA4XftdZ4nveEb1amGK31E3It\n/UOup7/I9fSPlVzLonIGWusPa637gLcD39Zafy//er3W+tsAnudNAZ8Evg8cBr7sed7RYs4rCIIg\n+EtRnoHneV8Hvj7H64PAh2b8/iQgzfMFQRBiigsrkA9EbUCCOBC1AQnjQNQGJIwDURuQIA4s9w1q\nKb04BEEQhGTjgmcgCIIgBIyIgSAIguBraamvSHM7f9FadwMXgRww6Xnevmgtcgut9eewRRHDnufd\nnn9tI/CXQCvQDWjP8xZuTSrMdy0fB34dGMkf9ql84YmwCFrrJuCL2MXBOeB/ep73meV+PmPpGUhz\nu0DIAfs9z7tLhGBFfB77eZzJ7wBPeZ63A/gR8LuhW+Umc11LgD/yPO/u/H8iBEvnBvCvPM/bA9wP\n/GZ+vlzW5zOWYoA0twsCRXzvd+zxPO8Z4Pyslx8FvpD/+QvA/D3EhWnmuZawcMNLYR48zxvyPO9g\n/ufLwFFsp4dlfT7jOjlIczv/McAPtNYvaq1/PWpjEkKN53nDYL+QQE3E9rjOJ7XWB7XW/0trndwm\nQAGitW4D7gSeB2qX8/mMqxgI/vOg53l3Ax/EupHviNqgBCJ12ivnz4B2z/PuBIaAP4rYHufQWlcA\nfwX887yHMPvzuODnM65isKzmdsLi5FeF43neKPA1bChOKI5hrXUtgNa6jreSn8Iy8Txv1PO8wmT1\nP4G3RWmPa2itS7FC8L89z/tG/uVlfT7jKgbTze201mXY5nZL3FFcmI3Wujz/1IDWeh3wPuBQtFY5\nieLmuPY3gY/nf/4Y8I3ZbxDm5aZrmZ+sCvwd5PO5XP4COOJ53n+b8dqyPp+xXYGcLy39b7xVWvoH\nEZvkLFrrrVhvwGDLib8k13N5aK3/L7a1ejUwDDyO7cv1FaAZ6MGW7s2996gwzTzX8t3YWHcOWwb5\nWCHeLSyM1vpB4CfYbYhN/r9PAT8DPJb4+YytGAiCIAjhEdcwkSAIghAiIgaCIAiCiIEgCIIgYiAI\ngiAgYiAIgiAgYiAIgiAgYiAIgiAgYiAIgiAA/x+aA6wAawFyUAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7fad0ed8a898>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"%matplotlib inline\n",
"x = np.arange(0, 20, .1)\n",
"y = np.sin(x)\n",
"plt.plot(x, y)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Simple usage of splrep/splev"
]
},
{
"cell_type": "code",
"execution_count": 24,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7fad0e85f7f0>]"
]
},
"execution_count": 24,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAEECAYAAADK0VhyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXl0XNd93z93QAAkQHADiW2wEgApkqJEbZRkeWHiRVaj\nRraT3sQ9TeK0J3Hr6LRpz+lpEjdHak7Tk/QkqeO6OXETN3XSpPZLE8eKnXhRZNqWbUmULFHcSYBY\nBzs3AFywze0fdwYaQgCIwbzl3jf3c46OgOGbd3+Y9+Z+3/1tVyilcDgcDkdxk4jaAIfD4XBEjxMD\nh8PhcDgxcDgcDocTA4fD4XDgxMDhcDgcODFwOBwOB7DBj5NIKT8HPAmMep53zwrHfBp4ArgOfMzz\nvDf8GNvhcDgchePXyuBPgMdX+kcp5RNAu+d5ncDHgT9c64mllEcKts4BuM/Sb9zn6S/u8/SP9XyW\nvoiB53kvAldWOeQp4E8zx74MbJVS1q7x9EcKs86Rw5GoDYgZR6I2IGYcidqAGHEk3zeEFTNIAgM5\nv6cyrzkcDofDAFwA2eFwOBz+BJDXQApoyvm9MfPa28j4uo5kf/c87xngmQBtKxo8zwP3WfqG+zz9\nxX2e/uF5HlLK3JeOep53dLX3+CkGIvPfcjwH/BLwRSnlI8BVz/NGlzswY/DRnJeeGRoa8tFMUHOz\npD/1DKK5g8RP/Yu3Xr95g/Qf/BdEczuJf/Lzvo5pAlVVVUxNTUVtxppJf+6/oaavkfjEJxGlpQCo\n9ALqf/931NVLJH75PyES0S1ubfs8gyD99S+hXj5K4pefRWzZ/tbrf/9XqB+8QOI//BaismpN53Kf\np380NDTged6z+bzHl2+SlPIvgO8De6SU/VLKn5dSflxK+YsAnuf9HdAjpewCPgt8wo9x14v6u7+E\nTZWIJRO+2FRB4hf/Peqlb6G6z0ZknQNAHX8F1XuexL/61UUhABCJEsTH/jXcuon63vMRWuhQY0Oo\nr/0/fY1yhAAg8cRPIDoPoL70ZxFZ58gXYUELa+XnykBdHif9G79M4tc/hajetewx6WPfRT33f0k8\n82nEhrA8acFjy5OXmp0h/czTJH7mE4j99y1/TP9F0p96hsR/+gyiamvIFmps+TyDQClF+vd+HXHw\nQRIf+NDyx1yfJv3r/4rEv/0NRFPbHc9ZzJ+n3zQ0NMDKnpplKboAsvrS/0G854kVhQBAPPhO2Lod\ndey7IVrmyKKefw6a21cUAgDRvBtx+N2or3wxRMsci5x6HSavIt77j1c8RFRuRvz4R0l/8Y+x4KGz\n6CkqMVBjw6iTryGe+MiqxwkhSHzwI6iv/7W7iUNGzc+hXvgKiR//p3c8Vjz+EdRLR1E3pkOwzJFL\n+htfQjz+EURJyarHiXc9DpfHoed8SJY51ktxicG3v4Z47L2IjRV3PvjA/SAScOLV4A1zLKJefREa\nmhHJ5jseK7ZXIw4+gPruN0KwzJFF9V+E4UHE4Xfd8VhRUoI48gTqW18NwTJHIRSNGKjZGdT3/wHx\n7g+u6XghBOLxD5N+/rmALXNkUUqhnv9bEqu4HpYi3v8U6oWvoBYWArTMkYv65pcR730SsaH0zgcD\n4rH3od48hppcrUmBI2qKRwxefRFaOxA19Wt+j3jgHdDXjbo8EaBljkV6u+DGNBx8YM1vES0dsH0n\nnHwtQMMcWdStG6jjLyMee/+a3yMqqxAPPIb6jlvBmUzxiMGL3ySxxlVBFlFahnjgHahj3wnIKkcu\n6pVvIx45gkis7odeinj4COoVd43CQL3+MnQeQFRtyet94l0fQL181MXgDKYoxEBdvQSpPrh77U+c\nWcTDR1A/+FYAVjlyUek06tUXEQ/d2Q+9FPHgY6gTr6FmbgVgmSMX9fJRxMPvyf+NrZ0wNwepXt9t\ncvhDcYjBD3+AuOeh24qX1kznfrh5AzXY479hjrfoOg2btyDqm+587BJE1VbYvQd1/JUADHNkUZNX\n4OJ5xL0P5/1eIYQW7WPfC8AyRxaVXmDhdz65rvcWhxi89j3EA4+t670ikUA89C7Uq+4mDhJ17Lvr\nWhVkEYff7VxFAaNe+z7ingcR5eXrer948J2oV190rqIg6bkA05PremvsxUBduwIDvXBg5QKmOyHu\nPYxyKaaBodJpvXp78J3rPoc49AicPYGamfHRMkcu6s1jiPseWf8JWjogvQADbpUdFOr4K4h7HlrX\ne+MvBsdfRhx8AFFatv6T7N4Ll8ddVlFQ9HdDxea8Mr2WIioqoaUdzr3po2GOLGrmFlw4A/sOrfsc\nQgjEfY+g3jzmo2WOXNTxVxD3Hl7Xe+MvBid/uK7AcS6ipARx4H63OggIdfI1RIHXCNAFaCdcimkg\nnDkObZ1adAtA3H0/6tQPfTLKkYsaH4Gpa9C2Z13vj7UYqPl5OHsCcWD9TzOL3POQe6IJCHXyh4iD\n9xd8HnHwQdSJV51POgDUiVcRBx8s/ESdB2CwF3XdtRDxG/XmMZ0os8627rEWA3rOw86at7XXXQ/i\nwP1w/iRq1vmk/URdn9Jpv513F36yhmZQCkYGCz+XYxGlFOrEa4h7ChcDUVqmBeHscR8sc+SiThZ2\njWItBurUD/Uk7gOicrOebNw+B76iTr0Oe+5eX9rvEoQQiLsfcO48vxnqh5ISqPVn23Jx4H7tvnX4\nhpqfg64zsPfgus8RczF4HXG3P2IAIO66B3X2hG/ncwCnX0cUkOm1FLH/kLtGPqPOnUDcdQ9C5NUe\nf0XE3VoMnDvPR3ovQE3DmneVW47YioG6Pq3dBe13+XZOsfcgymWr+Io6fwqx9x7/TrjnAHSddo3r\nfESdOwF7fXDjZamph4SAUX+3sy1m1Nk3EXcV9j2KrRjQdRra9qy5s+KaaN+ng1+3bvp3ziJGXR6H\nmzegIf+q45UQVVt147qBi76ds5hR6TScP4nYs373w1KEEIg9d6MunPLtnMWOOuPEYEXU+VOIPQd8\nPacoL9e57F2nfT1vsaLOn4Q9B3xzP2QRe+9GnTvp6zmLlqE+XQOyY6e/5+08AOfdNfIDNTMDfV3Q\nua+g88RXDC6cQviRobIEsfeg80n7xflTiD0BXaNz7hr5gTp3ElFAUHIlxJ4DqPNuZeAL3WegsXVt\nm3atQizFQN26qdMV2zp9P7fYe1A/0ToKRp07GYgY0HkAus6g0i5uUCg6XuC/GFCbhPk51MSo/+cu\nMlTXGURn4V6QWIoBF89B825E2foaaq1K6x5I9aHmZv0/dxGhrl7SDbWSLb6fW2zZBtt2uB44BaKU\ngq4ziM79vp9bCIHodKsDP1DdZxE+JMrEUgy0i8jfeEEWUV4O9U3aR+dYP11noGPfuqsl74To2Idy\nNSGFMT4CJRsQO3YFc/49Lm5QKCq9oItr2wuLF0BcxSCgp5ksov0uN9EUiLp4DrF7b3AD7N6rV4iO\ndaMu+vPEuRKifR/KXaPCGOqHLdvy3nluOWInBiq9oAsw1tmsaU2034XqcmJQCEGLgdi91000hdJ9\n1tc6nbeRbNHdgG+4PkXrRXX5J9ixEwOGB7VSbi5cKVdCtO+Di2ddBeU6UfNz2p8fQIB/kbpGmJ5C\nTV4NboyYo7rPBivYGzZA82798OZYHz4KduzEQF08hwhyVQCwYyckEuAyIdbHQC/sqis4FW41RCKh\nV4cX3QpuPahbN3WFcHN7oOOItr2oi+cDHSPOqO4ziI7C4wUQQzGg5zy0BeiLRmdC0H4XqvtMoOPE\nlcDjBRlEu3MVrZveC9DU5ksDwdVw7rz1o6YmYXpKr4J9IHZiENpE07YHel1G0bq4eC5YX3QGsfsu\nVLebaNaD6jmPCPihCtCB/p7zzuW6HvouQEu7bxl5sRIDdeumTodrag18LNHSgXK+znWhekJw5YF2\nE/V1u+KzdaB6u6C1I/BxxPZq2FCqv7eOvFC9FxAt/l2jWIkBfV26LNvP5nQr0dKhm9a57ph5oaYn\n9dZ8Pi1tV0NUboYtW113zPXQ14VoDTDAn8vuPageFzfIF9XbhfAxCSNWYqB6w7uBxaYK2F4NwwOh\njBcb+ru1LzqgYrOliOZ2/ZTrWDNq6hrcuA676kIZT7R06PvCkR+9Xfqh1CdiJQb0demuoiEhWjpQ\nrhI5L1Rft69L2zvS2uGqxfMl8z0KVbD7nBjkg7p6CRbmobrGt3PGSgxU/0VEwKlwt9Ha4YLI+dLX\nHYFgu4kmH7Rgh/g9ammH/m69d4JjbfRegNYOX9u/x0YM1M0bcGVC9w0KCdHS6VYGeaL6Q14ZNLfD\nwEUXRM4D1dsFLSHFC8hsSLSpgvTYcGhj2k4QLvHYiAEDPZBsQZSUhDdmU5vuYDo/F96YFqOuT8Pk\nNahtCG1MHUTeBiOp0Ma0nr4uRAiZRLfR3MGCqzdYM0Gs3mIjBqq/K9ylLSA2btI+u+HBUMe1lsXg\ncYiCjXMV5YOavAIzt2BnbajjipbdLPS4VO01M3ARmpwYLE/fxcBL55dDNO1Guf1214TqC1+wAe2T\ndu68tdHfowXb561I74Ro6WDBpZeuCXXtCszN6bY4PhIbMVD93eEGj7M0tbpNVNZKf0SC3diGGuwN\nfVwbUYM9iKbd4Q/c0s5CzwVXibwWBnr05l0+C3YsxEDNzMDECCSbQx9brwycGKwFNdCDaGoLf+DG\nVkj1uolmLQzolUHYiC3bobQULk+EPrZtqIGLgXyPYiEGDPVDTTKcyuOlNLXqbBU30ayKmp2BS2NQ\nH3zl8dvYuh0QcO1K+GNbRmSCDZQ07YZUbyRjW8VADwSweouFGOilbWskY+snmnK4PB7J+NaQ6ofa\nhkgEWwihN1JxrqJVUbMzui17FIINlDS3oVJ9kYxtE25lsBqDvdoVEBVNbTq671gRNdiDaIzmiRNA\nNLai3FPn6qT6oS6iFTZQ0rzbCfYdULdu6gfPAHp7xUIMIp9omtpQ/S5usCoDPaF0k12RxlY30dwB\nNXAx0u9Ronm3WxnciVQf1DfrXeJ8xnoxUErpL3lEfk4AXBD5jkQu2I2tqEE30azKQA80R3eNSpKt\nMDbsijhXQQ32IhpbAjm39WLA5QkoLdMl7REhMtkqjuXRgt0XrWDXN8NoCjU/H50NhhO5YJeV6SLO\nEVfEuSKpPki2BnJq+8VgsCfaeAFATT1cvaxTXB1v59IYlJVHK9jl5VC9C0ZdW4rlUErpmEFAE81a\nEckWt4JbBTXUjwgohd4Xx5OU8oPAp9Di8jnP8357yb+/B/gykI2y/rXnef/Zj7H1sqnVj1OtG1FS\novvtjAz42l88Ngz2QkBL23wQyVZ9vySjt8U4rlyC0lJE1ZZo7Whs1U+/juVJ9UFDMPdvwWIgpUwA\nnwHeCwwBx6SUX/Y87+ySQ7/jed6PFzre2xjogUMP+37afBENLahUX7gdOS1BP80YMAE3ZtJLH35P\n1JaYx1A/NIRftLkU0dhC+jvfiNoMI1GTV0GpTN2M//jhJjoMXPA8r8/zvDngC8BTyxwXSLMTNdQf\nWEAlLxqa9BfK8XaMmWhaXbbKCqihPjMEO9nqsr5WItUHDU2B9Y3yQwySQO7ej4OZ15byqJTyDSnl\nV6WU+30YV2cdTIxC7XLDhYtItqBSTgyWQw31IwwQAzfRrMJQv36giZrqGrhxXbc7d9xG0Cts/5NV\nl+c1oNnzvBtSyieAvwH2LHeglPIIcCT7u+d5VFVVLXvShYEeru+sZcuOat8NzpeFPfuZ/sIfrWir\nCZSVlYVun0ovcG10iKo9+3XL7whRlZVcu3WDSqFIbC7cNx7F5xkUUyMpNn3gKTZE+PeUlZWxZetW\npprb2HRlnA119ZHZYiI3xoYo2b2H8jVeIynlszm/HvU87+hqx/shBikg97GvMfPaIp7nTef8/PdS\nyj+QUu7wPO/y0pNlDD6a89IzU1NTyw6cvnAWVZdkpX8PE7WxEjU9yeTYKGJTRdTmLEtVVVXon5Ua\nSUHVVqbn5mEu+utEQzPTZ08h9t5d8Kmi+DyDQKXTpAd7ubFtJyLCvyf7eabrGrl+4TSJCNNcTWSh\nr5vEA48xu4ZrVFVVhed5z+Zzfj/cRMeADilli5SyDPhp4LncA6SUtTk/HwbEckKQN8P9iHoD3A+g\nNw+va3Rxg6UYEi/IoovPeqM2wywujcGmSkTF5qgt0TS26roUxyJKqcC/SwWLged5C8DTwDeAU8AX\nPM87I6X8uJTyFzOH/aSU8qSU8nV0CupPFTouaB+aEX7ODCLZom1yLGJMvCBLstUVCC5laCCS9u8r\nIZKuj9TbuDIBZRsRPrg3V8KXmIHneV8D9i557bM5P/8P4H/4MdZtDA0gnvRFV/yhodnlSC9lqB/u\nfiBqKxYRja2kf/BC1GYYhRrqM0uwG1v03uJKhb7jmrGk+gMXbGsrkN/KJIqm3e5yiGSzWxksIciK\nyXWRbIahflQ6HbUl5jDUr1t8G4KorIKNFdp95QDCEWxrxYDRYdixC1EaTbvdZWlocTGDHNTCAowN\nQ51BrryKzXqiuXIpalOMwThXHrxVIOjQpPoCj71ZKwamxQsA2F4Ns7Oo6cmoLTGDsWHYtkP3BTKJ\nhiYYHrjzcUWASi/oxnD1Zn2XREMzyl2jRdTQQOBFgdaKAQY+zQghXCVyLoZlEmURdY2oETfRADA+\nClXbIq8BeRt1jU6wM6j0gv4sAn74tVYM1HC/cU8z4CqRc9HuB/OuEfVNOoPGkWmJbE68IIuob0IN\nu1bWgI6NVm1FbAy2fslaMWBowKzAZJaGZhhyGUWAuSuDBjfRZDEyXgBasEcGdX59sZMKJ8BvpRiY\nmEmURTS4jKIs5k402gXhJhoymUTmXSNRuRnKyl2gH3Q35BCukZViwOgwVBuWSZQl2Qyp/qKfaNT8\nPIyPBLJxd8FUbdOtgKeuRW1J5Bgr2JBZHTh3XlgrbCvFQA2ZGS8A9EQjBFy7ErUl0TI2BNurEWWG\nZRKRE+gvcleRmp/PpP4aKNiAqG907jxADQ8iQpjvrBQDEzOJsgghtBui2PdxNTRekEUHKIv8qXN8\n2FjBBnR9SpFfI5VO6wer2obAx7JSDNSw4RNNXSOqyMVAux/My1JZxKUu6pWRoasCcIF+QMdMKioD\nzyQCS8WAoQFjVwaA/oKNFPnG60MDeoVkKHqiKW4xUKMpRF30G0OtSL0TbEYHQ9u8yzoxMGl3s5Vw\nOdKZicZgMaDexQwYSRn9PWLrDpifK+qKfjUSnmBbJwaMj8D2nWZmEmWpSxZ1zCBMP+e62b4Tbl5H\n3bwRtSWRYfrKQAiRcecV73eJkZSeT0LAPjEYHjTa/QBA9S6Yvoa6dTNqS6LhygRUbA7Fz7leFjcj\nKmY3RIgTzXop9kC/GhlEhBTXsU4M1GgKYfLSFhCJEqhpgNGhqE2JBtPdDxmKOXVRTU1COq1ToU2m\nvshXBqPhfZesEwMbnmaguDOKwvRzFkQxrwxGBqEuafzmMaK+qWibCqqZWzA9qT0NIWCdGIS5bCqI\nuiKuNRgdtEOwi7hNsv4emX+NinplMDoEu+q1pyEErBIDpZQ1KwPqkkV7E6uRFMLAvlFvo5hTF0N0\nPxTEzlqYuqqfkosMFfI1skoMyKaYBbgptF/o5W1xigGjlgj2rnq4ehk1OxO1JaGjXXnmC7ZIlOjr\nVIx1O8PhekHsEoNhO/ycgFb0sWG9MUURoWZuwVR4fs5CECUlsKuuOAP9tgg2RZxRFPI1skoMdF60\n+U8zgN7qccs2mCiyTb1Hh2BXXWh+zoKpa0QV2VOnmp/X92VNfdSmrI0ideeFnYhhlRhYEy/IUoTF\nZ8qiJ04AUdugA97FxMSo3pu6tCxqS9ZGEa4MlFL6wcrFDJZHjQwaX2OQS1Gml46YXwdyG3XJ4vNH\n2ybYxdjr6+plKC9HVFSGNqRVYsBIyvzq41yKsZR+xI600iyiNll8biJLgseL1DTA+AhqoYjibxF8\nj6wRAzU/B5fHdcDPEooxo8i6iaYuCWNDxbUz3Uh4nTD9QJSXw9btcGk0alNCI4rvkTViwPgI7NiJ\n2GBwg7ql1BdXzGDRz2nTyqCyCko2wOTVqE0JDWsqxHOpLTJ3XgR1IPaIwUjK6I04lqVqG6TTug9M\nMXD1MpSVISo2R21JfhRb3MCymAFk+kgV04NVBBXi1oiBjU8ziy14i+UmtixekEXUJlFFklGkrk/D\n3KzeK8Amim1lEEHmpDViYJufM4uoayyatDib6kBuo5hWBpnvkRWFmzmIuqROWy4C1OyMdltW14Y6\nrjViYPpGHCtSVCsDS/rdLEFPNMVRhWxDC/hlKabMvLEh2FmrK+RDxAoxUEoZv3n3Soi64kldtFaw\ni8kFYVvhZpZtO2B2FnVjOmpLgieia2SFGOgGdQqqtkZtSf7UJXXArhiwdGXArjq4PK7Tl2OObRXi\nWXT8rThEO6r4qB1ikJlkbPNzApmJZiL2E42am9XZRDvD9XP6gdhQCjt2wngR5LHbVgeSgy4QLAJX\n0WgKImgBb4UYWLOhzTK8NdGMRG1KsIwNw84axIYNUVuyPmqTse9RpNIL+j6saYjalPXhVgaBYoUY\nWOvnzFIMwS9bXUQZdHppzIPIE2NQtVVX9NpIXWPsM4r0Bl7RpGhbIQbWBiYzFENanDXbKK5EMTx1\nWhovyCKKYffAa1egtExXxoeMFWJgZfVxLsWQrWLLNoorUAxZX8q2jrJLqWmAidF4bxgV4ffIDjG4\nNKa3vrOU4lgZ2BuYBDIxg3hfI9vdrcWwYVSUnRbsEIPt1YhSixrULSUTM4hrZ0zt57R7omHrdpib\n1e0aYort7lYgs8qOsasowu+RHWJg8xMnZOojVKZeIoZMXYVEAjZvidqSdSOEiP/qYCSalEU/0Q3r\n4nuNotzAywoxsP1pZnGiietNnHmasbIOJIc4b3Sjbt6Am9dhe3XUphRG3FcGEQb5rRADq90PGXSA\nMp43sfWBySxxrhYfSUFtAyJhx1d+JeIcf1Nzc3DlEuyMZgMvK+4MYfnSFoi3C8LylMVFahviuzIY\ntbdw8zbivB9yxIWbVogB9fZPNHH2dcZlZSBivzKw/xqxbQfMzMSzYd1otG367RCDqm1RW1A4tTF+\norE9kyhLbRLGh1HpdNSW+E9MrpGOvzXE8rukhqMt3LRCDGwPTAJQUw+XxlDz81Fb4itqfg4uj1td\nB5JFlG+Eyi3674kZsUgrzSDqYrrKjrhw0xfnlJTyg8Cn0OLyOc/zfnuZYz4NPAFcBz7med4bfoxt\nC6K0VGdyTIzYnyqby/io/XUguWSfOi3svroSKp3WG6bUWtqgbil18cwoUiMpEu96PLLxC14ZSCkT\nwGeAx4EDwEellHctOeYJoN3zvE7g48AfFjqulcQxvXTUzk2HViKW2SpXJqBiM2JjRdSW+EMMr5FS\nKpOIEd13yQ830WHggud5fZ7nzQFfAJ5acsxTwJ8CeJ73MrBVShmfR681ope38XqiiUvweJE4Zn3F\nJXicQcQxo2jqGogEoiq6wk0/xCAJ5O74Pph5bbVjUsscE3/i2BkzJoHJLHFsWKdGUoj6+KzeqGmA\n8ZF4Nawz4Htk3E4kUsojwJHs757nUVUVfjvXIJhr6+TWse9E9veUlZX5PvbUxAgbf/THKI3JNVpo\n38v02PCaPqcgPs8guHF5lJLm3ZQbbuuaP8+qKq5t3U7lzeuUxORBZObqBAuNrVT4eI2klM/m/HrU\n87yjqx3vhxikgOac3xszry09pukOxwCQMfhozkvPTE1NFWykCagt20mn+onq76mqqvJ97IVUPze3\nbOdWXK5R+SbU5FUmJ8Z1dtEqBPF5BsHCQC+JfYeYNdzWfD5PVdPA9YvnEZX29sPKJd3XDdU1vt1P\nVVVVeJ73bD7v8cNNdAzokFK2SCnLgJ8GnltyzHPAzwJIKR8BrnqeVwQbzi5h63aYn0PFpGGdmp6E\n+Tn9d8UEkSjR+1aPDUdtin/ELGYAGXdejDa6MaEFfMFi4HneAvA08A3gFPAFz/POSCk/LqX8xcwx\nfwf0SCm7gM8Cnyh0XBsRQsSrnD4zycSiDiSXGMUN1Mwt3S23elfUpvhLXWO8Av1xiRl4nvc1YO+S\n1z675Pen/RjLdrKpi6JjX9SmFIwajVlgMoOoTeqU2TgwmoKaer3iiRGiLkn61RejNsMX3ircjKZB\nXRYrKpBjRZxqDWLofgBidY1UrK9RTAR7fAR27ERsiLZw04lByMQpdTHKLfqCJE7XiJheI7ZXx6dh\nnSGC7cQgbOLUGTMurauXUpeEsaF4bFMacb+boIhTwzpTHqqcGIRNtmBmwe6CGbWwoJe3NTHpd5OD\nqKyCkg0weTVqUwomdgVnOcRmBWdISxcnBiEjysp1KuaE5Zm1l0Zh63b998SRGDx1Lva7ieHKAIhN\nRpFbGRQzcWhLYYifMyh01pflAcorl6B8I6KiMmpLgiEuW8kakFYKTgwiQdQ1Wj/RmPI0Exhx2Iwo\nrjGdDHFoWKemJiGdNmIDLycGURCH1MXYTzRJ1OhQ1GYUROw6yi4lDg3rRgehzozCTScGESBisLyN\n/UQTg5hB7AW7vByqtsLEWNSmrBuTvkdODKIgBsvbuE801NTD5XFdHWopamQQURt9lkqg2B5ENiRe\nAE4MomHbDpidtbZgRt28AbduwrbqqE0JDLGhFHbs1Nt62opBE01Q2J5eqkYGI29Ql8WJQQRYXzAz\nku13E/Pbx+KWB2p2RtdJxGgv52WxfT9kg1bYMf82m4vNW2CqUXOeZoJE1Dfa+9Q5Ngw7axEl8WpQ\ntxT9PbLzGqn5eR3vqKmP2hTAiUF02FxrEPMag0UsXhkwMlg818jWmMHEKGzbgSgti9oSwIlBdGRa\nWduIGhk0ZmkbJLoexNZrFPM6kCzbq+HWTdSN61Fbkj+jKSPaUGRxYhARwuZaAwN2ZQqFuiQMD9rZ\nsM4gX3SQ6PibnasD0wTbiUFU1NpZMKPSaRgf1vbHnaqt+v8WblNqUv560FibUWSYK8+JQUSI8o12\nFsxcHoeKKsTGTVFbEjh6m1K9OrCJxQZ1Bj11Bkpdo5WxHbcycLyFjXsbFEHuei5Wxg0mr0JJCWLz\nlqgtCQc8hRrfAAAXG0lEQVSbVwb1TVFbsYgTgwgRtfbdxGrUrKeZwLExj90w90PQCAuzvtTUJCws\nwJboG9RlcWIQJTamlxZLWmkGKwV7eBBh0BNn4NQm7Yu/jZjToC6LE4MIERaml5pUPh8K9Rb6o0cG\ntd1FwmLDukvjUZuyZkz8HjkxiBIbA1/DZvk5A2dXHVyesKphnRo2b6IJHNvceYbFC8CJQbRsyxTM\n3LwRtSVrQt24Drdu6EKfIkE3rNul93u2hZEB4yaaoLGtLYV25ZnlbnViECEikdB9SWy5iTOBydg3\nqFuKReml6tZNXRdRvStqU8LFtvjbyKBR1cfgxCBybNoCUz/NmHUDh4FV6aWjKahpQCTi3aBuKTrQ\nb8n3aG5W70+9sy5qU27DiUHU2NSWwsCnmVCwaWVQjPECsGuTm9Eh2FWH2LAhaktuw4lB1Fi0BaYa\nHiiulMUMVq0MiiyTaBGbGtYZ2ujRiUHECJu2wBwu0okmk6liQ8M6NVycqzebGtaZunpzYhA1tQ0w\nPqwbwBmMmpvTfYkM2YgjVDZvAZGAqWtRW3JnRoqs4CwHaxrWGepudWIQMWLjJqio0hOtyYwNQ3WN\nTrUsMhYb1hnuzlMLCzoFtqYIOsouhyXxN2WoYDsxMAELJhqdu27e00xYiLpG7YIxmYlR2LJNV+QW\nI/XmbyWr0mljmz06MTAAGwpmdPC4eMWAhiYYHojaitUZLr5is1yEDTGDKxNQUYnYVBG1JW/DiYEJ\n2NDKejgFdUU80TQ0o4b6ozZjVUzsdxMqtQ0wNmx2wzqDA/xODAzAhs6YaqTIVwb1zRasDIo02yvD\n4oZRBjes0/ECM6+REwMTMLyU/i0/p5k3cSjs2Ak3b6JuTEdtyYoU/coAjP8uuZWBY3V27IIbU6hb\nhjasuzIBFZuN9HOGhRBCP3UPmbk6UEoVX0fZZTC9LYXJgu3EwAB0w7oGXaZuIsPFnUmURdQ3oUx1\nFU1ehZIEoqpItrpciXrDizgNrTEAJwbGYHLcwNSKydBJNhu7MmB4oKgD/FlErbkbRqnr0zAzY2wL\neCcGpmDyjlpFnrKYRa8MzMwoMjkwGSom1+wM9UNDk1FbXebixMAUDE5ddBNNhvomc1cGQwPGuh9C\nxeANo1SqD5FsidqMFXFiYAiioQVSZoqBC0xmqK6B61NmTjRD/UZPNGEhEgldb2Ciy3WoDxqao7Zi\nRZwYmEJtPVweR83ORG3JbajJq5BegC3bojYlckQioZ++DQsiK6Ug1QtODIBsRb95riKV6kcknRg4\n7oDYUJrZAtOwmzjVB8kWY/2cYSMams3LKJq6CgrYuj1qS8ygvlH75w1CKZVZGZgr2E4MDEI0NKMM\ncxVpP2dr1GaYQ4OBcYNUPySbnWBnEMlWVKovajNuxwLBdmJgEskW/SRuEpmVgUMj6puMC/SrVK+O\nOTk0yRbtNjMJCwTbiYFBiKR5GUWmZ0CEjondS1P9TrBz2VkL16eN2gJTpfqMF+yCdmSWUm4Hvgi0\nAL2A9DzvbdtBSSl7gWtAGpjzPO9wIePGlgazVgYqndYuEYODXqGzsxamrqJu3YSqqqitAfREk3jH\ne6M2wxhEIqGzdob6oGN/1OZohvqheXfUVqxKoSuDXwGe9zxvL/AC8KsrHJcGjnied58TglXYWatT\nF015ork0BpWViIrNUVtiDCJRktlRy4xA/6JgG5yyGAUi2YIaNOjBaqgfYfg1KlQMngI+n/n588CH\nVjhO+DBW7BGJRKawyRBXUarX6OyHqBD1zShTgsiXxmBTBaLSCfZtGBR/06m/ZtcYQOETdI3neaMA\nnueNADUrHKeAb0opj0kpf6HAMWONjhsYchOnXCHTsjQYJNhD/dDortFSRLIFZUoQ+fIElG9CbDa7\nieAdYwZSym8CtTkvCfTk/h+XOVytcJrHPM8bllLuQovCGc/zXlxhvCPAkezvnudRZYhvNgxu7d5D\nenyYigD+5rKysrw+y+ujKUrvf5SyIvr818Js+15mX/hK3p9nENyaGEG1drIpBtfIz88zvfcAU5/t\nZ/PmzZFn8Mx1nWKmeTebQ75GUspnc3496nne0dWOv6MYeJ73/lUGG5VS1nqeNyqlrAPGVjjHcOb/\n41LKLwGHgWXFIGNwrtHPTE1N3cnM2KCq60i/8iILAfzNVVVV5PNZLvR2sfC+p5gpos9/LajqWtJ9\n3czOzub1eQZB+uIFOHCI+Rhco3zvz1VJbEAlSpjq70Xs2OnPOddJuuss1CZDvVeqqqrwPO/ZfN5T\nqJvoOeBjmZ9/Dvjy0gOklBVSys2ZnyuBDwAnCxw3viSbjXBBqLk5mBh1PYmWo7oGrl8nPR39BKxS\nvc6VtxKNrWbEDTI1BqZTqBj8NvB+KeU54L3AbwFIKeullF/JHFMLvCilfB14Cfhbz/O+UeC48WXr\nDlhY0D2BomR0EHbWIkpLo7XDQEQiAY0tpAcuRmqHmp+HsWEn2CtgStzAhkwiKLDOwPO8y8D7lnl9\nGHgy83MPcKiQcYoJIYR+ikj1RdocTg32WXEDR4VItrDQdxEaI8wdHxuC7TsRZeXR2WAyyRY4dyJS\nE1R6AUYGdNKB4bh0TwMRyZboK5GH+lyWymo0trLQ3x2pCcoS90NUiGRr9LUG46NQtQ2x0fz9w50Y\nmIgBlch6ZeDEYCVEspWFgZ5ojXDxgtVpaIbRlHanRYVFvb2cGBiISDZH33Ux5VYGq9LYwsJAj64A\njghXB7I6orxc7zc8NhSZDWqw1+g9DHJxYmAiyVZI9Uc20aibN2B6EnbWRTK+DYiKzbqIaGI0OiNS\nvcZXtUZOQ0ukbeHVQA80tUc2fj44MTAQUbkZNlfpTJEoSPVBfZPOmnGsSEnzbhiMxlWkbt6Aa1d0\nnyTHiojGiNtZD1xENLVFN34euG+7qTTt1k8VEaAGexGNrZGMbRMlLR2o/ojSSwd69A50JSXRjG8J\nOr00GperujGtV9g1dqywnRgYimjeDQMRZav0dxvfbtcEStqiEwM10INoctfojjS2aeGMgoFeaGzV\nnW4twImBoYimtuhWBv0XEc12+DmjpKS1E6IqPBtwgr0maup1W/jr4VeLq4FuqwTbiYGpNO+GCJ46\n1fw8DPeDJX7OKEnsqoPZGdTkldDH1oJtz0QTFSKR0PdyFCu4/h6rvkdODExl+07dluJayBPNyADs\nqEGUbwx3XAsRQkBT+KKt5udgNOX2mlgjorkdFUGBoG2uPCcGhiKEiGR1oPqciygfRPPu8OMGQ/1Q\nXavz6B13prkd+sIVg0XBtqTGAJwYGI1oagv/iWbgovNF50MUgu1cRHkhWtrDF+xUH+yqs6pvlBMD\nk2luR/V1hTqk6ut2E00eiOZ2VNhB5H4n2HlR1whXJnRtRkio3i5Ea2do4/mBEwODEW2d0BueGKh0\nWhdRuYlm7dQl4dpVnVMeEqr3ArZNNFEiSkp0f6Aws/N6L0BrR3jj+YATA5PZVQ8zt8ILIo+moGor\notL+LRTDQiRKoLktNNFW83O6otbFdfJCu4pCfLByKwOHnwgh9NNFWBNNz3lsu4FNQLR2onovhDNY\nqg921SM2bgpnvLjQ2gk9IX2PZmZgLKV3WrMIJwaGE+pE03Medu8JZ6w40bYH1RPONVI9FxCWuR9M\nQLTtQfWcC2ewwR6oa0KUloUznk84MTCcMMVATzRODPJFtHZC7/lwBus9r59yHflRl4TpSdTUZOBD\n6ZiOfYLtxMB0Wjuh9wJKqUCHUXOzuvLYBY/zZ2ctzM+jrlwKfCjV24Voc4KdLyJRAi0d4Yh2b5eV\ngu3EwHDE9moo2QCXxoIdqP8i1DValRdtCkIIaNuj3WwBom7dhPERa3bOMg3tKgpeDGyNvTkxsIHd\ne1DdZwMdQi9t3RPnehGtHcG78/q6ddvqDaXBjhNTwhADNTUJk1esqjzO4sTAAkT7Pug+E+wgF89D\nm31PM6Yg2vYGL9hdp/W94FgfbXugJ2CX68Wz0LbHmrbVuTgxsADRsQ/VFawYqO4zbqIphPa90Nel\n6wACQnWfRXS4a7RexLYdUF4Oo8Htiay6ziDa7wrs/EHixMAGmtthbBh1K5hyenV5HGZndMaFY12I\nis2wqy6wPkUqnYbus+DEoCBExwHUhVOBnd/mhyonBhYgSkt1q+SLwfg71YXT0LFPB0Id60Z07A9u\nBTc8AJWbEVu3B3P+YqFzP1w4Hcip1fycfhjYvTeQ8weNEwNLEO13BTfRdJ1GdOwP5tzFRMc+VFdA\nE03XGeci8gHReSCwa8RAj+5UuqkimPMHjBMDSxAd+1ABBZHVhdOIzgOBnLuYEJ37oetMMAHKrjPO\nReQH9Y1w43ogNSHqwmlr4wXgxMAeOvbBxXN6W0ofUdenYWLMFZv5gNixC0rLAglQ6kwit3orFJFI\nQOf+QOIG6twJ2HuP7+cNCycGliA2b9EBSr8rKLvOwO49us2vo2BE537U+RO+nlONj8DMLWho8vW8\nxYoIIG6gFhbgwmnE3rt9PW+YODGwCLHvXtTZN309pzp7HLH3oK/nLGr23Qtn/L5GbyL23esC/D4h\nOu9GnT/p70n7u2HHTsSWbf6eN0ScGFiEuOtelN8Tzek3EPsP+XrOYkYL9nGdCuoXZ45rkXH4Q8tu\nuHbF17iBOvsm4i57XUTgxMAuOvfrwqaZGV9Op65egquXocVtlOIXYscu2LzFt121VDq9uDJw+INI\nlCDuugd1+g3fzqnFwO4VthMDixAbN0FTG/iUGqdOvwF33WNl6bzJiH2H/JtohvpgUwWiusaf8zk0\nB+6D06/7cio1Nwfd52CPvfECcGJgHWLfvf5NNKffQBxwLiK/EfvvRZ3x5xqp08fdqiAAxP77UKff\n8Medd+EkNDTpKnSLcWJgGeLgQ6g3jxV8HqUU6sxxxD4nBr6z5yBcPI+auVXwqdTJ1xAH7vfBKEcu\noto/d546fgxx72EfrIoWJwa20dION2+gRlKFnaf3AlRUInbV+WOXYxFRUam3Dz1VmBtCXZ/WeyS4\nAH8giAP3oU6+VtA5lFKo4684MXCEj0gkEPceRh1/paDzqNd/gLjvUZ+scixFHHoY9cZLBZ1DnXgV\n9h5ElG/0ySpHLuK+R1A//EFhJxnq1/+PwYZDTgwsRBw6jDr+8rrfr5RCvfYDxP1ODIJCHHoE9ear\nBVWMqzdeQhx62EerHLfRuR+uTOiivnWi3ngZcc9DsagBcWJgI3fdA4O9qKlr63v/0ADMz+o9YR2B\nIHbs1BXj62x7oOZm4fRxxD0P+WyZI4tIlOgVXAGrg7i4iMCJgZWI0jLEgftRr31vXe9Xr38fcd+j\nsXiaMZmC3BCnXtdbXFpc0WoD4oHHUD/8/rreq0ZSem9yy4vNsjgxsBTxjveivv9C3u9TSqFe/g7i\ngccCsMqRi3joXahXX9R56HmS/t4/IN7xowFY5biNvQdhbAh1aSzvt6qXvoU4/J7Y9PVyYmAr+w/B\n5QlUNoC1VrINulw75MARu+qgqQ31Rn7xHTV5Fc6dQDz4zoAsc2QRGzYgDr8H9eI383qfSqdRLx1F\nPPojAVkWPk4MLEWUlCAeOZL36kB99+uId33AuYhCQrzz/agXv5HXe9RLRxGHDlu7SYptiHd/EPXi\nN/ML9nedhmxHgJjgxMBixGPvRb30LR1sXAPp6UldIPOocz+EhbjvEejvXnPGilIK9b3nEY+9L2DL\nHFlEshl21kEexZzqha8i3vm+WD1UOTGwGFHfBM3tqBefX9Pxs//wVcTBBxFVWwK2zJFFlJbp+M7z\nz63tDcdfgZIS6/vc2IZ4zwdJf+urazpWjQyizp9EvPMDAVsVLk4MLCfx5E+hvvZXejPuVVA3ppn5\nyhcRT8qQLHNkEY9/BPXyt1GXxlc9TilF+m//L4l//NFYPXHagHjwMbg0hjpz/I7Hqq/9FeJHfkw3\njowRGwp5s5TyJ4FngX3AQ57n/XCF4z4IfAotPp/zPO+3CxnX8RZi916oS6K+9w+I93xwxePU1/+G\n0gceZaHe7ZYVNmLLNsS7H0d99YuIn3165QOPvwxKgSs0Cx2xoRTx4Z8h/VefJ/Frv6O3x1wGNTqE\neuMVEr/5hyFbGDyFrgxOAB8Gvr3SAVLKBPAZ4HHgAPBRKaW9u0YbSOIjP4f68p+v+OSpxoZQ3/57\nNv7Ez4VsmSOLePzDqDdeRnWfXfbf1fVp0l/4YxIf+Vm3KogI8cBjIATqpW8t++8qvUD6Tz6FePKn\nEJVVIVsXPAWJged55zzPuwCsdvceBi54ntfned4c8AXgqULGddyOaGlHvP9DpD/3u6j0wm3/pm7d\nIP2Z30R86J+RcE3pIkNUVpH42adJ/8//+rbKcaUU6c9/GnHoYcTdD0RkoUMkEiR+5hOov/wTVF/X\n2/5dfe2vYUMp4kefjMC64AkjZpAEBnJ+H8y85vAR8fiHoHwj6d//DZ2njt5IPf2Z30R07FvVheQI\nB3HoYcTDR0h/6llUqg8AdeM66nO/B5cnED/xsWgNdCCa20n8zC+R/sxvos6dADIrgi/9Geo7Xyfx\n8/9mRReS7dwxZiCl/CZQm/OSABTwSc/z/jYowxz5IRIlJJ7+ddSX/5z0f/yXsK0aJq8i3v8U4vEP\nO9eDIYgP/TOoriH9O5+Eqq1w7Qri8LtJ/Pv/gigtjdo8ByDuf5QEivTn/zuUbIDpa5BsJfHJ30VU\nbY3avMC4oxh4nvf+AsdIAc05vzdmXlsWKeUR4EjO+DQ0NBRoQhHx9K/o/1agqip+vs4oWdfn+dF/\nrv9zvA1j7s+Gn4QnfzJqKwpCSvlszq9HPc87utrxBWUTLWGlR89jQIeUsgUYBn4a+OhKJ8kYfDT7\nu5QSz/Oe9c3KIkZK+az7LP3DfZ7+4j5P/1jPZ1mQ80tK+SEp5QDwCPAVKeXfZ16vl1J+BcDzvAXg\naeAbwCngC57nnSlkXIfD4XD4S0ErA8/z/gb4m2VeHwaezPn9a8DeQsZyOBwOR3DYEBY/GrUBMeJo\n1AbEjKNRGxAzjkZtQIw4mu8bhFIqADscDofDYRM2rAwcDofDETBODBwOh8Pha2qpr7jmdv4ipewF\nrgFpYM7zvHjs4h0SUsrPoZMiRj3Puyfz2nbgi0AL0AtIz/OurXgSB7DiZ/kM8AtAdv/JX8sknjju\ngJSyEfhTdHFwGvgjz/M+ne/9aeTKwDW3C4Q0cMTzvPucEKyLP0Hfj7n8CvC853l7gReAXw3dKjtZ\n7rME+D3P8+7P/OeEYO3MA//O87wDwKPAL2Xmy7zuTyPFANfcLggE5l5v4/E870XgypKXnwI+n/n5\n88CHQjXKUlb4LGH1hpeOFfA8b8TzvDcyP08DZ9CdHvK6P02dHFxzO/9RwDellMeklL8QtTExocbz\nvFHQX0igJmJ7bOdpKeUbUso/llLGtwlQgEgpW4FDwEtAbT73p6li4PCfxzzPux/4R+hl5DujNiiG\nuDzt9fMHwG7P8w4BI8DvRWyPdUgpNwP/D/g3mRXC0vtx1fvTVDHIq7md485kqsLxPG8c+BLaFeco\njFEpZS2AlLKOt4KfjjzxPG/c87zsZPVHwENR2mMbUsoNaCH4M8/zvpx5Oa/701QxWGxuJ6UsQze3\nW+OO4o6lSCkrMk8NSCkrgQ8AJ6O1ykoEt/u1nwM+lvn554AvL32DY0Vu+ywzk1WWj+Duz3z5X8Bp\nz/N+P+e1vO5PYyuQM6mlv89bqaW/FbFJ1iKlbEOvBhQ6nfjP3eeZH1LKv0C3Vq8GRoFn0H25/hJo\nAvrQqXtXo7LRFlb4LH8E7etOo9MgP571dztWR0r5GPAd9DbEKvPfrwGvAB5rvD+NFQOHw+FwhIep\nbiKHw+FwhIgTA4fD4XA4MXA4HA6HEwOHw+Fw4MTA4XA4HDgxcDgcDgdODBwOh8OBEwOHw+FwAP8f\n2H64yXj7zb8AAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7fad0e930400>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"tck = splrep(x, y)\n",
"x_spl = x + 0.05 # Just to show it wors\n",
"y_spl = splev(x_spl, tck)\n",
"plt.plot(x_spl, y_spl)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### splrep with Weights"
]
},
{
"cell_type": "code",
"execution_count": 26,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7fad0e78cfd0>]"
]
},
"execution_count": 26,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAD/CAYAAAAT87ocAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHxZJREFUeJzt3XucVXW9//HXd80wILARQS4OMIAXQAElVFRUHDQFAyMt\nv0Edy0tJF8tOnWPaqbR+ntIu5rVTkhZZHfqaeUkP3hI0NYVUkKuA3Bnu12EQHWZ9f3+sPToijDOz\n18zaa8/7+XjMY2Ztv3t/P3zde733+q6b8d4jIiKtW5B0ASIikjyFgYiI5H8YWGvLk66hUGgs46Xx\njJfGMz5NGcu8DwOgPOkCCkh50gUUmPKkCygw5UkXUEDKG/uENISBiIg0M4WBiIhgUnBoad4XKCKS\nh0xjGhc3VxVxqqioSLqEgpDJZKisrEy6jIKh8YyXxjM+paWljX6OpolERERhICIiCgMREUFhICIi\nKAxERASFgYiIoDAQEREUBiIigsJARERQGIiICCkJA799a9IliIgUtHSEwex/JF2CiEhBS0cYzHou\n6RJERApaKsKAHdvwG9YmXYWISMFKRRiYk8/Av6ytAxGR5pKOMBhxFn7Ws6TgRjwiIqkUy81trLVj\ngVuJwuUe59zNB2hTDvwCaANsds6NbnAH/Y4GDKxcBv2PiaNkERGpI+ctA2ttANwJjAEGA5OstYP2\na3MocBcw3jk3BLi4MX0YYzCnjMK/PDPXckVE5ADimCYaASx1zq1yzlUD04AJ+7X5DPCAc24dgHNu\nS2M7MSPOwv/reXxYk3PBIiLyfnFME/UC1tRZXksUEHUNANpYa2cAHYHbnXP3NaYT07MXdO4Ki+fB\nccNyKlhERN4vln0GDexnOHA20AH4p7X2n865Zfs3zO5bKK9dds6RyWQA2DvqPMJXX6T9KWe2RM0F\np6Sk5N2xlNxpPOOl8YyXtfaGOosznXMz62sfRxisA8rqLPfOPlbXWmCLc24vsNda+xxwAvCBMMgW\nPLPOQ9dXVlYC4I8fQfjAVPbZKzAlbWMovXXJZDLUjqXkTuMZL41nfDKZDM65GxrznDjCYDZwtLW2\nL7AemAhM2q/Nw8Ad1toioC1wCnBLYzsynbtA36Nh3itw4sgcyxYRkVo570B2ztUAVwFPAguAac65\nRdbaydbaK7NtFgNPAK8DLwF3O+cWNqU/M2IUoY4qEhGJlUnBiVy+oqLivYU9VYTXfoHgx3djOmh+\nsTG0GR4vjWe8NJ7xKS0tBTCNeU4qzkCuy7TvgBkyHD9LVzIVEYlL6sIAwIw8B//PZ5IuQ0SkYKQy\nDDjuBNi+BV+xOulKREQKQirDwARFmFNH41/U1oGISBxSGQYAZuTZ+Jdm4mt0eQoRkVylNwyO6ANd\nDoeFc5IuRUQk9VIbBpDdkfzi35MuQ0Qk9dIdBiefiV/wGr5qd9KliIikWrrDoENHzOCP4GfrnAMR\nkVykOgxAU0UiInFIfRhw3DDYtgW/fs2HtxURkQNKfRiYoiLMqWfpnAMRkRykPgwgO1X00gydcyAi\n0kSFEQalZdC1O8x/JelSRERSqSDCAMCcOYbwuSeSLkNEJJUKJwxOOh2WLcJv25x0KSIiqVM4YdC2\nHWbEKPzzTyddiohI6hRMGACYUWPwzz+FD7UjWUSkMQorDPr0h85dYP6rSZciIpIqBRUGAObM87Qj\nWUSkkQovDE4+E5YuxG/fmnQpIiKpUXhh0O4QzEln4F/QjmQRkYYquDAA7UgWEWmswgyDvkdBx066\nC5qISAMVZBgAmFHakSwi0lCFGwYjRsEb87UjWUSkAQo3DNq1x5wyCv/c40mXIiKS9wo2DADM6HH4\n557AV1cnXYqISF4r7DA4og/06ot/5YWkSxERyWsFHQYAwdnj8DMeS7oMEZG8VvBhwPEnw87t+JVL\nk65ERCRvFXwYmKAIU34+/hltHYiIHEzBhwGAOeNc/NyX8ZU7ky5FRCQvtY4w6NgJ85HT8P94MulS\nRETyUqsIAwBz9jj8s9PxNbpekYjI/orjeBFr7VjgVqJwucc5d/NB2p0MvAh82jn31zj6bihTdhR0\n6QZzX4bhI1uyaxGRvJfzloG1NgDuBMYAg4FJ1tpBB2l3E5DYBYPM6HGE2pEsIvIBcUwTjQCWOudW\nOeeqgWnAhAO0+xrwF2BTDH02iRl+Gmxch1+zIqkSRETyUhxh0AtYU2d5bfaxd1lrS4FPOOf+BzAx\n9NkkprgN5uwL8E8+lFQJIiJ5qaV2IN8KfLvOcnKBcNYY/Ouz8du2JFWCiEjeiWMH8jqgrM5y7+xj\ndZ0ETLPWGuBw4HxrbbVz7pH9X8xaWw6U1y4758hkMjGUmZXJsOesMZjnn+SQz06O73VToKSkJN6x\nbOU0nvHSeMbLWntDncWZzrmZ9bU33vtcOywC3gDOAdYDs4BJzrlFB2n/W+BvjTiayFdUVORU4wde\ncMtGwhu/SfDjKZhD2sf62vksk8lQWVmZdBkFQ+MZL41nfEpLS6GRMzA5TxM552qAq4AngQXANOfc\nImvtZGvtlQd4Sm7pEwNzeA/MccPwzz+VdCkiInkh5y2DFhD7lgGAX7GU8Fc3EfzobkxRUeyvn4/0\nzSteGs94aTzjk8iWQVqZ/sfA4d11rwMREVpxGAAE512If+JBUrB1JCLSrFp1GDD0JHhnLyyZn3Ql\nIiKJatVhYIIAc+4EwiceTLoUEZFEteowADCnjoZVy/BrVyZdiohIYhQGJW0xH52An/6XpEsREUlM\nqw8DILot5sI5+A37nzgtItI6KAwAc0h7zOhx+Me1dSAirZPCIMucMx4/ZxZ+a2JX2BYRSYzCIMt0\nyGDOPA//eIvegE1EJC8oDOow507Az3oOv2Nb0qWIiLQohUEdplNnzGmj8U/p5jci0rooDPZjzrsQ\n//zT+MpdSZciItJiFAb7MV0Ox5x0Ov7pD9x3R0SkYCkMDsCM/ST+uen4qt1JlyIi0iIUBgdguvXE\nDDsV/6SuWSQirYPC4CDM+In4Zx/H79qRdCkiIs1OYXAQpms3zIhR+OkPJF2KiEizUxjUw4yz+H8+\ng9+2JelSRESalcKgHubQwzBnnIt/zCVdiohIs1IYfAgz9iL8Ky/gN29IuhQRkWajMPgQpmMnzNnj\n8H/736RLERFpNgqDBjAfnYCf/yp+/ZqkSxERaRYKgwYw7Ttgzv0E4cN/TLoUEZFmoTBoIHP2OFi2\nGL9iadKliIjETmHQQKZtO8zHJxHefw/e+6TLERGJlcKgEcwZH4U9VTDn5aRLERGJlcKgEUxQRPCp\nywj/8jv8vuqkyxERiY3CoJHMkOHQrQf+2ceTLkVEJDYKgyYILr4c/5jD79ElrkWkMCgMmsD06osZ\ndgr+sfuTLkVEJBYKgyYyEz6Lf/FpXaZCRAqCwqCJzKGHYc65AP/gfUmXIiKSM4VBDsy5F+LfXIRf\nsiDpUkREcqIwyIFp25bg4ssJ//QrfE1N0uWIiDSZwiBXJ54OnTrjZzyWdCUiIk1WHMeLWGvHArcS\nhcs9zrmb9/vvnwG+nV2sBL7snJsXR99JM8YQTJpM+JNr8SedgencJemSREQaLectA2ttANwJjAEG\nA5OstYP2a7YcGOWcOwG4EZiSa7/5xBzRO7oj2gO/S7oUEZEmiWPLYASw1Dm3CsBaOw2YACyubeCc\ne6lO+5eAXjH0m1fMOEt4/VfxS+ZjBgxJuhwRkUaJY59BL6DuXV/WUv/K/gvA9Bj6zSum3SEE9grC\nP/0av29f0uWIiDRKLPsMGspaOxq4DDijnjblQHntsnOOTCbT7LXFwZ81hqoX/k7xi0/TbtzFSZfz\nASUlJakZyzTQeMZL4xkva+0NdRZnOudm1tc+jjBYB5TVWe6dfWz/wo4H7gbGOue2H+zFsgXPrPPQ\n9ZWVlTGU2TK8vZy9N3+bd4achOlyeNLlvE8mkyFNY5nvNJ7x0njGJ5PJ4Jy7oTHPiWOaaDZwtLW2\nr7W2BJgIPFK3gbW2DHgAuMQ592YMfeYt07M3pnxcdO6BboIjIimRcxg452qAq4AngQXANOfcImvt\nZGvtldlm3wO6AL+01r5mrZ2Va7/5zHzsYti0Hl55IelSREQaxKTg26uvqKhIuoZG88sWEf7qZoIf\n3IHpkB/zoNoMj5fGM14az/iUlpYCmMY8R2cgNxNz9LGY4afh77836VJERD6UwqAZmYsuwS96Hb/w\ntaRLERGpl8KgGZl27Qku+Qrh1Dvxe6qSLkdE5KAUBs3MDDkRM2Q43t2TdCkiIgelMGgB5uLL8G/M\nw78+O+lSREQOSGHQAky79gSXfp3wvrvwu3clXY6IyAcoDFqIGTgUc+Lp+D/9WiejiUjeURi0IHPh\n5/BrV+Jfmpl0KSIi76MwaEGmbVuCK/8Df/+9+E3rky5HRORdCoMWZnr3j+598Juf61LXIpI3FAYJ\nMGePhw4Z/N+mJV2KiAigMEiEMYbgsqvxLz6NXzQ36XJERBQGSTGdOhNc/u+E99yC37E16XJEpJVT\nGCTIHHsCpvx8wrt/iq+pSbocEWnFFAYJMx+zUNIW/+B9SZciIq2YwiBhJggIrvgWfvZz+DkvJ12O\niLRSCoM8YDKdCK68hvD3d+LXr026HBFphRQGecIcNQhz0ecI7/pv/J7dSZcjIq2MwiCPBGecizlu\nGOFvbsGH2qEsIi1HYZBnjL0C3nkb/9Afki5FRFoRhUGeMcXFBJOvwf/rBcIX/p50OSLSSigM8pDJ\nHErwte/hH/idzlAWkRahMMhT5og+BJOvIZzyM/z6NUmXIyIFTmGQx8zAoZhPXUp4+w/xO7cnXY6I\nFDCFQZ4LRp6DOf0cwltvwO+pSrocESlQCoMUMOM+jTnmOMK7bsS/83bS5YhIAVIYpIAxBjPxi5jO\nXXVROxFpFgqDlDBBgLnsaqjZh//dbTopTURipTBIEVPchuBL1+G3b8X//k58GCZdkogUCIVBypi2\nbaNzEDatx//hlwoEEYmFwiCFTNt2BF//Pr5iNf6Pv1IgiEjOFAYpZdq1J7j6BnzFKvzUO7RTWURy\nojBIMXNIe4Jv/AC/YyvhlJ/i91UnXZKIpJTCIOVM23YEV30PakLCu36Ef1vnIYhI4ykMCoBp04Zg\n8jWYjp0Ib/kuvnJn0iWJSMoUx/Ei1tqxwK1E4XKPc+7mA7S5HTgfqAIudc7NiaNviZjiYrj8G/DQ\nHwhvuobg6usx3UuTLktEUiLnLQNrbQDcCYwBBgOTrLWD9mtzPnCUc+4YYDLwq1z7lQ8yxhBceAlm\nzIWEP7kO/+bipEsSkZSIY5poBLDUObfKOVcNTAMm7NdmAvB7AOfcy8Ch1toeMfQtBxCMGkvwuasI\n77yR8MVnki5HRFIgjmmiXkDdC+6vJQqI+tqsyz62MYb+5QDM8ScT/MePCO+6kXDtCswnL026JClw\nfu8e2L4NdmzFb98C27fC7krYsxu/ZzdUVcKeKnjnbaipgZp9sC/7EwTsbNMGX1QMxcVQ3AYOaQ8d\nO2E6doLan8O6Yrr1hG49IXMoxpik/9kFI5Z9BpKfTK8ygv/6eXRxu9t/QPiN68EUJV2WpJjf+xZU\nrMZvWAeb18Om9fiNFdHf+/ZB567RCvuwrtm/u0CvMoL2HaFDR2jfEUraRiv82hV/UTH4kI7t2rF7\nx/ZsQFTDW3tg967ogIjdu6ByF6xZTrhpPWzeEAVKtx6Y0jLoexSm7CgoOxLTvmPSw5RKcYTBOqCs\nznLv7GP7t+nzIW0AsNaWA+W1y845MplMDGW2UpkM/r9+xl53L7uvm0z7r32X4kFDk66qIJSUlBTs\ne9OHNYQbKqhZvZya1W8SrllBzarlhDu2UtSrjKLSMoKevSg66XSCI3oT9CjF5PhNvaSkhODQwxrc\nPtxdSbipgprVK6hZuYSaR6dRs3IZwaFdKBowmDZDh1M85ESCrt2aXFOaWWtvqLM40zk3s772xnuf\na4dFwBvAOcB6YBYwyTm3qE6bjwFfdc6Ns9aeCtzqnDu1gV34ioqKnGqUSLtlC6j6n5sxH/04ZsxF\nmEBHFucik8lQWVmZdBk5895HUzorl+BXLMWvWAKrlkGHDPTuh+ndD3plf3c/AlPUPFuXcYynD2tg\nwzr80oWwaC7+jdehQyfMsSdghp4Ix56AaVMSU8X5q7S0FKBRyZxzGMC7h5bexnuHlt5krZ0MeOfc\n3dk2dwJjiQ4tvcw592oDX15hEJNMJsOuVcsJp/wMiooJLvsGppV+a4pDWsPA79sHa5bjlyzAL1sE\nK5ZE8/f9B2D6D8D0Pwb6HRPN1beg5hhPH4awdgV+0Vz867NhzUrMkOGYE0fC8ScXbDAkFgbNTGEQ\nk9oPm6+pwT/xV/xTD2Muvhxz2mjtiGuCtISBf+dtWLEUv3Q+fsmCaOXftTtmwGA4+jhM/wFweI/E\n3wMtMZ5+13b8nFn4fz0Pq5djThyJOW00HHVs4v/+OCkMpF77f9j86uWE9/4CuvUkmDQZ0+XwBKtL\nn3wNA19dDcsX4xfOxS+ZB2tWQGkZ5pjB2QA4FtMh//Z1tPR4+m2b8S8/i3/x71BUjCn/GObUcky7\nQ1qshuaiMJB6HejD5qur8dPvx894DDN+Emb0+ZhARxw1RL6EwfumQhbOhTcXQ2mfaJ584FA4cmAq\nVnBJjaf3Hha/TjjjMViyADPybMy5n4iOiEophYHUq74Pm1+/hvC+u2DfPoKJX8QcObCFq0ufJMPA\nb96AXzQHFs7FvzEPOmailf+xw2DgkFQeXpkP4eq3bsY//TD+xWeiKaQxF2F6pO+yLgoDqdeHfdh8\nGOL/OQP/0H2YQcdjLvycpo7q0ZIrL1+5E7/49egImUVzofodzLEnwKATMMcej+mS/gMB8iEMavnK\nXfgZj+Jn/B9m2CmY8RNTdbCFwkDq1dAPm9+7Bz/9Afyzj2POGos570JMh/R902xuzbny8u+8DcsW\n4hfOiVb+mzfAMYPf+/Zf2qegdnhCfoVBLV+1G//kg9Fn4bTRmHG2xY+yagqFgdSrsR82v3UT/tE/\n4+e8hDnn45iPXoBp174ZK0yXOFdePqyB1cuz8/5zYMVS6NMvWvkfNwz6DYiuTFvA8jEMavmd26PP\nwisvYC6YiBk1ttnOt4iDwkDq1dQPm99YgX/kf/ELX4u2FM4ej+nUuRkqTJdcV1618/5+4RxYPA86\ndX5v5T9gCOaQ1hW8+RwGtfzalYTTpsDuXdEReAOHJF3SASkMpF45r7w2VeCffAg/+3nMyWdgRo/H\n9Cr78CcWqEZvae3cjl8yHxa/Hk39vPN2NO+fnfpJ89ErcUhDGED26KNXXySc9hvM0BMxn7o073bY\nKwykXnF92Pyu7fgZ0/HPPwndemLOOh8zfCSmTZsYqkyPD90hv31rtPJfMj/6vWtHNO8/cGj07b+0\nrODm/XORljCo5fdU4f86FT93FsGkKzHDRyZd0rsUBlKvuD9sft8+mDuL8Nnp753NecpZ0VmtreC6\nR3XH03sfXcHzzcXRjt835kHV7uzKfwhmwBDo3VfncNQjbWFQyy9ZQDj1DsxRAzETr8S075B0SQoD\nqV+zHv2ydTN+1nP4l2bAW3uiw/GGnRLNfRfgjk+/9y0O2biWPfNfiwJgxRtQ0hbTf+C7AUBpWasI\nxbikNQwA/Nt78e5e/IJXCS7/9+hM7wQpDKReLXLtF+9h/Rr8nJfxc16GjeuiQBg4FDNoKJT2Td0K\n0ldVRkf6rFke/V69HLZuoqjf0YT9jolO0DtyUKuf889VmsOglp87m/C+uzCjxmDG28S2BBUGUq8k\nPmx+x7ZoyuSNedHvPVVROPQ7GtPnSCjrj+nU8GvYNye/e1d0+eON66LfG9bC6uXRdE+fftl6j8SU\nHQmlZXQ6rEvqV175pBDCAKL3fDjlZ1BcTPCFb2Eyh7Z4DQoDqVc+fNj8ts3RlTNXv4lfsyJa2bYp\niVa23Uujq2ke3gMO7w5de0D7DrHsZPVhDVRVwc6tsH0rftsW2L4Ftm3Bb6qItmD27YMevTA9e0HP\nXpievaHPkdFO8gNszeTDeBaSQhpPX1ODf/iP+JdmEky+BnPUoBbtX2Eg9crHD5v3HrZtiS60tnk9\nbNmE37IRtmyErZugujq6XWKHTPanI6a4TXSrxKKi935qaqJbJVZX47O/easqul1i1W7YuwfatYfO\nXaLbMnbp9t4tGruXQs9e0XH+jQiefBzPNCvE8fRzZxNOvR3zyUsJTj+nxfptShgU3p49SRVjDHTt\nBl27HfCd66urYc/u6MbqVZWwpxJfvS+6GUvtTdVraqJAKG4DbdoQZH9zSIcoQDpmoi0MHckjLcyc\ncDLBf/6I8M4bCStWYT75+bx9H2rLoBUpxG9eSdJ4xquQx9NXVRL++ifRHQYn/2ezX9alKVsG6Tqs\nQ0QkhUyHDMHXr8cc1pXwZ9/F79qedEkfoDAQEWkBprgYc8lXMSeMILzp2/iN+TXjoTAQEWkhxhiC\nCyZizv8U4U+vw69cmnRJ71IYiIi0sODM8wj+7SuEt/8Qv2xh0uUACgMRkUSYYacQXPFNwrt+FF3F\nNmEKAxGRhJjBHyH40rWEd/8UP++VRGtRGIiIJMgMHEJw1XcJf3srfsFridWhMBARSZg5ahDBl68j\n/M3Po2t4JUBhICKSB8wxxxFMvobw1z9JZKeywkBEJE+YQcdHO5V/+WP86jdbtG+FgYhIHjGDP0Lw\n2S8T3vH/8JvWt1i/CgMRkTxjThyJGfdpwluvb7FLVygMRETyUFB+PubU0YS3/QD/1p7m76/ZexAR\nkSYxF0zE9BtAePdP8DU1zdqXwkBEJE8ZYzCfmQyhx0+bQnPeckBhICKSx0xREcHka/BvzMM/82iz\n9aMwEBHJc6Z9B4Kvfx8//QH83NnN0ofCQEQkBczhPQi+fC3h1Nvx69fG/voKAxGRlDBHDcJceAnh\nL/879iOMinN5srX2MODPQF9gJWCdczv3a9Mb+D3QAwiBKc6523PpV0SktQrOPI9w1TLCe24h+Mp3\nMEE83+lzfZVrgaedcwOBZ4DrDtBmH/BN59xg4DTgq9baQTn2KyLSapmJX4SqSvyjf47tNXMNgwnA\n1OzfU4FP7N/AObfBOTcn+/duYBHQK8d+RURaLVPchuBL1+L/8SR+fjz3Qcg1DLo75zZCtNIHutfX\n2FrbDxgGvJxjvyIirZo59DCCL3yL8He347dtyfn1PnSfgbX2KaL5/ndrADzw3QM0P+gZEdbajsBf\ngKuzWwgHa1cOlNcuO+fIZDIfVqY0QElJicYyRhrPeGk8m+Ck09i75iKq7/0FHb//C0xR0bv/yVp7\nQ52WM51zM+t7KZPLGW3W2kVAuXNuo7W2JzDDOXfsAdoVA48C051ztzWyG19RUdHkGuU9mUyGysrK\npMsoGBrPeGk8m8aHIeEdP8T07k/wyc8DUFpaCtEX9wbLdZroEeDS7N+fBx4+SLt7gYVNCAIREamH\nCQKCy7+Jn/VsTvdRzjUMbgbOtda+AZwD3ARgrT3CWvto9u/Tgc8CZ1trX7PWvmqtHZtjvyIikmUy\nnaKb4ky9A7+zaZe8zmmaqIVomigm2gyPl8YzXhrP3IUP/xG/fAl9fjoFWniaSERE8oQZPxHTq6xJ\nz83pDGQREckfpqgIY69o0nO1ZSAiIgoDERFRGIiICAoDERFBYSAiIigMREQEhYGIiKAwEBERFAYi\nIoLCQEREUBiIiAgKAxERQWEgIiIoDEREBIWBiIigMBAREVJy28ukCxARSaHCuu2ltfYHRP8o/eT4\no7HUeObzj8Yz9rFslLwPAxERaX4KAxERSUUYzEy6gAIyM+kCCszMpAsoMDOTLqCAzGzsE9KwA1lE\nRJpZGrYMRESkmSkMRESE4qQLOBhr7VjgVqLAusc5d3PCJaWatXYlsBMIgWrn3IhkK0oXa+09wHhg\no3Pu+OxjhwF/BvoCKwHrnNuZWJEpcZCxvB74IrAp2+w7zrnHEyoxVay1vYHfAz2IPt9TnHO3N/b9\nmZdbBtbaALgTGAMMBiZZawclW1XqhUC5c+4jCoIm+S3R+7Gua4GnnXMDgWeA61q8qnQ60FgC3OKc\nG579URA03D7gm865wcBpwFez68tGvT/zMgyAEcBS59wq51w1MA2YkHBNaWfI3//fec859zywfb+H\nJwBTs39PBT7RokWl1EHGEqL3qDSSc26Dc25O9u/dwCKgN418f+bryqEXsKbO8trsY9J0HnjKWjvb\nWvvFpIspEN2dcxsh+kAC3ROuJ+2ustbOsdb+xlp7aNLFpJG1th8wDHgJ6NGY92e+hoHE73Tn3HDg\nY0SbkWckXVAB0nHaTfdL4Ejn3DBgA3BLwvWkjrW2I/AX4OrsFsL+78d635/5GgbrgLI6y72zj0kT\nOefWZ39vBh4kmoqT3Gy01vYAsNb25L2dn9JIzrnNzrnaldUU4OQk60kba20xURDc55x7OPtwo96f\n+RoGs4GjrbV9rbUlwETgkYRrSi1rbfvstwastR2A84D5yVaVSrUXAqv1CHBp9u/PAw/v/wQ5qPeN\nZXZlVesi9P5srHuBhc652+o81qj3Z96egZw9tPQ23ju09KaES0ota21/oq0BT3Q48R81no1jrf0T\nUA50BTYC1wMPAfcDfYBVRIfu7UiqxrQ4yFiOJprrDokOg5xcO98t9bPWng48B8wj+ox74DvALMDR\nwPdn3oaBiIi0nHydJhIRkRakMBAREYWBiIgoDEREBIWBiIigMBARERQGIiKCwkBERID/D9vpwkC7\nuTsjAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7fad0e80fd30>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"tck = splrep(x, y, w=np.ones(len(x_spl)), s=len(x_spl) - np.sqrt(2*len(x_spl)))\n",
"y_spl = splev(x_spl, tck)\n",
"plt.plot(x_spl, y_spl)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"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.4.3+"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment