# 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.