Error Function, Gaussian Probability Integral, etc.
These subprograms compute values of the error function, the complementary
error function, and an exponentially scaled complementary error function.
This software is part of the MATH77/mathc90 collection here.