The introduction of epsilon factorized differential equations for Feynman Integrals in dimensional regularization, gave a massive boost to their computation in polylogaritmic cases. In this presentation I will demonstrate a method to obtain such a form 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. This method is a generalization of a similar method known to work for polylogarithmic Feynman integrals. I will demonstrate the method explicitly for a number of Feynman integral families with an elliptic highest sector. The presentation will largely be based on the recent paper arXiv:2110.07968.