From 3114ad297d40ecd47bd3c1fe7997bd76a81e8c75 Mon Sep 17 00:00:00 2001 From: Sarah Westcott Date: Tue, 13 Mar 2012 15:22:00 -0400 Subject: [PATCH] bug fix in linear algebra that caused an infinite loop --- linearalgebra.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/linearalgebra.cpp b/linearalgebra.cpp index ab75b7e..effb1ad 100644 --- 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; -- 2.39.2