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