monte_carlo_simulation_in_ptyhon

do_n_flips.py

def do_n_flips(n):
	return ''.join([str(random.getrandbits(1)) for i in range(n)])

monte_carlo_solve.py

def monte_carlo_solve(n, k, j):
	ret = 0.0
	for x in range(j):
		ret = ret + payoff(n, k)
	return ret / j

payoff.py

def payoff(n, k):
	if do_n_flips(n).find(''.rjust(k, '0')) > -1:
		return 1.0
	else:
		return 0.0