|   | 1 |  initial version  | 
PARI/GP has very advanced functionalities to deal with modular forms (any weight, \Gamma_0(N) and any character), see the documentation at https://pari.math.u-bordeaux.fr/dochtml/html/Modular_forms.html
For evaluating Eisenstein series, you can do
sage: mf = pari.mfinit([1,4],3)   # Eisenstein series weight 4 for SL(2,Z)
sage: f = pari.mfbasis(mf)[0]     # dimension is one so this is the unique Eisenstein series
sage: pari.mfeval(mf, f, I)
0.00606567871778629
|   | 2 |  No.2 Revision  | 
PARI/GP has very advanced functionalities to deal with modular forms (any weight, (\Gamma_0(N) , any weight and any character), see the documentation at https://pari.math.u-bordeaux.fr/dochtml/html/Modular_forms.html
For evaluating Eisenstein series, you can do
sage: mf = pari.mfinit([1,4],3)   # Eisenstein series weight 4 for SL(2,Z)
sage: f = pari.mfbasis(mf)[0]     # dimension is one so this is the unique Eisenstein series
sage: pari.mfeval(mf, f, I)
0.00606567871778629
 Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.
 
                
                Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.