pub fn norm_cdf(x: f64) -> f64
Standard normal CDF (Abramowitz & Stegun 7.1.26 via erf, |error| < 1.5e-7).