---

def hurwitz_zeta(z, q, epsilon):
    total_sum = 0.0
    n = 0
    while True:
        term = 1 / ((q + n) ** z)
        if term < epsilon:
            break
        total_sum += term
        n += 1
    return total_sum

# Пример использования:
# z_val, q_val, eps = 2.0, 1.0, 1e-6
# result = hurwitz_zeta(z_val, q_val, eps)
# print(f"Результат: {result}")