Stealing pages from the server...

Simpson Rule for Definite Integrals


Introduction

Simpson’s rules are numerous approximations for definite integrals in numerical analysis, named after English mathematician Thomas Simpson (1710−1761). In calculus, basically, there are two ways to approximate the value of an integral, Reimann sums and Trapezoidal sums. However, calculating the value of an integral, we need to compute the areas of a zillion rectangles or more to get a better result. Therefore, we use Simpson’s Rule, which is a way to approximate integrals without having to deal with lots of narrow rectangles.

Simpson’s 1/3 Rule

The most basic of these rules, called Simpson’s 1/3 rule, or just Simpson’s rule, reads

Introducing the step size this is also commonly written as

Because of the factor Simpson’s rule is also referred to as Simpson’s 1/3 rule.

Simpson’s 3/8 rule

Thomas Simpson proposed Simpson’s 3/8 rule, often known as Simpson’s second rule, as another approach for numerical integration. Rather than a quadratic interpolation, it uses a cubic interpolation. The 3/8 rule of Simpson is as follows:

Numerical Analysis

To obtain an approximation of the definite integral using Simpson’s Rule, we partition the interval into an even number n of subintervals, each of width is . If the function is continuous on , then

The coefficients in Simpson’s Rule have the following pattern: with points.

Example

The question is to use Simpson’s Rule with to approximate the integral .

It is easy to see that the width of each subinterval is and the endpoints . Calculate the function values at the points , which is , , , , .

Substitute all these values into the Simpson’s Rule formula:

The true solution for the integral is

Hence, the error in approximating the integral is

Conclusion

Simpson’s rule is a more accurate form of numerical integration than the Trapezoidal rule, and it should always be used before trying anything more complicated.

References

  1. https://web.stanford.edu/group/sisl/k12/optimization/MO-unit4-pdfs/4.2simpsonintegrals.pdf
  2. https://www.math24.net/simpsons-rule
  3. https://en.wikipedia.org/wiki/Simpson%27s_rule

Author: Yang Wang
Reprint policy: All articles in this blog are used except for special statements CC BY 4.0 reprint polocy. If reproduced, please indicate source Yang Wang !
  TOC