periods 20100;

var k, c, y, l, r, i, z;

varexo e;

parameters A, alpha, beta, delta, psi, sigma;

A     = 3.20;
alpha = 0.36; 
beta  = 0.99;
delta = 0.025; 
psi   = 0.95;  
sigma = 0.00712;



model;

A*c = (1-alpha)*exp(z)*(k(-1)^alpha)*(l^(-alpha));
1 = beta*((c/c(+1))*r(+1));
r = alpha* (k(-1)^(alpha-1))*exp(z)*(l^(1-alpha)) + 1 - delta;
y = (k(-1)^alpha)*exp(z)*(l^(1-alpha));
i = k - (1-delta)*k(-1);
i = (k(-1)^alpha)*exp(z)*(l^(1-alpha)) - c;
z = psi*z(-1) + e;
end;

initval;
k = 2;
c = 1.33;
l = 0.31;
z = 0;
e = 0;
end;

steady;

shocks;
var e = sigma^2;
end;



stoch_simul(hp_filter = 1600, order = 1);

statistic1 = 100*sqrt(diag(oo_.var(1:7,1:7)))./oo_.mean(1:7);

table('Relative standard deviations in %',strvcat('VARIABLE','REL. S.D.'),lgy_(1:7,:),statistic1,10,7,4)