%loop version for tabe5 (log linearized model, permanent shock, lump sum taxation)

var w, p, n, c, m, i, g, r, y 

varexo eta, e;

parameters phi, beta, gamma, sigma, Y_ss, i_ss, thetam, rhog;

set_param_value('beta',beta);
set_param_value('gamma',gamma);
set_param_value('phi',phi);
set_param_value('sigma',sigma);
set_param_value('Y_ss',Y_ss);
set_param_value('rhog',rhog)

i_ss=(1/beta)-1
thetam=sigma/(phi*(1-gamma)+sigma)


model(linear);
w-p=phi*n+sigma*c
m-p=sigma*c-(1/1+i_ss)i
sigma*c=sigma*c(+1)-i+(p(+1)-p)
w-p=0
y=(sigma/(phi*(1-gamma)+sigma))*g
c=(((-phi)*thetam)/sigma)*g
i=r+p(+1)-p
g=(rhog*g(-1)+eta)
r=(1-rhog)*phi*thetam*g
end;

initval;
p=0;
r=0;
i=0;
n=0;
y=0;
w=0;
c=0;
g=0;
m=0;
e=0;
eta=0;
end;
steady;

endval;
p=p_ss;
r=r_ss;
i=i_ss;
n=n_ss;
y=y_ss;
w=w_ss;
c=c_ss;
e=0.1;
eta=-e;
end;
steady;

check;

simul(periods=100);