A259490 Numbers k such that the denominator of the n-th convergent of the continued fraction expansion of Pi is prime.
2, 4, 9, 33, 595, 1127, 2003, 3611, 4356, 6926
Offset: 1
Programs
-
Mathematica
$MaxExtraPrecision = 25000; lst = {}; cf = ContinuedFraction[Pi, 10000]; Do[ If[ PrimeQ[ Denominator[ FromContinuedFraction[ Take[ cf, n]] ]], AppendTo[lst, n]], {n, Length[cf]}]; lst Position[Convergents[Pi,7000],?(PrimeQ[Denominator[#]]&)]//Flatten (* _Harvey P. Dale, Aug 12 2021 *)
Comments