function F = SMM_steady(ss, ss0, kns, lns, cns, Gs, thetaa, phhi, betta, deltta, alphha, gy);
options = optimset('fsolve','Display','iter','TolX',1e-6,'TolFun',1e-6);

w     =     ss(1);
L     =     ss(2);
C     =     ss(3);
K     =     ss(4);
Y     =     ss(5);



F(1) =  (1 - alphha)*kns^alphha - w;

F(2) = ((1/thetaa)*(1/cns)*w)^(1/(1+phhi)) - L;

F(3) = cns*L - C;

F(4) =  kns*L - K;

F(5) = K^(alphha)*L^(1-alphha) - Y;


end