function eppsilon_d = eps_steady_state_helper(eppsilon_d0, numb_banks, eppsilon_d_par, mmu_d, betta, thetta, gaamma_cb, gaamma_d, gaamma_m, b, i_cb)

%options = optimset('Display','Final','TolX',le-10,'TolFun', le-10);

eppsilon_d = fsolve( @(eppsilon_d) eppsilon_d_par*(numb_banks-1)/numb_banks+thetta/numb_banks-(gaamma_d/numb_banks)*(thetta-1/(b-1))*(((1/betta-mmu_d)*eppsilon_d/(eppsilon_d+1))/((gaamma_m+gaamma_d*((1/betta-mmu_d)*eppsilon_d/(eppsilon_d+1))^(thetta+1)+gaamma_cb*(1+i_cb)^(thetta+1))^(1/(thetta+1))))^(thetta+1)-eppsilon_d, eppsilon_d0 );
