diff --git a/matlab/dyn_ramsey_static.m b/matlab/dyn_ramsey_static.m
index 7427603e9..c2b619b3f 100644
--- a/matlab/dyn_ramsey_static.m
+++ b/matlab/dyn_ramsey_static.m
@@ -53,8 +53,12 @@ elseif options_.steadystate_flag
         %solve for instrument, using univariate solver, starting at initial value for instrument
         inst_val = csolve(nl_func,ys_init(k_inst),'',options_.solve_tolf,100); 
     else
-        %solve for instrument, using multivariate solver, starting at initial value for instrument
-        [inst_val,info1] = dynare_solve(nl_func,ys_init(k_inst),options_); 
+        %solve for instrument, using multivariate solver, starting at
+        %initial value for instrument
+        opt = options_;
+        opt.jacobian_flag = 0;
+        [inst_val,info1] = dynare_solve(nl_func,ys_init(k_inst), ...
+                                        opt);
     end
     ys_init(k_inst) = inst_val;
     exo_ss = [oo.exo_steady_state oo.exo_det_steady_state];
@@ -63,7 +67,9 @@ elseif options_.steadystate_flag
 else
     n_var = M.orig_endo_nbr;
     xx = oo.steady_state(1:n_var);
-    [xx,check] = dynare_solve(nl_func,xx,options_);
+    opt = options_;
+    opt.jacobian_flag = 0;
+    [xx,check] = dynare_solve(nl_func,xx,opt);
     [junk,junk,steady_state] = nl_func(xx);
 end
 
diff --git a/matlab/partial_information/PCL_resol.m b/matlab/partial_information/PCL_resol.m
index 154bfb883..ba0b1666f 100644
--- a/matlab/partial_information/PCL_resol.m
+++ b/matlab/partial_information/PCL_resol.m
@@ -89,7 +89,9 @@ else
             % nonlinear models
             if max(abs(feval(fh,dr.ys,[oo_.exo_steady_state; ...
                                     oo_.exo_det_steady_state], M_.params))) > options_.dynatol.f
-                [dr.ys,check1] = dynare_solve(fh,dr.ys,options_,...
+                opt = options_;
+                opt.jacobian_flag = 0;
+                [dr.ys,check1] = dynare_solve(fh,dr.ys,opt,...
                                               [oo_.exo_steady_state; ...
                                     oo_.exo_det_steady_state], M_.params);
             end
diff --git a/matlab/partial_information/dr1_PI.m b/matlab/partial_information/dr1_PI.m
index 4230828b4..0475d4bef 100644
--- a/matlab/partial_information/dr1_PI.m
+++ b/matlab/partial_information/dr1_PI.m
@@ -87,7 +87,9 @@ if options_.ramsey_policy && options_.ACES_solver == 0
     end
     old_solve_algo = options_.solve_algo;
     %  options_.solve_algo = 1;
-    oo_.steady_state = dynare_solve('ramsey_static',oo_.steady_state,options_,M_,options_,oo_,it_);
+    opt = options_;
+    opt.jacobian_flag = 0;
+    oo_.steady_state = dynare_solve('ramsey_static',oo_.steady_state,opt,M_,options_,oo_,it_);
     options_.solve_algo = old_solve_algo;
     [junk,junk,multbar] = ramsey_static(oo_.steady_state,M_,options_,oo_,it_);
     [jacobia_,M_] = ramsey_dynamic(oo_.steady_state,multbar,M_,options_,oo_,it_);
diff --git a/matlab/perfect-foresight-models/private/simulation_core.m b/matlab/perfect-foresight-models/private/simulation_core.m
index ecf5c8964..81674e4a2 100644
--- a/matlab/perfect-foresight-models/private/simulation_core.m
+++ b/matlab/perfect-foresight-models/private/simulation_core.m
@@ -79,7 +79,7 @@ else
                 illi = illi(:,2:3);
                 [i_cols_J1,junk,i_cols_1] = find(illi(:));
                 i_cols_T = nonzeros(M_.lead_lag_incidence(1:2,:)');
-                [y,info] = dynare_solve(@perfect_foresight_problem,z(:),1, ...
+                [y,info] = dynare_solve(@perfect_foresight_problem,z(:),options_, ...
                                  str2func([M_.fname '_dynamic']),y0,yT, ...
                                  oo_.exo_simul,M_.params,oo_.steady_state, ...
                                  options_.periods,M_.endo_nbr,i_cols, ...