OFFSET
0,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..500
FORMULA
a(n) = [x^n] Product_{k>=1} ((1 + x^k)/(1 - x^(4*k)))^n.
a(n) ~ c * d^n / sqrt(n), where d = 4.62579056836776492108784045382518984897... (see A192540) and c = 0.255113338880004277664416308115912337... - Vaclav Kotesovec, Dec 05 2017
MATHEMATICA
Table[SeriesCoefficient[Product[((1 + x^(2 k - 1))/(1 - x^(2 k)))^n, {k, 1, n}], {x, 0, n}], {n, 0, 26}]
Table[SeriesCoefficient[Product[((1 + x^k)/(1 - x^(4 k)))^n, {k, 1, n}], {x, 0, n}], {n, 0, 26}]
Table[SeriesCoefficient[(2 (-x)^(1/8)/EllipticTheta[2, 0, Sqrt[-x]])^n, {x, 0, n}], {n, 0, 26}]
Table[(-1)^n * 2^n * SeriesCoefficient[1/(QPochhammer[-1, x]*QPochhammer[x^2])^n, {x, 0, n}], {n, 0, 30}] (* Vaclav Kotesovec, Oct 07 2020 *)
(* Calculation of constants {d, c}: *) Chop[{1/r, 4/Sqrt[Pi*(77/2 - 4*s*(-r*s)^(7/8) * Derivative[0, 0, 2][EllipticTheta][2, 0, Sqrt[-r*s]])]} /. FindRoot[{s == (2*(-r*s)^(1/8))/EllipticTheta[2, 0, Sqrt[-r*s]], 7*I*r + 2*(-r*s)^(7/8)*Sqrt[r*s] * Derivative[0, 0, 1][EllipticTheta][2, 0, Sqrt[-r*s]] == 0}, {r, 1/5}, {s, 2}, WorkingPrecision -> 70]] (* Vaclav Kotesovec, Jan 17 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Dec 03 2017
STATUS
approved