12 lines
367 B
Text
12 lines
367 B
Text
pi[n_, prec_] :=
|
|
Module[{a = 1, g = N[1/Sqrt[2], prec], k, s = 0, p = 4},
|
|
For[k = 1, k < n, k++,
|
|
{a, g} = {N[(a + g)/2, prec], N[Sqrt[a g], prec]};
|
|
s += p (a^2 - g^2); p += p]; N[4 a^2/(1 - s), prec]]
|
|
|
|
|
|
pi[7, 100] - N[Pi, 100]
|
|
1.2026886537*10^-86
|
|
|
|
pi[7, 100]
|
|
3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628046852228654
|