function SS = steady(SS, delk, delj, c, b, alpha, rho, s, xi, beta)
SS = real(SS);
k = SS(1);
x = SS(2);

j = 0; %This is normalized to be 1;
p = log(delj*exp(j));
m = log(beta);
z=0; % no idio. shocks in steadystate;

CES = log(exp(k)-(1-delk)*exp(k));
i = log((((delk*exp(k)/c)^rho-(1-b)*((1-s)*exp(j))^rho)/b)^(1/rho));
CES1 = log(c*(1/rho)*(b*exp(i)^rho+(1-b)*((1-s)*exp(j))^rho)^(1/rho-1)*b*rho*exp(i)^(rho-1));
CES2 = log(c*(1/rho)*(b*exp(i)^rho+(1-b)*((1-s)*exp(j))^rho)^(1/rho-1)*(1-b)*rho*((1-s)*exp(j))^(rho-1)*(1-s));
qk = log(1/exp(CES1) );
y = log(exp(x+z)*(exp(k))^(xi*alpha)*(s*exp(j))^(xi*(1-alpha)));
SS(1) =exp(qk)-(exp(m)*(alpha*xi*exp(y)/exp(k)+(1-delk)*exp(qk)) );
SS(2) =1-exp(qk)*exp(CES2)-(exp(m)*((1-alpha)*xi*exp(y)/(exp(j))+1-delj));
SS = real(SS);
end