“…are the Bernoulli numbers, (1) reduces for r = 0 to the Euler-Maclaurin summation formula (note that L 1 = 1 2 ) or for s = 0 to the Gregory quadrature formula, respectively. In (1), the remainder term E rs n is of order O(h r+s+2 ) and vanishes for polynomials up to degree r + s. "Classical" representations of the remainder are known for r = 0 [4], for r even, s = 0 [1,7], and for r odd, s = 1 [2].…”