function [CSS, GSS, HSS, ISS, KSS, LSS, LAMDA_SS, MUK, MUX, QSS, USS, X_GSS, YSS, ZSS] = cycle_level_steadystate

global BETTA SIG ALPA DELTA0 DELTA1 MUA MUY MUI;
global DELTA2 THETA_C THETA_L RHO_Z RHO_G RHO_XG RHO_X RHO_A KAPA CHI;
%% steady state of growth rate : MUK and MUX
MUK = MUY/MUA;
MUX = MUY*MUA^(-ALPA/(ALPA-1));

%% steady state of endogenous variables other then growth rate and observable variable
X_GSS = (1/MUY)^(1/(1-RHO_XG));
ZSS = 1;
USS = 1;
QSS = 1;
KH = (QSS*(1/BETTA/MUA*(MUY)^(SIG)-(1-DELTA0))/(ZSS*USS*ALPA*(USS/MUK)^(ALPA-1)))^(1/(ALPA-1));            % steady value of k/h
YH = ZSS*(USS*KH/MUK)^(ALPA);                          % steady value of y/h
IK = (1-(1-DELTA0)/MUK);                           % steady value of i/k
CH = 0.8*YH - IK*KH;                                % steady value of c/h
F = ((1-THETA_C/MUY)^(-SIG)-THETA_C*BETTA*(MUY-THETA_C)^(-SIG))*ZSS*(1-ALPA)*(USS/MUK*KH)^(ALPA)*(1-THETA_L);
GSS = CHI*CH*((1-THETA_C/MUY)^(1-SIG)-THETA_L*BETTA*(MUY-THETA_C)^(1-SIG))+(1-THETA_L)*ZSS*(1-ALPA)*(USS/MUK*KH)^ALPA*((1-THETA_C/MUY)^(-SIG)-THETA_C*BETTA*(MUY-THETA_C)^(-SIG));
HSS = F/GSS;
KSS = KH*HSS;
YSS = YH*HSS;
CSS = CH*HSS;
LSS = 1 - HSS;
ISS = IK*KSS;
LAMDA_SS = (CSS-THETA_C*CSS/MUY)^(-SIG)*(LSS-THETA_L*LSS)^(CHI*(1-SIG))-THETA_C*BETTA*(CSS*MUY-THETA_C*CSS)^(-SIG)*(LSS-THETA_L*LSS)^(CHI*(1-SIG));
% Q_A = (ALPA*Y/K - IK)*(1-BETTA*(MU_Y)^(1-SIG));
GSS = 0.2*YSS*MUY^(1/(1-RHO_XG));