A046991 -id:A046991 - cp's OEIS frontend

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

0, 1, 1, 1, 17, 31, 691, 10922, 929569, 3202291, 221930581, 9444233042, 56963745931, 29435334228302, 2093660879252671, 344502690252804724, 129848163681107301953, 868320396104950823611, 209390615747646519456961, 28259319101491102261334882

Offset: 0

N. J. A. Sloane

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

L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 88.
CRC Standard Mathematical Tables and Formulae, 30th ed. 1996, p. 42.
Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 32, equation 32:6:3 at page 301.

T. D. Noe, Table of n, a(n) for n = 0..100

Cf. A046991, A002430, A050970.

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)!):
seq(A046990(n),n=0..19);  # Peter Luschny, Nov 16 2012

Join[{0},Numerator[Select[CoefficientList[Series[Log[1/Cos[x]],{x,0,40}], x],#!=0&]]] (* Harvey P. Dale, Jul 27 2011 *)
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 *)

a(n)=numerator(((-4)^n-(-16)^n)*bernfrac(2*n)/2/n/(2*n)!) \\ Charles R Greathouse IV, Nov 06 2013

{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 */

# uses[eulerian1 from A173018]
def A046990(n):
    def q(n):
        return add((-1)^k*eulerian1(n-1, k) for k in (0..n-1))
    return ((-1)^(n-1)*q(2*n)/factorial(2*n)).numer()
[A046990(n) for n in (0..19)]  # Peter Luschny, Nov 16 2012

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

A118817 Decimal expansion of Product_{n >= 1} cos(1/n).

Original entry on oeis.org

3, 8, 8, 5, 3, 6, 1, 5, 3, 3, 3, 5, 1, 7, 5, 8, 5, 9, 1, 8, 4, 3, 2, 9, 5, 7, 5, 6, 8, 7, 0, 3, 5, 9, 0, 5, 0, 1, 3, 9, 0, 0, 5, 2, 8, 5, 9, 7, 5, 1, 7, 9, 2, 1, 9, 1, 3, 1, 8, 4, 6, 1, 1, 9, 9, 8, 7, 9, 8, 7, 4, 9, 4, 3, 4, 6, 3, 3, 9, 3, 2, 7, 6, 8, 3, 8, 8, 4, 3, 1, 9, 7, 8, 1, 3, 8, 3, 4, 0, 8, 2, 2, 4, 1, 3

Offset: 0

Fredrik Johansson, May 23 2006

			0.38853615333517585918432957568703590501390...

Charles R Greathouse IV, Table of n, a(n) for n = 0..10000

Cf. A046990, A046991.
Cf. A051762, A085365, A246945, A230821, A249673.
Cf. A027641, A027642.

nn:= 120:
p:= product(cos(1/n), n=1..infinity):
f:= evalf(p, nn+10):
s:= convert(f, string):
seq(parse(s[n+1]), n=1..nn);  # Alois P. Heinz, Nov 04 2013

S = Series[Log[Cos[x]], {x, 0, 400}]; N[Exp[N[Sum[SeriesCoefficient[S, 2k] Zeta[2k], {k, 1, 200}], 70]], 50]
Block[{$MaxExtraPrecision = 1000}, Do[Print[N[1/Exp[Sum[(2^(2*n) - 1)*Zeta[2*n]^2/(n*Pi^(2*n)), {n, 1, m}]], 110]], {m, 250, 300}]] (* Vaclav Kotesovec, Sep 20 2014 *)

exp(-sumpos(n=1,-log(cos(1/n)))) \\ warning: requires 2.6.2 or greater; Charles R Greathouse IV, Nov 04 2013

T(n)=((-4)^n-(-16)^n)*bernfrac(2*n)/2/n/(2*n)!
lm=lambertw(2*log(Pi/2)*10^default(realprecision))/2/log(Pi/2); exp(-sum(n=1,lm,T(n)*zeta(2*n))) \\ Charles R Greathouse IV, Nov 06 2013

Equals exp(Sum_{n>=1} -c(n)*zeta(2*n)), where c(n) = A046990(n)/A046991(n).
Equals exp(-Sum_{n>=1} (2^(2*n)-1) * Zeta(2*n)^2 / (n*Pi^(2*n)) ). - Vaclav Kotesovec, Sep 20 2014
Equals exp(Sum_{k>=1} (-1)^k*2^(2*k-1)*(2^(2*k)-1)*B(2*k)*zeta(2*k)/(k*(2*k)!)), where B(k) is the k-th Bernoulli number. - Amiram Eldar, Jul 30 2023

Corrected offset and extended by Robert G. Wilson v, Nov 03 2013

cp's OEIS Frontend

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

Views

Author

Keywords

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Sage

Formula

A118817 Decimal expansion of Product_{n >= 1} cos(1/n).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

Extensions