Single scale Feynman integrals in quantum field theories obey difference or differential equations with respect to their discrete parameter $N$ or continuous parameter $x$. The analysis of these equations reveals to which order they factorize, which can be different in both cases. The simplest systems are the ones which factorize to first order. For them complete solution algorithms exist. The next interesting level is formed by those cases in which also irreducible second order systems emerge. We give a survey on the latter case. The solutions can be obtained as general $_2F_1$ solutions. The corresponding solutions of the associated inhomogeneous differential equations form so-called iterative non-iterative integrals. There are known conditions under which one may represent the solutions by complete elliptic integrals. In this case one may find representations in terms of meromorphic modular functions, out of which special cases allow representations in the framework of elliptic polylogarithms with generalized parameters. These are in general weighted by a power of $1/\eta(\tau)$, where $\eta(\tau)$ is Dedekind's $\eta$-function. Single scale elliptic solutions emerge in the $\rho$-parameter, which we use as an illustrative example. They also occur in the 3-loop QCD corrections to massive operator matrix elements and the massive 3-loop form factors.