erf
error function and related functions in pure standard Fortran 2008
...The calculation is based on calculating w(z), the Faddeyeva (plasma dispersion) function. Two algorithms are implemented: TOMS 680 and 916, http://www.netlib.org/toms/. The simple test program calculates values in the Abramowitz, Stegun.