login
Expansion of f(-q^2, -q^5)^3 / (f(-q^1, -q^6) * f(-q^3, -q^4)^2) in powers of q where f(, ) is Ramanujan's general theta function.
4

%I #20 Mar 12 2021 22:24:48

%S 1,1,-2,0,5,-4,-7,12,4,-22,7,29,-26,-28,52,14,-82,21,106,-85,-105,175,

%T 53,-268,70,326,-264,-301,505,142,-742,189,885,-698,-805,1323,374,

%U -1906,483,2205,-1732,-1946,3185,884,-4486,1120,5119,-3972,-4473,7229,2004

%N Expansion of f(-q^2, -q^5)^3 / (f(-q^1, -q^6) * f(-q^3, -q^4)^2) in powers of q where f(, ) is Ramanujan's general theta function.

%C Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

%H Seiichi Manyama, <a href="/A262933/b262933.txt">Table of n, a(n) for n = 0..10000</a> (terms 0..1000 from G. C. Greubel)

%H Michael Somos, <a href="/A010815/a010815.txt">Introduction to Ramanujan theta functions</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/RamanujanThetaFunctions.html">Ramanujan Theta Functions</a>

%F Euler transform of period 7 sequence [ 1, -3, 2, 2, -3, 1, 0, ...].

%F G.f.: T(q)/(T(q)-1), where T(q) = 1/q + 3 + 4*q + ... (cf. A108481). - _Seiichi Manyama_, Oct 11 2018

%e G.f. = 1 + q - 2*q^2 + 5*q^4 - 4*q^5 - 7*q^6 + 12*q^7 + 4*q^8 - 22*q^9 + ...

%t a[ n_] := SeriesCoefficient[ Product[ (1 - x^k)^{-1, 3, -2, -2, 3, -1, 0}[[Mod[k, 7, 1]]], {k, n}], {x, 0, n}];

%t (* alternative program *)

%t QP:= QPochhammer; a[n_]:= SeriesCoefficient[(QP[q^2,q^7]*QP[q^5,q^7])^3/ (QP[q,q^7]*QP[q^6,q^7]*QP[q^3,q^7]^2*QP[q^4,q^7]^2), {q, 0, n}];

%t Table[a[n], {n, 0, 50}] (* _G. C. Greubel_, Mar 18 2018 *)

%o (PARI) {a(n) = if( n<0, 0, polcoeff( prod(k=1, n, (1 - x^k + x * O(x^n))^[0, -1, 3, -2, -2, 3, -1][k%7 + 1]), n))};

%Y Cf. A108481.

%K sign

%O 0,3

%A _Michael Somos_, Oct 04 2015