login
Continued fraction for prime constant A051006.
7

%I #20 Jan 03 2016 16:04:21

%S 0,2,2,2,3,12,131,1,7,1,2,1,3,3,1,2,5,39,2,1,169,2,2,2,1,1,2,5,1,2,1,

%T 199,24,7,7,1,163,1,3,2,1,2,14,1,3,5,1,1,1,1,3,3,2,1,279,1,4,1,1,4,1,

%U 92,1,16,2,3,1,5,2,25,9,1,1,2,8,1,2,5,1,1,3,1,1,2,1,7,1,1,1,5,1,2,2,5,6,1

%N Continued fraction for prime constant A051006.

%H Harry J. Smith, <a href="/A051007/b051007.txt">Table of n, a(n) for n = 0..20000</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PrimeConstant.html">Prime Constant.</a>

%H G. Xiao, <a href="http://wims.unice.fr/~wims/en_tool~number~contfrac.en.html">Contfrac</a>

%H <a href="/index/Con#confC">Index entries for continued fractions for constants</a>

%e 0.414682509851111660248109622... = 0 + 1/(2 + 1/(2 + 1/(2 + 1/(3 + ...)))). - _Harry J. Smith_, Jun 15 2009

%t ContinuedFraction[ FromDigits[{{Table[ If[ PrimeQ[n], 1, 0], {n, 370}]}, 0}, 2], 95] (* _Robert G. Wilson v_, Jan 15 2005 *)

%o (PARI) { default(realprecision, 20080); x=0; m=67000; for (n=1, m, if (isprime(n), a=1, a=0); x=2*x+a; ); x=contfrac(x/2^m); for (n=0, 20000, write("b051007.txt", n, " ", x[n+1])); } \\ _Harry J. Smith_, Jun 15 2009

%Y Cf. A010051, A051006. Increasing partial quotients are in A102878.

%K nonn,cofr

%O 0,2

%A _Eric W. Weisstein_