function [ys,check] = SOE_steadystate(ys,exe);
 global M_

 
sigma = M_.params(1);
delta = M_.params(2);
rstar = M_.params(3);
alpha = M_.params(4);
omega = M_.params(5);
dbar = M_.params(6);
pssi = M_.params(7);
phi = M_.params(8);
rho = M_.params(9);
etatilde = M_.params(10);
sigmae = M_.params(11);
beta = M_.params(12);


a_t = 1;%technological factor
r_t = rstar; %interest rate
d_t = dbar; %debt
kh = ((1/beta-(1-delta))/alpha)^(1/(alpha-1)); %k/h
h_t = ((1-alpha)*kh^alpha)^(1/(omega -1)); 
k_t = kh * h_t; %capital
y_t = kh^alpha * h_t; %output
i_t = delta * k_t;%investment
c_t = y_t-delta*k_t-rstar*dbar;
tb_t = y_t - i_t - c_t;%trade balance
tby_t = tb_t/y_t;
ca_t = -r_t*d_t+tb_t;
la_t = ((c_t - h_t^omega/omega))^(-sigma); %marginal utility of wealth
cay_t = ca_t/y_t;






check = 0;
ys = [log(c_t);
log(k_t);
log(h_t);
d_t;
log(y_t);
log(i_t);
tb_t;
tby_t;
ca_t;
cay_t;
r_t;
log(a_t);
log(la_t);
kh];