%I #6 Feb 12 2024 13:29:38
%S 0,1,4,31,368,5857,116772,2796671,78190016,2499283841,89896028260,
%T 3593341846559,158017145220336,7581229628729569,394065923548717252,
%U 22060110489099436543,1323212563422417475328,84663543948545618984449,5755797775937679673467204,414332776323564390870654239
%N a(n) = b(n, 1/2) where b(n, x) = (Pi/4)*(Y(0, x)*J(n, x) - J(0, x)*Y(n, x)) and Y, J are Bessel functions.
%p b := (n, x) -> (Pi/4)*(BesselY(0, x)*BesselJ(n, x)-BesselJ(0, x)*BesselY(n, x)):
%p a := n -> simplify(b(n, 1/2)): seq(a(n), n = 0..19);
%t a = (Pi/4)*(BesselY[0, 1/2] * BesselJ[n, 1/2] - BesselJ[0, 1/2] * BesselY[n, 1/2]); Table[Round[a], {n, 0, 19}]
%Y Cf. A093985, A036243.
%K nonn
%O 0,3
%A _Peter Luschny_, Jan 30 2024
|