Skip to content

Instantly share code, notes, and snippets.

@carrotIndustries
Created September 16, 2015 19:47
Show Gist options
  • Save carrotIndustries/095ef01df03cd1e59512 to your computer and use it in GitHub Desktop.
Save carrotIndustries/095ef01df03cd1e59512 to your computer and use it in GitHub Desktop.
ts interpolation
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"%matplotlib inline\n",
"import scipy.interpolate\n",
"import scipy.signal"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"ts, xs = np.loadtxt(\"sinus_test.txt\", delimiter=\";\", unpack=True)"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7f629b055320>]"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEACAYAAABCl1qQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztfXmQnVd15+/2pu5WS91q7ZvlVbZs4wU7xkCWNquHGsyw\nVMEkk+AxWWrCQCphHINTNUgBMkxmCgOTImRmwlQIYGAIFaDY4pmhvUAwmyQsW5atfd/VavUi9Xbn\nj/M+v9dP733rPef77vfOr0ql16+/fr93vnvv+Z1z7vIZay0UCoVCoQhDW95fQKFQKBTFh4qFQqFQ\nKCKhYqFQKBSKSKhYKBQKhSISKhYKhUKhiISKhUKhUCgiwSYWxpgBY8zXjDE7jTHPGWPu5uJSKBQK\nBS86GD/7UwC+Y619hzGmA8BCRi6FQqFQMMJwbMozxvQD2Gqtvdr5hysUCoVCHFxlqKsAnDLG/C9j\nzC+MMf/DGNPLxKVQKBQKZnCJRQeAlwP4jLX25QDGAXyQiUuhUCgUzOCaszgM4LC19qeVn7+GOrEw\nxuihVAqFQpEC1lojzcmSWVhrjwM4ZIzZWHnrdQCebXBdaf99+MMfzv07qH1qXyvaV2bbrM0vxuZc\nDfU+AF80xnQB2APg3zJyKRQKhYIRbGJhrd0O4Fe4Pl+hUCgUctAd3EwYGhrK+yuwQu3zG2W2r8y2\n5QmWfRaxiI2xedbfFAqFwkcYY2DLMsGtUCgUinJBxUKhUCgUkVCxUCgUCkUkVCwUCoVCEQkVC4VC\noVBEQsVCoVAoFJFQsVAoFApFJFQsFAqFQhEJFQuFQqFQRELFQqFQKBSRULFQKBQKRSQKJxazs8DP\nfw5IHRt16RKwbZsMFwCMjwPPPSfHd/Ys8OKLcnwnTgAHD8rxHTxInFJ44QXg/Hk5vh07gIkJOb7t\n24GpKRkua2msz83J8En7lrKhcGLxiU8Ad94JfOMbMnwPPwzcfjvw9NMyfL/3e8BNNwF798rwvfWt\nwMaNwLlz/FzWAq98JXDrrcD0ND/f5CRw/fXAa18r4wCOHye+d72LnwsAnn0WeNnLgPe/X4ZveBi4\n7TbgIx+R4fvSl2isf/azMnyPPCLrW8qGwonFo48CDzwAfOEL/Fxzc8AXv0h8X/oSP9/YGPDtbwO/\n9VvAV7/Kz3fwIDmc++4DvvlNfr6f/ARYsAC44QZyPNz4/vdJnMbGyE5ufP3rwDveAfzoR8CZM/x8\njz4KvPvd1FckxPfznwfe8x7g7/+enwsg++6/n8agBD7/eeB3f1fGt5QRhRKLs2eB3buBD34Q+OEP\n+fl27QIWLgR+//eBxx/n5/vpT4Gbb6Zo/6mn+PmefBK45x7gjW+Use+pp4DXv544pex7wxuA17wG\neOIJfr4nngDe/GYSqCef5Od78kkKLDZskCmV/vCHwPveR+J79Cgv18wM2ffnfw5s3UrlYE6cP0/Z\n/EMPUd/UUlRyFEosfvlL4JZbgGuvpc505Ag/36230r9du/ijt+3bqeT1ilcAP/sZLxdA9t12G3DH\nHfSaG1u3kn133UXCyI1t24jvzjtlnOn27XQ/JfisrfJJ9Jdz54BjxyiYuesuyhI5sXcvsGwZsH49\njfcdO3j5fvYz6ivXXEP39vBhXr4yolBi8fzzwKZNgDHkwLlLC888Q+LU3U2dlnsieNs2smvtWpro\n5p5H2L6d7Nu0Cdi5k38icds2cm6bNpH4csLa+XzciwYmJ4EDB6jEFtxPThw4QFnv8uUy9/PFF2lu\nq70duPFGGouceP55upcAjQnuYGbXLporNEbmfpYRhRKLnTurHeiaa6gkxYnnnqMOBNAA4XY4u3aR\nfcbQ/9wdNrBv8WJgcBDYv5+Pa24O2LOHJoCvvpqywosX+fjOnqXVLStXVp03Z2nhxReBq64Curqo\nr3AHMi+8UB0LN9zA77x376YIHwCuu44/cKoVi2uu4V/wsXs38QDUR1UskqNQYrFrFzUkQA27Zw8v\n3/795AACvn37ePkOHACuvJJeb9zI22Gnp2n1zvr19PO11/Lad+IEsGgR0NsLdHaSnZztV3svV64k\ne0dGZPiuu46cG6c41fJdfz2JByfyEIvasS4tFtz3s4wolFgcPFgdINdey59ZHDwIXHEFvb7iCt79\nAZcu0QqaNWuqfIcO8fEdPUpOtLOzysdp34EDNBEbYMMGXvv276/yGUOiyMl34EC1r/T1UemSc0VU\n7f1ct44yNc4yYq1YSIy9/fspAwXof+7AcM+eqn0bNsjuBSoLCiUWhw/TwAD4G3R8nP4tX17lO3CA\nj+/QIRKK9nb6ef163km2WiEE5MVi3Tp+513Lxy0WBw/O5+O+n/v3VwOn7m6gvx84eZKPb+/eqvNe\nvRo4dYoWmXDhyBGauwOIlzOzsHa+fevW6QR3GhRGLEZHqQbd308/r1nDu3wvcKbG0M/cg7/e2UhE\nwpLOTdp5N7KvzHzcwcXRo9Wst7OTViodP87DZe18sVi5kkqIXDvHz52juaa+PvqZu2+WFYURiyNH\nSPED571iBTUy13LWRpE3Z2YhHXlL21cbCQP80VsjZ8othnlmapwOzlpaNrt6dfW9oPTFgdFRGueL\nF9PPbW003rnE6fjx+batXEkLJKSONSkLCiUWQaQBULmGswPVlrwAYOlSOoOHawVPIIYBuKObQ4fm\n861ZQw6BC8eOVSNTgN++Y8fm9xduvqNH5/OtXcsnhtZe7uA4g4vz5ymbWLiw+t7atXxicfjw/HsJ\nkK1c/bNeCNvb6WfufVxlQ2HEolEH4ixFnTxJEUYAY0icuOrC9XyBOE1O8vGtWlX9eeVK3gP36u1b\ns4Z3MJ44Qe0VYPVqPvusJftq+Vau5OsrIyO0qmzBgup7a9fyjYXjx+f3lYCPq/3qA0OAN5iRtq+s\nKIxYNOtAXAOk3tkAvA61ns8Ymlw/dYqHr965LV1KESTXpGU934oVfLYFfLXixCn0Fy5Q5N3bO5+P\ns6/U2hbwcd3P+sgbkBeL1av5xnoj+7j7ZxnBKhbGmHZjzFZjzLeirq2PhAHeVLHeuQG8YtGIj7PD\n1vO1t9PGPElxGhmhRQuuMT5O0X5t2YQ7K2zUVyT5li/n4zt27PKxxxnI1JfYAN4yVKPMgrO/lBXc\nmcUfAXgOQOT2pdOnaQVGLVasoPc50Mx5S4uFtIPjsO/iRSqnBSvZAKCjAxgY4NmLEGRpwWIIoOrc\nOPYiSPeVZpkFV19p5Lw5xeL06eqS9QCScxYAr31lBZtYGGPWAXgTgP8JwERc3lAsli2TFQvuaFHK\nAUxN0cmhAwPz3+cSi1OnLnfeAF803OheBksjOXZxhwk9xy7uRiVS7jJUo8yCa+w1GuucfJpZuAFn\nZvEIgAcBxIr1iiIWHM50Zoac2ODg/Pe5nOmpU/TZbXWtyyWGje4lwDcgi8DX00MT0BxPzWskhpxl\nqEZ8y5bxZhb1Y33pUr6x3ihT08wiOVjEwhjzLwGctNZuRYysApCNNubmqDxSnwpzicXp0yQUwe7t\nAFzRYpgz5bAvjI/DvkaDP+CTFCeu/tLIvr4+mv8ZH3fPd/YsOetaSGcWy5bxHZ9y5kzjErdmFsnQ\nwfS5rwJwnzHmTQC6ASw2xnzeWvs7tRdt3rz5pdcnTgxh2bKheR/ClVmcPUsbgoJzkwJwlk2aOVOO\no67ziLwbOW/p+8mZOQWH0NUiuJ/BgXgu+ertC5Z2nzo1f2LfBc6cuTzrHRigUub09OXjJCskqwjW\n0nhfsmT++z5lFsPDwxiWePRkBFjEwlr7MICHAcAY8xsA/kO9UABVsZiYAD7+8csHAVcq3MzZLF1K\nHUuKTzoSHhzkObAtD/tqd4sH4BSnV77y8vcHB3n6S1BGrEfg4BrZngVnz14uFm1t9N7p05dPDmdF\nI7EYHKQTG+bmLi+fZsHEBGX0PT3z3/cpsxgaGsLQ0NBLP2/ZsiWX7yG1zyJ0GjBIE+snSIM6putJ\nxNOnL0+7gWqHdY0zZxrzcUVTYfZxiWH94Ad4M8N65wbw2tfIeXPxNYqEAb7gqVEZCuApRU1NkQOv\nXTkH0Oq5RYvcL1CIGnv6eNX4YBcLa+3j1tr7wq5pFGkAFA10dlI67BLnzjUejFyDvxnfkiU84nTu\nnKwzbea8pe3j4mvmTLmCi7D2c83XrEwD8JRqAuddHxgCPMFFs77Z2Un+5cIFt3xlRiF2cDcTC4Cn\nAzVz3osWUdTj+vDCPMSiftksIC+GZeIrwv3k6C+jo+Q0u7ou/x1HJhM21pcudT/J3SyzAOh+crRf\nWeGFWLjusM0GozE8AzJs8I+MuN9IJu1MR0bkxVCSL8w+1/czOCusvsYe8Lm2L8yZDg66LwtJB4bN\nMguALzMsKwovFhwdtpmzCfhcO4BmfB0ddN6Q61S4KJE+Z1lIynkHy1WD47RrwXE/m9kG8NgX5kwH\nBty3X7P5NIBn+WxUZqFiER+FEIuwBuXosNJiEeUApOxbuJBKbK6PYZcUC2vDxdA138gICUWjFTpc\ngYVkJNxsPgaoZr4u0WhzaoDBQfdiESaGKhbJUAixGBlpXBMG6P0yZBaSHTaszMbhcCTLbJOT5Lib\nlWmkskJA9l4CfGUoycwiaqy73hEfJYY6ZxEfhRcLjuhGekAWpezFwTc7S6vVGpVpOjoomxkddcdX\nJOfdCmUojsyi2Vjv7+dZOquZhRsUQizOny9/ZiElTmFlmoDPpX1hZZqAz6V9Yc60t5fO4XJZZitz\nXwGiI2+OzKJ+j0UAjrEeZp9OcCdDIcQiqgP5Hi1KOoDJSSo3NSrTAO7tC7MNcG9fGB9HmS2Mr7+f\nsiaXz+yQnrOIKkNJZhYcfJpZuENhxKKsmUXYahpA1pkC8mLBwdds8APuM6cw+9rbqV1d9k/JrBCQ\nnwCOqiK4nrM4f755IKpzFslQCLEIa1DXYhFVpnGdyUSVaaSdt+sBWaTMApDNLAI+l+0XVWabnXVb\nZhsdlS0LSZehLlygzbaNoJlFMhRCLKImuCXLNP395Xamixf7bV+YMw34XM/JNOubgGx/4dg0GuZM\nA+ft8vwk6SqCioU75C4W1spGN1HOjWPwR5VNJMWi7GIonTlJ30/XmUyYM+3qon8un6EhvRpKxcId\ncheLsTGgu5uWWTZCHmLhcqlnVCTsu3NrdhRGAOn2k86cXPeXqODCdX8ZHW3uTAH3DjVszmLxYvIH\nrhYMXLpE/y9Y0Pj3KhbJkLtYRKX5rlNhaWcTZV8ezs1nZxqnLCS1rwPg6S/NsuyAz6V9YZE34Fbs\nZ2aoDNzX1/j3bW30XVzZd+FC84UlQPVe6jHl8ZC7WIRNbgOUdbS1VQ9YywrpGnRUhy2DM40SQ0nn\n5ntmEcc+abFwFX2fP0/fv9Hx5LV8rsQpyrauLjqq3JVvKTtyF4so5wa47UCjo/Gct6toQ3rwR9nn\n2pnG4ZN0bq6dd5z+4jq4iLLPFZ+1NB/RLNIH3JZqpMd61L0E3PfPMqPlxCKqA3V2UsQxMSHD59p5\nSzqbIvK5vp9jY3J81sbjc+Xcxsepnt9svhBwK75h8xUBXAeGKhbukLtYRJWhALerJOJGG64cQBRf\nXx8Jk6tJvSI6b58zC0m+iQkKVMKct8v7KR15R83HAG4n8OPY53o8lBm5i0WczGLxYnfPfIjbgVwO\nyLA0v62Nfu/q0bFRkalrZyoZCQOy4jQ3R/Xs3t5wPl+dWxy+RYvcjb04Yz2PwFAzi3gohFhERRsu\nGzTKuQHyA1LS4SxcSEsKXT06tmiZjEu+sTESima77wM+XyP9PDKLOGUoV3MkKhZukbtYxKljuoxu\nihi9STocY9wvT5SK9ONMyEo70zwCC1/FIk7JWbqKoGIRH7mLRdwylDqA+HxhzhSQtW/BAnLywQap\nLJiYiJ6QlS4L5VGG8lUs4lQRJAMZwH1ZtszIXSziRhuSzjSP6M3HuvfsbHRN3xh3AzKJc3Ox9Fna\n2cQpkeYhhq7si1qGHPBJZxY6wR0PuYtFnA5U9jJUHpmTC76xMZoDCavpA+4GZBzbgo1WLpY+F1Ho\nXfeVOM5bsm9Kj3UtQ8VH7mJRVGfqc6S/cGH4dZLOO+CTyiwAt2IYx7mNj7tZ+lzEMpRL5x03c/Ix\ncGoFqFg0gK914biRvitxknbeRRSn9na65y6WPrdCJiOdWUhmTmVHIcQizhyCj3XMuTkqh0jNkSRx\n3i7sixMpAv7al0ScXPFF9ZWFCyl7lMpk8hALl3MkOmfhDoUQC8kOJFkWmpighyxJRfpxnXcekb6W\nvdzwBSezugieirhST5fOFhe5i4VkHXN2lpZwRtX0y1w2AfwuQ0U5m4BPKtJ3zSfZX+JE3gsX0mNc\nZ2ay8xU1MFSxiAc2sTDGrDfG/MAY86wxZocx5v3118zM0E7iZo84DeAq2ghq+mFHJAd8Pg7+uHyS\nkSkgb18efL6WvaL4gk2crsZfnMBQM4tigjOzmAbwx9bamwDcDeC9xphNtRcEkVuU83YVbSRxplIT\nloB8JJyHWJTVmQZ8ZS17AbL9paeHAkgXx9GoWLgFm1hYa49ba7dVXo8B2AlgTe01SQajq84q7Uzj\n8Ek7m74+WfuknVsezlQyuPAxU4uzgRNwl8kkXVyiT8uLhsichTHmSgC3A3i69v24zmbBAmr8rEdG\n5OFMixgJl70MlUcmWub+4qL94i7rBtzcz/HxeItLOjv1aXlxEXLKjhsYY/oAfA3AH1UyjJfwyCOb\nMToKbN4MDA0NYWhoqMlnVDtQs4evx0HSSNHa6BKZKz5JZ+MyEr7iiujr8hDDvXuz88VdXZZHcOFb\nZhH3Xrrii2tbLV9U1pMXhoeHMTw8nPfX4BULY0wngH8A8AVr7T/W//5d79qMfftILKIQlKKWLUv/\nfeJ2oI4OijYuXoyefHfBJ13G8FWc4jqcPOw7c0aWz7eyVxLn7aL9kvKNjgKrVmXj5EJ9IL1ly5Zc\nvgfnaigD4G8BPGet/WSja9KofxYUtcP29VHanLVuKi0WSSJvSecmHennIU5Z+WZmgKmpeNG0iznD\nuCVnwE0ZMelYd/XwsTKDc87i1QD+DYB7jDFbK//urb0gjfpnQRI+Fw4nLl97O5XXsh5+V9Syia+Z\njKQ4xZ2QDfiy2hd3JSIgP/ZciFOcA0oDuOqfZQdbGcpa+xQixEi6AyWpm7rKLNauTcYXtWEwii+J\n85aak/FZnOKunsvqvMfHo5/KV8u3Z082PumsPunYkw4MNbOIRq47uJOkpmUuQwGymUywAuTiRRm+\nPCL9rHwzM7T6Lk6Zpsx9Bchnwll6zkIzi2jkKhZFjjZcOLik9rkqLcTlkxqQwamsPs3JjI0lK9NI\nOzcXfSVumaYVAjXNLKKRe2YhHW3EdabS0Zs0n+SA7OigOZksa9nn5uI9qwMof9v52lfKXEVoBXgl\nFr51oDyiRWmHIyW+Y2PJavpZM5k82i5JVijdVyT5NLMoJrwRizyct/I1R5KllwFflgGZxLaODnq8\napZMJo/FEGXOLIq8KU8zi3jwSix8m7NIGnn7FL0lWXrpki8usjrUNM7bt0xGMwuCZhbxkLtYFHFC\nFvAvepOO9JNEioCs8wZkMxkXq8uK3Ddd7uuIgzwCQ80souHVaiifnHdSPhc1felIv6jOG8hHnCQz\nmaxzMnEefBRAeixIL511tQ+o7Mg9syhq3TQr3+wsRZpxN9mV3Zn6mFkUtezV2UnzMpKZjPQy8qxj\nPakYahkqGi0lFpJ8cZ/K54qvyJFwwCdZ9vLNviK3n4uzy/IY63H3kWhmEQ/eiEXW1DTuI1wDZI2m\n0kTekpGwtDNthUxG2j6p/uni7LIiz09qZhEPuYrFzAzQ3R3vWh9X7xTZ2eQRCZedr+z9RYpv4UIq\nsc3OyvBpZhEPuYpFkZdeljlSBMqfWbSCfT5lvknKiMZUj4hJC80s3CNXsUjSWbu7sz3IXdrZFL3G\n7sK5lVl807RfmZciZ7mfSY5qCZClv0xNEWfcp2q62CfTCvBGLLI+yL3okb5vkXca51bmzKkVyohp\n7UtyVIsLvsC2uFWLri76blNT6fhaBd6IRXC9lFgEXGmjjTIPfh/4yj6h7lNmkTQLBWTHela+VoGK\nRRNkXcte9Ejft5q3D+LrC5+18plFkcc6oEd+xEHLiEUeHTbppi7NLIrFJz0nI7UJ8OLF6mGLSfiy\nZBZFFwvNLKKR+2qoJMgj2pDiCybjLl2S4St7jd2nslDwVD6pCeC0ztQXsUiyezuAZhbR8CqzyLIx\nz4cBIlkXzroCpOjO26c5oKS7/QH5SD9r4FT0wFAzi2h4JRY+daCi8wUnpaZ95kPRnXcW5zY9nWzD\nKOBX5J1HICNtX9yjPgJoZhEN78Qi7dHFPkQ3kmUvQNbBuTpVNy5c2JYk0lfn3Rw+THBrZhEN78RC\nnbdbPsk5mbm5dGvZg5p+3Gd1APLOLY++ouLkjk+P/IiGioXypULSSN+Y9A4n6bM6gPILfRbnnWYC\nOKt9Rc/qs5ZJWwG6GiqCzxcHIBlNpanpA+nbL41tCxbQ5H2aTKbsQu9Llu2Lfa0CzSyY+KTrtGn5\n0ohTmpo+kF4M07RdluNh0kTCWVaX+RBYZAmcfJmz0MwiHC0jFj7UodM6AB8i/Tz4JMVpwQISqDT7\nZMoeeftgn05wR6NlxMKHMlRWZ5o00i+z8wbyEcO0czJFHwutMMGtmUU42MTCGHOvMeZ5Y8yLxpiH\nGl2TtEHz2JSXhm9mJl2knzZ6y+K8y55ZlNW+7u7q0x8l+Fphglszi3CwiIUxph3AXwG4F8CNAP61\nMWZT/XVpGjTLPgvJSD/p6p2AL23knXQwZuVT5305JO3LsrpMugwlnTnpcR884Mos7gKw21q731o7\nDeDLAN5Sf5FUtDE1RY9ojPswlABZBn9a511mZyotTnnwlbWs19ND42hmRoZPM4vigUss1gI4VPPz\n4cp78yA9ISsZ6Rd98APp72eaSDELn3Skn8U+H9rPh0wmyyZOnbPgQQfT58ZaQLhly+aXXg8NDWFo\naCj0+qBBrU3m+H0YjHnw5TFHUlZnmpUvTSaax/0cGwMGBpLzJbWvdunz0qXx/85aYHy8XJnF8PAw\nhoeH8/4abGJxBMD6mp/Xg7KLedi8eXOiD+3ooIhjYiLZcc6+OBufavppy2wHDqTjW7EiHZ/Oybjj\nS2Pf3ByN1yxl2SRiMT5OFYv29mRcRc4s6gPpLVu25PI9uMpQPwNwnTHmSmNMF4B3Avimiw9OM0DK\nPvh9KZvkURYq85xTlv6Z9FTWgC9pfxkfT/787QCSY723lx4KNTub/G9bBSxiYa2dAfDvAXwfwHMA\nvmKt3enis9N0oKzOJumuXJ2QbV2+srdfUvvScgV8UmJhDFUrippdFAFcZShYa78L4LuuP1eyA3V1\nUTp76VKyyXhfBn8emZP00uCyi1OastDkZLIyboA0/SVt2wGyY72Wr78/3d+XHbnu4E6DxYuT77WQ\njm58WS1Udj5fyoiS4js2lq0sVNbMIi1fK8E7sZDuQGmjqTR8adeyZ4m8tWziji9NXwn2ACVdRg74\nMRbSBk6AikXR0DJi4UMqHKxll57Al5yT8SHSl8wMA9uS7gHKypcGPmQWaXZv1/LpnEVztIxYlHmA\nZJ2TuXgx2d+ldaY+lIWmpkg8k+72D/iK3leyONM0q6F8CdTS8rUSVCwKyFd0h1PmzMKXewnIZqF5\n8PlwP1sJKhYRaIUVIFLi1NubfE5maopW8KSJ9H1wpr7V9NNMqPtkn4pFc6hYlIAvjwGZRgzTnC+U\ntaZf9rKJZBbqS2aRZsNhWr5WQkuIRVZnWuYBmZQveBJcmkgfSN5+WWxL88wHn5x32cXJh0CtleCl\nWBR9n0WZ+bJwSfOleQ53mdsuK59OcLc2vBOLNE/L8ynST9phXUT6acpCaZH0fvokTkEbJHkOt09i\n4UuJVJfO8sA7sfBhgEh22CyRW8BXZOedxdkA8uKUhs+XyDttZuHT2NPMojlaRiyyDMgkA8Snmj5Q\n7kg/4EvSflnFSTJTS7Pjv+zOW5fO8qH0YmGtbId14WyK7Lx9cqYBX5HvZ9Y5mbSry9LAh8wi6w5u\nFYvmKL1YXLpEgyptpO9DGUOdt/Il4Uu7tDTN8TBjY9mz+iR8WobiQ+nFohUGv/IVh6/MZb3OTnpa\npdQEfkcHcU5OyvCpWITDO7Ho6aF183HXzpc98s5jAjjLhHoefLpgwB2fZNkLSGbf9DTN3/T08HO1\nIrwTi6Qns/o2GIsuTsoXzaftR5ibo8eqpnnQUhq+QOjT7PZPytWK8E4sgOQdSHowSkamWSb0Ar6k\nzi1tzTvgK6pzC/jKbp9UZjExQVF+e3t6vqRjPUvbpX2eTKvAS7FIsjEvawdKOqmXh7OR5FNxCkcr\nLIiQCtSk+dKsLmsleCkWkh1owQLqRFNT8a7P6kzV2eTL50IMpcU+rnO7dCn9Cb4BkjhT38QiKV+r\nofRikXXwJ+XzLRLWTGY+JJ131j1AAV9SZ5q2pg/IzhcCKhZFQunFIuvgT8OXpcP29CQ7KTUPZyqd\nybgoIybhk3Lek5PV5aESfC6caZIScJY9FgFULIqDlhCLrB0oicPJ6ryT1k1dRfpx52TyKNNI8QWR\nvtQCBd8CGYC+b9xTn6Wdt3QVodWgYhGTT8p5B3xS9nV1AW1t8TdalXkOYWKC6vkdHen5il6m8alv\nFp2v1eCtWMSNbvKYs/Apk0nKJzlnYS2t05es6Us6Ux/7ZpEzCxftp4cJNoe3YqEOgCCdOWW1L8mc\nTBDpZ12nL716R5IvSclydDT7WEgiFllPFwCSjz0XY12XzjZGS4iFT5E3UFwxnJsjB56lpp9kTkY6\na/Jt5VxSPlfONIl9OsFdHngpFkk35ZU50pe0b3ycMoO2jL0mLp8LZ5rkOdzSZQwfs9Ail6GkxbDV\nwCIWxpj/YozZaYzZboz5ujGm3+XnFzk1lXQAwdk7WaO3uA7OxeCX5kvyHG5fAwvJsaBi0brgyiz+\nCcBN1tpbAbwA4EMuP7zIA1KytDA2Roe0uYj045aFstoW8MWxz4XwJuFz0VcWLKCJ+Tg7/stcsgTk\n5yy0DMV/leLdAAAWZUlEQVQLFrGw1j5mrZ2r/Pg0gHUuP7+ocxazs7TRKsspm0n4XEX6ks40KZ8r\ncZKaI0mSyUiXoVzwJc0sJOcspDOLe+4Bjh3LxucTJOYsHgDwHZcfWNRoI4j0sxynkITPx8g74Ivj\nvPMoe7kQJ0mxTzIno2WoaCSZc9q2Lds5W74h9fYjY8xjAFY1+NXD1tpvVa75MwBT1tovNfqMzZs3\nv/R6aGgIQ0NDsbjj7rOYmQEuXswe6Sdxbq4i4VOn4vFJOlNpcfJZDKXEqXZ12ZIl/HxFDdSk+ax1\n1z+jMDw8jOHhYX6iCKQWC2vt68N+b4y5H8CbALy22TW1YpEESSL9LA9DScqXR+TtSpziiK9LPuky\nVFxx2rDBDZ9U2Svgu3AhWixc8AXPfJiejj7TysWcxcKFVNqdm4uem5NczBLsAcpyrldc1AfSW7Zs\n4SdtAK7VUPcCeBDAW6y1F11/ftCgUecZlT3yboWykK/3s4h8LpypMfGXrruwr60N6O2N7p8zM3Rk\nTW9vNj7Je+kbuOYs/huAPgCPGWO2GmM+4/LDOzroTKOJifDrWiESLrtz8zVTK2J/ySMTzTrBHfBF\n2Rf0Td+qCD4hw5FpzWGtvY7jc2sRRDdh8xF5OJsyT8i6vJ9x52R8m3BOwufrHFCczMLVHiAgnn15\nLOvWzMITJIk2JLiA8kf60ktZtf3i80lmTnFWRAUrA7Oc6xVAUiziri5TsfAIkmIhHekXuUxTduct\nLYZS/eXSJZrjc7HUM04Z6vx5oN/RuQ2SYz1YXSYlTj6h9GLhokGDXblRz3zwdcJZM5kqn49iKO3c\n4pShzp93xyeZWeTB5wu8Fouo6MbVYAx25UY5HF8jYekae1FXQ0na5zrSlwqcgHhlKOnMQlosWnGC\n22uxkGxQyQEZt27qc6Rf1DKUFF9wL7Ou3onL53IDWdxAzccyVMAXJzDUzMITSEdTcaJFl3VTSQdQ\n1DKNdE3fGHeRvlQWCsQrC7nOLMpehoojhioWnkA62ogzQFw6gP5+GnBhcGVf4CAvRmyfdGXfwEC0\nbda6s0/yXgLUVyT54tjnes6izGUo6fvpC7wVi7jRjasOK+1wBgaAkZHwa1ym+v39cqWFwJmG7cC/\neJGWXbo4TqG/P/peuoyE44iha+cWp69Il6F8nUOI2346Z+EJ4nSgkRF3zjSO8z5/nq5zgTjiNDLi\nji/Kvrk5d6WMzk46YyisVOPyXsYZ/K7bTrKvxOmbeZShfM0s4tzPkZHos7jKhlKLhfSAlHTeAZ+U\nGAa75V1ssgKiHarLe9nXR4fRzczI8MURJ2m+VihDSZb1XNrnC0otFj47gChnOjVF/7Iev17LF2af\nS9uA6PvpUgiDw+/CHByH8IaV2VzyxclkXM/JxBEL6TKUdGbhcjz4AK/FIqrDSg5Ia91GG1HONMia\nXCy9DPiiIn2XkVQUn8usUJqvu5v+D1swIF1mc+nc4kb6kpmFS3GSLgH7Aq/FokhlqMnJ6mm4Eny+\nO2/JMlQefHEyJ1d8QaQ/N9f8Gg6+MEiXoVzOIcSdn9QylCeI6kCzs3TqpdTqJA7nFiezcAVp+ySd\naVw+l4M/jji54mtvp3Jk2HhwLRZxavqSZSjJkvPsLC3O0KWzniCqA124QBObUU/Xios8nFvRIm9f\nnWkcPg7xlRZ7qf45MBCdyUiXoc6dc5dZRPWVYP7HlW/xBd6aG9WByuBMo5yN6zJUkZybz8404CuS\n2J87544vyGTCSlGSZaipKdqB72qxh3QJ2Bd4KxZRa719r+nnURYq0hyJtH2uxVc6uIhzP13uC1iy\npDlfsNjDVZkm6qw014s9pFcG+gJvxaKnhzpPs7XzZYhMJWvsRYuEpeeApNuPY8GAtH3nzjX+3cWL\n5LiDVWFZEZyV1iyTcW1bEIg2K7O5bjtf4K1YRD2kpOxlBenVSXmIr7QY+iy+YXwzM/S8ehePOA2w\nZElzsXA5XxGHz+V8BUBltijfomUozxAVbbhs0CCTafYApDI4myKtFvJdDMP4rJW9n8HubZcTsmHO\nm8OZDg4CZ88253Md6YfdT80sPETYxJfrBjUm3KFylRWa7QIuw5xMkcpCknzBZj1XZZooPg5nGhXp\nDw665RscDBcn1/ZF3U/NLDxDWJ2Wo0ElB2RXFx24NzHR+Pc6oZ4MUZG+6yOn40SmriZkgfD76XIl\nVIAwsTh7lkcswjIL14f6hbWfTnB7iLAOxJEqhg3IPKIbl3x9fSRMzRYM+F6GCruXY2NUZuzocMsn\n2VeknZu0WCxZIluGCms/LUN5iKVLgTNnGv8ujwHpOpOR5GtrC9+ZK+m8g7khV+vmAXlnKl3GiOJz\nHXnnkVmElb18bz8f4LVYhGUWZ8/634GiSguuHUAz+6x1PyB7eujYhEaH7QVr9H0u00iLUx5i2Iwv\njzKU5P3kGHs+wGuxCMsszpwBli1zyxe2tvzsWfo+LrF0afMBcuaMe75m9o2O0mSsi+dTBzCmuQPg\nsu38eRIoCT7ptpN2pmWfs5BuPx9QarFw3aDLljXnO33avTgtW0afW4/gQT4u180DzQcI1+BoZh/H\nvezoaB4tcvAtXUqf22g1G1dfadY3ObLsPOYsJMtQzfomwNN+PsBrsQiLNrjEolEHmpmh6Nt1hw0c\nTj0C21yWaQBZ5x3wNXJw0uLEwRdkYo2WdnNlMmHObflyt3xFyiwkAzVAMwvnMMZ8wBgzZ4xx3G2q\naJZZWCsrFkFk4+qRo1F8ZXCmYXyc4lQE+zhKpIsW0YF6k5OX/47jfpZdLJqJr7U8JWcfwCIWxpj1\nAF4P4ADH5wdoJhajozSB6upBRAHK7Gyi+MogTs0cwOnT8mLoms+Y5pkal1iMjDQus0mvhpLMLIL5\nO9e+xQdwZRafAPCnTJ/9EiQnSAGNhCX5ymBfVBnRNSTForMT6O1tvHqOc59FvThNTtJSa9fzd9J9\n0wc4FwtjzFsAHLbW/tL1Z9ejWWYhHZnmIU5cfKdOXf6+iq97Pkn7uO7nihXAyZPz35ub41md1NND\n2VN9me3MGZqP4Zq/qxcnFYuEMMY8Zox5psG/+wB8CMCHay938k0boK+Poor6tfplKCuE8ZVJnIpg\nXx6RvpR9c3N8NfZGYjEyQuPS5W74AI2Cw1OneISwt5cEqP64nVZdCQUAqZrUWvv6Ru8bY24GcBWA\n7Yakfh2Anxtj7rLWnqy/fvPmzS+9HhoawtDQUKLvUbtWf82a6vtcg7+/n57rPT1NaXgtn3QkvH69\ne77ly4vhvDmd6c6dsnx5Zxbnz9NO+Nr+6gqNxOLkSWDlSvdcAd+JE/P7PqfzDu5n7UkCeWQWw8PD\nGB4eliVtAKf6b63dAeClrmKM2QfgDmttw3UMtWKRFkFpqF4sODpQWxuJ05kzwKpV1fe5nM2SJdWN\nZLUrrc6cAW67zT1fEZxb2fi2bp3/3sQERfu9vTx89fZxOtNGYnHiBJ9YrFpFn18LCbHYsGE+n7RY\n1AfSW7Zskf0CFXDvs2hywLY7rFoFHD8+/z2u1BRoPCC5oo2ODloSWb+RjMu+ZhvJuPjyyCzqyxjT\n03SQIMfBcGHO23WNHWhcpslDLFas4OFbuTIfsZDiKzpYxcJae3WzrMIVVq++XCyOHaP3OdCoAx0/\nPj/T4Objsq/ZRjIuvka2zczQEknXm8iAxgsUjh8nLpcPBgqQR1+pX6DAGQlLZxYrV14+1rnFol58\nOX1L0eH1Dm6AGu7YsfnvcTboihWNxam2DOYSy5dfHk1x2lfPNz3N57wXLqSSzPh49b0TJ2iQut7g\nCFRr3rXgvpf1zpSTr1GWffw471jIuwzFmclIj72io5RiwTlA1q4FjhyZ/97Ro7x8R49Wf56ZoQl9\nrgFSz3fyJA0aDudtDIls7f08epRPeAcGaJfz2Fj1Pc7BX28bN1+jvnnkCN/9zCOzqHfeR46Q3RxY\ns2b+WABULLxGs8yCK9WvH5AzM5Sqcjrvw4erP588yRd5B3y19nHeSwBYt+5yPq7BaIws38AALU6o\nLetxi9PRo/PnnDjFtwhlKE776scewD8eiozSicXsLNUxuTpsvTM9eZJqwhzryhvxcUc2yucOxsjy\nBc8AGR2tvnf0KF/kvXLl5YEaZ1moURmKM7OoDyxmZykw5PItRUfpxOLUKVreyum8a1NTzvkKoHEk\nXJZIH8hfLDhLlo34uO2rL51wlqGWLqWyXm3mdOgQ9SEOrFo1f6xPT1NJVjIw5PQtRUfpxIJzcAD5\nOzfOtLsIfHncT26+2lKGBJ9U+xkDXHEFcKByXOilSxR5c/EtWULRfXAeVbCSjbskG5T1uNuu6PBe\nLBYtomWPwYmU+/YBV17JxxdkFkEH2r9//qYdDr7awb9vH3DVVbJ83Pezlm//flk+6fu5fz8vX21m\nMTVF44KrLATMF4vDh4mfK/I2hvrGvn30M3dg2NtLZ1IFy2e5+0rR4b1YGANcey2wZw/9zN2gPT10\n9k0wsbd3L3D11Xx8a9ZQtB08DpTbvnXrgIMHqz/v28dr3/r1VLoIwH0/16+v2metzP0M7BsZIQfO\nualrwwYSJIBsu+IKvsg74Avu54EDvIETQG0ViMXevfzOu7b9JPiKDO/FAiCx2L2bXnM7NwDYuBHY\ntYte79kDXHMNH1d3N6W+gQPg7rDr1lEdOFheys133XXACy+Q456dJcfD6XA2biQ+gAS/p4cmhiX4\nAmHi2L1dyxf0zd27aWxwYsOGamYhLRYvvkj9hxP17cftW4qM0onFnj386n/DDdUByR0JB3zPP08O\nldt5t7dXHfjFi7RggGu1CUCTpF1dVH8+dIhKJt3dfHzXXkuDfnpapq9s2kRtB8hEptdfP18suJ3p\nNddUnekLL/CL01VX0X0EZMQiGHuATH8pMkojFsEAeeYZ4Oabefmuv5460NQUdSCJDrtzZ3WNOffS\nvWCAPPss2cq9+iPg27EDuPFGXq7ubhK/vXtl+sq6dTQhOzpKfDfdxMt3/fXVTG3XLn7nfcstZBdA\n/7/sZbx8mzZRvwTITu6xVyv2O3bw95cioxRiceedwE9/SmuwJyepTsuJm26igfHcc5RVcJwgWosb\nb6SOum0bcPvtvGUMgAbIM89U+bgR8G3dKsN3441y9rW1kRju2CFj35IlNKe2bx/w858DL385L9/G\njTTRPDYmI4Z33AH84he08urZZ0msOLFpE7XdiROUaXP7liKjFGJx0000Cfytb5FwcDvTV78a+PGP\ngSeeoM7LjV/7NWB4GPjnfyb7JPgef7y8fL/+61U+ifb71V8FfvAD4Omn5e7n975HzpRbLDo6gFtv\nBb7yFQrUOOfvAFocMDgIfPnLFKgtWsTLd+utVB6V8i2FhrU2l39E7Q5vf7u1gLWf/KTTj22Ku++2\ntqPD2q9+lZ9rbs7aDRuI76mn+PkmJqxdtIj4du/m5zt61NrubmsXLLB2ZISfb/t2azs7rV261Nrp\naX6+732P7uWNN/JzWWvt5z5HfK97nQzfxz5mbXu7tb/92zJ8738/jfWHHpLhe9vbyL5PfUqGLwoV\n3ynus0uzF/GjH6Va/nveI8P3yU9SdPPWt/JzGQN89rPAj34EvOpV/Hw9PcDf/A3NkXBHigCt9nrk\nEZro7u/n57vlFmDLFioPSezGfcMbgA98gP6XwG/+JpW87r9fhu+976WVUA8+KMP34INUhvqTP5Hh\n+9jHyLc88IAMX1FhbP2TbqSIjbF5cSsUCoWvMMbAWiteECvFnIVCoVAoeKFioVAoFIpIqFgoFAqF\nIhIqFgqFQqGIhIqFQqFQKCKhYqFQKBSKSKhYKBQKhSISKhYKhUKhiISKhUKhUCgioWKhUCgUikio\nWCgUCoUiEioWCoVCoYgEi1gYY95njNlpjNlhjPnPHBwKhUKhkINzsTDG3APgPgC3WGtvBvBfXXP4\ngOHh4by/AivUPr9RZvvKbFue4Mgs/h2A/2StnQYAa+0pBo7Co+wdVu3zG2W2r8y25QkOsbgOwK8b\nY35sjBk2xgg8SFKhUCgUnEj1nDBjzGMAVjX41Z9VPnOJtfZuY8yvAPgqgKvTf0WFQqFQ5A3nT8oz\nxnwXwMettY9Xft4N4BXW2jN11+lj8hQKhSIF8nhSHscTiP8RwGsAPG6M2Qigq14ogHyMVSgUCkU6\ncIjF5wB8zhjzDIApAL/DwKFQKBQKQTgvQykUCoWihLDWhv4DcC+A5wG8COChJtd8uvL77QBuj/pb\nAIMAHgPwAoB/AjBQ87sPVa5/HsAbat7/GICDAC7UcS8A8JXK3/wYwIYomzyz734ApwBsrfx7wDf7\nAPQA+DaAnQB2gJZWl6b9IuxL3X5FsK3y/vcAbKvY9tcA2srSdhH2pW67ItlX8/tvAngm7diLMrYd\nwG4AVwLorNzQTXXXvAnAdyqvXwHgx1F/C+AvAfxp5fVDoAlxALixcl1n5e92o5r93AVagVXvTP8Q\nwGcqr98J4MsJGtMH+94N4NNJOmnR7AM509+oXNMJ4AkA95al/SLsS9V+RbGt8ru+Gs6vAXhnWdou\nwj7fx15bDdfbAHwRwC9r3kvUflH7LO4CsNtau9/SJrsvA3hL3TX3Afg7ALDWPg1gwBizKuJvX/qb\nyv//qvL6LQAetdZOW2v3Vwx+ReWzf2KtPd7gO9Z+1j8AeG2ETb7ZZyr/0qAQ9llrJ21ldVzls34B\nYG2Dz/Ky/SLsS9t+hbCt8tljAGCM6QTQBWCuwWd52XYR9vk+9u6q2NUH4I8BfLTOnkTtFyUWawEc\nqvn5MKqDIOqaNSF/u9Jae6Ly+gSAlZXXayrXhfE1/Y7W2hkA540xgxF/E/Xd41wjZZ8F8HZjzHZj\nzP82xqyLuD7Od49zDYt9xpgBAG8G8H/r+cvQfg3sS9t+hbLNGPP9yvWjoOh7Hr/vbdfEPt/H3prK\n64+Ajl2aaMYfp/2ixCLu7Hcc9TWNPs9SDhTGE/c7pIEP9n0LVEu8FVSr/LuI65N8dgAR+4wxHQAe\nBfCpSvSTFT7Yl7b9CmWbtfaNAFYD6EayDCLysyOQp32+jz1jjLkNwNXW2m/E5GqKKLE4AmB9zc/r\nMV+9Gl2zrnJNo/ePVF6fqKRbMMasBnAy5LOOIBxHAFxR+awOAP3W2rMRf9PsuxfOPmvt2UoqCgB/\nC+COsOsjvnve9v13ALustZ+u4y9L+11mX4b2K5ptsNZeAvANVEsiZWq7y+wrwdg7DOBuAHcaY/YB\neBLARmPM/6v5m/jtFzahAdqHsQc0YdKF6Emau1GdpGn6t6BJmocqrz+IyydpugBcVfl7U8fXaIL7\nryuv34Vkk2w+2Leq5vVbAfzIR/tA9dKvNbC3FO0XYl+q9iuKbQAWAlhd87lfAfCHZWm7CPtKMfZq\n+DZg/mqoRO0Xx+h/AWAXaMLkQ5X3/gDAH9Rc81eV328H8PKwv628Pwjg/6Dx8q+HK9c/D+CNNe//\nJai+NlP5/z9W3l8AOn8qWP51ZdwG9cS+vwAt6dsGqoNv9M0+UJQzB+BZ1C1DLEP7RdiXuv0KYtsK\nAD+pfP4zAD6F+UtnfW+7MPu8H3t13+dKzF8Nlaj9dFOeQqFQKCKhj1VVKBQKRSRULBQKhUIRCRUL\nhUKhUERCxUKhUCgUkVCxUCgUCkUkVCwUCoVCEQkVC4VCoVBEQsVCoVAoFJH4/yZFMfyRMFw3AAAA\nAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f629cf72e48>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.plot(ts,xs)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7f629b2af978>]"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ4AAAEACAYAAACDEBA8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xuc1VW9//HXWy7mJeVQyV0wwQTNIA3JUictIzpJmqVY\nRwU6cDKPJ33UEX3UL+x0Sj2ZN45KCYolImkpJl6QzmhmgFcEuciYo1wEE7EQb4N8fn+sL7Qdh70H\nZmZ/98y8n4/Hfsz3stb3u/b3McyHtb7roojAzMysXHbJuwBmZta+OPCYmVlZOfCYmVlZOfCYmVlZ\nOfCYmVlZOfCYmVlZlQw8koZLWiZphaTztpPmyuz8QklDSuWV1FXSHEnPSLpPUpeCc+dn6ZdJOq6B\ne82StKhgf1dJt2R55knquyMPwMzMyqto4JHUAZgEDAcGAaMkDayXZgTQPyIGAOOAaxqRdwIwJyIO\nAOZm+0gaBJycpR8OXC1pl4J7nQhsBAoHH40F1mf3vwy4eAefgZmZlVGpGs9QoCYiaiOiDpgBjKyX\n5nhgGkBEzAe6SOpeIu+2PNnPL2fbI4GbI6IuImqBmuw6SNoTOAf4MaCG7g/cBhzbiO9tZmY5KRV4\negErC/ZXZccak6ZnkbzdImJdtr0O6JZt98zSFebpmW3/F/Az4PXt3T8iNgN/k9S1xPcyM7OclAo8\njZ1PR6WToIauF2nOnmL3kaTBwIcj4o5G3svMzCpUxxLnVwN9Cvb78O4aSUNpemdpOjVwfHW2vU5S\n94hYK6kH8FKJaw0DDpP0XFbmfST9ISKOyfLsC6yR1BHYOyJeqf9FJHlSOjOznRARzfsf/ojY7of0\nR/5ZoB/QGXgSGFgvzQhgdrY9DJhXKi9wCXBetj0BuCjbHpSl6wzsl+VXvfv1BRYV7J8JXJNtnwLM\n2M53CWseP/zhD/MuQpvi59m8/DybV/a3s2is2NFP0RpPRGyWdBZwL9ABmBIRSyWNz85PjojZkkZI\nqgE2AaOL5c0ufREwU9JYoBb4WpZniaSZwBJgM3Bm9sUL1W+ymwL8StIKYH0WfMzMrEKVamojIu4G\n7q53bHK9/bMamzc7/grw2e3k+QnwkyLlqQUOKdh/iyxwmZlZ5fPMBbbDqqqq8i5Cm+Ln2bz8PCuf\n3tuS1TZJaqDVzszMipHU7J0LXOMxM7OycuAxM7OycuAxM7OycuAxM7OycuAxM7OycuAxM7OycuAx\nM7OycuAxM7OycuAxM7OycuAxM7OycuAxM7OycuAxM7OycuAxM6twr7wCGzbkXYrm48BjZlah7rwT\nunSBvn3hrrvyLk3z8bIIZmYV6uij4ZvfhG98A9SsCxM0Xkssi+DAY2ZWgWpq4IgjYNUq6Nw5v3J4\nPR4zs3bihhvg61/PN+i0FNd4zMwqzDvvQL9+6b3OIYfkWxbXeMzM2oF774V99sk/6LSUjnkXwMys\nvduwAa6+GpYuhcWL4fnn4Ze/zLtULcdNbWZmOdi8GVauhEWLYPJkePPN9E7nkEPgoINgt93yLmHS\nEk1trvGYmZXZmjWpx9rbb8PgwbD//jBxInzgA3mXrDwceMzMWtjGjXDBBbBwYWpOq6uDsWPh0kvz\nLlk+HHjMzFrYRRfBc8+lWs3AgdC9e34DQiuB3/GYmbWgVavgYx+DJ5+EPn3yLs2Oc3dqM7NWZMsW\nOPtsOPPM1hl0WkrJwCNpuKRlklZIOm87aa7Mzi+UNKRUXkldJc2R9Iyk+yR1KTh3fpZ+maTjCo7f\nI+lJSYslXSNpl+z4GZL+KumJ7DNmZx+GmVlzuv56qK2FCRPyLkllKRp4JHUAJgHDgUHAKEkD66UZ\nAfSPiAHAOOCaRuSdAMyJiAOAudk+kgYBJ2fphwNXS9taQk+KiMERcTDwIeCr2fEAbo6IIdln6s49\nCjOz5hOROg/8/Oewxx55l6aylOpcMBSoiYhaAEkzgJHA0oI0xwPTACJivqQukroD+xXJezxwdJZ/\nGlBNCj4jSUGkDqiVVAMcDsyLiNey63QCOgNbsvzKPmZmudqyBR5/HP7yF3jgAejQIc0wbe9Wqqmt\nF7CyYH9VdqwxaXoWydstItZl2+uAbtl2zyxdg/eTdG+W/u/ArdnhAL6SNfP9RlLvEt/JzKxZvfkm\nPPEEjBkDI0fCzTfDXnvBr3/dvnuvbU+pGk9ju4E15tGqoetFREgqdp8oSPt5SbsCNwHHAvcDdwLT\nI6JO0jhSDerYRpbbzGynTZ4Ml12WprjZf/8068D8+dDb//0tqlTgWQ0U9sXow7trJA2l6Z2l6dTA\n8dXZ9jpJ3SNiraQewEtFrrW6YJ+IeEvSHaRmufsj4pWC01OAS7b3ZSZOnLhtu6qqiqqqqu0lNTN7\njwi47z5YsQKefhp+9SuYMwcOPbTtLF9QXV1NdXV1i96j6DgeSR2B5aQaxBpgATAqIpYWpBkBnBUR\nIyQNAy6PiGHF8kq6BFgfERdLmgB0iYgJWeeC6aR3S71INZr+wO7AXhHxYnbdm4AHIuLqrQEsK8sJ\nwPci4ogGvovH8ZhZk9x5Z5px4CtfgUGD0rQ3hx6ad6laVtnnaouIzZLOAu4FOgBTssAxPjs/OSJm\nSxqRdQTYBIwulje79EXATEljgVrga1meJZJmAkuAzcCZWVPcHsAdWTPbLsAfgGuza50t6fgs/Xrg\njCY/FTOzBlx3XZqFYIwHbTSJZy4wM2uENWvSrNErV8Kee+ZdmvLxzAVmZjmZNg2++tX2FXRaigOP\nmVkJEXDjjXDGGXmXpG3w7NRmZttRVwfPPpuWon77bfjkJ/MuUdvgwGNmVs+aNfDZz6ag06sXDBkC\n117rwaDNxYHHzIw0HmfuXFi+PA0CPeWU1IOtrYzPqSTu1WZm7dKmTbBoUVon58knYeZMOPdcOOCA\nNAvB4MFprrX2ruzjeMzM2qKFC+Goo6B//xRgBg+G6uo05Y21PNd4zKxdiYBjjoGTT4Z/+7e8S1P5\nPI7HzKyJbrsNXnkF/vVf8y5J++WmNjNrNzZsSO9xpk3z+5s8uanNzNqN886DV19NyxlY47REU5sD\nj5m1C2+8AfvuC3/+c+pUYI3jdzxmZjtg7VoYPz4tP923L1RVOehUAtd4zKzN2Do2Z8WK9Jk1K3WV\nPu006NMH9tsPOvrN9g5xU1sTOPCYtW11dTB0aNo+8EAYMCAt1nbiiZ59oCk8gNTMrMD69fDNb8KS\nJfCXv6T51WbP9pxqlc6Bx8xarWnT0oDQ3/0uNaPttlveJbLGcFObmbVKEWmKm6uuSp0GrGW4V5uZ\nWeaxx1JngqOOyrsktqMceMysVbr++rQi6C7+K9bquKnNzFqdN99MC7Q9/ngan2Mtx73azKxdioA/\n/QmWLoWHH4aHHoIjjnDQaa1c4zGzihUBL7+ceq9deSV85jMwbFgKOgcf7Ik+y8E1HjNrF7ZsgTFj\n4Pbb05icffeFW2/9xwBRa90ceMysIvzmN/DEE/DII/Doo2nWgZoa+MAHPCC0rXFTm5nl7q674Nvf\nTr3UDjss1Wz22SfvUhl4rrYmceAxq0xvvZXe11x1FQwfnndprD4PIDWzNufyy9Okng467YdrPGaW\nm8WLU0+1+fPhwx/OuzTWkFxqPJKGS1omaYWk87aT5srs/EJJQ0rlldRV0hxJz0i6T1KXgnPnZ+mX\nSTqu4Pg9kp6UtFjSNZJ2yY7vKumWLM88Se7Zb9ZK/OQnMGGCg057UzTwSOoATAKGA4OAUZIG1ksz\nAugfEQOAccA1jcg7AZgTEQcAc7N9JA0CTs7SDweulrb1ZzkpIgZHxMHAh4CvZsfHAuuz+18GXLwz\nD8LMyuOdd9ISBnfdBXffDWPH5l0iK7dSNZ6hQE1E1EZEHTADGFkvzfHANICImA90kdS9RN5tebKf\nX862RwI3R0RdRNQCNcDh2bVfA5DUCegMbGngWrcBxzbuq5tZOUXA738P73tfWor6oovgf/4HunQp\nndfallLjeHoBKwv2V5EFghJpegE9i+TtFhHrsu11QLdsuycwr4FrASDpXuATwGzg1vr3j4jNkv4m\nqWtEvFLiu5lZC1u+HCZNSstRL1qUBobOmgVf+ELeJbM8lQo8jX0b35gXT2roehERkordJwrSfl7S\nrsBNpJrN/Y0sHwATJ07ctl1VVUWVF/Ewa3YbN6aBnwsXwsUXp1VBzz8fPvYx6NbNg0ErXXV1NdXV\n1S16j1KBZzXQp2C/D6kWUixN7yxNpwaOr86210nqHhFrJfUAXipyrdUF+0TEW5LuIDXL3Z+d3xdY\nI6kjsPf2ajuFgcfMmt8bb6Tpbfr0SWNzxo+Hs86Cjp4jpdWo/5/yCy+8sNnvUeodz6PAAEn9JHUm\nvfifVS/NLOA0AEnDgFezZrRieWcBp2fbpwO3Fxw/RVJnSfsBA4AFkvbIAhRZcPlnYGkD1zqJ1FnB\nzHLwu9/B4YfDU0/B9Onwne846Nh7Ff2VyN6ZnAXcC3QApkTEUknjs/OTI2K2pBGSaoBNwOhiebNL\nXwTMlDQWqAW+luVZImkmsATYDJyZNcXtAdyRNbPtAvwBuDa71hTgV5JWAOuBU5r+WMxsZ0ydCuPG\n5V0Kq3QeQGpmzaK2Ns2ztno17Lpr3qWx5uIpc8ysYt1wA5x6qoOOlebWVzNrkgh44QW47jq48868\nS2OtgQOPme2wN96As8+GJUvg6adht93gi1+EIUNK5zVz4DGzRlm3Li3StmQJzJ4NnTunudYGDvTa\nObZj3LnAzIp6+eVUqzn99DRG57DD0hidkSPT6qDWtnkhuCZw4DFrvNdegxNPTEtR19WlWs2JJ8L3\nvpd3yazcWiLwuKnNzN7j5z+H978/za/maW6subnGY2bvsnZtakp75BHYb7+8S2N58zgeM2tREXDu\nuTBmjIOOtRw3tZnZNnfckWaV/sUv8i6JtWWu8ZjZNpddBhMnwp575l0Sa8v8jsesHaurgxkzYMWK\ntGjbww+nZak7dcq7ZFYp3J26CRx4zJItW+D552Hx4rSMwfz5cNJJ0LcvHHUU9O+fdwmtkrg7tZk1\nyQ9+kGYb6NEDPvpRGDAA7r47DQw1KxcHHrM2bMsWuOqqNB5n8eI03c3zz0Pv3nmXzNozN7WZtWFz\n58IZZ8D3vw+DBqXxOf/0T3mXyloTv+NpAgcea4++8Q0YOjTNJG22Mxx4msCBx9qbV1+Ffv2gpgY+\n+MG8S2OtlWcuMLNGmzEDPvc5Bx2rPA48Zm3U1Klp6huzSuNebWZtxKZNUFubeq49+iisXp1qPGaV\nxoHHrA147DH49KfTeJyBA+GQQ9L4nI7+F24VyL+WZq3UrbfClClpddCXXkpr6Jx5Zt6lMivNgces\nlXjhBVi6FB5/PP387W9T4DnssLSEwS5+Y2uthLtTm7UCzz0HH/84HHpoakY76CA44ojUrGbWkjyO\npwkceKw1O/nkFGz+3//LuyTW3niSULN26IEHYN48uP76vEti1jwceMwq2Msvw6mnpok+d98979KY\nNY+SryMlDZe0TNIKSedtJ82V2fmFkoaUyiupq6Q5kp6RdJ+kLgXnzs/SL5N0XHZsN0l3SVoqabGk\nnxakP0PSXyU9kX08ZM7ajCuvhH/+ZzjxxLxLYtZ8igYeSR2AScBwYBAwStLAemlGAP0jYgAwDrim\nEXknAHMi4gBgbraPpEHAyVn64cDVkra2LV4SEQOBIcCnJA3Pjgdwc0QMyT5Td+5RmFWOiDQAdPJk\n+M538i6NWfMqVeMZCtRERG1E1AEzgJH10hwPTAOIiPlAF0ndS+Tdlif7+eVseyQpiNRFRC1QAxwe\nEW9ExAPZPeqAx4FeWR5lH7NW76mn4BOfSPOrDRoEJ5zgnmvW9pR6x9MLWFmwvwo4vBFpegE9i+Tt\nFhHrsu11QLdsuycwr4FrbZM1y30JuDw7FMBXJB0FPAOcExGrSnwvs4qwfj3MmpUCzuLF8Oc/p1VC\nR4+GffbJu3RmLaNU4Gls/+PG1DjU0PUiIiQVu8+2c5I6AjcDV2Q1IoA7gekRUSdpHKkGdWxDF5o4\nceK27aqqKqqqqhpRbLPmV1eX5lU7+2zYvBmOPTbNq3bdddC3b96ls/asurqa6urqFr1HqcCzGuhT\nsN+HVAsplqZ3lqZTA8dXZ9vrJHWPiLWSegAvFbnW6oL9XwDLI+LKrQci4pWC81OAS7b3ZQoDj1ke\ntmyBo4+G+fOhV680Nue229xjzSpH/f+UX3jhhc1+j1LveB4FBkjqJ6kz6cX/rHppZgGnAUgaBrya\nNaMVyzsLOD3bPh24veD4KZI6S9oPGAAsyK79Y2Av4JzCm2fvk7Y6HlhS8lub5eQPf4ANG+C119Js\nBL//vYOOtT9FazwRsVnSWcC9QAdgSkQslTQ+Oz85ImZLGiGpBtgEjC6WN7v0RcBMSWOBWuBrWZ4l\nkmaSgsdm4MysKa43cAGwFHg86+h2VdaD7WxJx2fp1wNnNMeDMWsJU6fC+PHQuXPeJTHLj6fMMSuT\nDRvSZJ7PPgsf+EDepTFrHC99bdaKzZgBn/+8g46ZA49ZGUSkudZGj867JGb5c1ObWQu56SZYuDCN\n0XnkkdTMNn8+dOiQd8nMGs/LIjSBA4+1tC1bYNWq9A7nwQfTIm3f+lbqMv2JT0CPHnmX0GzHeVkE\nswq1cSN885tQXQ0HHgj9+8PNN8OnPpV3ycwqjwOP2U6IgB/+EBYsSMtQv/RSmnngL3+BPfbIu3Rm\nlc2Bx2wnzJwJv/sdXHQRfOQj0K8fdPS/JrNG8Tsesx302mtpxujp0+HII/MujVnLcueCJnDgseby\nve/BmjWp15pZW+fOBWY5mzsXbrklvdsxs53jGo/ZDjj2WBg7Fk49Ne+SmJWHazxmZbZlSxr8uWIF\nPP44PP00nHRS3qUya90ceMwa8PrrKchcfjk89BAMG5bG58ya5ZmlzZrKgceswM9+BpMnw+rVcMAB\ncMgh8PDDadE2M2seDjzW7j34ICxalOZVmzkTHnggTXPjcTlmLcOzU1u79qc/wRe/mCbyPPjgNOXN\nxz7moGPWktyrzdq10aNh0KA0NsfM3ssDSJvAgcfq+/vfYd99Yfly6NYt79KYVSavQGrWjG65BT7z\nGQcds3Jz4LF264YbYMyYvEth1v74Faq1K5s3p2UM5s9Pg0KHD8+7RGbtjwOPtQtr1qQgs2IF9O2b\nuktfey106pR3yczaHwcea7PuvTfNOvDII/DYY/D1r6eazm675V0ys/bNvdqszdi0CZ57LvVSe/LJ\nVKM580z4+MfTp3dvULP2zTFr+zxJqNl2zJoFp5ySukcfeGBaFfSee+DQQ/MumZnV5xqPtXqbNqVB\noDfckLpHm1nz8TgeswZcfDF88pMOOmathWs81qotWgTHHJPWyunTJ+/SmLU9udR4JA2XtEzSCknn\nbSfNldn5hZKGlMorqaukOZKekXSfpC4F587P0i+TdFx2bDdJd0laKmmxpJ8WpN9V0i1ZnnmS+u7s\nw7DW59//HX7yEwcds9akaOCR1AGYBAwHBgGjJA2sl2YE0D8iBgDjgGsakXcCMCciDgDmZvtIGgSc\nnKUfDlwtbeuHdElEDASGAJ+StHXo31hgfXb/y4CLd+ZBWOvz1FNQUwNnnJF3ScxsR5Sq8QwFaiKi\nNiLqgBnAyHppjgemAUTEfKCLpO4l8m7Lk/38crY9Erg5IuoiohaoAQ6PiDci4oHsHnXA40CvBq51\nG3BsY7+8tT7/939w2mlpRdCjj4ZzzvEgULPWplR36l7AyoL9VcDhjUjTC+hZJG+3iFiXba8Dtk7T\n2BOY18C1tsma5b4EXF7//hGxWdLfJHWNiFdKfDdrBTZsgAUL0rucRYvgN7+BCy+EceNS89q+++Zd\nQjPbUaUCT2PfxjfmxZMaul5EhKRi99l2TlJH4GbgiqxGtEMmTpy4bbuqqoqqqqodvYSV0WuvpYGf\n/fqlJaiPPBLOPTct1GZmLaO6uprq6uoWvUepwLMaKHxt24dUCymWpneWplMDx1dn2+skdY+ItZJ6\nAC8Vudbqgv1fAMsj4sp6998XWJMFpr23V9spDDxWmTZsgLFj0+wDzzwDo0bBjTfmXSqz9qP+f8ov\nvPDCZr9HqXc8jwIDJPWT1Jn04n9WvTSzgNMAJA0DXs2a0YrlnQWcnm2fDtxecPwUSZ0l7QcMABZk\n1/4xsBdwTgP333qtk0idFayVuvFGeP11mD49LdTmoGPW9pQcxyPpC6T3KR2AKRHxU0njASJicpZm\na++1TcDoiHh8e3mz412BmaSaSi3wtYh4NTt3ATAG2Az8R0TcK6k38AKwFHg7K9pVETFV0q7Ar0i9\n3dYDpzTUDOdxPK3D4MHw85+nsTlmlj8vfd0EDjyV74kn4MQT4dlnYRfPqWFWETxljrVpU6emMTkO\nOmZtm2s8VhHeegt69YJHH0292MysMnhZBGszItIibUuWwOLF8OCD8OlPO+iYtQeu8VhZRaRlqKdP\nhyuugOOOS2N0hgxJ43TczGZWWVzjsVYrIs02MGMG7L477L8/3HILfOpTeZfMzMrNgcdazO23w2OP\nwdNPpx5re+8Nzz8PXbvmXTIzy5Ob2qxFzJ4N3/oWjBmTlqEeMgT694cOHfIumZntCI/jaQIHnvJ5\n6y04+GCYNAk+//m8S2NmTeFxPNYq/PjHcNBBDjpm1jDXeKxZLVwIn/tc+tmjR96lMbOmclNbEzjw\nlMeoUTB0aFqgzcxaPweeJnDgaTlvvgnPPZcGgo4bB7W1qQebmbV+HsdjFeXtt+F//xf+8z9hv/3S\niqCXXuqgY2bFOfDYDpk7F6ZMSctQ19SkKW4WLEjdpc3MGsNNbVbSX/8Ky5algaAXXJB6rX3ykzBw\nILzvfXmXzsxakt/xNIEDz85ZsSLNpTZ4cAo0xxwD3/hG3qUys3Jx4GkCB56dc8EFaUDopZfmXRIz\ny4MDTxM48Oy4d96Bvn3hnnvSTARm1v545gIrqzlzoGdPBx0za14OPLZdU6fC6NF5l8LM2ho3tdl7\nbNmSVgb99KfTYNAuXfIukZnlxQNIrcXU1cHZZ6d1c5Yvhz33hO9+10HHzJqfA0879vLL8PDD8Oyz\n8Mc/pvE6l14KAwbAPvvkXToza6vc1NYOvfgiPPoofP/78MEPwqBBqffaySenaW/MzLZyd+omaO+B\nJyJ1FLj//jSp52GHQVUVnHceqFl/pcysLfE7Httpt9yS1sj54x/T/GoONmaWF9d42oHXX4cDD4Sb\nboIjj8y7NGbWmriprQnac+A591xYuxamT8+7JGbW2uQyc4Gk4ZKWSVoh6bztpLkyO79Q0pBSeSV1\nlTRH0jOS7pPUpeDc+Vn6ZZKOKzj+35JekLSx3r3PkPRXSU9knzE7+hDasocfhltvhUmT8i6JmVlS\nNPBI6gBMAoYDg4BRkgbWSzMC6B8RA4BxwDWNyDsBmBMRBwBzs30kDQJOztIPB66Wtr2NuAMY2kAx\nA7g5IoZkn6k78P3bvMsvTwu1de2ad0nMzJJSnQuGAjURUQsgaQYwElhakOZ4YBpARMyX1EVSd2C/\nInmPB47O8k8DqknBZyQpiNQBtZJqgMOBeRGxILtO/TIq+xhpIOgtt6SxOUuXpvnWrrsu71KZmf1D\nqaa2XsDKgv1V2bHGpOlZJG+3iFiXba8DumXbPbN0xe5XXwBfyZr5fiOpd4n0bc4bb8Bjj8G0afCV\nr8BPf5qWMhgxAv70J9hrr7xLaGb2D6VqPI19G9+YGocaul5EhKRi9ylVhjuB6RFRJ2kcqQZ1bCPK\n0yb87Gfwgx+k2QY++lE44gi49to0q7SZWSUqFXhWA4Vj2fvw7hpJQ2l6Z2k6NXB8dba9TlL3iFgr\nqQfwUpFrraaIiHilYHcKcMn20k6cOHHbdlVVFVVVVcUuXbFuvx3uvDNN4LlgASxeDPvvn3epzKwt\nqK6uprq6ukXvUbQ7taSOwHJSDWINsAAYFRFLC9KMAM6KiBGShgGXR8SwYnklXQKsj4iLJU0AukTE\nhKxzwXTSu6VewP2kjgtRcL+NEfH+gv3uEbE22z4B+F5EHNHAd2kT3alffx1694b/+q9Uyxk0KO2b\nmbWEss9cEBGbJZ0F3At0AKZkgWN8dn5yRMyWNCLrCLAJGF0sb3bpi4CZksYCtcDXsjxLJM0ElgCb\ngTO3RossWI0CdpO0EvhlRPwIOFvS8Vn69cAZzfFgKtWtt6bmtG9/O++SmJntHA8gbWWOOgrOOQdO\nOCHvkphZe+CZC5qgLQSe5cvh6KNh5Uro1Cnv0phZe5DLzAVWOaZMgdNPd9Axs9bNs1NXsLffTrWb\nBQvSZ9o0+POf8y6VmVnTOPBUqNdfh4MPhnfegY9/HIYNg7vvho98JO+SmZk1jQNPBXn+efjud9NU\nN8uWwamnwo035l0qM7Pm5c4FOXrlFXjkkRRolixJ86qdcAKcdhp8+MOe6sbM8udebU1QaYGnrg6G\nDIE994RDD03NaoccksboeHVQM6sUXvq6DZk0CXr0gPvuc6Axs/bFNZ4cvPhimtDzoYfSktRmZpXK\n43jagAgYPRrGj3fQMbP2yU1tZXb//Wlszu9/n3dJzMzy4RpPGb39dlo/59xzoaNDvpm1Uw48LeyN\nN+BLX4JevWCPPWDDhjQ+x8ysvfL/u5tZBNxxB8yfD089BfPmpS7SDz0E/fq5B5uZmXu1NZOI1Fvt\nt79NzWljx6ZxOR/9aAo4u7huaWatkMfxVKjrr08Ls+2xR1qC+oYboJWuqm1m1uJc42miiLT89BVX\nwHHHNfvlzcxy5XE8FWj+/DSD9Oc+l3dJzMxaBweeJpo6NQ0IdacBM7PGcVNbE7z+OvTuDYsWpe7S\nZmZtjZvaKsxtt8HhhzvomJntCPdq2wEvvghTpqS1cxYuhGeeSWN2zMys8dzUVsKmTbB8ObzwAlx6\nKXTvDiNGpOWoP/IReN/7WqCwZmYVwgvBNcGOBp4IeOwxGD48rZvTrx8MGAA//jHsvnvLldPMrJJ4\nAGkLe/55+NGPUhPa4sWw995pFoIzzsi7ZGZmbYdrPJkIOPLItBz1iSempag/9KEyFtDMrAK5xtOC\nZs5M73MUVo2FAAAFqElEQVSuuMLzqpmZtSTXeICNG+Ggg+DXv4ajjipzwczMKlgu43gkDZe0TNIK\nSedtJ82V2fmFkoaUyiupq6Q5kp6RdJ+kLgXnzs/SL5N0XMHx/5b0gqSN9e69q6RbsjzzJPXd0Ydw\n/vnw2c866JiZlUPRwCOpAzAJGA4MAkZJGlgvzQigf0QMAMYB1zQi7wRgTkQcAMzN9pE0CDg5Sz8c\nuFraNhnNHcDQBoo5Flif3f8y4OJGf3vg5Zfhppvgkkt2JFf7Vl1dnXcR2hQ/z+bl51n5StV4hgI1\nEVEbEXXADGBkvTTHA9MAImI+0EVS9xJ5t+XJfn452x4J3BwRdRFRC9QAh2fXXhARaxsoY+G1bgOO\nLfGd2LgRHnwQrrsOxo+HE06AD36wVC7byv+wm5efZ/Py86x8pToX9AJWFuyvIgsEJdL0AnoWydst\nItZl2+uAbtl2T2BeA9dqVBkjYrOkv0nqGhGv1E+4eTM8+yz8y7+kXmwHHZSWNPj2t0vcwczMmk2p\nwNPYngeNefGkhq4XESGp2H2arffDXnulmQeOOQZ+8Qv3XjMzy0VEbPcDDAPuKdg/HzivXpprgVMK\n9peRajDbzZul6Z5t9wCWZdsTgAkFee4BDq93v4319u8BhmXbHYG/bue7hD/++OOPPzv+KRYnduZT\nqsbzKDBAUj9gDenF/6h6aWYBZwEzJA0DXo2IdZLWF8k7Czid1BHgdOD2guPTJf2c1IQ2AFhQooxb\nrzUPOInUWeE9mrs7oJmZ7ZyigSd7Z3IWcC/QAZgSEUsljc/OT46I2ZJGSKoBNgGji+XNLn0RMFPS\nWKAW+FqWZ4mkmcASYDNw5tbBN5IuIQWu3SStBH4ZET8CpgC/krQCWA+c0ixPxszMWkS7GUBqZmaV\noc2/Xm/MAFgDSbWSnpL0hKQF2bGdGeh7qKRF2bkr8vgueZA0VdI6SYsKjjXb82uOgdKtyXae50RJ\nq7Lf0SckfaHgnJ/ndkjqI+n/JD0tabGks7Pj+f1+NvdLo0r6kJr4aoB+QCfgSWBg3uWqxA/wHNC1\n3rFLgP/Mts8DLsq2B2XPslP2bGv4R+15ATA0254NDM/7u5Xp+R0JDAEWtcTzA84Ers62TwZm5P2d\nc3iePwTObSCtn2fxZ9kdGJxt7wksBwbm+fvZ1ms8jRkAa/9QvwPGDg30ldQDeH9EbO0QcmNBnjYt\nIv4IbKh3uDmf3w4PlG7NtvM8oeGhG36eRUTE2oh4Mtt+DVhK6ryV2+9nWw882xvcau8VwH2SHpX0\nr9mxYgN9VxXkLRw0XHh8Ne37eTfn83vXQGngb5K6tlC5K9lZSnNCTiloGvLzbKSsl/EQYD45/n62\n9cDjnhON96mIOBT4AvBtSUcWnoxUh/bz3El+fs3iGuDDwGDgReDSfIvTukjak1Qb+Y+IeNdky+X+\n/WzrgWc10Kdgvw/vjtiWiYgXs59/BX5HaqZcl827R1bNfilLXv+59iY919XZduHx1S1b8orWHM9v\nVUGefbNrdQT2jgamhWrLIuKlyADX8Y9Jg/08S5DUiRR0fhURW8dN5vb72dYDz7YBsJI6k156zcq5\nTBVH0u6S3p9t7wEcByziH4Nz4b0DfU+R1FnSfmQDfSNN4vp3SYdLEvAvBXnao+Z4fnc0cK3tDpRu\ny7I/jludQPodBT/PorLvPgVYEhGXF5zK7/cz7x4XZejR8QVSL44a4Py8y1OJH2A/Ui+WJ4HFW58T\n0BW4H3gGuA/oUpDnguyZLgM+X3D8UNIfhBrgyry/Wxmf4c2kGTreJrV1j27O5wfsCswEVpBm6eiX\n93cu8/McQ3qZ/RSwMPsj2c3Ps1HP8tPAluzf9xPZZ3iev58eQGpmZmXV1pvazMyswjjwmJlZWTnw\nmJlZWTnwmJlZWTnwmJlZWTnwmJlZWTnwmJlZWTnwmJlZWf1/5GyafEUGXF8AAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f629c6b9ef0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.plot(ts)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"interp = scipy.interpolate.interp1d(ts, xs)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"ts2 = np.linspace(min(ts), max(ts), len(ts))"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"xs2 = interp(ts2)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7f629b304438>]"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEACAYAAABCl1qQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsvWmMJHl63veLvI+IPOvsrr7mnj1mL3J5iISaECksCJkS\nBQESDcqwaMP+YEu0AQuWyA+aEWhZBA1QBm1KsC3bEC9LXtOUVqJEUZSa2uWSNLk7s9fcR1/VVVmV\nlVdE3kf4Q1R0b01VZf6v6OquzgcYYCo7q96Kyoj/ez3v+1i+77PEEkssscQS8xA7619giSWWWGKJ\nRx9LZ7HEEkssscRCLJ3FEkssscQSC7F0FkssscQSSyzE0lksscQSSyyxEEtnscQSSyyxxEJE5iws\nyypZlvV5y7LesCzrdcuyvjsqW0ssscQSS0SLRIQ/+38EftP3/b9gWVYCyEdoa4klllhiiQhhRTGU\nZ1lWEXjV9/2njP/wJZZYYoklHjqiKkNdA/Yty/o/LMv6qmVZ/6tlWbmIbC2xxBJLLBExonIWCeDT\nwC/6vv9poAv8jYhsLbHEEkssETGi6lncBe76vv9Hh19/ng85C8uylkupllhiiSUU4Pu+9bBtRpJZ\n+L6/C9yxLOu5w5d+EPjWCe87t//9rb/1t878d1he3/L6nsTrO8/X5vtnF2NHyYb6q8CvWJaVAt4D\n/kqEtpZYYokllogQkTkL3/e/BnxnVD9/iSWWWGKJh4flBHdEuH79+ln/CpFieX2PN87z9Z3naztL\nRDJnIWTYsvyzrL8tscQSSzyOsCwL/7w0uJdYYoklljhfWDqLJZZYYoklFmLpLJZYYokllliIpbNY\nYoklllhiIZbOYoklllhiiYVYOoslllhiiSUWYuksllhiiSWWWIils1hiiSWWWGIhls5iiSWWWGKJ\nhVg6iyWWWGKJJRZi6SyWWGKJJZZYiEfOWYzGU37l336V2ezh7I3qdIf8k9/9+kOxBbDbdPnCH7zx\n0Oy9f6/Jv/nquw/N3jc+2OUPXr/z0Oz9/uu3+eYHtYdm77f++G1u1VoPzd5vfPlb1Nu9h2bvH//u\na3T744diazbz+eXf+QqT6eyh2AvOllcf2tly3vDIOYs//3M/z49/8TP89C/904di7wf+u7/JX7zx\nCf7hb/3hQ7H3XT/zn/Ejv/URbnzt/Ydi71M/92f5oS88ywc7zchtTWczPv0/fw/f+0sv0R2MIrfX\n6PT53l99nu/4hT/1UA6Ar7+/y+f+xfN8x9/9schtAfzTL3+LH/3tj/EnfuavPRR7f+83bvCXbnyK\nz/3dv/1Q7P3V/+XX+Mtf+g5+/O/9g4di78//Dz/Pj3/x0w/tbDlveOScxb/d/1We9X6CX/rar0Ru\nazr1eW3yazzr/QS/cONXI7e31+xxO/PPudb5y/zsFz4fub0vv34LN/0GG60f4e/8+j+L3N7/+dt/\nRGyWJT94gV/4wu9Gbu9nf/1fUvK+l0msyz/9/WNCjMbxd/7fX+di+y9Qz36Zd+4eRG7v5/7Vr/GM\n+1d4O/lP6A2ij/Z/8ff+Ec95P8Hvd385clsAn3/rV3nW/Ql+886vPRR7v7P3KzzX/Ql+6bXon/Xz\niEfKWbx3r0E/+y5//z/8aXaSX4zc3m/+0RvEZhn+xg/+57w1iP5w+6V/9/s43Zf4Cx/7Ub5ycCNy\ne7/y73+PC6M/yfdv/RD//oMvRW7vN776RV5M/RAv2X+Kf/Gt6P+ev/P2l/mulR/kmnWdf/wH0d8v\nv7/9RX746T9Dtf89/O+/E729b3a+xH/yXT9GenCF//uLr0Vu7+bs9/iZ/+AnmSVcvvLOdqS2BqMJ\ne9kv8sv/6Su49ldodweR2rtVazLIvcvf/7Gf5l76xrIUpYBHyll84f/7Ok7vJX7gE9fwrQlffWcn\nUnu//Y2vsT79Tn70ez/OIP925KWTL77zNa6mP8OPfvenOEh8LVJbAH989zVeLH+Kz33ss9yZ/XHk\n9r5Z/xqfvvgJvv+Zz/BWK/rre7/3Kt//7Kf4zMZ38lrtK5Hb2+U1/vRLn+IF5zv4/Q+iPbynsxnt\n7Gv8ue/6FFfin+VffSPa67u122GcvcuPfM9HWBl+ls9/Odr75Xe//j6JUZXPvrBFpvcM//wPX4/U\n3ue/9BqF3ie4/olrYE159Z2H1+c6L3iknMUfvf8Wm6kXiMUsSoNP8i+/Gu2B8+qdN3i6+AJlJ0uy\nd4V//dU3I7X3+v4bfHzzRb7r+SvMkh7v3qtHau+97mt8z7VP8qc//QL97NuRNxJ3/Nf4wY99kj/5\nkY/QTET78M9mPq3Ma/zIZz/JZ596kZ3RW5HaO+j0GOVu8bnveIFPXPgI73Wivb7f+9ZNYhOHFy6v\n8FzlRd7cj/be/K3Xvk7O+xjpZILLuRd47U60f89//8ablKcvALARe4nf+Va0z/of33yLC4dnS6H/\ncf7NN74Rqb3ziEfKWby+/ybPlQ9voNQzfP3ue5Ha+8B9k09uBfaq/nN8+c1oWUO7kzf4rqdeJBaz\nyPde4N99/Z1I7bXi7/ADH3+BrdUCsXGJP3j9bmS2JtMZw/w7/OlPv8D1l55mkrlHw42OxfP+ThPf\nmvLxaxv84CdexM28SZTCize+9h6p3lXsbIrve+FF9onYWbzxHsXxcwB85vKL3BlE6yxevfk+K7Fn\nAHhh5QXea0frLF69/QaXsy8CcNV+nrfq0T4Lbx28xVPF5wG4mPoof3Qz2s/vPOKRchZ3+2/yqUvB\nB3q1+BTvHUTLGDqw3uQ7rwXOYjN7jTd2P4jUXjfzFt/7fHB91cRVvvJedPb6wwnT3Dafff4SAIXR\n83z57bcjs/fNmzWsUYGVYo5sOkGyf4kvffNWZPb+8K2bZAZXAfjY1XX82JAPdhuR2fvarVsU/SsA\n/MBLzzLMvR9p3fuNezdZSV4F4PtefI52PLrPDuDN2gdcsq8B8Jkrz1ObROss3m29xYtrwbP3/No1\ntr2bkdrbHrzFSxcCe09VrvLeQXT35nnFI+UsOvGbfPqppwD46IWn2e5H6yyGmdt857NXAbhausbt\nTnSHd8vrM0u6fPLpDQA2s1d57+BmZPZeffce8cEa+UwKgGriMm/eux2ZvT9+5xbZ0ZX7XzuzK3zj\nVnT2vnbrJiUCe7GYRXpwma+8G918x5u7t1hNBfY2KjbWNM0729Exot5r3GTLvgrAZ57dYpLdjrSM\neKdzk+dWA2fx2eev4MWj++wA6uPbfPRi8Pd86cpVDqbRBmot632+85ngbHlu7Sq1wdJZyOKRchaj\n9DafevoiAJ+8co2WFZ2zqDU9/NiQZ7fKALywfo3aKLob9ivvbJPoXyCRCHTWr5Wucbd7Mzp7790k\nN7p6/+vN/CU+aER3mH7jzi3K1gNnUU1c5q3d6A6ct2o3WU9fvf+1M7vEN29Hd303m7e45Dy4vvQw\nWud0r3uTZ1euAlB2slijIq/f2ovM3v7kAz5+KbD3yacvMMvtMhxNI7Pnxe7y0UtbAHz389fopqN7\n9nw/OFs++0Jwtnz04hXa3IzM3nnFI+Ms7u53wJpyabUIwMevXGSYuheZva+8vU1ycJFYLDi8X7p8\nhY4VXbTxtQ/ukp9u3f/6+fUr1Mc3I7P3re1bVOIPDrcr5S3uedH1LN7Zu8VG5oG9C/nLfNCM7u95\nu32LK8Wr979eSV3inVp017fbv8Uzqw+ur+hf5pt3onOGB7ObfGzr6v2vM6NLvPp+dM6pl7rJJ64E\nmUU+kyI2qPL1D3YjsTWb+Yyzd/nkYWD40rVN/FSTTi8a+uytWgdiMy5Wg7PlO565Sj+9zCxk8cg4\ni9fe3yY52Lp/eL94ZRU/1Y6MzvqNW9vkZw8O7088vckwGc3DAfDmvW3KiYv3v35x6yIu0VGD36vf\nYjP74HB7fuMS9VF0h83tzi2ulh/Yu1a+zE4vugeyNrjNs+uX73+9aW9xqxnd9TX8W3x068H1raYv\n8XYtOmfRS9zmU09fuv91kUu8fjea65vNfCaZHT751IP7MzPe4us3o7F3t94BLLZWCgAk4jHig3W+\nEZFz+sbNbZL9i1hWcLZ85MoqftKLzDmdVzwyzuL1O9vkpw9u1mQiRqy/Ftnen3d2t6nEH9h7bmsF\nP9Wh043GOd1s3GU9+8DeS09tMEpF5yx23F22ihfuf/2RrS268QgP09E9rlUfXN8zGxdojqNzvq6/\ny7MbG/e/vla5xE4vuusbJLb5+JUH13fBvsSddjSZzGzmM83s8bGr6/dfW01f4p29aJzTrb0mTHJU\niun7rxW5xJv3orm+V9+9S2rw4PAGSE82ePNuNPfLG9tHz5ZEwiI+WOPr7y1nLWTwyDiLt3e3KX/b\n4Q2QmWzyzdvRlKLudrZZST+wF4/FiA1Wef1WNDdQrbfNBftBJvPsxSp+qk3bi8Y5tcY1LlXW7n/9\n0csXGWWim8r1/D2urT043J67sE7Xiu5hHMT3eGbzwfU9tXqB9iyaeyU8vF+89OD6LhTWaQyi6SHc\n3mvDNE3Jzt5/bcPeYKcTjb1v3dolOVw/8tpK6iLv16O5X16/u32kJAvgsMG7u9E4i/f2tinGPnS2\nTDf41u3ogpnziEfGWdxq3mUte/QDLVibvLMTTfS9073Lpv2hG2i8wVvb0dxAB+NtLpe/LbqJx4kN\nVnnzTjQHgOvvcWXlwWH69IUyfqJLbxiNcwoO7wcHzvMX1xmlonMWk0yNFy8/sPfU+ho9az8SW9sH\nHZimqRYfHN6XKmu0p9Fc3xu3ayRHRw/vTWeNxjCae+XtezVys6P21uw1dtrR/D3fqR0tyQKUUxvc\nrEfz7N1ubrOaPuqcbNa5ub/MLGQQqbOwLCtuWdarlmV9YdF793s1Np2NI6+Vk5t8UI/GWTTHu1wq\nbR55Lc8m79WiuWFdv8bVlaPXlx5v8Pa9aOz1YzWe3njgLJKJGLHBCm/diWZqfJyu8fzWA3vPba3i\npxuMxuYZNXvNLlhTNsrO/deevbDGKBHN4fbWnT0Sw7Ujrz29vk6XaA7vd3drZKdHD++t8hrtSTT2\nPtir4cSO2ttwVjmIyDnda+9RzRz9e67nNtnuRPMs1Lq7bHzobCklNrjTXGYWMog6s/hJ4HVg4fRS\na3TARmH1yGur2U12vWichTutc6l61F4pvsHNg2jsDax9nlo/ai/vb0TmnCbpPV7YOnoApCarvL1t\n/gBoun2ID7m0Wrj/WiaVwBoVI5lFeOP2HvHh2n0yBMALl1aZZfYjmUV4d3ePzPTo4fbMhTWGiWgi\n0w/2a9jW0c/u6toanh/N4X2nWaOaPmrvYmmNzjga51vv1VnNHX0WLhY22O9F8yw0BnW2SkftVTPr\n7LhLZyGDyJyFZVlbwA8D/xtgLXg73rTOVnnlyGtr+VUag2gi4b5V59r6UXsrmU3uRXQDjVN1nts6\naq8Y3+BWw7xz8voj/KTL1Y3ykddz/hrv75o/AN66u098cPTwBkiO1nlr2/yB+u7OHpnJ0cPNzqaw\nxnYkwkQ39/ewraPO4sVLa0wze5FMcd9t7lFOHr2+ZzbXGCaicRY77i5r+aP2rqysRpY5tUZ11gtH\nn4Wt8gatSTTPXmda5/LKUXub9gb7/WUZSgZRZhY/D/x1QCjU61l1Lq9+6AMtrtAeR1Q2SdZ55sJR\ne+v2CvWueXuD0QQ/1eapzaOHdzm9ym7HvL037+wTG6yQiB/9eJ34GrcPzB8A79zbIzVZO/Z6brbO\n+7vm7b2/XyPP+rHXE6NV3r5r3hnebexRTBy9vkohC9M0d/Y7xu3tujVWckftPb+1xiQdzeF9MKxx\noXj073ltfY1BPBp7nek+F0tHn73LlbXInFOf44HhxeJ6pGy984hInIVlWX8G2PN9/1UEsgqAUeL4\nB7pVqeLNzJcxplOfWfqA57aqR17fLFVpDc3be3e7iTUskUrGj7y+mq9S75m3d9rhXU6tst02/0B+\nsLdH3j9uz4mtc+vAfPR256BGKXHcXna2xrsROKcdt0Y1e9xecrjOG7fNX99+v8am86EeQtkGaxr0\nawzDndbYKh2t6T9/cY1xRM6p59e5vHK0LHRlrcowFs36lFGizlMbH3JOK6t402i3Pp83JCL6ud8L\n/IhlWT8MZICCZVn/yPf9/+jb3/Tyyy/f///pTo1nLx79QK+srtC3zH+gt3Y7MM3e35sU4mKlijuN\nwlnUSY5Xjr2+XqjydtO8Hvf7tePsFoC13Bp7nvnI+/ZBjcIJh3cpGY1z2nH3qGaOX59trXGrbt7e\nfn+PZ8vPH3s9M13j3d0a8JxRe61JjUvlo9cXi1nEh2u8vb3PWjlv1F7X2udS5UOH93oJEj28/hA7\nmz7lO9UwTOxzbe3o83Bto8I4af7Z832YZY47iyurVfpW9GqHJnDjxg1u3Lhx1r9GNM7C9/2fAn4K\nwLKsPwn8Nx92FPDAWdTbPV75uZ9ltZg78u/X1quM4uY/0Le36yRHxw/vS9VobqAPanWy/nF7F8tV\n3PfM27t9sEchvnrs9TV7jZsd81rj9zp7VNLHnUUlV6Humdf+3u/tcfXbVn2EKCXXuNs07yxaoz0u\nlr7/2Os5VrjXNL/ptjurc6l6/H5JT9Z4Z6fG933sqlF7o1iDSyuVI6/F4xaxwSpv3d3nM89unfKd\napikjpeAn75QwU+3mExnx8qnOqi3e4BPtXD0bHnqQoVxIrotxSZx/fp1rl+/fv/rV1555Ux+j4c1\nZzG3C/jevQPiw5VjDdKnL1SZZerGm4g3a3XS0+MP45XVKuOE+cP7brNO3jrZOfV88/ZqnQOKqeP2\nNguruBPzmcV+t85K7ri91XyFZt/8A9kZNVh1KsdeL6VX2I+g5+RO62xVjl+fnSiz2zZ/fcNYk63q\n8evLUuXOvnl7k2STK+vlY68nx6u8Z5gQ4fVHkOhxea105PVMKgFjm9u1tlF7796rEx+uHJkWB3jm\nQoVZ+mApryqByJ2F7/u/6/v+j8x7z/u1OsnJ8YexWsjBLM5ey6yIznazTpbj9q5tVJmkoji86zjx\nE5zFSpVhBJlTo9+klDn+8G+WyvR885F+Z9RkJX/c3qpdoTOJIvJuslE4fphWsmVaA/PXN7CabFWP\nX18hWaHejeDwTjS5vHbcXj5WptYxe33T2Qw/1eLaxnF7Gcrc2Tdr753tOrFh9VhgCJAYVXlvx+zz\ncHPv5LOlkMuCH6fWjE6g67zhkZjgvrVfJ3dCmQYgPqry7rbZaHGnU8eJHbe3tVKERI/eYGzU3n6v\nTumESP/KWpVJ0vxh0xo2qWSPP/wXq2WGlvnD1Js0WXOO29soVejOzF9f32+yWT5ubyVfpjMyf33j\n+MnOopyt0OiZtzdLN7i2ftwZFhIV9g2X9e7ud2Bsk8scr0jnY2W2D8zae3/35P4dQHpa5daeWWdx\nu35yCRggPqrwvmHndJ7xSDiL7UYd+4TDGyA1XeG9XbPOYt+rU0getxeLWVjDMu8bVlxrDOpUs8ft\nPbVRwU83mU7NpsLu+ORI//JqmXHC/OHW85usl05wTuUKA6Ip01ysHLe3VijjTc1f3zR1cpmmmqvQ\nGhq+Vzp94JCa+yEU02UOembtfbDbID46fm0ATrzCnmfW3s29OpnZ8X4aQM6qcOfA7OF9r3VyCRgg\nOalyp/549C0eBTwSzmKnUz+xxg6Q9avca5q9ger9OpX0KTfQuMrNmll77VGd1fxxe9l0EiY57tbN\n1mm7sybrheMHwNWNMtOk+cN0aDW5eEKkv1WtMIpHUWNvcGn1pDJbhZ5v1t5oPIVkl4srhWP/tuaU\n6YxNH95NYsOTD+9ytkyzb/bzu73fIDU9nsUAFFJlDrpm7d1rNsjHTr4+O15lp2X277nTqVNMnPys\nZ2YVthtLZyGKR8JZ1LsHlDMn37B5q8pO2+zh3RodUM2dbC89q3K7btaeOz1g/YQaOwRlNtPO6bQy\nzVopD/ExbW9o1N4o3mRr5bi9K2sVJoad02zmM0s1uXpCpL9VLTOMmbV3q9bCGhVOZOhsFs2X2W7v\nN0hOTr5XVu2Kcee03WiS4RTnlDHfA9p3W9iJk+0Vk1X2XLPPwkG3QfGE/h0Emcy91rIMJYpHwlm0\nh23KmdKJ/5ZPlKh7Zlc4eOM2K/bJ9rJU2TacyfT9Nhulk+2lJhXjzml0SpkmFrOIDcvcrBlukqaa\nXDmhIXt1o3SfDmkKDbcPfuzEMs2llTLjuOnIu0l8fPJhc7FSYWC4B7TdaJKanWxvvVCma7jMttNq\nkI+d7JwquTJtwz2gRq+Fkzz5WSilqzQGhgPDQYtK7mR7dqJKrbPMLETxSDgLd9yikj/5Ay2mSjR7\nZss0vVmLtcIpN1C8zJ5hxsmQFpvl4on/lvGr3DVcp50kmyeWaQASk3IgdmMIQZnGY2v1+PUFdMg8\n23VzKzFu1prETqmxX90oM0sZbgDXm6SmJ9u7tGq+zLbTbJK1TrYXRZltr9PESZzsLFbsMt7E7N+z\nOWhROiUwrGartEZmn4XOsM3KaWdLskK9u8wsRPFIOIvetMWqc/JhWsqUaA3MZhYDv8XGKYe3kyxR\nd83aG8fbbK2cfMPmYiX2Ouac4bwyDUBqVuauQYbLvDINBIyTmzVzB9ztvSbJycnXtlo0X2a712yS\n8U9xTutlpobZbLVOEzt+8uF9sWK+zFbvNSilT76+jWKFrmHn1Bm2Tq0iVPNl3LHhKsK0xXrxdOdk\nOpM5z3gknEXfb58a6ZdzJdpDszfQKNbmQuWUVDhTotE3m8lMky0urZ2WCpfYd83Zm1emAchaZXYa\n5g6ceWUaCMpsJhkn243GqWWaKMpsu+0mudjJn92l1SJ+yjWq2VH3GhSSJ1/f5bWKcTZbs9+kkj3Z\nOW2Wygwwa8+btKieUgJedYr0Zmaf9f6szUbp5MCwki/RGZnfUnxe8Ug4i6HVOvUDXckX8SZmP9BJ\nosWl1dOcU5G2wUzmAZvGOfHfC6kija45e/PKNAB2rMKuwTrt3XqT5CllGoAMFe4emLO302qSO6VM\nA+bLbPtu89SGbCoZxxoVAhlUQzhtoBKCTGaWMhvpt0cNVu2T7W2tlBkZ7gH1Zi3WnJOfvfVikYFv\nNlALzpaT7a3YRbyxWXvnGY+EsxjHW2yeEumvFkr0fHOHqe/DLNnm8trJzqmaL+EavIHu7HWwxs6p\nZZpStkjTYCZzZ79J4pQyDYCTNEuHnFemAchZJWptc/b2Ok3yp5RpANKzitFBsoNek2Lq9OuLj8tG\ny2ytYZPqKZF+tZAFaxaITRlCb9Y8cXUKBIQB01Trvt9i7ZSy0Ga5xMgye3hPEi02yqc5pxL92dJZ\niOKRcBbTRJtLp9T0N8slBgadRdMdQGxKxTm5TLNiF+lOzdm7vd8iPj7ZMQFUciU6I3M37HajSfqU\nMg0E3PmmQTpkUKY53V4+UWTfM3d9+17z1DINBCsqTC73aw1Oj/QBktMSO01z19cZN1g5JdIPymwV\no2W2Pg3W7JOdxZX1En6qzXRmjs02irXYPOXwvlAtMo6brSJME+1TySUb5SIDls5CFGfuLGYzHz/d\nZmv1+NATwIVKiXHM3Ad6u9YmNiodWywWYq1QMpoK3ztok5yd/HBAmAqbe0B2WqezaSDgzreH5g6b\nutvCOaVMA+CkijS6Jss0jbmHdz5WZt81d33tUZPqCdPwIdJ+gVrL3PX1Zk3Wi/PLbLcNltmGscap\n9nKZYLnfTsMcm20cb3GxevLhfXGlyDRp7m/p+8HZcpq9C5Wi0bPlvOPMncVu04NJJphmPgFb1RKT\nhLnD9G69TXxyeqS/Xioyipmzt33QIu2f7izWCiV6BlPhvU7z1AlZgGqujDs2W6YpzCnTFDNF2gNz\nh037lL1XIexEiQODPaDT9l6FyFhFowSFPk0unDBQGSLlF6i1XGP2xrEOa8WTAzWA+KhktAc0S7VP\n7RdeXHEg2TVGGDhwg7PFzp1ytqwWmSSXDW5RnLmzuLvfJjanTHNprcgs1cI3tD5pu9EiNSfS3yyb\nzWT2Oi2y1unXt1ooMsDcDVtfUKap5Av0puYO7+agSXlOpF/OFHENltnccYvqKbx5CDKZlkHn1PPn\nR/q5eIEDz5y9kXU6ewcg7RfZN0i1niZc1ksnky8AEtMiOwdmrm8wmkCix0bFPtlWPAYjh3sHZpzh\nvUZQRTgNWytF/FTb2Nly3nHmzmL7oEVycvoHWrazEJsGvQYDqLXaZDjd3sVqkanBaKPRbZOPn25v\no2Q2FW72WxTSp9tbLRQZ+OYOt87o9AlZgFKuQHdiMPKedajkT4+EC+kCnaG56xtweo0dgp6MyaHR\nadxlbc7hnYkVqBt0Tn7SZaNyur2UX2Cnacbenb021qhIPHb6sROfFLmzb+b5u9eY3y8s5rJgTWm5\nZtffnFecubPYac6P9C3LIjYqcccQPbHuzY/0A+58x5goSrPfIp843d6FSolJwmDkPXIpZE5/+FcL\nBYYYjLwnHcq50w/var5otMw29D0q9unXV8oW8Mbmrm8S67BemuOcUkVaA7OR/kb59OvLxwscGCIM\nDMZDwKdaPF02NRMz15O5U28RH5/+rAMkp0XuHZixt9tqkZzOP1usUZFtQ/bOO87cWdTabTJzDm+A\n+LjE3QMz0caiSD+XScI0xX67a8Ree9g+dRcOBAwQk5lMd+xSnOMsNsoFRgYzmb7vUsmfbm/NcCYz\nslxWnNPtlQ1nMovKNCYzmdnMh5THZvV0e3aySKtvxt5u04WRQzJ5MtkDIBcrctA1Y2+nMf/wBkj7\n5thle502qTn9QoDEpMS9xtJZiODMncW+2zp1QjZEalbiXsPMgdoatLAT8+3FRiVja8M7oxbFOWWh\nzYoDyZ6xpl5/6lLKnn7YbFYKTBPmDu+h71Kdc3ivFYsMDXLnJzGX1eKczMkp0p+Zu75ZYn6Zppwz\nN9h10OnDNBns1DoFTqpA25BzqjVdYpPTrw2CTKZpylk055M94JAwYKgnU3dbCwPR5LTIztJZCOHM\nnUXdm1+mAUhj7gZqD1s4qfn2EtOiMefkTVqUMqfbC5p65uiJA9+bG+lvrRSZJQ1SIRdE+uulgtGe\nzCQ+P9Ke7GULAAAgAElEQVRfcQrGuPPjyRSS/WDn1Cmo5IrGMpla08Mazz+8S5kinaEhey2X+HSB\ns0gWjLHZap0WWWu+s8jFiuwb2s120G2TW2Av5ReN7mY7zzhzZ9Hstymk5n+gWatk7Abyxm1KcyJ9\ngOSsxK6hVLg3bc9tAAPEx0VjmcyiSH+lmIX4yJh07KJIPyizGVyUmHBZn1PTXy8VGBnqydSaXRjn\nT52+h8A59Q2V2WrNxYd3QBgwY2+/45KczbdXSJtzTnW3NbcEDEGZrW5oLqfZX1xFyFhF6gapz+cZ\nZ+4s2nNWFofIxYo0egYj/ez8zCLjF6m1zdjr+61TtTNCJKbmnNOimn4sZmGNCsbWhk8XHN5bK0V8\nQ5nMbOZDsst6+fRIf6NUZBI3V9OPLYj01wzuM9pruyQWOIuA+myo7OW6pPxFmUwB1xBhoNFrYSfn\nP3uFVIm2IcJAa7DYXi5WNEYYOO84c2fhjttUcvM/0HzSodU3w73uzdpzefoA2Zi5TbADq83KKevX\nQ6RmRXZbZpzTJOayUph/AMTGBe4Z4s7PFrB3ivkMAO2uPvV5v92FaXpuTX+zWjCWyewJlGnWi+ao\nz/XO4sPbZE/mwHNJWwucRdZcJtMadCikT2eWQcAuM7UJ1h215/YLIchkTG+ZPq84c2fRnbSo5Ocf\npk6qQGdgxln0/RbVBfbycXMrKsaxFmuFBZkM5jQtAp7+yUNPIZLTIrsGuPOj8RQSg7k1fcsCa1ww\nUmbbbbgLa/oXqwVj1Of99uIyzWalYIz6fOC5pKz5n51J6nOz65KNzb++ql00lsm4Q5dCer69ctYc\nYcBd0C+EwyHOpbMQwpk7i77fYcVe9IE6uCMzzmJkdVhdcHjbBjOZSbxzqtBSiGysQN1QKjxbMGQF\nkKJAra1vL6zpx+OnUy8h6MmY4M7XWi7xReydbEB9Pujob2atuy5J5h/eF6pFY+yyZtclsyDSXzc4\nxNnqu+QS8+2tOAVj1OfuZD6tG4LFmqacU2/qUsrOz2RKabOLPM8zztxZLGrIApSyDl1DddNxzGV1\nQZnGSTvGGCCLhqwAcgmHZlffOY0nM0j0WSudHukDpCmw39G/vt2mS2wy/zAFSM6K1Fr69kQifQjK\nbCYGrRqeS4b59i5Uze0zavZccnGRTMbMvdkZuNjJRWWvAiPLEFNv6lHKLXAWtmOszDaYzZ8BArPU\n5/OOM3cWAfVy/oFTzhfoTU1F+i6rxfn2ihkzmcxkNoHYiNXSyevQQ9hJh7aBTGa34cEkN5e9A5Az\ntDJCpKYPwT4jE5nMfscluaCmDxCfFIxQkRs9l8yCMk0yEYNxnp2G/ufXHrjkFhzeF6oFZoZ6Mu7I\nw04tymQKjGNmDu9FA5wQZDJDzDzrIoFoJWeuzHbecebOIji8FzFAzEUbs4TL+oIyTSlXoDs28PD3\nXRjZOM78Mo2TdugMDTgLAfYOHO4zMjBoVRegXgKkrYIRemLD8xY2ZAGSswI1Az2ZlkCkD2GZTd9e\nZ7g40l8v2ZDoG8lkuqPFPYTNSpGpoZ7M0Hep2PMDtZWCw8gyVXJ2WSnMt1e1Hfozc1t8zzPO3Fks\nWqcA5qKNINIfs7Yg0q/kHXoTfXt77aAhO2dvGmAuk9lvi0X6TqpgpKlXF6BeQkBPNNGTEanpQ5DJ\nmCAMdAYu+QWHNxxuZjVAffZGLs6CSD8et2DkGMlkupP50/4QZjKGSsDW4kh/reQYoz4vmgGCwDkN\n/aWzEMGZOws/4S1syK4WHUYGnEW778HIxrbnR/pVxzHS1NtreQvXKcAhPdFAJrPf9oQi/ULazMoI\nEeolgG1oM2uz75IViPQzMTM9GXe4uEwDQSaza2DZ3qK9XiHiEzPL7/pTd2EPIRjiHNMbjrTtTWLe\nwsN7veQwjZs5vKdxb2EgulZ0GBvKZM47InMWlmVdsizr31mW9S3Lsr5pWdZf+/B7BqMJxMenSpyG\nMBVt7DaDSP8Ukbz7MBVt7LUWD1nBYSZjQGNCNNIvZsxsZm12F9f0AeyUQ3vgadtr913yC9g7EPRk\nGgbKbJ5AmQYgQ5H9tr69nsDhDZCYmunJDHx37gZfOBziHJops00XrGqBgDDgJ81UEfzYcGG/cL3s\nMDbknM47oswsxsB/7fv+R4HvBv4Ly7Je/PY37BwENf1YbP7pvVEqME3of6C1pktcgL2zWnAYxwxE\n+h2XpL/YXpDJ6NtrCEb65ayZfUaiNf1C2sE1MCfjDl2hSD+fKNIwkMmIUD0BMjGHhqd/fYOZS3lB\nAxggNTNDfR7iUl3gLABik4KRuZxZ0mW9PP95KNlpsGZ0B3qZjDfyYGxTKCw4W8oOs6WzEEJkzsL3\n/V3f9187/H8PeAO48O3v2W26xKeLD9Mg2jBQFmq7xAXKNOulAhMDN9CB55JeQL2Ew0zGQJmtIRjp\nl3MFI4SBtmBNv5C18QyU2UQjfSdtZvmdSJkGIBu3jczlDHyxwztNkbqBMtvYWlzThyCT0XUW48kM\nkr2FA6OxmIU1doJAUgPtvgtDhwX9dC6sOPgpd6mWJ4CH0rOwLOsq8CngD7/99T2BISs4XBlhzXB7\netHGQccjJeQsHGYGMpmG55IWOLxN1U1bfZdsXCCTsQsMDEwBi1AvAcpZh64BwoBopF9IFXBHhso0\nApF+PuEE/TBNLNrgGyIbK1B39a9vEndZE3AWqVlRWwBpt+HBeDGtGyA2cQKtDQ2IkkvymRTAUi1P\nAKcv2TEEy7Js4PPATx5mGPfxD3/x5/Hf7fDyyy9z/fp1rl+/fuLPuB9tNFycXFX5dwnWKSx+ODYr\nDn6qg++zsL8xD62eR07EWRhq6nX6YpH+SqHAyIDGhDdyuVy8tPB9Fduhb2BOpjd1KeYWO8NitsAH\nnXe17Q1ZTPWEYHeZCerzWIC9A5BN2DR7+vYWbfANkbYc9jWd017LE6J1AySmDrWWprNoeUL9QgBr\nHOh+lwunKwaeJW7cuMGNGzfO+teI1llYlpUE/h/gl33f/40P//v3/bm/xJe+9AEvv/zywp8Vmzjs\nNjo8t6XuLBqeGPUynwlumrY3pOSo30BtgXUKENRNTfRkOgLUSwh6MhNLPxLuCUb6FdthiL490Zp+\nKevQn+jbC3j6ImUvh4Pegba9RfrbIfIJh44mYcD3ffykF4hvLUDGcmh6evZEBzgBkr5DvaP3POy3\nXRICVQSA+GEm89FrK1o2o8KHA+lXXnnlTH6PKNlQFvAPgdd93/97J71HlHoJkJwWtKONZk+sTANg\njQraXHaRISsInAUGlt+J1vTXSjaTuP5hKlrTr9oOYwM9GdGGbClnM5jpX98kJlamKWRsI2W2WXLx\nahgI2GW6mYw7CFT5CvbieDGXcLQzmb22S2Im9uylDDiLgBkoZi8x07f3JCDKnsWfAH4c+AHLsl49\n/O9z3/4Gka2XIZK+w54mPbE9EKNeQhBt6Dond+RipxffsLl0GvwYTc26qWikb6onI1rTXynYRthl\nI8Rq+lXbDGFgGveEyjTFjENPM5OZTGeQ6LFRWXy/OCkbT3OIc6chRiMHyCVs7Z5M3XVJCZA9IGCX\n6ep+H7hi/UIInNP+0lksRGRlKN/3v8QCZ9Tqu2QTYt4/TYEDV/PwHnpCkT4E0Ybu8rvu2OWKwOEG\n3O/JVIsZZXuikf56JQ/JLrOZv5C2PA8iu3cgaOCbYJeJTOQCVB3byBCnn1w8FwBQzjsMNFdG7LW6\nMMkKNYALGYebrQ+07Inob4dwDGQyza4nxAyEYE5Gd7FmoyceiKYwQ30+7zjTCe6ORKSfiTkcaH6g\n3sjFESjTgJm6aW/iUlywTiGEiUxmMJuvvx0im07ANE3D7WnZE63pb1TMZDIiE7kQUJEnMb1IOBgY\nHbJSzC18b8W2tYc4dxtie70Aillbe7Gm6MAogJO2cYd6f0/RAU4wI3bW7nlCM0AQ9GQaBrY+n3ec\nqbPwRp7w4Z2LOzR7epG+N3YpCJRp4DCT0XROg5lHWSDSh8O6aVvPXsDeEc1kbG16oij1cr3kQMrT\n7smIaHWE9nQzmVozGOoSybxWHP3ld7WW2MwRQDnnaPdkRDf4QigRoHd9wQCn2PWZYJcFe73E7GXj\nDi0D7LLzjrN1FmOxhixAPlHQXuMdRPpiN1DG0teYEB2ygiCT0a2bikb6APGpQ62ld+CIrG8AyKaT\nMEvQdNWlVYOa/mKtDgga+LOE3rXJRPorjsNYk11W73jC7J2KrT/xf+CKDYwCFHM2/ane9bUHYsxA\nOJz413UWgtP+EM7JLJ3FIpyps5Ap05gQJBrMXOFIPxvXZ4CMcKku0OoIkcbhQFNjQpS9AwGXfV+z\n7OUnxXj6EGQyOxqZjKhWBxA0iZN6mUytLV6mMTEnI6K/HaLq2NpDnA3BDb4QZjIP7/AuZvQzGalA\nNOnQMaTEeZ5xts5iKn54F9IF7ZURQ98TLtPkEw6tvuaKA8GGLJipm04Ti/W3QyR9m7qrHi0GNf0R\n1cL8RW0h4hOHPQ3nJKrVAZBJJWCaouGqS6vWOy4JQerletnWnpORoZGvFPQzmYBGLp7JDH09e92x\nJ3x4l3L6EgFSgWjKwTMwVHnecabOYuC7lPNiD2Qx42g7i5HgOgUIog3dVHgSXyzhGiIbL2g39XxB\nnj4cZjIa7LJgfYNYTR+Cnsy+Rk9GZqgLAnbZrsacjAzVM5iT0ctkZBrA60X9TEZ0YBSgatvaPRnR\n9esQsMt0BYn6U4+SwLQ/BOwyE3My5x1n6ixGAiuSQ5Ry+mu8RRenQbhfSO8GmiW8hap8IXTrpjLs\nHYB0TG9lRE0i0oeQXaYenYrqb4eIT2322ur2RPS3Q+QyQU+m5an3ZEQ3+ELQk5lq9mREB0YhZJfp\nHt7i/cKqrS8RIKK/HaKUNSN2dt5xps5ibHnCkXclX9CumwbsHYloQyOTEVXlC6HLAJFh70DALmv1\n1A8cmZo+hD0ZvUhflL0Dhz0ZjUym1fXICB7e8GBORhUyA6PrZRtSrlYm4w7FaeQmJv77vtiqFgjW\n0YwsvcBQdNofgkB0Ka26GGfqLET0t0MEU7l6N9BMoiGrm8m0eq6QKl8IJ+VoTeXKsHcgWKutk8ns\nt+UO77Qmu+xAcK9XiKTvcKDRk2n1xSN9gPjE1pqTkWkAm8hkZGjkJib+R74nfHivFvX1ZGSYgaYm\n/s87ztRZiOhvh9AVcr8f6ZfFIv2yZrQRqPLZwltrAwaIunOSjfSD/ULqh+mBoCpfiGzc1uKyN3su\nmZhYVgiQsmzqGj0ZmUgfwv1C6n9PmYFROOzJNNXt9aaL9bdDBJmMXk9G5vA2wS6bxMTPlqptRrb5\nvONMnYWfdIW2XsLhoJVGtNHqBfrbuZzY6V21C1p105qgVkeIYtbRmsqttz0SEoe3LgNEhr0Dh2Uv\njUymJcHeAf05GZkyDUBScxldV4K9A/rssoEEEzGTCjKZRkdjTibmslIQc/YbZYeZprTqLOEJMwNN\nZDJPAs7UWRCbULLFdiGtl/TWeMssTgP9qVyZFcmgzwAR1d8OEWxKVY9MZZZAQpDJuCN1ex1BVb4Q\n2Zit1ZOR4emDfibTn7qUJZxFYmazr9HAH0jU9EE/kxFdvw6HmUyyy2SilsmE+ttrZTGyx2rRYRrX\nF5M67zhbZyGgvx1io6In5B5E+uJljNWiXiZT73hSNf1KvqBVNxXV6gihywBpSVAvQT+TcYeekFZH\niFzCoa2h+y3D0wf9TEZUqyNEcuZoOSeZgVEIejI6mYzo+nWAZDwB0wy1ZlfJljcKqgiL9LdDmNKT\nOe84U2chU6apOnmID+kPJ0q29ltykf6a5n6hhqD+dogVR69uKsPTh8MVDhqZTGfgkhPcGAwBu0xn\nTkZUqyNEPqm3/K4vUaYByMb02GVDxDb4hkhha83JiG7wDaHTk5lMA/3t9bL4/RIbq0urHngujBzS\ngrplYSaju7vsvONsnYVEQzYWs2CcD4bBFFB35SJ9XQaIzDoFCDQftHoyfbnDW3cqV2YJJARlLx1p\nVVH97RCO5n4hmYFR0M9kRPW3Q2RiDs2u+ucnugQyRNK3lXsytWZXWH87RHyq3pPZa3rEJuLkklwm\n2MLc7KplMk8KztRZyBze8ED+UAUyi9MA1is5SPSZzmZK9mRU+UBf80Fm3TscTuVqZDKyNX1ddpmo\nVkeIQtqhO9aI9CWonhBkMh0dZyGx1wsOnYUGu2yWENvgGyLQfFD7e+42XGISJWAIymx7inMyNYn1\n6yF052SeBJyps5BpyIJetNGQLNOkUzEY55RTb5l1CqCv+eBKUi+DnowG9VKypl/J62Uyoqp8IXQ1\nH8YSVE8IMhlPo4E/TYjPAEGgXqeqwz2dBWUhEVW+EBnLVt5dtteWYwaCnp6MzPr1ELrssicBZ+os\nZMo0EEQbqlO57b5cpA8Qm6jvF5JR5YOwbtpj5qtlMrKR/ppmA19mCSRA1dEbfBJV5QtRzjkMNNZq\nywyMQpjJqF+fzAZfADvpKGcy9Y4H4xyZtPjjn9HQfJBZvx4ibdnK6nWyM0CgP/H/JOBsnYXEkBUE\nqbAqA6TdlyvTQOAsVKdyZSP9bCbIZBquWt1UNtJfL+tpPgwk1jeAPrtMZqgLoJS3tTQfZAZGQW9O\nJtzrtSq41wvATtnKu8t2D2nkMsjHHeVMpt4RX8oYIhtzaCg6pwNXbgYI9Cf+nwScqbOQKdNASE9U\n+0A7kpE+HGYyiqlwV2KdQghLgwEiM5ELAV3QT7n4igSQkS/XkF0t2Fo9GdmGbMAuU3/4ZQZGAco5\n9QZ+sNcrL6WHrtOT2Wt50mWhfNKmrUgYOPBcUsgFhtm4Q1uRXdbqedJVhDQO9aUO91ycqbOQGbKC\nQ0Eixbqp7DoF0KubBoe33A2ro8M9FNTfDpHPpACftjdSsjeOyUX662W9noyoKl+IqmMrT+X2BmOp\ngVHQY5cFDWC5e7OoMSezJ7nBF4Khyq5iT0Z2gBOCs0GVXSY7AwSHW5iXOtxzcabOQmbICkId7ocX\n6etsSu1P5co0oFc3leXpW5allcnIqPLBYSaT9PAVUxnZmr5OT2an4WKNHKlIv+o4DBUn/mX3egGU\ncrayDrfsBl8IqM+qczKtvtyqFgBHQ0+mMxBfvx4imJNZOot5OFNnIXt466zx7kns0w+hw2UfIq7V\nESKJ+soIWfYOQGxiU1N0FjKqfAB2Ltgv1OnJ7xeSVeWDMJNR++xqLflIv+rYTBSdhexqGAjYZaor\n+1Vq+oWMQ09xPUxn4En3CwsZ9S3M7siVD0QTDu3+smcxD2fqLGSGrEBv0KovIYYSIqvBABnhSdX0\nQU+9TramD6F6ndoDIqPKB2BZYI3UMhlZrQ4IGviq62H2W550pL+uMfGvwt4pa6jXyQ6MQtCTUXVO\nMuvXQwS7yxRLzmMXJyMXGNqaejJPAs7UWYjKHoYopNVXRsjob4fIJ9TpibLrFEBP80E20gf1Bn53\nMILYVKqmDxCbqmUyslodAIVcGqwZXl++J7Mvob8dIpj4V3O8Dc+TPrxXNXS4ZVT5QpTzjnLZyxu5\nOCm5v2cp5ygTBnoTT4oZCPp6Mk8CztRZyNb0dZp6I+TYOwB20laONiaSDVk4bOArZjKy7B0I9gup\nTOXuNjzpmj4c9mQUnJNKTT8Ws7DGttJ6GNkNvvBg0aXKfiHZvV6gR0VuS27whaAnozrxL0vrhrDM\npuacZPd6ATgZW2tO5knAmToL2cNbZ2WESqTvpB1lLvtMIdJX5bKrsHcg6MmoDD7tNuVr+hBqPshf\nn6wqX4jYxGGvJW9PdoMvgJNLARad3lDangp7J9Dh1igLSTsLWzmT6U08KVo36KnXyU77w+EWZo2J\n/ycBZ+ssJBuyFVu9qSejvx0iaOrJ2xtPg8NbVJUvRD6lVjdVYe+AuubDnsLuHQiGKlXYZSo1fYD4\n1FbaL9SUFFoKodqT6QzlB0Y3KwG7TAUqNPLVgnpPRnaAE/SUMVWqCEHZa9ngnofInIVlWZ+zLOtN\ny7LesSzrvz3pPbIN2ZWCeioso78dQjXaCFX5slm5w1u1bqoa6QcMEHl7QU1f3p7qCgdZVb4QSd+h\nrtDAV6npw2EmoyAQJKvKB1DMpyE2CbJKSXTH8mUhnUxm6MszA3Woz2PLpSqoyhcimJNZZhbzEImz\nsCwrDvxPwOeAjwA/ZlnWix9+n+zhvVpQkz+U1d8OoUpPDPS3xVX5QhQyavuF9lqu1Lr3EHnFnsyB\n65JWcBaqPRmVoS4IFlWqUJE7Q/maPgTsMpVMRnavFxyu7B85AVNMEn3JaX84nJNRbOCPLDmhJdBz\nTpO4J6zKF6Jq6yljPgmIKrP4LPCu7/s3fd8fA/8X8Gc//KZ1yZq+6hrvhhdE+pmM3OldsR2l/UK1\npvyWTYCS4qbUeseTnsiFYCpXZVNqQzHSz8UdOgpcdpWaPoRTuQ8n0odA80Glga/SAAb1if++Qk2/\n7GQgPg5mXiQhO8AJesqYs4TLhqSzWCnYjJfOYi6ichYXgTvf9vXdw9eOQLYhu6G4MmJXUn87RNVx\nlG6gvbZLYibnCCEoe6k08OuSWh0hiorqdbJaHSHspKO0X0hWfztE1lLLZLoKkT4EmYxKmU12g2+I\n+FRN6lR2gy+EmYytlMlMVSJ9J9jCPBjKbWGW1d8OsaqpJ/MkIBHRzxXiD/7tv/3K/f+/fv06169f\nn/v+MNrwfaQO/lrLJT6VP9zWimplr3pHrSGrWjdVjfQLGVtpKldWqyNEkMk0pb9PZSIXIJuwlaZy\nVXj6ELLL5O0NfDVnoTonI6vKFyI+cag1Xa6sl6S+b5aUp5HHY8EW5t1Gl6ub4t/rDoOljKL62yE2\nNJUxo8SNGze4cePGWf8akTmLbeDSt319iSC7OIKXX35Z6ofm08GgVcsdUS6khL9vT4GnD4EO9zQu\n//AfeJ70SmZQ57KrsnfKioNPKuwdCLjstdZt6e/zRi7r9pr096lqPvRnLmXJgVEI9gupZDIjhUgf\nQnaZ/P2pQiOHIJORLXsF+ttd1st5aXuxQ+ck4yxqzYAZGI/L2Vov2/iJQE8mZp0pSfQYPhxIv/LK\nK6e/OUJE9Vf5Y+BZy7KuWpaVAv4i8M90f2i4/E5W/vCg4ylF+mslWynaUBmygrCBrxjpK5SFynlb\nKZNxh3L62yGKGYeeAj2xN5XT3w5hp9TmZAYzeaonqLPLZDf4hkijVvaaJuRWtYRIzOQ1H/ZaXRhn\nSSYkT2+CIU5Z57TX8pSYgY4dh0kGd9CT/t4nBZE4C9/3J8B/CfwW8Drwj33ff8PEzw5TYRk0umpl\nmvVyDuJDxlO5pl6r5ykd3kEm8/DYO6rqdSrsHVDPZPpTT+nwLqRtJSqy7AbfEKrsMpUGMATqdSrr\nYWT1t0OkFHaX1ZpqhzdAYmZLb2HeU1jKCIe7yzS2MD8JiCzf8n3/X/q+/7zv+8/4vv/fm/q5KvTE\nhudKq/IBpNNBU082mmoP1Gr6G2VHiS6oyt5ZUaQiq7J3Ah1uhUhfsaavuh5mHHNZVYj0VedkVPZ6\nAeTiNi3JTGY0mUB8zEpRjkYOaiv7VZmBEDonuWdvv+OSVCCXgN4W5icBj1ZxTgBJX17zodX3lGr6\noKbD7Q5dbMnFaUDA4Ej0mc6mUt+nGukHg08KPP2Z2uFdcWylTCao6Suwy3JqQ5WqkX4wJyP/91TZ\n6wWHK/slezI7DRdGNvG4JDUQyMRt6Z7Mfkct0odAGbMhmTkduGr9QjgMRBXFx54EPHbOQiW6UdHf\nDqGiw+2N1Gr6+VwMRnnpFRy9iVpNX1W9bjCTn8iFcFPqw6vpq2Yys4T8tD8EPRnZtdpeX22DLwSZ\njCs5J1NTnPYHyMUc6UzmwHVJSkqqhgj0ZFSqCGrXp7q77EnBY+csVJbfBTV9tRs2OZPPZLpjtcM7\nrJvKNvD7ig3ZQL1OXod7pEi9DIYq1dg7KpF+sIxOIdJPqUX6JQUd7p0Dtb1ecKjDLVn2qinu9QK1\nTOZAcQYIAsKArHNq9lyyCiVnUN9d9qTgsXMWKvKHqpE+HO4Xkmzq9aaetCpfCJWpXFn97RDFXBbi\nI7yeXAM/UOVTa+DPFBr4U8VIf0Uhkwn1Lwr5tLS9isKmVBVVvhDFrHwms9dSLwvZSfmJ/2ZPXqvj\nvj0FPZl231PqF4Ja2etJwmPnLFSE3AP9bbXDW6XspRrpw6F6neSglSp7J6Ai29KZjIoqH8BGxcZP\nudI63Ko1/dWiLb2MLoz0VbDiOIwkMxmdSL+Ylc9kDlxXuabvpOUb+KpLGQHslDy7TEV/O0RWoSfz\nJOGxcxYqg1b9qfw+/RAqddOh71LOK5a9fIcDybqpKnsHQs0HuetTremXnBTM4ngDcR1urz8CaxYo\n30lirSTPLtOJ9FcU1sMcuGobfEFtd9mB55JRdBYqUqftgUtOsQSsQhjwRi52Ws1eLq42J/Ok4LFz\nFiqCRDqRfk5Bh1tFfzuEig73JOZK794JkZjKCQT5vi+tvx0i7MnI7BfaOQj2eqnU9DcVltHVWmp7\nveAwk5Ess9UVN/jC4aZUyUym2XPJKEb6KlKnroL+dgiVMlt37Cn1C+FQT0ZBfOxJwWPnLFQYJ8E+\nfcVoI2HTGcrdQKrrFCDQ4Zatm6pG+gAJX67s5fZGgKVU04dDKrIEl10n0i/bWYiP6Q/FezL1jpoq\nH6ixy1T3ekHALpNdq93qq5eFygor+72RS0HRWagMcXYV1q+HWOpwz8dj5yxKCoNWOpG+nZIve01i\nnrQqXwiVuqmfctkoK/ZkfLlMJlTlU0VCUr1Op6Yfi1kgqcOtor8dQkW9rtVTnwFS0eHWqelX8/IN\n/O7Eo6AY6as4p8HUpaRYAlbdwvyk4LFzFirSqpOYGnsHDqMNyRtIRZUvRF5yv1CnG2g+q0b6geaD\nxOGtwdMH+Qb+QcdTjvQB4pIrHHR4+uFnEH4mImhpzACtleUFgnTKQipSp/2pS0lhgBNCwoCks/DV\nmOmVDJYAACAASURBVIEQEAZUJv6fFDyezkKyqTeJy69IDlHMyqnXTWdTSPRZr8jt0w+RTzq4EgwQ\n3UhfdlOq6u6dECkcGhIrHA489UgfID6T03wIePrq9qyR3JxMZ6DuLDbL8j0Z1Wl/CJyF7MS/Tr+w\n6tiMLTl7I9QGOOGw7KWgJ/Ok4LFzFiqMk1lCXnwlRCkrVzdteF0Y5clm1P60sg183Ug/YICIP5DB\n7h11exlLbqhSJ9KHcCpX7vBW5ekDxCQFidyRi614eJedDMQm9IfiOty9iXpNf61kSy+6HPouFcWy\n0FpJXpBoHPOUqwhlxYn/JwWPnbNYLdpSy++ms5lWpF/Oy0UbtaaLNVG7WeGwgS9BF9SN9GXnVoKG\nrPr1yepw60b6KeSGKjsDT7mmD4fsMomeTHfsUlR0Fio63DplofWSw0wykxlZ6pH+ugL1eao47Q+h\nnsySDXUaHjtnIbvGu+F2YZxTjvQrkk29WstT3rIJoeaDuD1V/e0QsoNWOuwdCDalyvRkAv1tdeeU\nlsxk3JHaBt8QSV9uv1BXcYNviLgku2zge1QV9noBrJayEBtLreyfxDzlGaDNivw6GlWtDlDfXfak\n4LFzFhtluehmVzPSXynINdn2NSN9WQaIzu4dACdtSzXwVVX5QuSTjhQVWYe9A/I9GVX97RApyfUw\n/al6TR/kdbhHitP+AKlUsLJf5vqmcfUZoGIuD4k+3Z6YDvfMn0Gyx3pFXpUPYLUkPyfzJOGxcxZr\nRQdSLqORWLihs08/tCdzA9U7Hklf3TlVJRv4qsJOIYoZOSqyDk8fDtXrJBr47kidvQMBFVmmJ6Mb\n6actWyqTCWr6ej0ZGXbZyFKf9geIjeXEx2aKA5wAMSuGNckJU59b3aCKkM+pHWuPsg73o4DHzlmk\nEkmYJdlr9oXev9/xtCJ92bLXgacX6VccWyoV1o30ZXsynaFepF+Q5LIHe700nJPkepj+VE1/O0Q2\nLrcpdai4wTeErEDQRCPSB4hLTPyH+ttrJbVIH0JBIjF7O41g2t+SH/YHDneXJbtBhrLEMTx2zgKC\n6Ea0ThtM5Ko//OtlOR1uVVW+ELJ1Ux2ePhz2ZGbih4031CvTFDNyXPauBnsHDocqJXoyQ9+l6hSU\n7QVzMhKHt+VSVWTvQNCTkVl0OYurzwDBodSpYCaz3w70t1NJef3tEHEJZcy9lkdco+Rs5+MwztLp\nL3W4T8Jj6SyC6EbsBjpwXeUVyQArxQzEpgwnI6H366jyQcAAkSl7dQYdrTKNrA63N3YpZNQPU1ku\ne3/qUsmr25PVfNCN9POSa7WDGSD168tKLLr0fV9ZqyNE0MAXs7fb0KN1g5yejC4zcKnDPR+PpbNI\nSkQbza6nRfXMZCwYOjQ8sWix3dOL9Nck6YnuyKWQUj9sVhw5HW4dnj4caj5IZDLBXi91eyXJZXRj\nq6NV07fTttSczDTZ0Tq8swlbmPrc6Q3Bt3ByKWV7acTZZXuazMDQnug6Gp29XiFiE7mezJOEx9NZ\nSOwzamrs0w8ho8PdGXpK+tsh1ss5iI2YzMToibo1fVnNh/6so3V4y2YyQ6ujFemXJDOZadxlo6yZ\nyQj2ZGYzH5KelrOQEQgKpv0LyjV9kBMI0mUGwmGZTdDegeeRUpRwDZGYye0ue5LwWDqLjESdtt13\nyWs0ZEFOvc7T5OnncgE9sdkVi751a/qym1IH6EX6VUduqHJsuaxplGnKeVuKXTZNuGxoHN7FrENv\nIvbZNdw+TJPkMklle3baERYIqjVdYopLGUPISJ3WXf1IPxtzaAk+C01Pr+QMajLKTwoeT2chUaft\nDD3yGpE+HC6/E7yBPA1VPpCvmw40I/3NinOoXif2/hEuqxoN4DXJTanTuMu64kQuHGo++GKHzWzm\nQ8plXXGDL8j1ZHYP2Ts6CDIZseurtTrKG3xD5JO2sLPQnfaHQCJAdIizqbHBN0RQZltOcZ+Ex9JZ\nyAgS6ezTD5HyHeqCmUx/4lFWXKcQIj61hTOZge9StdUP72IuC7EJXk9sv5COVgeEmYz4wzhNdrQi\nfRkd7nq7B9M0mVRC2V5FYqhyt9khPlH/7EBuZX+945Ly9ezZSQdXUCCo0dWP9O2kQ1swc2r3XfIa\n0/4QBKIyczJPEh5LZyGzz6g79rQifThcGSHYI9HZshkiMXXYF3QWI032TqjDLTr4NIl3tKiXD1Y4\nLE5lgpp+V6umL7OMbrfpEhvrHaYVxxZe4x2o8undK6WcLbwepu52tGaA4HBORrCBr6O/HUJmiLOj\nsX49REZy4v9JwmPpLOyUeJ22N3UpaR7eMk29oe8pb9kMkfQd9gWd0zjW0Yr04ZABIuicdCZyAYp2\nCrDwBos1H+rtHkzSWjx9mWV0utP+IDcnoyO0FKKcE9+UeuC5ZCw9Z1jMig9Vdoaedr8wWKwpZs8d\nelozQHBY9pIUO3tS8Fg6i0Ja/IYdzFwqmmUhmanckaVXpgE5uuBUc8gKxNXrJtMZJHrKqnxw2JMZ\nibHL7jU62jz9DQnNh1qrox3pBz0ZsSztwO1orWoBOXZZs6e3fh3kpE47Q31ySVGizNbV7BeC/MT/\nk4TH0lmUsg49QWcx9D1l/e0Q+aR4k21sqavyhZBJhWdJl82KXrQoqsO91+rCJEsyoR7pQ5jJLD5Q\ng0hfs6ZvB5oPvcHinkzddUlplmnWJPYLHXguWc1If6XgMBZcq93ud7RmgCDoyYg28L2R3rQ/hM5J\n7Pp6E/X16yGCLczLBvdJiMRZWJb1c5ZlvWFZ1tcsy/p1y7KKJn++DHdeRzkrhJ0UFySaxtWFlkJk\n4w7t3uIbNoz018vqu3dAXId7p9EhpsnegUPNB4GylwmevowOd7DBV+/w3qjY+CnBmn5fb68XyC26\nbA9dHI0BTgioz6KZjIlIv2yLU5/7M0+LGQhyczJPGqLKLP418FHf9z8BvA38TZM/XEaHexxzWS3q\n3bAymg9BpK9JTxTU4d5teDDJkYjrfYxpyxGa66g1XeKa1EuAhC+2X2iv09Gu6YO4Dnejqye0BFDI\npQEft7e4J9PWXL8Oh+p1gpmMN9Ib4AQ5HW4doaX79iR0uHVU+UKUcnJ6Mk8SInEWvu//tu/fX934\nh8CWyZ9ftcXrtCYi/UJGbNBqMgkmcnV4+hBqPgg4i6ZrJNLPxsXmVvbbLsmZXmQK4abUxfZMNGQh\n2CUmMifT7HW0I/1YzArmZAQymWCvl971hewyEXjjDkWNvV5wuLtMUBe7P3O1aeQrBfF1NCP0+4Wl\nnC3ckyn/Vz/Aa+/taNl7nPAwehY/AfymyR8ow53XZe9AUDcViTb2W32YpkjG1Xn6EAgSiZS9TNT0\nQVy9zgR7BwJ2mUgm0+zq6W+HSAr2ZExE+iDOLjNR05fpyfQm+of3ekWcXTZEvywkU2YbxzxlSdUQ\nlby4nkw7+xqFbEbL3uME5VPNsqzfBjZO+Kef8n3/C4fv+Wlg5Pv+r570M15++eX7/3/9+nWuX78u\nZHut6DAWuIHGYzORvqh63W7TJaaxIjlEMeNwS6Bsst/WX6cAh3MrAplM3QB7B8TV65rdjjZPH8TV\n61wDNX0I2GUizkJ3VQsc1eG+tlme+96B31GWVA2xKcEuG6M3AwSHyphxsUxmGndZLek9f1VHrMw2\nmx1u8K3q35+LcOPGDW7cuBG5nUVQdha+7//QvH+3LOs/Bn4Y+FOnvefbnYUM1ssOMwFnsd/qw0w/\n0q/aYlz2oKavf9gUsw79usD1uWZq+nbKwRNwFk0DNX0I9guJZDKmIv10TEy9zh25XCle1raXECQM\n9KYu5ZyBMtuhDvciZ6Gr1QFQdh5kMot2Wo1jrnakv14O1tHMZhBbUAeZJfSWMkI4JyNAhugEVYRs\nWu9sEcGHA+lXXnklcpsnISo21OeAvw78Wd/3B6Z/vug+o51Gh9hI/2EUjTZqrQ5JA85ClJ7YMLA4\nDYKejMh+oVZfn6cPkE/YQlx2d+TipPX/nqKZTHfS0ZJUDRGwyxb/PQd+R3vaH8R1uEeW/uEd9GRs\ndgTmZKYJ/UjfTucg0afjzlevG499SHmsadoTVcbcOTDTL3ycEFXP4hcAG/hty7JetSzrF03+8FLO\nhmQPrzv/BtpttUlM9Vm7qwWx5Xf7nTYpX99exXYYCiy/C1T59A/TUsahK9CTaQ9c8gl9e3ZajIrs\njjraNX0IdomJqNf1p/qrWuBwPYxAJjP09WndIN6TmcQ72mQPCJQqRZzFLOGyrmlPVId7t9GFSUa7\nirBeElPG3G11jDADHydEkkP5vv9sFD83xLffQM6cJXrBRK7+4SZKTzzwOkbYO1VbLJPRlVQNURKc\nyu0M9davhyikHfb69xa+T1doKURecCp34Oup8oXIxsUymbHlagkthRDtyUwTelodIYItzPMP70B/\nu8dGRb+HFwoSvXDt9N/93kFbe68XHNXhjlmnx9LBBl+j42OPPB7LCW4Q0+He77RJa27ZBNioBFO5\ni5bfHXQ7ZGMGMpmiGF2wPehor1OAQPNBpCfjGYr0RTel9qYuZQOHt5MSy2SGvn5DFgJ2mch6mEnc\nZa2of32i7DI/2eGCgYZs0l+sVFlrdmGcI5nQP2ISAsqYu00zh3c+F+hwt7rdue/ba7dJG6giPE54\nbJ1FYrZ4CvjA65BB/wOtFtOAxXA6f9Cq2WsbKdOsC9ZNA/aO/sMvul/IVKQvymUfzDpUDJSFnLQt\nNFRpYikjBJmMyKbUaUJv/XoIkTkZtzcCyz8cGtRDSkCHe7tuJtKHUPd7vjOstdtGZoAe6MnMt7fv\nto1UER4nPLbOIjlbfMM2uh2ycf0PNJ0GhotT/c6wg500U/YSqZsGQkv69kQ3pQbr18008EUymaGB\nVS1wmMkIOKdJzGXdQKTvpB3coUikrz/tD5AVyGTuHbhYIwdLR1P1ECKaD/caZvqFILZYc79jJjCE\nQ2XMBVWLhtchZ6CK8DjhsXUWKRav8W722+Tj+h+oZYkJuXeGHQopfXtr5UCQaDydP2hlgqcPQdlr\nIsBlH8xcI5F+xXEYCiy/G8dcVgsGGviCPZlpQn+DLxyuh1mwX8hkpC/Sk9ltmGvIZuP2wkxmt9km\nNTPkLCx7oURA3TMX6ccFyl7NfpucgSrC44TH1llkBBgn7UHHyJAViOlwu+M2xayBBmk20OFe9ID0\nZx3tiVwIyl4icytDq2Mk0hfNZCYxfTYNHO4SE2CXmYr0ixl7YU9mpxFE+rGYfqQv0pOptVwSBmjd\nEOwuW5TJ7HXapA1F+jkBwkCz1yFnIDAESPn2wkymNejgJJeZxWOBrIC0qjtqUzDA04egR7Kobtqb\ndCgbcBb366YL6ImDmau9TgHCuZXOwveZYu+sFm2hFQ6zZEd7VQsczskw//o63WEQ6ef1I/2ywJzM\nzoG+VkcIEfW6/U7HyLQ/iImPHbgdsoaWTYsMcZrqF0JQtVhU4u4M2xTTS2fxWCAnIEgULE4zFG0I\n3EC9WYdK3oy9xHRxk21odVjVnMgFKOaD/TYtb/785CTusl7St3ehuljzYTbzg0jfAHtno1RkZLXn\nvieM9E2gahcYzOY7p722uUi/kiviTeZfX911yWiuXw9RTJfoDOfbO+i2ySfMPAv/f3tnHiNHdt/3\nT/XdPd3T1ffMdHMOnjtccrm70h6yDq+OlbSSLNmJA8d/JE4MIZaFOIYQ6IoAW4LjE7IBCYGFGIkC\nRTEsJbITxYpie3WsInmj3WgPXsslh+TMcI6+j+rp+6r8UU3O9PTFJV9x2VR9AYI91dX17V+/qt/3\n/X7v9967mQ2JFEEpYNAmcWZLo5+9kqAswiRhYsXCbRt/A5VbRSF186ClvcYNcNdRCI6Y9/FaYO2M\njyyaJoUZn4gxGQmp4WUzUxh5XttaIBq4fb5Y0ItqK2p7bA9BoVQD1cyUw3bbfFG/TNM82rltZRXM\nLfm2uQBmZJn6GHFKKWKWagEIeWQq7dF82gROMXyyw8tOczRfvqoIKfaAbsHAmA2JRGYRXJbx1WXl\nVhGf04gsJgI3s8dEpSPOeTvN4ysy6hSF9LyhO4A/ZpCtZVaIBsQ4OEtLZjMz3AFok6xKRIMCVrl1\nWKFl12rxh2Azo2BqiLHtQFimYxsthPFcAaug6p05n0zDNJovqRRwSKPXcrpZRLwyVXU0n4jl168j\nMCVTaY/mU2oK04LSNDezIVGpWUQW1NN3W8ePkVQ6ipDii0nCxIrFzezDXVeLhL1ibtgp6/gB56ap\nyKyAnjeAwzR+g6COrUAsJIbP2pHZzg13AFuZIjTd2Ky3t6XqdZiaMhvp4XybmQIWUT19nxssNWqN\n1tBzEgUFuyqGLxr00rKM7nmndwrCSi9nZC91abTzzlULeKxixCnolqmMEadiQ0EWlAL2utxjS58r\nbQW/Swyfx+5GGZO1EOlbJgUTKxY3UztflxQhyxuAVnEybpCtbVWIBgRGMiOqvUrVBpiahLwuIXwO\n1UtSGe7gNjMK5qa4h8PaktnKDnc427kCVkGllyaThNSYZjM93L6kUhA2IDsfllFthZFptmxZwW0R\nI05zN5FmU+oFpm1i+ELTXuqM5iu3xDlv302UPmvLr4t59mSnl2J99JhTXVKYFZACniRMrFj4psbf\nQC1zkVm/oFDYMXqMRBuQLRILiQuFR1V7baQVpIZXSOklgNMkk1JGRBZZcT19AFtHJp4fzpdUCjgE\n9fQBzE2ZjRFjMpmSwpRZDJ/sdgDSyIKBfLWAV5DzPhCS6VjH9PTrBWE59oh3fJqt0i4ScIvhC7q9\n1MdUs9VQCE2L4fO7ZErN/MhzWmaFWb+RhpoI+G9ij4m2pUhMQI4dtEhmVC17fqcGHQtu5+0PyII2\ngD+q2mszXcDcFOdMp8wymdII510QuxaOQxodyaSKCi6TOPusbS/b2eF8uXIBj1Ucn6kxOs1WqBWE\nVepFg9OotqI2rjQEpaZCYEpcJNOyjBaLmqoQFuS8w155bJqtIRWFOe+wx0d5zJhM21okGjQii4lA\n0OOhPmJiV7PVAYE7WfmnPJRHTLTaSItbCwdAdsgo9eE3bDyvCEvTAHhsXnKV4c40oRRwSOKcqcsk\nk94Z1dMvCCu9BLCrMonCcL58rSBsQBbA3PKyNUKcdhoKfpeY39NmNUNziu3s8Puz3CkQ8ogSCy9t\n6+g0VF0SNwA845Vpmsf39EVU6kG3YIDh90qzpRV7zAVuf0XdScLEikVo2kNrhFho69s7sVnErMIe\n9ExTbQ8PhbezYpcsDrh8FBvDH5BEQWyaxmuXKdRGOO+iInQtnCnL6EgmVxHb03dIMqnicL5iXcHn\nFJhma49Os5XbBWFpGgBz08vWiGq2qlpgxitoAD80BeY69ebw5WiaZnE5/ajfR2tMmq1jLRILCxIn\neXQkE8+WoDmF1SKm2GNSMLFiEZZHb+S+lS5ibonr6Ue8PqoMd97xXFHIqpfXEXSPDoXTRQWHoAFZ\nAJ/Ty05juLPJCnbeXptMvjrcPqWmIDtEptm8ZEojevrNgrAJlQAORkcy1Y5CeFqcfZb26IKBulQQ\nMicHwGKRkOqjS63bFoU5QT39aEBGtQ9/9iq1FliqRG5zl7zriAV8IyOZzYwi1LdMCiZWLGb8Hjoj\nNo7fyirCZsgCzPl9NEzDb6CUomAXNEMWIDLto6IO50uXCsKqaQCCUzI7zRFpmqqCR9AMWQCvw4sy\nIpIpNgr4BFXTAHisMtnycL5KWyHoFvd7Ok0yqeJwZ1pDXE8fwD6mYKBlVogFRRYMeIeOyXQ6KqpN\nISYop+/zOAHYqQ4uGNhK7yA1xayzBRALyrRHRDLbWXEr6k4SJlcsZC/YFW3v3QGIF/LYOmLqygEO\nBH00LcOdd7KYxyX5hfHN+mTq0nC+XFnBLTCnH/LIVDrDnZtSKwjt6fud8shZwKVWgaCgAVmAabtM\noTY6TRMR6LynzF6yI9JsDXNBWKUeaAUDqREFA21rgQMhkQUDwyOZfKkKqhmv+/bX2YLuWml1mWvp\nwc+D6J5+NDgNtp2hBQPJQvGnbuMjmGCxcFod0LGQzFcGvp8o5HEiTizmwz46tlE9/Rxuizi+OZ9v\nZHmiVk0j7uEPe73URky02mmKq5sHbWJXqTWmp+8RmGZzeEcWDNSlAhGBk6ymbTK5ynC+tkVsT39U\nwUCt0RK2xel12JGJ5waL03a2iNQQ60wtTd/QSGYzk8faEtdRs1nN0PBoE1EHICk4izApmFixADA3\n/GxkcgPfS+3kmTKJu4FC024wN7R9CAYgW8kzLWiGLMCBkI+2dbg4KfWCsBmy0I1kRohTuV0QmqYJ\neUYvUVGjQEQWGMm4RqfZmgKXToHRBQOqqs2+F9nTd1u8ZIak2bYyRaSGB4tZ3OPulLwkh8zL2cwo\nWATn9G0dmc3M4Odhu5DDoYp79qA7L2dIJJPayQn1LZOCiRYLS9PHdn5wg2bKOTwCe/raLGCZ9dRg\nvkItj88hMpKRUe0K7c7gULjUElc3D9ogYmvELGDRA7IzskxtRHliw1RgzidWnEYttte2FoTl2AF8\nLpnikIIBbbKe2p28JwbTdplCdTDfhuA5OQBui0xqSCSzlc1j64h1pg58Q8dkkorYFDB002xDlr/J\nlHJ4LIZYTBRsHT/bucGRRb6ax2sX26CWpm9o76bYyBNwieNzOSzQdA7dcKnSLghN08RC3pGDejWp\nQEQWGMn4vSPTbE1LgTmBOf3QtHdoJKMNyBaFrbMFEHB5h6bZNtMKUkNGwA6nNyA7ZApD0mxbWQVr\nW6xYTNtkcuXB4rSVy+FE7LM3ZZZJKEN6+oJTwAD2znBxyungWyYBEy0WLvxs5weLhVIX29MHsLV9\nbOUG37A77RxBt1g+U8PH5pA8bVVViAjs6c/6PWCtUG+0B77fNCnMCuzpRwPyyNr5jlXhQFigfb7h\ntfOJfAlaThw2MXNyAELTMpXOsDSN2KVTAHwuLzuNYfYVsCOWz+vwDi19Tiji0zQei2/omEy2nGPa\nKjiSkWSSQ8SpUM/hdxhiMVFwW3xDG7TYyhGYEtzbUH0kCoP5Kp084WmxfNaWj43sYL664J6+xWyC\nppvN9OBBvZZFbJpmPiyjWhVUtb+arVJrgqUmrG4eNLFomoakaVIFoYskgjbruDZksb3tXEHY/tTX\nEXLLlIek2ZLFAk6Bs+9Bq2ZThmyAlNoR77y9DplcZfCzkBecAgYtkhkmTsVWjoBLLN8kYKLFwmvz\nky4PjixKrTwRj+BQ2OQnVRx8w9akPDNewZGMKrM9JJJpWrLMhwJC+cxNeeAGSJ2OimrPsRgR93v6\nPMMX27uWLiDVxS2SCKPTbNcyOSxNsffKjE+mMSSSiRfywgdkgx4v1SGRTHqnwJSg/amvIzCimi1b\nyeET3PP2OYcXDCiNHH6BKWDQIplh83JK7RyRaSOymCj4nH6ylcFiUVXzzPnFPpBus49UabDzbphz\nLITF3kAuyUdiSN60bc9yaFasWNhag6vLNjNFaDnxuMQskgha7bypIbOa6P89VxNZLE2xts2Htd35\nmq3+goFr6SwOVTSfPHRezlYhw5RJLF/U76MqDX4WMuUcXqtYvvC0TLk9ZLywlhM6fgcQcvtQhix/\nU2rlCQvuGMoOmdyQNFu1k2fOb4jFRCE45aNQG9LTN+WYD4m+gXxkhohF25pnaUawOFl8xAek2XLF\nKkhtIr4poXxONchGLtt3/Go8i7kh1tkAWFo+1pL99q2lstg7YvkcNgtSY5r1ZL8D2MpncQt23odm\n/XTsuYF7WiSLWWRbUCjfUiRIw9zfdgDZagbZIVicfEGq0mC+YjNH2C1YLDwy5SGRTEXNERGcAg66\nZfLV4b5lQbBvmQRMtFhEPH52WoN7Ny1LnoWI2BsoMOUbeAPVG21Ua4n5iNhQ32vzkd7p57sSz2Kq\nB4SmaQDcpgBbuUzf8dVkBntbvFg4OyGupfv5NrMZXIh1pgCWRpDL2/18cSWDxyLWPp/HCR2LtqDl\nPmTKWXyCnfeh2SBte79tAPlahohb7O85HwxSNw/mK3dyzHjFOtNZefjyN1XyHAgKFqdpeehCni1L\nTrhvmQToJhaSJP1rSZI6kiS4AHoPZnw+yp3+Bu10VDr2HAdnxTZoyOOjUO/nW08qSA0PVovYn9Pn\nksmWB6RpklmsTfHO1GsLkij2O4CNbBan4DQNwJQUZGOAOG0XssKdN4CjE2Qt1c+XLmfx2cXzmRtB\nVrbSfcfztSyhKbF8cwEPmBta1LkPxVaWWVks31IkSMs2WCwq5ISnaaIB39Dlb7QUsNhnfU4OUGr3\nR056+ZZJgC5iIUnSAeBJYF2P61/HgYCfGv2RRSJXgY6FaZe4SU8As14fpdagHLv4AVKAkNs/MM22\nkdHHefsdQVLlAT3vQha3Wbw4TVuDbCv9zjRVEp+mAXARZCPbb1+ukiXoEs9nbwVZTfbzKc0sM9Ni\n289kkjDVglyJ9zu4kpoh5hdr32JE1saA2v37mjdMOeE9/aVIkIZ1sDi1rXkOzonlWwiGqEj992Y8\nW4a2HY9LzLpXkwS9Ios/AT6h07VvYD7kp2HpF4ur8Rzmhnjljwb8VNR+vvVUTviMVYA5bwil1X/D\nbubED5AChN1BctX+BzK5kxU+QArgtwdJlfr5spUMAad4vmlLkO3CYOcdEey8AZwEuZbp5yt1Msz5\nxPPZWkGuJvr56lKWBcGVcw67GanmYy3Z/zy0rHmhlXMAh2ZCdBypvjGgYqkJ1oo2T0ggDs6EaFj6\nn73VRA5z46dvvAJ0EAtJkj4EbKqqekb0tfdjcWbw+kmryTS2Vlg431I4RM3cfwNtZNNMERLONx8I\nU1JTfccTij5pmtnpIEqz39lkyln8Ojjv0FRooDjl61nCbvF8si1Icqefb6edZdYrns9jCrKV7+er\nSVnmg+IjGUcnxPoAcWpYMxya1WEMqBlkZauXr9XuoNrFF3v4p53QsZJSelc0uBLPYWrImE1iok/L\nPAAAFJVJREFUXdnRaIi2I83+aUDrqRy29k9fCgpuUSwkSXpakqSzA/59EPg08Nt7TxfyTQcgGvCC\ntcxOuXfHrvV0Gpcq3nkfmQvTsvU7781CCo9ZL3Hq59Mrxx71Byl1BqRpalmCLvF8M9NBCo0BzruV\nZdYr3rkFXUEylX6+ClliAT3EKTRQnBqWLIsR8XxuU7CvQKHd1ubIHJ7TZwxofyRzLalAc0robPjr\nMNfCrGz3dtZWtlNYGxHhXCHvFKCSLvQWKKyl0zh18C2TgFtqUVVVnxx0XJKkE8AScFrSFr6JAS9I\nkvSoqvZ3kT/72c/eeP3EE0/wxBNPvKbvYTaZMNWCrGynefjI3I3jm4UUbpN4570Y8aPailRqTVwO\n643jyZ0Ufpt4vsOzYZq2/kgmV80S88aE8y2Gg9RM/c6t2MwIz7GDVn650+63r6xmiPrF80U8QVby\nr/Ydb5jEp2kAAs4g6Uq/fW1bVhfn7bUGiSu97beRVqDl6rlfRcEt9afZXt1MYq2Ld94A9naIle0U\nb73/0I1ja6kUzo74Z0+SJMz1EJe20oT3lKhv5pN4dPAto/DMM8/wzDPP3FHOQRAq/6qqngNu3CmS\nJK0Cb1DVAYl+esXiVmFtRFiJJ3rEIlFM4bOKb1CLWROnS1tpHjy0y5eupgg5Z4XzHZoNoNoK1Jst\n7NbdplKaWR52nxLOd3AmSMMyKMeeZU4H560NIg5I0+jkvOd8QXZW+/latozwCY4AEXeQFxO9NR5Z\npQqmVrfnKhaDChQub2ewNMRHaaCJ0/4020o8ibOjj1hMEWIt1Su+G7kUHkkfPlsrxNVkmrecWLxx\nLF5M4bPqwzcM+zvSn/vc5+4o/3XoPc9i8DZ2AjGlznA5nuw5liqnCDr1UX9rM8zKPr58I0XEIz40\ntdvMSHUfVxO9FS7FdoqYT4dIZi5Ax57tG0SsmlIshsTbN0ycmtYMB2fEO7j5QJCy2stXrjVQrSUW\nBefYAWblIMVWL9+lrTSmWkj4HBmA0FSQXK2Xby2Z1WWODIDf2T8GtJZO6ua8vZYQG9lesdguppB1\n6BgCuNQQ6+levlQ5SUAn33K3Q1exUFX14LCoQhS8lghrmV7nnammmfPqk1d0diKspnozakozzUJA\nJ3FqaKH3XpRIcHROfCQjux3QsRHPlXqO160J7p8Xz3d4Lkh7X61+td5Etee574D49lsI9U8kO7ea\nxFQLC90Y6Dpi/iClfeJ0YTOOoyX+twRtDEjZNwa0mkngQp97MzzVPwa0kU/q1vMOOMJ9pdbJcpKg\nUx8+jynEVqGXL11NEpXvbGRxt2CiZ3ADBB0RNguJnmOFZoqFoD4PyLQpzHq613mXSXF4Vh8+RyfM\n1UTvDVu3xTk+P6MLn6Ue4sK1XfGt1PRz3mF5Ckxt0oXdrXHPrScx1ULa1paCcSwWomnv7Vhc2Ejg\naOrjvA9GQtT21eqvxON40IdvIRSi2O69N9ezcfyWuSGfuD3MySHyjV6+xE6SkEsfZxrxhEhVevmy\n1RSzHn2ePdkeYnufWBRaSRaDhlhMJGY9MyRKvQ6gTIpDM/rcQAFHmGu5Xr66JcWxAzqJkznM1eTu\nA1JrtFDtOY7P6xR6t6Kc39i68fcr11KYakFdnLfJJGGpznL6SnyXbz2BXSfnvRCRwdQksSdyWknE\ncevkvE8sztCwx3uOrWXj+Kz68B2PzVE2bfcc2yjEibj0EYul4BxFtZcvU00yO63PvTnnDZGv9zrv\nfCPFQkgfvog7RGKnl09P33K3Y+LFIuaLkKv3Ou+aJcHRmD4NOuOJEN/Zdd6NZoeOI83yvD5pr4Aj\nzHp2174L19KY6gFdnDeA1xTlUnxXLC5s6Oe8AVztGOeubd74W3Pe+kRNJpOEtRrlpSu79q1lEvgs\n+vAtRnxgbhDP7orTdjFB2KXP7/ngwRhN52bPmFOysk3Mqw/fsdkYFctmz7FCK8kBvz4978XQDEq7\nN4uw09HPec/7ZkhXe8W+bklyX8yILCYSB8O9N1Cj2abjSnBqKaoL3wFfmEx1l++V9RRS04vbIXZp\nkeuYc8+xVdztvb1yLY69oZ/zDjuirOV2nelKXD/nDSCbYlxK7HXe+vW8AVztKOf3iNOWEtfNeV8X\np5f3iFOqEic6rc/vGQtNA5K2pHwX+WacpZA+9p1ciNJwbPVsYFUiyVJIH2e6PBejYtnoOVazbXFy\nQZ/I6djsPPnOLl+zpXUMjy8YkcVE4uhchIpp13mfv5bAVAvgsovbe2Evjs3Ok1ev3fj7zPoGjvoB\nXbgADgbnSdV3+VbiCaZ0dN6x6RjbO7vOdC2TQNap5w3XxWmXb1OJE9ahDPk6fOYYK8k9zruaYE4n\n5w0w1Y5xbmPXvlwzzmJQP/us1Sinr+7at0OcY3P6ONOFWQ90LGzs2TCrat3g1KI+z8MDSwdoODdu\niFO13qLjTPLAkj72nZw/QNm8Kxbn15NIDVm3juHdjokXi4cOztN0XbtxA51e3cDeED9h7TpOLSxQ\nsezWzl/c3sSj6sd33+w8hc6uWFxObeC36Me3FIySru86m/XCBjMu/fg0cdrlS1a2iE3r8/CDJk7r\neyKnbHODpaA+USiAbI5yKbErFjtsc2RGP7Fwd2K8srlrX926zf0L+vBJEthqUV64pPEVy3U6jgwP\nHtZrzGka2tYb61GdvrqNqRrGYRM/4RDg4cMxWq5N2h1tw6zTq9dw1Od14ZoETLxYRIPT0LbfWAbg\n7PomsqRfT/8Nhw/Qcm7TarcBuJjYIGzXj+/UwjwV265YXMmtMj+9pBvf0dkoirrrbDZLqxwO6Md3\nMBgj3dh1pqnmKvfPHdSNL+qNsl3a5VOkVR5e0s++iDPWI041+xqP36cfn2yOspLU7Nup1OnY89y/\noF+O3aPGOLuu8f1kZRNzdU6XpT5AEyd7/QDPX9R6+6fXNnA29Xv2gl4XUtPNxU3Nt5xZv4ZPWtCN\n727HxIuFJIGzvsSzr6wCcCm5QcSpX09Y9tiRagHOrGnjCGv5DQ7osPTGdTx0OErbnqLW0Na/2qqs\ncjS0qBvfA4u9eeFMe5UH5vVzbsvRGIq6y1c0X+WNh/XjOxyKka5rzq3TUak71/iZ4/rxHfDG2Cpp\n9q0nC6jmBsdi+syoBpidinE1o9n3o/OrWCsHemb/i0bAuhs5vXR1HXdLX2fqZZ7Ta9rvefbaNXwm\n/cQCwF6b56UrGt/FxDVmnEZkMdEISEu8uKqJxaqywiHfYV35ppqLPHteS0VtVVc4MXdENy63y4K5\nFuG5C1rvNNde49SCfs7tkaMxOvbsjfLSkmWVR47ox/ezJ45Qdlyi01FpNNu0nJu8aVk/h/P40SPk\npEtAtzih7SQaFLu89V48OH+EeF3je/bCKo7KQV1mb1/HcvgIq0WN77lLV5juHBrzidvDvOcwl3Mr\nAJzfWidg1teZRuzzXIivAXAhcZWYe1FXPp9pnhcua8/6urLGos8Qi4lGdGqJC4mrAMQbF3kgekxX\nvrD1IC+sXgYgb7rIY4f05fO2j/CD86/S6ahUHFd4/Jh+zttmNeOoHOG7L1+iUKrRcWR4wxH9cvpH\nYgEk1crZ1STPX9zEXAvjdeu3scw7Th2hMbVKpdbk7y9cwVnX77cEeMfJZRTbBQBevHoVGX353nR4\nmVRHWyzxzNYV5uz6isXJmWXWK5p9r6ReZWn6Pl35DsvLrBReAeBq8QLHw8u68i26l3l5W+Pbbl7g\nVFRfvrsZ94RYHA8f53LxHAAF66u885S+N+xR+SRnkmeo1FrUnWu8/ZS+kcy8/STPr53lzGoCkDix\nqG+dd4hlnr10gW89dx5H+ahuOejr8NSW+e6ZCzx9+hxy87iuXLLbgaUS4wdnr/DDi+eYM5/Qle/R\nYwfoWItspBV+cu08S2597Xvng8eouLSOxcXsBQ779e3IvP3kMllJE4urpXM8sqDv7/nY0gk26+cB\nSKkXePSgvs771Oz9XClqfAXreZ58UN/2u5txT4jF+x96I1vqTzi3lqJjKfP4sr55zDctPshq9WX+\n+rnz2KoL+Kf1LaU7GTnJhdxZvvX8GXz1U7qmMQAOTi/z4tY5vnv+DLOmB3TlApizLfPsyjl+dOUl\nDrsf1J3P31nmO2fO8VLiZU6Gxa/euxdms4Srssy3/995XlVe5rEFfe1bnJExNT386Nw6640XeOfy\nw7ryveuhwzSdmyTzZXKWc7zz5P268n3w8RMojrOUKk0qzou8/1F9xeKtx06QVM/yylqGjqXCY/fp\nNz55t+OeEIsPPHYfTcc2n//mtwhUHhe+a9Z+/PITD5N3vMTXn/s+87xZVy6AX3z8Mdb5IX/76g85\n6n5Ed74PnHwb50rf5+83fsjDkcd053ti6a3838T3OZt/ljcv6W/fG0Nv4+mVZ7jaeJb33K+/fUcd\nb+YbL3yflO1ZfuFR/e2Ltt7Knz79bcquC/ziW/QVJ6fdimfnET7+n/+cjqnKE6f0q2QDWJ4PY2nK\n/OaXv4qjcphY0Ksr38+/+Th1+xZ/8M2/IlD5Gd07ancz7gmxsFstRKtP8ZX8h3nb3FO68x2JhvBU\nTvDfqx/jA8vv0Z3vA4/eD6j8SP1D/snj79ed78PvfjMl5yusOL/KR598n+58H3n3u9ie+jZp9/f4\nzZ8buK+WUPzqW9/LWdufUbdv8StP6u+8/9GD7+M7rd/G2gzwtgf0HbMAeHLpKb5e/Ff4y28i4nPr\nzvem4FN8NfdRltpP6bJ6734ct36AL2c+zKmp9+rO5bLbmKu9i6/mP8Lb5vTnu6uhqurr8k+jFodv\nPfeKev8nPqLGsztCrzsM/+Fvfqw+9KmPqZVa447w/e7X/k59y299Rm23O3eE7zf+/dfUD/7+H98R\nLlVV1V/+4y+p//wLX75jfO/+nd9TP/Gf/uqOcLXbHfXxz3xK/fxffu+O8BXLdfWBT/5L9b9894U7\nwreeLKjLH/+I+p0XL98Rvp9c3FaXP/4v1HOryTvC9+3nXlXv/8Sv3zHfMg5d33nHfbak7t+R/A5B\nkiT19eI2YMCAgUmFJEmoqnrH82H3RBrKgAEDBgzoC0MsDBgwYMDAWBhiYcCAAQMGxsIQCwMGDBgw\nMBaGWBgwYMCAgbEwxMKAAQMGDIyFIRYGDBgwYGAsDLEwYMCAAQNjYYiFAQMGDBgYC0MsDBgwYMDA\nWBhiYcCAAQMGxsIQCwMGDBgwMBa6iIUkSb8hSdIFSZLOSZL0h3pwGDBgwICBOwfhYiFJ0tuBDwIP\nqKp6Avi8aI5JwDPPPPN6fwVdYdg32biX7buXbXs9oUdk8evA76uq2gRQVTWtA8ddj3v9hjXsm2zc\ny/bdy7a9ntBDLI4Ab5Mk6ceSJD0jSdIbdeAwYMCAAQN3EJZb+ZAkSU8DMwPe+kz3mj5VVR+XJOkR\n4L8C+m7Ma8CAAQMGdIXwnfIkSfrfwB+oqvqD7t+XgcdUVc3uO8/YJs+AAQMGbgGvx055txRZjMH/\nAN4B/ECSpKOAbb9QwOtjrAEDBgwYuDXoIRZfBr4sSdJZoAH8Ux04DBgwYMDAHYTwNJQBAwYMGLgH\noarqyH/Ae4FXgRXgk0PO+WL3/dPAQ+M+C/iBp4FLwN8B8p73Pt09/1Xg3XuO/y5wDdjZx20Hvt79\nzI+BhXE2TZh9/wxIAy91//3qpNkHOIH/BVwAzqGVVt8z7TfGvltuv7vBtu7xvwFe7tr2JcB0r7Td\nGPtuue3uJvv2vP8/gbO3+uyNM9YMXAYWAWv3B13ed877gG93Xz8G/HjcZ4E/Aj7Rff1JtAFxgOPd\n86zdz11mN/p5FK0Ca78z/Sjwp93XvwR87TU05iTY9yvAF1/LTXq32YfmTH+2e44V+D/Ae++V9htj\n3y21391iW/c99x7ObwC/dK+03Rj7Jv3ZM+3h+gfAnwNn9hx7Te03bp7Fo8BlVVXXVG2S3deAD+07\n54PAVwBUVX0OkCVJmhnz2Ruf6f7/893XHwL+QlXVpqqqa12DH+te+3lVVRMDvuPea/0l8M4xNk2a\nfVL3363grrBPVdWq2q2O617rRSA64FoT2X5j7LvV9rsrbOteuwQgSZIVsAGdAdeayLYbY9+kP3uP\ndu1yAx8D/u0+e15T+40TiyiwsefvTXYfgnHnzI34bERV1WT3dRKIdF/Pdc8bxTf0O6qq2gIUSZL8\nYz4z7rvfzDl3yj4V+IeSJJ2WJOm/SZIUG3P+zXz3mzlHF/skSZKBnwO+u5//Xmi/AfbdavvdVbZJ\nkvS33fOLaL3vHv5Jb7sh9k36szfXff07aMsuVYbx30z7jROLmx39vhn1lQZdT9VioFE8N/sdbgWT\nYN9fo+UST6HlKr8y5vzXcu3ruCP2SZJkAf4C+EK393O7mAT7brX97irbVFV9DzALOHhtEcTYa4/B\n62nfpD97kiRJDwIHVVX95k1yDcU4sdgCDuz5+wC96jXonFj3nEHHt7qvk91wC0mSZoHUiGttMRpb\nwHz3WhbAq6pqbsxnhn33u84+VVVz3VAU4D8Cbxh1/pjv/nrb92fARVVVv7iP/15pvz77bqP97jbb\nUFW1DnyT3ZTIvdR2ffbdA8/eJvA48EZJklaBHwJHJUn63p7P3Hz7jRrQQJuHcQVtwMTG+EGax9kd\npBn6WbRBmk92X3+K/kEaG7DU/by0j2/QAPeXuq//Ma9tkG0S7JvZ8/oXgGcn0T60fOk3Bth7T7Tf\nCPtuqf3uFtuAKWB2z3W/Dnz0Xmm7MfbdE8/eHr4FequhXlP73YzRTwEX0QZMPt099mvAr+055991\n3z8NPDzqs93jfuA7DC7/+jfd818F3rPn+B+h5dda3f9/q3vcjrb+1PXyr8WbbdAJse/30Er6XkbL\ngx+dNPvQejkd4Dz7yhDvhfYbY98tt99dYlsYeL57/bPAF+gtnZ30thtl38Q/e/u+zyK91VCvqf2M\nSXkGDBgwYGAsjG1VDRgwYMDAWBhiYcCAAQMGxsIQCwMGDBgwMBaGWBgwYMCAgbEwxMKAAQMGDIyF\nIRYGDBgwYGAsDLEwYMCAAQNjYYiFAQMGDBgYi/8PwONICjX2bmcAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f629b304b70>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.plot(ts, xs)\n",
"plt.plot(ts2, xs2)\n",
"#plt.xlim(.000,.00002)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"(0, 1000000.0)"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEACAYAAAB78OvLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHJxJREFUeJzt3W+QHPV95/H3R1pJSEJI6MCS+E/qJBe6M9iQQhzxxWtD\niI6ygQcukCvhlJzqUrZyhStOYiAPjHgQDFTOGD9AVDnEFqoYSxViwBWCkTm27lwpUMJhEBYyEqCc\ntKAVEVr9QUjsst970L/RtkazfzUz3b37eVV1zW9+3T3T3dLOZ779m+lRRGBmZtYMU4reADMzmzgc\nKmZm1jQOFTMzaxqHipmZNY1DxczMmsahYmZmTTOqUJG0U9Krkl6WtDn1zZe0SdIbkp6VNC+3/J2S\ntkvaJum6XP8VkrakeQ/m+mdI2pD6X5B0YTN30szM2mO0lUoAnRHxmYi4MvXdAWyKiCXAc+k+kpYC\ntwBLgeXAQ5KU1lkLrIqIxcBiSctT/ypgX+p/ALjvFPfLzMwKMJbTX6q7fwOwLrXXATel9o3AYxHR\nFxE7gR3AMkmLgDkRsTkt92hunfxjPQ5cM4btMjOzkhhLpfKspH+R9N9T34KI6EntHmBBap8D7M6t\nuxs4t0F/d+on3e4CiIh+4ICk+WPZETMzK17HKJf7rYh4V9LZwCZJ2/IzIyIk+XovZmaT3KhCJSLe\nTbfvSfoJcCXQI2lhROxJp7b2psW7gfNzq59HVqF0p3Z9f22dC4B3JHUAcyPi/fw2OLTMzMYnIuqH\nL1pmxNNfkmZJmpPas4HrgC3AU8DKtNhK4InUfgpYIWm6pIuBxcDmiNgDHJS0LA3c3wo8mVun9lhf\nJhv4P0lEeIrgrrvuKnwbyjL5WPhY+FgMP7XbaCqVBcBP0ge4OoC/jYhnJf0LsFHSKmAncDNARGyV\ntBHYCvQDq2Nwz1YDPwRmAk9HxDOp/xFgvaTtwD5gRRP2zczM2mzEUImIt4FPN+h/H7h2iHXuAe5p\n0P8S8KkG/cdIoWRmZtXlb9RXUGdnZ9GbUBo+FoN8LAb5WBRHRZxzGw9JUZVtNTMrC0lEmQbqzczM\nRsuhYmZmTeNQMTOzpnGomJlZ0zhUzMysaRwqZmbWNA4VMzNrGoeKmZk1jUPFzMyaxqFiZmZN41Ax\nM7OmcaiYmVnTOFTMzKxpHCpWaZs2gS9ebVYevvS9VdrMmfD227BwYdFbYlZOvvS92SgdPZpNvb1F\nb4mZ1ThUrLJqYeJQMSsPh4pV1v79J96aWfEcKlZZrlTMysehYpXlSsWsfBwqVlmuVMzKx6FileVK\nxax8HCpWWfv3w1lnuVIxKxOHilVWby9cfLErFbMycahYZe3fn4WKKxWz8nCoWGW5UjErH4eKVZYr\nFbPycahYZdUqFYeKWXk4VKyy9u+Hiy6CAwdgYKDorTEzcKhYhfX2Zh8pnjULDh8uemvMDBwqVlED\nA3DwIMydC2ee6cF6s7JwqFglHTwIp58OU6fCvHkeVzErC4eKVdL+/VmYgCsVszJxqFgl9fZmYQKu\nVMzKZFShImmqpJcl/TTdny9pk6Q3JD0raV5u2TslbZe0TdJ1uf4rJG1J8x7M9c+QtCH1vyDpwmbu\noE1M+Upl3jxXKmZlMdpK5evAViDS/TuATRGxBHgu3UfSUuAWYCmwHHhIktI6a4FVEbEYWCxpeepf\nBexL/Q8A953aLtlkkK9UzjzTlYpZWYwYKpLOA64H/hqoBcQNwLrUXgfclNo3Ao9FRF9E7AR2AMsk\nLQLmRMTmtNyjuXXyj/U4cM2498YmDVcqZuU0mkrlAeDPgfzXyxZERE9q9wALUvscYHduud3AuQ36\nu1M/6XYXQET0AwckzR/DPtgk5ErFrJw6hpsp6YvA3oh4WVJno2UiIiRFo3nNtmbNmuPtzs5OOjsb\nbpJNAq5UzBrr6uqiq6ursOcfNlSAq4EbJF0PnAacIWk90CNpYUTsSae29qblu4Hzc+ufR1ahdKd2\nfX9tnQuAdyR1AHMj4v1GG5MPFZvcenvhk5/M2q5UzAbVv+G+++672/r8w57+ioi/iIjzI+JiYAXw\nvyLiVuApYGVabCXwRGo/BayQNF3SxcBiYHNE7AEOSlqWBu5vBZ7MrVN7rC+TDfybDau+UnGomJXD\nSJVKvdpprnuBjZJWATuBmwEiYqukjWSfFOsHVkdEbZ3VwA+BmcDTEfFM6n8EWC9pO7CPLLzMhlU/\npuLTX2bloMHX/HKTFFXZVmu9q6+G+++Hz34Wdu+Gq67Kbs3sRJKICI28ZHP4G/VWSa5UzMrJoWKV\nlB9TmTULPvoom8ysWA4Vq6R8pSJ5sN6sLBwqVjlHj8LHH8PMmYN9/lixWTk4VKxyalWKckOP/gKk\nWTk4VKxy8uMpNa5UzMrBoWKVkx9PqXGlYlYODhWrHFcqZuXlULHKGapScaiYFc+hYpUzVKXi019m\nxXOoWOW4UjErL4eKVY4rFbPycqhY5bhSMSsvh4pVTqNKxR8pNisHh4pVTqNKxR8pNisHh4pVjisV\ns/JyqFjlDDWmcuAA+HfczIrlULHKaVSpTJsGM2bA4cPFbJOZZRwqVikDA3DwIMyde/I8j6uYFc+h\nYpVy8CDMng0dHSfP87iKWfEcKlYpjcZTalypmBXPoWKV0mg8pcZfgDQrnkPFKmWkSsWnv8yK5VCx\nSnGlYlZuDhWrlOEqFQ/UmxXPoWKVMlyl4oF6s+I5VKxSXKmYlZtDxSrFlYpZuTlUrFJcqZiVm0PF\nKsWVilm5OVSsUlypmJWbQ8UqxZWKWbk5VKxShqtUTj8dPvwQ+vrau01mNsihYpUyXKUiDf5Yl5kV\nw6FilXH0KHz8McyaNfQyvlSLWbGGDRVJp0l6UdIvJb0maU3qny9pk6Q3JD0raV5unTslbZe0TdJ1\nuf4rJG1J8x7M9c+QtCH1vyDpwhbsp00Avb1ZaEhDL+PBerNiDRsqEXEU+HxEfBr4NLBc0jLgDmBT\nRCwBnkv3kbQUuAVYCiwHHpKOvwSsBVZFxGJgsaTlqX8VsC/1PwDc18wdtIljuPGUGg/WmxVrxNNf\nEXEkNacD04AAbgDWpf51wE2pfSPwWET0RcROYAewTNIiYE5EbE7LPZpbJ/9YjwPXjHtvbEIbbjyl\nxpWKWbFGDBVJUyT9EugBnk3BsCAietIiPcCC1D4H2J1bfTdwboP+7tRPut0FEBH9wAFJ88e3OzaR\nuVIxK7/RVCoD6fTXeWRVx3+smx9k1YtZS7lSMSu/jtEuGBEHJD0P/C7QI2lhROxJp7b2psW6gfNz\nq51HVqF0p3Z9f22dC4B3JHUAcyPi/UbbsGbNmuPtzs5OOjs7R7v5NgG4UjEbWVdXF11dXYU9v7JC\nY4iZ0llAf0T0SpoJ/Ay4F+gkG1y/T9IdwLyIuCMN1P8IuJLstNbPgX8fESHpReA2YDPwD8D3IuIZ\nSauBT0XE1yStAG6KiBUNtiWG21ab+P7yL+HwYfj2t4deZu1aeOUVePjh9m2XWZlJIiKG+cxkc41U\nqSwC1kmaSnaqbENEPC3pBWCjpFXATuBmgIjYKmkjsBXoB1bnkmA18ENgJvB0RDyT+h8B1kvaDuwD\nTgoUM8gqkLPPHn4ZVypmxRo2VCJiC3B5g/73gWuHWOce4J4G/S8Bn2rQf4wUSmbD2b8fFi8efhmP\nqZgVy9+ot8rwmIpZ+TlUrDL86S+z8nOoWGWMplLxtb/MiuVQscoYbaXS2wv+oKBZMRwqVhmjqVRm\nzICODjhyZPjlzKw1HCpWCQMD2e+kzJ078rIerDcrjkPFKuHQIZg9O6tCRuLBerPiOFSsEkYznlLj\nSsWsOA4Vq4TRjKfUuFIxK45DxSrBlYpZNThUrBJcqZhVg0PFKsGVilk1OFSsElypmFWDQ8UqYSyV\nii/VYlYch4pVwlgqlTPPdKViVhSHilWCKxWzanCoWCWMtVJxqJgVw6FilTDWSsWnv8yK4VCxSnCl\nYlYNDhWrhLFUKnPmwAcfwMcft3abzOxkDhWrhLFUKlOmwBlnZJfKN7P2cqhY6R07Bn19MGvW6Nfx\nuIpZMRwqVnq1KkUa/ToeVzErhkPFSm8s4yk1rlTMiuFQsdIby3hKjb8AaVYMh4qV3ngqFV+qxawY\nDhUrPVcqZtXhULHSc6ViVh0OFSs9Vypm1eFQsdIbb6XiUDFrP4eKld54KxWf/jJrP4eKlZ4rFbPq\ncKhY6blSMasOh4qV3ni/Ue9Kxaz9HCpWeuOpVGofKY5ozTaZWWMOFSu98VQqp52WXYDy6NHWbJOZ\nNTZiqEg6X9Lzkn4l6TVJt6X++ZI2SXpD0rOS5uXWuVPSdknbJF2X679C0pY078Fc/wxJG1L/C5Iu\nbPaOWjUNDGS/izLWUAF/AdKsCKOpVPqAP4mI/wBcBfyxpEuAO4BNEbEEeC7dR9JS4BZgKbAceEg6\nftHytcCqiFgMLJa0PPWvAval/geA+5qyd1Z5hw5lv6PS0TH2dT2uYtZ+I4ZKROyJiF+m9mHgdeBc\n4AZgXVpsHXBTat8IPBYRfRGxE9gBLJO0CJgTEZvTco/m1sk/1uPANaeyUzZxjGc8pcaViln7jWlM\nRdJFwGeAF4EFEdGTZvUAC1L7HGB3brXdZCFU39+d+km3uwAioh84IGn+WLbNJqbxjKfUuFIxa79R\nh4qk08mqiK9HxKH8vIgIwJ+zsaZzpWJWLaM6Uy1pGlmgrI+IJ1J3j6SFEbEnndram/q7gfNzq59H\nVqF0p3Z9f22dC4B3JHUAcyPi/frtWLNmzfF2Z2cnnZ2do9l8q7D9+8cfKq5UbDLq6uqiq6ursOcf\nMVTSIPsjwNaI+G5u1lPASrJB9ZXAE7n+H0n6DtlprcXA5ogISQclLQM2A7cC36t7rBeAL5MN/J8k\nHyo2OfT2jv/0ly/VYpNR/Rvuu+++u63PP5pK5beA3wdelfRy6rsTuBfYKGkVsBO4GSAitkraCGwF\n+oHV6fQYwGrgh8BM4OmIeCb1PwKsl7Qd2AesOMX9sgniVCuVd99t7vaY2fBGDJWI+AVDj71cO8Q6\n9wD3NOh/CfhUg/5jpFAyyzuVSmXePHj99eZuj5kNz9+ot1I7lUrFA/Vm7edQsVI71UrFYypm7eVQ\nsVJzpWJWLQ4VKzVXKmbV4lCxUnOlYlYtDhUrtVOpVM44I7sg5cBAc7fJzIbmULFSO5VKZepUmDMn\nu3S+mbWHQ8VK69gx6OvLLn0/Xh5XMWsvh4qVVu1iksd/jWccPK5i1l4OFSutUxlPqXGlYtZeDhUr\nrVMZT6lxqJi1l0PFSqsZlYpPf5m1l0PFSsuViln1OFSstFypmFWPQ8VKy5WKWfU4VKy0XKmYVY9D\nxUrLlYpZ9ThUrLRcqZhVj0PFSsuViln1OFSstFypmFWPQ8VKy5WKWfU4VKy0mlGpzJwJH38MR482\nZ5vMbHgOFSulgYHsd1BONVSkrNpxtWLWHg4VK6XDh7Mqo6Pj1B/Lp8DM2sehYqXUjPGUGg/Wm7WP\nQ8VKqRnjKTWuVMzax6FipeRKxayaHCpWSq5UzKrJoWKl5ErFrJocKlZKrlTMqsmhYqXkSsWsmhwq\nVkquVMyqyaFipdTMSmXePFcqZu3iULFSamal4su0mLWPQ8VKqdmVikPFrD0cKlZKza5UfPrLrD1G\nDBVJfyOpR9KWXN98SZskvSHpWUnzcvPulLRd0jZJ1+X6r5C0Jc17MNc/Q9KG1P+CpAubuYNWTc2s\nVObOhYMHsysfm1lrjaZS+QGwvK7vDmBTRCwBnkv3kbQUuAVYmtZ5SJLSOmuBVRGxGFgsqfaYq4B9\nqf8B4L5T2B+bIJpZqXR0wKxZcOhQcx7PzIY2YqhExP8B6k8e3ACsS+11wE2pfSPwWET0RcROYAew\nTNIiYE5EbE7LPZpbJ/9YjwPXjGM/bAL56KNsmj27eY/pcRWz9hjvmMqCiOhJ7R5gQWqfA+zOLbcb\nOLdBf3fqJ93uAoiIfuCApPnj3C6bAGpVyvEatwk8rmLWHqf8E0gREZKiGRszkjVr1hxvd3Z20tnZ\n2Y6ntTZr5nhKjSsVmyy6urro6uoq7PnHGyo9khZGxJ50amtv6u8Gzs8tdx5ZhdKd2vX9tXUuAN6R\n1AHMjYj3Gz1pPlRs4mrmeEqNKxWbLOrfcN99991tff7xnv56CliZ2iuBJ3L9KyRNl3QxsBjYHBF7\ngIOSlqWB+1uBJxs81pfJBv5tEnOlYlZdI1Yqkh4DPgecJWkX8C3gXmCjpFXATuBmgIjYKmkjsBXo\nB1ZHRO3U2Grgh8BM4OmIeCb1PwKsl7Qd2AesaM6uWVW1olLxpVrM2mPEUImIrwwx69ohlr8HuKdB\n/0vApxr0HyOFkhm0plLxpVrM2sPfqLfScaViVl0OFSsdVypm1eVQsdJpVaXiUDFrPYeKlU6rKhWf\n/jJrPYeKlY4rFbPqcqhY6bhSMasuh4qVjisVs+pyqFjptKJSmT178OrHZtY6DhUrlQg4cCD7Ya1m\nklytmLWDQ8VK5dAhmDkTpk1r/mP7C5BmredQsVJpxXhKjb8AadZ6DhUrlVaMp9S4UjFrPYeKlYor\nFbNqc6hYqbhSMas2h4qViisVs2pzqFiptLpScaiYtZZDxUql1ZWKT3+ZtZZDxUrFlYpZtTlUrFRc\nqZhVm0PFSsWVilm1OVSsVFpZqfgjxWat51CxUmllpeKPFJu1nkPFSqXVlUpvb3YlZDNrDYeKlUor\nK5Vp0+C00+Dw4dY8vpk5VKxEaj+iNXt2657D4ypmreVQsdKonfqSWvccHlcxay2HipXGm2+2bjyl\nxpWKWWs5VKxwBw7AN78JX/oS/Omftva5fuM34Fvfguef94C9WSs4VKww/f3w0EPwyU/C++/Dli3w\n1a+29jm//334gz+Ar30NrroKfvITGBho7XOaTSaKirxdkxRV2VYb2T/+Y1aVLFoE3/kOXHZZe59/\nYACefBK+/W04eBBuvx1+7/dg+vT2bodZq0kiIlo4Uln3fFV5oXaoTAyvvQZ/9mfw9tvwV38FX/xi\nawfmRxIBXV1w773wq1/BN74Bf/RHcPrpxW2TWTO1O1R8+svaYu/e7JTTF74A11+fhcuXvlRsoED2\n/J//PPzsZ/DUU/Dii3DxxXDXXfBv/1bstplVkUPFWuroUbj/fli6NPvi4bZtcNtt2RcRy+byy2HD\nBvinf4I9e2DJkmxbX3oJjhwpeuvMqsGnv6ypIuBf/xVefRVeeQV+8AO49NIsWJYsKXrrxubdd+G7\n383Gf7Zvh3POgUsuyQKyNl1yCcyZU/SWmg3NYypDcKiUzwcfZKexXnllMERefTUbj7jssixMrr8e\nfvu3i97SU9ffD2+9Ba+/Dlu3Dk7btsH8+ScGzeLFsGABnH12Nm+KzwdYgSZtqEhaDnwXmAr8dUTc\nVzffodJmEdkno3p6stNBPT3Zi2gtQHbvzt6pX3rpYIhceimcdVbRW94+AwNZZZYPmh07sjGk996D\nQ4eyYPnEJ7KQOfvswXa+74wzsoqnNp12WvHjTTYxTMpQkTQV+DVwLdAN/DPwlYh4PbeMQyXp6uqi\ns7NzzOt99FEWEgcODN7u3ZuFRW2qhUdtmj49e9ddm5YsGQyRJUugo6P5+zcW4z0W7dLXlw34v/fe\nYNDUt997L/v3OHRocOrvPzFk6qdZs2DmzBOnXbu6uOyyzuP3TzuNE9ozZmT/ntOnn9iePn3iBVjZ\n/1+0U7tDpeCXhOOuBHZExE4AST8GbgReH26liSQieyE5dmxwOnoUPvwwGyTO365f38X27Z0nzTty\nJHtxygdHfvr4Y5g7N3tXXJvOPhsWLswC4zd/M7ut3V+wIHvxKrOyv3hMm5Z9F2fRorGt19d3YsjU\nplr4HDky+P/jww+zL4/+4hdd9PZ2Hu/78MPBZY4ezd5UHDs2eOHOWruvL9vO+qCZNq3x1NExdH/9\nNHXq0P21qf7+aKYpU4a//dGPupg7t5MpU06cl58a9Q03SSffn2hh3AxlCZVzgV25+7uBZfUL/frX\n2YvvaKaBgexFtHabn+r7avf7+0+c+vpO7stPtT/Ivr7B9nB9+cBoNE2Zkv1h56fZs7N3mrV3prNm\nZadX5s/P2rX+T3wia+cDozbVgmTGDP8RVMW0adm/8fz5o19nzZpsGquIwaCpD5tGU+1vo35q9DfU\n339y/5Ejg/21qf5+/dTff/Lf9HC3u3bB5s3ZvtX68tNQfbXXjuGm2jIRg8FSHzrDtfPL16873Pyh\n7g/VV+tvt7KEyqjOa9W+1zCaqfZOJP/updE7nfz9od5tNeqfOTN7sZ427cR3dbX2UH35wKidkqhN\nU6eO7mCN98XDrBFp8P/gRNGOv5H8G9j620btWmjVr9foMYZrN7o/VF8EXH11a49DvbKMqVwFrImI\n5en+ncBAfrBeUvEbamZWQZNxoL6DbKD+GuAdYDN1A/VmZlZ+pTj9FRH9kv4H8DOyjxQ/4kAxM6ue\nUlQqZmY2MZT+u76SlkvaJmm7pNuL3p7xknS+pOcl/UrSa5JuS/3zJW2S9IakZyXNy61zZ9rvbZKu\ny/VfIWlLmvdgrn+GpA2p/wVJF+bmrUzP8Yak/9qu/R6OpKmSXpb003R/Uh4LSfMk/Z2k1yVtlbRs\nEh+LP0l/H1sk/Sht+6Q4FpL+RlKPpC25vkL3XdLFkl5M6/xY0shX7YuI0k5kp8J2ABcB04BfApcU\nvV3j3JeFwKdT+3SyMaRLgPuBb6b+24F7U3tp2t9paf93MFhZbgauTO2ngeWpvRp4KLVvAX6c2vOB\nN4F5aXoTmFeCY/IN4G+Bp9L9SXksgHXAf0vtDmDuZDwWZF8teAuYke5vAFZOlmMB/GfgM8CWXF9R\n+z43zdsI3Jzaa4GvjrgfRf4xjeIg/yfgmdz9O4A7it6uJu3bE2RXENgGLEh9C4FtqX0ncHtu+WeA\nq4BFwOu5/hXAw7lllqV2B/Bean8FWJtb52FgRcH7fx7wc+DzwE9T36Q7FmQB8laD/sl4LM4F/h9w\nZtrOnwK/M5mOBVlA5EOlsH0HBLwHTEn9V5F7PR5qKvvpr0Zfijy3oG1pGkkXkb0jeZHsP0xPmtUD\nLEjtc8j2t6a27/X93Qwek+PHKyL6gQOS/t0wj1WkB4A/B/I/5jsZj8XFwHuSfiDp/0r6vqTZTMJj\nERHdwP8kC5Z3gN6I2MQkPBY5Re77fLJ/g4EGjzWksofKhPsUgaTTgceBr0fEofy8yN4OTLh9rifp\ni8DeiHiZ7N3QSSbLsSB7x3g52WmJy4EPyCry4ybLsZB0JnAD2bv1c4DTJf1+fpnJciwaafO+j/t5\nyh4q3cD5ufvnc2KiVkoa5HocWB8RT6TuHkkL0/xFwN7UX7/v55Hte3dq1/fX1rkgPVYH2XnRfQ0e\nq+jjeDVwg6S3gceAL0haz+Q8FruB3RHxz+n+35GFzJ5JeCyuBd6OiH3pnfTfk50Cn4zHoqaov4lu\n4H1gnqQpucfqHnGLizp/Osrzix1kg0YXAdOp9kC9gEeBB+r67yedGyV7h1o/EDed7BTJmwwOxL1I\ndm00cfJA3NoYPJeaH4h7i2wQ7sxau+hjkrbtcwyOqUzKYwH8b2BJaq9Jx2HSHQuyC8u+BsxM+7AO\n+OPJdCw4eUyl0H0nG6i/JbUfpuoD9WlH/gvZJ6V2AHcWvT2nsB+fJRs/+CXwcpqWp3/QnwNvAM/m\n/yMDf5H2exvwu7n+K4Atad73cv0z0n+C7cALwEW5eX+Y+rcDK4s+Hrnt+hyDn/6alMcCuIzs5x5e\nIXt3PncSH4s1ZFcn30IWKtMmy7Egq9rfAT4iG/v4w6L3nSywXkz9G4BpI+2Hv/xoZmZNU/YxFTMz\nqxCHipmZNY1DxczMmsahYmZmTeNQMTOzpnGomJlZ0zhUzMysaRwqZmbWNP8fMaCAI10dEGsAAAAA\nSUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f629b271d68>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.plot(np.fft.rfftfreq(len(ts2), np.mean(np.gradient(ts2))), np.abs(np.fft.rfft(xs2)))\n",
"plt.xlim(0,1e6)"
]
},
{
"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