OFFSET
0,4
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
J. H. Conway and N. J. A. Sloane, Sphere Packings, Lattices and Groups, Springer-Verlag, p. 190, Equation (47).
FORMULA
Expansion of ( E_4(q) * 2 * (E_4(q^2) - E_4(q^4)) + E_4(q^2) * (32 * E_4(q^4) - 17 * E_4(q^2)) ) / 15 in powers of q. - Michael Somos, Nov 29 2007
EXAMPLE
1 + 7680*q^3 + 4320*q^4 + 276480*q^5 + 61440*q^6 + 2903040*q^7 + ...
MATHEMATICA
QP = QPochhammer; a[n_] := Module[{A, A1, A2, A4}, A = x*O[x]^n; A1 = QP[x+ A]^8; A2 = QP[x^2+A]^8; A4 = QP[x^4+A]^8; SeriesCoefficient[(A1*(A2^6 + x^2*32*A2^3*A4^3 + x^4*4096*A4^6) + x^3*3840*A4^4*(A1^2*A4 + A2^3)) / (A1*A2^2*A4^2), n]]; a[0] = 1; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Nov 30 2015, adapted from PARI *)
PROG
(PARI) {a(n) = local(A, A1, A2, A4); if( n<0, 0, A = x * O(x^n); A1 = eta(x + A)^8; A2 = eta(x^2 + A)^8; A4 = eta(x^4 + A)^8; polcoeff( ( A1 * (A2^6 + x^2 * 32 * A2^3 * A4^3 + x^4 * 4096 * A4^6) + x^3 * 3840 * A4^4 * ( A1^2 * A4 + A2^3 ) ) / (A1 * A2^2 * A4^2 ), n))} /* Michael Somos, Nov 29 2007 */
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved