We show that all Eichler integrals, and more generally all "generalized second order modular forms" can be expressed as linear combinations of corresponding generalized second order Eisenstein series with coefficients in classical modular forms. We determine the Fourier series expansions of generalized second order Eisenstein series in level one, and provide tail estimates via convexity bounds for additively twisted L-functions. As an application, we illustrate a bootstrapping procedure that yields numerical evaluations of, for instance, Eichler integrals from merely the associated cocycle. The proof of our main results rests on a filtration argument that is largely rooted in previous work on vector-valued modular forms, which we here formulate in classical terms.second order modular forms period polynomials Rankin-Selberg convolution MSC Primary: 11F11 MSC Secondary: 11F30, 11F75 I N recent years, there has been a surge of interest in the study of iterated Eichler-Shimura integrals, see for example [2][3][4]7]. This has in part been motivated by the fact that they are closely related to the string theoretic notion of holomorphic graph functions and modular graph functions. Indeed, the A-cycle graph functions and B -cycle graph functions introduced by Broedel-Schlotterer-Zerbini [1] can be expressed in terms of iterated Eisenstein integrals. Diamantis [7] showed that iterated Eichler-Shimura integrals are also closely related to higher order modular forms, first introduced by Goldfeld [12]. In a similar vein Mertens-Raum [16] showed that both higher order modular forms and iterated Eichler-Shimura integrals can be reinterpreted as components of vector-valued modular forms of higher depth arithmetic types. In this paper, we explore the problem of efficiently computing Eichler-Shimura integrals by means of the aforementioned theoretical frameworks, focusing on the special case of Eichler integrals of cusp forms of level one. Our hope is that this approach can eventually be generalized to provide alternative means of computing point evaluations and Fourier series expansions of modular graph functions, for example considered by D'Hoker-Duke in [6].We offer two theorems, that demonstrate the viability of our approach. They also prepare the route to a more general framework (see Remark 4.4). In Section 2, we introduce the notion of generalized second order modular forms, and show that Eichler integrals of cusp forms are examples. Generalized second order modular forms are holomorphic functions on H subject to an appropriate growth condition and whose modular deficits are cocycles taking values in polynomials whose coefficients are modular forms. Their name is motivated by the fact that they generalize second order modular forms. Indeed, the modular deficits of second order modular forms are cocycles with values in modular forms.