%CCRC_SS.M 
% Takes Kss and BETTA1 as given
function [RSTAR, SIG, OMEGA, BETTA1, ALFA, MU, z, PAIZ, KAPA, Kss,  QQss, Css, Hss, LAss, Ass, Yss, BETTAss] = ccrc_ss;

[RSTAR,SIG,OMEGA,BETTA1,ALFA,MU,z,PAIZ, Kss, KAPA] = cc_ss1;


WKss=1-MU*(1-(1/RSTAR)); %This is the adjustment due to the working capital constraint

Hss = ((1-ALFA) * WKss * Kss^ALFA)^(1/(ALFA+OMEGA-1));

Yss = Kss^ALFA * Hss^(1-ALFA);

Css = RSTAR^(1/BETTA1) - 1 +Hss^OMEGA/OMEGA;

COMP = Css-(Hss^OMEGA)/OMEGA; %This is a shorcut for the composite good

Ass = (Css-Yss*WKss) / (1/RSTAR - 1);

Uss = COMP^(1-SIG)/(1-SIG);

UCss = (Css-Hss^OMEGA/OMEGA)^(-SIG); 

%Marginal Utility of Consumption

LAss =  UCss;

BETTAss=(1+Css-Hss^OMEGA/OMEGA)^(-BETTA1);

QQss = (1/(RSTAR-1))*ALFA*((Hss/Kss)^(1-ALFA))*WKss;


%%%%%%%%%%%%%%
%CC_SS1.M 
% Takes Qss=1 and Ass/Yss=ALFA_a  given

function [RSTAR,SIG,OMEGA,BETTA1,ALFA,MU,z,PAIZ, Kss, KAPA] = cc_ss1;

KAPA = .10;
RSTAR = 1.04;
SIG = 2;
OMEGA = 1.455;
ALFA_a=0.5;
ALFA = 0.32;
MU = 0;%0.2;

%Productivity shock
ZH = 0.0258;
ZL = -ZH;
z = [ZL;ZH];

%Transition Probability matrix
PAIZHH = 0.71;
PAIZLL=PAIZHH;
PAIZHL=1-PAIZHH;
PAIZLH=1-PAIZLL;
PAIZ = [PAIZLL PAIZLH;PAIZHL PAIZHH];

WKss=1-MU*(1-(1/RSTAR)); %This is the adjustment due to the working capital constraint

Kss= ( ((RSTAR-1)/(ALFA*WKss)) * (((1-ALFA)*WKss)^((ALFA-1)/(ALFA+OMEGA-1))))^((ALFA+OMEGA-1)/((OMEGA-1)*(ALFA-1)));

Hss = ((1-ALFA) * WKss * Kss^ALFA)^(1/(ALFA+OMEGA-1));

Yss = Kss^ALFA * Hss^(1-ALFA);

BETTA1 = log(RSTAR) / log(Yss*WKss + ALFA_a*Yss*(1/RSTAR-1) - (Hss^OMEGA)/OMEGA + 1);