OFFSET
-1,2
LINKS
G. C. Greubel, Table of n, a(n) for n = -1..1000
D. Ford, J. McKay and S. P. Norton, More on replicable functions, Comm. Algebra 22, No. 13, 5175-5193 (1994).
FORMULA
Expansion of (T19A + 4)^(1/2), where T19A = A058549, in powers of q. - G. C. Greubel, Jun 24 2018
a(n) ~ exp(2*Pi*sqrt(2*n/19)) / (2^(3/4) * 19^(1/4) * n^(3/4)). - Vaclav Kotesovec, Jun 27 2018
EXAMPLE
T38a = 1/q + 2*q + q^3 + 3*q^5 + 4*q^7 + 7*q^9 + 8*q^11 + 13*q^13 + 15*q^15 + ...
MATHEMATICA
QP := QPochhammer; nmax = 100; G[q_] := 1/(QP[q, q^5]*QP[q^4, q^5]); H[q_] := 1/(QP[q^2, q^5]*QP[q^3, q^5]); T19A := -3 + (G[q]*G[q^19] + q^4*H[q]*H[q^19])^3/q; a:= CoefficientList[Series[(q*(T19A + 4) + O[q]^nmax)^(1/2), {q, 0, 60}], q]; Table[a[[n]], {n, 1, 50}] (* G. C. Greubel, Jun 24 2018 *)
PROG
(PARI) q='q+O('q^40); A = (eta(q)*eta(q^19)/(eta(q^2)*eta(q^38)))^2; B = - (eta(-q)*eta(-q^19)/(eta(q^2)*eta(q^38)))^2; T19A = (q^4/(A*B) - (A + B)/(4*q))^3; v=Vec((q^2 + T19A)^(1/2)); vector(#v\2, n, v[2*n-1]) \\ G. C. Greubel, Jun 25 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Nov 27 2000
EXTENSIONS
Terms a(12) onward added by G. C. Greubel, Jun 24 2018
STATUS
approved