erf(x) and the complimentary error function erfc(x) using
the integration routines from GSL.
Hint: you can implement mutually recursive functions:
double erf(double);
double erfc(double);
double erf(double x){
if(x < 0) return -erf(-x);
if(x > 3) return 1-erfc(x); /* more effective */
/* calculate erf(x) */
return result;
}
double erfc(double x){
if(x < 0) return 2-erfc(-x);
if(x < 3) return 1-erf(x); /* more effective */
/* calculate erfc(x) */
return result;
}