cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A036282 Write cosec x = 1/x + Sum_{n>=1} e_n * x^(2n-1)/(2n-1)!; sequence gives numerators of e_n.

Original entry on oeis.org

1, 7, 31, 127, 511, 1414477, 8191, 118518239, 5749691557, 91546277357, 162912981133, 1982765468311237, 22076500342261, 455371239541065869, 925118910976041358111, 16555640865486520478399, 1302480594081611886641, 904185845619475242495834469891
Offset: 1

Views

Author

Keywords

Comments

From Johannes W. Meijer, May 24 2009: (Start)
Absolute value of numerator of [2^(2n-1) - 1] * Bernoulli(2n)/n.
Equals the absolute values of the numerators of the LS1[ -2*m,n=1] matrix coefficients of A160487 for m = 1, 2, .. ,.
(End)

Examples

			cosec x
= x^(-1) + 1/6*x + 7/360*x^3 + 31/15120*x^5 + ...
= x^(-1) + 1/6 * x/1! + 7/60 * x^3/3! + 31/126 * x^5/5! + ...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, Tenth Printing, 1972, p. 75 (4.3.68).

Crossrefs

Cf. A160487.
Differs from A282898.

Programs

  • Maple
    a:= n-> (m-> numer(coeff(series(csc(x), x, m+1), x, m)*m!))(2*n-1):
    seq(a(n), n=1..20);  # Alois P. Heinz, Jun 21 2018
  • Mathematica
    a[n_] := Abs[ Numerator[ (2^(2*n-1)-1) * BernoulliB[2*n]/n ] ]; Table[a[n], {n, 1, 18}] (* Jean-François Alcover, May 31 2013, after Johannes W. Meijer *)
  • PARI
    a(n) = abs(numerator((2^(2*n-1)-1)*bernfrac(2*n)/n)); \\ Michel Marcus, Mar 01 2015

Extensions

Title corrected and offset changed by Johannes W. Meijer, May 21 2009