{ "cells": [ { "cell_type": "markdown", "id": "9a521ceb-3578-4995-9369-bdd4a2e483f6", "metadata": {}, "source": [ "# An Overiew of numpy.random\n", "\n", "***\n", "\n", "[Official numpy.random documentation](https://numpy.org/doc/stable/reference/random/index.html)\n", "\n", "***" ] }, { "cell_type": "code", "execution_count": 1, "id": "9a586b12-cbd7-4177-a10a-734b3c4487ba", "metadata": {}, "outputs": [], "source": [ "# Numerical arrays.\n", "import numpy as np\n", "\n", "# Plotting.\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "id": "bd16f9dc-97df-4fe8-8617-ff89147bfc8b", "metadata": {}, "outputs": [], "source": [ "# Change plot style.\n", "plt.style.use('fivethirtyeight')\n", "\n", "# Change default figure size.\n", "plt.rcParams['figure.figsize'] = [12, 8]" ] }, { "cell_type": "markdown", "id": "1ae54d4e-f19a-4b10-874e-f6f2605fc5b3", "metadata": {}, "source": [ "
\n", "\n", "## From the official documentation for np.random.normal\n", "\n", "https://numpy.org/doc/stable/reference/random/generated/numpy.random.Generator.normal.html#numpy.random.Generator.normal\n", "\n", "***" ] }, { "cell_type": "code", "execution_count": 3, "id": "a63a3795-d7ea-43e3-a05e-493ddf4ab69b", "metadata": {}, "outputs": [], "source": [ "mu, sigma = 0, 0.1 # mean and standard deviation\n", "s = np.random.default_rng().normal(mu, sigma, 1000)" ] }, { "cell_type": "code", "execution_count": 4, "id": "f4485a0b-f6c0-4e4c-b380-d055fa4bae80", "metadata": { "tags": [] }, "outputs": [ { "data": { "text/plain": [ "array([ 0.05096108, 0.06106065, -0.18797417, 0.07885399, -0.02960684,\n", " -0.0775354 , 0.16937639, 0.140479 , 0.09634012, -0.00776155])" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "s[:10]" ] }, { "cell_type": "code", "execution_count": 5, "id": "ca984002-db19-4f8d-904c-dc331cb16931", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.0026846616194203463" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "abs(mu - np.mean(s))" ] }, { "cell_type": "code", "execution_count": 6, "id": "5b84abe1-a3bd-4bb6-a597-60f93fa890de", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.002022387933364428" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "abs(sigma - np.std(s, ddof=1))" ] }, { "cell_type": "code", "execution_count": 7, "id": "50e9e6bf-5a39-4202-877f-e3244b21995c", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxcAAAHyCAYAAACZNjMkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABrSElEQVR4nO3dd3hUVeLG8e+dlkILHQKhFyGKIIICinREEUTB3lhdC4j+1o67FuyuZVUQLNiVXQWpKkWKdARBBEEg9JoEAqGEJFN/fwSiKJDCJGfK+3mePEtm7kze2WuSeXPuOcfKzMwMICIiIiIicoZspgOIiIiIiEhkULkQEREREZGgULkQEREREZGgULkQEREREZGgULkQEREREZGgULkQEREREZGgULkQEREREZGgULkQEREREZGgULkIspSUFNMRpBh03sKXzl140nkLTzpv4UvnLjyF43lTuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBQuRARERERkaBwmA4gIiJmJXy0q8SeO3NgrRJ7bhERCT0auRARERERkaBQuRARERERkaBQuRARERERkaAocrl47bXXSEhI4OGHHz7tcWvWrOGyyy6jRo0aNGvWjJdffplAIFDsoCIiIiIiEtqKNKF72bJlfPLJJyQnJ5/2uEOHDtGvXz/at2/P7NmzSUlJYfDgwcTHxzNkyJAzCiwiIgIlNxFdk9BFRIqv0CMXBw8e5O9//zvDhw8nISHhtMeOHTuW7OxsRo0aRfPmzenbty/3338/I0eO1OiFiIiIiEiEKnS5+L//+z/69u3LJZdcUuCxS5cupV27dsTFxeXf1rVrV/bs2cO2bduKl1REREREREJaoS6L+uSTT9i8eTPvvvtuoZ40PT2dxMTEE26rWrVq/n316tU76eNSUlIK9fyhLlJeR7TReQtfOndnKr5EnrVp1m52LNpITuXqYLP/5f4zP28lk1v/PZ2e/v8JXzp34SkUz1vjxo1PeV+B5SIlJYVnnnmGqVOn4nK5Cv1FLcs64fPjl0P9+fY/Ol3QcJGSkhIRryPa6LyFL527IFgQ3LkLNXIP8NrGz7l27xJYBgGnC3+teviTGuKv0xB/UkNSPBYNWrQ8sy8U5NzH6b+nU9P3W/jSuQtP4XjeCiwXS5cuJSMjg3bt2uXf5vP5WLRoER9++CG7d+8mJibmhMdUq1aN9PT0E27bt28f8PsIhoiIRBZbwM/du2by7JavqODLJssWQ1ylitj2pWLfugH71g35x7YA/JWqHiscjfAnNcSX1IBAjdonHeUQEZHwUGC5uPzyy2nVqtUJtw0ePJiGDRvywAMPnHQ0o23btjz99NPk5OQQGxsLwJw5c6hZsyZ169YNUnQREQkVrQ5vYeSGD2lzeDMAUyqfx/2Nb2HVPS3h6BFsOzdj374J2/ZN2HZswtq+Cdv+vdj274VfluQ/T8AVg79W/fwRDl9SQ/xJDaBMuVJ7LSW1ChVoJSoRiXwFlouEhIS/rA4VHx9PxYoVad68OQDDhg1j+fLlTJ48GYD+/fvz8ssvM2jQIB566CE2btzIG2+8wSOPPHLay6JERCS8lPMeZdiWcQzeNQM7AXbEVOL+RrcyuUprOP7zPr4s/iYt8Ddpkf+4lPXraVIhHtv2Tdi3b8S2YxO2HZuxZaRh37IO+5Z1J3wdf5Xq+GvnXVblq9MQf53wukxARCRaFGmfi1NJTU1ly5Yt+Z9XqFCBCRMm8NBDD9G5c2cSEhIYPHgw9957bzC+nIiImBYIcPXepfxn46ckujPxYuP12r0YVu9qshyxBT/eZiNQIwlfjSR8bTv9fnvWYWw7NuWNcuw49rFzC7Z9adj2pcHKRfmHvl2zC0OaDMRvFXk/WBERKSHFKhfffvvtCZ+PGjXqL8ckJyczderU4qUSEZGQVT87neEpH3Hp/lUALCnfiEFN/saqskG47LVMOfxntcR/Vsvfb/N5sdJ2HRvh2JxXPtau4K49synry+VvZ92FT/M0RERCQlBGLkREJPI5/V4e3PEt/9w2gTi/hwOOeB5vcB2ja3YmUJKjB3YHgcS6eBPrwoVd82767Wd8/36MG9MXEut3c1Pze/HY9CtNRMQ0jSWLiEiBOmb+xvKfhvLclq+I83v4vHoHktu+yvuJXUu2WJyCr1krep47lEx7PFfvW8a4X/9DjM9d6jlEROREKhciInJKVdyH+OC3d5i98jmaH93N+riadD/3cW5rNoh0VwWj2X6s0JjuLR9nn6Msl+9fyeTVrxLvyzGaSUQk2qlciIjIX1gBP3/bPYe1Sx/i1rT55FhOnqrXn1ZtXmROxWTT8fL9XK4+XVv+i1RnBbpmrmHqLy9TznvUdCwRkailciEiIidIPrKDH35+lvc2jKaSN4vvK57NuW1e4vl6/XDbnKbj/cWaskl0bvUEO2Iq0eHQBmb88iIVPUdMxxIRiUoqFyIiAkC8L4cXN/2Xn5b/kw6HNrDHlcANze6lV4vH2BRfw3S800qJr0mnlk+yObYqbQ5vZubK56nqPmg6lohI1FG5EBEReu9bweqlj/Dwjm+wB/y8ndids9v8m6+qt/t9M7wQty2uKp1aPcm6uJqcm7Wd2SufIzF3v+lYIiJRReVCRCSKWRnpjP31P0z89TXq5mawomw92p83jPub3MZBZxnT8Ypsd0wlurR6glVlkmh2dDdzfn6WOjl7TccSEYkaKhciIlHKOrCPuKfvot++nzhkj+UfjW6m3XnP8FP5hqajnZF0VwW6tfwny8o1oGFOOj/8/CyNjqaajiUiEhVULkREopHPS+w7z2I7dIAFFZpydttXGF770ojZ6Xq/sxw9zx3KovKNqZObwZyVz9Isa6fpWCIiEU/lQkQkCrkmfYp93S/4K1TimuT72R1TyXSkoDvkiKdXi8eYnZBMTXcms1c+R8vDW03HEhGJaA7TAUREpHTZ1/yEc/JnBCyL3HueIH1JyW2Gl/DRrlPcEw8LTnVf8GQ5YulzzkOMXfMGvfb/wve/PM/lLR5laflGJf61RUSikUYuRESiiJWZQcy7z2MFAnj63oKvWSvTkUpcjt3F1Wf/gwlVzqei9yjTf3mRizN/Mx1LRCQiqVyIiEQLv4+Yd5/HdvAA3matcPe9xXSiUuO2Obmu+X2Mqdaecr4cvl31b7rvX2U6lohIxFG5EBGJEs4pX+BYuwJ/+Yrk3v0viJDJ24Xls9m5rdk9fFCjE/F+NxNXv0bvfStMxxIRiSgqFyIiUcC2biWuCR/nzbO4658EEiqbjmSE37Jxd9PbGV6rBzEBL2PXvEH/9CWmY4mIRAyVCxGRSHcok9hRz2IF/Hh634jv7PNNJzIqYNn4R6NbeCWpN86Ajy/WjuCm1PmmY4mIRAStFiUiEgZOverS6VkBP9+seoWemRnMr9CUbgd74Cvmc0UUy2Jog+vIssfw9Nav+XjdO8T53byf2NV0MhGRsKaRCxGRCPbw9m/oeWAV+xxluanZ4IjZJC8oLIvn6l3Fow2uB2DUhg+5f8dUw6FERMKbyoWISITqkLmeZ7aMBeC2ZvewKzY651kU5LU6vRnS+Na8f2/6nKHbJpoNJCISxlQuREQiUCXPYT7/bQQO/LyS1JtplVuajhTSRtXqwR1N/44fi2e3jOWG1AWmI4mIhCWVCxGRSBMI8OG6d0nK3c/i8o15ov4A04nCwsc1OzGk8W0AvLrpCxI8WWYDiYiEIZULEZEI84+d39E742f2O8pwQ/N78dq0dkdhvZvYlXkVzqKa5xDPb/6f6TgiImFH5UJEJIJccDCFFzZ/CcDfzrqLHbFVDCcKM5bF4CYD8Vh27tozmwsPpphOJCISVlQuREQiRIIniy/WjsAZ8PFG7V58U6W16Uhh6bcytXkt6XIA3t7wIXa/z3AiEZHwoXIhIhIJAgFGr3+Pern7WFauAUMbXGc6UVh7vu6VbImtyrlZ2xmya7rpOCIiYUPlQkQkAty7azpX7vuJTHs81zcfgkfzLM5Itj2G+45N7n56yzhq52SYDSQiEiZULkREwlzrQ5v596YxANxx1p1sjatmOFFkmFq5JV9XaUNZfy7/2fip6TgiImFB5UJEJIyV9x5lzNrhuAI+RtTqwcSqbUxHiigPNLqZw/ZY+u37id77VpiOIyIS8lQuRETCVSDAe+vfp2FOOivK1uPRBtebThRxdsVW5ql6/QF4I+UT4n05hhOJiIQ2lQsRkTB11+5Z9N+7lEP2WK5vPoRcu8t0pIj0dq0e/Fy2LvVy9/GvrRNNxxERCWkqFyIiYejcw1t5bePnANzV9A42xdcwnChy+Wx2Bjf5G34s/rHzO5KP7DAdSUQkZKlciIiEmbLebP67djixAQ/v1uzC2GrtTEeKeEvLN+LdxK44Az7eTvkIK+A3HUlEJCSpXIiIhJNAgFEbPqBJdiq/lKnDg41uNp0oavyr/jWkOctz0cH13Jo6z3QcEZGQpHIhIhJG/rbnB65PX8wRWwzXJw8hR/MsSs1BZxkebHQTAC9v+i+V3YcNJxIRCT0qFyIiYeLsI9t5c+MnAAxq8jc2xCcaThR9/letPbMSkqnsPcJLm/9rOo6ISMhRuRARCQNlvDn8d+1w4vwePqxxCWNqXGQ6UnSyLO5tMpBcy8HA1LlcnPmb6UQiIiFF5UJEJAy8lfIxzY7uZk18Le5vfKvpOFEtJb4mL9fpA8CIDR/h9HsNJxIRCR0qFyIiIc6xYDq3ps3nqM3F9cn3kW2PMR0p6r1c5wpS4qqTfHQX/7dzquk4IiIhQ+VCRCSUZWcR89+3Abiv8a2sLVPbcCAByLW7GNJ4IABPbB1P3ey9hhOJiIQGlQsRkRDmnDYW68ghFpZvwsc1LjEdR/5gZqVz+F+1dsT73byV8jEEAqYjiYgYp3IhIhKqDmfimvYVAP9qcA1YluFA8mcPNbyRg/Y4Lt+/kr77fjIdR0TEOJULEZEQ5fpmDFbOUbzntGV+QjPTceQkUmMq5hU/4I2Nn1LWm204kYiIWQWWi/fff5/27duTlJREUlIS3bt3Z/r06ac8ftu2bSQkJPzlY+bMmUENLiISyaz96ThnTQDA3f8Ow2nkdN5N7Maycg1Iyt3PU1u/Nh1HRMSoAstFYmIiw4YNY+7cucyZM4eOHTty44038uuvv572cV9//TXr16/P/+jYsWPQQouIRDrX5M+wPB48bTrhr9fEdBw5Db9lY3CTv+HDYsjO6Zx7eKvpSCIixhRYLi6//HK6d+9OgwYNaNSoEU888QRly5Zl2bJlp31cpUqVqF69ev6Hy+UKWmgRkUhmpe3EMe87ApYN91UDTceRQlhRrj5v1+qBAz8jN3yIFfCbjiQiYkSR5lz4fD6+/vprsrKyaNu27WmPvfnmm2nUqBE9e/Zk0qRJZxRSRCSauCZ8jOXz4b2oJ4HEuqbjSCE9Vb8/u10JXHB4E3fsmWM6joiIEVZmZmaBa+etWbOGHj16kJOTQ5kyZXj//ffp2bPnSY/NyMhgzJgxXHjhhTgcDr777jtee+01Ro0axbXXXnvar5OSklK8VyEiEiFi03Zy1vvPELDbWXvPc3gSKgPQZkG84WRSGFen/8iXa9/igCOe5Lavku6qcML9yy46aiiZiEjwNG7c+JT3FapcuN1udu7cycGDB5k8eTKffPIJ33zzDc2bNy9UgAcffJDFixezaNGiwqcOUykpKaf9P1xCk85b+Iq0cxf7n8dxrFyEu/vVuG8akn97wke7DKaSQgsE+Gb1v7l0/yo+r96B25oNOuHuzIG1DAULjkj7fosmOnfhKRzPW6Eui3K5XDRo0IBWrVrx1FNPcc455zBy5MhCf5HWrVuzefPmYocUEYkGto1rcKxcRCAmFs8VN5qOI8VhWdzX+DaybU5uSltI5wNrTCcSESlVxdrnwu/343a7C3386tWrqV69enG+lIhIdAgEcI19HwBPj/4EKlQyHEiKa3NcdV6oeyUAIzZ8hMvvMRtIRKQUOQo64Omnn6ZHjx7UqlWLI0eOMG7cOBYsWMBXX+XtGjts2DCWL1/O5MmTARgzZgxOp5MWLVpgs9mYNm0ao0eP5umnny7RFyIiEs7sa5bjWLeSQHxZ3L1OPz9NQt9rSZdzQ9pCmh3dzcPbv+H5ev1MRxIRKRUFlou0tDTuvPNO0tPTKV++PMnJyYwbN46uXbsCkJqaypYtW054zKuvvsqOHTuw2+00bNiQESNGFDiZW0QkagUCuMbljVq4L78eypQzHEjOlNvm5N7GA5n1y/MM3TaJ/1Vrx6b4GqZjiYiUuALLxahRo4p0/w033MANN9xwZqlERKKIffkC7FvW469QCU/3q0zHkSCZW7E5n1a/mFvS5jM85WMua/Go6UgiIiWuWHMuREQkSPw+XF9/AICnz80QE2c4kATTow2vZ7+jDD0OrGbA3iWm44iIlDiVCxERgxyLZmLfvRV/lRp4OvU2HUeCbK+rAkMbXAfA6xs/h6NHDCcSESlZKhciIqZ4PbgmfASAu99AcDgNB5KS8GHNTiwu35ia7kxc4z80HUdEpESpXIiIGOL84Rts+1LxJ9bF276b6ThSQgKWjUFN/oYXG86ZE7B2bzMdSUSkxKhciIiYkJuNc/Knef+8+naw2Q0HkpK0umwdPqjZCSsQwPXtGNNxRERKjMqFiIgBzpkTsB08gK9+U3ytLzYdR0rBq3WuIGCz4Vj0PdbePabjiIiUCJULEZHSlnUY17f/BcDd/+9gWYYDSWnYElcN74Vdsfx+nFO/NB1HRKREqFyIiJQy19QvsbIO423WCl9ya9NxpBR5Ls/bB8o571uszAzDaUREgk/lQkSkFFkH9+OcMQ4Ad/87NGoRZfy16+NtfTGWx4Nz+jjTcUREgk7lQkSkFDmnfIGVm4O3ZXv8jZJNxxED3FfcCIBz9kQ4cshsGBGRIFO5EBEpJda+VJxzJhOwLNxX3246jhjir38W3rPbYOVk45w5wXQcEZGgUrkQESklrkmfYnk9eC/sir9OQ9NxxKDjoxeuGV9DzlHDaUREgkflQkSkFFi7t+GYP42A3Y67322m44hh/qbn4mt0NlbWIZxzppiOIyISNCoXIiKlwDX+I6yAH2/HywhUr206jphmWbj73ASAc9pX4M41HEhEJDhULkRESpht6wacy34g4HTi7nOL6TgSInwtLsBXpyG2zAwcC6aZjiMiEhQqFyIiJcw1bjQAnq79CFSqajiNhAzLwtM7b/TC9e1/wec1HEhE5MypXIiIlCDbul9wrF5KIDYed+8bTMeREONt0xF/jSRs+1JxLJltOo6IyBlTuRARKSmBADFf541auC+9BsolmM0jocdmx3181+5vxoDfbziQiMiZUbkQESkh9lU/Yt+wmkDZ8nguHWA6joQob/tu+CtVw757K/YVC03HERE5IyoXIiIlwe/Pn2vhvuImiCtjOJCELIcTz2XXAeD65nMIBAwHEhEpPpULEZES4Fg2F/v2jfgrVsHTpa/pOBLiPJdcjr98Rexb1mP/9SfTcUREik3lQkQk2HxeXOM/BMDd91ZwxRgOJCHPFYOnZ/+8f37zueEwIiLFp3IhIhJkjgXTsaXuwF+9Ft6Le5mOI2HC0/VKAvFlsK/7BduG1abjiIgUi8qFiEgwuXNxTfw475/9/gYOh9k8Ej7iyuDpdhUArm++MBxGRKR4VC5ERILIOWcytv178SU1xHtBZ9NxJMy4e1xNwBWL45cl2LalmI4jIlJkKhciIsGSfRTnlLy/OLv73w42/YiVIiqXgKfzFcCxfS9ERMKMfvOJiASJc8Y4bIcz8TVKxnduO9NxJEx5Lr2GgMOJY9kPWHu2m44jIlIkKhciIsFw5CCuqV8CkDvg72BZhgNJuApUqor3okuxAgFc3/7XdBwRkSJRuRARCQLXtLFY2Vl4z26D/6yWpuNImHNfdh0By4Zj0QysjDTTcURECk3lQkTkTOVm45w9GQD3lbcaDiORIFC9Ft4Lu2D5fDi/+5/pOCIihaZyISJyhhwLZmBlHcLXsBn+xmebjiMRwtP7BgCcc7/FOrjfcBoRkcLRAuwiImfC78c1YxwAN8Z0Y9xHuwwHkkjhr90A73kdcKxYiHP6ONzX3Gk6kohIgTRyISJyBuyrluTtxl25OhOqtDEdRyKMu/dNADhnTYSsw2bDiIgUgsqFiMgZcE7PG7XwdL8Kn81uOI1EGn/DZniTW2PlHMU5c4LpOCIiBVK5EBEpJtu2FBxrVxCIjcNzyeWm40iE8lyRN3rhmjEOco4aTiMicnoqFyIixeQ8NtfC0/EyiC9rOI1EKt9ZLfE1bI515BDOH741HUdE5LRULkREisHKzMCxeBYBy8LT/WrTcSSSWRbuPsfmXkz9Ejxuw4FERE5N5UJEpBicsyZi+bz4Wl9MoFqi6TgS4XzntsOX1BBb5j4cC6abjiMickoqFyIiReXOxTl7Ut4/e/Y3HEaigmXhueJGAFzf/hd8XsOBREROTuVCRKSIHAtnYB05hK9+U/yNzzEdR6KEt80l+KvXwrZ3N46lP5iOIyJyUioXIiJF8YdN8zw9rwHLMhxIoobNjvvyY7t2T/kc/H7DgURE/krlQkSkCOyrl2HbvQ1/pap421xiOo5EGW+HHvgrVcW+ayv2lYtMxxER+QuVCxGRInBOHwuAp9tV4HAYTiNRx+HE0+s6AFxTvoBAwHAgEZETqVyIiBSSbcdmHGt+IuCKxdOpt+k4EqU8l1yOv1wC9s2/YV+73HQcEZETFFgu3n//fdq3b09SUhJJSUl0796d6dNPvwzemjVruOyyy6hRowbNmjXj5ZdfJqC/rohImPt907xeUKac4TQStWJi8fQcAIBzyheGw4iInKjAcpGYmMiwYcOYO3cuc+bMoWPHjtx44438+uuvJz3+0KFD9OvXj2rVqjF79mxeeuklhg8fzogRI4IeXkSktFgH9+NY9H3epnk9tGmemOXp2pdAXBkcv/2MbeMa03FERPIVWC4uv/xyunfvToMGDWjUqBFPPPEEZcuWZdmyZSc9fuzYsWRnZzNq1CiaN29O3759uf/++xk5cqRGL0QkbDlnT8LyevC1ak+gem3TcSTaxZfF060fAK4pnxsOIyLyuyLNRvT5fEycOJGsrCzatm170mOWLl1Ku3btiIuLy7+ta9euPP/882zbto169eqdUWARkVLnzsUx6/imeQMMh5FwlvDRrqA9VxV3BzbbviJ+5WI6Dl/EvCHtg/bcIiLFVahysWbNGnr06EFOTg5lypTh888/Jzk5+aTHpqenk5iYeMJtVatWzb/vdOUiJSWlkLFDW6S8jmij8xa+SvrcVfp5PmUPZ3K0Rh3W2+LhlF8vvkRziPzRPld53q/Zhft3TeOR7VNISalaKl9XPyvDl85deArF89a4ceNT3leoctG4cWPmz5/PwYMHmTx5Mvfccw/ffPMNzZs3P+nx1p82lTp+OdSfby9K0HCRkpISEa8j2ui8ha8SP3eBAHEfPQ+Arc9NNG7S5NTHLgjeX6VFCuP1pMu4Z/f3DEhfQk65wQRqlOwle/pZGb507sJTOJ63Qi1F63K5aNCgAa1ateKpp57inHPOYeTIkSc9tlq1aqSnp59w2759+4DfRzBERMKF/defsO/aij+hMt4LOpuOI3KCXbGV+bTGxdgJ4Pp2jOk4IiLF2+fC7/fjdrtPel/btm1ZvHgxOTk5+bfNmTOHmjVrUrdu3eKlFBEx5MRN85yG04j81StJV+DHwrHoe6xDB0zHEZEoV2C5ePrpp1m0aBHbtm1jzZo1DBs2jAULFjBgQN6kxmHDhtGnT5/84/v3709cXByDBg1i7dq1TJ48mTfeeINBgwYVeFmUiEgosXZtxbF6KQFXDJ7OV5iOI3JSm+JrMKXyeVheD445U0zHEZEoV+Cci7S0NO68807S09MpX748ycnJjBs3jq5duwKQmprKli1b8o+vUKECEyZM4KGHHqJz584kJCQwePBg7r333pJ7FSIiJcA1PW/TPO9Fl0LZ8obTiJza27V70DdjOc7Zk/FcfgM4irQYpIhI0BT402fUqFFFvj85OZmpU6cWP5WIiGmHMnEsmg6AW5vmSYibnZCML7Ee9t1bcSyfh/eCLqYjiUiUKtacCxGRSOecMxnL48Hbsh2BmnVMxxE5PcvC0z1vUz3njPGGw4hINFO5EBH5M48b56wJef/UpnkSJrztuxOIL4N946/Ytm4wHUdEopTKhYjInziWzMJ28AC+Og3xNWtlOo5I4cTG4+l4OQDO7zV6ISJmqFyIiPxRIIBz2rHlZ3sOAK1yJ2HE0/VKApaF48dZcCjTdBwRiUIqFyIif2BfuwL7zs34K1TSpFgJO4FqifjOvRDL48E59xvTcUQkCqlciIj8Qf6meV2vBKfLbBiRYvB0z1vdzDlrIni9ZsOISNRRuRAROcbavQ3HL0sIOF14uvQp+AEiIciX3Bp/zTrYDuzDvmKB6TgiEmVULkREjnHN+BoAb4eeUC7BbBiR4rIs3N2vAsClid0iUspULkREAI4cxLHw2KZ5PfsbDiNyZrwdehCIK4N9wyps21JMxxGRKKJyISICOOdMwXLn4m1xAYHEuqbjiJyZ2Hg8F18KgHPmBMNhRCSaqFyIiHg9+W/AtGmeRApPt355y9Iu/h4OZ5qOIyJRQuVCRKKe48c52DIz8NWujy+5tek4IkERqF4bX4sLji1L+63pOCISJVQuRCS6BQI4p30FaNM8iTyebv0AcM6aBD4tSysiJU/lQkSimn3dSuzbN+IvXxHvhV1NxxEJKt/ZbfDXSMK2Px37ioWm44hIFFC5EJGo5px2bNO8Ln3BFWM4jUiQ2Wz5oxeumVqWVkRKnsqFiEQtK3UH9l8WE3A68XbtazqOSInwXNSTQGw89nW/YNu+yXQcEYlwKhciErWcM77GCgTwtutOoHxF03FESkZcmT8sS6vRCxEpWSoXIhKdjhzCOX8aAB5tmicR7vilUY7FM+HIQcNpRCSSqVyISFRyzv0Gy52D9+w2+Gs3MB1HpEQFaiThPactljsX59zvTMcRkQjmMB1ARKTUeb04v8+7PESjFhIpEj7addr7e3EJU1jKnklf0yTtIvxW4f++mDmw1pnGE5EooZELEYk6jmU/YDuwD19iPXzntDUdR6RUTKvUgpS46tTL3UfvfStMxxGRCKVyISLR5YRN8/pr0zyJGgHLxshaPQAYsmu64TQiEqlULkQkqtg2rMa+dQOBchXwtu9uOo5IqfqkRkeO2GLonLmW5CM7TMcRkQikciEiUcV1fNRCm+ZJFDrkiOeTGh0BGLxrhuE0IhKJVC5EJGpYabuw/7yQgMOZVy5EotDxS6NuTFtIgifLcBoRiTQqFyISNZyzJ+VtmndhFwIJlU3HETFifZlEvq94NmX8uQxM/cF0HBGJMCoXIhIdcnNwzstb39/T7SrDYUTMGlGrJwCDdn2PLeA3nEZEIonKhYhEBcfimVhHj+Br2Bx//aam44gYNbVySzbFVqN+zl4uy/jZdBwRiSAqFyIS+QIBnLMmAODpeqXZLCIhwP/HZWl3allaEQkelQsRiXi2lF+xb9+Ev1wC3radTMcRCQkf1+hIli2GrplraJa103QcEYkQKhciEvGcM/NGLbydeoPTZTiNSGg46CzDZzUuAvLmXoiIBIPKhYhENCszA8dPcwlYNjydrzAdRySkvH3s0qhbUudTQcvSikgQqFyISERz/PANls+H77wOBCpXNx1HJKT8VqY2sxKSKePP5bbUeabjiEgEULkQkcjl9eKcMwXQRG6RUxle+/iytDO0LK2InDGVCxGJWPYVC7Bl7sNfsw6+5ueZjiMSkr6r3IotsVVpmJNOr4yVpuOISJhzmA4gInIyCR/tKuSR8bDg5MfO+vlLLgH+r0xnRn68O//2zIG1zjygSITwWzZGJXbn35vHMHjXDL6toiIuIsWnkQsRiUhnH9nOJQd/47A9Nn9FHBE5uQ9rdiLLFkOPA6tpmrW74AeIiJyCyoWIRKS7d88E4LPqF3HYEW84jUhoy3SW4YvqHYC8uRciIsWlciEiEae89yg3pS4AYFSt7obTiISHkce+V25Jm09571HDaUQkXKlciEjEuSV1HmX9ucxJaM5vZWqbjiMSFn4tW4c5Cc0p58vhVi1LKyLFpHIhIhHFCvi5Z1feJVHHNwgTkcIZcex7ZtCuGVhallZEikHlQkQiStcDa2iavYcdMZWYUlmr3ogUxTeVz2NrTBUaZ6dx6f5VpuOISBhSuRCRiHJ8Mup7iV3x2eyG04iEF5/Nzju1ugEweNd0w2lEJBwVWC5ef/11OnfuTFJSEg0bNuTaa69l7dq1p33Mtm3bSEhI+MvHzJkzgxZcROTP6mbv5fKMn3Fbdj6o2dl0HJGw9GHNTmTbnFy6fxWNj+4xHUdEwkyB5WLBggXcfvvtTJ8+ncmTJ+NwOLjyyis5cOBAgU/+9ddfs379+vyPjh07BiW0iMjJ3LV7FnYCjK16IemuCqbjiISl/c5yWpZWRIqtwB26x48ff8Ln7777LnXq1GHJkiX06tXrtI+tVKkS1atXP7OEIiKFEONz87c9c4CCl58t/O7fItHp7Vo9uWPPD9yaOo8n6w8wHUdEwkiR51wcOXIEv99PQkJCgcfefPPNNGrUiJ49ezJp0qTi5BMRKZRr9i6hivcIy8vWY0n5RqbjiIS11WXrMLdCM8r7crgldb7pOCISRopcLh577DHOOecc2rZte8pjypYty7PPPstHH33E2LFj6dixIwMHDuTLL788o7AiIqcyaNf3wLFRC8synEYk/I2ofXxZ2u/Br2VpRaRwrMzMzEBhD3788ccZP34806ZNo169ekX6Qg8++CCLFy9m0aJFpzwmJSWlSM8pIpGrzYL4wh97aBOLVzxJhqMsddsNJ8fuKsFkItHB7veR8uM/qJObwcbr7+dww7NNRxKRENG4ceNT3lfgnIvjhg4dyvjx45kyZUqRiwVA69at+eKLL057zOmChouUlJSIeB3RRuctBC0o/LyIe46NWnxc8xIVC5Eg8dnsvJPYjRe2fEm9tT+Sc2k//awMYzp34Skcz1uhLot69NFHGTduHJMnT6ZJkybF+kKrV6/W5G4RCboq7kNcm74YPxbvJHYzHUckonxQszM5lhP7qh+x0rQQgogUrMBy8dBDDzFmzBhGjx5NQkICaWlppKWlceTIkfxjhg0bRp8+ffI/HzNmDGPHjmX9+vWkpKQwfPhwRo8ezZ133lkyr0JEotbf9vxATMDL1ErnsiWumuk4IhElw1WO/1VvhxUI4JythVlEpGAFXhY1evRoAPr27XvC7Y8++ihDhw4FIDU1lS1btpxw/6uvvsqOHTuw2+00bNiQESNGcO211wYrt4gItoCfu3bnbc45slYPw2lEItPIWj24LXUeznnfYZ2r/apE5PQKLBeZmZkFPsmoUaNO+PyGG27ghhtuKHYoEZHCuDxjBXVzM0iJq86MSueYjiMSkVaUq4+vYXPsm9ZS6del0FwTu0Xk1Iq8FK2ISKgYfGwi9zuJ3QlY+nEmUlI83foBUOWn2RAo9CKTIhKF9NtYRMJS06zddDvwK0dtLj6poUs1REqSt80l+MtXJD5tJ7YNq03HEZEQpnIhImHp7t15oxZfVO9AprOM4TQiEc7pwtupd94/Z04wHEZEQpnKhYiEnTLeHG5JnQ8c25FbREqcp/MVBCwbjuXzsPbvNR1HREKUyoWIhJ0b0xZQwZfNggpNWVW2ruk4IlEhUKkamWe1wvL5cP7wjek4IhKiVC5EJLwEAgw6NpF7ZKJGLURK077zOwPg+GEyeD2G04hIKFK5EJGw0vHgOs4+upM9rgQmVG1jOo5IVDlSpwm+2vWxHTyAY9k803FEJASpXIhIWBm0awYAo2t2xmMrcKseEQkmy8pfllYTu0XkZFQuRCRsJObu58q9P+Gx7LyX2NV0HJGo5G3fnUB8Gewbf8W2dYPpOCISYlQuRCRs3Ll7Ng78TKxyPntiKpqOIxKdYuLwXHwZAM5ZE81mEZGQo3IhImHB6fdyx+7ZAIzU8rMiRnm69gXAsXgmHDloOI2IhBKVCxEJC1ftXUoNz0FWl6nN/ApnmY4jEtUC1WvjbXEBlseNc95U03FEJISoXIhIWMhffrZWD7Asw2lEJH9i96yJ4PeZDSMiIUPlQkRC3rmHt9Lh0AYO2uMYU62D6TgiAvjOaYu/aiK2fanYf/nRdBwRCREqFyIS8u7ZnTdq8XGNS8hyxBpOIyIA2Gx4ul0JgPP78WaziEjIULkQkZCW4MnihrRFALxTq5vhNCLyR56LexFwxeBY8xPW7m2m44hICFC5EJGQdlvqXOL9bmZUPIeU+Jqm44jIH5Uph7dd3uptztmTDIcRkVCgciEiIcsK+Ll710xAy8+KhKr8S6PmT4Pso2bDiIhxKhciErJ67l9Fo5w0tsZU4bvKrUzHEZGT8NdphK9JC6ycozgWzTAdR0QMU7kQkZB1fPnZd2t1w2/px5VIqDq+LK1r5gQIBAynERGT9NtaREJSg+w0Lt3/CzmWkw9rdDIdR0ROw9v6YvwJlbHt3ob9t59NxxERg1QuRCQk3b1rJjYCfFntQjJc5UzHEZHTcTjwdO4DgHPmBMNhRMQklQsRCT25OdyWOheAUZrILRIWvJ16E7A7sK9YiJWRZjqOiBiiciEiIcexZBaVvFksLdeQn8o3NB1HRAohkFAZb9tOWAE/ztmTTccREUNULkQktAQC+ZdVaPlZkfByfGK3c+434M41nEZETFC5EJGQYktZjX37RvY6yzG26gWm44hIEfgbNsdXtwnW4YM4lv5gOo6IGKByISIhxTXjawDer9mFXLvLcBoRKRLL+n30QhO7RaKSyoWIhAwrIw378vkE7HZN5BYJU94LuxAoUx77lnXYNv1mOo6IlDKVCxEJGc6ZE7H8frxtOrEnpqLpOCJSHK4YPJdcDoBz5njDYUSktKlciEhoyM3OmwQKeHpcbTiMiJwJT5c+BCwLx9IfsA7uNx1HREqRyoWIhATHwhlYWYfxNWyGv2Fz03FE5AwEqtbE17I9lteDY+63puOISClSuRAR8wIBXN/nXT7h6dHfcBgRCYb8id2zJ4HPaziNiJQWlQsRMc7+60/Ydm/DX7EK3vMvMR1HRILA1/w8/DWTsB3Yh33FQtNxRKSUqFyIiHHO7/OWn/V0vRIcDrNhRCQ4bDY8XbUsrUi0UbkQEaOs1B04fllCwOnC06m36TgiEkSei3oSiI3DsW4ltp2bTccRkVKgciEiRjmPzbXwtusG5RLMhhGR4Iorg6dDTyBvqWkRiXwqFyJiTtZhnPOnAprILRKpjk/sdiycAVmHDacRkZKmciEixjjnT8XKzcHbrBX+pAam44hICQgk1sXb/Dwsdw7OhdNNxxGREqZyISJm+H35l0Rp1EIksuUvSztzIvj9ZsOISIlSuRARI+w/L8K2LxV/1UR8LS80HUdESpCvZTv8latjS9uJfc1PpuOISAlSuRARI5wzji0/270f2OyG04hIibI78HTpA2hZWpFIp3IhIqXOtn0TjnUrCcTG4bm4l+k4IlIKPB0vJ+B0Yv9lCVb6btNxRKSEqFyISKnL3zTv4l4QX9ZwGhEpFeUT8F7QBSsQwDlrouk0IlJCVC5EpHQdysSx+HsCloWn21Wm04hIKcqf2D3vO8jNMZxGREqCo6ADXn/9daZMmcLGjRtxuVycf/75PPXUUzRv3vy0j1uzZg0PP/wwK1asoGLFitx222088sgjWJYVtPAiYlbCR7uK/JjHt07gGY+Hbyu1pO9UCyj6c4hI6fr9ez0eFpzJ92w5FpZryAWHN/HAq+N4/Z83BSOeiISQAkcuFixYwO2338706dOZPHkyDoeDK6+8kgMHDpzyMYcOHaJfv35Uq1aN2bNn89JLLzF8+HBGjBgR1PAiEl6cfi93754JwFu1NddCJBq9XasHAIN2zYBAwHAaEQm2Akcuxo8ff8Ln7777LnXq1GHJkiX06nXyNwdjx44lOzubUaNGERcXR/PmzdmwYQMjR47k3nvv1eiFSJS6eu+PJLoz+TW+NrMqJpuOIyIGjKt2Aa9u+oKWWds5mrIaf5MWpiOJSBAVec7FkSNH8Pv9JCQknPKYpUuX0q5dO+Li4vJv69q1K3v27GHbtm3FCioi4e++nXm7846o3RP0RwaRqOS2ORldszOgZWlFIlGBIxd/9thjj3HOOefQtm3bUx6Tnp5OYmLiCbdVrVo1/7569eqd9HEpKSlFjROSIuV1RBudt+KIL/SRFxxMoe3hTWQ4yvJF9Q4lmElEQt17iV15ZPsUHMvmsW7FMrzlEkxHigr6PReeQvG8NW7c+JT3FalcPP744yxZsoRp06Zht59+06s/X/oUOHZd5ekuiTpd0HCRkpISEa8j2ui8FVMRJnYO2ZU3avFBzc5k22NKKpGIhIGdsZWZVKU1V+9bRtMtq3BffbvpSBFPv+fCUziet0JfFjV06FC+/vprJk+efMqRh+OqVatGenr6Cbft27cP+H0EQ0SiR62cDPqn/4gXGyNrdTcdR0RCwPFFHZyzJmlZWpEIUqhy8eijjzJu3DgmT55MkyZNCjy+bdu2LF68mJyc339YzJkzh5o1a1K3bt3ipxWRsHTP7pk48DO+aht2xlY2HUdEQsDCCk3wNWyGlXUIx4JppuOISJAUWC4eeughxowZw+jRo0lISCAtLY20tDSOHDmSf8ywYcPo06dP/uf9+/cnLi6OQYMGsXbtWiZPnswbb7zBoEGDtFKUSJSJ9bm5Y/dsAIbXvtRwGhEJGZaF+9JrAXBNHwt+n+FAIhIMBZaL0aNHc/jwYfr27UvTpk3zP4YPH55/TGpqKlu2bMn/vEKFCkyYMIE9e/bQuXNnHn74YQYPHsy9995bMq9CRELWDWkLqeI9wrJyDVhcPryuGxWRkuVrfRH+qjWxpe3CvmKR6TgiEgQFTujOzMws8ElGjRr1l9uSk5OZOnVqsUKJSIQIBBiyK+9yh+G1tPysiPyJ3YGn5wBiPn8L19QvyT7/YtOJROQMFXkpWhEJLwkfFX5Fp2DrnLmWc7J2sseVwNhqFxrLISKhKeGjXZTxnstWRzwVN/7KpW/8wJIKwRnhzBxYKyjPIyJFU+RN9ERECmvIzrxRi3cTu+Kx6W8ZIvJXWY5Y3k3sBsA/dnxrOI2InCmVCxEpEQ2y0+id8TO5loP3EruajiMiIWxErR64LTv99v1Eg+w003FE5AyoXIhIiRi8cwY2Avy3envSXRVMxxGREJYaU5Ex1TtgI8B9O7UsrUg4U7kQkaAr5z3KwNQfgGMTuUVECvCf2pcBMHDPXCp6jhRwtIiEKpULEQm6W1PnUd6Xw7wKZ/FLuXqm44hIGFhTNonpFVtQxp/LXbtnmY4jIsWkciEiQWUF/AzeNQOAt7RpnogUwetJlwNw787puPwew2lEpDhULkQkqC7LWEnj7DS2xlRhcpXWpuOISBiZVTGZlWXqUMNzkBvSFpqOIyLFoHIhIkF1fDLmyFo98Fv6ESMiRWBZ/Ccpb+7FP3Z8B4GA4UAiUlT6zS8iQZN8ZAddM9eQZYvhw5qdTMcRkTD0VbV27HRVJPnoLi7d/4vpOCJSRCoXIhI09+6aDsAnNS4m01nGcBoRCUcem4Phx+ZrPbDjO8NpRKSoVC5EJCgqeQ5zU9oCAEbU1vKzIlJ8o2t25rA9li6Za2h5eKvpOCJSBCoXIhIUd+yeQ5zfw7RKLdgQn2g6joiEsYPOMnxQszNwbO6FiIQNlQsROWMOv5d7dn8PwPBaWn5WRM7cW7V74sXGtemLqZ2TYTqOiBSSyoWInLEr9/1EUu5+1sXVZEalc0zHEZEIsD22KuOqXYADP0OOzecSkdCnciEiZ+z+Y8vPjqjdk4CWnxWRIPlP7bxlae/YPZty3qOG04hIYehdgIickfMPbaLdoRQy7fF8Vv1i03FEJIIsL9+AuRWaUcGXze17fjAdR0QKQeVCRM7IkJ15lyt8kNiZLEes4TQiEmleP7ap3n07p+Hwew2nEZGCqFyISLHVyD3AgL1L8GHxdq3upuOISAT6rnJLfotPpE5uBv33/mg6jogUQOVCRIrt7t2zcAV8TKpyPttjq5qOIyIRKGDZeKN2L+DYsrSBgOFEInI6KhciUiwxPjd37p4FwHBtmiciJejz6heR5ixP6yNb6ZS51nQcETkNlQsRKZbr0hdTzXOIn8vWZX6Fs0zHEZEIlmt3MbJWDwAe0KZ6IiFN5UJEii4Q4L5jy88Or9UTLMtwIBGJdO8kduOozcVl+1fSLGun6TgicgoqFyJSZJdn/My5WdvZ40rgy2rtTMcRkSiQ4SrHJzU6AvB/O6YaTiMip6JyISJFEwjw+LaJALyWdDm5dpfZPCISNd6s3Qs/FjelLaB6bqbpOCJyEioXIlIkXTLXcMHhTexzlOW9xC6m44hIFNkYX4PJVVoTE/AyaPf3puOIyEmoXIhIkQzdNgmAN5N6cdSuTfNEpHQd31Tv7l0zifflGE4jIn+mciEihdb+4Ho6Z64l0x7P28dWbhERKU2Lyjfhx3INqew9wi2p803HEZE/UbkQkUI7PmoxslZ3DjniDacRkahkWbyWdDmQN7HbFvAbDiQif6RyISKF0urwFnrt/4UsWwxvHtstV0TEhIlV27A5tiqNctLos2+56Tgi8gcqFyJSKMdHLd5N7EqGq5zhNCISzfyWLf+PHA/s+NZwGhH5I5ULESlQ86ydXLVvGTmWM38ypYiISR/XuIT9jjK0P5RCu4MbTMcRkWNULkSkQI9umwzAhzUvITWmouE0IiKQ5Yjl3cRuAPxjx3eG04jIcSoXInJaDY+mcl36IjyWnVfr9DYdR0Qk39u1upNrObhy3080yE4zHUdEULkQkQI8vGMKdgJ8Xv0itsdWNR1HRCRfakxFxlTvgI0A9++YajqOiKByISKnkZSzj1tS5+PD4t91rjAdR0TkL944NrF7YOpcKnkOG04jIioXInJKD+74FlfAx9hqF5ISX9N0HBGRv1hTNolplVoQ73dz1+5ZpuOIRD2VCxE5qeq5mdy+Zw4AL9XpYziNiMipvX5sU73BO2cQ43MbTiMS3VQuROSk/rHzO+L8HiZWOZ9fy9YxHUdE5JRmJySzskwdangOcn36ItNxRKKayoWI/EUlz+H8ywterNPXcBoRkQJYVv7oxQM7vsUK+A0HEoleKhci8hdDdk6nnC+H6RVbsLx8A9NxREQK9FW1C9npqkjzo7u5dP8q03FEopbKhYicoJz3KPfumg7Ai3U1aiEi4cFrczC89qWANtUTMUnlQkROcM+umVT0HmVehbNYkHCW6TgiIoX2fmIXDtlj6ZK5BtvWDabjiEQllQsRyRfny+X/duZtRPVC3SvNhhERKaJDjng+qNkZAOe0rwynEYlOhSoXCxcu5LrrrqNZs2YkJCTwxRdfnPb4bdu2kZCQ8JePmTNnBiW0iJSMO/bMoZrnEMvKNWBmxbNNxxERKbK3al+KFxuOH2dj7Us1HUck6hSqXGRlZdG8eXNeeukl4uLiCv3kX3/9NevXr8//6NixY7GDikjJcvk9PLT9GwCer3slWJbZQCIixbAjtgpfVmuH5ffjmvyZ6TgiUadQ5aJHjx48+eST9O3bF5ut8FdSVapUierVq+d/uFyuYgcVkZJ1a+o8arkP8EuZOnxbuZXpOCIixfZcvX4EbDYc86dipe40HUckqpTonIubb76ZRo0a0bNnTyZNmlSSX0pEzoDD7+WR7VMAeKluHwKWpmOJSPhKia+J96JL80YvJn5sOo5IVHGUxJOWLVuWZ599lgsvvBCHw8F3333HwIEDGTVqFNdee+0pH5eSklIScUpdpLyOaBO55y2+wCOuS19M/Zy9rI+ryddVLyiFTCIiJWt9i440XzAdx+JZbDznInKq1TIdybjI/T0X2ULxvDVu3PiU95VIuahcuTJDhgzJ/7xVq1bs37+fN99887Tl4nRBw0VKSkpEvI5oE9HnbcGu095tC/h5dPtkAF6ucwV+jVqISASo1+ZCvF364Jo5gcYrZpMz5BnTkYyK6N9zESwcz1upvYto3bo1mzdvLq0vJyKF1G/vMpod3c3WmCqMqd7BdBwRkaDxXHETAVcMjp/mYduy3nQckahQauVi9erVVK9evbS+nIgURiDA0G0TAfh3nSvw2kpkMFNExIhAQmU83foB4Br/oeE0ItGhUO8kjhw5kj/q4Pf72blzJ6tWraJixYokJSUxbNgwli9fzuTJeZdWjBkzBqfTSYsWLbDZbEybNo3Ro0fz9NNPl9gLEZGiuzzjZ1pmbWe3K4FPamipaBGJPO7LrsM5ezKOVT9i27Aaf5NzTEcSiWiFGrn4+eef6dixIx07diQ7O5sXX3yRjh078sILLwCQmprKli1bTnjMq6++SufOnenSpQtff/01I0aMYPDgwcF/BSJSPIEAjx8btXgtqTe5di0VLSIRqFwCnp4DAIj5ejQEAoYDiUS2Qo1cXHzxxWRmZp7y/lGjRp3w+Q033MANN9xwRsFEpGR1yVzDBYc3sddZjvcTO5uOIyJSYtyXDsA5czz2db9gX7scX/L5piOJRCwtCyMSpY6PWrxZuxdH7bFmw4iIlKT4srgvuw4A17gPNHohUoJULkSiUPuD6+mU+RuZ9nhG1upuOo6ISInzdL8Kf/mK2Df/hn3lYtNxRCKWyoVIFBq6bRIAb9fuwSFHwZvsiYiEvZg4PFfcCIDr6w/A7zccSCQyad1JkShz3uEt9Nr/C0dsMbxV61LTcURESkTCR3/dQDTG15p1MZVI2rGJO16ZwLhqFxbruTMHardvkVPRyIVIlHns2KjFu7W6keEqZziNiEjpybW7eL5u3r4XT28Zh93vM5xIJPKoXIhEkeZZO7lq3zJyLCf/qd3LdBwRkVL3cY2ObIqtxlnZe7ghfaHpOCIRR+VCJIocH7X4oGYnUmMqGk4jIlL6vDYHz9S7CoAnt47H6fcaTiQSWVQuRKJEw6OpXJu+GI9l59U6vU3HEREx5r/VO7A2PpH6OXsZuOcH03FEIorKhUiUeGT7FOwE+Kz6ReyIrWI6joiIMX7LxtP1+gPwz20TifW5DScSiRwqFyJRIClnH7ekzceHxb/r9DEdR0TEuAlV2/Bz2brUch/grt2zTMcRiRgqFyJR4MEd3+IM+Piq2oVsjK9hOo6IiHEBy8ZT9QcA8Oj2yZTx5hhOJBIZVC5EIlz13Ezu2D0HgJfq9DWcRkQkdHxXqSVLyjeimucQQ3ZNNx1HJCKoXIhEuH/s/I7YgIcJVc5nTdkk03FEREKHZfFE/WsAeHDHN1TwZBkOJBL+VC5EItmRg9y9ayYAL9a90mwWEZEQNKdiMrMTkqnoPcoDO78zHUck7KlciEQw17SxlPXnMq1SC1aUq286johISHry2NyL+3ZOo4r7kOE0IuFN5UIkQlkZ6TinjwXg+br9DKcREQldSyo05ttKLSnny+GR7VNMxxEJayoXIhHKNfY9LHcuX1W9gMUVmpiOIyIS0o6vHHXP7u+pmXvAcBqR8KVyIRKBbBvX4Fw8k4DTydAG15uOIyIS8laWq8e4qm2J83sYum2i6TgiYUvlQiTS+P3EfDECAM+l17ItrqrhQCIi4WFYvavxY3HHnjnUzd5rOo5IWFK5EIkwjsUzsW/+DX9CZdy9bzAdR0QkbPxWpjZfVO+AK+DjX9smmI4jEpZULkQiSW42rq/eA8A94O8QG284kIhIeHm23lV4LDu3pM6jydHdpuOIhB2VC5EI4vr2v9gy9+Gr3xRv+x6m44iIhJ3NcdX5qMYl2Anw5NbxpuOIhB2VC5EIYe1Lxfnd/wDIveFesOnbW0SkOF6oeyW5loPr0hdzzpHtpuOIhBW9+xCJEK6v3sPyuPFc0AV/k3NMxxERCVs7YyvzTmI3AJ7eMs5wGpHw4jAdQETOnG3Dapw/zibgdOG+9i7TcUREwt7Ldftwx5459M1YTptDm1hWvmH+fQkf7Sqxr5s5sFaJPbdIadDIhUi48/uJGXNs6dnLriNQubrhQCIi4S/dVYERtfPmrg3bMtZwGpHwoXIhEuYci2Zg37Ief0IV3JddZzqOiEjEeDWpNwftcfQ4sJqLM38zHUckLKhciISznKO4xr4PgPuaO7X0rIhIEB1wluU/SZcB8OyWsRAIGE4kEvpULkTCmOubMdgyM/A1aIa3XTfTcUREIs6btS9ln6MsFx1cT48Dq03HEQl5KhciYcrauwfntC8ByL1RS8+KiJSEw454XqlzBXBs7oVGL0ROS+9GRMJU3tKzHjztuuFvlGw6johIxBpZqzt7XAm0ObyZPvuWm44jEtJULkTCkG39KpxL5xBwxeAecKfpOCIiES3bHsOLdfoCMGzrWGwBv+FEIqFL5UIk3Pj9xHzxx6VnqxkOJCIS+UYndmZbTGXOydrJgPQlpuOIhCyVC5Ew41gwHfu2DfgrVcV92fWm44iIRAW3zclz9a4C4Omt43D5PYYTiYQmlQuRcJJ9FNe4Y0vPDrgTYmINBxIRiR6fVr+Y3+ITaZydxiPbp5iOIxKSVC5Ewojrmy+wHdyPr2FzLT0rIlLKfDY7g5v8DYCh2ybRNGu34UQioUflQiRMWHv34Jz+FXBs6VnLMpxIRCT6zEtoxgc1OhET8DJqwwdYmtwtcgKVC5EwEfO/UXlLz7bvjr9hc9NxRESi1mMNryfNWZ6OB9dxW+o803FEQorKhUgYsK1bieOneQRcsbgH/N10HBGRqHbAWZYHGt0MwMubxlDNfdBwIpHQoXIhEur8vvylZ92XX0+gkpaeFREx7ctq7ZhesQWVvFm8tvEz03FEQobKhUiIc8yfhn37RvyVquHpda3pOCIiAmBZDG4ykKM2F9enL6Znxi+mE4mEBJULkVCWnYVr3GgA3NfcpaVnRURCyNa4agyrdzUAI1I+It6XYziRiHkqFyIhzDX5c2yHDuBrlIz3wi6m44iIyJ+8UbsXP5etS/2cvTy5dbzpOCLGFapcLFy4kOuuu45mzZqRkJDAF198UeBj1qxZw2WXXUaNGjVo1qwZL7/8MoFA4IwDi0QLK20XzhnjAMi9cYiWnhURCUE+m517mtyOD4v/2zGVcw9vNR1JxKhClYusrCyaN2/OSy+9RFxcXIHHHzp0iH79+lGtWjVmz57NSy+9xPDhwxkxYsQZBxaJFjFfvoPl9eDp0BN/g7NMxxERkVP4qXxDRtTqiQM/72z4AJv2vpAoVqhy0aNHD5588kn69u2LzVbwQ8aOHUt2djajRo2iefPm9O3bl/vvv5+RI0dq9EKkEOy//Yxj+XwtPSsiEiaeqt+fHTGVaHN4M4N3zTAdR8SYEplzsXTpUtq1a3fCKEfXrl3Zs2cP27ZtK4kvKRI5/D5cY44tPdv7BgIVqxgOJCIiBTniiGNI44EAPLv5K5Jy9hlOJGKGoySeND09ncTExBNuq1q1av599erVO+njUlJSSiJOqYuU1xFtQuW8VV4xjzrbN+EuX4nfGrcmcMa54oOSS0RETu+bKufxdZU2XL1vGW+lfEK/sx8o8ny5kvxdFCq/56RoQvG8NW7c+JT3lUi5ALD+9M10/HKoP9/+R6cLGi5SUlIi4nVEm5A5b0ePEP/mFAD8N91Lo+bJZ/6cC3ad+XOIiEih/F/jW+l24FeuyFhBv33LmFC1bZEeX1K/i0Lm95wUSTietxK5LKpatWqkp6efcNu+fXnDg8dHMETkr1yTP8N2OBNfk3Pwtu1sOo6IiBTRnpiKPN7gOgDeTPmE8t6jhhOJlK4SKRdt27Zl8eLF5OT8vpnMnDlzqFmzJnXr1i2JLykS9qy0nThnfE3Assi94V4tPSsiEqbeS+zCovKNSXRn8vzmL03HESlVhSoXR44cYdWqVaxatQq/38/OnTtZtWoVO3bsAGDYsGH06dMn//j+/fsTFxfHoEGDWLt2LZMnT+aNN95g0KBBp70sSiSaxfxvFJbPi7dDT/z1m5qOIyIixRSwbNzT5HY8lp17ds+k3cENpiOJlJpClYuff/6Zjh070rFjR7Kzs3nxxRfp2LEjL7zwAgCpqals2bIl//gKFSowYcIE9uzZQ+fOnXn44YcZPHgw9957b8m8CpEwZ1+zHMeKhQRitPSsiEgkWFM2iVeSegMwav0HOP1ew4lESkehJnRffPHFZGZmnvL+UaNG/eW25ORkpk6dWuxgIlHD5/196dkrbiKQUNlwIBERCYYX6l7JgL1LOPvoTh7a8Q0v1r3SdCSRElcicy5EpPCc34/HvnML/io18PQcYDqOiIgESY7dxaAmtwPwz60TaXQ01XAikZKnciFikG1bCq6x7wOQe9N94IoxnEhERIJpTsVkPq7RkdiAh1EbPoBjS/OLRCqVCxFTcrOJHfUslteDp3MffK3am04kIiIl4JGGN7DXWY7OmWu5OW2+6TgiJUrlQsSQmC/exrZnO77EeuReP8h0HBERKSH7neV4sOFNALy68QuquA8ZTiRSclQuRAywL5uLc+43BJxOcu95AmJiTUcSEZESNKZ6B76veDaVvUd4ZdMXpuOIlBiVC5FSZmWkE/vRqwC4r70Hf52GhhOJiEiJsywGN/kb2TYnN6ctoNv+1aYTiZQIlQuR0uT3Efve81hZh/GeeyGebv1MJxIRkVKyOa46z9a9CoC3N3xInC/XcCKR4FO5EClFzm/GYF/3C/4Klci54zHQjvUiIlHl9aTLWFUmiYY56fxr2wTTcUSCTuVCpJTYNq7BNeEjAHL/PhTKJ5gNJCIipc5rc3B30zvwY/HAju8458h205FEgkrlQqQ0HD2St+ys34+717X4zmljOpGIiBiytHwjRtbqjjPg4531o7EF/KYjiQSNyoVISQsEiPn0DWz7UvHVbYK7/x2mE4mIiGFP1B/ATldFLji8ibt2zzQdRyRoVC5ESphj0fc4F88k4IolZ9AT4HCajiQiIoYddsRzX+PbAHh+85ck5u43G0gkSBymA4hEMittJzGf/geA3JvvI1Aj6ZTHJny0q7RiiYhICJhc9XwmVjmfK/f9xJspnzDg7H+YjiRyxjRyIVJSvF5iRz2HlZONp00nvBf3Mp1IRERCzP2NbuGQPZZ++36iz96fTMcROWMqFyIlxDX+Q+xb1uGvXJ3cgQ9q2VkREfmLXbGV+Vf9awF4K+VjyM4yG0jkDKlciJQA+9oVOL/7LwHLRs7d/4Qy5UxHEhGREPVOrW78WK4htd0HiPnyHdNxRM6IyoVIsB3OJObdF7ACATx9b8bfpIXpRCIiEsL8lo27m96B27LjnDMFx7yppiOJFJvKhUgwBQLEfvgKtsx9+BqfjbvPzaYTiYhIGFhdtg6Dm/wNgJhPXse2cY3hRCLFo3IhEkSO2ZNxrFhIIL4MOXf/C+xakE1ERArno5qdcHe/CsvrIfatJ7D27zUdSaTIVC5EgsS2cwsx/30bgNzbHiRQpYbhRCIiEm7c1w3C26wVtoP7iX3rX+DONR1JpEhULkSCwZ1LzKhnsTxuPBf3wntBF9OJREQkHDkc5Nz7NP6qNbFvWU/MR69CIGA6lUihqVyIBIHrq3ex79yMv3ptcm8aYjqOiIiEs7IVyLn/eQIxsTgXfY9z2lemE4kUmsqFyBmyr1yM6/vxBOwOcu55AmLjTUcSEZEw509qQM6djwPg+vJd7KuXGk4kUjiabSpyBqzMDGJHvwSAu/8d+Os3NZxIRETCWcJHu/7wWUOerHsVT24bT+5/nqZd62dIia9ZzGeOJ7NxMBKKnJ5GLkSKy+8n5r0XsQ4fxJvcGs+l15hOJCIiEebZev2YUOV8EnxHGf/r65TzHjUdSeS0VC5Eisk57Ssca34iUK4CuXc+DjZ9O4mISHAFLBsDz7qb1WVq0+zobj77bSS2gN90LJFT0rshkWKwbVmPa9xoAHLueJRAQmXDiUREJFIdccRx1dkPkuEoS++Mnxm2ZZzpSCKnpHIhUlQ5R4kd9SyWz4u7Wz98LdubTiQiIhFuS1w1rk8eghcbQ7dPYkD6YtORRE5K5UKkiGK+GIEtbSe+2vVxX3u36TgiIhIlZlc8m4ca3QjAB+veo+XhrWYDiZyEyoVIETh+nINz3ncEnC5y73kCXDGmI4mISBQZUasnH9a4hHi/m/G/vk5V90HTkUROoHIhUkjWvlRiPn4VgNzrB+Gv3cBwIhERiTqWxb1NBrK4fGPq5Gbw1Zo3cfq9plOJ5FO5ECkMr5fYd5/HOpqFt1UHvF36mk4kIiJRym1zMiD5fna6KnLxwfW8kfKJ6Ugi+VQuRAri9xHz3vPYN6zGn1CFnNsfBssynUpERKJYakxF+p/9D3IsJ3ftmc2du2aajiQCqFyInJ7fT8wHr+D8cQ6B2Hhy7n8OyiWYTiUiIsJP5RtyZ9M7AHhz46d0zPzNcCIRcJgOIBKyAgFcn7+Fc8E0Aq5Ysh98ifJzy8HcXaaTiYiIADCmxkWcm7WdB3d8y5dr3uSC1s+yPbaq6VgSxTRyIXIygQCuL9/BNWsiAaeTnP97Dn+TFqZTiYiI/MXQBtcxrVILqnoOM37168T7ckxHkiimciFyEq6JH+Oa+iUBu52ce4fhSz7fdCQREZGT8ls2bmp2LxviatAyazsfrHsPAgHTsSRKqVyI/Inzu//hmvgJActGzt1PaAduEREJeZnOMlx19gMcsscyYO+PDN0+yXQkiVIqFyJ/4Jw5gZgv3wEg945H8bXtZDaQiIhIIa0rU4ubmw3Gj8WzW8bSe99y05EkCqlciBzjmPcdMZ+9CUDOrf/Ae1FPw4lERESK5tsq5/FE/QEAfPrbSJpl7TScSKKNyoUIkLBmKTEfvgJA7vWDtUmeiIiErZfr9OHLqhdS3pfDhNWvU9FzxHQkiSIqFxL17MvnU2/iB1iBALlX/Q3PpQNMRxIRESk+y+KOs+5kRdl6NMpJY8za4dj9PtOpJEqoXEhUs69eSuzIZ7ACftyX34Cnz82mI4mIiJyxbHsMV5/9D9Kd5el+4Fde2vxf05EkSqhcSNSyrVtJ7Jv/wvJ6SG/TFfeAv4NlmY4lIiISFDtiq3BN8v24LTv/2DkV58wJpiNJFCh0uRg9ejQtWrSgevXqXHLJJSxatOiUx27bto2EhIS/fMycOTMooUXOlG3jGuL+MxTL48ZzyeXs6nGtioWIiEScBQlnMaTxbQDEfPYmrvEfaQ8MKVGFKhfjx4/nscce48EHH2TevHm0bduWAQMGsGPHjtM+7uuvv2b9+vX5Hx07dgxKaJEzYduWQtxrj2DlZONp143c2x5QsRARkYj1QWIX7mxyBwHLhmvSJ8R8/Dr4vKZjSYQqVLl4++23ueGGG7j11ltp2rQpr7zyCtWrV+fDDz887eMqVapE9erV8z9cLldQQosUl23nFuJeeQjraBbe8zuS+/fHwGY3HUtERKREfZjYmZz7niXgdOH8YQqxI54Cd67pWBKBCiwXbreblStX0qVLlxNu79KlCz/++ONpH3vzzTfTqFEjevbsyaRJ2ilSzLJSdxL77wexDh/E2+ICcu55AuwO07FERERKhe+8DmQ/8hqBMuVwrFhI3L8fhCOHTMeSCFPgO6uMjAx8Ph9Vq1Y94faqVauSnp5+0seULVuWZ599lgsvvBCHw8F3333HwIEDGTVqFNdee+0pv1ZKSkoR44emSHkdkcSZmUGTT/+N7dB+DtdtyqZetxDYsvWEYwp33uJLJJ+IiEhJS0lJASuW2JseouGYN3Cl/Ir9qbvYdP3/4alQyXQ8OYVQfF/ZuHHjU95X6D/bWn+6Jj0QCPzltuMqV67MkCFD8j9v1aoV+/fv58033zxtuThd0HCRkpISEa8jnCV8tOuEz2vmHmDOz//BlbOfReUb0yvpEbJ+jDWUTkRExIz89yeNG+M5qxn2Vx8hbtdWmn/xKtkPvUKgVj2j+eSvwvF9ZYGXRVWuXBm73f6XUYp9+/b9ZTTjdFq3bs3mzZuLnlDkDFRxH2L6Ly/QKCeN5WXr0bvFI2Q5VCxERCS6BSpVI/ufw/E1OQfb/r3EPz8E24bVpmNJBCiwXLhcLlq2bMmcOXNOuH3OnDlccMEFhf5Cq1evpnr16kVPKFJMCZ4spv3yIs2P7mZ1mdr0OvcxDjl0WZOIiAgAZcqR/fCreM/rgJV1mLh/P4h9xULTqSTMFWq1qMGDBzNmzBg+/fRT1q9fz6OPPkpqaioDBw4EYNiwYfTp0yf/+DFjxjB27FjWr19PSkoKw4cPZ/To0dx5550l8ypE/qSc9yjfrXqZllnbWR9Xk0tbDGW/s5zpWCIiIqHFFUPOvcPwdLoCy+Mm9q0ncPzwjelUEsYKNefiqquuYv/+/bzyyiukpaXRrFkzvvrqK+rUqQNAamoqW7ZsOeExr776Kjt27MBut9OwYUNGjBhx2vkWIsES58tl8upXaXt4E1tjqtDz3KGkxSSYjiUiIhKa7A5yb3uAQEIlXBM/IfajV8k9uB9Pn5u1D5QUmZWZmaltGoMoHCfeRBSPm0UPP0iPA6vZ5apIp1ZPsiWumulUIiIixmUOrFXgMY7Zk4j59E2sgB9P5z7k3nK/9oMyKBzfV2qRfzHizys6BUMFTxaf/TaSyw6sJt1Znh7nPq5iISIiUgTeLn0JlK9E7DvP4JwzGevQAXLu/he4YkxHkzBRqDkXIqHu/EObWLb8n1y2fyX7HWXoee5Q1pdJNB1LREQk7PjOv5jsh18jEF8Wx/L5xL3yMGQdNh1LwoTKhYS3QID7dkxl3s/DaJCzl+Vl69HuvGdYXbaO6WQiIiJhy9+0Bdn/fAt/xSrYN6wi7oX7sPbvNR1LwoDKhYStip4jjP/1dV7f9DmugI/htXpw8XlPsym+huloIiIiYc9fuwHZT7yNv2Yd7Du3EPfsYKzd20zHkhCnciFhqd3BDSz/6XH6ZKzggCOeq5P/j380vhW3zWk6moiISMQIVK7O0X8Nx9coGdv+dOKfG4Jt4xrTsSSEqVxIWLECfh7aPoU5Pz9LndwMlpZryPmtX2BS1Tamo4mIiESmshXIfuQ1vC3bYWUdIu7lB7CvXGQ6lYQolQsJG1Xch5i8+lVe2vw/HPh5vfZlXNLqSbbFVTUdTUREJLLFxJJz37N4Ol6G5c4l9s1/4Zj3nelUEoK0FK2EhYszf+PztW9Ty32ADEdZ/nbWXXxb5TzTsURERKKH3UHu3x4mkFAZ1+TPiP3g3+RmZuC54iZttif5VC4kpNkCfh7bNomntn6NnQALyzfhxub3sjO2suloIiIi0ceycF99O4EKlXB9/hYxX3+AlZmB+8Z7wa63laJyISGsem4mn/42kq6ZeRPHXqrTh6fq9cennUJFRESM8nTrh79CRWLfeR7XrInYU34ld+BD+BucZTqaGKZyISGpy4Ff+XTtSGp4DpLuLM9tze5hRqUWpmOJiIiErYSPdgX5GRtz0dmPMXvn+9i3byTumXvwdOuH++rbIa5MkL+WhAtN6JaQYvf7eHrLOKb98hI1PAf5IaEZrc9/QcVCREQkBC1IOIujL3yE+7LrwLJwfT+e+KG3Yl8+33Q0MUQjFxIyEnP38/nat+l4cB1+LJ6pexXP1euH31IHFhERCVkxcbivvRtvu27EfPQa9s2/EffWE3jP60DuzfcTqFTNdEIpRSoXEhJ6ZvzCx+tGUdVzmD2uBG5pNog5FZNNxxIREZFC8tdpRPYTI3DOmoRr3GgcKxZiX7sC99V34Ol2JWjOZFTQn4TFKIffywub/se3q/9NVc9hvq94Nq3Pf0HFQkREJBzZ7Hi6X8XRFz/G2/pirJxsYr4YTtwzg7BtSzGdTkqByoUYk5Szj9krn+ORHVPwYuOf9a/hshaPku6qYDqaiIiInIFApWrk3Pcs2fc/h79SVexb1hP39F24/jsScrNNx5MSpMuixIje+5bz4bp3qeTNYqerIjc2H8LChKamY4mIiEgQ+c67iKPNzsM1/gOc30/ANe0rHD/NJfeWf+A790LT8aQEqFxIqbL27sE1+TMm/vodAN9WasnfzrqbDFc5w8lERESkRMTF475xCN723fMmfG9LIe71x/C06YT7piEEErQxbiRRuZBSYWWk4Zr8OY7532H5fHgsO/+sfy3/SepFQKtBiYiIRDx//bPIfmoUzu/H4/r6Q5zLfsCxZhm5A+7E2+kKsOn9QCRQuZASZWWk45ryOY5532H5vAQsG5723Wnh7UlKfE3T8URERKQ02R14Lr0G7/kdifn0DRy/LCH2k//gWziD3IEP4q/dwHRCOUMqF1IirP3pOL8Zg3Put1heDwHLwtOuG+4+NxNIrEtK0HcJFRERkXARqFKDnH+8iP2nucR8Phz7xjXEPfl3PL2uw933FnDFmI4oxaRyIUFlHdiH85svcP7wze+l4oIuuK+8lUBiXdPxREREJFRYFr42nTiafD4xY9/HMWcyrm++wPHjHHJvewDf2eebTijFoHIhQWFlZuSNVPwwGcvjAcDTtjOevrfgr13fcDoREREJWfFlyb31H3g69CDmo1ex79xC3CsP5V3xcMNgAuUrmk4oRaByIWfEOrgf57f/xTl7EpbHDYD3/I64r7wNf5KumxQREZHC8TdKJnvY+zinfYlr4ic4F8/E8fMiPB164Ol6JYFa9UxHlEJQuZBisQ4d+L1UuHMB8La+GPeVt+Kv08hwOhEREQlLDgee3jfibdOJmM/exLF6Ka5ZE3HNmoj3rJZ4ul2Jr9VF4NBb2FClMyNFcygT19T/4Zw5EcudA4D3vA55IxV1GxsOJyIiIpEgUL0WOQ/9G9v2TThnT8KxaAaOdStxrFuJP6EKnk698XbqTaBiFdNR5U+szMzMgOkQkSQlJYXGjSPwTfbhTFxTv8Q5cwJW7rFS0bI97n634a/XpMhPl6DVokRERKQAmQNr5f3j6BGcC2fgnD0J2+5tAATsdrznXYy3a198Z7UEyzIXtISE4/tKjVzI6R05iGvqVzhnjsfKyQbAe+6FeaWi/lmGw4mIiEhUiC+Lp/tVeLr1w75uJc6ZE7CvWIBz2Q84l/2AL7Ee3q598XToAXFlTKeNaioXclJWRhrOH77BOeNrrJyjAHhbXJB3+VPDZobTiYiISFSyLHzNWuFr1gpr/16cP3yD44cp2Hdvxf7Zm7jGvoe3fQ88XftqQz5DVC4kn5W2C8dP83D8NA/75t/yb/ee0yavVDRKNphORERE5HeBSlVxXzUQd5+bcSyfj3P2ROzrfsE5exLO2ZPwNT0XT9e+eFtfDA6n6bhRQ+Uiylm7t+FYNjevUGzfmH97wBWL79wLcPccgL/x2QYTioiIiJyGw4H3gs54L+iMbedmHLMn41w4Hfv6X7Cv/wV/hUp4O/XG06k3gUrVTKeNeJrQHWQhP/EmEMC2Y1PeCMWyufmTogACsfF4W7XHe35HfOe0hZjYEouhCd0iIiJSUsp6s7kpbQH37Pqe5KN57zm82JhcpTWjanVnTkJzMv9W23DKgoX8+8qTULkIspD8jyAQwLZl/bERirnY0nf/fleZcnhbdcDbpiO+5PPB6SqVSCoXIiIiUuICAToeXMc9u77nyn0/4Qz4APgtPpEGl12OL7k1/nqNwWY3HPTkQvJ9ZQFULoIsZP4j8PuxbVyTVyiWz8eWkfb7XeUS8LW+GG+bS/KWbjOwEY3KhYiIiJSmmrkHuH3PHP6+eza13Afybw/El8V3Vkt8ya3xNj+PQM06IbOsbci8rywClYsgM/ofgc+Lff0q7D/NyysUmRn5d/kTquBt0xHv+R3xNznHeENXuRARERETHH4vl2Ws5MsqKdjXrsC2d/cJ9/sTquBrfh6+5PPwNT/P6DyNcCwXmtAd7rwe7GtX5M2hWLEA6/DB/Lv8VWrgPb8j3jaX4G/QDGw2g0FFREREzPPaHEyuej65A/sCYO3dg33NcuxrV2D/7WdsmfuwLZqBc9EMAPw1k/A2b51XOJq1gjLlTMYPeSoX4cTvx0rfhX3LBmxb12PbugH71g35+1AA+Gsk5RWK8zvm7ZwdIsN6IiIiIqEoULUm3k698XbqnTdPdecW7GuPlY11v2DbswPXnh0wayIBy8JftzG+5Nb4mrfG1/jsEl0AJxzpsqg/KMlLdfK3ry+sPxUJ+9b12LZtxMrO+suhvtr18R0foahVPywKhS6LEhERkVDn8Htpc3gzXQ6socuBNbQ7tAHXsUnhALmWg8UVGjM7IZlZFc/mp3IN8JXgpedFfj9pgMrFHxgrF34/VtpO7Fs3HBuNOHWR8CdUwV+vCb56TfDXb4K/XlMCCZVLLHdJUbkQERGRcBPny+Wig+vpeuBXOh9YQ6sj27Dx+1vpg/Y45iU0Y0W5eqyLT2RdfC1S4mqQYw/OapzhUC50WVRp+3OR2LIe27aUEy5tyj/0eJGo3xR/vSb46zUJyyIhIiIiEgmy7TF8X6kF31dqAUAlz2Euyfwtf2SjafYershYwRUZK/If48dia2wV1sXXYn18zfzSsS4+kQxX5M3f0MjFHwTzr+kxPje1c/eTlJtBnZx9vFPvQF6hOFWRqFgFf72meSMSUVAkNHIhIiIikaZ2TgYdM38j+ehOmh7dw1lZu2iUnYYD/0mP3+coy7oytVgXX/NY+UhkXXwi22Kr4Lf+uhCPRi4iVSBAVc8h6uRk5JeHOrkZJOX/bwY1PAdPfMz63/+ZXyT+OCJRoVLpvgYRERERCaqdsZUZU+OiE25z+r00zE6j6dHdnPWHj6ZHd1PFe4SLDq7nooPrT3hMts3Jhriax0Y68kY51sfXBHcVcMWU5ksqMo1c/MHxv6bH+tzHSkMGSbn7TigRx2+PDXhO+1xebOyMqcSO2Mpsj6nChviarChXn+Xl6pPuqlAaL+eMlWQ71siFiIiIRLVAgET3gbyykbX7hPLxx03+/sh92fW4r72rlIMWTaFHLkaPHs1bb71FWloaZ511Fi+++CLt27c/5fFr1qzh4YcfZsWKFVSsWJHbbruNRx55BCuEVzKau2IYjbNTqeY5VOCxGY6ybI+tzI6YymyPrcLOmEpsj6lyrExUZk9MxZMOZ4mIiIiIYFnsjqnE7phKzK549gl3lfMe5ayje04oHE2P7qZ+7fqGwhZeocrF+PHjeeyxx3jttde48MILGT16NAMGDGDJkiUkJSX95fhDhw7Rr18/2rdvz+zZs0lJSWHw4MHEx8czZMiQoL+IYKnmOUg1zyHclp0dMZXzRx2OF4YdsVXy/jemMlkOrWksIiIiIsF32BHPsvINWVa+4Qm3Z3YI/TkXhbosqmvXriQnJ/PWW2/l33beeefRt29fnnrqqb8c/8EHH/D000+zYcMG4uLiAHjllVf48MMPWbt2bciOXlw4YimZjnjSXBUIaNRBl0WJiIiIhJCImNDtdrtZuXLlX0YcunTpwo8//njSxyxdupR27drlFwvIKyjPP/8827Zto169emeWuoQsubet6QhRIxy+OURERESkaAr883xGRgY+n4+qVauecHvVqlVJT08/6WPS09NPevzx+0REREREJPIU+tqfP1/KFAgETnt508mOP9ntIiIiIiISGQosF5UrV8Zut/9lxGHfvn1/GZ04rlq1aic9HjjlY0REREREJLwVWC5cLhctW7Zkzpw5J9w+Z84cLrjggpM+pm3btixevJicnJwTjq9ZsyZ169Y9w8giIiIiIhKKCnVZ1ODBgxkzZgyffvop69ev59FHHyU1NZWBAwcCMGzYMPr06ZN/fP/+/YmLi2PQoEGsXbuWyZMn88YbbzBo0CBdFiUiIiIiEqEKVS6uuuoqXnzxRV555RUuvvhilixZwldffUWdOnUASE1NZcuWLfnHV6hQgQkTJrBnzx46d+7Mww8/zODBg7n33ntL5lUYkpuby8MPP0yDBg1ITEzkuuuuY9eu0y+xOnHiRDp16kSdOnVITEzkoosuYsyYMaWUWI4rzrn75JNP6NWrF/Xq1aNOnTr07t2bxYsXl1JigeKdt99++41bbrmFc889l4SEBF588cVSShvdRo8eTYsWLahevTqXXHIJixYtOu3xa9as4bLLLqNGjRo0a9aMl19+OX+unpSeopy3nJwc7rnnHtq3b0+VKlW4/PLLSzGp/FlRzt38+fO5/vrradq0KTVr1qR9+/Z89tlnpZhWjivKeVu3bh29e/emcePGVK9enXPPPZdnnnkGt9tdiokLVugJ3XfccQerV68mPT2duXPn0qFDh/z7Ro0axerVq084Pjk5malTp5KWlsb69et57LHHIm7UYujQoUyZMoUPPviA7777jsOHD3Pttdfi8/lO+ZiKFSvy0EMPMXPmTBYuXMiNN97IkCFDmDFjRikml+KcuwULFtCvXz8mTZrErFmzaNy4MVdffTWbNm0qxeTRrTjnLTs7mzp16vCvf/1Ll2WWkuMbrz744IPMmzePtm3bMmDAAHbs2HHS449vvFqtWjVmz57NSy+9xPDhwxkxYkQpJ49uRT1vPp+P2NhY7rzzTnr06FHKaeWPinruli5dSnJyMp988gmLFy/m9ttv5//+7/8YO3ZsKSePbkU9by6Xi+uvv57x48ezbNkyXnzxRT777DOee+65Uk5+eoXaRE/+6uDBgzRq1Ii3336ba665BoCdO3dyzjnnMG7cOLp27Vro5+rYsSNdu3Y96YaEEnzBOneBQICmTZvy4IMPctddd5VkZCE4561du3b06dOHoUOHlnTcqBYtG69GmqKetz96+OGHWbt2Ld9++21Jx5STOJNzd9xtt92Gz+fTCEYpCsZ5e/zxx1m2bBnff/99ScUsMm1DXUwrV67E4/HQpUuX/Ntq165N06ZNT7m54J8FAgHmzp3Lxo0bad++fUlFlT8JxrmDvA0mc3JySEhIKIGU8mfBOm9Sso5vvPrH8wTF23h1z549bNu2rUTzSp7inDcJDcE6d4cPH9bvs1IUjPO2efNmZs2adcLVRKGgwB265eTS09Ox2+1Urlz5hNtPt7ngcQcPHqR58+bk5uZit9t55ZVX6N69e0nGlT84k3P3R8899xxly5alV69ewY4oJxGs8yYlq7gbryYmJv7l+OP31atXr0Syyu+Kc94kNATj3E2bNo25c+cyffr0kogoJ3Em561Hjx788ssv5Obmcuutt/Lkk0+WZNQiU7n4k+eee45XX331tMdMmTLllPcVtLkgQLly5Zg/fz5Hjhxh7ty5+deCX3LJJcXKLHlK49wdN2rUKD7++GMmTpxI+fLli5RTTlSa501KjzZeDU9FPW8SOop77pYsWcLf//53Xn75ZVq3bl1S8eQUinPePvzwQ44cOcKvv/7Kk08+yRtvvMEDDzxQkjGLROXiT+65557867lPpXbt2ixbtgyfz0dGRgZVqlTJv2/fvn0FXuJks9lo0KABAC1atGDDhg289tprKhdnqDTOHeQVi+eff56xY8fqB3EQlNZ5k9KhjVfDU3HOm4SGMzl3ixcv5pprrmHo0KHcfvvtJRlT/uRMzlvt2rUBOOuss/D5fNx3333cd999OByh8bY+NFKEkMqVK//lsouTadmyJU6nkzlz5jBgwAAAdu3axfr160+5ueCp+P3+kFtGLByVxrkbMWIEL774Il999RXt2rULSu5oZ+J7TkrOHzdevfLKK/NvnzNnzgn7If1R27Ztefrpp8nJySE2Njb/eG28WnqKc94kNBT33C1cuJBrr72WRx99lEGDBpVCUvmjYH3P+f1+vF4vPp9P5SLcVahQgZtvvpknn3ySqlWrUrFiRf75z3+SnJxMp06d8o/r06cPrVu3zp/1/+qrr3L++edTr149cnNzmTFjBl9++SX//ve/Db2S6FPcc/fWW2/x7LPP8t5779GoUSPS0tIAiI2NpUKFCiZeSlQp7nlzu92sW7cOyFuXPz09nVWrVlG2bNn8EUQJrsGDB3PXXXfRunVrLrjgAj788MO/bLy6fPlyJk+eDORtvPryyy8zaNAgHnroITZu3Mgbb7zBI488oktySlFRzxvkrbvvdrvJyMggKyuLVatWAXmj8lJ6inru5s+fz7XXXsvtt9/ONddck//7zG63nzAyLCWrqOftf//7H7GxsTRv3hyXy8XPP//MM888Q9++fYmJiTH5Uk6gcnEGXnjhBex2OwMHDiQnJ4eOHTvyzjvvYLfb84/ZsmULtWrVyv88KyuLBx54gN27dxMbG0uTJk1455136N+/v4mXELWKc+7ef/99PB5P/jf9cddffz2jRo0qtezRrDjnbc+ePXTs2PGE+z/66CM6dOigZTNLyFVXXcX+/ft55ZVXSEtLo1mzZoXaePWhhx6ic+fOJCQkROTGq6GuqOcN+Mua/Me/1zIzM0sttxT93I0ZM4ajR48yfPhwhg8fnn97UlLSX/Ytk5JT1PPmcDh4/fXX2bx5M4FAgKSkJO64446QG3nSPhciIiIiIhIU2udCRERERESCQuVCRERERESCQuVCRERERESCQuVCRERERESCQuVCRERERESCQuVCRERERESCQuVCRERERESCQuVCRERERESCQuVCRERERESC4v8BxP76u8TdLW4AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "count, bins, ignored = plt.hist(s, 30, density=True)\n", "f = 1/(sigma * np.sqrt(2 * np.pi)) * np.exp( - (bins - mu)**2 / (2 * sigma**2) )\n", "plt.plot(bins, f, linewidth=2)\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "cfa55d59-9572-4fcb-9342-e237b60dc414", "metadata": {}, "source": [ "***\n", "\n", "## End" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.8" } }, "nbformat": 4, "nbformat_minor": 5 }