Purpose
To see the futility of law of large numbers

> n <- 10^6
> p <- 10^-6
> lambda <- n * p
> p.loss <- 1 - sum(dpois(1:10^6, lambda))
> p.loss
[1] 0.3678794
> p.loss.rec <- 1 - sum(dpois(1:5e+05, lambda))
> p.loss.rec
[1] 0.3678794
> p.win <- sum(dpois(1:round(2/3 * 10^6), lambda))
> 1/(2 * exp(1))
[1] 0.1839397