function[ys, check] = sims_steadystate(ys,exe);

global M_

beta=0.99;
alpha =1/3;
sigmae=0.01;
rho= 1;
theta = 1.40;
nu = 1.15;
phi_p = 37.7;
phi_w = 280;
rho_x = 0.9;

pi = 1;
omega = 1;
x = 1;
h = (((1-alpha)/(theta*nu))/(x^(rho-1)))^((alpha+rho*(1-alpha))^(-1));
y = x*h^(1-alpha);
c = y;
w_real = ((1-alpha)*x*h^(-alpha))/nu;
check=0;

ys = [log(x); log(c); log(y); log(w_real); log(h); pi;omega];