function [ys, check] = dsge11_steadystate(ys, exe)

%% Initialize parameters
alpha = 0.33;
betab = 0.945;
betae = 0.94;
betah = 0.99; 
delta = 0.025;
rho_a = 0.98; 
rho_g = 0.95; 
xi_h = 0.35;
xi_e = 0.027;
xi_b = 0.10;
eta = 0.95;
sigma_a= 0.01; 
sigma_g =0.01; 
chi = 15.02;
theta=0.96; 
theta_i = 4;
Abar = 1;
gbar = 2;


%% Calculate steady state equations
N= 0.33;
q = 1;
rd = (1/betah) - 1;
rl = (1/betab) - 1;
K = ((1/betae + delta - 1)/alpha)^(1/(alpha-1));
I = delta * K;
Y = Abar * K^alpha * N^(1-alpha);
gbar = 0.2 * Y;  
t = gbar;
W = (1-alpha) * (K/N)^alpha; 
R = 1/betae;
L = q * K * theta;
d = eta * L;
ce = Y - t - I - wn - rl * L; 
cb = rl * L - rd * d - t; 
ch = wn + rd * d - t; 
me = (xi_e * ce) / (1 - betae);
mb = (xi_b * cb) / (1 - betab);
mh = (xi_h * ch) / (1 - betah);
m = me + mh + mb;
tb = 0.1 * t;
th = 0.5 * t;
te = 0.4 * t;
Abar =1;
%% End of SS equations
check = 0; 

%% Assign output
ys = [Abar, q, rd, rl, K, I, Y, gbar, W, R, N, L, d, wn, ce, cb, ch, me, mb, mh, m, tb, th, te];

end
