pub fn error_f(x: f64) -> f64
Returns the error function erf(x) = (2/√π) ∫₀ˣ e⁻ᵗ² dt.
use cdflib::special::error_f; let y = error_f(0.8); assert!((y - 0.74210096).abs() < 1e-8);