e=[0.1:0.1:0.5];

for i=1:length(e)
x=e(i);

beta= 0.98;
deltak=0.072;
deltah=0.06;
sigma=0.92;
eta=0.66;
Psi=1.82;
alpha=0.32;
a=1.23;
g=0.2;
rho=0.5;
tauc=0.2;
taul=0.1;
phi=0.3;
save parameterfile.mat beta deltak rho sigma eta Psi alpha a g tauc taul deltah phi x;
dynare mod2 noclearall;
results_cell{i}=oo_.steady_state;
end