  We present a set of scalar master integrals (MIs) needed for a complete
treatment of massive two-loop corrections to Bhabha scattering in QED,
including integrals with arbitrary fermionic loops. The status of analytical
solutions for the MIs is reviewed and examples of some methods to solve MIs
analytically are worked out in more detail. Analytical results for the pole
terms in epsilon of so far unknown box MIs with five internal lines are given.
