function [ys,params,check] = LRRTANK_steadystate(ys,exo_ss,M_,options)
% read out parameters to access them with their name
NumberOfParameters = M_.param_nbr;
for ii = 1:NumberOfParameters
  paramname = M_.param_names{ii};
  eval([ paramname ' = M_.params(' int2str(ii) ');']);
end
% initialize indicator
check = 0;


%% Enter model equations here
Lcbar  = Lline/3;
Lncbar = Lline/3;
Mcbar = (NU-1)/NU;
Mbar  = BETA*exp(-GAM*MU);
Rbar  = PISS/Mbar;
Lbar  = LAMBDA*Lncbar+(1-LAMBDA)*Lcbar;
Ybar = Lbar;
Wbar = Mcbar*Ybar/Lbar;
Dbar = Ybar-Wbar*Lbar;
Cncbar = Wbar*Lncbar+TAUD*Dbar/LAMBDA;
Ccbar = Wbar*Lcbar+(1-TAUD)*Dbar/(1-LAMBDA);
% Parameters and Epstein Zin
CHI0C = Wbar*(Lline-Lcbar)^CHI*Ccbar^(-GAM);
CHI0NC = Wbar*(Lline-Lncbar)^CHI*Cncbar^(-GAM);

% Convert variables to log
ng = MU;
cc = log(Ccbar);
lc = log(Lcbar);
cnc = log(Cncbar);
lnc = log(Lncbar);
l = log(Lbar);
w = log(Wbar);
d = Dbar;
y = log(Ybar);
infl = log(PISS);
a = 0;
r = log(Rbar);
mc = log(Mcbar);
x = 0;
ccg = MU;
nccg = MU;
zg = MU;

%% end own model equations

params=NaN(NumberOfParameters,1);
for iter = 1:length(M_.params) %update parameters set in the file
  eval([ 'params(' num2str(iter) ') = ' M_.param_names{iter} ';' ])
end

NumberOfEndogenousVariables = M_.orig_endo_nbr; %auxiliary variables are set automatically
for ii = 1:NumberOfEndogenousVariables
  varname = M_.endo_names{ii};
  eval(['ys(' int2str(ii) ') = ' varname ';']);
end