Elliptic polylogarithms and Feynman parameter integrals

In this paper we study the calculation of multiloop Feynman integrals that cannot be expressed in terms of multiple polylogarithms. We show in detail how certain types of two- and three-point functions at two loops, which appear in the calculation of higher order corrections in QED, QCD and in the electroweak theory (EW), can naturally be expressed in terms of a recently introduced elliptic generalisation of multiple polylogarithms by direct integration over their Feynman parameter representation. Moreover, we show that in all examples that we considered a basis of pure Feynman integrals can be found.