Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- from random import choice
- import matplotlib.pyplot as plt
- %matplotlib inline
- from random import choice
- urn1 = ["W" for _ in range(8)]
- urn1 += ["B" for _ in range(3)]
- urn2 = ["W" for _ in range(7)]
- urn2 += ["B" for _ in range(8)]
- def experiment(N):
- cnt = 0
- M = 0
- for _ in range(N):
- balls = []
- urn1_copied = urn1.copy()
- urn2_copied = urn2.copy()
- for _ in range(2):
- b = choice(urn2_copied)
- urn2_copied.remove(b)
- urn1_copied.append(b)
- balls.append(b)
- ch = choice(urn1_copied)
- if ch == "W":
- M += 1
- if balls[0] != balls[1]:
- cnt += 1
- return cnt / M
- x = [i for i in range(1000,100000,1000)]
- y = [experiment(i) for i in x]
- plt.plot(x[3:],y[3:])
- plt.plot([1,100000],[0.5373,0.5373],color = "red")
- print(experiment(100000))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement