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

global M_

NumberOfParameters = M_.param_nbr;                            % Number of deep parameters.
for i = 1:NumberOfParameters                                  % Loop...
  paramname = deblank(M_.param_names(i,:));                   %    Get the name of parameter i. 
  eval([ paramname ' = M_.params(' int2str(i) ');']);         %    Get the value of parameter i.
end                                                           % End of the loop.  

check = 0;

Y_SS = ((eeta-1)/eeta)^(1/ttau)*gbar;
ppi_SS = ppibar;
b_SS = (eeta-1)/eeta*(((eeta-1)/eeta)^(1/ttau)*gbar)/(1-nbeta*ggamma);
c_SS = ((eeta-1)/eeta)^(1/ttau);
R_SS = ppibar*exp(oomega)/nbeta;
v_SS = 1;
a_SS = 0;
g_SS = gbar;
A_SS = exp(oomega);
H_SS = ((eeta-1)/eeta)^(1/ttau)*gbar;
Util_SS = c_SS^(1-ttau)/(1-ttau) - H_SS;
Welf_SS = Util_SS/(1-nbeta);
zzeta_SS = (gbar-1)/gbar;
llambda_SS = 1/((eeta-1)/eeta);
W_SS = 1/(1/(((eeta-1)/eeta)^(1/ttau))^ttau);
mc_SS = 1/(1/(((eeta-1)/eeta)^(1/ttau))^ttau);
P_SS = 1;
G_SS = ((gbar-1)/gbar)*((eeta-1)/eeta)^(1/ttau)*gbar;
Ystar_SS = (((eeta-1)/eeta)^(1/ttau))*(1/(1-((gbar-1)/gbar)));


Y = ((eeta-1)/eeta)^(1/ttau)*gbar;
ppi = ppibar;
b = (eeta-1)/eeta*(((eeta-1)/eeta)^(1/ttau)*gbar)/(1-nbeta*ggamma);
c = ((eeta-1)/eeta)^(1/ttau);
R = ppibar*exp(oomega)/nbeta;
v = 1;
a = 0;
g = gbar;
A = exp(oomega);
H = ((eeta-1)/eeta)^(1/ttau)*gbar;
Util = c^(1-ttau)/(1-ttau) - H;
Welf = Util/(1-nbeta);
zzeta = (gbar-1)/gbar;
llambda = 1/((eeta-1)/eeta);
W = 1/(1/(((eeta-1)/eeta)^(1/ttau))^ttau);
mc = 1/(1/(((eeta-1)/eeta)^(1/ttau))^ttau);
P = 1;
G = ((gbar-1)/gbar)*((eeta-1)/eeta)^(1/ttau)*gbar;
Ystar = (((eeta-1)/eeta)^(1/ttau))*(1/(1-((gbar-1)/gbar)));

ys = [Y
    ppi
    b
    c
    R
    v
    a
    g
    A
    H
    Util
    Welf
    zzeta
    llambda
    W
    mc
    P
    G
    Ystar];