MAT360 Section Summary: 4.6

Adaptive Quadrature

  1. Summary

    Objective: to approximate

    displaymath216

    to within tex2html_wrap_inline242 .

    The trick here is the usual one, of balancing errors, and trying for a better approximation by applying several methods and then taking an appropriate combination of them.

    For example, we might try a Simpson's rule on a single panel over [a,b] using a step-size of h = b-a, and then try Simpson's composite with h/2:

    displaymath217

    where tex2html_wrap_inline250 , and the error of the elemental Simpson's rule is

    displaymath218

    We now consider the error of the h/2 method, but used as a composite rule so that it will span [a,b]:

    displaymath219

    where tex2html_wrap_inline256 and

    displaymath220

    Now, provided tex2html_wrap_inline258 doesn't vary wildly, we can hope that tex2html_wrap_inline260 , so that the error of the h/2 composite method will be about a sixteenth of the error of the elemental rule:

    displaymath221

    Since the methods are both approximations to the same quantity, we can try to combine them to make a better approximation:

    displaymath222

    On the other hand, we can also to check to see if the difference is significant enough to justify using the smaller step size; if not, we can stick with the larger step size - maybe even make it larger! That is, we can adapt to realities ``on the ground'' (or on the interval, at any rate!).

    So

    displaymath223

    If we cavalierly assume that the approximation is exact, then substituting for I we find that

    displaymath224

    or

    displaymath225

    Hence

    displaymath226

    We can easily measure the quantity on the RHS: if

    displaymath227

    then we conclude that

    displaymath228

    If this condition is satisfied, then tex2html_wrap_inline266 is sufficiently close: otherwise, we divide and conquer: split the error in half ( tex2html_wrap_inline268 ), and give one half to each of the two panels of the h/2 method; iterate until satisfied.




Fri Nov 11 02:29:50 EST 2005