Return the fermi energy that correspond to the given particle quantity
       in a semicircular density of states of a bethe lattice in a multi
       orbital case that can be non-degenerate

def bethe_findfill_zeroT(particles, orbital_e, hopping):
    """Return the fermi energy that correspond to the given particle quantity
       in a semicircular density of states of a bethe lattice in a multi
       orbital case that can be non-degenerate"""

    assert 0. <= particles <= len(orbital_e)

    zero = lambda e: np.sum([bethe_filling_zeroT(e-e_m, t) \
                        for t, e_m in zip(hopping, orbital_e)]) - particles
    return fsolve(zero, 0)