We analyze the behavior of Euler-Maclaurin-based integration schemes with the intention of deriving accurate and economic estimations of the error. These schemes typically provide very high-precision results (hundreds or thousands of digits), in reasonable run time, even in cases where the integrand function has a blow-up singularity or infinite derivative at an endpoint. Heretofore, researchers using these schemes have relied mostly on ad hoc error estimation schemes to project the estimated error of the present iteration. In this paper, we seek to develop some more rigorous, yet highly usable schemes to estimate these errors.
David H. Bailey, Jonathan M. Borwein