login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A120083 Denominators of expansion for Debye function for n=1: D(1,x). 10

%I #15 May 01 2023 10:03:15

%S 1,4,36,1,3600,1,211680,1,10886400,1,526901760,1,16999766784000,1,

%T 1120863744000,1,181400588328960000,1,97072790126247936000,1,

%U 16860010916664115200000,1,324325300906011525120000,1

%N Denominators of expansion for Debye function for n=1: D(1,x).

%C Numerators are found under A120082.

%H G. C. Greubel, <a href="/A120083/b120083.txt">Table of n, a(n) for n = 0..447</a>

%F a(n) = denominator(r(n)), with r(n) = [x^n]( 1 - x/4 + Sum_{k >= 0}(B(2*k)/((2*k+1)*(2*k)!))*x^(2*k) ), |x|<2*pi. B(2*k) = A000367(k)/A002445(k) (Bernoulli numbers).

%F a(n) = denominator(B(n)/(n+1)!), n >= 0. See the comment on the e.g.f. D(1,x) in A120082. - _Wolfdieter Lang_, Jul 15 2013

%t Table[Denominator[BernoulliB[n]/(n+1)!], {n,0,50}] (* _G. C. Greubel_, May 01 2023 *)

%o (Magma) [Denominator(Bernoulli(n)/Factorial(n+1)): n in [0..50]]; // _G. C. Greubel_, May 01 2023

%o (SageMath)

%o def A120083(n): return denominator(bernoulli(n)/factorial(n+1))

%o [A120083(n) for n in range(51)] # _G. C. Greubel_, May 01 2023

%Y Cf. A120082.

%K nonn,easy,frac

%O 0,2

%A _Wolfdieter Lang_, Jul 20 2006

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 02:23 EDT 2024. Contains 371264 sequences. (Running on oeis4.)