OFFSET
1,2
COMMENTS
The first eight terms give Pi accurately to the first 75 decimal digits.
FORMULA
a(1) = 1, a(2) = 2, a(3) = 3, a(4) = 4; a(n) = ceiling((Pi - Sum_{i=1..n-1}(1/sqrt(a(i))))^-2).
EXAMPLE
Pi > 1/sqrt(1) + 1/sqrt(2) + 1/sqrt(3) + 1/sqrt(4) + 1/sqrt(8), but Pi < 1/sqrt(1) + 1/sqrt(2) + 1/sqrt(3) + 1/sqrt(4) + 1/sqrt(7).
MATHEMATICA
a[1] = 1; a[2] = 2; a[3] = 3; a[4] = 4; a[n_] := Ceiling[(Pi - Sum[1/Sqrt[a[i]], {i, 1, n - 1}])^-2]; Table[ a[n], {n, 1, 9}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Feb 13 2004
STATUS
approved