pub fn beta_cdf(a: f64, b: f64, x: f64) -> f64
The regularized incomplete beta function I_x(a, b), by its continued fraction (Numerical Recipes, section 6.4).