% RUN_MODEL_BasicRBC.M Solves Hansen's Basic RBC Model. 

% Solution of the model (make sure you have already executed the
% Mathematica code)

dirmodel='c:\users\prestrec\Codigos\BasicRBC\KPR\Hansen_IL\';

RHO = [0.95];
[H,Mke,dim,sstate]= kpr(dirmodel,RHO)
return
% Computation and plots of impulse response
% nombres1 = {'c1','c2','c3','s1','s2','cs1','cs2','ne1','ne2','nf1','nf2'};
nombres = {'c','h','k','\lambda','A','y','iv','LP'};

S = [0];
E = [0.1];
J = 40;

IRM = irf(H,Mke,dim,S,E,J);
plotimp(IRM,dim,nombres)

m=6;
flag='hp';
sde=0.00712;
[TSl,TSd,TS,msdev,mcorr,macor] = simul(H,Mke,dim,sstate,sde,nombres,flag,J,m)
plotfigs(TSd,dim,nombres)
plotcorrelogram(macor,dim,nombres)