diff --git a/matlab/dsge_simulated_theoretical_conditional_variance_decomposition.m b/matlab/dsge_simulated_theoretical_conditional_variance_decomposition.m
index 23cdfdff1..c85ca9b11 100644
--- a/matlab/dsge_simulated_theoretical_conditional_variance_decomposition.m
+++ b/matlab/dsge_simulated_theoretical_conditional_variance_decomposition.m
@@ -142,7 +142,7 @@ for file = 1:NumberOfDrawsFiles
             dr = temp.pdraws{linee,2};
         else
             M_=set_parameters_locally(M_,temp.pdraws{linee,1});
-            [dr,info,M_,options_,oo_] = resol(0,M_,options_,oo_);
+			[dr,info,M_,options_,oo_] = compute_decision_rules(M_,options_,oo_)
         end
         if first_call
             endo_nbr = M_.endo_nbr;
diff --git a/matlab/dsge_simulated_theoretical_correlation.m b/matlab/dsge_simulated_theoretical_correlation.m
index 50001e153..e77af705b 100644
--- a/matlab/dsge_simulated_theoretical_correlation.m
+++ b/matlab/dsge_simulated_theoretical_correlation.m
@@ -115,7 +115,7 @@ for file = 1:NumberOfDrawsFiles
             dr = temp.pdraws{linee,2};
         else
             M_=set_parameters_locally(M_,temp.pdraws{linee,1});
-            [dr,info,M_,options_,oo_] = resol(0,M_,options_,oo_);
+            [dr,info,M_,options_,oo_] = compute_decision_rules(M_,options_,oo_);
         end
         if ~options_.pruning
             tmp = th_autocovariances(dr,ivar,M_,options_,nodecomposition);
diff --git a/matlab/dsge_simulated_theoretical_covariance.m b/matlab/dsge_simulated_theoretical_covariance.m
index 86660c86d..d23fafc05 100644
--- a/matlab/dsge_simulated_theoretical_covariance.m
+++ b/matlab/dsge_simulated_theoretical_covariance.m
@@ -113,7 +113,7 @@ for file = 1:NumberOfDrawsFiles
             dr = temp.pdraws{linee,2};
         else
             M_=set_parameters_locally(M_,temp.pdraws{linee,1});
-            [dr,info,M_,options_,oo_] = resol(0,M_,options_,oo_);
+            [dr,info,M_,options_,oo_] = compute_decision_rules(M_,options_,oo_);
         end
         if ~options_.pruning
             tmp = th_autocovariances(dr,ivar,M_,options_,nodecomposition);
diff --git a/matlab/dsge_simulated_theoretical_variance_decomposition.m b/matlab/dsge_simulated_theoretical_variance_decomposition.m
index e17865aae..7d5574721 100644
--- a/matlab/dsge_simulated_theoretical_variance_decomposition.m
+++ b/matlab/dsge_simulated_theoretical_variance_decomposition.m
@@ -141,7 +141,7 @@ for file = 1:NumberOfDrawsFiles
             dr = temp.pdraws{linee,2};
         else
             M_=set_parameters_locally(M_,temp.pdraws{linee,1});
-            [dr,info,M_,options_,oo_] = resol(0,M_,options_,oo_);
+            [dr,info,M_,options_,oo_] = compute_decision_rules(M_,options_,oo_);
         end
         if file==1 && linee==1
             [tmp, stationary_vars] = th_autocovariances(dr,ivar,M_,options_,nodecomposition);
diff --git a/matlab/reversed_extended_path.m b/matlab/reversed_extended_path.m
index e65204767..748423d99 100644
--- a/matlab/reversed_extended_path.m
+++ b/matlab/reversed_extended_path.m
@@ -48,7 +48,7 @@ steady_;
 
 %  Compute the first order perturbation reduced form.
 old_options_order = options_.order; options_.order = 1;
-[dr,info,M_,options_,oo_] = resol(0,M_,options_,oo_);
+[dr,info,M_,options_,oo_] = compute_decision_rules(M_,options_,oo_);
 oo_.dr = dr;
 options_.order = old_options_order;
 
diff --git a/matlab/selec_posterior_draws.m b/matlab/selec_posterior_draws.m
index c42e1ced8..dd336f2fb 100644
--- a/matlab/selec_posterior_draws.m
+++ b/matlab/selec_posterior_draws.m
@@ -142,7 +142,7 @@ if info
             pdraws(linee,1) = {x2(SampleAddress(i,4),:)};
             if info-1
                 set_parameters(pdraws{linee,1});
-                [dr,info,M_,options_,oo_] = resol(0,M_,options_,oo_);
+                [dr,info,M_,options_,oo_] = compute_decision_rules(M_,options_,oo_);
                 pdraws(linee,2) = { dr };
             end
             old_mhfile = mhfile;
diff --git a/matlab/simulated_moment_uncertainty.m b/matlab/simulated_moment_uncertainty.m
index 8a892cc5a..152ac6b5c 100644
--- a/matlab/simulated_moment_uncertainty.m
+++ b/matlab/simulated_moment_uncertainty.m
@@ -71,7 +71,7 @@ else
     logged_steady_state_indicator=0;
 end
 
-[dr,info,M_,options_,oo_] = resol(0,M_,options_,oo_);
+[dr,info,M_,options_,oo_] = compute_decision_rules(M_,options_,oo_);
 oo_.dr=dr;
 if info(1)
     fprintf('\nsimulated_moment_uncertainty: model could not be solved')