import decimal from decimal import * getcontext().prec = 200 def my_factorial (number): facto = Decimal ('1') for n in range (1, round (number) + 1, 1): facto = Decimal (facto) * Decimal (n) return facto pi = Decimal ('3.14159265358979323846264338327950288419716939937510582097494459230781640628620899862803482534211706798214808651328230664709384460955058223172535940812848111745028410270193852110555964462294895493038196') a = [] s = Decimal ('0') zeta = [] for j in range (1, 101, 1): s = Decimal (s) + Decimal ('1') a.append ((Decimal ('-1') ** (Decimal (s) + Decimal ('1')) * Decimal (s) / my_factorial (Decimal ('2') * Decimal (s) + Decimal ('1')))) for i in range (1, round (s), 1): a[j - 1] = Decimal (a[j - 1]) + Decimal (a[i - 1]) * Decimal ('-1') ** (Decimal (s) + Decimal (i) + Decimal ('1')) / my_factorial (Decimal ('2') * Decimal (s) + Decimal ('1') - Decimal ('2') * Decimal (i)) zeta.append ((Decimal (a[j - 1]) * Decimal (pi) ** (Decimal ('2') * Decimal (s)))) print ("") print ("", 2 * j, "") print ("", '{:102.100f}'.format(zeta[j - 1]), "") print ("") print ("")