Fix provided implementation of erf and erfc functions.

This commit is contained in:
Xuejie Xiao
2012-11-19 15:13:18 -05:00
parent 1b2257c8ff
commit f7dd27a928
2 changed files with 9 additions and 4 deletions
+3 -2
View File
@@ -51,7 +51,7 @@ erf(double x)
term *= xsqr/j;
sum += term/(2*j+1);
++j;
} while (fabs(term)/sum > MATH_TOLERANCE);
} while (fabs(term/sum) > MATH_TOLERANCE);
return two_sqrtpi*sum;
}
@@ -64,7 +64,8 @@ erfc(double x)
double b = x;
double c = x;
double d = x*x+0.5;
double q1, q2;
double q1;
double q2 = b/d;
double n = 1.0;
double t;
if (fabs(x) < 2.2) {