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.