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.

A065478 Decimal expansion of Product_{p prime} (1 - p/(p^3-1)).

Original entry on oeis.org

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

Views

Author

N. J. A. Sloane, Nov 19 2001

Keywords

Examples

			0.57595996889294543964316337549249669...
		

References

  • Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, p. 106.
  • Paulo Ribenboim, My Numbers, My Friends: Popular Lectures on Number Theory, Springer-Verlag, NY, 2000, p. 16.

Crossrefs

Cf. A078079.

Programs

  • Mathematica
    $MaxExtraPrecision = 100; m0 = 200; dm = 200; digits = 101; Clear[f]; f[m_] := f[m] = (slog = Normal[Series[Log[1 - p/(p^3 - 1)], {p, Infinity, m}]]; Exp[slog] /. Power[p, n_] -> PrimeZetaP[-n] // N[#, digits+10]&); f[m = m0]; Print[m, " ", f[m]]; f[m = m + dm]; While[Print[m, " ", f[m]]; RealDigits[f[m], 10, digits+5] != RealDigits[f[m - dm], 10, digits+5], m = m + dm];RealDigits[f[m], 10, digits] // First (* Jean-François Alcover, Sep 15 2015 *)
  • PARI
    prodeulerrat(1 - p/(p^3-1)) \\ Amiram Eldar, Mar 17 2021