%I M1938 N0767 #33 Dec 09 2023 07:04:23
%S 0,0,0,1,2,9,49,306,2188,17810,162482,1642635,18231462,220420179,
%T 2883693795,40592133316,611765693528,9828843229764,167702100599524,
%U 3028466654021205,57708568527002410,1157199837194069405
%N Nearest integer to modified Bessel function K_n(2).
%D M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 429.
%D N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H Alois P. Heinz, <a href="/A000167/b000167.txt">Table of n, a(n) for n = 0..200</a>
%H M. Abramowitz and I. A. Stegun, eds., <a href="http://www.convertit.com/Go/ConvertIt/Reference/AMS55.ASP">Handbook of Mathematical Functions</a>, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
%H <a href="/index/Be#Bessel">Index entries for sequences related to Bessel functions or polynomials</a>
%F b(n) = (n-1)*b(n-1) + b(n-2) with b(n) = K_n(2). - _Christian Krause_, Dec 08 2023
%p Digits := 60: A000167 := proc(n) round( evalf ( BesselK( n,2 ) )); end;
%t Table[BesselK[n, 2] // Round, {n, 0, 21}] (* _Jean-François Alcover_, Mar 12 2014 *)
%o (PARI) a(n)=round(besselk(n,2)) \\ _Charles R Greathouse IV_, Jul 29 2016
%K nonn
%O 0,5
%A _N. J. A. Sloane_
%E More terms from _Herman P. Robinson_