product of the bethe lattice dos, fermi distribution an weighted
       by energy

def bethe_fermi_ene(energy, quasipart, shift, hopping, beta):
    """product of the bethe lattice dos, fermi distribution an weighted
       by energy"""
    return energy * bethe_fermi(energy, quasipart, shift, hopping, beta)