function nstar = rbc_steadystate_helper(nstar0,wstar,t_ratio,cstar_nstar,sig,thet,snakey)

options = optimset('Display', 'Final', 'TolX', 1e-10, 'TolFun', 1e-10);

nstar = fsolve(@(nstar) wstar*t_ratio*cstar_nstar^-sig - thet*((1-nstar)^-snakey)*(nstar)^sig, nstar0,options);
end