Factorizing one-loop contributions to two-loop Bhabha scattering and   automatization of Feynman diagram calculations
  In higher order calculations a number of new technical problems arise: one needs diagrams in arbitrary dimension in order to obtain their needed $\epsilon$-expansion, zero Gram determinants appear, renormalization produces diagrams with `dots' on the lines, i.e. higher order powers of scalar propagators. All these problems cannot be accessed by the `standard' Passarino-Veltman approach: there is not available what is needed for higher loops. We demonstrate our method of how to solve these problems. 
J. Fleischer, O.V. Tarasov, T. Riemann and A. Werthenbach