ep-Finite Basis of Master Integrals for the Integration-By-Parts Method
  It is shown that for every problem within dimensional regularization, using the Integration-By-Parts method, one is able to construct a set of master integrals such that each corresponding coefficient function is finite in the limit of dimension equal to four. We argue that the use of such a basis simplifies and stabilizes the numerical evaluation of the master integrals. As an example we explicitly construct the ep-finite basis for the set of all QED-like four-loop massive tadpoles. Using a semi-numerical approach based on Pade approximations we evaluate analytically the divergent and numerically the finite part of this set of master integrals.   The calculations confirm the recent results of Schr\"oder and Vuorinen. All the contributions found there by fitting the high precision numerical results have been confirmed by direct analytical calculation without using any numerical input. 
K.G. Chetyrkin, M. Faisst, C. Sturm and M.Tentyukov