OFFSET
-1,4
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, Commun. Algebra 22, No. 13, 5175-5193 (1994).
FORMULA
Expansion of 2 + ( eta(q) / eta(q^13) )^2 in powers of q. - Michael Somos, Jul 05 2012
G.f.: 2 + (1/x) * Product_{k>0} ((1 - x^k) / (1 - x^(13*k)))^2. - Michael Somos, Jul 05 2012
a(n) = A133099(n) unless n=0. - Michael Somos, Jul 05 2012
EXAMPLE
T13B = 1/q - q + 2*q^2 + q^3 + 2*q^4 - 2*q^5 - 2*q^7 - 2*q^8 + q^9 + 4*q^12 - ...
MATHEMATICA
QP = QPochhammer; s = 2*q+(QP[q]/QP[q^13])^2+O[q]^70; CoefficientList[s, q] (* Jean-François Alcover, Nov 13 2015, adapted from PARI *)
eta[q_]:= q^(1/24)*QPochhammer[q]; a[n_]:= SeriesCoefficient[2 + (eta[q] /eta[q^13])^2, {q, 0, n}]; Table[a[n], {n, -1, 50}] (* G. C. Greubel, Feb 18 2018 *)
PROG
(PARI) {a(n) = local(A); if( n<-1, 0, n++; A = x*O(x^n); polcoeff( (eta(x + A) / eta(x^13 + A))^2 + 2*x, n))} /* Michael Somos, Jul 05 2012 */
CROSSREFS
KEYWORD
sign
AUTHOR
N. J. A. Sloane, Nov 27 2000
STATUS
approved