A051148 Decimal expansion of BesselK(1,2)/BesselK(0,2).

1, 2, 2, 8, 0, 3, 6, 9, 2, 9, 8, 1, 8, 9, 0, 7, 9, 7, 5, 7, 4, 2, 6, 7, 2, 4, 5, 2, 0, 8, 1, 7, 1, 8, 9, 0, 4, 2, 3, 9, 3, 6, 6, 9, 5, 2, 6, 6, 0, 8, 5, 9, 6, 9, 4, 4, 6, 7, 0, 3, 8, 1, 4, 9, 1, 3, 7, 4, 7, 1, 9, 9, 1, 5, 7, 0, 3, 8, 9, 8, 3, 9, 7, 0, 5, 3, 8, 1, 5, 2, 4, 7, 3, 8, 2, 2, 9, 2, 1

Offset: 1

Author

N. J. A. Sloane, E. M. Rains

Keywords

• nonn
• cons
• easy

Examples

			1.228036929818907975742672452...
		

Links

• Harry J. Smith, Table of n, a(n) for n = 1..2000
• Index entries for sequences related to Bessel functions or polynomials

Crossrefs

Cf. A051149, A000994, A000995.

Programs

Mathematica

RealDigits[ BesselK[1, 2] / BesselK[0, 2], 10, 99] // First (* Jean-François Alcover, Mar 07 2013 *)

PARI

{ allocatemem(932245000); default(realprecision, 2080); x=besselk(1,2)/besselk(0,2); for (n=1, 2000, d=floor(x); x=(x-d)*10; write("b051148.txt", n, " ", d)); } \\ Harry J. Smith, Apr 29 2009

Extensions

Fixed my PARI program, had -n. - Harry J. Smith, May 19 2009