pub fn erf_Z(x: f64) -> f64
This routine computes the Gaussian probability density function Z(x) = (1/\sqrt{2\pi}) \exp(-x^2/2).