diff --git a/matlab/resol.m b/matlab/resol.m
index 2ad99cf2c..ec390ba0f 100644
--- a/matlab/resol.m
+++ b/matlab/resol.m
@@ -109,9 +109,11 @@ if info(1)
 end
 
 if options.loglinear
+    threshold = 1e-16;
     % Find variables with non positive steady state. Skip auxiliary
     % variables for lagges/leaded exogenous variables
-    idx = find(dr.ys(get_all_variables_but_lagged_leaded_exogenous(M)) < 1e-9);
+    idx = find(dr.ys(get_all_variables_but_lagged_leaded_exogenous(M))<threshold);
+    idx
     if length(idx)
         if options.debug
             variables_with_non_positive_steady_state = M.endo_names(idx,:);
@@ -132,7 +134,7 @@ if options.loglinear
             end
         end
         info(1)=26;
-        info(2)=sum(dr.ys(dr.ys<1e-9).^2);
+        info(2)=sum(dr.ys(dr.ys<threshold).^2);
         return
     end
 end