function u = u_steadystate_helper(u0,ddelta,nnu,phi_k,ggamma,bbeta_2)
options = optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);
u = fsolve(@(u) (bbeta_2/ggamma)*( ddelta*nnu*(u^nnu)*(1-phi_k*((((1+(ddelta*u^nnu +ggamma-1)*phi_k)-( (1+(ddelta*u^nnu +ggamma-1)*phi_k)^2 -2*phi_k*((ddelta*u^nnu +ggamma-1)+((ddelta*u^nnu +ggamma-1)^2)*phi_k/2))^(1/2))/phi_k)-(ddelta*u^nnu +ggamma-1))) +(phi_k/2)*((((1+(ddelta*u^nnu +ggamma-1)*phi_k)-( (1+(ddelta*u^nnu +ggamma-1)*phi_k)^2 -2*phi_k*((ddelta*u^nnu +ggamma-1)+((ddelta*u^nnu +ggamma-1)^2)*phi_k/2))^(1/2))/phi_k)^2 - (ddelta*u^nnu +ggamma-1)^2)+1-ddelta*u^nnu) - 1, u0, options);









