Little math problem(s) February 2019

[Albert Chan]

I tried best-of-21 simulation, the result match theoretical numbers well.

PHP Code:

from random import random

def simulation(n, best_of, total_win=0, random=random):
    for t in xrange(n):
        lst = [0,0]
        while max(lst) < best_of: lst[random() >= 0.6] += 1
        if lst[0] == best_of: total_win += 1
    return total_win / n 


for i in range(10): print simulation(10000, 21)
0.904
0.9069
0.9079
0.9043
0.9046
0.9042
0.8995
0.9003
0.9062
0.9037