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