function a=example3_steady_state_helper(alpha,beta,eta,mg,nu,psi,Rf,theta)
options=optimset('Display','off','TolX',1e-10,'TolFun',1e-10);
a=fsolve(@(a)(1/Rf-beta)*(-(1+alpha)*((1/beta-Rf)/exp(mg))^(1/eta)+1+(Rf-1)*a/Rf)^(-theta)-psi*a^(-nu),0.2,options);
