Add exercise 1.22
This commit is contained in:
parent
2a01de329a
commit
673fa5f790
1 changed files with 25 additions and 0 deletions
25
chapter-1/ex-22.py
Executable file
25
chapter-1/ex-22.py
Executable file
|
@ -0,0 +1,25 @@
|
||||||
|
#!/usr/bin/env python
|
||||||
|
|
||||||
|
import numpy as np
|
||||||
|
import random
|
||||||
|
|
||||||
|
random.seed()
|
||||||
|
def experiment(p, n, reps):
|
||||||
|
summed = 0
|
||||||
|
|
||||||
|
for r in range(reps):
|
||||||
|
heads = 0
|
||||||
|
for i in range(n):
|
||||||
|
if random.random() < p:
|
||||||
|
heads += 1
|
||||||
|
summed += heads
|
||||||
|
summed /= reps
|
||||||
|
return summed
|
||||||
|
|
||||||
|
print("p: 0.3, n: 10, reps = 8, result X: ", experiment(0.3, 10, 8))
|
||||||
|
print("p: 0.3, n: 100, reps = 8, result X: ", experiment(0.3, 100, 8))
|
||||||
|
print("p: 0.3, n: 1000, reps = 8, result X: ", experiment(0.3, 1000, 8))
|
||||||
|
|
||||||
|
print("p: 0.3, n: 10, reps = 1, result X: ", experiment(0.3, 10, 1))
|
||||||
|
print("p: 0.3, n: 100, reps = 1, result X: ", experiment(0.3, 100, 1))
|
||||||
|
print("p: 0.3, n: 1000, reps = 1, result X: ", experiment(0.3, 1000, 1))
|
Loading…
Add table
Reference in a new issue