Estimating Error w/ Trap and Simpson's Rule when Values are Equal

  • Thread starter Thread starter The Head
  • Start date Start date
  • Tags Tags
    Error
Click For Summary
SUMMARY

The discussion centers on estimating errors using Simpson's Rule and the Trapezoidal Rule for periodic functions. Participants confirm that the calculations for the area under the curve yield consistent results, specifically pi/2(y_0 + y_1 - y_2 + y_3). The conversation highlights the importance of understanding the signs of function values, particularly y2, and its implications on error estimation. The participants conclude that the error for the Trapezoidal estimate is less than that calculated by Simpson's Rule, emphasizing the need for clarity on error bounds in these methods.

PREREQUISITES
  • Understanding of Simpson's Rule and Trapezoidal Rule for numerical integration
  • Knowledge of error analysis in numerical methods
  • Familiarity with derivatives, specifically second and fourth derivatives
  • Ability to interpret graphical representations of functions
NEXT STEPS
  • Research error bounds for the Trapezoidal Rule and Simpson's Rule
  • Study the implications of the Intermediate Value Theorem in numerical integration
  • Explore the relationship between function derivatives and error estimates in numerical methods
  • Examine case studies where Simpson's Rule and Trapezoidal Rule are applied to periodic functions
USEFUL FOR

Students and educators in calculus, mathematicians focusing on numerical methods, and anyone involved in error analysis of integration techniques.

The Head
Messages
137
Reaction score
2
Homework Statement
I'm attaching the full problem with graph, but the problem has me find estimates of a periodic function using its graph by finding the area of a full period by using:
a) Simpson's Rule, starting at two different locations in the periodic function
b) Finding the average of these two estimates
c) Using the Trapezoidal Rule to estimate the area under one period
These are all very straight forward.

Then, part (d) has me suspicious. It tells me |max(f''(x))| AND |max(f''''(x))| on the interval =1 and asks me for the BEST estimate of |Trap estimate - integral (f(x))|
Relevant Equations
E<(b-a)^3/12n^2|max(f''(x))|
E<(b-a)^5/180n^4|max(f''''(x)|
I'm confident in calculating the beginning values. Only question I had was whether in this graph y2 is assumed positive or negative. Obviously f(pi)= negative value, but is y2 negative or is it = |f(pi)|. Assuming the latter, it all makes sense and when I sum the two Simpson's Rule estimates I get:
pi/2(y_0 + y_1 - y_2 +y_3)

I get the exact same thing for the trapezoidal rule. All of this makes sense. But the wording of "best estimate" in part (d) of |Trap estimate - integral (f(x))| has me suspicious that it's not just the error formula.

Further, I figure that since an average of the Simpson's Errors gives me a smaller value, the error for the trapezoidal estimate (since it's the same as the average of Simpson's) is less than this as well. This feels like a bit of an assumption, but I guess I can say:

T= 1/2(S1+S2), so there exists some estimate S3 that equals T (by intermediate value theorem). Thus, the error < that calculated by the Simpson's Rule estimate. Is that about as far as I can go?
 

Attachments

  • Screen Shot 2020-10-21 at 11.54.03 AM.png
    Screen Shot 2020-10-21 at 11.54.03 AM.png
    318.7 KB · Views: 189
Last edited:
Physics news on Phys.org
The Head said:
Homework Statement:: I'm attaching the full problem with graph, but the problem has me find estimates of a periodic function using its graph by finding the area of a full period by using:
a) Simpson's Rule, starting at two different locations in the periodic function
b) Finding the average of these two estimates
c) Using the Trapezoidal Rule to estimate the area under one period
These are all very straight forward.

Then, part (d) has me suspicious. It tells me |max(f''(x))| AND |max(f''''(x))| on the interval =1 and asks me for the BEST estimate of |Trap estimate - integral (f(x))|
Relevant Equations:: E<(b-a)^3/12n^2|max(f''(x))|
E<(b-a)^5/180n^4|max(f''''(x)|

I'm confident in calculating the beginning values. Only question I had was whether in this graph y2 is assumed positive or negative.
Obviously f(pi)= negative value, but is y2 negative or is it = |f(pi)|.
It should be obvious from the graph that y2 is negative.
The Head said:
Assuming the latter, it all makes sense and when I sum the two Simpson's Rule estimates I get:
pi/2(y_0 + y_1 - y_2 +y_3)

I get the exact same thing for the trapezoidal rule. All of this makes sense. But the wording of "best estimate" in part (d) of |Trap estimate - integral (f(x))| has me suspicious that it's not just the error formula.

Further, I figure that since an average of the Simpson's Errors gives me a smaller value, the error for the trapezoidal estimate (since it's the same as the average of Simpson's) is less than this as well. This feels like a bit of an assumption, but I guess I can say:

T= 1/2(S1+S2), so there exists some estimate S3 that equals T (by intermediate value theorem). Thus, the error < that calculated by the Simpson's Rule estimate. Is that about as far as I can go?
I think that what the problem is hinting at is a formula for the error in using the Trapezoid Rule. Has your textbook presented any discussion about the error bounds for this method? Here's some more information: https://en.wikipedia.org/wiki/Trapezoidal_rule, in the section titled Error analysis.
 
  • Like
Likes   Reactions: The Head
Mark44 said:
It should be obvious from the graph that y2 is negative.
I think that what the problem is hinting at is a formula for the error in using the Trapezoid Rule. Has your textbook presented any discussion about the error bounds for this method? Here's some more information: https://en.wikipedia.org/wiki/Trapezoidal_rule, in the section titled Error analysis.
But isn't y_2 showing a length? And why would they spend so much time going through Simpson's if that isn't applied here? It seems fishy they would give f''''(x), point to the idea that the two estimates are the same, and ask for "best estimate." Plus since the average of Simpson = the Trapezoidal, if they have the same estimate, then wouldn't both be within the same error rate of Simpson (which is smaller than the Trapezoidal).
 
The Head said:
But isn't y_2 showing a length?
##y_0, y_1, y_2,## and ##y_3## are coordinates, so can be positive or negative or zero. If your formula is using a positive value for ##y_2## your result is wrong.
The Head said:
And why would they spend so much time going through Simpson's if that isn't applied here? It seems fishy they would give f''''(x), point to the idea that the two estimates are the same, and ask for "best estimate." Plus since the average of Simpson = the Trapezoidal, if they have the same estimate, then wouldn't both be within the same error rate of Simpson (which is smaller than the Trapezoidal).
Simpson's Rule and the Trapezoid Rule have different error terms, one involving the 2nd derivative and the other involving the fourth derivative. For the particular function in this problem, it could well be that one technique will give a smaller error than the other, given how few subintervals are being used.
 

Similar threads

Replies
2
Views
6K
  • · Replies 1 ·
Replies
1
Views
7K
Replies
1
Views
3K
  • · Replies 5 ·
Replies
5
Views
3K
Replies
6
Views
2K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 1 ·
Replies
1
Views
5K
  • · Replies 6 ·
Replies
6
Views
6K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K