Skip to content

Instantly share code, notes, and snippets.

@qpleple
Created May 8, 2018 17:57
Show Gist options
  • Save qpleple/5e06b1d73f99ad81394d08a8f5caa3f0 to your computer and use it in GitHub Desktop.
Save qpleple/5e06b1d73f99ad81394d08a8f5caa3f0 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 16,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import cPickle\n",
"from utils import *\n",
"%matplotlib inline\n",
"from matplotlib.pyplot import *"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Loading a tube"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"all_tubes = cPickle.load(open('VIRAT_S_000201_02_000590_000623.pkl'))[0][1]\n",
"\n",
"nms_ids = nms3dt(all_tubes, 0.3)\n",
"tubes = [all_tubes[t] for t in nms_ids]\n",
"\n",
"boxes, score = tubes[21]"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Computing the super boxes"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"def sample_boxes(boxes, sampling_interval):\n",
" idx = range(len(boxes[:,0]))[::sampling_interval]\n",
" \n",
" frames = np.array([boxes[:,0][i] for i in idx])\n",
" box_sampled = np.array([np.average(boxes[i:i+sampling_interval,1:], axis=0) for i in idx])\n",
" \n",
" return np.concatenate((frames[:,np.newaxis], box_sampled), axis=1)"
]
},
{
"cell_type": "code",
"execution_count": 50,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"super_boxes = sample_boxes(boxes, 15)\n",
"x = (super_boxes[:,1] + super_boxes[:,3]) / 2\n",
"y = - (super_boxes[:,2] + super_boxes[:,4]) / 2"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# The track of the tube"
]
},
{
"cell_type": "code",
"execution_count": 87,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEKCAYAAADTgGjXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4leWd//H392QnJIQsJJAQEpYAIijIWvcVptqRWtux\nnVZbO7WbnU6ntdV2fl2m02rrtLbOdDqXVltb9af8WmutVam41gUkAioBgkASICQhZCEQQrZz//44\nT+IhnEAIOTk5J5/XdXEl536ec879XI/w8X7uzZxziIiIDIYv0hUQEZHopRAREZFBU4iIiMigKURE\nRGTQFCIiIjJoChERERk0hYiIiAyaQkRERAZNISIiIoMWH+kKhFt2drYrKiqKdDVERKLGm2++ecA5\nlzOQc2M+RIqKiigtLY10NUREooaZVQ303Ig+zjKzr5qZM7Ns7/VFZnbQzDZ5f74ddO4KMys3sx1m\ndmvkai0iIj0i1hIxs8nAFcDuPof+5py7qs+5ccAvgMuBvcB6M3vCObdlWCorIiIhRbIlchfwdWAg\nywgvBnY453Y55zqAR4Crw1k5ERE5uYiEiJldDVQ7594KcXiZmb1lZk+b2RyvLB/YE3TOXq+sv8+/\nycxKzay0vr5+6CouIiLHCNvjLDNbA+SFOPQt4JsEHmX1tQGY4pw7bGbvBx4HZpzqdzvn7gHuAVi4\ncKE2TBERCZOwhYhz7rJQ5WY2FygG3jIzgAJgg5ktds7VBr3/KTP7H6/TvRqYHPQxBV6ZiEhMe3xj\nNXeuLmdfcxuTMlK4ZflMVs7v90HMKZ9/uoa9Y9059w4woee1mVUCC51zB8wsD6hzzjkzW0zgcVsD\n0AzMMLNiAuFxHfCx4a67iMhwenxjNbc99g5tnd0AVDe3cdtj7wDwlUc3hexQTkmIC3l+uIJkpM1Y\nvxbYbGZvAXcD17mALuBmYDWwFVjlnCuLYD1FRMLuztXlvYHQo62zm3/pJ0B6jvd9fefq8jDVcARM\nNnTOFQX9/t/Af/dz3lPAU8NULRGRiNvX3DaiPieUkdYSERERz6SMlBH1OaEoRERERqhbls8kJSHu\nmLK+r08mJSGOW5bPHMpqHUMhIiIyQq2cn8/t18wl32tJxPuMH37wTGyA748z4/Zr5oZ1dJZCRERk\nBFs5P59Xb72En37kLLr8jpy0ZCruuPKk70vwGT/5yFlhDRBQiIiIRIX3z51IVmoiD7xeCUDlHVdy\n7rTMkOdmpCRw54fDHyAwAkZniYjIySUnxHHd4sn88sWd7Gk8wuTMMTz0mWWRrpZaIiIi0eIfl0wB\n4KF1fRc/jxyFiIhIlJiUkcIVZ+Tx6PrdHO0zqTBSFCIiIlHkhvcV0XSkkz+/tS/SVQEUIiIiUWXp\n1ExKcsfywOuVOBf5RcoVIiIiUcTMuH5ZEZurW9i4pznS1VGIiIhEmw/OzyctKZ4HXquMdFUUIiIi\n0SY1KZ5rFxbw1Ds17D90NKJ1UYiIiEShTyydQme345E39pz85DBSiIiIRKGpOWO5oCSHh9ZV0dnt\nj1g9FCIiIlHqhmVTqGtp569ldRGrg0JERCRKXTRzApmpCXzl0U0U3/oXzr3jeR7fWD2sddDaWSIi\nUerPb+2jpa2LLn9gvshw7Knel1oiIiJR6s7V5b0B0iPce6r3pRAREYlCfr+jup+908O5p3pfepwl\nIhJFGg63s6p0Lw+/UdXvOeHcU70vhYiIyAjnnKO0qokH11bx9Du1dHT7WVKcycUzJ7CqdA9HO98b\n4hvuPdX7UoiIiIxQh4528seN1Ty0djfldYdIS47nY0sK+fjSQqZPSANgQeF47lxdzr7mNiZlpHDL\n8pnD1qkOChERkRGnbN9BHly7mz9tquZIRzdz88fxow/N5QNnTWJM4rH/bK+cnz+sodGXQkREZJg9\nvrH6uNbDijPzePLtGh5cW8WmPc0kJ/j4wLxJfHzpFM6anBHpKvfLRsJ69OG0cOFCV1paGulqiIgA\ngQC57bF3aAvamTDOZyTGGW2dfqbmpPLxJVP40IICxo1JiEgdzexN59zCgZyrloiIyDC6c3X5MQEC\n0O13EGc8/JklLJuahZlFqHanTiEiIhJm7V3dlFY28dL2+n7ndhzt9PO+adnDXLPTpxAREQmD3Q1H\neGn7fl4sr+e1nQ20dXaTGOcjKd5He9fxq+4O59yOoaQQEREZAm0d3azd1cBL2+t5aXs9FQdaASjM\nHMOHFxZwYUkOS6dm8eyWuuP6RIZ7bsdQUoiIiAyCc46d9Yd5sTwQGusqGuno8pOc4GPZ1CxuWDaF\nC2dOoChrzDF9HD3DcSM5t2MoKURERAbo0NFOXtvptTbK3+vfmD5hLJ9YOoULS3JYXJxJckLcCT8n\n0nM7hlJEQsTMvgt8Bqj3ir7pnHvKO3Yb8GmgG/hn59xqr3wF8HMgDviVc+6O4a63iIwuzjm21hzy\nHlHtp7SyiS6/IzUxjnOnZ/OFi6dxwYwcJmeOiXRVIyaSLZG7nHP/GVxgZmcA1wFzgEnAGjMr8Q7/\nArgc2AusN7MnnHNbhrPCIhI7Qk34Wzk/n+YjHbyy4wAveY+p9h9qB2D2xHT+6fypXDQzhwWF40mM\n1yLoMPIeZ10NPOKcawcqzGwHsNg7tsM5twvAzB7xzlWIiMgp6zvhr7q5ja/9v7e469ly9jS14XeQ\nnhzP+SU5XFSSwwUlOeSmJ0e41iNTJEPkZjO7HigFvuqcawLygbVB5+z1ygD29Clf0t8Hm9lNwE0A\nhYWFQ1lnEYlSHV1+9jQdofJAK995YvNxE/66/I59B49y8yUzuLAkh7MKxhEfp9bGyYQtRMxsDZAX\n4tC3gF8C3wec9/MnwI1D9d3OuXuAeyCw7MlQfa6IjGzdfkd1UxsVDa1U1B+msuEIFQdaqWxoZW9T\nW2Bm+Al0dTv+9fKSE54jxwpbiDjnLhvIeWZ2L/Ck97IamBx0uMAr4wTlIhJD+uur6OH3O2pbjlJ5\noJVdB1qp9EKi4kAruxuP0Nn9XlCkJsZRnJPK3Pxx/P1ZkyjKSqU4J5UvPrSBmoNHj/vuaJ3wF0mR\nGp010TlX4738ILDZ+/0J4GEz+ymBjvUZwBuAATPMrJhAeFwHfGx4ay0i4Raqr+KW37/Fn9/aR0Kc\nj8qGQGAEb8KUFO+jKCuVGRPSuPyMPIqzx1CcPZai7DHkjE0KuQ7VN1bMiqkJf5EUqT6RH5vZ2QQe\nZ1UCnwVwzpWZ2SoCHeZdwBedc90AZnYzsJrAEN/7nXNlkai4iAytw+1dlNe2sLXmED98autxfRWd\n3Y7ntu1nak4qU7NTOW96NkXZgd+LslPJS0/G5zu1BQtjbcJfJGkpeBEZFn6/o6rxCNtqWthae8j7\n2cKextALEgYzoOKOK8NfSQG0FLyIDJP++i8OtnVSXnuIbbUtbK0JtDLKaw/1tjJ8BsXZqcwryOC6\nRYXMyktj1sR0PvK/r1HdrL6KaKKWiIgMSqjNlXwG6ckJNLd19pZljElgdl46syamMTsvndkT05mR\nOzbk0iChPjMlIY7br5mrR03DSC0REQmbupajvL6zgW/+8Z3j+i/8Dtq7/Hxjxaze0MhND925HYr6\nKqKPQkRETqiptYO1uxp4bWcDr+08wM761hOef7Szm89fNG3Q3xdLixOOBgoRETnGoaOdrK9s5LUd\ngeDYWtuCc4E5F4uLM7luUSHLpmVx0+9K2af+i1FPISIyioTqCF9xZh5vVjXx2s4DvLazgbf3HqTb\n70iM97Fwyni+enkJy6ZlM69gHAlBy4B8fbnmWog61kVGjf46wg3odhDnM86enMH7pmWxbFoWCwrH\nn3RfjJPNLpfopI51ETnGgcPtfPfPZSE7wscmxfFfH13AouJMxiad2j8J6r8QhYhIDHLOUbavhee3\n7ef5bft5a28z/T10aG3v5uJZE4a3ghIzFCIiMaK1vYtXdhzghW37eaF8P3Ut7ZjBvIIM/uXSEh5c\nV0W9t8FSMHWEy+lQiIhEsaqGVp7bGgiNdbsa6ej2k5YUzwUlOVw8awIXzcwhe2wSAFOyxqgjXIac\nQkQkinR0+SmtbAw8pirfzy5vzsa0nFRueN8ULp41gUVFmceMouqhiXwSDgoRkRGiv5FO9YfaebE8\n0Np4efsBDrd3kRjnY8nUTK5fOoVLZuVSmDVmQN+hjnAZahriKzIChBp+G+8zJo5LZk9TYJXb3PQk\nLpk1gYtnTuDc6dmknuJIKpGB0hBfkSji97uQ+2h0eTv4ffXyEi6eNYE5k9IHvAaVyHBRiIgMs44u\nP+9UH2R9ZSOllY2sr2ziYNCqt8G6uh1funTGMNdQZOAUIiJhdri9i427m1hf0cgblY1s2tPcu73r\n1OxUVszJY/WWWpqPHB8kGn4rI51CRGQQTrTcx4HD7ZRWNvJGRRPrKxvZUtNCt9/hM5gzaRwfWzyF\nRUXjWViUSU5aYPjtso1ZGn4rUUkd6yKnKFQneEKccU7hePYfbu8ddpsU7+PsyRksLs5kUVEmC6aM\nP+GyIlqHSkaKU+lYV4iInALnHEt++Bz7Q8z8NuCSWRNY5IXGmfnpJMWfeAFDkZFIo7NEhohzjqqG\nI7y+q4G1uxp4fWdDyADpcd8nFw1j7UQiTyEiEsQ5x96mNl7f2dAbHDUHAxsv5aQlsXRqFn/bXn/M\nHuI91Akuo5FCREa9vU1HWLurkdd3BkKjujkwuS8rNZGl07JYOjWLZVOzmJaTipmF7BNRJ7iMVgoR\niTkn66CuOei1NHY2sLaigT2NgdAYPyaBpVOz+OyFU1k6NYsZE8aGnNynNahE3qOOdYkpoVoJyfE+\nPrywgM5ux+u7GqhqOALAuJQElhRnsszbya9kQho+n2aEi6hjXUatO1eXH7d8yNEuP79bu5u05HiW\nFGfyiaVTWDYti9l56QoNkdOkEJGoV93cxrpdDazb1djbn9GXAZu+fQVxCg2RIaUQkajinGNPYxtr\nKwKhsa6igb3eKrfpyfEkx/s42uU/7n2TMlIUICJhoBCRiDpZJ7hzjooDrayraAy0Nioae4fcjh+T\nwOLiTG48t5glUzOZlZfOn9/ap5FTIsNIISIR07cTvLq5jdsee5u6ljZSkxJYu6uBNyoaeyf3ZY9N\nZElxFkumZrKkODB6qm+fhkZOiQwvjc6SiDn3juf77cOAwCZMwaHRM09DRMJLo7NkRDvc3sXrOxtO\nGCAvfu0ipmSNUWiIjHARCREz+y7wGaDeK/qmc+4pMysCtgLlXvla59znvPecA/wGSAGeAr7sYr0Z\nFSOcc5TXHeLF8npeKq+ntKqRzm6HAaFuYH5GCkXZqcNdTREZhEi2RO5yzv1niPKdzrmzQ5T/kkDw\nrCMQIiuAp8NYPzkNB4908sqOA7y0fT8vba+nriXQrzErL40bzyvmwpIc9jW18X/+VKZOcJEoFhWP\ns8xsIpDunFvrvf4tsBKFyLDrbzSV3+/YvO8gL5XX8+L2ejbubsLvAsNuz5+Rw4Uzc7hgRg5545KP\n+bz4OJ86wUWiWEQ61r3HWZ8EWoBS4KvOuSbvcVYZsN079m/Oub+Z2ULgDufcZd77zwe+4Zy7qp/P\nvwm4CaCwsPCcqqqqsF7PaNHfZkzz8sdR0XCExtYOzGBe/jguLAkEx1kFGcTH+SJYaxE5VSOiY93M\n1gB5IQ59i8Cjqe8TeCT+feAnwI1ADVDonGvw+kAeN7M5p/rdzrl7gHsgMDprcFcgfYVaUqSz27Fx\nTzNXn53PhSU5nD8jm6yxSRGqoYgMt7CFSE+r4WTM7F7gSe897UC79/ubZrYTKAGqgYKgtxV4ZTIM\nWo528vQ7Nf2OpnIO7vqHUN1YIhLrIjU6a6JzrsZ7+UFgs1eeAzQ657rNbCowA9jlnGs0sxYzW0qg\nY/164L8iUffRotvveGXHAf7w5l5Wl9XS3uUn3md0+Y9v2GkzJpHRK1Id6z82s7MJPM6qBD7rlV8A\n/LuZdQJ+4HPOuUbv2Bd4b4jv06hTPSzerTvE7zfs5fGN1dS1tDMuJYGPLJzMh84poKL+MN/842aN\nphKRXhEJEefcJ/op/wPwh36OlQJnhrNesa6/kVWNrR08samaP2yo5p3qg8T7jItm5vDdDxRwyewJ\nJMXHAXD25AzMTKOpRKSXlj0ZJUKNrEqMM0py0yivO0Rnt2POpHQ+tKCAvz97EtnqHBcZtUbE6CwZ\nWUKNrOrodmypaeHT5xXzoXMKmJWXHqHaiUi0UoiMAkc7u084supbV54xzDUSkVihEIlhuxuO8NC6\nKlaV7un3HI2sEpHToRCJMd1+xwvb9vO7tVW8/G49PjOuOCOX4uxUfv1qBW2d7+36p5FVInK6FCIx\nov5QO6tK9/Dwut1UN7eRm57Ely+dwXWLCnvXqyrJTdPIKhEZUgqRKOacY31lEw+ureLpzTV0djvO\nnZ7Fv105m8vOyCWhz5pVK+fnKzREZEgpRKLQoaOdPL6xmgfX7qa87hBpyfF8YmkR/7i0kGk5YyNd\nPREZRRQiUWRbbQsPrq3ijxuqae3o5sz8dH70obl84KxJjEnUrRSR4ad/eUaQUDPK/25uHs9sruXB\ntVWsr2wiKd7HVfMm8YllUzirYJy2jxWRiNKM9REi1IzyeJ+RFO+jtaObKVlj+PiSKVx7TgHjUxMj\nWFMRiXWasR6FQs0o7/I74v2O3964mPOmZ+PzqdUhIiOLQmQE6Or29zujvL3LzwUlOcNcIxGRgVGI\nRFBbRze/f3MP9/6tot9zNKNcREayk4aImX0JeNA51zQM9RkVmlo7+O3rVTzweiWNrR2cPTmDS2dP\n4JE3dmtGuYhElYG0RHKB9Wa2AbgfWO1ivTc+TPY0HuG+Vyp4dP0e2jq7uXTWBD574TQWFY3HzDir\nIEMzykUkqgxodJYFxpFeAXwKWAisAu5zzu0Mb/VO33CPzgo1THdG7ljueXkXT75dg8/g6rPzuemC\nqZTkpg1bvUREBmrIR2c555yZ1QK1QBcwHvi9mT3rnPv64KsaW/oO061ubuNfV23C72BsUjyfPq+Y\nT51bxMRx6ucQkdgwkD6RLwPXAweAXwG3OOc6zcwHvAsoRDyhhun6HaQnx/O3b1zCuJSECNVMRCQ8\nBtISyQSucc5VBRc65/xmdlV4qhWd9vUzTPfQ0S4FiIjEpJOGiHPuOyc4tnVoqxO9Orr8JCfEHdcS\nAQ3TFZHY5Tv5KXIyh9u7+PQD62nr7Ca+z6xyDdMVkVimyYan6cDhdj716/VsqWnhzmvnkRDn0zBd\nERk1FCKnoaqhlevvf4O6lqPce/05XDIrF0ChISKjhkJkkDZXH+STv36DLr/joX9ayjlTxke6SiIi\nw04hMgiv7jjAZ3/3JunJ8Txy02KmT9CkQREZnRQip+iJt/bx1VWbmJo9lgduXEzeuORIV0lEJGIU\nIqfg/lcq+Pcnt7C4KJN7r1/IuDGa+yEio5tC5ASC18FKTYrncHsXV5yRy90fnU9yQlykqyciEnEK\nkX70XQfrcHsXcT5jxZw8BYiIiEeTDfsRah2sbr/jJ89uj1CNRERGnoiFiJl9ycy2mVmZmf04qPw2\nM9thZuVmtjyofIVXtsPMbg13/fpbB6u/chGR0Sgij7PM7GLgauAs51y7mU3wys8ArgPmAJOANWZW\n4r3tF8DlwF4Cm2Q94ZzbEq46TspICbnvudbBEhF5T6RaIp8H7nDOtQM45/Z75VcDjzjn2p1zFcAO\nYLH3Z4dzbpdzrgN4xDs3bG5ZPpOEOK2DJSJyIpEKkRLgfDNbZ2Yvmdkirzwf2BN03l6vrL/ysFk5\nP5+lU7MwwID8jBRuv2auljQREQkStsdZZrYGyAtx6Fve92YCS4FFwCozmzqE330TcBNAYWHhoD+n\ntb2LRUWZrPrcsqGqmohITAlbiDjnLuvvmJl9HnjMBTZ4f8PM/EA2UA1MDjq1wCvjBOWhvvse4B4I\n7LE+mPp3dfsp29fCx5dOGczbRURGhUg9znocuBjA6zhPJLD97hPAdWaWZGbFwAzgDWA9MMPMis0s\nkUDn+xPhrOC7+w/T3uVnXsG4cH6NiEhUi9Rkw/uB+81sM9AB3OC1SsrMbBWwBegCvuic6wYws5uB\n1UAccL9zrixclXt8YzXffSLw8T/4y1ac0/LuIiKhWODf7ti1cOFCV1paOuDz+85Uh8CoLHWqi8ho\nYWZvOucWDuRczVjvI9RM9bbObu5cXR6hGomIjFwKkT40U11EZOAUIn30NyNdM9VFRI6nEOnjluUz\nSemzSq9mqouIhKal4Pvo6Ty/c3U51c1tGPCdD8xWp7qISAhqiYSwcn4+r956CU/cfC4OaO3wR7pK\nIiIjkkLkBOYVZLC4KJNfv1pBV7eCRESkL4XISfzT+cXsbWpjdVldpKsiIjLiKERO4tLZuRRljeFX\nr+yKdFVEREYchchJxPmMG88rZuPuZt6saop0dURERhSFyABce04B41ISuE+tERGRYyhEBmBMYjwf\nW1LIM5tr2dN4JNLVEREZMRQiA3TDsiKcgxU/f5niW//CuXc8z+Mb+93SRERkVNBkwwFau6sBnxmt\n7YHFGaub27jtsXcALRMvIqOXWiIDdOfqcrr7LJuv1X1FZLRTiAyQVvcVETmeQmSAtLqviMjxFCID\nFGp1X4CztAe7iIxiCpEBWjk/n9uvmUt+RgoGTMpIZuGUDJ7aXMvP17wb6eqJiESERmedgpXz848Z\nidXtd3z9929z15rtOBz/cllJBGsnIjL8FCKnIc5n/PjaeZjBz9a8i3PwlcsVJCIyeihETlOcz/jR\nh+ZhwM+fCzzWUpCIyGihEBkCvUFigSBxwFcum4GZRbpqIiJhpRAZIj6fccc18zCMu597F5zjK5eX\nKEhEJKYpRIaQz2fcfs1czODu53fggH9VkIhIDFOIDDGfz/jhBwNB8l/P78DvHF+7YqaCRERikkIk\nDHw+4wcr5wLGL17YSXntIbbWtLCv+SiTMlK4ZflMLdooIjFBIRImgSA5k8qGw6zZur+3XKv/ikgs\n0Yz1MPL5jKqG4zex0uq/IhIrFCJhVtN8NGS5Vv8VkVigEAmz/lb5jfMZz2+rw/XZo0REJJooRMIs\n1Oq/iXE+xo9J4MbflHL9/W+wve5QhGonInJ6IhYiZvYlM9tmZmVm9mOvrMjM2sxsk/fnf4POP8fM\n3jGzHWZ2t0XJmNm+q//mZ6Tw42vn8eqtl/Ltq87grT3NrPjZy/yfxzfT2NoR6eqKiJwSi8TjFDO7\nGPgWcKVzrt3MJjjn9ptZEfCkc+7MEO95A/hnYB3wFHC3c+7pk33XwoULXWlp6ZDWfyg1tXbwszXb\neXDdbsYkxvHlS2dw/bIiEuPVSBSRyDCzN51zCwdybqT+pfo8cIdzrh3AObf/RCeb2UQg3Tm31gVS\n77fAyvBXM/zGpybyvavP5Jkvn8+CwvH8x1+2csVdL/HsFvWXiMjIF6kQKQHON7N1ZvaSmS0KOlZs\nZhu98vO9snxgb9A5e72ymDEjN40HblzMrz+1iPg4H5/5bSkfv28dW2taIl01EZF+hW2yoZmtAfJC\nHPqW972ZwFJgEbDKzKYCNUChc67BzM4BHjezOYP47puAmwAKCwsHeQWRcfHMCZw3PZuH1+3mrjXb\nufLuv/EPiwr56hUlZI9NinT1RESOEak+kWeAHznnXvBe7wSWOufq+5z3IvA1oBp4wTk3yyv/KHCR\nc+6zJ/uukd4nciLNRzr4+XPv8rvXq0hOiONLl0znk+cWkRR//F7vIiJDJRr6RB4HLgYwsxIgEThg\nZjlmFueVTwVmALucczVAi5kt9UZlXQ/8KTJVHz4ZYxL5zgfmsPorF7CkOJPbn97G5T99mWc216i/\nRERGhEiFyP3AVDPbDDwC3OB1mF8AvG1mm4DfA59zzjV67/kC8CtgB7ATOOnIrFgxLWcs931yEb+9\ncTHJCT4+9+AGrrtnLZurD0a6aiIyykXkcdZwiubHWaF0dft5ZP0efvrsdpqOdPDhcwr42vKZTEhL\njnTVRCRGnMrjLIVIlDrY1sl/P/8uv3mtksQ4H1+4eDoT0pL42Zp32dfcpiXnRWTQFCJBYjVEelQc\naOWHT23l2S11GBB8N1MS4rj9mrkKEhE5JdHQsS5DpDg7lXuvX0j22ET6/u+AlpwXkXBTiMSIhsOh\n192qbm7j169WaOl5EQkL7WwYIyZlpFAdIijifcb3/ryF7/15C/MKxrF8Th7L5+QxfcLYCNRSRGKN\n+kRixOMbq7ntsXdo6+zuLevpE5lXMI7VZXU8U1bLW3uaAZg+YSzL5+SyfE4ec/PHESWLIovIMFDH\nepDREiIQCJI7V5efcHRWzcE2/lpWx+qyWtZVNNLtd+RnpHD5GbmsODOPRUWZxPkUKCKjmUIkyGgK\nkVPV1NrBmq11rC6r4+V36+no8pOZmsjls3NZfmYu507P1hIrIqOQQiSIQmRgWtu7eGl7Pc9sruWF\nbfs51N7F2KR4LpqZw4oz87ho5gTGJqkLTWQ0OJUQ0b8KAkBqUjzvnzuR98+dSHtXN6/tbOCvZbX8\ntayOJ9+uITHex/nTs1k+J4/LzsglMzUx0lUWkRFALRE5oW6/482qJp7ZXMvqslqqm9vwGSwuzmTF\nnDyumJPHpIyUSFdTRIaQHmcFUYgMHeccZftaWF0WCJTtdYcBOKtgHFfMyWPFmXlMywkMHR5IJ7+I\njEwKkSAKkfDZVX845NDh4qwxvPzuAdq7/L3nagkWkeihEAmiEBkewUOHX9vZEPKc/IwUXr31kmGu\nmYicKq2dJcNu4rgUbnhfEQ9/Zin9zTKpbm7jiw9v4L5XKti4u4mOoJaKiEQnjc6SIdffEiwpCT42\n7W7mL2/XAJAY72Nu/jgWFGawoHA8C6aMJzdd+6KIRBOFiAy5W5bP7HcJlpXz86lrOcqGqiY27G5i\nw+5mHni9inv/VgEEHnnNDwqVMyamkxivBrPISKU+EQmLUxmd1d7VzZZ9LWzY3cyG3U1srGpi38Gj\nACT1tFamjO9tsUxQa0UkrNSxHkQhEp1qDx4NtFS8Fsvm6hY6ugN9KPkZKceEymy1VkSGlEIkiEIk\nNrR3dVNMKlZyAAAJn0lEQVS2r4UNVU1s9FosNUGtlXkF41hQOJ75heNZMCVDe86LnAaFSBCFSOyq\nOdjGhqpmr2+libKg1krB+JRAv0phBgumBForCXFqrYgMhEIkiEJk9DjaGWitbPRCZUNVM7UtgdZK\ncoKPefkZzJ/iddoXjicnLan3vZphL/IehUgQhcjotq+5rTdQNuxuomzfQTq7A//NT84MtFbifMZf\n3q7RDHsRj1bxFfFMykhhUkYKV82bBPS0Vg72hsraXQ3UtbQf9762zm5+8JetXH5GLqlaAl+kX2qJ\nyKjmnGPqbU9xor8FhZljmJWXxqyJ6cz2fhZmjtEOkBKz1BIRGSAz63eGfWZqIp96XxHbag+xtbaF\nNVvr8Htpk5IQR0nuWGblpTNrYlrgZ14a47XPiowyChEZ9fqbYf/tq844pk/kaGc379YdZmttC9tq\nDrGttoVnt9bxaOme3nPy0pOZmZfGrIlpzPYCZmr2WM1jkZilEJFRrycoTjY6KzkhjrkF45hbMK63\nzDlH/eH23lDZVnOIrbWHeG3ngd4O/IQ4Y1rO2N5HYrPy0pg9MZ0JaUmY6ZGYRDf1iYiEQWe3n4oD\nrWytaWFb7SG2eT97JkgCjB+TEGi15KUz23skVpKbRkpiXARrLqI+EZGIS4jzUZKbRkluGlcHlR88\n0hlosdQGWi5baw7x6Po9vY/SzKA4K7U3XHoeixWMT8GnjnwZgRQiIsNo3JgElkzNYsnUrN4yv9+x\nu/FIb7BsqznE1poWnimrpedBQWpiHDPz0pgZ1GqZmZfGuJSECF2JSIAeZ4mMUEc6uthed7j3UVjP\no7GDbZ295+RnpDArL83rzA8MQS7OTiW+zxIvmpEvp2LEP84ys0eBmd7LDKDZOXe2d+w24NNAN/DP\nzrnVXvkK4OdAHPAr59wdw15xkWE0JjGesydncPbkjN4y5xy1LUe9fpb3Wi4vba+nyxt/nBjvY8aE\nsb3Djhta2/nNq5Uc9WbkVze3cdtj7wAoSOS0RbwlYmY/AQ465/7dzM4A/i+wGJgErAFKvFO3A5cD\ne4H1wEedc1tO9vlqicho0N7Vza761mNGiG2raWH/oeNn4/fQnvfSnxHfEulhgfGNHwF6/ku+GnjE\nOdcOVJjZDgKBArDDObfLe98j3rknDRGR0SApPo7ZE9OZPTEd5r9X3tjawYLvPxvyPftCTLAUOVWR\nngF1PlDnnHvXe50P7Ak6vtcr669cRE4gMzWR/IyUkMcm9VMucirCFiJmtsbMNof4Ezzi8aMEHl8N\n9XffZGalZlZaX18/1B8vElVuWT6TlIRj556kJMRxy/KZ/bxDZODC9jjLOXfZiY6bWTxwDXBOUHE1\nMDnodYFXxgnKQ333PcA9EOgTGXitRWLPQGfkiwxGJPtELgO2Oef2BpU9ATxsZj8l0LE+A3gDMGCG\nmRUTCI/rgI8Nc31FotbK+fkKDQmLSIbIdfR5lOWcKzOzVQQ6zLuALzrnugHM7GZgNYEhvvc758qG\nub4iItJHxIf4hpuG+IqInJpTGeIb6dFZIiISxRQiIiIyaAoREREZtJjvEzGzeqBqEG/NBg4McXVG\nGl1jbNA1xoaRdI1TnHM5Azkx5kNksMysdKAdS9FK1xgbdI2xIVqvUY+zRERk0BQiIiIyaAqR/t0T\n6QoMA11jbNA1xoaovEb1iYiIyKCpJSIiIoM2KkPEzJLN7A0ze8vMyszse155sZmtM7MdZvaomSV6\n5Une6x3e8aJI1n8gTnCNvzGzCjPb5P3p2ZbYzOxu7xrfNrMFkb2CgTOzODPbaGZPeq9j5j72CHGN\nMXUfzazSzN7xrqXUK8s0s2fN7F3v53ivPCqvEfq9zu+aWXXQvXx/0Pm3eddZbmbLI1fz/o3KEAHa\ngUucc2cBZwMrzGwp8CPgLufcdKCJwF7veD+bvPK7vPNGuv6uEeAW59zZ3p9NXtnfEVg1eQZwE/DL\nYa/x4H0Z2Br0OpbuY4++1wixdx8v9q6lZ5jrrcBzzrkZwHPea4jua4TjrxMC/7323MunALztwq8D\n5gArgP8xs7gQnxdRozJEXMBh72WC98cR2Kb39175A8BK7/ervdd4xy/1tvYdsU5wjf25Gvit9761\nQIaZTQx3PU+XmRUAVwK/8l4bMXQf4fhrPImovI/9CL5ffe9jrFzjifRuF+6cqwCCtwsfMUZliEDv\n44FNwH7gWWAn0Oyc6/JOCd6Ct3d7Xu/4QSBreGt86vpeo3NunXfoB95jgLvMLMkri9YtiH8GfB3w\ne6+ziLH7yPHX2COW7qMD/mpmb5rZTV5ZrnOuxvu9Fsj1fo/Wa4TQ1wlws3cv7+95bEeUXOeoDRHn\nXLdz7mwCuyQuBmZFuEpDru81mtmZwG0ErnURkAl8I4JVPC1mdhWw3zn3ZqTrEi4nuMaYuY+e85xz\nCwg8qvqimV0QfNAFhpHGwlDSUNf5S2AagcfONcBPIli/UzZqQ6SHc64ZeAFYRqBZ3LNRV/AWvL3b\n9nrHxwENw1zVQQu6xhXOuRrvMUA78Gveax6faGvikepc4O/NrBJ4hMBjrJ8TW/fxuGs0swdj7D7i\nnKv2fu4H/kjgeup6HlN5P/d7p0flNULo63TO1Xn/w+cH7iXK7uWoDBEzyzGzDO/3FOByAp2WLwDX\neqfdAPzJ+/0J7zXe8efdCJ9g0881bgv6S2kEnjFv9t7yBHC9N/JlKXAw6FHCiOScu805V+CcKyLQ\nAfm8c+4fiaH72M81fjyW7qOZpZpZWs/vwBUErif4fvW9j1F1jdD/dfbpz/kgx97L6ywwqrCY97YL\nH1EiuT1uJE0EHvBGOviAVc65J81sC/CImf0HsBG4zzv/PuB3ZrYDaCTwl3mk6+8anzezHAL71m8C\nPued/xTwfgKdd0eAT0WgzkPlG8TOfezPQzF0H3OBP3pjHOKBh51zz5jZemCVmX2awErcH/HOj8Zr\nhP6v83cWGKLtgErgs3Di7cJHEs1YFxGRQRuVj7NERGRoKERERGTQFCIiIjJoChERERk0hYiIiAya\nQkRERAZNISIiIoOmEBEZRma2yFtoL9mbwVzmrWkmEpU02VBkmHkz6ZOBFGCvc+72CFdJZNAUIiLD\nzAI7La4HjgLvG4lLWYgMlB5niQy/LGAskEagRSIStdQSERlmZvYEgWXdi4GJzrmbI1wlkUEbrav4\nikSEmV0PdDrnHvZWWH7NzC5xzj0f6bqJDIZaIiIiMmjqExERkUFTiIiIyKApREREZNAUIiIiMmgK\nERERGTSFiIiIDJpCREREBk0hIiIig/b/AWtBIVLcquxnAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x110498390>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plot(x, y, 'o-')\n",
"xlabel('x')\n",
"ylabel('y')\n",
"show()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Fitting the 3 models on the 4 first super centers\n",
"\n",
"See documentation of [Numpy Polyfit](https://docs.scipy.org/doc/numpy-1.14.0/reference/generated/numpy.polyfit.html)"
]
},
{
"cell_type": "code",
"execution_count": 71,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"p2 = np.poly1d(np.polyfit(x[:4], y[:4], 2))\n",
"p1 = np.poly1d(np.polyfit(x[:4], y[:4], 1))\n",
"p0 = np.mean(x[:4]), np.mean(y[:4])"
]
},
{
"cell_type": "code",
"execution_count": 81,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x10f720990>"
]
},
"execution_count": 81,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAfIAAAHVCAYAAAAKOHleAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Wd0VdXahuF77vQQSICE3nuvoQuCqIACioKgVBGQquIB\nj4ocFT2K4GdBEVR6kWJDQDqIHKSDoUhXAoRiCJ2E9Pn9SEAQhPSd8lxjZIzsVeZ+AsqbufZa8zXW\nWkRERCRrcjg7gIiIiKScCrmIiEgWpkIuIiKShamQi4iIZGEq5CIiIlmYCrmIiEgWpkIuIiKShamQ\ni4iIZGEq5CIiIlmYq7MDJJW/v78tVaqUs2OIiIhkiO3bt4dZawPudlyWKeSlSpVi27Ztzo4hIiKS\nIYwxR5NynC6ti4iIZGEq5CIiIlmYCrmIiEgWlmU+IxcRkdSLiYkhJCSEyMhIZ0eRRJ6enhQrVgw3\nN7cUna9CLiKSg4SEhJA7d25KlSqFMcbZcXI8ay1nz54lJCSE0qVLp2gMXVoXEclBIiMjyZ8/v4p4\nJmGMIX/+/Km6QqJCLiKSw6iIZy6p/ftQIRcREcnCVMhFROQfLfj1BE1Gr6H0yz/SZPQaFvx6wtmR\nUuV///sfVatWpVatWly9etXZcdKECrmIiNzWgl9P8Mp3uzlx4SoWOHHhKq98tzvLFHNrLfHx8Tdt\nmz17Nq+88gpBQUF4eXld3x4bG5vR8dKM7loXEcmh3lz0G3tPXvrH/b8eu0B03M2F8GpMHC99s4s5\nW47d9pwqRfLweruq/zhmeHg4TzzxBCEhIcTFxTFy5Eg6d+58fRluf39/tm3bxrBhw1i7di1vvPEG\nv//+O4cPHyYsLIyXXnqJvn37AjB27Fjmz59PVFQUHTp04M033yQ4OJhWrVrRoEEDtm/fzpIlSyhZ\nsiQAkyZNYv78+SxfvpylS5fSt29fRo4cSd68edm/fz8HDx7k0Ucf5fjx40RGRvL888/Tr18/AHx8\nfBgwYABLliyhcOHCvPPOO7z00kscO3aMjz76iPbt2xMXF8fLL7/M2rVriYqKYtCgQTz77LOcOnWK\nzp07c+nSJWJjY5kwYQJNmzZN1t/VnaiQi4jIbf29iN9te1IsW7aMIkWK8OOPPwJw8eLFu56za9cu\nNm3aRHh4OLVr1+bhhx9mz549HDp0iC1btmCtpX379qxbt44SJUpw6NAhpk+fTsOGDW8ap0+fPqxf\nv562bdvSsWNH1q5dy44dO9izZ8/1R7+mTJlCvnz5uHr1KvXq1ePxxx8nf/78hIeHc9999zF27Fg6\ndOjAa6+9xsqVK9m7dy89e/akffv2TJ48GV9fX7Zu3UpUVBRNmjThwQcf5LvvvqNVq1aMGDGCuLg4\nIiIiUvzndzsq5CIiOdSdZs4ATUav4cSFWz9HLurnxbxnG6XoPatXr86//vUv/v3vf9O2bdskzUwf\neeQRvLy88PLyokWLFmzZsoX169ezYsUKateuDcCVK1c4dOgQJUqUoGTJkrcU8X9Sv379m57fHjdu\nHN9//z0Ax48f59ChQ+TPnx93d3dat259/Wfw8PDAzc2N6tWrExwcDMCKFSvYtWsX33zzDZDwS8qh\nQ4eoV68evXv3JiYmhkcffZRatWol+c8rKfQZuYiI3NbwVhXxcnO5aZuXmwvDW1VM8ZgVKlRgx44d\nVK9enddee41Ro0YB4Orqev3z7L8/U/33x7OMMVhrr3/WHRQUxOHDh3nmmWcAyJUrV5Lz3Hjs2rVr\nWbVqFRs3bmTnzp3Url37ehY3N7frORwOBx4eHte/v/b5urWWTz755HqmI0eO8OCDD9KsWTPWrVtH\n0aJF6dWrFzNmzEhyvqRQIRcRkdt6tHZR3n2sOkX9vDAkzMTffaw6j9YumuIxT548ibe3N926dWP4\n8OHs2LEDSGhVvX37dgC+/fbbm8754YcfiIyM5OzZs6xdu5Z69erRqlUrpkyZwpUrVwA4ceIEoaGh\nKc4FCTPovHnz4u3tzf79+9m0aVOyzm/VqhUTJkwgJiYGgIMHDxIeHs7Ro0cpWLAgffv2pU+fPtd/\n5rSS6kvrxpghwCAgDvjRWvtS4vYawOdAHiAeqGetjTTG1AWmAV7AEuB5a61NbY6kWvDrCcYuP8DJ\nC1cp4ufF8FYVU/UfpYhIdvZo7aJp+m/k7t27GT58OA6HAzc3NyZMmADA66+/zjPPPMPIkSNp3rz5\nTefUqFGDFi1aEBYWxsiRIylSpAhFihRh3759NGqUcInfx8eHWbNm4eLi8ve3TLLWrVszceJEKleu\nTMWKFZN8ef6aPn36EBwcTJ06dbDWEhAQwIIFC1i7di1jx47Fzc0NHx+fNJ+Rm9TUUGNMC2AE8LC1\nNsoYU8BaG2qMcQV2AN2ttTuNMfmBC9baOGPMFuA5YDMJhXyctXbp3d4rMDDQbtu2LcVZ4a9HKa7G\nxF3f5uXmkurfMEVEsop9+/ZRuXJlZ8dIsjfeeAMfHx+GDRvm7Cjp6nZ/L8aY7dbawLudm9pL6wOA\n0dbaKABr7bXrGg8Cu6y1OxO3n00s4oWBPNbaTYmz8BnAo6nMkGRjlx+4qYhDwqMUY5cfyKgIIiIi\naSq1l9YrAE2NMf8FIoFh1tqtidutMWY5EADMtdaOAYoCITecH5K47baMMf2AfgAlSpRIZVQ4eZu7\nL++0XUREnOuNN95wdoRM766F3BizCih0m10jEs/PBzQE6gHzjTFlErffk7gtAlhtjNkO3P2BwRtY\na78AvoCES+vJOfd2ivh5XX+UoqQ5zVFb6Pp2ERGRrOiul9attfdba6vd5usHEmbU39kEW0i4qc0/\ncfs6a22YtTaChM/C6wAngGI3DF8scVuGuPYoRWvHFsa5fYobCY8M9GhUMqMiiIiIpKnUfka+AGgB\nYIypALgDYcByoLoxxjvxxrd7gb3W2lPAJWNMQ5PwQF4P4IdUZkiya49SHMzdkCejXyN/nlzkcnfh\n2x0hRERn3XV2RUQk50rtZ+RTgCnGmD1ANNAz8Sa288aYD4CtgAWWWGt/TDxnIH89frY08SvD3PQo\nRWw0R394k0e21eDV73bzYeda6tMrIiJZSqpm5NbaaGttt8RL7XWstWtu2DfLWls1cd9LN2zflrit\nrLV2cEY+Q36LM/spufcLxlYLYUHQSWZtvn0TABERyfzeeOMN3n///SQfHxwcTLVq1dIxUcbI2Wut\nF64Bg7fS0rckzaO38taivdQo6kvN4n7OTiYiIrcRGxuLq2vOLl1/pyVa85bC4TB8fL8Pnby3M3D2\nDs6HRzs7lYhIxpj6MPw6O+H7uJiE1zvnJbyOjkh4vSdxydTIiwmv9y5MeB1+NuH1gcRPSC//maS3\nDA4OplKlSnTt2pXKlSvTsWNHIiIiGDVqFPXq1aNatWr069ePaxdsmzdvzgsvvEBgYCAff/wxixYt\nokGDBtSuXZv777+fP//863137txJo0aNKF++PF9++SWQsAb68OHDqVatGtWrV2fevHm3zdS0aVPq\n1KlDnTp12LBhQzL/IJ1HhTyR74Z3ed19FpcuX2bo/CDi4513xV9EJLs7cOAAAwcOZN++feTJk4fP\nPvuMwYMHs3XrVvbs2cPVq1dZvHjx9eOjo6PZtm0b//rXv7jnnnvYtGkTv/76K126dGHMmDHXj9u1\naxdr1qxh48aNjBo1ipMnT/Ldd98RFBTEzp07WbVqFcOHD+fUqVM35SlQoAArV65kx44dzJs3j+ee\ney7D/ixSS9cnrmk/Dveoy7y03zJywR7G/3SYIS3LOzuViEj6evrHv753cbv5tbv3za89fW9+nSv/\nza9zF0zy2xYvXpwmTZoA0K1bN8aNG0fp0qUZM2YMERERnDt3jqpVq9KuXTsAOnfufP3ckJAQOnfu\nzKlTp4iOjr6pDek/tTx98skncXFxoWDBgtx7771s3bqVGjVqXD8vJiaGwYMHExQUhIuLCwcPHkzy\nz+JsmpFf45UX/ErQrUEJ3iizj89W7Wb9oTBnpxIRyZZu15p04MCBfPPNN+zevZu+ffve1M70xnaj\nQ4YMYfDgwezevZvPP//8puNuN25SfPjhhxQsWJCdO3eybds2oqOzzkesKuR/Y0L30vPk27yQZx3P\nzf2VUxe1fKuISFo7duwYGzduBOCrr77innvuAcDf358rV67wzTff/OO5Fy9epGjRhMeIp0+fftO+\n27U8bdq0KfPmzSMuLo4zZ86wbt066tevf8uYhQsXxuFwMHPmTOLibu7LkZmpkP9dwaqYnoto+fSb\nRMXEMWj2DqJj452dSkQkW6lYsSLjx4+ncuXKnD9/ngEDBtC3b1+qVatGq1atqFev3j+e+8Ybb9Cp\nUyfq1q2Lv7//TfuutTxt2LDh9ZanHTp0oEaNGtSsWZP77ruPMWPGUKjQzSuPDxw4kOnTp1OzZk32\n799/0xWAzC5VbUwzUlq0MU2uZdv28+W3S6jRuBWvt6uaoe8tIpIeMkMb0+DgYNq2bcuePXucmiMz\nSU0bU93sdgetg8fQzHsV9X4pQd2SeWlbo4izI4mIiNxEl9bv5IG3cOs6h4olCvPvb3ZxOPSKsxOJ\niGR5pUqV0mw8DamQ34lvUdzKNGF81zrUcg1m8Kwtaq4iIiKZigp5EhSOOcFMO4L7z83h1e92k1Xu\nKxARkexPhTwp/MvhaP8xPk0HqrmKiIhkKirkSVW7G/0eqMV9FfIxfdEadh6/4OxEIiIiObOQR8RE\nMGrjKM5FnkvWeQ6HYULeucx3e51/z1qn5ioiIing4+MDwMmTJ+nYsaOT02R9ObKQHzh/gIW/L6Tn\n0p6cunLq7ifcwOOegYTf8yp/XHFTcxURkVQoUqTIHVdwk6TJkc+R1y5Qm88f+Jwhq4fQf1V/vmv/\nHS4Ol6SdHFCR4vdXZKTPUcYtWM+na/x47v4K6RtYRCSdPL3s6Vu2tSrVii6VunA19ioDVw28Zf8j\n5R7h0XKPcj7yPC+uffGmfVNbT03ye9+4MMy0adNYuHAhERER/P7773To0OF6V7MVK1bw+uuvExUV\nRdmyZZk6der1Wb3k0Bk5QN2CdZnaeiqvNng16UX8Bt0qwk+5XiZ87Qf879CZdEgoIpKzBAUFMW/e\nPHbv3s28efM4fvw4YWFhvP3226xatYodO3YQGBjIBx984OyomUqOnJFfUzFfxevfzz8wn+K5i9Oo\nSKMknWv8SuDRoA8HdlXg67lBLB5yD0X8vNIrqohIurjTDNrL1euO+/N65k3WDPxuWrZsia+vLwBV\nqlTh6NGjXLhwgb17915veRodHU2jRkn7dzqnyLEz8hvFxMUw78A8Bq0exMqjK5N2kjG4Pfg6I3u2\nIyomjmGz1qu5iohIKnh4eFz/3sXFhdjYWKy1PPDAAwQFBREUFMTevXuZPHmyE1NmPirkgJuLG1Na\nTaFq/qoM+3kY3xxM+s0XZQN8+K7aRt4NHciHCzenY0oRkZynYcOG/PLLLxw+fBiA8PBwDh486ORU\nmYsKeSJfD18+f+BzGhdpzJsb32Ty7qT/xlexUVtOFWrOl1tCWbTzZDqmFBHJWQICApg2bRpPPvkk\nNWrUoFGjRuzfv9/ZsTIVtTH9m5j4GF5b/xpV81elR9UeST4vOjaeLl9s5I/T5/hmcAvKFdAdlSKS\n+WSGNqZyq9S0MdWM/G/cHG6Mbjr6ehE/dP4QsfF3b5Ti7upgQofifO14lYVT3lVzFRERyRAq5Ldh\njAEg7GoY3Zd258W1LxIVF3XX8woWKIJfyRoEXfJRcxUREckQKuR34O/lz5DaQ/jp+E/0X9mfK9F3\n6Ufu4krA07MJbNkxobnKpuAMySkiIjmXCvlddK3clXebvktQaBC9l/fm7NWzdz1ncIty/LvYbqot\n7cTuP3Tzm4iIpB8V8iRoW6YtH9/3MUcuHmHG3hl3Pd7hMHS7twYxrt4Mn7dNzVVERCTdqJAnUbNi\nzZj10CwG1x4McNfPv3NXfwiPXj8kNFeZ96uaq4iISLpQIU+Givkq4uZwI+xqGN2WdmPXmV13PL5m\niby82aYUnY+MYMnXkzIopYhI2rm4aBGH7mvJvspVOHRfSy4uWpThGfr06cPevXvveMyCBQvuekx2\npUKeApGxkZy7eo4+K/qw4eSGOx7bpUEpKvlEsHnXb2quIiJZysVFizg18j/EnjwJ1hJ78iSnRv4n\nw4v5pEmTqFKlyh2PUSGXZCmWuxgz2sygWO5iDFo9iBXBK/7xWOPmRcHn17DZvwPPzw3i5IWrGZhU\nRCTlQj/8CBsZedM2GxlJ6IcfpWrc4OBgKlWqRNeuXalcuTIdO3YkIiKC1atXU7t2bapXr07v3r2J\nikp47Ld58+ZcWxDMx8eHESNGULNmTRo2bMiff/7Jhg0bWLhwIcOHD6dWrVr8/vvvqcqX1aiQp1CA\ndwBTW02lun91hv087I7F3NvTkwnd6lIhZj+HJ3QmOuruz6SLiDhb7KlTydqeHAcOHGDgwIHs27eP\nPHny8MEHH9CrV6/rbUxjY2OZMGHCLeeFh4fTsGFDdu7cSbNmzfjyyy9p3Lgx7du3Z+zYsQQFBVG2\nbNlU58tKVMhT4dr67J0rdiaw0J1X0Ssb4MPL9d0oeXUf4xetz6CEIiIp51q4cLK2J0fx4sWvtybt\n1q0bq1evpnTp0lSoUAGAnj17sm7dulvOc3d3p23btgDUrVuX4ODgVGfJ6lTIU8nL1YsRDUeQzzMf\nMXExfHPwG+Lt7duZ1mrbn1l15/Hxtkg1VxGRTK/A0Bcwnp43bTOenhQY+kKqx762guY1fn5+STrP\nzc3t+rnXWp3mdCrkaWhZ8DLe3PgmI38Z+Y/rsw9/uBZ1i+fh2HcjCQlalcEJRUSSzrddOwq/NQrX\nIkXAGFyLFKHwW6Pwbdcu1WMfO3aMjRs3AvDVV18RGBhIcHDw9XalM2fO5N57703yeLlz5+by5cup\nzpUVqZCnobZl2jKo1iAW/r6QoWuHEhkbecsx7q4Oxj9RiYfMBjYumUl4lH6bFJHMy7ddO8qvWU3l\nfXspv2Z1mhRxgIoVKzJ+/HgqV67M+fPnGTp0KFOnTqVTp05Ur14dh8NB//79kzxely5dGDt2LLVr\n185xN7upjWk6mLt/Lu9sfoc6BevwyX2fkNs99y3HbNpzmCdn76d9zaJ81LnWLZeZRETSQ2ZoYxoc\nHEzbtm3Zs2ePU3NkJmpjmsl0qdSF0U1H88eFPzgVfvu7OxtWK8eL91dkc9Buds8cDvG3/1xdRETk\nTlTI08lDZR5i6eNLqZA34Q7Mi1EXbzlmUItyDCpygDK/z2Tfb9szOqKIiFOUKlVKs/E0pEKejnK5\n5QJg1t5ZPPbDYxw+f/im/Q6Hod0z/6GH5zj6/HhZzVVEJENklY9Uc4rU/n2okGeABoUbYLH0XNaT\noNCgm/b55fLgje6tOHM5illTPyH+lH5LFZH04+npydmzZ1XMMwlrLWfPnsXzb4/5JUeqb3YzxgwB\nBgFxwI/W2pcSt9cAPgfyAPFAPWttpDFmLVAYuLZW6YPW2tC7vU9WutntdkIuh9BvZT/CrobxYfMP\naVK0yU3756zfR/OVD3EpIJCKg792UkoRye5iYmIICQkhMvLWp2rEOTw9PSlWrBhubm43bU/qzW6p\nKuTGmBbACOBha22UMaaAtTbUGOMK7AC6W2t3GmPyAxestXGJhXyYtTZZVTmrF3KAsKth9F/ZnyMX\nj7DksSUUzFXw+j5rLe/OXMz0fXFM6t2EpuUDnJhUREScLaPuWh8AjLbWRgHcMLN+ENhlrd2ZuP2s\ntTYule+V5fl7+TOl9RTea/beTUUcElY5eqFLG0oWyMtLczZzfsN0J6UUEZGsJLWFvALQ1Biz2Rjz\nszGm3g3brTFmuTFmhzHmpb+dN9UYE2SMGWnu8AC1MaafMWabMWbbmTPZowVoHvc83F/yfgB+OfEL\nX+768vpnVd7urkzoVpe2sSvJu+I5YkKC7jSUiIgIrnc7wBizCih0m10jEs/PBzQE6gHzjTFlErff\nk7gtAlideIlgNdDVWnvCGJMb+BboDsy43Xtba78AvoCES+vJ/NkyvTXH1jD/4HzORZ5jeL3hOIyD\nsgE+1Hp8OI/PLUn1HW68UczZKUVEJDO7ayG31t7/T/uMMQOA72zClHKLMSYe8AdCgHXW2rDE45YA\ndYDV1toTieNeNsZ8BdTnHwp5djei4QjcXdyZtW8WF6IuMKrJKNwcbjxcsxjbjrVi6i/BtPA9zb01\nykHeks6OKyIimVBqL60vAFoAGGMqAO5AGLAcqG6M8U688e1eYK8xxtUY4594vBvQFsixz1s5jIOX\n6r3E4FqDWfzHYl746QWi4xKeJX+lTWXqF89FxTV9CP8+9Z2GREQke7rrjPwupgBTjDF7gGigZ+Ls\n/Lwx5gNgK2CBJdbaH40xuYDliUXcBVgFfJnKDFmaMYZnaz6Ln4cf+87tw82R8PiBu6uDj7s14JWP\nXyT8XHGmRsWSyyO1f10iIpLdqGlKJmKtxRjD8UvH8XLzwt/Ln/WHwug+ZTPtaxTmo8ZRmJKNnB1T\nREQygJqmZEHGGOJtPM/99Bw9lvYg5HII95T358X7K+C55yvM1NZwfIuzY4qISCaiQp7JOIyDNxq/\nwcWoi/RY2oOD5w8yqEU5zpV9lFdi+xFkyzs7ooiIZCIq5JlQzYCaTG89HYOh17Je7Arbydgu9Vnn\n04ZBX/3KhdATEB7m7JgiIpIJqJBnUuXylmPGQzPI55mPT4M+xdfLjQnd6nDucgQRX7TCftfP2RFF\nRCQT0G3QmVhRn6JMaz0NN4cbxhiqFc3DiHbVeXPh4zTPXZcnnR1QREScTjPyTM7fyx9fD1+i4qJ4\nduWzuPhuwLtmB17d7MK6g2fg4glnRxQRESdSIc8irLV4unryzpZ3KF3+F8oVyMUPcyZix9WCY5ud\nHU9ERJxEhTyL8HT15MPmH9K+bHu+3D2RurXX8b/4Kixwb0d0QFVnxxMRESdRIc9CXB2uvNXkLXpU\n6cGPR7+mbsOdDD3/OO+sPAqx0RAd7uyIIiKSwXSzWxbjMA6GBQ6jgHcBGhdpTOHYGKb+8juDjg0l\nwD8AOs+Cf+4MKyIi2Yxm5FmQMYaeVXtSPm95Xm5TkdJldzEurAKni7dWERcRyWFUyLO4kCtHOe/x\nDT8UPciTW3wJj4qFmEhnxxIRkQyiQp7FlfErw/j7x+PueZ5Qn//jg9njsR/XhFO7nB1NREQygAp5\nNtC4SGOmtp6Cl2c08+0c/perNHjldXYsERHJACrk2UT1gOrMeXgm7i6e9D3fkqDLucFaiI9zdjQR\nEUlHKuTZSPl85Vjy+CL8HXUZNGs7J759FhYPTSjoIiKSLamQZzMFc+dmQrc6hMXtpe2VzSy2erZc\nRCQ7UyHPhmoU8+OlFvcTGVGKVy5sYfb+rzQrFxHJplTIs6nejSvxQL4RxF6uwugto/lsaiPs2T+c\nHUtERNKYCnk2ZYxh9GN1KBbzLD6XqzHBJZz/Hf/J2bFERCSNaYnWbMzb3ZWJ3erzyKfRFHE9SIOq\nvZwdSURE0phm5Nlc2QAfxnSsycHgsry7ZB9//O89/v11WyJiIpwdTURE0oAKeQ7wUPXC9G5Smmkb\ng1l3eAPLIo7Sb2VfLkZddHY0ERFJJRXyHOKVhypRt2Q+PjzSl2G13mbv2X30WtaL0IhQZ0cTEZFU\nUCHPIdxcHIx/qg6ubh5MX+3HBw3f5eSFI/T4sSunrpxydjwREUkhFfIcpJCvJ+OerM0fZ66wef1F\nJoVdpLKnP3k9tS67iEhWpUKewzQp58+LD1RgwoFc7Kn3Ax+2m4OnqyeXoi+x64w6pomIZDUq5DnQ\nwObluK9SAf6zPIRfj52HP35mzNJ+9F7em5+P/+zseCIikgwq5DmQw2H44ImaFMzjyaDZO4je9AUv\nnjxKOd+yPP/T8yz6fZGzI4qISBKpkOdQft7ufNa1DmFXohkc0RffXsuY3HoKgQUDeXX9q8zaO8vZ\nEUVEJAlUyHOwGsX8eL19FVYcDueTX/4kl8OD8V4VaVm0GTP2zuBK9BVnRxQRkbvQEq053FP1S7A9\n+Dwfrz5Es1wh1PnpXd5v+yFhjUbi4+5DXHwcAC4OFycnFRGR29GMPIczxvDfDtWpUCA3z6yMJbT7\nWlzr9qRQrkJYa/nv5v/y8v9eJiYuxtlRRUTkNlTIBS93FyZ0q0NMnKXf0stEx8bDuSOYg8spnrs4\ny4KXMXjNYK3PLiKSCamQCwBlAnwY07EGQccv8M6SfbD8VVg8lKcrPsmoxqPYdGoTfVf05ULkBWdH\nFRGRG6iQy3XXm6tsCGZZ2Veh12Jw86RD+Q58cO8H7Du3j/6r+hNv450dVUREEulmN7nJKw9VYmfI\nBV5cfIJyg5tQDmD/ElqWa8nE+ycSGReJw+j3PxGRzEL/IstNrjVX8XJzof+sHVw9FgRzn4Stk6hf\nuD7NijUDYOHvC9l3dp+T04qIiAq53OLG5ir//sVin5oP9Z+9vj8yNpLPgj6j9/LebDu9zYlJRURE\nhVxu61pzlYU7TzLzbEVwcYWoy3DyVzxdPZnWehoB3gH0X9Wfn4795Oy4IiI5lgq5/KNrzVXeWrw3\nobnKwudgVkeIukKhXIWY3no65f3KM3TtUH44/IOz44qI5EipLuTGmCHGmP3GmN+MMWMSt3U1xgTd\n8BVvjKmVuK+uMWa3MeawMWacMcakNoOkj783V7nQ+FXoOBk8fADI65mXSa0mEVgokNCIUCenFRHJ\nmYy1NuUnG9MCGAE8bK2NMsYUsNaG/u2Y6sACa23ZxNdbgOeAzcASYJy1dund3iswMNBu26bPY51h\nV8gFOk7YSMOy+Znaqx4uDgN/7oWAiuBwITY+FhfjgjGGE1dOUCRXEfT7mYhI6hhjtltrA+92XGpn\n5AOA0dbaKIC/F/FETwJzE0MVBvJYazfZhN8gZgCPpjKDpLNrzVXWHTzDJ2sOwZmD8MW9sOETAFwd\nrhhjOB1+mk6LOvHWpreur9EuIiLpK7WFvALQ1Biz2RjzszGm3m2O6QzMSfy+KBByw76QxG23ZYzp\nZ4zZZoxX8TosAAAgAElEQVTZdubMmVRGldR4qn4JHqtdlI9XH+Lnc37w4NtQp8dNxxT0Lkjnip35\n+uDXDF83nOi4aCelFRHJOe66IIwxZhVQ6Da7RiSenw9oCNQD5htjyiTOtjHGNAAirLV7UhLOWvsF\n8AUkXFpPyRiSNq41V/nt5CVemBfE4ud6UNTbC+Lj4fJJ8C2GMYbn6zyPn4cf7297n8vRl/m4xcd4\nu3k7O76ISLZ11xm5tfZ+a22123z9QMKM+jubYAsQD/jfcHoX/pqNA5wAit3wuljiNskCbmyuMmj2\njoTmKkuHw6QH4Opfa7D3rNqTt5q8xdbTW5m8Z7ITE4uIZH+pXaJ1AdAC+MkYUwFwB8IAjDEO4Amg\n6bWDrbWnjDGXjDENSbjZrQfwSSozSAa61lxl4Owd/PfHvbxZvyf4VwRP35uOe7TcoxTzKUaNgBpO\nSioikjOk9jPyKUAZY8weEm5o62n/ug2+GXDcWvvH384ZCEwCDgO/A3e9Y10yl4eqF+aZe0ozfeNR\nFob6Q4N+YAxEnIMbnoIILBSIu4s7FyIvMGDVAIIvBjsvtIhINpWqQm6tjbbWdku81F7HWrvmhn1r\nrbUNb3POtsTjy1prB99Q+CULeblNJQJL5uXlb3dx6M/LcOEYjG8Amz+/5dgzV8+w9+xeei7ryd6z\ne52QVkQk+9LKbpIibi4OPn2qDt7uLgyYvYNwz0JQvSOUaX7LseXzlmd66+l4unjSe3lvtp7emuF5\nRUSyKxVySbFCvp6M65LQXOXl73/DtnoHClRK2BkdcdOxpXxLMaPNDAp5F6L/yv5sOLHBCYlFRLIf\nFXJJlcbl/PnXgxVZtPMkMzcdTdi4/iP4ojlEXrzp2IK5CjKt9TRalmhJ+bzlMz6siEg2pEIuqTbg\n3rK0vLG5StG6ULIRuHrecqyfpx9j7h1DgHcAsfGxrD662gmJRUSyDxVySTWHw/B/NzRXOVegAbT7\nGFw9IC72H8/77tB3vLD2BT7a/hG651FEJGVUyCVN+Hm781nXOoRdieaFeUHExVsID4NJLWHX/Nue\n83j5x+lUoROT90zmzY1van12EZEUUCGXNHNLcxWPPJC7MHjkvu3xLg4XRjYcSb8a/fj20LcM+3mY\n1mcXEUmm1K7sJnKTp+qXYHvweT5efYjaJfJy71Nz/9ppbcLCMTcwxjCk9hD8PPwYHzSe4EvBVMhb\nIYNTi4hkXZqRS5q61lylQoHcvDD3V05cuJqw47cFML0dxETe9rzuVbqzuMPi60U8Mvb2x4mIyM1U\nyCXN3ba5inFAfCzERPzjef5eCf125u6fS+fFnTkdfjqjIouIZFkq5JIurjVXCTp+gf/+uBeqtIde\nS8A7313PLedXjtCIULov7c4fF/++VL+IiNxIhVzSzU3NVXaeBIcjYcW37/rBoVX/eF5goUCmtJpC\ndFw0vZb24rew3zIwtYhI1qJCLunqluYqNh5C90LYwTueVzl/ZWa0mYG3mzd9VvTh7NWzGZRYRCRr\nMVllIY7AwEC7bds2Z8eQFDh9MZK2n/wPP293fhjUhFwu8eDqnqRzQyNC2XRqE+3Ltk/nlCIimYsx\nZru1NvBux2lGLunupuYq3+3Gurgl7Di1ExYOgTssBFPAu8D1Ir719Fa+P/R9RkQWEckyVMglQ9zY\nXGXGxsTmKid2wOE1cOlkksaYu38u/9nwH6bumZqOSUVEshYtCCMZZsC9Zdlx9Dxv/7iX6sV8qRP4\ndEIP839Y+e3vRjcdjcM4+GD7B5yPOs/QOkMxf1tgRkQkp9GMXDKMw2H44IlaFMzjyeDZOzgXHp1Q\nxK2F9R9CyPY7nu/m4sbopqPpXLEzU/dM5Y2Nb2h9dhHJ8VTIJUP5ersxoWtdwq5E8/zcXxOaq0Rd\ngm1TYffXdz3fxeHCiAYj6F+zP1FxUZqRi0iOp7vWxSm+2nyMV7/fzfMtyzP0gQpwJRRyBdyyFvud\nxNt4HMbB6fDT+Lj54OPuk46JRUQylu5al0ztyfrFeaxOUcatOcTaA6HgUyChiF8Jhf99kHC5/S4c\nxkFsfCz9V/bnmRXPcC7yXAYkFxHJXFTIxSmMMfz30epULJibofOC/mqusms+/DwGwg4laRxXhysv\nBr7I7xd+p+fSnpy6ciodU4uIZD4q5OI0Xu4ufNb1b81VGg2CAb9AQNJbmTYr1ozPH/ics1fPJqzP\nfkHrs4tIzqFCLk5VJsCH9zvd0FzFGMhfNmHngWXw594kjVO3YF2mtp5KbHwsY7aOScfEIiKZi54j\nF6drXa0wfe4pzaT1R6hbKh/taxaB6HD48UUoUhu6zE7SOBXzVWRmm5nkcs8FgLVWd7WLSLanGblk\nCv/+e3MV91zQ60foMDFZ4xTPU5x8nvmIiY/hhZ9eYOXRlemUWEQkc1Ahl0zBzcXBp0/VwdvdhQGz\ndxAeFQv5SicsGBMXk7Am+6ldSR4vMjaSc5HnGPbzML45+E06JhcRcS5dWpdM41pzlW6TN9Nt0mZC\nL0dy8kIkNXwjmONYiXfhmlC4RpLGyu2em88f+Jx//fwv3tz4JhejLvJM9WfS+ScQEcl4mpFLptK4\nnD9tqhXi1+MXOHEhEgvsvOhN0yvvssC1TcJBSVyW1dvNm3H3jeOh0g/x0Y6P+Hzn5+kXXETESVTI\nJdP59fiFW7adjXFn7PIDcOYAjK8PJ+68Lvs1bg433m36Lk9XfZpmxZqldVQREadTIZdM59SFyNtu\nP3nhKrh5g7c/eORJ8ngO4+DFwBepnL8yAIt+X0RUXFSaZBURcTYVcsl0ivh5/fN2v+LQexn4l0/Y\nePl0ssbef24/I9aPoP/K/lyJvpLaqCIiTqdCLpnO8FYV8XJzuWV71SK5sdb+1Vhlxwz4JBBC9yV5\n7Er5KjG66WiCQoPovbw3Z6+eTavYIiJOoUIumc6jtYvy7mPVKernhQGK+HnSsHQ+VuwN5dXvdxMb\nF59wYLkHILAX5C+frPEfKvMQ4+4bx5GLR+i5rCcnr5xM859BRCSjqI2pZAnWWt5fcYDxP/3Og1UK\nMu7J2njeOGuPvARhB6HYXTv+Xfdr6K8MWTOE1xu9zgMlH0iH1CIiKac2ppKtGGMY3qoSr7erwoq9\nf9JjyhYuXo3564BlL8PMx+DqrXe8/5PaBWqz5LEl14u4PjMXkaxIhVyylKeblObjLrX49dh5On++\nkdBLiXe43/8mdJwCXn7JGi+Pe8Ld75tPbabVt63YcHJDWkcWEUlXKuSS5TxSqyiTe9bj2LkIHp+4\ngSNh4eATAOXvTzjgyP/g9zXJGrOsX1kK5SrEoNWDWB68PB1Si4ikDxVyyZKaVQjgq74NCY+Ko+OE\nDewOuZiww1pYPQpWvQnx8Ukez9/Ln6mtp1LDvwbDfx7O/APz0ym5iEjaUiGXLKtWcT++7t8ITzcX\nunyxkV8OhyU8mvbkXOj6NTiS9593Hvc8THxgIk2LNeWtTW+x+dTmdEouIpJ2VMglSysb4MO3AxpT\nLK83T0/dyuJdJyFXfvApkDA7X/YKHEz6pXIvVy8+avERoxqPon6h+umYXEQkbaS6kBtjhhhj9htj\nfjPGjEnc1tUYE3TDV7wxplbivrXGmAM37CuQ2gySsxXy9WT+s42oWdyXIXN+ZcbG4IQd0eFwdAMc\n25Ss8dwcbnQo3wFjDMcuHeO9Le8REx9z9xNFRJwgVW1MjTEtgEeAmtbaqGtF2Vo7G5ideEx1YIG1\nNuiGU7taa/VQuKQZX283Zj7TgMFf/cp/fviNsMtRDH2gAubppeCWuORrXAy4uCVr3PUn1jNr3yxC\nLocw9t6xeLp6pkN6EZGUS+2MfAAw2lobBWCtDb3NMU8Cc1P5PiJ35enmwsRudXgisBjj1hzm1e/3\nEOfqlfC5+aWTMKEx7FucrDGfqvwUIxqM4OeQn+m/qj+Xoy+nU3oRkZRJbSGvADQ1xmw2xvxsjKl3\nm2M6A3P+tm1q4mX1kcZcWzhbJPVcXRy893gNBjYvy5wtxxg0eweRMXHg7gN+JSB3oWSP2aVSF95r\n9h47Q3fSe3lvzkWeS4fkIiIpc9dL68aYVcDt/vUbkXh+PqAhUA+Yb4wpYxPXfTXGNAAirLV7bjiv\nq7X2hDEmN/At0B2Y8Q/v3Q/oB1CiRIkk/1CSsxljeKl1Jfx9PBi1eC89p2zhy56B5On27V8HXQwB\n32JJHrNN6Tbkds/NtN+m4emiy+siknmkaq11Y8wy4D1r7U+Jr38HGlprzyS+/hA4Y6195x/O7wUE\nWmsH3+29tNa6pMQPQSf41/ydlC+Ym+lP16NAHk84tArmPglPzYOy9yVrPGstxhjCY8L5M+JPyviW\nSafkIpLTZdRa6wuAFolvWAFwB8ISXzuAJ7jh83FjjKsxxj/xezegLbAHkXTySK2iTO5Vj6Nnw3l8\n4gaCw8KhRENo0B+KN0z2eNc+CXpr01t0X9KdoNCgu5whIpK+UlvIpwBljDF7SCjYPe1fU/xmwHFr\n7R83HO8BLDfG7AKCgBPAl6nMIHJH9yauAnclMpaOEzewJywOHnwL3L0hNirZy7kCDK41GD8PP/qt\n7McvJ35Jh9QiIkmjNqaSYxwOvULPxK5pX3SvS+Ny/rB2NPz8HgzeBvnLJmu8sKthDFg1gMMXDvPu\nPe/SunTrdEouIjmR2piK/E25AgmrwBXx86TX1K0s2X0KmjwPnWcnu4hDwvrsU1pNoWZATd7f9j5X\nY6+mQ2oRkTvTjFxynIsRMTwzfSvbj51n1CPV6N6wZMKO03vg1E6o3TVZ40XGRhIaEUqJPCW49v+T\nnqoUkdTSjFzkH1xbBe6+igUYuWAPH648mFCAN3wCP72TsLRrMni6el4v4mO2jmHM1jHE26R3XhMR\nSY1ULdEqklV5ubvwefe6vPLdbj5efYiwK1GMavsxLhFnwD1Xisc1xjBz70wuRF1gVJNRuDmStySs\niEhyqZBLjuXq4mBMxxrk9/Fg4s+/cy48mg8718ITYN374JEHGvRL8njGGIYHDievR17G/TqOy9GX\nef/e97U+u4ikK11alxzNGMPLbSrx2sOVWbrnNL2mbuFSRCSc2AEntie0Qk3meH1r9GVkw5GsC1nH\nkDVDyCr3oYhI1qQZuQjQp2kZ/H08GPb1Trp8uZVpPSdSIE+uhIYrsdHg6p6s8Z6o+AR5PPLgZtx0\n45uIpCvNyEUSPVq7KJN6BnIkLJyOX2zn6PlIiLwEUx6EDZ8me7zWpVrTsmRLAFYdXUXI5ZC0jiwi\nokIucqPmFQvwVd8GXI6M4fEJG9hzJgYCKkH+cikeMyImgrc3vU2PpT04eP5gGqYVEVEhF7lF7RJ5\n+bp/Y9xdHHSZtI0NNd6Giomrtp0PTvZ43m7eTHpwEgZDr2W9tD67iKQpFXKR2yhXwIdvBzamsK8n\nvaYkrgJ3ejeMbwDbpyV/vLzlmPHQDPJ55qPvir6sP7E+7UOLSI6kQi7yDwr7evF1/0ZUL+bLoK92\nMPtIroQlXSu3T9F4RX2KMr31dEr7lmbXmV1pnFZEcioVcpE78PN2Z9YzDWhRsQAjftjHR3EdsV55\nIT4e9i1O9uNp+b3yM73NdAbUHADAuchz6RFbRHIQFXKRu7i2CtzjdYrx0apDjPxhD/G75sG8rnDk\n5+SP5+qFMYYTV07w6IJHmRA0Qc+ai0iK6TlykSRwc3Hwfqca+Od25/Of/+D8lXJ81GkWbqXvTfGY\nBb0L0rRYUz7b+Rnno87zcv2XcRj9bi0iyaNCLpJExhheaVOZAB8P3v5xH+ciAviibCy5r56EPd/C\nPUMTFpBJIleHK281eQs/Dz9m7J3BxaiLvH3P21qfXUSSRb/+iyRTn6Zl+LBzTbYGn6PLF5sI3zIT\nfvkILp1I9lgO42BY4DCer/M8S44sYdqeaWkfWESyNc3IRVKgQ+1i+Hm7M3DWDh6KbMxXnTtR1LdY\nisYyxtCneh9K+5amSZEmaZxURLI7zchFUqhFxQLM7tuAi5GxPDL7OL+dvAi/zoKlLyf7bnaAliVa\n4unqyaXoS/x73b85E3EmHVKLSHajQi6SCnVK5OWb/o1wdzF0/nwTJw8FQdgBiItO8ZhHLh7hp+M/\n0WNpD45fOp6GaUUkO1IhF0mlcgVy882AhFXgmu+8j2U1PgJXD4iNSnjePJlqBtRk8oOTuRxzmR7L\nenDg3IF0SC0i2YUKuUgaKOKXsApctaK+DJi7h682HoY5XWDJv1I0XvWA6sxoPQOHcfD0sqe1EpyI\n/CMVcpE04uftzuw+DWleIYBXf9jP5ujS2CK1UzxeGb8yzGwzk6r+VSngXSANk4pIdqJCLpKGvNxd\n+KJHII/VKUbnQy35z7E6xMVbOPcHxMcle7wiPkX48sEvKZSrEHHxcWw9vTUdUotIVqZCLpLG3Fwc\n/F+nmjzbrAwzNx3ltZmrsV+0gNVvpmrcOfvn0Ht5b2bvm51GSUUkO9Bz5CLpwBjDKw9VJr+PO+8s\n2U+lgk/yePWe+KRizE4VO7H19FZGbxnNhagLDKw5EJOMleREJHvSjFwkHfVrVpb/61STUWea0eXr\nk5y5FAm/fQ9xsckey8PFg/9r/n88Wu5RJu6cyDub3yHeJv+ueBHJXlTIRdLZ43WLMalHIIdDrzDq\ns8nwdS8IStnlcVeHK6Maj6JX1V58f/h7jlw8krZhRSTLMVmlfWJgYKDdtm2bs2OIpNj2o+fpPW0r\n9zp20r93X6oU9UvVeMcvHad4nuIAxMbH4urQJ2Ui2YkxZru1NvBux2lGLpJB6pZMWAVuq2sdOn+x\nmW17D8PK1yE2ZavAXSvi3x78lqeXPc3FqItpGVdEsggVcpEMVL5gbr4d0JiCvp7MmTOVuI0T4M89\nqRrT18OX387+Rq9lvQiNCE2jpCKSVaiQi2SwIn5efP1sI/4o/BDNrr7PVyH+qRrv/pL3M+H+CZy8\ncpIeS3tw7NKxNEoqIlmBCrmIE+TN5c7sPg0oX6ESr36/m4VfT8XO656wPnsKNCjcgMmtJhMeE06P\npT10mV0kB1EhF3ESb3dXvuwRyGO1i7IpaDcnjx4kLjoixeNV86/G9DbT6VejH74evmmYVEQyMxVy\nESdyc3Hwfqea5Grch3vPvcpzC44QFR0NMZEpGq+MbxmeqvwUALtmf8qeZk3YV7kKh+5rycVFi9Iy\nuohkEnpeRcTJHA7DiIer4O/jwbtL9/NEyDs0KRCDa/dvwSVl/4teXLQIRk/AJSZhwZjYkyc5NfI/\nAPi2a5dm2UXE+TQjF8kknr23LO93qsniS2WZ82dxwq4mv8nKNaEffoRbzM2rvtnISEI//Ci1MUUk\nk1EhF8lEOtYtRptuL/LfKw/TccIGTgYfgOjwZI8Te+pUsraLSNalQi6SydxXqSCz+zQgPDwcO60t\nl+b0SfYYroUL3357oUKpjScimYwKuUgmVLdkPmYPuJfPXLrz9B/N2fzH2WSdX2DoCxhPz5u2GU9P\nCrw4lLNXzxITF5OWcUXEiVTIRTKpCgVzM3DIcC7kqUj3KVsIWj4dIi8l6Vzfdu0o/NYoXIsUAWNw\nLVKEwm+NwuuhVvRZ0YfBawYTEZPyR91EJPNQ0xSRTO58eDT/nryYT8/24XC5p6nS/f9SNd73h77n\njY1vUC1/Nca3HI+fZ+qat4hI+sjQpinGmCHGmP3GmN+MMWMSt7kZY6YbY3YbY/YZY1654fjWxpgD\nxpjDxpiX0yKDSHaVN5c7H/Vvz/uFxvLob035dM0hUvMLeIfyHfjg3g/Yd24fvZb14s/wP9MwrYhk\ntFQXcmNMC+ARoKa1tirwfuKuToCHtbY6UBd41hhTyhjjAowH2gBVgCeNMVVSm0MkO/N2d2V4v6d5\nuHYpPl2xm02f9SU+4nyKx2tZsiUT75/I6YjTvL3p7TRMKiIZLS0WhBkAjLbWRgFYa6+1X7JALmOM\nK+AFRAOXgPrAYWvtHwDGmLkk/CKwNw2yiGRbbi4O/q9TTWrF76X2/gV8Mas2T/fuj4erS4rGq1+4\nPlNaTSHAKyCNk4pIRkqLS+sVgKbGmM3GmJ+NMfUSt38DhAOngGPA+9bac0BR4PgN54ckbruFMaaf\nMWabMWbbmTNn0iCqSNbmcBh6dHmS+U0WM/qPUjwzbRtXIlN+B3qV/FUI8A4gNj6WEetHsPX01jRM\nKyIZIUmF3Bizyhiz5zZfj5Awq88HNASGA/ONMYaEmXccUAQoDfzLGFMmOeGstV9YawOttYEBAZo1\niAAYY+jxYEPGdqxBxJHNHHm/BefOpG6hl0vRl9gTtof+K/vz07Gf0iipiGSEJBVya+391tpqt/n6\ngYQZ9Xc2wRYgHvAHngKWWWtjEi+3/wIEAieA4jcMXyxxm4gkQ6fA4oy8vygeMed5duovHD+X8sfJ\n8nnmY3rr6VTMV5Gha4fyw+Ef0jCpiKSntLi0vgBoAWCMqQC4A2EkXE6/L3F7LhJm7PuBrUB5Y0xp\nY4w70AVYmAY5RHKc2i0e53KvtRyIyM3jn/3CgaMnUzyWn6cfkx6cRL1C9Xjtl9eYu39uGiYVkfSS\nFoV8ClDGGLMHmAv0tAnPxowHfIwxv5FQvKdaa3dZa2OBwcByYB8w31r7WxrkEMmR6pYO4JsBjekf\nPxfH1AfZfiA4xWN5u3kzvuV42pdtT5X8ephEJCvQgjAi2cSZ3StZu3AGI64+yadP1uHBqmmzrvq6\nkHU0KdIEF0fK7o4XkZTJ0AVhRMT5Aqo/QMsXJlO5sC8jZq1hwbrtqR5z55mdDFo9iOHrhhMdF50G\nKUUkramQi2Qj+XK5M6dPfebkHke5Vb0Zv+ZgqlaBqxlQk2GBw1h5dCWDV2t9dpHMSIVcJJvx9nCj\n5JMf8HPJ5xi74hBvLtpLfHzKi3nPqj15u8nbbDm9hT4r+nAh8kIaphWR1FIhF8mG3Eo3ZsDTvXnm\nntIc3fQ9r89aQXRsfIrHe6TcI3zY/EMOnT/EltNb0jCpiKSWbnYTycZs5CUi/686KyKr8E2pN5jQ\nrS4+HilfmTk0IpQC3gUAiI6Lxt3FPa2iisjf6GY3EcF45sGr9yLiHvqADb+f5akvN3H2SlSKx7tW\nxLee3srD3z/M3rNqkSDibCrkItld4Ro81qgyXzxVk8dDP2XIZwtStQocQIBXAAZD7+W9tT67iJOp\nkIvkEC0LR9LVayOVIrbSceIG9p++lOKxSvmWYkabGRTyLkT/lf1ZfWx1GiYVkeRQIRfJKfKXxfW5\n7XTu/zoAT0zcwNbgcykerlCuQkxrPY1K+Srx4toXCQoNSqukIpIMKuQiOUkufyoWys0PXQow3/Ea\nIyYtYOXeP1M8nJ+nH18++CVD6wylRkCNNAwqIkmlQi6SAxXK5UK5fC5U8Pfk2ZnbmL/1eIrH8nbz\nple1XjiMg5NXTvLlri9TtQiNiCRPyp9DEZGsq2BVXAdu5L2YeC7O2s7b327kzJUoBjYvizEmxcMu\n/H0h44PGc+LKCUY2HKn12UUygAq5SE7lcJDLw8HU6r8RfvItHl7xH8KuRDHy4So4HCkr5s/WeJaY\n+Bi+2PUFF6Mu8l6z9/SsuUg606V1kRzOtfQ95Kn1CK0a1WHqL8EMnR+U4lXgjDEMqT2El+q9xKpj\nqxi4eiDhMeFpnFhEbqQZuUhOF1AB0+4jXrOWQt4wd/VmngmPZmK3uuRK4Spw3at0x8/Djzn75+jz\ncpF0phm5iAAJs+m+Fz5mcZ732HE4JNWrwLUr246ZbWbi4+7D1dirnA4/nYZpReQaFXIR+UuLEXg9\n+hEfdb+H/acv02niRkLOp3wVuGs3u73+y+t0W9KNPy7+kVZJRSSRCrmI/CV/WajcjgeqFOT7tgb/\nK/t4fMIGDpy+nKphn672NDHxMfRc2pPfwn5Lo7AiAirkInI78XFU2fEm0wp+jY23dErlKnCV81dm\nRpsZ5HLLRe/lvdl8anMahhXJ2VTIReRWDhd4ah7e3b7i24FN8PfxoNukzaxKxSpwJfOUZEabGRTx\nKcJrv7xGdFx0GgYWybl017qI3J5fCQCKW8viqmt4fX8Jnp1lefex6jwRWDxFQxbwLsC01tP4M+JP\nPV8ukkY0IxeRO7t6Hu+DC/lv1RAal83PS9/sYsLa31P8WJmvhy8V8lYAYNyOcUzePTkt04rkOJqR\ni8ideeeDvmtw98rL5DjLsPlBvLdsP2FXohjxUOUUrwIXFx/H8cvHWRa8jItRFxlad2iqlocVyalU\nyEXk7rzzAeB+9U8+vjyUWjX6MWr9Ec5eiWJMx5q4uyb/4p6Lw4XRTUfj6+HL1N+mciHqAv9p9B9c\nHfpnSSQ59H+MiCSLMYanm5ThauEAxi4/wLmIGCZ0rZOiVeBcHC6MaDCCvJ55mbhzIhGxEYxtNlYz\nc5FkUCEXkaTLXQj6/oQxhkElobBbBMN+PM5TkzYztVc98uVK/g1sxhgG1RqEn4cfeT3yqoiLJJNu\ndhOR5LlWaA8u57F1DzGnjQv7T12i48QNqVoFrmvlrjxU5iEA1p9Yz7nIlD+3LpKTqJCLSMoUrQvV\nHqNBg3uY+UwDzlyOouOEjRz8M3WrwF2KvsRLP79Ez6U9OXXlVBqFFcm+VMhFJGVy+UP7ceDhQ/3i\nPizs4EW8tXSauJH/b+++46Oo1j+Of55NT6ihSBdClSJV6SqKUi6CgFdRBFQEbIBcu171d70qdmkq\nIoLAFRRsKAqIihK69A5SpfcWCKnn98cuikIkpC2bfN+v176ye+bM7LOH2Xk4s2fOLN6W8d50gdAC\nDLtuGAfjD9Jtajc2H9H87CJ/R4lcRDJv5otU+OpmvrzjUqKjQuk6cgE/rM34LHD1LqnH6NajSU5N\npse0HqzcvzILgxXJXZTIRSTzmvaHDm9T6tLKTLq3MZWL56f3uMVMWrQ9w5usGl2VcW3GERUSxczt\nM7MwWJHcxTI6O1NOa9CggVu0aJG/wxCR89m7hvj9W+g1vyizNx7giTbV6HNVTIZHox9NOEqB0AKY\nGcKHIx8AACAASURBVCeTThIZEpnFAYtcnMxssXOuwfnqqUcuIlnrh+eJmPEEo+64nHaXl+Tlqet4\n8Zu1pKZmfEpXM2Nn3E5u/OJGPt3waRYHLBLYdB25iGStTu/BiQOEhkcwpEtdiuYLY+TsLRw8kcir\nN19OSFDG+g/R4dFUia7Cf+b9hyMJR+hZs6euORdBPXIRyWrhBaFIRQA884byXMwGHrmhCl8s3ck9\nYxZxMjE5Q5uNCI5gyLVDaFuhLYOXDOb1Ra+T6lKzMnKRgKRELiLZIyUJ1n+LbZjKg9dWZmCnWsT+\nup/b31/A4RMZuxd5iCeEgc0Hclu12xi7Ziwfrf0oi4MWCTwa7CYi2SfxBASFQVAwpKYyfe0++k5Y\nStnCEYzt2ZDShSIytFnnHJ/9+hltK7TV4DfJtTTYTUT8LzTKm8RPHYPRbWiVMotxd1/JvuMJ3Pzu\n3AzPAmdm3FzlZiJDIjmRdIKBCwZyPDFzM8qJBColchHJfp4gb1IPiaRhTBEm9mlMcmrmZ4EDWL5/\nORPXT6Tn9J4cjD+YRQGLBA4lchHJfqFRcMdncFk7AC7Ln8Dn9zWhcGQIXUcu4Md1GZ8FrkmpJgy5\ndghbjm6hx7Qe7IzbmVVRiwQEJXIRyRmnLxXbuQQG16bs7u/49L4mVC6en15jF/PZ4h0Z3nTzMs0Z\nccMIDp06RPep3dl0ZFMWBS1y8ct0Ijezvma2zsxWm9mrvrIQMxtjZivNbK2ZPXlG/a2+8mVmptFr\nInlN8cugblco34yi+cKY0LsRjWKieXjSct77OeMJuG7xuoxuNZoi4UUICwrLwoBFLm6ZGrVuZi2A\np4F/OOcSzKy4c26fmd0OtHfOdTGzSGANcI1zbquZbQUaOOcOXMh7adS6SC6Umgq/zSWhTGP+NXE5\n36zYTa/mFXiyzWV4PBmb7MU5h5mR6lL59fCvVI2umsVBi+SMnBq1fh/wsnMuAcA5t89X7oAoMwsG\nIoBE4Fgm30tEcpslH8KH/yBs92KGdqlLj8aX8n7sFh6ZtJyklIxN9nJ6trdxa8bR5ZsuTNs6LQsD\nFrn4ZDaRVwGam9kCM/vZzK7wlX8KnAB2A78BrzvnTg9NdcB3ZrbYzHr/3cbNrLeZLTKzRfv3789k\nqCJy0anTFW4aDmWuwOMx/q99DR6+vgqfL91J77EZnwUOoGPljlxe9HIe+/kxJq6fmIVBi1xczpvI\nzex7M1t1jkcHvHO1RwONgEeBieb97/CVQApQCqgAPGxmMb5NNnPO1QPaAA+Y2VVpvbdzboRzroFz\nrkGxYsUy9UFF5CIUHAZ1bvMOhDu6A1s6jr7XVealjrX4ecN+uo7M+CxwBUILMPz64TQv05z/zv8v\nI1aMIFAmwBK5EOe9aYpzrmVay8zsPuBz5/12LDSzVKAocDswzTmXBOwzszlAA2Czc26nb7v7zOwL\nvEl/VuY/iogEtLnDYNl4qNKG2xuWIzoqhH4fL+Of782ja8NyjIzdwq4j8ZQqFMGjrapyU93S591k\nRHAEg1oM4tk5zzJ8+XBuuPQGyhcsn/2fRSQHZXaw271AKefcs2ZWBfgBKAc8BlRzzt1lZlHAL0AX\nYBPgcc4d95XPAJ53zp33RywNdhPJ5VKS4NBmKPbH4LT5mw/S44MFJKT8+TgVERLEwE610pXMAVJd\nKhsOb6BadDXgjwFxIheznBrsNgqIMbNVwMdAD1/v/G0gn5mtxpvERzvnVgCXALPNbDmwEPgmPUlc\nRPKAoJA/kvjyT2D2WzSKKULByNCzqsYnpfDa9PXp3rTHPL8n8a83fU3/mf05lXwqS8IW8bdM3Y/c\nOZcI3HGO8jjgn+co3wzUzsx7ikgesOlHOL4LGvdl//GEc1bZdSQ+Q5uOT47np+0/0WdGH4ZdN4z8\nofkzE6mI32lmNxG5+Nz0Dtw+EYKCKVvw7B45wCUFwzO06Vuq3sKrV73KigMruHv63RyIv6ApLUQu\nOkrkInLx8QRBSAQkJzKx0FAeDv38rConTiUx+9eMJeHWFVoz7NphbDu2jR5TexCXGJfZiEX8Rolc\nRC5eniBKlCxLs7o1KV0oAgNKF4rgsdZVuaRgBN1GLeCVaesyNHlM09JNGXH9CDpW7ki+0HxZH7tI\nDsnUqPWcpFHrInmUc3/ccGXvaihSGYJDiU9M4fkpq5mwcDt1yxViSJe6lI2OzPDbrD64mqSUJOoU\nr5NFgYtkTk6NWhcRyV6nk/ipY/BhO5jyEAARoUEM7HQ5w26vy8a9cbQdHMuUFbsy9BbOOQYuGEjv\nGb2Zs3NOVkUukiOUyEUkMIQXgA5vQ1NvIic5EVJTaXd5Kb7t35yKxfPx4PilPPn5CuITUy5o02bG\noBaDKJe/HA/++CDTtuiqWAkcSuQiEjiqtYViVbzPZzwD/+sIKUmUjY5k0r2Nue+ainz8y3baD5vN\nuj0Xdp+mohFFGd16NLWL1eaxWY/xybpPsuEDiGQ9JXIRCUzFq0Oput6JZICQIA+Pt67G2Luv5Eh8\nEh2GzWHc/G0XNL96/tD8DG85nKvLXE3szlhSXcbuwCaSkzTYTUQC34FfIfZNaP0SRBTmQFwCD09c\nzs8b9tO6Rgle6Xw5BSND0r255NRkklOTCQ8O51jiMfKF5MNj6vdIztJgNxHJO3b8Apt+gGTvLHBF\n84Ux+s4reLrtZXy/di9th8SyaOuh82zkD8GeYMKDwzmVfIq7p93N07OfJik1KbuiF8kUJXIRCXx1\nbod+SyF/Ce/laks/wpNyil5XxfDZfU0I8hi3jpjP0B9+JSU1/Wchw4LCuKH8DUzZPIWHZj5EfHLG\npoUVyU5K5CKSO4RGef/uWgKT74dlHwFQu2whvunXjHaXl+SNGRvoOnI+e46m74YpZkbvy3vzTKNn\niN0Ry70z7uVY4oUNohPJbkrkIpK7lK4Pd0+H+nd5Xx/ZTv5QD4NurcNrN1/O8u1HaTN4Fj+s3Zvu\nTd5S9RZevdo7P/t/5v4nmwIXyRgNdhOR3CspHt5pBGWugM4jAdi4L46+E5aydvcx7mpanifaVCMs\nOChdm5u/ez7lC5SnRFSJ7IxaBNBgNxERCA6HFk9Dg57e16mpVCoWxRf3N+HOJuUZPWcrnd6Zy+b9\n6btpSqOSjSgRVYKU1BReWvASGw5vyMbgRdJHiVxEci8zuPwWuLSx9/W8ofDx7YSTyP+1r8H73Ruw\n80g87YbO5rPFO9K92f3x+/lh2w/cOe1Olu1blk3Bi6SPErmI5B3B4d5BccHee5lfX/0SpvZvTs3S\nBXl40nIGfLKMuITk826mRFQJxrYdS3R4NL2+60XsjtjsjlwkTUrkIpJ3NOwDnd739tTj9sG0pygZ\nnsKEXo0Y0LIKk5ftpN2QWFbuOHreTZXOV5oPW39IhYIV6PdjP83PLn6jRC4iecvpu6lt+hEWj4Zj\nOwnyGP1bVmZCr0YkJKfS6d05jIzdfN7pXYtGFOWDVh/QqFQjikUWy4HgRc6mUesiknfF7YN8xb3P\nf50BMddw+JTjsc9WMGPNXlpULcbr/6xNkXxh6d7k0n1LqVOsDnb6PwwiGaRR6yIi53M6ie9fDx/d\nDPOGUTgqlBHd6vN8hxrM2XSQNoNjmbvxQLo298ueX+g+tTsDFw7UDVckxyiRi4gUqwq3T4Qr+wBg\n8Yfp3uhSvry/KfnCg+n6wQJem76O5JS/T871L6lP9+rdmbBuAk/GPqn52SVHKJGLiABUaQWhkZCS\nDP/rDJ/dQ/VSBZjStxn/rF+Gt2du4tYR89lx+GSam/CYh0caPEL/ev35dsu39Puxn+Znl2ynRC4i\nciZPENS9A6q3ByAyNJhXb67N4C51WL/nOG0HxzJ15e40Vzcz7ql1D882fpa5u+Yye+fsnIpc8igN\ndhMR+TsrP4UN0+HGQWw7Dv0mLGX5jqPc3rAcz7arTnhI2tO7bj66mZiCMQCkulTd01wuiAa7iYhk\nheO74dhOCArj0iJRTLq3CX2ujmH8gt/oMGwOG/YeT3PV00l82b5ldP6qM9uPbc+pqCUPUSIXEfk7\nTfpCj68hKBgS4gid+wZPtqzAmLuv5OCJBNoPm834Bb/97TXnwZ5gDsQfoNvUbqw/tD4Hg5e8QIlc\nROR8PL7T5xumwcyXYPdyrq5SjG/7N+eK8tE89cVKHhy/lKPx5x6lXrNoTca0HkOwJ5i7pt3Fkr1L\ncjB4ye2UyEVE0qvWzfDAL1CuIQDFj61hTI/6PNGmGtNX76Ht4FgWbzt8zlVjCsUwrs04ikQUofeM\n3qw7tC4nI5dcTIlcRORCFK3k/Xt4G4xqhWfWK9x7dUUm3dsYM7jlvXm8PXMjKalnn2ovma8kY9qM\noXv17lQuVDmHA5fcSqPWRUQywjlYOQkqXgtRRSEpnmPJQTz1xSqmrNhN00pFeOuWOhQvEJ7mJvae\n2MucXXPoVLlTDgYugUKj1kVEstPpe51HFfUm9Ul3UeDb+xnapQ6vdK7F4m2HaTM4lpnr96W5iXFr\nxvHc3Od4e9nb571Bi0halMhFRDLLObi0MZS5AvN4uPWKckzp24xi+cO4a/QvvDBlDYnJZ0/v+lD9\nh+hYqSPDlw/nxQUvan52yZBgfwcgIhLwPB5o2v+P11tiqbT6C77s/R9emrGNkbO3sHDrIYZ0qUv5\nolG/Vwv2BPOfJv+hUFghRq8ezbGEY7zY7EVCgkL88CEkUKlHLiKS1XYuhi2zCA/28HyHmgy/oz7b\nDp7kH0Ni+XLpzj9VNTP+1eBfDKg/gI1HNxKfornZ5cJosJuISHZIPPnHTViWjmVnhZt5aNIqftl6\nmM71yvB8hxpEhf35pGhCSgJhQWEkpCRwKvkUBcMK+il4uRhosJuIiD+FRnr/bpgGUwZQ+uA8JvRq\nRL/rKvP50h3cOHQ2q3Ye/dMqYUFhADwz+xnunHYn+06mPVBO5DQlchGR7HRZO+g5A6q0IjjIw7/q\nBTG+Z0NOJCbT6Z25jJq95awR652rdGZX3C66T+3Ob8d+81PgEiiUyEVEslvZK71/4/bB+y1ovGUI\nU/tfRfPKRXl+yhp6jV3EoROJv1dvWLIhH7T6gBNJJ+g2tRtrD671U+ASCJTIRURySlQxuPYZqNeD\n6KhQRnary3M3VmfWhgO0GTyLeZsO/l61ZtGajGkzhhBPCAN+GkBS6rnncRfRYDcREX/55mFIOsWq\nBi/S7+NlbDl4gr4tKtHvusoEB3n7WbvjdnPo1CFqFK3h52Alp2mwm4jIxcw5iCwCUUWoWaYQX/dt\nRud6ZRjy40Zue38+u454L0Mrma/k70l8xIoRfL3pa39GLRehTCdyM+trZuvMbLWZveorCzWz0Wa2\n0syWm9k1Z9Sv7yvfaGZDzMwyG4OISMAxgxZPwfXPAxB19FdeL/wFQ26uzppdx2gzOJbpq/f8Xj05\nNZmFexby1Oyn+N+a//krarkIZSqRm1kLoANQ2zlXA3jdt6gXgHOuFnA98IaZnX6vd33LK/serTMT\ng4hIrrBhOiwZR/tqUXzTrznloiPpM24xz3y5ilNJKQR7gnnnundoWa4lr/zyCkOXDtX87AJkvkd+\nH/Cycy4BwDl3+qLH6sCPZ5QdARqYWUmggHNuvvPugWOBmzIZg4hI4Gv2EDywEPIVp3yRSD5vcYhe\nzS5l3Pxt3PT2HDbuO05oUCivX/06nSt3ZsSKEby04CV/Ry0Xgcwm8ipAczNbYGY/m9kVvvLlQHsz\nCzazCkB9oCxQGthxxvo7fGXnZGa9zWyRmS3av39/JkMVEbnIRRXx/t08k5BPu/F0ubWMvusK9h9P\n4Mahc/jkl9/wmIfnGj9Hz5o9qV6kun/jlYvCeW+aYmbfAyXOsehp3/rRQCPgCmCimcUAo4DLgEXA\nNmAukHKhwTnnRgAjwDtq/ULXFxEJSDEtoMsEqNKaFh4P03pV46Gvd/D4ZyuJ/fUAL3WqxUP1H/q9\n+pK9S6gWXY3IkEg/Bi3+ct5E7pxrmdYyM7sP+Nx3mnyhmaUCRZ1z+4EBZ9SbC2wADgNlzthEGeDP\ndxAQEcnrzKBaW+/zU8coNqEt4y5rz7uV7uLNGRtYvuMIQ7rUpW65whw+dZh7v7+XyoUq8/Z1b1Mo\nvJB/Y5ccl9lT618CLQDMrAoQChwws0gzi/KVXw8kO+fWOOd2A8fMrJFvtHp3YHImYxARyb1Co6DB\nnXhqdOCBFpWY2Lshqanwz+HzePenTRQMLcTA5gNZd2gdd067k70n9vo7YslhmU3ko4AYM1sFfAz0\n8PXOiwNLzGwt8DjQ7Yx17gdGAhuBTcDUTMYgIpJ7eYKg+cO/T/Naf/uH/Fj9G1pXL8or09bRY/RC\nahVuyrst32XPyT10n9qdrUe3+jdmyVGZSuTOuUTn3B3OuZrOuXrOudMj1bc656o65y5zzrV0zm07\nY51FvvoVnXMPOl0/ISKSfvFHCEs8zNDbG/BSx1os3HKItoNjOXW8Ah+0+oD45Hgmb9KJzrxEU7SK\niASa1FTweODYbvbPHUO3NVeybl88fa6K4Y5mhSid/xI85iEpNYkQT4i/o5UM0hStIiK5lcd36F71\nKcUWDWLyHeXo2rAc783azIPjNrPj0Cl2x+2m4+SOzPxtpn9jlWynRC4iEqia9IX75xJWvBIvdqzF\nhNawZf9x/jEklp/XH6VAaAEG/DSAyRt1qj03UyIXEQlk0THev9t/ofFPt/Pzdb9RpUR+npi0iRIn\n+1O/eAP+PeffjFk9xr9xSrY573XkIiISAErXh/bDKFzrZj5pEsbb01cwKHYnFYreSsPqkby+6HXy\nh+anU+VO/o5Usph65CIiuYHHA/W6QUgEwS6Z/tv7MbvuTOJOGbPmtKJ5ke60LJfm/F4SwJTIRURy\no4rXUbrW1Uzt35ymFYvx7ezqPPzJBvYcO8o7y94hMSXR3xFKFtGpdRGR3CY4FFo+B0ARYFTdTSxn\nIbetb8GNH8RyKnoUS/ctZXCLwZqfPRdQj1xEJJezPSuow3om9mlGVHIdEnb9kwW7F9Jzek8Onzrs\n7/Akk5TIRURyu9YD4Y7PqVUumil9avNCYTi5vSurD6zjjm97sOfEHn9HKJmgRC4ikheEhAOQb81E\nbtn7FsOb1CN51z38dug4P6zb4efgJDP0G7mISF7SsA+UvZKWpesxpVYjXv5fQZ6etJcNO1bR46oC\nxBQu5+8I5QKpRy4ikpeYQel6AMS4Hbx3/CHerbyI8evHctPkzny57mc/BygXSolcRCSvKlIJu+5Z\n2tzWl5db9YCkgvx7fn+e//ETAuWGWqJELiKSdwUFQ9N+EBlN58ur811oPsomRzDxtxe5ZfwQjp9K\n8neEkg5K5CIiAilJlChWhUk1ulI2/HLWJo2mzduTWb79iL8jk/PQ/chFRORPElMS+eand9k7bwNv\nnmzDI60v455mMXg85u/Q8hTdj1xERDIkNCiUjom76FVgAXWrbuaNJS9x54cLOBCX4O/Q5ByUyEVE\n5GytXybonu+4unYwoYUXsj3+OVoP+oHZvx7wd2TyF0rkIiJyNjOIjOb+OvfzWMkWHMi3i8jiQ+j2\n4SxembaOpJRUf0coPkrkIiLyt7q1HMRLlW7nWPgeSlcbw6Sf53PLe/PYfuikv0MTlMhFROR8PB5u\nbPokg1oM4pZKVzA3379pvPdj2g6OZcqKXf6OLs/TFK0iIpIu15S9hmtKNoaTqTQpfiU/xcbx4Pgl\nzNl4gGfb1SAiNMjfIeZJSuQiIpJ+wWG4a59h8JRbOB69hzc8Jdm+JIr2W+9k6O11qVaigL8jzHN0\nal1ERC6ImfHG1W8QFRzJqxEbqFjDw5H4JDoMm8O4+ds0vWsOUyIXEZELVq5AOca2HUepQhV4Lnkp\nT96cxD/LHOXXr9/gvnGLOXpS07vmFCVyERHJkOKRxfmw9YdUL1Kd77ZP5vlyS3giagoL122m7ZBY\nFm095O8Q8wRN0SoiIplyMsl7GVpkcASnDm1k/cli9B2/hNLHltDk2g7c36ISQZre9YJpilYREckR\nkSGRRIZEcjI5nrvnP8cPe0cx7YaDTAj5L4t+mMQdIxew99gpf4eZaymRi4hIlggPDqd6keqMXj2a\nl08sJ6n9UNp17Mqy7Ufo+NZ0fli7198h5kq6/ExERLKExzw83fBpCocXZvjy4Rwtey2vXn0r9UsE\nEzX6akb87wZmN3qAJ9pUIyxY15xnFfXIRUQky5gZD9R5gCeufIIft//If+f9l5hLClGkbgcuqdmC\n0XO20umduWzeH+fvUHMNDXYTEZFsMW3LNGoUrUHZ/GV/L5uxZi9rJ/0fiSlQof1TdG5Q9m+2kLdp\nsJuIiPhV6wqtKZu/LKkulaFLh7IrbhfXX1acXpXjaBi5m4c/XcGAT5YRl5Ds71ADmhK5iIhkq53H\ndzJh7QS6Te3GpqObibhtLE3+9TEDWlZhwbIVDH7jeVZuP+LvMAOWErmIiGSrsgXKMrr1aFJSU+gx\nrQcrD6wiKDSc/i0r82m9lQxIfI97h09hZOxmTe+aAUrkIiKS7apGV2Vcm3HkC8lHz+96Mm/XPABK\ndX6Z5B5TqV61Gi98s5Z/j5jEwbgEP0cbWJTIRUQkR5QtUJZxbcYRUzDmj0JPEAUq1GNEt/q83+w4\nL+7uxStvvc7cTQf8F2iA0ah1ERHJUakuFY95+5GbjmyiYqGK3gVJp9j7/RC6r67LhoOJPHh1efpf\nfxnBQXmzz6lR6yIiclE6ncTn755Px8kdGblypPe38ZBwLmnzGF/0a0GXusW5fm5Xxr71ODsOn/Rz\nxBc3JXIREfGL+pfUp21MWwYvGcwbi974faBbZGgwAztcRuFLa7L4eGHaDo5l6srdfo724pXpKVrN\nrC/wAJACfOOce8zMQoH3gAZAKtDfOfeTr/5PQEkg3reJG5xz+zIbh4iIBJYQTwgvNXuJgqEFGbNm\nDEcSjvB/Tf6PYE8whOWjbM//8fjBk+z4eCkzJgwiYX4+Wvd4kvBQzS5+pkz1yM2sBdABqO2cqwG8\n7lvUC8A5Vwu4HnjDzM58r67OuTq+h5K4iEge5TEPT1z5BPfXuZ/JmyYza8esPy0vVySSSX0a0/uS\n9RT/7VtuGjabDXuP+ynai1NmT63fB7zsnEsAOCMpVwd+PKPsCN7euYiIyJ+YGffVvo+P2n7EteWu\nBfjT9eShwR6q9f2c1Fs/4sDJJG4f9h0/fzVW15z7ZDaRVwGam9kCM/vZzK7wlS8H2ptZsJlVAOoD\nZ06oO9rMlpnZM2aW5t3mzay3mS0ys0X79+/PZKgiInIxu7zY5QCsPriau6ffzcH4g38s9HhoVjOG\nb/s356nCP9J0cX+eG/MNR+OT/BTtxeO8idzMvjezVed4dMD7G3s00Ah4FJjoS8yjgB3AImAQMBfv\nb+jgPa1eC2jue3RL672dcyOccw2ccw2KFSuWiY8pIiKB4uipo6w6sIoe03qwK27Xn5YVzx/OTX3f\n4tu67zB+g4e2g2NZsXadnyK9OJw3kTvnWjrnap7jMRlvsv7ceS3EO7CtqHMu2Tk3wPcbeAegELDB\nt72dvr/HgfHAldn14UREJPA0Kd2E9294n0OnDnnnZz+y6U/LPSGhtL/pNibd25iqbhNVP27K9E/e\nITU1b55qz+yp9S+BFgBmVgUIBQ6YWaSZRfnKrweSnXNrfKfai/rKQ4B2wKpMxiAiIrlMneJ1+LD1\nhzjn6DGtB1uObjmrTt1yhRn0wM3Miu7Mo0uL0G3UAvYdjT/H1nK3zCbyUUCMma0CPgZ6OO/og+LA\nEjNbCzzOH6fPw4DpZrYCWAbsBN7PZAwiIpILVSlchTFtxtC6fGvK5C9zzjoFChSmZb/3eLpzY5Zs\nO8jat9qx4es3czhS/9IUrSIiEhAOxB9g5f6VtCjX4pzLN+3cw77R3fk6viaRje/hsdbVCA0O3HnP\nNEWriIjkKu8ue5f+M/szcf3Ecy6vWLoEdR/7luAr7mLk7C28MGQY++eMgwDpsGaUpscREZGA8MgV\nj7Dn5B7+O/+/HE04yj217uGvVzCHhwbz/E21aFKpGJ5Ph3D4u53MDW9Gh/oV/BR19lOPXEREAkJE\ncASDWgyiXUw7hiwdwmuLXiPVpZ6zbuuaJajR7zNeLz6Q/pPW8Ngni4j/ddY56wY6JXIREQkYIZ4Q\nXmz2IndcdgexO2KJS4pLs27p6Hy8c+8/6HddZQqt/ICIj25k46oFORhtztBgNxERCTjOOY4nHadA\naAGSUpJIcSmEB4enWX/B+p18PfF9Jp5qyJNtq3Fn/WgsvGAORnzhNNhNRERyLTOjQGgBAJ6Z+wx9\nZvTheGLaN1NpWLU0//rXUzSvXJQPpvzEiVdrEPfLhJwKN1spkYuISEC7psw1rDiwgrum3cWB+ANp\n1ouOCmVkjwbc16oeU5Mb0PU7x7xNB9OsHyiUyEVEJKC1rtCaYdcO47fjv9Fjag92xu1Ms66Z0fWa\n2lzWZwzHw0px+8j5LBt+DymLxuRgxFlLiVxERAJe09JNGXH9CI4kHOH+7+8nJTXlb+vXLF2Qr/s2\n49Y6xYjbuYYvZs5j15HAnN5Vg91ERCTX2Hh4I8cSj1HvknrpXufLJdt59ssVWFAI714XTJOip6Ba\n22yMMn002E1ERPKcSoUr/Z7Ex60Zx+yds8+7zk31yvJVv2soFx3J/umvcfSzfpw6mfbAuYuNErmI\niOQ6iSmJfL3pa/r+2JepW6aet375olF8dl8T1lz5Mp3iHuOmEUvZuPcY7Ln4b9CpRC4iIrlOaFAo\nH7T6gNrFavP4rMf5ZN0n518n2MOTN9bm33fexP7jCYx++0VShzfHbf8lByLOOCVyERHJlfKH5md4\ny+FcXeZqXljwAu8tfy9d67WoWpyp/Zuzt8wNvJDUlb6zPBw7lQTJidkcccYokYuISK4VHhzOmy3e\npH3F9hQKK5Tu9YoXCOe9e66lSMuHmLp6L7cN/oaEwQ1gxaRsjDZjlMhFRCRXC/GE8ELTF7i1YJVf\n2wAACc1JREFU2q0ArD24lqTUpPOuF+QxHmhRiYl9GkOq46ejxZm0LZLU1Ivrai8lchERyfVO3+50\n/8n99JjWg4dmPkR8cvquG69/aWHGP9SOr6q+yqNzoMfohcT98Bqs/DQ7Q043JXIREckzikUW45EG\njxC7I5Z7Z9zLscRj6VqvYEQIw26vy0sda7F4y342xn7KnqXnHw2fEzQhjIiI5DnTt07nidgniCkY\nw3vXv0fRiKLpXnfD3uP0++gXtu47Qo+rLqN2vmNMi53H18crU6pQBI+2qspNdUtnOkZNCCMiIpKG\nVuVb8fZ1b7P9+HbGrx1/QetWuSQ/X/a9ms4NK/PerM3EzXiRFxJfJj8n2Hkknic/X8mXS9Oe7z2r\nqUcuIiJ51uYjmylXoBzBnmCcc7//lp5e9Z7/jpMn46hu21jiqvxeXrpQBHOeuDZTsalHLiIich4x\nhWII9gSz7+Q+bp1yK0v3Lb2g9Q+fTOIUYX9K4kCO3oBFiVxERPK85NRk4pPj6f1db2btmJXu9UoV\nirig8uygRC4iInleqXyl+LD1h1QoWIH+P/bnm83fpGu9R1tVJSIk6E9lESFBPNqqanaEeU5K5CIi\nIkCRiCKMajWKOsXr8GTsk0zfOv2869xUtzQDO9WidKEIDO9v4wM71cqSUevppcFuIiIiZ0hISWDI\nkiH0vrw3BcMK+i0ODXYTERHJgLCgMB694lEKhhUkISWB8WvHk+pS/R1WmpTIRURE0vDd1u8YuHAg\nT8Q+QVLK+edn94dgfwcgIiJysbqx4o3sO7mPQUsGcTzxOG9e8yYRwTk3Ij091CMXERH5Gz1r9eS5\nxs8xd9dcen/Xm6MJR/0d0p8okYuIiJzHzVVu5rWrXmNH3A4Oxh/0dzh/olPrIiIi6XBD+RtoXqY5\nEcEROOc4nHCY6PBof4elHrmIiEh6nf59fMzqMXSa3In1h9b7OSIlchERkQt2VZmrCPYEc9e0u1iy\nd4lfY1EiFxERuUAxhWIY12YcRSKK0HvGhc3PntWUyEVERDKgZL6SjGkzhoqFKjJg5gD2ndznlzg0\n2E1ERCSDosOj+eCGD1iybwnFI4v7JQb1yEVERDIhX2g+ripzld/eX4lcREQkgCmRi4iIBDAlchER\nkQCW6URuZp+Y2TLfY6uZLTtj2ZNmttHM1ptZqzPKW/vKNprZE5mNQUREJK/K9Kh159ytp5+b2RvA\nUd/z6kAXoAZQCvjezKr4qr4NXA/sAH4xs6+cc2syG4uIiEhek2WXn5mZAbcA1/qKOgAfO+cSgC1m\nthG40rdso3Nus2+9j311lchFREQuUFb+Rt4c2Ouc+9X3ujSw/YzlO3xlaZWLiIjIBUpXj9zMvgdK\nnGPR0865yb7ntwETsiow3/v2BnoDlCtXLis3LSIikiukK5E751r+3XIzCwY6AfXPKN4JlD3jdRlf\nGX9T/tf3HQGMAGjQoIFLT6wiIiJ5SVadWm8JrHPO7Tij7Cugi5mFmVkFoDKwEPgFqGxmFcwsFO+A\nuK+yKA4REZE8JasGu3XhL6fVnXOrzWwi3kFsycADzrkUADN7EJgOBAGjnHOrsygOERGRPMWcC4wz\n1g0aNHCLFi3ydxgiIiI5wswWO+canK+eZnYTEREJYErkIiIiAUyJXEREJIApkYuIiAQwJXIREZEA\npkQuIiISwALm8jMz2w9sy8CqRYEDWRxOoFObnE1tcja1ydnUJmdTm5wtq9rkUudcsfNVCphEnlFm\ntig91+HlJWqTs6lNzqY2OZva5Gxqk7PldJvo1LqIiEgAUyIXEREJYHkhkY/wdwAXIbXJ2dQmZ1Ob\nnE1tcja1ydlytE1y/W/kIiIiuVle6JGLiIjkWkrkIiIiASygE7mZhZvZQjNbbmarzew/f1k+xMzi\nzngdZmafmNlGM1tgZuVzOubsllabmNeLZrbBzNaaWb8zyof42mSFmdXz7yfIen/TJteZ2RIzW2Zm\ns82skq881+8np5lZkJktNbMpvtcVfJ95o68NQn3leblNPjKz9Wa2ysxGmVmIrzzXf3dO+2ubnFGe\n546xp51jP/HbMTagEzmQAFzrnKsN1AFam1kjADNrABT+S/2ewGHnXCXgLeCVnAw2h6TVJncCZYFq\nzrnLgI999dsAlX2P3sC7OR5x9kurTd4Fujrn6gDjgX/76ueF/eS0/sDaM16/Arzl++yH8bYF5O02\n+QioBtQCIoB7fOV54btz2l/bJC8fY0/7a5vciZ+OsQGdyJ3X6f8NhvgezsyCgNeAx/6ySgdgjO/5\np8B1ZmY5EmwOSatNgPuA551zqb56+3x1OgBjfevNBwqZWcmcjjs7/U2bOKCAr7wgsMv3PNfvJwBm\nVgb4BzDS99qAa/F+ZvC2wU2+53myTQCcc9/69iEHLATK+Bbl+u8OnLtN8vIxFs7dJvjxGBvQiRx+\nP72xDNgHzHDOLQAeBL5yzu3+S/XSwHYA51wycBQokpPx5oQ02qQicKuZLTKzqWZW2Vf99zbx2eEr\ny1XSaJN7gG/NbAfQDXjZVz1P7CfAILwH4lTf6yLAEd9nhj/vC3m1TX7nO6XeDZjmK8oT3x3O3SZ5\n+hjLudvEb8fYgE/kzrkU36nRMsCVZnYV8E9gqH8j859ztElNIAw45Zs28H1glD9jzGlptMkAoK1z\nrgwwGnjTnzHmJDNrB+xzzi32dywXi3S0yTvALOdcbA6G5VfnahMzK0UePsb+zX7it2NscE69UXZz\nzh0xs5lAC6ASsNF3RifSzDb6frPZifc3jB1mFoz3dOpBf8Wc3c5ok9Z4/xf4uW/RF3gTF/zRJqeV\n8ZXlSme0SRugtq9nDvAJf/S08sJ+0hRob2ZtgXC8PzEMxnvaL9jXmzpzX8iTbWJm/3PO3WFmzwHF\ngD5n1M8L351z7Ser8Y47yavH2HPuJ/jxGBvQPXIzK2ZmhXzPI4DrgcXOuRLOufLOufLASd8OBvAV\n0MP3/GbgR5fLZsRJo03WAV/i/U8OwNXABt/zr4DuvpGVjYCj5zhdFtDSaJO1QEEzq+KrdroM8sB+\n4px70jlXxvcd6YL3M3YFZuL9zOBtg8m+53m1Te4ws3uAVsBtp3//9Mn135002qRwXj7GprWf4Mdj\nbKD3yEsCY3wDLzzAROfclL+p/wEwzsw2Aofw/iPkNudsEzObDXxkZgOAOP4Yefst0BbYCJwE7vJD\nzNktrTbpBXxmZql4R2jf7aufF/aTtDwOfGxmLwBL8bYF5O02GY73FsrzfD3Qz51zz5M3vjsXKi/v\nJy/jp2OspmgVEREJYAF9al1ERCSvUyIXEREJYErkIiIiAUyJXEREJIApkYuIiAQwJXIREZEApkQu\nIiISwP4fuupACyCEC8AAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x10fe3bed0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"figure(figsize=(8, 8))\n",
"plot(x[:4], y[:4], 'o-', label=\"super frames\")\n",
"plot(x[:4], p2(x[:4]), ':', label=\"parabola\")\n",
"plot(x[:4], p1(x[:4]), '--', label=\"line\")\n",
"plot(p0[0], p0[1], 'o', label=\"point\")\n",
"legend()"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 2",
"language": "python",
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.13"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment