login
Numerators of Taylor series for log(1/cos(x)). Also from log(cos(x)).
10

%I #36 Mar 24 2020 03:30:19

%S 0,1,1,1,17,31,691,10922,929569,3202291,221930581,9444233042,

%T 56963745931,29435334228302,2093660879252671,344502690252804724,

%U 129848163681107301953,868320396104950823611,209390615747646519456961,28259319101491102261334882

%N Numerators of Taylor series for log(1/cos(x)). Also from log(cos(x)).

%D L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 88.

%D CRC Standard Mathematical Tables and Formulae, 30th ed. 1996, p. 42.

%H T. D. Noe, <a href="/A046990/b046990.txt">Table of n, a(n) for n = 0..100</a>

%F Let q(n) = Sum_{k=0..n-1} (-1)^k*A201637(n-1,k) then a(n) = numerator((-1)^(n-1)*q(2*n)/(2*n)!). - _Peter Luschny_, Nov 16 2012

%e log(1/cos(x)) = 1/2*x^2+1/12*x^4+1/45*x^6+17/2520*x^8+31/14175*x^10+...

%e log(cos(x)) = -(1/2*x^2+1/12*x^4+1/45*x^6+17/2520*x^8+31/14175*x^10+...).

%p q:= proc(n) add((-1)^k*combinat[eulerian1](n-1,k), k=0..n-1) end: A046990:= n -> numer((-1)^(n-1)*q(2*n)/(2*n)!):

%p seq(A046990(n),n=0..19); # _Peter Luschny_, Nov 16 2012

%t Join[{0},Numerator[Select[CoefficientList[Series[Log[1/Cos[x]],{x,0,40}], x],#!=0&]]] (* _Harvey P. Dale_, Jul 27 2011 *)

%t a[n_] := Numerator[((-4)^n-(-16)^n)*BernoulliB[2*n]/2/n/(2*n)!]; a[0] = 0; Table[a[n], {n, 0, 19}] (* _Jean-François Alcover_, Feb 11 2014, after _Charles R Greathouse IV_ *)

%o (Sage) # uses[eulerian1 from A173018]

%o def A046990(n):

%o def q(n):

%o return add((-1)^k*eulerian1(n-1, k) for k in (0..n-1))

%o return ((-1)^(n-1)*q(2*n)/factorial(2*n)).numer()

%o [A046990(n) for n in (0..19)] # _Peter Luschny_, Nov 16 2012

%o (PARI) a(n)=numerator(((-4)^n-(-16)^n)*bernfrac(2*n)/2/n/(2*n)!) \\ _Charles R Greathouse IV_, Nov 06 2013

%o (PARI) {a(n) = if( n<1, 0, my(m = 2*n); numerator( polcoeff( -log(cos(x + x * O(x^m))), m)))}; /* _Michael Somos_, Jun 03 2019 */

%Y Cf. A046991, A002430, A050970.

%K nonn,easy,frac,nice

%O 0,5

%A _N. J. A. Sloane_