|
%I #16 Mar 12 2021 22:24:46
%S 1,1,0,1,0,0,-1,-2,0,1,4,0,0,-1,0,-3,-8,0,4,14,0,1,-4,0,-6,-23,0,5,40,
%T 0,1,-10,0,-10,-60,0,11,98,0,4,-24,0,-19,-140,0,17,224,0,4,-54,0,-31,
%U -304,0,31,478,0,9,-112,0,-50,-627,0,46,968,0,11,-224
%N Expansion of q^(1/2) * eta(q^2)^2 * eta(q^6)^2 / (eta(q) * eta(q^9)^3) in powers of q.
%C Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
%H G. C. Greubel, <a href="/A182036/b182036.txt">Table of n, a(n) for n = 0..1000</a>
%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 18 sequence [1, -1, 1, -1, 1, -3, 1, -1, 4, -1, 1, -3, 1, -1, 1, -1, 1, 0, ...].
%F a(3*n) = A132179(n). a(3*n + 2) = 0.
%F Expansion of psi(x) * f(-x^6)^2 / f(-x^9)^3 in powers of x where psi(), f() are Ramanujan theta functions. - _Michael Somos_, Aug 10 2017
%e G.f. = 1 + x + x^3 - x^6 - 2*x^7 + x^9 + 4*x^10 - x^13 - 3*x^15 - 8*x^16 + ...
%e G.f. = 1/q + q + q^5 - q^11 - 2*q^13 + q^17 + 4*q^19 - q^25 - 3*q^29 - 8*q^31 + ...
%t eta[x_] := x^(1/24)*QPochhammer[x]; A182036[n_] := SeriesCoefficient[q^(1/2)*(eta[q^2]* eta[q^6])^2/(eta[q]*eta[q^9]^3), {q,0,n}]; Table[A182036[n], {n,0,50}] (* _G. C. Greubel_, Aug 09 2017 *)
%t a[ n_] := SeriesCoefficient[ 1/2 x^(-1/8) EllipticTheta[ 2, 0, x^(1/2)] QPochhammer[ x^6]^2 / QPochhammer[ x^9]^3, {x, 0, n}]; (* _Michael Somos_, Aug 10 2017 *)
%o (PARI) {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^2 * eta(x^6 + A)^2 / (eta(x + A) * eta(x^9 + A)^3), n))};
%Y Cf. A132179.
%K sign
%O 0,8
%A _Michael Somos_, Apr 07 2012
|
