Break loop whem sum==0 to avoid zero division; ref #4345

This commit is contained in:
Yukihiro "Matz" Matsumoto
2019-03-28 14:02:30 +09:00
parent a26856f883
commit 3b92b642eb
+1
View File
@@ -120,6 +120,7 @@ erf(double x)
term *= xsqr/j;
sum += term/(2*j+1);
++j;
if (sum == 0) break;
} while (fabs(term/sum) > DBL_EPSILON);
return two_sqrtpi*sum;
}