{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Levenberg-Marquardt Fitting" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Author: Maurice Wilson\n", "\n", "Objective:\n", "Use Levenber-Marquardt Statistics to determine the best 4 fitting parameters for simulated data of a normal distribution." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline \n", "from scipy.optimize import leastsq" ] }, { "cell_type": "code", "execution_count": 19, "metadata": { "collapsed": true }, "outputs": [], "source": [ "x_vals = np.linspace(0,10, 300)\n", "\n", "A = 10.\n", "B=5.\n", "C=2.\n", "E=4.\n", "\n", "A_g = 20\n", "B_g = 7\n", "C_g = 5.\n", "E_g = 9\n", "\n", "sig = 0.25\n", "\n", "sigmas = np.ones(len(x_vals))*sig" ] }, { "cell_type": "code", "execution_count": 20, "metadata": { "collapsed": false }, "outputs": [], "source": [ "noise = sigmas*np.random.randn(len(x_vals)) \n", "Data = A*np.exp(-(1.0/2.0)*((B-x_vals)/C)**2)+E + noise" ] }, { "cell_type": "code", "execution_count": 21, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 21, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAl8AAAG2CAYAAABBHOPKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3X+c3FV96P/XGcLqoqiEvRIpNiD+SAVcN+G2WOzd9N7u\nbq/15jak39YCbcTWXO6lKHGgaZqCKTE3pbLCldpyQwukcoOtjdj1XpvJljapwVZrSGMMlwsi5ouX\ngiyhtdbYFefcPz4zuzOzM7szuzu/X8/HYx6Z+cznM3N2svPZ9+ec93mfEGNEkiRJjZFqdgMkSZK6\nicGXJElSAxl8SZIkNZDBlyRJUgMZfEmSJDWQwZckSVID1S34CiHcHUJ4NoRwtGT7tSGE/x1C+EoI\n4ZZ6vb8kSVIrqmfP1z3ATxZuCCH8OLAGeEuM8ULg1jq+vyRJUsupW/AVY/wc8ELJ5v8M7Igxfi+3\nz3P1en9JkqRW1OicrzcA/yaE8DchhP0hhIsb/P6SJElNtaQJ73dGjPGSEMK/Bv4YeF3pTiEE1zyS\nJEltI8YYqt230T1f3wA+BRBj/FsgG0I4s9yOMUZv87x98IMfbHob2vnm5+dn5+fXnjc/Pz+7Zt1q\n1ejg69PAvwUIIbwR6IkxPt/gNkiSJDVN3YYdQwj3A4PAmSGEp4CbgLuBu3PlJyaBX6zX+0uSJLWi\nugVfMcafr/DUL9TrPZVYvXp1s5vQ1vz85s/PbmH8/BbGz2/+/OwaK8xnrLLeQgixFdslSZJUKoRA\nbOGEe0mSpK5m8CVJktRABl+SJEkNZPAlSZLUQAZfkiRJDWTwJUmS1EAGX5IkSQ1k8CVJktRABl+S\nJEkNZPAlSZLUQAZfkiRJDWTwJUmS1EAGX5IkSQ1k8CVJktRABl+SJEkNZPAlSZLUQAZfkiRJDWTw\nJUmS1EAGX5IkSQ1k8CVJktRABl+SJEkNZPAlSZLUQAZfkiRJDWTwJUmS1EAGX5IkSQ1k8CVJktRA\nBl+SJEkNZPAlSZLUQAZfkiRJDWTwJUmS1EAGX5IkSQ1k8CVJktRABl+SJEkNZPAlSZLUQAZfkiRJ\nDWTwJUmS1EAGX5IkSQ1k8CWp7WQyGYaH1zE8vI5MJtPs5khSTUKMsdltmCGEEFuxXZKaL5PJsHbt\nek6evAWA3t5NPPDALkZGRprcMkndKoRAjDFUvX8rBjkGX5IqGR5ex/j4GmB9bssuhobG2LdvTzOb\nJamL1Rp8OewoqeUVDjNOTDy/oOMdppTUbPZ8SWpppcOMPT3XAacyOflhYO5hx+njrwQeIpV6nJtv\n3siWLVuK9hkd3QlAOr3BIUxJNXHYUVJHKTfMODBwF319ZwFzB0vJ8ecB9wFJAJdKbeSzn72fkZER\nc8gkLVitwdeSejZGkuqhr++sGnO8HiIJvJIALpuF0dGdjIyMMDq6Mxd4rQcynDx5Hpdffg27d3/M\nAExSXZjzJamlpdMb6O3dBOwCdtHbu4l0esOcx+XzvCYmniWExyruc+jQkfwjkgDsak6cuJG1a9eb\nHyapLhx2lNTy5srJKn0eKBpKXLLkGr7//VOJ8XYgGVrcsuVatm+/I5cLtgt4PXA1zqKUVKuWGXYM\nIdwN/BTwzRjjRSXPpYEPA30xxhP1aoOkzjAyMlIx4JqYeJZjxx6bSsA/eHA9K1asyAVey4CdvPji\nRZx//rd43evGAEind5UMNw4Bv9zIH0lSF6vnsOM9wE+WbgwhvJbkTHe8ju8tqc3UUg4inyQ/Pn4e\nhw8fzwVe64Gkt+v48W8AR3Pb1gBX8+STf086vYF9+/aUyeUaAX6OVGojtQ5vSlKt6hZ8xRg/B7xQ\n5qmPAL9ar/eV1H6mg6k1jI+vmTPfKum1upJkBuM5M54/44yXAHcznWS/nmz2tqmhSYDBwZUlwdZ9\n3HxzmqGhMYaGxuYsX1EuULSemKRqNHS2YwjhPwLfiDF+OYSqh0YldbjiIUA4efIol19+DatW9c9S\nSiI/g3EZ03laSR2wp546FfiBiu+XyWTYvv0Ostn3AHeSSj3Oz/7sOzlw4GFg9vIVpaUpDh5czwMP\n7AIou90Zk5JKNSz4CiGcBvw6yZDj1OZK+2/dunXq/urVq1m9enW9miappWSAXZw4cSvj48VBTGGe\nVwiPkczLGSHpvdrK0qXPsXx5P4cPX0VpUJYMIyZBUmmwl81ez8c/fjfZ7G3A7IHTzECRqR61ctsN\nvqTOs3//fvbv3z/v4xvZ83U+cC5wJNfrdQ5wKITwwzHGb5buXBh8Seps6fQGDh5cz8mTAHcCt1Iu\nuJk5g/E68hOje3ufZPfuXQVDi8VB2e7ds/VCPZQLvGYPnKZLU6zJbwHu5NCh51i+fObwp6TOVNop\n9Ju/+Zs1Hd+w4CvGeBQ4K/84hPAksMrZjpJGRkZ44IEkcDp06DlOlDkrlPY4vfgiuUr30zMYASYm\nniWV2kg2mxyXD8oKA6niYA9Sqcen9q+keJmi60kS+ncBt3LiBHz729fR03MDk5P5953uaZOkQvUs\nNXE/MAicGUJ4CrgpxnhPwS4W8pI0JV9OYjrISbbng5jCZPm8wkr3xblYR0ml0vT3X8iOHTN7vAqD\nPYDBwY1s375pxnsWKl+aYrqHbnKyfDA4PLwu99g1IyUl6hZ8xRh/fo7nX1ev95bUnvI5XStWvB64\nh76+M0mnp4Onwt6q0gBp8+YdJXlcF9HXN1Yx4CmtHXbxxRcXBGPXMjq6k9HRnRWCphHKzbKsHAya\ngC9pmms7Smqawsr0g4MrcxXnyy9wXdpbVRiUZTIZjhz5Sk3vVxpUzex5Kw6aSocqe3oeBSoPM1ZK\nzDf4kmTwJakpSoOcBx9Mk82OMluwUq7SPST7ZbPvBjZNbUulNpJO31/x/Sr1RJULmjZv3kZf31lF\nPXKDg9ezZ8+fcfz4NpYvP6fs8KYklWPwJakpZpZ7uHOBr3gRSQL8TuBp+vvfXBQMzb8n6ihHjjxC\nNvteoHRdyFtyr7VpxlGlPWUm4EvKM/iS1CIuLZmlWH2wMh3o3AKsobd3Ezt2zC/QmX6toySFXB+b\nUYbiIx/ZNmtRWGDW3DVJ3S3E2HqTDkMIsRXbJWnxlA4D5nuU8lXmBwdX5ob1nskN622eNXiZLZ+r\n0vtVSoDfvn07N900mgu67gSuZrpg6y6WLt3GiRM35rZlgCtJZj4mFfbh1KmFvmd7H0mdIYRAjLHq\npXsMviQ1TaWAKZPJsGbNu5icXMJ0UHMDY2MfX1AQM1eAljc8vI7x8TWUC65mDjuWBmdvozRYGxoa\nm5oFKanzGHxJantJ8PM09Qpi5grCpoOvZSQ5ZI/Q2ztBb+8rp3rhgFxR2CO5XrD8vl8EPlSXdktq\nTbUGX6l6NkaSKslkMgwPr2N4eB2ZTKaKI45y6NCRGvav/L5r165nfHwN4+NrWLt2/YzXS6c35IYP\nryRZSuinOHlykhMnbuTw4atYuzYJrPbt28Pu3R8r2ffngPeRJP/vyuWubZh3eyV1oBhjy92SZknq\nVHv37o29vWdFuDfCvbG396y4d+/eoud7el4VoS+3TzrCKyruX4uhoctyrxNzt3vj0NBlM/YbGBgs\n2G/2Y4r3jRHScenS8+PQ0GXzbqek9pGLW6qOc+z5klQ3lXq3iss+JOslXn75NVP7jYyMMDb2CQYG\n3sTSpds4/fRPAh8t2P+WsssNLaa+vjMXsO9FrFrVz759e0y0lzSDpSYk1UW5oqb52YyHDh0hGaKD\nJKF9FydO3Mr4eHHx03zgkuRgLU67qq2/VbzfeSRDieWPmU9Nr2qT/yV1HhPuJdVF8YxBgOtJpe7O\nlW84CtxF0ps1s5RDaYJ6LWUiqlFt4FO6/FG+DEalUhbVBlOL/fNIai5nO0pqCTODr9ISDNezdOmn\nAQpqZkGl2YGd0FOU/xmmZ0g6I1LqBLUGXw47SqqL0qG4VOrxqer1iYtYtepJ0ukNuV6gZGulIbtK\n6zq2qtJgESjo7Xq6ZO/pmZztGlhKqp49X5LqpnTYbrow6VFSqXvp77+wqGYWtG+vVqHSYcWenut4\nyUtewT/9083MLNxaOATrEKTUjhx2lNSyMpkMmzdvyy1UfRvQmcFG+Qr5r2chw66SWpdFViU1zVyF\nU0dGRujrO6tgoerGlI1orp0kPVxbgU1MF1+9jw984KpmNkxSk5jzJWlRlCst0Wk9WtUqzncrzO96\nPfAbnH56YNOm/PqQVwLXT+2Rz3nrhAkGkspz2FHSopg5u7HyrMVuKLOQD54mJp7l6NGv8OKLp1K4\nSPgFF7yZw4evYnpocitLlz7H7t0fA+iKz0jqFM52lNTSRkZGeOCBXQW9Op0ZVBTOzly5cnVBoAWT\nk3D8+LbCvYFnWLVqjJGREYaH1xWsAAAnTyYTEjrxc5K6kcGXpEVRS5X3disbsVDllipavnwZJ09u\nqqkqvqTO4LCjpEVjnlJ5lYZaoXyJjW4ZmpU6haUmJLWkbg/Mav35u/3zktqJwZeklmNPjqROZp0v\nSS1ndHRnQQJ5N9T2mt1c9dAkdTaDL0k1Kxc8GFBUJ98LOD6+hvHxNaxdu97PS+oyDjtKqknxEGKy\nRuN5553NU089w+Tkh4FkLcMLLuinr+/MMotKd/ewY7X10CS1D4cdJdXV9BDiMuA+stlRnnjiZbnA\naz2wjMnJJRw+fNVUzw7AAw8kQcbQ0FjXBl6zma3n0F5FqbNY50tSVfKz7w4dOgKsIVmzMJ/HNVaw\nZ34tw+ICofv27THgonw9tMHBaysuzeSyTVLnMfiSNKfiAOA84H3Amwv22ABcmbv/dOnhKlCuwn/x\nhITiivazPSepPRl8SZpTaQAA8PKX7+Y739lINps87ul5kQsuuAc4hWPHbmByMtme79kZHl4HWLMK\nZlb4rzTzM5PJFPQ0TpuYeNbPU2pjBl+S5uEi3va2t5FObyjowflEUYX2/PbBwWvZvv0Oh81mMT0U\neRR4iFTqcc4++5253sYrgeun9u3puY5jx05lcvK9gJ+n1I6c7ShpTgspkursvups376dm24aJZu9\nDYBUKk02O0ryuWWArSxd+hzLl59TtEi3n6fUfLXOdrTnS1JFhT1YW7Zcy4EDSWJ9Oj134DUzQV+z\nOXDg4VzglQRb2WxvwbMjwDOsWjVW/mBJbcXgS1JZM2fZbcoFYA9PBWTlhhln1vXKJ+gnens3kU7v\natwP0nYyJAFY8XBj4edWOlvSz1NqLw47Spohk8lw+eXXcOLEjUwPb11PKnX31LBYfugRZhZQXbFi\nRcnQ2PUsXfppVq3qN0G8gulg9zzgakqHG3fv/ljFYNfPU2ouhx0lLUhxEFDooYJhsemSB8n94lII\nx49vKzn2IlatetK8pFnkS1AkQe/UVvLDjYUBVulsSUntxeBLUpHiCvbTpSVSqcenykrMZfnyZZw8\nucmhsRqNjIywe/fHcsFvss3PTuo8Bl+SKhgBdpEf9vrABzayfft0QNXTcx0TE/3Ai/T0FNf12rEj\nCRYKC4naU1OdckVY/eykzmLOl6Qis5WVyOcaTUw8y7Fjj1VcSNtgQVI3qTXny+BL0gxzJXRbu6s5\nTLSXWpMJ95IWzITu1lNahNXK9lL7sudLUs0WUvFetctkMrzjHVcUVLwHy3dIrcNhR0kN4RBY4yTD\nvE9TXP/rSuBWwJw7qdkMviSpwyTB13nAfcAtwJ1UCsTshZQar9bgK1XPxkiSFi6d3kBv730kQdad\nwGMFz+4kCbzWA8lQcL5HUlJrMviSRCaTYXh4HcPD68hkMs1ujkrka38NDT3J0NDZfOhDH6C3dxNJ\nHbanm908STWq67BjCOFu4KeAb8YYL8pt+zDwTmASeAK4Ksb4jyXHOewoNUg1yfPmd7WeTCbD5s3b\nePzxr/HP/3ySGG8HHHaUmqHVSk3cA9wB/GHBtn3AphhjNoTwW8Bm4Nfq3A5JFUwvJzRzzcZyBVUt\ncdA6Hn30q7n/u6OkUmn6+y9kxw7/b6RWV9dhxxjj54AXSraNxxjzK8R9ATinnm2QlKhlaHFi4lnW\nrl3P+PgaDh/+fi7wMqeolRQHzbeSzY7S13emgZfUBppdZPU9wP1NboPU8UqHFgt7r9LpDRw8WLyQ\nM6wo+MM+1qxmS1JHalrwFULYAkzGGHeXe37r1q1T91evXs3q1asb0zCpA1UaWsxXsi9dyLm4Z2sD\nySy7RG/vJtLpXY1rvMoqFzT7/yI1xv79+9m/f/+8j697na8QwrnAZ/IJ97lt7wbeC/y7GON3yxxj\nwr20iKpdi9GFs9tLpYkQTpCQGqvliqyWBl8hhJ8ERoHBGONEhWMMvqRFVO2MxsJ9DLjaU/H/41FS\nqXtzifib/T+U6qSlgq8Qwv3AINAHPAt8kGR2Yw9wIrfbX8cY/0vJcQZf0iKbqzek2t4xtbbp/8dl\nJP+Xrr8p1VtLlZqIMf58mc131/M9JZWXz+9St9hJEnjNzPOT1FxWuJcE5JewyVdN35VL4N7Q7Gap\nRtP/j1a+l1qVC2tLmmKidmfIV78/cuQRstnbAIcdpXpqqZyv+TL4kqSFM5iWGsPgS5IkqYFqDb7M\n+ZKkDlHLElKSmseeL0nqANXUcpNUH/Z8SZpTvodk5cq3s3LlantKOsDmzTsKlpAqXgDdHjGptTR7\nYW1JDTbdQ3Il8FfArQAcOPAuK9q3qUwmw5EjXynZepRDh46wcuXbi5aKKlxUXVJzOOwodZnpCuhj\nQL6ifYZk8ewkEHPIqr0k/6fnAfeRFFY9CtwFfBS4E7gaVy6Q6sdhR0lTqh9u2kkSeM0cslK7uIik\nQO4Y8Eckgdd64OyKRzgcKTWHw45ShypNwM4PN6XTGzh4MD/seH1ub6uht7Pp/9NbgDWkUgfIZvPP\nbiDp1TwKPEQq9TiDgxsr/n7Y2ynVn8OOUoeabaHsfAX0xx//GiGcyqtf/UqeeuqZqbwghx3bT2FB\n1cHBlWzffsdUYLVkyTV8//unEuPtQPL/u2LFCg4fvgqHI6WFa6mFtSW1rkcf/erUH+cXX9zETTe9\nnwMHxgBIpw282k3pwukXX3zxVDA2MXFxUaB18iQcP76tGc2UhMGX1LGmh6IAjpJK3cvXvnY2K1eu\n5vjxbxSUJUj+GB84YK9HJykMxoaH1+W2Zkjy+57mjDNewsmTm3K/H+QWUt/VjKZKXceEe6lDjYyM\n8MADuxgYuItU6m6y2XfzxBN/z+HDV3HixL9qdvPUQOn0Bnp6riPJ/VoDXM1TTz3Hli3XMjQ0xtDQ\nmMPMUgPZ8yV1oML8H1hCNnsbySy4/IzGZSR/iBP2enS2kZERLrigv2jocXLS3k6pWQy+pA5TOost\nlUqX2WsEWM/SpdtYtarfHK8u0Nd3ZsGjDHAnhw49RyaT8f9eajBnO0odZnqW4zKS/J5HCOEZYvwl\nkjpQFlLtRsUrG/h7IC2mWmc7GnxJHWZmtXMI4Vd461tXAS8CS1xCqEtlMhkuv/waTpy4EUtMSIvH\nUhNSl0unN/Dgg1eQzY6S/wMbI/T1+Qe2242MjLBqVT/j41A483Fi4pTmNkzqMs52lDrMyMgI/f0X\nNrsZalHlZj4eO/aYywtJDeSwo9SBSpPuzetRoZUrV1vdXlpEDjtKmqrxlS834WxGFSqe+Sip0ez5\nkqQuU65ndMuWazlw4GEAJ2NINXK2oyRpTrMtxO0wtVQbgy9JUk2ma8OZAybNR63Bl7MdpQ6SyWQY\nHl7H8PA6Z69JUosy+JI6RD6PZ3x8DePja1i7dr0BmKqSTm+gt3cTSeX7XfT0XMfExPMG8VKdONtR\n6hCjoztzOTvJ0NHJk8k283ZUjRUrXs/x49s444zTeOqpU3OlKODgwfXmf0mLzOBLkrpY6czHf/iH\ndNHqCAbx0uIz+JI6RDq9gYMH13PyZPK4t3cT6fSu5jZKLa+0xzSbvbO5DZK6gMGX1CEsrKrFcSmp\n1Eay2eSRQby0+Cw1IUldzIKr0sJZ50vqQoUFM/1jqVr5+yMtjMGX1GVcRFuSmsvgS+oyVieXpOay\nwr0kSVILM/iS2lxpdfJkdtqGZjdLbcjlqaTGsNSE1AbKJURnMhk2b97G8ePPcPbZrwZGeeGF77B8\n+YrmNlZtqTR38ODB9c56lOrEnC+pxU3/UbwSeIhU6nF+4Rfeyf33/ymTk0uAW4GjwF3ARwGT7lW7\nmbmD15NK3U02exvg75Q0G3O+pA6TVCC/ErgPuJRs9g3s2jXG5OQKksBrPfAkSeC1Hkh6L/I9ZdL8\nPJQLvPydkhabw45SW3gIyAdgtwAuAaPFVbo8VSr1+FSVe0mLy+BLanHp9AYefPCK3B/C/Bp8y4Cf\nBa7P7XUe8L6pY1wSRrUqXZ5qcHAj27dvygVjR0ml7mVi4kIymYxDj9ICmfMltYHt27fzG78xCuSH\ngQCu5+Uv301Pz2ksX34O69YNmRytRZWf1HHkyCO5IcgkCOvvv5AdOzb7OyblWGRV6lDbt2/npptG\nTYBWQ00n4i8jCfxdSUEqVWvw5bCj1Ca2bNnCxRdfXFBywj98aqSdTA97w8mTyWQQfwel2s052zGE\n8PYy2y6tT3MkzWZkZIR9+/aQTm9gdHSnxTBVd9NFfJ9udlOkjlFNqYk7ymz7nbkOCiHcHUJ4NoRw\ntGDb0hDCeAjhsRDCvhDCq2pprKTpul/j42sYH1/D2rXrDcBUN/lE/IGBU0ilNuJKCtLCVcz5CiG8\nDfhRYCPwESA/lnk6sDbG2D/rC4fwY8C3gT+MMV6U2/bbwESM8bdDCJuAM2KMv1bmWHO+JMpXtnch\nbTVLud9HSYtbZLWHJNA6Jffvy3O3bwE/M9cLxxg/B7xQsnkNyWUTuX9/utqGSt3GHi61Goe9pcVR\nMeE+xngAOBBCuDfG+PVFer+zYozP5u4/C5y1SK8rdZyksv3MBOfSYpjW9FIjlVsD0lmPUm2qme34\nkhDCXcC5BfvHGOO/XcgbxxhjCKHi2OLWrVun7q9evZrVq1cv5O2ktjMx8XzZ7aXFMJ31qHqoNMRY\n6aLA30F1k/3797N///55Hz9nna8QwpeB3wMeBr6f2xxjjIfmfPEQzgU+U5Dz9SiwOsb4TAjhNcBf\nxhhXlDnOnC91tUwmw5o17ypYOBt6em5gbOzj/pFT3ZX2bhXW9DLnUJqpHnW+vhdj/L0FtKnQGNNV\n+tYDn16k15U6yujoTiYnbycpbLkTeJoLLnijgZcaYrberXLD3oOD1zI8vA4wEV+qRjXB12dCCNcA\nnwL+Jb8xxnhitoNCCPcDg0BfCOEp4Cbgt4A/DiH8EvB1ksXpJFU0krvtoq9vrNmNkcqsAXkt27ff\nYQ6YVINqhh2/DszYKcZ4Xp3a5LCjut5swz5SvdXy++cwpFSHYccY47kLapGkmplUr2by90+qr2p6\nvl4GfAD4wRjje0MIbwDeFGP8n3VrlD1f6lIWsVSrqvS7aS+tVHvPVzXB1x8Dh4BfjDFekAvGPj9X\nhfuFMPhSpyv3h8w/YmpVc/1uetGgbleP4OtQjHFVCOFwjHEgt+2IwZc0P5X+kI2O7jR3Ri3JvC5p\ndou5vFDev4QQegve4HwKZj1Kqs30NP5lwBgnT57H5s3bmt0sSVKDVBN8bQX2AueEEHYDfwFsqmej\npM53lKQXYQ1wNUeOPMLg4Ep6ezeRLHt6PalUmomJ5107T02XTm8o+N3clVvSakOzmyW1rTmHHQFC\nCH3AJbmHfxNjnKhroxx2VAfLZDK84x1XkM2OUjqMk05vYPPmbRw58gjZ7G2AuV9qDaV5XYB5XlJO\nPYYdAX4AOAV4CfBvQgiXzadxkpJp/P39F1Z8rq/vrFzgtR5IcsPyf+SkZhkZGWHfvj1TeV5r165n\nfHwN4+NrWLt2vT20Ug3mrPMVQrgHuAg4BmQLnvpUvRoldbodOzbnku6Tx8kwzq7mNkqqkotrSwtT\nzfJCPwJc4DigtHjKFbGEZFbZxMSz9PTcwORksq+BmdqVJSik8qopNbEL+O0Y47HGNMmcL3Wf0vIT\nPT3XccEF/fT1nekfLbWcamrSWbdO3aQedb5WA2PAM0yXmIgxxrfMt5FzNsrgS13GOkpqN3P1avk7\nrW6y6Gs7An8AXAl8heKcL0lSl8oHW6OjO6eCMHu1pOpUE3x9M8Y4VveWSF0snd7AwYMm4Kt9lA4r\nHjy4vmhY0d9pqbJqhh1/F3gV8BkglwJMjDHWbbajw47qRiYnq51UM6zo77S6RT2GHU8jCbqGS7Zb\nakJaRCMjI/5xUkcw6JJmV1WF+0az50uSWlul2YyAsxzVdRa9wn0I4U0hhAdDCMdyj98SQviNhTRS\nUrFMJsPw8DqGh9dZKVxtIV+rbmhojKGhsakAq7gAqys0SOVUM+x4F3ADcGfu8VHgfuBD9WqU1E3m\nSlyWWpVD5dL8VJXzFWP8QghJb1qMMYYQvlffZkndw6Va1Emc5SjNrZrg67kQwuvzD0IIPwP8ff2a\nJElqN5lMhs2bt3H8+DOcffargVFeeOE7LF++otlNk1pONaUmzgd2Am8D/gF4Ergixvj1ujXKhHt1\nEZdhUbvLZDKsWfMuJieXALeSZKfcBXwU8HdanW9RlxcKIZwC3BJjvD6E8HIgFWP81iK0c/ZGGXyp\nw1Saep/fPjHxLLDEtRzVlpKaX08DV5MMn68DXFpI3WNR63zFGL8fQnh7SKKhby+8eVL3me7ZuhJ4\niAcfvIKbb97IxRdfbI+XJHWhaoYd7wTOBj4JfCe32Qr3UpWSXoHzgPuAJNBKpTbS3/8WDh++CnsH\n1O6qHXYELL6qjlSPCvcvBU4A/7ZkuxXupao9RBJ4JYFWNgvHj29raoukxTIyMsLY2CdyCffbWL78\nHNat+1UOHEiWBU6nd/GlL32Jm24aJZu9DbCkirqbFe6lOstkMrzjHVeQzY4Cy0jmrzzN+ed/i6ef\nft5hR3W84u+APb3qPIve8xVCuKdkUwSIMb6nxrZJXWlkZISbb97IjTf+CjG+lGRYBp566gZuuun9\nRb0DBl7qRKOjO8lm39DsZkgto5phx/9FLuACeoG1wNN1a5HUgbZs2cKePeNFOV6Tk3DggFf+6nwT\nE88DlwLPDCnIAAAgAElEQVSbpralUhtJp+9vWpukZpoz+Iox/knh4xDCbpIEFklzKCwxIXWjTCbD\nsWNHgGMkFx53EsJj3Hxz2p5eda1qer5KvRH4V4vdEKnTlBZP7em5jp6eG5icTJ532RV1g9HRnUxO\n3s50viO89a0/xJYtW5raLqmZqsn5+jbTw44ReJbCvmNJZZWu2Tg5CQMDd9HXZ46XutFI7rZr6juQ\nV6kIsdSpqhl2fHkjGiJ1g76+s8zxUleZa6Ht0h5iS1CoG1RTZPVS4EiM8dshhF8ABoD/FmM8XrdG\nWWpCHcA1G6XEbD1bSRFilyJSe6tHkdU7gf4QQj/wAeAPgD8EBufXRKk7jIyM8MADuwr+6Bh4qTuN\njIz4uy8VqCb4ejHGmA0h/DTwsRjj74cQrPElVcE/OtJMhT1hg4MrOXhwU8VhSakTVTPs+FfAXuAq\n4MeA54C/izFeVLdGOeyoFmeCsDQ/xQvN/y/gm7zmNa9i2bLX0td3pt8ntaVahx2rCb5eA1wOfDHG\n+LkQwg8CPx5jrNulicGXWpm5XNL8TS80/wckgy/Jig9Llryfiy56qwGY2tKiB1/NYPClVmaCsDR/\nyfcnv0jK1STfowxwJflAzAsatZtag69UFS/4thDC34YQvh1C+F4IIRtC+NbCmil1tkwmw/DwOoaH\n15HJZJrdHKllpNMbSKUeL9m6kyTwWg8kvcquDKFOVk3C/e8A7wL+GLgY+EXgTfVslNTKytUtGhy8\nluHhdUCSQLx9+x3WLZLKmF5o/reI8frcVpcLVnepJufrUIxxVQjhyzHGt+S2/V2M8a11a5TDjmpx\nmUyGzZu3cfz4M5xxxmk89dQzTE7+Ismyp48BH8FhSamy8t+hDwMOO6r91KPO1z+HEF4CHAkh/Dbw\nDFD1G0id6tFHv8rJk7dw4sSdJB3C9wG3kJTGkzSb0jIsxTOIDbzU2arp+TqXZD3HHmAj8Argd2OM\nX61bo+z5UoubTrpfBlxDstb81bnH24BjwO2AV/FSrSzlonaz6D1fMcavhxBOA5bFGLcupHFSZzlK\nssb8lcA9BY9vyd3fyMDAW9ixw8BLqpZrPaobVDPbcQ1wmGQuMCGEgRDC2OxHSZ0tmbF1L0mgdSvJ\nylt35x6vz217D8ePf4PR0Z3OeJSqNDq6Mxd4OfNRnWvO4AvYCvwI8AJAjPEw8LqFvGkIYXMI4VgI\n4WgIYXcup0xqGyMjI/T3X1iwZQvwAwWPM8AuTpy4kfHxNaxdu94ATJIEVJdw/70Y4z+EUDSUmZ3v\nG+ZyyN4L/FCM8V9CCH9EUsrCxbzUVnbs2JwbHkke9/R8A7iByUlIku7zdYvg5Mnkit6hE6myTCbD\nxMSzpFIbyeb+yrjWozpRNT1fx0IIVwBLQghvCCHcAXx+Ae/5LeB7wGkhhCXAacD/XcDrSXVXrmjq\nyMgIDzyQlJEYGhpjbOwTjI19nKGhMZYufa7JLZbaSz7X6/Dh95LNvodUKs3AwD1T+V4WLlYnqWa2\n48tIxlSGc5sywLYY43fn/aYhbABGgZNAJsb4CyXPO9tRLWM+azm6/qNUm9mW7fL7pFZXj9mO/wz8\neu62YCGE84HrgHOBfwQ+GUK4Isb4Pwr327p169T91atXs3r16sV4e6lmxQnA1Q0h5nvFrFskzd/E\nxLMMD6/j0KEjNX8HpXrav38/+/fvn/fxFYOvEMJngEj5gqoxxrhmnu95MfD5GOPzuff5FPCjQMXg\nS2pHpUUkJVVWumxXT891HDt2KpOT78Xlh9RqSjuFfvM3f7Om42fr+boE+AZwP/CF3LZ8ILaQMcFH\ngRtDCL3Ad4GfAL64gNeT6qrcWo4mAEuLq7S3eGKin8OHryLp7VpGUk8v4XdQ7a5izlcuGX4I+Hng\nIuB/AffHGI8t+E1D+FWSb1QWeBj45Rjj9wqeN+dLLcWK21JjzcwBu56lSz/NqlX9fgfVcmrN+Zoz\n4T73oi8hCcJuBbbGGH9n/k2solEGX2pRBmFSY5hkr3ayqMFXCOGlwE+R1OE6FxgD7o4x1rU0hMGX\nmmGuwMo/BlJjebGjdrFowVcI4ePABcBngT+KMR5dnCZW0SiDLzVYNYHVbFPhJUndq9bga7Yiq1cA\nbwDeD3w+hPBPBbdvLbShUiupZj25iYnnm9I2qdtZYFWdpuJsxxhjNdXvpa6QyWQ4duwIcP3Utp6e\nG0inP968RkldoLRX+uDB9Q73q+1Vs7aj1PHmKicxOrqTycnbSaa87wSe5oIL3ugfAKnO5lPkWGp1\n9m5JzFynsfKV9QiwB7iavr6zGtxKSTBd+d5hSLWrqkpNNJoJ92o1znSUmqP0u9fTcx1wKpOTHwb8\nLqo11KXOV6MZfKkVlE5zB5z2LjVI4fdvcHAlBw48DCQTX6Yr32eArSxd+hy7d3/M76SaxuBLWgTb\nt2/npptGyWZvA7y6lhpptp7m6ZIvy0gCMHuj1XwGX9ICZTIZ3vGOK8hmR7Gml9R4s9XUmw7MzgOu\nLruP1GiLWedL6kqjozvJZt/Q7GZIKpEfilyx4vWcfvrTzW6ONG+WmpDKuhTYNPUoldpIOn1/85oj\ndZFypV8GB68tGopcsuQaQriOZJDkKKnUvUxMXEgmk3HoUS3Pni+J4grag4Mr6e29D7gSuJNUKs3N\nN6c9oUsNUq70y4EDDxfU+1rGiy/2EuMvAb8F3EU2O8rhw1exdu16y0+o5Znzpa5XLrl3y5Zrp2ZX\nObNRar7iPLB1QLn7YO6XmqHWnC+HHdX1ylXQPnDAk7fUSoqHIsvle2XIrz7xta99i+HhdVPHefGk\nVuOwo7qKC/RK7alwKHJg4BR6em4AdgHnAf+FJE1gDXApTzzxDcbH1zA+vmZqGNLvvlqJw47qGpVq\nBwFWr5faTGER1q997Ws88cR1VBqGHBi4i0cf/arfcdWNw45SBZUW6N23bw8PPLCroHq9J2Wp1Y2M\njEx9T4eH1/HEE5X3PX78GRfnVktx2FFdY2Li+YrPjYyMsG/fHtLpDYyO7nRoQmoj6fQGens3MT0M\n+b7c/V309m5i+fJzmto+qZQ9X+oKmUyGY8eOANdPbevpuYF0+uNF+xQOPx48uN6hCakN5PPBpteC\n/FUOHBgDkp5syKcWJPv39m6a2i41gzlf6grF68ElM6IGBk7h4YcPltmnOFekr+8swFlTUjspzAlL\npzcAFD32u6zFZM6XNKuR3G0XfX1jc+x7lCNHHiGbfS9gT5jULir1Yls+Rq3Cni91hdlmOk4PVaxk\n+/Y7pvZJpdIuri21odkW5pbqwZ4vqYzSnJDiPJD81XG+sn3SIzYxcSGHDzenvZJqlx9qPHToCEm5\nienCqxMTpzS3cVIBe77Utea6Oq7UW+awo9R6ir+vR4HfA04DbgWSCTZjYx/3+6u6sOdLWiTless8\ncUutqbSOX9Lrdf3U48lJa3updRh8qWsNDq7kwQc3ks0mj8tNPy8s5Cipnby82Q2QKjL4UlfKZDJs\n334H2ex7gDtJpR5ny5aNBlpSmypeeBt6eh4FbmByMnlsbS+1EnO+1JWcDSV1Hmt7qVnM+ZIKlJ6M\nPflKnatcmkC573y584LnCjWSPV/qWLPNVnQmo9Sdyn33t2y5tqjGn+cD1cqeL4nkBHv55dcUzX46\nebJ4ttOKFa/n+PFtLF9+Djt2eKKVOlVhr9bExPMzzgsf+ci2Wc8V0mIz+FJHyWQybN68Lbcs0Jsq\n7lN45Xvy5KZGNlFSA5V+31OpdJNbJBl8qYNs376dm24azQVdt5Esor1+6vn8bKfSekBe5Uqdq/T7\nns0eJZUqLjHzgQ9cy/btm6ZmSjozUvVm8KWOkMlkuOmm28hmbwPyC2YnC2jDVpYufY7du5Ohxfzw\ng6RudBH9/W+mry85T+SLJ1988cUWVFbDmHCvjpCUjngauJrpHq/ihNoDBx4GZi6gbXKt1LlKlx1K\npe6lv/9CduzYDFiKQouj1oR7gy+1tcKFdE+c+GngPpKg6yjw3zn99KW8+tWv5KmnnmFy8sPAzGDM\nk67U2YpzQW8DoKfnOuDUovNC4UWYpSdUC4MvdY2ZC+neBbwXeIgQjnHKKT28+OIocCdJj5gFVaVu\nNbOw8tuodF6wFI1qZakJdY3Nm3eULKQLS5d+mlWr+pmYeAmHD1+Ve26s0ktIEski3Hdy6NBzUz1e\nTspRPRl8qS1lMhmOHPlKydaLWLXqSfbt28Pw8LqC7RuAK6ceOZNJ6j6V1348SjIx51ZOnIC1a9ez\nYsWKJrZU3cDgS21pdHQn2ey7gekaXanURtLp+4FyJ9oXueCCe+jrO9OZTFIXGhkZ4YEHdhXkcX0C\ngMsvv4YTJ26lsJcL7qK3d+7SE+aFab4MvtTGLiK5Yt0JPE1//5unTn7lTrSeGKXuVm7tx1Wr+hkf\nL96vr+8sHnjgxllLT5TmhR08uN68MFXNhHu1jcKrTMtFSFoM802un5nA70SebmbCvTrSzKvMfLmI\n4kKJklSLmb3knktUf/Z8qS14lSmplViOQoXs+ZIkqc7sMdNC2POltuBVpqRGq3U2o7Mfu1dbVLgP\nIbwK+H3gAiAC74kx/k3B8wZfmsETm6RGqXTBB+XXg/QCsbu1S/C1CzgQY7w7hLAEeFmM8R8Lnjf4\n6hCLHTAZgEmql8Lzy8TE8wWrZADsYmDgLh599KtlAyzzUrtby+d8hRBeCfxYjHE9QIzxReAfZz9K\n7Wix6+BYV0fSYim9kAOKzi+pVHrGMcePP1Ow7FCGkyfP4/LLr2H37o81qtnqFDHGht6AtwJfAO4B\nHiZZDfm0kn2i2t/Q0GUR7o0Qc7d749DQZbMes3fv3jg0dFkcGros7t27t2j70qXnV3y9SsdJUqm9\ne/fG3t6zcueTe2Nv71lxYGCw4PyyN8KKCK+MkI5wSUylzoznn//W3D57IxQf/6EPfWjGa3ou6h65\nuKXqWKgZsx2XACuBX4kx/m0I4Xbg14CbCnfaunXr1P3Vq1ezevXqBjZRzVCuZ2vLlmvZs+fPOHLk\nEbLZN1V9nD1ikiopt3D28ePbcs9mcttvAT5D0j/wUbJZeOqp6+jpuYHJyfNzz08ff+DAmLMfu8j+\n/fvZv3///F+glkhtMW7AMuDJgsdvB/5nyT51iEvVaOWuLme7EizuKdsb4ZwIr4hwScWrzXyPV609\nbJK6V7lzxsDApbnz1SUFz5Xfb7ZeeHUnauz5Si00+ptHsPcM8FQI4Y25TT8BHGt0O1R/+To4Q0Nj\nDA2N1dAblQHeBXwH+Chwdv4VSdZyvJOlS7fZuyVpXtLpDfT2biI5n+yit3cTO3bcyAMP7GLp0udm\nPbav7yx27/7YjOPzeWNSNZo127GfpNRED/AEcFV0tmPXmx4+PK9g69UknaX5YYCZU7id4i2pVpVm\nThefT46SH3aE4nOLM69VqC1KTczF4KszVXOyymQyXH75NZw48a+AS4H7SIKuo8DdDAy8hR07Ns84\n1hOhpMVSeD4ZHFzJgQMPA55bVJnBl1pSLQULM5kMa9a8i8nJJSQ9Xg8RwmNs2/YBtmzZ0qSfQFI3\n8wJPszH4UksqV4CwtGBhT891XHBBP319ZzI4uJI9e/6M48efYfnyc2b0dnkilNQopjZoLi1fZFXK\nKy1YODm5JFdRGg4erHxys7SEpEYqV5pidHSn5xzNm8GXGiKd3sDBg+s5eTJ53Nu7ieXLV3DiRH6P\nncCtVHNy80QoSWpnDS81oe6ULzsxMHAXS5duY8WKFaxbN1QwXfvpZjdRksoqV5qitLREJpNheHgd\nw8PryGQyTWmn2oc9X2qofI7XiRPw6KOb2LLlWg4cGGNi4hSOHbuByclkv+Tktqvsa5TrRau0ryQt\nVP7isVL1+kamQpjv2hlMuFfDlEu6HxoaY9++PUBtJxVPQJIardJ5Z65z22K+v4n/rcmEezXVQoKi\nkZGRqvevZV9JWqj59G4t5HxY7ljzXTuHwZcWzVwnJ4cLJbWr2QKfcue2wcFr5z0UWelcqs5h8KVF\nM9dV2Vx5E5LUjsqd2xbSS1XpWC9gO4ezHdVQIyMj7Nu3h3R6A6OjO50ZJKktzDXjsfTc9jd/87eL\n3oZ8kDc0NMbQ0Jj5Xm3MhHstmtKu8sKK9ZUXrk2u3pJZj66fJql1lcvDKl0Hcvv2Ozh58kqSBbl7\nSOoXQk/PDYyNfRyYuaRaufcxsb69uLyQmiaTybB58zaOH3+GM844jaeeeobJyQ8DxSePmTODrieV\nupts9rYZ+0pSqyoNklKpNNnsKDAGrAGWkRSQfpqBgVPYsePGqoMqZ3S3F4MvNUXlk9B0gLV06adZ\ntaqfiYnnc8sI5Z97G3A19Z6mLUmLaeaFZP5clg++iteyPX78GU6cuBHPdZ3HUhNqqPzV2aFDR4oS\nRLPZOwv3AnZx4sStjI8nw5E9PdMFVVOpx8lmG91ySVpsl5JKbSSbfQ9w/dTWnp7rOHbsVCYnz29e\n09RSDL40L/khxiNHHskNF5YuD5Q/CQHcSeG6jZOTMDBwF319YwAMDm5k+/ZNzuCR1FZmzj68jy1b\n0hw48DATE28C7qGv70wmJvpzvf3LmO718lzXzQy+VLPpIcbzgNtITibLgCun9pk+CY1x6NBzBQto\nJ/r6zirqar/44ostQSGprVQqn7NlS/F+w8Pr8keQzJbcytKlz7F7t+e6bmXOl2o2nedQmtcwndc1\n1+xGE+oldYvZzoEm1ncGc77UQBso7kK/r+yVnMVVJXWz/DkwPxt8+fIVQGMX5FZrseeri8x1hVXt\nFVjxCeMoqdS99PdfyI4dmz1pSFKB/Hl1YuJZjh17rKj8ztlnv4YnnriOamc/2kvWumrt+SLG2HK3\npFlaTHv37o29vWdFuDfCvbG396y4d+/eqp8v93pDQ5fFoaHLZt1PkrpV8Xn1kty/MXdLRzijZNu9\ncWjosipea+5ztBorF7dUHec47NihSq+Q5lpnrNZ1yEZGRrzqkqRZFJ9Xx0qefQh4D7BpaksqtZF0\n+v4qXmv6HJ1/DuwNaycGXx2oXB7BihUrmtwqSepmGyicEZ7UN7yaZPZjUgW/v//NRcFT4UX0xMTz\nM15xYuJZc8balAtrd6DiK6T8F/PFWReFnWvR2FplMhmGh9fNunB2NftIUrsqPq8+Q0/PiwwM3MPQ\n0Bg337wx99wzwBp6e59kx44bp47NX0SPj69hfHwNx44doafnBgrP0bCk4Fy/jJMnz+Pyy6/xfNoO\nahmjbNQNc74WZGjosrJ5BHPlae3duzcODFwaly49Pw4MDMYPfehDc+5f7vlqchPMX5DUDWY775Z7\nLr9t6dLzZ5zHBwYuLdp/YGAwt8/eCJ5Pm4kac76aHmiVbZTB14LMN7ApPi4d4RVVJ+j39LwqDgwM\nxqGhywpOCJWTSCsFiJLUrcon6O+NcFmES+LAwKVF+/b0vCpCX5lkfs+njVZr8OWwYwfK15QZGhpj\naGisKAdgtqG+4uHKJ4GPUjh0mc89mLnvMiYnl3D48FWMj6/hyJGvNOYHlaQOUnxe3Qq8jyRPbA1w\nNceOPTZ13h4d3cnk5O3AfcBzzWmw5s2E+w5VOhtx5lqM80vOLFxIOzkhQJIsOr12YzZ7tGBdx/Lr\nl81cE801ziRp2gjwgyQLdE+viztzFvoI8DGKl3fzfNrqDL46SKUCfOXXYpxZTqI4IDqP5Kor0du7\nicHBawtm1hQ+X7qo9kX09795auHschXtrXovScVKL0pTqb+fuoida9+enhe54IJkIe/ZzqcWam0N\nVrjvELOtHVZ+LcYM04u7fgxgqgozLKGv70wGB1dy4MDDwHStsOR1itdyXL582YzKzU53lqTqFVbC\nLzwHb99+R8V1cWsNpFxnt36scN+lZktgn34uPyMmnUvSTEe4JIZwelyy5Mw5E/Rnew8r3kvS/Mw2\nSWqh59b88QMDl8bTT/9BE/PrBCvcq9R09/QtJHkB9wBXkSRq3kKMd/Lii1dTOBy5efO2GVdUs+Vp\nWfFekuZnOtF+GbCTkyfPY/PmbVPn1fnk5QIFPWdXAn8FvL4ezdc8GHx1iEqBUf6LuGLF64EkH2Bi\n4iIOH34IqLTsxdFcYv57geLEfPO0JKkejpIsNZQMCR45spFMJlNz4FU4rPjgg2my2VGSc/ytJMHd\n+qn9TcxvHoOvDlEuMALKju8DvOMdV+QSOTPAs8B1U6+VSt2bmxHpGmKSVG/p9AYefPCKXKCUnzU+\n+/q65RSXqsiQzfaW7DFCUiE/n++b/D0YHl431Q7P641h8NVBSrunh4fXlV2Idd++Pdx880ZuvPFX\niPGlJFdERwnhA7z1rRcBF3L4cP5VMsCdfO5zj/Hnf/4XxHg74BpikrRYRkZG6O8vPO8uVIbkvP92\nkgvrXyIpWZHo7X1yKvBybcjmMPjqUlu2bGHPnnEOH76KfHAW40X09Y2RTm/IfSGPklwlree7330M\nuJ1KZSokSdUrnam4Y8fm3Hk3eX4+Q4LT6SfnkeT33kcSeD0EfJfzz7+d173udVMpI5Uu0D2v158V\n7jvYXItl9/WdWfa4/BDm0qWfJukVexJ4YwNaLEmdr3TR7LVrk+Cn0sok1Zo+dz9HEnDdQnIO/2vg\nd3jd617Hvn17DK5agD1fHWyuBPm5Zi+uWtXP+Hh+70tJkkETqdRG0un7G/BTSFJnKc7NKk4JWWhx\n1JGREXbv/lhBXm9lg4MrefDB2VcjUX0YfHWY0mnGhUVSa60yPx2cXUl++BHuJJV6nJtvTnv1JEkN\nUDqLca7crJGREW6+eSM33VQ5sMpkMmzffgfZ7HvIn9e3bNnoeb1BrHDfQYq/oEeBu0gWx55/JeNy\nVZedESNJ81drpfnpVUryZSKS4cl9+/bM+T6VestWrlxdlPNb7WuqvFor3Nvz1UGKu7LXkQReC0uk\ntHiqJC2uxaiZODHx7JwlIiqdvzOZDEeOfGUeLddiMfiSJKnBarmwnbmI9nUcO3Yqk5MzC2FXY3R0\nJ9nsuzGPt3kMvjpI8Rf0POB9U8+ZSClJ7am0p2xior9oyHB+hbAvIsnl3Qk8TX//mx3laCBzvjpM\nLQn3kqT2Uy4HbGDgLh599KtV5ZHVmnOmudWa82XwJUlSGykXPK1YsaKmBPpqS1eoOibcS5LUwcol\n7Ofv1/IaBlzNY/AlSVKbKRc8VSqardbTtGHHEMIpwJeAb8QY/0PJcw47SpJUA4cSm6dtcr5CCB8A\nVgGnxxjXlDxn8CVJktpCrcFXUxbWDiGcA7wD+H2g6sZKkqT5yWQyDA+vY3h4HZlMptnN6WrNyvm6\nDbgBeEWT3l+SpK5R6/qQqq+GB18hhHcC34wxHg4hrK6039atW6fur169mtWrK+4qSZJmUbz83PyX\nnFNi//797N+/f97HNzznK4TwX4FfAF4EXkrS+7UnxviLBfuY8yVJ0iKZ7+Lcqk7bJNwDhBAGgeud\n7ShJUv1Y1b6+2rHIqlGWJEl1VK4wq4FX87i8kCRJ0gK0RakJSZKkbmXwJUmS1EAGX5IkSQ1k8CVJ\nktRABl+SJEkNZPAlSZLUQAZfkiRJDWTwJUmS1EAGX5IkSQ1k8CVJktRABl+SJEkNZPAlSZLUQAZf\nkiRJDWTwJUmS1EAGX5IkSQ1k8CVJktRABl+SJEkNZPAlSZLUQAZfkiRJDWTwJUmS1EAGX5IkSQ1k\n8CVJktRABl+SJEkNZPAlSZLUQAZfkiRJDWTwJUmS1EAGX5IkSQ1k8CVJktRABl+SJEkNZPAlSZLU\nQAZfkiRJDWTwJUmS1EAGX5IkSQ1k8CVJktRABl+SJEkNZPAlSZLUQAZfkiRJDWTwJUmS1EAGX5Ik\nSQ1k8CVJktRABl+SJEkNZPAlSZLUQAZfkiRJDWTwJUmS1EAGX5IkSQ1k8CVJktRABl+SJEkN1JTg\nK4Tw2hDCX4YQjoUQvhJCeF8z2tGp9u/f3+wmtDU/v/nzs1sYP7+F8fObPz+7xmpWz9f3gI0xxguA\nS4BrQgg/1KS2dBy/RAvj5zd/fnYL4+e3MH5+8+dn11hNCb5ijM/EGP8ud//bwP8Gzm5GWyRJkhqp\n6TlfIYRzgQHgC81tiSRJUv2FGGPz3jyElwP7gQ/FGD9dsL15jZIkSapRjDFUu++SejZkNiGEU4E9\nwH2FgRfU9gNIkiS1k6b0fIUQArALeD7GuLHhDZAkSWqSZgVfbwf+CvgykG/A5hjj3oY3RpIkqYGa\nmvMlSZLUbZo+27GSEMKHQwj/O4RwJITwqRDCK5vdplYXQvjJEMKjIYTHQwibmt2edmLh38URQjgl\nhHA4hPCZZrel3YQQXhVC+JPcee+REMIlzW5TuwghbM59d4+GEHaHEF7S7Da1shDC3SGEZ0MIRwu2\nLQ0hjIcQHgsh7AshvKqZbWxlFT6/mmKWlg2+gH3ABTHGfuAxYHOT29PSQginAL8D/CTwZuDnLVxb\nEwv/Lo73A48wnU6g6v034LMxxh8C3kJS/1BzyJUrei+wMsZ4EXAK8K5mtqkN3EPyt6LQrwHjMcY3\nAg/mHqu8cp9fTTFLywZfMcbxGGM29/ALwDnNbE8b+GHgqzHGr8cYvwd8AviPTW5T27Dw78KFEM4B\n3gH8PuCM5RrkrpJ/LMZ4N0CM8cUY4z82uVnt4lskF0+nhRCWAKcB/7e5TWptMcbPAS+UbF5DMhGO\n3L8/3dBGtZFyn1+tMUvLBl8l3gN8ttmNaHE/ADxV8PgbuW2qkYV/5+024AYgO9eOmuE84LkQwj0h\nhIdDCHeFEE5rdqPaQYzxBDAK/P/A08A/xBj/vLmtaktnxRifzd1/FjirmY1pc3PGLE0NvnLjy0fL\n3P5DwT5bgMkY4+4mNrUdOMyzCHKFf/8EeH+uB0xVCCG8E/hmjPEw9nrNxxJgJfC7McaVwD/jsE9V\nQgjnA9cB55L0Vr88hHBFUxvV5mIyE8+/KfNQbczStCKrADHGodmeDyG8m2QY4981pEHt7f8Cry14\n/J6Nf5oAAAb3SURBVFqS3i9VabbCv5rTjwJrQgjvAF4KvCKE8Icxxl9scrvaxTeAb8QY/zb3+E8w\n+KrWxcDnY4zPA4QQPkXy+/g/mtqq9vNsCGFZjPGZEMJrgG82u0HtppaYpWWHHUMIP0kyhPEfY4zf\nbXZ72sCXgDeEEM4NIfQAPweMNblNbSNX+PcPgEdijLc3uz3tJsb46zHG18YYzyNJdv4LA6/qxRif\nAZ4KIbwxt+kngGNNbFI7eRS4JITQm/se/wTJpA/VZgxYn7u/HvACtAa1xiwtW+crhPA40AOcyG36\n6xjjf2lik1peCOHfA7eTzPb5gxjjjiY3qW1Y+HfxhBAGgXSMcU2z29JOQgj9JJMVeoAngKtMuq9O\nCOFXSQKGLPAw8Mu5iUcqI4RwPzAI9JHkd90E/Cnwx8APAl8HfjbG+A/NamMrK/P5fZBkdmPVMUvL\nBl+SJEmdqGWHHSVJkjqRwZckSVIDGXxJkiQ1kMGXJElSAxl8SZIkNZDBlyRJUgMZfEmaIYTw0yGE\nbAjhTbPsc24I4WiF5+4KIazI3f/1erWzEUIIZ4cQPlnH1/96CGFpvV5fUuuxzpekGUIIfwT0Ag/H\nGLeWeX4JcA7wmRjjRXO81j/FGE8vsz3A1Dpyc7UnFWPsyAW7QwhPAqtyC0RL6gL2fEkqkltc/EeA\nXyFZpiq/fXUI4XMhhD8FvkKyEsCSEMJ9IYRHQgifDCH05vbdH0JYFUL4LaA3hHA4hPDxEMLyEML/\nCSHsAo4Crw0h/G4I4W9DCF8JIWwteL+vhxB+K4RwCPi13L/5595Q+Di3bUkI4Yu5CvuEEHaEED5U\n5ud7WQjhz0MIh0IIXw4hrMlt/9chhCMhhJfk9vlKCOHNhT18IYQLQghfyP08R0IIry957f8UQvjt\ngsfvDiHckbv/6RDCl3Kv+94y7SrqSQwhXB9C+GDu/vkhhD/LHf9X+R7JEML/F0I4GkL4uxDCgVn+\nWyW1khijN2/evE3dgCuAO3P3/wpYmbu/Gvg2sDz3+FyS5Vzelnv8ByTLCgH8ZcFx/1Tw2ucC3wd+\nuGDbGbl/T8kdd2Hu8ZPA9QX7/QXQn7v/X4FryrT9zSTr+v0EyTIzS8rscwpweu5+H/B4wXPbgA8D\nvwNsKmjz0dz9O4DLc/eXAC8tee3S1/ss8KMlP2cvSeB5RsHPubTwfXLb08BNufsPAq/P3f8R4MHc\n/S8Dr8ndf0Wzf3e8efNW3c2eL0mlfh7I5zh9Mvc474sxxuMFj5+KMf517v59wNureP3jMcYvFjz+\nuVwv1sPABSQBVN4fFdz/feCqEEIK+Flgd+kLxxgfybXjMyRrI75Y5v1TwI4QwhFgHDg7hPDq3HM3\nA8PAxcBvlzn288Cv59YSPDeWLKAbY5wAvhZC+JEQwpnAihjj53NPvz+E8HfAXwOvBd5Q5vVLhRDC\ny4AfBT4ZQjgM3Aksyz3/ELArhPDLJMGgpDbgl1XSlFzi948DF4YQIkkvUQRuyO3yzyWHFOZrhZLH\nlUy9RgjhPJIenotjjP8YQrgHeGm5fYE9JAvY/gXwpRjjCxVe/yLgBeCs3Hv8CEnAAskCwmeS9FCt\njDF+P5dzlX/PPuBlJD93L/Cdoh82xvtDCH8DvBP4bAjhP8UY/7Lk/T9BEhw+Cnwq14bVwL8DLokx\nfjeE8JclPyfAixSngvSSfJ4p4IUY40DpDxpj/M8hhB8Gfgo4FEIwd0xqA/Z8SSr0M8AfxhjPjTGe\nF2P8QeDJEMKPVdj/B0MIl+TuXw58rsw+38sl6JfzCpIA61shhLOAf1+pYTHGfwEywO8B95TbJ4Rw\nGfAqYBC4I4TwyhjjF2KMA7nbZ3Lv+c1c4PXjwPKCl/jvwG+Q9KrdUub1XxdjfDLGeAfwpySBXqkH\ngJ8m6TH8RMHP+UIu8FoBXFLmuGeBV4cQloYQXkIS4BFj/CeS/4OfybUhhBDekrt/fozxizHGDwLP\nkUyCkNTiDL4kFXoXSfBQaA9JIBEp7tmKwP8BrgkhPAK8kiQwKrUT+HII4eOlrxFjPAIcJukl+n/t\n3DFKQ0EURuFz92AtuAfdQBZgaSOx0SXYp4zYuQUhks5WrBUhCIEkxMJN2IrFtZgRJDwIRJgYOF85\njxnudJeZf94IeFpT3x0lZ/a4+iEi9oAhcJGZ75Tc1k3HGiPgMCJmQB94K9PjDPjMzDFwBRzVE6uf\nvQKc1MD8lHJFeru6eGZ+UHJn+5n5WocfKI8TlrXGl455X5Rrz0nd3/LX51PgvF5bLoDjOn5dHw3M\ngefMnHXsV9I/468mJO2MiLikhOUH265FkjZl5kvSToiIe+AA6G27Fkn6C0++JEmSGjLzJUmS1JDN\nlyRJUkM2X5IkSQ3ZfEmSJDVk8yVJktTQN0FrDqP5yG5jAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "\n", "plt.figure(figsize=(10,7))\n", "plt.scatter(x_vals, Data, marker='o')\n", "plt.ylabel('Measurement')\n", "plt.xlabel('Arbitrary x-axis values')\n" ] }, { "cell_type": "code", "execution_count": 22, "metadata": { "collapsed": false }, "outputs": [], "source": [ "eps_data = sig\n", "def residual(vars, x_vals, Data, eps_data):\n", " amp = vars[0]\n", " centroid = vars[1]\n", " width = vars[2]\n", " offset = vars[3]\n", "\n", " model = amp*np.exp(-(1.0/2.0)*((centroid-x_vals)/width)**2)+offset \n", "\n", " return (Data-model)/eps_data" ] }, { "cell_type": "code", "execution_count": 27, "metadata": { "collapsed": false }, "outputs": [], "source": [ "def mwgaussfit(x_vals, Data, eps_data):\n", " vars = [A_g, B_g, C_g, E_g]\n", " out = leastsq(residual, vars, args=(x_vals, Data, eps_data))\n", " print('Best Fitting Parameters: \\nA = '+str(out[0][0])+'\\nB = '+str(out[0][1])+'\\nC = '+str(out[0][2])+'\\nD = '+str(out[0][3]))\n", " return out" ] }, { "cell_type": "code", "execution_count": 28, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Best Fitting Parameters: \n", "A = 10.0332472863\n", "B = 5.00660239344\n", "C = 1.99338121722\n", "D = 4.0142970844\n" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 28, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlAAAAGnCAYAAAB4nWl5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XlclOX6x/HPMwo6mpKESraoaWW5cNQ226A8MOUx1xbT\niqwkWyx1NDLSyDXLySyzol8mLdqqRXUOj2RJ6fHUSc1Mj7nvW4hroqjz/P4YQFYFBGaG+b5fL17C\nM8/MXIzMzDX3fd3XbViWhYiIiIiUns3bAYiIiIj4GyVQIiIiImWkBEpERESkjJRAiYiIiJSREigR\nERGRMlICJSIiIlJGp0ygDMOYbhjGLsMwlhc6PsgwjP8ZhvG7YRgTKzdEEREREd9yuhGod4Fb8h8w\nDOMmoBvQzrKsNsCkSopNRERExCedMoGyLOtHYG+hw48AEyzLOpZzzp+VFJuIiIiITypPDdTFwI2G\nYfzHMIz5hmFcUdFBiYiIiPiymuW8TgPLsq4xDONK4BPgosInGYahPWJERETEb1iWZZT23PKMQG0F\nZufc0X8Bt2EY55QQiL7K+fXcc895PQZ//tLjp8dOj59/funx02Pnra+yKk8C9QVwM4BhGJcAwZZl\n7SnH7YiIiIj4pVNO4RmGMQuIBM4xDGMLMAqYDkzPaW2QDdxX6VGKiIiI+JBTJlCWZd1dwkX3VkIs\nkk9UVJS3Q/BrevzKT4/dmdHjd2b0+JWfHruqZZRn3q9UN2wYVmXdtoiIiEhFMgwDqwxF5OVZhSci\nIlLtGEap3zvFz1XEAI8SKBERkRyaOan+KipR1mbCIiIiImWkBEpERESkjJRAiYiIiJSREigRERGR\nMlICJSIi4kcsy6JevXps3LjR26EENK3CExER8WHNmjVj9+7d1KhRA/CsIluzZg3h4eEA3H///Vxw\nwQWMGTPGm2EGHI1AiYiI+DDDMPj66685ePAgBw8e5MCBA3nJk3iPEigRERE/Y7PZWLduHUlJScyc\nOZMXX3yRevXq0b17d2+HFjA0hSciIuLjimvwaRgGcXFxLFq0iAsuuIDRo0d7IbLApQRKRESkNCpq\nq5cydju3LIsePXpQs6bnLbu4TYPVQb3qKYESEREpDS8lKYZh8OWXX3LzzTfnHbPZVIHjbfofEBER\n8WPaBNk7lECJiIj4scaNG7N+/XpvhxFwlECJiIj4mfyjTg8++CArV66kQYMG9OrVy4tRBRajsgrP\nDMOwVNQmIiL+wjAMFWMHgJL+n3OOl3o+VCNQIiIiImWkBEpERESkjJRAiYiIiJSREigRERGRMlIC\nJSIiIlJGSqBEREREykgJlIiIiEgZKYESERERKSMlUCIiItXUhAkTGDBgQIWfezo2m63aby+jTuQi\nIiL4fifyGTNm4HK5WL9+PfXr16dnz55MmDCBkJAQb4dWhM1mY+3atVx00UVFLouKiuKnn34iKCgI\nwzC4+OKLueOOOxgyZAjBwcFnfPuno07kIiIiPuDw4cP06zeAs89uwoUXXs6XX35Z4ffhcrl4+umn\ncblcHDhwgP/85z9s2rSJ6Ohojh07Vux1Tpw4UeFxVATDMHj99dc5cOAAO3fuxOVy8dFHH9GlS5cy\n3Y63k10lUCIiIqfw+++/c+WVNxMe3pLu3fuyZ8+eApc/+OAgZs/ey/79i9iyZSp33x3HL7/8UuCc\n48eP8+233zJ79mx27dpVpvs/cOAAiYmJTJ06lZiYGGrUqEHTpk355JNP2LhxIx988AEAiYmJ3H77\n7dx7772EhIQwY8YMEhMTuffee/Nu67333qNp06aEhYUxduxYmjVrxnfffZd3/dxzN27ciM1myzu/\nYcOGjB8/Pu92fv75Zzp16kSDBg1o0qQJgwYNKjGRK05u8mO324mMjCQlJYVFixbxzTffnPb2b7zx\nRgAiIiKoV68en376Kfv27aNr1640atSI0NBQbrvtNrZt21amx7mslECJiIiUICMjgxtuiGHx4rvY\ntesb/vWvBtxyS+8Cox/ffPM1R45MAZoCN5OdfT+maeZdnp2dzY033krPnsPp3/9dLr64HUuWLCl1\nDP/+9785cuQIvXr1KnC8bt26dOnShbS0tLxjKSkp3HHHHezfv59+/fphGCdnpFauXMljjz3GrFmz\n2LFjB/v372f79u15l+c/N9fChQtZvXo18+bNY/To0fzxxx8A1KxZkylTprBnzx4WLVrEvHnzmDZt\nWql/p8L3dcEFF3DFFVfw448/nvb2f/jhBwB+++03Dh48yB133IHb7ebBBx9k8+bNbN68GbvdzuOP\nP17qeMpDCZSIiEgJFi5ciNv9NyzrYeBSjh17jd9++7XAKFTduvWBDXk/BwdvKFCX9M477/DrrzYO\nHfqFAwe+4uDBScTGlv7NPSMjg7CwMGy2om/Z4eHhZGRk5P187bXX0q1bNwBq165dINH77LPP6Nat\nG9deey1BQUGMHj26QCJT3JTYc889R61atWjXrh0RERH8+uuvAHTo0IGrrroKm81G06ZNiYuLIz09\nvdS/U3GaNGnC3r17y3X7oaGh9OzZk9q1a3PWWWfxzDPPnHE8p6MESkREpAR169bF7d4JuHOO7MXt\nzsZut+edM2XKeOrUuQPDSKB27T40bryC++67L+/yjRu3kJV1PVAj58gNbNu2pdQxhIWFkZGRgdvt\nLnLZjh07aNiwYd7P559/fom3s3379gKX2+12zjnnnFPed3h4eN73derU4a+//gJg9erVdO3alXPP\nPZeQkBASEhKKTG2W1datWwkNDS3X7R8+fJiHH36YZs2aERISQmRkJPv376/UOiklUCIiIiWIjIzk\nsstCsNt7AJOoW7czjz/+BHXr1s0758477+Dbb2czalQwL7zQiWXLFlG/fv28y6+99mrq1p0JeBKx\noKApXHXVVaWOoVOnTtSqVYvPP/+8wPFDhw6RmppK586d844VNw2Xq0mTJmzdujXv56ysrHInPY88\n8giXX345a9euZf/+/YwbN67YBK+0tmzZwpIlS7jhhhvKdfsul4vVq1fz888/s3//ftLT07Esq1IT\nqJqVdssiIiJ+LigoiB9++BfTpr3B+vWbue66ePr06VPkvE6dOtGpU6dib6N79+4MHryMiRObY7MF\n0bp1BO+/P7vUMYSEhPDcc88xaNAg6tevz80338y2bdt49NFHueCCCwoUiZ9K79696dSpE4sWLaJj\nx44kJiaWO8E4dOgQ9erVo06dOqxatYo33niDRo0alfr6ufd7+PBh/vvf/zJkyBCuvvrqvJV4p7v9\nxo0bs27durw2BocOHcJutxMSEkJmZibPP/98uX6vstAIlIiIyCnUrl2boUOHMHXqZO6+++5TjvKU\nZOzYURw8uJedOzexePEPBabdSmP48OGMHz+eYcOGERISwjXXXEPTpk2ZN28eQUFBgGf0qXBs+Y+1\nbt2a1157jT59+tCkSRPq1atHo0aNqFWrVrHXP9XvOWnSJGbOnEn9+vWJi4ujT58+pb4uwOOPP079\n+vUJDw9nyJAh3HHHHaSmppb69hMTE4mNjaVBgwZ89tlnDB48mKysLMLCwrj22mu59dZby/X/VBZq\npCkiIoLvN9KsaIcOHaJBgwasXbuWpk2bejucKqNGmiIiIlImX331FYcPH+avv/5i2LBhtGvXLqCS\np4qkBEpEvMI0TWJiehMT07tAzxwRqTwpKSmcd955nHfeeaxbt46PPvrI2yH5LU3hiUiVM02Tnj1j\nycqaCIDdHs+cOck4HA4vRyaBLNCm8AJVRU3hKYESkSoXE9ObtLRuQGzOkWSio1OYO/fzU11NpFIp\ngQoMqoESEb+Sf8ouI6PsvWc05ScivkQjUCJS6QpP2QUHDwaCyM5+CTj9FN7J698DLMRmW8Po0UNI\nSEgocI7LlQSA0xmn6UApM41ABQZN4YmI3yhuyq59+7cJC2sMnD7h8Vy/OfAB4EnCbLYh/POfs3A4\nHKqpkgpR2X2DxHdURAKlTuQi4hVhYY3LWPO0EE/y5EnC3G5wuZJwOBy4XEk5yVMsYJKV1Zy+fR9j\n5szXlURJqelDv5SFaqBEpNI5nXHY7fFAMpCM3R6P0xl32uvl1j1lZOzCMFaXeM7ixctyf8KTRA0k\nM3MkPXvGql5KRCqFpvBEpEqcrkap8OVAgWm5mjUf48SJICzrFcAzTZeQMIhx417LqY1KBloCA9Hq\nPhEpqwqdwjMMYzrwD2C3ZVltC13mBF4CwizLyixPsCISOBwOR4lJU0bGLlasWJ1XVL5gQSytWrXK\nSZ7CgSSOH29LixYHuOiiFACczuRCU3fRwENV+SuJSAA73RTeu8AthQ8ahnEBnlerTZURlIj4p7K0\nGsgt/E5La87SpZtykqdYwDPqtGnTVmB5zrFuwEA2bNiB0xnH3LmfF1Pb5ADuwmYbQlmnCkVEyuqU\nI1CWZf1oGEazYi56GXgK+LISYhIRP1R4JdyCBbGnXAnnGT26B8/KuuZFLj8/pCbtMt/iGm6hBT/S\nhO00codQ784+0PoyqFePCfv3EWE8wmJrJYtoiWH/gIQEJ+npJ0epTtUaobgpRbVDEJHSKPMqPMMw\nugNbLcv6TUs+RSRXwek0yMpaTt++j9GxY8QpEpHclXXhQCyN2UdvfuEuYxYdNpzgV+rwEwf5iRh2\ncC67+S9XtfqBqZPGsXjePD4c8xIXWpfRm2m04RCbQ8JJ+/Az6jY+n4FPP16KvlIFkz2gTEmgiASu\nMiVQhmHUAZ7BM32Xd7ik8xMTE/O+j4qKIioqqmzRiYifMoFkMjMnkZZWMBHJX/dkGKuxLOhAGE4u\n41acpAXbmX1+K2LWP8lRLsSTkN0N7MFuf4sxo5Ph2msZkegi7dhr5CZsZ/EE0bum03fn+cz43zy+\n+fFbFrz+Ctc//HCR6Iome+SNOhV3XAmUSPUzf/585s+fX/4bsCzrlF9AM2B5zvdtgV3AhpyvY8BG\noFEx17NEJHCkpqZadntjC2ZYcE3Ov1bO1wwrOrpXoXNmWH+rYbdSqWltJtQayl1W49oNrdTUVCs6\nule+66dacI0VGtrCSk1Nzbu/gudYBe6zAXuseG63MoNrWVbv3pa1YkWBOENDWxR7++3bRxYbt4hU\nfzl5y2nzotyvMo1AWZa1HGic+7NhGBuAjpZW4YkEPIfDwZw5npVxixf/SWYxrwq5Iz+hdGU8z9Dj\nRBAzzj+XVy5txwnbMZKd7wOQkbELm20Ibrfnenb7BmbOLDiV5nTGsWBBLFlZnp9ttjV55+8llIl0\nZeX12aRccw1ERsL99/PttdfSs98jObVXw/AUqScDk8jMhEOHBhMcPJzs7Nz7jcfpTK6Mh0tE/N2p\nsitgFrAdOApsAfoXunw9EFrCdaskYxQR31N4pMlub5w3snQbT1rbONd6hSess3m9wAhPwes5LZvt\nHKt9+8gCI0+F7yc6upcVHd3LGjt2bLH3aVmWZe3YYVn33WftrmW3/s6wfCNP5xcZcWrf/rq820xN\nTS1wHyXFISL+jzKOQKmRpohUqPw1TlCTsLBzPEXk11/Pth49ODLve2KteBZySZE96zp0iGLp0v6U\ntxFm/hV0kZEdSE9fApxcTRff8QYGLVnBx/TnGcaTTRSnarypPfZEAof2whORKlU4afF0Bi+UcFx0\nEVxzDeddcQXfzv6MOtPeJ5pVBdoMmKbJsmW/l+n+Cq/uy23WWdIqu5vHP8s1Pe7l1SPpLKQVdwVl\nstUoecqupGJzJVAiogRKRMqtcKIyb54Tt9tF/oTjh6dH49i2BkaPhocf5u+Gwd979ChyWy5XEm73\n/UB83jGbbQhO56wS76+kNgPFJT4jRowhLKwxjS67hLFWDR756yi/7tzDoPBwvtozhqZNz2fCBI0u\niUjpKIESkXIrnKi43W8WuPxB0hm6cjHM/x46dSrFLbbFU9SdBGwnIuLyAglN+UeElrNs2Urc7gGA\nZ6Rpd8Ig7hs9ibfWbCSbvnyR9UWRaxUuVFdRuYjkUgIlIhXoupzVcxaj+JJY4ytWvpHEDaVInk4m\nKxOBbtjt8UyYUL5k5eRtLcfTrHM1bvdk8ideL788hszsV1hDR1K5hbOz2hRo/Ame5KxVq5bAuzm1\nXBqhEpEcZak4L8sXWoUnUu0Vt9pu7Jgx1ucXtrRW1wuxJsfHW+3bX5fXY+l0q9hOt+KtpNV9xRk7\ndqxlszUosS9V/l5QzZhhrcFmJdDLgnet4OCzreDghqW6HxGpHtAqPBGpSgWKuocOwGGa8OOPzHv6\nabr0G0B2dk1gEgDBwcNJSXn/jEZxSrtXXUxMb9LSuuEZdTKBe/LisNvjSUgYlK/g/U0a04dveZtP\nuJMx/ItTrc4TkeqnrKvwlECJSMWwLHA64YcfIC2NmLseIi1tO5WViJwukTqZQIXjqalaid2egd0e\nklMwPgIgp/HnMjIzR9KIWqQzgHex8SKvVkrcIuKbyppA2SozGBGp3kzTJCamNzExvVkTGwvffw9p\nadCgQQnXWM7ixcuIiemNaZpndL89e8aSltaNtLRu9OwZW+T2nM44goMH4xl56gb8g6ysbDIzR7J0\naX969vQkR3Pnfs7Mma8THDyY3QyiM2MYgI1BPIynoD05p3g8rtzxikg1VJb5vrJ8oRookWotfz3S\nAO631hk1rO9mzixweXDw2RaE5XUWh/oVUldUdB+84vesK7i33amvk//cC9lobeEsq3/dxupALhIg\nKGMNlEagROSU8o8y5R/lyW0p0J0QEjGJsWK5/fGReec5HA5SUj6ifftLCQ0dQ716n0LetJhntV3u\nFFxlCQs7p1znbqYpXRnOy9n7mPvsk1p5JyJFlSXbKssXGoES8XvFrrIbO9aKju7lWVlHorWbMKsD\nr+UbaSp+dKm0o0bljev0q/ZOPQJW3G3+d/x4y2rUyLJWriwxDu2TJ1I9UMYRKCVQIlKiokmPM681\nQDgPW5swrF48VmybgMLJUVlaEJRGaZOXwhsOn65NQpHLp0+3rIsvtqzMzEr9fUTEu8qaQGkVnoiU\nqGArAIBOwEBqcyfpRJJCMG+E7gQgM3Mkp1u1VtoWBD5n8GBYtQq++Qbz228LrNzTSj2R6kGbCYtI\nhSm8lYnNtga3G15jEOu5iHHcQnTHr3A643L2qPOcV9KWJ7mb/fqL3ITP5nbz4Z5dHLjrLnr+c0FO\n76jthc4+ucLQr5JDESkXjUCJyCnlHzWKjOzAxudfZOixulzJ7WTZZhIR0aZATyXws9GlEhTeuDg8\n6An+ffwvhlkDmc1UCjbnXA68jadI3pNAFrfJsYj4LjXSFJHK8/vvZF9/PX3PvYg5qzfm7C9XPROG\n4jqZX0E437CJq1nGRpoDwwgN9WxCrOk8Ef+mRpoickZKalvAoUNw++0ET5nCgQua59uct2paEnhX\nEjCJX5jEBCw+pjNBvIPd/gFDh/b3dnAi4gVKoEQkT4kdvi0LHn4Yrr8eYmNPf0PVgNMZh90ej6cb\n+cl6p1dox0524AoanLefXmZmD2AYhTuXl5iMiojf0xSeiOQpuuouZyqqzz9gyhT4z3/Abi9SH1Qd\np/DgZP1XRsYuli//nePHg4BJNOAQS3mCyRddxpT1wzk5zZdIaOifzJz5OkBAPEYi1YVW4YlIhQo/\n/BfEx8P8+WC3A57VdHPmJOcrGq+eiUH+VYMdOkSxdGl/IJa9wN3s4IsNL/AJe9nhORvYSceOKTgc\nDmJieuckT55kNCvLU2RfHR8nkUCkBEpE8hRuW1C39lNM2dcQRoyA1q0LnOtvLQnOVOFtYRZxMV80\nbsL03fHc6j4bMEps3yAi1Y+m8ESkgPxtC167sAGXrl0L330HtsAumSxu2vKLT/6PawYP5v9qhZB6\nXrMC7RsCZZpTpLpQGwMRqRi//QadO8Mvv0DTpoAfdxKvIMX+/suXw803F3icTnm+iPgkJVAicuaO\nHoWrroIhQ+D++wGNqJzSxIkwdy6kpQX8SJ2Iv1IfKBE5c+PGQfPmBVoWuFxJ+YqiA6H306kVaFHQ\nti0cPgzTpnk7LBGpIkqgRAJUcT2KTNNkwLXR7HthIvPvvBOMUn8YCyhF+mXd/gA/PvQQJCbC+vXe\nDk9EqoCm8EQCUMHpuOXYbDNo3rwJWzfv4LtjZ/Me1/JucAqtW0cQFnYOTmccoL5GuUrsl9X5Kvj+\ne/jXv5R8ivgZTeGJyGmdnI4LBz7A7Xaxbl1dHjx2CycIJ4k+ZGfXZOnS/nkdyQHmzPEkCtHRKQGb\nPJ3K3DZtWL9wEeMjri7SeVxdyUWqGcuyKuXLc9Mi4ktSU1Ot6OheVmhoCwtmWNAr51/LOo9brN3U\ns1qxssBxz9cMKzq6l7fD9xmpqamW3d445zGaYdntja2xY8dadntj62qetbZxtnVu7YZWampqiefn\nXiYiviEnbyl1nqMRKJEAkb9ux7N32xOc3OPNYhp7mcoJVvEz+fd+k6JyO7HnH41LT19CVtZEfmIM\nc+hL4pHL84rsVYAvUv2oE7lIgCj4Ju5x1lkzOXx4CN3dv9CSLfQLCqJ9m3eBGqxYMZzsbM95dns8\nkZGDiInpDainERTtxJ4/IXqG8aykGf/b2w7TNFm8eBnQrcD1MzJ26fEU8WNKoEQCVls6derE8Mfu\no12fu3mp9RV8Ni6hQCft3KQgMnIQ48a9lldAvmBBrGqgCjm5Dc5yDrCQIcZxJq9ZxuU97uPAkXuB\nYXnnBgcPZsWKILKzBwB6PEX8kVbhiQSIEhthLlgAa9fCrFklXrfEVWdzP6/8wP3IuHHjGDXKhds9\nGbD4F3Gk0YuX+QgwgURCQ/+kadPz8zYm9tDjKeJtZV2FpxEokWou/0hSQsIg0tNTAHA6k3G0aAH9\n+sGvv57yusVNQUlR6elLcpKnWMDkSUJZyDfMZAc7cQA76dgxxctRikhFUAIlUo0VHnVasCA+J4la\ngsuVRPv9W2n01FNw/vlF9m2D/H2fmuMpOvew2+NxOpOr+tfxIyYQy2ruYTrTmMAd9GdAgcfNM93n\nOVuPp4j/0RSeSDVlmiZ9+z5GZuZITk4VDcNmm47bPZmuLMVlTGVjyhdYQUFFpvdatWpVaJppGKGh\nX9CxY4SKnktwMmFtDgwEYjmLOaziLh6o14ihn75TbI2ZHk8R79NmwiJS7Bu5RydgILW4i5VczkB6\nQ7Rn65HCNU6hoWMKJV+q0ymN4hLXfjzMqHqfcMm+PdpsWMRHqQZKRAp1Gj/ZtsBmW4PbDUN5mV/5\nG2m0IZri925r2jScrKx4TTOVkcPhYObM13MSWM+x2bW/4LULmsD06fDQQ94NUEQqhBIokWrNASST\nu/pr6NAhvD12GEOPHOFqRhEcPJiMjAjgOMHBBfs+TZjgSZZOTjNpmX1p5TbaPPnYvUeDhg2hSxe4\n/XY4+2wvRygiZ0pTeCLVUIktCxwOtnTpwk8r1jD+nMasWLGa7OyXAE9vovybBytZqgQPPghhYTBx\norcjEZFCVAMlIkAJRcrLl0PnzvDHH8Tc9ZB6O1W17dvJbtWKByOuZZe9rhJVER+iGigRAYpuNYJl\ngdMJI0dCgwbeCyyAjXv3XY4fPMYtC/ZzD3erA7mIH9NyEJFAkZoKmzbBwIGAZ1TKbo/HUyOVnFMk\nHufVEKsz0zQZNWoyL/EKUWzmSi4jK+se+vZ9jJiY3pim6e0QRaQMNIUnEgiOH4d27Ty1N7fdlndY\nvYiqjmc7nO3AQB7gOLG8QiQ7ABegGjQRb1MNlIgUlZQEH30E8+aBUerXB6lAngSqOfABNiawlMd5\njgf5glfxdC6/B5gEFCz6F5GqoQRKRAo6fBguuQTmzIErr/R2NAHr5MrIe4CFRLOCqdShNVs4Th88\new2qoF/EW8qaQKkGSqSaME2TmJjeRetpXnsNrrlGyZOX5faGio7eQHR0EyLHxrPFtpcHGABs93Z4\nIlJGpx2BMgxjOvAPYLdlWW1zjr0EdAWygXVAf8uy9he6nkagRKpIiX2frrrKM/r044+Ymzap3snH\nLHr1VS4aNpy2waFkHD6CZb0CaApPxBsqfArPMIwbgEPAe/kSqGhgnmVZbsMwXgCwLOvpQtdTAiVS\nRTz1NUWngN5sUJNVC3/m2UbnFWiaqTdo32CaJllduvNvd3de4gJsthlERLRhwoQR+r8RqWIVPoVn\nWdaPwN5Cx9Isy3Ln/PgTcH6ZohSRcilxmq44OzYT8smnxG17gqVLT+QkT7GAZ6QqdzRKvMflSmKE\n+zmG8T0hPIvb7SIs7BwlTyJ+oCJqoB4A/lkBtyMip5A7TZeW1o20tG707Bmbl0QV19Mpbsc2puNg\nG0OAJt4MXU5hFU34mq4My1mBJyL+4Yw6kRuGkQBkW5Y1s7jLExMT876PiooiKirqTO5OJKC5XEk5\nNU6eabqsLM+x3I7j+TevHdlnPO0GPkJzuuZcOw7PMnkPT9PM5Kr9BaQIpzOOBQtiScx6iqWM4v9q\n23E6P/B2WCIBYf78+cyfP7/c1y9VGwPDMJoBX+XWQOUcux8YAHS2LOtIMddRDZRIBSqpzqnwUnfT\nNKkV+xC/22w492Rps2Afl9vMdOAfy+j4t9Y0/fLLAsdBRf8iVaFS+kAVTqAMw7gFT/vcSMuyMkq4\njhIokQpU4kq7fG+spmnyfPe+fHbUzcVM5HhwvJImf7F7N1x2GSxejPnHH/n+r5eruFykClTGKrxZ\nQCQQBuwCngNGAMFAZs5piyzLerTQ9ZRAiVSw041KxMT0ZnDaWr4hjmk8hhoy+plRo2DzZmK2H8wZ\nbQzHM+JYctIsIhWjrAnUaWugLMu6u5jD08sUlYhUiNx6p5Jctm8PrdlGDwZUYVRSYZxOuPhiLmwZ\nkXMgCU/yVLTuTUS8S53IRaqRkScO8mJQNseYRe5qPKczztthSWmFhMDw4YyteTRnVaU6lIv4KiVQ\nItXFggWE7d1L9zmziI5OITo6RdM9/ujRRwlfvZq0lxNp374GNtsQ8renUEIs4hu0mbBIddG5M/Tr\nBw884O1I5Ey5XLBoEXz2mVbjiVSRSlmFV85AlECJVJX58+Ghh+B//4OgIG9HI2fqr7+gZUswTWjX\nztvRiAQTyouGAAAgAElEQVSECt/KRUR8nGXBc895VnApefJ7pmkS0/M+3qrfkF2PPnr6K4iIVyiB\nEvF3330HO3dC377ejkTOUP7teoasfhz3wkX8e9o0b4clIsVQAiXip0zTJCa6F0u79yThqI2YLned\nfoNh8WkjRkzI264nizgmchdW4ligjBtJi0ilO6O98ETEO3JHKq7Pup7aHOaFv57CvclGenofdR73\nU6ZpsmzZ7wWOvUVD4jN2c3ervzF7w/a8bXkWLIjVCksRL1MRuYgf8uyLdxuLGMpk+vAJ0wATz4bB\nkwB1rfY3nv/T5sAHeJpnLgfeZhC30Zl/0oPJnG4fRBEpPxWRi1QzJU3dRLOC+hznU67MOZKEJ3mK\nBTz7qOUufxd/0RZPz6cU4GPgVZL4PzqSRQc2FnsNTe2JeIem8ER8WOENhHOnbpzOOOrOu41x7mgs\nnsLzWUhdq/2Z0xnHggW5/9fdsNnScbvhKLV5gQd4jvF05yCwEJttDZGRQ0r8+9Coo0jl0xSeiA/z\nTOt0o8jUTcIT/NW3L1GNmrFq7QYMI4hGjULYsmVnXp2MpvD8T/6mmZGRHRg37jWysiZSi2zW8zD/\noC6/MhXw/P+2atWKpUv7o6k9kTNX4ZsJi4gPGjOG9X36sOKND/NGH44fj2fUqCdJT08BwOlU8uRv\nCm8WfcUVV+QlVLNWNmfEtnO4K9/Gwps2jfFKnCKiBErEp52c1gFYjs02g5CVIWz/cxt/X7Ihb8k7\neN5Q09M1+lCd5E+ounXuztvb5tOKt1lFKrCdBg1qkZUVn/P3Qc5eecneC1gkgKiIXMSHORwO5sxJ\npn37t7HZpuN238/927YwOvsudu9t5O3wpAo99tSjTK2RzTMMAroBA9my5U8SEgZp82gRL9AIlIiP\nyl8PAzVxuyfTgWTaUY9eJAHz8bQt8NDoQ/XmcDgY17ojT/62jIu4nvW0IDtbo44i3qIESsQHFV5d\nZbM5AXiWVbzErWRTC3AAsYSGjqFjxwjVPAWA2o0bM42beZoXiON24E0WL/4T0zT1fy9SxbQKT8QH\nnVx9F46nv9NK2rINEzctqE0WLkAr7QKNaZo82ONelh05SHvsbGEyoL8DkYpQ1lV4SqBEfFDRrtQw\ni4fY2eRC3mt8LlBT27UEKNM0+aNHH4wjHXmCb3OOqn2ByJlSAiVSDZimSZcu/XC7XUAsl7KKH7iK\nh266iZTvvvR2eOJld0b+gzd/SOdy3mYXnwHbad++BkuWLPB2aCJ+S1u5iFQDDoeDiIg2eT+PYAKv\nEs2RmipbFHjwmSeYaTvGMB4kd0XeihWrtZWLSBXSCJSIj8otJA/PcvIzo2lT207yF+9ryk4AuKVN\nJ2auWM4lbGQPYWgaT+TMaARKpJrI7QE19bxpmM3PU/IkBbibNOFTrmQwr3g7FJGApBEoEV+2fTu0\naQNr1sA553g7GvEhpmkypHs/fjh6lJZMJNs+moSEQaSnLwHQAgORMlIRuUh1MmwYnDgBkyd7OxLx\nQaZpYusfx9Y6Z7G9f9+8zYdBrQ1EykoJlEh1kZkJLVvCsmVwwQXejkZ81cqVcNNN3NbmGr7+rhe5\neyOqJkqkbFQDJeLnTNMkJqY3yVddz9Yrr1TyJKd2+eVw7bU4tm3ydiQiAUUJlIgPyV15tzDNwS3r\nttE1fbGWpsvpPf00D+7dQb3aTwHJQDLBwYPJyNhDTExv/Q2JVAJN4Yn4kNwtXJ5gHzfyA7fTTdMw\nUip7/vY3Xtl3mGkH3TRoUIctW3aSnf0SoHookdLQFJ6InwviOE5cTGCEt0MRP2GaJrH/20SfTX+x\nNzOBDRu25yRPsYBnU2qXK8nbYYpUK0qgRHyI0xnH/UFD+YP6LGYFdns8Tmect8MSH+dyJfFN9mSO\nEM5thOJ2X+ztkESqPe0LIeJDHNHRXBcewuhGDYgOTcHp1LSLlJbBBEYwggmkcB022xDcbs8lnkQ8\n2bvhiVQzqoES8SWzZ8MLL8BPP4FR6ql4CXC5iw+OZk1gJc/wRPBRbhzlVFNNkTJQHygRP2WmptL0\nrr583PwSrpn4vN7wpExM08TlSuKWbRu5rzaELV7s7ZBE/IoSKBE/ZJomU7r1wZVtpzXjqG0foVVT\nUj7Z2dCiBXz5JXTo4O1oRPyGVuGJ+CGXK4mh2U2YyHgs+mvVlJRfcDAMHeqZChaRSqMESsQHXLI/\nk0vYyUz6ejsUqQ4GDIDvv/dsQi0ilUIJlIgPGBl0jClBxzjGLCBZ7Quk3EzTJKZXLO/XD2PrE094\nOxyRaksJlEgVyd3jLv/WGqZp0uvyDhj//g9zz7+QFi1chIaOoVWrVl6OVvxR7mq8tLRuDF7/OHVS\nTV6Njy/ydyciZ05F5CJVIPeNLSvrHmAhNtsa7r23K7Nmfcmb2VmspwtjuQh4G3gV0PYbUna5WwF5\nOpDDy3TEzQqG8RagvymRU1ERuYgPcrmScpKnD4DrcLsvJjk5hUbZzelOTabyDrABT/Kk7TekYriw\n6E9NGnAb+psSqVhKoESqzEIgN4kaCFyKkx1M5wb20cC7oUm14HTGYbfHA8lAMjtsm/mCDjzOVG+H\nJlLtaCsXkSrgdMYxb16/nK01JgKxnIOd+7iLNvyE5w2vOXCy6Ffbb0hZORwO5sxJzhtliowcwpQx\nk0k7+isu9nDE9iEZGW0wTVPTeCJnSDVQIlVk3LhxPPusC5gMxPI8owjnC5xnZRAcXIemTc+nd+9o\nbb8hFco0TWr06cfX+w8zxXoDWI7NNoOIiDZMmDBCf2MiOdSJXMSHjRs3jlGjXNRxj2c9T3FzrWAm\nffmh3sSkUj129c3E/7yMliRzjIfwjIKqqFwkPyVQIj7ONE02DhrOpQcyOZr8jt68pNLFxPRmeNof\nfEhtkhlE7io9SCY6OoW5cz/3ZngiPkGr8ER8nCMqiof/2kPwyBG4XEnqzyOVzumM4+Xg7cSzAgO3\nt8MRqRZOOQJlGMZ04B/Absuy2uYcCwU+BpoCG4E7LcvaV8x1NQIlUpy33+bPpCSarthCVpamUqRq\nmKmpNOl9B89nWXxuvQ7o704kvwqdwjMM4wbgEPBevgTqRSDDsqwXDcOIBxpYlvV0MddVAiWCZ8ou\nd1WUc/CDOJ58EmfIuby8+EE0lSJVas4c9j/9NHdc2BoMQwsVRPKp0Ck8y7J+BPYWOtwNz5prcv7t\nUaYIRQJI/q010tK68X6PvuwNCmJ5gzBvhyaBqHt3QmrUYO7Tj+F0xmkKWeQMlKcPVGPLsnblfL8L\naFyB8YhUK54O5J6+T2Ax9NhzTA4OwTnsYRYsjCUry3Oeej5JlbDZID6ePcOG0fOPHXlTyAsWxGoq\nT6SMzqiRpmVZlmEYJc7TJSYm5n0fFRVFVFTUmdydiN/JyNiT930McwnmOD81PJfRhRoeOp1685KK\nV2D6OHe67u67OTHwEVofGc4vOVPIWVmeZF9/gxJI5s+fz/z588t9/dO2MTAMoxnwVb4aqFVAlGVZ\nOw3DOBf43rKsIlvHqwZKAp1pmnTr1ofs7JrAJL7jBd6rsY0+33yqNyqpdCc3sC66UGFaqwjC/wii\nN7/knK0aPJGqaGOQwsnK11jgi3Lchki153IlkZ39CvABVzOD5qxnRds2Sp6kShScPi64kfDFE0dz\nHUtpxQQgGbs9nsjIDsTE9FZNlEgpnXIKzzCMWUAkEGYYxhZgFPAC8IlhGA+S08agsoMU8W8OnuYN\nJnEnZzc87O1gRIju3p01993DK9++iav1FURGDmLcuNdUEyVSBupELlJJcqdQmmc9wTwm0rp2bWZ+\n8Z7elKRKnGoKD4C9e6FlS1i6lJiHhpCW1g211ZBAVtYpvDMqIheRkjlyCsVt/eP4Z93zmTn1ZSVP\nUmUcp1uo0KABPPAAuFxeilDEv2kESqSC5V/5lHBPdyKHDIF16+Dss70cmQS6Iqvy2raFNm347o03\n6Nr/SXXGl4CmzYRFqkhxS8QLT5u8XuMRbu3VleaffOLNUEVKntKbPRsaN8a87rqiLQ9EAogSKJEq\nUNKbkcuVlFdLEsaf/EFzHr3xRj5K/6d3A5aAFxPTu/g6p2kToVMnWL8e6tXzZogiXlUVbQxEAt7J\nJeLhQApZWc0ZMWJMgXOeZAqfcCWZtexeiVGkVFq2hM6d4a23vB2JiF9RAiVSbsvxfJrvBgxk2bKV\nREZ2wG6Ppx5vMhAXLuNXMjL2qK+OeJ3TGYfdHo9nC9PknO2D4jwXPv00TJ4MR496M0QRv6IpPJFy\nME2TLl364Xa7KDwl4nTG8b8HHiZ8+3bu5m1ARbniGwrX7QF5P3+4bwsNH3oI4uK8Fp+IN6kGSqSK\ndOgQxdKl/SlSU5LyIXtCGtA5ewTLGFXwMvXVER9RuI6vc/AQUs6pTZ3Nm6GmOtxI4FENlEgVmTBh\nRPFTIu+9x7p6ISyjqbdDFClR4a1e5mVPZn1WNnyuJF+kNJRAiZRTbqPC6OgUoqNTmDMnGU6cYNvg\noSSFNiQ4eDjF1puI+KiPm10CEyZAvtkD0zS1R55IMTSFJ1JBTNPkg9vuYuCxhlxPAsHBQ2jdOoKw\nsHPUV0d8TrGtOGbPwPHUUzBxItx66+m3gxGpRrSVi4iXuCa9xcRjIYxkMtCV7GyDsDDVPYlvKnGr\nl8xMzyjUrbcWmuaDrCzP1J8SKBFN4YlUmCv27KImbv5JF2+HIlIqDoejwGo80zThzjth61ZYuNDL\n0Yn4No1AiVSQYccP4gw6gHXsfYCcuqdkL0clUrLCU3QLFsR6puiGD4cXXsDpfJwFC2LJyvKcr79p\nkZNUAyVSERYtgr59mfv660x65R1A+4mJ7ytxe5eUD6F5c5g7F3P7du2RJwFBNVAi3jB2LMTHE9Ol\nCzFdNIUnfq52bVZ36cLWmH/ganulEieRYmgESuRMLV0KXbvCunVQu7a3oxEptZJW2QHc1+NeVhzJ\n4ipGsdPu0uo7qfbUiVykqt1+O1x3Hebll2uqQ/xO4e1dHA5H3tTeWFYTSiaPco066Uu1pyk8kaq0\nYgX8+CPf3ntv8cW4SqLExzkcjhL/TqfwJH9wKc/TqoqjEvF9GoESORP33AOtWxPz/S/FF+PqE7v4\nofxTe6/yPsdrLuLyr2frA4FUa9oLT6SqrF0Lqanw2GPejkSkQpimSYcO19O372M0adKIFi1cJIWs\n5QErm6D9+70dnohP0RSeSHm98IIneapfH6czTv1yxK+Zpkm3bn3Izq4JTCIzcznwNvAqH/Mu+/o9\nwLGQEI1CieTQFJ5IKeUvtk24tweRgwfD6tWYv/yCy5VERsYuoKb2vhO/5Ckc3w4MxDMV3RvwTEs3\nYwO/0JYBUTcx+/uvvBqnSGVREblIJThZE3IPsJBuaXEE3XAlB3/5RZutSrW3keak0J4em9d5OxQR\nn6ERKJFS8Hw6bw58QGNGsJIEWhs1Ofdvf2Pp0v6oeFz8XeEpPDg5hQfQppaTxbWP88P06bz45oeA\n2nVI9aIicpFKsxCYiJOtfMAD7LSmsGnTVm8HJVIhHA4HKSkf0b79pYSGjqF9+18YO/YpoqNTiI5O\nYdKXH7K6aVO+u70faWndSEvrRs+esZ4NiEUCkEagRErBNE26dOlHA3ciqxlBBNezlX20aHGA7dv3\naApPqj3TNHHeehfzLIuL2M5h6qIRV6lONAIlUgkcDgejRw9hMEP5jGNspQ8wkC1b/iQhYVDep3Ql\nT1JduVxJrLAu4wcuYSBvejscEa/TCJRIae3fz95zwuh4YjwbGJ5zUJ/AJTB06BDF0qVX0I53+Rdu\nWvAi2bZ4/vnPWfrQINWCVuGJVLDc9gV9NqyiUf1z2LC3kbdDEqlSpmmyYsUyYAW/0Z//Mp0BDKbx\n6GeUPEnA0giUyCnkti8wskaznuHEBFmsMmqTnf0SoJonCQy5mwtDOJBER1bzddBqwg8egFq1vB2e\nSIVQDZRIBXK5ksjKmsij7OcHYvjt2Gu0bn2Jap4kQDmAz1nMMLbUbwAzZuRdYpomMTG9iYnprZV5\nEhA0hSdyGnU4ihMXf+dbYDFhYY1V8yQBpbititwJ8Z7tjB54APO77wo0lF2wIFYfLqTa0xSeyCmY\npskPXW+n3fHL6cOjmrKTgJV/K6O8BpqdO8M99xAz6+ucKT41lBX/pSJykQrkuO46os6qyeOX1Sf6\nrBScTiVPEpgcDkfRv/2RI2HAAGxN23gnKBEvUgIlcipTp1LL4eDtjz7ydiQiPsU0TVyT3sKVsZdH\nImrwgz2+wBSf05ns3QBFKpmm8CRgFTslkd/Bg9CiBcyfD5dfXvUBivio/Jtr38QnvMU2bg5vSsNz\nLyQs7BztkSd+SavwREoh9w3glHt6TZ3qqfFQ8iRSgGd16j3AO3xPFlu5hL/vvJnly38lI2MPLleS\nVuJJtacRKAlIJ/valFD0euAAtGwJ6elw2WXeClPEJ3meP9tzfhrI9VzEe9zBpRznGC5APdLE/2gE\nSqQiTJ3K9jZtiHnyWfW1ESnE6YzDZluT9/MCbmANFvdzG54PJZ6WBrlT5CLVkYrIJSAV19cmMnIQ\nMTG9qXP8GB/+NI8uJ+wsO+r5NK2+NiIn5W6uPXLkC1jWMABG0ZBPSCGZo2Sj7uRS/WkKTwKWaZqM\nGDGGTZt20qBBHbZs2Ul29n0k8DGt2M29JKG+NiIlK/wcem39H3xt3ck0/q4pPPE76gMlUgarVq0l\nK2simZlvAvdRn/d4kiNch+qeRE6ncG+oRVOn8vyw4Wy67iCDnlLyJNWbEigJWLn73Hk2SP0TWMgT\nXM+/OMAajgCD885VXxuR0+v0+OPw7bd8HRWFiafYHEpoEyLi5zSFJwHLs5KoOfABcA9n8w6rOcp1\n1GENLmA5MJ327dsxYcIIvQGIlMayZRy96SbOywpiz5EXAa3IE/+gVXgipeRZSTQDmAhMYhgRfMmJ\nnOQpFpgEPMCmTVvV10aktCIi+CmoDv2PRKEVeVKdlTuBMgxjhGEYKwzDWG4YxkzDMLTsQvyKw+Eg\nIsKzh1cjdjGQ5YymWb4zTCCZzMyRJTfbFJEi3m9xGcNI5SwOejsUkUpTrgTKMIxmwACgg2VZbYEa\nQJ+KC0ukakyYMAK7PZ4E7uF9OrIreDfBwcOBZCARzyiUPkWLlJZpmiw+ksV3ZDGI/kByTg1hnLdD\nE6lQ5R2BOgAcA+oYhlETqANsq7CoRCqBaZrExPQu0BjT4XDwrzdf5P6gH/j3jUGkpHxESsr7REen\nEBr6p5cjFvEvuVskLV06gETuYQizuaHt23n1T8U9B0X8VbmLyA3DiANcQBZgWpZ1b6HLVUQuPuPk\n5qcTgUJFrQ88AE2awNixpb+OiBRReIuk/+NG7M1303f9Kj2fxOdVSR8owzBa4Fnj3QzYD3xqGEY/\ny7I+zH9eYmJi3vdRUVFERUWV5+5EztjJlgWeF/asLM8xR9Om8NVXsGZNkes4HA7mzEnOm7ZzOvVi\nL1IWifRg+eZ4+kR2Ie331cU/B/WcEi+ZP38+8+fPL/f1y9sH6grg35Zl7QEwDGM2cC1QYgIl4pNG\njoRhw+Dss4u9uHCjQBEpWeEtknYHj+Hd40FE/nCEj2no3eBECik8sPP888+X6frlmsIzDCMCT7J0\nJXAEmAH8bFnW6/nO0RSe+Izipg/mvfgsncaP94w+1a3r5QhFqgfTNPNGbTMy9rBp6R38wXN04kXW\nEo9nYYam8MT3lHUK70xqoJ7CMxbrBpYAD1mWdSzf5UqgxKfkf2F3OuNwvPIKdO0Kjz3m5chEqqfc\nmqgRbKMdv3E35xMa+gUdO0aoO7n4nCpLoEoRiBIo8UmmaTI3YRzDVyzmt08/JqZrV2+HJFIt5Y78\nGlmjWUM8t9ey8dyXM5U4iU9SAiUBrcgoU6EXatM06dnjPuYfqcsUoplj/1LTCCKVKPc52XXLevrW\nsRG2eLG3QxIplhIoCVilWSYdE9ObBmmNeYqfuJL/YuHp+TR37ufeClskMBw7BpddBklJcPPN3o5G\npAjthScBq2CrguI7h+/f/Sfj+YzhvISlP3+RKmN+9x3j64SyquftmKmp3g5H5IzpHUQChmmaXP/7\nL6xiH9+zBUgmOHi4tpgQqWS5o8PPLn+Ewwfqk9y9jzqRi98rbx8oEZ9TuAeNZ/+t5LzL33xhKm+e\nCKIzk4EUYDutW1+i+ieRSpZ/dHg4F/JWdl+efOkNPffEr2kESqqN3M7h0dEpREenFKl/umvjar7i\nb6zgYeBzYCBhYY29Fq9IIPqOzqziXG5etUz74olfUxG5BIbNm8lu04bLj9Vi3RE18hOpSoUXeLQL\nGkTasUNcxmtkcpaei+ITtApPJJ/cJdTDf/+FiyKvZe3995+yzYGIVJz8bUUiIzuQnr4E8HQof3Bp\nDY7RliHcCiQSGvonM2e+ruekeI0SKJEc48aNY9QoF23dg0nFRUTtWrz3xft6gRapAqdqKxIT05ul\naVGsZCTXUZM1uIqcI1LVlECJ4Hnx7tKlH273JNL4gDn0ZBpnqeeTSBXJ3cbF01YEIDnv+ZebXA3K\nCuYaGtKLxUXOEalq6gMlgmfVj9t9Md1ZQjg7eYuHvR2SiHByWq9Vq5ZMPwvas4kbSfd2WCJlpgRK\nqq1grsbFNAbThRN8iM02RD2fRKqI0xmH3R4PJAPJ2O3xREZ2oGfPWNLSurF06QD2HcnkGbJwcT8G\nTmw2JxkZe7QqT/yCEiipNkzTzFsWHRnZgeE13+F3mjGPH7HZnIwe7VRthUgVKa6tSHr6kny7BYRz\n/LidWQzkBDvpyzTcbhdLl/anZ89YJVHi81QDJdVC4YLV5rWH8bvtMI+2v57tdc7SijsRH1CwLqo3\n4Pm+E1F8wm+0YjN/cRaqhRJvUBG5BKTCBavTuYE6zXZz14Y/vBuYiOQp+EHnTWAgucnUe+xgKxfw\nDMeB7bRocYCLLmoFqOWIVA0VkUu1l3+qrrhh/iv4Lw5+Z1bzS70QnYiUJP+0Xvv2NQgOHo6nRqo5\nT/ErA/iUllwFXMe6dVtJS+tGWlq3vCm90z33RaqSRqDEr5TUWwbIOf4CCxnHe0G76fnVJ/rUKuLD\n8jfaXL9+PT3WXcpNHKQrtcmd3vNIpn37t1m1am2xfaVEKoKm8KRaO11vmSXDRtJ70xo2fDwTx623\nejNUESmDmJjezE/rwm+8xDBC+YaHyf88Dw0dQ2bmSIp77otUBE3hSbWWkbGnxMsc113HiH072Pvc\ns7gm/5+G+UX8iNMZR017Ak9wG1NYSS0Gkb8FQtOm53s7RJECano7AJHSMk2TFSuWAcPyjgUHD8fp\nfN/zw/PPs61VK24a+VLeMP+CBbEa5hfxA7n1US5XEjt/rUNqRGvGGykAOJ35p+k959vt8XnHRbxB\nU3jiN05O34UDScB22revwZIlC2D5cujcmTsvv5JP0++kcO1EWFhjQKt5RPzChg1w5ZWkv/IK496b\nA5DXBFebgUtlKesUnkagxA85cr6SCQtLAbcbHnkERo9m3+y0QucuZ9mylbjdAwCNSIn4hebNWXvL\nLfx5/0OknXgLOPncVc2T+ArVQInfKG5rCKczjuXDh7Nq+Upu+XwukZEdCpxjs83A7Z6MZ0TKs3ov\n9xOsiPiuITsO0uFEff5OE/TcFV+kESjxG/lrJMBTFxF04ACNXn6FWxnF0m+b8cPCeBISBpGe7qmd\nyMhow9Kl3oxaRMoit7XB4l9X8Bj38iaxtOVKsthNRkYNb4cnkkc1UOLX/nl+c9Zua8mT5E7dFVza\nXFLfKE3hifiegs/X5cAbzOI464khgdsJDh5OSsr7ev5KpVAbAwkc//kPV2XsZCS9SjyluA1N9eIr\n4ptcrqR8mw1PAi5iMC/yED/Rhg5kZ7+kaTzxGZrCE/90/DgMHMiP3W7j0OcJ4K4DFL+02eFwKGkS\n8UtnsYuzGckY3uJhrmeAtwMSyaMESvzTlCnsMQz6fZWO2/0A8CY22xoSEoYoWRLxU05nHAsWnOz1\nFBy8ChjO29kTuY8/GRT0JF2cn3o1RpFcqoES/7NuHVx9Nfe1uoL3F96NtnYQqT7y74+Xv/dTs0P7\nmbryF4L/9z8491xvhijVlPpASbVT4AV16AAcL74II0aw0/y3lyMTkYpW3JR73s/PPANPPgmffFIk\n0XI4HMUeE6ksGoESn1Z4Fd0jQYOY0PxcQlaswJw3TyvsRAJJVha0bcuSe+/l+olvFHjuJyQMYty4\n1/R6IOWmESipNkzTpG/fx/JW5ZzLdhKPwVOhF/BWTc+fbqtWLdm0aQxNm57PhAl6sRSprnJHlyLO\nbsKwceOpdWwyWTnT91lZ8PLLY/Kt4PMcc7mS9JoglUYJlPgc0zQZMWJMzhYsl+YctZjGo7zJTWyo\nV7PIyFRWVrz3AhaRSpX/+Z4GNOdnXmYWD/Cot0OTAKYESnzKuHHjGDXKlZM4TcazcXAsd/ITl/Bf\n7q99nI+d7xXqF6NPmyLVWeHn+1MsZjmvcytD+RcR2O3xDB06iHHj4vNW8BXX0kSkIqmRpvgM0zQZ\nNWpyzt51TXKOOgjnZV7l/3iyvsHHX7ynJEkkwP1FR15oeRnJtd6kR9RnzJmTTEJCgprmSpVSEbn4\njJiY3qSlbQcGkjvyBC+QwiusqLmWE4kjSE9fAkBkZAcVjIoEiMJbvNhsM4iIaMMXTc7CduwYD1hn\nAVp5J2emrEXkSqDE6/I2D128jMzMHsAHgOeFsj9TedLmpk+zS1i/dTfZ2S8BJ1fd5CZUeuEUqd4K\n1kZOBiA06EkWHzvIIwwmlXZFPkiprYGUhRIo8StFNw99GxgALKQpy/kvh7mZ0fzON3hGptQ0UyRQ\neUapu3HydaATkUTyIe8TwTL28E3e64I2EpeyUhsD8SsjRkwoUBwKEBr6BVd0aMfY/6xi0qHb+Z1n\ngaVei1FEfFc6l/Eh/Xib7vTiBIsXZ+SNPGmhiVQmFZGL15imybJlvxc62paOHSMwu9xAnRo2JnFr\nzqwInaUAABsRSURBVPE4YBiQDCTnrLCJq9J4RcS7nM447PZ4cl8HgoNXERw8nJFk05yf6E9rMjNH\n0rNnLBkZe7wdrlRzGoESr3G5knC77wdO9nCy2YaQ2H0MJCayZ8rL1HoknqwsT54fHHyc1q3fJSzs\nHJxODcWLBBqHw8GcOcn56po+AqBv38folzma+UwmnWdYnzUReBu7/fRtDVQnJeWlGijxmpP1DOFA\nErCdTu3g39n7YORI6NtXL24iclq5ryVPsI8+fMQNDODm6G9wOuNO+fqhOinJT0Xk4tPyJ0TFtSL4\nI6oDF4SFwXvveTNMEfEjuYnQkawJpOJicc2NdPj689MmQkWL0rU4JZCpiFx8VuFPewsW5LYiSAHg\nxasHcMGsWfDxx94MU0T8TP6pvXeONiX59x3UDgrydlhSzWkESqrMKT/tbdkCV1wBX30FV13lzTBF\nxN99+y3ExsKSJdC4cYmnaQpP8ivrCJRW4Yn3ZWfDnXeC06nkSUT+v717j4+yuvM4/jkhpESLSAgF\nWUApWwRZG4JWqegmRZMI26JCtYIC4oV166JgoEiplEIsVaSirl2KyMULFlaLpdXOGKlBbat1JUsR\npV7KTV3ASL2sjUac3/4xM8lMMknmmUkymeT7fr3yMpnnmec5mZc5/J5zfud3knfeeXDllTBlCgQC\njZ4WHrnS9i+SCI1ASZtp9Gnvt7+Fv/4VHnsMMhTTi0gLOHoUxozhtcGD+fe3PwTiW4iihSudV5sl\nkTvnjgdWA8MBA640s+cjjiuAkgYadE4ffQTf+x689BL07Jni1olIR1Lx4IMMmzKNi5nHs5xc+9AG\nxAySNKXXubVlALUe2GZma5xzmcCxZvZBxHEFUB1ESz+Rha/3Dx9/xM93/YmsrVvhtNNaoqki0slF\n9ldVVe/Ru/J01vAwX+NF/pdy8vPvZffuN2IGSVqV17m1ySo851wP4BwzmwZgZkeBD5p+l6Sjhivn\npiX1RBa+HtWL+SNlzO0aYFxVFXq+ExGv6j/cAVH9VUZGKQGm8zN68gjfppCr2LfvYMQWL36qqwcx\nefJ1bNhwT4p+C0lbZub5CxgBvACsBbYT3AH2mHrnmKS/oqIJBusMLPS1zoqKJjT5Hp/PZ0VFE6yo\naIL5fL6o13NyBhustQeZbA9wmcHa2us19j4Rkfp8Pp9lZ/cJ9U/rLDu7j+XnF0T0Vz6DoQY9zHGj\nbaan3cMXbPDgEaFzfAbR7y8rK2twTfVFnUcobok7Fkq0DlQmMBL4dzN70Tm3ArgJWBh50qJFi2q/\nLywspLCwMMHbSbqINWK1YMFMHn30t+zY8QqBwMnM4bcM5Q3O5jlgU6PvU+6BiDQm1mbB+/YtCR31\nh16/Ffg1xmqm8RP+xGK27/0LB7LmUlMzOHS87v3btm2pt1WM+qCOrKKigoqKisQv4CXasrrRpb7A\nnoifzwZ+U++c1g8XpdXFespr6oksesTKZ9Df4DiDUQbrrIQye5sM68/yqOslMtIlIp1XrD4jP390\nqL8aFXGs7rxh7LLDdLfJJ+eFRsPV50gdPI5AJbRm3MwOAgecc0NCL50H7Eo8jJP2KvE6KX7gUuDv\nwF1AP77CQdZzF5dwG2/xX+TkLNEok4gkpLR0BtnZ84D1wHqys+exdOnNbN68npycd2O+51VOYQZX\nsGLvazx6148avD+cRyUSj2RW4eURLGOQBbwJTDetwuv06qbiBkW8ei3d6c7zXMIdTGE1hQ2WB2v5\nsIh41dgK4ej+ZCfBNN27gGDf8j/fGcuQnTspv/lmlt1zf4P3S+ekzYSlVcVT0sDv9zN58nUcOdIb\nGE0mD/AbvsTrZDKTfeTnf5WlS+fH3BldBexEpCXU37h827btQKhvKS6GK66Ajz+GTZtUwFcABVDS\nihobJYKGRen8fj/jx19KTU0X7uUE+vIWF+FYVFbKggULUvY7iEjnFRlUzZl5BcW33gqFhVBWltqG\nSbugAEpaTawic/WL0mVlzWL48Dxyc3tRUDCSHves4pzDB/nu8DNYeNsPokaVNOIkIm0l1gPgb9be\nyZj582HhwuCIlHRqbVJIUySsflG6mppMKiunA9Cv4nr+s2c3svfv4ff9+kW9T2ULRKQtxSp78JP7\nNjHmiSeCo1B9+sDYsSlto6QXTfxK3GKtejnxxP4RZ6wCbgemUcCJ3PaZccOXT4V6wRPU78yCgVR4\nNEpEpM0MHQqbN8PUqfDCC6lujaQRBVASt3BJg/z8e8nJWcLQoUOZOLEoIqh6B4CRvMQmLuFS/o29\n3XuktM0iIhD7AbC2bMHXvw5r1/Lp2LFMH11McfFE/H5/KpsraUBTeOJZOOfpyBHYvXseCxbMZNu2\nLVRVdeHoy7P5zWcBZjCV57PXs7l0fcxrlJbO4LnnplFdHfw52JnFPldEJFnhB8DGqoz7u3Zl8/99\nzsI/VDKaBVzUimkFyv/sGJRELp40uVv5/v1Uf+1r/EfvAZT3O7HZjkGdiIi0tcb6nXDfdhPvcDkP\nUsi15BdVBPu2Fr6/at61T0oil6QlFNgcPgxFRWTPn8/cWbOYG8d9SkpK1GmISJuJZ/HKT5jPMfyd\ncpax5LO8pB70Yr03VjL78uWr1BemIQVQEqW5DibW1Nv8q+6EMWNg0iSYNStVTRcRaVJTwUtk37aQ\nwXTPfI+7X69k5IVTOPTJMsDbauHG+lLpOBRASZTmno7q5xHMv+pOvrFkCUyYAD/8YaqaLSKSlPp9\n27AbH+XFGdfxyCfHcD4T+ZgvehotaqwvVf5nx6EASjwLT709vXEjJ109gwd796P3WWdR4uKeOhYR\naXPNBS/hvi089fbC32q4nS/j43z+hcf5sAXa0Fwyu6QPJZFLlPrDzpGVxSPn/5/euJHeky5ns41l\nIRPIzr4ptBovYr8pdQoi0s7Eykuqv2/eLbfcTXX15cC9OLpyN3l8jT2M7/oh63/9ENBw+6pY91Gy\neHrRVi6SFL/fz/z5S9i37yA9ex7DgQMHqakJzv/XdgDDhnFg2HDW/P1cFrEZcMAcMjLWEAjcEX2u\nOgsRacfqBzoZGaUEAsuBLcB4oC/wc5bxPN/K/pC3197HN6ffEFdgpJXG6UUBlCSs8Y4kXLJgDmf2\n2MSvPq3igV79mPv2zRHHvg5cS8zyBiIi7VTD0izhviwcQIVfX8fP+95MweFDfCNwG/9LeMGM+rqO\nQmUMxLPwU9JLL+2ISnoMBFZGnsVI7mPzB8b3uIxfvPsIWVlzqakJHs3IeJ1AoM2bLiLSwkaTkTGb\nQOBKYE7tq1lZs5l5pCuzAn35A2WMpYTdDEtdMyXlFEB1YuHpuh07XglNvb1T74xwRwJjuJWH+ZwZ\n3M+vuBBqziY//15yc7cAUFAwm1tumaeVJSKSVhomlj/IggWlbNu2naqqk4G15Ob2oqoqj8rK6dxG\nXw5yCRWMYgIzqcxerb6uk1IA1UnVTdcNAu4gOOrUF7i89pxwR9Jtw11MefU1Lra5PMOFtcdzc/tE\nDVuffvrpWlkiImmlsVVxCxZEn1dcPDH8Du5nE4eYxa/cUg7M/QH56us6JeVAdVJ18/715/nnkJPz\nGKedlkfp7Kspee452LiRZ2+6iZLrf6AVJSLSKcVaVbf11gV8felSXr3gAm544xA4p2TxNKYcKPFo\nBnXBU3DUacOG9ZQUFMAVV8CBA/DHP3JO795sHjBAI0wi0imFR6rCq5RPPHEoHw4ZwjNLl9Jz+tVc\nYmfxXaZ6qlYu6U0jUGmmuWWx8S6bjX6a2klGxjry8v6JpUvnUzJsWLCy+Fe+AmvXQrdurfo7iYi0\nZ+F+tarqELt2vRZV2qVfvxM49Oa/8gA+evEeE7mEEU1sQqzSBu2X1xEozKxVvoKXlpbk8/ksO7uP\nwTqDdZad3cd8Pl/cx2Ndr6hoghUVTag7b+tWs759zZYtMwsEWvtXEhFp16L71VGh/1roq9Sgp8E6\nc3xuS1hge+hl1545Jo5rNd9HS9sKxS3xxzleTvZ0YQVQSasf4BQVTaj3x7vOioom1J7f3PEmBQJm\nt99u1qeP2VNPtdJvJCKSXqL71fp97KhQEFUXFH3HHWuf9uhhtnJlg4fQxvromA+z0ua8BlDKgWqn\nYu3kPXTo0Na52ZEjcPXVsH8/vPACnHhi69xHRCStzSBypXKw/t21wHpgFfAOr40YQdbD98HFF8Mz\nz/DUxRdz288eAKCq6r0GV6yqOtSgr1cOVXrISHUDJLbonbzDf1xHyc6eR/CPdX2o1tKM2veUls5o\n8nhMzz4LI0bASSfB738fFTz5/X6KiydSXDwRv98f8+3xnCMikq6i+9WDZGUdJT9/LUVFW1i8eHbo\n2EFgPNnZe1i69GY4+WR4/nneqqqi/4Rvc7A8n/Ly8ezatYOsrLlE9tGQGdHX96W6ehCTJ1+n/jQd\neBmu8vKFpvCSkuhQr8/ns/z80ZaTM9jy8wusrKws9vlHj5otWmSf9Oxp3x9xVoPj8czVaz5fRDqD\npvrdWMfCr+XkDLYpXGOHybU53GYZrLH8/NFR5+fnF4T6UF/UVKD607aHcqA6hkSDk+j3lRoc1/Aa\nr75qduaZVjVihH25W27t8ays4y0/v8CKiiZE/FE3nk+VVM6ViEgHFCvp/CTW2TZ6WQXdbdwpp0Wd\nm5V1vEFujAR19adtzWsApSm8dipcc6SoaAtFRVui5sSbmjaLnvrbA9xFeBrw0+ql7J85B845B6ZN\nY3LuIP76ye2Eh45rajKprJxOefl4dux4uU1/XxGRjiC6D14EXM9e5vANlvE4Y1n3ynZ2zZoFgQDL\nl6+ipmYF8CDwbiqbLQlQEnk7VlJSEpVI2HDvuvgTDofzMqspo1vVRzxzx+2UPbSFl7b/GbggdMYq\nIBxMQSCws3YfPIi9t13DPaS0/52ISJ0SYCAwhwDTWMZ0fJzKprU/hcpK+h/9QsR59xC9lZb603bP\ny3CVly80hZeQxuba64aFmx7mrT+Fdyzd7TbOt8N0t5ldu1vZ4sWNTPE1vG79uXov7RUR6Yzqp19k\nZPRq0LeWnHeR2YoV9ulxx9mizGOtK6sbpFF4ruEnSUM5UOmrqbynunyjyLwjn8Eoy8kZbD6fr/aP\nKj9/tOWP+Gdb9NUz7f0ePezJEwbaxf88rpFaUqWhhPPRlpXVWwmMIiIJiuqDQ4FQWVlZ4/mse/fa\noTPOsP3HfNG+P+KsBPJc1Ve3JAVQaayppOy6Y+GVGqWhxMNSg1HmXHfLzAw+6eTxI/tdRpZ9NHCg\n2dNPx30PPdWIiCSmqcCm2b718cfNhgwxGzfObPfumNcOB2bduw9Usnkr8RpAKQcqTdTlG91KcJ58\nLTCdYPLhrZitpM/RSyjjac7Hx48ClzIv53WO//Hd8OO7a/dcaipvqX7OlYiIxKcuebwvsIrq6kHM\nn7+ktl9tsm8dNw7OOw/uvhtGj2bf6NHMff8o73+hGwUFI7nllruprr4ceAb4x7b5haRZ2ky4Half\nfTw7ex6bNweDm/BGlpBJbm4vqqreo7LyU+BachnHXEZzFe+wkhu4lXl8xGIyMtbUJpuHr1VSUqLN\nLEVEWlhx8UTKywcRfqgFyMiYzRNPPOypj31640Z2XXYlkz53/CdjWO6e4327A9gCjCcYoE2rvUdk\n3y7J8bqZsEag2pFw6YK64CYYPDUWVE0ZO4kbbRPXMJONHEseXXibIcBmMjLWhYKn4Kq66mpqr6vg\nSUSkZZWWzmDr1ssIBJZTt5o52N966WeX3reJ8s9/xjIKWcQ1/MU+4E62cA81fAAEV+ytBxaRk/Mu\nGzYE/z0oLp5Y2w71621DAVQ7U3+ot7h4YkRNkWAgdP8tK3jojOHs7fYJD1Q/ST7dOcBPgJ04dyMj\nRpwK/BOVleGr+IGVPPvsazz11O8wWwFozyURkZZSUlJCXl5kv5uc/ezmSv7MUM7lJn7Nm3Th52xl\nBR/yLseRnb2nNnjSXnop4iVhyssXSiJvEZFJ36eyw9Zzln2Q2dVs9myzffsarRhel9AYmWyeo+RD\nEZEWUj85vCVWyEWXrCmtXTR0EiPsHjLtbxldbEv/QfbsypVmph0hWhKqRN6xzLl+Opd3vYHfMQwf\nBbye+WdefHgD/PSnMHAgubm9Yr4vPB2Yk/MYwQKZe4Ahbdl0EZEOK5yzWl4+nvLy8Vx0UXCWoLEd\nJOJV13e/C/yeYK7T7eylkutYzbVnF/Otq6Zy9qJFcO65nHX4HTIItPSvJ/HwEm15+UIjUMl5802z\nefPMvvQlqxoxwpaceoaNPfdCzxv6RtePCj/NhAu89VS5AhGRBCQy8uOlVIzP54tZhLP2Hp9+avbQ\nQ3ZgwADbQ4bN49vWlxWqC5UEVMYg/b0xdSq5v9jIUycM5N1Jk9j8ygEgdnJgrMTzyHPqyhZcTjDx\ncBqwkoyM11m8uFTz5CIibaD+KuvmcpVKSkpYvHg2Cxc2sqVWVhb+Xr24qKqGU5jEv1LOn/glD81b\nqH69jaiMQTvj9/u57sLLeeuTH/MpfwHuJbghcOLLVcNlCyLLIGilhohI4horO9NYvxosczCe8IIg\nCE71Pfnko83ep7GV0yNHFlJZOb32ml1Yw5iix5u9psSmMgZpbvnyVbz5SXhT34kEg6foUgSJzKkr\nWBIRaTnNjf7Ho6rqULPlBxrrv/1+Pzt2vBz12ud08XR/SY4CKBERkQR4eTitvwtEVtYsdu3qSk3N\nNYD38gPLl68iELgCmFf7WkbGbEpLH/byK0gSFEC1M9F/ZIOA62uPRc1/i4hI2qg/YlVVlRc1/ZZY\nseNTCea2rgLeIS/vFM02tCHlQLVDkXPeBQUj2bZtO6AKsyIiHUWsnKj8/HvZvfuNuPKqvOZgSfO8\n5kApgBIREWljsQKgoUOHRo1KNZdorn1NW5aSyEVERNq5WEno4e+9XENBU+oogBIREUmBWAFQZKK5\n8l7bt6Sm8JxzXYD/Bt4ys2/VO6YpPBEREQ80LZc6bZoD5Zy7ETgN6G5m4+sdUwAlIiIiacFrAJXw\nZsLOuf7AOGA1EPcNRUREJDF+v5/i4okUF0/E7/enujmdWjI5UHcAc4HjWqgtIiIi0giv++lJ60oo\ngHLOfRM4bGaVzrnCxs5btGhR7feFhYUUFjZ6qoiIiDRh+fJVoeApue29JKiiooKKioqE359QDpRz\n7sfAFOAo0I3gKNSjZjY14hzlQImIiLSQRDcklvi0eSFN51wBMEer8ERERFqPqo+3rlQV0lSkJCIi\n0opiFd9U8JQ62spFREREOr02K2MgIiIi0lkpgBIRERHxSAGUiIiIiEcKoEREREQ8UgAlIiIi4pEC\nKBERERGPFECJiIiIeKQASkRERMQjBVAiIiIiHimAEhEREfFIAZSIiIiIRwqgRERERDxSACUiIiLi\nkQIoEREREY8UQImIiIh4pABKRERExCMFUCIiIiIeKYASERER8UgBlIiIiIhHCqBEREREPFIAJSIi\nIuKRAigRERERjxRAiYiIiHikAEpERETEIwVQIiIiIh4pgBIRERHxSAGUiIiIiEcKoEREREQ8UgAl\nIiIi4pECKBERERGPFECJiIiIeKQASkRERMQjBVAiIiIiHimAEhEREfFIAZSIiIiIRwqgRERERDxS\nACUiIiLikQIoEREREY8UQImIiIh4pABKRERExCMFUCIiIiIeKYASERER8UgBlIiIiIhHCqBERERE\nPFIAJSIiIuKRAigRERERjxRAiYiIiHiUcADlnBvgnHvaObfLOfeyc+76lmxYZ1dRUZHqJqQ1fX6J\n02eXHH1+ydHnlzh9dm0rmRGoz4DZZjYcGAVc55wb1jLNEv0hJEefX+L02SVHn19y9PklTp9d20o4\ngDKzg2b2P6Hv/w94FejXUg0TERERaa9aJAfKOXcSkA+80BLXExEREWnPnJkldwHnvghUAGVm9ljE\n68ldWERERKQNmZmL99zMZG7knOsKPAo8GBk8eW2EiIiISDpJeATKOeeA9cB7Zja7RVslIiIi0o4l\nE0CdDTwD/BkIX2S+mflaqG0iIiIi7VLSOVAiIiIinU2rViJ3zi1zzr3qnNvhnPulc65Ha96vI3DO\nne+c2+2ce905Ny/V7UknKu7aMpxzXZxzlc65X6e6LenGOXe8c+6RUL/3inNuVKrblC6cc/NDf7s7\nnXMbnHNfSHWb2jPn3Brn3CHn3M6I13Kcc+XOudecc086545PZRvbs0Y+P08xS2tv5fIkMNzM8oDX\ngPmtfL+05pzrAvwHcD5wCjBJxUk9UXHXlnED8Ap1U/MSvzuBJ8xsGPBVgvXxpBmhUjjXACPN7FSg\nC3BpKtuUBtYS/Lci0k1AuZkNAbaGfpbYYn1+nmKWVg2gzKzczAKhH18A+rfm/TqAM4A3zGyvmX0G\n/AK4IMVtShsq7po851x/YBywGtBKWg9CT6vnmNkaADM7amYfpLhZ6eJDgg9AxzjnMoFjgLdT26T2\nzcyeBf5W7+XxBBd3EfrvhW3aqDQS6/PzGrO05WbCVwJPtOH90tE/AAcifn4r9Jp4pOKuCbsDmAsE\nmjtRGhgEvOucW+uc2+6cu9c5d0yqG5UOzOwIsBzYD7wDvG9mT6W2VWmpj5kdCn1/COiTysakuWZj\nlqQDqNB8684YX9+KOGcBUGNmG5K9XwenKZMWECru+ghwQ2gkSuLgnPsmcNjMKtHoUyIygZHAz8xs\nJPAxmkKJi3NuMDALOIngqPEXnXOXpbRRac6CK8T0b0oC4o1ZkiqkCWBmRc005AqCUwLnJnuvTuBt\nYEDEzwMIjkJJnJoq7irNOgsY75wbB3QDjnPO3W9mU1PcrnTxFvCWmb0Y+vkRFEDF63TgD2b2HoBz\n7pcE/398KKWtSj+HnHN9zeygc+4E4HCqG5RuvMQsrb0K73yC0wEXmNknrXmvDuK/ga84505yzmUB\n3wG2pLhNaSNU3PU+4BUzW5Hq9qQbM/u+mQ0ws0EEE3h/p+ApfmZ2EDjgnBsSeuk8YFcKm5ROdgOj\nnHPZob/j8wguZBBvtgDTQt9PA/QQ6YHXmKVV60A5514HsoAjoZf+aGbfbbUbdgDOubHACoKrUO4z\ns6UpblLaUHHXluOcKwBKzWx8qtuSTpxzeQQT8LOAN4HpSiSPj3PuewT/0Q8A24GrQ4tpJAbn3MNA\nAZBLMN9pIfArYBMwENgLXGJm76eqje1ZjM/vhwRX3cUds6iQpoiIiIhHbbkKT0RERKRDUAAlIiIi\n4pECKBERERGPFECJiIiIeKQASkRERMQjBVAiIiIiHimAEhEREfHo/wE6/K8gF5ZQ1gAAAABJRU5E\nrkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(10,7))\n", "out = mwgaussfit(x_vals, Data, eps_data) \n", "new_D = out[0][0]*np.exp(-(1.0/2.0)*((out[0][1]-x_vals)/out[0][2])**2)+out[0][3]\n", "plt.plot(x_vals, new_D, 'r', label='Fit')\n", "plt.scatter(x_vals, Data, label='Original Data')\n", "plt.legend()" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [ "\n", " " ] } ], "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.10" } }, "nbformat": 4, "nbformat_minor": 0 }