""" Berechnung von Zeta(4) (c) Jürgen Meier www.3d-meier.de 05.07.2024 """ import math p = 0 k = 15 for n in range(1, k + 1): a = (1.0* math.factorial(n)**2)/((n**4) * math.factorial(2*n)) p = p + a print(36.0*p/17)