Computes the error function, which is typically used to compute areas under
the Gaussian probability distribution. See the manual for more information.
Also see the function `gaussian-cdf.'

This implementation follows Numerical Recipes in C, section 6.2