function [ys,params,check] = probset1sol_steadystate(ys,exo,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;

%% Steady state equations
% omitting the time notation

Z = 1;

V = 1;

K = (1/alph)*(1/(bet) + del - 1)^(1 - alph);

C = (1/alph)*(1/(bet) + del - 1)^(alph/(alph - 1)) - del/(alph)*(1/(bet) + del - 1)^(1 - alph);

Y = (1/alph)*(1/(bet) + del - 1)^(alph/(alph - 1));

%sol: K* = 0.238; C* 0.64

%BK condition (not necessary, stable and instable eigenvalues (equilibrium) - >1 are unstable

for iter = 1:length(M_.params)                                               %update parameters set in the file
  eval([ 'params(' num2str(iter) ',1) = ' 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

end

