%I #32 Mar 16 2020 14:24:10
%S 1,-1,1,-1,-67,-1,-64397,-113249,-3679787,-810304169,-6040635661561,
%T -428305999661,-16827172241810597,-5620292762592913,
%U -1550760014054450957,-4168373361283100017,-8551022502876237590534947
%N Numerators of coefficients in expansion of sqrt(sin(x)/x) (even powers only).
%H G. C. Greubel, <a href="/A008991/b008991.txt">Table of n, a(n) for n = 0..255</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/SincFunction.html">The Sinc Function</a>.
%F sum(n>=0, A008991(n)/A008992(n) ) = A117017 - _Johannes W. Meijer_, Feb 10 2013
%p A008991 := n -> numer(coeff(taylor(sqrt(sin(x)/x), x, 2*n+2), x, 2*n)): seq(A008991(n), n=0..16); # _Johannes W. Meijer_, Feb 10 2013
%t Numerator[CoefficientList[Series[Sqrt[Sin[x]/x], {x, 0, 50}], x][[1 ;; -1 ;; 2]]] (* _G. C. Greubel_, Jul 21 2018 *)
%o (Sage)
%o length = 16; T = taylor(sqrt(sin(x)/x),x,0,2*length+2)
%o [T.coefficient(x, 2*n).numerator() for n in (0..length)]
%o # _Peter Luschny_, Dec 13 2012
%Y Denominators are in A008992.
%Y Appears in A220002 and A220466.
%K sign
%O 0,5
%A _N. J. A. Sloane_
|