The $\varepsilon$-form of a system of differential equations for Feynman integrals has led to tremendeous progress in our abilities to compute Feynman integrals, as long as they fall into the class of multiple polylogarithms. It is therefore of current interest, if these methods extend beyond the case of multiple polylogarithms. In this talk I discuss Feynman integrals, which are associated to elliptic curves and their differential equations. I show for non-trivial examples how the system of differential equations can be brought into an $\varepsilon$-form. Single-scale and multi-scale cases are discussed.