from random import random from math import ceil, log def exponential(p): y = random() x = ceil( -log(1-y) / p ) return(x) p = 1/6 for i in range(0,5): y = exponential(p) print(i, y)