|
%I #3 Mar 30 2012 17:34:40
%S 1,1,4,9,7,0,5,0,8,8,1,2,6,6,9,0,2,1,5,7,4,0,6,3,3,4,4,2,8,9,5,9,1,6,
%T 5,9,1,7,3,3,4,7,0,4,1,9,1,2,3,3,9,8,8,2,3,7,0,0,2,3,4,7,2,6,6,8,9,9,
%U 8,5,5,8,6,6,4,4,0,6,1,3,4,8,6,2,1,8,5,5,8,6,5,2,8,2,9,9,8,7,3,8,1
%N Decimal expansion of the product 1/(1 - 1/A000668(n)^2), n= 1 to infinity.
%e equals 1.149705...
%t (*A000668 Mersenne primes *);
%t $MaxExtraPrecision = 200;
%t a = Flatten[Table[If[PrimeQ[2^n - 1], 2^n - 1, {}], {n, 0, 2000}]];
%t c = N[Product[1/(1 - 1/a[[n]]^2), {n, 1, Length[a]}], 100];
%t Table[Mod[Floor[10^n*c], 10], {n, 0, 100}]
%Y Cf. A000668
%K nonn,cons
%O 1,3
%A _Roger L. Bagula_, Apr 17 2010
%E Invalid claim concerning A065443 removed, keyword:cons inserted, offset corrected - The Assoc. Eds. of the OEIS, Apr 23 2010
|
