next up previous
Next: Alternative Derivations of Euler's Up: Numerical O.D.E.'s Previous: Euler's Method

Convergence and Error Analysis of Euler's Method

  When iterating Euler's formula, equation (5.11), we obtain a sequence of values tex2html_wrap_inline3688 all of which are built upon the previous calculated values. In general, it is unlikely that tex2html_wrap_inline3644 and tex2html_wrap_inline3692 would completely agree, because Euler's method produces only an approximation to the true solution. Thus, in computing tex2html_wrap_inline3648 from the already slightly-off point tex2html_wrap_inline3644 , we're only compounding the effects of errors.

    The error that occurs in a single step only as a result of the numerical formula is referred to as the local truncation error or the discretization error.

    For Euler's method, this truncation error can be examined using Taylor's theorem. Using Taylor's theorem, write

  equation1017

for some tex2html_wrap_inline3698 . Assuming that y(t) satisfies the differential equation, y'(t) = f(t,y(t)) so that in fact we have

displaymath3680

for some tex2html_wrap_inline3698 . Since we're just addressing the truncation error, which is introduced in a single step, we assert tex2html_wrap_inline3706 , or

displaymath3681

    The term

  equation1030

is the truncation error associated with a single step of Euler's method.

    To analyze the global truncation error of Euler's method, also referred to as the cumulative error, subtract

displaymath3708

from equation (5.12), which yields

  equation1039

    Equation (5.14) shows that the error in computing tex2html_wrap_inline3710 consists of two parts: (1) the truncation error tex2html_wrap_inline3712 which is incurred in stepping from tex2html_wrap_inline3714 to tex2html_wrap_inline3716 , and (2) the contribution from the propagated error

  equation1049
    To further the analysis, the term inside of the brackets above can be viewed as simply a function of the single variable found in the second argument. I.e. let
displaymath3682

Then, the Mean Value Theorem can be called upon to produce the result

eqnarray1052

where tex2html_wrap_inline3718 belongs to the interval bounded by tex2html_wrap_inline3720 and tex2html_wrap_inline3620 . Placing this result in Equation (5.14) yields

  eqnarray1057

for some tex2html_wrap_inline3724 and some tex2html_wrap_inline3718 in the interval bounded by tex2html_wrap_inline3720 and tex2html_wrap_inline3620 .

    This result can be used to rigorously produce a bound on the global error. Although the derivation is omitted, the result is summarized in the following theorem:

  theorem1069

Note the first term on the right-hand side of this expression deals with the error in the initial value. If, in the ideal situation, the initial value used for Euler's method is equal to the true initial value, we have

  corollary1081

    Note that a strong generalization can be made now between the global error and the truncation error. As shown by equation (5.13), the truncation error is tex2html_wrap_inline3738 . In contrast, the global error of Euler's method as seen from the above corollary is tex2html_wrap_inline3740 .

example1087


next up previous
Next: Alternative Derivations of Euler's Up: Numerical O.D.E.'s Previous: Euler's Method
Paul Gray

Wed Oct 28 11:42:13 EST 1998