OFFSET
0,2
COMMENTS
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..2000
Michael Somos, Introduction to Ramanujan theta functions
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions
FORMULA
Expansion of psi(-q) * psi(q^9) / (psi(q) * psi(-q^9)) in powers of q where psi() is a Ramanujan theta function.
Expansion of eta(q)^2 * eta(q^4) * eta(q^18)^3 / (eta(q^2)^3 * eta(q^9)^2 * eta(q^36)) in powers of q.
Euler transform of period 36 sequence [ -2, 1, -2, 0, -2, 1, -2, 0, 0, 1, -2, 0, -2, 1, -2, 0, -2, 0, -2, 0, -2, 1, -2, 0, -2, 1, 0, 0, -2, 1, -2, 0, -2, 1, -2, 0, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (36 t)) = g(t) where q = exp(2 Pi i t) and g() is the g.f. of A128143.
a(3*n) = A261320(n). a(3*n + 1) = -2 * A261325(n). a(3*n + 2) = 2 * A260057(n). - Michael Somos, Nov 08 2015
a(n) ~ (-1)^n * exp(2*Pi*sqrt(n)/3) / (2*sqrt(3)*n^(3/4)). - Vaclav Kotesovec, Nov 16 2017
EXAMPLE
G.f. = 1 - 2*x + 2*x^2 - 4*x^3 + 6*x^4 - 8*x^5 + 12*x^6 - 16*x^7 + 22*x^8 + ...
MATHEMATICA
a[ n_] := SeriesCoefficient[ QPochhammer[ q, q^2] QPochhammer[ q, -q] QPochhammer[ -q^9, q^18] QPochhammer[ -q^9, q^9], {q, 0, n}];
PROG
(PARI) {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^2 * eta(x^4 + A) * eta(x^18 + A)^3 / (eta(x^2 + A)^3 * eta(x^9 + A)^2 * eta(x^36 + A)), n))};
CROSSREFS
KEYWORD
sign
AUTHOR
Michael Somos, Aug 13 2015
STATUS
approved