function [ys,check] = guosicm_steadystate(ys,exe)
global beta sigma alpha delta v fai eta share dynatol_
klratio=((1/beta-1+delta)/alpha)^(1/(alpha-1));
k=((1-alpha)*(klratio^(alpha-1)-delta)^(-sigma)*klratio^(alpha+1)/fai)^(1/(1+sigma));
l=k/klratio;
c=(klratio^(alpha-1)-delta)*k;
x=delta*k;
y=k^alpha*l^(1-alpha);
p=(alpha*k^alpha*l^(1-alpha)-x)/(1/beta-1);
s=share;

ys=[c c k k l l p p s s x x y y 0 0]';
%if max(abs(HW2cm_fff(ys))) > dynatol_
 % error('Bad steady state')
%end
check=0;