Fix provided implementation of erf and erfc functions.

1 files changed, 3 insertions, 2 deletions

diff --git a/src/math.c b/src/math.c
index 9aae87acd..3e4550946 100644
--- a/src/math.c
+++ b/src/math.c
@@ -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) {