%I #10 Mar 14 2021 05:20:35
%S 0,1,7,2,3,1,1,1,7,1,1,6,1,5,1,1,3,2,1,1,3,2,1,1,1,1,13,1,16,1,1,16,1,
%T 80,1,1,1,1,7,5,1,4,1,33,3,8,1,8,1,16,11,1,2,6,1,19,1,12,5,11,1,7,5,1,
%U 1,1,2,5,1,4,1,3,4,4,4,1,11,1,2,5,4,12,3,1,4,1,3,1,168,1,4,1,1
%N Continued fraction expansion of Product_{p prime} (1 - 1/(p^2*(p+1))).
%t digits = 93;
%t $MaxExtraPrecision = 4 digits;
%t terms = 4 digits;
%t LR = Join[{0, 0, 0}, LinearRecurrence[{-2, -1, 1, 1}, {-3, 4, -5, 3}, terms + 10]];
%t r[n_Integer] := LR[[n]];
%t c = Exp[NSum[r[n] PrimeZetaP[n - 1]/(n - 1), {n, 4, terms}, NSumTerms -> terms, WorkingPrecision -> digits + 10]];
%t ContinuedFraction[c][[;; digits]] (* _Jean-François Alcover_, Aug 01 2019 *)
%o (PARI) contfrac(prodeulerrat(1 - 1/(p^2*(p+1)))) \\ _Amiram Eldar_, Mar 14 2021
%Y Cf. A065465.
%K nonn,cofr
%O 1,3
%A _Benoit Cloitre_, Dec 02 2002
|