OFFSET
0,2
COMMENTS
Notice the unusual pattern in the signs of the coefficients of the series e(q)/e(-q), which appears to consist of a group of either 7 or 8 consecutive positive coefficients followed by a group of either 7 or 8 consecutive negative coefficients. Empirically, sign(a(n)) - sign(a(n+15)) is frequently zero. Similar behavior is seen in the coefficients of the power series e(-q) (A392020), 1/e(-q) (A392021) and e(q)*e(-q) (A392022).
FORMULA
e(q)/e(-q) = 1 + 2*q + 2*q^2 + 4*q^3 + 6*q^4 + 6*q^5 + 8*q^6 + 2*q^7 - 6*q^8 - 26*q^9 - 62*q^10 - 96*q^11 - 156*q^12 - 150*q^13 - 106*q^14 + 150*q^15 + 698*q^16 + 1590*q^17 + 3142*q^18 + 4630*q^19 + 6234*q^20 + 5436*q^21 + 528*q^22 - ...
MAPLE
# uses 2 X 2 matrices to determine rational function approximations to e(q)/e(-q) - see Morier-Genoud and Ovsienko, Section 2.4
# calculate the partial quotients of the simple continued fraction expansion of e
c := proc(n) if n = 1 then 2 elif irem(n, 3) = 0 then (2/3)*n else 1 end if end proc:
# define the q-analog of the integer n
f(n, q) := simplify( (q^n - 1)/(q - 1) ):
with(LinearAlgebra):
Mat(n, q) := Matrix( [[f(c(2*n-1), q), q^(c(2*n-1))], [1, 0]] ).Matrix([[q*f(c(2*n), q), 1], [q^(c(2 n)), 0]] ):
Mat_prod := proc(n, q) option remember; if n = 1 then Mat(1, q) else simplify( Mat_prod(n-1, q).Mat(n, q) ) end if end proc:
# series expansion of e(q)/e(-q)
S := series( simplify( Mat_prod(14, q)[1, 1])/simplify(Mat_prod(14, q)[2, 1]) * simplify(Mat_prod(14, -q)[2, 1])/simplify(Mat_prod(14, -q)[1, 1]), q, 101 ):
seq(coeftayl(S, q = 0, n), n = 0..100);
CROSSREFS
KEYWORD
sign,easy
AUTHOR
Peter Bala, Jan 06 2026
STATUS
approved
