Analytic continuation of harmonic sums

We present a method for calculating any (nested) harmonic sum to arbitrary accuracy for all complex values of the argument. The method utilizes the relation between harmonic sums and (derivatives of) Hurwitz zeta functions, which allows a harmonic sum to be calculated as an expansion valid for large values of its argument. A program for implementing this method is also provided. (C) 2009 Elsevier B.V. All rights reserved.