# hw14pi - Pablo Blanco
# NOT OK to post


JesusG:=proc(N) local i,j:
 13-add((mul(2*j+1,j=0..i-1)/mul(2*j,j=1..i))^5 * (20*i*(50+41*(i-1)) +13)/(2^(10*i))  ,i=1..N)
end:


# Comparisons:
# JesusG(10):  0.00005167 from true value
# JesusG(100): 0.00005167 ... from true value ?


# opting out of this one
PiComp:=proc(N):
 evalf(sqrt(128/JesusG(N)),10*N^2):
end:


# Remark: Since I'm sure Jesus's computation should improve... I assume it exceeds Maple rounding limits?