OFFSET
0,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
D. Ford, J. McKay and S. P. Norton, More on replicable functions, Comm. Algebra 22, No. 13, 5175-5193 (1994).
FORMULA
Expansion of A - q/A, where A = q^(1/2)*(eta(q^4)*eta(q^5)/( eta(q)* eta(q^20))), in powers of q. - G. C. Greubel, Jun 26 2018
a(n) ~ exp(Pi*sqrt(2*n/5)) / (2^(5/4) * 5^(1/4) * n^(3/4)). - Vaclav Kotesovec, Jun 27 2018
Empirical: Sum_{n>=0} a(n) / exp(n*Pi) = (8/25) * exp(-Pi / 2) * 2^(2/5) * 5^(3/4) * Gamma(9/10)^2 * Gamma(7/10)^2 * (5+5^(1/2))^2 * (1/4*5^(1/2)-1/4)^2 * (1/4*5^(1/2)+1/4)^2 / Gamma(4/5)^4 = A388522. - Simon Plouffe, Sep 17 2025
EXAMPLE
T40a = 1/q +3*q^3 +4*q^5 +4*q^7 +4*q^9 +7*q^11 +12*q^13 +...
MATHEMATICA
eta[q_] := q^(1/24)*QPochhammer[q]; A:= q^(1/2)*(eta[q^4]*eta[q^5]/( eta[q]*eta[q^20])); a := CoefficientList[Series[A - q/A, {q, 0, 60}], q]; Table[a[[n]], {n, 1, 50}] (* G. C. Greubel, Jun 26 2018 *)
PROG
(PARI) q='q+O('q^50); A = eta(q^4)*eta(q^5)/(eta(q)*eta(q^20)); Vec(A - q/A) \\ G. C. Greubel, Jun 26 2018
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Michael Somos, Aug 28 2005
STATUS
approved
