Function normal_cdf 

fn normal_cdf(z: f64) -> f64

Standard normal CDF via Abramowitz–Stegun 7.1.26 on erf, accurate to ~1.5e-7 — orders of magnitude finer than any decision read from a p.