n = numerical_eigenforms(1,24)
n.ap(2)
[8388609.0, -4016.3511717179363, 5096.35117171844]
n.ap(3)
[177148.0, 252.00000000001896]
eisenstein_series_qexp(24, 5, normalization='integral')
236364091 + 131040*q + 1099243323360*q^2 + 12336522153621120*q^3 + 9221121336284413920*q^4 + O(q^5)