# verval.py - radioactief verval nabootsen met dobbelstenen
import math
import random

random.seed(159)               # vaste startwaarde: elke keer dezelfde 'toevallige' worpen
N0 = 1000                      # aantal dobbelstenen ('kernen') bij de start
p = 1 / 6                      # kans dat een steen een zes gooit ('vervalt') per worp

over = N0
print("worp  over  model")
for worp in range(1, 16):
    vervallen = sum(1 for _ in range(over) if random.randint(1, 6) == 6)
    over = over - vervallen
    model = N0 * (1 - p) ** worp
    print(f"{worp:4d} {over:5d} {model:6.1f}")

T_half = math.log(0.5) / math.log(1 - p)
print(f"halveringstijd volgens het model: {T_half:.2f} worpen")