pub fn incomplete_beta(a: f64, b: f64, x: f64) -> f64
Regularised incomplete beta I_x(a, b) by its continued fraction (Lentz’s method), using the symmetry I_x(a, b) = 1 − I_{1−x}(b, a) where the fraction converges faster.