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