% computes the steady state 


function [ys,check] = model_FA_NL_ZLB1_steadystate(ys,exe)

global M_ lgy_


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

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

lambda     	=	lambdabar;
c          	=	cbar	;
R          	=	Rbar	;
pip        	=	pipbar	;
rreal      	=	rrealbar;
lh         	=	lhbar	;
w          	=	wbar	;
piw        	=	piwbar	;
y          	=	ybar	;
l          	=	lbar	;
k          	=	kbar	;
s          	=	sbar	;
we         	=	webar	;
z          	=	zbar	;
Q1         	=	Q1bar	;
Q2         	=	Q2bar	;
Pstar      	=	Pstarbar;
Q1w        	=	Q1wbar	;
Q2w        	=	Q2wbar	;
Wstar      	=	Wstarbar;
lomegatilde	=	lomebar;
clag       	=	cbar	;
Fome       	=	Fomebar	;
Fnome      	=	Fnomebar;
omegatilde 	=	omebar	;
SIG        	=	SIGbar	;
muG        	=	muGbar	;
ffome      	=	ffomebar;
deltau     	=	deltaubar;
deltauu    	=	deltauubar;
adji       	=	adjibar	;
adjii      	=	adjiibar;
nw         	=	nwbar	;
i          	=	ibar	;
ktilde     	=	ktildebar;
b          	=	bbar	;
xx         	=	xxbar	;
q          	=	qbar	;
rk         	=	rkbar	;
rr         	=	rrbar	;
LAM        	=	LAMbar	;
v          	=	vbar	;
ce         	=	cebar	;
SIGp       	=	SIGpbar	;
muGp       	=	muGpbar	;
gammat     	=	gammabar;
pref       	=	prefbar	;
em         	=	embar	;
techno     	=	technobar;
chil       	=	chilbar	;
g          	=	gbar	;
sigome     	=	sigmaWbar;
a          	=	abar	;

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);
for kk = 1:NumberOfEndogenousVariables
  varname = deblank(M_.endo_names(kk,:));
  eval(['ys(' int2str(kk) ') = ' varname ';'])

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