function  [Cbar,Kbar]=ckss(alph,bet,gam,delt)
%function  [Cbar,Kbar]=ckss(alph,bet,gam,delt)
Kbar=(alph*bet/(1-bet*delt))^(1/(1-alph));
%Cbar=Kbar*((1-bet*delt*(1-alph))/(alph*bet)-1);
Cbar=Kbar^alph-(1-delt)*Kbar;