Contents

Our approach to computing numerical values of these sums involves the compound application of the Euler-Maclaurin summation formula (see [2, p. 289,]), which can be stated as follows. Suppose has at least

We will briefly present a method for computing . See [4] for more details. Let and . By the Euler-Maclaurin summation formula, Since for all

Let **c** be a large integer, and let . Applying the Euler-Maclaurin summation formula (1) again,
we can write

This formula suggests the following computational scheme. First, explicitly evaluate the sum for , using a numeric working precision of 150 digits. Secondly, perform the symbolic integration and differentiation steps indicated in formula (4). Finally, evaluate the resulting expression, again using a working precision of 150 digits. The final result should be equal to to approximately 135 significant digits.

The difficulty and cost of performing the symbolic integration and differentiation operations indicated in (4) can be greatly reduced by approximating as follows: first, expand , the numerator of , into a sum of individual terms; next, write as ; next, expand using the binomial theorem to 18 terms; next, multiply together the resulting numerator and denominator expressions; finally, omit all terms whose exponent of is greater than 18. The result is a linear sum of terms of the form for modest-sized integers

We have performed many computations of this type. The integration and differentiation operations required in (4) can be handled using a symbolic mathematics package, such as Maple [11] or Mathematica [17]. The explicit summation of the first

Whatever software is used, this explicit summation is an expensive operation. For example, the evaluation of to terms, using the MPFUN package with 150-digit precision arithmetic, requires 20 hours on a ``Crimson'' workstation manufactured by Silicon Graphics, Inc. Thus while such runs can be made, clearly this is pressing the limits of current workstation technology. Fortunately, it is possible to perform such computations on a highly parallel computer system. The details of this parallel algorithm are given in [4].

Contents