login
Numerators of expansion for Debye function (D(1,x)) A120082 with 1's instead of 0's.
0

%I #12 Dec 08 2023 16:04:43

%S 1,-1,1,1,-1,1,1,1,-1,1,1,1,-691,1,1,1,-3617,1,43867,1,-174611,1,

%T 77683,1,-236364091,1,657931,1,-3392780147,1,1723168255201,1,

%U -7709321041217,1,151628697551,1,-26315271553053477373

%N Numerators of expansion for Debye function (D(1,x)) A120082 with 1's instead of 0's.

%C (Bernoulli numbers numerators) A027641(n)/a(n) is an integer sequence.

%C Note 1's in denominators A120083. A120082 and A027641 are analogous.

%t swfact[n_] := n!/Floor[n/2]!^2; a[n_] := 2*Zeta[n]*swfact[n]/(2*Pi)^n*If[Mod[n, 4] == 0, -1, 1]; a[0]=1; a[1]=-1; a[_?OddQ] = 0; Numerator[ Table[a[n], {n, 0, 36}]] /. (0 -> 1) (* _Jean-François Alcover_, Aug 09 2012 *)

%K uned,sign

%O 0,13

%A _Paul Curtz_, Aug 20 2008

%E Typo a(14)=7 instead of 1 fixed by _Jean-François Alcover_, Aug 09 2012