function [ys,check] = hansen_stoch_exog_estim_steadystate(ys,exe)
  global M_
  
  alpha = M_.params(1); 
  beta = M_.params(2); 
  delta = M_.params(3); 
  gamma = M_.params(4); 
  eta = M_.params(5);
   

  check = 0;

    dmu = exp(eta);
    mst = (dmu)^(-1);
    a=0;
    c = ((1-alpha)/gamma)*exp(a)^(-1/(1-alpha))*(1/alpha*1/dmu^(1-alpha)*(dmu/beta+delta-1))^(alpha/(alpha-1));
    k = c*(mst-delta*mst+exp(a)*(c*gamma*(1-alpha)^(-1)*exp(a)^(-1))^((alpha-1)/alpha)*mst^(alpha)-1)^(-1);
    h = k*(c*gamma*(1-alpha)^(-1)*exp(a)^(-1))^(-1/(alpha));
    y = k^alpha*h^(1-alpha)*exp(-alpha*(eta))*exp(a);
    Y_obs = 1;
  
  ys =[
c
dmu
h     
k 
y
Y_obs
    ];

