Bernoulli poisson.py

From Werner KRAUTH

import random, math
p = 0.7
lam = - math.log(1.0 - p)
N_trial = 1000000
count_a = 0
for iter in range(N_trial):
    t = -math.log(random.uniform(0., 1.0)) / lam
    if t < 1.0: count_a += 1
print(count_a / N_trial)