|
%I M1810 N0716 #41 Jun 01 2022 01:48:14
%S 0,1,2,7,44,361,3654,44207,622552,10005041,180713290,3624270839,
%T 79914671748,1921576392793,50040900884366,1403066801155039,
%U 42142044935535536,1349948504738292193,45940391206037470098
%N Nearest integer to modified Bessel function K_n(1).
%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 Vincenzo Librandi, <a href="/A000155/b000155.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>
%p Digits := 60: A000155 := proc(n) round( evalf ( BesselK( n, 1) )); end;
%p a := proc(n) options remember;if n<3 then n elif n=3 then 7 else a(n-4)+2*(n-3)*a(n-3)+2*(n-1)*a(n-1) fi end: # _Mark van Hoeij_, Nov 21 2011
%p series(hypergeom([1,1],[],2*x/(1+x^2))*x/(1+x^2), x=0, 20); # _Mark van Hoeij_, Nov 21 2011
%t Table[Round[BesselK[n, 1]], {n, 0, 18}] (* _Ray Chandler_, Nov 28 2006 *)
%o (PARI) a(n)=round(besselk(n,1)) \\ _Charles R Greathouse IV_, May 11 2016
%K nonn
%O 0,3
%A _N. J. A. Sloane_
%E Extended by _Ray Chandler_, Nov 28 2006
%E More accurate definition by _Sean A. Irvine_, Mar 28 2010
|
