%I #4 Sep 14 2014 03:42:52
%S -10,7,-3,1,0,1,5,31,222,1807,16485,166657,1849712,22363201,292571325,
%T 4118361751,62067997590,997206323191,17014575491837,307259565176257,
%U 5854946313840720,117406185841990657,2471384848995644517,54487872863746170031,1255692460715157555230
%N a(n) = 2*(K(n,2)*I(4,2) - (-1)^n*I(n,2)*K(4,2)) where I(n,x) and K(n,x) are Bessel functions.
%F a(n) = a(-n).
%F a(n+4) = A058309(n).
%p a := n -> 2*(BesselK(n,2)*BesselI(4,2)-(-1)^n*BesselI(n,2)* BesselK(4,2)); seq(round(evalf(a(n), 99)), n=0..24);
%Y Cf. A058309.
%K sign
%O 0,1
%A _Peter Luschny_, Sep 12 2014