function eps=allowtax_steady_state_helper(alpha,beta,delta,Thetabar,Pibar,Hbar,rb)
options=optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);
eps=fsolve(@(eps)beta*((alpha/(alpha*Pibar-rb*Hbar))*(delta*Pibar*eps/(1-(eps/2)^2)-beta*(eps/2+(1-(eps/2)^2)/eps)*Thetabar*Hbar^2)*(eps/2+(1-(eps/2)^2)/eps)+1-delta+((eps/2+(1-(eps/2)^2)/eps)-1)*(eps/2+(1-(eps/2)^2)/eps)*beta*Thetabar*Hbar^2/Pibar)-1,1.5,options);