|
%I #24 Mar 12 2021 22:24:46
%S 1,-3,0,8,-9,0,17,-27,0,46,-57,0,98,-126,0,198,-243,0,371,-465,0,692,
%T -828,0,1205,-1458,0,2082,-2463,0,3463,-4104,0,5678,-6642,0,9085,
%U -10623,0,14370,-16632,0,22273,-25758,0,34178,-39246,0,51674,-59220,0,77362
%N Expansion of q^(1/4) * (eta(q) / eta(q^3))^3 in powers of q.
%C Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
%H Seiichi Manyama, <a href="/A199659/b199659.txt">Table of n, a(n) for n = 0..10000</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 Expansion of (f(-x) / f(-x^3))^3 in powers of x where f() is a Ramanujan theta function.
%F Euler transform of period 3 sequence [ -3, -3, 0, ...].
%F G.f. is a period 1 Fourier series which satisfies f(-1 / (48 t)) = 3^(3/2) / f(t) where q = exp(2 Pi i t).
%F G.f.: (Product_{k>0} (1 - x^k) / (1 - x^(3*k)))^3.
%F Convolution cube of A137569. Convolution square is A007262. Convolution fourth power is A030182.
%F a(3*n + 2) = 0.
%F a(0) = 1, a(n) = -(3/n)*Sum_{k=1..n} A046913(k)*a(n-k) for n > 0. - _Seiichi Manyama_, Apr 28 2017
%e 1 - 3*x + 8*x^3 - 9*x^4 + 17*x^6 - 27*x^7 + 46*x^9 - 57*x^10 + 98*x^12 + ...
%e 1/q - 3*q^3 + 8*q^11 - 9*q^15 + 17*q^23 - 27*q^27 + 46*q^35 - 57*q^39 + ...
%t QP = QPochhammer; s = (QP[q]/QP[q^3])^3 + O[q]^60; CoefficientList[s, q] (* _Jean-François Alcover_, Nov 25 2015 *)
%o (PARI) {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x + A) / eta(x^3 + A))^3, n))}
%Y Cf. A007262, A030182, A137569.
%K sign
%O 0,2
%A _Michael Somos_, Nov 08 2011
|
