Computes the Hurwitz zeta function \zeta(s,q)
for s > 1, q > 0.
