diff --git a/models/pshock.mod b/models/pshock.mod
index e442e17fe282e67ff24ad0cff5fd6eaac9b33f19..49481aeca17a49c038d5dae0fb20823328a304b3 100644
--- a/models/pshock.mod
+++ b/models/pshock.mod
@@ -46,6 +46,8 @@
 
 @#include "modeqs.mod"
 
+for i=1:@{ITERATIONS}
+
 load_params_and_steady_state('eagle_steady.txt');
 
 initval;
@@ -58,3 +60,5 @@ end;
 
 perfect_foresight_setup(periods=@{PERIODS});
 perfect_foresight_solver(no_homotopy, markowitz=@{MARKOWITZ_VALUE}, stack_solve_algo = @{STACK_SOLVE_ALGO_VALUE}, solve_algo = @{STEADY_SOLVE_ALGO_VALUE});
+
+end;
diff --git a/models/tshock.mod b/models/tshock.mod
index 794d51733fd43069631700f7f4023ad45dae2671..71b706d86c1d72f144394b311d23b0cdcf5fb906 100644
--- a/models/tshock.mod
+++ b/models/tshock.mod
@@ -46,6 +46,8 @@
 
 @#include "modeqs.mod"
 
+for i=1:@{ITERATIONS}
+
 load_params_and_steady_state('eagle_steady.txt');
 
 initval;
@@ -86,3 +88,5 @@ end;
 
 perfect_foresight_setup(periods=@{PERIODS_VALUE});
 perfect_foresight_solver(no_homotopy, markowitz=@{MARKOWITZ_VALUE}, stack_solve_algo = @{STACK_SOLVE_ALGO_VALUE}, solve_algo = @{STEADY_SOLVE_ALGO_VALUE});
+
+end;