Document Type

Conference Proceeding

Publication Date



The Euler-Maclaurin summation formula is frequently used to efficiently estimate sums of infinite series of the form $\sum_{j=1}^{\infty}f(j)$. The purpose of this article is to describe a modification of this numerical technique designed to simplify and reduce the computational effort required to obtain an acceptable estimate of the sum. The modified formula is obtained by replacing $f\left( x\right) $ with an easily constructed polynomial like interpolating function $a\left( x\right) $ designed to simplify the calculation of the integral and derivatives associated with Euler-Maclaurin. This approach provides a more tractable algorithm which can be written as a matrix equation. Examples are provided to demonstrate that the accuracy of the new algorithm compares favorably with that of the traditional formula. The paper concludes with a brief discussion of a method for approximating the error incurred when replacing the exact value of the sum of the original series with the estimate.


To view the content in your browser, please download Adobe Reader or, alternately,
you may Download the file to your hard drive.

NOTE: The latest versions of Adobe Reader do not support viewing PDF files within Firefox on Mac OS and if you are using a modern (Intel) Mac, there is no official plugin for viewing PDF files within the browser window.