Advertisement
Guest User

Untitled

a guest
Oct 23rd, 2019
90
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.89 KB | None | 0 0
  1. from random import choice
  2. import matplotlib.pyplot as plt
  3. %matplotlib inline
  4. from random import choice
  5. urn1 = ["W" for _ in range(8)]
  6. urn1 += ["B" for _ in range(3)]
  7. urn2 = ["W" for _ in range(7)]
  8. urn2 += ["B" for _ in range(8)]
  9. def experiment(N):
  10. cnt = 0
  11. M = 0
  12. for _ in range(N):
  13. balls = []
  14. urn1_copied = urn1.copy()
  15. urn2_copied = urn2.copy()
  16. for _ in range(2):
  17. b = choice(urn2_copied)
  18. urn2_copied.remove(b)
  19. urn1_copied.append(b)
  20. balls.append(b)
  21. ch = choice(urn1_copied)
  22. if ch == "W":
  23. M += 1
  24. if balls[0] != balls[1]:
  25. cnt += 1
  26. return cnt / M
  27. x = [i for i in range(1000,100000,1000)]
  28. y = [experiment(i) for i in x]
  29. plt.plot(x[3:],y[3:])
  30. plt.plot([1,100000],[0.5373,0.5373],color = "red")
  31. print(experiment(100000))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement