Simpson's rule
From Conservapedia
This article/section deals with mathematical concepts appropriate for a student in late high school or early university. |
In Calculus, Simpson's Rule is a method of estimating integrals by using parabolas or other higher order polynomial expressions. It is named after its discoverer, Thomas Simpson^{[1]}.
Contents |
Simpson's Rule
Specifically, Simpson's Rule states:
Example 1
Given:
We can use Simpson's Rule to get a quick approximation of the actual value of the definite integral.
- a = 1
- b = 4
- f(x) = − x^{3} + 4x^{2} − x + 1
- f(1) = − 1^{3} + 4(1)^{2} − 1 + 1 = 3
- f(4) = − 4^{3} + 4(4)^{2} − 4 + 1 = − 3
Now we can evaluate the integral:
In the picture on the right the "error area" is the yellow and dark red area. In cases where the approximation does not give a good answer, we increase the number of subintervals.
If we solve this using integration, we find that Simpson's Rule gave the exact value for the definite integral. However, this is rarely the case, as shown in the next example.
Example 2
Given:
Directly solving this integral would require an understanding of partial fractions. If we want to quickly attain an approximation of the integral, we can use Simpson's Rule.
- a = 4
- b = 12
Example 2 |
---|
Explanation |
This shows the parabolic approximation we are using. Note how much more error area that there is compared to the previous example. |
Now we can solve:
As the function given does not resemble a parabola, the answer we get is completely off of the real answer, which is:
This is why we will sometimes split the interval of integration into multiple pieces. When we increase the number of subintervals, the approximation approaches the actual value (see example below). When we do this, we use the Composite Simpson's Rule.
Composite Simpson's Rule
Simpson's Rule can be applied to segments of an integral, effectively making it a Riemann Integral.
If the interval [a,b] is split into N segments, with N as an even number, there is a formula that can be used to solve for the integral.
Note that:
- x_{n} = x + nh
As an extension of Simpson's Rule, we can find that:
The above formula can be rewritten as:
- The formulas above do not work unless N is even.
Example 2 Continued
Example 2 Continued |
---|
Click on the image to get a better view. |
As our estimation was quite off, we should use the Composite Rule to get a better approximation.
For this example, we will use 8 subintervals.
- N = 8
- x_{0} = 4,x_{1} = 5,x_{2} = 6,x_{3} = 7,x_{4} = 8,x_{5} = 9,x_{6} = 10,x_{7} = 11,x_{8} = 12
This new answer is much closer than the original; it is off by about .05. In the picture on the right, note how much less "error" area that there is compared to when we only used one parabola. The graph represents the first parabola we use for Example 2 Continued, evaluating from 4 to 6. By reducing the interval that each individual parabola approximates, our total approximation became more accurate.
If the number of subintervals is increased, it can be shown that the estimated value approaches the real value:
N | Estimate |
---|---|
2 | 9.34179 |
4 | 8.26138 |
8 | 7.99289 |
10 | 7.96783 |
100 | 7.94161 |
1000 | 7.94160 |
After we reach 106 subintervals, the approximation is correct for five digits.
Other Versions of Simpson's Rule
There are a few other rules that Simpson proved. For example, Simpson's 3/8 Rule^{[2]} shows that cubic functions can be used to approximate integrals.