OFFSET
0,3
COMMENTS
Compare to: exp( Sum_{n>=1} x^(2*n)/(2*n) ) = 1/sqrt(1-x^2).
EXAMPLE
E.g.f.: A(x) = 1 + x^2/2! + 12*x^4/4! + 370*x^6/6! + 21490*x^8/8! +...
where
log(A(x)) = x^2/2 + (1+1/2)*x^4/4 + (1+1/2+1/3)*x^6/6 + (1+1/2+1/3+1/4)*x^8/8 + (1+1/2+1/3+1/4+1/5)*x^10/10 + (1+1/2+1/3+1/4+1/5+1/6)*x^12/12 +...
Explicitly,
log(A(x)) = x^2/2! + 9*x^4/4! + 220*x^6/6! + 10500*x^8/8! + 828576*x^10/10! + 97796160*x^12/12! + 16145775360*x^14/14! + 3554072121600*x^16/16! +...
PROG
(PARI) {H(n)=sum(k=1, n, 1/k)}
{a(n)=local(A=1); A=exp(sum(k=1, n\2+1, H(k)*x^(2*k)/(2*k))+x*O(x^n)); n!*polcoeff(A, n)}
for(n=0, 25, print1(a(2*n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jan 13 2014
STATUS
approved