bug fix in linear algebra that caused an infinite loop

diff --git a/linearalgebra.cpp b/linearalgebra.cpp
index ab75b7e94bed1e90f44928a2e6e156220027f493..effb1ad44e04abd9c1bdd8523ca3b05f4253de94 100644 (file)
--- a/linearalgebra.cpp
+++ b/linearalgebra.cpp
@@ -222,7 +222,7 @@ double LinearAlgebra::betacf(const double a, const double b, const double x) {
         if (fabs(d) < FPMIN) d=FPMIN;
         d=1.0/d;
         h=d;
-        for (m1=1;m1<=MAXIT;m++) {
+        for (m1=1;m1<=MAXIT;m1++) {
             m2=2*m1;
             aa=m1*(b-m1)*x/((qam+m2)*(a+m2));
             d=1.0+aa*d;