login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A006524
Egyptian fraction for 1/ Pi.
(Formerly M3509)
26
4, 15, 609, 845029, 1010073215739, 1300459886313272270974271, 1939680952094609786557359582286462958434022504402
OFFSET
1,1
REFERENCES
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
MATHEMATICA
a = {}; k = N[1/Pi, 1000]; Do[s = Ceiling[1/k]; AppendTo[a, s]; k = k - 1/s, {n, 1, 10}]; a (* Artur Jasinski, Sep 22 2008 *)
PROG
(PARI) x = 1/Pi;
f(x, k) = if(k<1, x, f(x, k - 1) - 1/n(x, k));
n(x, k) = ceil(1/f(x, k - 1));
for(k = 1, 7, print1(n(x, k), ", ")) \\ Indranil Ghosh, Mar 27 2017
CROSSREFS
Sequence in context: A000881 A354462 A109923 * A299683 A341598 A295231
KEYWORD
nonn,frac
EXTENSIONS
More terms from Herman P. Robinson
Offset changed to 1 by Indranil Ghosh, Mar 27 2017
STATUS
approved