%Last updated by Dennis : 21-August-2024 
function [ys, params, check] = model2_steadystate(ys_init, exo_ss, M_, options_)



paramfile_model2

nparams = size(M_.param_names,1);
for icount = 1:nparams
    param_name = M_.param_names{icount}; M_.params(icount) = eval(param_name);
end


check=0;
%% Steady state

load ss_values
ys = [ACSS 
AHSS 
ZSS 
tauSS
uCSS 
uHSS 
pi_CSS 
upsilon_wCSS 
upsilon_wHSS 
upsilon_pSS 
Y_FSS
Q_starSS
R_starSS
barphiSS
PhiSS
thetaSS
R_polSS 
R_DSS 
R_govSS                                  
R_pvtSS 
R_BSS 
LambdaSS 
Lambda_primeSS 
Q_govSS 
Q_pvtSS
R_SSS 
p_FSS 
p_DSS 
p_D_starSS 
p2SS
p_KCSS 
M2SS 
M1SS 
MU_CSS 
MU_HSS 
qSS 
p_KHSS 
CSS 
C_FSS 
C_DSS 
N_HSS 
N_CSS 
MU_NCSS 
MU_NHSS 
w_CSS 
w_HSS 
MU_NC_primeSS 
MU_NH_primeSS 
MU_C_primeSS 
N_H_primeSS 
N_C_primeSS 
C_primeSS 
C_F_primeSS
C_D_primeSS
N2_CSS 
N2_HSS 
K_CSS
K_HSS
K_LSS
w2_CSS 
w2_HSS  
w2_C_starSS 
w2_H_starSS 
f1_CSS 
f1_HSS
f2_CSS 
f2_HSS
Y2SS 
H2SS 
YSS 
p_LSS 
I_C_hatSS 
I_H_hatSS
I_CSS 
I_HSS 
GSS 
XSS 
lambda3SS
MU_H_primeSS
JSS 
HSS 
H_primeSS
b_primeSS 
GDPSS 
b_pvtSS 
b_gov_maSS
b_pvt_maSS 
sSS 
b_pvt_fiSS 
b_gov_fiSS 
Pi_maSS 
TSS 
b_starSS 
nSS 
dSS 
phiSS 
OmegaSS 
lambda1SS 
divSS
x1SS
x2SS
GDP_obsSS 
pi_C_obsSS
p_F_obsSS 
q_obsSS 
R_pol_obsSS 
C_obsSS 
b_star_obsSS 
w_C_obsSS
w_H_obsSS 
I_C_obsSS
I_H_obsSS
N_H_obsSS
N_C_obsSS 
NX_obsSS 
];



params = M_.params;
  % Check for issues in steady-state calculation
    if any(isnan(ys)) || any(isinf(ys))
        check = 1; % Indicate that there was a problem
    end

