function [rss,wss,css,mpss,kss,hss,yss]=findss(betax,deltax,theta,barg,B, D)
rss = 1/betax -1 + deltax;
wss = (1-theta)*(rss/theta)^(theta/(theta-1));
css = -wss/B;
mpss = D*barg*css/(barg-betax);
kss = css/((rss*(1-theta)/(wss*theta))^(1-theta)-deltax);
hss = rss*(1-theta)/(wss*theta)*kss;
yss = kss*(rss*(1-theta)/(wss*theta))^(1-theta);
end