login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A009057 Expansion of e.g.f. x*cos(sinh(x)) (odd powers only). 1
1, -3, -15, 21, 1953, 32461, 77649, -18557595, -894805183, -21405607651, 421226202033, 90171090228021, 6418365899545825, 230309031988304109, -11126456208720437487, -3093396343171182148731, -342554566975833968005503 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

LINKS

G. C. Greubel, Table of n, a(n) for n = 0..250

FORMULA

a(n) = (2*n+1)*Sum_{m=1..2*n} (Sum_{i=0..m} (-1)^(i+m/2)*(m-2*i)^(2*n) *binomial(m, i))/(2^m*m!), n>0, a(0)=1. - Vladimir Kruchinin, Jun 29 2011

MATHEMATICA

With[{nn=40}, Take[CoefficientList[Series[x Cos[Sinh[x]], {x, 0, nn}], x] Range[ 0, nn-1]!, {2, -1, 2}]] (* Harvey P. Dale, May 22 2017 *)

PROG

(Maxima)

a(n):=if n=0 then 1 else (sum((sum((-1)^(i+m/2)*(m-2*i)^(2*n)*binomial(m, i), i, 0, m))/(2^m*m!), m, 1, 2*n))*(2*n+1); /* Vladimir Kruchinin, Jun 28 2011 */

(PARI) x='x+O('x^60); v=Vec(serlaplace(x*cos(sinh(x)))); vector((#v-1)\2 , n, v[2*n-1]) \\ G. C. Greubel, Jul 23 2018

CROSSREFS

Sequence in context: A087674 A212846 A276804 * A289712 A083795 A083793

Adjacent sequences:  A009054 A009055 A009056 * A009058 A009059 A009060

KEYWORD

sign

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 October 16 13:20 EDT 2019. Contains 328082 sequences. (Running on oeis4.)