Merge pull request #563 from xxuejie/master

Fix manual implementations of erf and erfc functions

2 files changed, 11 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) {
diff --git a/test/t/math.rb b/test/t/math.rb
index d71e44fc9..780b805d2 100644
--- a/test/t/math.rb
+++ b/test/t/math.rb
@@ -113,5 +113,13 @@ if Object.const_defined?(:Math)
   assert('Math.erfc 1') do
     check_float(Math.erfc(1), 0.157299207050285)
   end
+
+  assert('Math.erf -1') do
+    check_float(Math.erf(-1), -0.8427007929497148)
+  end
+
+  assert('Math.erfc -1') do
+    check_float(Math.erfc(-1), 1.8427007929497148)
+  end
 end