%I #21 Mar 22 2022 21:27:54
%S 48,768,30720,2211840,247726080,39636172800,8561413324800,
%T 2397195730944000,843812897292288000,364527171630268416000,
%U 189554129247739576320000,116765343616607579013120000
%N Denominators of coefficients of expansion of Bessel function J_3(x).
%D Bronstein-Semendjajew, Taschenbuch der Mathematik, 7th german ed. 1965, ch. 4.4.7
%H T. D. Noe, <a href="/A014401/b014401.txt">Table of n, a(n) for n=0..50</a>
%H <a href="/index/Be#Bessel">Index entries for sequences related to Bessel functions or polynomials</a>
%F a(n) = 2^(2n+k) * n! * (n+k)! here for k=3, i.e., Bessel's J3(x).
%F D-finite with recurrence: a(n) - (4*n^2 + 4*n*k)*a(n-1) = 0, a(0) = 2^k*k!, here for k=3. - _Georg Fischer_, Mar 22 2022
%e a(1) = 768 = 32*24, J3(x) = x^3/48 - x^5/768 + x^7/30720 - x^9/2211840 +- ...
%p k:=3: f:= gfun:-rectoproc({a(n)-(4*n^2 + 4*n*k)*a(n-1), a(0)=2^k*k!}, a(n), remember): map(f, [$0..16]); # _Georg Fischer_, Mar 22 2022
%t Denominator[Take[CoefficientList[Series[BesselJ[3,x],{x,0,30}],x],{4,-1,2}]] (* _Harvey P. Dale_, Dec 10 2011 *)
%Y J0: A002454, J1: A002474, J2: A002506.
%K nonn
%O 0,1
%A _N. J. A. Sloane_