|
EXAMPLE
|
a(1) = ceiling(r(1)) = ceiling(1/tau) = ceiling(0.618...) = 3;
a(2) = ceiling(r(2)/(1 - r(1)/1) = 48;
a(3) = ceiling(r(3)/(1 - r(1)/1 - r(2)/2) = 5215.
The first 3 terms of the series r(1)/a(1) + ... + r(n)/a(n) + ... are
0.961..., 0.997..., 0.99999997...
|
|
MATHEMATICA
|
$MaxExtraPrecision = Infinity; z = 16;
r[k_] := N[2/Log[k + 1], 1000]; f[x_, 0] = x;
n[x_, k_] := n[x, k] = Ceiling[r[k]/f[x, k - 1]]
f[x_, k_] := f[x, k] = f[x, k - 1] - r[k]/n[x, k]
x = 1; Table[n[x, k], {k, 1, z}]
N[Sum[r[k]/n[x, k], {k, 1, 18}], 200]
|