A051149 Continued fraction for BesselK(1,2)/BesselK(0,2).

1, 4, 2, 1, 1, 2, 8, 1, 6, 1, 4, 2, 1, 8, 1, 1, 7, 1, 6, 4, 2, 1, 6, 1, 1, 2, 1, 1, 1, 2, 5, 1, 3, 3, 1, 3, 1, 8, 1, 1, 1, 1, 1, 2, 3, 8, 3, 29, 1, 7, 1, 57, 1, 121, 2, 14, 2, 8, 1, 1, 16, 1, 3, 1, 5, 1, 5, 1, 4, 17, 1, 5, 6, 1, 3, 2, 9, 7, 1, 4, 4, 1, 1, 16, 3, 5, 2, 1, 2, 1

Offset: 0

Author

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

Keywords

• nonn
• cofr

Examples

			1.228036929818907975742672452...
BesselK(1,2)/BesselK(0,2) = 1 + 1/(4 + 1/(2 + 1/(1 + 1/(1 + ...)))). - _Harry J. Smith_, Apr 29 2009
		

Links

• Harry J. Smith, Table of n, a(n) for n = 0..2000
• G. Xiao, Contfrac
• Index entries for sequences related to Bessel functions or polynomials
• Index entries for continued fractions for constants

Crossrefs

Cf. A051148, A000994, A000995.

Programs

Mathematica

ContinuedFraction[BesselK[1,2]/BesselK[0,2],90] (* Harvey P. Dale, Sep 07 2018 *)

PARI

{ allocatemem(932245000); default(realprecision, 2100); x=contfrac(besselk(1,2)/besselk(0,2)); for (n=1, 2001, write("b051149.txt", n-1, " ", x[n])); } \\ Harry J. Smith, Apr 29 2009