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 c
terms, as indicated in (4), could be performed by utilizing the
multiple precision facility in the Maple or Mathematica packages.
However, it was found that the MPFUN multiple precision package and
translator developed by one of us [3] was significantly
faster for this purpose.