diff --git a/matlab/dsge_simulated_theoretical_conditional_variance_decomposition.m b/matlab/dsge_simulated_theoretical_conditional_variance_decomposition.m
index 23cdfdff190133667aaf380385b2687096435a02..c85ca9b1182d66d08777a155fe281b53b535af76 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 50001e153b4a8aff5bd39ff696dc7e8aaa381619..e77af705b2d23d14ec710381ffd87819860a6f44 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 86660c86d994103d2f8b69985e59a99d3791f1a9..d23fafc05b182072ecd731d996891cbfe44e8970 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 e17865aaeb5d9e7ad49ae38e236074f32b79a96a..7d5574721b602c79623e2b3c71ade5ccea8cc7b1 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/identification_analysis.m b/matlab/identification_analysis.m
index ff773195cf7810d42f5618469f271a4b3b0f7140..5f2b47660e73354c3306b91669e7a95998704c1d 100644
--- a/matlab/identification_analysis.m
+++ b/matlab/identification_analysis.m
@@ -133,7 +133,7 @@ no_identification_minimal     = options_ident.no_identification_minimal;
 no_identification_spectrum    = options_ident.no_identification_spectrum;
 
 %Compute linear approximation and fill dr structure
-[oo_.dr,info,M_,options_,oo_] = resol(0,M_,options_,oo_);
+[oo_.dr,info,M_,options_,oo_] = compute_decision_rules(M_,options_,oo_);
 
 if info(1) == 0 %no errors in solution
     % Compute parameter Jacobians for identification analysis
diff --git a/matlab/reversed_extended_path.m b/matlab/reversed_extended_path.m
index e65204767aa4b2926895cdca02bd538cc733a7bf..748423d992c879d4c05e120da2c9a44f2bc962c6 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 c42e1ced80fdc769dfe10a731cc0e2a45483c7e9..dd336f2fb4fabc6d03ba0c0b24e77a37f8a7e41c 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 8a892cc5ae82d32dfd6160a3cea5012e81756d7e..152ac6b5cb13cb39be2baf5db8c881719c7e9a8d 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')