Automatic Computation Of Three-Loop Two-Point Functions In Large   Momentum Expansion
  We discuss the calculation of two-point three-loop functions with an arbitrary number of massive propagators and one large external momentum. The relevant subdiagrams are generated automatically. The resulting massless two-point integrals and massive tadpoles are transformed on-line to FORM-expressions ready to be used by existing FORM packages which calculate them analytically. As an example we compute the quartic mass corrections to the photon polarization function. 
K.G. Chetyrkin, R. Harlander, J.H. Kuehn and M. Steinhauser