pub fn erf(x: f64) -> f64
Formula and coefficients from https://en.wikipedia.org/wiki/Error_function#Numerical_approximations (Abramowitz and Stegun), Max error 3e-7.
3e-7