diff --git a/matlab/perfect-foresight-models/perfect_foresight_solver.m b/matlab/perfect-foresight-models/perfect_foresight_solver.m
index 0ffe5187687372e4aa771b866e23389bd7e0d684..d303d59db43642f0c0b8b50a814577849ea8ce2b 100644
--- a/matlab/perfect-foresight-models/perfect_foresight_solver.m
+++ b/matlab/perfect-foresight-models/perfect_foresight_solver.m
@@ -234,7 +234,7 @@ else
     initial_period = options_.initial_period;
 end
 
-ts = dseries(transpose(oo_.endo_simul), initial_period, M_.endo_names);
+ts = dseries([transpose(oo_.endo_simul(1:M_.orig_endo_nbr,:)), oo_.exo_simul], initial_period, [M_.endo_names(1:M_.orig_endo_nbr); M_.exo_names]);
 assignin('base', 'Simulated_time_series', ts);
 if oo_.deterministic_simulation.status
     oo_.gui.ran_perfect_foresight = true;