pub fn normal_cdf(x: f64) -> f64
Standard normal CDF (error function approximation).
Returns Φ(x) = P(Z ≤ x) for Z ~ N(0,1).