pub fn std_normal_quantile(p: f64) -> f64
Φ⁻¹(p): a rational approximation (Abramowitz and Stegun 26.2.23, error below 4.5e-4), polished by Newton’s method on the exact CDF.