fn std_normal_cdf(x: f64) -> f64
Standard normal CDF using the error function. pnorm(x) = 0.5 * erfc(-x / sqrt(2))