diff --git a/matlab/backward/simul_backward_nonlinear_model.m b/matlab/backward/simul_backward_nonlinear_model.m
index c13b02045bd67639f677e16ee6c3e48911d51800..fd5f8af90bc1e56db0f0308c23b0c1c7d36d3d8f 100644
--- a/matlab/backward/simul_backward_nonlinear_model.m
+++ b/matlab/backward/simul_backward_nonlinear_model.m
@@ -65,7 +65,7 @@ if nargin<6
     DynareOutput.exo_simul = zeros(sample_size,number_of_shocks);
     DynareOutput.exo_simul(:,positive_var_indx) = DynareOutput.bnlms.shocks;
     if isfield(DynareModel,'exo_histval') && ~ isempty(DynareModel.exo_histval)
-        DynareOutput.exo_simul = [DynareModel.exo_histval; DynareOutput.exo_simul];
+        DynareOutput.exo_simul = [transpose(DynareModel.exo_histval); DynareOutput.exo_simul];
     else
         DynareOutput.exo_simul = [zeros(1,number_of_shocks); DynareOutput.exo_simul];
     end