The erf functions return the error function of x, where erf(x),
    erff(x), and erfl(x) equal 2/sqrt(pi) times the area under the
    curve e**(-t**2) between 0 and x.

    The erfc functions return (1.0 - erf(x)). The erfc function can
    result in an underflow as x gets large.