login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A009453 E.g.f. sin(log(1 + tanh(x))). 1
0, 1, -1, -1, 8, -14, -46, 314, 48, -7244, 13484, 233804, -1002112, -10692824, 79853864, 661278344, -7758425472, -52788251504, 928828095344, 5226295555184, -135732896444032, -617858881186784, 23892314499643424, 83518025870961824 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,5

COMMENTS

Related to A102573: letting T(q,r) be the coefficient of n^(r+1) in the polynomial 2^(q-n)/n times Sum_{k=0..n} binomial(n,k)*k^q, then A009453(x) equals (-1)^(x+1) times the real part of Sum_{k=0..x-1} T(x,k)*i^k, where i is the imaginary unit. See Mathematica code below. [John M. Campbell, Nov 17 2011]

LINKS

Vincenzo Librandi, Table of n, a(n) for n = 0..200

FORMULA

a(n) = Sum_{m=0..(n-1)/2} (Sum_{r=2*m+1..n} (stirling1(r,2*m+1) * Sum_{k=r..n} binomial(k-1,r-1)*k!*2^(n-k)*stirling2(n,k)*(-1)^(r+k+m)/r!)). - Vladimir Kruchinin, Jun 21 2011

MATHEMATICA

poly[q_] := 2^(q - n)/n*FunctionExpand[Sum[Binomial[n, k]*k^q, {k, 0, n}]]; T[q_, r_] := First[Take[CoefficientList[poly[q], n], {r + 1, r + 1}]]; Table[Re[Sum[T[x, k]*I^k, {k, 0, x - 1}]]*(-1)^(x + 1), {x, 0, 23}] (* John M. Campbell, Nov 17 2011 *)

With[{nn = 30}, CoefficientList[Series[Sin[Log[1 + Tanh[x]]], {x, 0, nn}], x] Range[0, nn]!] (* Vincenzo Librandi, Feb 09 2014 *)

PROG

(Maxima) a(n):=sum(sum((stirling1(r, 2*m+1)*sum(binomial(k-1, r-1)*k!*2^(n-k)*stirling2(n, k)*(-1)^(r+k+m), k, r, n))/r!, r, 2*m+1, n), m, 0, (n-1)/2); /* Vladimir Kruchinin, Jun 21 2011 */

CROSSREFS

Sequence in context: A059676 A242198 A111050 * A332658 A198307 A218555

Adjacent sequences:  A009450 A009451 A009452 * A009454 A009455 A009456

KEYWORD

sign,easy

AUTHOR

R. H. Hardin

EXTENSIONS

Extended with signs by Olivier Gérard, Mar 15 1997

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 10 08:09 EDT 2021. Contains 342845 sequences. (Running on oeis4.)