%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
// Modelito con DIEGO MARADONA
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


// Name of the endogenous variables
var  c, Ny, Nr, R, N, alphad, alphay, RatioN;  


// Name of the exogenous variable
varexo exd, exy;                                    
                                             

// Name of the model parameters  
parameters  gamma, beta, X, delta, sigmaexd, phid, phiy, sigmaexy;               
    

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
// Calibration
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

gamma       = 0.2;
beta        = 0.98;
X           = 1;
delta       = 0.02;

sigmaexd    = 0.01;
phid        = 0.8;

sigmaexy    = 0.01;
phiy        = 1;


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
// Description of the model
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

model;


X/(1-Ny-Nr)=c^(-1)*(1-gamma)*(alphay)*(R)^(gamma)*Ny^(-gamma);

X/(1-Ny-Nr)=beta*(c(+1)^(-1)*(alphay(+1))*gamma*R(+1)^(gamma-1)*Ny(+1)^(1-gamma));

c=(alphay)*R^(gamma)*Ny^(1-gamma);

R=Nr(-1)+(1-delta-alphad)*R(-1);

N=Ny+Nr;
                       
alphad = phid*alphad(-1)+exd;

alphay = phiy*alphay(-1)+exy;
   
RatioN=Nr/Ny;

end;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Nyss = (1-gamma)/(X+1-gamma+beta*gamma*delta);
Nrss = ((beta*gamma*delta)/(1-gamma))*Nyss;
Rss=Nrss/delta;
css=Rss^(gamma)*Nyss^(1-gamma);
Nss = Nyss + Nrss;

initval;

Ny = Nyss;
Nr = Nrss;
R=Rss;
c=css;
N=Nss;
RatioN=Nrss/Nyss;

alphad=0;
alphay=1;


end;

resid;
steady; 
check;


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

shocks;
var exy; stderr sigmaexy;
end;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
stoch_simul(order=1, periods=1000, irf=20) c, Ny, Nr, R, N, alphad, alphay RatioN;