fn regularized_incomplete_beta(x: f64, a: f64, b: f64) -> f64
Regularized incomplete beta function via continued fraction (Lentz’s method).
I_x(a, b) = B_x(a, b) / B(a, b)