pub fn normal_ppf(p: f64) -> f64
Inverse normal CDF (probit function) via rational approximation.
Reference: Peter Acklam’s approximation (max error ~1.15e-9).