pub fn gamma_cdf(a: f64, x: f64) -> f64
The regularized lower incomplete gamma function P(a, x): its series below a + 1, its continued fraction above (Numerical Recipes, 6.2).