From 6bbfa2fa6d69a405f8a37bea9bc09233103e3c2a Mon Sep 17 00:00:00 2001
From: Michel Juillard <michel.juillard@mjui.fr>
Date: Tue, 30 Sep 2014 08:13:31 +0200
Subject: [PATCH] removed 1e6 in expression forcing a non-singular Jacobian

---
 matlab/solve1.m | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/matlab/solve1.m b/matlab/solve1.m
index 17bb0a3ab..5fecdd356 100644
--- a/matlab/solve1.m
+++ b/matlab/solve1.m
@@ -101,7 +101,7 @@ for its = 1:maxit
     end
     if rcond_fjac < sqrt(eps)
         fjac2=fjac'*fjac;
-        p=-(fjac2+1e6*sqrt(nn*eps)*max(sum(abs(fjac2)))*eye(nn))\(fjac'*fvec);
+        p=-(fjac2+sqrt(nn*eps)*max(sum(abs(fjac2)))*eye(nn))\(fjac'*fvec);
     else
         p = -fjac\fvec ;
     end