next up previous
Next: Convergence and Error Analysis Up: Numerical O.D.E.'s Previous: Notation

Euler's Method

  Euler's method is a very simple numerical method for solving the first-order initial value problem (5.1). Euler's method can be intuitively derived from several different avenues and, for this reason, is often considered to be the cornerstone method on which the framework to formulate and discuss subsequent methods is built.

    To derive Euler's method, recall how the derivative is computed. By definition, the derivative of tex2html_wrap_inline3620 is

displaymath3622

From the differential equation, this limit is functionally given as f(t,y), i.e.

displaymath3626

    One way of deriving Euler's method comes about by considering the quotient in this limit to be a good approximation when h is small,

eqnarray948

    Euler's method is derived from the above equation by asserting equality in this approximation and replacing h with the step size tex2html_wrap_inline3614 , y with the approximation tex2html_wrap_inline3636 .

displaymath3638

Rearranging terms yields the recursion relationship

  equation955

which is Euler's method. Equation (5.11) builds the approximation at each step using the previous values of time and tex2html_wrap_inline3636 . For the initial approximation, tex2html_wrap_inline3642 is used or some close approximation. With this initialization, tex2html_wrap_inline3644 can be computed, then tex2html_wrap_inline3644 may be used to compute tex2html_wrap_inline3648 , and so on.

   example960

    The next section analyzes the error associated with Euler's method. As a motivation for the aspects of error analysis, the preceding example is presented with a smaller step size in t. Note the differences in the error in the computations of Example 5.4 where tex2html_wrap_inline2330 and Example 5.5 where this step size is halved.

   example988
 


Paul Gray

Wed Oct 28 11:42:13 EST 1998