1 / 6

Simpson’s Rule

Simpson’s Rule. Nicole Typaldos. Basic Simpson’s Rule:. NOTE: P (x) ≈f (x) then ∫P (x) dx ≈ ∫f (x) dx Using two equal subintervals say, [ a, a+h ] and [ a+h , a+2h] with the interpolating polynomial P(x) of degree two, on the integrated interval. We have: With error term:.

Download Presentation

Simpson’s Rule

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Simpson’s Rule Nicole Typaldos

  2. Basic Simpson’s Rule: • NOTE: P(x)≈f(x) then ∫P(x)dx ≈ ∫f(x)dx • Using two equal subintervals say, [a, a+h] and [a+h, a+2h] with the interpolating polynomial P(x) of degree two, on the integrated interval. • We have: • With error term:

  3. Questions about approximating: • What if the graph oscillates greatly in some subintervals? • What if the oscillation of the interpolating polynomial P(x) is greater than the function we are trying to approximate?

  4. Adaptive Simpson’s Rule: • With the basic Simpson’s Rule those subintervals with large magnitude first derivatives are not very accurate and grossly underestimate the amplitude and number of oscillations of the function. • Let ε be equal to the total error over the entire interval.

  5. Adding more points • Let ε be equal to the total error over the entire interval. • Then with the 1st Simpson’s we have 2 subintervals: • With ≈ errors: ε/2ε/2 • If the error in say “A”≤ ε/2 then that subinterval is left alone and the next is evaluated. • If the error in “B”> ε/2, then “B” is broken up into further subintervals until the total error of “B” ≤ ε/2 . A B

  6. Adaptive Simpson’s Rule: • For example: • We want to approximate a function using Simpson’s Adaptive Rule y=exp(t).*sin(t.*cos(exp(t))) using the error of 1*e^(-5) • To the example using MATLAB.

More Related