cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A278389 Decimal expansion of Sum_{k>=1} (-1)^(k+1)/(k*prime(k)).

Original entry on oeis.org

3, 7, 4, 4, 8, 5, 1, 8, 7, 9, 7, 4, 7, 4, 6, 1, 6, 3, 2, 1, 7, 0, 9, 4, 0, 8, 6
Offset: 0

Views

Author

Jon E. Schoenfield, Nov 20 2016

Keywords

Comments

Alternating sum of the reciprocals of the products of k and the k-th prime.
From Jon E. Schoenfield, Jan 15 2021: (Start)
The second Mathematica program appears to compute partial sums through k=10^10 and includes a comment that it is "good for the first 27 digits". For the (10^10-1)st, (10^10)th and (10^10+1)st partial sums, I get
.
k prime(k) S(k) = k-th partial sum
=========== ============ ====================================
9999999999 252097800611 0.3744851879747461632172924014592...
10000000000 252097800623 0.3744851879747461632168957300096...
10000000001 252097800629 0.3744851879747461632172924014591...
.
Since the sum is alternating, successive partial sums oscillate around a central curve, and taking the mean of the k-th and (k+1)st partial sums gives an estimate of the value around which the sum appears to be converging. Those means for k = 9999999999 and 10000000000 are
.
k (S(k) + S(k+1))/2
=========== =====================================
9999999999 0.37448518797474616321709406573443...
10000000000 0.37448518797474616321709406573440...
.
Of course, even though these results happen to agree through their first 31 significant digits, they are certainly nowhere near sufficient to establish the first 31 significant digits of the infinite sum (since the sequence of mean partial sums tends to meander). However, these results are clearly lower than the value given in the Data (i.e., 0.374485187974746163217094086), so (if correct) they would appear to indicate either that the last two terms in the Data are incorrect or that those terms were obtained using partial sums significantly beyond the (10^10)th. (End)

Examples

			0.374485187974746163217094086...
		

Crossrefs

See the following for alternating sums of reciprocals of primes, composites, and related expressions: A078437 (primes), A242301 (primes^2), A242302 (primes^3), A242303 (primes^4), A242304 (primes^5), A269229 (composites), A275110 (composites excluding prime powers), A275712 (nonprimes), A276494 (composites^2).

Programs

  • Mathematica
    RealDigits[N[-Sum[(-1)^k/(k*Prime[k]), {k, 1, 8*10^6}], 30]][[1]] (* G. C. Greubel, Nov 22 2016 *)
    s = 0; k = 1; p = 2; While[k =< 10^10, s = N[s - (-1)^k/(k*p), 48]; k++; p = NextPrime@p]; RealDigits[s, 10, 20] (* good for the first 27 digits *) (* Robert G. Wilson v, Mar 07 2019 *)

Formula

Sum_{k>=1} (-1)^(k+1)/(k*prime(k)) = 1/(1*2) - 1/(2*3) + 1/(3*5) - 1/(4*7) + 1/(5*11) - ... .

Extensions

Edited and a(21)-a(26) from Robert G. Wilson v, Mar 07 2019

A285308 Decimal expansion of the sum of the alternating series of reciprocals of cubed composite numbers.

Original entry on oeis.org

0, 1, 2, 2, 1, 0, 3, 4, 4, 6, 8, 7, 8, 2, 2, 3
Offset: 0

Views

Author

Terry D. Grant, Apr 16 2017

Keywords

Examples

			Equals 1/(4^3) - 1/(6^3) + 1/(8^3) - 1/(9^3) + ... =
0.0122103446878223...
		

Crossrefs

Programs

  • PARI
    composite(n) = {
    local(c,x);
    c=1;
    x=1;
    while(c <= n,
    x++;
    if(!isprime(x),c++);
    );
    return(x)
    }
    \p 35000
    sumalt(k=1, (-1)^(k+1)/composite(k)^3)
    
  • PARI
    upto(n) = my(c=4,s=1,t=0.); while(c<=n, t+=s/c^3; c++; while(isprime(c), c++); s=-s); t
    upto(n) = my(prev=3,s=1,t=0.); forprime(p=5,n,for(i=prev+1,p-1,t+=s/i^3;s=-s;prev=p));s=-s; t+sum(i=prev+1,n,s=-s;s/i^3) \\ David A. Corneth, Apr 29 2017

Formula

Equals Sum_{k>=1} (-1)^(k+1)/(A002808(k)^3).
Showing 1-2 of 2 results.