ζ(1) = +∞, при s>1 ряд для ζ(s) сходится, но только достаточно медленно вплоть где-то до s=3. Можно попробовать ускорить сходимость используя тождество Эйлера:
ζ(s) = 1/Π(p∈Primes)(1-1/p^s).
Ещё можно рассматривать ζ-функцию Римана как аналитическую. А Вам что вообще нужно решить? Вот график на всякий случай:

У графика две асимптоты: вертикальная х=1 и горизонтальная у=1.