function H=benchmarkcb_steady_state_helper(tau_h,tau_k,iota,theta,chi_s,sigma_s,sigma_c,nsigma_s,xi,delta,beta)
options=optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);
H=fsolve(@(H) 1- (1-tau_h)*(1-theta)*iota/((1-chi_s)*(1-iota)*((1-H)^(-1/sigma_s))*((chi_s*((1/(xi*((((1/beta)+delta-1)/((1-tau_k)*theta))^(1-theta))))^nsigma_s)+(1-chi_s)*((1-H)^nsigma_s))^((sigma_c-sigma_s)/(sigma_c*sigma_s-sigma_c))))+delta*((((1/beta)+delta-1)/((1-tau_k)*theta))^(1/theta)),0.2);