%I #37 Jan 21 2019 07:45:11
%S 8,4,8,9,6,9,0,3,4,0,4,3
%N Decimal expansion of Sum_{k>=1} 1/(k*prime(k)).
%C From _Robert Price_, Jul 14 2010: (Start)
%C This series converges very slowly. I could not find any transform that converges faster, so I did this by brute force using 256 bits of precision.
%C After k=596765000000 terms (p(k)=17581469834441) the partial sum is 0.848 969 034 043 245 206 069 544 346 415 327 714...
%C The next two digits are either 29 or 30. (End)
%C The table in the Example section shows, for increasing values of j, the results of computing the partial sum s(j) = Sum_{k=1..j} 1/(k*prime(k)) and adding to it an approximate value for the tail (i.e., the sum for all the terms k > j). See the Links entry for an explanation of the method used in approximating the size of the tail of the summation beyond the j-th prime. - _Jon E. Schoenfield_, Jan 20 2019
%H Jon E. Schoenfield, <a href="/A124012/a124012.txt">Notes on approximating the size of the summation's tail beyond the j-th prime</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PrimeNumberTheorem.html">Prime Number Theorem</a>
%e 0.848969034043...
%e From _Jon E. Schoenfield_, Jan 14 2019: (Start)
%e We can obtain prime(2^d) for d = 0..57 from the b-file for A033844. Given the above result from _Robert Price_, and letting j_RP = 596765000000, the partial sum through
%e prime(j_RP) = 17581469834441
%e is
%e s(j_RP) = Sum_{k=1..j_RP} 1/(k*prime(k))
%e = 0.848969034043245206069544346415327714...;
%e adding to this actual partial sum s(j_RP) the approximate tail value
%e t(j_RP) =
%e h'(prime(j_RP), prime(2^40))
%e + (Sum_{d=41..57} h'(prime(2^(d-1)), prime(2^d)))
%e + lim_{x->infinity} h(prime(2^57), x)
%e (see the Links entry for an explanation) gives the result 0.84896903404330021273712255895762255... (which seems likely to be correct to at least 20 significant digits).
%e The table below gives, for j = 2^16, 2^17, ..., 2^32, and j_RP, the actual partial sum s(j) and the sum s(j) + t(j) where t(j) is the approximate tail value beyond prime(j).
%e .
%e j s(j) s(j) + t(j)
%e ==== ====================== ======================
%e 2^16 0.84896790758922908159 0.84896903393397518971
%e 2^17 0.84896850050492294891 0.84896903400552099072
%e 2^18 0.84896878057566843770 0.84896903404214147367
%e 2^19 0.84896891330602605081 0.84896903404317536927
%e 2^20 0.84896897639243509768 0.84896903404350431035
%e 2^21 0.84896900645590169648 0.84896903404376063663
%e 2^22 0.84896902081581006534 0.84896903404343742139
%e 2^23 0.84896902768965496764 0.84896903404337393698
%e 2^24 0.84896903098637626311 0.84896903404331189996
%e 2^25 0.84896903257029535468 0.84896903404329806633
%e 2^26 0.84896903333252861584 0.84896903404330030271
%e 2^27 0.84896903369988697984 0.84896903404330084536
%e 2^28 0.84896903387717904236 0.84896903404330042023
%e 2^29 0.84896903396285181513 0.84896903404330024036
%e 2^30 0.84896903400430044877 0.84896903404330021861
%e 2^31 0.84896903402437548991 0.84896903404330021472
%e 2^32 0.84896903403410856545 0.84896903404330021655
%e ... ... ...
%e j_RP 0.84896903404324520607 0.84896903404330021274
%e (End)
%Y Cf. A033286, A085548, A209329, A210473.
%K cons,more,nonn
%O 0,1
%A _Pierre CAMI_, Nov 02 2006
%E Offset and leading zero corrected by _R. J. Mathar_, Jan 31 2009
%E Four more terms (4,0,4,3) from _Robert Price_, Jul 14 2010
%E Title and example edited by _M. F. Hasler_, Jan 13 2015