%TVCCI_SIMU.M
%Simulates equilibrium dynamics in an economy with an individual time-varying collateral constraint. For a description of the underlying economic model, see ``Individual Versus Aggregate Collateral Constraints and the Overborrowing Syndrome,'' by Martin Uribe, April 2006.
%(c) Martin Uribe, April 27, 2006 

function [dm, A, Z, R, LA, H, Y, C, BETTA,Q, u] = tvcci_simu(zL, zH, DL, DH, aL, aH, a1L, a1H, n, RSTAR, SIG, OMEGA, BETTA1, ALFA, PAIZ, KAPA, Kss);



rand('state',sum(100*clock)); %sets a different state each time.

%Length of simulation
T = 5000;
Tcut = 1000;%number of periods thrown out in computing the test

%Usable periods
Tleft = T-Tcut;

A0 = 5.10;%initial debt position 

b=0;h=0;aa=0;%initial conditions for construction of den Haan and Marcet (RES, 1994) accuracy test

%initial state
x = rand<0.5; 
sL = x; %low state
sH = 1-x; %high state

for t=1:T

%External debt
A(t,1) = A0;

%Update state of nature
x = rand<PAIZ(1,1);
sL = sL * x + (1-sL) * (1-x);
sH = sH * x + (1-sH) * (1-x);

%Update productivity shock 
Z(t,1) = sL * zL + sH * zH;

%Update parameter vector
a = sL * aL + sH * aH;
a1 = sL * a1L + sH * a1H;

%Assume  country premium is nil
R(t,1) = RSTAR;

%Regressors
X(t,1:n) = chebyreg2(A(t,1),DL,DH,n);

%Marginal utility of consumption
LA(t,1) = exp(X(t,:)*a) .* R(t,1);

%Price of land
Q(t,1) = exp(X(t,:)*a1);

%Hours
H(t,1) = ( (1-ALFA) * exp(Z(t)) * Kss^(ALFA) ) .^(1/(OMEGA+ALFA-1));

%Output
Y(t,1) = exp(Z(t)) * Kss^ALFA * H(t).^(1-ALFA);

%Consumption
C(t,1) = LA(t).^(-1/SIG) + H(t)^OMEGA/OMEGA;

%Next period's external debt
Ap(t,1) = RSTAR .* (A(t) + C(t) - Y(t));

%Ensure debt is above lower bound
Ap(t,1) = max(Ap(t),DL);

%Ensure debt is below upper bound and that it doesn't violate the collateral constraint
Ap(t,1) = min([Ap(t) DH KAPA*Q(t)*Kss]);

%consumption
C(t) = Ap(t,1) / RSTAR - A(t)  + Y(t);

%Marginal utility of consumption
LA(t,1) = (C(t) - H(t)^OMEGA/OMEGA)^(-SIG);

%Subjective discount  factor
BETTA(t,1) = (1 + C(t) - H(t)^OMEGA/OMEGA)^(-BETTA1);

%Interest Rate
R(t,1) = LA(t) / exp(X(t,:)*a); 

%Lagrange multiplier on collateral constraint
XI(t,1) = 1/RSTAR - 1/R(t,1); 

%Update debt initial condition
A0 = Ap(t);

%Prepare for computing denHaan-Marcet (RES, 1994) statistic
if t>Tcut;
tt = t-Tcut;

%Error in bond Euler equation
u(tt,1) = LA(t-1)-R(t-1).* BETTA(t-1).* LA(t);

%Error in land Euler equation
u1(tt,1) = LA(t-1)*Q(t-1) - BETTA(t-1) / (1-KAPA*XI(t-1)) * LA(t) * (Q(t) + exp(Z(t)) * ALFA * (H(t)/Kss)^(1-ALFA));

%Instruments
h = [A(t-5:t-1)' 1 Z(t-5:t-1)'];

%Orthogonality conditions
UH = kron([u(tt) u1(tt)],h);

%Mean of orthogonality conditions
b = b + UH/Tleft  ;
%average of var/cov of orthogonality conditions
aa = aa + UH'*UH/Tleft;
end %if t>Tcut
end %for t=1:T

%denHaan-Marcet (RES, 1994) statistic
dm = Tleft * b*inv(aa)*b';