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

if isfield(M_,'param_nbr') == 1
NumberOfParameters = M_.param_nbr;
for i = 1:NumberOfParameters
  paramname = deblank(M_.param_names(i,:));
  eval([ paramname ' = M_.params(' int2str(i) ');']);
end
check = 0;
end

BETTA =0.99;
LAMBDAYSTAR=0.5;
PHIy=0.3;
PHIrstar=0.5;
PHIpi=0.25;
RSTARSS=0.001;
YSS=1;
LSS=0.33;
ETTA=1;
GAMA=0.5;

%express y
y=((l)*(GAMA*(LSS)^(GAMA-1)*LSS))/YSS;
%express c
y=c;
%express rstar
c=-BETTA+c(+1)-rstar(-1);
%express l
c=-(ETTA-1)*l;
%express pistar
rstar=PHIrstar*RSTARSS+(1-PHIrstar)*((1+PHIpi)*pistar+PHIy*(y-y(-1)));
%express yh
yh=l;
%express pistar
pistar=BETTA*pistar(+1);
%express yhstar
yhstar=LAMBDAYSTAR*yhstar(-1)+vyhstar;
%express bstar
exp(yf)=exp(yhstar)+exp(bstar)-exp(bstar);


%%
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;
end
ys = zeros(NumberOfEndogenousVariables,1);
if options_.linear
    
else
for i = 1:NumberOfEndogenousVariables
  varname = deblank(M_.endo_names(i,:));
  eval(['ys(' int2str(i) ') = ' varname ';']);
end
end

else
ys=zeros(length(lgy_),1);
for i = 1:length(lgy_)
    ys(i) = eval(lgy_(i,:));
end
check = 0;
end