Updated Lecture 2 and Lecture 3

Week2/Lecture2.ipynb

+%% Cell type:markdown id:4e98befd-a710-4193-9393-d52d8a7fbcc3 tags:
+
+## Poisson distribution
+$$ P(k; \lambda) = \frac{\lambda^k}{k!} \mbox{e}^{-\lambda}$$
+
+%% Cell type:code id:9c2b256c-333a-4b6c-8241-523b8aeb919f tags:
+
+``` python
+from scipy.stats import poisson, binom
+```
+
+%% Cell type:code id:ccf39158-28be-4ba4-a898-a5b14be924e1 tags:
+
+``` python
+poisson.pmf(2,3.5)
+```
+
+%%%% Output: execute_result
+
+    0.18495897346170082
+
+%% Cell type:code id:6643ab02-3c24-40b8-90e2-c0ffab308f16 tags:
+
+``` python
+poisson.pmf(0,3.5)+poisson.pmf(1,3.5)+poisson.pmf(2,3.5)
+```
+
+%%%% Output: execute_result
+
+    0.3208471988621341
+
+%% Cell type:code id:19fc9d45-64f9-4e0f-adab-b044c59a083e tags:
+
+``` python
+poisson.cdf(2,3.5)
+```
+
+%%%% Output: execute_result
+
+    0.32084719886213414
+
+%% Cell type:markdown id:9d2578c1-302c-4ea4-8308-8c80d5ab9288 tags:
+
+## Binomial distribution
+$$ B(k;N,p) = \left( \begin{array}{c} N \\ k \\ \end{array} \right) p^k (1-p)^{N-k} $$
+
+%% Cell type:code id:a9d78b0f-49ca-44c8-964b-ab1d593ce540 tags:
+
+``` python
+binom.pmf(2,10,0.3)
+```
+
+%%%% Output: execute_result
+
+    0.2334744405
+
+%% Cell type:code id:4c029640-fec5-4aeb-9644-7e965b0dcbee tags:
+
+``` python
+binom.pmf(0,10,0.3)+binom.pmf(1,10,0.3)+binom.pmf(2,10,0.3)
+```
+
+%%%% Output: execute_result
+
+    0.3827827864000002
+
+%% Cell type:code id:fa781bea-3265-4915-b299-0870a76b8053 tags:
+
+``` python
+binom.cdf(2,10,0.3)
+```
+
+%%%% Output: execute_result
+
+    0.3827827864
+
+%% Cell type:code id:a88c9312-41d1-44ba-b600-ca8fcbbb5c75 tags:
+
+``` python
+```

Week3/Lecture3.ipynb

