{ "metadata": { "name": "", "signature": "sha256:16740d4797a0b8049545d9fde770792babee50a593ce58576923df2dd397b21a" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "%matplotlib inline\n", "\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "from pymacula import MaculaModel" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "markdown", "metadata": {}, "source": [ "#Generating a random model\n", "`MaculaModel` generates a model with 3 random spots (by default), using default stellar and spot parameters." ] }, { "cell_type": "code", "collapsed": false, "input": [ "model = MaculaModel() # default\n", "ts = np.arange(0,500,0.05)\n", "plt.plot(ts, model(ts))" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 5, "text": [ "[]" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEFCAYAAADwhtBaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXm4HGWd7z99zklycrKHEJIQkkBCgKwYZHGIbBFEAQUc\n74hMZhC5wIMiRkGYK4uijDN4ryAYuOIgICCKDAxcR8FBkHUQAwQTloQtJCQhC9nOyckJZ6n7x6/e\ndHWlqmvvru7+fZ6nnuquevutpbvfb/2W931BURRFURRFURRFURRFURRFURRFURRFURRFURRFURRF\nURRFURRFSYuvAouALuC2gLILgLXAVuBWoL9j30jgAaADWAGc4frsPOB1YDvwGDAh4XkriqIoFeY0\n4LPATZQXjE8C7wMHAcOBx4EfOPbfYy9twJHAFmCavW+U/f5ziMhcC/x3alegKIqiVJTvUV4wfgl8\n3/H+WMTaABgE7ASmOPbfQVFQzgWeduxrAzqBqQnOV1EURUmRpghlCwH7pwEvO97/FdgLGIE0/D3A\nm479LwPT7dfTXZ/ttMvOiHB+iqIoSoZEEQwrYP9gJHZh2Gavh9j7trnKt9v78Nm/zd6uKIqi5ICW\nCGWDLIwOYKjj/TB73e6xz+zf5ijjtb/d4zhBwqUoiqLsTlAbHkiaFsYrwMGO97OBdcBmYDkiTlNc\n+19xfHa2Y98gYLJjf9hzUBRFUbxJ3H6GEYxmoBVp8JuBAfbazS+ALyNZUiOAKygGybcD9wNXIwHt\nucApwJ32/geQeMXp9rGuAhYjQrMbU6ZY7LOPxe9/b2FZlV86Oy26uiz6+rz3+20Pu2zaZDFsmIV8\nv6XLgAEWHR3F771S17xpk8Uf/mDxy19a3HijxVVXWXzlKxZ/93cW8+ZZzJ5tMWGCbI97jBNOKL3W\nlhaLceMsDj/cYtEi/89ldS9ef13Owet7cC5tbRZbtsQ7Rl+fxYsvWrz0ksU778h97umJf85p3ItL\nLvG/1sGD5bv+whcs3n03/jG6uy3ef99i6VKLJ56w+Pd/t7jlFosf/cji7bfT+f6i3ItXX5XrmzKl\nfLmdO6Vcc3Pw//yCC6TsT35SvtwPfiDlvvWt8uXWrZNyI0cGX8/IkVJ2w4bSe1EJvgP0uZYrkX4S\n7cB4R9kFSGqt6YfRz7FvBKX9ML7gOs484DUk4O3XD8P+5VoWWNY++1jW9u2WL4sXW9a991pWX59/\nmaj87GeWVSjI8ZuaLKutzbJGjrSsYcMsa+BAy2puln3DhlnWPffEO8Zvfyt1tLRY1je/aVmXXmpZ\n3/62teu6Fy+2dv2DK0Ffn2VNm1Y8frmlf3/L2rEj+jE2bSrWMXWqZQ0fXlrv//yf/p/N6l7ce2/x\n+F/8opzD179uWf/rf8n3cdllxf2//W28Y/zgB973ccQIyzrtNMvq6YlWXxr3Yt48OYdBgyzr6KMt\na8YMyxo71rL69Ss9x7PPjlf/P/9z8X/itZx4YqLT30WUe7FokRz7Ix8JLtvSImV37ixf7stflnK3\n3FK+3E03SbnzzitfbvlyKbfffsHnOGmSlH3rLXlfbDsbCwuwXn/dsubMkRtyxRXeN6y93bL22EPK\n/PGPwTc4LJ/4RFEsghrPj3403jH+9V/l8xdeWLrd/JEfeaSygrFhQ/GaPv95y7rgAsu68krL+vGP\nLevuu+V8XnhBGhiwrJdein6Mxx6Tzx5xRHFbV5dl3XCDbD/lFP/PZnUv7ryzKBZ+fOUrUuZHP4p3\njIkT5fOjRlnWhAmWNXRo6W/o5Zej1ZfGvfjYx+TYTz9dur2vz7I2brSsf/s32T9pUvS6OzqKDe7I\nkZZ14IGWNXeuZX32s7IGy5oyJdHp78LciyefDC777LNy7MMPDy47bJiU3bSpfLn586Xc7beXL/eL\nX0i5v//78uWWLpVyBx0UfI4zZxYfLi0rPcGIEvTODQccADfcAHPnwrXXwpe+BPvuW1rm5pvhgw/k\n9ZNPwnHHJT9uXx88/7y8XrUKRo+Gri5ZWlqgXz/o3x86OmDkSPjrX6G7W7ZH4a23ZD3V1Qtlr71k\nvW4dFTUzt9mpCZMmwb33+pebNw8eegjefBMOPti/nBcbNsh6772L2wYMgI98RF5v3Oj/2azuxc6d\nxfPwY7/9ZP3229Hr37ED3n1Xfjvvvw/NtqO3t1eue8mS8tftRRr34sMPZd2/f+n2QgH22APOOgsu\nuABWrIDt22HQoPB1v/wy9PTA7NmweHHpvvXr5Te+ZUuSsy/S3m5x+OFw1FHwpz/B0Uf7l/W7Zi8G\nDYKtW+XaR4wIrrPc78e53/ze0jjHsHVGJUrQO1cceSSceabckLPOgs7O4r6tW0VIDKYBTsry5dJ4\n7r03jBsnf/TBg2HUKBg+XH5I/frJj2jsWPmC338/+nHMH2bkyNLto0fLev36ZNcRla12svSwYeXL\njRsn67Vry5crd4zhw0u3D7YTq7dvj15nUsIIhhG4ONf8zjuynjSpKBYgr80D0Natu30sc4IapuZm\nmDxZXr/xRrS6X3xR1occsvs+892HEYz/+i+YMgUeecR7v2XB+efDq6/K+6VLy9cXVTAg+DcZtk4V\njApx7bUwZoxYEFOnwnnnwY9/DF/8YumT2Zo16RzPWBeHHx5c1jSeq1dHP45f4+m0MCqJsTCGuhOf\nXSQRDNNIuK+5rU3WO3ZErzMpYQRjzz1lHdUSAHlCh92tYyje623u3kkVIEzDdMABsl62LFrdL7wg\n6zlzdt/Xvz80NYkF0tNTvp4bbpAHwRNP9G64f/xjuPvu4vt2rwR9B1Ea49ZWWafVwJv9QfV1d4er\nD1QwPBk3Dv74R5g5UxrmW26Br38dfvc7GDIE7rlHyqUtGIcdFu7c4h7b74m+WhZGWMEYO1bWSawq\n9zUbwXBakJUijGCMGiXrOIJhXKZGdJxUUyjDNHTGwjBWUliMheElGBC+oevqKr7+5jdL9/30p7Bg\ngbyebo8lESS8URpjU8bcJz+qaWGEPceo1GQMw8m0aeILff55eO45MUEHDoSvfEVEA2Dz5nSO9ec/\nyzqMYIwZI+s41oBf42kEo9IWRliXlHGhxbnfflbVwIGyrkfBKHdfwz51ZkGYhmniRFkbKykMXV3w\nyitiRcya5V1mwAARyZ07y8dGnPflpz8Vi+Tkk+HXv4Zf/Uq2/5//I27jiy5K18KoBcHIysKoecEA\n+QEecYQsTszTWRp+4B07JGDX1OTtf3VjGr44x65Vl1QUH7SbIJdUXgVjjz1kvXGjJEU0RbDZ/b5n\n5zHTfkIMQ5iGadIkWb/7bvh6lyyRgP60af5iELXxXLAAFi6EW2+VxdSxcCF8+cvidYhSnwpGeepC\nMPxobZWbazKZjO8xDi+8IGbrrFnBDSekIxjuJ08jGHkNeptrjmNh+AmG018ctUFOShjBGDBAGr/t\n22UxVm0Yyt3XrP7wYYhiYUQRjCB3FERvPP/+7yXp5cYbYeVK+X9+9avF88vS3WPcWEnrjHqOYbIu\ns3rgqGvBKBTkz7hhg/w5kwjGs8/K+m/+Jlx50whEFYzubnmabmoqZggZnDEMy5Lrc2NZEvzfsEFM\n86CUvjCEtTBMimESC8PdeBYK4pbasUOWKCmcSQkjGCD3Zft2+a6jCIbfNTuPmXfBWLHC/7foJgvB\n6N8fZsyAn/3Mu1xY155aGOGo6aB3GJK4SZyY9L25c8OVN41A1OM6nzrdf8LWVmmcurv9n+Lfflv+\nPP/xH/DMM9GO7UdUl1SaMQyonlsqrGCY7zpqRlPeLYxy1z1smHxXO3aEj9+EEQzzUOcMasc9R+f+\ntBp3Z5mgOs13p4JRQ6QhGGvXwhNPSADt058O9xljHURt5ILcP0FuKZO2CGKip0FYl5QzFbS3N9ox\n/FxSUAx8VzpjKIqFAdkIRqVjGL29shQKpX1DvIgS+O7ulo6sUOyM6UXYhi5sYxzVwgjj7jFlwopQ\nNUVNBSMiSQTjD3+QPhef/KT8iT7zmfI9O52Yp+KoHc7KNZwQnCnl7AORVnZYWAujubnokgnKSnET\nRjDyamGY+xLV/Rgm6F1pC8OZXhrkZorSafHVV6XB23//8r+jtJ+2q5myGrUfRlB9cfphBMVZolLX\nMQxIJhiXXQYvvSSvhw6Fa64J/9m4bpSwFoafYDj7QKTVSzishQFy3e3tYg34iZ6bnh4ZTqVQ2D1u\nA9WzMMI+IcZ1SeUxhhGl4YxiRYdxR0F2glENl1TagpGHfhh1b2HEjSV0dBTHuvnJTyQl8MADw3/e\nBGcrLRjOp720BCOshQFFH3SUxt3UP2yYdxZUI7qkqtUPI0qjFOWhKKxgpJ2BVAtB71oSjLq3MOJm\nK8kAf+Kn/cpXoh83roUR5JIKimFkYWFEEYw4jXvQNdeKYMR1SeUphhFHMMJ8L1EFI2xAuRoD+6Ut\naiYm0t1dPuMsD4LRMBZG1D+z8f+HjVm4ycolFRTDqLZLqhEFI45Lqlz6NIRvlNImSqMUNrZkWcXB\n/2bOLF82TENnWcX7EhSkroUYRqFQKhppnGPYwHxUVDB8yKtgVNolZVnxLIygtEgnQdecd8GI45Jy\nuuG8niizekIMIgvBeP99ud6RI4sPPH6EuW6nWAQF5qvl7untLXY0Dco2C1unWhgVIKlguIcYD0tW\nWVLlXFI9PcV5JSAdwdi5U/6gAwaE6wSYhYVRrYH4snRJlQt4Q/UEI2y6KoR3m3V0yHrEiHQa+Cz6\nI6TdGEepL4s6VTBiUi0LwzmkRZQ+CUlcUhs2iEXgrisJUawLqK5LauvWdIcDDysYcQQt6HuuBQsj\njBvFuT9MH4esGs604g0Qzt0TRXid5epBMEZSOhf3GT7lBgDXAauBTcBCSoPqByFzdW8B3gBOdX3+\nHHt7O/B7YGyYCwiiWoJRKMRrSMJaGF6CYdxRZrjsNBrPqIIRtqeukzQEo7tbhtw+4IBS0UxC2oFV\nJ40kGGk3dGG/F2d9ac41kYWFEUaEonQurKZgLAS6gNHAmcDNwDSPcpcBc4DpwFT79eX2vhbgQeAh\nYARwLnAXsL+9/xjgGuAziEC9A9wT9WK8qJZgQLw4RlBDMmSINMqdnUVT32AEw6T/pmFhmDrCCkYW\njWcYwVi1SuaXeP/94jwTSamEYPiJZD0JRh4sjErHMML24YlSZxZzdkQlSDAGAacDVwCdwDNIwz/f\no+zJwI2IBbERuAE42953IGIxXI9MRP64Xdd8x2d/A7wGdAPfA44CPOYii0Y1BSNOX4ygxrNQKM61\n4e5ha97vt58E2rq6kv9goloYcRrPNCyMVauKr9Pq4Z6lYOQ1hlFLjXG1zzGtjCZnuUq74qISJBhT\ngR7gTce2lxErwgtnSKsJGA/4jd/Z5KjH8vgswIyA8wskrmBs2iTrSlsYQY0nFKf0dM92ZgRj7Nj4\n1+2mVgTDGeyvBcFIyyXV1wc33RR9bm0/ojwZV8slFaW+FtspbsbIyuM5OsulJULVsjAGA25PeDve\nIvAwcBEwChgDfA0RgjZgGbAeuAToB5yAWBBtjs9+HpgJDASudHw2Ec4c+Si+7aRZUpCNSwrEggAZ\nmdaJl2AkjWNUQjCC3DNh0jed2WiNJBg//7l0LI0yCkE5asElFSWgXCik3xhHiTdU28KotGB0AO6m\nYhgiGm6uAV4CFgNPI4HyHmAd4mY6FTgJWAssAO4F3rM/+0fgO8C/I/GLd+xjvIcHhUJh1xJEv37S\n4PT2RktxTTOGEfa4lpVMMEynvVq1MJLEMJz3OI3YjXkibWoqPqX6Ya45Tt8TP5EM2/HKzNPS1xf+\n2OXIi7un0vGBvFgY6Z5jgT/9KVw7GZYgwViOBKynOLbNBpZ6lO0CLkTcUFOQTKlFjv1LkOD2KOBT\nwGTgecf+mxAX2Bjgfvu4XseJTJzGsxpB7x07in0eyk325CcY79nyOm5c/OEq3NSKSyptwYiSiVNN\nCyPFtqDkeHnOkspHY1wbMYy0CRKM7UjjfTXiHpoLnALc6VF2nL0UgCOQDKmrHPtnAq12PRcDewG3\n2/sGIPGKAjABuAUJkHv+9S3L2rWEoVYEI+wQHJMny/r110u3v/WWrPfbL30LI8ywIKCCEZa0gt5t\nDqdtGunEteCSitoYh+lgmHbKaj4Ew2LWrPDtZBjCpNVegMQV1iOpsOcj2UwTELfReLvcZCTzqQO4\nDbgUeNRRz3xgDeKiOhY4HnFVgQjJ3XZ9f7bruSLmNe1G1Maztzfa+El+xBWMoGHBZ8wQN8lrrxVT\na7dskUB9W5tkUdVTDCNMfxanYKTR/yROj+dqWBjOtiDqHCReqIURXF8jd9wLM1rtZuA0j+0rKQ1+\nP0X5NNhv2YsXWxFXVyZEbTydf+Yw48D4ETWtNuip0zBwoAzi9tJLsGgRHHNMMUtmv/2Kc5lD8qft\nrPth9PUFHyPvFoZpQHp6wtcfJBhhYxjO+/zBB+G/Jz9qIYZRC41xtUVNhwZJQNTGMw13FGTnkgIR\nCYAHHpC1GT561ixZ10oMo6NDRGPwYP/gct4Fw5x3HMHws6qam2WxrPLpoM5AexodFmvJJVULQe9G\n67hXF8QVjCQptRA9SyqsSwrgi1+U9a23yvSXzzwj7w89VNa1kiUV5przLhhxLIww1mSYP71TMDZu\nDH98P9QlFb6+PAe9dXjzBERtPNPotAfRLYywLimAj34U/vZvpaE85BC4005DOPHE0jqSNp7mnLIS\njDAdFfMuGEksjFoWjLC9ietVMKp1jtUcS6ruZ9yD+nRJAdxxh/x47rFH3fr0p4sduNIKehs3xx57\nhCsfVzDKXXO9CcbOnbK0tBSvzYso/nyonoUR1ChVs+Ne2Dq14144GkIwovrz0xaMLFxSpv5f/hIu\nuED6ZHzuc8V9aVkYxtrKWjAaycJwfs/l+lHk3cKotksqz/EBFYwaploWRlZZUm7mzpXFSRpB764u\nOfd+/bynEfUiDzGMNNNqwzRKTU3S8FtWcZa1coS1JKMKhnM8rbjUkmCk1RibXv2FQnVmx8uiTudD\nTIrdMDSG4UW1LYwkfT8MaVgYTndU2B7FWVgYzsmo/IbAqKaFAcU/aJjRQcM+GFTTJZXm4IPV7rgX\nVKfTugjzO6+FoHfYMbSiooLhQdoWRljBCNN4hiWNGEbU+AVkE8MoFIqi4WdlOO/xhx9GG9fJi7iC\nEcYtlZWFkYZgRIkPhI1h5MXC8PtNRq2vGh33LCta0DtMnXFQwfAgrbTaqC6pPFsYYcnCwoBgt5QR\nDOM6S2plRBWMKKm1YWNVYRqmWsiSihIfqEZjXG1rIEydpi9OS0uwyzNsnXFQwfCg2i6pNCwM03C2\nt0ebU9xJJQQj7DWXGx6kp0f+FIVCcQrbSgtGFAujFlxSWcQw8tpxrxYEI2p9YeqMgwqGB9V2SaVh\nYTQ3y3SuEH+MITNc+ujR4T9TDQvD3N9Bg9LLDqtFl9QHHyQf5rwRg95xXVLd3f4B5aiiFnQvVTAq\niLMRCZMxUK0sqTRdUs564sYxVq+W9d57h/9MFjEMiC4YSTOlKiEYQSIZRTD69y8dNNOLJUuKadJ+\nZDk0SL0IRqEQ/H1X+xzD1BmHhhCMAQNk6ekpn5ppSKundxQLo68v+jAcQSR92jaCMX58+XJOsnJJ\n1ZOFkZZLyrKK52lE3c8t9corMs7YEUeUP2aUhsmkoPb1lbds0nZJZRXDCBtMjlJnns8xDg0hGBC+\nIentlYbbOeJrXKLEMNrbpQEYMiR4hrewJO2LUUkLIy3BSGvQxSzTatNySTkbJeM29BOM3/1O1kFz\nf0dN3QxjZdSbSyqLOtXCyBlhBcP59Bc2G8GP/v3lKcwEZcMcN42AtyHp0/a778o6roURxv3XyDGM\nsC4pv8bYuKMGDIBRo+S1n2CEjaMl8ef7kZVLKux3EzSBkgpGeFQwXKSVUgvyBBY2jpE3wejslOFG\nWlpgypTg8gbnsNxBjadlZRfDyHNabVouKXOOra1FwfDr7e0c+jyM+yhNwYjiTmlulv9NX59/dl+e\nGuO0gtRROheGJYsRa1UwXKQV8DaEdUtlKRhxAsCvvSYN+v77R/uRQni3VGen/BEGDgxumPMuGNXI\nkjIWRmtrMZ3YZLa5MfO9Q/nvJUvBiNpDOc9P70GNcR7OMQsLoyHGkoLwDYkJeKdhYUD4wHcWghHH\nn9/RATfdBE89Je8POyz6cQcMEDHYubP8GFRRxLkeBSNplpTTJTVunLxes8a77MqVxdc7dviPkptF\nn4QoQ4OYOnfulDq9zrPaHffC1FkL5xiHMBbGSOABZK7uFcAZPuUGANcBq4FNwEJKBekg4DFgC/AG\ncKrr86cCrwDb7PVnw1xAWKotGGFdUmml1DrritJ4fuc7cOml8NvfyvvTvCbnDSBoKAZDloKR57Ta\nLFxSUQXDj7T7DzjrzGvP7DjuHs2S8mch0AWMBs4EbgameZS7DJgDTAem2q8vt/e1AA8CDwEjgHOB\nu4D97f2jgbuBbwBDgUuAXwKjol6QH9USjDy4pKIIxhNPFF+ffz6cckr044Z1SUURjHI9vWvJwrCs\nbFxS5QSjs7M0GF5unK1qu6Sc5YIaz2r19M6iznqxMAYBpwNXAJ3AM0jDP9+j7MnAjYgFsRG4ATjb\n3ncgMBa4HrCAx+26TD1TEAvmEfv974DtwOSoF+RHtS2MIMFIc1gQQ9Sn7b4+ydcHuQ833xwvUywL\nwTAWhpelZrbVQlptZ6cEcwcODP7zp+WSWrWq9H0YCyOLLKm0Bs3Lw9N72J7ZeRa1OAQ1B1OBHuBN\nx7aXESvCC+fgwE3AeGBImWPPsF//1T7OyUAz4p7qsrenQrUFoxaypDZulMZk5MhkQf8sBaPWLYwo\nrscoLqmxY+X12rW7Z0GZ9GhDpQUj7QbeXHfaAhT2u45SZx7cZpUc3nwwElNw0o63CDwMXIS4kcYA\nX0OsiTZgGbAecTX1A04AjgJMSKsDOA/4NSIUd9vvPX/ahUJh1xKWagtGLQS916+Xtcm4iUtYwYjS\noz7vghE2rTbK8C9BmThOl1Rrq/xme3p274vx9tul72tdMGrh6T0PVtAddxSAAmedFb6dDCJIMDqQ\nmIKTYYhouLkGeAlYDDyNBMp7gHVAN2I1nASsBRYA9wIm2W8OcAvwcURQjgZuBWZHupoyaAwjGCMY\nUQYb9EItDH+iuB6juKTA3y3lFgy/78Wyig1/2NEGqumSakTBiJrinjZBgrEcCVg7u27NBpZ6lO0C\nLkTcUFOQTKlFjv1LgGMQC+RTSHzieXvfPOA54EX7/SLgz8AnvE7KsqxdS1ga0cKI2niaTl977pns\nuNUSjLY2ud/NzVIuiSleCy4pp4UBRcEwQ7oY3ILhF/SOOvMcZGth+NUZN2U1rQmUoDb6YSxYYAEW\nP/pRenO0BgnGduB+4GrEtTQXOAW406PsOHspAEcgGVJXOfbPBFrtei4G9gJut/e9jFgXxqL4iP3+\n5SgXU45qC0Y1Yxhhg971YGEUCkVXXJLU2qwtjLRjGAATJ8r6nXdKyxnBMEO8+AlGkoazltNqa6Gn\nd60EvQEuQGIN65FU2POB14AJiGvKjDQ0Gcl86gBuAy4FHnXUMx9Yg7iojgWOR1xVAH8ArkXEqR24\nD3FxOT+fiFpxSaXZD8MZwwhjjNW6hWHEOY1MqVoQDLdLygzh8tZbxTKWVXw/bVrp59xEfXJ3llWX\nVClpd9zLoq9IHMJ4KjcDXt23VlIa/H4K2LdMPd+yFz9+aC+ZEEYwLCu9oc0N1XRJ9e8vT59dXWLh\nmHPxox4sDEgnjpFVWq35nsNcc1Bj7HZJTbaT0N905DRu2CCW1tChRQsjyDWTZnqpc0yxNAQji7mt\n8xTDCLJYaqHjXl0QZhKl9nb5cbe1RUuxK0cYl1SUzlxRMfWZhqoc9WJhGOswaLKgcrgb4yCixjDS\nCHq7XVJGMJwWxlI72jh9evEeBrmkovz2g3z5TusibFyk3HU7g/JpzW2dtmA4RS1qttmHH3q3T7Xk\nkqoLWlvlBnZ3+/9h0nrCdhLGJdXeLrnzbW3RniDCEKXxNOmYoxL2r4+aVhvG/RdFMMwc5M4RWqOS\nVVptmoLhdkkZwXj77WJfjCVLZD1zZrFcFoJRqSfjPFgDUP66e3ul0TcjN4fBObeI129IBaMKBLkq\n1q2TddJ+CE7CuKRMw2YaujQxjX+YxjOt8wgjGJaVzdAgkI5g1IKF4T7HIUPkt7tzJ6xYIduMhTFj\nRrFcJYPeaTd0eehkF1Rn3BTYctZaFplccWgowQhqSLIUjHIuKXM+SZ/svYjSeFbSwti+XRrXMEOb\nQ3kLo6ND1kPsiFpSwejrix5kzINLCmDOHFm/+GLpeubMYjm/7yVO0DtIMNIO1sY5x0pPoBTHqgpb\np1oYFSRoVrIsBCOMSypLCyNs42lZlbUwos474jeWVG+viEihUCyTVDCc7qiwfveoQe8sXFIAhxwi\n6xdekGD34sVyboceWl0Lo95cUmFELWocNEzsRoPeFaQaghHGJZXWk70XYRtPE/AfNCh5wD8LwTDC\n6xYMpzvKBEFNTCQNwQhLWAsjiqsrqksKioLx3HPw7LNiLc2ZI/cnrGBkEcOoN5dUGPdRmoKRFwuj\nYSZQgvwKRiUsDL9rNqQpWlkIRmur/El37pQGzzR+xh3lnKgpLQsjbPwCwge9ozwpxnFJHX20BFqf\neqr4XR5/vKzDBr3rzSVV6Y57cc7RWT7PgtFQFoZJF62GYJSLYZjzqWbQO03RykIwCoWiG8eZtNBu\nj2qWpmB4uXqCCGthROmTEHbwQed5jhgBxxwjrrr77pNtX/iCrIMsjCw67uXBJeUUtUqkrNazhdFQ\ngmEaT9PfwE21YxjVdEnlXTDAO8utnIURtx9GHAsjrGCkaWH4ubf+6Z+Kr087TTKknOWyGFMpzy4p\nZ8qq13lmFcPIwsKII7xpDm+uLikH79lj5+69d3rHdLqkLMs7iJqlhRFWMCrtkjICWm7ObzdeFoY7\nQwqqa2EE/TmzdkkBzJsHjz4Kr74KX/pScXs1gt55cEmZ8t3dUqf7s1nFMNIUjLwMDdKQFoaXYHR3\ni2AUCrBiyJwOAAAgAElEQVTPPukds39/WXp7/RvQPFgYxuqqlGBE7ecA3r3WvSwMY7Vs3iz33c11\n18Fxx/n3fq81C8NL2ObNgwsvLL0vWQa9g+ID1XRJZVFnpbOk1CVVBcoJxnvvSUbJuHHRf4xBBI2e\nmoegt5lDwQyRnYSsBKOcheFsGFtapKxleYvCN74Bjz8Od3qNuUw8C6MaQe+o97AaQe+8dIrLqjFO\ny8XlLK+CkRNM0NsrhmGmsZw0Kf3jBglGli4pk2Lq97RtqAXB8LIwvILeEM6yCgr+ZmlhhJmgKK5L\nyo8sgt6VdknVQkBZLYw6oZyFYYZSqIZgZGlhBD1tG2pJMLwsDPdIvH6CYcqD/zg/WWZJxbEwgkar\nDXue1Qx6p+WSShpQ9rr2WrKCtONeBRk8WG5iZ+fuaa7GwjCT0KRJOcHYsUOWfv2iBYCjYAZTNFlg\nhi1b4MQT4Re/qA3BMC4pp/D5DRdurvn990u3r11bfO2XuVZrMYy0LIxa6rhXbXdPudiNWhh1QqFQ\nTJl1N57VsjBMeunIkeGHoYjK2LGydjaWAL/5DTzyCPzjPxYF08yZkISsBMOIgjNd1m/EW79rjiIY\naWdJRZ0borlZeq/39nq7E7MSDHVJha8vixiGV52aJVUlTMqse97jalkYpsHLwh1l8Gs8ne6anTvF\nujDnmoQogmHGfwrDmDGydl6HU3Cd+F2zsaSg1D3ldW5p9/Q2jX5zc/S5Icr1KA7beAYFvbOIYeTN\nJdVIMQwdrTYF/ASjWhaGabTTmhLWC+NmcjaWACtXlr4/8MB0jhdGMMyos1EaZS8RiGphOO9BFhZG\nOcFIcwA5y8qHhRG2p3e9uaRq/RzjEkYwRgIPIHN1rwDO8Ck3ALgOWA1sAhZS2jHwIOAxYAvwBnCq\nY9+ZyFzeZtkO9AEfCXcZ4fESjN5eWLVKXk+YkPYRw1kYWQqGX+PpFoy5c9M5XlYuqTQEw/k+TQuj\n0oLR0yNp4M3N4TKuoDpB77RHWc2LS6paMYxaCHovBLqA0UjDfjMwzaPcZcAcYDow1X59ub2vBXgQ\neAgYAZwL3AXsb++/G5kf3CwXAG8BL0W9oCC8BGPNGvkDjhkTrZEISznBMBk/ac7l7SZIMAYOFMvq\nnHPSOV4lBMOMCeQ3B3ueLYywDTz4/+njnKPzeyk3plKeg9614JJqZAtjEHA6cAXQCTyDNPzzPcqe\nDNyIWBAbgRuAs+19BwJjgesBC3jcrsurHoCzgF+EvIZIeAmGcUdlEb+AcIKR9lzeToxLyt14mrjN\nu+/CO++k18M9K8FoaxNh/fDDYl8ak7zgnofcTzDM8C+QjYVRLuid5tN2nHMsFMp/N7UU9I7aGPtN\notTbK5ZaU1P46VSdx08jtuSus5xQxomrpTmWVJBgTAV6gDcd215GrAgvnKG8JmA8YjH4HXuGx/aJ\nwMepoGBk2WkPwglGGsFmP0zj6Xy63rpVzmfgwPSHJMlKMAAOOEDWr74qnfa2bpU63EkDJhtu/frS\nDCPzXYO/hWG2m4Ejw5CVS8rP9RGnr4izvFccI0nQu9JDg6TVGOfFGgiqM8737Rx0MS2CBGMw4G7m\n2vEWgYeBi4BRwBjga4g10QYsA9YDlwD9gBOAowCvHJl/AJ4E3vXYB0ChUNi1RMUIholZQLYBbygv\nGGZbJSyM1auLrghnVlja6bxBU2JCfMEwI6++8kpR9MeP3/0a+veXvhh9fUUro7c3nIVh+ui4OwOW\nI0yWVJoWhmnco2SZQfk4Rl6GBqnkwH5ZNO5pWxjOBIdoglGgu7tA6XN8MoIEowNwP/sOQ0TDzTVI\nzGEx8DQSKO8B1gHdSJD7JGAtsAC4F3jPo55/AO4Id/rRmThRzM9Vq4pfbB5cUllaGEOHyhP4jh3F\njmxZphE7f/hevnKILxgHHyzrJ54QNxr49x3Zd19Zv/22rNesKbU2/CwMIxhpWxhR+mAY0nRJOct7\nWRi14JJKO4Zh6ov6JJ6FqPnV6XzQiOI2y4IgwViOBKynOLbNBpZ6lO0CLkTcUFOQTKlFjv1LgGMQ\nC+RTwGTgeVcdRyKxjvvKnZRlWbuWqPTvLw1JX1+xIVm2TNZTp0auLhTGeqiWhQEwebKs33pL1ibg\nnUVWWFNT+Ml/oj4hn3yyrP/zP+EPf5DXM7wcm+x+zUYkzffsZ2HkySUVJBhRn2LDCEYWQe+8jlYb\n93dYSQsj7ndtWRZ77WUhjp50CBKM7cD9wNWIa2kucArgNc7nOHspAEcgGVJXOfbPBFrtei4G9gJu\nd9Xxj4hYlJluKDmmwVi+XNZGMIx/PG2qbWFAsfF8045GZW1VBcUx4j4hT5oEn/2sWAHXXy/bzDzW\nbsw1mwcDIxjT7By/jg5vCyiJhZHnoLezfKUsjLxnScW9j1n29E7rHOOcQxBh0movQGIN65FU2POB\n14AJiGvKOAQmI5lPHcBtwKXAo4565gNrEBfVscDxiKvK0Ap8ngzdUYb97WTe5culp/CGDeKvTnPi\nJCdGDLwG/6uUhTHFthHN0/Zrr8k6K5EMEow4HfcMN91UFLrBg4tWh5v99pO1uebXX5f1tGnyR/Kb\noyRODKPSabVxG5G0g95BHfeyckmlZRHEvY9OC9r90JG2hREnQ8pdZ1qE+eluBk7z2L6S0uD3U8C+\nZer5lr340YX00cgc06N56VIJnoI0nFmN5WQ6lW3atPuse5W2MIxV9eqrsp7m1aMmBcoJRtwgnmHc\nOHjuObjnHvj4x/2zvNxWlfOaBw+W76OjY/c/YhyXVNZBb3eDnNTCqHTQOy2XVFKLIC3BaG6WpbdX\nvnPn9WVlYcT5r1TDwqg7DjtM1v/93/DnP8vrQw/N7nj9+0sD1du7u1uqUhbGzJmyXrxYnqBXrJAn\n3SlTyn4sNuUEo7tbRKOlJdrTtpMxY2DBAvjoR/3LTLeTv5cskT+1saqMYIB34DuroHfeXVJZ5Pqn\n7ZLKi2A4P2OsZUNWMYw8WBgNKRizZolJ+8Yb8MADsu3ww7M9pt/8DJWyMGbMkB/w8uXwxz9Kg21c\nM1lQTjCS/AGiMHKkuKW6uuD55yVW1dws1qQRDK/Ady0IRlw3RTnBiOMmzDpLyu3uievKzKIxNm4x\nt2CkbWHkySXVkILRrx8cdZS8fuYZWc+bl+0xvQSjr684Y9wQv+6NKdG/P8yeLa9/+ENZH3lkdsfL\ng2BA0QK5/nq533PmiBCY+ISXYBirI04Mo5aD3klG6U3TJWWGdbes3Yd1T2phuH+PcbOknJ9x38us\nRE1dUlXk7LOLr+fNyya91ImXYDjFohL51SeeKOunnpJ1liKZF8Ew1/ib38j66KNlXQ0LI4t+GGn2\n9DYNXZTGM4ue3lC5NNgsLIw49xGyOcdK9/SuWz7/ebjySjjtNPjpT7M/npdgmPhF1u4ow5lnypMb\nSKD4pJOyO1YYwYjzVBeV004rfcr64hdl7RfDsKyiYMRpOGslhpGWkIe1MKrtz6+GYER54AD/BAd1\nSeWAQgG++124//5iNk2WmEwpp2BUYuBBJ1Onwq23wvHHS4ZRlk/4ebEw9twTbrhBxpa65BL4iD1g\nvp+F0dkpbpDW1mh/tnoIeiexMNIMejvLpy0Yfi6pWrAw8uCSipmjokTFWBjO6UUrFfB2ctZZsmRN\nXgQD4LzzZHHiF8OIK+K10g/DTzD6+uI17uZaent3TxmH/AiGedpPK97g/EyeXVJqYdQoXi6pSlsY\nlSSMYMR5YkoLPwujEoKR5mi1aQmG0+0RpT+Sc0TUPPd6NoLhdkHmycIIStONc45R3WJBqGBUCC/B\nqMR83tUiTxaGF1kJRqWypNJOq43byEHlBMOy4lsEpuE08SlDngQj6BzjPGBFyfQLgwpGhTC9kc3E\nP1CZ+byrRa0IhvuJ0wzfEnUGxKyD3mn53v2+lzQycbwEI04/DGd5p2CYaWnjdPjMUjCc4tvXF1/M\nszhHtTBqFK9Z7xrVwkhiYqdF2haGyT7r65PFizhptX65/mm7pPJmYXil/6bRcGZtYTjrizrUkLEG\n/NxmcSwMFYwaxWvWu0a3MKoZw0g76F0olAaAvYhjYaTtpvATjKwsjLiCYa7brzGOivm+sxaMJMLr\n913H6UjqrjMtVDAqxB57yB9r8+bij6pRLYw8uaTcgpFExIMC31kIRloWRl4Fw/m0naWFEaeB9wpS\nJ7XUWlrk9+N0xcUZPdmgMYwapampaGUYt1SjWhh5ckm5zX8zI6H5rqIQFPiOk1ab9pNxli6pNCc8\n8rruLATDXHccazdtCwO8r1stjAbFxDGMW6rRLYw8ptUawRgzJnqdWVoYaaWDVjrondTCcDacabh7\n3PfRvDe/hyhkIRhe5xlnuH13fWmhglFB3HGMRrUw8uCS8othGOsvjmAEZUql6ZKK+9RZK0FvrwBw\nkk52AwZInOnDD0u/H/P950Uw1MJQduG0MCwL1q+X93vuWb1zyopacUm5BcNM5xpnMMpKxjCyEowk\nvYndgtHbKxljTU3RB9f0uu4kDWehUPycs4E333+cOr0y2NKyMNK67rRjGDo0SAUxjdCKFTJS7fbt\n8gOp5NAglaJWXFLt7TII5A9/KEOfr14t5zVpUvQ6a1kwTMMZZ5h9PwsjrnUB3hZGEvcRyL3s6JB7\naa4zDQsjrcbdnKOznqR1VsPCGAk8gMzVvQI4w6fcAOA6YDWwCVhIqSAdBDwGbAHeAE51fb4NuAnY\nYJd5IswF1BJmdrs33yy6pcaNy25q2GqSd5eUc9rc//2/4fvfh9NPl22zZsUbbj4o6B2nH4Zfbn7c\nRsSrkYPiUPtxGs4sBcN5nkkadyj/9B6nTrfoQPL5berBJbUQmW97NHAmcDPgNRP0ZcAcYDow1X59\nub2vBXgQeAiZt/tc4C5gf8fnbwGGAwfaZb4e7VLyj5dgxMnGqQXy7pIaMED+1D09cN99pfv+x/+I\nV2ctWBimITMNmyFvFobXk3YWgpGkTuMZcE67nHTKglq3MAYBpwNXAJ3AM0jDP9+j7MnAjYh1sBG4\nATDTFB0IjAWuByzgcbuu+Y79pyBC8oFd5qU4F5RnzDDqb78Nq1bJaxPXqDfy7pKC4nAtZq7v666D\nn/1M5gqPQ9igd5S0Wq9GrrdX7muhEF10vRo5SNZw1pqFkZYImXtpOntC/iyMSvfDmAr0AG86tr2M\nWBFeOJ0rTcB4wO/WNTnqOQx4F7gacUn9FRGqumLQIBGI7m54wna47bNPdc8pK/JuYUBpssHw4XDR\nRXDOOfFnP8zCwnC6kMz81s4GJKo70zRk27aVzpedN5dUFhaGWyw//FDOubk53jma0QCc4ptUMGrd\nwhgMuJ5FaMdbBB4GLgJGAWOAryGWQhuwDFgPXAL0A04AjrL3gQjLDMQ6GQt8FbgDsTx2o1Ao7Fpq\njem2RP7qV7I+4IDqnUuW5D2GAaWCcdBByWNJWQhGc7PcS8sq3rckDUhLizQilpVeY1xpCyPuU/OI\nEbLevLm0vsGD43335VxScQXDiJCxWnp744tQoVBg1KgCpc/xyQgSjA7A7Y0bhoiGm2sQN9Ji4Gkk\nUN4DrAO6kSD3ScBaYAFwL/Ce/dkddpnv2595EnFbnRDpamqAI46QtUm/U8GoHsYlBSIYSQnb0zvq\nPMvuxjNpJo5p6JxxjEaIYZgRiN2CkfQ+elkYcWMYRtRMp17n2GZxLd80CRKM5UjAeopj22xgqUfZ\nLuBCxFqYgmRKLXLsXwIcg1ggnwImA8/b+/5qr91SaOGBZVm7llrDCAbID8BMGVpvhHFJVTuG4Uyd\nneaVxhGRLCwM2L3xTCoYTreUIQ2XlHtokCSC4Ux7NiQVDNMYmyHszdpsj0prq1z7zp3F33RSl5TJ\n3jOiFne4fSi2kxMnptdOBgnGduB+JLbQBsxFgtN3epQdZy8F4AgkQ+oqx/6ZQKtdz8XAXsDt9r4n\ngJXAPyECdSQiLo9Eu5z8c9xxxaFAPv7x+D/+vFMLFsaMGcXX0/2ichHIoqc37P4km6Rx96oP8ueS\ncjeckL5Lyoy0EHdonkJh93uZdNplt4WRVNSSftZNmLTaC4CBSAziLuB84DVgAuKaGm+Xm4xkPnUA\ntwGXAo866pkPrEFcVMcCxyNuKBA31GeBTyNxjJ/a5ZfHu6z80toK994L8+fDzTdX+2yyoxYE45hj\nZAiQ8ePhqKOS1xdkYcTphwG7NyKmwYvbEHgJhmmY4kwX7NfTO4lguK8Zitcddygdt0vK1J1kaB73\nvTQTpMUdvcEtlGYdx8IwpCkYYRL8NgOneWxfSWnw+ylg3zL1fMte/HgV+JsQ51PzHHecLPWM30xx\nzm3VdkmNGgXLlsl5pHEuYV1SUWeLcz8ZZyEYSZ62/R4OkgiGaSC3bCkOL5K0gXe7pNIQDHOemzZJ\n2nxSwfCzMPIiGDqWlJIJtWBhgDSeaQlXVjEMt2AkdVOYz23cKGvLSiYY5ntMUzD69ZM4QF9f0QWX\nlmCYa01DMEzHWzPKsREMZ0JFFNwWxrp1st5rr3j1gQqGUgPUimCkSVZZUmlbGO7pgrdtE5EbPDie\nePqNT5VEMGD3p+2kozub6169WtZJG3coHYF6+3bJfhwwIH58yVgm778vQm7ONUkH39Gj43/WjQqG\nkgm1kCWVNlkFvf0EI66bwj0vS9Lgr9f825CeYGzeLHV3dooox22MTSfZVaukMTajLSTpPOu8l86h\n8eP26Rk+XOJI27fL9+Iccy4uaXYOVsFQMsE5V4A7+7neLYysBSOtJ23TGCUdZj/Iwoh6vQZzfRs3\nlloDSRrjQYMk22rrVli5UrbHGcre4LyXb70lr/fbL359UEz3fued4jkmEYzx44PLhEUFQ8mElhZp\nSCyrtLeus8dyvVkYWQmGefI3jWZSN8W+dmrKq6/K+t13ZR234cxinnAoPhmvXJlO414oFOtcsSLZ\n3CcGIw6vvQZvvCGvzZhxcTHfz/LlsNTu8Zakn9D++weXCYsKhpIZXsM/9/SIaLS05KPnappkJRim\nQTONprEM9t47Wj2GGTOkM+Bbb4lYrFgh2ydOjFefX9DbCEbcyYRMw7liRTqCATB7tqzvv18sq+HD\nk7lsDjlE1i+9VBwfLmln3MMPl/VPflI8x7jfDcCBngMsxUMFQ8kMr1nt6tUdBeHnw4iaVmtcFCtW\npBMIbWmBU06R13fcAa+8Iq/jPon6xTCSzOIHpa4Z4+5J0nCCdJYF+N73iu+TjCG2xx4wc6ZY0WaY\n/OOPT3aOn/60rJ97TtannJKfOXNUMJTM8Jp7oREEI20LY9w4scbWrhU3SkeHxDWS5Oafc46sFy6E\n3/1OXh96aLy6snJJmfG9Fi2Cv/xFXh98cLy6DGeeWZyvvbkZvvGNZPUB/Mu/FAP73/xmchfQrFky\nodewYXDYYfDP/5z8HNNCp2hVMsPLwjDxjLSHXc4DWWVJtbSIW+GVV+CWW2TbrFnJnjqPO0784iaO\nse++MkVtHIIEI65Las4ceehYvlyWQgGOPjpeXYbhw+HFF+Hhh8WdNGtWsvpALIK1a8WiiusmdPPt\nb8uSN9TCUDLDy8KoZ8HIysIAOPZYWV97ben7uDQ1wZ13wsc+Jk+x99wj2+KQlYXRrx9cfHHx/Ze/\nnE6K6Nix8KUvpSMWhpEj0xOLPKMWhpIZjWZhZCkYX/0q3Hab5OfvuSece268c3QyZw48+2zyeoKC\n3kncj1dcIeLY05PculCSo4KhZIaXhWGG5q5nwUi7pzfIvClLl4o//6ij0u29m5Ssgt4gbigTqFaq\njwqGkhlqYRTp65MF4qcTT5pUOodHXsgqhqHkD41hKJnRaDGMckFvZ0ptXlIk0yKrGIaSP1QwlMxQ\nC6NIEndU3lHBaBxUMJTMaDQLo9EFI4ugt5IvVDCUzPCal7kRBMMr6F3PghEU9NYYRv2ggqFkhnPW\nNEMjCEajWhjqkqp/wgrGSOABZL7uFcAZPuUGANcBq4FNwEJKM7EOAh5D5u1+AzjVsW8S0IfME26W\nHPZ1VMJiRll1zstcz4JRLuhdz4LRr58E8ru7obe3uF0Fo/4IKxgLgS5gNHAmcDPgNeDuZcAcYDow\n1X59ub2vBXgQeAgYAZwL3AW4R14ZiswVPgS4JuT5KTnEzGdg5m+A+haMRrUwCgVvK0MFo/4IIxiD\ngNOBK4BO4Bmk4Z/vUfZk4EbEgtgI3ACcbe87EBgLXA9YwON2Xe561E1WJxjBaBQLo1EFA2RiIiid\n+8S81hhG/RCmcZ4K9ABvOra9jFgRXjizzJuA8Yi14Hf8Ga5t7wKrgJ8DMSeNVPKAUzDMrHsqGJU7\nn0rilUJtevXHnVJVyR9hBGMwsM21rR1vEXgYuAgYBYwBvoZYE23AMmA9cAnQDzgBOAowzx8bgI8C\nE4BD7Prv9jqhQqGwa1HyS2urCEN3d7EhadShQRpFMExGXG+vPBwUCvX5XdcCWbSTYQSjA4krOBmG\niIaba4CXgMXA00igvAdYB3QjQe6TgLXAAuBe4D37s9uBF5HA93rgq4ioDAp9NUrucLulttmPHkPd\nv6g6QC2M3R8MBg+uv57tjUwYwViOBKynOLbNBpZ6lO0CLkTcUFOQTKlFjv1LgGMQC+RTwGTg+ajn\naFnWrkXJN+5MqXoWjEbNkoLdp+M1a3VHVY8s2skwgrEduB+4GnEtzQVOAe70KDvOXgrAEUiG1FWO\n/TOBVruei4G9gNvtfYcBB9jntAcSMH8cb0tGqRHcmVL1LBhqYahg1DthM5IuQGIN65FU2POB15B4\nQztiUYBYDM8gbqzbgEuBRx31zAfWIC6qY4HjEVcVwH7A75F4yRJgB/79PZQaYdQoWa9fL2sjGMOG\nVed8skQFQwWj3gk7vPlm4DSP7SspDX4/Bexbpp5v2YsXv7IXpY4ws5CtWSPrRrAwNOitglGvaJ8H\nJVOMYKxeLZkz9dyQqIWhFka9o4KhZMq4cbJevbrYiAwZEn/+6DxjxKCRLQwVjPqmDv+2Sp5wWhhb\nt8rrenRHAfTvL+tGFAz3UPbmux7i12VXqUlUMJRMcQqGSa0dMaJ655MlRjA+/HD3ffUuGGZkYvMd\nm6y4PXSshrpCBUPJFGfQe+1aeT16dPXOJ0saWTD23FPWGzfKWgWjPlHBUDJl4EARje5u+MtfZJsK\nRv1hBGPDBlkbwTD9cJT6QAVDyZwDDpD100/L2jQu9UYjC4bpb+MWDLUw6gsVDCVzjGA8+aSsVTDq\nD7eFYVxTKhj1hQqGkjlGMHbulPXEidU7lyxpZMEYMkTm9u7slJTaVatk+/jx5T+n1BYqGErmHHxw\n6fspU7zL1TqNLBiFAkyaJK+XLZOsuEIB9tmnqqelpIwKhpI5hx9e+v7AA6tzHlnTyIIBRUvy0Ueh\nrw/Gji3eE6U+UMFQMqe1FU49VV4femj9Zs40N8tTdW+vLE4aQTDMg8DPfy7rWbOqdy5KNqhgKBXh\nppvgu9+Fu+6q9plkR6Hg39u7EQTjuONkvXy5rI88snrnomSDCoZSEcaOhSuvhKlTq30m2eLnlmoE\nwTj22KJbqq0N/uEfqns+SvqEHd5cUZQQNLJg9O8Pjz0Gv/61iMeECdU+IyVtVDAUJUUaWTBARide\nsKDaZ6FkhbqkFCVFGl0wlPpGBUNRUkQFQ6lnwgjGSOABZJ7uFfjPsz0AuA5YDWwCFlLq8joIeAzY\nArwBnOpTz5VAH3BciHNTlFyhgqHUM2EEYyHQBYwGzgRuBqZ5lLsMmANMB6bary+397UADwIPASOA\nc4G7gP1ddUwG/hZYE+UiFCUvqGAo9UyQYAwCTgeuADqBZ5CGf75H2ZOBGxELYiNwA3C2ve9AYCxw\nPWABj9t1uev5CXAp4DFnmaLkHxUMpZ4JEoypQA/wpmPby4gV4UXBVfd4wG+SxiZghuP95xFL5vcB\n56QouUUFQ6lnggRjMLDNta0dbxF4GLgIGAWMAb6GWBNtwDJgPXAJ0A84ATgKGGh/dghwjf35QAqF\nwq5FUfKECoaSF7JoJ4MEowMY6to2DBENN9cALwGLgaeRQHkPsA5xMZ0KnASsBRYA9wLv2Z/9DnAn\nsNJRn6qBUnOoYCj1TJBgLEcC1s4BqWcDSz3KdgEXIm6oKUim1CLH/iXAMYgF8ikkwP28ve84xCJZ\nay/7IIJyiddJWZa1a1GUPDFggKzdgmHe6+itSqXIop0M6um9HbgfuBo4B8l8OgX4mEfZcfZ6LXA4\nkiF1tmP/TCSdtgm4ANgLuN3eN89xLgXgL4gV8nDoK1GUHOBnYZjJo4ygKEotEiat9gIk1rAeSYU9\nH3gNmIC4psycWpORzKcO4DYk2+lRRz3zkXTZdcCxwPEUs6E22fWvt/f3ApsRwVKUmkEFQ6lnwowl\ntRk4zWP7SkqD308B+5ap51v2EoZy9ShKblHBUOoZHRpEUVLETzA0hqHUAyoYipIiamEo9YwKhqKk\niAqGUs+oYChKihjB6Ooq3a6CodQDKhiKkiID7bELnIJhWdpxT6kPVDAUJUWMYHR2FrcZ91S/ftCk\n/zilhtGfr6KkiBGMHTuK29QdpdQLKhiKkiLlBENTapVaRwVDUVLESzCMS0otDKXWUcFQlBRRl5RS\nz6hgKEqKqGAo9YwKhqKkiMYwlHpGBUNRUqRcWq1aGEqto4KhKCnS1iZrdUkp9YgKhqKkiMYwlHpG\nBUNRUkRjGEo9o4KhKCmi/TCUeiaMYIwEHkCmXl0BnOFTbgBwHbAamXJ1IaUz+h0EPAZsQeb2PtWx\nbxqwyP7cFmSq17khr0FRcoOXYJiBCFUwlFonjGAsBLqA0cCZwM1IA+/mMmAOMB2Yar++3N7XAjwI\nPASMAM5F5gff396/Gvg8sIe9/1fAfZGvRlGqjBlgsKenOEKtEQ8TEFeUWiVIMAYBpwNXAJ3Ik/+D\nwHyPsicDNyIWwkbgBuBse9+BwFjgesACHrfrMvVsBd6x9zUDfcDaOBekKNWkUNjdyjAptioYSq3T\nEnIs+tsAAAgrSURBVLB/KtADvOnY9jJwjE/5guN1EzAeGOJTtgmY4dq2BRGpNcBxAeemKLlk0CDY\nvl2WoUNVMJT6IcjCGAxsc21rx1sEHgYuAkYBY4CvIRZDG7AMWA9cAvQDTgCOAga66hgODENcUr+h\nVIAUpSYYOlTW7e2yNpbGQPevXVFqjCDB6ACGurYNQ0TDzTXAS8Bi4GkkUN4DrAO6kSD3SYiraQFw\nL/CeRz2dSDxkKjDT66QKhcKuRVHyhhGMbfajlloYSjXIop0MEozliNtqimPbbGCpR9ku4ELEDTUF\nyXha5Ni/BHFljQI+BUwGnvc5brN9bp0++xUltxjB2LpV1ioYSr0QJBjbgfuBqxHX0lzgFOBOj7Lj\n7KUAHIFkSF3l2D8TaLXruRjYC7jd3vcJ4GBEKIYCP0LcWM7YyS4sy9q1KEreGDZM1m4LQ11SSiXJ\nop0Mk1Z7ARJrWI+kwp4PvAZMQFxT4+1yk5HMpw7gNuBS4FFHPfORYPY64FjgeMRVBRK7uAcJei8D\n9gQ+E/OaFKWquF1Smlar1AtBWVIAm4HTPLavpDT4/RSwb5l6vmUvXtyH9rtQ6gQ/C0MFQ6l1dGgQ\nRUkZjWEo9YoKhqKkjF+WlMYwlFpHBUNRUsa4pNTCUOoNFQxFSZk99pD1xo2yNpbGUHePJkWpMVQw\nFCVl9tpL1u+/L2tjaRjLQ1FqFRUMRUkZIxjr1oFlqWAo9YMKhqKkzJgxsl63TubC6O6W2fZaW6t7\nXoqSFBUMRUmZoUNlsqSODlhrD9Kv1oVSD6hgKErKFApFt9SyZbJWwVDqARUMRcmASZNk/eKLsh4+\nvGqnoiipoYKhKBmwvz358NNPy9pYHIpSy6hgKEoGTJ0q6yeflLUJhCtKLaOCoSgZcMABsja9vMeO\nrd65KEpaqGAoSgYcfnjp+wkTqnMeipImKhiKkgFjxsAUxzyV06dX71wUJS1UMBQlI844Q9YjRsAh\nh1T3XBQlDcJMoKQoSgwuv1zE4ogjpCOfotQ6hWqfQEQsQOfyVhRFiUChsKupT9Tmh3VJjQQeQObr\nXgGc4VNuAHAdsBrYBCyk1Io5CHgMmbv7DeBUx74jgP8CPkDmD78X0GRERVGUnBBWMBYCXcBo4Ezg\nZmCaR7nLgDnAdGCq/fpye18L8CDwEDACOBe4C7C7ODEc+L/ARHtpB26LdDWKoihKZoQxTwYh1sJ0\n4E172x3AGuCfXGX/AvwrcJ/9/gz7/QRgBvDfwBBH+UeAPwNXehx3DvAnwDntjLqkbIyJqfdC74UT\nvRdF9F4UqaRLairQQ1EsAF5GBMQL5wk1AeMpFQn38Wf47DsKWBri/BRFUZQKEEYwBgPbXNva8RaB\nh4GLgFFI/OFriFXQBixDYhOXAP2AExBRGOhRzyzgCrusoiiKkgPCpNV2UOoWAhiGiIaba5BYxGIk\n5vFvwMHAOnv/qcCNwKWI++peu5yTKcDvELF5xuuEHOZVw6P3oojeiyJ6L4rovUiPMBbGckRYHP1W\nmY23u6gLuBBxQ01BYh+LHPuXAMcgFsingMnA8479E5FMqauBuz3q129eURQlHhVrP+8Bfom4luYi\nabEHeZQbZy8FJE12JfAJx/6ZQKtdz8XAW4h7CmBv+/030z99RVEUpVKMoLQfxhfs7RMQ19R4+/3H\ngXeA7cBr7N5f41rE6mgH/hPYz7HvKqDP3mcWd+xEURRFURRFURRFURRFUSpI2OFJ6oGvIskCXeze\n230e8Dri9nsMcQs6+Vdgo738S7anmTn9gVuR73sb8BJwomN/I90LkJER1iL34m3g2459jXYvDPsj\n/5M7Hdsa7V78CdhB0ZX/mmNfo92LXdxjL23AkUjg3Wt4knrgNOCzwE2UCsYo5Lo/hzSm1yK95w3n\nIT8Ok3zwir2tVmlDYlvmR34S0lhOQO7FVhrnXoB0lm21Xx8AvA98ksa8F4Y/AE8Cv7DfN9p/BOBx\n4GyP7Y14LwAZnmQnpam9dwA/qM7pVIzvUSoY5wJPO963AZ1Ib3yAZ4FzHPu/ROkPpB54GTgdvRcH\nAKuQIXQa9V58Afg18lBhLIxGvBePA1/22J76vaiVCZSiDk9SL7jzpqcj123oRO6JuQ/TXPv/Sn3d\no72Q38JSGvde3IS4F15BOsq+SGPei6HAd4EFlP5PGvFegDw8b0AE4mh7W+r3olYEI8rwJPWEe9S0\nQex+H7ZRvA+DEdeEc9/gbE6t4vRDOnPejnQmbdR7cQFyHZ8Avg8cRmPei+8hI0msQf4n5r/SiPfi\nUmBfxK10C/D/kC4Lqd+LWplxL8rwJPWE28IIug/u/cPsbbVOE+Jy6EISAqBx7wVI4/gn4DdI8kej\n3YuDkWDuR+z3BYr/lUa7F1A6WsYvkN/Ep8ngXtSKhRFleJJ6wm1hvIJct2EQMrzKK479Bzv218M9\nKiCZUnsiwbtee3sj3gs3/Si6pxrpXhwNTEJGkliLjA7xOeAFGu9elKOh70XY4UnqgWYkG+YHyBPD\nAHubyXo43d5/LRK4MpwHvIqYpnsjP4hzK3bW2fB/kUDcINf2RrsXeyJB3kHIb+GTiDvhUBrvXgxE\nJnMbjcS1fohYW3vQePdiGPJbaEUeqs9ErIQpNN69KMFveJJ65DvIMCnOxUwyNQ/Js+7EP6/6A3up\n9bzqici1d1I6ZIzpg9NI92IU4obajDQCzwOfcexvpHvh5iqKabXQWPdiFPJb2Ib8Np5Frt/QSPdC\nURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURSlDvn/FR0Z1q4Ag6UA\nAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 5 }, { "cell_type": "markdown", "metadata": {}, "source": [ "You can also change the default parameters; e.g. longer rotation period, different numbers of spots, different maximum spot size, etc.:" ] }, { "cell_type": "code", "collapsed": false, "input": [ "model = MaculaModel(nspots=5, star={'Peq':50}, spots={'alpha_max':10})\n", "plt.plot(ts, model(ts))" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 6, "text": [ "[]" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEFCAYAAADwhtBaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXmcFNW59789M+yrMIAC7ggqxAVxN6jBuBsVEqMxajRe\nNVyX+MZE87rGhJh7/XhdiVndjVFfxeW6RsWNqJEoLgiIC4uo7DgMwwAz0+8fTx2ruqe6u9bu6urn\n+/n0p3qqqrtP11Sf33mW8xxQFEVRFEVRFEVRFEVRFEVRFEVRFEVRFEVRFEVRFEVRFEVRFEVRouJc\nYCbQCtxe4twLgS+Ar4C/Al0dxwYA04BmYAFwUt5rJwBzgXXAC8BWIdutKIqilJnjgWOB31NcMA4D\nvgR2AvoD04FrHMfvsx49gf2BNcDO1rFG6+9JiMj8N/BaZN9AURRFKSu/prhg/A34jePvgxFrA6AX\nsAEY4Th+J7agnAW86jjWE2gBRoZor6IoihIhdT7OzZQ4vjPwjuPvd4EhwGZIx98GfOQ4/g4w2no+\nOu+1Lda5Y3y0T1EURYkRP4KRLXG8NxK7MDRZ2z7Wsaa889daxyhwvMnaryiKoiSABh/nlrIwmoG+\njr/7Wdu1LsfM8SbHOW7H17p8TinhUhRFUTpTqg8vSZQWxmxgN8ffuwJLgdXAh4g4jcg7Ptvx2l0d\nx3oB2zuOe22DoiiK4k7o/tOLhVEPdLHOrQe6IfGI9rzz7gLuAO5FsqUuxw6SrwMeBq4GzgTGAscA\n+1rHpwHXAhOBJ4ErgVmI0HQim1XdyGRksKDXQq+FE70WNnotbMy1CIsXC+NyJAh9MfBDYD1wKTJP\nYi0w3DrvGSQddjoyz+JjpOM3TAZ6AMuAe4BzgDnWsRVISu0UYBUwDjgx2FdKL/PmwZ/+BL/8pb1v\nrZvTTlGUquPTT2HyZNhrL/jVryCJOheN7JSPLNTeiGHGDLjkEnjVmXhs/esaG7Pcey8cemhFmpYI\ndCRpo9fCppquxQsvwHHH5Q4A33wTxo2L5v0dFkaoPt9P0FupAL/+NVxxhTzv0weOOgrGjAHI8r//\nC6+/DkcfDc88AwcfXMmWVo5q6BDKhV4Lm2q5Fh99BMcfL2Jx3HGwZImIxYwZ0QlGVPgJeitlZs0a\nuOoqeX755XIj3XcfXHqpPGbMgHPPhU2b4Ac/gNWrK9rcRPDBB3DbbfDKK8k06RXFSTYLZ54JTU0i\nGg89BCecIMc++aSybXNDBSPBvPEGdHTAvvvC1VeLheGkrg5uvBEOOAC+/BIuu6wy7UwCHR1w4YUw\nejT8+McwfjxMnAjr11e6ZYpSmEcfhZdegoED4a9/ld/0ttvKsQULKto0V1QwEsxH1rz4b3yj8Dl1\ndfCHP0AmA3/+MyxaVJ62JY3Jk+GGG6BLFzHr+/eHRx4RC0xRkkg2C1OmyPMrr4TNNpPnW1llVxcv\nrky7iqGCkWC+sCpxDR1a/LzRo+H73xfX1A03xN+upHHXXfDHP0L37vDsszBtGrz4ovx9223w/POV\nbqGidOa112DmTLEuzjzT3j9okGxXrKhMu4qhgpFgjGBsvnnpc//P/5HtXXfBhg3xtSlpzJ8PP/mJ\nPL/lFjjoIHm+664S9wHJMNN4hpI0br5ZtmefDT162PsbG2WrghERSTTV4sAIxhZblD533DjYbTdY\nuVJcMbVANgtnnQUtLXDSSXDGGbnHf/pTGDxYRnEzZlSmjYriRnOzxC9ABMNJr17QrZvE31pagn9G\nNgvt+dOrQ1KVgnHqqbUxYjRZT2bEUYxMBk4/XZ4/8EB8bUoSd9whrqdBg2S0lj+ZtWdP+I//kOe3\n3FLu1ilKYR59VAThgAPsmIUhkwlvZTQ3i6t6yBA70zIKqlIwXnwR7ryz0q2In3XrZNurl7fzJ06U\n7VNP2a9NK8uWwc9+Js9vuEH8wG6cfbYkBjz8sKQpK0oS+NvfZHtS/rqjFmEF4/HHYc4c8Tj86lfB\n3sONqhQMkHkIaffVG3O0Z09v5w8fDnvvLSOXp5+Or11J4Kc/FQvssMMK/+gAttwSDjxQEgKMC0BR\nKsmaNZKcUV8P3/ue+zlmABRUMGbOlO2ECXD44cHew42qFIwxY+Dzz+Heeyvdknjxa2GAbWX87/9G\n356k8I9/yATGnj3h1ls7u6LyMROh7r8//rYpSimefx7a2sQdZTKi8jGCEXQy7rx5sj33XPE4REVV\nCsYvfiFbk2WQVvxaGGDXlHruuXTGeTZssOdWXHGFPcmpGJMmyWjuH/+Ar74qfb6ixInpwIuN/M2c\njFWrgn3Gl1/KdtiwYK8vRFUKxgknyAWdNQvefbfSrYkPIxh+LIxddpFRy2ef2aOMNPE//wMffgg7\n7igzu70waJDMlm9rEyFVlEqRzdru4mKCMWCAbIMKhp+UfD9UpWB06wYnWsXP77qrsm2Ji40bpYNr\naJDZy16pq4NDDpHn//hHPG2rFGvX2jNjb74Zunb1/tojjpBtlOa5ovhl9mypCbf55jJXqBBhBKOj\nA5YuleeDB/t/fTGqUjAATj5Ztg8/nE7XSxB3lMEIRtpG09OnS1xn333t7+gVIxhPP53O+0WpDpzW\nRbHYWxjBWLlS5l8MGCCD6yipWsHYZx9JPfv0U0kfSxtBAt4GM9t5xox0dY7Gxbbnnv5fu9tuMqpb\nskRGebVCNiuDqiOPlCDrlCnQ2lrpVtUuXuIXEE4wTPwiancUVLFg1NfL2hAgOcdpI4yFse22crOs\nXCn+/rTw8ceyHTGi+HluZDKw337yPM1xLyft7TIPZdIk6ahmzJCKxgcdpCs1VoLmZim773QbF0IF\nIwaOPlq2afRLBwl4GzIZ2H9/ef7Pf0bXpkpjKvFus02w1++4o2znzo2kOYnnooukgnGPHjK58ZFH\nYOutpWz+j36ULuuzGpg+XeYD7bVX4YmmBhWMGPjWt6RzfP319JnZxiUVxMIAezSdphpK5sfjpVSK\nG6YEw+efR9OeJPP44yISDQ3w5JNwwQVw7LES1+rbV9xUDz1U6VbWFiZ+YeJpxVDBiIEBA2StiA0b\nZNSUJsK4pMC2MNIoGObH5JchQ2RrflBpZckSu67YNdfYMS0Qd97vfifPf/ELGfEq8ZPNeo9fQK5g\n+LUEzf3tpWipX6paMEDKPoCsWpUmwrikAHbfXdaDmDs3eC530ggrGGbEZVIO00h7O/zwhxK/Ouww\nu+y9k7POglGjJGEkrWnpSWP+fLneAwfCHnuUPr9HD/n9btzov2KtWhhFSKtghHVJde1q53m/8040\nbaokHR12mQQzC9YvJid92bJo2pRErrlGinMOGSIFOutcfuH19bLCG8B112ksoxwY6+LQQ+X6eyGo\nW6rSgjEAmAY0AwuAQqXeugHXA0uAVcBUoMFxfCfgBWANMB84znFsG6ADWOt4XOrlC4wfL9vXX5eJ\nbmkhrIUB9tKu770Xvj2VpqlJRKNvX/HLB6F/f9mmtTzISy/ZQnDXXbYLzo3vfldiQXPm2EsBK/Fh\nBMNL/MJgBMNvPSkjGMX+/0HxIhhTgVZgMHAycCuws8t5lwBjgdHASOv5ZdaxBuBR4DFgM+As4B5g\nh7z36Av0sR5TvHyBQYMkjbSlJV3zMcJaGGALxgcfhG9PpQnrjgIRGxDxSduoevly+MEPRFR/+Uu7\nplghunSRbB2onTTjStHSIlYfiJvQK9VoYfQCJgKXAy3ADKTjP8Xl3KOBmxELYgVwE2DWQNsR2AK4\nAcgC0633yn+fQC4yc+P/619BXp1Mwga9IV1ZQU1NsjWdfhAaGuR6ZrOSE58WOjpkUbHPP5fJeVdf\n7e11O1jDtU8/ja9tSeO55yRbbNw4OO+88vw2XnxREnPGjfNXqiOIYGzaJOfX1YUbXBWiVAc9EmgD\nnEbrO4gV4YZzsnsdMByxFgp99pi8fQuBxcBtQIlMZZs0C0YYl5QZYZhCZNVMFNcDoF8/2RoBSgPX\nXispmwMHStl3ry47k0WThvujFNmsuOu+/W147DH4979lFcaxY+N32T75pGz9uKMgmGAsXy7bQYO8\nx0r8UEowegP5P621uIvA08AFQCOwOXA+Yk30BOYBy4CfA12AQ4HxgFn6fDkwDtgK2MN6/4KrXWQy\nma8fYJeKePPNEt+miojCJWUEIw1ppEYwevQofl4pnG6pNDBjhiwmBhK3GD7c+2trSTDuuUcsr7o6\nWYHuhRfg4IMlY+7II+2ONmqc6bTlEAyTAThkSOd+MgpKCUYzEldw0g8RjXymAG8Ds4BXkUB5G7AU\n2IQEuY8CvgAuBB4APrNeuw54Cwl8LwPORUTF03hy7Fi5Ed59V1abSwNRjKhN0Gvp0ur32UfhogPb\nwkhD4Purr+D735dU2p//XDo+P9RCmjHA4sX2Gip/+pOso3LwwTLy33dfWQrgvPPi+ez58+GTT6Tz\nN54QrwQRDJMBGHWVWkMpwfgQCVg7q/fsCrzvcm4rcB7ihhqBZErNdBx/DzgIsUCOALYHSjmRXNuX\nzWa/foB0qmPGyA/n7bdLvGOVEEUH2aOHpNdu2lT9QhqVYKTJwnjgAZmkt8cedtl3P5j05DSvdd7R\nIZMYm5okdnHGGfax7t3F8ujZU1ZjfPbZ6D//scdke9hh/l1EYS2M/H4yCkoJxjrgYeBqxLV0AHAM\ncLfLuUOtRwbYB8mQutJx/BtAd+t9LgKGAHdYx/YCRlntGYgEzKfjbsm4MnasbJMw52DlSslUCSNe\nUbikID0+exWMzhgX7Kmn+lszxZAma6sQU6fKkqiDBol1ke+d2W47sTgALrlEBCZKTPkVs3SyH8II\nRqUsDIDJSKxhGZIKew4wB4k3rEUsChCLYQbixroduBhwrshwCvA54qI6GPg24qoC2A54ComXvAes\np/B8D1d22UW2SUgRvOIKKb9wilsumUeiDvJWe6dgLKSwgmFeX+0WF9jVe0eODPb6tM9LmTcPLr5Y\nnv/xj4U70fPOg6FDZYAXZX2tzz6T+WE9eviPX0B4CyMOvORTrAaOd9m/iNzg9ytAsRWWf2E93Pi7\n9QiMEYxKWxhr1shIBmTdhXXrgnX6UVsY1d4pRGVhmNf7LbeQRD75RLbbbRfs9ebeWLNGYlwRxkYr\nTlubWF7r18v2eLcezKJnT/i//1fiHLfcAt/7XjRtmDZNtocfHqwPCBPDiEswqr40iMEIxnvvRW9W\n+uHuu3NnnJuS3H7RIG8uUQtGGiwMk9kTtMhct27ix29rS8f1cPK730ma/ZZbwo03lj7/1FMl3vfy\ny9G5K421MmlSsNdXq0uqKhg0SH44zc2VnYj0t7/l/r1kSbD3icollRaffVSCYdJyq93C2LBBrNCG\nBujdO/j7pGVA4eSttyR1FuD2223XWzH69IGdrfoVUazIuHChlGrp3t1et8cvSXRJpUYwoPJxjI8/\nFp9lr162CRxUMNQllYu6pHIx9YUGDAjnSirngGLjRkllvfVWmWgYR2n11laxFtraJDYxYYL31xrX\nXlCvgJM775Tt8cfbv0G/9O4tA4J162SA4AUVDB9UujrrfffJ9rjj7LIcQUuLR2Vh9LGiTNW+JKcK\nRi5R1NYC+/6K+3pMmyar/R11FEyeLEHgbbeFJ56I9nMuv1wshJEj7XU/vBLVRMaODrjjDnlu1iUJ\nQiZjpz57KUDY0ZE70zsOUiUYptje+26zRMqA8VmedFL4HPeoLIy0+OxVMHKJSjDM9TD3WxzccIOk\nlX75pcyX+vGPYaedxPo++mhvMQYvzJwp5drr6yWW6PdeiUowXn5Z3OLDh8uqoGHw45ZatUrmovXv\nL/GpOEiVYJi1noO6gcKwZAnMmiU36YQJtt/Ub2likIyVqEphpKWDVMHIpVosjHvugQsvlOfXXSfW\n/1/+Iskpv/mN7P/pTyXWEJY77pDfzn/+p/9Z1WC7ccKul3LbbbI9/fTw9Zz8CEacZc0NqRKMoUNl\nW4nqrKbA2IQJEugyghHEwjDrk3fvHv6GS0sHqYKRixmIeAnoFiNOC+Opp2yXzHXXyep/ZkGn+nqp\ngTV1qvx9/vn2fR+U11+XbdCspKjmpbzyimyjSM/1syaGGSgPGxb+cwuRKsFwmpTlTq01vtijjpJt\nGMGIyh3lfI9q7yBVMHIx90iYDCmIz8J4/XVZpKmtTSbPuS0VCxLPGDNGshtnznQ/xysLF8p21Khg\nrw/zmzV89RUsWCAuoZ12Cv4+Bj8WhgqGT3r0kKyPTZvKWx9nwwapsw/RCEZUAW9ITwepgpFLVPdI\nHBbGu+9KULulRSyMa64pfr7Jbgyz8t+GDbBihVguQQO+UVgYRrS22y74ypBOVDBixgSby5lG+sYb\n8oMbM8YuLx0mnTWqztH5HmnpIDWmI0RlhUZtYcybJ6v9rVkj2YJu9Zvy2daqD7FgQfDPNf77LbZw\nX8fcC86Z70Exaa1RrXY30FoVyEv59c+s2t8qGD6Iwqz0i1l+8eCD7X3mhxhk5JZGl9TGjeFerxZG\nLkm8Hq+9BvvvL53mIYfA3//ubZRtgrRh1qQwccugs94hmr4j6nkQfjK3jIXhZ00Uv6RWMIJkJwXl\npZdke+CB9r4wgpE2l9RFF4mrMExhtyR2kJXE3Fdh75Ew96khm5Xg9cEHS6Xmo46CRx7xntrZ2Cjb\nMIJhXDbmvYLg9AoErQgetWCYRB4/gqEWhg/KbWFs2AD//Kc8Hz/e3m+CkbVuYbz5pmTIbNgg6ZRB\nibo0iM5LyX19mPvjoYekcN+GDRLEfuQRf0JmYg5hBCOKQVbXrpKZ2N4e/HrEZWF4yfxUl1QAyi0Y\nb74p6YCjR+cG28yN29zsf7SSFgsjm83Njpk7N9j7dHTkphqHIW2CkQQLw1RlveIKsTT8BnujEIyo\nYzpBr4exdEzsISxeLYzWVjvor/MwfFBuwTDxi4MOyt1fXy8meTbrP788LUHvRx+FV1+1Tf2FC4PV\nDzKde48ewQOahkoWH8xmpYLqNdfYk8yCEnUlgDDXwxTrC7LmA0Sz8l9Uv5mwgmFK8PTpU/w8rwwa\nJH3JihXF44Cm4OrWW4efu1WM1AlGGFdQEMwkHWf8whD05kuDSyqbhauvlue/+Y3c+Nms3Ph+iVJA\nu3aVrJ1Nm8T1UC4+/FBclnvvLWsvnH66lJAIStQuqTC/l/nzZbvjjsFeH0W9s6gtrqC/l+Zm2UYl\nGHV1dsaVyQRzwyymtf320XxuwfbE+/blpxy1cQwdHZJSC7Dffp2PBxWvOFxS5RJQwzPPyApmQ4bA\nmWfawchKC0YmU3631L/+BfvuK9aWMyj7wgvB3zOqoHfYWmMtLfLo1i14VVbTuQZx3xqS4pIyghF2\nQqUTL3EMIxgjRkT3uW6kTjCi8Ml6Zf58yagYOtQ90OSMY/ghSgujUj773/5Wtj/7mcQdjJ+60oIB\n5b0mixfDd74jvu1jjpF7xtQaCjNRLSnrg5j/Z2Nj8DLrDQ1yj3R0BG9HUiyuqF1SYMcxitXIUwsj\nIOUUDGNd7L13tG2J0sLo1q38Lph33xVXXd++cM45si8pFgaUTzA6OuCHP5TMmW99S7KJ+vePZqJa\nUiwMp2CEIaxbKuo047AuqSgtDHO/FFsYzgw+VDB8Uq76/iCuBihcGTOsYETlgjGZReWyMv74R9me\ncordEZiskZUr/b9ftQrGbbdJnGLwYHjgAejSRfabLJYg4mlIooURhrCCkbSgd5SCYdxMxSxSs2jc\n6NHRfa4bqRWMJFgYQWMYUbqkoLwumHXrpKQ1wNln2/vDdAim3dV0PZqapOgeyHoPzjTLMNaWQS2M\nXJIiGFEHvcEWDJNckM+KFeKu6tVLLQzflEswWlultn8mA3vsUbwtfmMYUbqkwO4gw5aP9sIDD0hn\nue++9oJWEK5DqEYL48YbJW7xzW/C97+fe8wsq7pqlVRz9YtzvZRKp9VGtS5H0lxSQfqPjo7o2uGk\nlIVhVhjdZZfwaeel8PL2A4BpQDOwADipwHndgOuBJcAqYCrgnMKzE/ACsAaYDxyX9/oJwFxgnXXe\nVl6+QD7lygqaNUviAjvvbK+LnE8SXFJQXgvjgQdke8YZufvNNaoFwVizRma3A/z6152DwfX1dgcb\nxEVn4lENDbabKyhhr4UZDBX6DXglKRZGGAFtaREx79kz2rkQW28t/+vFi93bNWOGbAsNXKPEi2BM\nBVqBwcDJwK3Azi7nXQKMBUYDI63nl1nHGoBHgceAzYCzgHuAHazjjcBDwKXW8ZnA/b6/DeWLYZSK\nX4DtkqpklhSUTzC++gqef15GOccem3ssSRZG3MvWXn+9XItvfct9fg6Em6wW5Sg2rIVh7u2wbQn6\nWzEkwSUVR8AbZFAwerSI0dtvdz7+/POyDbscrBdKCUYvYCJwOdACzEA6/lNczj0auBmxIFYANwFm\nnLkjsAVwA5AFplvvZd5nIvA+IhobgauAXRHh8UW5XFKl4hdh2hKXSypuwXjiCRn9jh/feU2CJAlG\nnNdj1SoRDLAnLrqRlOvRpYuMhtvags3Cj2ohp7AingSXVBwBb4MZmJqBqmHpUrEw6usLD06ipJRg\njATaAKf37B3EinDDaXzXAcOBQuGfOsf7jLbe19BifeaYEu3rRLkEw4uFkYSZ3lA+wXj4Ydkef3zn\nY1F0kGHXwjDEeT2uu06+42GHSanvQoS5HlH7ycN01lFZGGH/J0myMKIMeBtMP2MqYxt+/3txTx59\ndPg4khdKCUZvoClv31rcReBp4ALEvbQ5cD5iTfQE5gHLgJ8DXYBDgfHWMRBLJv9zmqzP70Qmk/n6\nkU85XFIrV0oAqkcPWTSpELVkYaxfL2s4Q/UIRtT3yIoVEuwG+NWvip+bFAsDwl2PqC2MSk/cS6qF\nccQREgt7+mlxY2azcPfd9mqGbkvgFusng1KqrmQzkB/O6oeIRj5TgP7ALCTm8RdgN8Aq+MtxiMvq\nYuBN4AHrPL+fUxQz5yDOjKA335Tt2LHFg461FPR+9llp9557wpZbdj4eRQeZdAG99lr5Xx95ZHFX\nJURjYSQhphOVhRFWMKJOMw7SjjgtjGHDJEbx/PNSMWDVKvjgAzl24YW5SyvESSkL40NEVJwVSnZF\n4g35tALnIW6oEUimlHNZ9/eAgxAL5Ahge8B45GZb72voZR2f7daobDb79SMfp082SMqiF7zEL6C2\nXFLF3FGQzBF1lNdj0SK4+WZ5Xsq6gGRej0paGElzSYURjDgsDJABSY8eUpPsgw+kxtRf/2pn5OVT\nrJ8MSikLYx3wMHA1cCaS+XQMsK/LuVbFE74A9kYypJzJld9A0mnrgMnAEOAO69g04Fok+P0kcCVi\nqXzo58sYuneXG7m1NZ5/nrEw9tyz+Hm1YmFs2gSPPSbPJ050PydMFkzSr0c2CxdcIO93wgkwblzp\n1zgL7vklSVljSbAwNm2SwWF9ffg04zDtiNMlBbD77pIl9cQTMkHvyCPDf1+/eFnqZDJwGxKDWAGc\nA8xB5knMRuZXfIZYBHch6beLENfTc473OQURnS7Ay8C3AZOXsQKYBNyCpNu+DpwY9EvFKRhmTQOI\nx8Job5eVy5wlPcISd2mQF18Uv+rOO8OoUe7nhPENO9fDiIKoBePWW2WVud694X/+x9trkjjzvVot\nDKc7Kqy7Psw8rjhdUoZRowr/xsqBF8FYDbg5GhaRG/x+Bdi2yPv8wnoU4nlEfEITZxxj4UJZGWzg\nQNhmm+LnBikNEuXNb4jbwjArrhVyR0HuD7Gjw9+M1CRaGF99JX7kv/8dLr1U9v3xj96Xx0ySS6ra\nLYwor0cYl1TcFkYS8LmYYnUQp2AYd9Ree5Xu0IOUBomjtECcgtHRYQtGIXcUiLugZ0/5Ia5f7+/7\nJU0w5s6VhAfn66dMgR/8wPt7RGFhJCFrLAlZUklZP6YcFkalSV0tKYhXMLzMvzAEccNUm2C8/rqs\nBLb11uJjLUbQ2lpJKz547732aw86CJ58UlbR80OYUX0aLYwoXFJqYcRPKi2MOIvteQ14QzDBiDqF\nFOIVDKc7qpTF1bu3uPOCJgEkJYbx3nuyve8+ODFgpC3MiDopFkZ7u/zGnKsYBiUpLqmkptUmhVRb\nGFF3kO3tMNNKFI5LMKJOqYX4BCObtdNpi7mjDGGr9ybFwgi7hjUkp4MM05Y4gs1hg95hcYqn32zU\nuNNqk0CqBSNqC2PuXLk5t9lGFsUphfOH6PXmqyaX1LvvwiefyLVwW9M8n7BrnCdFMJYtk61ZCCkI\nSbQw/F6PqOIXzjZUWkDr66X/yGb99x+14JJSwfCBiV94sS4g9+bz+mOsJsF46CHZHnust3LOaYhh\ntLfbJcnDLBqUJguj0q6guK6H34GNuqSqlLgFw0vA2+DXLRWnYER5PbJZePBBef6973l7TdIsjCCd\n06pV8t379w83aSpMG6K2MIK6g+LITqq0S8r5Pn7/N2phVClxCYZZqKTUhD0nfjvJagl6z54tLrqB\nA+Hgg729JoiF0dFh/x+jmsgYpnNavly2+eXbg7YhCSPqoOIVZTJCUlxSzvfx2xa1MKqUOEbUq1fD\n++9D167eXVLgv5OslqC3sS6OP15WA/NCEAvDOZpOwkRGs4Z1VIIRpA1JsTCidBUmMc3YryWsFkaV\nEkeW1D//Ka6IvfbyN9JNkksqytpJ990nz726oyCYhRF1/ALCXQ/TKYRdkjRNFkaUCzmZulB+SIpL\nSrOkqpQ4XFKvvirbAw7w97o0CsbLL0tqqSm57JUgFkbUnSNEV7coijYESd9MioUR5f/GOZejku1w\nvo+f+3TTJqkBV1cX3f8liahgeOSVV2T7zW/6e10aBeMvf5Ht6ad7d0dBMAsj6kl7YN8fGzZIjMQP\nUf1/GhrEvdnRARs3+nttGi0M5/tUuh1BLAxninGE6xUlDhUMD6xfLzO8MxnY162wexHSJhhLl9rx\nizPOKH5uPkFKnMdhYTgrAfu9R6L8/wTtIJNmYVR6xnlcS9b6aUctxC8gpYIRddD75ZdlFLjbbrDZ\nZv5e61cwkp4ldf31MjI/9ljYtlhtYheCzHyPI4YB4SerVVIw1MIoTzv83Ke1EL+AlApG1BbGM8/I\n9vDD/b+Nf5bpAAAgAElEQVQ2qIURl88+zOJbX34JU6fKc1PS2w9JsTCgugUjLgsjaDsqXQQxCS6p\nWkiphZQLRlQ+eyMYhx3m/7V+O8k4XFINDfLo6JDgXBCyWVk7uLlZ1hT2k1psCBP0jjqQmATBSMrI\nPmnB5qS4pNTC6EyqBSMKC2PxYlk/t3dv//ELSEYMA8K7pf7wB1ksqHt3uPHGYO+RlLRaSIZgBOkg\n29tF9DMZCZpHQdJcQZVuRxgLQwWjColSMB5/XLYTJgT7gaZBMFatgp/+VJ7/+c/+YxeGpKTVQvUK\nhlNAKz2RMWrrLymCEaQdKhhVTJSCYTKCJk0K9vokBL0hnJvOBP333x9++MPgbQiTVquCkXtulC66\noHNCktBRQzJcUpolVcVElSX15Zfw0ktiWXznO8HeIwlBbwhnYcyZI9t99gnXhiTGMCrZOYW1MKKi\nS5dgMa6kCEaSXFIa9K5CorIw7r9fRlyHHQb9+gV7D7+CYUYqUd94YQRjyRLZbrlluDZoDCOXpFgY\nQdsSV5ZUpQVDg96FSbVghMmSymbh1lvl+Y9+FPx9/AhGR0d8N14UgjFsWLg2OK+FV9eHuqRySZKA\nJsHCyGajt8o16F0YL4IxAJgGNAMLgJMKnNcNuB5YAqwCppK7Zvhw4HFgJfAFcDNglt3ZBugA1joe\nATL9hSgsjBdegHnzYOjQ4O4o8CcYzpuuLmIpD+Om+/xz2YYVjIYG6NYtt2R5KeISjKA5/3EIRpBO\nOgkWRhIEY8MGEY2uXf2VqYm6HSoYNlOBVmAwcDJwK7Czy3mXAGOB0cBI6/lljuM3ASuALYDdgAOB\nyXnv0RfoYz2meP0S+URhYVx/vWzPPjvcjehnHkZTk2zDVkJ1I4yFsWaNbP3Ocncj6PogSRhRQ3ot\njCSIV5DrEUfMT11ShSklGL2AicDlQAswA3gUOMXl3KMRq2ENIgw3Ac5qQ6OB+4GNwFLgaWufn/Z4\nIqyFsWkTPPecPD/nnHBt8WNhxBW/gGRUaHW+h9c4Rpon7iUphhEkCSAJFkYcWYVBXFJx/naTRKkO\neiTQBnzk2PcOnTt6gzMzvA5xQ5lL+AzwA6AHMAw4Angq7/ULgcXAbcDAEm0rSNgsqQ8+EFN3hx1g\n8OCgrRD8CEZSLYwoO0m/M9+TFtOp9EzvuC2MahWMOCwMdUl1ppRg9Aaa8vatxRYBJ08DFwCNwObA\n+UAWMP/Kq4Ax1vstBt5ErBWA5cA4YCtgD+v97y3UqEwm8/XDjbCzmhcvlu322wd7vZMgFkaaBSNo\nmnESBGPjRlncp74+mlnWSbQw/FyPJGRJxeGSClIkM4mCUaqfDEIpwWhG4gpO+iGikc8U4G1gFvAq\nEihvQ9xPGcTCeBARkEYkmP5f1mvXAW8hge9lwLnAoYhLzDddusgs2E2bpJSCX774QrZbbBHk03Nx\n/ghKZQYZCyNJLqlNm+QRVSdZzRaGUzij+A0mMYbhtS3ZbDJiGHG4pNTCKEwpwfgQyXQa4di3K/C+\ny7mtwHmIG2oEkik10zrWiFgOtwCbrGN3AEcGaV82m/364YZz9a4gbqkoBaO+XmIq2WzpzimJLqmo\nO0m/o7ekCkYUBAmuJsXCMIOxhgYZoEVBkNhBHC6pbt3kXt+wwfuAM4mCUaqfDEIpwVgHPAxcjVgG\nBwDHAHe7nDvUemSAfZAMqSutYyuQVNqfIKm0/YHTkHgIwF7AKKs9A5GA+XTcLRlPhMmUilIwwHsn\nmcSgd9SdZFosjCioZgsjKbGDOFxSmYx/8dKZ3jaTkUD1MuAe4BxgDhJvWItYFADbI1lUzcDtwMWA\nlWtEFsm2OgYRj/nABuBC6/h2SAC8CXgPWE/h+R6eCGNhrFol28bGMC2w8ZpKmnQLIwrUwrAJ4itP\nSppxUgQjrtprftuSRAsjDrzMMFgNHO+yfxG5we9XgGJ1TN8ACq2I/XfrERlhLIy4OoZSo2q1MDoT\nt2AEGc1WUjCiXjzJUK0WRlwC6ud/09YmA9O6OrvfSSupLA0C4SyMqM1cYzE05eeb5VELFoafiXsb\nN4qvvKEhurUfDNVuYVQ6hhGHayyMS6qSFoZzUBNhQlIiSb1gBLEwojZzTeHCr74qfl4tCIafiXtx\nmvnVKhhJ6CCd50UpXElKM/aTkFAr7ihIsWAkySWlgmHjx8KIaw4GVL9gVDqmE4crKMi6HHEJqJ+g\ntwpGCgjjkor6x+BVMEywfcCAaD7XSVIEQy0MmyDVe801SYqFEaVgGNejWYbWC0mwuMzvOo6BXtJI\nrWAkycIwN5IKhj8LI+2C0aWLPDo6JOe/Em0w+E0CiLuKsNd2xG1heLlPze866Jo51URqBUMtjFyS\nIhhqYeSSlFIpfqvVJm0+SCUtDONKVsGoYoJ2kHEsyOJFMLJZWzCiKCOeT1IEIykWRpBy3nF01kHn\npSTFwqhksNl5XlxptV4GNmphpICgLqnW1ugXZPEiGOvWid+2R494crmDlnxPgoURdecIuQJayQBr\nUAsjrhF1JYPezvertEvKzIVa66HWhApGCgjqkorDxPUiGKtXyzYOdxQkz8KotEuqrs6e21HJ+IHf\nEhRxu6Q0hiF4nTsFKhipIKiFEccN6EUw4oxfQHIEw8+IOs60WvB/TSptYXR0JGfiXlIEI64YhhEM\ntTBySa1ghLUwovwhVLNgRO0W8mNhmB+rCobgvDejXvM9CRP3grQjSRaGptVWMdVqYcQR8IbcGIaf\naseVnHsQ98it2gQjSVljSXGNxRX0NjEML4KhWVIpIOiIulIWxtKlsg27JGwh6uqkzj/4s7qi7iS7\ndpW5B21tUiuqGGvWyDZuwajkaDaIiy6OJICkVGetZgtDBaOKCeqSitPCaGoqPKr+8kvZbr55dJ+b\nT7XNPTA/xP79o/tsJ9VmYZSrVIoXCzQJM86dKfAqGOUhtYIR1CUVh4XRpYu8X3t7Yd+9EYyoFm1y\nIymC4TWOoS6pXOJMMzZL8Gaz3rLGkuCScq76F3U1YxUMd1IrGEmyMAAGDpTtihXux80qf7VgYXid\nvBe3SyroHIgoBxNJcUmBPxddElxSccUvQLOkCpFawQgb9I76Jhw0SLbLl7sfr0WXlFcLIy6XlJ/J\nWVB5CyPuNOMgnXUlXVJxCqgfC0OD3ikgbNA76puwlGBEvY64G0kRDK8WRtwjNyMYXlJ829vFVZPJ\nRJtKmhSXVNC2JMHCiCsJoK5O2tHWVvi81lZ5dOmS/tX2IMWCkZRSGIZigtHWJhZGJpMsCyPuTrJU\nRx23S8p0dl4sDKflGeWqaklySfkZVcctGJW+HpmMNwt05UrZDhyY/tX2IMWCkaS0WiguGIsWSec8\nbJid+hoHQYO8UXeSXiyMDRvk0dAQ/eQwgx8LIwlltON2Sfnx2yfhesTlDTB4EVCnYNQCqReMarAw\nPvlEttttF+1n5pOErCDn+xXrqJ3xi7hGbn5iGEnoION2SXmdrLZpk4h5XV30bpgg/5M4gt6gguFG\nagUjSWm1UN2CEfWI1ouFEbc7ytmOahGMuGM6Xi0M530RtZgnQcQNKhid8SIYA4BpQDOwADipwHnd\ngOuBJcAqYCrgLBA+HHgcWAl8AdwM1DuOTwDmAuuAF4CtPH4HV5JSbM9QTDDmz5dtUgWjEhaGqd4b\nV4YUVJ9LKs4138G7hRFnR50kwfAbw6gFvAjGVKAVGAycDNwK7Oxy3iXAWGA0MNJ6fpnj+E3ACmAL\nYDfgQGCydawReAi4FNgMmAnc7++r5BI26B21hWGC2Z9/3vnYrFmy3WWXaD8zn6QIhpeJe0ZYjdDG\ngVoYuXgNesdZ0ypJgmGusxm8uKGCkUsvYCJwOdACzAAeBU5xOfdoxGpYgwjDTcAZjuOjERHYCCwF\nnsYWnonA+4hobASuAnZFhCcQSUur3Xpr2S5cmLs/m4W335bnu+8e7Wfm4zeuU8lO0kxwjFMwkmBh\n+MkKitvC8OqSSopgxNkOgMZG2RpRcEMFI5eRQBvwkWPfO0jn74bTo1mHuKGsW4BngB8APYBhwBGI\naGC93zuO17ZYnzmmRPsK0tAgQbn29uJ51PnEZWEMGSJWz6pVuT+Gjz6Sm66xEbbcMtrPzCdpFkax\nTsFYGOZHGwdJGM36WbI2aS6pSguGOSeu62FEQAXDppRg9Abyb5+12CLg5GngAsS9tDlwPpAFTNd7\nFSIATcBi4E3EWgGxZPI/p8n6/E5kMpmvH4Vwzh3wY2XEZWFkMrCVFZVxWhnPPCPbCRPiz+NOimCY\nNT9MSXc3ymFhJMEl5WxDpcu9+7Uw4lo6t65OrOBSAz0jbH3ceqMI8CIY5j5NomB46Sf9UkowmoF8\n/e6HiEY+U4C3gVnAq0igvA1xP2UQC+NBREAakWD6fwX4HM8kZWazwbilFiyw9z34oGwPPzz6z8vH\nb+ZY3HW1iv0QyxHDSIJLqls3ebS1lXYVJsXCiNMV5HXCnPN4XNfDi0uqHCV9kkQpwfgQyXQa4di3\nKxJvyKcVOA9xQ41AMqVmWscagT2AW4BN1rE7gCOt47Ot9zX0Ara39ncim81+/ShGkMB3XGm1ACOt\niMzs2fb25ZelE5o4MfrPy6eaLIykuaTinCTmZb0U5/FKWxhxTyD0+n8pl4VRqGAolKekT1C89pN+\nKCUY64CHgasRy+AA4Bjgbpdzh1qPDLAPkiF1pXVsBZJK+xMklbY/cBp23OIRxF01EehuvW4WIliB\nSZqFMXasbN96S7a//KVsTzutPMs7JkUw/Jj65XJJlfpNVbrQXUdH/EvW+s2SqmQ6K8RvcZW6T9vb\n1cJwYzISqF4G3AOcA8xB5kmsRSwKEItgBuJeuh24GHjOOpZFxOAYRDzmAxuAC63jy4FJiFtrFTAO\nODH41xL8WhjOBVnisDD22EO2r70msYvHH5cf/+WXR/9ZbiRFMJwWRqkFpeJagRDEFeR19b84BcOL\nheEUi/r6wueFwatLKu41rKvFJbV8uQh5Y2P063EklYbSp7AaON5l/yJyg9+vANsWeZ83gG8WOf48\nsJOH9njGbwcZ54IsAGPGiOm6eLEds7j00vKNTpIiGN27iyC3tMhoNd+l0NEBS5bI8+HDO78+Svr0\nsTPXitXxqrSFUY4S2l5dUuWqIlxKuCod9DbuqKFD4/n8JJLa0iAQrtheHNTXw1ln2X9/85tw0UXx\nfJYbSREMKP5jXL5cxHvgwPgKDxq8rv5XaQsjbvcL+O+o4xaMSlsYffvK4HHtWncL1EzCTWL8Ii68\nWBhVi1+XVNzVL0Esiu7d5bMuukhuyHKRJMEYMEAsrVWrYJttco999pls456XAnZnUyrgXGkLoxy1\ntfr0kSyl5mZx0xW6N8tlYVQ66J3JiEv088/FRbpVXrGiRYtkG7cVnCRSLRhJ6iANXbrAJZfE9/7F\nSNL1KGZhLF4s23L8EDfbTLbFyj9AvCN8LxZGOSaI1dVJ7a7Vq+VRKOEgCYKRzdrH4xIMkHvw889l\nEJMvGB9/LNvtt4/v85NGql1SQS2MuFxSlSaJguGWsmjmqZTDwvAqGHF2Tkmqimrev1jKcxIEw0zs\nM/NY4sIMWozV60QFI2X4WdQeymNhVBI/he6c58VxPYzf160Y49y5st1xx+g/Nx+TsVVJC8OLW6zc\nglEs5TkJRRDLYV2ACkY+qRYMP+sDQ/xB70rjZ6IaxJtvP2yYbE02lJN582Q7alT0n5uPsTCKjagh\n3gCr6XiLdZDGEotzIiMkQzBMSXsTt3Ej7viFoZBgtLXZyxKMGEHNkGrB8DuiLkfQu5I4BcPL5M84\nLQzzQ8wXjGwWPvhAnpdTMLxaGHG6pKrNwogrO8lLJQDz/zL/v7gw96mJqxnmzhW32LbbxpuIkDRU\nMByk3cLo0kXiOu3t3uIYcQpGIQtj4UJYtkw6LlN7K068CEZbm1yvurp47g0vFkZSBKO9XSxPZ82n\nqPEiGOaYOTcujPVgrF6DqdZgqjfUCioYDtJuYYD3yVlQGcF47TXZ7rNP/NV7wZtLyukvj6NNRjCK\nuWDKVRXVdMCFBMNpadXF1Hv4EYy4r8dO1lTiefNyq+dOny7bvfeO9/OThgqGg7RbGOB9cpbTConj\nejgFw/lDfOop2X6zWE2ACPES9K50zSLnsUpnScUdv4BkWRi9e8s8oY0bZe0akN/G09ZKPuWoMp0k\nVDAcqIVh40wxjmMk2aOH5LVv2mRnm6xfL/W1AI53K0YTA15cUnELhglkF6uKWq4id6XEqxzClSTB\nANjVqqNtrN/HHpP/x4gRUu6nllDBcJD2tFpIxqpqhtHWuo2m3Pvtt4tbZs897VLwceNFMOJO4dxs\nM3F1rV7tvmhQa6t0kA0N8WdJmcl6y5a5Hy9HtlbfvjJIWbtWBhRulFMwDjlEtk8+KYJpyvmcd155\n3KZJQgXDQS24pJJSwhpswXjvPekcfvtb+fvnP4/vM/Mp5bOH+C2M+vririDnmgtxxQ0MppCeW7oz\nlEcw6upKC3k5BeOYY6RNDz0kQe5PPhGr4yc/if+zk0aqBcPPeslQGy4pr3MxymFh7LWXbF94Aa68\nUjqpPfcsz2JShgEDpDNYvbpwifNyTBIzHbBZOMqJ6bxN3CdOnLElt9Trcs0HKeWWKqdgbL01nHOO\nXI9FiyQQ/sQTknVYa6S6lpRaGJ3xa2HEKRjf/ra4WV5+WR51dXDrrfGt9+BGfb0UmPvyS1i61L0c\nSTkqxRaLY5jZ8OUoo92njzzWrhURze+QkyIYxiIsh2AA3HQTjB8vg4rvfjf+KspJJdUWhga9O+PX\nwojzWvTvD//xH/bfv/2tvchUOTFlSkxgOZ9yCIaJHbgJhrEwyrXuQrFZ+OUSDPP+hWIpS5fKNs4F\ntpzU18P3vw+nnFK7YgFqYeRQCxaG16B3OSwMgBtuENfU4MFw5JGlz48Dk3lUSDDM/IhKWRgmi2zb\nYsuTRciwYTKTeckS+MY3co+VSzCMOJr4jZNsNtlraacZFQwHtZAl5TWttlzXomtX+NGP4v2MUpQS\njHJMEjMWhhk5OzE1i8qVOVaoHAbYI/5yCYZbccqmJknB7tUr/lpSSi4145KqdO2kpJA0CyMJmFGq\n22gWyhNgNbETt076ww9lu8MO8X2+E1MOw3yuk3ItnVtMMNS6qBypFgyzNndHB2zYUPr8WhCMpFkY\nScBYGJUUDLM4j1nFzbBuneyrry+fS8qUw5gzJ3d/NmtXbY07Y0sFI5mkWjDAn1uqFjpJr0uS1pKF\nUcwFA5UVjLfekgHPLrvI4KccmHVIzLokhpUrZeDVr1/894UKRjJRwXBQC4LhtZx3LVwLg1kAxwSX\n8ym3YDjdp2+8IVszZ6UcjBghFs2nn+b+box1UY6lc4tNIPzkE9nmrwWvxI8XwRgATAOagQXASQXO\n6wZcDywBVgFTyQ2qNwNrHY824Cbr2DZAR97xSz1/iyJ4FYxstjY6SS91eqC2LIzttpPtp5/KaD6f\ncghG374ycm9pyU0lNUXuDjggvs/Op1s3mcnc0WELFkjpeSjP0rmDB0u24sqVnQc3RthraeGipOBF\nMKYCrcBg4GTgVmBnl/MuAcYCo4GR1vPLHMd7A32sx+bAeuCBvPfo6zhnitcvUQyvgrFxo1ShNHGP\ntOKlMirUhngaeveGIUPE3ZI/ou3oKN9iPc5SKSDC8dJLMtovd8rx+PGyfekle59xUZVjYau6Ovtz\n8teiMFVja2lp1KRQSjB6AROBy4EWYAbwKHCKy7lHAzcDa4AViPVwRoH3/S6wFHjVZ3t841UwaqWD\n7NdPfoxNTYULu0FtWRhQ2C311VdiffbrJ4OJODFVUWfNku3NN0sxwqOOKt+MZsO3viXbadNsF5kR\nDBMUjxu3WEo2a2dvqYVRfkp10CMR19FHjn3vIFaEG87ajXXAcMRayOc04C6X/QuBxcBtQCRZ7yoY\nuXgp7Aa1cz0MZo6DqZxrMLWd4l6HAmD33WU7Y4bMvbj2Wvn7F7+I/7PzOeww+c7vvQf/+pfse/tt\n2e7s5l+IASMY779v7/vsM7G8Bgwoj2tMyaWUYPQG8jP21+IuAk8DFwCNiMvpfCAL5M+b3hoYD9zp\n2LccGAdsBexhvf+9hRqVyWS+fpTCjJDNiLkQtdRBeolj1JqFYZbaNEtvGsqZkWMW43nkEZgwQVxk\np50G++8f/2fn07Ur/PjH8vzii2Vw8e67sn/cuPK0waxm98or9j4TUxk3rvZKi/vFTz/plVKC0YzE\nFZz0Q0QjnynA28AsxNU0DbFO8ueungK8glgThnXAW0jgexlwLnAo4hILhddie7UkGF7iGLV0PcAW\njH//O3e/mf1dDsHYcks49FB5vnixjLBvuCH+zy3EJZfI4OKll6TGVzYL++1XvlpK++8v8ZuZM+2S\nJE88IdsDDyxPG5RcSgnGh0imk9NbuCvwvsu5rcB5iBtqBJIpNdPlvFPJtS58ty+bzX79KIVZSrLU\nvINa6iC9WBjlKLiXJHbbTdx1s2fbRSih/Dn/t90GJ50Ep54qZd/79y/P57qx2WZw990ykv/0U9lX\nzjUgevcW11hHB/zudxLLuP9+OVbOEvjVip9+0iulBGMd8DBwNeJaOgA4Brjb5dyh1iMD7INkSF2Z\nd85+1jkP5u3fCxhltWcgEjCfjrsl4wu1MDrjxcIox9rNSaJXL4khtLXlukCMYMS9NKph2DD429/g\nzjuTMTHtyCNlpbljjpGYyve+V97Pv/hiEazrrpNg+/r1MGmSHd9QyouXrKTJQA/EVXQPcA4wB4k3\nrEUsCoDtkSyqZuB24GLgubz3OhV4CBEiJ9sBTyHxkveQlNtC8z18oRZGZ4qt7gbiejDXq1YsDLDd\nQc8+a+/TWcUSW3nsMVmatNxxg/HjZeneYcOge3c44QT4y1/K2wbFxkui4GrgeJf9i8gNfr8ClKp2\nc06B/X+3HpGjgtEZIxhuq7uBjOLa2+UHmuY5Kfkceihcc41MlrvuOtmns4orz2mnyUOpPKkvDeLV\nJVULiycZiq01ALXnjjLst5/47T/4wE6vNZPENOdfUWpAMNTC6Ixxr7gVdoPadEeBWFOTJsnz++6T\n1OKlS6VURjnW01aUpKOCYVFLglGsEijUroUBkqEEIhhmotqoUZJBpSi1Tup/Bpol1ZlSgmGuVS0K\nxoEHynyITz6BK66QfeWsFKsoSSb1gqEWRmcGDoQuXWT2bmtr5+O16pICmSh2wQXy/MUXZXvQQZVq\njaIkCxUMi1oSjLq64suS1rJLCuCss+zV7TbfHI49trLtUZSkkHrB6NFDqoy2tkoJ80LU2sxm45Yy\ni+I4qXXB6NNHymHceKNYGbVST0tRSpF6wchkvMUxak0wzAjazDNwUssxDMOWW8L555dn7QdFqRZS\nLxhg1+NZs6bwObUmGGZewUcfdT5myp5Xso6RoijJoyYEw0vtJBUMG1MZtLGxfO1RFCX51IRgmI7P\ndIRuqGDYqGAoiuJGTQjGoEGyLVQ7CWpXMObPt5fgNKhgKIriRk0IRikLo71d0mozmdpIqwUR0cZG\nyYhavDj3mAqGoihuqGAAa61VN/r0qZ0SEJmMvYa0KYFhUMFQFMWNmugeTcdXyCVVa+4og5tgtLTI\no2tXnX+gKEouNSEYJoZRyMJQwbD3mUyygQPLv1iOoijJpiYEo5RLqlYFY489ZPv663bg2xQkrOUV\n5hRFcaemBENdUrmMGAFDhsCyZfDhh7LPlArZcsvKtUtRlGRSE4JhRstLlnROIYXaFYxMRtZMBnj5\nZdkawRg+3P01iqLULjUhGP36SQZUS4t7eRCzrxZrJ6lgKIrilZoQDLBdLPlzDiA30FtrGMGYPl2s\nrwUL5G91SSmKkk/NCIYZMbsJxqpVsh0woHztSQpjxkip8yVLJFvqgw9k/047VbZdiqIkDy+CMQCY\nBjQDC4CTCpzXDbgeWAKsAqYCDY7jzcBax6MNuMlxfAIwF1gHvABs5fE7eMKMmN3Wf6hlC6OuDo45\nRp4/9JCUCslkYMcdK9suRVGShxfBmAq0AoOBk4FbgZ1dzrsEGAuMBkZazy9zHO8N9LEemwPrgQes\nY43AQ8ClwGbATOB+f1+lOMVcUrVsYQB85zuy/e1vYdMmGDkSevasbJsURUkepQSjFzARuBxoAWYA\njwKnuJx7NHAzsAZYgVgPZxR43+8CS4FXrb8nAu8jorERuArYFRGeSNhmG9l+/HHnY7VsYQAccggM\nG2b/PWFC5dqiKEpyKSUYIxHXkbMI9juIFeGGc25wHTAcsSjyOQ24y/H3aOt9DS3WZ44p0T7PGBfL\nnDmdj9W6hdG1K1x9tf188uTKtkdRlGTSUOJ4byB/YdO1uIvA08AFwHTrfc8HskBP6zWGrYHxwOmO\nfb2A/Gl1TdbndyLjqFmRdZtY4YIRjHnzoKMjt8hgrVsYAGecIaVC+vWD7bardGsURQlLJobaPqUs\njGYgfzpbP3IFwDAFeBuYhbiapiHWydK8804BXgEWBvycQPTrJxP4WlthoeOTOzpswahVC8Ow++4q\nFoqiFKaUYHyIWAsjHPt2ReIN+bQC5yFuqBFIptRMl/NOBe7M2zfbel9DL2B7a38nstns1w8/7GyF\n6t93tH7FCmhrE7Ho3t3X2ymKoiSWoP1kMUoJxjrgYeBqxLV0AHAMcLfLuUOtRwbYB8mQujLvnP2s\ncx7M2z8NiVdMBLpbr5uFCFZkjBsn2zfesPeZYntDh0b5SYqiKOnDS1rtZKAHsAy4BzgHmIPMk1iL\nWBQgFsEMxL10O3Ax8Fzee52KZEKty9u/ApiEuLVWAeOAE/19ldLss49sX3/d3qeCoSiK4o1SQW+A\n1cDxLvsXkRv8fgXYtsR7nVPk2PNArPOLjWC88YYsy1pfr4KhKIrilZopDQKw+eYS1G1uhplWdMXM\n/Nb1HxRFUYpTU4IBcMQRsn3ySdl+ZM0wGTHC/XxFURRFqDnBOPJI2T7xhGzNwkEjI5tTriiKkk6q\nbTIX+GAAAAX/SURBVNXmLHifrOdGS4uswLd+vVgXY8fKAkpLl8LgwZG1U1EUJTE4JvGF6vNrzsLo\n2RMmTZLnF18sYjF8OAwaVNl2KYqiJJ2aEwyQMhgg5bwBDjhASnoriqIohalJwTjoIBEJwylutXcV\nRVGUHKptXB06hmH47DOYMgV22AEuvFAtDEVR0ktUMYxq6yYjEwxFUZRaQYPeiqIoSllRwVAURVE8\noYKhKIqieEIFQ1EURfGECoaiKIriCRUMRVEUxRMqGIqiKIonVDAURVEUT6hgKIqiKJ5QwVAURVE8\noYKhKIqieEIFQ1EURfGEF8EYAEwDmoEFwEkFzusGXA8sAVYBU4GGvHNOBOZY7/URYIqMbwN0AGsd\nj0u9fQVFURSlHOR36G5MBVqBwcDuwBPAO8AHeeddAowFRlvv+zhwGXCVdfzbwO+AE4B/AVvQuXJi\nX6yKtIqiKEqyKFXqthdiLYxGLAKAO4HPgV/mnfsm8F/A/7P+Psn6eyvr738CfwZud/mcbYBPgC5A\ne5H2aHlzC1OuWK+FXgsnei1s9FrYlKu8+UigDVssQKyL0YXalffew4E+QD2wB2KlzAcWAzcD3fNe\nv9A6dhswsHTzFUVRlHJRSjB6A015+9YiIpDP08AFQCOwOXA+YhH0BIYg1sMkJG6xG+Leusx67XJg\nHGKN7GG9/73+voqiKIpSSXYH1uXtuwh4zOXc7ojV8BlikVwCbLCObYYEtZ2rZ08E3irwuUOs83vl\n7c/qQx/60Ic+Aj9CUcrC+BAJYI9w7NsVeN/l3FbgPMQNNQKJfcy0jq1GhCRs+6ptSVlFUZSkUJb+\n8z7gb4hr6QBgDbCTy3lDrUcG2AdYBBziOP4rJDtqEGJxvGLtA9gLGIUIxEDgfuD5iL+HoiiKEjOb\nkTsP40Rr/1ZIPGO49fc3gU8RF9YcOs/XaEBSdFcDXwA3AF2tYyciWVLNSAbWHUiAXFEURVEURVEU\nRVEURVEqiNcSJdXOuUiyQCudJzlOAOYibr8XsCdFGv4LWGE9fhdvM8tCV+CvyP+7CXgbONxxvNau\nxz2IO7cJceE6y+fU2rUw7ID8Vu527Ku1a/EisB67rNIcx7FauxZfc5/16AnsjwTfd65oi+LheOBY\n4PfkCkYj8p0nIR3pfwOvOY6fjdwYJvlgtrWvmukJXIl9kx+FdJZbIdfjK2rreozGnuw6CvgSOIza\nvBaGZ4GXgbusv2vxdzIdOMNlfy1eC0DmY2wgN733TuCayjSnLPyaXME4C3jV8XdPoAWZjQ9SeuVM\nx/HTyb050sI7yByeWr8eo5CqCGOp3WtxIpJReSW2hVGL12I68GOX/ZFfi2opb+63REkayM+ZHo18\nZ0MLcj3MNdg57/i7pO/6DEHuhfep3evxe8S9MBuYgkx+rcVr0RdJy7+Q3N9KLV4LkMHzckQgDrT2\nRX4tqkUw/JQoSQv5szJ70fkaNGFfg96IW8J5rHc8TasIXZByMXcgE0pr9XpMRr7HIcBvkDlMtXgt\nfg38BUnDd85irsVrcTGwLeJW+hNSKXw7YrgWXsqbJ4FmZEThpB8iGmkl38IodQ3yj/ez9qWBOsTl\n0IokBUBtX48sEuh8EEn+qLVrsRsSzN3d+juD/XuptWsBMiHacBdyTxxJDNeiWiwMPyVK0kK+hTEb\n+c6GXsD21n5zfDfH8bRcnwySKTUICd6Z8ve1ej2cdMF2T9XStTgQWRJhEZI19jPk3vg3tXctilHT\n18JriZJqpx7JhLkGGS10s/aZjIeJ1vH/RoJWhrORRa2GAsOQm+GssrU6Pv6ABOLyC1HW2vUYhAR5\neyH3w2GIO2FPau9a9EAqQQxG4lrXItbWQGrvWvRD7oXuyKD6ZMRKGEHtXYscCpUoSRtXIZV6nY8r\nrGMTkBzrFgrnVK+0HmnIqd4a+f4t5C7fa+bg1NL1aETcUKuRTuBfwHccx2vpWuRzJXZaLdTWtWhE\n7oUm5N74J/L9DbV0LRRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRF\nUVLI/wdu3uuHqHAWjgAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }