% Steady state of Chapter_3_DSGE
function [ys,check] = Chapter_3_DSGE_steadystate(ys,exe)
  global M_
  
  alpha = M_.params(1); 
  beta = M_.params(2); 
  gamma = M_.params(3); 
  mst = M_.params(4); 
  rho = M_.params(5); 
  psi = M_.params(6); 
  delta = M_.params(7); 

  check = 0;
  
  dA = exp(gamma);
  gst = 1/dA;
  m = mst;
  
  khst = ( (1-gst*beta*(1-delta)) / (alpha*gst^alpha*beta) )^(1/(alpha-1));
  xist = ( ((khst*gst)^alpha - (1-gst*(1-delta))*khst)/mst )^(-1);
  nust = psi*mst^2/( (1-alpha)*(1-psi)*beta*gst^alpha*khst^alpha );
  n  = xist/(nust+xist);
  P  = xist + nust;
  k  = khst*n;
  l  = psi*mst*n/( (1-psi)*(1-n) );
  c  = mst/P;
  d  = l - mst + 1;
  y  = k^alpha*n^(1-alpha)*gst^alpha;
  R  = mst/beta;
  W  = l/n;
  ist  = y-c;
  q  = 1 - d;
  e = 1;
  
  P_obs = 1;
  Y_obs = 1;
  
  ys =[
m     
P     
c     
e     
W     
R     
k     
d     
n     
l     
Y_obs 
P_obs 
y     
dA          ];