double erf(double);