What function to use to model probabilistic bitcoin mining?

I would like to calculate profitability of various mining strategies. For that, i need to simulate many random mining intervals, that reflect probabilistic distribution of real block mining. What would be such formula that could use random number generator?

Block arrival times follow the Poisson Distribution with λ = 10 minutes. See this SO question on how to do it in pure python