OFFSET
-1,2
LINKS
G. C. Greubel, Table of n, a(n) for n = -1..1000
FORMULA
Expansion of (eta(q) * eta(q^5) / (eta(q^3) * eta(q^15)))^2 in powers of q.
Euler transform of period 15 sequence [ -2, -2, 0, -2, -4, 0, -2, -2, 0, -4, -2, 0, -2, -2, 0, ...].
G.f. A(x) satisfies 0 = f(A(x), A(x^2)) where f(u, v) = u^3 + v^3 - 4*u*v * (u + v) - u^2*v^2 - 9*u*v.
a(n) = A058589(n) unless n=0.
EXAMPLE
G.f. = 1/q - 2 - q + 4*q^2 - 3*q^3 - 2*q^4 + 11*q^5 - 6*q^6 - 11*q^7 + ...
MATHEMATICA
a[ n_] := SeriesCoefficient[ q^(-1) (QPochhammer[ q] QPochhammer[ q^5] / (QPochhammer[q^3] QPochhammer[q^15]))^2, {q, 0, n}] (* Michael Somos, Oct 23 2013 *)
PROG
(PARI) {a(n) = local(A); if( n<-1, 0, n++; A = x*O(x^n); polcoeff( (eta(x + A) * eta(x^5 + A) / (eta(x^3 + A) * eta(x^15 + A)))^2, n))}
(PARI) {a(n) = local(A, u, v); if( n<-1, 0, A = 1/x; for(k=0, n, u = A + x * O(x^k); v = subst(u, x, x^2); A += x^k * polcoeff(u^3 + v^3 - 4*u*v*(u + v) - u^2*v^2 - 9*u*v, k-5) / 2); polcoeff(A, n))}
CROSSREFS
KEYWORD
sign
AUTHOR
Michael Somos, Mar 17 2004
EXTENSIONS
Edited for better readability and more coherence. - Michael Somos, Oct 23 2013
STATUS
approved