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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A009002 Expansion of (1+x)/cos(x). 1
1, 1, 1, 3, 5, 25, 61, 427, 1385, 12465, 50521, 555731, 2702765, 35135945, 199360981, 2990414715, 19391512145, 329655706465, 2404879675441, 45692713833379, 370371188237525, 7777794952988025, 69348874393137901, 1595024111042171723, 15514534163557086905 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,4

LINKS

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

FORMULA

E.g.f.: (1-x^2)/U(0) where U(k)= 1 - x/(1 - x/(x + (2*k+1)*(2*k+2)/U(k+1)) ; (continued fraction, 3-step). - Sergei N. Gladkovskii, Oct 17 2012

a(n) ~ n! * 2^n * (Pi + 2 + (-1)^(n+1) * (Pi - 2)) / Pi^(n+1). - Vaclav Kotesovec, Jan 22 2015

MAPLE

seq(coeff(series(factorial(n)*(1+x)/cos(x), x, n+1), x, n), n=0..25); # Muniru A Asiru, Jul 21 2018

MATHEMATICA

With[{nn=40}, Take[CoefficientList[Series[(1 + x)/Cos[x], {x, 0, nn}], x] Range[0, nn]!]] (* Vincenzo Librandi, Oct 24 2012 *)

PROG

(PARI) x='x+O('x^30); Vec(serlaplace((1+x)/cos(x))) \\ G. C. Greubel, Jul 21 2018

(MAGMA) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!((1+x)/Cos(x))); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Jul 21 2018

CROSSREFS

Sequence in context: A290509 A208800 A249935 * A119882 A276968 A074701

Adjacent sequences:  A008999 A009000 A009001 * A009003 A009004 A009005

KEYWORD

nonn

AUTHOR

R. H. Hardin

EXTENSIONS

Extended and formatted by Olivier Gérard, Mar 15 1997

More terms from Ralf Stephan, Mar 08 2004

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 November 14 07:19 EST 2019. Contains 329111 sequences. (Running on oeis4.)