pub fn erf<F>(x: F) -> Fwhere F: Float,
The error function, 2/√π times the integral of e^-t² from 0 to x.