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.

Showing 1-3 of 3 results.

A120086 Numerators of expansion of Debye function for n=4: D(4,x).

Original entry on oeis.org

1, -2, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -691, 0, 1, 0, -3617, 0, 43867, 0, -174611, 0, 77683, 0, -236364091, 0, 657931, 0, -3392780147, 0, 1723168255201, 0, -7709321041217, 0, 151628697551, 0, -26315271553053477373
Offset: 0

Views

Author

Wolfdieter Lang, Jul 20 2006

Keywords

Comments

Denominators are found under A120087.
See the W. Lang link under A120080 for more details on the general case D(n,x), n= 1, 2, ... D(4,x) is the e.g.f. of the rational sequence {4*B(n)/(n+4)}, n >= 0. See A227573/A227574. - Wolfdieter Lang, Jul 17 2013

Examples

			Rationals r(n): [1, -2/5, 1/18, 0, -1/1440, 0, 1/75600, 0, -1/3628800, 0, 1/167650560, 0, -691/5230697472000, ...].
		

Crossrefs

Cf. A060054. [From R. J. Mathar, Aug 07 2008]
Cf. A000367/A002445, A027641/A027642, A120097, A227573/A227574 (D(4,x) as e.g.f.). - Wolfdieter Lang, Jul 17 2013

Programs

  • Magma
    [Numerator(4*(n+1)*(n+2)*(n+3)*Bernoulli(n)/Factorial(n+4)): n in [0..50]]; // G. C. Greubel, May 02 2023
    
  • Mathematica
    r[n_]:= 4*BernoulliB[n]/((n+4)*n!); Table[r[n]//Numerator, {n,0,36}] (* Jean-François Alcover, Jun 21 2013 *)
  • SageMath
    [numerator(4*(n+1)*(n+2)*(n+3)*bernoulli(n)/factorial(n+4)) for n in range(51)] # G. C. Greubel, May 02 2023

Formula

a(n) = numerator(r(n)), with r(n) = [x^n](1 - 4*x/(2*(4+1)) + 2*Sum_{k >= 0} (B(2*k)/((k+2)*(2*k)!))*x^(2*k) ), |x| < 2*Pi. B(2*k) = A000367(k)/A002445(k) (Bernoulli numbers).
a(n) = numerator(4*B(n)/((n+4)*n!)), n >= 0, with the Bernoulli numbers B(n) = A027641(n)/A027642(n). From D(4,x) read as o.g.f. - Wolfdieter Lang, Jul 17 2013

A120087 Denominators of expansion of Debye function for n=4: D(4,x).

Original entry on oeis.org

1, 5, 18, 1, 1440, 1, 75600, 1, 3628800, 1, 167650560, 1, 5230697472000, 1, 336259123200, 1, 53353114214400000, 1, 28100018194440192000, 1, 4817145976189747200000, 1, 91657150256046735360000, 1, 11856768957122205686169600000, 1, 1396008903899788738560000000, 1
Offset: 0

Views

Author

Wolfdieter Lang, Jul 20 2006

Keywords

Comments

From Wolfdieter Lang, Jul 17 2013: (Start)
The numerators are given in A120086.
See the link under A120080 for D(n,4) as e.g.f. of 4*B(n)/(n+4) = A227573(n)/A227574(n), n>= 0. (End)

Examples

			Rationals r(n): [1, -2/5, 1/18, 0, -1/1440, 0, 1/75600, 0, -1/3628800, 0, 1/167650560, 0, -691/5230697472000, ...].
		

Crossrefs

Programs

  • Magma
    [Denominator(4*(n+1)*(n+2)*(n+3)*Bernoulli(n)/Factorial(n+4)): n in [0..50]]; // G. C. Greubel, May 02 2023
    
  • Mathematica
    Table[Denominator[4*BernoulliB[n]/((n+4)*n!)], {n,0,50}] (* G. C. Greubel, May 02 2023 *)
  • SageMath
    [denominator(4*(n+1)*(n+2)*(n+3)*bernoulli(n)/factorial(n+4)) for n in range(51)] # G. C. Greubel, May 02 2023

Formula

a(n) = denominator(r(n)), with r(n) = [x^n](1 - 2*x/5 + 2*Sum_{k >= 0}(B(2*k)/((k+2)*(2*k)!))*x^(2*k) ), |x| < 2*Pi. B(2*k) = A000367(k)/A002445(k) (Bernoulli numbers).
a(n) = denominator(4*B(n)/((n+4)*n!)), n >= 0, with the Bernoulli numbers B(n) = A027641(n)/A027642(n). From D(4,x) read as o.g.f. - Wolfdieter Lang, Jul 17 2013

A120084 Numerators of expansion for Debye function for n=2: D(2,x).

Original entry on oeis.org

1, -1, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -691, 0, 1, 0, -3617, 0, 43867, 0, -174611, 0, 77683, 0, -236364091, 0, 657931, 0, -3392780147, 0, 1723168255201, 0, -7709321041217, 0, 151628697551, 0, -26315271553053477373, 0, 154210205991661, 0, -261082718496449122051
Offset: 0

Views

Author

Wolfdieter Lang, Jul 20 2006

Keywords

Comments

Denominators are found under A120085.
This sequence appears to coincide with A120082.

Examples

			Rationals r(n): [1, -1/3, 1/24, 0, -1/2160, 0, 1/120960, 0, -1/6048000, 0, 1/287400960, ...].
		

Crossrefs

Programs

  • Magma
    [Numerator(2*(n+1)*Bernoulli(n)/Factorial(n+2)): n in [0..50]]; // G. C. Greubel, May 02 2023
    
  • Mathematica
    max = 38; Numerator[CoefficientList[Integrate[Normal[Series[(2*(t^2/(Exp[t]-1)))/x^2, {t, 0, max}]], {t, 0, x}], x]] (* Jean-François Alcover, Oct 04 2011 *)
    Table[Numerator[2*(n+1)*BernoulliB[n]/(n+2)!], {n,0,50}] (* G. C. Greubel, May 02 2023 *)
  • SageMath
    [numerator(2*(n+1)*bernoulli(n)/factorial(n+2)) for n in range(51)] # G. C. Greubel, May 02 2023

Formula

a(n) = numerator(r(n)), with r(n) = [x^n]( 1 - x/3 + Sum_{k >= 1} (B(2*k)/((k+1)*(2*k)!))*x^(2*k) ), |x|<2*Pi. B(2*k) = A000367(k)/A002445(k) (Bernoulli numbers).
a(n) = numerator(2*B(n)/((n+2)*n!)), n >= 0. See the comment on the e.g.f. D(2,x) in A120085. - Wolfdieter Lang, Dec 03 2022
Showing 1-3 of 3 results.