// baseline rbc model 

var c r w n k i y a;

varexo epsilon ;

parameters bet del khi eta alp rho nss;
bet=.985;
del=.025;
nss=.33;
eta=1;
alp=.33;
rho=.95;
khi=(1-alp)*(1-nss)^eta/nss*(1/bet-1+del)/(1/bet-1+del-del*alp);


model;
1/c = bet*(r(+1)+1-del)/c(+1);
w = khi*c/(1-n)^eta;
k = (1-del)*k(-1)+i;
y = a*k(-1)^alp*n^(1-alp);
log(a) = rho*log(a(-1))+epsilon;
w = (1-alp)*y/n;
r = alp*y/k(-1);
y = c+i;
end;

steady_state_model;
a = 1;
r = 1/bet-1+del;
n = nss;
k = (alp/r)^(1/(1-alp))*n;
y = k^alp*n^(1-alp);
w = (1-alp)*y/n;
i = del*k;
c = y-i;
end;

steady;

resid;
check;

shocks;
var epsilon; stderr .009;
end;


//options_.loglinear=1;
stoch_simul(order=1,periods=1000,irf=30,hp_filter=1600);
//stoch_simul(order=2,periods=1000);



break
stoch_simul(irf=30);
M=nan(5,30);
ii=1;
alphavec=nan(5,1);
for alphai=0.18:.04:.35
	alp=alphai;
	khi=(1-alp)*(1-nss)^eta/nss*(1/bet-1+del)/(1/bet-1+del-del*alp);
	stoch_simul(order=1,irf=30);
	M(ii,:)=oo_.irfs.n_epsilon;
    alphavec(ii)= alphai;
    ii=ii+1;
end
close all
figure
plot(M'),grid on,legend(num2str(alphavec))
title('N')