[Home]Euler-Maclaurin formula

HomePage | Recent Changes | Preferences

Difference (from prior major revision) (author diff)

Changed: 7c7
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

Changed: 10c10
1
0

Changed: 12,13c12,13
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

Changed: 15c15
S - I = ∑ ---- (f(2k-1)(n) - f(2k-1)(1)) + R
S - I = ∑ ---- (f(2k-1)(n) - f(2k-1)(0)) + R

Changed: 22c22
(2π)2p 1
(2π)2p 0

Changed: 28c28
n-1
n

Changed: 30c30
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.

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
          0
can 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 f(k) at the end points of the interval 0 and n. For any natural number p, we have
               p   B2k      
    S - I  =  ∑   ---- (f(2k-1)(n) - f(2k-1)(0))   +   R
              k=1 (2k)!                
Here, B2 = 1/6, B4 = -1/30, B6 = 1/42, B8 = -1/30 ... are the [Bernoulli numbers]?.

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) = x3, we can choose p = 2 to obtain after simplification

     n 
    ∑  i3  =  1/4 n4 - 1/2 n3 + 1/4 n3.
    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:


HomePage | Recent Changes | Preferences
This page is read-only | View other revisions
Last edited August 10, 2001 9:35 pm by AxelBoldt (diff)
Search: