|
%I #11 Aug 11 2020 10:00:01
%S 0,1,3,7,3,4,5,6,5,3,2,0,0,5,4,9,1,9,2,4,3,1,4,7,1,7,5,4,0,8,6,2,5,4,
%T 5,1,6,5,0,4,4,6,8,7,7,0,2,3,0,8,3,6,5,6,6,2,9,7,2,5,2,5,5,7,0,1,5,3,
%U 1,9,1,1,3,3,8,5,0,0,8,1,0,9,3,2,2,2,8,9,1,3,7,0,7,3,9,7,9,3
%N Decimal expansion of the convergent to the sum of (1/9)^p where p ranges over the set of prime numbers.
%F Equals 8 * Sum_{k>=1} pi(k)/9^(k+1), where pi(k) = A000720(k). - _Amiram Eldar_, Aug 11 2020
%e 0.0137345...
%o (PARI) /* Sum of 1/m^p for primes p */ sumnp(n,m) = { local(s=0,a,j); for(x=1,n, s+=1./m^prime(x); ); a=Vec(Str(s)); for(j=3,n, print1(eval(a[j])",") ) }
%Y Cf. A000720, A132799 (base 8), A010051 (base 10).
%K cons,nonn
%O 0,3
%A _Cino Hilliard_, Nov 17 2007
%E Offset corrected _R. J. Mathar_, Jan 26 2009
|
