Skip to content

Instantly share code, notes, and snippets.

@tmiyama
Created April 4, 2013 16:10
Show Gist options
  • Save tmiyama/5311742 to your computer and use it in GitHub Desktop.
Save tmiyama/5311742 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"metadata": {
"name": "chi2"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"<h1> Learn $\\chi^2$ (chi-square) distribution </h1>"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pylab inline"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"\n",
"Welcome to pylab, a matplotlib-based Python environment [backend: module://IPython.zmq.pylab.backend_inline].\n",
"For more information, type 'help(pylab)'.\n"
]
}
],
"prompt_number": 1
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from scipy.stats import chi2,norm\n",
"import numpy as np"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 43
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"<h2> If $v$ is normally distributed random variable, then $v^2$ is $\\chi^2$ distributed with one degree of freedom </h2>"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"df=1\n",
"rv=chi2(df)\n",
"x=np.linspace(0,10)\n",
"chi2pdf=rv.pdf(x)\u0000\n",
"plot(x,chi2pdf),ylim([0,0.5])"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "pyout",
"prompt_number": 139,
"text": [
"([<matplotlib.lines.Line2D at 0xea14f4c>], (0, 0.5))"
]
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAD9CAYAAABDaefJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAG9dJREFUeJzt3XtwVGWexvFvQ5pbEu4Qhu5ghMSkA3LRxAsWawZRlJUM\nCjUTdmdwmMCmVByZomota6vWYM2wZHFG2cGysq7X1Q2MlyWZEaIbMLUokqCoMICCSKQTLoPIRa5J\nenr/eDc3SDoX0jnd5zyfqlPdnZyc/tGlT978znve4woGg0FERMQWelldgIiIdB+FuoiIjSjURURs\nRKEuImIjCnURERtRqIuI2Ei7oV5aWkpaWhopKSkUFBRc8f3y8nIGDRrElClTmDJlCr/+9a/DUqiI\niLQvJtQ3A4EAS5YsoaysDI/HQ2ZmJtnZ2fh8vhb73X777ZSUlIS1UBERaV/IkXplZSXJyckkJSXh\ndrvJycmhuLj4iv10/ZKISGQIOVKvqakhMTGx8bXX66WioqLFPi6Xi61btzJp0iQ8Hg9PPfUU6enp\nV+wjIiKd19lBc8iRekfC+IYbbsDv9/P555/zyCOPMGfOnDYL6+y2bl2Q++/v/M9F8vbEE09YXkOk\nbPos9Fnoswi9dUXIUPd4PPj9/sbXfr8fr9fbYp/4+HgGDBgAwD333ENdXR3fffddl4q5XGIiNHt7\nERFpR8hQz8jIYP/+/VRVVVFbW8u6devIzs5usc+xY8caf6NUVlYSDAYZOnRotxQ3ZoxCXUSkM0L2\n1GNiYlizZg0zZ84kEAiQm5uLz+ejsLAQgLy8PN58802ee+45YmJiGDBgAGvXru224kaNghMn4NIl\n6Nu32w5rqaysLKtLiBj6LJros2iiz+LquIJdbdx05k1cri73h665Bt5/H8aO7eaiREQiXFeyM+Kv\nKFULRkSk4yI+1HWyVESk46Ii1A8dsroKEZHoEBWhrpG6iEjHRHyoq6cuItJxER/qar+IiHRcVIS6\nRuoiIh0T8aE+bJi5+OjsWasrERGJfBEf6i6XRusiIh0V8aEO6quLiHRU1IS6RuoiIu2LilDXtEYR\nkY6JilBX+0VEpGOiJtQ1UhcRaV9UhLraLyIiHRPx66mDmaM+ciScO2emOIqIOIEt11MHiIszdz46\nccLqSkREIltUhDqoBSMi0hFRE+o6WSoi0r6oCnVNaxQRCS2qQl0jdRGR0KIm1NVTFxFpX9SEukbq\nIiLti6pQV09dRCS0qLj4CMyNMgYOhPPnoXfvbipMRCSC2fbiIzAXHw0dCkePWl2JiEjkippQB7Vg\nRETaE3WhrpOlIiJti6pQ17RGEZHQoirUNVIXEQkt6kJdPXURkbZFVair/SIiElpUhbraLyIioUXN\nxUcAgQAMGABnzph56yIidmbri4/AXEk6ejRUV1tdiYhIZGo31EtLS0lLSyMlJYWCgoI299u+fTsx\nMTG8/fbb3Vrg5dSCERFpW8hQDwQCLFmyhNLSUvbs2UNRURF79+5tdb/HHnuMu+++u1vaLKEo1EVE\n2hYy1CsrK0lOTiYpKQm3201OTg7FxcVX7Pf73/+eefPmMWLEiLAV2kDTGkVE2hYT6ps1NTUkJiY2\nvvZ6vVRUVFyxT3FxMZs3b2b79u24XK5Wj5Wfn9/4PCsri6ysrC4VPGYM7NzZpR8VEYlo5eXllJeX\nX9UxQoZ6WwHd3NKlS1m5cmXjWdq22i/NQ/1qJCbCO+90y6FERCLK5QPe5cuXd/oYIUPd4/Hgb9bA\n9vv9eL3eFvt88skn5OTkAPDtt9+yceNG3G432dnZnS6mI9RTFxFpW8h56vX19aSmprJp0yZGjx7N\nTTfdRFFRET6fr9X9Fy5cyOzZs7n//vtbvkk3zVMH+O47GDsWTp3qlsOJiESsrmRnyJF6TEwMa9as\nYebMmQQCAXJzc/H5fBQWFgKQl5fX9Wq7aMgQqKuD77+H+Pgef3sRkYgWVVeUNkhLg7ffhvT0bjuk\niEjEsf0VpQ00rVFEpHVRGeparVFEpHVRGeqaASMi0jqFuoiIjURlqI8Zo566iEhrojLUNVIXEWld\nVE5pvHABhg2D06fB7e62w4qIRBTHTGns3x+uuQa++MLqSkREIktUhjrA5Mnw2WdWVyEiElkU6iIi\nNqJQFxGxkag8UQpw7JhZ++Xbb6EDy76LiEQdx5woBUhIgD59oLra6kpERCJH1IY6qAUjInK5qA71\nSZMU6iIizUV1qGukLiLSkkJdRMRGonb2C0AgAIMGweHDMHBgtx9eRMRSjpr9AtC7N0yYADt3Wl2J\niEhkiOpQB7VgRESaU6iLiNiIQl1ExEai+kQpwLlzMGKE1lYXEftx3IlSgNhYcyekL7+0uhIREetF\nfaiDWjAiIg0U6iIiNqJQFxGxEduE+uefQ/hP+YqIRDZbhPqoUdCrl1kuQETEyWwR6i6XWjAiImCT\nUAeFuogI2CjUdcMMEREbhbpG6iIiNlgmoEF9vVlb/ehRiI8P61uJiPQIRy4T0CAmBsaP19rqIuJs\ntgl1UAtGRKTdUC8tLSUtLY2UlBQKCgqu+H5xcTGTJk1iypQp3HjjjWzevDkshXaEQl1EnC5kTz0Q\nCJCamkpZWRkej4fMzEyKiorw+XyN+5w7d47Y2FgAdu3axX333cdXX33V8k16oKcOsHUrLF0KlZVh\nfysRkbDr9p56ZWUlycnJJCUl4Xa7ycnJobi4uMU+DYEOcPbsWYYPH96pArrT9dfD7t3mpKmIiBPF\nhPpmTU0NiYmJja+9Xi8VFRVX7Ld+/Xoef/xxjhw5wnvvvdfqsfLz8xufZ2VlkZWV1bWKQ4iPB48H\n9u2D9PRuP7yISFiVl5dTXl5+VccIGeoul6tDB5kzZw5z5sxhy5Yt/OxnP+PLVu5Y0TzUw6nhIiSF\nuohEm8sHvMuXL+/0MUK2XzweD36/v/G13+/H6/W2uf+0adOor6/nxIkTnS6ku+hkqYg4WchQz8jI\nYP/+/VRVVVFbW8u6devIzs5usc+BAwcaG/k7duwAYNiwYWEqt32TJ8P/lyEi4jgh2y8xMTGsWbOG\nmTNnEggEyM3NxefzUVhYCEBeXh5vvfUWr776Km63m7i4ONauXdsjhbdl6lTIyYGLF6FfP0tLERHp\ncbZZJqC5W26BFStg+vQee0sRkW7n6GUCmrvzTvif/7G6ChGRnmfLUL/rLmhjZqWIiK3Zsv1SVwfD\nh8NXX8GIET32tiIi3Urtl//ndsPtt8OmTVZXIiLSs2wZ6qC+uog4k+1DvQe7PiIilrNtqKemmsdW\nViwQEbEt24a6y6UWjIg4j21DHczURoW6iDiJLac0Nvj2Wxg3zjy63T3+9iIiV0VTGi8zfDgkJ8O2\nbVZXIiLSM2wd6qCrS0XEWWwf6jpZKiJOYuueOsClS2apgG++gSFDLClBRKRL1FNvRd++cNttsHmz\n1ZWIiISf7UMd1IIREedwRKhrvrqIOIUjQn38eLhwAQ4csLoSEZHwckSoa8kAEXEKR4Q6mFDXfHUR\nsTvbT2lscPQo+Hxw/DjExFhaiohIh2hKYwijRsGYMfDxx1ZXIiISPo4JdTAtmHfftboKEZHwcVSo\nz5kDf/iD7oYkIvblqFC/7TYztfHTT62uREQkPBwV6i4XLFgAr75qdSUiIuHhmNkvDQ4cgKlTobpa\nN84Qkcim2S8dMG4cpKTohKmI2JPjQh3UghER+3Jc+wXg5Em49lo4eFBrrItI5FL7pYOGDDErN77x\nhtWViIh0L0eGOqgFIyL25Mj2C0BdHXi98OGHkJxsdTUiIldS+6UT3G6YPx/+8z+trkREpPs4dqQO\nsGMHzJ1r5q73cuyvNxGJVBqpd9KUKRAba1owIiJ20G6ol5aWkpaWRkpKCgUFBVd8//XXX2fSpElM\nnDiR2267jZ07d4al0HDQsgEiYjch2y+BQIDU1FTKysrweDxkZmZSVFSEz+dr3Oejjz4iPT2dQYMG\nUVpaSn5+Ptu2bWv5JhHafgGoqYHrrzeP/ftbXY2ISJNub79UVlaSnJxMUlISbrebnJwciouLW+xz\n6623MmjQIABuvvlmqqurO1m2tTweyMiAkhKrKxERuXohb+xWU1NDYmJi42uv10tFRUWb+7/wwgvM\nmjWr1e/l5+c3Ps/KyiIrK6tzlYZRQwvmJz+xuhIRcbLy8nLKy8uv6hghQ93lcnX4QO+//z4vvvgi\nH7Zx1rF5qEea++6DRx4x9zEdNcrqakTEqS4f8C5fvrzTxwjZfvF4PPj9/sbXfr8fr9d7xX47d+5k\n8eLFlJSUMCQKF1OJjTV3RXr5ZasrERG5OiFPlNbX15OamsqmTZsYPXo0N9100xUnSg8dOsT06dN5\n7bXXuOWWW1p/kwg+Udpg506YORO+/lonTEUkMnQlO0O2X2JiYlizZg0zZ84kEAiQm5uLz+ejsLAQ\ngLy8PJ588klOnjzJgw8+CIDb7aaysrKL/wTrTJwImZnw0kvw0ENWVyMi0jWOvqL0ctu2QU4O7N+v\nuyKJiPV0RelVuuUWc2ekoiKrKxER6RqN1C+zaRMsWQK7d2s9GBGxlkbq3WD6dBg4EP77v62uRESk\n8xTql3G54J/+CVasgCj540JEpJFCvRX33gu1tfDee1ZXIiLSOQr1VvTqBY8/Dr/5jdWViIh0jkK9\nDT/+sVm5ccsWqysREek4hXobYmLgscfgX/7F6kpERDpOUxpDuHTJzFv/4x/NXZJERHqSpjR2s759\nYdkyMxNGRCQaaKTejrNnYexYKC+H9HSrqxERJ9FIPQzi4sy89Ucf1bx1EYl8CvUOePhhOHYM3njD\n6kpEREJT+6WDPvgA5s+HPXsgPt7qakTECbqSnQr1Tvj5z2HECFi1yupKRMQJFOphduwYTJhgTpqO\nH291NSJidzpRGmYJCfDEE2ZpXhv8jhIRG1Kod9KDD8Lp07qRhohEJrVfuuCjj2DePNi716y9LiIS\nDuqp96BFi8wsmKeftroSEbErhXoPOn7cnCwtK4OJE62uRkTsSCdKe9CIEfDkk+bCpL/+1epqREQM\nhfpVWLwYAgH43e+srkRExFD75SpVVcFNN0FJCdxyi9XViIidqP1igaQk+Pd/h5wc+O47q6sREafT\nSL2bLF0KBw/C+vXgclldjYjYgUbqFvrXf4XDh2H1aqsrEREn00i9G339temrv/MOZGZaXY2IRDuN\n1C02diw89xz85Cdw6pTV1YiIE2mkHgZLlsCRI/Dmm+qvi0jXaaQeIZ56ypw0ffZZqysREaeJsboA\nO+rXD/7wB5g6Fa67Du66y+qKRMQpNFIPk+RkeOst+Pu/h+3bra5GRJxCoR5G06bBf/wHZGfDvn1W\nVyMiTqD2S5j96EdmRceZM+HDD2H0aKsrEhE7a3ekXlpaSlpaGikpKRQUFFzx/S+++IJbb72Vfv36\n8dvf/jYsRUa7RYvM4l93362pjiISXiGnNAYCAVJTUykrK8Pj8ZCZmUlRURE+n69xn+PHj/PNN9+w\nfv16hgwZwrJly658E4dNaWxNMAiPPgqffQbvvgv9+1tdkYhEum6f0lhZWUlycjJJSUm43W5ycnIo\nLi5usc+IESPIyMjA7XZ3vmIHcbngmWfA44G/+zuor7e6IhGxo5A99ZqaGhITExtfe71eKioquvRG\n+fn5jc+zsrLIysrq0nGiWa9e8Mor8Ld/C//wD/D889C7t9VViUikKC8vp7y8/KqOETLUXd14OWTz\nUHeyPn3g7bdhzhyznMDrr0PfvlZXJSKR4PIB7/Llyzt9jJDtF4/Hg9/vb3zt9/vxer2dfhNpKT4e\nNmwwLZlZs+DMGasrEhG7CBnqGRkZ7N+/n6qqKmpra1m3bh3Z2dmt7uv0E6Gd1bcvrF0LKSnwwx/C\nX/5idUUiYgftLui1ceNGli5dSiAQIDc3l8cff5zCwkIA8vLyOHr0KJmZmZw5c4ZevXoRHx/Pnj17\niIuLa3oTzX5pUzAI+flQVATvvWfupCQiAl3LTq3SGCHWrIGCAtOWuf56q6sRkUjQlezUFaURYskS\nGD4cZswwS/ZOm2Z1RSISjbT2SwTJyYHXXoN588ycdv1xIyKdpfZLBKqqMsF+7bXwwgswcKDVFYmI\nFXSTDJtISoIPPjDtmIwM2LXL6opEJFoo1CNUv37mfqf//M8wfbq5ElVEpD1qv0SB3bth7lxz8vTf\n/k2LgYk4hdovNjV+vLl70tmzMGWKWZddRKQ1GqlHmbfegkceMSP3FSvMkgMiYk8aqTvA3Lnw5z/D\nuXPmIqV337W6IhGJJBqpR7H33oO8PPibv4Hf/Q6GDbO6IhHpThqpO8xdd5npjkOGmFH7Sy9BIGB1\nVSJiJY3UbaKyEpYtg9OnYdUqc6NrEYluWtDL4YJBKC6Gf/xHcwHTqlUwaZLVVYlIV6n94nAul7mj\n0u7d8KMfmdH6z38O1dVWVyYiPUWhbkNuNzz8MOzbZ250PWmSWQWyqsrqykQk3BTqNjZwIPzmN2bk\nHhcHN94ICxaY1yJiTwp1Bxg1ClauhAMHIC0N7rjDtGm2bbO6MhHpbjpR6kDnz5vpj6tWmROqDz9s\nevB9+lhdmYg0p9kv0il1dWbZgcJC+OILWLgQFi8267iLiPU0+0U6xe02d1t6/32zXbgAmZkwaxaU\nlEB9vdUVikhnaaQuLVy4AG+8YUbvBw+a0J8/39ysw+WyujoRZ1H7RbrVF19AURH813+Z1/Pnm83n\ns7YuEadQqEtYBIPwyScm3NeuhYQE+PGPYfZss9a7RvAi4aFQl7ALBOB//9ecYP3jH6F3bxPus2eb\n1SI1g0ak+yjUpUcFg2aVyJISE/D79pmVI++5x8yFT0y0ukKR6KZQF0sdOQLvvGPWed+8GYYONeE+\nYwb88IfmtYh0nEJdIsZf/wo7d0JZGWzaZO6rmpJiWjRTp8Jtt8Ho0VZXKRLZFOoSsWproaLChPuH\nH8LWrWZtmoaAnzrVnHR1u62uVCRyKNQlagSDpgffEPBbt8I338CECWbhsYZNQS9OplCXqPb99/DZ\nZ2b6ZMNWVWWC/frrzTZhgnlMSNBUSrE/hbrYztmzpje/a5fZ/vxn89irlwn4CRMgNbVp83rN90Ts\nQKEujhAMwtGjJtx374Yvv2zaTp82J2Svu85sY8eabdw4c2K2d2+rqxfpOIW6ON6ZM6ZX/+WXsH+/\nWb/m66/NduIEXHONCfmkJBgzpuU2erT69xJZFOoiIVy4YHr0Bw6Yk7KHDrXcjh0zvXqv1wT86NHm\ndoDNnyckwJAh6udLz1Coi1yFujqoqTHb4cNmu/z5sWPml8PIkSbgR40yjwkJMGIEDB9+5WNsrNX/\nMolWCvUoUF5eTlZWltVlRIRo/SwuXjTh3rAdPWoev/3WbMePNz0eP25+ZujQtrfBg+Hw4XKmTs1i\n8GAYNIjGxwEDnPdXQbT+dxEOXcnOmPZ2KC0tZenSpQQCARYtWsRjjz12xT6//OUv2bhxIwMGDODl\nl19mypQpnSrCSfQfbJNo/Sz69TO9+Wuu6dj+58/Dd9+13E6caHr+zTfwwQflbN+exalT5mTvyZPm\n/EBtLcTHm4AfONBs8fEtt7i4lo+xsW1vAwaYRdci+RdFtP53ESlChnogEGDJkiWUlZXh8XjIzMwk\nOzsbX7MFtTds2MBXX33F/v37qaio4MEHH2Sb7mgs0mjAALN5vW3vk59vtsvV1Zn5+2fONG2nT5up\nnt9/3/TY8Mvh7Fk4d67p8fLtwgVzR6sBA6B//6ba+vdv2vr1a/nYvz/07WteN2wNr/v2NVufPk3P\nm7/u06dpa/7a7dbU03AJGeqVlZUkJyeTlJQEQE5ODsXFxS1CvaSkhAceeACAm2++mVOnTnHs2DES\nEhLCV7WIQ7jdTW2a7lJfb8L9/HmzNTy/cMG0li5cuPL5pUvm9Zkz8Je/mOcXL5qvX7pk/qJoeN6w\n1dU1fb22tmlr+F7v3ubf53Y3Bb3bbX4hrV1rnsfENH294Xnzx8uf9+7d8rG1rzV/bG/r1evK5x15\n7OzmcrX+vCtChnpNTQ2JzdZP9Xq9VFRUtLtPdXX1FaHuiuS/93rY8uXLrS4hYuizaOK0zyIQMNvF\ni1d+7+RJZ30W3SlkqHc0iC9v5F/+czpJKiLSM0IO8D0eD36/v/G13+/He1lj8PJ9qqur8Xg83Vym\niIh0RMhQz8jIYP/+/VRVVVFbW8u6devIzs5usU92djavvvoqANu2bWPw4MHqp4uIWCRk+yUmJoY1\na9Ywc+ZMAoEAubm5+Hw+CgsLAcjLy2PWrFls2LCB5ORkYmNjeemll3qkcBERaUUwzDZu3BhMTU0N\nJicnB1euXBnut4tYhw4dCmZlZQXT09OD48ePD65evdrqkixXX18fnDx5cvDee++1uhRLnTx5Mjh3\n7txgWlpa0OfzBT/66COrS7LMihUrgunp6cEJEyYE58+fH7x48aLVJfWYhQsXBkeOHBmcMGFC49dO\nnDgRnDFjRjAlJSV45513Bk+ePNnuccI6U7RhnntpaSl79uyhqKiIvXv3hvMtI5bb7ebpp59m9+7d\nbNu2jWeffdaxn0WD1atXk56e7viZUY8++iizZs1i79697Ny5s8WUYSepqqri+eefZ8eOHezatYtA\nIMDatWutLqvHLFy4kNLS0hZfW7lyJXfeeSf79u3jjjvuYOXKle0eJ6yh3nyeu9vtbpzn7kSjRo1i\n8uTJAMTFxeHz+Th8+LDFVVmnurqaDRs2sGjRIkfPjjp9+jRbtmzhF7/4BWBanoMGDbK4KmsMHDgQ\nt9vN+fPnqa+v5/z5846adDFt2jSGDBnS4mvNrwN64IEHWL9+fbvHCWuotzaHvaamJpxvGRWqqqr4\n9NNPufnmm60uxTK/+tWvWLVqFb0cflnhwYMHGTFiBAsXLuSGG25g8eLFnD9/3uqyLDF06FCWLVvG\nmDFjGD16NIMHD2bGjBlWl2Wp5hdyJiQkcOzYsXZ/Jqz/Rzn9z+rWnD17lnnz5rF69Wri4uKsLscS\nf/rTnxg5ciRTpkxx9CgdoL6+nh07dvDQQw+xY8cOYmNjO/Qnth0dOHCAZ555hqqqKg4fPszZs2d5\n/fXXrS4rYrhcrg5lalhDvSPz3J2krq6OuXPn8tOf/pQ5c+ZYXY5ltm7dSklJCddeey3z589n8+bN\nLFiwwOqyLOH1evF6vWRmZgIwb948duzYYXFV1vj444+ZOnUqw4YNIyYmhvvvv5+tW7daXZalEhIS\nOHr0KABHjhxh5MiR7f5MWEO9I/PcnSIYDJKbm0t6ejpLly61uhxLrVixAr/fz8GDB1m7di3Tp09v\nvNbBaUaNGkViYiL79u0DoKysjPHjx1tclTXS0tLYtm0bFy5cIBgMUlZWRnp6utVlWSo7O5tXXnkF\ngFdeeaVjg8FwTc9psGHDhuB1110XHDduXHDFihXhfruItWXLlqDL5QpOmjQpOHny5ODkyZODGzdu\ntLosy5WXlwdnz55tdRmW+uyzz4IZGRnBiRMnBu+7777gqVOnrC7JMgUFBY1TGhcsWBCsra21uqQe\nk5OTE/zBD34QdLvdQa/XG3zxxReDJ06cCN5xxx2dmtLYIzfJEBGRnuHsqQciIjajUBcRsRGFuoiI\njSjURURsRKEuImIjCnURERv5P2zBHLSWHEybAAAAAElFTkSuQmCC\n"
}
],
"prompt_number": 139
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"v=norm.rvs(size=10000)\n",
"v2=v*v\n",
"h,bin_edges=np.histogram(v2,bins=100,density=True)\n",
"bin_center=(bin_edges[1:]+bin_edges[:-1])/2.\n",
"bar(bin_center,h,width=0.2);plot(x,chi2pdf,'r-');xlim([0,10]);ylim([0,0.5])\u0000"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "pyout",
"prompt_number": 151,
"text": [
"(0, 0.5)"
]
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAD9CAYAAABDaefJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAH7RJREFUeJzt3Xt4VPW97/H3hAy33Ei4BJkJBEkkCchFklJB2xQEWlpS\nRNqGtgdLg5vHSivWs7fH9rQ7eFo37LZb2eBxc9it1WIDrVhCK0SMmBZFEhUVBCoRCUyCpAIhFy65\nTNb5IxcyyWRyIZM1l8/redaTrFm/mfmaRz755bt+s5bFMAwDEREJCCFmFyAiIn1HoS4iEkAU6iIi\nAUShLiISQBTqIiIBRKEuIhJAugz1vLw8kpKSSExMZP369R2OFxQUEBUVxfTp05k+fTo/+9nPvFKo\niIh0LdTTQafTyerVq8nPz8dms5GWlkZGRgbJycku4z7/+c+za9curxYqIiJd8zhTLyoqIiEhgfj4\neKxWK5mZmeTm5nYYp88viYj4Bo8z9bKyMuLi4lr37XY7hYWFLmMsFgsHDhxg6tSp2Gw2fvnLX5KS\nktJhjIiI9FxPJ80eZ+rdCePbbrsNh8PB+++/z/e//30WL17cRWFGu639Y9f3v8Z2XmCJyzHDMPx6\n+9d//VfTa/CVTT8L/Sz0s/C89YbHULfZbDgcjtZ9h8OB3W53GRMREcHQoUMB+NKXvkR9fT0XL17s\nVTHtOYgjDkfXA0VEBOgi1FNTUykuLqakpIS6ujq2b99ORkaGy5jy8vLW3yhFRUUYhkFMTEyfFHeG\nsQp1EZEe8NhTDw0NZdOmTSxYsACn00lWVhbJycls3rwZgFWrVvHCCy/w9NNPExoaytChQ9m2bVuf\nFXeO0QznAgOppY5Bffa6ZkpPTze7BJ+hn8V1+llcp5/FjbEYvW3c9ORNLBYMw2ju0bd/u/aPue6X\nMI4v8BqnuBmw9LrPJCLib1qysyd8/hOlasGIiHSfz4e6TpaKiHSfX4T6WM6YXYaIiF/wi1DXTF1E\npHt8PtTVUxcR6T6fD3W1X0REus8vQl0zdRGR7vH5UL/AcAZRSxg1ZpciIuLzfD7UwaLZuohIN/lB\nqKuvLiLSXX4T6pqpi4h0zeMFvXxDDGeoII7fAq7XeI+IiKaqqm8u8ysiEgj8YKZegYMtjOU7tL/B\nRnV1hamViYj4Gj8IdbVfRES6yy9CXZ8qFRHpHr8I9eszdV1LXUTEE78I9cuEU8sghnPB7FJERHya\nX4Q6qAUjItIdfhPqOlkqItI1vwp1fapURMQzvwp1zdRFRDzzm1BXT11EpGt+E+qaqYuIdM2vQl09\ndRERzyyGYXj9Ez0WiwXDMJovxtX+7do/5n5/ILVUEclQrtDIgNZj/VC+iIgpWrKzJ/xmpl7HIC4S\nw2jOmV2KiIjP8ptQB7VgRES64nehrpOlIiKd86tQ17JGERHP/CrUNVMXEfHM70JdPXURkc75Vair\n/SIi4plfhbraLyIinvlVqJcTSzQVDKTW7FJERHySX4V6IwM4yxjslJpdioiIT+oy1PPy8khKSiIx\nMZH169d3Ou6tt94iNDSUF198sU8LbE8tGBGRznkMdafTyerVq8nLy+PYsWPk5ORw/Phxt+MeeeQR\nvvjFL3r9WiwKdRGRznkM9aKiIhISEoiPj8dqtZKZmUlubm6HcRs3bmTp0qWMHDnSa4W20LJGEZHO\nhXo6WFZWRlxcXOu+3W6nsLCww5jc3Fz27dvHW2+91Xwlxo6ys7NbvgPSm7eeO8NYpnC4V88VEfFl\nBQUFFBQU3NBreAz1zgK6rTVr1rBu3brWS0R21n7Jzs5m7dq1NIV67zmI48u8dEOvISLii9LT00lP\nT2/db8rMnvEY6jabDYfjev/a4XBgt9tdxrzzzjtkZmYCcP78efbs2YPVaiUjI6PHxXSHeuoiIp3z\neJOMhoYGJk6cyKuvvsqYMWP4zGc+Q05ODsnJyW7Hr1ixgkWLFrFkyRLXN+mDm2S0iOYiHzOC6A6v\n0yQiIpqqqoud/SeJiPiN3twkw+NMPTQ0lE2bNrFgwQKcTidZWVkkJyezefNmAFatWtX7anupgmis\nGIRTRQ0RHY5XV3fdMhIRCVR+czu7to5jYQlHOU6Ku3fTLe5EJCAE9O3s2nKAljWKiLjhl6F+BnSy\nVETEDb8MdQcKdRERdxTqIiIBxC9D/QzqqYuIuOOXoa6ZuoiIe34Z6qcBO6WEUm92KSIiPsUvQ/0a\ncJpxJPF3s0sREfEpfhnqAO8xjWm8Z3YZIiI+RaEuIhJAFOoiIgEkAEJd13kREWnht6H+D2KpYyB2\nSs0uRUTEZ/htqINaMCIi7fl1qL/PVIW6iEgbfh3qmqmLiLhSqIuIBBC/vPNRy2MhOKkkijGcpZrI\n1mO685GIBIKgufNRi0YG8AGTmcJhs0sREfEJfh3qoBaMiEhbCnURkQCiUBcRCSB+faIUYCiX+ZSR\nRFFJA1Z0olREAkXQnSgFuEIYDuKYyIdmlyIiYjq/D3VQC0ZEpIVCXUQkgCjURUQCSMCE+lTeR9dW\nF5FgFxChfo7RNBLCGM6aXYqIiKkCItTBohaMiAgBE+rqq4uIQACFum6YISICoWYX0DdieI8K1jbv\nNX1ytUlERDRVVRfNKUtEpJ/5/WUCWvYHUE8lUYzmHDVEuBzTZQNExB8F5WUCWjgJ5SiTdG11EQlq\nARPqoJOlIiJdhnpeXh5JSUkkJiayfv36Dsdzc3OZOnUq06dPZ8aMGezbt88rhXaHQl1Egp3HnrrT\n6WTixInk5+djs9lIS0sjJyeH5OTk1jGXL18mLCwMgCNHjnD33Xfz0Ucfub5JP/TUweB2DvAka5hJ\nkcsx9dRFxB/1eU+9qKiIhIQE4uPjsVqtZGZmkpub6zKmJdABampqGDFiRI8K6EtHuJVJHGUADabV\nICJiJo9LGsvKyoiLi2vdt9vtFBYWdhi3c+dOHn30UT755BP27t3r9rWys7NbvgPSm7e+VUMEZdi4\nhRMcJ6XPX19ExJsKCgooKCi4odfw2H7ZsWMHeXl5bNmyBYCtW7dSWFjIxo0b3Y7fv38/K1eu5MMP\nXW9Y0V/tF4A/8DX+xN3k8M3WY2q/iIg/6vP2i81mw+FwtO47HA7sdnun4++8804aGhq4cOFCj4ro\nSzpZKiLBzGOop6amUlxcTElJCXV1dWzfvp2MjAyXMSdPnmz9TXLo0CEAhg8f7qVyu/Ye07iNQ6a9\nv4iImTz21ENDQ9m0aRMLFizA6XSSlZVFcnIymzdvBmDVqlXs2LGD5557DqvVSnh4ONu2beuXwjtz\ngFlsI5NBRFPLJUCXDRCR4BEwlwlou/8mn+VHFPKa25tmqMcuIv4hqC8T0NYrzGOe2UWIiJggIEN9\nL/OZb3YRIiImCMj2Syj1nGcgCfyD84zs8Fy1X0TEH6j90qwBK38F5vKq2aWIiPSrgAx1gFeAebxi\ndhkiIv0qCEJdrRYRCR4BG+otFyqYyIcex4mIBJKADXVoWdqoFoyIBI+ADvW9zFeoi0hQCcgljS2P\nDedTTjKBEZynAWvr41rSKCL+QEsa27nACD4igc9y0OxSRET6RUCHOrR8utT9jTtERAJNwIe6TpaK\nSDAJ6J46GAyklk8ZyThOc4lo1FMXEX+hnrobdQziDWYzh31mlyIi4nUBH+qgFoyIBI+gCHWtVxeR\nYBEUoX6USQzhKjdz0uxSRES8KihCHSxqwYhIUAiSUG/qq2u9uogEuoBf0tgilnMcJ5mRXKJBSxpF\nxA9oSaMH5YzmDGNJNbsQEREvCppQh6YWzAKzixAR8aKgCvWdLObrAGq/iEiACqpQf4PZDAF4912z\nSxER8YogCvUYIITngCdnzMBisbRukZExZhcnItIngmb1S8v+zZzkALOwU6obZ4iIT9Pql274mAkU\nk8gCXja7FBGRPhd0oQ7wHMtZznNmlyEi0ueCrv0CMIwKTjGe8ZzSNdZFxGep/dJNl4hmL/P5Gn80\nuxQRkT4VlKEOasGISGAK2lB/mQUkUswEPjK7FBGRPhO0od6AlRwq+R8kAmjduogEhKA8UdpiOhZ2\nEM8ETmK4/H7TiVMRMZ9OlPbQu8BlwpjNG2aXIiLSJ7oM9by8PJKSkkhMTGT9+vUdjj///PNMnTqV\nKVOmMHv2bA4fPuyVQr1FJ0xFJJB4bL84nU4mTpxIfn4+NpuNtLQ0cnJySE5Obh3z5ptvkpKSQlRU\nFHl5eWRnZ3Pw4EHXN/HR9gtYGEMpR7gVG2Vca7rcF2q/iIgv6PP2S1FREQkJCcTHx2O1WsnMzCQ3\nN9dlzO23305UVBQAM2fOpLS0tIdlm+ssNt4mlQx2mV2KiMgNC/V0sKysjLi4uNZ9u91OYWFhp+N/\n/etfs3DhQrfHsrOzW74D0ps339DSgvkD3zC7FBEJYgUFBRQUFNzQa3hsv+zYsYO8vDy2bNkCwNat\nWyksLGTjxo0dxr722ms88MADvPHGG0RHR7u+iQ+3X8BgKJcpxU4yBuVUdvhvA4iIiKaq6qLbYyIi\n3tDn7RebzYbD4Wjddzgc2O32DuMOHz7Mfffdx65duzoEuj+4Qhg7Wcx3qKQp+Dtu1dUVZpYoItIt\nHkM9NTWV4uJiSkpKqKurY/v27WRkZLiMOXPmDEuWLGHr1q0kJCR4tVhveoKHeBAYzNV2R2JomtHr\nA0oi4vs89tRDQ0PZtGkTCxYswOl0kpWVRXJyMps3bwZg1apVPPbYY1RUVHD//fcDYLVaKSoq8n7l\nfewIU3gLWMEzPM332hypoGPbBqqrLf1VmohItwX1J0rbPzYTC9sYRyLFLndFchfqWvYoIt6mT5Te\noELgJBNYRo7ZpYiI9Ipm6u3GzCGfTaxmEkebrwejmbqImEMz9T6wjzlUEcnd/MnsUkREekyh3oGF\nn/NjfsTjuJ+hi4j4LoW6G3/hKwykjvnsNbsUEZEeUai7YRDCv/EoP+bnZpciItIjCvVO/IGvY6OM\nO8wuRESkBxTqnXASynoe4VGzCxER6QEtafQwZiC1nGQwizjEe0zv8FwtaRQRb9KSxj5WxyB+Bc0r\nYUREfJ9CvQtbgM/zV5I5ZnYpIiJdUqh34TLwc37MBh5E69ZFxNcp1LvhKR4glnK+xh/NLkVExCOd\nKO3mmNm8Tg7LSOEYNUSgE6Ui4m06UepFb3AHrzKXnzISdzfN0I0zRMQXaKbegzGjKOcDRpPOBxxj\nEh1p9i4ifUczdS/7B7GsBTaxGp00FRFfpFDvoaeBKCp1Iw0R8Ulqv/RizGc5wAssJZnjVBPpckzt\nFxHpK71pvyjUezlmCyupJoIf8oTLMYW6iPQVhXo/jhnBpxxlEneRzxGmtB5TqItIX9GJ0n50npH8\nlMd4igew0Gh2OSIigEL9hmzhPgbg5If8h9mliIgAar/c8JhxRFFEFRlAYbtnRkREU1V1ERGR3lD7\nxQSnqeKf+BPbGEc0F2gK/KaturrC5OpEJNhopn7DY5r2n2AN4znFYnbSchkBnTgVkRuhmbqJ/oV/\nZwxneZANZpciIkFMM/U+mqkDjOdjDvJZvsxLvE0amqmLyI3QTN1kp7iZ+3ma7XyDKC6ZXY6IBCHN\n1Ptwpt5iI6u5iU9YyouaqYtIr2mm7iP+J79kPKd4oHk/MjKmw7XXdQ12EfGGULMLCES1DObr/IED\nJMLevc1LG93/tq2utrh9XESkN9R+8UL7pcUdWNg/YgRp58/zdqfXX9fJVBFxT+0XH/M6wH//N7uA\nRE6YXI2IBAOFurd99av8BHiZBdzEWbOrEZEA12Wo5+XlkZSURGJiIuvXr+9w/O9//zu33347gwcP\n5le/+pVXivQ/MbR8qtRisfBrmi7+lccXtdRRRLzKY0/d6XQyceJE8vPzsdlspKWlkZOTQ3JycuuY\nTz/9lNOnT7Nz506io6N5+OGHO75J0PXU3Y1pZAMPMo33WMDLXGNI6zH11EXEnT7vqRcVFZGQkEB8\nfDxWq5XMzExyc3NdxowcOZLU1FSsVmvPKw4qFtbwJGXY+D3fZAANZhckIgHI45LGsrIy4uLiWvft\ndjuFhe0vMNs92dnZLd8B6c1bcDEI4V6e5SW+zP/jn7iPLbq9hoi0KigooKCg4IZew2OoN7VL+kZ2\ndjZr166lKdSDVz0DWcKL7CSG7TzDt3D9Oesa7CLBKz09nfT09Nb9pszsGY/tF5vNhsPhaN13OBzY\n7fYev4m4qiGChTRgsJTdzCGCSnQNdhHpCx5DPTU1leLiYkpKSqirq2P79u1kZGS4HauTfT1TB2Sy\njWISeY0vMJJ/mF2SiASALj9RumfPHtasWYPT6SQrK4tHH32UzZs3A7Bq1SrOnTtHWloaVVVVhISE\nEBERwbFjxwgPD7/+Jlr94mGMQTbZLONnzKeR03SkloxIcOrN6hddJsD0UG/yABYewc5CdvMBt3YY\nq7+ERIKPLhPgx54C/plfkM9d3MF+s8sRET+lUPch28nk22zlBZbyIE/S2YXCREQ6o/aLj7Rf2u6P\no4QXWMopDpNFPdXtnqkeu0hwUPslQJwmnjt4nfPU8zaJTOYwLSdVtexRRDxRqPuoWgbzPeAxfso+\n5rCcZ80uSUT8gNovPth+af9YCkfZwT3s505+wH9yjaFaDSMSBNR+CVDHmEQabxFODe8ynVlmFyQi\nPkszdT+YqV8XwxIq2AjsAH4E1DQf0clTkcCjmXrAq+BFDCZzgTBWcIRxzKfpk7vV1RVYLJbWLTIy\nBoDIyBiXx92NEZHAoZm6X83UXffnsZfNLOBvLOeH/AcXGe4ytvOfuesYEfFNmqkHmVeYz61ABdEc\n4Va+wzOE4DS7LBExkULdz10GHuJJFrOT7/Ib3mU683m5i2e53kNVLRmRwKH2ix+3Xzo+ZvBVIvh3\nLlMC/DNwuPlxz6/rekwtGRHfoPZL0LOQy2UmUUcum3iZWJ4BbJSaXZiI9BOFegBqwMr/5QFu4QRl\nwPtMZSOrGUeJ2aWJiJcp1ANYNZH8b2ASR6khnHeYwbMsJ8XswkTEaxTqQaCc0TzKOiZwkr+TxKvA\nn1jMTA6aXZqI9DGFehCpZBj/xo+4GXiFeWwjk9dIZynhWLHguiLGqg8sifghrX4JqNUvPXtOKPXc\nww5WsYwkRvMMK9jCfZQwvpPnNT1Xq2NE+odWv0iPNGBlO5nMAb7AawzhKm+RxkssZBEwgAazSxSR\nHtJMPYhn6u72B3OVr/FHVnEv47mJbWSSwzLeZj5wqd3rhIKb4NfFxUT6hmbqcsOuMYTfsZw7gDns\no5oIfs83OcEl1vITkjhG0y8Ag6ZAN9ps0UDHi4upFy/SfzRT10y9G2MMZhDCN3mITLZRTix/4Ov8\nmR9xlEZaTrDqk6oifas3M3WFuk+Epv/UF4KTz/E37mEHi3gKJ+P5M4v4M4v4G/OoV6iL9BmFuilj\ngru+W3mfDHaxiD9zC0Xs5evs4Uu8ylxKiXMZq1AX6RmFuiljVF+L0Vj4MluYz17msI+LxPAqc8nn\nLl5jKRcV6iI9olA3ZYzqc7dvoZEpHOYu8pnLq8wmj8jbboPPfQ5mzYLZs2HMGESkcwp1U8aovs73\nY4AKAKzATGA2IcymkVlAFXAAeKP565nwYVysrkBEmijUTRmj+nr3HINbOMFs3mAWB5jFAcZxnLDP\nfAZmzLi+TZoEVivtRUbGUO3mF0DbNfLdGSPiyxTqpoxRfTf2nOuz+XBgGjCjeUsLCSFp0CDeqa3j\n/UYnR4APgCNAObh57abXb/lfuvP7s+qkrfgHhbopY1Rf37/u9aAPA6YAt/Jf3MoRJvMBt3KERi7y\nAel8wGQ+ZGLrVso4GhXqEiAU6qaMUX39/7oGownhVl5mEkfbRPqHRPEJYVOnwi238LM//pGP+TUf\nczMnmcBZxtDIAFpCvbP2DHhu0aitI/1FoW7KGNXnG+/dNLuPAG4BJgKJwHiWczMfczMfM5yznAY+\nBkqAM8AZtnKGsZxhLGcZQwPW1tfv7J+G/gKQ/tKbUA/1Ui0i/awCMKgG3mnemoL/2dYRg7EQzzEm\ncJJxnGYsq5nCX5oj/QyxlFNOLKXYOQvw/e+Dzda09HLMmKbvY2PdvPf1dlFT4F+n2bv0N83U/Wo2\n6m/1+dfPJpR6bIzARhVjgDGArfnrGCDOEkJiZAS1lZX8AzvlxHKO0ZTzEuU8yqeM5Dwj2n2N53If\n/hPrbctI/JNm6n6hAEg3uQZfUWB2AS4asHKaKk67ba0AhgXj0iWiLBZieZ1YyomlnNG8RCxDGcdp\nZvAOI/mUEbzKSJyMBK5aLFwELgKVA0K5Y9FXICaGDb97nnP1tVwCTgADCeESjVTSdJHjSgZwBae7\nQtyWV11tcfu4vykoKCA9Pd3sMvxWlzP1vLw81qxZg9PpZOXKlTzyyCMdxvzgBz9gz549DB06lN/+\n9rdMnz7d9U00U2+zn928+Wp9/fne2cBaH64P2rZWXHX/dYZwhRjGEEMlMc2vOLz5awz/QhSV5PM6\nKznKMGYSRSXRVBBJOQMZQDURVBJFFZFUcYRqvkg1Ec3b89RQSzVQA1QDl5s3Y0g4L7/+VwgLu74N\nHQoDB0K7NpGZ9NdH5/p8pu50Olm9ejX5+fnYbDbS0tLIyMggOTm5dczu3bv56KOPKC4uprCwkPvv\nv5+DB3VDYwkUTb16Vz0LxKsMpYxKyty+zvrm77PZwVFwuRm4hVCuEkE1kc2RHslUolhNODVEUE04\nvyGCx4ihhnFUE04NYWwjnHrCrtbw7owZhEHzZmHYQCs0NFDd2MhV4ErzdrXNdg24ioVrGK2P1QLX\nmh+71rofQi2N1AJ1zY/VMoBanK37dQygzmW/6bF6nO1+Gm33spu36399dAx+9zdogZ79IgjEXyge\nQ72oqIiEhATi4+MByMzMJDc31yXUd+3axb333gvAzJkzuXTpEuXl5cS6PaEkIj3RgJUKYqig7U1G\nvtzm+/uAn7R71u9w/9eGAXV1DACGAEM5x1CuMISrDGUSQ/gbg7nGEK4yhK8ymN82f3+VQdQymEeJ\n5H8ximsM5hqD+S8GsYxB1DKQOgZRyyBeYRAzGUQtVuoZyFEGMZaB1LVug6hhINDAAOqxUs816hjR\n/L2VX1HD/fwf6mmkASiyWNgDNPA56rHSQCj1vEwDi5u/b3qsge04qaOhuoJNFgsNWGjAwElT/Dd9\nDcHZ/LpOIAtwsgEnAzpu1d+B3/8eBgxo2kJCWr/e841vcvlqDU6gEXAygEacrfuNzfuNrfvuHru+\nb7R5zGgzpjc8hnpZWRlxcdcvn2q32yksLOxyTGlpaYdQv74qwN0sp/1jXe372piePmdtH753d8b4\n08+mv8cES33XOWlq1dQwut2Rz7Xb/46bZ69rt5/jZkxhu/0zHippOWdwweXIf7od/7d2+zs7ed0W\n7toW7qLywU5f4flvfauL92jR/tyHu3MhvRnTcx5Dvf3yrM607/m0f57W7oqI9A+P9yi12Ww4HI7W\nfYfDgd1u9zimtLQUm83Wx2WKiEh3eAz11NRUiouLKSkpoa6uju3bt5ORkeEyJiMjg+eeew6AgwcP\nMmzYMPXTRURM4rH9EhoayqZNm1iwYAFOp5OsrCySk5PZvHkzAKtWrWLhwoXs3r2bhIQEwsLCeOaZ\nZ/qlcBERccPwsj179hgTJ040EhISjHXr1nn77XzWmTNnjPT0dCMlJcWYNGmSsWHDBrNLMl1DQ4Mx\nbdo04ytf+YrZpZiqoqLCuOeee4ykpCQjOTnZePPNN80uyTSPP/64kZKSYkyePNlYtmyZce3aNbNL\n6jcrVqwwRo0aZUyePLn1sQsXLhh33XWXkZiYaMybN8+oqKjo8nU8tl9uVMs697y8PI4dO0ZOTg7H\njx/35lv6LKvVyhNPPMHRo0c5ePAgTz31VND+LFps2LCBlJSUbp+QD1QPPvggCxcu5Pjx4xw+fNhl\nyXAwKSkpYcuWLRw6dIgjR47gdDrZtm2b2WX1mxUrVpCXl+fy2Lp165g3bx4nTpxg7ty5rFvXfvVR\nR14N9bbr3K1Wa+s692A0evRopk2bBkB4eDjJycmcPXvW5KrMU1payu7du1m5cmVQr46qrKxk//79\nfPe73wWaWp5RUVEmV2WOyMhIrFYrV65coaGhgStXrgTVoos777yT6Ohol8fafg7o3nvvZefOrpZx\nejnU3a1hLysr8+Zb+oWSkhLeffddZs6caXYppnnooYf4xS9+QUiIV/8X9HmnTp1i5MiRrFixgttu\nu4377ruPK1eumF2WKWJiYnj44YcZO3YsY8aMYdiwYdx1111ml2Wqth/kjI2Npby8vMvnePVfVLD/\nWe1OTU0NS5cuZcOGDYSHh5tdjin+8pe/MGrUKKZPnx7Us3SAhoYGDh06xPe+9z0OHTpEWFhYt/7E\nDkQnT57kySefpKSkhLNnz1JTU8Pzzz9vdlk+w2KxdCtTvRrq3VnnHkzq6+u55557+Pa3v83ixYvN\nLsc0Bw4cYNeuXYwfP55ly5axb98+li9fbnZZprDb7djtdtLS0gBYunQphw4dMrkqc7z99tvMmjWL\n4cOHExoaypIlSzhw4IDZZZkqNjaWc+fOAfDJJ58watSoLp/j1VDvzjr3YGEYBllZWaSkpLBmzRqz\nyzHV448/jsPh4NSpU2zbto05c+a0ftYh2IwePZq4uDhOnDgBQH5+PpMmTTK5KnMkJSVx8OBBrl69\nimEY5Ofnk5KSYnZZpsrIyODZZ58F4Nlnn+3eZNBby3Na7N6927jllluMCRMmGI8//ri3385n7d+/\n37BYLMbUqVONadOmGdOmTTP27NljdlmmKygoMBYtWmR2GaZ67733jNTUVGPKlCnG3XffbVy6dMns\nkkyzfv361iWNy5cvN+rq6swuqd9kZmYaN910k2G1Wg273W785je/MS5cuGDMnTu3R0sa++XORyIi\n0j+Ce+mBiEiAUaiLiAQQhbqISABRqIuIBBCFuohIAFGoi4gEkP8P76hwT95k9a8AAAAASUVORK5C\nYII=\n"
}
],
"prompt_number": 151
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"<h2> If $v = v_1 + v_2 i$ where $v_1$ and $v_2$ are normally distributed random variable, then $|v|^2$ is $\\chi^2$ distributed with two degree of freedom </h2>"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"df=2\n",
"rv=chi2(df)\n",
"x=np.linspace(0,10)\n",
"chi2pdf=rv.pdf(x)\n",
"plot(x,chi2pdf),ylim([0,0.5])"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "pyout",
"prompt_number": 161,
"text": [
"([<matplotlib.lines.Line2D at 0xed6becc>], (0, 0.5))"
]
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAD9CAYAAABDaefJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHn9JREFUeJzt3Xt41NWdx/F3ICNyE8IdZgJBEslE5GYiClVSELFUstws\nYVdBDDa1C0prt162uwWfloK2q2i6GllvFAsoVkJdiDalUbZAgqQCCtUoRCYBETGgATQwzP5xmoRA\nGHKbnLl8Xs9znpnJ/PL7fZ1HvnPy/Z1LlM/n8yEiImGhle0ARESk+Sipi4iEESV1EZEwoqQuIhJG\nlNRFRMKIkrqISBi5aFLPzc0lMTGRhIQElixZct77+fn5dOrUiWHDhjFs2DB+8YtfBCRQERG5uGh/\nb3q9XubOnUteXh5Op5OUlBTS0tJwu921jhs9ejTr1q0LaKAiInJxfnvqhYWFxMfHExcXh8PhID09\nnZycnPOO0/wlEZHg4LenXlZWRmxsbPVrl8tFQUFBrWOioqLYvHkzQ4YMwel08utf/5qkpKTzjhER\nkYZraKfZb0+9Psl4+PDheDweduzYwbx585g0aVKdx3372z58PrWf//zn1mMIlqbPQp+FPgv/rTH8\nJnWn04nH46l+7fF4cLlctY7p2LEj7dq1A+A73/kOp06d4osvvjjvXNu3w5EjjYpRRETqyW9ST05O\npri4mJKSEiorK1m9ejVpaWm1jjl06FD1N0phYSE+n48uXbqcd64bb4Q//rEZIxcRkfP4ralHR0eT\nlZXF+PHj8Xq9ZGRk4Ha7yc7OBiAzM5M1a9bw1FNPER0dTbt27Vi1alWd55o0CV59Fe64o9n/G0JK\namqq7RCChj6LGvosauizaJooX2MLNw25SFQUX3zhIy4ODhyA9u0DfUURkdAXFRXV4Np6i80ojYmB\na66BN95oqSuKiESeFl0mYPJkeO21lryiiEhkabHyi8/no6wMrroKDh0ChyPQVxURCW1BXX4BcDoh\nIQHeeqslryoiEjlafJVGlWBERAKnRcsvAB98AGPGgMcDrbTwr4jIBQV9+QVg4EDo1Am2bWvpK4uI\nhD8rfeVJk1SCEREJBCtJffJkWLvWxpVFRMKblaSenAzHj8OePTauLiISvqwk9agolWBERALB2vgT\nDW0UEWl+LT6kscqpU9CrF7z7Lpy1uZKIiPxDSAxprOJwwC23QB1bnoqISCNZnf6jEoyISPOyVn4B\nOHECeveGvXuha9dARyEiElpCqvwC0K6d2eZOY9ZFRJqH9dVXZsyAlSttRyEiEh6sll8ATp6EPn3M\nRKRevQIdiYhI6Ai58gtA27aQlgYvv2w7EhGR0Gc9qYMpwfz+97ajEBEJfdbLL2AmIjmdsHUrXH55\noKMREQkNIVl+ATMR6dZbdcNURKSpgiKpQ00JJvB/N4iIhK+gSeojR0JFBezaZTsSEZHQFTRJvVUr\nSE9XCUZEpCmC4kZplR074J/+CfbtM2uui4hEspC9UVpl8GCzdMCWLbYjEREJTUGV1KOi4J//WWPW\nRUQaK6jKLwAffQSjRkFZGURHBzgwEZEgFvLlF4D4eOjXD/78Z9uRiIiEnqBL6mBKMBoFIyLScEFX\nfgE4cACuvNI8tm0bwMBERIJYWJRfwCzFO3w4rF9vOxIRkdASlEkdtHmGiEhjXDSp5+bmkpiYSEJC\nAkuWLLngcdu2bSM6Opo//OEPzRLY1Knwpz/B0aPNcjoRkYjgN6l7vV7mzp1Lbm4uu3fvZuXKlezZ\ns6fO4+6//35uvvnmBtd/LiQmBsaNg1WrmuV0IiIRwW9SLywsJD4+nri4OBwOB+np6eTk5Jx33JNP\nPsm0adPo3r17swaXkQHPPtuspxQRCWt+p/eUlZURGxtb/drlclFQUHDeMTk5OWzcuJFt27YRdYFF\nWxYsWFD9PDU1ldTU1IsGd9NN8P3vw86dZgkBEZFwlp+fT35+fpPO4TepXyhBn23+/PksXry4eujN\nhcovZyf1+mrdGu64w/TWly5t8K+LiISUczu8CxcubPA5/CZ1p9OJx+Opfu3xeHC5XLWO2b59O+np\n6QB8/vnnbNiwAYfDQVpaWoODqcvs2TBiBDzyCLRp0yynFBEJW34nH50+fZqBAwfy5z//mT59+nDN\nNdewcuVK3G53ncfPnj2biRMnMmXKlNoXacQA+rONHWvKMNOnN/oUIiIhp9knH0VHR5OVlcX48eNJ\nSkpi+vTpuN1usrOzyc7OblKwDaEbpiIi9ROUywSc6+RJcLmgqMgs9iUiEgnCZpmAc7Vtaxb5ev55\n25GIiAS3kOipA7z7rtnqbu9eMypGRCTchW1PHWDoUOjWTeusi4j4EzJJHXTDVETkYkKm/AJmca+4\nOLPlXbduTY9LRCSYhXX5BaBzZ7jlFlixwnYkIiLBKaSSOtSUYAL/94WISOgJuaQ+ejScOAHbttmO\nREQk+IRcUm/VCu68UzdMRUTqElI3SquUlcFVV8Enn0DHjs12WhGRoBL2N0qrOJ1mka8XX7QdiYhI\ncAnJnjrApk1w112we7cpyYiIhJuI6akDfOtbcOmlZnNqERExQjapR0XBPffAE0/YjkREJHiEbPkF\nzJK8/frBX/8KCQnNfnoREasiqvwCZknejAz47W9tRyIiEhxCuqcOsH+/WcFRwxtFJNxEXE8doG9f\nDW8UEakS8j11MMMb58yBPXs0vFFEwkdE9tTBDG9s2xbefNN2JCIidoVFUq8a3vjkk7YjERGxKyzK\nL6DhjSISfiK2/AKm/DJnDmRl2Y5ERMSesOmpA3g8MGSIhjeKSHiI6J46QGyshjeKSGQLq546mOGN\ns2fDBx9A69YtckkRkYCI+J46mOGNPXvCmjW2IxERaXlhl9SjouDBB+FXv9Lm1CISecIuqQN897sm\noW/YYDsSEZGWFZZJPSoKHngAFi2yHYmISMsKy6QOcOut8Omn5sapiEikCNukHh0NP/2peusiElnC\nbkjj2b75Bi6/HF5/HYYNa/HLi4g0iYY0nqNNG/jxj81IGBGRSBDWPXWAigro398s9HXFFVZCEBFp\nFPXU69ChA/zrv8KSJbYjEREJvIsm9dzcXBITE0lISGBJHZkxJyeHIUOGMGzYMK6++mo2btwYkECb\nYt48eO01s+CXiEg481t+8Xq9DBw4kLy8PJxOJykpKaxcuRK32119zPHjx2nfvj0Au3btYvLkyXz0\n0Ue1L2Kx/FLlvvvA64XHH7cahohIvTV7+aWwsJD4+Hji4uJwOBykp6eTk5NT65iqhA5QUVFBt27d\nGhRAS/nxj2H5cjh82HYkIiKBE+3vzbKyMmJjY6tfu1wuCgoKzjtu7dq1PPjggxw8eJA3L7BR6IIF\nC6qfp6amkpqa2riIG8nphGnTTE/9l79s0UuLiNRLfn4++fn5TTqH3/LLq6++Sm5uLsuWLQNgxYoV\nFBQU8OQFNgPdtGkTc+bM4YMPPqh9kSAovwDs2wfJybBnD/ToYTsaERH/mr384nQ68Zx1d9Hj8eBy\nuS54/PXXX8/p06c5cuRIg4JoKf37w4wZGrcuIuHLb1JPTk6muLiYkpISKisrWb16NWlpabWO+fjj\nj6u/SYqKigDo2rVrgMJtup/9zNTW9++3HYmISPPzW1OPjo4mKyuL8ePH4/V6ycjIwO12k52dDUBm\nZiavvvoqy5cvx+Fw0KFDB1atWtUigTdWr16QmQkPPwz/8z+2oxERaV5hP6O0LuXlZnbp//0fDBxo\nOxoRkbppRmk9xcSYIY7/+Z+2IxERaV4R2VMHOH4c4uNh/Xqt4CgiwUk99QZo3x7+/d9NExEJFxGb\n1AHuusuMWdfuSCISLiI6qbdpAwsWwEMPmY2qRURCXUQndYDbboMjRyA313YkIiJNF/FJvXVr+MUv\nTG/9zBnb0YiINE3EJ3WAyZPNRtUvv2w7EhGRponYIY3nevttuP12c+O0XTvb0YiIaEhjk9xwA1x7\nLTzyiO1IREQaTz31s+zfbyYiFRVBv362oxGRSKeeehP17Qv33AP/9m+2IxERaRz11M9x4gQkJcEL\nL0ALb84kIlKLeurNoF07ePRRuPdeOH3adjQiIg2jpF6HadPMSo7/2MVPRCRkqPxyATt3wrhxZohj\nly62oxGRSNSY3Kmk7scPf2hmnF5gn20RkYBSUm9mR46A2w0bN8KgQbajEZFIoxulzaxrV7M70r33\nahVHEQkNSuoX8YMfwGefwSuv2I5EROTiVH6phy1bYOpUeO893TQVkZajmnoA3XMPfPUVPP+87UhE\nJFIoqQdQRYW5WfrMM3DTTbajEZFIoBulAdShAzz9NGRmmgQvIhKM1FNvoJkzzaiYxx6zHYmIhDuV\nX1rAkSOmDLN2LYwYYTsaEQlnKr+0gKpe+pw5UFlpOxoRkdqU1Bth+nSIi4PFi21HIiJSm8ovjVRa\nanZJeusts/66iEhzU/mlBblc8PDDpgzj9dqORkTEUFJvgsxMaNPGbKohIhIMVH5pIo8HkpPh9dch\nJcV2NCISTlR+sSA2FrKy4F/+RZOSRMQ+9dSbyZ13QlQUPPus7UhEJFyop27RE0/A22/DmjW2IxGR\nSKaeejMqLISJE+Gdd0xZRkSkKQLSU8/NzSUxMZGEhASWLFly3vsvvfQSQ4YMYfDgwYwaNYqdO3c2\nKIBwcs01ZpekmTM1zFFE7PDbU/d6vQwcOJC8vDycTicpKSmsXLkSt9tdfcyWLVtISkqiU6dO5Obm\nsmDBArZu3Vr7IhHSUweTzMeMge98Bx54wHY0IhLKmr2nXlhYSHx8PHFxcTgcDtLT08nJyal1zHXX\nXUenTp0AGDFiBKWlpQ0MO7y0bg2/+x3813+ZMoyISEuK9vdmWVkZsWcVh10uFwUFBRc8/tlnn2XC\nhAl1vrdgwYLq56mpqaSmpjYs0hDSt68Z5jhjhkns//jOExHxKz8/n/z8/Cadw29Sj4qKqveJ/vKX\nv/Dcc8/x17/+tc73z07qkeB73zPrwsyaBX/4A7TSOCMRuYhzO7wLFy5s8Dn8phqn04nH46l+7fF4\ncLlc5x23c+dO7rrrLtatW0dMTEyDgwhXjz0Gn30GixbZjkREIoXfpJ6cnExxcTElJSVUVlayevVq\n0tLSah2zf/9+pkyZwooVK4iPjw9osKHmkkvMuPWnnoING2xHIyKRwG/5JTo6mqysLMaPH4/X6yUj\nIwO32012djYAmZmZPPzww5SXl3P33XcD4HA4KCwsDHzkIaJPH1i9GqZOhc2bYcAA2xGJSDjT5KMW\nkpUFy5aZxN6+ve1oRCQUaI/SIObzmZumXi+sWGHWiRER8UdrvwSxqCjIzobdu2HpUtvRiEi4Uk+9\nhe3bB9deCy+/DKNH245GRIKZeuohoH9/M+N0+nQoLrYdjYiEGyV1C266yexvOmECHD5sOxoRCScq\nv1j00EPwl7/Axo3Qtq3taEQk2Gj0S4g5cwZuvx1OnoRXXjGLgYmIVFFNPcS0agXPPQfl5fCTn9iO\nRkTCgZK6ZW3amAW/3nwTHn/cdjQiEur8LhMgLSMmBtavh1GjzLK9U6bYjkhEQpWSepDo1w/WrYPx\n46FXLxg50nZEIhKKVH4JIsOHmzHskydDUZHtaEQkFCmpB5mbb4annzZj2Hftsh2NiIQalV+C0OTJ\nUFlpSjEbN0Jiou2IRCRUKKkHqenT4Ztv4MYbIT8ftP+IiNSHknoQmzkTvv4axo41+53GxdmOSESC\nnZJ6kPv+902PvSqx17FFrIhINSX1EDBvnknsY8aYxN67t+2IRCRYKamHiJ/8xNw8veEG+NOfVIoR\nkbopqYeQhx6Cjh3h+uvhjTcgKcl2RCISbJTUQ8y8eWZZgTFj4I9/hJQU2xGJSDDR5KMQdNttsGwZ\nfPe7Zj12EZEqSuohauJEs8/p9OmQk2M7GhEJFiq/hLDUVLO648SJcOyYGdcuIpFNST3EJSebpQTG\nj4cDB+D++yEqynZUImKLtrMLE2VlkJYGgwbBM8+YzTdEJLRpO7sI5nTC22/DV1+Z9WI+/9x2RCJi\ng5J6GGnfHtasMePYR4yA3bttRyQiLU3llzC1fLmZhfq735l6u4iEHpVfpNrMmWZD6zvugKws0Heq\nSGRQTz3M7d1rbqBefTX893+bEo2IhAb11OU8l18OBQXm+YgRsGeP3XhEJLCU1CNA+/bwwgvwox+Z\nVR5//3vbEYlIoKj8EmF27IBp08ymG48/DpdeajsiEbkQlV/kooYMge3b4cgRGDkSPv7YdkQi0pwu\nmtRzc3NJTEwkISGBJUuWnPf+3//+d6677jouvfRSfvOb3wQkSGlel11mFgObPRuuvRZefFGjY0TC\nhd/yi9frZeDAgeTl5eF0OklJSWHlypW43e7qYw4fPswnn3zC2rVriYmJ4b777jv/Iiq/BK0dO+D2\n280N1WeegR49bEckIlWavfxSWFhIfHw8cXFxOBwO0tPTyTlnndfu3buTnJyMw+FoeMRi3ZAhsG0b\nuN3m+Wuv2Y5IRJrC7yqNZWVlxMbGVr92uVwUVI2Pa6AFCxZUP09NTSU1NbVR55Hm16YN/OpXZgnf\nmTNh7VpYuhQ6d7YdmUhkyc/PJz8/v0nn8JvUo5pxDdezk7oEp5Ej4d134ac/hcGD4dlnYdw421GJ\nRI5zO7wLFy5s8Dn8ll+cTicej6f6tcfjweVyNfgiEjo6dDAzT5ctgzlzTL39s89sRyUi9eU3qScn\nJ1NcXExJSQmVlZWsXr2atLS0Oo/VjdDwMn48vP8+9OpVs0b7mTO2oxKRi7no5KMNGzYwf/58vF4v\nGRkZPPjgg2RnZwOQmZnJp59+SkpKCl9++SWtWrWiY8eO7N69mw4dOtRcRKNfQtrOnZCZaXZUevpp\nU5oRkcBrTO7UjFKplzNnTEnmZz8z49t//nMtDiYSaJpRKgHTqpXprb/3ntkLdeBAM2lJJRmR4KKe\nujTK5s1w333wzTfwm9/At79tOyKR8KPyi7Qonw9eeQUeeMDcTH30UdODF5HmofKLtKioKPje98wa\n7TfcAN/6FsydC4cP245MJHIpqUuTtWlj9kPds8fU3hMTTe/9889tRyYSeZTUpdl06wZPPGFmpR47\nZkoxDz1klvkVkZahpC7NLjYWnnoKiopMb/2KK+A//gPKy21HJhL+lNQlYPr1MzNR33nHDINMSIAH\nH4SDB21HJhK+lNQl4Pr3N4uDFRZCRQVceSVkZMDu3bYjEwk/SurSYi6/HJ58EoqLIS4OxoyBW26B\nt97SzksizUXj1MWakydh+XIzealzZ5g3D269VZthi1TR5CMJSWfOwOuvmyV/t2+HO+6AH/wABgyw\nHZmIXZp8JCGpVStIS4PcXNiyxfzs2mvh5pshJwdOn7Ybn0goUU9dgtLXX5slCJ56Cjwes83erFlm\neKRIpFD5RcLSrl1mRciXXjI3WGfNgunTISbGdmQigaWkLmHt9Gl44w2T4N980+zONHOm2Uf1kkts\nRyfS/JTUJWKUl8Pq1bBihVlzJi3NLC42dqwSvIQPJXWJSKWl8Oqr8PLL8Pe/w6RJJsGPGQMOh+3o\nRBpPSV0inscDa9aYXvyHH5oRNGlp5rFzZ9vRiTSMkrrIWQ4cMOPf162Dt9+GlBST4CdONLNbRYKd\nkrrIBRw/Dnl5JsH/7/+akTM33WTa6NHQoYPtCEXOp6QuUg9nzpg1399807Rt2yA5uSbJDx0KrVvb\njlJESV2kUSoqzKJiVUn+0CGzNd/o0ZCaqiQv9iipizSDTz81Sb6qlZXBqFEmyY8aBVdfrUXHpGUo\nqYsEwGefmRutb78NmzebcfGDBsF119W02FizEbdIc1JSF2kBJ06Y3Zy2bDFt82YzHj452bSrrzat\nZ0/bkUqoU1IXscDng337zLLB77xjHrdvNyNqkpNh+HAYMsS0vn3Vo5f6U1IXCRI+H+zda5J8URHs\n3Ak7dphe/uDBJsEPHgxXXQVuN3TqZDtiCUZK6iJB7vDhmgS/Ywe8955Z2qBLF0hKqt0SE6FrV9sR\ni01K6iIh6MwZ+OQTsxH37t3mRuz775tkHx1t1pA/tw0YoAlTkUBJXSSM+HymZ//hh+e3vXuhY0ez\n3MGAAeaxqsXFQZ8+5gtBQpuSukiE8PnMePqPPzYJfu/emuclJfD559C7t7kx269fTXO5alrnzrpp\nG+yU1EUEgMpKs2Ll/v2mtFPVSkvNZKrSUrPpSFWCdzpN7753b9POft6une3/msilpC4i9fbllzUJ\nvrQUDh407cCBmucHD0KbNtCjhxl3f27r3t20bt3MY5cuWlKhOSmpi0iz8vng6FEzq/bQofPb4cOm\n1FP1ePSoKet062Zaly5mBE/VY9XzmJja7bLL9GVQFyX1EJCfn09qaqrtMIKCPosa4fJZeL3wxRcm\nyR85YtoXX5z/vLzctKNHzeNXX5kbvzEx0Lp1Pk5nKp06cV677LLarWPHmseOHc2aPOF0n6AxufOi\n98dzc3OZP38+Xq+XOXPmcP/99593zD333MOGDRto164dL7zwAsOGDWtQEJEkXP7xNgd9FjXC5bNo\n3bqmJNMQXq8pB5WXwyOP5JOensqxYybpHztm2mefmZvBX35Z0776quZ5RQWcOmWGenboYJJ81fP2\n7Wvaua/btau7tW1b06pet2kT/F8afpO61+tl7ty55OXl4XQ6SUlJIS0tDbfbXX3M+vXr+eijjygu\nLqagoIC7776brVu3BjxwEQkfrVvXlGJ69TJLHjfGqVNmQ5SKCpPwv/rKvK5qFRW1n5eXm1m+57bj\nx+HkydrtxAlz/ksvrWlt29Z+XdXatDGt6vnZP7vkkgs/ntsaw29SLywsJD4+nri4OADS09PJycmp\nldTXrVvHrFmzABgxYgRHjx7l0KFD9NRqRiLSwhwOU9MP1H60Xi98/XXd7eRJ8/jNNxd+rKw0xx07\nZl5XtVOnzHtnt2++aVyMfpN6WVkZsbGx1a9dLhcFBQUXPaa0tPS8pB4V7H+ztKCFCxfaDiFo6LOo\noc+ihj6LxvOb1OubiM8t5J/7e7pJKiLSMlr5e9PpdOLxeKpfezweXC6X32NKS0txOp3NHKaIiNSH\n36SenJxMcXExJSUlVFZWsnr1atLS0modk5aWxvLlywHYunUrnTt3Vj1dRMQSv+WX6OhosrKyGD9+\nPF6vl4yMDNxuN9nZ2QBkZmYyYcIE1q9fT3x8PO3bt+f5559vkcBFRKQOvgDbsGGDb+DAgb74+Hjf\n4sWLA325oLV//35famqqLykpyXfllVf6li5dajsk606fPu0bOnSo75ZbbrEdilXl5eW+qVOn+hIT\nE31ut9u3ZcsW2yFZs2jRIl9SUpJv0KBBvhkzZvi+/vpr2yG1mNmzZ/t69OjhGzRoUPXPjhw54rvx\nxht9CQkJvnHjxvnKy8sveh6/5Zemqhrnnpuby+7du1m5ciV79uwJ5CWDlsPh4LHHHuP9999n69at\n/Pa3v43Yz6LK0qVLSUpKiviRUffeey8TJkxgz5497Ny5s9aQ4UhSUlLCsmXLKCoqYteuXXi9Xlat\nWmU7rBYze/ZscnNza/1s8eLFjBs3jg8//JCxY8eyePHii54noEn97HHuDoejepx7JOrVqxdDhw4F\noEOHDrjdbg4cOGA5KntKS0tZv349c+bMiejRUceOHWPTpk3ceeedgCl5dorQve0uu+wyHA4HJ06c\n4PTp05w4cSKiBl1cf/31xMTE1PrZ2fOAZs2axdq1ay96noAm9brGsJeVlQXykiGhpKSEv/3tb4wY\nMcJ2KNb86Ec/4tFHH6VVq4D+Lxj09u3bR/fu3Zk9ezbDhw/nrrvu4sSJE7bDsqJLly7cd9999O3b\nlz59+tC5c2duvPFG22FZdfZEzp49e3Lo0KGL/k5A/0VF+p/VdamoqGDatGksXbqUDhG6H9nrr79O\njx49GDZsWET30gFOnz5NUVERP/zhDykqKqJ9+/b1+hM7HH388cc8/vjjlJSUcODAASoqKnjppZds\nhxU0oqKi6pVTA5rU6zPOPZKcOnWKqVOncttttzFp0iTb4VizefNm1q1bR//+/ZkxYwYbN25k5syZ\ntsOywuVy4XK5SElJAWDatGkUFRVZjsqOd955h5EjR9K1a1eio6OZMmUKmzdvth2WVT179uTTTz8F\n4ODBg/To0eOivxPQpF6fce6RwufzkZGRQVJSEvPnz7cdjlWLFi3C4/Gwb98+Vq1axZgxY6rnOkSa\nXr16ERsby4cffghAXl4eV155peWo7EhMTGTr1q2cPHkSn89HXl4eSUlJtsOyKi0tjRdffBGAF198\nsX6dwUANz6myfv163xVXXOEbMGCAb9GiRYG+XNDatGmTLyoqyjdkyBDf0KFDfUOHDvVt2LDBdljW\n5efn+yZOnGg7DKveffddX3Jysm/w4MG+yZMn+44ePWo7JGuWLFlSPaRx5syZvsrKStshtZj09HRf\n7969fQ6Hw+dyuXzPPfec78iRI76xY8c2aEhji2ySISIiLSOyhx6IiIQZJXURkTCipC4iEkaU1EVE\nwoiSuohIGFFSFxEJI/8PI2WVMPSfEvkAAAAASUVORK5CYII=\n"
}
],
"prompt_number": 161
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"v1=norm.rvs(size=10000);v2=norm.rvs(size=10000)\n",
"v=v1+v2*1j\n",
"vabs=np.abs(v)**2\n",
"h,bin_edges=np.histogram(vabs,bins=100,density=True)\n",
"bin_center=(bin_edges[1:]+bin_edges[:-1])/2.\n",
"bar(bin_center,h,width=0.2);plot(x,chi2pdf,'r-');xlim([0,10]);ylim([0,0.5])"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "pyout",
"prompt_number": 167,
"text": [
"(0, 0.5)"
]
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAD9CAYAAABDaefJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VPXZ//H3QIYtCZCwykw0YAKZgGwmIiiSAoqlmgKi\nxj6K0kBxQcXapxbrY9GfRbC1iuBjU8QFsYCFSrA/iDZglAokLApUECIQmYSlLAECAUKG8/yREBIy\n2cgkZ5bP67rOFebMd+5zm8vc+eY+33OOxTAMAxER8QtNzE5AREQ8R0VdRMSPqKiLiPgRFXURET+i\noi4i4kdU1EVE/EiNRT0tLY2YmBiio6OZOXNmpfczMjJo06YN/fr1o1+/frz00ksNkqiIiNQsqLo3\nXS4XkydPJj09HZvNRnx8PImJiTgcjgrjhgwZwvLlyxs0URERqVm1M/WsrCyioqKIjIzEarWSlJRE\nampqpXG6fklExDtUO1PPy8sjIiKi7LXdbiczM7PCGIvFwtq1a+nTpw82m40//vGPxMbGVhojIiJ1\nV9dJc7Uz9doU4/79++N0OtmyZQuPP/44o0aNcp/Yj36EYRgBv/3ud78zPQdv2fS90PdC34vqtytR\nbVG32Ww4nc6y106nE7vdXmFMaGgorVq1AuDHP/4x58+f59ixY5WDbdoER49eUZIiIlI71Rb1uLg4\nsrOzycnJoaioiMWLF5OYmFhhzKFDh8p+o2RlZWEYBuHh4ZWDDR8On3ziucxFRKSSanvqQUFBzJkz\nhxEjRuByuUhOTsbhcJCSkgLApEmTWLJkCW+99RZBQUG0atWKRYsWuQ82ahQsXQoPPeTp/wafkpCQ\nYHYKXkPfi0v0vbhE34v6sRhX2ripy0EsFoxjxyAyEvbvh+Dghj6kiIjPs1gsde6tN94VpWFhcMMN\n8OmnjXZIEZFA07i3CRg9Gj7+uFEPKSISSBqv/WIYkJcH110Hhw6B1drQhxUR8Wne3X4BsNkgOhq+\n+KJRDysiEiga/y6NasGIiDSYxm2/AOzcCUOHgtMJTXTnXxGRqnh/+wWgRw9o0wY2bGj0Q4uI+Dtz\npsqjRqkFIyLSAMwp6qNHw7JlphxaRMSfmVPU4+Lg9GnYscOUw4uI+KtGP1HaunU4BQX5vAEcAF4u\nNy40NIyTJ93c4VFEJAD5xInSgoJ8wOBjVjGaOMAo20reExGRK9XoM/WSB28YBHGeg3SmL9+Qy8Wn\nK9X9t5KIiL/yiZn6RcVY+Qd38FMqP/NURESujKlX/3zMaEajpY0iIp5iWvsFoCWFHOAqurGHY7RD\n7RcRkUt8qv0CcIZWpDOcUWjNuoiIJ5h+85WF3Md9LARKnmtqsVgqbK1bu3neqYiIuGVq+wWgBWfY\nTxccHOcQ7lJRS0ZEApPPtV8AztKS5SRyj9mJiIj4AdOLOpS0YH5mdhIiIn7A9PYLQBDnyaMZN7Kb\nvXS7/NNqv4hIQPLJ9guUXIj0Nyg9YSoiIlfKK4o6wELgZ/yVy2fxIiJSe15T1NcCIZziOraZnYqI\niM/ymqJuAItIUgtGRKQevOJEaekoevMNqfyUruwFLGX7daJURAKRz54ovWgrvSmkFQNZZ3YqIiI+\nyauKOlj4Kz8rPWEqIiJ15VXtFzC4lu/5ipuwkYeLINR+EZFA5fPtF4DdRPED1zCMVWanIiLic7yu\nqAP8lZ9Ve+dG3b1RRMQ9r2u/AFzFfr6lJ104ztkqL0ZSW0ZE/JtftF8ADtCFzfRnpNmJiIj4GK8s\n6nDx4RkiIlIXNRb1tLQ0YmJiiI6OZubMmVWO27BhA0FBQfz973/3SGJLuYtbgTYc90g8EZFAUG1R\nd7lcTJ48mbS0NLZv387ChQvZsWOH23HPPPMMt99+u8f63McJ459AEos8Ek9EJBBUW9SzsrKIiooi\nMjISq9VKUlISqamplcbNnj2bsWPH0qFDB48mNw9IZp5HY4qI+LOg6t7My8sjIiKi7LXdbiczM7PS\nmNTUVFavXs2GDRtKV7dUNm3atHKvMoCEGpP7DPgLB7mOrWyjd43jRUR8WUZGBhkZGfWKUW1Rr6pA\nlzdlyhRmzJhRtvSmqvbLxaL+wgsvUJuCDnABeI+HSGYeU5hVq8+IiPiqhIQEEhISyl6X1Mu6qbao\n22w2nE5n2Wun04ndbq8wZtOmTSQlJQFw5MgRVq5cidVqJTExsc7JuPMu48lkAL/mFYpo7pGYIiL+\nqtqLj4qLi+nRowerVq2iS5cu3HDDDSxcuBCHw+F2/Pjx47nzzjsZM2ZMxYPU8eKjy/elM4y/8As+\n4t4K7+niIxHxZx6/+CgoKIg5c+YwYsQIYmNjuffee3E4HKSkpJCSklKvZOtiHsk6YSoiUgteeZuA\ny/e14Ay52OnPZvZxTdl7mqmLiD/zm9sEXO4sLfkrP2M875qdioiIV/OJmTpAn9JH3XVjDxdoimbq\nIuLv/HamDrCFvhyhve6zLiJSDZ8p6qATpiIiNfGZ9guU3Nwrh0ii+J6jdFD7RUT8ml+3XwBO0JZ/\ncAf3s8DsVEREvJJPzdQBEvicNxhOby64PVZoaBgnTx6rb8oiIqbz+5k6wBcMoRUXiCeTkoJfcSso\nyDc1PxERM/lcUTdowjvolrwiIu74XPsFoAsWthHGNfzAKUIrfUYnUEXEHwRE+wVgP7CKYTzI+2an\nIiLiVXyyqAO8wRNMZg6WKk6YiogEIp8t6v/iZs7Sglv5p9mpiIh4DZ8t6mDhDZ7gCd4wOxEREa/h\nkydKy9+S9weu4Sa+4nuiy97TiVIR8QcBc6L0orO0ZB7JPMabZqciIuIVfHqmDhDBPr6hb7nljZqp\ni4h/CLiZOoCTq7W8UUSklM/P1AFuZg1vMwEHOzBoqpm6iPiFgJypQ8nyxjO05DY+MzsVERFT+UVR\nv7i88XFmm52IiIip/KL9ApRb3niYbLVfRMQPBGz7BUqWN77NBCabnYiIiIn8ZqYOYMfJFq4m/ORJ\nCL387o0iIr4loGfqALlEsArgfS1vFJHA5FczdYCbsbDm2mth505o2vSKcxYRMVvAz9QB/gXQqRMs\nWWJ2KiIijc7vijoAU6fCyy+DVsGISIDxz6L+k5+UFPSVK83ORESkUflnUbdY4De/genTzc5ERKRR\nBZmdgGeFAyUnF5oC3wHjLRb+RRBQXGl0aGgYJ08ea9QMRUQakp/N1PMpWRVj4MLgFVJ4ltspKehG\npa2gIN+8VEVEGoCfLWmsuL8Z59hDN+5gP99UMV53dBQRb6UljZcpojl/4pdMNTsREZFG4tczdYBg\nTrGXUG5iJ9l0rzReM3UR8VaaqbtxmhDeBJ5hptmpiIg0uBqLelpaGjExMURHRzNzZuXCmJqaSp8+\nfejXrx/XX389q1evbpBE62M2MJqPseM0OxURkQZVbfvF5XLRo0cP0tPTsdlsxMfHs3DhQhwOR9mY\n06dPExwcDMC2bdsYPXo033//fcWDmNh+ubj/j/ySprh4itcr7Ff7RUS8lcfbL1lZWURFRREZGYnV\naiUpKYnU1NQKYy4WdIBTp07Rvn37OiXQWP7ELxnHfNpz2OxUREQaTLUXH+Xl5REREVH22m63k5mZ\nWWncsmXLmDp1KgcOHOCzz9w/J3TatGnlXmUACXXPth72Y2MJY5nC6zzH7xv12CIitZGRkUFGRka9\nYlTbflm6dClpaWnMnTsXgAULFpCZmcns2e6fBbpmzRomTJjAzp07Kx7EC9ovYBDJXjYSh4MdHKYj\nar+IiDfzePvFZrPhdF46ueh0OrHb7VWOHzx4MMXFxRw9erROSTSWHLqykPuYystmpyIi0iCqLepx\ncXFkZ2eTk5NDUVERixcvJjExscKY3bt3l/0m2bx5MwDt2rVroHTr7yWeYxzziWCf2amIiHhctT31\noKAg5syZw4gRI3C5XCQnJ+NwOEhJSQFg0qRJLF26lPnz52O1WgkJCWHRokWNkviVOkRnUpjE87zI\nRLOTERHxML+/otTd/rbks4vu3MwRdqqnLiJeSleU1tJxwvgTp3iRkm/a5Vvr1uFmpygickUCsqgD\nvMFZbqEzfdmMbskrIv4iYIt6IfB7fsvv+a3ZqYiIeEzAFnWAuUzEwQ5uZo3ZqYiIeERAF/UimjON\naUznWdyfYBUR8S0BXdQBFnA/7TjK7aSZnYqISL0FfFG/QFOe4yWm8ywWLpidjohIvQR8UQf4mNEU\nE8Q9fGR2KiIi9RKQFx+52zeYL/mAB3CwgzME60ZfImI6XXxUD2u4hfXcyK95xexURESumGbq5USw\nj6/pR3+O8YNm6iJiMs3U68nJ1bzBE/zB7ERERK6QZuqXaUkh2wkm8vPPISHBTSwRkcahmboHnKEV\n/w3w5JNQXGx2OiIidaKi7sYSgLAwKH2Mn4iIr1D7pZJwIJ/rgH8CDuDiPRtDQ8M4efKYm8+IiHje\nlbRfVNSrGfsmj+KiKU8wu+w9rV8Xkcaiou7hoh7OUXbgYCir+ZZeqKiLSGPSiVIPO0Y7XuR5ZvEk\nuoujiPgCFfUa/JmH6ch/uJu/mZ2KiEiN1H6pxdgbWcdS7qIXBzim9ouINBK1XxrIegayhLG8anYi\nIiI10Ey9ljGCOcW/CSXy00/httvcfE5ExLO0+qUBizrACCykRUbCtm0QEuJ2jIiIp6j90sA+BRg8\nGP7nf8xORUTELc3Uax2j5D3jyBHo1QuWLYMBA6oYJyJSf5qpN6hwACzt25N08CDbbryRZhYLFouF\n1q3DTc5NRKSEinqt5VMygzdYzAVyuIPf8AJgUFCQX8NnRUQah9ovVxjDRi5f048hfMEOeur2ASLi\ncWq/NKI87DzPi7zNBH0TRcRraKZejxgWLrCK5nxKMTMvG63b9IpIfWmm3sgMmjCOYp6iI3FkcbHn\nrj67iJhFRb2ecoHJzOFD/otgTpmdjogEOLVfPBR7Hj/HwMIE5pXt18lTEakPtV9M9ARvcAtfclfJ\nE05FREyhmboHY8eTxSfcSRwbyeVqzdRFpF4aZKaelpZGTEwM0dHRzJx5+RoP+PDDD+nTpw+9e/fm\npptuYuvWrXVKwJ9s4AZm8STzGac/gUTEFNXO1F0uFz169CA9PR2bzUZ8fDwLFy7E4XCUjVm3bh2x\nsbG0adOGtLQ0pk2bxvr16yseJEBm6gBNcLGaoazkS2Zopi4i9eDxmXpWVhZRUVFERkZitVpJSkoi\nNTW1wpiBAwfSpk0bAAYMGEBubm4d0/YvF2jKA2zhl0Bc6b1hym+6T4yINKSg6t7My8sjIiKi7LXd\nbiczM7PK8fPmzWPkyJFu35s2bVq5VxlAQu2z9DFOTjCZxSzkt8SxkZO0KXuvoMBiYmYi4s0yMjLI\nyMioV4xq2y9Lly4lLS2NuXPnArBgwQIyMzOZPXt2pbGff/45jz32GF999RVhYWEVDxJA7Zfy++bw\nGDbyGMPfMcr+KNJSRxGpHY+3X2w2G06ns+y10+nEbrdXGrd161YmTpzI8uXLKxX0QPYUr9GR//As\n081ORUQCRLVFPS4ujuzsbHJycigqKmLx4sUkJiZWGLNv3z7GjBnDggULiIqKatBkfc15mjGWJTzC\nW9zOSrPTEZEAUG1PPSgoiDlz5jBixAhcLhfJyck4HA5SUlIAmDRpEi+++CL5+fk88sgjAFitVrKy\nsho+cx9xgC7cy2KWcheDWMsesxMSEb+mi48aKfZjzGEicxnEVk6rpy4itaDbBHixN3mMLfThLwAq\n6iLSQFTUG42FSXxCLDClSROtXxeRBqGi3ojOcpwx7GEqHbmFDHT/dRHxNBX1RpZDVx7gAxZzL1Fk\nm52OiPgZFXUT/JPbeJ4XWcFI2nPY7HRExI+oqJtkLr/gb9zNchJpwRmz0xERP6EljSbGtnCBD3iA\nlpzhbj7GpVUxIlKOljT6GIMm/Jx3CCOfP5qdjIj4BRV1kxXRnDF8zW3AFN2qV0TqSUXdCxznBCPJ\n4VfYGM1StNRRRK5Utfd+kcazj2tIZDmfMoKDdGYdg8xOSUR8kGbqXuRr+vMAH/Axo+nHZrPTEREf\npKLuZT7ldh7mz6xgJL3YZnY6IuJj1H7xQssYTTOK+JQ+DOXiMtBLQkPDOHnymDnJiYhX0zp1L479\nABamYyOBDHYTVWG8Hokn4v+0Tt3PfAC8yPOsYhjXkGN2OiLiA9R+8XJz+QXNOccqhjGEL8ij8jNi\nRUQuUlH3AXN4nOacYzVDGcIXHDQ7IRHxWirqPuJVfkUzivgSO7dS+eQp6ASqiOhEqc/FnoyFX2Nn\nBJ+yg9hK43UCVcR/6ERpAJgDTOVlVjOUODaYnY6IeBkVdR/0Ifczkbn8f35CAp+X7i258dflNwSz\nWKxu9ulmYSL+Su0XH449hAw+4h4mMpfljKpj7JL31K4R8V5X0n7RiVIf9gUJjGQFn3AnbShZ1y4i\ngU3tFx+3iTiGspqXgGeYQdWzchEJBCrqfuA7HAwE7uZvvMdDNOOc2SmJiElU1P3EfuAWviSUAtIZ\nTjuOmJ2SiJhARd2PFBLMWJawhsFkMgAH281OSUQamU6U+hmDJvyW6XxHDBkk8ADwmdlJiUij0Uzd\nT33AOMbwd94DHmMOOoEqEhi0Tt3PY3fFwnJ6sonreZT/pZDgCp/ROnUR76XbBEgle4EBZAKQyQBi\n2FH6jvsrUHWVqYhvU1EPAIUE8xDv8RpP8SW3cB9/BfIpmdmX38IoKMjXLQVEfJjaLwEWuzdbWMJY\nVvE9UzjDOVrUIkbJe2rViDQutV+kRlvpw/Vsoh2wlkF0Y7fZKYmIB9VY1NPS0oiJiSE6OpqZM2dW\nev+7775j4MCBtGjRgldffbVBkhTPKqA19wDvMp713Mg43kerY0T8Q7XtF5fLRY8ePUhPT8dmsxEf\nH8/ChQtxOBxlYw4fPswPP/zAsmXLCAsL4+mnn658ELVfvDZ2b7bwAQ+wh278glQOq/0i4jU83n7J\nysoiKiqKyMhIrFYrSUlJpKamVhjToUMH4uLisFqtdc9YTLeVPsSzgR042AKM4mOzUxKReqj2itK8\nvDwiIiLKXtvtdjIzM6/oQNOmTSv3KgNIuKI44nlFNOdZXuYTZjCf/2YUy3iSWZygrdmpiQSUjIwM\nMjIy6hWj2qLu7uHGV+piUX/hhRdQQfdO64C+fMMr/Jqt9CaZeaRzq9lpiQSMhIQEEhISyl6X1Mu6\nqbb9YrPZcDqdZa+dTid2u73OBxHfcZoQHuN/mchc3mYC83mADqUzdq1fF/F+1Rb1uLg4srOzycnJ\noaioiMWLF5OYmOh2rE6i+ZfPGEFPvuUgnfk3J5hIChZcXH7BUkFBvrmJikgFNV58tHLlSqZMmYLL\n5SI5OZmpU6eSkpICwKRJkzh48CDx8fGcPHmSJk2aEBoayvbt2wkJCbl0EK1+8enY12EhhRsxsPAw\nf2YbvSt8Rr/QRRrGlax+0RWlil2rsRZcTGQuL/Ec7zKeF/gdhURQcruBykJDwzh58pjb90SkdnRF\nqTQYgyb8hUn04t90YT876cE48t22ZNSWETGPirrUyX/oxAMs4B4+4mFgI3Ek8LnZaYlIKRV1uSLr\nGMQgYAa/YR7JpJJId3aanZZIwFNRl3r5G/fgYAdfcgv/4mZmM5n2HK40rnXrcLdLIrUsUsSzVNSl\n3opozqv8Cgc7uEATvqMTLwPtyxXukh575d67+u8inqWiLh5zlPY8yRv0xaAND7OTcH7PVMI5YnZq\nIgFDRV08Lhd4lLfoz2bac4RddOdFoG0Vyx9FxHNU1KXB7OMaJvEX4thIFyCbaKYzlc4cMDs1Eb+l\noi4NLoeuTABuIIsQTvEtPXmbZBxsNzs1Eb+joi6NZi/deILZRJNNDpGsZiifcAe3AOhWAyIeoaIu\nje4Y7XiJ/6Ere/mEO5kLMGAAfPABnD0LVLcE0qplkSLV0L1fFNvk2OFYyOcO4FHgeuA94M/Anjoe\nUzcWE3+je7+ID8rHwOATDH6MwUCygV+xHljJ7SSSSlOKzU5SxGdopq7YXhm7ORbuZj6P8BYROJnP\nON7nQbLpUWUczdTF32imLn7jHLCAB7iJtYxkBc05xxcMYS0wiT9rzbtIFTRTV2yfid2UYkZg5UHu\n5jY+41NGMJ9x/JNbOU9zzdTF72imLn7NRRArgHv5iK7s5XN+xFRe5gBX8Q7AypVQVGRyliLm0kxd\nsX0+to1c7iKCWYMGwXffwahRcM89hI9NIv/UcTcxgqCKk696YpN4Ez3OzkuKjGI3duxwLj5Wzw6M\nBe4FugNp3MdyEknjdk7QtobYJe+pjSPeQu0XCVCXbuubi8HrGAzEoBeQQQL3s4B9XE06w3iCWXQ1\nOVuRhqSZumIHROxWnGY46SSynJ/wDvnE8Bm38Rm38QVDOE1I2Wc0UxdvofaLDxUZxTYvtgULfdlU\nWtI/I54NbCSu9NWzbCouhqZN3cQTaVwq6j5aZBTb3NjBnGIIX5QW+TdwhIXBzTfDkCGQkAB9+6rI\niylU1P2kyCi2WbFLTrh2AoaU22xYaDvyxyVF/qab4PrroUULN58X8SwVdb8rMoptfuxwOpDPLcAt\nwCDAAfwbWFduc5b7hJZFiqeoqHtNIVBsf47dkkLiCGYgMxjIOgaxlvNY2UgcG4ljE8+z4uBB6NTJ\nzTEuad06vMqHbusXg4CKeh33K7Zieyq2QVf2cj2biGMj1/MKw9u2hZAQiIuD/v259+WZrD9zmn2V\n4lT146dVOKKiXsf9iq3YDRG75D3jwgXYswc2boTNm1n5yiv04SpaUchWerOFPmxlDttYxw4cnKRN\n5Rgq6gHvSop6UAPlIhKgSp7AZGni7rq+/bTnML3ZSh+2cDPwMI8Rw3ccI5ztxJbbgKNHoV27sk+r\nXSO1oZm6Yiu2ybEtXOAafigr5w6eoydFxFByh5pdFbaP2EV3dnNtuQumSmJpZu9/1H7xsh9WxVbs\n+sU26MBhurOrdEumO6Ppzi66sYcCQtlDN3azmT0UsQfYA+QA+2mKC5eb2JrV+xIVdZ/5YVVsxa5v\nDIPOHORadtONwXTjd6X/3kMkObRnPwe4hn1czQ9cU7bl8ktyOUUucPn9K1XsvY+Kul/8sCq2Ytd/\nvxULEXzP1ewrV9J/wM672IjFTi5BFJOLnVzs5PEv9lPEAeAAsL/06wEgSMXeNCrqXvIDpdiK7Qux\nQzmJjbzSsn4bV/F7ruIAXdjPVRwo3XI4B/wHOFS2WTiEwSHgcOl2pPTr+ZC2HK/iZK7UnYq6D/1A\nKbZi+0rsthyjI/+hE4dKt3voxHN04hAdOEx7jpR9bctRgtq1g/btS7bw8JIVPBe/Xvx3WBi3/HQM\nuYUF5AMngQulR1Qb6BIVdZ/4gfocSGig2L72PckAftRAsT21v7FiZ3Dp/wtfyru8cJqQTzjQAWhX\nuoXThHZcKP13yb6w0q0tkYSRTygFFBBKPgWs4gLRwInLtrPNW/L/Xv8TtG59aQsNvfQ1NLTknjwW\ni5vcfHNJaIOsU09LS2PKlCm4XC4mTJjAM888U2nME088wcqVK2nVqhXvvfce/fr1q1MSgSWDykU9\nUGWYnYAXycD3/7/I5wIGRyhpx1xS3S+MvQA0wUVrThJGOGeYRA+SaMMJ2nKcNpygDc/Q5twZUh55\nhNZQtoWW+3cIYAVOlW6nmzShe79+JVf2BgfzdkE+p3mIU4RwmuCyrZBWFBY8AkuWQKtWl7aWLS9t\nF183b17lLw13zPhFUm1Rd7lcTJ48mfT0dGw2G/Hx8SQmJuJwOMrGrFixgu+//57s7GwyMzN55JFH\nWL9+vccTFRH/dYGmHCesdEVOZw5W+gX3JLX5SyKI8wRzmhAiCb1wgtBNmwiGsi2Em8vKeQinCCOf\nVrxJK2DJ3XfTCmhVOrYl0MrShG6dO8GZM1BYCOfPU2gYnAXOAmdKv16+nSvdZgHnmMRZWnCO5pyj\nOUU04xzTKSrI52GLhaLSsUVAEU0o4kLpv69MtUU9KyuLqKgoIiMjAUhKSiI1NbVCUV++fDkPPvgg\nAAMGDOD48eMcOnSITjXczEhExNOKsXKCtpzgBO7bRsluPvWam7EA4WDkw4EDZXuaAC2AFhyhBWfL\nbb1oyRpacJbmnCv39b9oTp+y180ooiVnaMMZmvN4aZk/h5XzNKOIZiyhGbfTjCKac47BfFXn70G1\nRT0vL4+IiIiy13a7nczMzBrH5ObmVirqlgp/slT154u7/XUZ6wuxXyjdGiK2r35PFLvEC1Xs90Rs\nT+9vyNie+Bmpan9dY1R0ASgECmnv5t3BVXzq0Sr2z65if1qtcqlKtUXdUsve0eWN/Ms/p8uXRUQa\nh7u7DpWx2Ww4nZdu/+90OrHb7dWOyc3NxWazeThNERGpjWqLelxcHNnZ2eTk5FBUVMTixYtJTEys\nMCYxMZH58+cDsH79etq2bat+uoiISaptvwQFBTFnzhxGjBiBy+UiOTkZh8NBSkoKAJMmTWLkyJGs\nWLGCqKgogoODeffddxslcRERccNoYCtXrjR69OhhREVFGTNmzGjow3mtffv2GQkJCUZsbKzRs2dP\nY9asWWanZLri4mKjb9++xh133GF2KqbKz8837rrrLiMmJsZwOBzGunXrzE7JNNOnTzdiY2ONXr16\nGffdd59x9uxZs1NqNOPHjzc6duxo9OrVq2zf0aNHjeHDhxvR0dHGrbfeauTn59cYp9r2S31dXOee\nlpbG9u3bWbhwITt27GjIQ3otq9XKa6+9xrfffsv69et58803A/Z7cdGsWbOIjY2t9Ql5f/Xkk08y\ncuRIduzYwdatWyssGQ4kOTk5zJ07l82bN7Nt2zZcLheLFi0yO61GM378eNLSKq58mTFjBrfeeiu7\ndu1i2LBhzJgxo8Y4DVrUy69zt1qtZevcA1Hnzp3p27cvACEhITgcDvbv329yVubJzc1lxYoVTJgw\nIaBXR504cYI1a9bw85//HChpebZpc/mj7QJD69atsVqtFBYWUlxcTGFhYUAtuhg8eDBhYWEV9pW/\nDujBBx9k2bJlNcZp0KLubg17Xl5eQx7SJ+Tk5PD1118zYMAAs1MxzVNPPcUf/vAHmrh97Fvg2Lt3\nLx06dGA3ir1JAAACMklEQVT8+PH079+fiRMnUlhYaHZapggPD+fpp5/m6quvpkuXLrRt25bhw4eb\nnZapyl/I2alTJw4dOlTjZxr0JyrQ/6x259SpU4wdO5ZZs2YREhJS8wf80D/+8Q86duxIv379AnqW\nDlBcXMzmzZt59NFH2bx5M8HBwbX6E9sf7d69m9dff52cnBz279/PqVOn+PDDD81Oy2tYLJZa1dQG\nLeq1WeceSM6fP89dd93F/fffz6hRo8xOxzRr165l+fLldO3alfvuu4/Vq1czbtw4s9Myhd1ux263\nEx8fD8DYsWPZvHmzyVmZY+PGjQwaNIh27doRFBTEmDFjWLt2rdlpmapTp04cPHgQgAMHDtCxY8ca\nP9OgRb0269wDhWEYJCcnExsby5QpU8xOx1TTp0/H6XSyd+9eFi1axNChQ8uudQg0nTt3JiIigl27\ndgGQnp5Oz549Tc7KHDExMaxfv54zZ85gGAbp6enExsaanZapEhMTef/99wF4//33azcZbKjlORet\nWLHC6N69u3Httdca06dPb+jDea01a9YYFovF6NOnj9G3b1+jb9++xsqVK81Oy3QZGRnGnXfeaXYa\npvrmm2+MuLg4o3fv3sbo0aON48ePm52SaWbOnFm2pHHcuHFGUVGR2Sk1mqSkJOOqq64yrFarYbfb\njXfeecc4evSoMWzYsDotaWyUh2SIiEjjCOylByIifkZFXUTEj6ioi4j4ERV1ERE/oqIuIuJHVNRF\nRPzI/wE3VTs9UmMbFgAAAABJRU5ErkJggg==\n"
}
],
"prompt_number": 167
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment