function [ys check]=rbc_steadystate(junk,ys)
global beta gamma theta eta delta abar rho 
check=0;

  kappa = eta/beta - 1 + delta;

  lambda = eta - 1 + delta;

  hss = ((1-theta)/gamma)/(1-theta*lambda/kappa);

  yss = (abar^(1/(1-theta)))*((theta/kappa)^(theta/(1-theta)))*hss;

  kss = (theta/kappa)*yss;

  iss = (theta*lambda/kappa)*yss;

  css = yss - iss;

  ass=abar;

a=ass;
c=css;
h=hss;
i=iss;
k=kss;
y=yss;
yobs=y; 
cobs=c;
hobs=h;
uy=0 ;
uc=0;
uh=0;
ys=[a;c;cobs;h;hobs;i;k;uc;uh;uy;y;yobs];