binomial <- rbinom(100000, 1000,
                   1 / 1000)

poisson <- rpois(100000, 1)

compare_histograms(binomial,
                   poisson)

$$X \sim \text{Poisson}(\lambda)$$

$$E[X]=\lambda$$

$$\text{Var}(X)=\lambda$$