OFFSET
0,3
LINKS
Philippe Flajolet, Xavier Gourdon, and Philippe Dumas, Mellin transforms and asymptotics: harmonic sums, Special volume on mathematical analysis of algorithms. Theoret. Comput. Sci. 144 (1995), no. 1-2, 3-58.
FORMULA
a(n) = numerator(b(n)), where b(n) = Sum_{k=1..n} b(n-k)*(k/(2*n)-1)/(2*k)!, with b(0)=1. - Tani Akinari, Sep 17 2023
a(n) = numerator((-1)^n*A186491(n)/(4^n*(2*n)!)). - Andrew Howroyd, Sep 19 2023
EXAMPLE
1/sqrt(cosh(x)) = 1 - (1/4)*x^2 + (7/96)*x^4 - (139/5760)*x^6 + (5473/645120)*x^8 - (51103/16588800)*x^10 + ...
MAPLE
a:= n-> numer(coeff(series(1/sqrt(cosh(x)), x, 2*n+1), x, 2*n)):
seq(a(n), n=0..19); # Alois P. Heinz, Sep 19 2023
PROG
(Maxima) b[n]:=if n=0 then 1 else sum(b[n-k]*(k/n/2-1)/(2*k)!, k, 1, n)$ a[n]:=num(b[n])$
makelist(a[n], n, 0, 20); /* Tani Akinari, Sep 17 2023 */
CROSSREFS
KEYWORD
sign,frac
AUTHOR
N. J. A. Sloane, May 05 2011
STATUS
approved