Bernoulli.py
From Werner KRAUTH
(Difference between revisions)
Revision as of 05:29, 6 June 2024
import random
p = 0.7
count_a = 0
N_trial = 1000000
for iter in range(N_trial):
Upsilon = random.uniform(0.0, 1.0)
if Upsilon < p: count_a += 1
print(count_a / N_trial)
