pub fn regularized_incomplete_beta(a: f64, b: f64, x: f64) -> f64
Regularized incomplete beta function I_x(a, b).
Uses a continued fraction expansion (Lentz’s method).