{ "cells": [ { "cell_type": "markdown", "metadata": { "zanadu": {} }, "source": [ "Illustration of the Glivenko-Cantelli Theorem\n", "---" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "zanadu": { "code_type": "" } }, "outputs": [], "source": [ "import pandas as pd\n", "from pandas import DataFrame, Series\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline \n", "import numpy.random as rd\n", "import seaborn as sns\n", "from scipy.stats import cumfreq\n", "from scipy.stats import norm\n", "sns.set();" ] }, { "cell_type": "markdown", "metadata": { "zanadu": {} }, "source": [ "### Generate the random sample" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "zanadu": { "code_type": "" } }, "outputs": [], "source": [ "n = 100000 ## size of the random sample\n", "nbBins = 200\n", "xx = rd.normal(0.,1., n)\n", "data = pd.DataFrame(xx, columns=['X'])\n", "vals = data['X'].values\n", "uu = np.linspace(np.min(data['X']), np.max(data['X']), nbBins) ## bins" ] }, { "cell_type": "markdown", "metadata": { "zanadu": {} }, "source": [ "### Plot the empirical CDF" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "zanadu": { "code_type": "" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAEXCAYAAAC59m+aAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3deXhU5cH+8W92EggCIRQCuCIPIpZFBREQFEQJIrvihjvaYl2K60+raEut1bbW5VXfQlVaUDGA6K9YEQGRRUFUFoHHyiZLiCEsZieZmfePmdAhBDJJJjmz3J/r8vKcmTMz94zDzeNZnonxeDyIiEjkiHU6gIiIBJeKXUQkwqjYRUQijIpdRCTCqNhFRCKMil1EJMLEOx1ApKEZY9oDnwNdrbX7fLedCUwDWgIFwHhr7WbffbcA9wMJwELgbmttmTEmBZgKdMc7SHrIWvue7zG9gJeAJsAe4HprbbbvvkeAG/H++fsn8KS1VucdS9BoxC5RxRgzHlgKZFS6awbwqrW2M/AEkGWMiTHGdAGeBPoDBmgG3Od7zGSgwFp7FnAp8LIxpp0xJhHIAu713ZeF9y8NjDGZwFXAuUAX4GJgbD29XYlSGrFLWDLGDACmAFvxFmQCcAdwNXBRpc1LrbW9jDEZwAjgMsD6PVdboBPwNoC19kNjzCt4R+JDgPettbm+bV8DXgD+CIwErvU95gdjzMd4S/sL4Cdr7XLfS0wDnjfGpPkeM9NaW+h7vteB64FZwflkRFTsEt56AROttd8YYyYBv7fW9j/extbaPcAoAGOM/13tgT3WWrffbbuAdr77tldxe8XjdlZx3x7/2621h40xuUBb32M+Oc7ziQSFil3C2Q5r7Te+5a+Am4wxL3CcEfsJnicWqLyPOwZwVXFfxe1VPe54jwn0+USCQsUu4azYb9kDxFhr767F8/wAtDHGxPgdxMzAO5r+gaP3x1fcjt99OX73fVP5McaYBCAN2F3N84kEhQ6eStSz1u4Cvse7fx5jzGWAG1gPvA9caYxpZYyJASYA7/keOs+3jjGmHXA58P/x7mNPM8Zc6NvuFmCltfag7zHXGWMaG2OSgJv8nk8kKDRiF/G6BvibMeYxoAQY69vnvs4Y8xSwCO8B2i+AZ3yPeQJ4xRjzLRAHPGCt3QJgjBkFvGSMaQzkAeMBrLUfGGPOAVYBiXiLfnoDvUeJEjGatldEJLJoV4yISIQJeFeMMaYpsAK4wlq7vdJ93fBegdcU78Ufd1pry4OYU0REAhTQiN13efQyoONxNvkncJe1tiPe07duD048ERGpqUB3xdwOTMR74cVRjDGnAMnW2s99N72BLpEWEXFMQLtirLW3wTFX61XIALL91rMJ/Eq6JOB832N0kYaISGDigDbAaqC08p3BON2xqivp3MfZtrLzgc+CkEFEJBr1w7ub/CjBKPZdeP/mqNCaKnbZHEc2wIEDhbjd4XHaZVpaE/LyCpyOETL0eRxLn8mxjveZHDx4AGs3sXXrFrZs2cLWrd+zd292Fc/wX7GxsZzUrDnNm7XgpGbNSE5JJiW5CUnJjSgojeGUjJYkpzRmf76LM05uyZ59JZxxchrx8QnszC2m06ktiY9PYFt2AZ3PaEV8fDzf/XCQzqelERMXx+ZtBzjrtBYAbNq2P+DlmmwLsPPHQs5ok0pyYlyNP8/Y2BiaN28MR+8tOaLOxW6t3WGMKTHG9PHNZncD8GGAD3cBuN2esCl2IKyyNgR9HsfSZ3Ist9vDgQMHWLNmFatXf8Hq1V+wdeuWY7ZLTEykfftTyGjXnvjkNLp26UibjLbkFcZzQfczaXpSMzZs2U/XjukArP0u97jLZ3dMp/y7XE7v4L3twOFc0tK9yzG5bhKTkgEoLQM3seDyUFRSRrnLu9OhJss1fVxxaTnlLjdud53OOq9yF3ati90YMx943Fr7JXAd3qv2muKdjOmF2j6viESWXbt2Mnv2Uj744F+sXfs1/hdFNmqUTOezz+Gk9FPo3asHHc40HChtQo+zWgPHFnXzFmmOvIdwU6Nit9ae6rec6be8FugZvFgiEs4KCgpYsOBD5s2bw9q1Xx+5PTExkVPO6Ezfvn1IbdWRKwb3JSEh4agC/+m7XKdiRwzNFSMiQbNjxzb++c83+eCDeZSUeCffTE5Opkv3C8nMHEKvC/ry3a5CunZMZ+13uSQkJDicODKp2EWkzr79dj1Tp77GkiWfHNnVcob5OWPGjqX/gEvJPuji9IxU39aFzgWNEip2Eam17du38tJLf2Xhwo8ASEhIYPDlwxgz7gYOlTU9snuFg/kOpow+KnYRqbH8/HxeeukvZGW9g8vlIjEpiZFjruHs84fSr6f3Qsa12lfuGBW7iATM4/GwYMGHPPvs0+zbl0tcXBy9BwzlvvvuJT39ZyrzEKFiF5GA7N+fx+TJj7J06RIAupzTjfsefIx8V3PSfeeGS2hQsYtItZYv/4zHH3+EvLx9NEltytDRtzHhtvHExsZqlB6CVOwiclzl5eW88MKfmT797wB073E+D/9mCtmH4oiN1e/0hCoVu4hU6aeffuLhh3/NihXLiIuLZ8jIm7j3nonExcWRfUij9FCmYheRY2zfvpW77/klP+zYTrPmLXhyyp/wpJxMXFzNJ6yShqf/lxKRo6xfv47x46/hhx3byWh/Oq9MncnPu53rdCypAY3YReSIVas+5957f0lRUREX9unP8BseoHWbDKdjSQ1pxC4iACxZsoi77ppAUVERl142lCef/jNJjZKdjiW1oBG7iLBk6adMmnQ3Llc5fQcO5+HfPKmzXsKYil0kyq1e/QUPPuAt9auuGU/vweNV6mFO//VEotj69eu4555fcLi0lAsvHsadd00iJibG6VhSRxqxi0SpHTu2c9ddt1NUVMSgwUPJHHe3Sj1CaMQuEoUOHTrIr+6+k0OHDnFhn/489NhT2v0SQfRfUiTKlJWVMWnSPb7z1M/gsSefIT5ev2QUSVTsIlHE4/Hw9NNP8eWXX9AirSW33/s7klNSnI4lQaZ97CJRZPacLObMeZfEpCSmPPMCpfGtnI4k9UAjdpEosXnzJv74zG8BmPTg43Tq3MXhRFJfVOwiUSA/P58HHriXw4cP07v/UAYPGeZ0JKlHKnaRCOfxeHjyycfYuXMHHc7sxMjrJjodSeqZil0kws2eO5uFCz8iJaUxk6c8R2JiktORpJ6p2EUi2K5dO/nTc78H4N4HHqNtu5MdTiQNQcUuEqFcLhe/+c0jFBcV0e38/gwanOl0JGkgKnaRCPWPf7zB119/SYu0lowZf4+mC4giOo9dJAJt2fI9L7/8PAD3PzyZlNSTHE4kDUkjdpEI43a7eeq3j1NWVsbQYaPo3ecipyNJA1Oxi0SYOXPeZe03X5HatDl3TPy103HEASp2kQiyb18uzz//HAAjr5tIatOmDicSJ6jYRSLIs88+TUFBPr1696N7zwFOxxGHqNhFIsSy5cv46KP5NGqUzL33/z+dBRPFVOwiEaCsrIxnn50CwPhb7qB1m7YOJxInBXS6ozHmWuAxIAF43lr7cqX7ewCvAYnATuB6a+3BIGcVkeOYNestdmzfRsuftWX0Vdc7HUccVu2I3RjTFpgC9AW6AROMMZ0rbfZX4HFrbVfAAvcHO6iIVG3//v288sqLAIy85hckJiY6nEicFsiumEHAImvtfmttIZAFjKm0TRxQcfg9BSgOXkQROZH/+Z+/UlCQT88L+tC56wVOx5EQEMiumAwg2289G+hZaZtfAwuMMc8DhUCvmoRIS2tSk80dl56e6nSEkKLP41gN9ZmsXvMNs2fPIi4+nnsmPUJJTDIJCd4/1qlNGtVqubaPO9Ey0CCv01DvJxivA5CSkkR6i+D/NGEgxR4LePzWYwB3xYoxJhmYBgyy1q4yxvwamA4MDTREXl4Bbren+g1DQHp6Krm5+U7HCBn6PI7VkJ/JH55+Go/Hw4hRV5PR9mTWfpdLWVk5APkFJbVaru3jTrQMNMjrNNT7CcbrABQVlZLrch3nv+7xxcbGnHBAHMiumF1AG7/11sAev/UuQLG1dpVv/TVgQM1iikhNrVr1OZ+vXEaj5MbccPMdTseREBJIsS8EBhpj0o0xKcBo4N9+938PtDfGGN/6cGB1cGOKiD+Px3PkCtOBmVdz0knNHE4koaTaYrfW7gYeBRYD3wAzfbtc5htjzrPWHgBuAmYZY9YBtwA312Nmkaj38ccfsXHjBtJapnPR4FFOx5EQE9B57NbamcDMSrdl+i1/CHwY3GgiUpWysjJefPEvANx4yy9ISkp2OJGEGl15KhJm5sydzc6dO2h/8qlkXjHC6TgSglTsImHk8OHDTJv6CgA33z6RuHj9Vo4cS8UuEkbmzs3ixx9zaNPuNPpffKnTcSREqdhFwkRpaSnTpr0GwOXDxxMbqz++UjV9M0TCxNy57/Ljjzmc3qEj55zb1+k4EsJU7CJhoKSkhGnT/heAm2/9pUbrckL6doiEgaw575Kb+yNndjyLPhdd7HQcCXEqdpEQV1ZWxj+m/x2AG26eoF9Gkmqp2EVC3L///S9y9mbTqs3J9Omn0bpUT8UuEsIOl7uZ9ve/Ad45YbRvXQKhb4lICFu0aCHbt22h1c9ac27vgU7HkTChYhcJUR6Ph+lvTgVg7LjxxMcnOJxIwoWKXSRErVmzmo3frqNxk6YMvVIzOErgVOwiIarivPV+l44kOTn4P58mkUvFLhKCNm/eyMqVy2iUnEy/gZrBUWpGxS4Sgv7uOxNm2PCxNG7S1OE0Em5U7CIhZsfOXSxc+BHx8fGMHXe903EkDKnYRULMzJnTcbvdXHLpENJbtXY6joQhFbtICCkoKOCDebMBGH3VdQ6nkXClYhcJIfPmzaGoqJAzzM/paDo7HUfClIpdJES4XC5mvvUPAPoPHu1wGglnKnaREPHpp4vZvWsnbTLa0qV7b6fjSBhTsYuEiBkz3gRg1NhriY2NcziNhDMVu0gI2LBxI2vWrCYlpTFDrhjpdBwJcyp2kRBQMVrPHDaSxo2bOJxGwp2KXcRh+/blsnDBfGJiYhg55lqn40gEULGLOGzWrLcoKyujS/cLyWjbzuk4EgFU7CIOKis7zOzZswC46FJNzSvBoWIXcdCiRQvJy9vHqaedQYdOXZ2OIxFCxS7ioHfemQnA8FFXExMT43AaiRQqdhGHbLaWr776kuSUFAZfPszpOBJBVOwiDnn7nRkADL58GCmNGzucRiKJil3EAfn5+Xz04QeAdzeMSDCp2EUc8K9/zaO4uJgOnbpy2ukdnI4jESY+kI2MMdcCjwEJwPPW2pcr3W+A14DmwF5gnLX2QJCzikQEj8fDO++8BUDfS4Y7nEYiUbUjdmNMW2AK0BfoBkww5r8TRRtjYoD3gT9Ya7sCXwMP109ckfC3evUXbNu2hZYtW3FOjz5Ox5EIFMiumEHAImvtfmttIZAFjPG7vwdQaK39t2/998DLiEiVZs3ynuJ4xfDRxMUH9D/NIjUSyLcqA8j2W88GevqtdwD2GmOmAd2BTcCvgpZQJILk5OSwePEnxMXFM/TK0ezWDkupB4EUeyzg8VuPAdyVnmMAcJG19ktjzG+BPwM3BRoiLS28ZrNLT091OkJI0edxrON9JtP+/houl4tLBl1Gm4wMfio7QEKC949hapNGJ1yu7v5gLdfHcwMR9X6C8ToAKSlJpLdIIdgCKfZdQD+/9dbAHr/1vcB/rLVf+tbfwru7JmB5eQW43Z7qNwwB6emp5ObmOx0jZOjzONbxPhOXy8Us30HTzGGjKSsrJ7+ghLKycoBql2uybV2W6+O5gYh6P8F4HYCiolJyXS5qKjY25oQD4kD2sS8EBhpj0o0xKcBo4N9+968A0o0xFRNdDAPW1DipSIRbvvwzcnL20rJVBt3P7Vn9A0Rqqdpit9buBh4FFgPfADOttauMMfONMedZa4uBkcDfjDHfApcAk+oztEg4mjPHO4vjBf0ziY3VJSRSfwI6JG+tnQnMrHRbpt/yFxx9QFVE/OTk5LB06RLi4uLp1fcyp+NIhNOwQaQBvPdeFm63m779LyH1pBZOx5EIp2IXqWcul4u5c73nEwwbPqaarUXqTsUuUs8+W76MvXuzaduuvQ6aSoNQsYvUs6x33wZg6LDROmgqDULfMpF6lJOTw4rlnxIbF8dlQzXhlzQMFbtIPZo3bzZut5uf9+hLixZpTseRKKFiF6knLpeL2XPeBaD3gKEOp5FoomIXqScrViwjZ282GRntOPOs7k7HkSiiYhepJ7Nne680vWL4GB00lQalb5tIPcjJyeGzz7xXmuqgqTQ0FbtIPZj73mxcLhcXDRiog6bS4FTsIkHmvdLUe9B06JWjHU4j0UjFLhJkS5YsIWdvNmnpbXSlqThCxS4SZDNmzAC8pzjqoKk4Qd86kSDKyclh4cKFmp5XHKViFwmiefO8B037XnSxpucVx6jYRYLkqOl5R4x1OI1EMxW7SJCsXLmc7Ow9tNP0vOIwFbtIkLyb9Q4AV44cq4Om4ih9+0SCICcnh2WfLSE2Lo7Lh45wOo5EORW7SBBUHDQ9p3sfWqS1dDqORDkVu0gd+R80vfDiKxxOI6JiF6mzioOmGW01Pa+EBhW7SB1l+Q6a6jdNJVToWyhSB/7T815+hQ6aSmhQsYvUwfvvz/FOz9v/Ek3PKyFDxS5SSy6Xi9mzvdPzDhs+xuE0Iv+lYheppZUrl7N37x7S0tvQVVeaSghRsYvUUsVB0979NT2vhBZ9G0VqYc/e/x407anpeSXEqNhFamHu3Kwj0/M2babpeSW0qNhFasjlcvH+PO+VplfooKmEIBW7SA2tXLn8yG+a9jivl9NxRI6hYhepoazZswAdNJXQpW+lSA38+GMOny1drIOmEtICKnZjzLXGmI3GmP8YYyaeYLuhxphtwYsnElrmzZujg6YS8qotdmNMW2AK0BfoBkwwxnSuYrufAc8BMcEOKRIKXC4Xc+Z4rzTVQVMJZYGM2AcBi6y1+621hUAWUNW3eirwZDDDiYSSpZ8tPTI9rw6aSiiLD2CbDCDbbz0bOOr6aWPM3cBXwOe1CZGW1qQ2D3NMenqq0xFCSrR8HllZbwMwYvQ4kpISSW3SiISE+CP/Bkht0gigytvruhzs52vI16n4TCLl/QTjdQBSUpJIb5FCsAVS7LGAx289BnBXrBhjugCjgYFAu9qEyMsrwO32VL9hCEhPTyU3N9/pGCEjWj6PnTt/YMXypcTHJzB4yJWUlZWTX1By1L8B8gtKAKq8va7LwX6+hnydis8kUt5PMF4HoKiolFyXi5qKjY054YA4kF0xu4A2fuutgT1+62N9938JzAcyjDGf1TipSAh799238Xg8dO91MSed1MzpOCInFMiIfSEw2RiTDhTiHZ1PqLjTWvsE8ASAMeZUYIm1tl/wo4o4o6SkhHnzZgPQd+Bwh9OIVK/aEbu1djfwKLAY+AaYaa1dZYyZb4w5r74DijhtwYIPOXToEKbT2Zxyeien44hUK5ARO9bamcDMSrdlVrHdduDUYAQTCRWzZr0FwPBRVzucRCQwuvJU5ATWbVjPhg3rSG3alIsH6UpTCQ8qdpETeOftGQBcnjmCRo2SHU4jEhgVu8hxHDp0kI8//hCA4aOucjiNSOBU7CLHMW/eHA6XltKpy/m0bXey03FEAqZiF6mCy+Xibd9umL4Dr3Q4jUjNqNhFqvDpp4vZs2c3GRnt6NxV88JIeFGxi1Rhxow3ARg59lpiY+McTiNSMyp2kUo2bNzImjWrSUlpzJArRjgdR6TGVOwilVSM1jOHjaRx4/CaeVQEVOwiR8nL28fCBfOJiYlh5JhrnY4jUisqdhE/7777NmVlZZzdrTcZbWs1C7WI41TsIj6HDx8+Mi9M/8GjHU4jUnsqdhGfjz6az/79eZzeoSMdOnV1Oo5IranYRQCPx8OMGdMBGHPV9cTE6DfZJXyp2EWAL1avYvPmjTRr3oKBlw5xOo5InajYRYA3Xv8bACNGXU1iUpLDaUTqRsUuUW/z5k18vnIZiYmNGDHmGqfjiNSZil2i3htvTAWg94Ch+qFqiQgqdolqu3btZMGCD4mLi2fAZWOcjiMSFCp2iWrTp7+O2+1m4OAhNE9r5XQckaBQsUvU2r8/j3nzZgMw7rqbHU4jEjwqdola/5zxD0pLS+nTdwCnnd7B6TgiQaNil6hUWFjAu7NmAjDuhlscTiMSXCp2iUqzZ88iP/8nTu/YhS7ndHM6jkhQqdgl6pSUlDB9+usADMwc53AakeBTsUvUmTNnFvv25dLhzE507nqB03FEgk7FLlGlpKSE133TB9x46y802ZdEJBW7RJWs2bPIzc3lzI5n0affAKfjiNQLFbtEjaNH63dqtC4RS8UuUWPOnFnk7cul3SkduLDvAKfjiNQbFbtEheLiIqb93Ttav3zEjRqtS0RTsUtUmDFjOnn7cul0VhfO7tbb6Tgi9UrFLhHv4MEDR6bmvf0X92i0LhFPxS4Rb+q0/6WgoIDze11Ij/N6OR1HpN7FB7KRMeZa4DEgAXjeWvtypfuHA08CMcA24GZr7YEgZxWpsb17s5n1zgwAbrvzbofTiDSMakfsxpi2wBSgL9ANmGCM6ex3f1PgFWCotbYrsA6YXC9pRWro1Vdf4vDhw3TvOYCO//3aikS0QHbFDAIWWWv3W2sLgSzA/6dmEoCJ1trdvvV1wMnBjSlSc5s3b2LevDnExcWTOVozOEr0CGRXTAaQ7beeDfSsWLHW5gFzAYwxycDDwItBzChSYx6Ph2ef/T0ej4cRo68m/WdtnY4k0mACKfZYwOO3HgO4K29kjDkJb8Gvtda+WZMQaWlNarK549LTU52OEFJC8fOYP38+a9asplmz5tx6x11kH3CTkOD9uqc2aRTU5cq3AVXeHuzXaaj3E4zlis8kUt5PMF4HICUlifQWKQRbIMW+C+jnt94a2OO/gTGmDfARsAi4r6Yh8vIKcLs91W8YAtLTU8nNzXc6RsgIxc+jtLSUyZOfBOCW2yaSnNyY/J25lJWVA5BfUBLU5cq3AVXeHuzXaaj3E4zlis8kUt5PMF4HoKiolFyXi5qKjY054YA4kGJfCEw2xqQDhcBoYELFncaYOOADYJa19nc1TigSZDNmvMnu3bto0+40Mq8c5XQckQZXbbFba3cbYx4FFgOJwFRr7SpjzHzgcaA90AOIN8ZUHFT90lp7W32FFjmeXXuy+dvUVwEYMe5O4uIDOqNXJKIE9K231s4EZla6LdO3+CW60ElCxB+f+S3FRUX06z8Q0+U8p+OIOEKFLBFj8eJPWPrpIpIaJfOr+x5yOo6IY1TsEhGKigr5wzPeQzyZo28hvVVrhxOJOEfFLhHhlVdeJGdvNh07dabfwOFOxxFxlIpdwt769WuZMWM6sbGx/PrBx4mNjXM6koijVOwS1oqLi3n0sYdwu91cdc14TCfNByOiYpew9uKLf+GHHdtpnXEKN9820ek4IiFBxS5ha/Xqz5k5czpxcXFcN+FhEpOSnI4kEhJU7BKW8vPzeeKJRwG4/sbbaX9qR4cTiYQOFbuEHY/HwxOTH2XPnt10NJ25/qbbnY4kElJU7BJ23nlnJos+WUBSoxQe/+0fiY9PcDqSSEhRsUtY2bTpW/70pz8AMO6WSbRtp990EalMxS5h46efDvHgg/dRVlbGlSOvonvPAU5HEglJKnYJC+Xl5dz/wH3s3PkDZ3bsxMS7H3A6kkjIUrFLWHjuuadZ9cUKmjRtxlNPP69TG0VOQMUuIW/WrLd4++0ZJCQkcOuvnqJ1mwynI4mENP0KgYS0pUuX8Ixv1sZJDz9B6zPOdjiRSOjTiF1C1po1q3nggXtwuVxcf+NtXDbkSqcjiYQFFbuEpE2bvuXue35BaWkpvfsP5ZYJv3I6kkjY0K4YCTnff/8dv/zl7RQWFND/ksEMv/4eYmJinI4lEjY0YpeQsnHjBm699QYOHNhPzwv68OgTT2t+dZEaUrFLyPj66zVMmHAThw4donPXXvz26edJSNB0ASI1pV0xEhI+/XQRDz00iZKSYi4eeBlDr7lP56qL1JJG7OIoj8fDm29O4957J1JSUkzmFSN4dPIfNLGXSB1oxC6OOXz4MFOmTGbevDkAZI66hfsn6UCpSF2p2MURO3Zs48GHJmE3byQpqRGP/OZ3NGvXXaUuEgTaFSMN7oMP3mPcuNHYzRtJS2/DX195g/6XDHY6lkjE0IhdGkxOTg7PPPM7Fi36GICBlw5h0Mg7MZ1OcziZSGRRsUu9c7lcZGW9wwsv/InCwkKSU1K4576HGZw5nHX/2ed0PJGIo2KXeuPxeFixYhl//etzfPedBaBvv4sZOGICA3prMi+R+qJil3rx9ddf8eqrL/LFFysBaPWz1gwdcwc3XDtKo3SReqZil6Bxu90sXbqE11//G2vXfg1AampTrr/xNkaMvoZNO37SWS8iDUDFLnWWk5PD++/PYe7cLPbs2Q1AatOmjBg1jk7nDeHCHmc4nFAkuqjYpVb2789j0aKFLF36CcuWLcPtdgPQunUGF1w8ggm3jic5JYW13+U6nFQk+qjYJSDl5eV8++16Vq5czsqVy1m/fu2RMo+PT6D/gEFcMXw0XXv0ZMOW/SSnpDicWCR6qdjlGB6Ph+zsPWzYsJ4NG9ayfv06Nm36lpKSkiPbxMcncEHvvlwy6DLOu6Af23PK6Nox3cHUIlIhoGI3xlwLPAYkAM9ba1+udH83YCrQFFgK3GmtLQ9yVgkil8vFvn37yMnJZu/evezcuYNt27aydesWtm/fSlFR0TGPOfmUUznv/N6c3/NCzul+Lt/vLua8s9tQVlYOOdrlIhIqqi12Y0xbYApwLlAKrDDGLLbWbvTb7J/Abdbaz40x04DbgVfqI7B4R9Tl5eWUl5dTXFxMcXGR759iiouLKSryrhcWFnDgwEEOHTrAwYMHOXDgAAcPHmD//jxyc3+kvPz4f/c2a9acTmedTafOXUhsejKZl/bhpGbNWftdrt/IvLhh3rCI1EggI/ZBwCJr7X4AY0wWMAZ4yrd+CpBsrf3ct/0bwJMEVuxxALGxNT8Fbgx3ROQAAATISURBVMeObcyYMZ3Dh0vxeLy3eSoWOHr92H8f2aLS46rbDuLjYykrc53w+Su2PzZH5ZzHvp7L5S1sl8vl9++ySut1+5+huLhYWrduTbNmzUlrmU5aWkvSW7UmJqkF53XrREbb9uzKc3HWaS0A2LRtP2lp3uWURgnEx8X6LcfgccdWcfvxl2uybV2WnXudmn8mof1+6r6cnBQfUe8nGK9T8ZnUpv/8HlPlz4vFVC7DyowxjwCNrbWP+dZvA3paayf41nsDz1pr+/rWOwDzrbUdA8jXF/gsgO1ERORY/YBllW8MZMQey9FD0BjAXYP7T2S1L1g24ArwMSIi0S4OaIO3Q48RSLHvwlu+FVoDeyrd3+YE959IKVX8bSMiItXacrw7ApmPfSEw0BiTboxJAUYD/66401q7AygxxvTx3XQD8GEdwoqISB1UW+zW2t3Ao8Bi4BtgprV2lTFmvjHmPN9m1wF/McZsBpoAL9RXYBERObFqD56KiEh40U/jiYhEGBW7iEiEUbGLiEQYFbuISITR7I61ZIzpDnxurU1yOovTfKe6/gVIBPKAW3ynwUad6ibMizbGmCeAq3yr/7LWPuhknlBijHkOaGmtvSnYz60Rey34zud/EW+RCczAOwlcN99yVJ7u6jdhXl+gGzDBGNPZ2VTOMcYMAgYD3fF+HucaY0Y6myo0GGMGAjfW1/Or2GvnT8DzTocIBcaYJOAxa+06303rgJMdjOSkIxPmWWsLgYoJ86JVNjDJWnvYWlsGbCJ6vxtHGGNa4B0A/L6+XkO7YmrIGHMlkGKtzTLGOB3HcdbaUrzTNmOMiQUmA+85mclBGXjLrEI20NOhLI6z1n5bsWyMORPvLpk+x39E1HgN70Wf7evrBVTsx2GMGYt3v7G/zXh/TGRQwydy3vE+E2vtIGNMIvAm3u9UvY1EQlxdJsSLWMaYs4F/AQ9Ya//jdB4n+WbH3Wmt/cQYc1N9vY6uPK0B33+UR4B8301dgbVAP2tt/nEfGOGMMU2A9/EeOL3eN4qPOsaYG/F+F27zrf8GiLHWPuVsMuf4DqzPBu611r7tdB6nGWM+xjtpYjnQAu8ULG9aa+8L5uuo2OvAGOOx1tZ8lvwIY4x5D/gR708iRu0I1XfwdBne3S+FwApggrV2laPBHGKMaQ98BVxtrV3kdJ5Q4xuxD6iPs2K0K0bqxHfa53BgI/CV77jDHmttpqPBHGCt3W2MqZgwLxGYGq2l7nM/0Aj4s9/xqFetta86Fyk6aMQuIhJhdLqjiEiEUbGLiEQYFbuISIRRsYuIRBgVu4hIhFGxi/gxxvQwxhzy+z1fjDEtjTFbjDFDncwmEiid7ihSiTHmDuBhoAfeq4wXAJ9Ya6c4GkwkQCp2kSoYY/6B93LvLcDpwGhrrf6wSFhQsYtUwRjTGPgG7w9mdLHWFjgcSSRg2scuUjUDpALNgHMdziJSIxqxi1RijGkJfIl3Js9GeKch7m6t3etoMJEAqdhF/Bhj4vAeLN1orf2V77bX8e5nv8Ra63Iyn0ggtCtG5GjPAo2BSX63TcQ7d3a0/oCIhBmN2EVEIoxG7CIiEUbFLiISYVTsIiIRRsUuIhJhVOwiIhFGxS4iEmFU7CIiEUbFLiISYf4PBr2+jy9SLbUAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "sns.distplot(data['X'], bins=uu, hist_kws={'cumulative': True, 'weights': np.full(n, 1./n)}, kde=False, kde_kws={'cumulative': True});\n", "plt.ylim(0., 1.)\n", "plt.plot(uu, norm.cdf(uu), 'k', linewidth=2)\n", "plt.title(\"n=%s\" %n)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": { "zanadu": {} }, "source": [ "### Empirical density" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "zanadu": { "code_type": "" } }, "outputs": [], "source": [ "def gaussianDensity(x): return np.exp(-0.5*x*x)/np.sqrt(2.*np.pi) ## Gaussian density" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "zanadu": { "code_type": "" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD7CAYAAABpJS8eAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3deXAc53nn8W/PgXtwD0kAvAny5SFe4iEflCybjC1LKTtaSfZGWktOSlF5HW82zm4UpyRt4lTkStnejRJbUqzVbqxYZrKObEl2TMmSrMOUKJHgfYEvb4I4CIA4CAyAObv3jx4oA/DAgJhBz/F8qliFfrsb82sCeNB45+33NSzLQgghRO5zOR1ACCHE9JCCL4QQeUIKvhBC5Akp+EIIkSek4AshRJ7wOB3gKgqBDUAHEHM4ixBCZAs3UAc0AaHxOzO14G8AtjsdQgghstTNwLvjGzO14HcA9PUNYZrZ/ZxATU0ZPT0Bp2NMG7ne3CbXm9lcLoOqqlKI19DxMrXgxwBM08r6gg/kxDVMhlxvbpPrzQpX7AqXN22FECJPSMEXQog8IQVfCCHyRFJ9+Eqpe4FHAS/whNb6yascdwfwfa31gvh2JfBjYCHQDXxBa30hFcGFEEJMzoR3+EqpBuBxYBOwBnhIKbX8CsfNBL4LGAnNfw1s11ovA/438HepCC2EEGLykunS2QK8qbXu1VoPAS8Ad1/huGeBb45ruwP7Dh/gn4HPKqW81xtWCCHE9UumS6eesWM6O4CNiQcopf4I2At8cLVztdZRpdQA4AfakwlXU1OWzGEZz+/3OR1hWuXy9VqmyYVXX6P7nd9gxUzaAU9pCbO/cDcVKy77wzcn5fLX90py6XqTKfguIHEgqgGYoxtKqRuAu4DNwOxx5xpX2DZJUk9PIFvHwH7I7/fR3T3odIxpk8vXG+nro/Mfn2X46JHL9vUfOEjVZz5LzefvxOXN3T9ic/nreyXZdr0ul3HNG+VkCn4r9mO6o2Yx9g79Huy5G3YDBUC9Umq71vpmoC1+fKtSygP4gJ5JXYEQGWDo8CE6nvkHzOEhogXFvFG7nk6X/YO1ItzOuu6D9L26jaHDh2j4o6/jra52OLEQl0um4L8B/KVSyg8MYd/NPzS6U2v9F8BfACil5gNvx4s9wDbgfuBbwBex38CNpCy9ENMg0t1Nxw+ewhwZobN2Pj8pW8/K1Qv4VEMFFwdCNJ/p4VhBHb/d+R6VrefpeOZp5vzpNzDcbqejCzHGhG/aaq3bgEeAt4D9wFat9S6l1Dal1PoJTn8M+IhS6gjwVeAPpxpYiOlkRaO0P/M05sgILdULeK7qFj5322ru+7Ri4/JZNM6u5Na1DazdvJGXl/4Og+5igidP0POLl52OLsRljAxdxHw+cEb68LNPrl1v9ws/oe/VbQwXlvFM/R18dMNC5sywu3JWL/Fzun2QwUAQgJhp0n/oKDds/2cAZv+3hyldlltv5Oba13ci2Xa9CX34C4Czl+2f7kBCZIuhI4fpe3UblmHws9qPs3hx3YfF/krcLhe3/e5m2ldswgDOPPkUAz19RJMepiBEeknBF+IKrFiMrq0/AmBHzRoi9fNZ3Vg74XnhqAm3/BYXK+spCAbY98yPCEWi6Y4rRFKk4AtxBYO7PiDS2clQcQW7/avYtLoOl2v8KOMrM9weop+5Cwvwn95PT+sVpyYXYtpJwRdiHCsWo+fffg7A274V3LFpIaVFkxtbb/pnMbhgBR7L5PBzP8HMzPfKRJ6Rgi/EOKN395cKfAwtWc3Nq+qv6/MEP7IZC5jffoQ9TcdTG1KI6yAFX4gEiXf371at4ku3LU+6K2e8aM1MRhpX4sGk/aWfEzPl3VvhLCn4QiTo/8C+u+/z+qj46MeorSpmKiODAxs+iYXB4q5mPnhfpy6oENdBCr4QcZZl0ferbQDsrFlF/Ywympo7iU7hzjxaM5PiG9fjwaRt2ytEZIymcJAUfCHigmfOEG1vY8hdhGftBooKklofaEK+zVsAaLx4nLf3tDAUijIUisr4fDHtpOALEXdp+zsAHK9qRC2oSdnntWbPJ1I1g7JYEP3r92hq7qSpuVPG54tpJwVfCMAMBrn0wfsAxFbfhNuVuh8NwzAYXmFPO7XwwlEuXgqm7HMLMRlS8IUABpp2YUTCdJbNZNbSBSn//CNqDZbLzcLhdlpOtqb88wuRDCn4QgDtr78JQNFHNuEyrm8Y5rWYxaUEFy7DhUXFif2EI7GUv4YQE5GCL/JeqL0NT/tZIi4vSz57a9peZ3iZ3a2zsv8Ep9svpe11hLgaKfgi753Z9joA0eVr8BQVp+11QnMWEfVVUhkNENTH0vY6QlyNFHyR1yzTJLxvNwALb/+t9L6Yy8XI0hsBaOjQ9AdC6X09IcaRgi/yWvfhYxSFAoRLKyhbvDjtrzeyeCUAS4ZaOKA70/56QiRK6skSpdS9wKOAF3hCa/3kuP13At8E3EAT8JDWOqyUegD4G2D0O/uXWutHUhVeiKmImnD89d8wAyhacyPD4diUplFI6jWrZxCp8lPc10377gPwsYXpfUEhEkx4h6+UagAeBzYBa4CHlFLLE/aXAt8HfktrvQIoAr4c370e+BOt9Zr4Pyn2ImMMBcMUnTwMQOesxVOeRiFZwUUrAPCdbZZuHTGtkunS2QK8qbXu1VoPAS8Ad4/ujLfN11p3KqVKgBlAX3z3BuABpdQhpdTzSqmqFOcX4rqd2X2Y8kiAcFEZkVlzpu11R+IFf8lQC3uPSbeOmD7JdOnUA4lL9nQAGxMP0FpHlFKfBZ4H2oDXEo79LrAD+Bb2XwL3JRsuvhhv1vP7fU5HmFbZcr2XmpooAcxlq/D5SgDwej34yorGHHelttF2YMy+pM4vXYBZWUNpfw+H9h7C/9kVKbqi6ZEtX99UyaXrTabgu4DEnk0DuOzvXq31K0CNUupbwNPAvVrrO0f3K6W+DZyaTLiengBmujtV0yzbVr2fqmy5XtM0KTl9BIDAvGWEA/Z0B5FIlMHA2KkPrtQ22g6M2Zf0+QuX49u7He/xQ7S0fobiwtRM1JZu2fL1TZVsu16Xy7jmjXIyXTqtQF3C9iygfXRDKVWtlPp0wv4fA6uUUhVKqa8ntBuAzBYlMkLLQU1FaIBwQQnh+vnT/vrBRTcAsGTwHPpc77S/vshPyRT8N4DNSil/vI/+LuDVhP0G8LxSam58+x7gXSAAPKyUuine/jXgxdTEFmJq2n6zA4CRhcsghROlJSsyowFXVTW+2Ahndx+a9tcX+WnC73StdRvwCPAWsB/YqrXepZTappRar7XuAR4C/k0pdQBQwJ9prWPAF4CnlVLNwDrg4XRdiBCT4T1pd+dEl6x0JoBhULJ2nZ3hyAGZI19Mi6Q6DrXWW4Gt49puT/j4JeClK5y3HbhxihmFSKnOs21UDvcS8xQQakj9zJjJci9fCW++Tl1fC2/uOY+vpIANy2biyZL+fJF95ElbkXfObt8FgGfxUnA7V1y9CxuJeQrwh/vpabngWA6RP6Tgi7wTPWZ355SvWe1oDsPjITy3EQD3aZlMTaSfFHyRV2LhMJXdZwEoXuFQ/32C0LwlANR0nc76Icgi80nBF3mlfe8hvGaUSM0s3FXVTschGC/4c4c66OsfcjiNyHVS8EVe6WraC0DZqlUOJ7GZZRWEqmdSYEWJnDrhdByR46Tgi7ziOtUMQO26tQ4n+Xfh+QqAkvNS8EV6ScEXeSPc1UVpoJeop4CSxvTPfZ+s0X78GRfPYlnSjy/SRwq+yBsdu/YAEJ23GMOTOWPdw7PmEvUUUB2+RNfpVqfjiBwmBV/kvKgJQ6EofXv3AVC6ajVDoWjaFztJmtvNcMMiALqa9jgcRuQyKfgi54UiUZoOt+FtPQ1Aa/nsaVvsJFmxhXY/fuyEjMcX6SMFX+QF74XzeM0ol4qrsHwVTse5THiO/QBWeXcLVizmcBqRq6Tgi/xw2h4BMzhrvrM5riJWXsVwSSWFsTAXm2W0jkgPKfgiL3jPn7Q/WLDE2SDXMNxgL2jesXu/w0lErpKCL3KeGQxS3teBiYFrYaPTca5ugT1UNCL9+CJNpOCLnBc6eRyXZdLnmwFFJU7HuaronEVYQEnXecxw2Ok4IgdJwRc5r/+gvaJUYJZzc98nwywuZbhyJm4rRuC4djqOyEFS8EXOCx6zp1MwFmZu//2H4hm79x50OIjIRUk9bqiUuhd4FPACT2itnxy3/07gm4AbaAIe0lqH4+vcPg/MADRwn9Y6kML8QlxTdGCAgp4LRA037vmZfYcPULlqJZG92xlpPup0FJGDJrzDV0o1AI8Dm4A1wENKqeUJ+0uB7wO/pbVeARQBX47vfgp4Smu9FNgNPJbS9EJMYOiYXTgvVtRheAscTjOxmpXLiRkuvN1txIZkumSRWsl06WwB3tRa92qth4AXgLtHd8bb5mutO5VSJdh3831KKS9wS/x4gB8C96QyvBATuRjvGhmpX+hwkuS4i4q4VFmHAQxrGa0jUiuZLp16oCNhuwPYmHiA1jqilPosdvdNG/AaUAsMaK2jCefNnky4mpqyyRyesfx+n9MRplUmXe/RE8fwAL6VKyksKxqzz+v14LvOttF2YMy+yZ4/vr2kpBDPkmWws43wiWP4P3PrhNc43TLp6zsdcul6kyn4LiBxmikDuGwSEq31K0CNUupbwNPAn447jyuddy09PYGsX/bN7/fR3T3odIxpk0nXG+npwXOpl7C7gHDNLMKB4Nj9kSiD19k22g6M2TfZ88e3Dw+HKFVLYecbdO87RE2G/F+OyqSv73TItut1uYxr3ign06XTCtQlbM8C2kc3lFLVSqlPJ+z/MbAK6AIqlFLueHtd4nlCpNvAUbv/fnjmXHBlx4A0w2VQe4MiYrhxXbzAwMU+opkzx5vIcsn8FLwBbFZK+eN99HcBrybsN4Dn4yNywO6nf1drHQG2A1+Mt98PvJKa2EJMrHOfPf6+eIlyOEnyQpEYR1oG6CydCcDRt3cRikQnOEuI5ExY8LXWbcAjwFvAfmCr1nqXUmqbUmq91roHeAj4N6XUAUABfxY//avYo3qOAjdjD+0UYlpE4xOmzVq32uEkk3epdg4ABW1nHE4icklS4/C11luBrePabk/4+CXgpSucdw64dWoRhZi8SG8vRYE+Ip4CiufNg+MXnY40KZGGBXB2J57zp52OInJIdnRsCjFJnfvs4ZiR+vkYbvcER2eeggULiRhuivq7iAWy501Dkdmk4IucNNp/X71yhcNJrk95RQkdxX4AwiePO5xG5Aop+CInGS2nAJixZqXDSa6PyzDoqbIfWwmdkIIvUkMKvsg5g53dlA33E/MUUDh3ntNxrttI3XwAgjJzpkgRKfgip0RN0Nt32xtzFzIctcjWZ/fcc+cTNVzE2luJBWTOQTF1UvBFTglFonQfOAzA4Mx5NDV3EjWz88ml2hof7YV2P/7ICbnLF1MnBV/kFNOyqOo5D0BkduZPh3wtxYUeOnz2Q+4ykZpIBSn4Iqe0nGyjOjxAzO0lUlvvdJwpG5ppvwcxIv34IgWk4Iuc0r7XHo4ZmjUXsnD8/XiuOfOJYRA6f57Y8LDTcUSWk4IvckrklD2dQjTLu3NG+WeU01FUC5ZF8NRJp+OILCcFX+SMUCRG+UW7/z5cP9/ZMClSUVbIhfhEasPSrSOmSAq+yBknT13AH+rDdLkJz5zUWjsZy2UYROfYq3VJP76YKin4Ime07j6EAYRnNIDH63SclClfthQLCJ49gxkOOx1HZDEp+CJnBONzzkQacqP/ftSC+TPoKqiCWIzg6VNOxxFZTAq+yAnDwUjO9d+PmjernLYSux9/RObVEVMgBV/kBH2qi1nBi1iGQbhu7sQnZBGvx0Wo3v6rZeS4FHxx/aTgi5zQtv8obiy8DXOwCoqcjpNy5UvtZRpHTp3AisqSh+L6JLXilVLqXuzlCb3AE1rrJ8ft/zzwTez1bc8Av6e17lNKPQD8DdAZP/SXWutHUhVeiFGR+HKGhYuXOJwkPRapOfR4y6kJDxBsaaF44UKnI4ksNOEdvlKqAXgc2ASswV6jdnnC/nLgaeAOrfVq4CDwl/Hd64E/0Vqvif+TYi9SLhozKeu2++8LGxc7nCY9lsypoLV4tB9fhmeK65NMl84W4E2tda/Wegh4Abg7Yb8X+MP4YudgF/zRTtQNwANKqUNKqeeVUlWpCi7EqPMd/dSPdAFQsCg3C35RgYdQvcyrI6YmmS6deqAjYbsD2Di6obXuAV4EUEoVA98Avpdw7HeBHcC3gO8D9yUbrqamLNlDM5rf73M6wrSa7uvd+/YeyqwY3rp6fDP9+Pq7xuz3ej34yi7v179Se7Jto+3AmH2TPT+ZYwsKvVhuF5UrVsCxNxg+cQLTMCgpKcBXUnDZ5003+X7OXskUfBeQuISEAVw2wbhSqgK78B/QWj8HoLW+M2H/t4FJDSLu6QlgZuvqFXF+v4/u7vxZhNqJ6+3ad5AyoLhxMcPDIQYDwTH7I5HoZW1Xa0+2bbQdGLNvsucnc2xgOMSB490ECkoZ8JRQPjzEjjf2sHbTKoJDocs+bzrJ93Nmc7mMa94oJ9Ol0wrUJWzPAtoTD1BK1QHbsbtzHoy3VSilvp5wmAHI8AKRct72swCULMnNN2xH+auKP+zHL4hfsxCTkUzBfwPYrJTyK6VKgLuAV0d3KqXcwC+An2it/1hrPXpLHgAeVkrdFN/+GvGuHyFSZXgkTO0lu8exOEdH6Ixyu1wM1NhzBEnBF9djwi4drXWbUuoR4C2gAHhWa71LKbUN+B/AHOBGwKOUGn0zd7fW+kGl1BeAp+N9+8eB+9NyFSJvnTt0nGIzjOWrwFNTSzgcczpSWplzFsDZ7RS0n8Wysru7U0y/pMbha623AlvHtd0e/3A3V/lLQWu9HfuXgRApFzWh88AR6oGCxsUMh2NZu2B5ssrmzWHkvQKKhwaI9VyEhrqJTxIiTp60FVkrFIkSji940l1en9ULlieruqKY9vi8OqPXLkSypOCLrGVZFrX99viBXJsw7WpchsFA7RwAQiel4IvJkYIvslbf+XbKosNEvEVEq/1Ox5k25hx7IjVZAUtMlhR8kbW69h8GYHjmXDDy51vZt2ghYcMDF7uIDgw4HUdkkfz5KRE5Jxjv0rDm5NaCJxMpKyuiq9T+i0bmxxeTIQVfZK2i+Fj0XFvhKhmD8X78IX3M4SQim0jBF1kpfKkf30g/UZeHiL/e6TjTzppnT4986agUfJE8KfgiK3XsPQTAYE0DuN0Op5l+pY2NxDCgs43YyIjTcUSWkIIvslLf4aMARGbnX3cOgLe4mD7fDAzLIijj8UWSpOCLrGSdOw2AMW+Rw0mcY86Nd+scaXY4icgWUvBF1okND1HS34lpuIjMmuN0HMdUrVgBQL8UfJEkKfgi6wwe0xjAcE09eLxOx3HM7A2rsABXRwuBwWGGQlGiuT2zhJgiKfgi64w+cOVemJvLGSbLLCyir7QWl2Vy4J3dNDV3EorIkhPi6qTgi6wTPGk/bFSzcrnDSZw3NCO+fHT8PQ0hrkUKvsgqZihEQXc7JgYVy5Y6Hcdxrvn2m9bu82ccTiKygRR8kVVGTp3EZZkMVc7AVVzsdBzHuRc0AlDe2wax3F78RUydFHyRVS4eOgKAMa/R4SSZwSwpY6C4Eo8ZxdPV5nQckeGSWvFKKXUv8CjgBZ7QWj85bv/ngW9iL1R+Bvg9rXWfUmou8DwwA9DAfVrrQArzizwTaD6GF6iW/vsPDc+cS/nZfswzJ4ENTscRGWzCO3ylVAPwOLAJWAM8pJRanrC/HHgauENrvRo4CPxlfPdTwFNa66XYSyE+ltL0Iq+YkQjujhYAGtaudDhN5rDiD5+5W6UfX1xbMl06W4A3tda9Wush4AXg7oT9XuAPtdajf08eBOYqpbzALfHjAX4I3JOS1CIvhc6ewRWLMlBag7ei3Ok4GcM13+7equhpxcrxJR7F1CTTpVMPdCRsdwAbRze01j3AiwBKqWLgG8D3gFpgQGsdTThv9mTC1dSUTebwjOX3+5yOMK3Sdb2nXz9pf7CgEb/fh9U7jK+saMwxXq8nqbbJHHut84Ex+6b79X1lRVA2i6HiCkpHLuG6eAH/jSsuOzaV5Ps5eyVT8F2AlbBtAJfdRiilKrAL/wGt9XPxriBr3GGTuv3o6QlgmuM/RXbx+310dw86HWPapPN623fuwwCKFyu6uwcZDkUZDATHHBOJJNc2mWOvdT4wZt90v/5ou8s/l9KWQ7R9sJeSOXMvOzZV5Ps5s7lcxjVvlJPp0mkF6hK2ZwHtiQcopeqA7djdOQ/Gm7uACqXU6Ny1dePPEyJZVjSKFR9rPu+mNQ6nyTyj/fiBZplXR1xdMgX/DWCzUsqvlCoB7gJeHd0ZL+i/AH6itf5jrbUFoLWOYP8S+GL80PuBV1IZXuSHqAl9x0/gikXpL6nG8JUzFIqS5X/8pZQRfwCroP2s9OOLq5qwS0dr3aaUegR4CygAntVa71JKbQP+BzAHuBHwKKVG38zdrbV+EPgq8JxS6lGgBfjddFyEyG2hSJQT23dTBfRVz6apuROA1Uv8zgbLIGZFNYGicsqCAwTPt1A8b77TkUQGSmocvtZ6K7B1XNvt8Q93c5W/FLTW54Bbp5BPCABcLacAMOfm7/z3ExmcMZeylsN07TvEPCn44grkSVuR8axolNLuVgAKFi9xOE3msuK/DC8dOeJwEpGppOCLjBduOYvHjNJXVInbJ+Pvr8YVn1fH3XpG+vHFFUnBFxlvKD7yZMCfvuGGucCsqGa4uBxvJETwfIvTcUQGkoIvMt7AEXvB8pj0308oOsf+P7p44LDDSUQmkoIvMpoVjeJqPQtAYaP030+kfIU9zVX/IenHF5eTgi8yWvDMadyxCP1FVRjSfz8hf3xSOXfraSyZH1+MIwVfZLSe/QcBGJg5z+Ek2aGgppbhkgo8kRDBc+ecjiMyjBR8kdH6D9l90daC/F6wfDJi8+z/q9FflkKMkoIvMpYZCuHpaMHEwCv990mrWmV364z+shRilBR8kbGGjmtclslgxUwoKnE6TlYwXAbVq+P9+G1nCQRGiMqQfBEnBV9krK49BwAwFsndfbJCkRj720foLanBbUY5+PZuQpHoxCeKvCAFX2SsoWZ7/P2MtascTpJ9BmbOB8B19oSzQURGkYIvMlIsEKCwp4OY4aLqhmVOx8k61jx7mgVPy0mHk4hMIgVfZKRLR45gAMP+ObgKCp2Ok3UKGpcQw6CsrwNzZMTpOCJDSMEXGakz3n9fvEzu7q+HUVRET9kMXJZF6KR2Oo7IEFLwRUaKnDgGwOyNax1Okr0G4/34AZkuWcRJwRcZJ9LbQ9FgLxF3Ab7GRqfjZC0z/rDa8NGjDicRmSKpFa+UUvcCjwJe4Amt9ZNXOe6fgDe11j+Mbz8A/A3QGT/kl1rrR6YaWuSuqAltO/cAEGxYwHDUkrVrr1PRokbCv/ZQ0NNJpLcXb3W105GEwya8w1dKNQCPA5uANcBDSqnl446pV0r9Arh73OnrgT/RWq+J/5NiL64pFInS9r5d8AfrFtLU3ElUFvO4Lm6vl67KBgCGjhxyOI3IBMl06WzBvmvv1VoPAS9weWG/D3gZ+Mm49g3AA0qpQ0qp55VSVVNOLHKaFYtR0XkGANdiecN2qoZn211ivXv3O5xEZIJkunTqgY6E7Q5gY+IBWuvvACilNo07twP4LrAD+BbwfexfDkmpqSlL9tCM5vf7nI4wraZyvUM7T1EQCzNYXEnpbPvu1Ov14CsruuzYK7Un25aq84Ex+6b79Sc61rV2LRz6NeHjx6itLsFwuy/7PJMl38/ZK5mC7wISe1ENIKm/sbXWd45+rJT6NnBqMuF6egKYWd6B6/f76O4edDrGtJnq9Z7fvhMXdncOgSAAkUiUwfjHia7Unmxbqs4Hxuyb7tef8NjyKgYKyykPDXC+6SDFi6b2Jrh8P2c2l8u45o1yMl06rUBdwvYsoH2ik5RSFUqpryc0GYBM6iGuaaTZHkIo3TmpMxLv1hk4KNMl57tkCv4bwGallF8pVQLcBbyaxHkB4GGl1E3x7a8BL15fTJEPYoEAJRfbiBkurHmyfm2q+FbZcxH17jvgcBLhtAkLvta6DXgEeAvYD2zVWu9SSm1TSq2/xnkx4AvA00qpZmAd8HBqYotc1L13PwbQX9WAJdMppMzcj9xIDANXewuxoSGn4wgHJTUOX2u9Fdg6ru32Kxz35XHb24Ebp5BP5JELTfsoASILZDrkVCosK+FSdQPVva0MHT1C+YaNE58kcpI8aSsygmVZuE7bc74YjUsdTpN7ipbdAMCFXXsdTiKcJAVfZIShc+coCgUIF5YSq62b+AQxKQtutu/qw82HsazsHvkmrp8UfJERWrZ/AIC1eDkYhsNpck/5ogUEC8soCAYYOXfW6TjCIVLwRUYYOWiPIJnxkQ0OJ8k9hstgOBwjttieEeXs2+8zFIrKWrd5SAq+cFyorw9fXwcxl4eS5SucjpNzQpEYTc2dBOYoAAIH9tHU3Clr3eYhKfjCcWfeiXfnzG/EVSjDMdMlNm8RUZeHysFujMFLTscRDpCCLxzXv9ceOSLdOWnm8XJp5gIAos0ye2Y+koIvHBUNh/FdsGfHrFgrj2yk3RK7y8x7qtnhIMIJUvCFo06/uxuvGSU6owFvlcyenW6RhfYzDrW954mFLp+QTeQ2KfjCUV07mwCoXi9399PBLPUxWFWHx4pxfscep+OIaSYFXzgmGIlR0mI/XVu4cjVDoagsZzgNYo328Mze3VLw840UfOGY5h0HKIsMEy4q40CgSJYznCbhRXbB97UeJxaLOZxGTCcp+MIxF9+3h2OGGlfI07XTKFozk5HSSkqiQU7tkLl18okUfOGIcCRK5bmjAIQWr3Q4TZ4xDMJL7P/zznffdziMmE5S8IUjjrx/kIpIgHBRKeG6eU7HyTvh+C/ZsnPNxKLSrZMvpOALR3S9Z3fnhBtvAJd8G063iL+eaHkVpdERTuzY53QcMU3kJ01Mu6GRCJUt9oM/wcYbHE6TpwwD33r7yeauHdKtkwUvt2oAABKVSURBVC+SWvFKKXUv8CjgBZ7QWj95leP+CXhTa/3D+PZc4HlgBqCB+7TWgRTkFlns8I6DVEUGsUp9hOvnOx0nb/k2bGTkzdfwnW3m0lAIj8dNodeDR24Dc9aEX1qlVAPwOLAJWAM8pJRaPu6YeqXUL4C7x53+FPCU1nopsBt4LCWpRVa7+IHdnVN64zrpznGQWT+HYEkFpdFh3vn5uzKDZh5I5qdtC/Zde6/Wegh4gcsL+33Ay8BPRhuUUl7glvjxAD8E7plqYJHd+gMhalqPAVCydp3DafKbYRiEF9tdau5jBx1OI6ZDMl069UBHwnYHMGYVZK31dwCUUpsSmmuBAa11NOG82ZMJV1NTNpnDM5bf73M6wrS61vU2vbOf6sggLp+PylUr8Z3sGbPf6/XgKyuasG0yx6bzfGDMvmzL7169Dg68R13ncaJFXkpKCvFXl1x2bCL5fs5eyRR8F5D4wLsBJPM45PjzSPK8D/X0BDCz/Fl7v99Hd/eg0zGmzUTX2/nmW1QA5Rs2MhKKMhgYO4FXJJJc22SOTef5wJh92Zb/km8G7rJKygL96Ka9DC/YQvc1nr6V7+fM5nIZ17xRTqZLpxVIXFV6FtCexHldQIVSyh3frkvyPJGjLnQP0nDBnjun/GObJjhaTAvDILTMnriuqFmeus11yRT8N4DNSim/UqoEuAt4daKTtNYRYDvwxXjT/cAr1xtUZL+Dr71LaSyIe2YdhfPmOx1HxAWXrgWgvvsUocCww2lEOk1Y8LXWbcAjwFvAfmCr1nqXUmqbUmr9BKd/FXtUz1HgZuyhnSIPmaZFbN8uAKo+/nEMmTsnY8Qqqhn0z8FrxTj75nan44g0SmocvtZ6K7B1XNvtVzjuy+O2zwG3Xn88kQuiJuw/fI55l85hYeC+caNMhZxhIivWwdvnCe/+AO65w+k4Ik1kELRIu1Akit72Nh7LJNSwgH0XIjIVcoYJLV5JzOWmuuc8g+2dTscRaSIFX6TdcDBCXesRAEaWycpWmcgqLGJgjgLgxLY3HE4j0kUKvki7ve8fYXawm5jHS3Dh8olPEI6IrbTfkrMONGFZ0t+Wi6Tgi7QyLYuh7e8AdreBVVDocCJxNeG5jURKfPhG+mlvkhk0c5EUfJFWh3UHjd322PuhlR9xOI24Jpebko/fAkDHq685HEakgxR8kVb6lbcoMsOEZswmMqPB6ThiAjWf/BSm4aKsRRPqueh0HJFiUvBF2nT0DDHz5B4AhlfJ3X02cFdUEF1yAy4sTr084fOVIstIwRdp8/5rO5kV6oWSUkZkoZOsYLgMZt72GQCiu3cQGAoSldGzOUMKvkiLgeEwRtN7AJR9bBN4vA4nEskIRWIciVYwWFZDUXiY3S++LnPk5xAp+CIt3n5XowbPgGFQevOtTscRk2EYjMS74Ar37XA4jEglKfgi5YLhKANv/RqPZVK6chWeWr/TkcQkxW5YR9hdQM2lDoZOHHc6jkgRKfgipQaHw7z1/ilW9RwFoHjzbTJnThayCgrpXWo/iNX+s5ccTiNSRQq+SKnBoTAdr71uD8Wsn8/BSLnMmZOlrI23EDXclLQcZ6Slxek4IgWk4IuUen9fC2suHgYgsO4TDqcRU2GV+uhetAaAUz/5mcNpRCpIwRcpEzNNTv/bryiLBQnX1hGau9jpSGKqPvZJTMOF59gBQp0XnE4jpkgKvkiZpiMdLG/bD0Bg/SdAFjnJelZ5FcNL1uDC4swL0pef7aTgi5QwLQv989eojAaIVNYSXLjC6UgiRebe9TtYAPt3MdLe4XQcMQVJrXillLoXe3lCL/CE1vrJcfvXAM8C5cBvgK9oraNKqQeAvwFGV1T4pdb6kVSFF5lj/5E2VrbsBGDwps3gknuJXFHcUM/I0rWUHNvH/qf+L43/9Y+cjiSu04Q/lUqpBuBxYBOwBnuN2vGTmj8PfE1rvQQwgD+It68H/kRrvSb+T4p9DrIsi7M/fZmyWJDChYsINq50OpJIsaVfvo+o4cY8tJfhkyedjiOuUzK3YVuAN7XWvVrrIeAF4O7RnUqpeUCx1vqDeNMPgXviH28AHlBKHVJKPa+UqkpddJEpjhw6y5JWe/50/xd/V/ruc4zhMoj5Khm5cRMAp5/7EZGYPFyRjZLp0qkHEjvuOoCNE+yfnfDxd4EdwLeA7wP3JRuupqYs2UMzmt/vczpCWnX89KfMs6JUbdxA8RKFT3eN2e/1evCVFV123pXak23LlPOBMfuyLX8yx5oY6PP9mDdvYeTAToo7zhHYv5clt9162efMRbn085tMwXcBib/ODcBMZr/W+s7RRqXUt4FTkwnX0xPAzPLHNP1+H93dg07HSJvmnYeZ03YEy3BR8bn/AMBgIDjmmEgkelnb1dqTbcuU82Hs9WZb/smd72J43Scp3vkKrT/+MRWrVuLyFlz2eXNJtv38ulzGNW+Uk+nSaQXqErZnAe0T7VdKVSilvp7QbgAy7V4OCYWj9G99DhcWJTffSqTKTygSczqWSKOyWz7JpeJKCgd6af/Zi07HEZOUTMF/A9islPIrpUqAu4APV0bQWp8Dgkqpj8ebvgS8AgSAh5VSN8XbvwbId0gOOfD8T6kZ6iZY7OP00o/T1NxJNCbTKOQ0t5tLn7oTCwi88StC52XKhWwyYcHXWrcBjwBvAfuBrVrrXUqpbUqp9fHD7gP+Vil1DCgD/l5rHQO+ADytlGoG1gEPp+MixPQLtHdQ+r697mngU7+DVXB5H7HITZ4FjXQvXo/LMjn7zDNYMfmrLlskNQ5fa70V2Dqu7faEjw8w9o3c0fbtwI1TzCgyjGVZHH/qB5RZMXrn30B4wVKnI4lptuzBL3HusWOUdbRy8Vev4r/9DqcjiSTI0zFi0k6/+AvKLpwlXFBMePPnnI4jHFBUVkrxXfcC0PPyiwRbzjmcSCRDCr6YlEF9nMgr9pwqVb/7JcziUocTCacs2/wxLsxfjSsWpeX73yM2POR0JDEBKfgiadGBAc59/3u4LJPoxluo3HhZL57IE4bLYCgUZelXHuRicQ30XqTlmWfkgawMJwVfJMUyTY797d9TMDLIQM1slv3+/U5HEg4KRWI0NXdytHWAi5/5jwRdXiKHD9Dzyi+djiauQQq+mJBlWRz/wbMUnD9J0FvMDX/6dQxPUu/3izxQPrsevc5+0/bSyz9lcE+Tw4nE1UjBFxM68aN/xtizg6jhpvYP/jMxXwVDoaisVSs+VLtxPUfn34RhWbT94B8Ybj7qdCRxBVLwxTUd/9eX4DevYWLQufkejrtqaGrutB+ykrVqRZxhGFTc/tu0LliLy4xx7u+eYPj0pGZSEdNACr64qiP//CL8yh6R0/WJz2MslWmPxdUZLhfr//grdM9Zjjsa5vR3vsPFw81OxxIJpOCLy1imyZ7v/W+8v34ZgILP3Y25coPDqUQ28BR4WPXf/wvD85dSEAnS9Xffpfn17U7HEnFS8MUYI8Mh3v+r7+A78B6m4aLsvt+j5tO3OR1LZIlQJMaeE730334vvY1r8VgxjP/3f3nvH/+VmHQBOk4KvvjQxRNnOPSNP6e2tZmIu4Ce376f41WN0lcvJs/lJviZu+hb90lcWPjf+yW/+Ytvc6mn3+lkeU0KvsA0Tfb8+EW6vv3XVA73MlJWRd89DxGZt9jpaCKbGQYjH91C1f2/j+nx0tBxjBOPPcahd/YQlXsIR8hg6jx3du9ROrY+T02/vcSBsfYm+jd8Bqug0OFkIld41n+E7lg55a/8C+W9F7B+9D2aPriR9V95AG9FhdPx8ooU/Dx1/uR5Tm39V+paDlIDRItKqf2P91G8fiNtxzqdjidyTKzKT98X/zPh91/Hd2AHNSf2oh8+TMmn72D+b9+Gq1BuMKaDFPw8c3z3Edpf/gWzOjT1WJiGi4GVH2H4pk/RVVjMaumvF+ni9jC06bMEl91IxXuvUNRygugrL3Ls169S/olPMuuzt+EpL3c6ZU6Tgp8HRnr60K++RXD3TqoHO6kHLMPAu2Y9bcs3Ea3yOx1R5JFYzUwWfePP6N+zl46fvUh5XwfDr2/j5K9/hVutYNatt1C6alXOr5frBCn4OciKxehqPknbzt1E9VHKe9spwaIEML0FlHx0ExVbPo2rupZz0n0jHBCOmpwoboB7/5D+Myfx7HyHWT1nsZoP0tF8kJi3ABYsoWrtaqrXrMZb68cwDKdjZ72kCr5S6l7gUcALPKG1fnLc/jXAs0A58BvgK1rrqFJqLvA8MAPQwH1a60AK8+e92MgIg+fO033qHIGz54idP0dxbwdeM8ro2vWm4SK6cCmXFq0kvGAZlrcAumKsrpTuG+Eww8CzcDEsXMz5vj5KTx3CfXgvVYFuOH6YgeOHGfh/PyZUUEywtgH3nHlUzJ/LDLWQsvo6mcRvkib831JKNQCPY69JGwJ2KKXe0lonzo70PPCg1voDpdT/Af4AeBp4CnhKa/0vSqnHgMeAP0v1ReQKKxbDioQxgyHMUBBzJEhseAhzeIjwQIDAxV5GevsI9/Vj9vXgHrxEYWTkw/MTlyIJ+6ooUMuoWruaIrUMikrokrt5kcG8VVWseeALWOY9jHR20b5zL2F9hOL2sxSGRyhsPwntJ4nthA7AxGCksJRwSQWmrxJ3eTmFVZUUV1dSXOGjtKqc4spyPEXFGIWFuAoLMbxeDFf+jkZP5tfjFuBNrXUvgFLqBeBu4K/i2/OAYq31B/Hjfwh8Uyn1LHAL8DsJ7e+QXMF3A7hck/8TzjQt9p3oZnAkMtpA1f7teAb7PjzGuuyDcawr7LhSG1ZCe/zj0W3LPsdlWFgxEzAxTAvDsjBME8wYhmVimCaGGcVlxnBZE99xF8f/4QGqyogZFUTLqzBqaimeVUfpgnmcipRiFpUwDPQDXAiybEEJJUXeMZ/L43Zdd9vVjzWmeP5UX396zy8u9BCLeq95bCbnn+z5o9ebrkwx06L5TC/ggZUbYeVG5iyo4sTBM7gvtkPnBYz+i5QE+ikavmT/LBCGYJf9r8v+PCYwGP932WsYbiyXG9PlBrcby3CBy4XlcoPLAAy7zTA47nFhmvZ7Xox2KRkGBgYW9qRxGAZul0F5acGYbqer9kAl0TVVumoNZTeum/C48RJqpvuKL21dsZD9O6XUnwOlWutH49sPAhu11g/Ftz8KfEdrvSm+3QhsAz4BNGmtZ8fbPcCw1jqZd2I2ATIBhxBCXJ+bgXfHNyZzh+9i7L2wgf0LdKL949sZd961NGEH7gBiSZ4jhBD5zg3UYdfQyyRT8Fuxi++oWUD7uP11V9jfBVQopdxa61j8mMTzriXEFX47CSGEmNBVFyJI5t2LN4DNSim/UqoEuAt4dXSn1vocEFRKfTze9CXgFa11BLtb5ovx9vuBV64jvBBCiBSYsOBrrduAR4C3gP3AVq31LqXUNqXU+vhh9wF/q5Q6BpQBfx9v/yrwkFLqKPZfCY+m+gKEEEIkZ8I3bYUQQuSG/B2QKoQQeUYKvhBC5Akp+EIIkSek4AshRJ6QmYemiVJqLfCB1jqnV3qID8/9W6AA6AF+Pz50N6dMNKFgrlFK/QXwhfjmL7XWDzuZZ7oopb4L1Gqtv+x0llSQO/xpEH9+4XvYRTDX/Rh7Ir018Y//foLjs07ChIKbgDXYQ4+XO5sqfZRSW4BPA2uxr3edUupOZ1Oln1JqM/CA0zlSSQr+9PifwBNOh0g3pVQh8KjW+mC86SAw18FI6fLhhIJa6yFgdELBXNUB/DetdTj+QGUzufl1/ZBSqhr7l/q3nM6SStKlk2ZKqc8BJVrrF5RSTsdJK611CHuqbJRSLuAvgZeczJQm9dhFcFQHsNGhLGmntT4y+rFSajF2187Hr35GTvgB9gOnc5wOkkpS8FNEKXUPdt91omPYi8Jsmf5E6XW169Vab1FKFQDPYX9/5dQdUtxEEwrmJKXUCuCXwJ9qrU84nSdd4jMCn9da/1op9WWn86SSPGmbRvFvnD/n36flXg0cAG7WWl9pqu6sp5QqA36O/Ybtf4rf9ecUpdQD2F/DB+PbjwGG1vqvnE2WPvE3438K/LHW+l+czpNOSqnXsSd7jALV2NPFPKe1/rqjwVJACv40UkpZWuucXphTKfUS9kypX9Fa5+Rdb/xN23exu3GGgB3AQ1rrXY4GSxOl1BxgL/BFrfWbTueZTvE7/FtzZZSOdOmIlIkPPf08cBTYG3/Pol1rfbujwVJMa92mlBqdULAAeDZXi33cfweKgP+V8D7UP2it/8G5SOJ6yB2+EELkCRmWKYQQeUIKvhBC5Akp+EIIkSek4AshRJ6Qgi+EEHlCCr4QQuQJKfhCCJEnpOALIUSe+P9Ldh4YdgzDfQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "xx = np.linspace(-5., 5., 100)\n", "sns.distplot(data['X'].values, kde=True, kde_kws={'cumulative': False})\n", "plt.plot(xx, gaussianDensity(xx), 'r', linewidth=2)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": { "zanadu": {} }, "source": [ "# Convergence of the mean" ] }, { "cell_type": "code", "execution_count": 16, "metadata": { "zanadu": { "code_type": "" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA7EAAAGECAYAAADgC6duAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeZgcV3ku8LeX6ZnRaLEtC2zjgE2AA4YQEwjLNQ7cQHITJ0AWssElyTXYOOw4bMFA2GO84xUvyMa7LduyLe+LJGvf9+1oG2k0o5E0Gs2+9FZ1/6g61VWnq6qru2emZ9Tv73n8WNNLdXV1d536zvedc2KmaYKIiIiIiIhoKojXegeIiIiIiIiIomIQS0RERERERFMGg1giIiIiIiKaMhjEEhERERER0ZTBIJaIiIiIiIimDAaxRERERERENGUka70DVB+EEAkAXwfwGVjfuxSABQB+JKVM13LfpgohxJ8DuBPAUQAfkVKOuO67E8BvpJTrhRCLAdwspXysNnsaTgjxUwB7pZT3+ty3CcBHpZS9FWx3MSbx+yYimmhse6vHtrfkdhdjHN63EOIsAI9JKf/XWG6XTh4MYmmi3AbgVAAfk1L2CSFaADwA4C4An6vpnk0d/wzgTinlz33u+zMAt0/w/lRESvmjkPvOn8h9ISI6ybHtrR7b3hqQUh4GwACWAjGIpXEnhDgHwGcBnCml7AcAKeWQEOIyABfYj5kF4BYA5wMwATwP4PtSypwQYhTAlQD+HMCZAK6SUt4mhFgB4Fop5eP2Nn5lb/u7QojPA/gSrJL5bgBfkVLuEkLcA+A0AL8P4BkA1wC42/67G8ARANuklD8WQrwDwK8BzAaQAHCjlHKuEOKjAH4BYD+AdwFoAPBFKeVyIcR0ADfZ7ysH4EkAV9iP+RWAj9jb2gjga+p4uI5VA4DrAHwMQB7AagDfBHAZgL8BMCKEmCWl/LbrOb8AcBaAB4QQ/2rf/CkhxLcBnAHgFQCXSCkNIcT/svejxd7+T6SUz/h8ZmHv/X8AtAEQAIbsz+Zr9t+PSym/aT/uVwAOAng7gBEA/y6l3Gl/BtuklNcIIdIAngLwh7C+I2sBzJFSHhdC/BeAf7OP4x4A/27/+zYAb7X3bQDAZ6SUUn8PrveyGMB6AB8E8DoAd9jH5SP2cfhHKeVW+zv4awB/YH9erwL4tv0dvBjAF2FlMU4DcKX9Hfx3AH8LwLD3aRjAv0kpdwbtDxHRRGDby7YXNWp7hRBvA7ACwFlSyoxdEdAG4OMAZgG4CkAjrO/Vy1LKz9vf16UAdgI4x96Hl6WU04UQr4fVWfB6+9gehNV2HxNCHABwj/3ZvRHAvVLKH9r7cTGA/7SP+XFY7fMhIcQnAPwAVps+DOBbUsqVfu+FJi+OiaWJ8F4A2/VGQ0p5RDWCAG6E1ZD9AYD3wTqxfsu+rxHAcbuk5NMArhdCNMEq7/l/gFMy9X8B3CWE+Aisk9+FUsr3wDpZzne99DQp5TullN+1X3e7lPIdAP4Bdq+fECIJ4DEA35NSvhdWA/gtIcQH7W18AFYj/h5YDfEv7dt/CqAJwDtgXRRcYD/3e7AagfdKKf8QwGFYDZDuB7AaxT+0/4sDuFpKeTWApwFc725E7eN4hb29z0opV9s3z7DfyzsA/CWAC4QQp9r7+jkp5R8B+BSA24QQb3RvL8J7/2NYQdz5APoB/BeAvwLwRwC+bJcAAdbneJOU8t32697n835TABZIKYWUcp1rHz4Jq+H8kJTyXQBaAXzFfi+9UsoPSSnfBqvh/YrPdnXnSCkvgPUduQrAYinl+wC8AOCr9mOuB7Defs/vAXA6gMvti6NLAFxkf97/ZG9D+QiAr9r7uRrWZ01EVGtse9n21qTtlVLuBrAdwCftm/4cQKvdwft1WOXsHwBwHoBPCiHeaz/ubAA/s1+j07XJfwawUkr5IQBvhhV4uisJpkspL4R17L8lhDhXCPGHsAL6v7CPxdMArhBCvBXW90a16ZcCeMKuUqAphJlYmggGSneY/CWAC6SUJoC0EOI3AL6BQmPzlP3/DbAa1hYAjwC4RghxBqyT+G4p5R4hxCUA3gJghRBCbf9UIcRp9r+XuV73Ivu5kFJ2CiHUmI63weohnuvaRjOs4GYngINSyk2uffp3+98fB3C5lDIPq+fvIwAghLgKwCkA/szeXgrAsYDjcIWUMms/7yZYPcrlesTeh2EhxB5YGcgPwer1fNL1nkwA74bVQ6qUeu+tUsqN9u37APRJKTMAjgsh+mH1tgPAZinlUvvfcwHcIoSY7bOvS31u+ziAeVLKHgCQUl6u7hBC7BdCfBXWZ/xRAFF6T59w7S9gBa/q74/a//5rAO+3MwmA9Z4hpRwUQvw1gL+yG7/zAUx3bXu9lLLd/vcGAH8XYX+IiMYb2162vbVse++C9fk8BqvT40779n8DcJEQ4vuwssXNsNrUblgdDkXblVL+WghxoRDicljZYNVprDxlP65DCHEM1rH4CIAXpZSH7PtusN/Hl2B9Hq+6jrNhv6/NJd4TTSIMYmkirAbwDiHEDCnlgLpRCPEGWKWdn4bV0Jqu58RhlQEpIwAgpTTtk05MSjkshJgHa8KKD8E6YQJWCc59dm8vhBBxWD2sPfb9g67t5gDEXH/nXdvok65xInY5Sx+sstQR13NM1zZy7vchhPg9WD2GCQBfl1I+b98+HVavsS5R4jhElfXZvwSAnXbvp9q/swB0+exD2HvXJwPJwl/O9W91fPI+jxv0uU0/jqfAuhD5S1i9pjcDeBDACQDnBry+m2ef1YWKJgHgH+yeYvWaphDibFiN6h2wLsIegxXwKkHfBSKiWmLby7a3lm3vPADX2SXSH0Ghw2EJgC2wOpMfhZVdV/uZllLmtO2okvX3wwrKF8H6bNzfH7/vhf5emgG8CdZxflVK+U+u+34PVladphCWE9O4k9bg/Adg9S7OBAD7/7cC6JbWTH8vAviKECImhGiEdbJ8OcLm74TVq3cBAFUe9SKAfxFCnGn/fRms8Y1+ngXweXufZsMa32gCkLDGwPxf+77fA7ANVnlWmFcA/JsQIm6/j8dg9wba7y9lN+x3whrfonsBwH8IIRrsx30Z0Y5DDqUb3FUA3iqE+BP7PZ0Pa7zLG7THVfredecLId5t//tSACtk9JkPXwHwd+r7AuDHAC4H8H8A3COl/K29n5+A1SCNhRcBfNP1HXwaVrnU+2BdbPwcwEuwA1i7jI6IaFJi28u2FzVse6WUowAehjVe9XG78+MUWGXR35VSPgGrfPgtpbZlv/4NUsr7YGXS/yzCcxYB+Ljr+/hFWCXurwL4cyHE2wFACHERrKC6ucT2aJJhEEsT5UsAdsAqM9oEq4d4B4Av2Pd/DVbZzVb7PwlrAodQUsr1sHoYH7NPmJBSvgRrHMTLQogtsHqL/84ul9J9E8DbhRBbYTXEBwEM2yU6nwLwBXsbLwH4oZRyeYld+gmADKySlI0AnrNP1D8DcMC+bQesXsL/9Hn+z2FNcLEJVvlQA6zxI6U8AeB+YS0F4EtK2QXg7wFcLYTYDGuczOeklAe0x1X63nVHAPzCPrZ/gzJmwpRSPgdrLM9y+/lnwJqk4xoAX7T3aymscrK3lLlfQb4Gq1RuK6wGbSusBu8lAO2wvpM7YU0c0TWGr0tENF7Y9rLtrWXbeyesDOpd9vZ7YXUibBBCbIM1Znl5hG39FFYJ+xZYHczLSj1HSrkVwLcBvGAf978AcJmUcges4P5h+/afAfiklNIvM02TWMw0/c4tRPXBHhuxUUq50u69XQrgv1XpEVVGWDMk3mxPDEFERORg2zs+2PZSPeGYWKp3OwDcZJeGpmBNaMBGlIiIaPyw7SWiqjATS0RERERERFNGVZlYIcRnYK2t1QBrwPUt2v3nw6qDnwlrNrLLpLWA9pn27WfBmj3us/rYACIiIiIiIiJdxRM72VO0/wLAh2Gtm3ipEOI87WH3A/iKvWhxDMAl9u33wVpk+T32v39V6X4QERERERFR/agmE/txAAullCcAQFgLVX8a1gxiEEK8CUCzlHKV/fh7APzEXlvsD2FNjw1Ys6AFTcHupxHW9Nyd8F/3ioiIqBwJAGcCWIvitRgpGrbNREQ0lkLb5mqC2LNgNVZKJ6xptMPuPxvA7wNoA3CtEOJCWFOBf6WM1/1jWLPYERERjaULYS3dQOVj20xEROPBt22uJoiNw1qYWokBMCLcnwTwHlhTqV8uhPgCgN8B+GjE1+0s/RAiIqKysX2pXCcA9PQMwTAqnzBy9uzp6O7mco2l8DhFx2MVDY9TNDxO0YzFcYrHYzj11BYgoG2uJohthxUZK2cAOKzdf6bP/UcADEgpn7FvfxDAjWW8bh4AursHq2oo58yZga6ugYqfXy94nKLjsYqGxykaHqdoxuI4xeMxzJ49HWAZbDXyAGAYZlVts9oGlcbjFB2PVTQ8TtHwOEUzhsfJt22ueGInAK8A+JgQYo4QYhqAvwfwgrpTSnkQwKgQ4gL7ps8BeF5KuQ9AuxDiL+3bPwFgfRX7QURERERERHWi4iBWStkB4AoAiwBsAvCglHKNEOI5IcT77Id9FsD1QohdAKajkHH9OwDfFUJsA/B1ABdXuh9ERERERERUP6paJ1ZK+SCscmD3bRe5/r0Z3sme1O0S0cfAEhEREREREQGorpyYiIiIiIiIaEIxiCUiIiIiIqIpg0EsERERERERTRkMYomIiIiIiGjKYBBLREREREREUwaDWCIiIiIiIpoyGMQSERERERHRlMEgloiIiIiIiKYMBrG2TDaPW57YiuO9I7XeFSIiorq0elsnFqw4UOvdICKiSY5BrG3Lvm6s392FRxburfWuEBER1aWV2zrx2qaOWu8GERFNcgxiiYiIaFKIIVbrXSAioimAQazGrPUOEBER1TGTDTEREZXAINYWY+cvERFRTbEtJiKiKBjEEhER0aRhMhVLREQlMIjVsPEkIiKqjVgsxmE9RERUEoNYB2uYiIiIaikWAyenICKikhjEEhER0aTBGJaIiEphEEtERESTQiwW4/TERERUEoNYIiIimhRYTUxERFEwiCUiIqLJgYlYIiKKgEEsERERTQqcYpGIiKJgEGvjAutERES1FYvFuNQdERGVxCBWw7aTiIioNtifTEREUTCItbHhJCIiqjGOiSUioggYxBIREdGkEIvFODsxERGVxCCWiIiIJgWrKophLBERhWMQS0RERJMDy4mJiCgCBrFEREQ0KcTAcmIiIiqNQSwRERFNCrEYWE1MREQlMYhVOD0xERFRzZmMYomIqAQGsURERDQpxGIxZmKJiKgkBrEakzNKEBER1QSriYmIKAoGsbYY64mJiIhqKsbZiYmIKAIGsURERDSJMIolIqJwDGKJiIhoUojFYszEEhFRSQxiiYiIaFKIcWQPERFFwCBWww5gIiKi2mEmloiISmEQq7D3l4iIqKZisRjXiSUiopIYxBIREdGkEANYEkVERCUxiCUiIqLJIcYYloiISqvrIPZozzD6hzIAWE1MRERUa1yznYiIoqjrIPa/bl+Fb926HAB7fomIiGpNzU5scnYnIiIKUddBLADk8mwoiYiIJgOVh2XLTEREYZK13oHJggVMRER0shBCzASwAsBfSykPBDzmXgALpZT32H+/EcD9AF4HQAL4rJRycEJ2WHFSsWDDTEREgeo+E6tjBRMREU1lQogPAFgG4G0B958lhFgA4NPaXbcCuFVK+XYA6wD8cFx31EchhmVjTEREwRjE2mLs8SUiopPDJQC+DOBwwP2fBfAUgEfVDUKIBgB/AuAx+6Z7APzD+O2iP6ecmDEsERGFYDkxERHRSURK+QUAEEIE3X+1ff+HXTefDqBfSpmz/+4EcPY47qY/digTEVEEDGI1LGEiIqI6FEfxfEpGuRuZPXt6VTuhltg5/fTpaEgmqtrWyW7OnBm13oUpg8cqGh6naHicohnv48Qg1sHuXyIiqlvHAMwSQiSklHkAZyK4HDlQd/cgDKPyzmA1tKera4BBbIg5c2agq2ug1rsxJfBYRcPjFA2PUzRjcZzi8Vhox2hVQawQ4jMAfgCgAcANUspbtPvPB3AXgJkAlgC4zFWqBCHEewCsklI2VrMfREREVDkpZVYIsRTAPwF4EMC/Ani+VvvDMbFERBSm4omdhBBvAPALAB8GcD6AS4UQ52kPux/AV6SUb4OV6rzE9fxpAG4CkKp0H4iIiKg0IcRzQoj3lXjYl2C15TsAXAirk3pCxexULGNYIiIKU00m9uOw1pc7AQBCiMdgTdf/U/vvNwFollKush9/D4CfALjN/vtaADcAuKCKfSAiIiIfUspzXP++yOf+f9f+Pgjgo+O9X2GcgT2MYomIKEQ1S+ycBWv2QkWfyTDwfiHEJwFMk1I+BiIiIiJwnVgiIoqmmkysPpNhDN6ZDH3vF0KcAatE6eNVvHbVMyC6zZkzA7OODwMAUg1Jzjqm4fGIjscqGh6naHicouFxOpnY5cSMYYmIKEQ1QWw7rDEzyhnwzmTYDmt2Q/3+vwYwG8AStYadEGITgAullJGnsap2BkT3RU9X1wD6+0cAAJlMjrOOuXAWtuh4rKLhcYqGxymaiZgBkSZOjAsFEBFRBNUEsa8A+LEQYg6AIQB/D+BSdaeU8qAQYlQIcYGUcjmAzwF4Xkp5F6wZiwEAQghTSnl+Ffsxptj5S0REVBtOOTEbYyIiClHxmFgpZQeAKwAsArAJwINSyjXaDIifBXC9EGIXgOkAbqx2h4mIiOhkpVKxjGKJiChYVevESikfhLWenPu2i1z/3gzg/SW2MSmKh9jrS0REVFuFiZ2IiIiCVTM7MREREdGYcfKwjGKJiCgEg1gHW0wiIqKamhS1WURENNkxiLWx15eIiKi2Ys4SO2yUiYgoGINYHRtOIiIiIiKiSYtBrI2hKxERUW1xYiciIoqCQazCFpOIiKimnCGxbJOJiCgEg1gN200iIqIasVOxbIuJiCgMg1ibySaTiIiopmJcY4eIiCJgEEtERESTghPD1nQviIhosmMQq7DFJCIiqi1VTsw2mYiIQjCItbG9JCIiqq1Y6YcQERExiNWx95eIiKg2nCV22BgTEVEIBrE2tpdERES1xlwsERGVxiDWwSiWiIiolgqZ2NruBxERTW4MYomIiGhSKMxOzCiWiIiCMYi1sdeXiIiotmJcY4eIiCJgEEtERESThL3ETo33goiIJjcGsURERDQpOGNia7sbREQ0yTGItbHBJCIiqq1COTFbZSIiCsYg1sY16YiIiGqN5cRERFQag1gNg1kiIqLa4MROREQUBYNYIiIimhQYwxIRURQMYomIiGhysFOxrIoiIqIwDGJtbC+JiIhqK1b6IURERAxiFZPFS0RERDXlLLHDJpmIiEIwiCUiIqJJIcbZiYmIKAIGsQpbTCIiotriOrFERBRB3Qax+qQRbC6JiIhqyyknru1uEBHRJFfHQWz02xdtaMdl1yzmbIlERETjSJUTM4olIqIwyVrvwKQR0mDe/9JumLAC3BinTiQiIhofjGGJiCiC+s3EFv1dusnkDMZERETjpzAklu0tEREFq9sgNqie2O9WdRvbVCIiovHDaiciIoqiboPYokwsA1QiIqIas5fYYZtMREQh6jeIraCBZKNKREQ0fpiJJSKiKOo2iK0Mo1giIqLx4oyJZXtLREQh6jiI1daJjdBeGmxTiYiIxk0sxnJiIiIqrW6D2MAGMqzlZKNKREQ0flhOTEREEdRvEFv0N5fYISIiqqXCEjs13Q0iIprk6jaILRKhwWSjSkRENH6ccmJ2GhMRUYi6DWKDFlJns0lERFRjbIyJiChE3QaxegMZpb1kJpaIiGj8qCV22NwSEVGYug1igxrI8Dkl2KwSERGNlxgYxRIRUWn1G8RWUE7MJXaIiIjGkRPDssElIqJgdRvE6oKCWiIiIpoYnJ2YiIiiYBCrCWs3GegSERGNnxjXiSUiogjqNojV49FIEzuNy54QERERUBgTy05jIiIKU79BbMkbojyJiIiIxgwzsUREFEHdBrGBA25CAlXGsEREROPHWWKHDS4REYVIVvNkIcRnAPwAQAOAG6SUt2j3nw/gLgAzASwBcJmUMieEuADA9QBSALoBXCylPFjNvlQrUnvJVpWIiGjcOOXENd4PIiKa3CrOxAoh3gDgFwA+DOB8AJcKIc7THnY/gK9IKd8Gq0joEvv2BwB8QUp5vv3vGyvdj0oVlxOXbjK5xA4REdE44vTEREQUQTWZ2I8DWCilPAEAQojHAHwawE/tv98EoFlKucp+/D0AfiKEmAvgB1LKLfbtWwB8tYr9qEhQ+ziSzqF/KIOZLamJ3SEiIqIxIoSYCWAFgL+WUh7Q7guqkvo3AFcCOGo/9Fkp5RUTt9euGHYiX5SIiKacasbEngWg0/V3J4CzS90vpUxLKe8HACFEHMCPATxZxX5URJ/5UP3VcXwI37hpWaTnEBERTTZCiA8AWAbgbQEPCaqSeh+Ay6WU59v/TWgACwCxGMuJiYiotGoysXF425kYACPq/UKIFIDf2fvwy3JffPbs6eU+xaNvMO38e86cGZje0ui5f86cGUXPOe206ZhzanNVrzsV+R0L8sdjFQ2PUzQ8TtHwOBW5BMCXAdyn3xFUJQXgNgB/DOCtQojvA9gM4KtSyp4J2WMdo1giIgpRTRDbDuBC199nADis3X+m3/1CiOkAnoY1qdOnpJTZcl+8u3sQRhWDVFPNhXLhY8f6MeAKagGgq2vA9zWRy1X8mlPRnDkzfI8FFeOxiobHKRoep2jG4jjF47GqO0YnEynlFwBACOF3d1gVVSeAa2CVIf8SwM0APjtuO0pERFShaoLYVwD8WAgxB8AQgL8HcKm6U0p5UAgxKoS4QEq5HMDnADxv330/gL2wxuEYqAF3ZXDUUNg8SbqG248NoncwjXe9eXatd4WIiCZWYJWUlPJv1Y1CiKsA7Ct349V2BnQNdgMAZs5qZoa9BB6f6HisouFxiobHKZrxPk4VB7FSyg4hxBUAFsFaKucuKeUaIcRzAH4kpVwHqwf3TnuCiQ0AbhRCvAfApwDsALDB7ik+LKW8qMr3UhZPQBo9ij0p/GjuGgDA3O/9aY33hIiIJphvlZQQYhas5e6ut2+PASi79KjaKim1xE5f7zArEUKwUiM6HqtoeJyi4XGKZiKqpKpaJ1ZK+SCAB7XbLnL9ezOA92tP24jCBIS142pjDdOMNGlTTVLGREREYySkSmoQwHeEECuklKsBfAXA/IneP3tep3HtMz7WM4xlW4/gby8815lIioiIppZqZiee0vQGMlKDydmJiYhoChJCPCeEeJ/952cBXC+E2AVgOoAbpZR5AP8I4DYhxE4A7wXwnQnfURXEjmNze/MTW/HMigM41jsyfi9CRETjqqpM7MkiamPJELZ81z6yCRf8wRn44Hln1HpXiIjqipTyHNe/S1VJQUq5FMAfTcjOBSjkRaO1uIZp4oZ5m/F/3v9GvPOc0yI+x/p/Nsf6KiKiqap+M7GuyNU0zWjtJaPYsm1vPYE7nt5R690gIqIpQJX3dvenIw3zGU3nsW3/Cdw6f2vk10gmrNfI5RnEEhFNVXUbxLqZqDyGXbn9CLa3nhjjPSIiIqpf970o8cKatsiPL6f8uCFhXfrkcuyZJiKaquo2iPU0eFHLiX1ayWdXHsTCDe1js1MnGYNjiImIqAzueZZ2HeyN8Izy25mGpHXpk83ly34uERFNDnUbxLoZplnxLBKmaSKXn5rBWjXLIEyG7RMR0cklVubiBaqZKWeS4aSdic1wTCwR0ZRVt0GsO2b98vVL0DeUKfkcv5jMMMwpO64mb4zvfjOIJSKispS54k0lFT8MYomIpr76DWK1EqSjPRGm2vdpLE0TyE7RIHa8M8h5BrFERFSGcldtNStoZ5J2OXE6w3JiIqKpqm6DWH0YTZRZEP0eYZgm8lM0iK02yLzkqkW4Yd7mwPujHFMiIiIl5qoLjlIiXEkz1mDPTpwpMSbWME3cuWA7DhzpL/9FiIhoXNVtEKu3e3q85RuA+d1kmshO0RkOqw1i84aJLfu6ce+LEt+7feWYb5+IiOpLOWNbgUJbXU6fqSonTmfDg9gT/aNYuf0obnkifPmejXu68M2blnHdWSKiCVS3QWwpEWNYGOb4jy0dL2OVQV68sQPHfMqxOSaWiIjGUyXtjBPElignTsTtpXhKvMaDL+9B31AGvYPpsveFiIgqU7dBrJ5p1SeH8Jsswi87a5imp/fVME30TZGGbLwzpYxhiYioHLEyU7GVTOyknlFqYqd43NqX0oGydX+5WWQiIqpc3QaxOr0djNq7a9qzE4+kc+gZSOP5VQfxzZuX43hvhImiasD9vsY7iJ2qGWoiIqqNcuPASpox1Q6WKidWHdelrgecZX7K3nsiIqpUstY7UCulxsD6Z2KLt2OYgJE38dN71uJozwjeevYsAMCJgTROP6V5zPZ3rLhnUh7vCalYTkxERGVxxYH7OvrwxJJ9eP/bX4+B4Qzecc5pRQ+vpJ1Rz8mUDGKt/0ft8GUmloho4tRvEKuNcC0KYn3iO/056nl5wyxaokff3v7D/TjnzBmI17iVcweuY52JNQzTKb8aj+0TEdHJzd1CDo3m8MyKg3hmxUEAwA1f+zBmTkt5Hl/JLPiqkzqdDe/IVdsu1ZY5GVvOyE9ENGFYTmzT26hyMrG5EhnN7a0n8PN71+HVde3V7OKYyLrWhq0myPTr/dbXy2UMS0RE5QgbExvUBpcraiZWPa5UtlfdyzaPiGji1G8QW0E5sR/DNJHLhz/2mD0+9nD3UPT9Gye5nLucuPIW12+8qx7Ms5yYiIjGSiJeHOBWVE5sRgxiI76Gulww2eYREU2Yug1i9aZGz0r6NUa+y+5ECHZVA1jrUmIAyBnucuLKx8T6Be76GnkMYimq430jFZUFEtHJpdxmspIS3nyZEzuVfAWWExMRTbj6DWKLxsDqmVif5/g0ZVHiQNWw1TqITWfz+P7tq5y/qykn9nuuHsRyTCxFcfDIAL5z20os3thR610hohoLLycutCk9A2n8eO4anP/fiOYAACAASURBVOgvf0k71d6XqqKKGpM65cRs84iIJkwdB7Hev/XFzH0bowozsSqrG6vgaBuGiSMnhks+7oZ5m3HDvM2hj+nsHvK8haqCWJ9xwEXlxOyVpgjU91se6q3xnhBRrYV19bqblM7uIbQdG0R712DZr6G2U3LpnKhL7antsckjIpowdRvE6oozsVEndvJvtdwBovpnJZnYp5a14vt3rCoZyG7Z140t+7pDH9OQTHj3sYoxsVO5nNgwTby09hDSmfBSMpoY6mcxSb8uRDSRQppJ9ylCH9dazulDtc+lOlrLHeIwWds8IqKTEYNYm56V9A1ifZfY8d+euzEzqygn3ne4D4A1ZrBauaJy3yrGxPo8V5+deLKWE6/bdQwPv7oH85fur/WuEAq/C46JJaJYaBTr6hy2m5tMiWVy/Kj2vfTSOdG257fEzvBoDnc9swPDo9my94+IiEqr2yBWv2COMrOu3qCF9eLmPJnYysuJkwnrSblc9Rf4+nusrpy4+Ll6kDye5cSyrQfbD5yo6LlDI9ZFRamZKWliqL4dxrBEFNbX6z5FqDY8nbPO4+V0EUddOidqG1YoJy48/oU1bVix7QheXV/7pfWIiE5G9RvEan9HmdipaBshJcfu7VUzO3GDCmJLrEUbRVEQW1U5sU8mtkQ58db93egbylT8mm6/enAjrn14U0WBcsbez2Sybr/+k0qMmVgiisB9ioi6TI4f1TaNWSZW/d/VBKr90ofxEBHR2OBVvC3aEjv65E/F2/ErU6pmTKwKtPRS3Uro41hzhoEr71+Py65ZXPa2fGcnDsn0GoaJ6x/djGsf3lj2a4U50DlQ9nNUsN3AIHZSYCaWiJSw2YndqiknzkdcEid6J2nx9lRnaaqB7QwR0Xio37Or1jbps+1Gabz8MkeFHl6j6La4z0LtpahMrJ7lrIRfkLm7vc9pbMvhl8UNy8Sq144y07JucCRbtO1TpqcAADsqKClWWWR1bJVFGzsq2r96YZomFm3sGPvyOLOwfSKqb+GzExfPNVFqrVff7UQsJ456SlKbcW+vkImt38ssIqLxVLdnV71tijSxU4QxsU4m1hXkqQmhKlkmVmViK2modUUTO01gObF6fFA2etOe45i3eK/ntj3tvcgbBr7266W49pFNnvvUWOHRCmYY9svEGoaJ+16U+OV968veXr3YsPs47ntR4oGXdwOwLtJGMzkAVqdNpSXvhZlCx2Y/iWgKi9hOOuXEufJnJ1bPHasxserF/TKxemcpERGNjbo9u+pZHz2IzeVM9A16F1HXn+O75I6h/m8WPS5qmZRbMmE9p5JgTaeCjP930dsBlB4PtHrHUXT1+M+K7PfcsHViVQAdlI2+8fEteH5Vm/N3a2c//uf+DZi/pBUAsFtbQ1QFopVkqLM+mVi1ncERziQZpNueIbu50Rrj9d3bV+JL1y0BAHzr1hX46q+XVrRd9V1iJpaIwmYn9oyJtc8b6UwFlUSRx8RGnNjJp5w4a3c8T9ZZ+omIprq6DWJ1ekNz74sS37x5uZNpAnwmg/LN1vqMiXUmdqp8/9z7oYu6Np0K3s49cyaA8CV2DMPE7U9vx3du9g9MomRi8z7lxFHHBffbE0AdOua/kL3qfa9krLAKqGOuD0Rtj4KpagD1GfYNFibp6hvMVLzurvoe1vpSL53J45s3LcP21spmvSai6oXPTlzcOVzZxE72/0uuExtxg045ceEmlYllEEtEND4YxNr0QLC9ywqe3BfmelPkm4n1ndipsjGxG3Z3YV+HtU5sWCY2aiCnyoebGhKev/2oIPV4b/RMrL4fnjGxJTKxznPUckQlZqxVk3lkKwg+1b6YPvt3slu+tROLN3ZUVPqbrmAClSjU9zAs62GaJmRbz7i8vtJxfAh9Qxk8sWTfuL4OEQUL7ef0mZ04XUEbELmcOGIAqh7lLSdmJpaIaDwla70DtaJfL0dqaPQxsT7P8Zu6X/XOljs78c1PbHX+PZoOCWIjBmAqyEylEkX7qCt1PPyCoKJ1Yj1jYqNlo7NZA42pBOLx4P3IG4ZzezXlxHmzOIitIlk+6Zmmid8+uxMAMOeUZrzz3NPKer7KxObG+KKsUE4c/JhX1rXjoVf34Bv/8G68+/dPH9PXL+yH9R1IxNm3RzQZuU8Rqn0Z3yV2Igax9uM8HaN2p58+aSQREY2Nur1aiz7Wxf1vfUxs8ORPntmJtexiJcLKiaNm1dTjUsk4YjH/fYy6Td/ZifXZj91jYvPRMrHOwvX2sfILUt1LKlQUxNrP8csUn8xRrDubX8mFn3qO3llRLfeY2BP9o76TmB040g8AGBiufMxy3jCwctuRwDJ69Z1OVFP3T0RVCWsnPUGss05s+R2QhqvjLOxawH2mCCs9dtaHdz0mzXJiIqJxVb9B7Bg8ya9t8itTchrJgEbQMExc9eAGbA9ZLia0nDgkqOgZSDuTFangI5mIIxGPewJRvbdYX1NWl/MJBHI573P8gsRSAUJGG3fpVyrtXhKoqkysz5jdSnz+yoV4+NU9FT8/qlXbjwSWdwNWR0FY+duQa9KqSrKpaddEJX5LXVRKffcM05og6nptJmoAGLErEZpSlRePPLviIO58ZgfWyy7f+9V3OpFgEDsVbNzdhdc2ddR6N2iMRV1iR53CnE6wMl7D8GwnfBiDEiWj6j6tZllOTEQ0ruo2iI3a4rnbt+IxsSETO+WLG9ugtmxwNItdbb24/antgfsRFsSGZU3/85bl+PqN1uRM2byJGKxAMpGIeYM4LQANm/QJKH+dWLWPpbLRqlddxbp+Wb+sK1NXURCbLc7EquA5bGZMP4Zp5edfWnuo7P0oh2mauGPBDvz83nWBj7n06sX45f3WEkGGYeKJJfucCbIA63umVDIm1p0Bz/t8tpVSGXv129nd3lf0GFWJ0GSXwldiw24reA3qSMk5mdj6PS1OJUu3dGLRBgaxJ52oY2KrCA49QWzIdtzNYFjHrpOJ9bQphv08lhMTEY0HXq2V4LdUjt99Ss5nrI1TuuQTOb+87hC+feuKwO0pYeXEpQI5td+5vIFEIo5YLIZkPOYJRPWGttQasn69y/oMv77rxJbIxDqlvvZT/TKkaXcmtqIJiop7yCvNxI6kgz+XsaT2tb9EOe3+w1bZ7Y6DJ/DMioO490Xp3Dc0UtjXUp+vYZpF3wl3ma/7vkzEjoQXVrdh3a5jRberfQnbzkjImHA/fUMZLN182PnbNE202TNdB33WqiMjWUUmdvHGDqzZebTi51N0OcNg1vwkFLrEjvvfZVSApDN5XPfIJnR2DwHwnvvDMqWmz5CYsMd5JnZS7UwV67ETEVGwug1i/QJKP0ZI2aRfk5bzGQdjOJmm4sc/9MoeJ3DLhzTKYYFq1GxkLm+gIWldICTiMU+2VW+gw8pNj/eN4IGXdxfdPqJli70TJ4WPN1S3q0BJBcD+S/morGnwe396WSv+6/aV2NfRFxiMedb0c4JnE8Oj4YHi0s2Hcd9LVnCoSrXH+1K63N58deHkPj5D7kxsiUz7bfO34dKrF3tu8waxxRmHUh5dtBe3PrkNn/jPp9A/XMgQq99K2BI9I3YnTqklMZRb5m/F3c/vQnffKACg7Whhqaag/VXvr5oxsfe+KPEbn4qKbM7AjY9twVUPbqhoRu1q9Q1lcPMTW0t+tyeD1TuOer4fQfJ5k1nzk1DoEjs+5cRRbD9wAttaT2DeImvmcdOnk9n/9Qr/Ds3E2v9/bWMH2o8NIpvLc4kdIqJxVrdXAFE7cT1BrL6NkGVmPBlc18Q1oa9V4WzB0Sd2MpFMWB95IhH3BKrFmdjgbT61rNX3dv0C2TMmNu8d66qMpHNIZ/NORkVfliBsYqeW5obAIPbJZa042jOCX9y3Hlc/tNFznx4o66/z/TtWhQZLdz+/yyljHLQzo82N4zvRd6kxyjr1aPfhHnSNiS2VHVhvl966v7NBmVi/7Goph7uGCvtiB9RhWe1R+76oF4Q9/Vbwqj5H9X4A+E4cBRS+Vx3Hhyp6T2GO941g097j2NXWixMD6THdtp+rH9qIH89d4/y951AvNuzuClx3ebIYHs3h9qe3F/1m/eQNs6qsOU1t5ZQTq0qmRp+Z+cM2EzUTq+xu78OP5q7xTEBXamgOFdvW2o2Nu/3nLiAiUuo2iI3K01AWTezkU05sZxxznpl/fbZV6rXKuC9yJjZnFILYeMwz3jSrBTV6OXRrZ7+TKZsxLeW7/aFRbxDiKSe2j4s+JvYHd63G13+9FEk7o6ICCXVs/S4cVKA7rSkZ6b3vae9D72AhcFDvIx8QxPYPZ9E3WDoTBBQCw6hBbDqTrygTV/YyDT5fF8/EThG35z4u7kyp+/kPVTCpVaNrbKsKqIdDgliV5ffrOAKsMclb93c7f6uHqYvQ/Yf78IbTWwAEz8ysgtvO7mHc+uS2KG8jMvdxDFsua6zsPNjjlE8DhfHQYSXbOw6cqPl6yeq30eHq5AiSzxucSfokFDo7sbs5LrOcGCiMqTfMwrCBsI4xo8wgVvFUmuTd7UyeY2QjuO6RzbjJtcQgEZGfug1iozaA3hjWDLyvcJs9NsZnduJSL1npuq1llRO7glhPZk3bhruh7R1M42e/W4fbn7bKJGf6BLHTGpMY0YJYv8l/9Oq/noE0MjnDCWAyWpZUn3DKeoydiW1KRh7L6g7g/LLl+sV9t53JK6XcIPY/rnsNP7xrTekHakplYvXvs/quui8HB0dyTiY8akbT3TGRyead4DNofyqZqTjKmr/ZEqV5D7+6B9c/urloP1TnTC5nYMa0BsQApIPKibVy5msf3oiOrtKZy+VbO7Fy+5HQi1N3CXPY+PbxoioGgkqp244O4JqHN+HRRXvL2m42l8fFVy7E9Y9uxm+eqj7wLyeIzhkmEom6bcJOWnoMe+bsac6/3b/+ok7kkFOPal8aG+wg1ihUJUUtJ446bAIA+of8Z4L/4jWv4ft3rIq8HSIiCsYrgBLMkExs2AW735jYqOP5/IRmYqNm1fIGkslCObEniLUzxwPDGeTyhqf3+ES/lcWUh3oAAKmG4q/NjJaUZ8wl4D/mVM+cNGqzzapgMmxMrMrEtjQVlxM/v+ogFm/sQEtTUnuO9TjTNH1nJ9a3o8ZSlqLKxqY1Rp8191jIMjlBSo1hLTpO9ltzZzWGRrOYNT3l//gA6uJv2ZZO9A9nMd0+rkHPD/qO67d7lncqoyzQb/vuoFDvMFKdMzn7ojWVSpTMxCrbD/Rg96Hekvv022d34s4FO0LLod2TnoXNNA4A+w73obWzv+TrlkN1tuiTrykqc3SkOzgDapom9nb0ec57qmJh6/5urNl5rOoxt+5zmbv83Y81JpaZ2JON/omee+ZM/OP/fov1R8RZhXV9Q97KGsMsBLFh5b7u8005k/gNeDKx3u0fj9i2ENWD+Uv2Y/6S/bXeDZqi6jaIjdr8hY2JjTqG1VnHror5HSrJxOpBdj5vIhkvTOzkKQ+1t/H1G5fhlie2eoKmEwNWo6uyb35ZuJnTGjCczmlLFxTuP3JiGEDx7MT6BYvfpEs61SPuV048b/E+3PuiLMqSZ9xrnNq3qWNqGKbnogOwxjC6/eapbbj4yoVF+6IC91RDeUu/DI9mywpmS2Vii5Y3CignVtnIqGNsVVAy97mdAIBpTQ0AgsfUBt2u75/7wrFUEGtqF667D/Xi81cudErEe1xjTNXFqnqOCrZzeauUvjEZx44DJ5xZnN38gtuwEudyHuvO9AeNyQWs790v7l2Pn/0ueCklXS5vYMHy1qJOJDcVEAadK9S5LB4yUdKBIwP45X3rsetgT+G1tc9OjbldvKkjdN3rIO796yrx+8gbBjOxJyMtFRuLAbNnNQHQZyfWnxe8SbXUmHs+hAa7Q/c7t61Ee0DFReVBrPV7a2lKVjyx05Z93Vi5/UhFzwWs4QHfuGlZWZUf62VXRb9bnWla5+lq1xCfKNsPnMATDKRqYsGKA1iw4sCYbOvV9e245KpFVSWMaGqp3yuAiN9x7+zE2iZCtuHJNNn/1suRy5nwIW+YgUFzUFZMbzw9mdh4DKMBs81u3tft2X8VJKjX8RufOXNaCqbpLcl07++ijdZESPrETvrJJuOz/I1OXei2NDcglzcKpaOuC2D9gkMFEe7HqP17dNFePLvyoOfx3f3eyXfW7PSf6EddrES5UHEHLz/93Tp87zcrSz5HcR9zvxN0cRCrxiAXbhsczaKlqQGJRBz5vIGB4Qyue2STZy1ZJWV/T4a1EnGV4Q7K/gcFx/r+5Ty/j/DfgTsANEwTC1YcgAkrqAIKlQJAYRZi9XEUvrMmEokYUg0JtHcN4ef3rsPFVy5Eu2vcqF+ZsXr/B470+x53M+AiV/+tugPkoEzs1v3d+M5t0b8TytItnZi/tBUv+6xVrPbZycQGBNDqVBSW2FQdBIe7h5HN5fHyukNYqy0npI7/vS9IXPvwJuf2geEMlm/tLPle3J910L4q7k45Onnon2gs5lp0x/WzKudC1TeIdXWALNl02Pd57pcop0OrfziDhmQcLc0Nvm3D2l3HSgZ4N8zbjDsX7Ij8mrp5i611wju7hyM/55b5Wz2/20q9ur4dVz6wAVv3Vx8QT4RrH96EZ8YokKLaeeDl3cgbhWo7OvnVbRAbeXZi12+haImdkI247wvKxOoTIYXJ5Q184apFvj2zQdkVPSvmmdgp4c3EZvOGNolF4d+q1NAp0fRplGe0WGWq7mxQlMBOPyYZn1Jfnbq4bXGVtm7c04UvXrO45HPcx0otAbRsS/HFdZRyYvdyPH77axgm1uw86hxX97jcYz3llRS7P4/XtAuunoG0Z3Il0zSdzPruQ71YteOI/fo5TG9uQDIRQ94wsWhDB7a1nsAr69uLXq/BFcR61/u1/h20HE5Qx4wekLg7Xkp9T9zBYd4w0WdnYNU+qkoBAE52tjAmVsvEauXr+11lu35B09BoDvsP9+On96zDc1pHB+ANSN0Bv55t9U7s5P+7P1BmCbFpmhhJ5yDbrMxoS3ND0WNUY14oJ/Z+PplsHmt2HnU+g7B1nFVG52jPMFZsO4KHXtmD+UtbPY9pOzbg+9ybHt+K3z6705M19+OZSCxrYOv+bud1t7V248vXv4YV26zfK9eJPTnpY2LdX8nQMbEhVAdMOpOHYVrdyapD13pN/+9RUCdVKQNDGcyY1mAtZefTSXfbk9vGPcBTQ6H0juOJsKvNGoYxOBJtgkSanB5duBdb9h2v9W4AsKoELr5yYVHVnJ90DZaxo9qo3yC2gnVii7YRmol1Z/y8M+4qQyXGfPlZveNo0W1BWbFs3pt1tCZ2UuXE3jGx+bzhCVbcwciAa5KKF1a3YYWdUfn4+852bp85zbqAdl/IBwV2YX+rMXt+z73izlW45/ld6BlMW73cdmlrNmdgQ8B0/O97++us7WaLM7Grth/FHU9vR4PPGN+giZ08ZbB508ne+QViizd14DdPbceSzVbQ6V52oVzuoO++F6VnP372u7WeTHHe/qwBKwi742mrN39wJIuW5gYkE3ErqLMv4vwy+ao8ejid8wRqKhjyy96q1/YTmon1+54EXDyaJtBrj8NUs/z2uDKx6rFOJlbNFm5n7VJJbxDblEoglzewXh7zLfMdTuecTpw97X1F97sb1GtcGQw92xolExsWQPp5ee0hfPn6Jdi0x7rI8PvNqPdUmNjJ+9obdnfhN09tx2F7LGzYGFO138d6RnzPQ4D1Hfeb1E6N8dUv6PXvkft7cvj4EK5/dDN++6xVyr525zGMpPOYv6TV3hbXia0HsVjMCTK9Qwuib2PENXGg+n66l2cK+hpVmokdGMliRnMKyUQcezv6cPDIQNHv82jPcOja2Epn9xAWLG8tmblNZ/P4+b3rnKESqpO23PPKWOgbss7JI2XMxL5oQ3uksZH7Dvfh4isX4kTEyRfLMVXKnyu1ee9x7CyjXPyFNW24Yd6Wil6rdzA9pmW9z6w8AAA4GiEBUKqKh04evAIowXNBpt8XcWKnwEzsSPmzlPoFHEGZWPcPOZs3rHJi9+zEeibWvc8+mVjAKr092jOCRDyGaa4ZedWyO54g1uf4eMbMmmZRAOMsEO/z3M7uYWzY3YUDnQN40xkznEArmzN8A6H3v+N1+JePvRVAoWdOD/hX7TiKxmTxeNbjfSO+DZp7hsq8YbjKnwsdFcu2dCKXN5ysYL8ddA349EpHnWBJv/hX353h0awT1BX2sXgZB8M0MTSaxfTmpLW8kmvNYH1maqCQ5RwazTqZsHgshr/7kzcDAO58xr/MLej96BlAb2dA8XMOdA6g7aiV1XNfCOUN0wmk01lrv3oG005WXgVa+hJN+YCsXVMqgc17u3HL/G2+42RHRrPOjKZ+Y8v6Azom9Me6339QEKsH0WFjZwFgtV3Kq7atfs/u760aMjDgysSaZiGbrY5lr50h1S9409k8XtvUgXQ273QatHcNQrb1Fk2eBljfPfd7V6+rfp9p13GQbT34xk3LPB1Q7nPZcNp6rir5VpNsjbjWDGYm9uSjZ0WtINb6t/uUXM7sxO7vn/p9uMuJg7KV7jYxKBPr186lM3k0pRJIxGPoHczgJ/esLfo9P/TKHnz39tLDB657ZBPmL20tWbm1t70P+w/347HFez37XovATFVc6BNqhbnvpd2RxkYuttdo397qDcZO9I9i3uK9VQVOJ/vyR79+bAuujlguXs2xON43gstvXl40TMuto2uwaLhSGNW++V2v6Gq9VBxNnLoNYqOe59wX2EXLmEQMYnMBjYnfRXGpBienjQ98cU0b9vpkiDbs7sK3bl3h/J3NWTMOu8uJ3RfWuZyp7bMrE+tTvpFIxDwN/ww7E+tuaP0CS09w73OSXLShA7sP9QaWEzc3JnDw6ADefOZM5yJEvTddMhFHo51lfXJpK7bs6/Y9ufnNtpzJGr6zo3rKkQ2zEHTb+7ty2xHMfW4nXlzTBvforr0dfbjukc3QRV4eSTseat/8LhIyOaNo3d/RdB6mac3onExYJW4qK+/XWKnv4fBozlmj9ZJPnIfXndocup9BmVh9VtxS5cQ/v3cdfnz3WgDei0e/jObQaA4zW1JINcSdx+rjpHN5azkWfT9Ms1Dy5pdpGRrNOd9T1Ygu3NCOi69ciJw9rtiPfsGqvictTcnAiVb0rORgibIp/bCpgNV9bNMZq0NDHZds1sDGPcfxzZuXY+eBE85nq147Fos5zzdNEz+euwa/e0Fix4ETzn6f6E/DhDVrbPH7NjydDt+4cZmns+z2p7Zh5TarvP3wcSv76y5Xy+aLP99c3kT/cAZHe0bQkLQ+Y8O0qg2SzMSedPRwMh7zn7OpqKonpO1UlR+ZTN453yTdQWxAttK9xcAg1uf8peafcHey+HVK9Q9lSrf5ee+5LIgaVnHqjEbPfvmdX+c+uxO3zB+fdVhzecOZpyCoYicqv2MTsz8r/fO++7mdeH5VG/Z3VD6zu9+SflPVoWPlBYm64YB2NwrVibFlr38psmGa+OFv1+CahzdG3qZqf8ImMFTKWQ6LpjZeAZQQuhC6fZ9fJ65fVlM/H/ttO+c81v913cHMM8sP4JGFe7HN1SOpnqdKWJVszvBM7KRf/OUM75hYTyZ2qPikkYzHncYEKGRiwya30W8LajA6jg8FHveu3lFkcwbOOXOGky3M5g3fQCyZiDvZ2v6hDFZs6/S9ENB74U+baV0EdPePwjRNvLSmzbkv4ynBNp2gSO2vylrrpcN+WT6gOEMZRL2/3z/LChxUEOsXaGey+aIOAhWoORM7udZJ9CtHV8dpeDTnNB7NjQnPhZ+fwNmJtUalnCV23N8p93FVQc5oOoemVBLNqaRzm/oqZ12Z2GQ8XtQY5w0ztORtOJ1zPiMVJKq1VLM5I7BEXC8TzGTziMdimNaUxGg2j50He4p+H/rvbFDbdlrPsGuHTb2m+zt1x4Ltnll+07k8jtt/L9t6xDm2ffb3dvWOo7j06sUAgI6uIad0K53JF2WQ3zCnpeh9Z7J5jGhBuvv127uGcOczO5A3DGdt5UFXRYr796mqDXKG4ZQOnnvmTJj2/jATe5LSPtJYrBDFuocB6UFM6GoBqiMqm3fKkN1jYgMzsfZrNDcmMBxwnvBtx3PWmuzucvegQMBvnHjW1dmmdi2oMiOby2PplsPOsIrTZjZ59t2vbVy2tRPrpf8QnFIWLG/Ftv3dvve1Hxv0BBDVBrF+x1Z9VvpdhWC/8lLSkykT+99z1+Dqh6IHiTp3u3s4ZOk1P/GAjgalMGGi/xwKftT3P8pcMqWqmADrmuBgGa9Pk1PdBrFRS2zcQWPQ7MT6ODtALyf2HxPrd4JWP76g3XOfZHe19RTdr16jV2sYszkrsHEvsePZrlaS6/53cCa28PcMe1KZkXQOC5a3Qrb1BJQT+78XADj/Ladbr503vOvz+pjV0lgIYgPKia2LiMLMloMjWd8GTs/Avf7UaQCsyZ16BzN4eOFe576Mnol1zXYJFD4390WRCf9sL2CVVAY1nH1DGWeCKRX0ffQ9bwBQGE+tBzqA1QupB6Yr7OzXdNeY2MJFTnBnykgm55SRNqWSnnFkfkqVE//Nhedaj7OP15Z9x7Fln/8FkeJuTN2flQqqRtI5TGtMoCmVcAJu9fsezeTx7MoDGEnnkUzEimYgzhtm6Fi34dGc01GhXi/ryr7rv43XndLseWxn9xDyhoFszkCqIY6mVBKyrRdXP7QRjy3e53luv7Ytfdv/ce1r+J/71zt/6+P61Xt3X0R2dA1hwfIDzt/ZnIEGu2NnT3uvc2z9LjbX7CqMs87kjKJM1BtOn170nHQ2X/Q4v+VLdrX1Ot9Rd8bZE8TaPe75vOns3xmnWcd3JJ3jOrF1woph1ZjYwu3lrBbgTEiXzTvnvYYomVj7PNXS1BCYifXruHMysXF3Jtb/3LjPp4PTH519TgAAIABJREFUOyzHfn7AMIR5i/bh7ud2YckWq+NadQ45mdiIy6lFvSaav7QV1z1aXFX0yrpD+NHcNdjqCnDV2NijJ4adsfvl8O10VgGS1uarDi11LZA3/Du3o75e31AGz6w4MCXHyar3ffCof5AWJdAfdXXaRJno0k39XoPGrZdTZg5Y3021P0Fzybg/uygVbo+/th8/uWctOn0C9INHBiKtEU+1V7dBbFTuBkBdNI5mcs4sh0Bh/KDboWODzoWXHuAofj3HTlAUcOJ0749fb5N6mhqP6Ww3ZyBnmM66inoGI5f3LuHjPvn7ZQuTibgnUFMzo/YMpDF/aSuuemij7/vLB7wGUAj0cnnTd0ysW0MiXjqITcYRi8WcS/3B4axv1lHv2Xv9aYUgVs9oesqJ84YTMOS1DLresR90AfSbp7bjwZd3F73GvEV78c2bluHbt1kl4epYnTLdyhKrDOyAXyY2ZxQFpk/bgUyLa0yss+6vz+erjlM2ZzjBUVOqOBN76SfP8/wdPLGT9V1969mneN6PmjQi5fMbAlRmzz0DcOH9dnQNon84g5FMHk2NSTT5ZGKfXXkAj79mTRaSSMSLfjP5fHFw5jY8mnOyyGlt2/m8gT5tPPLXPv1u67HZPHoG0rjiztV4ZOFeZLJ5pJJxNKUSTvZWn4ysqJzY9dmq80FrZ+GiRP+JFIJs73t0f78z2bzzPrr7RguZWO21RzM5HDwygNl2VkdlYptcszsHZmLti41//tO3ALAqK3QDQxlnP3a392GpXTnizcRa7z+XN5z9U7/NwZEsDNMsWRlAU48+JjbuysS6FWViA9oMwzSd+9LZwtwP3tmJ/fdFnc6mNSU954mX1x7CttbuwNdVkyiWKicG/NdDHvapaAoahtBpr8GuSngLQZx9fo8wA1bvYDpSVVBYWamaqdy97+qc9l93rMKNj0ebIMgI6ExXVL+AftzVuaBvKIMf/XY1vnj1a05VCQB87ddL8fK64mXI3NzXJEs2H8YTS/bjeJkBXLm6+0bxqs8KAdUoFcT5Vddt3NOFT/znU07nqbecuLzOgMJcFP6/yYEyg9jhdK5wHRdQTuxuL6OUP0t7Bm2/zO5P7lmLKx/YUNY+Um1UdQUghPiMEGKHEGKPEOLLPvefL4RYJ4TYLYS4SwiRtG9/oxBiiRBilxDiKSFEcZf+OIvat+aegEadM7903RJcfsuyQibWJ8s2OJLFz+9dZ2/DeqB+0vVrXNQFbmA5sesk6zdBjHqePuGMGjeaSARkYrWJnfQA83R7sXklEY95LjZSDXGkknHIQ1Z2eFZLyvcE1t03gs32OImiINbOaOv74qchGXeVI2Z9x9cmk973ODCSdU7uarIewCpHfee5p+G8c0519h2wTpxprTFwX4jkjOJyYrXb7rUN1VIoMQBXXfahov3cedCbUV+65TCeX93muU0FnLOmW/umTuRB5cRBPdBqiZ2cq5faL7BXga01UY/KxBYHsXpZujOBSiaPBctbnddQF0hqMjA9O3DW6cUBEQBcdu1rzhhKwJuh2LjnOK64YxVG0jk0p5Jobkw4y9eo35q73DeZiAWUEwcHsYZpOh0FRQGwYaJ3MO18DwE4Qd5oJu9kErfu60YmZ6AhmUBTKukck2PaxWv/cAbvFXPwt/bkWe4MpR4sn+gfLRobpPZPvxhVFyUNyTgyOaPwPmKFSbP0LE/fYAbD6Sxm26X1mZz1PXj9adOQSsYxc1qDU33h3YdCp4cqazzsE8SOZPKe43n387tw8ZULnfH9qWQcg/ZnPZrJO8uRnHFqIYgFwmdTpqmpeJ1YuM6lhdv9qnX8Akr3uSbjycS6ZicuscROS1OD59yzYMUBrNhqnZf82rmMvZydux0LCmL1qinAe57LO0Gs//P185ehtUVhmdiV248gm8vj8puXR1qT1q+9Ua/pngEaAGbPbELvYCZSaaebu+0KG/6jd76qc8F62YX2riFteJQ1x8VDr+xBGHdb2GpnyMvd/3Ld+PgWPPDy7rKzk2H8OiTcx8PvtdQ641+/cRmeWLLP870qd8kadRyDOpb0qiMgvBLgiGut46AJUd2VS1H2V3X29gykcfGVC/HSmjaWF09BFQexQog3APgFgA8DOB/ApUKI87SH3Q/gK1LKt8Fqhy6xb78VwK1SyrcDWAfgh5XuR8UiRrGeYMr1z5F0Hg/YGbQGn3JiAE4PXtCYWL9A7aYntmDV9iOecWJuuRKZWMMInsFYjQsEULQ0xasb2j0X/HrDrK9BmUjEPeXEiXgMzY1JJ1N09pzpvu8vlzfx68e2wDTNosmHGpxMbOkgNpmIOYH18b5R37Vr9YBr0B3EujJKJoA/ePNsfOlv3oUL/uAM/Pkf/x4aUwmMZvJFWa3ATKxRmAwHsJZscPeIj6TzaG5MotlnRlf94qTz+LB2fw57O6yLezUx06CrnFjPYlqZWP8gtqXZGhO780CPM64pZ48pVvvunjU66yojbWosLifWM/qqM+H51Qcxf2krltnLManj1tyYsPcxj0cWFi4ogjKxgFUSpT5Lvdd0yB6z29RoBYgjmTxM0/QtLUwm4rjsU+/CTLuTAigdxALeqgZPR49homcw7WQkz3/L6U5AOzSS9YxRzTjlxN5zhfre5Ozv0htfNx0f+yNr6Sr37/FYj/c78a1bVziZFyWtlTsrnfYFwGkzGpHN5j1DFoImyegdTGN4NIeZ0xsRQyE4bU4lMOeUZpw6s8n3u5zJ5p0LcCeI7R4uetyotnSTsn53F2IAmpuSnrKxdbuOIdUQxyn2pDXq2HBM7MlHjyc9sxOHjIkF/NtUdS6c1pj0TManT+y0Xh4rWrZFvYQ7E2sYpvX71oaSuI2m80gm457fYlB2SK+aAvwzsUHBlD55j9OhqjKxIUHsnQt24PlVVodp0DJ1bu4g1h10uCuC1HnoHeecirxhOm1XVO62a3Ak68xSrxvRzh8qiPUbHhLUAaDL2m3nSDrnrCE+3pMEqWM6HGHCIt2mPcfx/KriGYD9vmvu21ZuO+L5/HYf6kVXb+G7/8yKg552UZ/TohTVwRB0HecXRIfNjaGqeRqS8cA2a8eBQjIgkzWwZd9x3PbktsBtqt/mDnvJoYcX7sVP7lkb+PhKGIY5YeXo62UXrnu0MPP0Ay/tdubvOJlVk4n9OICFUsoTUsohAI8B+LS6UwjxJgDNUspV9k33APgHIUQDgD+xH+/cXsV+VCTqOrEvrimUn+jPOWKX8QRdgKuGVw9wFL2H9N2/Pxsn+tO4Y8EO3BRQeuPO3rqzJwnXQHq/E3bWLjF1MrHaxV/fYAbXzyuMc8nlTcRc7605lcC3/vl85/6klolNxONo0jJS+kXGmbOnefdHu9hOxq0xRHmjeOkdXUMyjlktKTQk4+jqHfHNxOpl3tmc4VwYNzV4g4lpjUlMa2rA5//qPDQ3JtHUYI2v1Bsw97HNuMqY9ZLxGGJOljaTtRrF5sakb/mj/nl19XkzdDc9vtWZqCuZiKGluQGDw1lksnm0HRvA9GneDgYrE1t8/P7mwnMxc1oKSfsYb7dP+qOZPC69erFTduv+XDI5w9m/5lTCKUdX9PejOhPUcVPHWzWg6jsi23o9vy39gkSnvoe+FyjpPKY1JtHUaH1mQQF8Mh7DH7/9dfjlJR90btPLif0+H3cW1F3KlM8b6B1IY86sZvzsCx/AFz/1TjQ3JjF7ZhMWrDiAB+1e/2zOsMuJE0VBrOpVVsensSGB5sYE4rGY54LRfYFx7SPFSyQ0pRLOxFN6EKsufk+d0Yjd7X14wZXlD5p4pXcwg6HRHFqakkilEk6ZcHNjEp/88Ln4qw++Cc2p4iDWRKGXXWVxj54oDmJHMlZZczwWw1vPnuW5ryEZR2NDouhiJZMtTAalet05O/HJSF9ix3Wb67TmVyXrLh9V1Dl6mt3pMmIHfXo58S3ztznVU85rqImdUsnCmssjWZgonLf92p50No+GhDeIDczEDhb/Bv2Wqis7E+saGxrmkM+Y9SDuYNWvbNi9n+e96VTEYzHsclUaldoXwBt0X/fIJvz47rUwDBMrtx3Bz3631jkv6NUjqm3Sj0feMHxLsY+eGMaq7Uc8t+XyBm6dvw1fvn6J857GY83R1s5+3PjYFuTyhnPtFpTlDnPj41swb/G+wNnwAeC6Rzfhpse3eK5lXt3QjvauQoXMlQ9sQLfWgeNZFaDMTGyuVCbWp93xy7ofPDKAn/1uLfZ19CHVEMebXj8jcEzs8q2dzuoJ2ZyBG+ZtwdpdxwLHkqtjpO/LWE7u9YWrFgUuRzjWbpm/Fdv2n3A+t1c3tOOF1W2hnSMLlrdiexnrBk9GxVch0Z0FoNP1dyeA95e4/2wApwPol1LmtNsn1CsvLMDuV+8t+Tj3aMXrN03Db1tS2N3mHfB9qDHhO8NpMhnH53ffjX0dfcjlTRxdncKSeYVArncw47nAO7G+Ccfti9XdRVuzJOIxfH7j7dZjDvUWJhKKx2AYJr4k58IwgX1a7+dPt7Sg4/gQjqxqwov3N+Foz4hvGZNyz9pG9AykEY/HkM+bODytAfLVFhyzn3cwlcCOl1M4esIKuC7bczcOHhlwGrDDKxpgohDENDcmcawh7jTYl8q5yOa8s8MdX9toHZPVKcRisdDJBL69ayaSiThaO/txcGkC2ZxRdLK6c0MzHpvR6Pm8bljXhO6+URxMJTyPv2lLC+52ZZtbO/uxK5XACy2NzjqVAPA/m6Y5pS3f3zodh+z7EvEYPr1zLg53DaK7bxT3rGuyxvINZnBsdQo5w0Q2Z+BL62cUfX8A4PM77nL+vf9wv+eE7v4ufG33XLQdHcTOZByP3xVH70DaKg92Nfy/2jQNw+lcUQnqk7tPwVNzrYl23CU5+5Nx5HIG9iwEFtwzE6ZZmE35QEMcG5tT6BkYxWV77rb2x7X/v3QdAwD42dbpaGlK4njfKLr7RnHHuiY8MasJJ/rT6OodwTd234097b3YEwNMV1vRmozjrNNbYBgmYjGg7ajPhVUM2Gd/H3UnNjTbswVncOnqmdjjs+zU7Rub8Yg9prirZwQ9A2nctKEZfa7xmYlE8fYPLStkYb6887dOmdm3d8zAwSMDkDMaserJwtJDh7uHMDCUheqGisWBXakkTNMKNt2zkX5111z78zOwr6Mfd2ychkemp7C3ow+3rm3EPFe1gfo9+J0bGhrigAn865LpSGfz6OjylvDG4kBnU0PRbMp7Yv4T4vxqQzO6ekewf0Yj+oYyOLS0AUOjVkfMJldnlN93+Y41jegdSOPyPXdbk2OYKPqO3r2m0V67OIeulpTnt56Ix5BMxn0vPr4nZ2Fvex9usn/HG1r+Fp+5SC8AoqlMz8TGPZnYgrCZT6e5qgRUkNnS3IDjfaNO8OWe2En95vWAUr1EYyrhnJNVB4qzRnjAfjQk457y1KBZ0N2Z2IHhTNHkaCoYHRrNwjTNojHDeseeYZhYtKHdub1Uh7DeTri1HR3AWae3OJ177okEewbSaGmy2kx3eagKGGdMS+Hs17V42vhczkSiUAjjyx1EqIzdSCaHF9a04dCxQafa6+V1h3DOGTPwoXed4bxvoDgoGR7NIRVTFWiFY/ffc9cgkzPwwXee4do/A5u0ZWHKDeCiuPu5nWjvGsKhY4NO8B00030U+zr68JY3zHJWY3AH3tvsoRj/8jHv+xgcsb5PevCqeMqJI2Ril2w+jJnTUjj/rac7v5WgTKxfEOv3Pb3nhV04eGQAh44N4uw50zG9ucGzvxdfuRAfe+/Z+Kc/fQs6jg/hLz7wRrywus3z/nuH0nh9alrRtgu/Z+9xd2/f7/fmZ+mWwzh9ZhPecc5pzm3q/LRq+1Fc+ol3ltzGWOkdTHs6y9fv7sKF7z7L97FqHPvc7/3phOzbeKgmiI3D26bEABgR7tdvh/a8SGbPnvBhtEgk4kj5ZB/i8TiA4hNdQ9J6vOk8LuZ5fiKR1R7vX5Y8a3oj+uyGzgSQSiVh2CWTsXjMMzYo2ZDwnwwhbl3kqvfgzsSeNqsJmWze00DF43G7jCsGwJpAJZUqZBLj8dj/Z++94+uo7rTxZ+ptuurVkots2de23DE2LuBCDb3EhLKQ0FsgJBBgN4X8ks0bEkIgu0lIQgLZZMmS5d1s8iaEkNAJzWCKwcA1xR1bkmXJ6rfO74+Zc+bMmTP3Xl1JxrbmySfIupo7c6ad8y3P9/naWTkJ1j6Z6yDBcZcd28OsByQLoH09ZLP+SZI81SIJyFgCuuopYqCqiut+kWvD19IFdNV5b2QZgHsc7IRGWgxJsi0eRf4uyxKkrEQ+hGGYWXBNU1zXhpwPYNJGc4ky6JoKTZWt8yBUVOfOJC5Lbh9Dsa4tlz2lfUGBD7kee4ZhGpX8s0ugcRltRZEcz3zWMKDrKiQ56RhD1mqPUlkaREfXILKGgVLLwQTMqNZOzpGV4L0o6poCwzD/rngwI9jnoaEmYgZpFNnBkFBkGRmrV6lm0QHZ98nxPEgSDMP97AR1Fb39LO3OvjY8O0BVzTZQhmX4aZo5RnMcBt2vqG6ZhabIGBhK46OPe2i2uzSiYyCRRjptlhGILp0X0ymTNecXTVPMd0WSTCElVTwHmtcDALkH1vOiWEGHgK4izUbQrW0V2X1NJFnyrHUNWMcmww5brb1qaqLiMfk45JCrJpadN71oep09QwgHbfuAzI8VJQFsQy9VzGVZF15zLjlGwFpX0xm7f/hQDjox2T+7X1HfZ0WW0N2XoMbyTT9+HumMgXNWTaXbkPnn4ac+RHdvElv29MDIGrjnpjXCsWeyBn7zNzvUlS+zRK4Hj77BFL75q1dx8UkxLJhejW/96hU0T7BZE929CTTVmNdZlIlVFQmlEd0xF6YyWQQgtnMIRPWcA0NpIePtvj+/g4poADMnV9DrwFNVB4bSyMrmMdl5xe7xztTgZrIoL9EdwYyxoBMTTY59PUO0Y0SuTGz/UMplLwGmVsne/UO45+E3kc4Y+LcvHI2u3oQwe0yytcTRG0qm8dj6HZ6U08FkBpoqQ5bdWhIi/OrR9wCYDlG+mliRwy56B4lycDpjoLYiBE2Rsb095dj3Ext24lNLJwEAaqzuAIlUhnZh2N+XpB0nRMfj63Pbu2wm3GAig6FkGpWlQVqmJLJNH/iLfe4EI+nRC5jisBOqw67Sv3zo6k2ghEnIeAWp8pXsHSoYiRO7E8DRzO/1AD7m/t4g+Hs7gLJYLKbE4/GMtY2zqWkB6OzsG9FNOPaE0/DUlqphfefSk2dh5bwGXHrHk47P502rErYJmVwXxTmrpuKehzciaxhYMbcel51iZw3+9soOPPSEXRd46cmzcP9f3nXs44LjpqMiGsCP/9fk9quKhHu/vAZ9gync8MPnUFaiY39fEromI5nK4q4bVqKnL4mv37/esZ+LT4rh13+N45xVU3HKsin476c+oLTCn960Cs+/vQe/eSxOt1+zqBHr32lDSVhH274BrFnUiItOiOF/nvkQj7y4DdObyrCstR6/fiwOTZVx782r8aPfv0XramZMLEffYApt+waQyRqY01yJ+qowHn/VqcLHmp6nLJuM5zbuxqLp1QjqKv7K9Gfl8W9fPAahgIr//FscL25qQziguiKKovs1a3IF3t3WRX8SfPPSJWiqtQ2f7z74GgwAqxdMwM8ZwYszVjbjj//YQn+f0WoKQQ0k0vi/3z0Nd/3nq3j69V047ogmdPUmsIGpM5o3rQo3rpuPK+982mVY3HPTKgQ0BT/6/VvY3dmP/sGUS5wLAO69ZQ0efHwz1r/ThgUt1Xj+7T24cd08qvILmLThnR39eJVpkQIAP7Um2B///i3HuLygKhKCuoq5U6uweUc37rx2uXldmev5lYuPwLd/bbd9ufbMOVg8sxb//j8b8fr7ezFzUjluuWARHn76A/z9lR346ZfX4PofPof+wRQWTq/GWcdMxdd/uR6KLOHeW9bQ/WyIt9NnviIaQE1ZEFv29Hoam1ef0Yq2fQP43+e24JuXL8VXf/Gya5vLTpmFFXPNKSmZyuDqu57BOaum4rH1O6gBMbnezK5e8qmZmDGpHP/8s5cQ0BRqANy4bj7usWj3F50Yw28ei9NzJugdSOKrv3jZsVA31URQXRbC9KYyPMy01vnXK5aioSqCrXt68M1fvYrrz5mLhdNr8K3/eBXl0QCuPHU2tu7pwX1/fsdVA8viiBk19J4GLJbBv16+FI+t347nNu7GpLoS9A+mPaPuPBqqwtjdOYCLT4zhiQ07UVcZxsYP9+L4xROxbk0L3Y59FiJBFf1DadRZlK7vXLUMt9z7AvbuH8JRs+vw0jttjvGmLNXhNQsbqQEEmIZZTXnI8X4GdQWXnjwLi2fW4so7n8aCliq8Gu/AoiWzAAAdHSMT5JBl6RMJjPoQgAvASfDqEyv+euf+IUxk5nJSglNh1VMTo44Nnng5esRQJhoKbG9owhTwqjlVFckxX4mU5KvLTFZU/1AaJSGN7svLoXnq9V15nVKeOpsvEytSqgVMxylrGGjvGsQb7+9FZ08CnT32mtLWNYg5dB9uJ1ZTFZSENCrWBnhfZzYgIWJgDCbSQhVnwOw1OnNyhSddu38oDcNKEIhq6FlnI5U20FAVcTixYyHsRDKmbV2D1LEWtTIETMbdHQ++hlmTK3DTeQuc7fu4Nnkvv9OGB/++GUtm1br2Q5z2UivwN5TM4KV39ri2IyAlUDCMgpSrWXjVxP7pha2Y3lgmfL4f37ADZx8zleutbB83qCvQVJv1yLYAIoGLcqvELJXOQlNNJ1ZUc/7g3zbT+8o7eSyL6ed/2oSNH3bivGOn4+2POvH2ln345a1rHNt7BdNEtOdte3rx0jt7cO6alpwZ3s79Q7j9/vU4dlETLjxhhud2LEgpXndfgt5jwF2GMJRM474/vYNTlk3x3Neujj587Zfr8dWLF2PqhNKCjv9JYSQFRY8DODYWi9XEYrEwgHMA/JX8MR6PbwMwFIvFVlgfXQTg0Xg8ngLwHIDPWJ9fDODREYyjKBRaE1vId/gI4XlrWzB1QimS6Qyeen0XSiMadE3OK+zE18sBZmaPfanJZEUWqpOWTMLiWA1OXjrZHKMhztqQF16R3bQaTZVpNJDgqdd2QZIkqrxM6t9oJpaheJF9sXtIZ7LoG0hSQSi+JQ8L8rkiS5R2WEhNLGBmqQcTaWGxv6j10T6LyskvZmFOpCZg9Rwlk/dXLj4CgDi6FgqolI5GIpapTNalkEeUeUW9VkmmPZ3JQtcUzwlOliXUlYfQP5RGV18CDVVhzJtW7ThGb39K2DaHQMqT5ab7CmpU2Im/PgR8TSIxGIlRs68ngUQyg1QqS5kGhMZXEQ3QiCF/v1n16DuvWY5bL1xEI6ABTaEK0gQhq8UOAHyLq2sjYO8521OQpU3Re6TKNPLNGjFsfQmh45dyY4mGdVx68izHZ129CbGwU8bA/v4knt9oGhOkVjsSUtE7mMLd//0Gvvvb13M6sPwYiRGoqzKtl55cF3XQmPOBiEGFgyp0TUbPQBLpjIFo2Hmuy1rr6L9JvWpb1yBikyoc27H18GUlOoaSprBTUFPoNScgNbEszl3TQgMF4aBKny9fnfjwA39HZVkS9on1CmLv5TQFaCbWcmKJUcvOw17BMXII8l4mrZIFwH7nvGrONFVxOG0io7bKKhfYb2Vj7W3FWZySkP2uZLOG0MHixany9Yn1ctKIA9HVm3DUtUaCKsIB1dE6S0QnVhUJJUHNYUR7ObHs/C8S19nXmxAGdQFgIGF+7nUPB4ZSVLU+nTHQxonksWs6K3B4tqUST2wAU/CwsHZF9/1pk5OOm3R2DCDXnAT5AXfgIpnKYO/+QVra8+62LodSP+AWOiQ9ej/e69YhILYJ6XAwlMzkXBNM1X8FuqYMuy6YiHZmsgbuefhN3PXQ69jXM4T/ffYjfO+/Xhfaa4++tB3PvmlXIPK1zZqqIBRQMZRIwzAMx3tH1oPSEh26aiZ0iPq4SETqidfsZAr//G9vtwOiJDn10BPv4+0tJi2bPGebd3RjR3ufZ5mAyFH/1aPv4bH1O8TlUjDL2G776Yv03Xr53TbXNolkBlt2u3tLk7WwqzfhOOdBJqjV0T2I/3nmI7z+/l789nGbrcFnzN+0zvvPL2wVjpPH5h3dBfXmHQsU7cTG4/FdAL4C4CkAbwD4bTweXx+Lxf4Si8UWW5tdCODuWCz2HoASAP9mfX4tTDXjd2Bmc79a7DiKRTGCYV7f4WnA5dEAGirDZk/GVAZVZUGURXQYhoG/vrwdP3z4TVx6x5MuZzMYcDuxsiSmLxBDtbI0iGvPmksNacMwhBMOqY/hW+yoikkbFkUo+wdT9NyIqixZ+A04nU/AKWueTGXRN2hOguR7XgYnUSWWZQmqLCOdza9OTPYVsZwrdqEk4yION+vMdvUOCcfBO2lBK5tFrmV1mZldIgsmi1BAsejd9rVPpbNuRT/rsCLxIHJ/Uums6eTlsM1rrV6ZO9r7XI56NKKjdzCZM1qfy8FlEQmqSKYzGBhKORwNdvyKIjmGSgwmMom2dw/imh88g8c37ESp5VCRZy0UUB20FxaserRs0aPJfQ3qCr722cVoZNryhHQVMyaaPWi9hBxYh5vsayiRcRhQQfq8yvSZJ9cCcGZTyCLFO2EAXM5e/1AaZZEAdbQJ0tks7v7dG3RR1a3jR4Ia9nYP4gOG3k1EK0QICsagaQqOO2IiLj91Fi46MYYvnjvf8/sspjGR13BQRUBTsNfKgkQ5EbErTmvFEbEaui1Ba7NZG0QU2lus/sAAUB4JUGGngK7Qa0ugyJKrbRnbxqiiJIAOa7+80JiPQx9udWJxH1ezT7D7DzzbgNTEEieW1F+z76IXXd/gMrHJVMaVifXKmmpcJlbkMFSU2GrbLHOjbzAlnBtLQva8sq9nSBhU7eQCXiLhqUJA1rIaubeaAAAgAElEQVSu3iGHxkBVWRCNNRHs6ujD2x91or1rAD39STpf2plY2XUOT762S5hxzGcAE4Xi8hJ3QS1x5LxqV/uH0jSIn0pn8c8/e8nhjPY7MrFZpNJZtE6pwFpLJZ60rLv7d2/gq79YnzcY+M7WfXhxUxve2247/tf84Bl8l+k7SoLWbV2DdO3nKbZ3//ebuOXeF9HeNYBIUEVNeRAb4jaDKp3JurJspJ1ZQHfPi6SuldiKQ8l0zjpcUluuqTK2t/Vie1svXn2vHd94YL0nTRgw3xlif/QOpLDxw05s2tqFPzAMNlFvVsD5LvFBH7OtogIDpuPJ7oPYG2VhHbqmmLaoNWnkqvkWId/2xLm+48HXcPv96x3Pc1dvgmqEiOYFElR+cZPJfOQd9dff70B79yBlNIr28acXtuJb//Eqtu6xghtb92Hjh530HeruSzqyz2zG+vb719OexKwd7FUX/OHH3sri6UwWb2/pxBvv78UdD76Gp1/fBcB8Bg+UIjMwMjox4vH4bwH8lvvsZObfb8Ip9kQ+3wZg9UiOPVKM5iXmDS5NMevcEimzbYauWnWABhz1B0PcAyxS+zQzsc6FOmsYtB6HZCyIo5u1BIR4UEVGy+hT6E+nw8fCgJ1lJgs+FcMwDHpMsi8SuQ0HVPQOJJE1DPo91aqtEEFTZCRgKpUqViY21yQJ2LWnogyhpspIpDL0ut1x1TL0DabwjQfWI5nKIqArjjpiRZZcmZ+grmIomaFR2KBu1gaKom7EwGbbN6TSWdei2mHVWogCBuT+pDNZoWHGgtA1ewdStNbj7GOmYsPmDuiq7Kk4S8C3DfJCOGCKEfUOplBfaWfSfv7l1bj6+08jmc7SfsFk0iIOYc9AEmsXNWIgkcZLm8xoIqmfItuK+s4S8PcDsCfdoK6gsjSImZMraMQyGFDQVGO2p2GjrCzYY5H3qnfQea3IcyXBZC0Qam5ZScDMfjMGInFoRcEn3tkDzJYTfG/LB/++GdsZYSyaiQ2qNNN71jFT8cq7bZg5qQKPbxCf20UnzMDSWbV45MVt2GoJqehWL+Xlc0wKdWtzJS751Ew8wFB3RWhtrsS2tl6kMwbCAQ26plB6HZ8BB+zrytZsxSaVO7aZ2mA7xuUlOtq7B5E1gDpdQZir9drV0Y8pzPYTqiOYM9UWzKgsDVCjlmeQ+LARi8VKAbwA4NR4PL6V+9sCAL8AUArgWQBXx+PxdCwWmwSzNV4tgDiAC+PxeOHytaMAnoXC99wmyGYN6KqCdMa5ju7t5pxY652LhDSoikwDKywrwqvukQo7OTKxKfrvbNbwdGL5mlgRHbbCUvDuG0w5+kb3DaVQEQ249s1mYvd09iObEjmxzvPn288VamDaTmwCfYMpzJtWhVBAxZlHN+Nv63fgpXfa8IP/NksrFFlCQ1UYOzv6bSdWkV2t+f768nZsb+vFzectpJ9df8+zWDq7DrlAnLOmmhJ09znVVMna6XUPB4ZS0Lk5hl3H2YxeOmM6sdGwTu26ZCqDFzftoWr+j63fjvOOnc58P41vPLAepyybjFULGqlDtW1PLxZOr6G2zIdWRtUwDLpNd2+CZq75ex3fYYrm7e4cQF1l2NQ9YGxG8m/Z0isA7KAhn8mXYN9P0t+7L4+QVFdfAuUlAWQN81y+8cArlLK6vy9Jg0IDQ2k88tJW+j2zNaF9L0pCGvoHU3jdcszCARUDQylK+2XBPpu8g66rMrWRBxMZR6cAmomNkEys3eqNBAz29yVw809ewNVnzEEueNG6CXghNZaFcO8f38YHO/fjrutWOO7nC2/vxvI5DTQ58zerJ+/caVVY0GIz6Qj1/u0tdnmiaRPatgsJnj312i5ccnIp7nzI2a2guzeB/da9qSkP4uV32lBeomPdmhbHNWVn2X09QwhoMn77+Ps4d00Ldlk2Se9ACql0RqjX85vH4nhu427KsspkDXz48X58+9cb0FRTgmVz6nDxqbmv9WjAD2MPA16TP58N01TZrFFNm5m8gEUP5bNj/MsgohPLsltcaE/nAK1DJN+xF3mxMEK/lUGkGVjigHLZVPe5cZlY61wNwEUnJrSMaESnkzT9ntU+R3wMm+KsKqYq6aBAEl8EkdgB2R+5XxXRACbWltCosKY4RW7CQdVlOAWtPrFk4if3VEQdIxNrJuvMxPL3gYgOiNqCDCQYJ1aVcyViaX9c9lxPXT4Ft3/uSETDOnoHUjl7AyYKzMSSPqDdfUlXtpEELlRFBns6mUzWymBnUREN4MrTWjFrskktJTXHxNggzv/VZ7Tim5c6Y10B0btgXRQSGGFp/GR8uiACTSDqcUsW8lOXT8G1Z86hzwx5HkiWMKQrCGiKw0Dsz5GJLQ27nb3YxHLXO84LaQU0OxMLmAbKcUc04ZuXLaV9V0WIhnUcEaultVaAO7gGuPs9A+a5n7umhS6mJSENLY1l1jhUxz75DDNgv//sdSDnf9EJM7B0dp3jfkbDOoaSGSSSaQQ0dybWgO3MB3UF/3r5Usd7XlkapEFIv0+sGLFYbCmAfwDwKqg6eHu4c5Al0MWGXYJNUTX3PNHWxdOJSV9YCaURW92Uze6zmVi2rykxOsm7lEpnHIGvRCrj7cSqMqbU26oPfIYUYDOxSRrkNH9PUSeBBctQats3IAyquoSeOLsjX6kOAVm/OnuGkEhlMG1CKa46vRV1FWHUV4YdWSRFkfDp1dMA2PRFUSYWcArndHSb9cBPvrYr51hIdoy/Jo01EYbFJA7O7utNuBIG7NgdmdiMuW5rqgxVMW2WZDqLdqt2taEqTB1Fgsc37MDe/UP4j7/GTQfVGiuhjPYyQWXDUmQn63PvYJJmyvoGU/jry9tNlh5zD3d29KGuIoxQQHWMmzhpDdVu0SLiiDVZPcx13dZ1COoKArpCS0a8sHf/ECJBFQFmnSW2KLsO/uG5j2ivYQDYtGWfw9adOqEU9VVhep0lyaR1Vwqeb/b95p1YTZVpwHgwkXZkavdbTABdM/9vXmNzDKSE7MnXdiGTNfD062Yg+JxVU11rD5BfJZp/51j1beJMP7Z+u2Ne+MWf38X+voSrPpfN6meyWdqbmC0f2tM54BSftNY8EaUYMOe/7r4EArqCshKzC8CjL2+ngov02MxY9vUM4R8bd+MfG3fjzy9sdbRfYtkehmHQ1qIkW0yeI1mWaKBuZ0cf3vzArRM0Fhi/TmwxdGKIHVmdi1JoqgxdVZBMZZFIZaBrMmTJvbgMDKUdDxVPNQRIRshpqG3eabe1cGViDXcmVlUk+iCqHJ2Yz6by0LiaWJKJNQw3nfj4xSb9hkyc7DmpiuRZE0v2KUlmTexbH3XS7F0+5HJi+bWaOrGq7DAGRE5I0Jr0E6kMdNWs59VVcSslmonNZCllJ5XJOmjdp6+YgotOjAEQ18TaC7EBzaJ4e0FTFboPPmNZGtas+kXnM3AUU7uYr3E5OR9yXQYTaerQEpB7TjKxBGkmM0EcHuJsEuebGEfkmVoyq84hqiU6L8B+VokjyAaPyHPG9/91jJl7xhVZpmNtaSzF4pm1dssq65TINdBUGZGQ6qg3I4aC6L1ljeubz1uAC46bbtbtCp410ffI4jqproTeD94pXTi9GjzINpoqC1UNRUbl/JYqnLR0Ep3bdE3BvGnVUGQJ0bDmMGL4+l/ADqDxzwgArFnUhKtOd7YXIMbYkEUnJplY8kw11UToOEVzRhXjzA9XuXEc4QoA10Egmniw93Dn4cjEMp9nDUOoWNveNYCsYaB3IIlL73gSz75h9diWZRrkA5xBY7bE4v/8ZgN9FwwY5txPW5dkHUbuC2/vwR+e2yIch6bIuPHc+TiXEULjUcbQibsZleBuKwvGYyiVofPB3v2DBQV7+YBmoX0wyTxNvs8GwFgG1LlrWvDDG47G3KmmUCbJFKkeTiy7LhLxNpF+BQsSFGed2PISHaVhndo2XuJDj7y4DRu5tjmsYe7IxKazSKUz9H7qlqhfV28CZSU6KqMBdPclsL2tl9aqshTf7W19dibWYovsYxyVnv4k/XtDVRiDiQx9phPJDGXqvc/YeEQsLxxU6bizhkFFQSdU2fYW+x0A+MpFi3HKssk0m0zOKagr1BnxQiKZQTioOYKYZD5mW6LxNaU/+cPbDlG+2ooQDeCzYyP3kn0P2QALL1CmqbZ+wmAy7aITl0YC1vnJjvaR8e3duP+Rd/Enq76T2AGVpUFhECxfv16eAvzOVkaAkAh27Rtw7SeRyricWNae2Nnej2Qq61pjv37/evyC6TVLnHtR3+jmhlLs2TdgPq8R3cHu3Nc75JhTWEe5rWuQ2j67OvowkEjT5AN7nR95cRv+5ecvYVdHn4sSnkxlHOt1qYCNNhYYEZ34UEYxwk4wxFFMfgFTrawdYD4ARKiHr1UdTKQpPQMQ0xJFdGLWCSHfIc/O7s5+PMXRKSNBjUatqLAToRNzmVmvcyPGt8pkOYnzQqgOy+c0YPmcBvzmb7bKMRmfouSgEzOZ2Fw1broquxYqUSRN5zKxBNGQht3W8djefjydEbAdk96BFB2fpspCQQIyCSeSGXqdU+msg960cHoNNT5EFFpCicpkswXV+QU0k0bnqokN6+gbTDkygVef0Yols2wn1qt2iFzf+soQtuzudTiSvKPPOuKlYZ1GZjMZg14jcm9OXjYZm7buo5MiiZiLnnf2/HhQYSehE6t4fk80ZsB83ohBSu4NfWaIE2s9G7qmIBzQHJS/3sEUArqStx3U7CmVmG31jxOxLVjQTKxl/JE6X8Bt6J26fApef99pnJGAmtdxyD2RJDvqHeLa1miqjOMWN2FOc6XLiBHRpAnIM3eS1e7AC6GAQhfgoK4iFFAwe0oFjl88EbFJ5ZAlCc9bAiYi0ZnKUnshzke9H6+Ix+OXA0AsFhP9+aDu4f7www9j8xM/pr/f91oIIV3B9rY+fPOtEupA7ezoQzpjCGvgL3n7PqTSWWxv68PmJ8zPvvl2FJ09Q3SO/vamKM2i7HpedRhll73/ACSY2ah9PUP49lsl2N7Wh6+/VYKO7kF6zNutfatWr20Wd240gzFDyQw9DvveAcCdb5Vg195+3LdehyQ5DcuPXwhiX+8QjKz53hoAtmSyyGQMZLMG/v21EHRVdvWE5vHAK0H88YEA+odSUGQZqiJh6+7cat6Xvf8AevqTjmzdT9+M4EHr/e8bTNHj/uKNMB62DO/3md71129+AMlUhpY3EMiyhMve/DkAYGdHv6W94aaWsvhANtuy/eb1MNr2DSAYUCDVRfHu3n4kUhl8+LTZG1zXZEcNYDRs9sbe+qzsWI9ve7sEO6xM6Q83BCkF/WcbgujqTeDdkIYX/jeMD3ftx9ZnTYFDwzDwrqqgP5HGX39j7is2qRybd3QjFDAdzFs2mq3biLP5uVfLMDCUpnToa96+D4CE7W29+DikORydLcw1uPkl3VEW9B+vh5BMmWrul20oQyJpX9eOVwKOMhcW177/ADr3D2Hv/iEQq+yLm+/Htj29Qvp1SVhDeYmOne3meDtfDZpsAxK4kQAYwB0bgpQZtGffgKuOlO1l3v1aCAOJNPoGUo5+4bssmjF773/4tITf/TSMoG4mDD5mxMN+9noYAU3G9rY+fPWtEgwl0+iw7tu258xyr8teLsHOjn4MDKVgGOZaM5jI4L3H7fGQZ+GHGyPYs28gJ2tNhK9sCKKuMoLNlojWh0/ZNvxH1vl99LTk2u+Nb5tzCGvz73lJx19+HUY6Y6C9awC9AylUlgZcQo6bnwCefbgUmWwW3X1J9PQn8aEi4bJXyxy92vda3938hESTHcTGuf7t+/Dx3n6h0NuPX9YR0BS0dw3S3vHdFSG0dw3ixjcjdN3/YNd+ZDIGbnwzjN2ceNh9rwahqTL2WHNGnXEygOXDurbFYPw6sUVmYkUPQGNNxOFgkZpYwJzsA6oCSXI3jB5IpKEoEq47bS4CAuVSwO7RyYJVHiORH+JQsq1WCEIBe4EebiZW5+jEJGuaNdjvOsfHNpFfPqceL21qw8q5DXhnq7OWhYA4xqawk70v0gONjkVTXE6sqCa2vjKMtq5BV2aIZAZLQpqjNlHkCJN70TOQpPcyqCtClVjiAP3zT56nTg5PJy5nIshCmX+HsFN+wzyom/eUD6BEw5o5ATERP9758VIavOeGlejqTeD3z3yELbt7HZMw7+gvm1OPR1/ajqCu4tYLFuK1zR146MkPkMnafRRJFH56Uzl+/uU19Lu88yQCceZY55lE+WiUnHHYyLOo53AS+aydoki0rpU6sdbfaM01m4kNqtjZYZcH9g+mhPRags+eFENF1EkBZs+LLycgxwHETizP+NA1BecfN91hxJLgmZcTS46vawo1xMm5U+qkRaUj2XGb4qx61jADZpb27s+vEGZrWbA0zmhYgyRJjho5wK7bEgUNK/1M7EhxSPVwVxQZqmqXs5A+xWaA17ktcSbN9Ykr89Fkx1zomCu4jL9qMRnsecU8pmwZoIoiOW0BgUFBej6zjzD73pFtVKu8hZ/1zR7NptChriuQICE5kHT04M7F2CFIZbLY0d7nKXgngq6rkLhMkq7bfbZ1JpurMf23WX0EXVeESvjZrNn/mqWE5qM4Z2mg3547dF2FqipmP1PNFFeMBDWHE1tVFhKW1zivm+T43DAAxXrOZNlkRWWyBgKWUu9+pgYyC/PWR0lGWHI+H+msM12yo73PrrMOKHStVDgRMF7XQtNMQaOsYUDTFAww5xgOakInVpLINXK+BwGrD7no1TYMsxc9PS6fOLDGnjHMfff2J12ZSR7BgBmA6RtIocLqC0/OmRyDnHsmY9D1rLLMuXZq7LsvSzCY+5ZIZlBRGoCuq6a9S2yQsI7BhHm85sYy7Onsp+PVdbVgaj0ATJ9Ujve3d2Nv95AjAMDugzxnIseY1KaHgioNCmSsd2HPnl7qbJaWuJ1YwMzsZ7MGFYci7xEL8t10xkAkJFtzIOmrax6vtsLsbkGuUSiomvoXQZtlqSoyohEd7V2D+HhvP0rCGibVRe0uHIKAEz8XlZea2fex7uE+bp3YomAYtIUIi6kTynDPDStx7Q+eBeBuEaFrsrkAcRmwgaE0ZEmi6p4iyIJMLNucnEzqXlRdwMxcdVmZMr7Fjq1S7JEl9WixA2YK4cfHGgrTJpTRBtAsxUS0vSxLjgmXd74Cmow+ThtD5MSeu7YFy+bUY/ZkZ5uPiCWKMbkuSpXdvPZBsn09/awTq2Ig4Yx6f+WiI6i40C7GwTEzsUxEmKFVsY5AbUUIvQNJZ02sIgvVOFmQLKbOZR6Jgi2bVdA4S4+dgEh/YXJ+DVUqLjoxhvJoALGJ5ZSCFOKypuesmoZTjpqMcFBFOKjiuMUT8dCTHyCdMeixRfWXLEI5qLWyLOFzn5rpcOKIQUkz4+TZZPbD0om/fcVSpDMGbrd6JvPPuCpL2G9dd+KsEooT+b2hOow3PjDf1XKuhiedMXJmVlctaHR9RjL84aDtxJ5w5EQq9EAWgtYplbj67HmY32ILGvEBi4Am4/jFEx2fEUc3oImvbVVZECcumYhFM2rwnf98zXGuRByEF3EgjnFJDocdMOeBMgEFkuCkJZPQ1jXgYAl4qVPncoQr8wSEfOTFQd3Dfd26dfj1S/b9P29tC6Y0lOKOB1/DTectQKvFarjzv15HKpPFNqZ/9NypVXjro05cfGIM4aCKn/5xE0ke4RuXHIln3vgYT1kqmt+/djlu/skLAEwaHltj9v0bVqI0rOO/n/wAT762E7dfciS+ct/LuPK02fjlI++iqbbEUQsnwq0XLERsUgXauwZw289M5vac5kraqgMwW7c99Pj70DUFmio7+s1fcvJMPPLiNrR3DWLNwkZoqoxn3vgYqXTW7Ds/fwIUGI62JAR1lWHMaa6kSqTDxb23rcVfX97uEKL82mcXo9kSXCN9rQHgxnXzaJu3z9/9rBmclyX89JY1tB/3ynkN+MdGe5zXn7cAr2/eC/m1nWiqiThq8HiwwkVf/9xifPNXr6KuMozvXHkUvT8/+uIxuPLOp3HW0c14+o2Paa3hVy46At/+zQbXPi89eRbu/8u7AICVcxvwj7fMsZ24ZCIef3UnTlwyCZ9ePQ1f/+XLCOoqtu7pwaoFjaivDOPBv9t5xqvPaMVP/7gJt124CHc8+BrOWTUVj760HfNbqvDipjacs2oqPt47gA3xdovJZa/L5x83Hf/1uEkJntZYSvURrjh1Nl55rx2ZrEFb5lx75hx09gzhd09+gB/ceAwe37ADf3huCxZOr8Y1Z87B+zv3487/et1xjuGAih998Rj87ZUdlHoMmP3iv/fb1/De9m5HVhAA5kyrwrmrpuHr1pr52ZNiePmdNrzHZPsAYM7USnzp3AW49I4nMYX5vKkmgvqqCO1RP6U+ii+fvxCqImN3Zz96B1O4yxIiuuX8hXhu48c4et4EfI8bO39NyLWeOqEUt9z7IpbPqccLb+9BJbP9GSubccbKZjz/1m788hHz3t7+uSPx//3qFdqv/ge/e4O+f9+6bAm+9sv1ruPyKI3opi3SVI7P3/Ms/bzeY3vWpjp+8UT8/VVzbS8v0VHSl8Q5q6bif575CIAp0vmdq5bR+7F6YSNOOWoyvnzvC459qoqEbNbdCudHt6zGFd97GoD5nvz7jUfjurvNMR63uAkSJHp8Inp53toWvLutC29+2AlZkrBibj3e/LATqxdMwP97fisAU9zxguOm4yv3vUyPNaGuBOmppo07Y6LJQGCxemEjIkEVj7y4DQCwdLnJAhrrHu5+GHsYyBp2JvbC4229DFlyOpG6qjgMTl1TzJpYLhM7mEjnpY5KsttJNIvYVdx13QrqLOdyeoJWoTvAZGIVLhPLRK2XzKrFd648Cnddt8LOdlE6Mam9FX/XPIab5mluJx6kyjjUrOgRH53mHTbRsQHTqVkyq84VGSL1rJPro46JW0QnJk77/r4kQ6l2H795QqkwgNA3mEQma+DU5VPwnauOcmTT7RriifjmpUsQDqiMOrEBVZWxav4E1z5ZiOpCASA2scK1Lb8NW0vJO7iAOWlfePwMh3MYDjivkSxJjusmyxLKSnQ8sWEnpQCJ6pUd55CDTgwAx8yf4FBFZttCsWNnnzE2eNRQFcFEptZWVBNLQM71vGOn48rTZlPn+eh55n3Ys2+A0qFZiOqpc4GMlQ2cTKh21zRpqoxTVjQ7xqhpvBPrvn40E+txbSVJwmfWTqeCMzLTC5o4HrxjSBz7WZziMI981N5z17bg+nPmIRqxn4uohxObi7ZcXhJw1eP7KBwHew93Ho55nLHhslmzXpVdb6osqnnvYIpS1lkROnaucLTY4TILhBqcNQxIskTn0P395rxeU+YtskZA+2KrzjmJBanT7RtMuVRRI0GNrkOKIlGdBmLIPv/mx0IHFjADXBcePyPne5QLhmHQoDu5+ixjiWXRsNdRVZyBRl1TcP9ta3Hikc5g2/cfstuKVeUQrAPM9ZqArDkGk9VKprOUvqtrCu66bgXuuX4lLvnUTKqIz4Ol8bIiRclU1syO0RIiBR/s2o90xkB5ie6qU/7zC9ugKhKm1Eeha7JJJU6kUV8ZRnVZEI+8uA0vbtqDJbPrHGsZAEed6KQ6+xybaktww6fn4cQl9jUL6gpdowYTaezZN4Cq0gCuP2ceVEV2rE2ke4HG6CMQ3HrBQmt/5r74lkXVZSEH0yUS1FwMQk2VHSJkLNataaFK9BXRAL7+uSMRCphteibVRenYADOgesVpra7gMAFPUdZV+xq8wPXLBezniL3nk+ujuOm8Bfj82fMA2C1uAKfd5xVMBYAlM2uxoKU6r71CMKnWvJf1lWFHW7zuviRWzm2gteOAKexkGAZ6B1OYP60KF50wgyqWA8BVp7dClkxqco2gxR57jaJhDaGAStfr8pKAo6SOtMuJhDTKUgwFFNRVhtHTn3TUbldEAy67eHtbH46Z3wBFloQBvHe27MMzb9gxz3yMrNHCuHVii+1jRJwf1mCTJKeCsKlky2ViJcnRRxUAjVjmgiy51Yl7+pIoi+gOkQORI/Wvly/FVae3IqArdNx2LSxXE8ucz5lHT0VdZRgV0QCTiXXSic2aWFj75DJcik37ZA2QfLWDkuQcRyaTxQ+/cDSuO8uU6RY5ESJ4UR6JoVJfGaYGeySoYnpTmWtbUne3vz/JUDTdDgtpC0RQVxFCU02EtiSZUB2mbXDY7wBmZljXzImZ0okzZp/YU5dPwT3Xr/Q8R6rQyzs2ukIn5Jpyc1LnHf1/OiGGq88wxXZyvQbs4ifKVvP47Ikz0TeYQtzqj8e2gxAhF51YBHLdWAMJcGZiRUINBCJ1YsA00sj3ApqCo1rtGGt9ZRjnrW3BVae3Yk5zJXjkE2pynYPVA5V18POJmhDwdGKhE5unJpZAVUwF7FBAoe+oLWrlvE4r5zXgrutW4OKTZubcZ773m8CRifUwsnPRtGVZQkXU/LvfJ7ZwHCo93HmYfaLNf7PkzKxhQJb4ecpsozOYSNPafPLeq5YjSMDOnXyfWCJ+ZAoY2u8VcXiIoU8MtcUCNhXv0AHmesBvU2KJ8fX0O/txhwKqo2+1aP3xAivAWAyyhoFkymz3Vmo5Oiyzhp33HIFqJmDAgqe0EuiqLAwisyBrj6rIqC4L4thFTfj82XPNvwWIgn6C7g8w78vR8ye4An8ErCpse9cAJMl0ZAjVlHyPpfVGgpqjMwBg1mV/Zu106JqCaEinytglIQ0NVREaSLnw+BkuoUnWia1mHEdia7HzfUC3RY0GEmnssdruiEA+DxAKPrPuTak3HUxyz1jmzBWnzsY158xzrPXhoOoqP1o6uw579w85eu0S6KpM27CJymWqy0LUmSd2ilfwc+/+IcffNNUZhJrWWIrbLlzE7Nu8hqTly6IZ5jvZOqWSnhPrrLLvWhnnzCuyLURKbCNiLGQAACAASURBVINcbEdWxIgEzk21ZOfz1zyh1HHNk+ks+ofS2NczhOrykOlLSBKWz6nH+cdNx9LZdbjitNkATLEoHoRJAthrJvELyiK6Y14jNdSRoEbPPRRQ6Xrczuy/oiQgDNA3VpegJKwhkco45jVZktDePegIDok6NIwFxi2duJi5nW3izDoG7CIL2O1YCHSrJpbPxKasPpssYhPL0TdkiyaIhJ329yddURlRbcyE6ggmVEeoFLY5bo+aWEGmEDAfxFBAoQ+sLexkv9S8KBRZQHnjXOSwl1o1ceTv7OKXzmRREtJwRKwWd1x1FAwDeaX4vY4DABccNx2NNRG0NJbRBeXL5y90REEJqsvs65vPMWCPd/GJMbwS76D0KHY/BHxrojAjnZ9OEzqx5FhMrj6j1ZE9tq+xe0z/ctER2NczhN8/+xE6uodcWXpVkekCmqsfL/sMF5JxJAtBV28CuiYLx8aCpyjng8Q5seSnw4kVOHZETIVvbUT7zgbUnAvUCUtsoaKzj5mKSXVR3PPwm65jF4r6ijDqq8K2MmeBjpirH7XAMCTb8A4vD0mSrObx9vjJs8C/PoosC9t9uPaZszGUDTZC6xUBzxc0qSwNorMn4Qs75UE8Hp/C/PuQ6OHOQ5KYZ4srQ5VlyaEhoGsywgEFg4m07fha3+EdQXbe5tuzpCzGVdYwIMFmKxBV1kaLPXHGymYsb62Hqkp41aL1EZD3k8xJtRUhN2tJlVFeEkBPfxKKLKO+MowBqywlFLAzT3wW2QuqIiOdscUBvUT88iGbNbUTdFVBZTSAnr6kY64IOTLabraVq+2gxxwXDqoIeDiagOmQBalTZ3YJuPAEmwVHSoRI2x4+iGl2FXALQrKaEZ09CZSX6FAVm+5LxkuCFp9aOgnLWuuhazIViyJYZgU9o2GNZiiDukqzcOUlpmjOZ9a24I4HX6PfY8si2KQHCRCw831AU2xBs/Y+bNvTi1OXTxFeM6rjIMjEks/Iu1FbHqIqy8vm1CMc1NDfa2emzUys8xma3lSGf2zc7Wo1RM6DMG28arBv/9yR2NbWm1PokiCoq9Qx4hX3mxtKKb0dsJ1YXVPw3auXCdW92fVG10z17L7BlK23YVGBA5opitQzkMr53pFnYUJ1BD0W5ZokQCIhzWWPVEYDlH1UXqKjuy+JXR19GExkHIyEy0+dbX+n1H0epH89i1LrupO2XZJkJ25aGsto67BISKXXvqY8RK8JK1pZEQ0I7YuG6jBKwzrtE0zeu8rSgOt5KJYFMlyM3zB2scJOAsqdxNGJSYsdApKJJQsKiSICbofr1gsX4XOfsjMeEsTCTjwNL5e+CTuxkwmDjJ/0iWUnEnZMpx09FV+9eLGtRCzbTizrfLLgnQzHyTC4+/MrcPslS+iEKsvOrCZbHF9bEfac7L507nzccr4tDOOV2aouD+GcVdMgW0qH/HmzCAdVSp8i0V2vzCF7/jUVIceCzUduAfa62U7YwFAaWcNAJmu4xLcAsw3NMiZDSJV6BedaXxnG7CmV1HgTBTjsNkQGmhtK0TpFREO2nxuediQCuU5dfYmc9ByCXIuXCDyd2L4vzog1DxpscakTm98PD8OZPnX5FMxkaLWhAgxLHl/97GKcubKZ/q6qMk44ciJWL3TX0LLgr5fovpJ5p5CkqJcTW4hYjHMgw9ucfTaiIfFzlSuoANiZMJ9OfPhDMr1YAO4WO5IkUdEnwKYcsqJBxAhXZMlBCWSfcy86sWEYJnvCOgYRXGmqLcHPbl6FNQsbEdAVYVmLygRzb/rMAnzts4tdz6sqy6go0WEYZtC2ijFYQwGVOjJ8FtkLtMzEmuvytVPzQjZr0ol1zXSyw0HVZePYx2QCA1ypB4H3OutUP2fxz/+0CF+9+Aj7nATObrUlHkNqiUWlGaL9861OKqIBy4lNCb9z9qqpCOgma+U7Vy5z2GjEuYwyKv3BgEKdWDLHzphYjusZ28+rDIY67ZpzXSP7efTl7QCAo+ezZe1mKczaRY12yQpxZhV3i5zlcxqwYm49zl3bgktOnokVc50VnmSODuoKls62Oxu0TqmgzLK2fW5Ksa4peTNwAV1xaF3wzwarESNiTpCAaktjGTRVxq0XLMTS2XUOGnRNeUhoB7LBeEmS8K3Ll+LbVyx1/V3XZLq/XAwI8j6ztOFlrfVYPqceF58Qc80LlaVByLKEn928CteeZT4L7+80nUuRrQjYTimLcgFVl1z3s1dNxeS6KOZOraLnsHKe/ayEgxp9/lubK+m97maoyaIAAAA0VEZo1rnS0e7KvT2f3R4rjN9MbJHyxHbjdGcqnV0QFVlyTEoBqyaWOGWsESei4LGLhSxLwgWSj3LkMjzZsdiZWNlxfHZxZccUDmqOOh7yNwMGdZx5Sp/tZDgXAl7sg9AqyNH4+ia+n51X1mUOU2PAn4sXaNYpx7ZE1Xmilan1zsTaY2YjWKoi5awLoJnYoIqdHbb0ObkGkiThuCOaMG9aleu7VOQoR0sZgzol7r8RBzCbNQU7RGDrunMJ9hCQ806msohU5Hdih+0sccchP4N5MrG2oiiXibWep+FmU9lASzGZWFVxBrk0VcZ5x07P/8UCpiyS1SiE2qtrisMJP215M/79fzZ61pGNFth3nM8u82hpdFP9ATs67asTH/5gfFgHPZbUxLKZ2GmNpVj/bhsGExmXQqiiyJ6BSJcTy9CJSYcAVZHQYWUryiJ6XqYJa0S3WqUIotIb1gBkFVlDukqdQVWRc5ZKEBDxQzLXFSuvRejEuqZgWWu9q5aXnbvZdZGwsngHwtOJDahCJ5OIQZr7dwaTWVRbJTOvv98BSbIz5CwCmuLq2el2YoMoLzFoyzIy/pvPW4D27kHHPBMOqlgyqxa/evQ9xz6iYY0mOUK66iojApyUbEmS8OnV09DSWOboOMEHaQGnYOHOjj7Mn1blYnmdYNUd//oxs5lOi1UmpanutaC1uZI+k0fPm0C1HwiuOqMVDz/1ASpLgzjrmKk4+ajJGEpmUBLSKE14e5u7LjKgylALKD1iwdp1VaVBXHvmHHz/oTfw7rYux7NF3oWbz1uA3oEUPb/YpArEJhVWakao1uetNXs3l0V0lEV0u2NCQEV3XxK6qqCqNIite3od92HVggkYGErjFUu4ynQIkyiLMO9waZBmUj9ixOIA2wHXrP0DoIJnLL2cRVRgQ5ZFbOo63c5yYhuqIrj9kiMBAKevmILm+igWTK9GfWUY723rQkNVGDMmluP5t/bgyJm1wr7RPPNqQUs13vhgLypKA0hYc+WU+lIq+MWzKdj3d6wxbp3YYsC22GEdJd4WlySu9kZ1Ss0HdYWqwomMMNaJlSSxU8arj+ayWR2LDCfsRJ1a1inPsTPWmPASVyGUYz5ySpzYshIdlzCRTDD7YRdH3gjxqqtxjbEA54iMJde5EnrEfMuJJA6LBKdx4BRtsts4VEQDwowSX0tcVhLA/v4kNaTYBf8CRkCMBbnmooWdoLE6gq17eoVUYGI45ArmVEQDWDi9Gqcsm+K5DQvWqBO1LSI4ZUUzXtg4bNFTaiBoXJCEPb9cmVj+XpBrOFxHVJYk+gwMV9iJHput9SkwI03qR6c3lWHl3AbhNhrNxOZ/B/hM7LxpVbjvljUFjWW0kOtd/cUtazyzvJNqo1Z2YviZcB+HFmTWi/WoiZ01uQLXnDkHJSGNZmITKbfD6BWI5A25FJOJtcsYFAwk0qbjmYNeT+YGkePGv5eqKjv2xQrRhQIKXfNUWSqoJpbM68NlufDIGgYSFp148cxaLJ5Z67kt67BStgxPJxY4UoBpm+SiEwN22YmoRMIMJpg04IaqsDCwQNYEQrUGgK5ep2hQRTSAxuqI7cRa12/2lErMhhtBXcW8aVU4islSsomFYEChv7MOAa/Yf/JRkwEAGz909vsG3JlY9vdcrB3iXM60HLtC1xcWrVMq0XqJXW0QCtgU1FJFw4TqCB4XKF9rmpJXC4MH+6zeea3ZV5TU1bLPPAnYN1RF0OCO6xeEOc2VuPv6lXT/PEhWXddkVFkBpR5GcO2zJ81EJpulTuwFx03H65v3YsnsOvzpha1oqnEGUfgWPqxdRMbQ3j2I+dOqMKlOHDxmA/ORoNVaURD4mTvVrdmhKjIWWrXBMyaW0wz4yrkNWDKrDgFNcZwfAWHe6ZrZaeWaM1sxkMhAliSsXdiIsrCOM49uxsIZ1WjvGnSULObqtjIWGLdOLHm0blw3D5u2dFEZ6txfMoTCTiKDkeWx65rs4G2rivlgeAk7OR1kt7ATv43XGAgcmVhCJyaGvaDXa65Mpl1nZC/ufIaL1mtyn5Nrt2RmHZXkB5yZ2AxjTGT4TGyerEtFNOAQbMgFsvAqOa7bp46ahEdf2k7rLlhFYLbGho+uz5hYjtryEBUW4MHTsCtKAkhnsvivxzc7xpYLhbAo/+nEGJbPqRcKQNiZWG8nVlVkXH/OvPwHsuBFMeNx9dnzcM7RzZ5/9wLJnqt8JtaDlkVw/nHT8ZvH4gjo4prY4TqxkiXmlc4YwxZ24o8NFC7sVF4SwPevXW6q8+YVScv/gBy7qKmgWtdPCrnOccmsWsybVjUssRsfhyYk2OubMxNrPiNHzqyDpsqU4RQKqOjpH3AFQETiSN+45Ejc+8dNLtEU4uxkGQFDXZUxmDCzRbnW23JrHRIxh0R0YtagnsrU+Gmq4qn0r1rzD2DSF0krtAB1YkdGs//Hxt1IpDJ5HUzAOdd40Yk9BdgkMd2XBb1nglOSJAllER179w9hskDbArCvRSigoHfAvK980KIiGsAcxgnIxxABgBvXzXf8zorRhaya2POOnY4ls+wAgFeZjVikj2WsKJAlCd+9ehm27ekVsrMITjhyIu770ztUsLLQ4H+hkCQJn1o6ibax4cc8XHaMIks4dlETlsy2r1O9Jc7E1mkW44zzIM+LG3YGHQBlIPztlR2YPcXpHLLnV10ewoUnzIBhGDhn1VSs4ILLi2M1eG9bA56zsq1eQqfnrJ6Wc80mbabKSgLoH0q7tr32zDkuRmIuSJLNFmUd68bqCGrKQzT7+283HA3AnIvKrADRUa31VPxyelM5pjeV01ZQxy+eiHVrphU8jtHAuLUAyGJYFgnk7clJkGXoxOyDLK45ZLn8imMbTTWpQQOJtLDPIZu1lT0ysbyBVyidmNB9iC1AI6esUFWOfZHo1CnLJlNHihd2IpMmnyWkLTz487F+lbieZbx7pXpEcwm+/rkjsXe/WPqdxxc+PQ8vbWoTFs0TrFvdgnOOmUavtacTKzmdoVmTK3DH1cs890vOgjxDxJF43pKNL8QAIWPK1aw7oCmYNcUdnTPPwcrE5j1S4dAci+7o0zxpJjaHsJMoM33M/Ak4RtCyiBhWhdTvur4ry0hnMkVnAh3zwTAW58o87Shsmnz+fR3Ptb0YKYbDDr/twkUuCufwjiUVReX2ceiBfVdYJ9YwTDrxSUsnObYPBeyMKQu+JhYwW5tEwxra9jk+xp9e2Irq8pDZSo4RfQFMwzUXbr1gId78oFOsZs+tfQpXcsKrcpO5IZ119qQOBVT0DqRw9LwGXHLyLOrEsmrGI8HvnjT7wxbaEYCA2DN8YM7LpjAMsQPHgpxTJiNerQhjatkccedO8swENAW9SAm3qS0POYR18lHFRYhy6s2SJFGKL4EXQ0nEIOLVXwGTcupFOyVYMqsOS2bZGeJCg6TDAcsYYEHekSNiNZg2QVwKwkPixLoAs6f4u1u7MGtKBf7w3BbHvscSLJV7cn00Ly3WbnEpCRlrmqrgkpNnYcmsOgwl3WrNDVVme5t8JTyaKiORyqAsouPjvf2utXYk95j1ZY6e1+AQsswXYCIgtn1laWDEc89wMY6tALtesBDqHUFG4IjlMxjNVjP272wLHnFNLPtvcSaW/yxXYoadIInhztNpC83EBnWVvtik2bE7K2z+5Os1vepQifKkhNyZwXwRPlLfUAhqK8I4nRHX8QI7VrtXrri2suBCdkInJplYLhtWyCRAvptLXTgXiIF3RgHXoFCwE2k+w6QYZC1Jf2LYlUZ0NFSFHRH44dTZkmtYjBOrKhISqeG3CRJhNI0MW2F4ZJmYsQYr7OHDRy5IEsMA4ujEkmCtInRiVuGalKqI6MQiNs77O/fjrofewMzJ5XQ9q6+MoKN7KG9f09qKMI4/Utz+hF8r+RIaVy2p9Xs6nXU4xWHLieUd5cAI6cRhpt0bkLtcRQQS0C40wGQYRl7nhNwzL9vg3DUteOHtPbTGkwf5Xi5hrPrKcM77UAjYAITXsbzui1hVf3Tm8LFwKry0PoiNdt1Zc4V/LxS6puBWq30OdWLHwBknMO9dP71vhQoGFsJUAOD5bN7+uSML+j5xYhfPrMW727qwekEjFTMDRu/asD2ZhwPCCim2vGokGMdOLIMC5woDBqW4shOD12RDKAABLhOrKjIt1OezmIDTcRK12DH3zZ0CW0cLp0qfs4+bM4NHjuWgRxf4AoscYMAU9gHcL5Z9TH7s5k8DZsQ5H/jo5oEC3yuXgGb1CxA/AphaYkXsxBaygBLjpVhnRZKkUS++N4W5TJpboRG84SDLZWIDmoJvX3GUa7toWMNyj6i8CHydUiGgVORhiliIMJpGBhlXISIwPnwcCjBr0JlFwkLWEAdvwwEVQ8mMY/0hc20h2VGCRCqDbNZeW5fOrsVbH3UW3WNedCyy7y+fv1DY1oKqDGeyDkeHOOj8e27XfzqPc+sFC/Hd376ed3yT66O09Rfg7awAJhV7kOsFSoLkhTJUDMPuZ+oFcs+8WEcnLZ3kysazIN/LNScSddma8iA6uoeKomOTGljS1scLqxc2YgbXm34sgr4EY5GJLSbwO1KMZYbvytNbsf7dNqQzWbz0TpswOCZCvlZ2eb9f4H0n93BSbQnuv22tK6CTS+BzOBC1mywEhKGarzXeWGDcOrHsOlSwH2DYTpaSpyYWMPs2dfclrcbJ9ueaKudUEeWFnQqqiWV+//GXjnG8HCzNhdCGKbVXIHpTqBNrCxQ5JxdCE+QnT5r99brghrsOlseBVD3jQRZTTZXx+bPnUkl0onzIKtQVAnIPeUOhkMn6jJXNUBRpWM7agYCmKkhn0mOyKPPCTl74oVXHkQ9DlgFWFJ2YGGsHWSZ2WWs99uwbwKkFinGNBk5aOhnx7d2Ym6NOy4ePYiExr4cBs8VcWUSHYakT8yDzNOn5CdjrnjAT67HeybJkqvBbx1gyqw67Ovo9RXWWzKpFR7e7d2YhxyK0XZ5ZQ2itqXSWC0Zbgk+cs2Ur4zvnFK+WGSy+8Ol56BtMUSe2pjyIs4+Z6rm9yOAdrs5AZWkAOtMSZiDhplyS4HGxrCOaic2xJhF7aVlrPf7f81uLWr+IExsKKDmzqBefGPM8/lhgNGpJebDP1w2fnodZkyqQyRZfHlIIRiszLUJZRMfxiyficUsbp9DkQKG28kjB63/wxx2tTGyxJTrEiS1WI2QkGL9OrPVTYqO8FljRBP47InViL3zh0/Px1/XbUVEacNXAkZudX51YEr5Q7oiu/W/NagpOwDpJxPluba7A4pm1WLd6Gj2O6Pi54KVOTKTPVy1wLvYiKjbAKB7DyFnj+UmDLKKaKjtEm1oayxDQFZx7XAFtUuAUsgLcBkchUeBQQMW61S0FHe9Agqomj0H9Cu3tO0oT9qDVLDyXkrIXbGNt5MbHaDqxqiIf8OeisTqC712z/IAe08f4gdnCzvz35h3d+NHv38LVZ7TSPrE8SDaAXUvIuidUDPZYyxVZoi12yHfXrfF+t64+Y07+c8ljN/BrL9GgIGJuuiYjmcrS4DN//sSe4M9TpL0BmMb7NWfOgSxJaGkqo7W1gEl1DAeHF+Aj489FK7zpvAXY1dGPaEjDwhnV2LbHVNMtK9GFTiyxlfIFuL1AnoMjYrXYtLULc6dWUSEaHmesbMbS2XWulkKFgNCJC+nny2MsM7GjLezEo6EqbGW5D332D7lW+Wzg799wNF59e/eBGBIAoKmmBO1dg457ee+XVuGaHzwDYOQ2xPevXT4ijYp0gQmGscC4dWKJFyvBnYlVFFO0xfUVw6DRJk+1PQaT66O46vRW8ziMwyfLtihJPtEmrzWPd35zOaGsE6tS4QUF156Zf9HNBeqEcWOpKQ8JM6arFzTitc0dWMn1JKNjN3ILFX3SmFhbgiNm1LjqSEsjOu790irU1ETR0eHuncZDZHj9nyuPwr/8/CUAY0ubGWsQqt0nmYktFIOjkIkdjRqQQ/l++/Ax1mDnyx3tfQCA+PZus8WO4NURORH8O8aqu3oFpE0nVuwoFwuytpdFdHznKncpBI/FsRpcf/ZczG8x1fyDmoJkKkv3ww+dii1a63x5iY7uvqRwjrn4xBjmTq2ijjLgvBbFzEspkpHJwVBpnVKJVkZwkIgkVpeFsLtzwLU9FXYq0jYgweeZkytw/21r8bdXdlAn9nvXLHOVhhXjwJJxqopUVDaK3C++RctoYKwdi7F0wAFTJGpDvCP/hqMAcq3yvfKxyZWoDB84SvVlp8zCstY6R+9hlh4/UlpzPsHIfEgL2kMeKIxbJ9ZgvFj+gVVlCV6NWkgmVlTLmgs0mksUVXMUkPM1sSK4Fy/v77A1f6Mp+GLTiQvbZ1VZUFjDSGBgbIv3RwpVkXHd2SMTLACcNcAE9UwbnEPZqSF2xljQozKjnIkdsjKxRQk7EcXqUXFiDwwlyYePQxFmoNnJ+klnssh60IlFWQl2Prr3plUFaUAosoRs1ii83KgAkPHLefq+kt7kkiTRPo+A6Rz2DKQYJ9b8eeXps7Hpo310TSFryLevOAqpdFZ4jktm1boyrez1LGaeTabMOVVUG/fVixcLP589pQInLZmE44+ciJt+/Lzr7+Q6FUsnzlDHWnH8BEzHebQgSRKiYb2oEhNJkvC1zy52KQ/fuG5+weJBXiDvw1jVsY61E3vdWXNHVIc+HJD35mATRgwFVBwRc/dqlmDakdoBUG7OhYtOjOHhpz7AxNrRD8Lkw6FrLY8SJOa/BLwTQZTWsgaQppnY4T3k5KUgkR6aiRXsh1cnFu4vB53Y69jmdqMYVfagEw8X1KkzzGbSZxXRQ/RQBD8x861jDkXYmdixoxMfHJlYU3BtNIIuY1nr48PHoQ6z5McE21osa4gdUFEQMMpkTQKa4mAyeWZiFRmGxzGKRSHG8c+/vNqzPzevnkrGdtTselx26my75IKxM0ojusNpn1RntvMQrTNOFtjwz5tkVUW1dVMnlDqCtQSKLOPctS0o91D3z9diJx+OXWwKQZKyEVLqMhazbmU04HjWhoPmhlLXWjRvWhVik4bX5oiHqsi48PgZ+JeLjhjRfrww1k4scODWSHKYQ2VJpi2tPuHER3NDKW65YFFRralGivGbiWXmQ36N4jMjrc0V9Et2TezwHhryUtAC7QCJLrq3dWZixftzt7U58G8dUXAbqRNrw0BpRMdpK5rxv89tyamMeCiDTMh8cFFXZaTS2UM6E8v25Btt8H1iR2t/kWHWfQHm+x/Ucwt4+PDhY+SQJbjakmWyhifVV8SSiuYIVHkGiq3uArxmxkigFEBXzDX/Ewpha3MlNn7YiakTSh1/J/Mvb8Ow+7zi1NkIBVShwTlSO4JkYovRCvCaSwO6gtlTKnDCkd4KxLlw2vIpuOT0Odi716SiBzXL7B2DqfuK01uHzdI7EDj2iKZR3+cXPj0P699tP2DiRgcCXm0gD1bcuG4+ntiwc1x3Ixi3Tiyj7OT6E1/vShYAA4w4kSCDetuFizwdLzJBk8WF1NKRSZ9FIZnT4WRixwpkCIXUB+cCqclho7ffueqoopyLQwFet0rXFPQPpQ9peinJxB4KdOLmhii27O4tSoRKUaSilfx8+PBROESZ2JffaTN/F0yVovkhVxutfMJOo2nPjjTgS6i1R86sxanHtCAx4Cx8IvOvq5c5c9yArnjWwA0zNu9CKkcmtljIkoSbz1s4on2wdhQx+EczOEFQWz569OSDHfNbqmmt9uECEgQ62OjEXpg9pRKzp4h70I4XjFsrjNTEsvU2BJ5ZTsOWklZlGbdduAgf7NpPt5sxsdzzeHYm1lmXIXRiWUqPtaisWz0Nuqbgwb9vdo5J8J0DBVvYaWTHXrd6GqZNKKNtBgA4CtgPN5Bb526ncOhmYAmyYyjsRDBa1JmbPrMAXb2JorKpqix9Io29ffgYb5Ak0EmzEAYSm3UM6AoSyUzOkgGvtVOWrUzsaNKJZVEBU+EI0j6wMkojOjo4J5Ywu3ixRb5Pfb7xFQtCJy52bpQlCScsGds+8MWoB/sYHyBBoEPEh/WB8ezEEl0nSSDsxE/yjBAPm4mdMbE8p+PKQuYysSRSKZK1dmRirYN/6qjJaOuylfu8mqZ74V8vX4r2rsGCxlooKJ14hE6FpipYOrtuNIZ0SICIN/DZgdbmSrR37TogNSZjBfJejUWLHYLRcvbDQW3YLSQI5rdUCwNQPnz4GF04MrH8OidY9lhHl/wrV2bQKzuazmRHPRM70pkrQJ1Y8aBIcD5XxjenEztqdOLiTMtf3LpmRMcvBGR99R0VHzwCFtX8cC1lOxwxbp1YFvxcxi8Q5DfDMKjS3XApn3xNLJnkEymBE8usMexEyy4w7oh07uNPqI5gQvXoKoeRYY5eTez4wFnHTMXUCaWYPdkp2HD+sdOxdmHjiOXOP0kcSpnYkeDEJcXVZ/nw4WN4kGVv9kpPf9K1PeukkSBxruybV/Yxlc6OeosdG8Xtk83EikCEnXINOdd6PdK1PJkafTrxaONQDhL7GFssnFGNi06YgRVzGz7pofgoEAfvTDPGcCyFeai57CKWtoSdhhuxdKsTW3TidO6aWM9/HwTCTiVhDVWlATRUHb7U37GAqshCuXRVkdFYU/IJjGj0QDOxY6hSp6qHR9DkB59fMaIG4z58/RohqAAAH7NJREFUjAdIkkTXYL5X6L5edzM8NsBMtg/mcFwUZu1cOa8B/9i4G4DpxGaN0W2xM9JGIaVhHbomezqbtKavyEzsSB32686eg7+/suOgbpVHa2IPj2XExyhCliSsWTT6Ilg+xg7j1okF5b5LOXuuEkgw602yWQOKLA17snf1ic0h7MTuW5LFjis/xk9iQg7qKu68dsWBP7CPgx5jqZY3XGXwgxXlJYFPegg+fBz0kGHnLbOcE9vVM+TaXuSk5ZqPyLqqKhIuPXkW48RmLDrxKKoTW8cqtg3LsUc0YcH0ak/7I8vYNZ5jyMEiG2lN7MLpNVg4vSb/hp8gSDZ7SkNpni19+PBxsGPcOrFUnFjwN2GUUzK/lUxniqr3k/hMrE6c2NyZGEfPWOaXg6HFjg8fXhiLSPySWbVY/277qO/Xhw8fBy8kSaILNZ+JFSnYs05sbUUI7V2DVIdABL7nKsFQMoO+wdSoUmMrS4O46MQYFhSp6hoKqGjKwdahWh859pHLVhgPpUGqYopyNtaMbnmVDx8+DjzGrxPLeLHurKYoE2vK7afS2aIa+rprYs195Csglxx0Yvvz4Qo7+fBxILA4VoNX4x1jorR85WmtuOyUWaO+Xx8+fBy8YHxYRyZ2cn0UN5wzz7U9Syc++5ipqIwG0Zwj60bWUuLAXX1GK/7w3Bbs2TeAHe19mDmpMPHGQrFmYeOo7o+FUUAmNhfGSzC8UEFOHz58HNwYt04sgUT/Y0NEqSFzezKVLSrLxGdiNVXBVae3YnpTWZ7vicflqokdhQjqty5firRfo+djBLjitFZceHxqTIIqsixBln1RDh8+Dnfc8Ol5+Okf3kYynbUyse6a2KWz6lARdVPyWbV8XVXQkmeNJY4b+blkVh12dw7gj//YAgDoG0yP7GQOIOya2OK+/0m06vPhw4ePYjGOnVh7MeSbXnvZ32YmNlNUlknmMrEACmor48zEejuxo+EzNI6yerGP8QdNlVHm13r68OFjBFjQUo3Gmgi27O6FzLTYYTOxAY+yHjYTW0gXAZKBZddadp3u7Bnd1nRjiUJqYnPB92F9+PBxKOHwUEgpAjadWHI5gEJhJ8nswZZMZ4tSXiWLSi5lQBE8a2ILoED78OHDhw8fhzJYOnGGabGjeygO52pFJ9xesA3rxA4mDp1+0BNrzXrZiiIDieRajBdasQ8fPg5tjNtMLCuAwE/XYidWAkhN7EiEnYaZxS28xc6wh+TDhw8fPnwclGB1K8gizWZivZxYNqCrFBA0jgRNM2gwYdOGx6Km/0DgzKObsWB6NSbXR4v6vuwhcuXDhw8fByMOzZl6VGDRbgAXF1cUhJRgLqrJdKbImljz53AzsQ46MfNVX9jJhw8fPnwc7pCs/wHOmlgvOjGLQjKxpNUVu+9MZqQdXT8ZKLKMaRNy1wDnAgmUjweVYh8+fBz6GLeZWBZedOK1ixqxbnWLtZFFJ05lEQ3lVhQWQUZxmViHsFMOmpRP//Hhw4cPH4cbGF0nZya2gLKeXD1RCUQ1/ANDKfrvFXPqCxjl4QGvdkM+fPjwcTBi3DqxLFWJn64l2XY4SZN00mInmc4WRTXiW+wUCtY59RJ5Yvfvw4cPHz58HI5wCDvp+Z1YtQCZ3vISd1CaHOXC42fg2COaCh7fwYxvXb4UiWTu+l7Jz8T68OHjEML4dWKtnxIkFxWXzN8GyyiyPksVTScuTtjJyzl1ZWL9RceHDx8+fByGIOtn2pGJLYBOXEAmVtSr/cQjJwEGsGrBhGGM8uBGId0HiBnhmxM+fPg4FDB+nVgi7CR504mzBtuGx6qJTWU9BSVygQSEh08nFq8mki/s5MOHDx8+xgFsOrHdx7yQdbiQjKIosBzQFZy+srnwAR4myNI+s75B4cOHj4Mf41jYyYZruiYfMJlY0mInVTSd2KIoj6DFDgt+cfaFnXz48OHDx+EIcZ/YQpxY38QZDoIWRXvF3IZPeCQ+fPjwkR/jNhPLeqhuOrH5O0sntmtiM9CLarFj/hytTKxbndj504cPHz58+DhU4dAHltzqxIWU5hRCJwaAGRPLHcyr8YpQQMW9X1pVVBtBHz58+DjQGLdOrE0ndis7UTqxw9EFMpksDAPQClBF5CEXWRPrpTqsCBzvmZPKcfziicMemw8fPnz48HGwgs/Enrp8MkKB0aETA8BtFy4qdmiHHQoRzPLhw4ePgwHj14ll/s07ipSBxAVmEymzHmckfWJH0mKHhahP7C0X+AuxDx8+fPg4zGAtdyQTe/Yx0wr62nCDxj58+PDh49DB+J3hrVSsyEmUKJ3YSTl+cdMeAEU6sUX3iS2MTuzDhw8fPnwcNnCU85jIGsawSmb8VjE+fPjwcfhi/DqxFiS4hZ1sdWL7M1ZQohg6MVl4x6rFjg8fPnz48HG4wewgYK3JWe8SGxEKrYn14cOHDx+HHsatE+toAcstipPqSgAAsUnl9LMMk5UtRtiJLLzDzcR6ZVyHs5D78OHDhw8fhzoyWWNYLCR/nfThw4ePwxfjtyaW+KSS5Mp2NjeU4p7rVzqaoDszscXXxKrDjAx7be13DvDhw4cPH4cryPqrKjLTJ9YY1vrrt57z4cOHj8MX49aJpTWxgj9JkuRwYAGnE5vODF+KXyoyE+u1CPt0Yh8+fPjwcbjiitNmY0O8HU01JRgYSgEwa2L9AK4PHz58+AB8OjEgFUY5YvvThYPD9/1tdeLRka/3hZ18+PDhw8fhipKQhlULGq3f7PXOpwj78OHDhw9gBJnYWCw2CcB/AqgFEAdwYTwe7+O20QH8EsBiAIMALojH4+/FYrESAPcDmAlzdfp2PB5/qNixjAQS/Q/zWY418rjFTWidUjns49Ca2FESmvAXch8+fPjwMR7ALneFUISXtdbhxU1tYzgiHz58+PDxSWMkmdifAPhJPB6fCeBVAF8TbHMDgP54PD4LwI0AfmV9fhuA7fF4fB6AYwH8IBaL1Y1gLMOGwcr3805sju9NqY8WdbzW5kqsWdSIUGB0GNx+rY8PHz58+BgPYJe7QlhIV5zWivtvWzuGI/Lhw4cPH580inJiY7GYBuAYAP/X+uhXANYJNj0FwIMAEI/HnwVQY2VwnwHwb9bn7QD2AagvZizFwgDpEyvRHq6FQCmyIGdSXRQXnRDznU8fPnz48OFjGJAcdOJPcCA+fPjw4eOgQbFpwWoAPfF4PG39vhtAk2C7CdbfwG4Xj8f/Tj6IxWLnAggA2FTkWIpDjkxsLj6xL6jkw4cPHz58HEAMMxPrw4cPHz4Of+R1YmOx2DoAd3Mfvw9nq1UAyAq+LnPbSex21r5/COAkxiEuCFVVJcPZ3AXj3XYAQE11Ccq6Bh1/m9xUjopoUPi9yoowamqKoxQXA69jHQxj8OGGf60Kg3+dCoN/nQqDf50Ob7Buq68H4cOHDx8+gAKc2Hg8/jCAh9nPLDpxZywWU+LxeAZAA4CPBV/faf3tQ+v3erJdLBa7HsCXAZwQj8ffGu7AOzv7HG1vikVnZx96e4YAmAnY7161DOmhFDosSX8e/f0JdHT0jvi4hcLrWAdqDDU10QN6vocy/GtVGPzrVBj861QYRuM6ybI04sCoj7GDoybWd2J9+PDhwweKrImNx+MpAM8B+Iz10cUAHhVs+hfrb4jFYisBDMXj8e2xWOxMAF8EsKIYB3Y0YAj8X0WWUF0eyvk9n07sw4cPHz58HEjY667kr8E+fPjw4QMjaLED4FoA/xGLxb4KYDuA8wEgFotdDWBCPB7/OoB/B/CzWCy2Cf9/e3cfZElVHmD8uTO7CMIqAoN8qRUivooKq4LrBygloIWIwYCgGAEjgkFLjZpES9T4RUApBRVE0WQxgFqSmICgIqBhiRUJiRoD8hYSRRcwbi0aWBDcZW7+6L4zvcPsbs/Hnbl9+/lV7dbtPt33nnv2zLz79jl9Gh4EXlue+wFgG+DyiOi930mZeeMc6jNDvSy2M7HY0nSJ7VQmsZIkLZyNR2IXrx6SpMEx6yQ2M28HDppm//mV1w8AJ0xzzL6z/dz50ktYO53NPxd2qtHRuTyVSJIkzZYLO0mSYG7PiR0avSR2kEZin/oHOyzI50iSNMi8J1aSNNVcphM3WjVhndlzYhcmgP75K/dlvE5WLUnSENv4ObEmsZKkFiexvXtiOx0m1ozoPuypQQ+3UNOJR0Y6jMwguZYkaSht9JzYxauGJGlwtDYcTNwTS2dyoYgaA59LFvl+nJ23sHqyJEnDxOfESpKmau1I7ES+2pn4q5bFXp34QyetoOs0Y0lSS3QqiWvHJFaSRIuT2J4Ok0v210kNF3t14qVLWjt4LklqqQ5FjHZxYkkStHo68UYrO9W22COxkiS1Thl6fcSOJAkciS2eEzuDLHbJqAFUkjS4IuI44DRgKXB2Zp47pfww4Mxy88fAKZm5LiJeCPwj8Muy7AeZ+boFqvZmdejQpet0YkkS0OqR2N6rDjOJiaMujShJGlARsTvwEeAAYDlwckTsXSnfHrgQeFVm7gP8CDi9LN4POCszl5d/BiKBhclnxToSK0mCNiexlUfszCQkjjoSK0kaXIcA12bm3Zl5H3ApcHSlfC/g9sy8udz+OnBk+Xp/4MUR8V8RcVlEPG7Bal2TqxNLkqDFSexGZhAUvSdWkjTAdgPuqmzfBexR2b4VeFxE7FtuHwPsUr7+LfCpcoT2SuDLfa5rbRMjsYZgSRJtvie2sq7TTIKi9+NIkgbYCBsvtt8BxnsbmfnbiDge+FxEjAAXAL8vy95YOe78iDgjIh6dmf9X98N33HG7udafsbFlD9tXxN4uW2+9dNryNrId6rOt6rGd6rGd6ul3O7U2ie1FeHNSSdIQWQ0cWNneBbiztxERo8DqzFxRbu8P3FYmtO8GzsjMhyrnb5jJh69du47x8dk/y3xsbBlr1ty7yfIN6x/abHlbbKmdNMm2qsd2qsd2qmc+2mlkpLPZC6OtnU7cW9ipQ8fRVUnSsLgaODgixiLikcBRwDcr5V3gqojYPSI6wNuBr2TmOPCK8njK0drvl/fVLrpelO44n1iSRIuT2EoW62isJGkoZOYdwHuA7wA/BC7JzBsi4sqI2K9MVk+hSGwT+A3wsfL0E4C3RcRNwOuAkxb8C2yK98RKkipaO524p8PMnhMrSdIgy8xLgEum7Htp5fUVwBXTnHcT8Ly+V3AWenHaR+xIkqDFI7EbrXphTJQkaXBNjMQasCVJbU5ie7OJZ/qgWEmStKBGTGIlSRXtTWIrY7FOJ5YkaZD1phMvcjUkSQOhveGgW1nt0BxWkqSB1QvTjsRKkqDNSSxMREVjoiRJg6sXp33EjiQJWpzEbrSwk9OJJUkaeI7ESpKgzUlstzuRvBoTJUkaXJ0yUJvESpKg5c+JnUks3H2nbTcavZUkSQvLHFaSBC1PYnvqXNn90EkrFqAmkiRpql6YHvWeWEkSrZ5OXNkwJkqSNLAmVic2iZUk0eoktju52qHzkyRJGlxlnHYkVpIELU5iC53K35IkaSCV06dc2EmSBC1OYrvdynPnjImSJA2s8fIWIKcTS5KgxUksTI7AOp1YkqTBNV5msU4nliRBi5NY13WSJKkZumXU9qKzJAnanMR2u5Wh2EWtiiRJ2ozedGJHYiVJ0OIkFqDTW9jJK7uSJA2sbpnFek+sJAlanMR2u0yMwBoSJUkaXC7sJEmqam0SC9WFnRa1GpIkaTO6XRd2kiRNam0S260s7dRxLFaSpIHVW53Y58RKkqDFSSw+J1aSpEboXXYeae//WiRJFa0NB+Vi/cXfZrGSJA0874mVJEGbk9hu10nEkiQ1yKhDsZIkWpzEgtOJJUlqEgdiJUnQ5iR2cl0npxNLktQATieWJEGLk9hKDuu0YkmSGsBH7EiSoM1JbLc7MQLrQKwkSYPPR+xIkqDNSWzltdOJJUkafE4nliRBi5NYSZLULE4nliRBm5PY7uQ0YmOiJEmDz5lTkiRocRLbpbqgk0FRkqRB50isJAnanMR2J4divbArSdLg855YSRLAktmeGBGPBy4CdgYSeE1mrptyzFbAF4D9gN8Bx2XmLZXyJcAq4LOZuXK2dZmtXig0iZUkafC5OrEkCeY2EnsecF5mPhm4EXjvNMe8BbgvM58CvA1YOaX8fcCT5lCHedFxOrEkSQPP6cSSJJhlEhsRS4EXAJeWu1YCr5zm0MOBiwEy8zpgrBzBJSKeB+wLXD6bOsxVt3JTrBd2JUkafE4nliTB7EdidwLuycwN5fZdwB7THLdbWUb1uIh4FPAJ4ORZfv6cbf2IUbbeqphNbRIrSdLgM4mVJEGNe2Ij4pUUCWfVrRQL/FaNT3P6yJTjOuVx5wKnZ+b/RkT92lbsuON2szqv51WHBi97/p6M7bQt6zc8NLF/bGzZnN53GNkm9dlW9dhO9dhO9dhO7WESK0mCGklsZn4V+Gp1XzmdeG1EjGbmQ8CuwJ3TnL66LLut3N6FYjT2YODpEfEB4PHAiyJifWZeXLfia9euY3x8ah5d39jYMh6470HWrLmXDQ9N5t9r1tw76/ccRmNjy2yTmmyremynemyneuajnUZGOnO+MKqFMerUKUkSs1ydODPXR8Qq4FjgEuB44BvTHHplWXZ9RBwAPJCZt1NMMwYgIlYC351JAjvfjImSJA0+R2IlSTC31YlPBU6OiJuBA4HTACLijRHxwfKYTwGPiIibgE8Cr51LZfulYxYrSdLAM4mVJMEcnhNbjqgeNM3+8yuvHwBO2ML7nDjbOswXQ6IkSYPPR+xIkmBuI7FDw5FYSZIGnzmsJAlMYiVJUkM4nViSBCaxkiSpIZxOLEkCk1hJktQQ3v4jSQKTWEmS1BBOJ5YkgUmsJElqiBFHYiVJmMRKkiRJkhrEJFaSJEmS1BgmsZIkSZKkxjCJlSRJkiQ1hkmsJEmSJKkxTGIlSZIkSY1hEitJkiRJagyTWEmSJElSYyxZ7ApIkqT5ExHHAacBS4GzM/PcKeWHAWeWmz8GTsnMdRGxPXAxsCewBjgmM3+1cDWXJKkeR2IlSRoSEbE78BHgAGA5cHJE7F0p3x64EHhVZu4D/Ag4vSz+MLAqM58CXACcs5B1lySpLpNYSZKGxyHAtZl5d2beB1wKHF0p3wu4PTNvLre/DhxZvj6cYiQW4EvAYRGxdAHqLEnSjJjESpI0PHYD7qps3wXsUdm+FXhcROxbbh8D7DL13MzcANwDjPW1tpIkzYL3xFY8/rHbLXYVJEmaixGgW9nuAOO9jcz8bUQcD3wuIkYopg3/vnIsmzq3jh13nHscHRtb9rB9j9x6Cfc/sGHasrayLeqzreqxneqxnerpdzuZxJY++dYD2WqJA9OSpEZbDRxY2d4FuLO3ERGjwOrMXFFu7w/cVhbfUR6/OiKWAMuAtTP58LVr1zE+3t3ygZswNraMNWvufdj+s059HuPj3WnL2mhT7aSHs63qsZ3qsZ3qmY92GhnpbPbCqElsabttvO1HktR4VwN/HRFjwH3AUcDJlfIucFVErKBIbt8OfKUsuxI4nmKhp2MpFnlav1AV35ytt/K/K5KkSQ49SpI0JDLzDuA9wHeAHwKXZOYNEXFlROyXmePAKcA3gQR+A3ysPP29wHMi4ibgVOBNC/4FJEmqwUubkiQNkcy8BLhkyr6XVl5fAVwxzXl3Ay/vewUlSZojR2IlSZIkSY1hEitJkiRJagyTWEmSJElSY5jESpIkSZIawyRWkiRJktQYJrGSJEmSpMYwiZUkSZIkNYZJrCRJkiSpMUxiJUmSJEmNsWSxKzALowAjI505v9F8vEcb2E712Vb12E712E71zLWdKuePzrky7WVsXmC2U322VT22Uz22Uz39js2dbrc7pw9YBAcAqxa7EpKkoXMgcP1iV6KhjM2SpH6YNjY3MYl9BLA/cBfw0CLXRZLUfKPArsC/Aw8ucl2aytgsSZpPm43NTUxiJUmSJEkt5cJOkiRJkqTGMImVJEmSJDWGSawkSZIkqTFMYiVJkiRJjWESK0mSJElqDJNYSZIkSVJjmMRKkiRJkhpjyWJXYDFExHHAacBS4OzMPHeRq9R3EfF+4Jhy84rM/MuIOAT4OLAN8JXMPK08djnweeBRwHXAGzNzQ0Q8HrgI2BlI4DWZuS4itgcuBvYE1gDHZOavFvDrzbuIOAvYKTNPnK/2iIitgC8A+wG/A47LzFsW/MvNk4g4Ang/sC1wVWa+1T71cBHxJ8C7y81vZOY77VOTIuJRwPeAl2Xmz/vdh4ahzYaVsdnYvCXG5i0zNtdjbN68JsTm1o3ERsTuwEeAA4DlwMkRsffi1qq/yo73YuAZFN/5WRHxauBvgT8CngLsHxGHladcBLw5M58EdIA3lPvPA87LzCcDNwLvLfd/GFiVmU8BLgDO6f+36p+IOBg4obJrvtrjLcB95f63ASv7+T36KSL2BM4HjgT2AZ5Z9h/7VEVEPBL4JPBCYF/gwPLn0T4FRMQK4HrgSeX2NvS/DzW6zYaVsdnYvCXG5i0zNtdjbN68psTm1iWxwCHAtZl5d2beB1wKHL3Ideq3u4B3ZObvM3M98BOKjnlrZv4sMzdQdMJXRsQTgG0y89/Kc1eW+5cCL6Bor4n95evDKa6qAHwJOKw8vnEiYgeK/0idXm7PZ3tM7M/M64Cx8kpVE72C4krc6rJPHQvcj31qqlGK37PbUowuLQXWY5/qeQPwJuDOcvvZ9L8PNb3NhpWx2di8Scbm2ozN9RibN68RsbmNSexuFIGj5y5gj0Wqy4LIzJt6HSwi9qKYujTO9O2wqfbZCbin7LzV/VTPKcvvAcb68mX677PAe4DflNvz2R7D1PeeCIxGxGUR8UPgVDb9/VrbpzLzXoqrj7cAq4GfA7/HPgVAZp6UmasquxaiDzW6zYZY6/5djM0zYmyux9hcg7F585oSm9uYxI4A3cp2hyJoDL2IeCrwbeAvgP9h+nbYVPtM3Q+T7daZsr+RbRoRJwG/zMxrKrvnsz2Gqe8toRg5eT3wXGAFxf0N9qmKiNgH+FPgCRS/oB+imD5on5pe3b5imw2f1v67GJs3z9g8I8bmGozNMzaQsbmNSexqYNfK9i5MDpcPrYh4PnAN8K7MvJBNt8Om9v8aeHREjJb7d2Wy3e4ojyMilgDLgLX9+SZ9dSzw4vLq5QeBlwMnMX/tMUx971fA1Zm5JjN/B3yNInDapzb2EuCazPx1Zj5IMaXmIOxTm7IQv5eGrc2GRSv/XYzNtRib6zM212NsnpmBjM1tTGKvBg6OiLHyxu6jgG8ucp36KiIeB/wTxUpfXy53f78oiieWnew4itXZbgceKAMrwGvL/euBVRTBBOB44Bvl6yvLbcryVeXxjZKZh2bm0zJzOfA+4LLMfB3z1x4T+yPiAOCBzPxFv79Xn3wdeElEbF/2n8Mo7n2wT23sR8AhEbFtRHSAI4B/wT61KQvxe2nY2mxYGJsLxuYpjM0zYmyux9g8MwMZm1v3iJ3MvCMi3gN8B9gK+Hxm3rDI1eq3dwJbAx+PiN6+84ETgX8oy65k8gbs1wAXRLG89n9SrOAGxb0VF0bEacAvgFeX+98LrIyIm4DflucPk/lqj08Bny33P0jxw95Imfn9iPgoxep1Symmwn2G4v4S+1QpM6+KiGcA/0GxaMQNwBkUV8ftU1Nk5gMRcSL97UND1WbDwthsbJ4FfydMYWyux9g8M4Mamzvd7tQpy5IkSZIkDaY2TieWJEmSJDWUSawkSZIkqTFMYiVJkiRJjWESK0mSJElqDJNYSZIkSVJjmMRKQygiuhGx0zy8z34RcemWj5QkSZtjbJbmT+ueEyupvsy8ETh6seshSZIKxmbJJFbqq4jYDvg7YC9gnOLB2qeUxZ8AngMsAzrASZn5rxGxErgfeDrwWOAyYC1wBLBLedy15XG/A5YDOwNXAW/JzPVT6vB6iodOj5Tv8+bMvKVmPV8AfDoznxYR3yrrA7AdsCcQwO3AmcALgVHgB2U97pltu0mS1C/GZqn5nE4s9dcrgGWZuRzYv9y3J7AC2A14bmbuDVwIvKty3jOBF1EEqncA6zLzecA5U45bARwK7F3+OaVSRkS8EDgBODAznwF8FPjaDOo5ITNfUpavAO4A3p2Zt5b12QA8KzP3Be4Ezthy00iStCiMzVLDORIr9df1wOkR8V3g28DZmflT4KcRcRpwSkT8IXAQcG/lvMvLq7a/ioj7gG+W+28DdqgctzIz1wFExBeBI4FPV8oPB54IfC8ievseExE7ZObdW6pnROxR/TIRMQJcBPwkM88sd78M2B44tPyMrYBf12wfSZIWmrFZajhHYqU+ysyfUQSqvwEeBVwdEUdExOHAFeVh/wycTzFtqefBKW+1nultqLweAR6aUj4K/H1mLi+v1D4T2A/4TZ16TvN55wDbAm+a8hlvrXzGs/FeHUnSgDI2S81nEiv1UUT8GcX9LFdl5l8B36IIVodSXNH9DHAjxVXa0Vl8xLER8YiI2JpiatLlU8q/Bbw6InYtt98IXDODelaPeRfwXOCYzKwG5G8Bb46IrcqrwRdQBFxJkgaOsVlqPqcTS/31RYrpSDdHxP3AL4BPUizC8KWI+DHFz+FVwFFloJmJ+4FVwGOASymC3YTMvCoizgS+HRHjwD3AH2dmt2Y99wWIiN2A04EErqvU833Ah4CzKBaNGAV+SHGvkCRJg8jYLDVcp9ud+vMiqQnKFRD/OzPPWuy6SJIkY7O0UJxOLEmSJElqDEdiJUmSJEmN4UisJEmSJKkxTGIlSZIkSY1hEitJkiRJagyTWEmSJElSY5jESpIkSZIawyRWkiRJktQY/w/iNE/4hJgLcgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "nn = 200*np.arange(1,500)\n", "means = [np.mean(rd.normal(0.,1., n)) for n in nn]\n", "var2s = [np.var(rd.normal(0.,1., n)) for n in nn]\n", "\n", "plt.figure(figsize=(16, 6))\n", "plt.subplot(1,2,1)\n", "plt.plot(nn, means, 'b')\n", "plt.plot(nn, [0. for _ in nn], 'k')\n", "plt.title(\"Convergence of the empirical mean\")\n", "plt.xlabel('sample size')\n", "\n", "plt.subplot(1,2,2)\n", "plt.plot(nn, var2s, 'b')\n", "plt.plot(nn, [1. for _ in nn], 'k')\n", "plt.title(\"Convergence of the empirical variance\")\n", "plt.xlabel('sample size')\n", "plt.show()" ] } ], "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.7.3" }, "zanadu": { "authors": [ "Antoine Jacquier" ], "bibliography_data": {}, "category": "1", "clearance": "Private", "group_name": "ImperialMScMathFin", "md5": "040cba41cfebd4356c99f2c286825736", "notebook_id": "CA90D37D-032C-49B1-A80E-71AA78783EB0", "python_name": null, "reviewer_id": null, "status": "Approved" } }, "nbformat": 4, "nbformat_minor": 1 }