✐ ✐
✐ ✐
Chapter 16 Differential Equations
A vast number of mathematical models in various areas of science and engineering involve differ- ential equations. This chapter provides a starting point for a journey into the branch of scientific computing that is concerned with the simulation of differential problems. We shall concentrate mostly on developing methods and concepts for solving initial value problems for ordinary differential equations: this is the simplest class (although, as you will see, it can be far from being simple), yet a very important one. It is the only class of differential problems for which, in our opinion, up-to-date numerical methods can be learned in an orderly and reasonably complete fashion within a first course text. Section 16.1 prepares the setting for the numerical treatment that follows in Sections 16.2Ð16.6 and beyond. It also contains a synopsis of what follows in this chapter. Numerical methods for solving boundary value problems for ordinary differential equations receive a quick review in Section 16.7. More fundamental difficulties arise here, so this section is marked as advanced. Most mathematical models that give rise to differential equations in practice involve partial differential equations, where there is more than one independent variable. The numerical treatment of partial differential equations is a vast and complex subject that relies directly on many of the methods introduced in various parts of this text. An orderly development belongs in a more advanced text, though, and our own description in Section 16.8 is downright anecdotal, relying in part on examples introduced earlier.
16.1 Initial value ordinary differential equations Consider the problem of finding a function y(t) that satisfies the ordinary differential equation (ODE) dy = f (t, y), a ≤ t ≤ b. dt = dy The function f (t, y) is given, and we denote the derivative of the sought solution by y dt and refer to t as the independent variable. Previous chapters have dealt with the question of how to numerically approximate, differen- tiate, or integrate an explicitly known function. Here, similarly, f is given and the sought result is
Downloaded 09/14/14 to 142.103.160.110. Redistribution subject SIAM license or copyright; see http://www.siam.org/journals/ojsa.php different from f but relates to it. The main difference though is that f depends on the unknown y to be recovered. We would like to be able to compute the function y(t), possibly for all t in the interval [a,b], given the ODE which characterizes the relationship between y and some of its derivatives.
481
✐ ✐
✐ ✐ ✐ ✐
✐ ✐
482 Chapter 16. Differential Equations
Example 16.1. The function f (t, y) =−y + t defined for t ≥ 0 and any real y gives the ODE