How does simpsons rule work

  • Thread starter brandy
  • Start date
  • #1
161
0
can anyone tell me how it works, like visually, what does it do. Im curious.
iv heard its a three point approximation whereas the trapazoidal rule is a two point approximation. But im not entirely sure i understand.
 

Answers and Replies

  • #2
HallsofIvy
Science Advisor
Homework Helper
41,833
963
The trapezoid method approximates the curve, between successive points, by a straight line while Simpson's rule uses three points at a time and approximates the curve by the parabola passing through the three points.

Suppose your function is f(x) and you take three points, [itex](x_0, f(x_0), (x_1, f(x_1), (x_2, f(x_2)), equally spaced on the x-axis. For simplicity, take the x-values to be [itex]x_0= x_1- h[/itex], [itex]x_1[/itex], and [itex]x_2= x_1+ h[/itex] where h is the distance between succesive x values so the points themselves are [itex](x_1- h, f(x_0)), (x_1, f(x_1)), (x_1+ h,f(x_2))[/itex].

The equation of a parabola through [itex](x_1, f(x_1))[/itex] can be written [itex]y= a(x- x_1)^2+ b(x- x_1)+ f(x_1)[/itex]. Taking [itex]x= x_0= x_1- h[/itex] we get [itex]y= f(x_0)= ah^2+ bh+ f(x_1)[/itex]. Taking [itex]x= x_2= x_1+ h[/itex] we get [itex]y= f(x_2)= ah^2- bh+ f(x_1)[/itex].

Adding those two equations, [itex]f(x_2)+ f(x_0)= 2ah^2+ 2f(x_1)[/itex] so that
[tex]ah^2= \frac{f(x_2)+ f(x_0)}{2}- f(x_1}= \frac{f(x_2)- 2f(x_1)+ f(x_0)}{2}[/tex]]
and
[tex]a= \frac{f(x_0)+ 2f(x_1)+ f(x_2)}{2h^2}[/itex]

Subtracting those two equations, [itex]f(x_2)- f(x_0)= 2bh[/itex] so that
[tex]b= \frac{f(x_2)- f(x_0)}{2h}[/tex]

Now, [itex]\int_{x_0}^{x_2}a(x-x_1)^2+ b(x-x_1)+ f(x_1) dx= \int_{x_1-h}^{x_1+h} a(x-x_1)^2+ b(x-x_1)+ f(x_1) dx[/itex][itex]= (1/3)ah^3+ (1/2)bh^2+ f(x_1)h- (1/3)a(-h)^3- b(-h)^2- f(x_1)(-h)= (2/3)ah^3+ 2f(x_1)h[itex] so b really isn't important. The integral of the parabola through [itex](x_1-h, f(x_0)), (x_1, f(x_1)), (x_1+h,f(x_2))[/itex] is
[tex](2/3)ah^3+ 2f(x_1)h= (2/3)\frac{f(x_0)- 2f(x_1)+ f(x_2)}{2h^2}h^3+ 2f(x_1)h[/tex]
[tex]= \frac{2(f(x_0)- 2f(x_1)+ f(x_2)}{3}h+ 2f(x_1)h= h\frac{2f(x_0)- 4f(x_1)+ 2f(x_2)}{3}[/itex]

That is the formula for exactly 3 data points. If we have more, we can separte into groups of three and do that for each such interval. Bacause the overlap at the ends, we must have an odd number of points to do that.

For example, if we had only an even number of points, [itex]x_0, x_1, x_2, x_3[/itex], our first parabola would be over [itex]x_0, x_1, x_2[/itex] but, even including [itex]x_3[/itex] as our first endpoint in the next interval we would have only [itex]x_2, x_3[/itex], not enough for a parabola. But with 5 points, we could have [itex]x_0, x_1, x_2[/itex], [itex]x_2, x_3, x_4[/itex]. Similarly, with 6 points, an even number, we would have [itex]x_0, x_1, x_2[/itex], [itex]x_2, x_3, x_4[/itex], [itex]x_4, x_4[/itex] and not be able to complete the last parabola but with 7, an odd number, we would have [itex]x_0, x_1, x_2[/itex], [itex]x_2, x_3, x_4[/itex], [itex]x_4, x_5, x_6[/itex].

Now see what happens when we "attach" parabolas: on [itex]x_0, x_1, x_2[/itex] we have
[tex]h\frac{2f(x_0)- 4f(x_1)+ 2f(x_2)}{3}[/itex]
while on [itex]x_2, x_3, x_4[/itex]
[tex]h\frac{2f(x_2)- 4f(x_3)+ 2f(x_4)}{3}[/itex]

Those add to give
[tex]h\frac{2f(x_0)- 4f(x_1)+ 4f(x_2)- 4f(x_3)+ 2f(x_4)}{3}[/itex]
so you see why we have "4" multiplying everything except the first and last terms. Taking the integral to be from a to b and using n intervals, h= (b-a)/n so the formula becomes
[tex](b-a)\frac{2f(x_0)- 4f(x_1)+ 4f(x_2)- 4f(x_3)+ 2f(x_4)}{3n}[/itex]
Simpson's rule.
 

Related Threads on How does simpsons rule work

  • Last Post
Replies
1
Views
15K
Replies
31
Views
3K
Replies
1
Views
525
Replies
1
Views
9K
  • Last Post
Replies
1
Views
2K
  • Last Post
Replies
4
Views
3K
  • Last Post
3
Replies
65
Views
8K
  • Last Post
Replies
2
Views
2K
  • Last Post
Replies
6
Views
8K
  • Last Post
Replies
5
Views
793
Top