quantile_gamma_poisson function quantile_gamma_poisson(p:Real, k:Real, θ:Real) -> Integer Quantile of a gamma-Poisson distribution. p: The cumulative probability. k: Shape. θ: Scale. Return: the quantile.