+%% Cell type:code id: tags:
+
+``` python
+import numpy as np
+```
+
+%% Cell type:code id: tags:
+
+``` python
+import matplotlib.pyplot as plt
+```
+
+%% Cell type:code id: tags:
+
+``` python
+T0=1
+sigT=0.2
+```
+
+%% Cell type:markdown id: tags:
+
+### error propagation
+
+%% Cell type:code id: tags:
+
+``` python
+f0=1/T0
+sigf=1/T0**2 * sigT
+print("error propagation: freq ",f0," +/- ",sigf)
+```
+
+%%%% Output: stream
+
+    error propagation: freq  1.0  +/-  0.2
+
+%% Cell type:markdown id: tags:
+
+### simulation
+
+%% Cell type:code id: tags:
+
+``` python
+T=np.random.normal(T0,sigT,1000)
+```
+
+%% Cell type:code id: tags:
+
+``` python
+plt.hist(T)
+```
+
+%%%% Output: execute_result
+
+    (array([ 13.,  52., 116., 158., 215., 197., 139.,  80.,  19.,  11.]),
+     array([0.48111168, 0.59232677, 0.70354186, 0.81475695, 0.92597204,
+            1.03718713, 1.14840222, 1.25961732, 1.37083241, 1.4820475 ,
+            1.59326259]),
+     <BarContainer object of 10 artists>)
+
+%%%% Output: display_data
+
+    ![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAANGElEQVR4nO3db4xl9V3H8ffHpTRqGwvuQNYFHGxWLSSCdUVi1WwlkX8Ptk2KAU1LGpLVSE1NfMDSB2JiSLYPrMYoNtgSMLEgsbRgwCpBK5pK6WIosCB2hRVWCLuUxtaaYHb5+mAOOO7OMHfm3jt359v3K5nMveeeO+f7Y8N7z5ydeydVhSSpl++a9QCSpMkz7pLUkHGXpIaMuyQ1ZNwlqaGTZj0AwObNm2t+fn7WY0jShvLII4+8XFVzSz12QsR9fn6evXv3znoMSdpQkvz7co95WUaSGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaOiFeoSqtZH73vTM79oE9l8/s2NJaeeYuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWpoxbgnOTPJ3yV5Ksm+JB8dtp+a5P4kXxs+n7LoOdcn2Z/k6SQXT3MBkqTjjXLmfgT4zap6F3AhcG2Sc4DdwANVtQ14YLjP8NiVwLnAJcBNSTZNY3hJ0tJWjHtVvVhV/zzc/hbwFLAV2AncNux2G/C+4fZO4I6qerWqngX2AxdMeG5J0ptY1TX3JPPAjwNfBk6vqhdh4S8A4LRht63A84uednDYduzX2pVkb5K9hw8fXsPokqTljBz3JG8DPgv8RlV98812XWJbHbeh6uaq2l5V2+fm5kYdQ5I0gpHinuQtLIT9z6rqrmHzS0m2DI9vAQ4N2w8CZy56+hnAC5MZV5I0ihV/WUeSAJ8GnqqqTyx66B7gamDP8PnuRds/k+QTwA8A24CHJzm0tJ5m9YtC/CUhGscov4npPcAHgceTPDps+xgLUb8zyTXAc8AVAFW1L8mdwJMs/KTNtVV1dNKDS5KWt2Lcq+ofWfo6OsBFyzznRuDGMeaSJI3BV6hKUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLU0EmzHkAby/zue2c9gqQReOYuSQ0Zd0lqyLhLUkMrxj3JLUkOJXli0bbfTvIfSR4dPi5b9Nj1SfYneTrJxdMaXJK0vFHO3G8FLlli++9V1fnDx30ASc4BrgTOHZ5zU5JNkxpWkjSaFeNeVQ8Cr4z49XYCd1TVq1X1LLAfuGCM+SRJazDONfePJHlsuGxzyrBtK/D8on0ODtuOk2RXkr1J9h4+fHiMMSRJx1pr3P8YeCdwPvAi8LvD9iyxby31Barq5qraXlXb5+bm1jiGJGkpa4p7Vb1UVUer6jXgT/i/Sy8HgTMX7XoG8MJ4I0qSVmtNcU+yZdHd9wOv/yTNPcCVSd6a5GxgG/DweCNKklZrxbcfSHI7sAPYnOQgcAOwI8n5LFxyOQD8CkBV7UtyJ/AkcAS4tqqOTmVySdKyVox7VV21xOZPv8n+NwI3jjOUJGk8vkJVkhoy7pLUkHGXpIaMuyQ1ZNwlqSF/E5N0gprlb706sOfymR1bk+GZuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQyfNegCt3vzue2c9gqQTnGfuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNbRi3JPckuRQkicWbTs1yf1JvjZ8PmXRY9cn2Z/k6SQXT2twSdLyRjlzvxW45Jhtu4EHqmob8MBwnyTnAFcC5w7PuSnJpolNK0kayYpxr6oHgVeO2bwTuG24fRvwvkXb76iqV6vqWWA/cMFkRpUkjWqt19xPr6oXAYbPpw3btwLPL9rv4LDtOEl2JdmbZO/hw4fXOIYkaSmT/gfVLLGtltqxqm6uqu1VtX1ubm7CY0jSd7a1xv2lJFsAhs+Hhu0HgTMX7XcG8MLax5MkrcVa434PcPVw+2rg7kXbr0zy1iRnA9uAh8cbUZK0Wiv+mr0ktwM7gM1JDgI3AHuAO5NcAzwHXAFQVfuS3Ak8CRwBrq2qo1OaXZK0jBXjXlVXLfPQRcvsfyNw4zhDSZLG4ytUJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktTQim/5K+k7z/zue2dy3AN7Lp/JcTvyzF2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDJ43z5CQHgG8BR4EjVbU9yanAnwPzwAHgF6vqG+ONKUlajUmcub+3qs6vqu3D/d3AA1W1DXhguC9JWkdjnbkvYyewY7h9G/BF4LopHGfm5nffO+sRJGlJ4565F/A3SR5JsmvYdnpVvQgwfD5tzGNIklZp3DP391TVC0lOA+5P8i+jPnH4y2AXwFlnnTXmGJKkxcY6c6+qF4bPh4DPARcALyXZAjB8PrTMc2+uqu1VtX1ubm6cMSRJx1hz3JN8b5K3v34b+AXgCeAe4Opht6uBu8cdUpK0OuNcljkd+FyS17/OZ6rqC0m+AtyZ5BrgOeCK8ceUJK3GmuNeVc8A5y2x/evAReMMJUkaj69QlaSGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1NBJsx5Akl43v/vemRz3wJ7LZ3LcafLMXZIaMu6S1JBxl6SGjLskNWTcJamhFj8tM6t/YZekE5Vn7pLUUIszd0kaxyy/+5/Wz9h75i5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1NLW4J7kkydNJ9ifZPa3jSJKON5W4J9kE/BFwKXAOcFWSc6ZxLEnS8aZ15n4BsL+qnqmq/wHuAHZO6ViSpGNM6+0HtgLPL7p/EPipxTsk2QXsGu7+V5KnpzTLetgMvDzrIabEtW08XdcFDdeWj79xcy1r+8HlHphW3LPEtvp/d6puBm6e0vHXVZK9VbV91nNMg2vbeLquC1zbakzrssxB4MxF988AXpjSsSRJx5hW3L8CbEtydpKTgSuBe6Z0LEnSMaZyWaaqjiT5CPDXwCbglqraN41jnSBaXF5ahmvbeLquC1zbyFJVK+8lSdpQfIWqJDVk3CWpIeO+CqO8pUKSHUkeTbIvyd+v94xrsdK6knxfkr9M8tVhXR+exZxrkeSWJIeSPLHM40nyB8PaH0vy7vWecS1GWNcvD+t5LMmXkpy33jOu1UprW7TfTyY5muQD6zXbuEZZ28QaUlV+jPDBwj8M/xvwQ8DJwFeBc47Z5x3Ak8BZw/3TZj33hNb1MeDjw+054BXg5FnPPuL6fg54N/DEMo9fBvwVC6/NuBD48qxnntC6fho4Zbh96UZZ1yhrG/bZBPwtcB/wgVnPPME/t4k1xDP30Y3ylgq/BNxVVc8BVNWhdZ5xLUZZVwFvTxLgbSzE/cj6jrk2VfUgC/MuZyfwp7XgIeAdSbasz3Rrt9K6qupLVfWN4e5DLLzWZEMY4c8M4NeBzwIb4f+xN4ywtok1xLiPbqm3VNh6zD4/DJyS5ItJHknyoXWbbu1GWdcfAu9i4YVojwMfrarX1me8qRtl/RvdNSx8d9JCkq3A+4FPznqWKZhYQ6b19gMdrfiWCiz89/wJ4CLgu4F/SvJQVf3rtIcbwyjruhh4FPh54J3A/Un+oaq+OeXZ1sMo69+wkryXhbj/zKxnmaDfB66rqqML30y2MrGGGPfRjfKWCgeBl6vq28C3kzwInAecyHEfZV0fBvbUwkXA/UmeBX4UeHh9Rpyqtm+VkeTHgE8Bl1bV12c9zwRtB+4Ywr4ZuCzJkar6/EynmoyJNcTLMqMb5S0V7gZ+NslJSb6HhXfCfGqd51ytUdb1HAtnEiQ5HfgR4Jl1nXJ67gE+NPzUzIXAf1bVi7MealxJzgLuAj54gn/nuGpVdXZVzVfVPPAXwK81CTtMsCGeuY+olnlLhSS/Ojz+yap6KskXgMeA14BPVdWb/jjXrI2yLuB3gFuTPM7CZYzrqmpDvO1qktuBHcDmJAeBG4C3wBtru4+Fn5jZD/w3C9+lnPBGWNdvAd8P3DSc4R6pDfJuiiOsbcNaaW2TbIhvPyBJDXlZRpIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrofwFjeCFv4/gf5wAAAABJRU5ErkJggg==)
+
+%% Cell type:code id: tags:
+
+``` python
+f=1/T
+plt.hist(f)
+```
+
+%%%% Output: execute_result
+
+    (array([ 78., 259., 293., 174., 103.,  59.,  15.,  12.,   5.,   2.]),
+     array([0.62764293, 0.77273059, 0.91781824, 1.0629059 , 1.20799356,
+            1.35308122, 1.49816887, 1.64325653, 1.78834419, 1.93343184,
+            2.0785195 ]),
+     <BarContainer object of 10 artists>)
+
+%%%% Output: display_data
+
+    ![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAQQklEQVR4nO3df6zddX3H8edLQOYGmTAupJa6y0zdhGVWd9eRsS0o22DyRzGRpW7BxpDUZbhg4h8W/pguS5OaTF0Wh6YKsSYqayaObjA3ZDpmVPBCKlAq2kkH1zb0+mMTt4Sl5b0/7hc5tuf2nnvPvfecfnw+kpvzPZ/v93u+r9v2++r3fs/3fG+qCklSW1406gCSpOVnuUtSgyx3SWqQ5S5JDbLcJalBp486AMB5551Xk5OTo44hSaeUBx988DtVNdFv3liU++TkJNPT06OOIUmnlCT/Od88T8tIUoMWLPckP5XkgSRfS7IvyZ934+cmuSfJN7vHc3rWuSnJgSSPJ7lyJb8BSdKJBjlyfxZ4fVW9GtgAXJXkUmAbcG9VrQfu7Z6T5GJgM3AJcBVwS5LTViC7JGkeC5Z7zflh9/SM7quATcCubnwXcE03vQm4vaqeraongAPAxuUMLUk6uYHOuSc5Lcle4AhwT1XdD1xQVYcBusfzu8XXAk/1rD7TjR3/mluTTCeZnp2dHeJbkCQdb6Byr6pjVbUBuBDYmOSXT7J4+r1En9fcWVVTVTU1MdH3Sh5J0hIt6mqZqvov4AvMnUt/OskagO7xSLfYDLCuZ7ULgUPDBpUkDW6Qq2Umkry0m34J8DvA14E9wJZusS3And30HmBzkjOTXASsBx5Y5tySpJMY5ENMa4Bd3RUvLwJ2V9U/JvkysDvJ9cCTwLUAVbUvyW7gMeAocENVHVuZ+JKkfjIOv6xjamqq/ITq4Ca33TWybR/ccfXIti3pxyV5sKqm+s3zE6qS1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBlnuktQgy12SGmS5S1KDLHdJapDlLkkNstwlqUGWuyQ1yHKXpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNej0UQc4lU1uu2vUESSpL4/cJalBlrskNWjBck+yLsnnk+xPsi/Jjd34e5J8O8ne7usNPevclORAkseTXLmS34Ak6USDnHM/Cryzqh5KcjbwYJJ7unkfqKq/7F04ycXAZuAS4GXA55K8sqqOLWdwSdL8Fjxyr6rDVfVQN/0MsB9Ye5JVNgG3V9WzVfUEcADYuBxhJUmDWdQ59ySTwGuA+7uhtyd5OMltSc7pxtYCT/WsNkOf/wySbE0ynWR6dnZ28cklSfMauNyTnAV8GnhHVf0A+BDwCmADcBh43/OL9lm9Thio2llVU1U1NTExsdjckqSTGKjck5zBXLF/oqruAKiqp6vqWFU9B3yEF069zADrela/EDi0fJElSQsZ5GqZALcC+6vq/T3ja3oWeyPwaDe9B9ic5MwkFwHrgQeWL7IkaSGDXC1zGXAd8EiSvd3YzcCbk2xg7pTLQeBtAFW1L8lu4DHmrrS5wStlJGl1LVjuVfVF+p9Hv/sk62wHtg+RS5I0BD+hKkkNstwlqUGWuyQ1yHKXpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBlnuktQgy12SGmS5S1KDLHdJapDlLkkNstwlqUELlnuSdUk+n2R/kn1JbuzGz01yT5Jvdo/n9KxzU5IDSR5PcuVKfgOSpBMNcuR+FHhnVb0KuBS4IcnFwDbg3qpaD9zbPaebtxm4BLgKuCXJaSsRXpLU34LlXlWHq+qhbvoZYD+wFtgE7OoW2wVc001vAm6vqmer6gngALBxmXNLkk5iUefck0wCrwHuBy6oqsMw9x8AcH632FrgqZ7VZrqx419ra5LpJNOzs7NLiC5Jms/A5Z7kLODTwDuq6gcnW7TPWJ0wULWzqqaqampiYmLQGJKkAQxU7knOYK7YP1FVd3TDTydZ081fAxzpxmeAdT2rXwgcWp64kqRBDHK1TIBbgf1V9f6eWXuALd30FuDOnvHNSc5MchGwHnhg+SJLkhZy+gDLXAZcBzySZG83djOwA9id5HrgSeBagKral2Q38BhzV9rcUFXHlju4JGl+C5Z7VX2R/ufRAa6YZ53twPYhckmShuAnVCWpQZa7JDXIcpekBlnuktQgy12SGmS5S1KDLHdJapDlLkkNstwlqUGWuyQ1yHKXpAZZ7pLUIMtdkho0yC1/pR+Z3HbXSLZ7cMfVI9mudKryyF2SGmS5S1KDLHdJapDlLkkNstwlqUGWuyQ1yHKXpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBC5Z7ktuSHEnyaM/Ye5J8O8ne7usNPfNuSnIgyeNJrlyp4JKk+Q1y5P4x4Ko+4x+oqg3d190ASS4GNgOXdOvckuS05QorSRrMguVeVfcB3xvw9TYBt1fVs1X1BHAA2DhEPknSEgxzzv3tSR7uTtuc042tBZ7qWWamGztBkq1JppNMz87ODhFDknS8pZb7h4BXABuAw8D7uvH0Wbb6vUBV7ayqqaqampiYWGIMSVI/Syr3qnq6qo5V1XPAR3jh1MsMsK5n0QuBQ8NFlCQt1pLKPcmanqdvBJ6/kmYPsDnJmUkuAtYDDwwXUZK0WAv+guwknwIuB85LMgO8G7g8yQbmTrkcBN4GUFX7kuwGHgOOAjdU1bEVSS5JmteC5V5Vb+4zfOtJlt8ObB8mlCRpOH5CVZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBlnuktQgy12SGmS5S1KDLHdJapDlLkkNstwlqUGWuyQ1yHKXpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMWLPcktyU5kuTRnrFzk9yT5Jvd4zk9825KciDJ40muXKngkqT5nT7AMh8DPgh8vGdsG3BvVe1Isq17/q4kFwObgUuAlwGfS/LKqjq2vLH1k2Zy210j2/bBHVePbNvSUi145F5V9wHfO254E7Crm94FXNMzfntVPVtVTwAHgI3LE1WSNKilnnO/oKoOA3SP53fja4Gnepab6cZOkGRrkukk07Ozs0uMIUnqZ7nfUE2fseq3YFXtrKqpqpqamJhY5hiS9JNtqeX+dJI1AN3jkW58BljXs9yFwKGlx5MkLcVSy30PsKWb3gLc2TO+OcmZSS4C1gMPDBdRkrRYC14tk+RTwOXAeUlmgHcDO4DdSa4HngSuBaiqfUl2A48BR4EbvFJGklbfguVeVW+eZ9YV8yy/Hdg+TChJ0nD8hKokNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBg3ym5jG3ih/S48kjSOP3CWpQZa7JDXIcpekBlnuktQgy12SGmS5S1KDLHdJapDlLkkNstwlqUGWuyQ1yHKXpAZZ7pLUoCZuHCatpFHdmO7gjqtHsl21wSN3SWqQ5S5JDRrqtEySg8AzwDHgaFVNJTkX+FtgEjgI/EFVfX+4mJKkxViOI/fXVdWGqprqnm8D7q2q9cC93XNJ0ipaidMym4Bd3fQu4JoV2IYk6SSGLfcC/iXJg0m2dmMXVNVhgO7x/H4rJtmaZDrJ9Ozs7JAxJEm9hr0U8rKqOpTkfOCeJF8fdMWq2gnsBJiamqohc0iSegx15F5Vh7rHI8BngI3A00nWAHSPR4YNKUlanCWXe5KfSXL289PA7wGPAnuALd1iW4A7hw0pSVqcYU7LXAB8Jsnzr/PJqvpskq8Cu5NcDzwJXDt8TEnSYiy53KvqW8Cr+4x/F7himFCSpOH4CVVJapDlLkkNstwlqUGWuyQ1yHKXpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQcP8gmxJK2hy210j2/bBHVePbNtaHh65S1KDPHKXdIJR/dTgTwzLxyN3SWqQ5S5JDbLcJalBlrskNcg3VCWNDd/IXT4rduSe5Kokjyc5kGTbSm1HknSiFTlyT3Ia8DfA7wIzwFeT7Kmqx1Zie5I0jBY/MLZSR+4bgQNV9a2q+j/gdmDTCm1LknSclTrnvhZ4quf5DPDrvQsk2Qps7Z7+MMnjA7zuecB3liXhyjLn8jLn8joVcp4KGWEZcua9Q23/5+ebsVLlnj5j9WNPqnYCOxf1osl0VU0NE2w1mHN5mXN5nQo5T4WMMN45V+q0zAywruf5hcChFdqWJOk4K1XuXwXWJ7koyYuBzcCeFdqWJOk4K3JapqqOJnk78M/AacBtVbVvGV56UadxRsicy8ucy+tUyHkqZIQxzpmqWngpSdIpxdsPSFKDLHdJatDYlfsgty1IcnmSvUn2Jfm31c7YZThpziQ/m+Qfknyty/nWEeW8LcmRJI/OMz9J/rr7Ph5O8trVztjlWCjnH3X5Hk7ypSSvXu2MXY6T5uxZ7teSHEvyptXKdtz2F8w5JvvRQn/vI9+PkqxL8vkk+7sMN/ZZZiz2ox9TVWPzxdybr/8B/ALwYuBrwMXHLfNS4DHg5d3z88c0583Ae7vpCeB7wItHkPW3gdcCj84z/w3APzH32YRLgftH9He/UM7fAM7ppn9/XHP2/Pv4V+Bu4E3jmHMc9qMBc458PwLWAK/tps8GvtFnfx+L/aj3a9yO3Ae5bcEfAndU1ZMAVXVklTPCYDkLODtJgLOY+0d5dHVjQlXd1217PpuAj9ecrwAvTbJmddK9YKGcVfWlqvp+9/QrzH12YtUN8OcJ8KfAp4FR/NsEBso5DvvRIDlHvh9V1eGqeqibfgbYz9yn8HuNxX7Ua9zKvd9tC47/Q3wlcE6SLyR5MMlbVi3dCwbJ+UHgVcx9eOsR4Maqem514i3KIN/LuLmeuaOksZNkLfBG4MOjzrKAcdiPBjFW+1GSSeA1wP3HzRq7/Wjc7ue+4G0LmMv8q8AVwEuALyf5SlV9Y6XD9Rgk55XAXuD1wCuAe5L8e1X9YIWzLdYg38vYSPI65sr9N0edZR5/Bbyrqo7NHWyOrXHYjwYxNvtRkrOY+4nsHX22P3b70bgduQ9y24IZ4LNV9T9V9R3gPmC131wbJOdbmfuxt6rqAPAE8EurlG8xTplbRST5FeCjwKaq+u6o88xjCrg9yUHgTcAtSa4ZaaL+xmE/GsRY7EdJzmCu2D9RVXf0WWTs9qNxK/dBbltwJ/BbSU5P8tPM3W1y/xjmfJK5oyKSXAD8IvCtVU05mD3AW7p3+y8F/ruqDo861PGSvBy4A7huDI8uf6SqLqqqyaqaBP4O+JOq+vvRpuprHPajQYx8P+rO998K7K+q98+z2NjtR2N1WqbmuW1Bkj/u5n+4qvYn+SzwMPAc8NGqOullaaPICfwF8LEkjzD3I9u7uiOkVZXkU8DlwHlJZoB3A2f05LybuXf6DwD/y9yR0qobIOefAT/H3JEwwNEawd34Bsg5FhbKOQ770SA5GY/96DLgOuCRJHu7sZuBl/fkHIv9qJe3H5CkBo3baRlJ0jKw3CWpQZa7JDXIcpekBlnuktQgy12SGmS5S1KD/h/fH36MP5XHewAAAABJRU5ErkJggg==)
+
+%% Cell type:code id: tags:
+
+``` python
+print("simulation: freq ",np.mean(f)," +/- ",np.std(f))
+```
+
+%%%% Output: stream
+
+    simulation: freq  1.0321737435939056  +/-  0.2231584633106588
+
+%% Cell type:code id: tags:
+
+``` python
+```