function [ys,check] = rbc0_steadystate(ys,exe)
global M_

%% DO NOT CHANGE THIS PART.
%%
%% Here we load the values of the deep parameters in a loop.
%%
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;
%%
%% END OF THE FIRST MODEL INDEPENDENT BLOCK.

beta=.99;
gama=2;
tau=0;
phi=1;
delta=.025;
alpha=.33;
kesi=10;
theta=0.75;
phi2=.9;
phi1=1.5;
rho=.9;
sigma=.002;

r=1/beta+delta-1;
rf=1/beta-1;
p=1;
ps=p;
%lp=0;
%lps=0;
mc=(kesi-1)/kesi;
a=1;
%lmc=0;

yk=kesi/(kesi-1)*(1/beta+delta-1)/alpha;
yl=yk^(alpha/(alpha-1));
w=(1-alpha)*(kesi-1)/kesi*yl;
cy=1-delta*1/yk;
y=(1-delta*1/yk)/((1/phi*yl^tau*w)^(1/gama));

y=y^(1/(-1-tau/gama));
c=cy*y;
k=y/yk;
l=y/yl;
S=c^(-gama)*p^kesi*y*mc/(1-beta*theta);
D=c^(-gama)*p^kesi*y/(1-beta*theta);
e=0;

%% DO NOT CHANGE THIS PART.
%%
%% Here we define the steady state values of the endogenous variables of
%% the model.
%%
NumberOfEndogenousVariables = M_.endo_nbr;                    % Number of endogenous variables.
ys = zeros(NumberOfEndogenousVariables,1);                    % Initialization of ys (steady state).
for iiii = 1:NumberOfEndogenousVariables                         % Loop...
  varname = deblank(M_.endo_names(iiii,:));                      %    Get the name of endogenous variable i.                     
  eval(['ys(' int2str(iiii) ') = ' varname ';']);                %    Get the steady state value of this variable.
end                                                           % End of the loop.
%%
%% END OF THE SECOND MODEL INDEPENDENT BLOCK.