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