|
MATHEMATICA
|
q0 = 1/16^(n/2); t0 = 4/(8*n + 1); x0 = -2/(8*n + 4); y0 = -1/(8*n + 5); z0 = -1/(8*n + 6); f[n_] = q0^2*(x0 + y0 + z0 + t0) g[n_, 1] := If[Mod[Prime[n], 12] - 1 == 0, f[n], 0] g[n_, 2] := If[Mod[Prime[n], 12] - 5 == 0, f[n], 0] g[n_, 3] := If[Mod[Prime[n], 12] - 7 == 0, f[n], 0] g[n_, 4] := If[Mod[Prime[n], 12] - 11 == 0, f[n], 0] PiN[m_] := Sum[g[n, m], {n, 1, 10000}] Pi5 = Sum[f[n], {n, 0, 2}] aout = Flatten[{1/N[Pi5, 100], Table[1/N[PiN[n], 100], {n, 1, 4}]}] Floor[aout ]
|