""" Berechnung von Zete(7) (c) Jürgen Meier www.3d-meier.de 05.07.2024 """ import math p = 0 k = 1 for n in range(1, k + 1): p = p + 1.0/((n**7)*(math.exp(2*math.pi*n)-1)) print(19.0*math.pi**7/56700 - 2*p)