If n is a natural number and f(x) is a smooth (meaning: sufficiently often differentiable?) function defined for all real numbers x between 1 and n, then the integral |

If n is a natural number and f(x) is a smooth (meaning: sufficiently often differentiable?) function defined for all real numbers x between 0 and n, then the integral |

1 |

0 |

S = f(1)/2 + f(2) + ... + f(n-1) + f(n)/2 (see [trapezoidal rule]?). The Euler-Maclaurin formula provides expressions for the difference between the sum and the integral in terms of the higher derivatives f^{(k)} at the end points of the interval 1 and n. For any natural number p, we have |

S = f(0)/2 + f(1) + ... + f(n-1) + f(n)/2 (see [trapezoidal rule]?). The Euler-Maclaurin formula provides expressions for the difference between the sum and the integral in terms of the higher derivatives f^{(k)} at the end points of the interval 0 and n. For any natural number p, we have |

S - I = ∑ ---- (f^{(2k-1)}(n) - f^{(2k-1)}(1)) + R |

S - I = ∑ ---- (f^{(2k-1)}(n) - f^{(2k-1)}(0)) + R |

(2π)^{2p} 1 |

(2π)^{2p} 0 |

n-1 |

n |

i=1 |

i=0 |

The **Euler-Maclaurin formula** provides a powerful connection
between integrals (see calculus) and sums. It can be used to approximate
integrals by finite sums, or conversely to evaluate finite sums and infinite series using integrals and the machinery of calculus. The formula was
discovered independently by Leonhard Euler and [Colin MacLaurin]?
around 1735. Euler needed it to compute slowly converging infinite series while MacLaurin? used it to calculate integrals.
*f*^{(k)} at the end points of the interval 0 and *n*. For any natural number *p*, we have
*B*_{2} = 1/6, *B*_{4} = -1/30, *B*_{6} = 1/42, *B*_{8} = -1/30 ... are the [Bernoulli numbers]?.

If *n* is a natural number and *f*(*x*) is a smooth (meaning: sufficiently often differentiable?) function defined for all real numbers *x* between 0 and *n*, then the integral

n I = ∫ f(x) dx 0can be approximated by the sum

S = f(0)/2 + f(1) + ... + f(n-1) + f(n)/2(see [trapezoidal rule]?). The Euler-Maclaurin formula provides expressions for the difference between the sum and the integral in terms of the higher derivatives

p BHere,_{2k}S - I = ∑ ---- (f^{(2k-1)}(n) - f^{(2k-1)}(0)) + R k=1 (2k)!

*R* is an error term which is normally small if *p* is large enough and can be estimated as

2 n |R| ≤ ----- ∫ |f^{(2p+1)}(x)| dx (2π)^{2p}0

By employing the [substitution rule]?, one can adapt this formula also to functions *f* which are defined on some other interval of the real line.

If *f* is a polynomial and *p* is big enough, then the remainder term vanishes. For instance, if *f*(*x*) = *x*^{3}, we can choose *p* = 2 to obtain after simplification

n ∑ i^{3}= 1/4 n^{4}- 1/2 n^{3}+ 1/4 n^{3}. i=0

With the function *f*(*x*) = log(*x*), the Euler-Maclaurin formula can be used to derive precise error estimates for Stirling's formula (see factorial).

**External links:**