%I #19 Dec 08 2023 23:37:00
%S 1,-2,1,0,-1,0,1,0,-1,0,5,0,-691,0,7,0,-3617,0,43867,0,-174611,0,
%T 854513,0,-236364091,0,8553103,0,-23749461029,0,8615841276005,0,
%U -7709321041217,0,2577687858367,0,-26315271553053477373,0,2929993913841559,0,-261082718496449122051
%N Numerators of rationals with e.g.f. D(4,x), a Debye function.
%C The denominators are given in A227574.
%C For general remarks on the e.g.f.s D(n,x), the Debye function with index n = 1, 2, 3, ... see the W. Lang link under A120080.
%C D(4,x) := (4/x^4)*int(t^4/(exp(x) - 1), t=0..x) is the e.g.f. of the rationals r(4,n) = 4*B(n)/(n+4), n >= 0, with the Bernoulli numbers B(n) = A027641(n)/A027642(n).
%C See the Abramowitz-Stegun reference for the integral appearing in
%C D(4,x) and a series expansion valid for |x| < 2*Pi.
%C Initially coincides with A227570, A176327, A164555 and A027641 for n <> 1. - _R. J. Mathar_, Jul 19 2013
%C Differs from these sequences for n = 1328, 2660, 2828, 2880... - _Andrey Zabolotskiy_, Dec 08 2023
%D M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, Tenth Printing, 1972, pp. 998, equ. 27.1.1 for n=4, with a factor (x^4)/4 extracted.
%F a(n) = numerator(4*B(n)/(n+4)), n >= 0, with the Bernoulli numbers B(n).
%e The rationals r(4,n), n=0..15 are: 1, -2/5, 1/9, 0, -1/60, 0, 1/105, 0, -1/90, 0, 5/231, 0, -691/10920, 0, 7/27, 0.
%t A227573[n_]:=Numerator[4BernoulliB[n]/(n+4)];
%t Array[A227573,50,0] (* _Paolo Xausa_, Dec 08 2023 *)
%o (Sage)
%o print([(bernoulli(n)*4/(n+4)).numerator() for n in range(30)]) # _Andrey Zabolotskiy_, Dec 08 2023
%Y Cf. A227570, A227574, A027641/A027642, A120086/A120087 (D(4,x) as o.g.f.).
%K sign,easy,frac
%O 0,2
%A _Wolfdieter Lang_, Jul 17 2013