OFFSET
1,2
FORMULA
a(n) = ((2*n+1)!*sum(m=0..n, binomial(n-1/2,n-m)/(2*m+1)!*sum(k=1..2*m+1, (-1)^(n-m+k+1)*k!*2^(2*m+1-k)*Stirling2(2*m+1,k)))). - Vladimir Kruchinin, Jun 17 2011
E.g.f.: tanh(x/sqrt(1+x^2)) = (x/sqrt(1+x^2))*G(0) where G(k)= 1 - x^2/(x^2 + (1+x^2)*(2*k+1)*(2*k+3)/G(k+1)); (continued fraction, 2-step). - Sergei N. Gladkovskii, Aug 06 2012
a(n) ~ (2*n-1)! * (-1)^(n+1) * 16 * (4+Pi^2)^(n-3/2) / Pi^(2*n). - Vaclav Kotesovec, Feb 02 2015
EXAMPLE
tanh(sin(arctan(x))) = x - (5/3!)*x^3 + (121/5!)*x^5 - (6677/7!)*x^7 + (651985/9!)*x^9 - ...
MATHEMATICA
nn = 20; Table[(CoefficientList[Series[Tanh[x/Sqrt[1 + x^2]], {x, 0, 2*nn+1}], x] * Range[0, 2*nn+1]!)[[n]], {n, 2, 2*nn, 2}] (* Vaclav Kotesovec, Feb 02 2015 *)
PROG
(Maxima)
a(n):=((2*n+1)!*sum(binomial(n-1/2, n-m)/(2*m+1)!*sum((-1)^(n-m+k+1)*k!*2^(2*m+1-k)*stirling2(2*m+1, k), k, 1, 2*m+1), m, 0, n)); /* Vladimir Kruchinin, Jun 17 2011 */
CROSSREFS
KEYWORD
sign
AUTHOR
Patrick Demichel (patrick.demichel(AT)hp.com)
STATUS
approved