function n_ss=cc_steady_state_helper(eta,beta,delta,theta)
options=optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);
n_ss=fsolve(@(n_ss)(eta*(1-beta*(1-delta))/(beta*(1-delta))+(theta-1)*(n_ss-1)/(theta*n_ss-theta+1)*n_ss^((2-theta)/(theta-1))),1,options)


