Overview of Chapter Seven

The goals of this brief chapter are to give an introduction to numerical methods other than Euler's method and to discuss error analysis in more detail. In particular, we introduce the concept of the order of a numerical method. Section 7.1 considers the errors involved in Euler's method. Sections 7.2 and 7.3 present improved Euler's method and Runge-Kutta, and Section 7.4 deals with the problem of roundoff errors arising from finite arithmetic.

Since the entire chapter deals almost exclusively with first-order equations, most of this material can be covered immediately after Section 1.4 (Euler's method). However, Sections 7.1 and 7.4 are somewhat more sophisticated than the rest of the text, so it may be appropriate to skim Section 7.1 and jump directly to Section 7.2.

We expect that the students will use some type of computer technology as they solve most of the exercises in this chapter.

DETools

Both EulersMethod and EulersMethodForSystems provide four different step sizes for Euler's method as well as Runge-Kutta with one step size. These tools may help illustrate some of the ideas discussed in this chapter, but alone they are probably not sufficient. Actually it is surprisingly easy to do most of the necessary computation using a spreadsheet.

7.1 Numerical Error in Euler's Method

This section is a detailed discussion of the errors involved in Euler's method including an explanation of the fact that Euler's method is a first-order method. Precise error bounds are derived in Exercises 11 and 12, but the text contains a discussion of the reasons why estimates of this type are often far too conservative.

Comments on selected exercises

Exercises 1-5 study the errors involved in Euler's method as a function of step size.

Exercise 6 is a good exercise to see if the students have a practical understanding of the fact that Euler's method is a first-order method.

Exercises 7-9 use a convenient method to estimate the error in an approximation.

Exercises 10-12 are theoretical. In particular, Exercise 11 is a long step-by-step derivation of precise error bounds.

7.2 Improving Euler's Method

Improved Euler's method is introduced as an example of a second-order method. One goal of the section is to develop the analogy between numerical methods to approximate solutions to differential equations and numerical methods to approximate integrals. The other goal is to contrast the accuracy of a second-order method with that of a first-order method (Euler's method).

Comments on selected exercises

Exercises 1-4 of this section involve the initial-value problems in Exercises 1-4 of Section 1.4. Similarly, Exercises 5-8 involve the initial-value problems in Exercises 6-9 of Section 1.4. These exercises compare the accuracy of improved Euler's method with that of Euler's method. The computations are possible by hand, but they would be tedious.

Exercises 9-13 give students practical experience with the fact that improved Euler's method is a second-order method.

Exercises 14-17 compare the accuracy of Euler's method versus improved Euler's method.

7.3 The Runge-Kutta Method

In this section the (fourth-order) Runge-Kutta algorithm is derived. As this is a method of choice for many applications, implementations of the method are given for the TI calculator, in Mathematica , and in the C programming language. We also use vector notation to discuss Runge-Kutta for first-order systems.

Comments on the exercises

In Exercises 1-5, we revisit initial-value problems that were studied in Sections 1.4 and 7.2. Using the results, we can compare the accuracy of Runge-Kutta with that of Euler's and improved Euler's methods.

Exercise 6 gives students practice with the fact that Runge-Kutta is a fourth-order method.

In Exercises 7 and 8, we illustrate how Runge-Kutta can be used to approximate solutions to first-order systems as well as second-order equations.

7.4 The Effects of Finite Arithmetic

The goal of this short section is to point out that, in practice, using a ridiculously small step size may not yield more accurate numerical approximations. The exercises study this phenomenon for particular examples.

Review Exercises

Exercises 1-5 are true/false problems. We always expect our students to justify their answers.

Exercise 6 compares the results of Euler's method, improved Euler's method, and Runge-Kutta for a first-order equation.

Exercise 7 compares the results of Euler's method, improved Euler's method, and Runge-Kutta for a competing species system.


Comments on the Labs

Lab 7.1

This lab relates the geometry of an equation with the accuracy of the numerical approximation. It addresses the question of why some methods are more accurate than expected on certain equations.

Lab 7.2

In this lab, numerical methods are used to approximate solutions of Newton's equations. Parts of the lab could be used as detailed exercises or projects much earlier in the text.