The Hurwitz zeta function is defined by \zeta(s,q) = \sum_0^\infty (k+q)^{-s}.
This routine computes the Hurwitz zeta function \zeta(s,q) for s > 1, q > 0.