<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 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)