Acceleration of generalized hypergeometric functions through precise remainder asymptotics, preprint arxiv.org/1102.3003

Abstract. We express the asymptotics of the remainders of the partial sums {sn} of the generalized hypergeometric function q+1Fq α1,...,αq+1 β1,...,βq ∣∣ ∣ z) through an inverse power series znnλ ∑ ck nk, where the exponent λ and the asymptotic coefficients {ck} may be recursively computed to any desired order from the hypergeometric parameters and argument. From this we derive a new series acceleration technique that can be applied to any such function, even with complex parameters and at the branch point z = 1. For moderate parameters (up to approximately ten) a C implementation at fixed precision is very effec-tive at computing these functions; for larger parameters an implementation in higher than machine precision would be needed. Even for larger parameters, however, our C implementation is able to correctly determine whether or not it has converged; and when it converges, its estimate of its error is accurate. 1