{ "cells": [ { "cell_type": "markdown", "metadata": { "toc": "true" }, "source": [ "# Table of Contents\n", "
" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# 課題解答例\n", "\n", "## Gaussian(正規分布)へのフィット\n", "\n", "正規分布で知られる,ガウス関数\n", "$$\n", "f(x)= \\frac{1}{\\sqrt{2\\pi\\sigma}}\n", "\\exp \\left(\\frac{- (x-\\mu)^2}{2\\sigma^2} \\right)\n", "$$\n", "フィットをやってみましょう.\n", "\n", "例えば,平均値($\\mu$)が60点,偏差値($\\sigma$)が15点として,ピークの人数が20人としましょう." ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAD8CAYAAABw1c+bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VPXZ//H3PVkhQFgSMKxhiSEJu2GVIqICgksrWsHdqlTFrZtibftof/XRat21RSu2LtQVrRubUhGtoAl7ICyBsIQtYUuAELLdvz9m8IkxkZCZyZnlfl3XXMw5c+bM56vhw+Q7Z84RVcUYY0z4cDkdwBhjTNOy4jfGmDBjxW+MMWHGit8YY8KMFb8xxoQZK35jjAkzVvzGGBNmrPiNMSbMWPEbY0yYiXQ6QF0SEhI0OTnZ6RjGGBM0li1btk9VExuybUAWf3JyMtnZ2U7HMMaYoCEi2xq6rU31GGNMmLHiN8aYMGPFb4wxYcaK3xhjwowVvzHGhJmTFr+IdBGRz0RknYisFZE7PevbisgnIrLJ82ebep4/XkQ2iEieiEz39QCMMcacmoa8468EfqWq6cAwYJqIpAPTgYWqmgIs9Cx/h4hEAM8B5wPpwBTPc40xxjjkpMWvqrtVdbnn/mEgF+gEXAy87NnsZeDHdTx9CJCnqltUtRx4w/M8Y8xJHNx8gK/ueIPPL32GvA/WodV2mVTjG6f0BS4RSQYGAl8DHVR1t+ehPUCHOp7SCdhRY7kAGFrPvqcCUwG6du16KrGMCSlLfvk2Lf7+OOlHvmEE1e6Vs6Egoht5Ay5l6IL/R7O2zZwNaYJagz/cFZEWwGzgLlUtqfmYuq/Y7tXbEVV9QVUzVTUzMbFB3zo2JqRotbJo1B8Y/sRPia04zBejfkfO35dQ8EU+i698np2JAxi97DHyuo2hKGev03FNEGvQO34RicJd+rNU9V3P6r0ikqSqu0UkCSis46k7gS41ljt71hljaig7VMby/tczevsbfJFyPUOXzyClRfS3j3ce6f6FeOnd79Lv0as4MGAom97+iJSf9HEutAlaDTmqR4CZQK6qPl7joQ+Aaz33rwXer+PpWUCKiHQXkWhgsud5xhiP6spq1px+CSO2v8Gi8Q8zcv1MomuUfk3DHrmErS8vJlLLaX/JmWxbmNfEaU0oaMhUz5nA1cAYEVnpuU0AHgbOE5FNwLmeZUSko4jMAVDVSuA2YD7uD4XfUtW1fhiHMUFr8U+eYHDRXBZP/iuj596DuOQHt0+/JpPKz7+iWiI4etEUyo+UN1FSEyrEPT0fWDIzM9XOzmnCQe5ry+h59XCWd7yQoTveOWnp17T07ncZ9ugkFg29h9FLH/ZjShMMRGSZqmY2ZFv75q4xDjmy5wixP5vCfld7Uhf//ZRKH9zTPot7T2XU14+w/NGFfkppQpEVvzEOWX72r+hWkcfex16jTc+2jdrHGZ8/Tn50Kh2nX03xtkM+TmhClRW/MQ7Y9F4OI9f/ncWD7mLAXaMbvZ+49nGUz3yN06p3s+KKR30X0IQ0K35jHHDg1t9xmJb0e+M+r/eVdtUZfNV1MoO/epLC1Xt8kM6EOit+Y5pYzotLGbrnfVac8xvaprTzyT47vvhHYjhO7lUP+mR/JrRZ8RvThLRaqfjNvRRJezJfu8tn+00+L4Wv0m5k+Jrn2bE432f7NaHJit+YJrT8z58w8NAi1l3yO1qc1sKn+0555fdUEcG26//Hp/s1oceK35gmotVK7IO/oyCiG8Nemurz/SdlduLrwbczYstrbP4o1+f7N6HDit+YJpLz9yVkHM1iy2X3ENMqxi+v0efl31BONDunP+OX/ZvQYMVvTBM5/L/PUEw8g5642m+vkZCWSFavKQxa+wrF24v99jomuFnxG9ME9izfxeDt77Bi0M98PrdfW8L9t9OCo6y88x9+fR0TvKz4jWkC63/xPBFU0f3RaX5/rbQrB7G65Qi6ffwc1ZXVfn89E3ys+I3xs/Ij5aR/+TzZ7SfQbUzPJnnNI9fdTnJFHssenNckr2eCixW/MX6WdffbtK/ei+uO25vsNQc/PIk9riR41j7kNd9nxW+Mn7V69Vnyo05n0D3nNdlrRjWPYv1ZNzN43zzy529sstc1waEhV+B6SUQKRSSnxro3a1yUZauIrKznuVtFZI1nOzvBvgk7eR+so++RpWybcDOuyKZ9n5X+5FQqiWDbH19u0tc1ga8hP4n/BMbXXKGql6vqAFUdgPtavO/W9USPsz3bNugCAcaEkoKHXqWSCDL+dEWTv3b7fqexInEsvb5+zT7kNd9x0uJX1cXAgboe81yP96fA6z7OZUzQq66sJiVrFisSx5HYp4MjGSouv5rOVdtZ/exiR17fBCZvf/f8EbBXVTfV87gCn4rIMhHx/XfUjQlgq57+nE5VO6iY7L8vbJ3MgP+5mBJacvivrzqWwQQeb4t/Cj/8bn+kZzrofGCaiIyqb0MRmSoi2SKSXVRU5GUsY5x3ZMarlNCSgfdf7FiG5gnNWdVrEv03vc2xA8ccy2ECS6OLX0QigUuAN+vbRlV3ev4sBN4DhvzAti+oaqaqZiYmJjY2ljEBoXRfKQM2vc2qlEtp1raZo1la3Hw1rTjMivvfdzSHCRzevOM/F1ivqgV1PSgicSLS8sR9YCyQU9e2xoSalQ+8T0uO0PLWa5yOQv87R7MzogtRb9h0j3FryOGcrwNLgFQRKRCRGzwPTabWNI+IdBSROZ7FDsCXIrIK+Ab4WFXta4QmLES9+SoFEV3pd1u9s5tNxhXpYtPgKxlYNJ+inL1OxzEBoCFH9UxR1SRVjVLVzqo607P+OlWdUWvbXao6wXN/i6r299wyVNWuCWfCwr7cIgYWLSBv8BVNfux+fTrfezWRVLHu/recjmICQGD8VBoTQnIffp9Iqki663Kno3yr10XpbIrJIP7T2U5HMQHAit8YH4udM5ttkT04/bL+Tkf5jp1DJ9G3+AuK1hY6HcU4zIrfGB86lH+Q/vsWkj9wEuISp+N8R8fbJxFBNbkP/dvpKMZhVvzG+FDOQx8STQUJP5/kdJTvSbmkL1ujehE39x2noxiHWfEb40NRH85mZ0QXMq6v9ysrjhGXsPWMS+l34DMObq7zLCwmTFjxG+Mjh3cdpv+e+eT1vSTgpnlOaH/LJKKoZO1DHzgdxTjIit8YH1nz5znEcpzWNwTeNM8JaVedQUFEN6I/tKN7wpkVvzE+IrPfYa/rNPpMHeF0lHqJS8gbMIn+hQsoKShxOo5xiBW/MT5w7MAx+u2cw4a0nxARHeF0nB/U7qZJxFDOmoc+cjqKcYgVvzE+sObJhcRRStxVP3E6ykll3DCMva7TkA9tnj9cWfEb4wNlb3/IYVrQ97aznI5yUq5IFxt7TSRjxzwqSiucjmMcYMVvjJe0WknZ+BE5ncYR3SLa6TgNEn3JBcRTTM6ML52OYhxgxW+Ml9a/voKk6l1UnX+h01EarM8vzqOMGIpnfeh0FOMAK35jvLT3xQ+pRuj9ywlOR2mwuPZxrEkcQ/KaD9FqdTqOaWJW/MZ4qf03H5LTcjgJacF15bhj51xIckUe+fM2OB3FNDErfmO8sGf5LtJLl3Fg+AVORzllPe+YCMD2v9phneGmIVfgeklECkUkp8a6+0Vkp4is9Nzq/B1XRMaLyAYRyROR6b4Mbkwg2PjkxwB0ujl45vdP6DS8Kxti+9P6S5vnDzcNecf/T2B8HeufUNUBntuc2g+KSATwHHA+kA5MEZF0b8IaE2hiFnzIjshkel2c4XSURtmdeSF9iv9rJ20LMw259OJioDE/FUOAPM8lGMuBN4CLG7EfYwLSsQPH6Lv3U7akXxiwJ2U7mcSfXei+JONjc52OYpqQN3P8t4vIas9UUJs6Hu8E7KixXOBZZ0xIyHl2Ec05RovJwTe/f0La1ZkUujrAnO/90m5CWGOL/29AD2AAsBt4zNsgIjJVRLJFJLuoqMjb3Rnjd6Wz51JKMzJuGeV0lEZzRbrYmDyO3tvnU1Ve5XQc00QaVfyquldVq1S1Gvg77mmd2nYCXWosd/asq2+fL6hqpqpmJiYG12FxJjx1y53L2sSziW0d63QUr7gmjKed7if31Wyno5gm0qjiF5GkGos/AXLq2CwLSBGR7iISDUwG7KxQJiRsW5hHckUepWed73QUr/W+YyxVuNj32jyno5gm0pDDOV8HlgCpIlIgIjcAj4jIGhFZDZwN/MKzbUcRmQOgqpXAbcB8IBd4S1XX+mkcxjSprTPcJZl8S/AXf9uUdqxrMYSELPuAN1xEnmwDVZ1Sx+qZ9Wy7C5hQY3kOYJ8amZDT/PO5bI3qRfKYnk5H8Yn9meMZtegB9m/YR7vUBKfjGD+zb+4ac4rKDpWRUfQZ29KC/93+CYnXnI8LZf3TC5yOYpqAFb8xp2jt3xbTnGM0nxQ6xZ92dSb7JAGdZ/P84cCK35hTdOTtuZQRQ8atgX/RlYZyRbrY0HUsqfnzqK6sdjqO8TMrfmNOUZe1c8lpN5rmCc2djuJb559Pohax/l/LnU5i/MyK35hTsGNxPj3KN3B0VOhM85yQesc4qhEKX7aje0KdFb8xpyB/xnwAutw4zuEkvpeQlsj65mfQNmu+01GMn1nxG3MKohctoCCiK93HpzodxS8KB4wl/fBSircXOx3F+JEVvzENVFFaQdruhWxJGRe0Z+M8mTaTxxFJFbnP/cfpKMaPrPiNaaDcl78hnhKiJox1OorfpF0/jMO0oPwjO54/lFnxG9NAB16fTxUu0m47x+kofhPdIpp1p40heeN8uwh7CLPiN6aBEpYvYF2LIbTuXtflJ0JH2Y/G0rUyn+2fbXY6ivETK35jGuDg5gOkHc1i/xmhdzRPbck/d49x6/N2dE+osuI3pgHWP7eQCKppNyV05/dP6Hp2T7ZHdif2C5vnD1VW/MY0QMWcBRQTT9q1dV1zKLSIS9h6+jjS9/yHitIKp+MYP7DiN+YktFrpmTef3I7nEBl70jOZh4ToC8bSkiOsm7nE6SjGD6z4jTmJ/Hkb6FS1g/KzQ39+/4S0aWOoJIKDb9g8fyhqyBW4XhKRQhHJqbHuURFZLyKrReQ9EWldz3O3eq7UtVJE7IKeJijtmOme6+4+9TyHkzSd+K7xrGs5jMRVnzgdxfhBQ97x/xMYX2vdJ0AfVe0HbATu/YHnn62qA1Q1s3ERjXFWsy8XkB+VQpdR3Z2O0qQOZI4l7Wg2BzbtdzqK8bGTFr+qLgYO1Fq3wHNNXYClQGc/ZDPGceVHykkvXMT23qF/NE9tCVeMdV+V668LnY5ifMwXc/w/A+o7j6sCn4rIMhGZ+kM7EZGpIpItItlFRUU+iGWM99bNXEILjhJ7QfhM85zQ+6pMiomncq5N94Qar4pfRO4DKoFZ9WwyUlUHAOcD00RkVH37UtUXVDVTVTMTExO9iWWMzxx6awGVRJB269lOR2lykbGR5HY8hx55C+z0DSGm0cUvItcBFwBXqmqdPxWqutPzZyHwHhD6B0GbkNJ+5QLWthpOq86tnI7iiPLRY+lctZ38+RudjmJ8qFHFLyLjgbuBi1S1tJ5t4kSk5Yn7wFggp65tjQlEBzbtp3fpMg6eEX7TPCck3+Qe+46XbLonlDTkcM7XgSVAqogUiMgNwLNAS+ATz6GaMzzbdhSROZ6ndgC+FJFVwDfAx6o6zy+jMMYP1v91IS6UhCvC74PdE7qO7sG2yJ52+oYQc9KvIarqlDpWz6xn213ABM/9LUB/r9IZ46CqOQs4JK3pfVV4H4m8LXUsA9e+SkVpBVHNo5yOY3zAvrlrTB20Wum++ZOwOk1DfWIu9Jy+4aWlTkcxPmLFb0wd8udvpHPVdirOCt/5/RN633K2+/QNb9p0T6iw4jemDidO03Diw81w5j59w1ASV9h5e0KFFb8xdWj2xQK2RvWi6+geTkcJCAfOcJ++4eDmAyff2AQ8K35janGfpuEztqWG79E8tbWb4jl9w3N2+oZQYMVvTC3fnqbhQiv+E9KuGUwx8VTMsXn+UGDFb0wth95aQAWRYXmahvrY6RtCixW/MbW0XzmfdWF8mob6lJ89zn36hnkbnI5ivGTFb0wN+3KL6F26nIODbZqnthMXojlxxJMJXlb8xtSw8W/u0zQkXmnFX1uXUd3Jj0qh2ZdW/MHOit+YGqrmLeCgtKH3lWc4HSUgbe89lozCzzhectzpKMYLVvzGeGi10mvzAtZ3OpeI6Ain4wSk2AvHEkcpuS8tcTqK8YIVvzEemz9cR1L1TirH2DRPfdJuGU0FkRx6077FG8ys+I3xKPiHe+76xIeY5vtadW7FulbD6bDKij+YWfEb49Hiy3lsju5N5zO7OR0loB0cMo60YysoytnrdBTTSA25EMtLIlIoIjk11rUVkU9EZJPnzzb1PHe8iGwQkTwRme7L4Mb40rEDx8jYv5iCjPFORwl4Ha51/zfa+JxdlStYNeQd/z+B2n8bpgMLVTUFWOhZ/g4RiQCew32h9XRgioike5XWGD9Z+9fPaUYZcZOs+E8mdfJAiiQR5tkF9YLVSYtfVRcDtU/JdzHwsuf+y8CP63jqECBPVbeoajnwhud5xgSco7PncYxYMm4Z5XSUgOeKdLEheRyp2+ZTXVntdBzTCI2d4++gqrs99/fgvr5ubZ2AHTWWCzzrjAk4XdbNIydhNM3aNnM6SlCQ8eNJ0H2s/9dyp6OYRvD6w11VVcDrszaJyFQRyRaR7KKiIm93Z0yDFXy5lR7lGygdOc7pKEHj9GnuI58KX7HpnmDU2OLfKyJJAJ4/C+vYZifQpcZyZ8+6OqnqC6qaqaqZiYmJjYxlzKnb8jf3oYldbrL5/YZKzGjPuuZn0CbLDusMRo0t/g+Aaz33rwXer2ObLCBFRLqLSDQw2fM8YwJK9GfzKIjoRvfxqU5HCSqFA8eTUbKE4m2HnI5iTlFDDud8HVgCpIpIgYjcADwMnCcim4BzPcuISEcRmQOgqpXAbcB8IBd4S1XX+mcYxjRO+ZFy0ncvZMvp4xGXOB0nqLS9YjyRVLHuWbsqV7CJPNkGqjqlnofOqWPbXcCEGstzgDmNTmeMn62buYQBHCb6QpvfP1Vp1w2leForKj6cB49OcjqOOQX2zV0T1g69Mc99ta1pY5yOEnSimkeR2/FceuXNs6tyBRkrfhPWOq78mJzWI4nvGu90lKBUcd4EOlYVsOndNU5HMafAit+ErV1f7+D0sjUUj5hw8o1NnVLuOB+AXS/ajG4wseI3YSvvaXdZdf75RIeTBK/TBnUkt9lA2iz52Oko5hRY8ZuwFbPwY3ZEJtPzgjSnowS1vWdMpE/JVxzcXPvMLiZQWfGbsFR2qIw+exeypfcEO4zTS+2unkAE1eQ+ZdfiDRZW/CYs5Tz3OXGU0vxSm+bxVvp1Q9gnCVR/ZNM9wcKK34Slo2/P4Rix9Ln9bKejBL2I6AjWJ4+n99Z5VJVXOR3HNIAVvwk7Wq0kr/2YnMQxdjZOH3FNnECC7iP3lSyno5gGsOI3YSd//ka6VW6mdIxN8/hK2l3jqMLFvpdtuicYWPGbsLP9efdhnD1vs+P3faVNz7bktBpBh2VW/MHAit+EnfjPPyQvJp3OI5OdjhJSDg6fSNqxFezOrvfs6yZAWPGbsHJw8wH6HlpMwSC7Cqivdb71IgA2PWZnXw90VvwmrKx7bC6RVJFwgxW/r/W8II2tUb1o9qkVf6Cz4jdhxfXh++x1nUb6tYOdjhJyxCVs7XcR/fb9h8O7Djsdx/wAK34TNo6XHKdPwVw2nH4hrkj70feHNtdeTAzl5PzFrsUbyBr90y8iqSKyssatRETuqrXNaBEprrHNH7yPbEzjrHlmES05QrPLbZrHXzJuGsF+aUfVv226J5Cd9Apc9VHVDcAAABGJwH0h9ffq2PQLVb2gsa9jjK+Uvv4+R2lOnzvsoiv+EhkbSW73iWTkf0hFaQVRzaOcjmTq4Kvfd88BNqvqNh/tzxif0molJfcD1iSNs2/r+lnkpItpowfJmfGl01FMPXxV/JOB1+t5bISIrBaRuSKS4aPXM+aUrP/XcpKqd1I50aZ5/K3PL8dSRgwlr9l0T6DyuvhFJBq4CHi7joeXA11VtR/wDPDvH9jPVBHJFpHsoqIib2MZ8x17X3ifKlyk/dpO0+BvLU5rwZrEc+i+5n27Fm+A8sU7/vOB5aq6t/YDqlqiqkc89+cAUSKSUNdOVPUFVc1U1czExEQfxDLm/3T65j3WxI+kXWqdP37Gx46N/zFdK/PZ+M5qp6OYOvii+KdQzzSPiJwmIuK5P8Tzevt98JrGNNiWOetJOZ5D8bmXOh0lbKTf+2OqcLH7mXecjmLq4FXxi0gccB7wbo11N4vIzZ7FS4EcEVkFPA1MVlX73c80qe1PzAbg9OmXOJwkfCSkJbKqzWi6fvO2TfcEIK+KX1WPqmo7VS2usW6Gqs7w3H9WVTNUtb+qDlPVr7wNbMypSvrybVa3HEFSZieno4SVI+MupUf5BvLeX+t0FFOLfX3RhLStn2witWwVB8bYNE9T633vT6hG2PmUTfcEGit+E9K2Puae5km5Z5LDScJP+36nsTr+R3RaasUfaKz4TUhr/8U75MQNodPwrk5HCUvF515KyvG1bP4o1+kopgYrfhOydizOJ710GfvOsmkep5z4QH3Hk7MdTmJqsuI3IWvzI+6y6fEbm+ZxSlJmJ/cH6/+16Z5AYsVvQlbiorfIbTaIrqN7OB0lrB0YcympZavYMneD01GMhxW/CUlbP9lExtEs9o6Z4nSUsJf6u59SjbD94X85HcV4WPGbkLT1wVlUI6Teb8XvtKTMTqxsczbdv5plX+YKEFb8JuRotZL81SxWtR5tX9oKEEcvvpJulZtZ+49vnI5isOI3IWjdy1kkV+Rx5OIrnY5iPPo9MIkyYtj/9Cynoxis+E0I2vfULI4TTb8H7GieQBHfNZ4VnS4gfc2bVJZVOh0n7Fnxm5BSWVZJ+uo3WNHxAuK7tXY6jqlBrrySRC1k5V8+dTpK2LPiNyFl5WMLSdRCuNKmeQLNwPsmcEhaU/aSTfc4zYrfhJSyl2ZxSFoz4LcTnI5iaolpFcPq1MsYkP8eRwuPOh0nrFnxm5BRUlDCwC2zWZ16GbGtY52OY+oQf+tVtOAoK35np3BwkhW/CRkrp79BHKW0/fUNTkcx9eg37UdsjepFyzdfdDpKWPP2ClxbRWSNiKwUkew6HhcReVpE8kRktYgM8ub1jPkh7f79Ipti+pBx/RCno5h6iEvYes6N9C/5wk7h4CBfvOM/W1UHqGpmHY+dD6R4blOBv/ng9Yz5no3vrCbjaBa7JtyIuMTpOOYHpP/5WiqJYPv9LzkdJWz5e6rnYuAVdVsKtBaRJD+/pglDu/80k+NE0/fPVzkdxZxE+36nkZ10IRlZ/6SitMLpOGHJ2+JX4FMRWSYiU+t4vBOwo8ZygWfd94jIVBHJFpHsoqIiL2OZcFJ2qIx+q19lWddLaJvSzuk4pgEipt5Iohay7IGPnI4Slrwt/pGqOgD3lM40ERnV2B2p6guqmqmqmYmJiV7GMuFk+e/fo40eJPa2G52OYhpo4PRx7HZ1IuIf9iGvE7wqflXd6fmzEHgPqP2p2k6gS43lzp51xvhM7KwX2R7ZnQG/ONvpKKaBImMj2TDiegYVzWPX1ztO/gTjU40ufhGJE5GWJ+4DY4GcWpt9AFzjObpnGFCsqrsbndaYWvI+WMegg/9hy9k34oq0o5ODSa+H3IfdbvzlDIeThB9v/qZ0AL4UkVXAN8DHqjpPRG4WkZs928wBtgB5wN+BW71Ka0wtu6Y/TRkx9H2mro+YTCDrPDKZrKSL6LvkeY4dOOZ0nLDS6OJX1S2q2t9zy1DVBz3rZ6jqDM99VdVpqtpTVfuq6veO9TemsQ7lHyQz9xWyUq6kXWqC03FMI8T8+g7a6X6yf/W601HCiv1ubILWyttepDnH6PC/dzodxTTSgLtGsyG2H+3feMquztWErPhNUKosq6TX/GdZ0Xo0p1/az+k4ppHEJRRefgepZatZ9fTnTscJG1b8Jihl/f59Oldtp/xme7cf7DIfv4L90o6yR55yOkrYsOI3Qan5C0+xPbI7mQ9c6HQU46VmbZuxZvjPGbL7fXYsznc6Tliw4jdBZ83zX7lP8jXxdiKiI5yOY3wg9albqSSSLdP+4nSUsGDFb4JO2e//xD5JYPALdghnqEjK7MTS3tcxNGcme5bvcjpOyLPiN0El97VlDC6aS855vySufZzTcYwPJf9tOpFUsv4me9fvb1b8JqgU3/0nDklrBs2c5nQU42NdR/dgaY8rGbJ8BkVrC52OE9Ks+E3Q2Dh7DcN2/5uVo+6kVedWTscxftDx2d8SSxlrb3zC6SghzYrfBI19v3iQw7Sg/8w7nI5i/KTH+aks7fJTzlj6LAc3H3A6Tsiy4jdBYePsNQzb8RbLhk6jTc+2TscxfpTwxH205Airrra5fn+x4jdBofjmuymRePrPutvpKMbPTp/Ul/8mX8nQJU/YKZv9xIrfBLxlD3/C4H3zWDnxd/ZuP0x0e+1BBGXLFb9zOkpIsuI3Aa2qvIoWD/yGHZHJDJ91m9NxTBPpfGY3lg6+gxFbXmXDmyudjhNyrPhNQFsy7TVSy1ax45aHiGkV43Qc04QGvv1bDkkbjt7yaztzp49Z8ZuAdbTwKD1fuo+1cYMZ/uTlTscxTSy+W2vW/PgPDDq4kOz/N9fpOCHFm0svdhGRz0RknYisFZHvnSZRREaLSLGIrPTc/uBdXBNOsibeT1L1TqoeeRxxidNxjAOGv3IL+VGn0+FPt1G6r9TpOCHDm3f8lcCvVDUdGAZME5H0Orb7QlUHeG5/9OL1TBjJnbWcH2U/zuLeN9Hv1pFOxzEOiW4RTfGjL9C1Mp9vJtzvdJyQ4c2lF3er6nLP/cNALtDJV8FM+Kosq4SbbmK/K5F+c/7sdBzjsAF3nsUXp9/AyKzHWf/6CqfjhASfzPGLSDIwEPi6jodHiMhqEZkrIhk/sI+pIpItItlFRUW+iGWC1JeXPUXaseVsvvMZWndv43QcEwD6znuUA65U05oRAAAMSUlEQVQEqm+8yf3GwHjF6+IXkRbAbOAuVS2p9fByoKuq9gOeAf5d335U9QVVzVTVzMTERG9jmSC1bWEegz/6A193uJBhf7nU6TgmQLTu3oa8O54mvXQZX06y8/h4y6viF5Eo3KU/S1Xfrf24qpao6hHP/TlAlIgkePOaJnQdLznO0YsmUy4xdHn/OftA13zH8McuY2nSTxgx5z7WvZLtdJyg5s1RPQLMBHJV9fF6tjnNsx0iMsTzevsb+5omtC0dPZ300mVsvOclOg7t4nQcE2DEJaR+8SKFEUnE3XA5xduLnY4UtLx5x38mcDUwpsbhmhNE5GYRudmzzaVAjoisAp4GJquqfRPDfM83v/+Qs1Y8yef9bmfoQz92Oo4JUG16tuXAM/+iU+U21o78uX2xq5EkEHs4MzNTs7PtV7lwUfDfbTT/0SAKY7uSvGsJsa1jnY5kAtyisf/L6E/uY/EVMxg16+dOxwkIIrJMVTMbsq19c9c4qnh7McfOmUiEVhLz3ptW+qZBRs2ZTna7cQz/120sf+RTp+MEHSt+45iK0gryBl5K8vENbH7kXbqPO93pSCZIuCJdpCx/k60xvel5zyQ2vZfjdKSgYsVvHKHVytIBN3PGgU/5+oa/M+g35zgdyQSZ+K7xNP/sY4654mh22UT2rtztdKSgYcVvmpxWK58PvZsfbXqJRaP+wMgXr3M6kglSnYZ35eArH9G6aj/Fw8fZRdobyIrfNCmtVhYPuovR2X/h877TOOuz+52OZIJc2pWD2PjI+3Quy6Nk0Gh7598AVvymyVRXVvNFn1s4a9XTLDrjl4xa+Yx9Scv4xKDfnMPGJ+fSoXw7pUPOYndWgdORApoVv2kSR/Yc4evknzIq93kWDb+Xs775i5W+8akBd55F/owFtK3YS/XwEeTOWu50pIBlxW/8bvuiLezqPoIhO99j0cWPc9aXD1rpG7/o+/MR7HnjcxQh+aoz+er2152OFJCs+I1fZT0whxZjBpN4vICVD81j9L9/YaVv/Cr18gHErMpiQ/wQRjx7BYsyf8XxkuNOxwooVvzGL4q3HeKLlOsZfP9E9sV0ouTTLM6Yfp7TsUyYSMxoT8auT/m87zRGL3ucHe0Hse7lLKdjBQwrfuNTWq0svftdSnv0YXjeqywa8Vu67c2i25ieTkczYSaqeRRnrX6WrAfmEFdRzOnXDWfR0HsoKah99vjwY8VvfGbN81+xpvVIhj06iSNRbdj48lJG//dBYlrFOB3NhLHBfzif5ltyWJJyLaO/eYTjXXvx+WXPUn6k3OlojrHiN17RamX5I5+S1X4ifW8+k/ZH81l89Qt0P7SC9GsadL4oY/wuvltrfrRxJutezqIgPoOz3rmd3W3S+PzSZziy54jT8ZqcFb9plAOb9rP4ihlsiuvPoHvOI3lfNovO/RMtd29i1Cs3ERkb6XREY74n/ZpMBuz/D1n/8xElse05a/YdVCV1ZlHmr9j4zmqn4zUZOy2zabCitYVseGY+0e+9xcDCeURRyYbYfhROuYvBj0+xM2uaoJPz4lJK/vgkQ3a8QyRVbIrpw85RU+g49QJSLukbVEegncppmb0qfhEZDzwFRAAvqurDtR4Xz+MTgFLgOlU96bcqrPgDQ+HqPeS/vpRj/1lChzWfknbM/b9ut6sTGwZNIenXV3L6Zf2D6i+HMXXZl1vE2gfeps2cWfQ7/BUAe1xJ5CWfR/WZI0n6yXB6TEwjIjrC4aT1a5LiF5EIYCNwHlAAZAFTVHVdjW0mALfjLv6hwFOqOvRk+7bibzpV5VXsXbGLom/yKVmxmaqcXOK25tDpwBo6Vrm/9l5OFLmthnJwyHg6XDue1MkDcUXaLKEJTbuzd5L31wVEfDqf1IJPaafuq8UepgVbW/TlYOc+aFoGzfv1om1mDzqOSKZZ22YOp2664h8O3K+q4zzL9wKo6kM1tnkeWKSqr3uWNwCjVfUHz6JkxX9y1ZXVVJRWUH6knIqj5ZQfPk55SRnlxccoP1RK+aFSKg4dpWJfMZUHSqg+WAz79xNxaD/RJfuIO7ybtmW7SKzaQxSV3+63jBi2xfZmX4cMKgZk0m7CMFJ+OtCmcUxY0mpl28I8CmYvpeqrr4nfnkOXkpxv/zE4Yb+0Y39MR0rikihr1Z6q+HZo23ZI2zZEtGlFVEI8UW1aEBXfnJi2cUTHNyO6Vaz71jKGqOZRRLeI9uo3ilMpfm8+gesE7KixXID7Xf3JtukE+OX0eeviMomuOnZKzxGvPuOo+7lSY/3/7V+/Xe9e514WtI7lagTFRTUurSKCqm//jKSSKCqIoJoY4FQOlDxONAdcCRyObkdJXBKbk9LZmJiEq0cycX26kzC4O53OTCY1NpLUxvznMCbEiEtIPi+F5PNScF9i3P2PQVFuEXuXbKF4xRYqNmzBtWcXMft3EXd4N+2LNxGfv594Tv37Antdp9Ghyv9nFw2YQy9EZCowFaBr166N2sfB9r1xVTTiq9nS+Dlqrfe5NdZ7tnFvK9+uUxHPY/LdZXGhLheIC1wu1BUBERFoRARERkFkpPsWHQ3R0Uh0FNK8Ga7msbiaxRAZH0dU6zii28TRrEMrmp/Wipad42me0Jwkl5DU6NEaY8QlJGa0JzGjPTCs3u3Kj5RzeGcJR3YWU7q7mPKDR92/hRcfpfrIMaqPHae6tAw9VgYVFe5bbCwdmmAM3hT/TqBLjeXOnnWnug0AqvoC8AK4p3oaE+jM/Nca8zRjjPG56BbRtEtNoF1qgtNRvsebT+iygBQR6S4i0cBk4INa23wAXCNuw4Dik83vG2OM8a9Gv+NX1UoRuQ2Yj/twzpdUda2I3Ox5fAYwB/cRPXm4D+e83vvIxhhjvOHVHL+qzsFd7jXXzahxX4Fp3ryGMcYY37KDsY0xJsxY8RtjTJix4jfGmDBjxW+MMWHGit8YY8JMQJ6WWUSKgG2n8JQEYJ+f4gSqcBwzhOe4w3HMEJ7j9mbM3VQ1sSEbBmTxnyoRyW7oyYlCRTiOGcJz3OE4ZgjPcTfVmG2qxxhjwowVvzHGhJlQKf4XnA7ggHAcM4TnuMNxzBCe426SMYfEHL8xxpiGC5V3/MYYYxooqItfRMaLyAYRyROR6U7n8RcR6SIin4nIOhFZKyJ3eta3FZFPRGST5882Tmf1NRGJEJEVIvKRZzkcxtxaRN4RkfUikisiw0N93CLyC8/Pdo6IvC4isaE4ZhF5SUQKRSSnxrp6xyki93r6bYOIjPNVjqAtfs/F3p8DzgfSgSkiku5sKr+pBH6lqum4L/kzzTPW6cBCVU0BFnqWQ82dQG6N5XAY81PAPFXtDfTHPf6QHbeIdALuADJVtQ/u07xPJjTH/E9gfK11dY7T83d8MpDhec5fPb3ntaAtfmAIkKeqW1S1HHgDuNjhTH6hqrtVdbnn/mHcRdAJ93hf9mz2MvBjZxL6h4h0BiYCL9ZYHepjjgdGATMBVLVcVQ8R4uPGfYr4ZiISCTQHdhGCY1bVxcCBWqvrG+fFwBuqelxV83Ff12SIL3IEc/HXdyH3kCYiycBA4GugQ40rmu2BJrlcZ1N6ErgbqK6xLtTH3B0oAv7hmeJ6UUTiCOFxq+pO4C/AdmA37iv1LSCEx1xLfeP0W8cFc/GHHRFpAcwG7lLVkpqPeS56EzKHaInIBUChqi6rb5tQG7NHJDAI+JuqDgSOUmuKI9TG7ZnTvhj3P3odgTgRuarmNqE25vo01TiDufgbfCH3UCAiUbhLf5aqvutZvVdEkjyPJwGFTuXzgzOBi0RkK+5pvDEi8hqhPWZwv6srUNWvPcvv4P6HIJTHfS6Qr6pFqloBvAuMILTHXFN94/RbxwVz8TfkYu8hQUQE95xvrqo+XuOhD4BrPfevBd5v6mz+oqr3qmpnVU3G/f/2P6p6FSE8ZgBV3QPsEJFUz6pzgHWE9ri3A8NEpLnnZ/0c3J9jhfKYa6pvnB8Ak0UkRkS6AynANz55RVUN2hvuC7lvBDYD9zmdx4/jHIn717/VwErPbQLQDvdRAJuAT4G2Tmf10/hHAx957of8mIEBQLbn//e/gTahPm7gAWA9kAO8CsSE4piB13F/jlGB+7e7G35onMB9nn7bAJzvqxz2zV1jjAkzwTzVY4wxphGs+I0xJsxY8RtjTJix4jfGmDBjxW+MMWHGit8YY8KMFb8xxoQZK35jjAkz/x+8gnCHEzihyAAAAABJRU5ErkJggg==\n", "text/plain": [ "