OFFSET
0,2
MAPLE
assume(x > 0): egf := (log(2*x - 1) - log(x - 1) - log(x + 1))/(x - 1):
ser := series(egf, x, 24): seq(n!*coeff(ser, x, n), n = 0..20);
# Recurrence:
a := proc(n) option remember; if n < 4 then return [0, 2, 6, 34][n+1] fi;
(2*(n-4)^4 + 14*(n-4)^3 + 34*(n-4)^2 + 34*n - 124)*a(n-4) - (3*(n-4)^3 + 19*(n-4)^2 + 38*n - 128)*a(n-3) - ((n-4)^2 + 7*n - 16)*a(n-2) + (3*n - 2)*a(n-1) end:
seq(a(n), n = 0..20);
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Apr 09 2024
STATUS
approved