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