Computes the error function erf(x) = (2/\sqrt(\pi)) \int_0^x dt \exp(-t^2).
err contains an estimate of the absolute error in the value z.
This function is from the GNU Scientific Library, see http://www.gnu.org/software/gsl/ for documentation.
Package: gsl