// Modelo Macroeconomia Estocastica
 
var y1, y2, c, k, k1, k2, z, r;
varexo ze;

parameters beta, delta, a1, a2, pz;

a1 = 0.36;
a2 = 0.5;
beta  = 0.99;
delta = 0.05;
pz = 0.8;

model;
beta*(c/c(+1))*(r(+1)+(1-delta))=1;
y1 = exp(z)*(k1(-1)^a1);
y2 = exp(z)*(k2(-1)^a2);
r = a1*exp(z)*(k1(-1)^(a1-1));
r = a2*exp(z)*(k2(-1)^(a2-1));
k(-1) = k1(-1) + k2(-1);
k = y1+y2-c+(1-delta)*k(-1);
z = pz*z(-1) + ze;
end;

initval;
y1 = 0.373321;
y2 = 0.240949;
c = 0.644976;
k1 = 0.064765;
k2 = 0.0580565;
k = 0.122822;
r = 2.07513;
z = 0;
ze = 0;
end;

steady;

shocks;
var ze = 0.0016;
end;

stoch_simul(periods = 20000, irf=100) k1, k2, k, y1, y2, r;