pub fn erf_e(x: f64) -> Result<Result, Value>
This routine computes the error function erf(x), where erf(x) = (2/\sqrt(\pi)) \int_0^x dt \exp(-t^2).