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

%% DO NOT CHANGE THIS PART.
%%
%% Here we load the values of the deep parameters in a loop.
%%
if isfield(M_,'param_nbr') == 1
    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
%%
%% END OF THE FIRST MODEL INDEPENDENT BLOCK.


%% THIS BLOCK IS MODEL SPECIFIC.
%%
%% Here the user has to define the steady state.
%%

Rh_1=1.0/delta_1;
Rh_2=1.0/delta_2;

Rf_1=Rh_1;
Rf_2=Rh_2;

ret_1=(Rh_1-1.0)+dep_bar_1;
ret_2=(Rh_2-1.0)+dep_bar_2;

K_1=(ret_1/alpha)^(1.0/(alpha-1.0));
K_2=(ret_2/alpha)^(1.0/(alpha-1.0));

kh_1=K_1;
kf_1=0.0;
kh_1=K_2;
kf_2=0.0;

a_1 = kh_1/(kh_1+kf_1);
a_2 = kh_2/(kh_2+kf_2);

wage_1 = (1.0-alpha)*((K_1/L_bar)^alpha);
wage_2 = (1.0-alpha)*((K_2/L_bar)^alpha);

tau_l_1=(G_1)/wage_1;
tau_l_2=(G_2)/wage_2;

cash_1=(1.0-tau_l_1)*wage_1*L_bar+Rh_1*kh_1+Rf_1*kf_1;
cash_2=(1.0-tau_l_2)*wage_2*L_bar+Rh_2*kf_2+Rf_2*kh_2;

c_1 = cash_1-kh_1-kf_1;
c_2 = cash_2-kh_2-kf_2;

Y_1= (K_1^alpha)*(L_bar^(1-alpha));
Y_2= (K_2^alpha)*(L_bar^(1-alpha));

dep_1=dep_bar_1;
dep_2=dep_bar_2;

%%
%% END OF THE MODEL SPECIFIC BLOCK.


%% DO NOT CHANGE THIS PART.
%%
%% Here we define the steady state values of the endogenous variables of
%% the model.
%%

for iter = 1:length(M_.params)
  eval([ 'M_.params(' num2str(iter) ') = ' M_.param_names(iter,:) ';' ])
end

if isfield(M_,'param_nbr') == 1

if isfield(M_,'orig_endo_nbr') == 1
NumberOfEndogenousVariables = M_.orig_endo_nbr;
else
NumberOfEndogenousVariables = M_.endo_nbr;                   % Number of endogenous variables.
end
ys = zeros(NumberOfEndogenousVariables,1);                   % Initialization of ys (steady state).
for i = 1:NumberOfEndogenousVariables                        % Loop...
  varname = deblank(M_.endo_names(i,:));                     % Get the name of endogenous variable i
  eval(['ys(' int2str(i) ') = ' varname ';']);               % Get the steady state value of this variable.
end                                                          % End of the loop.
else
ys=zeros(length(lgy_),1);
for i = 1:length(lgy_)
    ys(i) = eval(lgy_(i,:));
end
check = 0;
end
end

%%
%% END OF THE SECOND MODEL INDEPENDENT BLOCK.