A229502 Expansion of q * f(-q) * f(-q^4) * f(-q^16) * f(q, -q^3) in powers of q where f() is a Ramanujan theta function.
1, 0, -2, -2, 0, 2, 2, 4, -1, -2, 0, 0, -2, -4, 2, -4, -2, 2, 6, 4, 0, 2, -6, -4, 3, 2, 0, 4, 6, 0, -8, 0, -2, 0, -4, -2, 0, -6, 2, -4, 0, 4, 0, -4, 2, 12, 8, 8, 3, -6, 4, 0, 0, -8, 2, 0, -6, -6, 0, -4, -18, 0, 2, 8, 2, 0, -10, 4, 0, 4, 10, 0, 4, 6, 0, 0, 4
Offset: 1
Keywords
Examples
G.f. = q - 2*q^3 - 2*q^4 + 2*q^6 + 2*q^7 + 4*q^8 - q^9 - 2*q^10 - 2*q^13 + ...
Links
- G. C. Greubel, Table of n, a(n) for n = 1..2500 (terms 1..77 from Michael Somos)
- Michael Somos, Introduction to Ramanujan theta functions
- Eric Weisstein's World of Mathematics, Ramanujan Theta Functions
Programs
-
Magma
Basis( CuspForms( Gamma1(16), 2), 78)[1];
-
Mathematica
a[ n_] := SeriesCoefficient[ q QPochhammer[ q^3, -q^4] QPochhammer[ -q, -q^4] QPochhammer[ -q^4] QPochhammer[ q] QPochhammer[ q^4] QPochhammer[ q^16], {q, 0, n}]; -
PARI
{a(n) = my(A, m); if( n<1, 0, n--; A = x * O(x^n); polcoeff( eta(x + A) * eta(x^4 + A) * eta(x^16 + A) * sum( k=0, n, if( issquare( 8*k + 1, &m), (-1)^((m\2 + 2) \ 4) * x^k, 0), A), n))}; -
Sage
CuspForms( Gamma1(16), 2, prec=78).0;
Comments