|
MATHEMATICA
|
$MaxExtraPrecision = Infinity; z = 9; f[n_] := 1/(2 n + 1); g[n_] := (1/n - 1)/2; r = 2; s = 0; a[1] = NestWhile[# + 1 &, 1, ! (s += f[#]) > r &]; p = Sum[f[n], {n, 1, a[1]}]; a[2] = Floor[g[p - r]]; a[n_] := Floor[g[((-1)^n) (p - r - Sum[((-1)^k) f[a[k]], {k, 2, n - 1}])]];
Table[a[k], {k, 1, z}]
|