From 38827443f99ec052f519dd58ee2875d21a8cbebf Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?S=C3=A9bastien=20Villemot?= <sebastien.villemot@ens.fr>
Date: Thu, 27 Oct 2011 16:45:17 +0200
Subject: [PATCH] Fix simult_.m when order=3 but k_order_solver has not been
 specified

---
 matlab/simult_.m | 6 ++++++
 1 file changed, 6 insertions(+)

diff --git a/matlab/simult_.m b/matlab/simult_.m
index 46823f87e..1bcc12725 100644
--- a/matlab/simult_.m
+++ b/matlab/simult_.m
@@ -39,6 +39,12 @@ iter = size(ex_,1);
 y_ = zeros(size(y0,1),iter+M_.maximum_lag);
 y_(:,1) = y0;
 
+% stoch_simul sets k_order_solver=1 if order=3, but does so only locally, so we
+% have to do it here also
+if options_.order == 3
+    options_.k_order_solver = 1;
+end
+
 if ~options_.k_order_solver
     if iorder==1
         y_(:,1) = y_(:,1)-dr.ys;