function [ys,check]=benchmark_steadystate(ys,exe)
check=0;

global M_ Ahbar Afbar KhKf KhNh KfNf alpha beta gamma delta dzeta eta thetah thetaf kappa mu nu rhok rhop rhom sigma tauh tauf upsilon phi phistar null nullpih nullpif nullyh nullyf chi psihbar psifbar xi tetah tetaf omegabar C1 C2 C3 C4 C5 C6 C7 C8 C9 F1 F2 F3 F4 F5;


%Evaluate SS values

NhKh=1/KhNh;
NfKf=1/KfNf;
Rh=1/beta;
rh=Rh-1;
Rf=1/beta;
rf=Rf-1;
psih=psihbar;
psif=psifbar;
Ah=Ahbar;
Af=Afbar;
lambdah=rh;
lambdaf=rf;
Vh=1;
Vf=1;
Zh=Rh-(1-delta);
Zf=Rf-(1-delta);
Kh=(Zh*C1-delta)^((-eta)/(eta+null))*Zh^((-(alpha+null))/((1-alpha)*(eta+null)))*(phi*C4*C6+(1-phistar)*C5*C7*KhKf^(null/eta))^(eta/(eta+null));
Sh=Kh;
Kf=(Zh*C2-delta)^((-eta)/(eta+null))*Zh^((-(alpha+null))/((1-alpha)*(eta+null)))*(phistar*C5*C7+(1-phi)*C4*C6*KhKf^((-null)/eta))^(eta/(eta+null));
Sf=Kf;
Ch=C4*C6*Zh^(-(alpha+null)/((1-alpha)*eta))*Kh^(-null/eta);
Chh=phi*Ch;
Chf=(1-phi)*Ch;
Cf=C5*C7*Zh^(-(alpha+null)/((1-alpha)*eta))*Kf^(-null/eta);
Cff=phistar*Cf;
Cfh=(1-phistar)*Cf;
Lh=C4^(eta/(null+alpha))*Kh^(alpha/(null+alpha))*Ch^(-eta/(null+alpha));
Lf=C5^(eta/(null+alpha))*Kf^(alpha/(null+alpha))*Cf^(-eta/(null+alpha));
Yh=Ah*Kh^alpha*Lh^(1-alpha);
Yf=Af*Kf^alpha*Lf^(1-alpha);
Gh=thetah*Yh;
Gf=thetaf*Yf;
Ih=delta*Kh;
If=delta*Kf;
Wh=chi*Lh^null*Ch^eta;
Wf=chi*Lf^null*Cf^eta;
Th=Gh;
Tf=Gf;
Nh=NhKh*Kh;
Nf=NfKf*Kf;
Vhs=1;
Vfs=1;
Qh=1;
Qf=1;
MCh=(kappa-1)/kappa;
MCf=(kappa-1)/kappa;
Dhh=Kh-Nh;
Dhf=(beta/(1-beta))*(Ch-Wh*Lh-(1/beta)*Kh+Nh+Th+Dhh);
Dff=Kf-Nf-Dhf;
Ph=1;
Pf=1;
PH=1;
PF=1;
pih=1;
pif=1;
piH=1;
piF=1;
Ex=1;
leverh=KhNh;
leverf=KfNf;

N_Endo=M_.endo_nbr;                                %Number of endogenous variables.
ys=zeros(N_Endo,1);                                %Initialization of ys (steady state).
for kk=1:N_Endo                                    %Loop...
varname = deblank(M_.endo_names(kk,:));            %Get the name of endogenous variable i.                     
eval(['ys(' int2str(kk) ') = ' varname ';']);      %Get the steady state value of this variable.
end                                                %End of the loop.

