  In this paper we develop and demonstrate a method to obtain epsilon
factorized differential equations for elliptic Feynman integrals. This method
works by choosing an integral basis with the property that the period matrix
obtained by integrating the basis over a complete set of integration cycles is
diagonal. The method is a generalization of a similar method known to work for
polylogarithmic Feynman integrals. We demonstrate the method explicitly for a
number of Feynman integral families with an elliptic highest sector.

High Energy Physics - Theory