pub fn erfcinv<F>(q: F) -> Fwhere F: Float,
The x with erfc(x) = q. 0 gives ∞, 2 gives -∞, and q outside 0 to 2 gives NaN.