# Aurora Hiveley, 3/10/25, Assignment 14pi

Help := proc(): print(`JesusG(N)`): end:
with(combinat):

Digits := 100:

JesusG := proc(N) local n:
evalf(add( (-1)^n * binomial(2*n,n)^5 * (820*n^2 + 180*n + 13)/(2^(20*n)) , n=0..N)):
end:

# abs(JesusG(10) - 128/Pi^2);
# output: 1E-32
# abs(JesusG(100) - 128/Pi^2);
# output: 0.