simulate_gamma_poisson function simulate_gamma_poisson(k:Real, θ:Real) -> Integer Simulate a gamma-Poisson distribution. k: Shape. θ: Scale.