Skip to content

Instantly share code, notes, and snippets.

@jeffhussmann
Last active December 13, 2015 21:58
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save jeffhussmann/4980841 to your computer and use it in GitHub Desktop.
Save jeffhussmann/4980841 to your computer and use it in GitHub Desktop.
Class on Monday, Feb 18 2013.
{
"metadata": {
"name": "jhussmann_sum_of_uniforms"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": "####Sampling"
},
{
"cell_type": "markdown",
"metadata": {},
"source": "How well does the sum of 12 Uniform[0, 1] variates minus 6 approximate a Normal(0, 1) deviate?"
},
{
"cell_type": "code",
"collapsed": false,
"input": "import numpy as np\nimport matplotlib.pyplot as plt\n\ndef uniform_sums(size):\n def sum_of_12():\n return sum(np.random.uniform(size=12)) - 6\n \n return [sum_of_12() for i in range(size)]\n \nnum_draws = 100000\n\nuniforms = uniform_sums(size=num_draws)\nnormal = np.random.normal(size=num_draws)\n\nbins = np.linspace(-6, 6)\nplt.hist(uniforms, label='sums of 12', bins=bins, histtype='step')\nplt.hist(normal, label='normal', bins=bins, histtype='step')\nplt.legend()",
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "pyout",
"prompt_number": 1,
"text": "<matplotlib.legend.Legend at 0xb15a7cc>"
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAYEAAAD9CAYAAABazssqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XtYlGX+P/D3wAwiCogHGJ0Rp2AIRlFxE9HW/dEaiiRk\nqQiW4Gn7rpWr1W7mtush18DKXU2XbfMil/TbolutaGvkTxNzLSmPlZRMBspwSh1QUGQ43L8/yPmJ\nMyqHYR7geb+ua64L7vu5n/k84jzvec4KIYQAERHJkovUBRARkXQYAkREMsYQICKSMYYAEZGMMQSI\niGSMIUBEJGN3DIF58+bBz88PoaGh1jaz2YyoqCgEBQVh4sSJqKystPalpKRAr9cjODgYe/futbYf\nO3YMoaGh0Ov1WLx4sbW9trYWM2fOhF6vR0REBM6dO+fIZSMioru4YwjMnTsX2dnZzdpSU1MRFRWF\n/Px8TJgwAampqQCAvLw8bN++HXl5ecjOzsZTTz2FG5cgLFy4EOnp6TAajTAajdZ5pqeno1+/fjAa\njXj22WexdOnSjlhGIiK6jTuGwPjx4+Hj49OsbdeuXUhOTgYAJCcnY+fOnQCArKwsJCYmQqVSQafT\nITAwELm5uSgtLUVVVRXCw8MBAElJSdYxN89r2rRp2L9/v2OXjoiI7kjZ2gHl5eXw8/MDAPj5+aG8\nvBwAUFJSgoiICOt0Wq0WxcXFUKlU0Gq11naNRoPi4mIAQHFxMQYPHtxUiFIJb29vmM1m9O3b1zq9\nQqFow2IREVFLbgjRrgPDCoXCKStpIUS3fa1YsULyGrh8XD65LZsclq+lWh0Cfn5+KCsrAwCUlpbC\n19cXQNM3/KKiIut0JpMJWq0WGo0GJpPJpv3GmPPnzwMA6uvrcfny5WZbAURE1LFaHQJxcXHIyMgA\nAGRkZGDq1KnW9szMTFgsFhQUFMBoNCI8PBxqtRpeXl7Izc2FEAJbt27FI488YjOv9957DxMmTHDU\nchERUUuIO0hISBADBw4UKpVKaLVa8fbbb4tLly6JCRMmCL1eL6KiokRFRYV1+jVr1oiAgABx3333\niezsbGv70aNHxbBhw0RAQIBYtGiRtf369etixowZIjAwUIwZM0YUFBTY1HCXEru8AwcOSF1Ch+Ly\ndV3dedmE6P7L19J1p+KniTsthULRqv1bRETU8nVnq88OIpKz/x6pxYwlRyBg/8P10oLhWLSg6x7X\n6tu3LyoqKqQug1rBx8cHZrO5zeMZAkSt8J/vPsalh+ZhlHaYTd+3pefwr/NJWIRVElTmGBUVFdzy\n7mLae4YmQ4CoFRrRgP5Xx+PIr/9t0/fQy3/C9frrElRF1Ha8gRwRkYwxBIiIZIwhQORAFlSjrLrM\n5lVeXc597TJTU1OD2NhY9OnTBzNnzpS6nNviMQEiB/GGP464bMLINzNt+iqvVyJzeiamBk+VoLL2\n6dsX6MgThnx8gHac3NJpvffee/jxxx9hNpvh4mL7ffubb77B888/j+PHj+PSpUtobGy09lksFixc\nuBD79++H2WxGQEAAUlJSEB0d7fA6GQJEt6hrqMNbx96CpcFi0/flla9uOy7MJQkhliT86be2fYnv\nJ6KmrsaRZTpNRQXQkRsx3fUekefOnUNQUJDdAAAANzc3JCQk4Omnn7beeeGG+vp6+Pv749NPP4W/\nvz/+85//ID4+Hl9//TWGDBni0Dq5O4joFgWVBXjpk5dw/sp5m1cv1z7wr0iWukTZWrt2LbRaLby8\nvBAcHIwDBw4AAObMmYM//vGP1ulycnKsdygGAJ1Oh9dffx3Dhw+Hp6cn5s+fj/LyckyePBne3t6I\nioqyPiDr+vXreOKJJ9C/f3/4+PggPDwcP/74o916vv32W0RGRsLHxwfDhg3D7t27AQArVqzA6tWr\nsX37dnh6emLLli02Y4OCgjB37lwYDAabPg8PD6xYsQL+/v4AgIcffhj33HMPjh8/3sZ/udvjlgCR\nHb69fPGXSX+xaf/3NeCdyxIURDhz5gz++te/4ujRo1Cr1Th//jzq6+sB3P2OxgqFAh988AH279+P\nuro6hIWF4cSJE9iyZQuCg4MRExODN954A8uXL0dGRgauXLkCk8mEHj164OTJk+jZs6fNPOvq6hAb\nG4sFCxZg3759OHToEB555BEcPXoUq1atgouLC86ePYt33nmn3cteXl6O/Px8DB06tN3zuhW3BIio\nS3B1dUVtbS1Onz6Nuro6+Pv7495777X23+3A+6JFizBgwAAMGjQI48ePx9ixYzFixAj06NEDjz76\nKE6cOAGgaTfNpUuXYDQaoVAoEBYWBk9PT5v5HTlyBFevXsWLL74IpVKJBx98EFOmTME///lPaz2O\nOBmgrq4Ojz/+OObMmYOgoKB2z+9W3BIgukVdHXD5MrBpk23fqVPOr4eaBAYGYv369Vi5ciVOnz6N\nSZMm4c9//jMGDhzYovE3HoYFAD179mz2u7u7O6qrqwEAs2fPRlFRERISElBZWYknnngCa9asgVLZ\nfHVZUlLSbJcTAAwZMsT60CxHaGxsxOzZs+Hu7o5N9v5DOgC3BIhuYTIBFy8C331n++rRA3jiCfvj\n+vQBXnkFcHOzfe3YDnz1tXOXoztKTEzEoUOHcO7cOSgUCutzyXv16oVr165Zp7vxzJM7ud23dKVS\nieXLl+P06dP47LPP8OGHH9rdpTNo0CAUFRU1m8+5c+eaPUmxPYQQmD9/Pi5cuID3338frq6uDpnv\nrbglQGSHqyuwaX3rxjz9NPDkk/b7Al8EqqvaX5ec5efnw2Qy4YEHHkCPHj3g7u5uXQGPHDkS69at\nwx/+8AfU1tZi/fpW/vFucuDAAfTv3x8GgwGenp5QqVR2V8ARERHw8PDAq6++iueeew6HDx/Ghx9+\niJUrV7b4va5fvw6LpekstNraWgBAjx49AAALFy7Ed999h3379lnbOgK3BIgcRKGwvxXg5ta1T4P0\n8Wmqv6NePj4tq6O2thbLli3DgAEDMHDgQFy8eBEpKSkAmnbhjBgxAjqdDtHR0UhISLjrjdVu7r/5\nwHJ5eTlmzJgBb29vGAwGREZGYvbs2TbjVSoVdu/ejY8++ggDBgzAM888g61bt1r329/tYHVhYSE8\nPDwwbNgwKBQK9OzZEyEhIQCatijeeustnDp1Cmq1Gp6envD09LQeb3AkPk+A6BYfH81H7D+nwLIu\n32HzHPJ8IuKC4rDxfxIdNs+OwM9b13O7v1lL/5bcEiAikjGGABGRjDEEiIhkjCFARCRjDAEiIhlj\nCBARyRhDgIhIxnjFMMlWTV2N3WcGXK3npb0kHwwBki2/VwehvrEBCjS/qrNRAMqr4Q5/P7MZyMuz\nbVepgMDArn1VsZy5uLjg+++/b3ZH066EIUCyVVVbjejj16BUqGz6QkMd+159+gA5B4ETW237zp4F\nTp4EfrpjQKfTd21fVFzvuOdL+rj7wLy0Gz5fsotgCJCsvf8e4OHe8e9jMABxU4FEO+EyfHjT7as7\nq4rrFRArOu5WEopVztsEqq+vt7kltNzxwDCRE/Rx74Mn/v0E3Fa72by+eaQXiq8WSl1ip6fT6bBu\n3TqMGDECffr0QUJCgvXOm5s3b4Zer0e/fv3wyCOPoLS01DrOxcUFaWlp0Ov1uO+++3Dw4EFotVq8\n9tpr8PX1xaBBg7Bz507s2bMHQUFB6NevH1JTU63jv/jiC4wdOxY+Pj4YNGgQFi1ahLrOnNqtJTq5\nLlAidVH4o1JcrbE45b0aGhtEbX2t3VePZ4eKf316yil13I29zxtWduxnsKXz1+l0YsyYMaK0tFSY\nzWYREhIi3nzzTbF//37Rv39/ceLECVFbWysWLVokfvGLX1jHKRQKMXHiRFFRUSGuX78uDhw4IJRK\npVi9erWor68XmzdvFv369ROzZs0S1dXV4vTp06Jnz56isLBQCCHEsWPHRG5urmhoaBCFhYUiJCRE\nrF+/vtn8z54969h/lFa43TqypetObgkQOYGLwgVurm52XwpukLfYb37zG6jVavj4+CA2NhYnT57E\nu+++i/nz52PkyJFwc3NDSkoKPv/8c5w/f946btmyZejTp4/1vvwqlQovvfQSXF1dMXPmTJjNZixZ\nsgS9evWCwWCAwWDAyZMnAQCjRo1CeHg4XFxcMGTIEDz55JM4ePCgJMvfEfi/j4i6DLVabf3Zw8MD\n1dXVKCkpgb+/v7W9V69e6NevX7PHPN76GMh+/fpZ7/V/4yHytz5+8urVqwCaHmYzZcoUDBw4EN7e\n3njppZdw6dIlxy+cRBgCRNSlDRo0COfOnbP+fvXqVVy6dAkajcbadrcHzNzJwoULYTAY8P333+Py\n5ctYs2YNGhsb21VzZ8IQIKIuSfz0wJTExERs2bIFp06dQm1tLX7/+98jIiKi2dZBe1RXV8PT0xMe\nHh747rvv8Le//c0h8+0seK4UEd2Rj7tPh57G6ePewudL3uLG4xsnTJiA1atXY9q0aaioqMADDzyA\nzMzMZtPZG3un32/2+uuv48knn8Srr76KsLAwJCQk4MCBAy0a2xXw8ZIkW4rlKlz9/TV4uNteLOZM\nPZ8bjq2PbsP08cMlrQPg560rkuzxkikpKRg6dChCQ0Mxa9Ys1NbWwmw2IyoqCkFBQZg4cSIqKyub\nTa/X6xEcHIy9e/da248dO4bQ0FDo9XosXry4reUQEVEbtCkECgsLsXnzZhw/fhxff/01GhoakJmZ\nidTUVERFRSE/Px8TJkywXnCRl5eH7du3Iy8vD9nZ2XjqqaesCbVw4UKkp6fDaDTCaDQiOzvbcUtH\nRER31KYQ8PLygkqlwrVr11BfX49r165h0KBB2LVrF5KTkwEAycnJ2LlzJwAgKysLiYmJUKlU0Ol0\nCAwMRG5uLkpLS1FVVYXw8KabdSUlJVnHEBFRx2vTgeG+ffvi+eefh7+/P3r27IlJkyYhKioK5eXl\n1nNt/fz8UF5eDgAoKSlBRESEdbxWq0VxcTFUKhW0Wq21XaPRNDu394aVK1daf46MjERkZGRbyiai\nu/Dx8enyBzrlxsen6cB6Tk4OcnJyWj2+TSFw9uxZrF+/HoWFhfD29saMGTOwbdu2ZtPcOHLvCDeH\nABF1HLOZd/Psqm79grxq1aoWjWvT7qCjR49i3Lhx6NevH5RKJR577DF8/vnnUKvVKCsrAwCUlpbC\n19cXQNM3/KKiIut4k8kErVYLjUYDk8nUrP3mCzyIiKhjtSkEgoODceTIEdTU1EAIgX379sFgMCA2\nNhYZGRkAgIyMDEydOhUAEBcXh8zMTFgsFhQUFMBoNCI8PBxqtRpeXl7Izc2FEAJbt261jiEioo7X\npt1BI0aMQFJSEu6//364uLhg1KhRePLJJ1FVVYX4+Hikp6dDp9Nhx44dAACDwYD4+HgYDAYolUqk\npaVZdxWlpaVhzpw5qKmpQUxMDKKjox23dEREdEe8WIxkixeLUXfW4ReLERFR18d7B1G3VlQE/HSu\nAhHZwRCgbu3RR4GaGsDDw7ZPEQO4ujq/JqLOhCFA3Vp9PfC//wuMHGnbp1oNuHCHKMkcPwJERDLG\nLQHq1i4E/hl/OHYCPgW2fQ2NDc4v6DYaRIPdehQKBVwU/K5GHYchQN2aefA7GNJ7OiLuHWLTN8Mw\nAypXaU8PBQDX2gGYdeB+zMqx7evt1hvmF8xwdeHBC+oYDAHq9n6hnoKZI+wcFOgk7j20H9u2AcPt\nXCagfFkJAV4nQx2H25lERDLGLQEiifXqBfz854DSzqexYRHQ2Ah+XaMOwxAgktiBA8C1a/b7+m34\nKQSIOghDgEhi7u5NLyIpcCOTiEjGGAJERDLGECAikjGGABGRjDEEiIhkjCFARCRjDAEiIhljCBAR\nyRhDgIhIxhgCREQyxhAgIpIxhgARkYwxBIiIZIwhQEQkYwwBIiIZYwgQEckYQ4CISMYYAkREMsYQ\nICKSMYYAEZGMMQSIiGSMIUBEJGMMASIiGVNKXQAR3VlNDSAabNtVKkDJTzC1U5u3BCorKzF9+nSE\nhITAYDAgNzcXZrMZUVFRCAoKwsSJE1FZWWmdPiUlBXq9HsHBwdi7d6+1/dixYwgNDYVer8fixYvb\ntzRE3dDAgUDfvs1fffoAY8dKXRl1B20OgcWLFyMmJgbffvstvvrqKwQHByM1NRVRUVHIz8/HhAkT\nkJqaCgDIy8vD9u3bkZeXh+zsbDz11FMQQgAAFi5ciPT0dBiNRhiNRmRnZztmyUg2KmoqcKrslN2X\ncK2Rurx2cXUFqqubtgZufn31FXDlitTVUXfQpo3Jy5cv49ChQ8jIyGiaiVIJb29v7Nq1CwcPHgQA\nJCcnIzIyEqmpqcjKykJiYiJUKhV0Oh0CAwORm5uLIUOGoKqqCuHh4QCApKQk7Ny5E9HR0Q5aPJKD\nJR8vwYGCA/Dp6WPT51LvCR83XwmqIuoa2hQCBQUFGDBgAObOnYtTp07hZz/7GdavX4/y8nL4+fkB\nAPz8/FBeXg4AKCkpQUREhHW8VqtFcXExVCoVtFqttV2j0aC4uNjm/VauXGn9OTIyEpGRkW0pm7qp\n2ro6/D48FY/pZ9n0Pfgg4DvH+TU50icFn8BV4dqsragMqPUaCMAgTVHU6eTk5CAnJ6fV49oUAvX1\n9Th+/Dg2bdqE0aNHY8mSJdZdPzcoFAooFIq2zN7GzSFAdKsvvwSycoDlhbZ9KhXgY7uB0GVMN0zH\n2sNrbdovXbmK0vFXAOQ5vyjqlG79grxq1aoWjWtTCGi1Wmi1WowePRoAMH36dKSkpECtVqOsrAxq\ntRqlpaXw9W3aDNdoNCgqKrKON5lM0Gq10Gg0MJlMzdo1Gk1bSiIZa2gA5s0D/vprqStxvMzpmXbb\ns788g7iCOCdXQ91Rmw4Mq9VqDB48GPn5+QCAffv2YejQoYiNjbUeJ8jIyMDUqVMBAHFxccjMzITF\nYkFBQQGMRiPCw8OhVqvh5eWF3NxcCCGwdetW6xgiIup4bT7LeOPGjXj88cdhsVgQEBCALVu2oKGh\nAfHx8UhPT4dOp8OOHTsAAAaDAfHx8TAYDFAqlUhLS7PuKkpLS8OcOXNQU1ODmJgYHhQmInIihbhx\nrmYnpVAo0MlLJInpnp+Fh/VT8Ndf2x4Y7q6yvzyDuMw4WNadkboU6qRauu7kbSOIiGSMIUBEJGMM\nASIiGWMIEBHJGEOAiEjGGAJERDLGECAikjGGABGRjDEEiIhkjCFARCRjDAEiIhljCBARyRhDgIhI\nxhgCREQyxhAgIpIxhgARkYwxBIiIZIwhQEQkYwwBIiIZYwgQEckYQ4CISMaUUhdARG0jBHD6tP2+\nQYMAHx/n1kNdE0OAqAvy9gZcXYD4eNu+K1eAceOA7dudXxd1PQwBoi6ob19giM7+lsAHHwDbtjm9\nJOqiGALUJRiNwO7d9vsuX3ZuLZ1F5fVKvHXsLZv2E2agoucoAPc7vyjqchgC1CXs2AFkZQE//7lt\nn1oN6PXOr0lKGi8NZhhm4GjJUZu+3MsFuOCbDeAD5xdGXQ5DgLoEY+PHuDh5Hb4JtO2rKP8Kar84\n5xclod5uvbEpZpPdvhe2fIBtJu4PopZhCFCXUIyj6Cn64bfj5tr0KaDAA/4PSFAVUdfHEKAuow/u\nxcSAiVKXQdSt8GIxIiIZYwgQEckYQ4CISMYYAkREMsYQICKSMYYAEZGMMQSIiGSszSHQ0NCAsLAw\nxMbGAgDMZjOioqIQFBSEiRMnorKy0jptSkoK9Ho9goODsXfvXmv7sWPHEBoaCr1ej8WLF7djMYiI\nqC3aHAIbNmyAwWCAQqEAAKSmpiIqKgr5+fmYMGECUlNTAQB5eXnYvn078vLykJ2djaeeegpCCADA\nwoULkZ6eDqPRCKPRiOzsbAcsEhERtVSbQsBkMmHPnj1YsGCBdYW+a9cuJCcnAwCSk5Oxc+dOAEBW\nVhYSExOhUqmg0+kQGBiI3NxclJaWoqqqCuHh4QCApKQk6xgiInKONt024tlnn8Vrr72GK1euWNvK\ny8vh5+cHAPDz80N5eTkAoKSkBBEREdbptFotiouLoVKpoNVqre0ajQbFxcV232/lypXWnyMjIxEZ\nGdmWsomIuq2cnBzk5OS0elyrQ+DDDz+Er68vwsLCbvuGCoXCupvIEW4OASIisnXrF+RVq1a1aFyr\nQ+Czzz7Drl27sGfPHly/fh1XrlzB7Nmz4efnh7KyMqjVapSWlsLX1xdA0zf8oqIi63iTyQStVguN\nRgOTydSsXaPRtLYcIiJqh1YfE3jllVdQVFSEgoICZGZm4pe//CW2bt2KuLg4ZGRkAAAyMjIwdepU\nAEBcXBwyMzNhsVhQUFAAo9GI8PBwqNVqeHl5ITc3F0IIbN261TqGiIico923kr6x2+fFF19EfHw8\n0tPTodPpsGPHDgCAwWBAfHw8DAYDlEol0tLSrGPS0tIwZ84c1NTUICYmBtHR0e0th4iIWkEhbpze\n00kpFAp08hLJCaJWr8G1ums4/PIaqUvp9F7Y8gG2fbUNJX/h4yXlrKXrTl4xTEQkYwwBIiIZYwgQ\nEckYQ4CISMYYAkREMsYQICKSsXZfJ0DkKBcvAv/4B2DvrLbCQsCXF5QTORxDgDqN//4XePNN4NFH\nbfvUfkBgkPNr6op6ufZBqfdOuL7sarf/tYf+jOfG8fkd1IQhQJ1GVb0ZHpHvwn+67abAt2c/xyD1\nCAmq6nomBf0Sb6+os7tF9eN9f8IuXMRz45xfF3VODAHqNPKufoqzvn9Gvvlhm757fO7Bw3rbdrIV\nEQGcK7C/FTDhZRdYGpxcEHVqDAHqVLxqhmPj5I1Sl0EkGzw7iIhIxhgCREQyxhAgIpIxhgARkYwx\nBIiIZIwhQEQkYwwBIiIZYwgQEckYQ4CISMYYAkREMsYQICKSMYYAEZGMMQSIiGSMIUBEJGMMASIi\nGWMIEBHJGEOAiEjGGAJERDLGECAikjGGABGRjPFB8+RU1+uvI+W/KbA0WGz6ciryJaiISN4YAuRU\n3xQV4ZV9mzCo8Lc2fdeu3Y+Q3hESVCUfCgDffQskJNjpUwDLlwMhIU4viyTEECCnqjAD4lpfvLtw\nmd3+e+5xckEyM2oU4NUXmNrftm/DBuDUKYaA3DAEyOlcFMADD0hdhTx5eSlwxuUD7FKdtekrGatA\nce0fAQQ7vzCSTJsODBcVFeHBBx/E0KFDMWzYMLzxxhsAALPZjKioKAQFBWHixImorKy0jklJSYFe\nr0dwcDD27t1rbT927BhCQ0Oh1+uxePHidi4OEd3J3JFz8dL4lzAlaIrN66rqB5yrPSl1ieRsog1K\nS0vFiRMnhBBCVFVViaCgIJGXlyd+97vfibVr1wohhEhNTRVLly4VQghx+vRpMWLECGGxWERBQYEI\nCAgQjY2NQgghRo8eLXJzc4UQQkyePFl89NFHzd6rjSVSJ7X3aL5QPRcodRlkx+BnZ4pFf/+n1GWQ\ng7R03dmmLQG1Wo2RI0cCAHr37o2QkBAUFxdj165dSE5OBgAkJydj586dAICsrCwkJiZCpVJBp9Mh\nMDAQubm5KC0tRVVVFcLDwwEASUlJ1jFERNTx2n1MoLCwECdOnMCYMWNQXl4OPz8/AICfnx/Ky8sB\nACUlJYiI+P9nfWi1WhQXF0OlUkGr1VrbNRoNiouLbd5j5cqV1p8jIyMRGRnZ3rKJiLqVnJwc5OTk\ntHpcu0Kguroa06ZNw4YNG+Dp6dmsT6FQQKFQtGf2VjeHABER2br1C/KqVataNK7NVwzX1dVh2rRp\nmD17NqZOnQqg6dt/WVkZAKC0tBS+vr4Amr7hFxUVWceaTCZotVpoNBqYTKZm7RqNpq0lERFRK7Up\nBIQQmD9/PgwGA5YsWWJtj4uLQ0ZGBgAgIyPDGg5xcXHIzMyExWJBQUEBjEYjwsPDoVar4eXlhdzc\nXAghsHXrVusYIiLqeG3aHXT48GFs27YNw4cPR1hYGICmU0BffPFFxMfHIz09HTqdDjt27AAAGAwG\nxMfHw2AwQKlUIi0tzbqrKC0tDXPmzEFNTQ1iYmIQHR3toEUjIqK7Ufx0KlGnpVAo0MlLpFb4v8eM\nePjdGFjWGaUuhW7h/1wCpgZPxRtP2rmnBHU5LV138i6iREQyxhAgIpIxhgARkYwxBIiIZIwhQEQk\nY7yVNHWI4GDgp+sGm6nzBBQ8+aTT+vpr4O9/t21XqYAnngDc3JxfE3UshgB1iB9+AEwm25XG2Qpg\n5m5paqI78/cHXAqA48dt+/71L2D0aCA01Pl1UcdiCFCHEF7nUVx3Bapb2svqzsFBt5QiB7t3iBv2\nNTyHEnfbe85c81DhQs3HAAY6vzDqUAwB6hD1/xOCWf8eAlcX28NOEVo+R7gz2jh5I0qqSuz2jVwb\nh0rLRTAEuh+GAHUM1zp8Of8kevfkTuSuwtvdG97u3nb7FI3uTq6GnIVnBxERyRhDgIhIxhgCREQy\nxhAgIpIxhgARkYwxBIiIZIwhQEQkYwwBIiIZYwgQEckYrximNlu4+2l8VvSZ/U5FIxTgTYKIOjs+\naJ7arPfzYbDkvABlZbBNn5d7bxR/pYerqwSFkcN5PD8cdVf6wqXOy6avh/BBxdv/gKsrQ78zaem6\nk1sC1GaNjUDGumAkRoZJXQp1sE8WvQPjj+ft9iXtmYqGxrfhysTvkhgCRHRXEbqRiNCNtNuX9B9u\nAXRlPDBMRCRjDAEiIhljCBARyRhDgIhIxnhgmO5o2f5l+PvRv9vtq/G6AndXDydXRESOxBCgOyo0\nn8fS+1Mw5Z4ZNn3TH1NCN9X2vHGSmUYl7nt9jN2LA4f1H4Vdv7L/JYI6B4YA3dGXXwJZn/bGhvN9\nbfpUKqCvbTPJTHLNV/ix7IpN+wVLIT7u/0cJKqLWYAjQHdU3AHPmAGm/lroS6qz+8fp9dtv3fOGJ\nR7c7uRhqNR4YJiKSMW4JEFGHEYp6XLh6wW6ft7s33FzdnFwR3YohQEQdwtPNG/WoweBUg01fg6IW\no70ewWe/3SpBZXQzhgBh4yc78PLhl+z2mT3KoFRMdXJF1B2MCx2I9x8ohsVi27ft2Ps4Xv2u84si\nGwwBieUm9ayZAAAGQElEQVTk5CAyMlLSGnYf/h74YQJi+jxv06dQKPDML+9p87w7w/J1pO68fO1d\nNldX4NFH7fcdvQoc/6bNs3aI7vy3a41OEQLZ2dlYsmQJGhoasGDBAixdulTqkpzGWf8RD377Debs\neBoCtvcXL79+Hj/TzkbGar3D37e7f9C68/J15LL1cHNFieogevz6F3b7fzXsWWx65jYJ4iDd+W/X\nGpKHQENDA5555hns27cPGo0Go0ePRlxcHEJCQqQurcsRQuCTgk9gabDd/v7Hwc9RdqEeM3xSbAcq\ngbkPD3NChURN/jBzMkaF/NvuQ09W/Psf+O+F3ThiGmjTp4ACYQPDeEDZgSQPgS+++AKBgYHQ6XQA\ngISEBGRlZTEE7mD3f39AifmyTXvhpWKs/SEenub/Y9NXVwfchwS8s9r+Ny8iZ3JX9cBjPxtvt+/j\nT83Y/G0qxq5aYtvZ14hH703GtIjRNl0NogHj/COgdLFdrfVS9cKAXgPaXXd3JPnjJd977z18/PHH\n2Lx5MwBg27ZtyM3NxcaNG5sKVPCBFUREbdElHi95t5U8ny9MRNRxJL9iWKPRoKioyPp7UVERtFqt\nhBUREcmH5CFw//33w2g0orCwEBaLBdu3b0dcXJzUZRERyYLku4OUSiU2bdqESZMmoaGhAfPnz+dB\nYSIiJ5F8SwAAJk+ejDNnzuD777/HsmXL7E6zceNGhISEYNiwYd32OoJ169bBxcUFZrNZ6lIc6ne/\n+x1CQkIwYsQIPPbYY7h82fbMpq4oOzsbwcHB0Ov1WLt2rdTlOFRRUREefPBBDB06FMOGDcMbb7wh\ndUkO19DQgLCwMMTGxkpdisNVVlZi+vTpCAkJgcFgwJEjR24/segCPvnkE/HQQw8Ji8UihBDixx9/\nlLgixzt//ryYNGmS0Ol04tKlS1KX41B79+4VDQ0NQgghli5dKpYuXSpxRe1XX18vAgICREFBgbBY\nLGLEiBEiLy9P6rIcprS0VJw4cUIIIURVVZUICgrqVssnhBDr1q0Ts2bNErGxsVKX4nBJSUkiPT1d\nCCFEXV2dqKysvO20nWJL4G7+9re/YdmyZVCpVACAAQO63/m+zz33HF599VWpy+gQUVFRcHFp+q82\nZswYmEwmiStqv5uvb1GpVNbrW7oLtVqNkSNHAgB69+6NkJAQlJSUSFyV45hMJuzZswcLFizodmcg\nXr58GYcOHcK8efMANO1y9/b2vu30XSIEjEYjPv30U0RERCAyMhJHjx6VuiSHysrKglarxfDhw6Uu\npcO9/fbbiImJkbqMdisuLsbgwYOtv2u1WhQXF0tYUccpLCzEiRMnMGbMGKlLcZhnn30Wr732mvXL\nSXdSUFCAAQMGYO7cuRg1ahR+9atf4dq1a7edXvIDwzdERUWhrKzMpn3NmjWor69HRUUFjhw5gi+/\n/BLx8fH44YcfJKiy7e60fCkpKdi7d6+1rSt+M7nd8r3yyivWfa5r1qyBm5sbZs2a5ezyHE4uFzFW\nV1dj+vTp2LBhA3r37i11OQ7x4YcfwtfXF2FhYcjJyZG6HIerr6/H8ePHsWnTJowePRpLlixBamoq\nXn75ZfsDnLOHqn2io6NFTk6O9feAgABx8eJFCStynK+//lr4+voKnU4ndDqdUCqVYsiQIaK8vFzq\n0hxqy5YtYty4caKmpkbqUhzi888/F5MmTbL+/sorr4jU1FQJK3I8i8UiJk6cKP7yl79IXYpDLVu2\nTGi1WqHT6YRarRYeHh5i9uzZUpflMKWlpUKn01l/P3TokHj44YdvO32XCIE333xTLF++XAghxJkz\nZ8TgwYMlrqjjdMcDwx999JEwGAziwoULUpfiMHV1deLee+8VBQUFora2ttsdGG5sbBSzZ88WS5Ys\nkbqUDpWTkyOmTJkidRkON378eHHmzBkhhBArVqwQL7zwwm2n7TS7g+5k3rx5mDdvHkJDQ+Hm5oZ3\n3nlH6pI6THfczbBo0SJYLBZERUUBAMaOHYu0tDSJq2qf7n59y+HDh7Ft2zYMHz4cYWFhAICUlBRE\nR0dLXJnjdcfP3MaNG/H444/DYrEgICAAW7Zsue20kt9AjoiIpNP9Do0TEVGLMQSIiGSMIUBEJGMM\nASIiGWMIEBHJGEOAiEjG/h8TOXQwnIjQAwAAAABJRU5ErkJggg==\n"
}
],
"prompt_number": 1
},
{
"cell_type": "markdown",
"metadata": {},
"source": "####Characterizing p.d.f.s"
},
{
"cell_type": "markdown",
"metadata": {},
"source": "If $X_i \\sim U[0, 1] \\ \\ $ and $Y_i = X_i - 1/2\\ $, then $Y_i \\sim U[-1/2, 1/2]\\ \\ $\nand $$\\varphi(Y) = \\int_{-1/2}^{1/2} e^{i t y} dy = \\int_{-1/2}^{1/2} \\cos(ty) dy = \\frac{2}{t} \\sin\\frac{t}{2}.$$\n\nso if $Y = \\sum_{i = 1}^{12} Y_i$ for independent $Y_i$, then $\\varphi(Y) = \\left( \\frac{2}{t} \\sin\\frac{t}{2} \\right)^{12}.$\n\nUse the inversion formula to produce a way to numerically evaluate the p.d.f. of Y: \n\n$$p_Y(y) = \\int_{-\\infty}^\\infty e^{-ity} \\left( \\frac{2}{t} \\sin\\frac{t}{2} \\right)^{12} dt = \\int_{-\\infty}^\\infty \\cos tx \\left( \\frac{2}{t} \\sin\\frac{t}{2} \\right)^{12} dt$$ "
},
{
"cell_type": "code",
"collapsed": false,
"input": "import scipy.integrate\nimport scipy.stats\n\ndef inverse_transform(x):\n def integrand(t):\n return 1 / (2 * np.pi) * (2 * np.sin(t/2) / t)**12 * np.cos(t*x)\n \n p, _ = scipy.integrate.quad(integrand, -np.infty, np.infty)\n return p\n\nsum_of_uniform_pdf = np.vectorize(inverse_transform)\nnorm_pdf = scipy.stats.norm().pdf\n\nxs, step= np.linspace(-6, 6, 500, retstep=True)\nps_uniform = sum_of_uniform_pdf(xs)\nps_normal = norm_pdf(xs)",
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 2
},
{
"cell_type": "code",
"collapsed": false,
"input": "ps_uniform.sum() * step",
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "pyout",
"prompt_number": 3,
"text": "0.99999999913984783"
}
],
"prompt_number": 3
},
{
"cell_type": "code",
"collapsed": false,
"input": "plt.plot(xs, ps_uniform, label=\"sum of uniforms\")\nplt.plot(xs, ps_normal, label=\"N(0, 1)\")\nplt.legend()\nplt.title(\"Both pdfs\")",
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "pyout",
"prompt_number": 4,
"text": "<matplotlib.text.Text at 0xb4d976c>"
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEICAYAAABWJCMKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlcVOX+B/DPYRhA9kVEYJCRRUFkS1BR0HEB1BRBu4aZ\nS6KSpTfbs1LRFrW08kqZmZmmGZYlpkguOAok4E9RSlxQQRbBFRBEGGY4vz/MSQR0mBk4s3zf9zVX\nzszzzPkcGuY785zlYViWZUEIIUTvGHAdgBBCCDeoABBCiJ6iAkAIIXqKCgAhhOgpKgCEEKKnqAAQ\nQoieogJASBtmzJiBRYsWdchzi8ViuLi4yJfPnz+PgIAAWFpaIjExsUPWScijqAAQnSEUCmFqagoL\nCwvY2tpi7NixKC0tVajv999/j7CwsGb3MQwDhmE6ImoLn3zyCUaMGIE7d+5g3rx5nbJOQqgAEJ3B\nMAz27NmDmpoalJeXw8HBAfPnz+c6lkKuXLmCPn36cB2D6BkqAEQnGRsbY+LEicjPz5ffV11djWnT\npqFbt24QCoX46KOPwLIszp49i7lz5+LYsWPybw8P3L59G2PHjoWlpSUGDhyIy5cvt7q+oqIiGBgY\nYMOGDXB2doaTkxNWr14tf/zevXuYMWMGbG1t4ePjg+PHj8sfGz58OMRiMebNmwdLS0sUFBQgJSUF\nPj4+sLS0hEAgaPZchKiLIdcBCFGnB1c2qaurQ1JSEkJCQuSPzZ8/HzU1NSgsLMTNmzcREREBR0dH\nzJw5E19//TW+/fZbpKenN3uun376CampqQgMDMT06dPx3nvvYfv27W2uXywW4+LFi7h06RKGDx+O\ngIAAjBgxAkuXLkVhYSEuX76M2tpajBo1Sj68lJaWhmHDhmHq1KmYOXMmAGDIkCH45ZdfMHjwYFRX\nV7dZeAhRBX0DIDqDZVlER0fDxsYG1tbWOHToEN544w0AgEwmQ1JSEpYvXw4zMzO4urri9ddfxw8/\n/CDv+yiGYTBhwgQEBQWBx+NhypQpOHXq1GMzLFmyBF26dEHfvn3xwgsvyIvFzz//jPfeew/W1tYQ\nCAR45ZVXWqzz4WUjIyOcOXMGd+7cgZWVFQIDA1X63RDSGioARGcwDIPk5GRUVlaioaEBa9euxdCh\nQ3H9+nXcvHkTjY2NcHV1lbfv0aMHysrKHvucDg4O8p+7dOmC2trax7Z/+MieHj16oLy8HABw9erV\nFo+1lv+BnTt3IiUlBUKhECKRCFlZWY9dLyHKoAJAdBLDMIiJiQGPx0NGRga6du0KPp+PoqIieZvi\n4mIIBAJ5e3UoLi5u9rOTkxMAwNHRscVjjxMUFIRdu3bhxo0biI6OxqRJk9SSj5CHUQEgOuXBMArL\nsvJvA97e3uDxeJg0aRLee+891NbW4sqVK/j888/x/PPPA7j/Sb+0tBSNjY0tnqs9PvzwQ9y7dw9n\nzpzB999/j2effRYAMGnSJCxfvhxVVVUoLS3F2rVr28ze2NiIbdu2obq6GjweDxYWFuDxeO3OQsiT\nUAEgOmXcuHGwsLCAlZUVFi1ahC1btsDb2xsAsHbtWpiZmcHNzQ1hYWGYMmUKXnjhBQDAiBEj4OPj\ng+7du6Nbt24AWj8P4EnfFIYOHQoPDw+MHDkSb775JkaOHAng/r4BV1dX9OzZE6NGjcK0adMe+9xb\nt25Fz549YWVlhW+++Qbbtm1T7RdDSCsYmhCGENUVFRXBzc0NUqkUBgb0uYpoB5VeqampqfDy8oKn\npydWrlzZZrvjx4/D0NAQO3fubHdfQgghHUPpAiCTyTBv3jykpqYiPz8f27dvx9mzZ1tt9/bbb2PU\nqFHt7kuINumsy0YQoi5KF4CcnBx4eHhAKBSCz+cjNjYWycnJLdqtXbsWzzzzDOzt7dvdlxBtIRQK\nIZPJaPiHaBWlzwQuKytrdlyzQCBAdnZ2izbJyclIS0vD8ePH5Z+QFOkL0CcqQghRliK7d5X+uKLI\nm/OCBQuwYsUKMAwDlmXlgdrzxv6gny7elixZwnkG2j7Fb01NTfjm/75Bt0+7wWmcE2J/iYXH/zzg\nv84fmcWZnOej/3a0fQ9uilL6G4CzszNKSkrkyyUlJfKTah44ceIEYmNjAQA3b97Evn37wOfzFepL\niCaRNkkRvycepypOIXVKKpJrk5EwMQFNbBN+PfsrYpJi8OGwDzG732yuoxKiMKULQFBQEAoKClBU\nVAQnJyckJSW1uEjWwxeweuGFFzBu3DhERUVBKpU+sS8hmoJlWby450WUVJdgz8Sj2LjODJs2JSMp\nCXB1NcDUqc/g6PRAjPhhKCyNLfFs32e5jkyIQpQuAIaGhkhMTERkZCRkMhni4uLg7e2N9evXAwDi\n4+Pb3VffiEQiriN0KF3Zvq//72tkl2VjhUcWnuprhogIYP58EUaPBs6cAdauBZq+cse3iSmYum8E\nfLr5oG+3vlzHVomu/Ldri65vn6I0+kSwB/sOCOFKwa0ChGwMQYIgC8sWeOCnn4Dhw5u3aWoCPv0U\n+PJLYN7G77Dt0hrkzMqBsaExN6GJ3lP0vZMKACFtYFkWI7aMgJ9xFLbNX4ADB4CAgLbbf/UVsGo1\ni95LxiOs50C8G/Zuu9Zna2uLyspKFVMTfWJjY4Pbt2+3uJ8KACEq2nFmBz46sgI3Pj6ODet5ePrp\nJ/d55RUgv/wycoP6I29uHpwsnBReH73eSXu19ZqhAkCIChpljfD5ygeCvC/hZxaOL75QrF9DAzB4\nMNDtuXfg6HEDG6M2KrxOer2T9lK1ANBpi4S0YvPpzTCVuqBYPBIff6x4P2NjYMMG4PgXb+G3s7tQ\nWFnYcSEJUREVAEIe0cQ24dPMT1GVvBiff8bA1LR9/QMDgedibCG8+SJWZK7omJCEqAEVAEIekVKQ\ngvo75nDjDcHYsco9x+LFQNH2V/HTXztwrfaaegMSuczMTHh6esLCwgK7d+9W63PPnTsXH374oXx5\n3bp1cHBwgKWlpc7srKd9AIQ8Yvj3I5C35QXs+fh5DByo/PMsWwb8UDUH08a7YNHQRU9sT6/39hsx\nYgSio6Mxf/78Dl1PY2MjrKyskJOTg759NeccD9oHQIganao4hVOl5+CDSSq9+QPAq68CN/fOR2L2\nOkhkEvUEJM0UFxejT58+Hb6eiooK1NfXK3XCanuvz9OZqAAQ8pDEnC/BP/Uy3n3bSOXnsrAA5k70\nBb+6N3ad26WGdNxauXIlBAIBLC0t4eXlhcOHDwMAZsyYgUWL/v2GIxaLm13tVygUYtWqVfDz84OF\nhQXi4uJw7do1jB49GlZWVggPD0dVVVWb692wYQM8PT1hZ2eH8ePHo7y8HADg7u6Oy5cvY9y4cbC0\ntGw2n/MDBgYGzS5J83BWsVgMgUCAzz77DA4ODnBycsL333/fom1BQQG8vLwAANbW1vJpPv/8808E\nBwfD2toa/fv3x7Fjx+R9RSIR3n//fQwePBjm5ua4fPkyDAwMsG7dOnh6esLS0hKLFy/GpUuXEBIS\nAmtra8TGxsq34ebNmxg7dixsbGxgZ2eHIUOGdEwRYTWYhscjOuau5C5r/oEN2zekjG1qUs9zXr3K\nsqYDf2BHfDf6iW01+fV+7tw51sXFhS0vL2dZlmWvXLnCXrp0iWVZlp0xYwa7aNEiedvDhw+zAoFA\nviwUCtmQkBD2+vXrbFlZGdutWzc2MDCQPXXqFFtfX88OHz6cXbp0aavrPXToENu1a1c2NzeXbWho\nYOfPn88OGTKk2XMfOnSozdwMw8hzPpr18OHDrKGhIbtkyRJWKpWyKSkprKmpKVtVVdWibVFREcsw\nDCuTyViWZdlbt26x1tbW7NatW1mZTMZu376dtbGxYW/fvs2yLMsOHTqUdXV1ZfPz81mZTMZKJBKW\nYRg2OjqarampYc+cOcMaGRmxw4YNYwsLC9nq6mq2T58+7ObNm1mWZdl33nmHffHFF1mpVMpKpVI2\nIyOj1e1r6zWj6GuJvgEQ8o/fzv4G4xsD8c7LTlDXVBSOjkC0Vwz+LD6GqzVXVX4+hlH9pgwej4eG\nhgacOXMGjY2N6NGjB9zc3OSPs0/4dDp//nzY29vDyckJYWFhCAkJgb+/P4yNjRETE4Pc3NxW+23b\ntg1xcXEICAiAkZERli9fjmPHjqG4uFi5DXkkK5/Px+LFi8Hj8TB69GiYm5vj/PnzLdo+un179+5F\n7969MWXKFBgYGCA2NhZeXl7yHdEMw2DGjBnw9vaGgYEB+Hw+AOCtt96Cubk5+vTpA19fX4wePRpC\noRCWlpYYPXq0/PdgZGSE8vJyFBUVgcfjYfDgwUpv7+NQASDkH+uyvockZwYmTlTv87482wz8ixPx\nw+mtKj8Xy6p+U4aHhwe++OILJCQkwMHBAZMnT5YPxSjCwcFB/nOXLl2aLZuYmKC2trbVfuXl5XB1\ndZUvm5mZwc7ODmVlZUpsRUt2dnbNZnEzNTVtM8vDrl69ih49ejS7z9XVFVev/lvkHx4Ge+Bxv4cu\nXbrI1/3mm2/Cw8MDERERcHd377B506kAEAKguLoYJ6+exIyQKJiYqPe5Q0IA2+IZWHfse43dGaiI\nyZMnIz09HVeuXAHDMHj77bcB3H9Trqurk7erqKh44nMp+ntwcnJCUVGRfPnu3bu4desWnJ2dFepv\namraLFt5eblaZhp0dnbGlStXmt135cqVZrlUWY+5uTlWrVqFS5cuYffu3fjss8+Qlpam9PO1hQoA\nIQC+z/0ByJ+EubPV/O6P+8Mur8QMxu2qRhy/elztz98ZLly4gLS0NDQ0NMDY2BgmJibg8XgAgICA\nAKSkpKCyshIVFRX4QtHrZihg8uTJ2LRpE06fPo2Ghga8++67GDhwYItP320JCAjAtm3bIJPJkJqa\niqNHjyq87scVqTFjxuDChQvYvn07pFIpkpKScO7cOYx96MQRRYrcw20e/nnPnj24ePEiWJaFpaUl\neDye/PetTlQACAHwXfYOeNx7Dh01LcW0aQwkJ2Ox5cTPHbOCDtbQ0ICFCxfC3t4ejo6OuHnzJpYv\nXw4AmDp1Kvz9/SEUCjFq1CjExsY+8dPvw48zDNNm+xEjRuCDDz7AxIkT4eTkhMLCQvz0008K516z\nZg1+//132NjY4Mcff0RMTEybOVrL+GjOB2xtbbFnzx6sXr0aXbt2xapVq7Bnzx7Y2tq2+dytraut\n38PFixcRHh4OCwsLDBo0CC+//DKGDh2q4FYrTuUTwVJTU7FgwQLIZDLMmjVL/rXwgeTkZCxevBgG\nBgYwMDDAp59+iuH/XFD9wc4PHo8HPp+PnJyc5uHoxBjSCS7cugC/z4fiK/dSzHxB/Z+yHhj9wmnk\nuI3HzfcL23wzoNc7aQ9VTwRT6bgzqVTKuru7s4WFhaxEImH9/f3Z/Pz8Zm1qa2vlP+fl5bHu7u7y\nZaFQyN66davN51cxHiEKWbT/I5Y//mW2urpj17N7dxNr8qYne7zseKuP0+udtFdbrxlFX0sqDQHl\n5OTAw8MDQqEQfD4fsbGxSE5ObtbGzMxM/nNtbS26du36aAFSJQIhKtt8/BcMtHwGlpYdu57ISAbM\n2WewMUs7h4GI7lF6TmAAKCsra3aok0AgQHZ2dot2u3btwsKFC1FeXo79+/fL72cYBiNHjgSPx0N8\nfDxmz57dom9CQoL8Z5FIRHN5ErW6dPsSKu6W4ZNxYR2+LiMjYFSPZ7Dj7//gq5gVajkahRDg/lnN\nYrG43f1UKgCKvoCjo6MRHR2N9PR0TJ06VX6iRWZmJhwdHXHjxg2Eh4fDy8sLYWHN/xAfLgCEqNum\nnF+AcxMQ9X7Hjf0/bN7EQOxNBXIrcvGU41Odsk6i+x79cLx06VKF+qk0BOTs7IySkhL5cklJCQQC\nQZvtw8LCIJVKcevWLQCAo6MjAMDe3h4xMTEtdgIT0tG2n/gdoV3Ho0uXzlnf0KEM+JeisSnz985Z\nISGPoVIBCAoKQkFBAYqKiiCRSJCUlISoqKhmbS5duiQf5z958iSA+2ff1dXVoaamBsD9kzv2798P\nX19fVeIQ0i636m7hSn0eXhoj6rR18nhAhOtY/Ja/p9PWSUhbVBoCMjQ0RGJiIiIjIyGTyRAXFwdv\nb2+sX78eABAfH4+dO3diy5Yt4PP5MDc3lx/DW1FRgQkTJgAApFIppkyZgoiICBU3hxDF7Tz9B5gr\nIox5V/0nfz3Oi2NCsVt8EeU15XC0cOzUdRPyMJoQhuitQaunoO7sEJz6Nr5T19vYCJjPiMVHM8Px\nxog4+f30eiftRRPCEKIEWZMMJ6r+wNQBYzp93Xw+EGw1FltzdHMYaPLkyS0OB+8IiYmJeOeddzp8\nPbqMCgDRS0cvZ0Fa6Yxp0S2v2NgZZolG4UxdGuql9ZysXxlCoRAODg7NLq727bffYtiwYfLlvLw8\n5OXlYfz48fL7fvzxR7i6usLc3BwxMTHtmk93zpw58PLyAo/Hw+bNm5s9Nnv2bGzbtg03btxQYav0\nGxUAope+TtsLp9qnYW/Pzfr/83RXsNd8sffMEW4CKKmpqQlr1qxp8/H169fj+eefly+fOXMGL774\nIrZt24Zr167B1NQUL730ksLrCwgIwFdffYWnnnqqxWHnxsbGGD16NLZs2dL+DSEAqAAQPXWweA+i\nvJ/mbP1mZoAHRuG7I/uf3FhDMAyDN954A6tWrUJ1dXWrbVJTU5tdtGzbtm2IiopCaGgozMzM8MEH\nH+DXX3/F3bt3FVrnSy+9hOHDh8OkjWt0i0Qi7N27t/0bQwBQASB66EplCSqlVzE/WsVZ31UU1Tcc\nf1ZoTwEA7h/6LRKJsGrVqhaP3b17F4WFhejdu7f8vvz8fPj7+8uX3dzcYGxsjAsXLqglj5eXF06f\nPq2W59JHKh0GSog22pC2H+bXw+HVu3PO/m1L/NggrNpQhrLqcjhbKXY4KLNU9ctHsEuUP9KIYRgs\nW7YMgwcPxiuvvNLssQcTu1tYWMjvq62thZWVVbN2lpaW8nOAVGVhYdHmtxHyZFQAiN5J/vsgBnYb\nyXUMuLvxYH5jODaKD2Dx+GkK9VHlzVtdfHx8MHbsWKxYsQLeD02gYG1tDQCoqamBnZ0dgPszWz36\nBl1dXd2sSKiipqamRYEhiqMhIKJXmtgmnGs4hKmDuS8AABBsF45fTx/gOka7LV26FBs2bGg2N6+Z\nmRnc3d2bTaru4+PTbIjm0qVLkEgk6NWrl1pynD17FgEBAWp5Ln1EBYDolcyCvyCrs8J/wl2f3LgT\nTB8cgfz6A1p3Api7uzueffbZFkcEjRkzBkeO/Htk05QpU/D7778jIyMDd+/exaJFizBx4kT5ZeIT\nEhKaHUb6qMbGRtTX16OpqQkSiQT19fXNfldHjhzB6NGj1bx1+oMKANEr36YdgkAyUu0TvytrUkRP\nSOvMkVHwF9dR2m3x4sWoq6trdnjmnDlzsG3bNvlynz598PXXX2PKlClwcHDAvXv38NVXX8kfLykp\nQWhoaJvrCA8Ph6mpKbKysjBnzhyYmpoiPT0dAFBfX499+/Zh+vTpHbB1+oEuBUH0iss7YzDSLg6b\n3pzIdRQ515fnIszHHdteflMnXu9TpkzBpEmTmp0M1pbAwECkpaXBxsam3etJTExEaWkpVqxYoUxM\nnaDqpSCoABC90SCVoMuSrsj8TxFCAmyf3KGTzPnsN+y9/jWurtxPr3fSLnQtIEIUtOPPLPCre2Og\nv+a8+QPA3NHDUM77k+sYRA9RASB6Y+ufB+FtMhKaNhNjoLc1jO54P7khIWpGBYDojZybBzHeVzMO\n/3xUb2MR1xGIHqITwYheuF1biyrjPMSPGcR1lFaN9hbhrOn/aKJ40i7K7Dx/mErfAFJTU+Hl5QVP\nT0+sXLmyxePJycnw9/dHYGAg+vXrh7S0NIX7EqJOW9KyYFYTACf7Tpr8t51mR4aicYEhau7dA8uy\ndKObQrfbt2+r9sJjlSSVSll3d3e2sLCQlUgkrL+/P5ufn9+sTW1trfznvLw81t3dXeG+LMuyKsQj\npJlhyxazQe+8w3WMx+ry32D2q71HuI5BdICi751KfwPIycmBh4cHhEIh+Hw+YmNjW8wC9OBsP+D+\nRaG6du2qcF9C1Ol0ZTpGeYdxHeOx+piKsPOEmOsYRI8ovQ+grKwMLi7/zqYkEAiQnZ3dot2uXbuw\ncOFClJeXY//+/e3qC9w/VfwBkUgEkUikbGSip+41NOJ2l+OYMUIzx/8fGNtXhP/93yoAi7mOQrSM\nWCyGWCxudz+lC4CiO6uio6MRHR2N9PR0TJ06FefOnWvXeh4uAIQoY/uRkzCpc4O7szXXUR5rVkQo\nluZPQmVNPWwsNORaFUQrPPrheOnSpQr1U3oIyNnZGSUlJfLlkpISCASCNtuHhYVBKpXi9u3bEAgE\n7epLiCp+O5EOT2PNHv4BAIG9Jczq+uD7AzlcRyF6QukCEBQUhIKCAhQVFUEikSApKQlRUVHN2ly6\ndAn390cAJ0+eBADY2dkp1JcQdTl+LR3DPTS/AACAr4UIv+WKuY5B9ITSQ0CGhoZITExEZGQkZDIZ\n4uLi4O3tjfXr1wMA4uPjsXPnTmzZsgV8Ph/m5ub46aefHtuXEHVrlDbhukkGZgxbx3UUhUT5ibD8\nCO0HIJ2DLgZHdNrOo2cweXcUJKsucR1FIder78DhEydUvHYTDna0H4Aohy4GRwiApKx0CA20Y/gH\nALpZWcKyoQ82prZ+VBwh6kQFgOi0Y1fTMbSn9hQAAPC1HIK9f6dzHYPoASoARGc1NQFXeRmYOkS7\nCsCYvqH4+04G1zGIHqACQHRW2oliMPx6hPXx5DpKu0wVDcYdq2OorJJxHYXoOCoARGdtS0+HoClU\n666w6WJrD1OpE7YdzOM6CtFxVACIzjp6JR2DBNo1/PNAb9NQJOfSMBDpWFQAiE5iWaAY6YgdpJ0F\nIKJ3GE7eoh3BpGNRASA6KefvW5CZl2J0oD/XUZQydWgoKi0yUFtL58GQjkMFgOikLYcz0L1xIPg8\n7Zz0ro9jT/ANGfx6+DLXUYgOowJAdNLhy+no3107h3+A+2dyuvND8ctx2g9AOg4VAKJzWBa41JiB\nicGhXEdRicgtDMev0X4A0nGoABCdk19wFxKbvzBxwACuo6hkSlgorhlnoL6e6yREV1EBIDpnS1o2\n7Br9YWqkmRPAK2pgT18YWFYgNf0611GIjqICQHTO/vPpCLTT3vH/B3gGPPRgQpB0LJPrKERHUQEg\nOuf8vXREB2p/AQCAwS6h+LOEdgSTjkEFgOiUwuJG3LPLRuygwVxHUYvJg8NQykuHRMJ1EqKLVC4A\nqamp8PLygqenJ1auXNni8W3btsHf3x9+fn4YPHgw8vL+vb6JUCiEn58fAgMD0b9/f1WjEIKtB3Nh\nKesJOzMbrqOoxbBewYD9GWRk3+U6CtFBKp0lI5PJMG/ePBw8eBDOzs4IDg5GVFRUs+kd3dzccPTo\nUVhZWSE1NRVz5sxBVlYWgPvHOovFYtja2qq2FYT8I+XvdPja68bwDwB04XdBN9YfP6ZnYXjYCK7j\nEB2j0jeAnJwceHh4QCgUgs/nIzY2FsnJyc3ahISEwMrKCgAwYMAAlJaWNnucpnwk6vT3nQw83Ve7\nj/9/1ADHMBwtpP0ARP1U+gZQVlYGFxcX+bJAIEB2dttT2W3cuBFjxoyRLzMMg5EjR4LH4yE+Ph6z\nZ89u0SchIUH+s0gkgkgkUiUy0WEVFSzu2mXgudC1XEdRq4nBodhzYg1kMoDH4zoN0URisRhisbjd\n/VQqAO25zvrhw4fx3XffITPz30PaMjMz4ejoiBs3biA8PBxeXl4IC2v+9f3hAkDI42w/eA4mPHO4\n2gi4jqJWT/sOBuv0HE6cakT/fnyu4xAN9OiH46VLlyrUT6UhIGdnZ5SUlMiXS0pKIBC0/OPLy8vD\n7NmzsXv3btjY/LtzztHREQBgb2+PmJgY5OTkqBKH6Lndp9Lhbao74/8P2HaxhRVckSQ+zXUUomNU\nKgBBQUEoKChAUVERJBIJkpKSEBUV1axNcXExJkyYgK1bt8LDw0N+f11dHWpqagAAd+/exf79++Hr\n66tKHKLnTt1OR6S37hUAAAi0C8WBC3RdIKJeKg0BGRoaIjExEZGRkZDJZIiLi4O3tzfWr18PAIiP\nj8eyZctQWVmJuXPnAgD4fD5ycnJQUVGBCRMmAACkUimmTJmCiIgIFTeH6Ktbt4Bq63RMHvQe11E6\nxPjAMLx58lew7KvQshkuiQZjWA0+DIdhGDpKiCjk2x0leDmvH+o/uKZ1cwAr4krVFbit6I+/plWg\nTx/d2z6iXoq+d9KZwEQn/HoiHR5G2jcBvKJ6WPWACd8IOw9f5DoK0SFUAIhOOH4tA8M9dOv4/4cx\nDIO+FmHY+zftByDqQwWAaL3qauCWufZOAK+o0T6h+Ks6AzQqStSFCgDReqlHbsPA5goG9AjkOkqH\niukXigaHDBQVcZ2E6AoqAETr7cjKhCtvAAwNtHMCeEX5OvSFgcUN/H74GtdRiI6gAkC03rGydAwR\n6vbwDwAYMAboZTIIybl0XSCiHlQAiFarrQWuGafj2RDdLwAAMKJXKE7eoAJA1IMKANFqhzPqgO55\nGOKm3RPAK2pCUChqbDNQXs51EqILqAAQrZaUmY3ujC9M+aZcR+kUAwTBQNd8HDhSy3UUogOoABCt\nln4lA4ME+jH8AwAmhiZw4QdiZ3YW11GIDqACQLTWvXtAqUE6ngnWnwIAAEOEoci6SvsBiOqoABCt\n9WeWFBBkYWQv3ZgAXlETg8JwyywDt29znYRoOyoARGvtOHoKtrwesDO14zpKpwoTDgIE2RCnN3Id\nhWg5KgBEa6VdSkewg34N/wCATRcbdOUJ8WvmKa6jEC1HBYBoJYkEKJSlI+Yp/SsAADDQKRRHr9B+\nAKIaKgBEK+XksGBcMzBKR2cAe5KYfmG4ysvAP5PqEaIUlQpAamoqvLy84OnpiZUrV7Z4fNu2bfD3\n94efnx8GDx6MvLw8hfsS8ji/Hj0PU74pXKxcuI7CiWHuoTDomY4//6RLgxLlKV0AZDIZ5s2bh9TU\nVOTn52Nn8BbbAAAeIklEQVT79u04e/ZsszZubm44evQo8vLysGjRIsyZM0fhvoQ8zh/njuIpu6Fc\nx+BMD6se6MI3wa70Aq6jEC2mdAHIycmBh4cHhEIh+Hw+YmNjkZyc3KxNSEgIrKysAAADBgxAaWmp\nwn0JaYtUClyQHEF0wBCuo3Cqn30YDhXQfgCiPKWvn1tWVgYXl3+/fgsEAmRnZ7fZfuPGjRgzZky7\n+yYkJMh/FolEEIlEykYmOuLECRZwPYIxPglcR+HUOP9QvHU4HfX1M2FiwnUawiWxWAyxWNzufkoX\ngPbMvXr48GF89913yMzMbHffhwsAIQCwS1wEI2MZPGw9uI7CqZGeoeC5fYbjx4Ew/dwXTv7x6Ifj\npUuXKtRP6SEgZ2dnlJSUyJdLSkogEAhatMvLy8Ps2bOxe/du2NjYtKsvIa1JyT+KAJshOjsBvKJ8\nuvmAMb2FvUcquI5CtJTSBSAoKAgFBQUoKiqCRCJBUlISoqKimrUpLi7GhAkTsHXrVnh4eLSrLyGt\nkUqBs3ePYpy/fo//A/cniOlrNQipZzK5jkK0lNJDQIaGhkhMTERkZCRkMhni4uLg7e2N9evXAwDi\n4+OxbNkyVFZWYu7cuQAAPp+PnJycNvsS8iS5uQDT8wjG9n2V6ygaYbRPKFYcTodUOhGGuj0jJukA\nDMuyGnsgMcMw0OB4hAPvryzDqjo/1CXcgAFD5zFmFmdi5GfzcXTKSQQHc52GaApF3zvpL4holb1/\np8PPKoze/P8R7BwMmVUB/jhSxXUUooXor4hoDakUOFN7FGN9afz/ASOeEbwsBuL3vHSuoxAtRAWA\naI3cXMCg51GM8aEC8LAx3iKcrj6MpiaukxBtQwWAaI29h2+CtShBQPcArqNolCi/YYBQjDNnuE5C\ntA0VAKI1fs9LRx/LQTA0oMNdHhbkFIQm6wKkiiu5jkK0DBUAohWkUuDvOzT80xojnhF6m4Vg16mj\nXEchWoYKANEKubmAgdtRjOlDBaA1o71FyK0Ug46aJu1BBYBohX2HqyCzPo8gpyCuo2ikmIBhkLnQ\nfgDSPlQAiFbYffoovM1DYGxozHUUjRTkFATW+hJ2H7jNdRSiRagAEI0nlQJ/1R7C2L7DuY6isfg8\nPrwtQvDbSdoPQBRHBYBovNxcgHFPw3jfEVxH0Whj+4qQd0cMqZTrJERbUAEgGi/50DWwFiV4yvEp\nrqNotHF9h8HATYwTJ7hOQrQFFQCi8ZJPH0aA9VA6/v8J+jn2Q5PVZew5RPsBiGKoABCNdu8ecK4h\nDTGBNP7/JHweH32tBiH59BGuoxAtQQWAaLTMTIDneQhj+9D4vyLG+Ypw7p4Y9fVcJyHagAoA0Wg7\nDxXBsEstfOx9uI6iFUZ7DYNhrzQcO8Z1EqINVC4Aqamp8PLygqenJ1auXNni8XPnziEkJAQmJiZY\nvXp1s8eEQiH8/PwQGBiI/v37qxqF6KB959Iw0GG43s//q6ggpyCw5mVITivnOgrRAirtVZPJZJg3\nbx4OHjwIZ2dnBAcHIyoqqtn0jnZ2dli7di127drVoj/DMBCLxbC1tVUlBtFRt24BV43T8E4/Gv9X\nFM+Ah+Cuw7FXfABfYBrXcYiGU+kbQE5ODjw8PCAUCsHn8xEbG4vk5ORmbezt7REUFAQ+n9/qc9CU\nj6Qthw6xMHBPQ4Qnjf+3x6R+ESg0OIA7d7hOQjSdSt8AysrK4OLiIl8WCATIzs5WuD/DMBg5ciR4\nPB7i4+Mxe/bsFm0SEhLkP4tEIohEIlUiEy3ys/gcTLsboad1T66jaJWnvcLxuudiHEprQkw07ebT\nB2KxGGKxuN39VCoAqo7LZmZmwtHRETdu3EB4eDi8vLwQFhbWrM3DBYDol7SigxANGEHj/+3U06Yn\nLI0tsD3tL8RE+3Mdh3SCRz8cL126VKF+Kn08cHZ2RklJiXy5pKQEAoFA4f6Ojo4A7g8TxcTEICcn\nR5U4RIdcvgzc7f4Hnu0XyXUUrTSyZwQOXj5Al4cmj6VSAQgKCkJBQQGKioogkUiQlJSEqKioVts+\nOtZfV1eHmpoaAMDdu3exf/9++Pr6qhKH6JCU/Q1ocjmKcPeRXEfRSpOCwlHnuB8FBVwnIZpMpSEg\nQ0NDJCYmIjIyEjKZDHFxcfD29sb69esBAPHx8aioqEBwcDDu3LkDAwMDrFmzBvn5+bh+/TomTJgA\nAJBKpZgyZQoiIiJU3yKiE3Ycy0APbx/YdqEjxJQxvOcwNDlNxe+p9/B6ry5cxyEaimE1+DAchmHo\nKCE91NgIWEx4C/PiTbFqbALXcbSW96ehsDy5BNnbw7mOQjqZou+ddIgA0TjHjgGMZyqeCRjFdRSt\nNt43HLk1++myEKRNVACIxtmx7yoYyzIEOwVzHUWrjfeJAN9rP9LTuU5CNBUVAKJxkv/ej5DuI8Ez\n4HEdRasFOwejyaIYO/fTZSFI66gAEI1SWgpct0zFc/1p+EdVhgaGCHWMwO9n93EdhWgoKgBEo6Ts\nkwFuBzHKg44IU4fng8fips1ePHS6DiFyVACIRtl2JBsOps5wtnTmOopOGNNrFNieB7EntYHrKEQD\nUQEgGkMiAbIqd+MZ39ZPJiTtZ29mD6GZD344epTrKEQDUQEgGiMjAzDw3o3JgVQA1OnZwLE4cWcP\n7t3jOgnRNFQAiMbYmnIRhhaV6OfUj+soOuU//mNh4P07Dh6kkypJc1QAiEZgWSD57O+IcB0HA4Ze\nlurk280XJqZSbNl3jusoRMPQXxrRCH/9BdS57Mb0kHFcR9E5DMPgac+x+OPyHjQ1cZ2GaBIqAEQj\nbP+tEk3dT2CkG83+1RGeCxoLmcceHD/OdRKiSagAEI3w04l9COoqginflOsoOml4z+GQdT2N7buv\ncx2FaBAqAIRzxcXAVYvdmDaAjv7pKCaGJghzHIWf/9rFdRSiQagAEM79vOseWPdURHvT+H9HmjV4\nIm7a78SlS1wnIZpCpQKQmpoKLy8veHp6YuXKlS0eP3fuHEJCQmBiYoLVq1e3qy/RH5sz/kBvy0A4\nmDtwHUWnPd1rNCDIwpafK7mOQjSE0gVAJpNh3rx5SE1NRX5+PrZv346zZ882a2NnZ4e1a9fijTfe\naHdfoh9u3wbOGvyMmQP/w3UUnWduZI7+9iOwJWs311GIhlC6AOTk5MDDwwNCoRB8Ph+xsbFITk5u\n1sbe3h5BQUHg8/nt7kv0w45f68H02ovJARO4jqIXZodOxFXrnbh8meskRBMoPSdwWVkZXFxc5MsC\ngQDZ2dlq75uQkCD/WSQSQSQSKZWXaKavD/wBz6AAdDfvznUUvTDeayxmu76ErTtqsPgdC67jEDUR\ni8UQi8Xt7qd0AWAYRtmu7er7cAEguuXaNSCf+RkrB9HwT2exMrFCYNdQbNq/F4sRy3UcoiaPfjhe\nunSpQv2UHgJydnZGyUMXGS8pKYFAIOjwvkR3/PRLPdBrLyb7T+Q6il6JHzwJV21/oqOBiPIFICgo\nCAUFBSgqKoJEIkFSUhKiolo/jvvR2enb05form8O7UdvK38a/ulkE31iAKEY3++4xXUUwjGlh4AM\nDQ2RmJiIyMhIyGQyxMXFwdvbG+vXrwcAxMfHo6KiAsHBwbhz5w4MDAywZs0a5Ofnw9zcvNW+RH+U\nlgIXu2zD6sE0DNHZLI0tEeowGt8f2IEPMJfrOIRDDPvox3MNwjBMi28PRDd89FklllX3RMU7hbDp\nYsN1HL3z+7m9eGbNxzjxcib69uU6DVE3Rd876UxgwokNf+7AgK7h9ObPkVGeETB0KMCaLXQ8qD6j\nAkA63ZkzQLn9FrwRPp3rKHqLz+Mjptez+Cl/G6RSrtMQrlABIJ3us+8LYOhwEaM9I7mOotfmD3ke\nUu+tOHCAhln1FRUA0qkaG4Gfzm3BJK/nwOfxn9yBdJj+zv1hYdmEz3ZkcR2FcIQKAOlUe/Y2Qeaz\nBa+IaPiHawzD4KWBs3Hk7jeoruY6DeECFQDSqT755RC6WdrA38Gf6ygEwEuDZgBev2HT9iquoxAO\nUAEgnaaiAjjBrMNrQ15U6VIiRH26mXVDiP0orD74A+iIa/1DBYB0msTNZWB6ihEXPIXrKOQhS55+\nEddcvkZ2NlUAfUMFgHQKmQz4MutbjOsZCwtjugqlJhnWcyisbWRI+C6T6yikk1EBIJ1iT4oUd3tv\nwJIxdOkBTcMwDOYPmoO0O1/j5k2u05DORAWAdIolP/2KntY94evgy3UU0oqXBk0H02svvthYznUU\n0omoAJAOd/48izNWq/Dh0288uTHhhJ2pHaJ6PofE7ETIZFynIZ2FCgDpcAu/ToeFfRUm9h3HdRTy\nGMvHv4q7Xt/g5113uY5COgkVANKhrl8H9txajXeGvgYDhl5umszD1gNPdQ3Dwh2b6JBQPUF/kaRD\nLVl7FobCLMwLm8Z1FKKAVRPeQKnLZxAfbeQ6CukEVABIh6mtBTZd/gAv93sVpnxTruMQBYQJB8Hd\nTogF323jOgrpBCoXgNTUVHh5ecHT0xMrV65stc1///tfeHp6wt/fH7m5ufL7hUIh/Pz8EBgYiP79\n+6sahWiYZevOAm4HsXjUy1xHIe2wduISnLH7ECdy6TrRuk6lAiCTyTBv3jykpqYiPz8f27dvx9mz\nZ5u1SUlJwcWLF1FQUIBvvvkGc+f+exw4wzAQi8XIzc1FTk6OKlGIhqmtBdae/gDx/q/SiV9aJrzX\nUPS0dcGcxK1cRyEdTKUCkJOTAw8PDwiFQvD5fMTGxiI5OblZm927d2P69PtXfhwwYACqqqpw7do1\n+eM05aNuen/NGbA9D+LDcfO4jkKU8OV/EnDa+gMcy5FwHYV0IKUnhQeAsrIyuLi4yJcFAgGys7Of\n2KasrAwODg5gGAYjR44Ej8dDfHw8Zs+e3WIdCQkJ8p9FIhFEIpEqkUknqKoCvrr4Ft6IeZc+/Wup\niN5D4dW1N2auW4ez/V/hOg55ArFYDLFY3O5+KhUARa/o2Nan/IyMDDg5OeHGjRsIDw+Hl5cXwsLC\nmrV5uAAQ7fDy6oMwcT6PhKd/4zoKUcHWaZ+i39ph2HNwGsaOpLmbNdmjH46XLl2qUD+VhoCcnZ1R\nUlIiXy4pKYFAIHhsm9LSUjg7OwMAnJycAAD29vaIiYmh/QA64NJlGZIqX8eqyJUw4hlxHYeoIMDZ\nB8McYzDz+4/o7GAdpVIBCAoKQkFBAYqKiiCRSJCUlISoqKhmbaKiorBlyxYAQFZWFqytreHg4IC6\nujrU1NQAAO7evYv9+/fD15euE6PtJq38Bs72lpgdOoHrKEQNts5chkrXzfhwfT7XUUgHUGkIyNDQ\nEImJiYiMjIRMJkNcXBy8vb2xfv16AEB8fDzGjBmDlJQUeHh4wMzMDJs2bQIAVFRUYMKE+28SUqkU\nU6ZMQUREhIqbQ7j064GrOGWzGNkviGnCFx3R3cIBrz+1FB/tj8f8yiOwtaFTh3QJw2rwYTgMw9BR\nQlqivh7o9vIkhD/liZ0vf8R1HKJGsiYZnBYNhlfdLBz5fBbXcYgCFH3vpHJO1GLqR7+hyeEkfpj9\nPtdRiJrxDHjYOWM9Mozfxc9/lHIdh6gRFQCisgNZV7GzYS5+/M8WmBp14ToO6QChnv6Y7PYKpv82\nHbV3m7iOQ9SECgBRSe3dJsRsfgHjnOIRFTiI6zikA22e9Q7MrRsQvvgzrqMQNaECQFQyfPEKmFjd\nwc/zaOhH1/EMeDj40g/I4X+CD7fQ/MG6gAoAUdrbG/fipOGXyPjvLzAy5HMdh3QCvx49sTrseyz5\naxL+/LuM6zhERVQAiFJ+EZ/FpxdewIaIHfBycuY6DulEC54eg6e7zcPIbyag/GYd13GICugwUNJu\nx8+XIGRDKF7y/gD/i6OJXvRRUxML73dnoLL+Jq6s3IUuxvQNUJMo+t5JBYC0y/nSG/D/fAhG2MzG\n3vdf4zoO4VBdfSOEb0+AhZEVzi3fDL4hj+tI5B90HgBRuzMlpQj4YgiCujyLPe/Rm7++MzXhI3/p\nDtxqqECfRVMgkdI0ktqGCgBRyJ/nL+CptUMQxJuJ9A8SQFd6IADQ1boLLiTswc3qOni8H4Oqulqu\nI5F2oAJAnmjjkVSEbQrFGPP3cXTFm/TmT5rpZmuCy8t3oumOA3osHYz8q0VcRyIKogJA2tQoa8Tk\n9QmYs3cm3nD5Fb8tnklv/qRVNlZ8XP7iW/Spj4P/2oH48tAuriMRBdBOYNKqE8X5GLthJqoqrPHz\n899h7FAnriMRLcCywDtf/olVl6disJMIv730KexMbbmOpXdoJzBRSlV9FSZteBP91w2F4NY0FC9P\noTd/ojCGAVbOG4T0507h7F8mcPrIG8tS1kHWRDPKaCL6BkAAAJX3KrEw+Qt89/eXMCkajy8nfoyp\nMQ5cxyJarLEReHPVaXx5+RVYdL+OpeHv4MXBk8Hn0TkDHY3OAyBPxLIsjpedQMLvG3Cg7GfwL0Xj\nteB3sXieB4xoNkeiJjdusHjx04PYfWsljJ0uYKbfi3gzYipcrFy4jqazOmUIKDU1FV5eXvD09MTK\nlStbbfPf//4Xnp6e8Pf3R25ubrv6EvVrYpuQXZqDWT++h24Jvhj0+STkHnbFpz3/QtXm7/Dha/Tm\nT9TL3p7Bzk/CUfLRQUzl/4aNv5TA7ZMA9PpwBJam/g8Xb1/kOqLeUvobgEwmQ+/evXHw4EE4Ozsj\nODgY27dvh7e3t7xNSkoKEhMTkZKSguzsbLzyyivIyspSqC9A3wDU4frd6zh19W/sOZUN8eUMnL/7\nJ5rudId52XiM84zGa8/2h7+fAR3dQzpNYyOw9496rE3dh8zreyBzS4GpkSn8bUMwovdAjPYdgD7d\nvGFuZM51VK2l6Hun0nMC5+TkwMPDA0KhEAAQGxuL5OTkZm/iu3fvxvTp0wEAAwYMQFVVFSoqKlBY\nWPjEvuTxpE1SVNdXo6q+ChXVlbhQXo6L10tReKsEhbdKUXb3Cq435UPaJAN7zQfWdcHwsYjDh/7f\nYeIMB7i5cb0FRF/x+UD0WBNEj42BVBqD/zvRhB2Hz+LQqWx8kpWFZXbfAnYXYMLaorthbwgtesOt\nqwA97bvDrZsjPBwc4WTVDZbGljDjm9H80ypQugCUlZXBxeXfMTyBQIDs7OwntikrK8PVq1ef2PcB\noUgEFvcrmZWrK6yErgAgr24PHvu32snvwYO7Htzz733soy0f3NXyPqaNPmzz1m1lafH/j+bGv1W6\nCVI0MQ2QoQEypgEySCBFA6TsPzc0QMqrgdSwCk0GdTBotAJ7zxqot4aRxBFmUgGsGBc4WwzHiG4u\n6NejD/r16g4/PwZmZq3+egnhlKEhMHCAAQYO8AHgA2Ambt4Eck814fiFYuSVncfFSxewL68Md5ou\n4J5hOaQm5WDMrwNGNWAN68GTWsBQZgnDJgvwWXPwGWMYMkbg4d+bfJnhwwAGMDBgwMj/vf/zg3/x\n0PLDP98vNE8uNq23aHmvYnWLUWCNQFVREaqKrijyhM0oXQAUrbqqDuH0HDbs/vr++TUwYOS/S/l9\nDNNqm4d/dY9rc/+hB/f/8yjzSOuH2hg8tO2PrqvFeph/n7fNx/751wC8f168xvIXcRcjY5gZG8PM\nxBhmJkawNbOEnZk1bM3NYWlhgK5dATMzRV9MhGi+rl2B8JEGCB8pBCAEENnscZkMqK4Gbt8G7tRK\ncau2Brdra1BZdwdVdbWoa5DgnkQCiUyCxqZ//31wa2JZNDU1oYllIWtqAsuyuP+/JuCff5sv/3vf\nE7XSRKF+rT7Vvx9in8Ra6AZr4YOv9SyuHDmiUD+lC4CzszNKSkrkyyUlJRAIBI9tU1paCoFAgMbG\nxif2feDwkiXKRiSE6CAeD7C1vX+7/xZm88+NPMAsXapQO6WPAgoKCkJBQQGKioogkUiQlJSEqKio\nZm2ioqKwZcsWAEBWVhasra3h4OCgUF9CCCEdS+lvAIaGhkhMTERkZCRkMhni4uLg7e2N9evXAwDi\n4+MxZswYpKSkwMPDA2ZmZti0adNj+xJCCOk8dCIYIYToGLoWECGEkMeiAkAIIXqKCgAhhOgpKgCE\nEKKnqAAQQoieogJACCF6igoAIYToKSoAhBCip6gAEEKInqICQAgheooKACGE6CkqAIQQoqeoABBC\niJ6iAkAIIXqKCgAhhOgpKgAcEovFXEfoULq8fbq8bQBtn75QugDcvn0b4eHh6NWrFyIiIlBVVdVq\nu9TUVHh5ecHT0xMrV66U35+QkACBQIDAwEAEBgYiNTVV2ShaS9dfhLq8fbq8bQBtn75QugCsWLEC\n4eHhuHDhAkaMGIEVK1a0aCOTyTBv3jykpqYiPz8f27dvx9mzZwHcn7HmtddeQ25uLnJzczFq1Cjl\nt4IQQki7KV0Adu/ejenTpwMApk+fjl27drVok5OTAw8PDwiFQvD5fMTGxiI5OVn+OE33SAghHGKV\nZG1tLf+5qamp2fIDP//8Mztr1iz58g8//MDOmzePZVmWTUhIYF1dXVk/Pz925syZbGVlZYv+AOhG\nN7rRjW5K3BRhiMcIDw9HRUVFi/s/+uijZssMw4BhmBbtWrvvgblz52Lx4sUAgEWLFuH111/Hxo0b\nm7Vh6RsCIYR0mMcWgAMHDrT5mIODAyoqKtC9e3eUl5ejW7duLdo4OzujpKREvlxSUgKBQAAAzdrP\nmjUL48aNa3d4QgghylN6H0BUVBQ2b94MANi8eTOio6NbtAkKCkJBQQGKioogkUiQlJSEqKgoAEB5\nebm83W+//QZfX19loxBCCFECwyo5znL79m1MmjQJxcXFEAqF2LFjB6ytrXH16lXMnj0be/fuBQDs\n27cPCxYsgEwmQ1xcHBYuXAgAmDZtGk6dOgWGYdCzZ0+sX78eDg4O6tsyQgghj9euPb8c+d///sd6\neXmxPj4+7FtvvcV1nA6xatUqlmEY9tatW1xHUas33niD9fLyYv38/NiYmBi2qqqK60gq27dvH9u7\nd2/Ww8ODXbFiBddx1Kq4uJgViURsnz59WB8fH3bNmjVcR+oQUqmUDQgIYMeOHct1FLWrrKxkJ06c\nyHp5ebHe3t7ssWPH2myr8QUgLS2NHTlyJCuRSFiWZdnr169znEj9iouL2cjISFYoFOpcAdi/fz8r\nk8lYlmXZt99+m3377bc5TqQaqVTKuru7s4WFhaxEImH9/f3Z/Px8rmOpTXl5OZubm8uyLMvW1NSw\nvXr10qnte2D16tXsc889x44bN47rKGo3bdo0duPGjSzLsmxjY+NjP3Rp/KUg1q1bh4ULF4LP5wMA\n7O3tOU6kfq+99ho++eQTrmN0iPDwcBgY3H+ZDRgwAKWlpRwnUs2Tzm3Rdt27d0dAQAAAwNzcHN7e\n3rh69SrHqdSrtLQUKSkpmDVrls4daVhdXY309HTMnDkTAGBoaAgrK6s222t8ASgoKMDRo0cxcOBA\niEQi/N///R/XkdQqOTkZAoEAfn5+XEfpcN999x3GjBnDdQyVlJWVwcXFRb4sEAhQVlbGYaKOU1RU\nhNzcXAwYMIDrKGr16quv4tNPP5V/MNElhYWFsLe3xwsvvICnnnoKs2fPRl1dXZvtH3sYaGd53PkG\nUqkUlZWVyMrKwvHjxzFp0iRcvnyZg5TKe9z2LV++HPv375ffp42fSNravo8//lh+eO9HH30EIyMj\nPPfcc50dT60ed26LLqmtrcUzzzyDNWvWwNzcnOs4arNnzx5069YNgYGBOnk9IKlUipMnTyIxMRHB\nwcFYsGABVqxYgWXLlrXeoXNGpZQ3atQoViwWy5fd3d3ZmzdvcphIff766y+2W7durFAoZIVCIWto\naMi6urqy165d4zqaWm3atIkdNGgQe+/ePa6jqOzYsWNsZGSkfPnjjz/WuR3BEomEjYiIYD///HOu\no6jdwoULWYFAwAqFQrZ79+6sqakpO3XqVK5jqU15eTkrFArly+np6ezTTz/dZnuNLwBff/01u3jx\nYpZlWfb8+fOsi4sLx4k6ji7uBN63bx/bp08f9saNG1xHUYvGxkbWzc2NLSwsZBsaGnRuJ3BTUxM7\ndepUdsGCBVxH6XBisVgnjwIKCwtjz58/z7Isyy5ZsuSxR05qxBDQ48ycORMzZ86Er68vjIyMsGXL\nFq4jdRhdHF6YP38+JBIJwsPDAQAhISH46quvOE6lPENDQyQmJiIyMlJ+bou3tzfXsdQmMzMTW7du\nhZ+fHwIDAwEAy5cv19mr9eri39zatWsxZcoUSCQSuLu7Y9OmTW22VfpEMEIIIdpN93aDE0IIUQgV\nAEII0VNUAAghRE9RASCEED1FBYAQQvQUFQBCCNFT/w+tofxUhYEgGwAAAABJRU5ErkJggg==\n"
}
],
"prompt_number": 4
},
{
"cell_type": "code",
"collapsed": false,
"input": "plt.plot(xs, ps_uniform - ps_normal)\nplt.title(\"difference\")",
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "pyout",
"prompt_number": 5,
"text": "<matplotlib.text.Text at 0xb56796c>"
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEICAYAAAC0+DhzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XtYVOW+B/DvIKh4RU1AGXTkJqCCpkaW7c2WwNT0aHYh\nS6ncnqzMrF3b3B13ek4qure7m09ZnTIvbbeeUnGnkGihmRc08VKooKDc8YJ4RYbLe/54hZSBAWbW\nzFozfD/PMw/O8L5r/WYc1m+t97Z0QggBIiKi27ioHQAREWkPkwMREZlgciAiIhNMDkREZILJgYiI\nTDA5EBGRCSYHanGeeeYZzJ07F7t370ZwcHDt6ydPnsTAgQPRqVMnLFu2DDdv3sTYsWPh4eGBJ554\nQsWIiezPVe0AiOxNp9NBp9Nh+PDhOHHiRO3rS5YsQVRUFA4fPgwAWL16Nc6dO4eSkhK4uPA8iloW\nfuOpRapv7ufZs2cRGhp6x/OgoCCLEkNlZaVV8RGpjcmBnF5aWhruvvtudOrUCbGxsbh58yYAICUl\nBb6+vgCAESNGICUlBTNmzEDHjh0xadIk/M///A/WrVuHjh07YsWKFQCAL774AqGhoejatSseeugh\n5OTk1O7HxcUFH330EQIDA9G3b18AwLfffouBAweiS5cuuP/++3Hs2LHa8gaDAUuXLkV4eDg8PDwQ\nGxuL8vLy2t8nJCRg4MCB6Ny5MwICAvDdd98BAC5fvoypU6eiZ8+e0Ov1mDt3Lqqrq237IVLLI4ic\nWHl5uejVq5d47733RGVlpfj666+Fm5ubmDt3rkhJSRF6vb62bGRkpPj8889rn8+bN09Mnjy59vmm\nTZtEQECAOHHihKiqqhLvvPOOuO+++2p/r9PpRExMjLh06ZK4efOmOHTokPD09BSpqamiurparFy5\nUhgMBmE0GoUQQhgMBhERESEKCwtFSUmJCAkJEcuXLxdCCLF//37RuXNnsX37diGEEPn5+eLEiRNC\nCCHGjx8vpk+fLm7cuCHOnTsn7rnnHvHJJ5/Y7kOkFolXDuTU9u3bh8rKSrzyyito1aoVJk6ciKFD\nhzZYXtzW3CSEuOP58uXLMWfOHPTt2xcuLi6YM2cODh8+jNzc3Noyc+bMgYeHB9q0aYNPP/0Uzz//\nPIYOHQqdTocpU6agTZs22LdvX235mTNnwtvbG126dMHYsWNr+zs+//xzTJ06FVFRUQCAnj17om/f\nviguLkZiYiLeffdduLu7o3v37pg1axb+9a9/KfaZEQFsViInV1BQAB8fnzte6927d719DoDsrG7I\n2bNn8corr6BLly7o0qULunXrBgDIz8+vLVPTTFVTfunSpbXlu3Tpgry8PBQUFNSW8fb2rv23u7s7\nrl+/DgDIy8uDv79/vTFUVFSgR48etducPn06zp8/b+5jIGo2jlYip9ajR487Dt6APMAGBAQ0Wrdu\noujVqxfmzp2LJ598skl1evXqhbfeegt/+ctfmhm1TDKnTp2q9/U2bdrg4sWLHEFFNsVvFzm1++67\nD66urvjggw9QUVGBDRs24MCBAwDqH7FUt1npdtOnT8fChQuRnp4OQHYM/9///V+D+542bRqWL1+O\n1NRUCCFw/fp1bNmyBdeuXWuwTs0+p06dihUrVuD7779HdXU18vPzcfLkSfTo0QMxMTF47bXXcPXq\nVVRXV+P06dPYtWtX0z8UoiZgciCn5ubmhg0bNuDLL79Et27dsH79ekycOBHAb/Mdbnf787q/Hz9+\nPGbPno3Y2Fh07twZAwYMqB1BVLcuAAwePBifffYZZsyYga5duyIwMBCrVq1qsOnq9v0NHToUK1as\nwKuvvgoPDw9ERkbWjoxatWoVjEZj7aipxx57DEVFRVZ8SkSmdKKhxlciImqxrL5ySEpKQnBwMAID\nA7F48eJ6y8ycOROBgYEIDw9HWlpak+suXboULi4uKCkpsTZMIiJqBquSQ1VVFWbMmIGkpCSkp6dj\n7dq1OH78+B1ltm7dilOnTiEzMxOffvopXnjhhSbVzc3NRXJyMnr37m1NiEREZAGrkkNqaioCAgJg\nMBjg5uaG2NhYJCQk3FFm8+bNiIuLAwBERESgtLQURUVFjdZ97bXXsGTJEmvCIyIiC1k1lDU/P/+O\ncd16vR779+9vtEx+fj4KCgoarJuQkAC9Xo+wsLAG921uPDoREdWvqd3MVl05NPUA3Zw+77KyMixc\nuBDz589vtH7NDFZne7z99tuqx8D3x/fH9+d8j+aw6srBx8fnjqUDcnNzodfrzZbJy8uDXq9HRUVF\nvXVPnz6NM2fOIDw8vLb84MGDkZqaCk9PT2vCJSKiJrLqymHIkCHIzMzEmTNnYDQasW7dOowbN+6O\nMuPGjcOqVasAyHVuPDw84OXl1WDd/v37o7i4GNnZ2cjOzoZer8ehQ4eYGIiI7MiqKwdXV1csW7YM\nI0eORFVVFaZOnYqQkBB88sknAIDnn38eo0ePxtatWxEQEID27dvXLn3cUN26WmLfQmRkpNoh2BTf\nn2Pj+2sZHHYSnE6na3YbGhFRS9ac4yaXzyAiIhNMDkREZILJgYiITDA5EBGRCSYHIiIyweRAREQm\nmByIiMgEkwMREZlgciAiIhNMDkREZILJgYiITDA5EBGRCSYHckhGI3DkCHDunNqROI7iYuDoUaCi\nQu1IyBEwOZDD2bgR6N0bePJJICgImDGDBzxzjEbghReA4GAgNhYwGIDNm9WOirSOyYEcyldfAS+/\nLBNEejqQnS0fzz4LcAV3U0IAkycD+fnAmTPyM1u/Hpg+Xf4kagjv50AOIyMDuO8+ICUF6N//t9fL\nyoBhw4BZs4BnnlErOm368kvgww+Bn34C2rb97fUjR4CoKCA1FfDzUy08srPmHDeZHMghVFUBw4cD\nTz0lm5Hq+vlnYMwY4PRpoH17+8enRdeuAX37Ahs2ABERpr//xz+ATZtksnVhG0KLwJv9kNNZs0Ye\nwF58sf7fDx4sk8dnn9k3Li37xz+A3/++/sQAAK+8Aty8CaxbZ9+4yDHwyoE0z2iUnalffgn87ncN\nlzt4EJgwQV49tG5tt/A0yWgEevWSVwXBwQ2X275dJtz0dMDVqjvKkyPglQM5lZUrgcBA84kBAIYM\nkaOXNmywT1xalpAAhIaaTwyA7Hfw8QFWr7ZPXOQ4mBxI04QA3n0XePPNppV/5hnZBNXSLV8uRyQ1\nRqeTn+1773G0F92JyYE0bft22dwRGdm08hMmALt3A+fP2zQsTcvLAw4fBsaPb1r5mBjZDLVrl23j\nIsfC5ECa9uGHwMyZ8gy3KTp0AEaPbtlj+DdtAsaObXq/i04nP+MPPrBtXORYmBxIs4qK5Nnsk082\nr15sLPDNN7aJyRFs3CivoJrj6aeBHTuACxdsExM5HiYH0qy1a2XTSHPnLURFAQcOAFeu2CYuLbtw\nQY7aiolpXr2OHeU8EQ5rpRpMDqRZq1YBU6Y0v1779nLG9I4dysekdUlJMjm6uze/7pQp8jMnApgc\nSKOOHpVnwU3tiK5r9GggMVHRkBzC9u3Nv2qoERUF5OQAJ04oGxM5JiYH0qTVq2U7uKXLOoweDWzd\n2rKGZwohk0NUlGX1XV3l8iSc80AAkwNpUHW17G94+mnLtxEYKEfhnD6tXFxal5Ehk2lAgOXbePpp\n4J//bFlJlerH5ECac+CA7CDt18/ybeh0ckZ1Sxq7v2MH8OCDTR/2W5/wcKBVKzlPglo2JgfSnA0b\ngEcesX47LS05/PADMGKEddvQ6eRnzyVIiMmBNEUIOUdh4kTrt9WSkoMQwJ49wP33W7+tiRNb9jwR\nkpgcSFOOHZP3bhg0yPptBQcDV68CubnWb0vrcnJkX43BYP22hg6Vc0SOH7d+W+S4mBxIU2qalKxp\nN6+h08l7PPz0k/Xb0ro9e+TcDiU+NxcX+X+wcaP12yLHxeRAmmLJ0g/mRETIW2E6u717ZXJQCvsd\nyOrkkJSUhODgYAQGBmLx4sX1lpk5cyYCAwMRHh6OtLS0Ruu+8cYbCAkJQXh4OB555BFcvnzZ2jDJ\nAeTkAPn5yh7k7rlHjn5ydnv2yPtrK2X4cCArCygsVG6b5GCEFSorK4W/v7/Izs4WRqNRhIeHi/T0\n9DvKbNmyRYwaNUoIIcS+fftEREREo3W3bdsmqqqqhBBCzJ49W8yePdtk31aGThr08cdCPPWUstss\nLRWifXshKiqU3a6WlJUJ4e4uxI0bym73iSeE+N//VXabpK7mHDetunJITU1FQEAADAYD3NzcEBsb\ni4SEhDvKbN68GXFxcQCAiIgIlJaWoqioyGzd6OhouNyaGhsREYG8vDxrwiQHsWWLXPxNSZ07A76+\nwK+/KrtdLTl2TE76s2Q9JXMefhj49ltlt0mOw6q7xubn58PX17f2uV6vx/79+xstk5+fj4KCgkbr\nAsAXX3yBJxtYs3nevHm1/46MjESkpQvxkOrKyoCdO+UtQZVW07QUHq78trUgLQ24+27lt/vQQ8BL\nLwHl5UCbNspvn2wvJSUFKSkpFtW1Kjnomjg0Qlg4F3/BggVo3bo1Jk2aVO/vb08O5Nh27pQH765d\nld/20KGyU/qPf1R+21qQlqbM0N+67rpLzlLftQuIjlZ++2R7dU+a58+f3+S6VjUr+fj4IPe2QeS5\nubnQ6/Vmy+Tl5UGv1zda98svv8TWrVvx1VdfWRMiOYgtW+RiebYwaJBzLwdx6JBtkgPApqUWzZrO\njYqKCuHn5yeys7NFeXl5ox3Se/fure2QNlc3MTFRhIaGivPnzze4bytDJw2prhaiTx8hjhyxzfYv\nXxaiXTshKitts301VVTI93b5sm22f/iwEH5+8v+IHF9zjptWNSu5urpi2bJlGDlyJKqqqjB16lSE\nhITgk08+AQA8//zzGD16NLZu3YqAgAC0b98eK1asMFsXAF5++WUYjUZE37qWHTZsGD766CNrQiUN\nO3ECqKgABgywzfY7dQK8vIBTp4C+fW2zD7WcPAn07Cnfoy2EhQFGo9xPcLBt9kHapLuVTRyOTqez\nuC+DtOXvf5cH7uXLbbePCROASZOAxx6z3T7UsGYNsHkzsH697fYxfbocDfWnP9luH2QfzTlucoY0\nqc4WQ1jrCgsDjhyx7T7UYKvO6NuNGcN+h5aIyYFUdfkycPCg9UtNNyY8XN561NnYIzmMGCH/j0pL\nbbsf0hYmB1JVcrJcZrp9e9vuJzzc+a4chJCjsGydHNq3l8tpJCfbdj+kLUwOpCp7NCkBQJ8+QEmJ\nc5395ucDbm6ys93WxoyR/1fUcjA5kGqqq4HERPskBxcXORrKmZqW0tOtu5Vqc4wZI/+vqqvtsz9S\nH5MDqebnn4EuXQA/P/vsLyzM+ZJDaKh99tWnD9Ctm+x7oJaByYFUY8tZ0fVxtn6HX3+135UDwKal\nlobJgVSzdat9mpRq9Osnz7adhT2vHAAmh5aGk+BIFcXFcsZtcTHQurV99nnunNznxYvK3E5TTULI\nRQozM+UCefZQUQF4esqk1KOHffZJyuIkONK8xETgwQftlxgAoHt3mRTOnbPfPm2lsFCOVLJXYgDk\n/qKj5f8dOT8mB1KFvYaw3k6nA0JC5FpOjs7eTUo12LTUcjA5kN1VVADbtwOjRtl/3yEhwPHj9t+v\n0uzdGV1j1Chgxw65GB85NyYHsrvdu+VCbvaYvFVXcLBzJAe1rhw8PeXKtj/+aP99k30xOZDdqdGk\nVMNZrhzUSg4Am5ZaCiYHsju1k4Oj9zkIoV6zEsDk0FIwOZBdZWXJ9Y3uvlud/ffqBVy4AFy7ps7+\nlXD+vPzZvbs6+x80CLhyRd6Dg5wXkwPZ1ZYtslPTRaVvXqtWQFCQY189ZGTIdn+15mq4uMiZ7bx6\ncG5MDmRXajYp1XD0foeMDJng1MSmJefH5EB2c+UKsGePnEilJkfvd9BCcoiOBvbudezmOTKPyYHs\nJjFR3jSmUyd143D04axaSA4dOwIREXLOAzknJgeym40bgUceUTsKOcfCkTtTtZAcAGDsWGDzZrWj\nIFvhwntkF+XlgLc3cPKknEilpqtXZSxXr6rXMW6pqiqgQwc54srWt1ZtTE6OHHVWs84TaR8X3iPN\n2bFD3olN7cQAyCaRjh3lQc3R5ObKxfbUTgyAHBbs5wfs3Kl2JGQLTA5kFxs3AhMmqB3FbwID5XLX\njkYrTUo1Jk4EvvlG7SjIFpgcyOaqqoCEBG0lh4AAx+x30GJy2LhR/h+Tc2FyIJv76SfAxwcwGNSO\n5De8clBGQIBcQHHPHrUjIaUxOZDNbdigrasGgFcOSnr0UTYtOSMmB7Kpqipg3TrgiSfUjuROvHJQ\nTk2/Q3W12pGQkpgcyKZSUmSTUt++akdyp4AA4PRpucKpoygvBwoKtNU8B8ilwzt2BPbvVzsSUhKT\nA9nU2rXAk0+qHYWpmuGsBQVqR9J0Z88Cer025xRMmgSsWaN2FKQkJgeymfJy2d+gtSalGo7W75CV\nJecVaNHTTwPr1/P2oc6EyYFsJjERCA+XZ7ta5Gj9DllZgL+/2lHUz2CQCxomJqodCSmFyYFsZtUq\n2dygVY525XD6tHavHABg8mRg9Wq1oyClMDmQTRQVAT/8AMTGqh1JwxzxykHLyeGxx4DkZKCkRO1I\nSAlMDmQTX34phzh27Kh2JA1ztCsHrScHDw95hzh2TDsHrspKiquulmPxv/pKrvmvVVeuAD16yBvW\nqHXLzaYSQt4HIy8P6NxZ7WgatmsX8PzzQHq69j/Tlsiuq7ImJSUhODgYgYGBWLx4cb1lZs6cicDA\nQISHhyMtLa3RuiUlJYiOjkZQUBBiYmJQWlpqbZhkRykpQLt2wD33qB2JeZ06yeWvHWF11gsXgNat\ntZ0YAOCBB+R9ulNS1I6ErGVVcqiqqsKMGTOQlJSE9PR0rF27Fsfr3GJr69atOHXqFDIzM/Hpp5/i\nhRdeaLRufHw8oqOjkZGRgaioKMTHx1sTJtnZxx8D//mfjnHmGBDgGP0OWm9SqqHTAS+8IL8D5Nis\nSg6pqakICAiAwWCAm5sbYmNjkZCQcEeZzZs3Iy4uDgAQERGB0tJSFBUVma17e524uDhs2rTJmjDJ\njk6flmeNzzyjdiRN4yh3hXOU5ADIUUvJyY41wZBMuVpTOT8/H76+vrXP9Xo99teZQ19fmfz8fBQU\nFDRYt7i4GF5eXgAALy8vFBcX17v/QYPmAZBnK97ekejZMxI6HWofNb9T6nlDZVxcgLZtTR/u7r/9\n7NIF6NoV6NZN3qjFEc6qLfGPf8irhg4d1I6kafz8gOxstaNonCMlh06dZIJ47z1gyRK1o7ENIWSf\n1fnzwKVLwI0bwPXrdz5u3gQqK+Wjquq3f9/+qKqS27q9G6Dm33V/Nvd3AFBUlIKiohSL3qNVyUHX\nxCNcUzpAhBD1bk+n0zW4n5dfnnfHB3L7o+5rzX3enDrV1fKLUFoqf5aVyZ81jxs35Bfo4kU5zK+y\nUiaKHj2A3r3lw2CQ6w+FhQE9ezpm8rhwAfjnP4E6LYua1qePY0zcysoC7r1X7Sia7o035ATI2bPl\nCZGjEUIuV3LsmPzss7Pl4+xZ4Nw5+V1v0wbo3l2e+LVvb/po21YuddKqFeDqKh9t2/7271at5KPm\nVrW3/83ffiJqze+AyFsP6bnn5jf5M7AqOfj4+CA3N7f2eW5uLvR1psPWLZOXlwe9Xo+KigqT1318\nfADIq4WioiJ4e3ujsLAQng3cW/K556yJXj03b8pEUVAgv2xnz8rmmC1bgCNH5NlEWJjs0I2MBO6/\nX56Nad2yZXL4qre32pE0nZ+f/OPXuqwsbU8orMvXV34XPvgAmN/045FqLl8Gdu+WTaKpqfLvsH17\neWvbgAB5EvG738kTOS8veavWtm3Vjrr5mnXMFFaoqKgQfn5+Ijs7W5SXl4vw8HCRnp5+R5ktW7aI\nUaNGCSGE2Lt3r4iIiGi07htvvCHi4+OFEEIsWrRIzJ4922TfVoauaUVFQnz3nRDz5gkRGSlE+/ZC\n3HuvEPHxQpw8qXZ09btwQYhu3YTIzFQ7kuYpKBCie3e1o2icr68QWVlqR9E8mZnyO1FaqnYk9Tt+\nXIiFC4W45x4hOnQQ4g9/EGL+fCG2bRPi3Dm1o7ON5hw3rT7Cbt26VQQFBQl/f3+xcOFCIYQQy5cv\nF8uXL68t89JLLwl/f38RFhYmfv75Z7N1hRDi4sWLIioqSgQGBoro6Ghx6dIl08CdODnUVVYmRHKy\nEC++KESPHkL06yfE4sUyiWjF668L8fzzakfRfNXVQri7C3H1qtqRNOzmTSFatxaiokLtSJpvyhQh\n3npL7Sh+U1goxKJFQoSECOHjI8RLLwmxfbv8jFuC5hw3OQnOwVRXA3v3Al98IVc8HTECePFF+VOt\nfoqsLGDoUODoUXnvBkfTrx/wr3/JJgQtysiQM48dYVRVXbm5wMCBQFoa0KuXOjEIIZuLPvhA/nz0\nUeDZZ2UfjksLWyPCrpPgyL5cXGQfxOefy76K6Ghgxgx5cF6/Xp0bvb/yCvD6646ZGADZnqzlfgdH\nGqlUl6+v/H7+6U/q7P/HH4E//EHO2h41CsjJAT77DLjvvpaXGJqLH48D69QJmD4d+PVX4K9/Bd5/\nHwgOliOG7HXLxm++kZPI1PrjV4LWO6UdOTkAwJtvyg7ejRvtt8+9e+WJU1ycnHOTni6HWGt5rS+t\nYXJwAi4uwLhxwE8/AZ98Anz4obyU37zZtrfBLCgAXnpJLrLXurXt9mNrTA625e4ur3RnzAAamLKk\nmAMHZBNcbCzw+OPAyZMyObhaNS6zZWJycDIjRgB79gALFgD/9V/AsGHA998rvx+jEXjqKblUgiON\nv68Pk4PtPfCAHEY5aZJtmj4PHwb+4z+ACROAhx+W/TTTpmnzlqqOgsnBCel0wNix8g9m5kzZ3vrg\ng8rdAF4I2QneqZNMQI5O67OknSE5AMC8eXLS10svKXdF+8svcj7FqFHyxOjUKfndbNNGme23ZEwO\nTszFRZ6ppafL+zg/+qg8uzp2zPJtVlcDL78sRyatWSP/2B1dnz4yOdirn6Y5hHCe5NCqFfD118DB\ng3IAgzWf9/HjsukoKkpeHZ8+LQdGOOLENK1icmgB3NzkJXZmppxx/eCDwJNPyj/S5igtlQnm8GG5\nsJqzdO61by+vgoqK1I7E1MWLsr3cw0PtSJTRqRPw3XfyKjY2Vs5Mbo6ff5bfwd//Xq4icOqUTDTt\n2tkm3paMyaEFadsWePVV+Qc1aJC8HL//fnmvZ3O3zDAaZadz//5yuOqOHdq/r0BzabXfwVmuGm7X\nrRuwfbtcX2zAAHnf6YqKhsuXlsrv3/Dh8sp3+HB5pfeXvzjPCYoWcRJcC1ZZCSQkyOTwww8yYQwa\nJJexdneXiwUeOQIkJck/4nfekZfwzuipp4CRI4EpU9SO5E5r18ohoOvXqx2JbezcKfsi0tNlv8GA\nAfIq6fp12VR06JC8Uv3DH2SH9pgx7GS2RnOOm0wOBAC4elWODT9yRJ6tlpXJ1Sb79ZMdfc529lrX\n3Lmy+ebtt9WO5E4LFsj/G2e/31VWFrBtG3DihFwK291dfuf695cL3rm7qx2hc2jOcZOjfwmAvDyP\niZGPlsjPT5u3tszK0vZ9uJXi5ycndJJ2sM+BCNpdQiMrC/D3VzsKaomYHIjADmmiutjnQAQ5a7d9\ne9kJr5X2baNRNvddv87lH0gZXJWVqJlatZJLSp89q3Ykvzl7FtDrmRhIHUwORLdorWmJTUqkJiYH\nolu01inN5EBqYnIguoVXDkS/YXIguoXJgeg3TA5Et2ht6e7Tp5kcSD0cykp0y+XLcnTQlSvynhhq\nEkIubnj2rFzGhEgJHMpKZIHOneWibhcuqB2JXKq7VSsmBlIPkwPRbbTS78D+BlIbkwPRbfz8ZFu/\n2rimEqmNyYHoNlrplOaVA6mNyYHoNmxWIpKYHIhu4+/P5EAEMDkQ3UFLfQ5MDqQmznMguk1lpVy6\n+8oVoE0bdWKoWar72jXeL5mUxXkORBZydZUT4dRcuvvsWcDHh4mB1MXkQFSH2v0ObFIiLWByIKpD\n7X4HJgfSAiYHojrUHs7K5EBawORAVAeTAxGTA5EJtfscuFQ3aQGHshLVcfmyHC109ar9l+7mUt1k\nS3YZylpSUoLo6GgEBQUhJiYGpaWl9ZZLSkpCcHAwAgMDsXjx4kbrJycnY8iQIQgLC8OQIUPwww8/\nWBoikUU6d5ZzHM6ft/++L1yQw2mZGEhtFieH+Ph4REdHIyMjA1FRUYiPjzcpU1VVhRkzZiApKQnp\n6elYu3Ytjh8/brZ+9+7d8e233+Lo0aNYuXIlJk+ebGmIRBZTq9+BTUqkFRYnh82bNyMuLg4AEBcX\nh02bNpmUSU1NRUBAAAwGA9zc3BAbG4uEhASz9QcOHAhvb28AQGhoKMrKylBRUWFpmEQWUavfgUt1\nk1a4WlqxuLgYXl5eAAAvLy8UFxeblMnPz4evr2/tc71ej/379ze5/jfffIPBgwfDrYGpovPmzav9\nd2RkJCIjIy19O0R3UGuuw+nTTA6knJSUFKSkpFhU12xyiI6ORlFRkcnrCxYsuOO5TqeDrp6eu7qv\nCSEaLFf39V9//RVvvvkmkpOTG4zv9uRApCQ/P+Cnn+y/39OngeHD7b9fck51T5rnz5/f5Lpmk4O5\nA7OXlxeKiorg7e2NwsJCeHp6mpTx8fFBbm5u7fO8vDz4+Pg0Wj8vLw+PPPIIVq9ejT59+jT5zRAp\nxc8PWL3a/vvNygJutbYSqcriPodx48Zh5cqVAICVK1di/PjxJmWGDBmCzMxMnDlzBkajEevWrcO4\ncePM1i8tLcWYMWOwePFiDBs2zNLwiKyiVp8Dm5VIKyye51BSUoLHH38cOTk5MBgMWL9+PTw8PFBQ\nUIBp06Zhy5YtAIDExETMmjULVVVVmDp1KubMmWO2/jvvvIP4+HgEBgbW7is5ORl33XXXnYFzngPZ\nUFUV0K6dnPPQtq199llWJoewXr8OtGpln31Sy9Kc4yYnwRE1IDAQ+Pe/geBg++wvPR2YMAE4edI+\n+6OWh/cGG62BAAANmklEQVRzIFKAvec6sEmJtITJgagB9u534AQ40hImB6IG2HuuAyfAkZYwORA1\nQI1mJV45kFYwORA1gH0O1JJxtBJRA65cAXr0AK5ds/3S3dXVcuhsSYn8SWQLHK1EpIBOneSBup5l\nvxSXny/nODAxkFYwORCZYa+mJTYpkdYwORCZYa8RS7xvNGkNkwORGQEBwKlTtt9PZqackU2kFUwO\nRGYEBckDt61lZsp9EWkFkwORGUFBQEaG7feTkcErB9IWDmUlMqOkBDAY5OqsthrOWl0NdOggR0V1\n7GibfRABHMpKpJiuXYHWrYFz52y3j/x8wMODiYG0hcmBqBG2blpikxJpEZMDUSMCA23bKc3OaNIi\nJgeiRtjjyoHJgbSGyYGoEYGBbFailofJgagRtp7rwGYl0iIOZSVqxLVrgKen/Omi8OlUZaUcxnr5\nMtCmjbLbJqqLQ1mJFNShg1wxNS9P+W2fOSOXBWdiIK1hciBqAlt1SrMzmrSKyYGoCYKCgJMnld8u\nO6NJq5gciJogJAQ4flz57aanA6Ghym+XyFpMDkRNEBoqD+RKY3IgrWJyIGqCfv2UTw5CyG3266fs\ndomUwORA1AQ9ewJlZcDFi8pts7gYaNUK6N5duW0SKYXJgagJdDrlm5Z+/ZVNSqRdTA5ETaR0cmB/\nA2kZkwNREzE5UEvC5EDUREp3SrMzmrSMyYGoiUJDZT+BEoRgnwNpG5MDURP5+gJXrwKXLlm/rfPn\n5b2jvbys3xaRLTA5EDWRTqfcTOma/gadzvptEdkCkwNRM/Trp0zTEpuUSOssTg4lJSWIjo5GUFAQ\nYmJiUFpaWm+5pKQkBAcHIzAwEIsXL25y/ZycHHTo0AFLly61NEQixYWFAUePWr+do0eB8HDrt0Nk\nKxYnh/j4eERHRyMjIwNRUVGIj483KVNVVYUZM2YgKSkJ6enpWLt2LY7fuiZvrP5rr72GMWPGWBoe\nkU0MHAikpVm/nbQ0uS0irbI4OWzevBlxcXEAgLi4OGzatMmkTGpqKgICAmAwGODm5obY2FgkJCQ0\nWn/Tpk3w8/NDKK+7SWPCw+VZf3W15duorJTNSmFhysVFpDRXSysWFxfD69ZQCy8vLxQXF5uUyc/P\nh6+vb+1zvV6P/fv3m61/7do1LFmyBNu3b8ff/vY3szHMmzev9t+RkZGIjIy09O0QNUnXrvKRlQUE\nBFi2jYwMwMcH6NhR2diI6kpJSUFKSopFdc0mh+joaBQVFZm8vmDBgjue63Q66OoZdlH3NSFEg+Vq\nXp83bx5effVVtGvXrtF7nd6eHIjspaZpydLkwCYlspe6J83z589vcl2zySE5ObnB33l5eaGoqAje\n3t4oLCyEp6enSRkfHx/k5ubWPs/Ly4OPj4/Z+qmpqfjmm2/w5z//GaWlpXBxcYG7uztefPHFJr8p\nIlsaOBA4fBh47DHL6h8+zORA2mdxn8O4ceOwcuVKAMDKlSsxfvx4kzJDhgxBZmYmzpw5A6PRiHXr\n1mHcuHFm6+/atQvZ2dnIzs7GrFmz8NZbbzExkKYMGmRdpzSvHMgRWJwc3nzzTSQnJyMoKAjff/89\n3nzzTQBAQUFB7SgjV1dXLFu2DCNHjkRoaCieeOIJhISEmK1PpHVDhgAHDsglMJqruho4eFBug0jL\ndKKxhn2N0ul0jfZJENmKXg/s3An4+zev3vHjwJgxskObyN6ac9zkDGkiC9x7L3Br4F2z7N8v6xJp\nHZMDkQXuvRfYt6/59fbtY3Igx8DkQGQBJgdyduxzILLAjRvAXXcBFy8C7u5Nq3Ptmlyiu6QEaNPG\ntvER1Yd9DkQ21q4dMGBA8/od9uwB7r6biYEcA5MDkYWiooDt25tefvt2WYfIETA5EFnowQeBHTua\nXn7HDlmHyBGwz4HIQjdvAt27A3l5QOfO5stevAj06QNcuAC0bm2f+IjqYp8DkR20bQtERMjJcI35\n4Qdg+HAmBnIcTA5EVnjoIeDf/2683LffyrJEjoLNSkRWOHsWGDwYKCwE3NzqL1NeDvToARw7Ju/j\nQKQWNisR2Unv3vK+Dt9/33CZ5GSgf38mBnIsTA5EVnr8cWDduoZ/v369LEPkSNisRGSlggJ5ZZCd\nbTpqqaRErtx64oScHU2kJjYrEdlRz55ATAzw2Wemv/vkE2DsWCYGcjy8ciBSwC+/yNnPJ04AXbrI\n1y5cAEJCgN27gb591Y2PCOCVA5Hd9e8PxMYCzz4LVFTIxzPPyAcTAzkiXjkQKeTmTWDiRDm8VQg5\nimn9ei60R9rRnOMmkwORgoQAfvwR0OnkjGidTu2IiH7D5EBERCbY50BERFZhciAiIhNMDkREZILJ\ngYiITDA5EBGRCSYHIiIyweRAREQmmByIiMgEkwMREZlgciAiIhNMDkREZILJgYiITDA5EBGRCSYH\nIiIyweRAREQmmBw0KCUlRe0QbIrvz7Hx/bUMFieHkpISREdHIygoCDExMSgtLa23XFJSEoKDgxEY\nGIjFixc3qf7Ro0cxbNgw9O/fH2FhYSgvL7c0TIfk7F9Ovj/HxvfXMlicHOLj4xEdHY2MjAxERUUh\nPj7epExVVRVmzJiBpKQkpKenY+3atTh+/LjZ+pWVlZg8eTI+/fRT/PLLL9i5cyfc3NwsDZOIiCxg\ncXLYvHkz4uLiAABxcXHYtGmTSZnU1FQEBATAYDDAzc0NsbGxSEhIMFt/27ZtCAsLw4ABAwAAXbp0\ngYsLW7+IiOzJ4ntId+nSBZcuXQIACCHQtWvX2uc1vv76a3z33Xf47LPPAABr1qzB/v378eGHHzZY\n/7333sOhQ4dw7tw5nD9/HrGxsXjjjTdMA+ed24mImq2ph3xXc7+Mjo5GUVGRyesLFiy447lOp6v3\nYF33NSFEg+VqXq+srMTu3btx8OBBuLu7IyoqCoMHD8aIESNMtkVERLZhNjkkJyc3+DsvLy8UFRXB\n29sbhYWF8PT0NCnj4+OD3Nzc2ud5eXnw8fExW9/X1xe/+93v0LVrVwDA6NGjcejQIZPkQEREtmNx\nY/64ceOwcuVKAMDKlSsxfvx4kzJDhgxBZmYmzpw5A6PRiHXr1mHcuHFm68fExODYsWMoKytDZWUl\ndu7ciX79+lkaJhERWcDiPoeSkhI8/vjjyMnJgcFgwPr16+Hh4YGCggJMmzYNW7ZsAQAkJiZi1qxZ\nqKqqwtSpUzFnzhyz9QHgq6++wqJFi6DT6TBmzJh6R0IREZENCQf3wQcfiODgYNGvXz/x5z//We1w\nFPf3v/9d6HQ6cfHiRbVDUdTrr78ugoODRVhYmJgwYYIoLS1VOyRFJCYmir59+4qAgAARHx+vdjiK\nysnJEZGRkSI0NFT069dPvP/++2qHpLjKykoxcOBA8fDDD6sdiuIuXbokJk6cKIKDg0VISIjYu3ev\n2fIOnRy+//578eCDDwqj0SiEEOLcuXMqR6SsnJwcMXLkSGEwGJwuOWzbtk1UVVUJIYSYPXu2mD17\ntsoRWa+yslL4+/uL7OxsYTQaRXh4uEhPT1c7LMUUFhaKtLQ0IYQQV69eFUFBQU71/oQQYunSpWLS\npEli7NixaoeiuClTpojPP/9cCCFERUVFoydkDj2B4OOPP8acOXNqJ8l1795d5YiU9dprr2HJkiVq\nh2ET0dHRtfNXIiIikJeXp3JE1jM3r8cZeHt7Y+DAgQCADh06ICQkBAUFBSpHpZy8vDxs3boVf/zj\nH51uNOTly5fx448/4rnnngMAuLq6onPnzmbrOHRyyMzMxK5du3DvvfciMjISBw8eVDskxSQkJECv\n1yMsLEztUGzuiy++wOjRo9UOw2r5+fnw9fWtfa7X65Gfn69iRLZz5swZpKWlISIiQu1QFPPqq6/i\nb3/7m1NOus3Ozkb37t3x7LPP4u6778a0adNw48YNs3XMDmXVAnNzLSorK3Hp0iXs27cPBw4cwOOP\nP46srCwVorSMufe2aNEibNu2rfY1RzyTaej9LVy4EGPHjgUg32vr1q0xadIke4enuJYyMfPatWt4\n9NFH8f7776NDhw5qh6OIb7/9Fp6enhg0aJBTrq1UWVmJQ4cOYdmyZRg6dChmzZqF+Ph4/Pd//3fD\nlWzf0mU7Dz30kEhJSal97u/vLy5cuKBiRMo4duyY8PT0FAaDQRgMBuHq6ip69+4tiouL1Q5NUStW\nrBD33XefKCsrUzsURezdu1eMHDmy9vnChQudrlPaaDSKmJgY8e6776odiqLmzJkj9Hq9MBgMwtvb\nW7Rr105MnjxZ7bAUU1hYKAwGQ+3zH3/8UYwZM8ZsHYdODsuXLxd//etfhRBCnDx5Uvj6+qockW04\nY4d0YmKiCA0NFefPn1c7FMVUVFQIPz8/kZ2dLcrLy52uQ7q6ulpMnjxZzJo1S+1QbColJcUpRys9\n8MAD4uTJk0IIId5+++1GR3dqvlnJnOeeew7PPfccBgwYgNatW2PVqlVqh2QTzthc8fLLL8NoNCI6\nOhoAMGzYMHz00UcqR2UdV1dXLFu2DCNHjqyd1xMSEqJ2WIr56aefsGbNGoSFhWHQoEEAgEWLFuGh\nhx5SOTLlOePf3IcffoinnnoKRqMR/v7+WLFihdnyFk+CIyIi5+V83fJERGQ1JgciIjLB5EBERCaY\nHIiIyASTAxERmWByICIiE/8PT8HIi1f0Ly0AAAAASUVORK5CYII=\n"
}
],
"prompt_number": 5
},
{
"cell_type": "code",
"collapsed": false,
"input": "plt.plot(xs, (ps_uniform - ps_normal) / ps_normal)\nplt.title(\"relative difference\")",
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "pyout",
"prompt_number": 6,
"text": "<matplotlib.text.Text at 0xb6e324c>"
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEICAYAAABRSj9aAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlclFX7P/APCJkLohiLCjgBKuCKu/ak4zLikuaC5pdU\ncnuyfpVYltkmWipq5qNZtpmh5lqhJkSiiZWK5pK7aQqx44KoCAoM5/fHSRQZEGYG7lk+79drXjAz\n577nGpdrzlz3WWyEEAJERGSxbJUOgIiIqhYTPRGRhWOiJyKycEz0REQWjomeiMjCMdETEVk4JnpS\njFqtxsqVK/U6NikpCQ4ODqjq0cGJiYmwtbVFUVERAGDgwIFYs2ZN8fPvvPMOnJ2d0bhxYwBAZGQk\nPDw84ODggGPHjlVpbEQVxURPirGxsYGNjU2F2qpUKvzyyy/F9z09PXHz5s0KH28s0dHRGDt2LAD5\nYfPRRx/h7NmzSEtLAwBMnz4dn376KW7evIm2bdtWa2xEZWGipypRWFho1PPZ2NhUee+9spKSktCw\nYUM0bNgQACCEQFJSEvz9/fU6n1arNWZ4RMWY6MloVCoVFi5ciDZt2sDBwQFFRUWIj49H9+7d0aBB\nA7Rr1w579uzReeyFCxfQu3dvPPbYY3B2dsaYMWNw/fp1AMDYsWORlJSEwYMHw8HBAR9++GGJksrG\njRvRqVOnEudbsmQJnn76aQDAnTt3MH36dDRt2hRubm544YUXcPv2bZ1xFBUVYfr06XB2doa3tzei\noqJKPH+33LRr1y7069cPaWlpcHBwQHBwMOrVqwetVou2bduiWbNmAIC0tDSMGDECLi4u8PLywscf\nf1x8rrCwMAQFBWHs2LFwdHREREQErl+/jokTJ6Jx48Zwd3fHu+++W1w2+uabb/Cf//wHr7/+Opyc\nnODl5YWYmJji82VlZWH8+PFo0qQJnJycMGzYsOLntm/fjnbt2qFBgwZ44okncOLEiQr9nZKFEERG\n0rRpUxEQECBSUlLE7du3RUpKimjYsKH46aefhBBCxMbGioYNG4orV64IIYRQq9Vi5cqVQggh/v77\nb7Fz506Rn58vLl++LHr06CFCQ0OLz61SqcSuXbuK7yckJAgbGxuh1WrFrVu3hIODgzh//nzx8x07\ndhQbN24UQggRGhoqnn76aXHt2jVx8+ZNMXjwYDFz5kyd72HFihXC19dXpKSkiKysLKFWq4Wtra3Q\narWlYo6LixPu7u4ljrexsREXLlwQQgih1WpF+/btxfvvvy8KCgrExYsXhZeXl/j555+FEELMmjVL\n2Nvbi61btwohhMjLyxNDhw4VU6ZMEbm5ueLSpUuic+fO4vPPPxdCCLFq1Sphb28vvvrqK1FUVCRW\nrFghGjduXPzaAwcOFKNHjxbZ2dmioKBA/Prrr0IIIY4cOSJcXFzEwYMHRVFRkYiIiBAqlUrcuXOn\ngn+zZO6Y6MloVCqVWLVqVfH98PBwMXbs2BJtAgMDRUREhBCiZNJ8UGRkpAgICChx7rISvRBCjBkz\nRsyZM0cIIcS5c+eEg4ODyMvLE0VFRaJOnTrFyVcIIfbt2ycef/xxna/bq1ev4sQqhBA7duwo8Tr3\nx7x79+5yE318fLzw9PQs8fy8efPE+PHjhRAy0ffs2bP4uYyMDFGzZk2Rl5dX/Ni6detEr169hBAy\n0fv4+BQ/d+vWLWFjYyMyMzNFWlqasLW1FdnZ2aXe05QpU8S7775b4rEWLVqIPXv26PwzIMtjp/Q3\nCrIsHh4exb//888/2Lx5M3788cfixwoLC9G7d+9Sx2VmZmLq1Kn4/fffcfPmTRQVFcHJyanCrxsc\nHIzXXnsN7777LtatW4dhw4bh0UcfxaVLl5Cbm4sOHToUtxVCFJdDHpSenl7iPXh6elY4hgf9888/\nSEtLQ4MGDYof02q16NGjR/F9d3f3Eu0LCgrQqFGj4seKiopKxODm5lb8e+3atQEAOTk5uHLlCpyc\nnODo6KgzjtWrV5coGxUUFCA9PV3v90bmhYmejOr+UTCenp4YO3Ysvvjii4ce99Zbb6FGjRo4efIk\n6tevjy1btuDll1/WeV5d+vbti8uXL+PYsWPYsGED/ve//wEAHnvsMdSqVQunT58ukUDL0qhRIyQl\nJRXfv//3yvL09MTjjz+Oc+fO6Xz+wVFHHh4eqFmzJq5evQpb28pdPvPw8EBWVhauX79eKtl7enri\n7bffxltvvVX5N0EWgRdjqcqMGTMGP/74I3bs2AGtVovbt28jLi4Oqamppdrm5OSgTp06qFevHlJT\nU7Fo0aISz7u6uuLChQtlvpa9vT1GjhyJ6dOn49q1a9BoNAAAW1tbTJ48GaGhobh8+TIAIDU1FTt2\n7NB5nlGjRmHZsmVITU3FtWvXEB4eru/bR+fOneHg4ICFCxciLy8PWq0WJ0+exKFDhwCg1CiiRo0a\noV+/fnj11VeLv9VcuHABv/7660Nfq1GjRhgwYABefPFFZGdno6CgoPi4yZMn47PPPsPBgwchhMCt\nW7cQFRWFnJwcvd8bmRcmeqoy7u7u2Lp1K+bNmwcXFxd4enpi8eLFOodJzpo1C0eOHIGjoyMGDx6M\nESNGlOjtzpw5Ex988AEaNGiAjz76CEDpXn5wcDB27dqFkSNHlugRL1iwAD4+PujatSscHR2h0WjK\n7GVPnjwZgYGBaNu2LTp27Fgqjgc9+Nz9921tbbF9+3b8+eef8PLygrOzM/773//ixo0bxW0fPH71\n6tXIz8+Hv78/nJycMHLkSGRkZJTZ/v77a9asgb29PXx9feHq6oply5YBADp06IAvv/wSL730Epyc\nnNCsWTOsXr26zPdElsdG6PpfR0REFsPgHn1MTAx8fX3RrFkzLFiwoNTz3377Ldq2bYs2bdrgiSee\nwPHjxw19SSIiqgSDevRarRYtWrTAzp070aRJE3Tq1Anr16+Hn59fcZv9+/fD398fjo6OiImJQVhY\nGOLj440SPBERPZxBPfqDBw/Cx8cHKpUK9vb2GD16NLZu3VqiTbdu3YpHAXTp0gUpKSmGvCQREVWS\nQcMrU1NTS4w5dnd3x4EDB8psv3LlSgwcOLDU49W9MBURkaWoSFHGoB59ZRL07t278fXXX+us4wMy\nWEu9zZo1S/EY+N74/vj+LO9WUQb16Js0aYLk5OTi+8nJySVm+t11/PhxTJ48GTExMSVmCRIRUdUz\nqEffsWNHnD9/HomJicjPz8fGjRsxZMiQEm2SkpIwfPhwrF27Fj4+PgYFS0RElWdQj97Ozg7Lly9H\nYGAgtFotJk6cCD8/P3z++ecAgOeffx5z5szBtWvX8MILLwCQMxgPHjxoeORmRK1WKx1ClbHk9wbw\n/Zk7S39/FWUSE6ZMcVMJIiJTV9HcySUQiIgsHBM9EZGFY6InIrJwTPRERBaOiZ6IyMIx0RMRWTgm\neiIiC8dET0Rk4ZjoiYgsHBM9EZGFY6InIrJwTPRERBaOiZ6IyMIx0RMRWTgmeiIiC8dET0Rk4Zjo\niYgsHBM9EZGFM2jPWCIhgIsXgfPngcREIDkZuHEDyMmRNyEAe3vgkUeAWrWAxx6TN2dnwNUVUKkA\nT0/5PFWeEMDVq/LPPSkJSEkBrl0Drl8HsrOBmzeBoiLAxkbe6ta993fQqBHg6wu0aCEfJ8tl8J6x\nMTExCA0NhVarxaRJkzBjxoxSbV555RX89NNPqF27Nr755hsEBASUDIJ7xpqNoiLg6FEgOhrYvRs4\ncgRwcJAJo2lTmbTr15eJo04doEYNID8fKCgAbt2SSenyZeDKFSA9XX44pKUBLi7A448DzZoB/v7y\n5ucnz2fL753IywP++gs4cwY4fVr+PHMGSEgAHn0U8PCQN3d3wMlJ/h3Ury//bmxt5QeCEPLD98oV\neUtNBc6eBc6dk3/OPXoAPXsCAwbIc5Dpq2juNCjRa7VatGjRAjt37kSTJk3QqVMnrF+/Hn5+fsVt\noqOjsXz5ckRHR+PAgQOYOnUq4uPj9QqWlHP+PLByJbB6NVCvHjBwINC3L9Cxo0zShigslD3RhAT5\nOqdP37tlZ8sPET+/kh8AXl6AnQV+H71xo2Qyv/szLQ3w9r73/u//czC0N67VAidOAL/+Cvzyi/wA\n/89/gOBgYORIftsyZdWS6Pfv34/Zs2cjJiYGABAeHg4AePPNN4vbTJkyBb169cIzzzwDAPD19cWe\nPXvg6upa6WCp+v3+OzB7NnD8OBASAkyYIBNvdbl+XfY6T58GTp2615NNTy/d+/f3l4+ZQ2K6cqV0\nMj99WpZdfH1LJ3Rv7+r7YLt5E/jxR+Drr2VML70EvPAC0KBB9bw+VVxFc6dB/3RSU1Ph4eFRfN/d\n3R0HDhx4aJuUlJQSiR4AwsLCin9Xq9VQq9WGhEYG+vNP4PXXgQsXgHfeAaKilEmgjo5Aly7ydr9b\nt0qWMtatkz+TkmTdv3lzWQpSqe79VKnk+aqDVit74Rcvym8qd38mJMhSSX5+yWTev7/plKocHGRv\nPjhYfsB/9JGs48+YIZN+zZrKxmfN4uLiEBcXV+njDEr0NjY2FWr34CeOruPuT/SknFu3gFmzZInm\n/fdlD97eXumoSqtTB2jfXt7ud+eOLP+cOyfr/xcvynLE3SRraysvAru43Pvp7HzvmkKdOvd+t7O7\nV9u++0+4sFDWuW/elLe7v1+6BGRkAJmZ8mdWljyvl5f8oPHyAvr0kT99fOSF0Ar+91FUmzbAN9/I\nD9Q33gC+/BL46itZ2qHq92AnePbs2RU6zqBE36RJEyQnJxffT05Ohru7e7ltUlJS0KRJE0NelqrI\n/v2yF/fEE8DJk4bX3pVQsybQqpW8PUgIWQq6dOneLTNTXhy+elV+G7h1SybvW7dkrxy4N2LFxkZe\nXHZwkB8GDg7y5uoKtG4tf7q5yZ/OzpZ1DcHPT5ZzfvgBeOYZYNQoYMEC8yiTEQBhgIKCAuHl5SUS\nEhLEnTt3RNu2bcXp06dLtImKihIDBgwQQgixf/9+0aVLl1LnMTAMMlBRkRCLFgnh4iLEli1KR0Om\n7upVIZ5+WohOnYRISFA6GutW0dxpUJ/Dzs4Oy5cvR2BgILRaLSZOnAg/Pz98/vnnAIDnn38eAwcO\nRHR0NHx8fFCnTh2sWrXKCB9PZCx37gCTJ8uv5gcPyiGSROVxcgIiI4ElS4Bu3eTvXbsqHRWVx+Bx\n9EYJgqNuFHH9OvD00/I/7tq1QO3aSkdE5iYqCnjuOeCLL4Bhw5SOxvpUNHdyKoqVunYN0GjkiI/N\nm5nkST+DBgE//wy8+CKwfr3S0VBZLOhyEVVUdrYcAdKzpxw6Zw6jP8h0tW8PxMYC/frJmdPPPqt0\nRPQglm6szO3bcsx269bAsmVM8mQ8p08DvXvLiVYDByodjXWolpmxxsJEXz2KioDRo+Uwww0b5FBB\nImOKjweGDAG2bpUXaqlqsUZPpUyfLseOr1nDJE9Vo2tXOcFqxAg5L4FMAxO9lVi3Tk54iYyUqx0S\nVZWBA4FXXwWGDgVyc5WOhgCWbqzCyZNAr17Arl1ySjtRVRNCLoInhPwGSVWDpRsCIJe9HTECWLyY\nSZ6qj40N8NlnwOHDTPSmgD16Czd2rBwj/+9kZaJqdeyY3Ldg/365mBsZV7UsU0ym7Ycf5CiIP/9U\nOhKyVm3bytVQ/+//gL17uQiaUtijt1CXLsn/ZN9/D3TvrnQ0ZM2EAAYPBjp1kkmfjIfj6K2YEMDw\n4XKziH83/SJSVEoKEBAA7Nkjl90g4+DFWCv2/fdy440K7klAVOXc3YE5c4BJk+6t80/Vh4newuTk\nANOmAStWcMs3Mi3PPy93+FqxQulIrA9LNxbmjTfkVnarVysdCVFpZ8/KbQhPnJDbKZJhWKO3QqdO\nAWq1nCD1wN7rRCbjjTeAK1fk4mdkGCZ6KyOEXHp4+HDgpZeUjoaobDduyIEC27cDHTooHY1548VY\nKxMdLUs2U6YoHQlR+erVA95/H5g6VXZQqOrpneizsrKg0WjQvHlz9OvXD9nZ2aXaJCcno1evXmjZ\nsiVatWqFZcuWGRQs6abVAjNmAAsWAHacAkdmYPx44NYtYNMmpSOxDnon+vDwcGg0Gpw7dw59+vRB\nuI4B2/b29liyZAlOnTqF+Ph4fPLJJzhz5oxBAVNpERFy39ennlI6EqKKqVFD7m729ttAQYHS0Vg+\nvRP9tm3bEBISAgAICQnBli1bSrVxc3NDu3btAAB169aFn58f0tLS9H1J0iE3F3jvPWDhQu4WReal\nVy/g8ceBVauUjsTy6f1FPzMzE67/Du1wdXVFZmZmue0TExNx9OhRdOnSRefzYWFhxb+r1Wqo1Wp9\nQ7MqH38sN3vo2lXpSIgq74MPgKAgYNw47pNQEXFxcYiLi6v0ceWOutFoNMjIyCj1+Ny5cxESEoJr\n164VP+bk5ISsrCyd58nJyYFarcY777yDoUOHlg6Co270kpMDeHkBcXGcVk7ma8gQOWJs6lSlIzE/\nRlm9MjY2tsznXF1dkZGRATc3N6Snp8PFxUVnu4KCAowYMQJjxozRmeRJf598Ir/+MsmTOXv/fblh\n/aRJQJ06SkdjmfSu0Q8ZMgQREREAgIiICJ1JXAiBiRMnwt/fH6GhofpHSaXk5MiLWe++q3QkRIZp\n2xbo0QP49FOlI7Fcek+YysrKwqhRo5CUlASVSoVNmzahfv36SEtLw+TJkxEVFYXff/8dPXr0QJs2\nbWDz75XC+fPno3///iWDYOmm0j78EDhwANi8WelIiAx37BgwYABw8SJr9ZXBmbEWLDdX1uZjY4HW\nrZWOhsg4nnpKrlv//PNKR2I+ODPWgq1aJUfZMMmTJZk5U076KyxUOhLLw0RvZrRaWZt/4w2lIyEy\nrieeADw8OFu2KjDRm5kffpArU3J7QLJEM2cC8+dzDRxjY6I3I0IAixYBr7+udCREVSMwUG5OsmOH\n0pFYFiZ6M/Lbb0B2tpxgQmSJbGyA0FDgf/9TOhLLwlE3ZmTIEGDQII5KIMt2+zagUgG7dwN+fkpH\nY9o4vNLCXLggR9okJQG1aikdDVHVCguT+yt89pnSkZg2JnoLM326rF0uXKh0JERVLzMT8PUF/v4b\naNhQ6WhMFxO9BcnNBZo2BQ4elMu6ElmD8eOBZs2At95SOhLTxQlTFmT9elm2YZIna/LKK7J0o9Uq\nHYn5Y6I3cUIAy5dzw2+yPgEBQOPGcj9kMgwTvYnbv1/uranRKB0JUfWbMoUXZI2BNXoTFxwMdO4s\nxxYTWZvcXMDTE/jjD5YudeHFWAuQkSHHESckAPXrKx0NkTKmTZNDiufNUzoS08OLsRbgm2+AESOY\n5Mm6Pf888PXXQH6+0pGYLyZ6EyUE8NVXwOTJSkdCpCxfX/nNNjJS6UjMFxO9idqzR35d7dxZ6UiI\nlPfCC7woawgmehP15ZeyN//vDoxEVm3oUOD0aeDcOaUjMU+8GGuCsrLkVoEXLwJOTkpHQ2QaXnsN\nqFmTF2XvV+UXY7OysqDRaNC8eXP069cP2dnZZbbVarUICAjA4MGD9X05q/Ltt8DAgUzyRPcbPx6I\niOBMWX3onejDw8Oh0Whw7tw59OnTB+Hh4WW2Xbp0Kfz9/WHDOsRDCSHLNpMmKR0JkWlp1Qpo0oSb\nkuhD70S/bds2hISEAABCQkKwZcsWne1SUlIQHR2NSZMmsTxTAX/8IWfCqtVKR0JkesaPB1atUjoK\n82On74GZmZlwdXUFALi6uiIzM1Nnu2nTpmHRokW4ceNGuecLCwsr/l2tVkNtpZlu5Upg4kS5JDER\nlTR6tNxXNivLOkubcXFxiIuLq/Rx5SZ6jUaDjIyMUo/PnTu3xH0bGxudZZnt27fDxcUFAQEBDw3u\n/kRvre7cAb77Djh6VOlIiExTgwbAgAHAunXWudDfg53g2bNnV+i4chN9bGxsmc+5uroiIyMDbm5u\nSE9Ph4uLS6k2+/btw7Zt2xAdHY3bt2/jxo0bGDduHFavXl2h4KxNdDTQpo1c24OIdJswAZgxwzoT\nvb70LhAMGTIEERERAICIiAgMHTq0VJt58+YhOTkZCQkJ2LBhA3r37s0kX441a4AxY5SOgsi09e4N\nXLkCHDumdCTmQ+9E/+abbyI2NhbNmzfHL7/8gjfffBMAkJaWhkGDBuk8hqNuypaVBezaJde2IaKy\n1agBhITwomxlcMKUifjiCyA2Fti8WelIiEzfhQtAt25ASgrwyCNKR6Mcrl5pZtasAcaOVToKIvPg\n7S0XO4uJUToS88BEbwISEoAzZ4D+/ZWOhMh8jBkDrF2rdBTmgYneBKxbB4waZd1fQYkqa+RI4Oef\ngXJWX6F/MdErTAjZK2HZhqhyGjQA+vQBvv9e6UhMHxO9wg4fljvndO2qdCRE5mfsWJZvKoKJXmFr\n18paI0eeElXewIHA8eNAUpLSkZg2JnoFFRYCGzZwkhSRvmrWBIKCgPXrlY7EtDHRKyg2FlCpgGbN\nlI6EyHyNGSOHJ1v5VJxyMdEr6G7Zhoj098QTQE4Ol0QoD2fGKuTmTcDDAzh/HnB2VjoaIvP29tty\n9dcPP1Q6kurFmbEmLjISePJJJnkiYxgzRs5H4TaDujHRK4RlGyLj8fMDGjcGdu9WOhLTxESvgLQ0\nuWXgkCFKR0JkObgkQtmY6BWwYQMwbBhQq5bSkRBZjtGjga1bgdxcpSMxPUz0CuAGI0TG5+YGdOki\nkz2VxERfzU6eBC5fBnr2VDoSIsvD8o1uTPTV7NtvgeBguUsOERnX0KHA3r2yM0X3MNFXo6Iimei5\nUiVR1ahbF3jqKWDjRqUjMS1M9NXo11/l0qqtWysdCZHlYvmmNL0TfVZWFjQaDZo3b45+/fohu4zV\n/7OzsxEUFAQ/Pz/4+/sjPj5e72DNHcfOE1W9vn2BxEQ565wkvRN9eHg4NBoNzp07hz59+iA8PFxn\nu6lTp2LgwIE4c+YMjh8/Dj8/P72DNWd5ecAPP8j6PBFVHTs7OdTy22+VjsR06L3Wja+vL/bs2QNX\nV1dkZGRArVbj7NmzJdpcv34dAQEBuHjxYvlBWMFaN5s3A59/DuzcqXQkRJbv0CGZ7M+ft+y9Hiqa\nO+30fYHMzEy4uroCAFxdXZGZmVmqTUJCApydnTF+/HgcO3YMHTp0wNKlS1G7du1SbcPCwop/V6vV\nUKvV+oZmkli2Iao+HTrInv2BA5a1e1tcXBzi4uIqfVy5PXqNRoOMjIxSj8+dOxchISG4du1a8WNO\nTk7Iysoq0e7QoUPo1q0b9u3bh06dOiE0NBT16tXDnDlzSgZh4T36K1cAb28gORmoV0/paIiswwcf\nABkZwPLlSkdSdYzSo4+NjS3zubslGzc3N6Snp8PFxaVUG3d3d7i7u6NTp04AgKCgoDJr+ZZs0ya5\n5RmTPFH1CQ6WM2WXLAHs7ZWORll6X4wdMmQIIiIiAAAREREYOnRoqTZubm7w8PDAuXPnAAA7d+5E\ny5Yt9X1Js8WyDVH18/ICWrQAfv5Z6UiUp/fF2KysLIwaNQpJSUlQqVTYtGkT6tevj7S0NEyePBlR\nUVEAgGPHjmHSpEnIz8+Ht7c3Vq1aBUdHx5JBWHDp5u+/ge7dgdRU9iqIqtuKFcCePXIhQUtU0dzJ\nHaaq2KxZQHY2sHSp0pEQWZ+rV2XP3lKvj3GHKRMghFypctw4pSMhsk4NGwJqtdzRzZox0VehvXuB\nRx8F2rdXOhIi68UlEZjoq9Td3rwlT9ggMnVPPQUcPix3drNWTPRV5PZt4LvvgGefVToSIutWq5bc\n0W39eqUjUQ4TfRXZvh1o1w7w8FA6EiKy9vINE30VWb2a684TmYqePeUM9ZMnlY5EGUz0VeDyZbn2\n/IgRSkdCRABgaytnylrripZM9FVgwwZ5AcjBQelIiOiuu+UbrVbpSKofE30VWLOGZRsiU9O6NeDm\nZp1LhTPRG9nJk3K5gz59lI6EiB40YQLw9ddKR1H9uASCkU2bBtSuDcydq3QkRPSg7GxApQIuXJCz\nZs0dl0BQwJ07sgY4YYLSkRCRLvXrA4MGAevWKR1J9WKiN6KtW2Ud0Ntb6UiIqCzWWL5hojeilSuB\nSZOUjoKIytOrF3DtGnD0qNKRVB8meiP55x+5IfGwYUpHQkTlsbUFxo+3rl49L8YaSViYnHlnyftT\nElmKf/6RG4inpMgVZs0VL8ZWI61W9g5YtiEyD02byuXDv/9e6UiqBxO9EcTGAs7OchEzIjIPL74I\nfPqp0lFUDyZ6I/j0U2DKFKWjIKLKeOopICkJOHZM6Uiqnt6JPisrCxqNBs2bN0e/fv2QnZ2ts938\n+fPRsmVLtG7dGsHBwbhz547ewZqihAS5k1RwsNKREFFl2NkB//2v3EDc0umd6MPDw6HRaHDu3Dn0\n6dMH4eHhpdokJibiyy+/xJEjR3DixAlotVpssLDt2D/7DHjuOaBOHaUjIaLKmjQJ2LgRuHFD6Uiq\nlt6Jftu2bQgJCQEAhISEYMuWLaXa1KtXD/b29sjNzUVhYSFyc3PRpEkT/aM1MXl58iLsCy8oHQkR\n6aNRI0CjkQsRWjI7fQ/MzMyEq6srAMDV1RWZmZml2jg5OeG1116Dp6cnatWqhcDAQPTt21fn+cLC\nwop/V6vVUKvV+oZWbTZuBDp1Anx8lI6EiPT14ov3bqa+v3NcXBzi4uIqfVy54+g1Gg0yMjJKPT53\n7lyEhITg2rVrxY85OTkhKyurRLsLFy5g8ODB+O233+Do6IiRI0ciKCgIzz6wkao5jqMXQib5OXOA\ngQOVjoaI9CUEEBAAzJ8PDBigdDSVU9HcWW6PPjY2tsznXF1dkZGRATc3N6Snp8PFxaVUm0OHDqF7\n9+5o+O8yccOHD8e+fftKJXpz9PvvwPXrQP/+SkdCRIawsQFeew1YvNj8En1F6V2jHzJkCCIiIgAA\nERERGDp0aKk2vr6+iI+PR15eHoQQ2LlzJ/z9/fWP1oQsWABMny6nUxOReXvmGeDsWeDPP5WOpGro\nvQRCVlYWRo0ahaSkJKhUKmzatAn169dHWloaJk+ejKioKADAwoULERERAVtbW7Rv3x5fffUV7O3t\nSwZhZqUy4vy/AAAQpElEQVSbkyflBZyEBPOePk1E9yxcCJw4YV4XZiuaO7nWjR7GjQP8/ICZM5WO\nhIiMJTsb8PKSE6g8PJSOpmKY6KtIUpJc6uDCBaBBA6WjISJjevVV+fOjj5SNo6KY6KvItGlAjRrA\nhx8qHQkRGVtqqtw86MwZ4N/R4yaNib4KpKUBrVrJGn3jxkpHQ0RV4ZVXAHt7OQrH1DHRV4H/9/+A\nWrXYmyeyZHc7dKdPA25uSkdTPiZ6I0tMlBsVnD0rlyQmIss1dar8uXSpsnE8DBO9kT33nLwS//77\nSkdCRFXt0iXA3x/Yvx9o1kzpaMrGRG9Ef/wBDBkie/OOjkpHQ0TVITwcOHgQ+OEHpSMpGxO9kRQV\nAU88AUyeDEyYoHQ0RFRdbt8GfH2B1auBHj2UjkY37hlrJN9+CxQWytINEVmPRx+VS5289BKQn690\nNIZhj74cV64AbdrIr25duyodDRFVNyHkloPduwNvv610NKWxdGMEzzwjL8ByOCWR9UpKAtq3lyvW\n+voqHU1JLN0YaNMm4PhxjrIhsnaenkBYGDB2LGCuW16zR6/D33/LC7A//gh07qx0NESkNCGAoUMB\nlcq0xtazR6+nnBxg2DD5Cc4kT0SA3Jzkm2+AbduA775TOprKY4/+PlotMGoUUL8+8NVXpr9/JBFV\nr8OH5a5yW7bIb/1Ks+ge/aVLcpU5YyoqkuPkb9wAPvmESZ6ISuvQAVi7Fhg+HDh1yrjnvnNHrppZ\nVGTc8wJmmOh//FFe+W7dGnjvPVk7M1RhIfD883LHqC1buGsUEZUtMBBYskTuMnf4sHHOeeqUXHKh\nb1+gXz+goMA4573LrBJ9bi4waRIQFQWcOyfHtxt6YeTmTbm8wT//ANu3A3XqGCdWIrJcwcHAp5/K\nzcR37DDsXOnp8sPjvfeA5GS538UnnxgnzrvMKtGvWwd06QJ06wY89phM+IsWATEx+p3v9Gk5EcLD\nQ56rXj3jxktElmvoUOD77+Ws+Q8+0K/kUlgIjB4tO7AhIYCtLTBnDvDxx8Yt4eid6Ddv3oyWLVui\nRo0aOHLkSJntYmJi4Ovri2bNmmHBggX6vhwAYMOGkuvNNG0KrF8v/6CTkip+nqIi+WncsycQGgp8\n9pncaICIqDKefBI4dAj4+WdZyrl4sXLHv/suULOm/HlX586yfHzokBEDFXo6c+aM+Ouvv4RarRaH\nDx/W2aawsFB4e3uLhIQEkZ+fL9q2bStOnz5dql1Fwrh5U4i6dYW4caP0cwsXCtG5sxC5uQ+P+/hx\nIbp1E6JrVyHOnn14eyKihykokHmoYUMhPvpIiMLChx+zdq0Qnp5CXLpU+rnp04WYPfvh56hoCte7\nR+/r64vmzZuX2+bgwYPw8fGBSqWCvb09Ro8eja1bt+r1er/9Jq94OziUfm76dLlm9NChQF6e7uMz\nM+VmAn36yK9Ie/cCLVroFQoRUQl2dsDrr8v163/8EWjXDvjpp7IHi0RFyf2no6J0b2TUrx8QG2vE\n+Ix3qtJSU1Ph4eFRfN/d3R0HDhzQ2TYsLKz4d7VaDbVaXeL5gwdlbV6Xu5MZnnsO6NQJWLZMlmW0\nWuDPP+VwqLVr5RTmkycBFxfD3hcRkS7NmgG7dsmJVdOmAbNnAy++KC/aPvaYXChx6VJg1SrZplUr\n3efp3Bk4elTW8O3uy9JxcXGIi4urdFzlJnqNRoOMjIxSj8+bNw+DBw9+6MltKjEY/f5Er8sffwDj\nx5f9vJ0dsGYNsHkz8PLLwF9/yce8vYGgIODYMXnRlYioKtnYAE8/LVe93LYNiIiQ+00XFMhrgSNG\nAAcOAO7uZZ/D0RFo0kQOGGnT5t7jD3aCZ8+eXaGYyk30sQZ+d2jSpAmSk5OL7ycnJ8O9vHdXjkOH\n5AXU8tjYyJmto0bJ3nxhobzQQURU3WrUkMupDBsmSzh5efIiq20FC+adOsm8d3+i15dRhleKMgpR\nHTt2xPnz55GYmIj8/Hxs3LgRQ4YMqfT5r16Vf0iV6ZHXqMEkT0SmwcYGqF274kkekJNCjTX7Vu9E\nHxkZCQ8PD8THx2PQoEEYMGAAACAtLQ2DBg0CANjZ2WH58uUIDAyEv78/nnnmGfj5+VX6tc6ckbPG\nuCwBEVkLf39ZujEGs1jU7IsvgPh44OuvqzEoIiIF/f23XBIhMbHsNha1qNmZM4AeXwSIiMzW44/L\nBRxzcgw/l1kk+osXAR8fpaMgIqo+NWrI2f/l9egryiwSfWKi3NmFiMiaqFRM9EREFs1qEn12thyD\nWr++0pEQEVWvpk3lEuqGMvlEf7c3z6GVRGRtrKZHz7INEVkrq0n0//zDRE9E1kmlsqLSTdOmSkdB\nRFT9XFzkdqe3bhl2HrNI9OzRE5E1srUFPD0N79WbfKJn6YaIrJkx6vQmn+hZuiEia2aMIZYmnehz\nc+WtYUOlIyEiUoa7O5CSYtg5TDrRp6cDjRtzDD0RWa8mTYC0NMPOYdKJPi1NJnoiImvVuDGQmmrY\nOZjoiYhMmEX16AsKSj/GRE9E1s6ievQ//FD6MSZ6IrJ2DRvKPbNzc/U/h8kkel27qKSmMtETkXWz\nsQEaNTKsfKN3ot+8eTNatmyJGjVq4MiRIzrbJCcno1evXmjZsiVatWqFZcuWlXk+XYmePXoiIsPr\n9Hon+tatWyMyMhI9evQos429vT2WLFmCU6dOIT4+Hp988gnOnDmjs62utRyY6ImIdNfpi4oqfryd\nvi/s6+v70DZubm5wc3MDANStWxd+fn5IS0uDn46dvh/s0QvBRE9EBOju0VemZq93oq+sxMREHD16\nFF26dNH5/K5dYQgLk7+r1Wq0b68GADg4VE98RESm6m6PPi4uDnFxcQB0l7vLUm6i12g0yMjIKPX4\nvHnzMHjw4Aq/SE5ODoKCgrB06VLUrVtXZ5tWre4legA4e5azYomIANmjP3xYdoLVajUA4MIFYPHi\n2RU6vtxEHxsba3CABQUFGDFiBMaMGYOhQ4eW2e7BTyeWbYiIJF01+sr06I0yvFIIUebjEydOhL+/\nP0JDQ8s9x4NBp6bKTzEiImunq0ZfLYk+MjISHh4eiI+Px6BBgzBgwAAAQFpaGgYNGgQA2Lt3L9au\nXYvdu3cjICAAAQEBiImJ0Xm+B0fdMNETEUmNG8tEf3+fujK7TtmIsrrj1cjGxgYdOggcOnTvsZdf\nBnx8gKlTlYuLiMhU1K8PXLwIODnJ+z/8AIwYYVNmReV+Jjszlj16IqJ7HqzTV3uN3hh0lW54MZaI\nSHqwTl+Z0o3JJHr26ImIymYRPfqcnHsXGrRa4NIluZAPERGV7tGbZaKvUUMuxQkAmZlAgwbAI48o\nGxMRkal4sEd/40bFjzWZRN+wIXD1qvydZRsiopIe7NFfuVLxY5noiYjMwIM9+rv5siKY6ImIzMCD\nPXqzTPSPPcZET0RUFldX4PJloLBQ3jfLRN+w4b2aExM9EVFJ9vayQ5yeLu+bfY3+4kXA21vZeIiI\nTI2XF5CQIHv1Zjm88v7Szd9/M9ETET3I21vmx6wsOQS9okwm0d8t3eTkANevc/kDIqIH+fjIDUeu\nXJE5s6JMJtG7uwNJSfJNPP44YGsykRERmYa7PfqkpMpdxzSZdOrvD5w6BZw5A7RooXQ0RESmp0UL\nmSNPnQJatqz4cSaT6F1d5c/vvwe6d1c2FiIiU9S2LZCYCOzda6aJ3sZGBv7dd8CTTyodDRGR6bG3\nBzp2BCIjK5foy90cvLpNny5r9e3bKx0JEZFpeuEFWavv0KHix5hMjx4AhgwB1q2Tn1qWJC4uTukQ\nqowlvzeA78/cWeL7GzkS+PJLoFatih+jd6LfvHkzWrZsiRo1auDIkSPlttVqtQgICMDgwYP1fTmz\nZon/2O6y5PcG8P2ZO0t/fxWld6Jv3bo1IiMj0aNHj4e2Xbp0Kfz9/WFjY6PvyxERkZ70TvS+vr5o\n3rz5Q9ulpKQgOjoakyZNqtBu5UREZGTCQGq1Whw+fLjM54OCgsSRI0dEXFyceOqpp3S2AcAbb7zx\nxpset4ood9SNRqNBRkZGqcfnzZtXoXr79u3b4eLigoCAgHJrZYI9fSKiKlNuoo+NjTXo5Pv27cO2\nbdsQHR2N27dv48aNGxg3bhxWr15t0HmJiKjijDK8sqwe+bx585CcnIyEhARs2LABvXv3ZpInIqpm\neif6yMhIeHh4ID4+HoMGDcKAAQMAAGlpaRg0aJDOYzjqhoio+umd6IcNG4bk5GTk5eUhIyMDP/30\nEwCgcePGiIqKKtW+Z8+e2LZtW7nn/Pjjj+Hn54dWrVphxowZ+oZm0hYvXgxbW1tkZWUpHYpRvf76\n6/Dz80Pbtm0xfPhwXL9+XemQDBYTEwNfX180a9YMCxYsUDoco0pOTkavXr3QsmVLtGrVCsuWLVM6\npCphyXN4srOzERQUBD8/P/j7+yM+Pr7sxnoMtKkSv/zyi+jbt6/Iz88XQghx6dIlhSMyvqSkJBEY\nGChUKpW4evWq0uEY1Y4dO4RWqxVCCDFjxgwxY8YMhSMyTGFhofD29hYJCQkiPz9ftG3bVpw+fVrp\nsIwmPT1dHD16VAghxM2bN0Xz5s0t6v3dtXjxYhEcHCwGDx6sdChGN27cOLFy5UohhBAFBQUiOzu7\nzLYmswTCihUrMHPmTNj/u/6Bs7OzwhEZ36uvvoqFCxcqHUaV0Gg0sP13E4EuXbogJSVF4YgMc/Dg\nQfj4+EClUsHe3h6jR4/G1q1blQ7LaNzc3NCuXTsAQN26deHn54e0tDSFozIuS57Dc/36dfz222+Y\nMGECAMDOzg6Ojo5ltjeZRH/+/Hn8+uuv6Nq1K9RqNQ4dOqR0SEa1detWuLu7o02bNkqHUuW+/vpr\nDBw4UOkwDJKamgoPD4/i++7u7khNTVUwoqqTmJiIo0ePokuXLkqHYlTTpk3DokWLijsgliQhIQHO\nzs4YP3482rdvj8mTJyM3N7fM9tW6emVZ4/Lnzp2LwsJCXLt2DfHx8fjjjz8watQoXLx4sTrDM1h5\n72/+/PnYsWNH8WPm2MOoyLyKuXPn4pFHHkFwcHB1h2dU1jJwICcnB0FBQVi6dCnq1q2rdDhGU9E5\nPOaqsLAQR44cwfLly9GpUyeEhoYiPDwcc+bM0X1A9VSTHq5///4iLi6u+L63t7e4cuWKghEZz4kT\nJ4SLi4tQqVRCpVIJOzs70bRpU5GZmal0aEa1atUq0b17d5GXl6d0KAbbv3+/CAwMLL4/b948ER4e\nrmBExpefny/69esnlixZonQoRjdz5kzh7u4uVCqVcHNzE7Vr1xZjx45VOiyjSU9PFyqVqvj+b7/9\nJgYNGlRme5NJ9J999pl47733hBBC/PXXX8LDw0PhiKqOJV6M/emnn4S/v7+4fPmy0qEYRUFBgfDy\n8hIJCQnizp07FncxtqioSIwdO1aEhoYqHUqVK2/5FXP25JNPir/++ksIIcSsWbPEG2+8UWZbk9l4\nZMKECZgwYQJat26NRx55xKInVlliWeDll19Gfn4+NBoNAKBbt2749NNPFY5Kf3Z2dli+fDkCAwOh\n1WoxceJE+Pn5KR2W0ezduxdr165FmzZtEBAQAACYP38++vfvr3BkVcMS/899/PHHePbZZ5Gfnw9v\nb2+sWrWqzLY2QphhsZiIiCrM8i5HExFRCUz0REQWjomeiMjCMdETEVk4JnoiIgvHRE9EZOH+Pzf3\nzFpxvrvtAAAAAElFTkSuQmCC\n"
}
],
"prompt_number": 6
},
{
"cell_type": "markdown",
"metadata": {},
"source": "####Bonus"
},
{
"cell_type": "markdown",
"metadata": {},
"source": "What is the expected number of independent Uniform[0, 1] deviates you need to sum to exceed 1? (I learned about this from http://nbviewer.ipython.org/4731599/.)"
},
{
"cell_type": "code",
"collapsed": false,
"input": "def how_many_to_exceed(threshold):\n total = 0\n count = 0\n while total < threshold:\n total += np.random.uniform()\n count += 1\n return count\n\nnum_draws = 100000\ndraws = [how_many_to_exceed(1) for i in range(num_draws)]\nrunning_average = np.cumsum(draws) / (np.arange(num_draws) + 1.)\n\nplt.plot(running_average, label='Average draws')\nplt.axhline(np.e, color='red', linestyle='--', label='e = {:.8f}'.format(np.e))\nplt.ylim(1, 4)\nplt.legend()\nplt.gca().set_xscale('log')",
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEFCAYAAAAG45eHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlcVPX+P/DXIKQhyKKBCKQmWOACSEneXPCSmuu1cq0U\nlxL5XrO81c9bXUO7ptJyTeteU1MDvaTlbaHErWzIJaTSLLdcSmQRFREXhFjm/fvjkzOODDDIDEPH\n1/PxmIee+XzmnM98OLzOZz6cOUcnIgIiItIkJ0c3gIiI7IchT0SkYQx5IiINY8gTEWkYQ56ISMMY\n8kREGmZVyFdWViIiIgJDhw61WD59+nQEBwcjLCwMe/futWkDiYjoxlkV8osWLUJoaCh0Ol2VsrS0\nNBw7dgxHjx7FsmXLEB8fb/NGEhHRjak15HNycpCWlobHH38clr43lZqaitjYWABAVFQUioqKcPr0\nadu3lIiI6sy5tgozZszAa6+9hosXL1osz83NRWBgoHE5ICAAOTk58PX1NT5n6RMAERHVrr4XJahx\nJP/555/Dx8cHERERNW7o+jJLoS4iNnskJCTYtH5N5ZbKanvu+vKaytgX7IuboS/qssy+MC3bQo0h\nv2vXLqSmpqJ9+/YYO3Ystm3bhvHjx5vV8ff3R3Z2tnE5JycH/v7+NmlcdaKjo21av6ZyS2W1PXd9\neV3bWxfsixtfN/vC+vr17Yu6LtvSTd8XYiW9Xi9Dhgyp8vyGDRtk4MCBIiLyzTffSFRUVJU6ddiM\n5iUkJDi6CY0G+8KEfWHCvjCxRXbWOid/ravTMEuXLgUAxMXFYdCgQUhLS0NQUBCaN2+OVatW2fYo\npDH2HLH80bAvTNgXJuwL29L9frSw70Z0OpvNLxER3SxskZ11GskTkXW8vb1x/vx5RzeD/iC8vLxQ\nWFhol3VzJE9kB9znqS6q219ssR/x2jVERBrGkCci0jCGPBGRhjHkiUiT9Hq92SVXblYMeaKbWHR0\nNLy9vVFWVuboppCdMOSJblInTpxAZmYmfHx8kJqaavP1V1RU2HydttKY22ZrDHmim1RycjLuv/9+\njBs3DklJSQCA3377DZ6enjhw4ICx3tmzZ+Hq6oqCggIA6sKF4eHh8PLywn333YeffvrJWLddu3Z4\n9dVX0bVrV7i7u6OyshILFixAUFAQWrRogU6dOuGTTz4x1jcYDHjmmWdw22234Y477sDbb78NJycn\nGAwGAMCFCxcwefJktGnTBgEBAZg1a5ax7HolJSWYMGECvL290alTJ3z77bdm5XVtW9u2bbFnzx4A\nwH//+184OTnh0KFDAIAVK1bgwQcfBABkZmbi7rvvhoeHB1q3bo1nnnnmxn4g9lLvCyNYoYE2Q9Ro\n/BH2+Q4dOsiaNWvkyJEj4uLiImfOnBERkUmTJsmLL75orPf2228br0+1Z88e8fHxkczMTDEYDJKU\nlCTt2rWTsrIyERFp27atRERESE5OjpSWloqIyIcffiinTp0SEZF169ZJ8+bNJT8/X0RElixZIqGh\noZKbmyvnz5+XmJgYcXJyksrKShERGT58uEydOlWuXLkiZ86cke7du8vSpUstvp+ZM2dK79695fz5\n85KdnS2dOnWSwMBAY3ld2zZ+/Hh54403RETkiSeekKCgIFmyZImIiIwbN07efPNNERG59957Zc2a\nNSIiUlxcLBkZGXX+WVS3v9hiP2LIE9mBNfs8UP/Hjdq+fbs0a9ZMLl68KCIiYWFhsnDhQhER+eKL\nL6RDhw7Gun/6059k9erVIiIydepUmTVrltm67rzzTvn6669FRKRdu3ayatWqGrcdHh4uqampIiLS\nt29fWbZsmbHsiy++EJ1OJ5WVlZKfny9NmzaVkpISY3lKSor07dvX4nrvuOMO2bx5s3F52bJlEhAQ\nYFy2tm2ffvqpiIisWLFChg0bJiIiISEhsmLFChkzZoyIqAPG3r17RUSkd+/ekpCQIGfPnq1x3TWx\nZ8hzuobIQWwR8zcqKSkJ/fv3h7u7OwBg5MiRximb6OhoXLlyBZmZmThx4gT27dtnnJrIysrCG2+8\nAS8vL+MjJycHeXl5xnVff0ZLcnIyIiIijPX3799vnPo5depUlZsOXZWVlYXy8nL4+fkZXzt16lSc\nPXvW4nvKy8szW9ftt99epY41bTt37hwAoHfv3ti+fTvy8/NRWVmJkSNHYufOncjKysKFCxcQHh4O\nQE3dHDlyBCEhIejevTs2bNhQU9c3OF67hugmU1JSgg8++AAGgwF+fn4A1Fx8UVERfvzxR3Tt2hWj\nRo3C+++/Dx8fHwwdOhTNmzcHoILzxRdfxAsvvFDt+q+9aVBWVhamTJmCbdu2oUePHtDpdGY3IfLz\n8zO7H8W1/w8MDETTpk1x7tw5ODnVPh718/PDyZMnERISAgA4efJkvdoWFBQEV1dXvPXWW+jTpw/c\n3d3RunVrLFu2DL169TKuJygoCCkpKQCA//3vfxgxYgQKCwtx66231trmhsCRPNFN5pNPPoGzszMO\nHTqEffv2Yd++fTh06BB69eqF5ORkAMAjjzyCtWvXIiUlBY888ojxtU888QTeeecdZGZmQkRQXFyM\nDRs24PLlyxa3VVxcDJ1Oh1atWsFgMGDVqlXYv3+/sXzUqFFYtGgR8vLyUFRUhMTERGMQ+/n5oX//\n/vjb3/6GS5cuwWAw4Pjx4/j6668tbmvUqFGYP38+ioqKkJOTg7feeqvGfqitbQDQp08fvP322+jT\npw8A9Snn2mUAWLNmjfHThYeHB3Q6nVUHpYbSeFpCRA0iOTkZkyZNQkBAAHx8fODj4wNfX19MmzYN\nKSkpMBgM6N69O9zc3HDq1CkMHDjQ+NrIyEgsX74c06ZNg7e3N4KDg5GcnFztfZxDQ0PxzDPPoEeP\nHmjdujX279+Pnj17GsufeOIJ9O/fH127dkVkZCQGDx6MJk2aGEMyOTkZZWVlCA0Nhbe3N0aOHIn8\n/HyL20pISEDbtm3Rvn17PPDAAxg/fnyN95eurW2ACvnLly+jd+/eFpcBYPPmzejcuTPc3d0xY8YM\nrF27Fk2bNq3lp9BweBVKIjvgPn9jNm7ciPj4eJw4ccLRTWlQvAolEWlSaWkp0tLSUFFRgdzcXMyZ\nMwcPPfSQo5ulKRzJE9kB93nrlJSUoE+fPjh8+DBuvfVWDBkyBIsWLYKbm5ujm9ag7DmSZ8gT2QH3\neaoLTtcQEdENYcgTEWkYQ56ISMMY8kREGsaQJyLSMIY8EZGGMeSJyC6effZZdOzYES1atEBISAhW\nr15dbd158+bB3d3d+HB1dUWTJk1QWFgIAPjggw/wpz/9Cc2bN0ffvn2rvH7Hjh2455574OHhgQ4d\nOmD58uXGsqSkJONNPQIDAzFz5kxUVlYay3NycjB06FC0bNkSfn5+ePLJJ43l5eXlGDFiBNq3bw8n\nJyekp6ebbbeiogJPPvkk/Pz80LJlSwwbNszsipx9+/aFj4+PsQ+ubVeDqffFiq3QQJshajS4z4sk\nJCTIzz//LCIiu3fvFi8vL9m1a5dVr509e7bExMQYl7/44gv58MMP5eWXX5bo6GizuhUVFdKqVSvj\ndem//fZbcXNzkx9//FFE1I1JduzYIeXl5ZKbmyuRkZGyYMEC4+sffPBBmTBhgvz222+Sn58vXbp0\nkcWLF4uISFlZmSxatEh27Nghfn5+kp6ebrbtRYsWSVhYmJw5c0ZKS0tl/Pjx8tBDDxnLf/zxR+MN\nVXbv3i1NmzaVw4cPV3m/1e0vttiPOJInugnl5eXh4Ycfho+PD+64445ar9h4I2bPno2OHTsCALp3\n745evXrhm2++qfV1IoKkpCTExsYan4uJicGIESOMl0a+1unTp3Hu3DmMGzcOAHD33XcjJCQEBw8e\nBABMnToV9913H5ydndGmTRs8+uij2Llzp/H1Bw4cwOjRo3HLLbfA19cXDzzwgPH2hy4uLpg+fTru\nu+8+NGnSpMq2Dxw4gAEDBuC2225D06ZNMWrUKLNbJ3bp0gUuLi7GZTc3N7Ro0aLWPrClGkO+tLQU\nUVFRCA8PR2hoKJ5//vkqdfR6PTw8PBAREYGIiAjMnTvXbo0lovozGAwYOnQoIiIikJeXhy+//BJv\nvvkmtmzZYrH+ggULzG4Scu3D29vbqm2WlJTg22+/RefOnWutu337dpw9exYPP/ywVetu06YNunbt\nipUrV6KyshK7du1CVlZWlStKXpWenm7WjgEDBiAlJQUlJSXIzc3Fxo0bza68WZP+/ftj48aNOHXq\nFK5cuYL//ve/GDRokFmdIUOG4NZbb0V0dDRWrlxp8UBlV7UN9YuLi0VEpLy8XKKiomT79u1m5V99\n9ZUMHTq0xnVYsRkiTbFqn09IsHzDp4QE6+pXV68WGRkZcvvtt5s9N2/ePJk4ceINrc8a48ePN94n\ntjaTJk2qti3Lly+vMl0joqZCWrVqJc7OzuLs7CzvvvuuxdevWLFCAgMD5dy5c8bnzp07JxEREeLs\n7Cw6na7abQcEBFSZrhFR702n04mzs7N069ZNCgsLq9SpqKiQDz/8ULy8vCQrK6tKeXX7iy2ys9bp\nGldXVwBAWVkZKisrLR65hdfoIKq72bMt39Nv9mzr6ldXrxZZWVnIy8szG5HPnz8fZ86cucE3UrPn\nnnsOBw8exAcffFBr3StXrmD9+vVmUzW1yc3NxZAhQ5CSkoLy8nIcOHAAiYmJSEtLM6v3ySef4IUX\nXsDGjRuNOSYiGDBgAEaOHIkrV66goKAAhYWFmDlzplXbfvbZZ3Hp0iUUFhaiuLgYDz74oMVPAU2a\nNMGIESMQFRWFjz/+2Or3Zgu13v7PYDCgW7duOH78OOLj4xEaGmpWrtPpsGvXLoSFhcHf3x+vv/56\nlTqAmp+7Kjo6GtHR0fVuPBHV3e2334727dvjyJEjVtWfN28e5s+fb7FMp9Ph4sWL1b42ISEBmzdv\nRnp6ulVXlvz444/RsmVLszsvXb+96+3atQsBAQHo168fAKBjx44YPHgwNm7caJw62bRpE6ZMmYK0\ntDR06tTJ+NqCggJ8//332LZtG1xcXODt7Y0JEyZg1qxZSExMrLW9mzZtwvz58+Hp6QkAmDZtGl56\n6SUUFhZaHBCXl5cbb6VoiV6vh16vr3W7dWLtkL+oqEiioqLkq6++Mnv+4sWLximdtLQ0CQ4OrvLa\nOmyGSBMa8z5fWVkp3bp1k8TERLly5YpUVFTITz/9JN9++61NtzNv3jwJDg6W/Px8q1/Tr18/SbAw\nDVVZWSklJSWyZMkS6d27t5SWlhrPWjl48KC4urrKtm3bxGAwyLFjxyQoKEiWL18uIiJffvmleHt7\nV5lqFhExGAzSpk0bSUxMlIqKCjl//rwMHz5cHn30UWOd0tJSKSkpkYCAANmyZYuUlJQYy8aOHSsP\nP/ywXLhwQcrKyuSVV16RgIAAERE5fPiwpKWlyZUrV6SsrExWr14tHh4eDT5dU6c1vPzyy/Laa6/V\nWKddu3Zm810ijXuHJ7KHxr7P5+XlydixY6V169bi5eUlPXr0kC+//NKm29DpdNKsWTNxc3MzPubP\nn28sd3Nzkx07dhiXc3JyxMXFRY4fP15lXatWrRKdTmf2uHbuPCkpSUJCQsTd3V0CAgLk73//u7Gs\nb9++4uLiYtaOQYMGGcszMjKkZ8+e4unpKa1atZLRo0fLmTNnjOVt27YVnU4nTk5Oxn+vBnV+fr6M\nHDlSWrVqJZ6entKrVy/jwfLQoUMSFRUl7u7u4u3tLX369DF7v9eyZ8jXeD35goICODs7w9PTEyUl\nJRgwYAASEhIQExNjrHP69Gn4+PhAp9MhMzMTo0aNqnLrLl5bm2423OepLux5Pfka5+RPnTqF2NhY\nGAwGGAwGjBs3DjExMVi6dCkAIC4uDuvXr8eSJUvg7OwMV1dXrF27tl4NIiIi2+GdoYjsgPs81QXv\nDEVERDek1lMoiajuvLy8LJ7uR2SJl5eX3dbN6RoiokaK0zVERFQjhjwRkYYx5ImINIwhT0SkYQx5\nIiINY8gTEWkYQ56ISMMY8kREGsaQJyLSMIY8EZGGMeSJiDSMIU9EpGEMeSIiDWPIExFpGEOeiEjD\nGPJERBrGkCci0jCGPBGRhjHkiYg0jCFPRKRhDHkiIg1jyBMRaRhDnohIwxjyREQaxpAnItKwGkO+\ntLQUUVFRCA8PR2hoKJ5//nmL9aZPn47g4GCEhYVh7969dmkoERHVnXNNhc2aNcNXX30FV1dXVFRU\noGfPntixYwd69uxprJOWloZjx47h6NGj2L17N+Lj45GRkWH3hhMRUe1qna5xdXUFAJSVlaGyshLe\n3t5m5ampqYiNjQUAREVFoaioCKdPn7ZDU4mIqK5qHMkDgMFgQLdu3XD8+HHEx8cjNDTUrDw3NxeB\ngYHG5YCAAOTk5MDX19es3uzZs43/j46ORnR0dP1aTkSkMXq9Hnq93qbrrDXknZyc8MMPP+DChQsY\nMGAA9Hp9lYAWEbNlnU5XZT3XhjwREVV1/QB4zpw59V6n1WfXeHh4YPDgwfjuu+/Mnvf390d2drZx\nOScnB/7+/vVuGBER1V+NIV9QUICioiIAQElJCbZu3YqIiAizOsOGDUNycjIAICMjA56enlWmaoiI\nyDFqnK45deoUYmNjYTAYYDAYMG7cOMTExGDp0qUAgLi4OAwaNAhpaWkICgpC8+bNsWrVqgZpOBER\n1U4n10+o22MjOl2VeXuqqrwcyM0F2rVzdEuIqDGwRXbW+odXso/Ll4EffwT27lWPH34ADh4EDAZg\n82agTx9Ht5CItICXNWhg69YBd94J+PgATz0F7NsHREYCb70FnD0L/PvfwCuvOLqVRKQVnK6xg9JS\n4ORJoGNH8+dFgJAQYMECYMgQwNnC56iyMiAoCFi/HujevWHaS0SNky2ykyN5Gzl8GHjtNaB/f+C2\n24DwcOCnn8zrZGSooP/LXywHPADccgvw//4fR/NEZBsMeRu4eBHo0QPIygL++lcgJweYOxd46SXz\neu+9B0yYAFj4rpiZyZOBzMyqBwkiorridI0NrFwJfPYZ8PHHpudKSoDgYPXcPfcAV64AAQEquK35\nrtirr6o/yL7/vv3aTUSNG6drGomrI/Rr3Xor8OKLwKxZavnjj4GoKOsCHgDi44EvvgCOHLFlS4no\nZsOQr6djx4CffwYGDapaNnmyCunt2y0fCGri7g5Mm6b+SHuzOHIE6NoV+P0L1ERkA5yuqaeXXgIu\nXQIWLrRcnpQEvPEGkJen5uqbNbN+3YWFasrn++//+F+Q+uUXNfV08aLqk7VrgesvRHr1bxWenuoM\no3/9S01xnTunvlMwYADg5wc4cWhCNwlbZCdDvh4MBqB9eyA1FQgLs1ynogLo3BmIiVHnwNfVrFnq\ndMykpPq11dH+8x/1R+lr/fqr6eD11lvA9OlAero6A2n6dBXsgAr2U6fU/1esACZNarBmEzkU5+Qd\nTK8HvL2rD3hAnSr52WfAyy/f2Daeew7YtEl9acqRKiuB//1Pfbq40df/3/8BRUXq8g1//jMwfLg6\nUO7bB2zZAnz0EdC7t/q273ffqb9J/PyzOhhs2ABMnKimwPz81Kj/55/VH7iJqHocydfD+PHq26pP\nPWXf7SxerII+Lc2+26mOiArfw4eBMWPUqDsrC2jbtmrdX35Ro/Prp1Rmz1bruXp57MpKVc9gUFNZ\ngLpuT5s2NbflzTeBGTNMy3feqaazmje/wTdH1IhxJO9AFy+qaZpHHrH/tqZOVaPWr76qWrZjhzrd\n0p727FHTJfv3Aykp6tu67dqpPyhfa/duoEMHNT2Vk6Oe++ADYPRo4N13gZYtTXWbNFGXeMjLU8H+\n9tu1BzwAPP20+kZxVpapX9zc1B+2b3bFxepASmRGGkADbcaizEyRgQNFXnnFtutdsUJk+HDbrrMm\n778vcvfdIpWV5s+PHi3SpInIgQOWX5eXV/U11sjNFRkwQOT4cZE//1nk+efV80uWqL5cvlykRw8R\ng8H0mhkzRMaOFWnRQgQQOXxYJDBQ/R8Q+fLLqtvZt0/k8uW6t++q48dF3N1N29i4UaSw8MbXZwuV\nlSLl5TXXubbf6mPvXtN7v/rw9a36XOvWIqmpIhcv2ma71DBskZ2aDfkDB0QeekikTRuRf/1L/btl\ni+3W36uXyCef2G59tamsFImMFFm71vRcSYmIh4fIzJkigwZVfc3rr4u4uKh/6+rtt1U43HqrSEKC\nSGmpeXlFhUiXLqoPSktVaLVvr0K7vFwkPl5t+8knRcrKqr7e1qZONQ+1xYur1jl/3r5tKC01b8Nz\nz6nH1eUFC0Qef9y8ziOPiOh0IrfdppbvvVfk/vtFVq0SOXLEfP1nz4oMGSIyfbo6mMXEmK9n1SqR\njz5S+0nTpuoA+9hjIosWVQ19b2+RRx9VP6dly0QWLlTP33OPqU7HjurgPmOGyLhxIgcP2rf/DAZ1\n0Nq6Vf38Pv1U5K9/VYOD5GSRNWtEXnhBZMMGkXXr1EBr507bHTAboz9WyFt6JCRYfkFCgk3q7+iX\nIMXFqsq2bSJ+fmqEaqv1V8yyb/uvr//llyJ33CFS/g/L9Y+NM9UvKFAj3C1bRBJd7dOezz9Xo8bZ\nsFz/1NSG6x+DQR0I//UvkfVdLNdPQIJx8fPPRUaM+P3gY6P2XF1/aqrIa69ds5pq+ufa9pht1sr6\nTZuKJCbWrf0//STyUdfq13/nneqTYWCgSEBAze1Z6JEgbduKeHmpA3qLFuqTnK3fb6tWav1t24qE\nhtbQHs8EcXJSiy4uav+PihJZ2baa/TMuQY4fFykqMh0oDAaRy89arr/lTwkybpxI//5qgJOYKDJr\nlsj3Qy3Xn98sQSZOVAemF14QmTdPZNIkkdVBluvr+ybIhAnqoDp5snoPf6yQt7P8fDVq9PYW+cc/\n1A/uenPmiERHq1HojdqyReT223//5XKAYcPUziIiEhurRmkiagTXpYvpvb33nsiDD6r/Jyaq996y\n5e8HuRps2aJGiu7utU+jGAxqBHh1pNiQ01e1+dvfTL8/Dzyg/nV1tZyFx4+r9/rVV6oPBwxQn2Te\neENk9WqRM2dM6zUY1C+sh0fV9Tz+uPl0SEmJyNGjpuXKyqqjzppGoatXq4HJ1fU/+qh6PitLJDu7\n3l1UZ8eOicyerUb3998v0rmziI+PeR9ERqrQj4wU6d5dTZW6uKjlYcPUFN4HH4hMmKBG5jNmiOze\nLZKTo36H6zoqNxhEvv5aHVwSE1XbBgxQn9x79FC/q87OIoMHi4wapT799O5teT+4/nH77aafs5OT\n+mQDqE/NkZEi/fqZpsb69ze9/zvvNNW9+nBxUX3Wtas6aEVFmcpatVJTovffLxIcrA6uMTFXf/YM\neTl/Xv3SeXuLPPWUyOnT1detqFCd99JLdd/OhQsiU6aoH/zmzTfe3vo6ftx0IPPxMYW2waB23nfe\nUcvDh4skJZm/9vnnRcaPV9MA3buLfPihyJUrpvKiIjUaa9tW/ULWRVGRGD81NQYXL6qf0wsvqF+k\nWbNMB8CKCpHvvlN/r2jWzLpfeB8fNboaMsT8+SlTVPjZeyqI7KuyUmT/fpFffmlc0z83dcgXF4vM\nn6+OgpMmqdGNNfLz1VF+61brt7V1qwq+xx+3/Amhob38svqbwPXv+YcfVBjt2qXC+tw58/KLF9V7\nb9lS5Omn1cGiWTNVX0TNcw4cKLJnj5ob1YqaPrmVlanRaffuIps2mZ5LShJJT1ej+Gv/eAyI/Pvf\njSsISLtuypA3GERWrlQfZUaOFDl0qO7ruDo/n5dXc72SEvWHn8BAUwA0dk8+KXLLLSLvvmu5fOdO\nNYoVUXXmzRMJD1fvtX9/kf/8p+Ha+kd07ScfInuzRXb+ob4MdeCAujpjaam6RMA999z4ul5+WZ13\n/sUX6pzt6/3yCzByJHDHHeocbw+PG99WQyopAbKzq96VqjoiQN++wIkTQJcu6tx1V1e7NpGIrHTT\nfBmquBj4+9/VBa3GjAG++aZ+AQ+oywA7OVm+3EBqKnDvvUBsrPoyzx8l4AF1iWNrAx5QlwdYtkwd\nPD/5hAFPpDWNfiSfmqouVtWzJ/D660Dr1rZrV36+uixBUhJw//3qYmIvvqiulrhunbrbExGRo9hi\nJF/NnUYd7+RJFe6HDqkrD8bE2H4brVsDq1cDjz2mLiI2Y4Yaye7ZA7RqZfvtERE1tEY3XVNerm6I\n3a2bGmX/+KN9Av6qP/8ZiIsDuncH+vVTFwFjwBORVjSq6ZqdO9VFp/z91QWrgoLs3TLFYFB/aG2o\n7RERWUMzNw0pKABmzgQ2b1Z3WBoxwnSXICKim9Uf/uwagwFYuRLo1End0/TgQXXaIgOeiMg2agz5\n7Oxs9O3bF506dULnzp2xePHiKnX0ej08PDwQERGBiIgIzJ0716oN79+v7gC0dKm6IcabbwItWtzY\nmyAiIstqPLvGxcUFCxcuRHh4OC5fvozIyEj069cPISEhZvX69OmD1NRUqzZYXKzOTV+5EvjnP4En\nnrD8ZSQiIqq/GkfyrVu3Rnh4OADAzc0NISEhyLt6r7ZrWDtn9OmnQGiouhvQ/v3qj6wMeCIi+7H6\nPPkTJ05g7969iIqKMntep9Nh165dCAsLg7+/P15//XWEhoZWef1dd81GQQEweDAwcWI0fH2j6914\nIiIt0ev10Ov1Nl2nVWfXXL58GdHR0fjHP/6B4cOHm5VdunQJTZo0gaurKzZu3IinnnoKR44cMd+I\nTod//lPw3HNA06Y2bT8RkWY1yCmU5eXlGDJkCAYOHIinn3661hW2b98e33//Pby9vW3aUCKim43d\nT6EUEUyePBmhoaHVBvzp06eNjcjMzISImAU8ERE5To1z8jt37sSaNWvQtWtXREREAADmzZuHkydP\nAgDi4uKwfv16LFmyBM7OznB1dcXatWvt32oiIrJKo/jGKxERVfWH/8YrERHZF0OeiEjDGPJERBrG\nkCci0jCGPBGRhjHkiYg0jCFPRKRhDHkiIg1jyBMRaRhDnohIwxjyREQaxpAnItIwhjwRkYYx5ImI\nNIwhT0RvxBllAAAJGklEQVSkYQx5IiINY8gTEWkYQ56ISMMY8kREGsaQJyLSMIY8EZGGMeSJiDSM\nIU9EpGEMeSIiDWPIExFpGEOeiEjDGPJERBrGkCci0rAaQz47Oxt9+/ZFp06d0LlzZyxevNhivenT\npyM4OBhhYWHYu3evXRpKRER151xToYuLCxYuXIjw8HBcvnwZkZGR6NevH0JCQox10tLScOzYMRw9\nehS7d+9GfHw8MjIy7N5wIiKqXY0j+datWyM8PBwA4ObmhpCQEOTl5ZnVSU1NRWxsLAAgKioKRUVF\nOH36tJ2aS0REdVHjSP5aJ06cwN69exEVFWX2fG5uLgIDA43LAQEByMnJga+vr1m92bNnG/8fHR2N\n6OjoG2sxEZFG6fV66PV6m67TqpC/fPkyRowYgUWLFsHNza1KuYiYLet0uip1rg15IiKq6voB8Jw5\nc+q9zlrPrikvL8fDDz+Mxx57DMOHD69S7u/vj+zsbONyTk4O/P39690wIiKqvxpDXkQwefJkhIaG\n4umnn7ZYZ9iwYUhOTgYAZGRkwNPTs8pUDREROYZOrp9rucaOHTvQu3dvdO3a1TgFM2/ePJw8eRIA\nEBcXBwCYNm0aNm3ahObNm2PVqlXo1q2b+UZ0uipTOkREVDNbZGeNIW8rDHkiorqzRXbyG69ERBrG\nkCci0jCGPBGRhjHkiYg0jCFPRKRhDHkiIg1jyBMRaRhDnohIwxjyREQaxpAnItIwhjwRkYYx5ImI\nNIwhT0SkYQx5IiINY8gTEWkYQ56ISMMY8kREGsaQJyLSMIY8EZGGMeSJiDSMIU9EpGEMeSIiDWPI\nExFpGEOeiEjDGPJERBrGkCci0jCGPBGRhjHkiYg0jCFPRKRhtYb8pEmT4Ovriy5dulgs1+v18PDw\nQEREBCIiIjB37lybN5KIiG6Mc20VJk6ciCeffBLjx4+vtk6fPn2Qmppq04YREVH91TqS79WrF7y8\nvGqsIyI2axAREdlOrSP52uh0OuzatQthYWHw9/fH66+/jtDQ0Cr1Zs+ebfx/dHQ0oqOj67tpIiJN\n0ev10Ov1Nl2nTqwYhp84cQJDhw7FTz/9VKXs0qVLaNKkCVxdXbFx40Y89dRTOHLkiPlGdDqO9omI\n6sgW2Vnvs2vc3d3h6uoKABg4cCDKy8tRWFhY39USEZEN1DvkT58+bTzSZGZmQkTg7e1d74YREVH9\n1TonP3bsWKSnp6OgoACBgYGYM2cOysvLAQBxcXFYv349lixZAmdnZ7i6umLt2rV2bzQREVnHqjn5\nem+Ec/JERHXWKObkiYio8WLIExFpGEOeiEjDGPJERBrGkCci0jCGPBGRhjHkiYg0jCFPRKRhDHki\nIg1jyBMRaRhDnohIwxjyREQaxpAnItIwhjwRkYYx5ImINIwhT0SkYQx5IiINY8gTEWkYQ56ISMMY\n8kREGsaQJyLSMIY8EZGGMeSJiDSMIU9EpGEMeSIiDWPIExFpGEOeiEjDGPINTK/XO7oJjQb7woR9\nYcK+sK0aQ37SpEnw9fVFly5dqq0zffp0BAcHIywsDHv37rV5A7WGO7AJ+8KEfWHCvrCtGkN+4sSJ\n2LRpU7XlaWlpOHbsGI4ePYply5YhPj7e5g20pK47QW31ayq3VFbbc9eX23OnZV/c+LrZF9bXr29f\n1HXZlm72vqgx5Hv16gUvL69qy1NTUxEbGwsAiIqKQlFREU6fPm3TBlrS2H5o1z93M+/A1z/HvtBX\nW34z9QVD/saX601q8euvv0rnzp0tlg0ZMkR27txpXI6JiZHvvvuuSj0AfPDBBx983MCjvpxRTyrD\nTXQ6Xa11iIioYdTr7Bp/f39kZ2cbl3NycuDv71/vRhERkW3UK+SHDRuG5ORkAEBGRgY8PT3h6+tr\nk4YREVH91ThdM3bsWKSnp6OgoACBgYGYM2cOysvLAQBxcXEYNGgQ0tLSEBQUhObNm2PVqlUN0mgi\nIrKOTjhhTkSkWfzGKxGRhjkk5IuLixEbG4spU6YgJSXFEU1oNH799Vc8/vjjGDlypKOb4nCffvop\npkyZgjFjxmDr1q2Obo5DHT58GPHx8Rg1ahRWrFjh6OY4XHFxMe655x5s2LDB0U1xKL1ej169eiE+\nPh7p6elWvcYhIf/RRx9h1KhRWLZsGVJTUx3RhEajffv2ePfddx3djEbhL3/5C5YtW4Z33nkH69at\nc3RzHOquu+7CkiVLsHbtWmzevNnRzXG4V199FaNHj3Z0MxzOyckJ7u7u+O233xAQEGDda+zcJoty\nc3MRGBgIAGjSpIkjmkCN2Ny5czFt2jRHN8PhPvvsMwwePBhjxoxxdFMcauvWrQgNDcVtt93m6KY4\nXK9evZCWloYFCxYgISHBqtfYLOSru5jZpk2bcNdddyE4OBiJiYkAgICAAOP59QaDwVZNaDTq0hda\nV5e+EBHMnDkTAwcORHh4uCOaa1d13S+GDh2KjRs3IikpqaGband16Yv09HRkZGQgJSUFy5cv19yX\nK+vSF1e/bOrp6YnffvvNug3U+zuzv/v6669lz549ZpdAqKiokA4dOsivv/4qZWVlEhYWJgcPHpTi\n4mKZOHGixMfHS0pKiq2a0GjUpS/OnTsncXFxEhQUJAsWLHBgq+2jLn2xePFiiYyMlKlTp8o777zj\nwFbbR136Qq/Xy/Tp02XKlCmycOFCB7baPurSF1e99957smHDBkc0167q0hcfffSRxMXFyejRoyU9\nPd2q9dss5EWqXudm165dMmDAAOPy/PnzZf78+bbcZKPFvjBhX5iwL0zYFyb27Au7zslfO/cOqGma\n3Nxce26y0WJfmLAvTNgXJuwLE1v2hV1D3tLFym5W7AsT9oUJ+8KEfWFiy76wa8hffwGz7Oxsq0/7\n0Rr2hQn7woR9YcK+MLFlX9g15O+++24cPXoUJ06cQFlZGdatW4dhw4bZc5ONFvvChH1hwr4wYV+Y\n2LQv6v0Xg9+NGTNG/Pz85JZbbpGAgABZuXKliIikpaVJx44dpUOHDjJv3jxbba5RY1+YsC9M2Bcm\n7AsTe/cFL1BGRKRhvEAZEZGGMeSJiDSMIU9EpGEMeSIiDWPIExFpGEOeiEjDGPJERBrGkCci0jCG\nPBGRhv1/8YeFUQcDoO0AAAAASUVORK5CYII=\n"
}
],
"prompt_number": 7
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment