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

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A009843 Expansion of x/cos(x). 12
1, 3, 25, 427, 12465, 555731, 35135945, 2990414715, 329655706465, 45692713833379, 7777794952988025, 1595024111042171723, 387863354088927172625, 110350957750914345093747 (list; graph; refs; listen; history; internal format)
OFFSET

0,2

COMMENTS

Related to the formulae sum(k>0,sin(kx)/k^(2n+1))=(-1)^(n+1)/2*x^(2n+1)/(2n+1)!*sum(i=0,2n,(2Pi/x)^i*B(i)*C(2n+1,i)) and if x=Pi/2 sum(k>0,(-1)^(k+1)/k^(2n+1))=(-1)^n*E(2n)*Pi^(2n+1)/2^(2n+2)/(2n)! - Benoit Cloitre (benoit7848c(AT)orange.fr), May 01 2002

Expanding x/cosh(x) gives alt. signed values at odd positions.

LINKS

Peter Luschny, The lost Bernoulli numbers.

FORMULA

a(n)=(2n+1)*A000364(n)=sum(i=0, 2n, B(i)*C(2n+1, i)*4^i)=(2n+1)*E(2n) where B(i) are the Bernoulli numbers, C(2n, i) the binomial numbers and E(2n) the Euler numbers. - Benoit Cloitre (benoit7848c(AT)orange.fr), May 01 2002

Recurrence: a(n) = -(-1)^n*Sum[i=0..n-1, (-1)^i*a(i)*C(2n+1, 2i+1) ]. - Ralf Stephan, Feb 24 2005

a(n) = 4^n |E_{2n}(1/2)+E_{2n}(1)| (2n+1) for n > 0; E_{n}(x) Euler polynomial. [Peter Luschny, Nov 25 2010]

a(n) = (2*n+1)! * [x^(2*n+1)] x/cos(x).

E.g.f.: x / cos(x) =x+x^3/Q(0); Q(k)=8k+2-x^2/(1+(2k+1)*(2k+2)/Q(k+1)); (continued fraction). - Sergei N. Gladkovskii, Nov 15 2011

MAPLE

seq((2*i+1)!*coeff(series(x/cos(x), x, 32), x, 2*i+1), i=0..13);

A009843 := n -> (-1)^n*(2*n+1)*euler(2*n): [Peter Luschny]

MATHEMATICA

x/Cos[ x ] (* Odd Part *)

PROG

(PARI) for(n=0, 25, print1(sum(i=0, 2*n+1, binomial(2*n+1, i)*bernfrac(i)*4^i), ", "))

# (Python) The objective of this implementation is efficiency.

# n -> [a(0), a(1), ..., a(n)] for n > 0.

def A009843_list(n):

....S = [0 for i in range(0, n+1)]

....S[0] = 1

....for k in range(1, n+1):

........S[k] = k*S[k-1]

....for k in range(1, n+1):

........for j in range(k, n+1):

............S[j] = (j-k)*S[j-1]+(j-k+1)*S[j]

........S[k] = (2*k+1)*S[k]

....return S

print(A009843_list(10)) # - Peter Luschny, Aug 09 2011

CROSSREFS

Bisection of A009391, A009392, A065619, A083008.

Cf. A099028.

Sequence in context: A143925 A074708 A160143 * A182962 A136173 A003024

Adjacent sequences:  A009840 A009841 A009842 * A009844 A009845 A009846

KEYWORD

nonn

AUTHOR

R. H. Hardin (rhhardin(AT)att.net)

EXTENSIONS

Extended and signs tested Mar 15 1997 by Olivier Gerard.

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
Recent Additions | More pages | Superseeker | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified February 17 06:27 EST 2012. Contains 205998 sequences.