Overview of Chapter Two

Chapter 2 marks a radical departure from the traditional ordering of topics. Rather than following first-order equations with second-order, linear equations, we go directly to autonomous first-order systems. The emphasis in this chapter is on interpreting first-order systems and their solutions and on qualitative techniques. One reason for going directly to systems is that it leads to a better qualitative understanding of autonomous second-order equations (achieving a major goal of the standard approach). It opens up the course to a wealth of topics, both theoretical and applied, associated with the modern theory of dynamical systems. A second, more technical reason is that the numerical treatment of second-order equations requires the use of first-order systems. Consequently, if numerical techniques are to have an equal role in the course, we must be able to discuss systems and vector fields at this point.

This chapter is organized in much the same way as the first half of Chapter 1. It begins with a section on modeling in which we return to the predator-prey models mentioned in Section 1.1. We also discuss the undamped harmonic oscillator. The next three sections are devoted to introducing our three-pronged approach in terms of the study of systems of differential equations. Section 2.2 introduces geometric and qualitative concepts, Section 2.3 presents some analytic techniques that apply to special systems, and Section 2.4 describes Euler's method as it is applied to systems. Unlike Chapter 1, we do not feel the need for a separate section on existence and uniqueness theory, so the Existence and Uniqueness Theorem is presented at the end of Section 2.4. Section 2.5 begins the qualitative discussion of the Lorenz system as an introduction to three-dimensional systems. This last section is the first part of a discussion of three-dimensional systems that reappears in Sections 3.8, 5.5, and 8.5. These four sections can be safely omitted if necessary. (At Boston University, we often omit these sections due to time constraints.)

Although all of Sections 2.1-2.4 contain material that is consistently used throughout the remainder of the text, it is conceivable that one could skip to Chapter 3 immediately after Section 2.2.

The use of technology to draw vector fields, direction fields, and phase portraits is essential. In order to begin to get a feel for what pictures to expect, students must see many examples accurately drawn. HPGSystemSolver often fits the bill. The CD also contains other special purpose tools that are useful at various times (see the comments below).

2.1 Modeling via Systems

This modeling section is similar in purpose to Section 1.1. Our goal is to introduce the relationship between various graphical representations of a system of differential equations, its solutions, and the interpretation of the solutions in terms of the model. The graphical representations include the phase portrait and graphs of the component functions of the solutions.

The transition to systems is natural for students. Relating the solution curve in the phase plane to the graphs of the components is considered difficult by some but is generally mastered after sufficient effort. (More practice to help develop this skill is provided in Section 2.2.) It is important for students to realize that both types of graphs are necessary because neither graph alone contains all of the information about a solution. A good analogy is trying to understand a slinky from its shadows (see Figure 2.16). Bring a slinky to class. It can be used both to demonstrate a mass-spring system and to illustrate Figure 2.16.

The section is based on predator-prey models and the undamped harmonic oscillator. We focus on these two models because we feel that these are two of the best to use at this point. In the case of the predator-prey equations, we delay an explanation of how we obtain the graphs of solutions until Section 2.4. For the undamped harmonic oscillator we simply guess solutions. At this point, we do not have enough theory to discuss the general solution. Many of the exercises focus on interpretation---both of the equations in a system and of the behavior of the solutions.

DETools

In addition to the general purpose HPGSystemSolver, the special purpose tools PredatorPrey and MassSpring provide easy-to-run classroom demonstrations that are very relevant to this section.

Comments on selected exercises

Exercises 1 and 15 involve the interpretation of the parameters in a system while Exercises 9-14 involve the interpretation of the equations.

Exercises 2-6 require an analysis of a predator-prey system similar to the one carried out in the section.

Exercises 7, 8, 16, and 17 give practice in going from the phase plane to the graphs of component functions and in the interpretation of solutions. Exercise 17 is particularly good for assigning an essay. (This predator-prey phenomenon really does occur.)

In Exercises 9-14 and 18, modifications are made to a predator-prey model. Modifying an existing model is easier than developing models from scratch, but it is still challenging. In Exercise 18, there is more than one reasonable answer.

Exercises 19 and 20 provide practice with solutions of equations for undamped harmonic oscillators. Both of these exercises prepare the students for ideas that are discussed in more detail in Sections 2.2 and 2.3.

Exercises 21-24 are modeling problems that involve Hooke's law applied to mass-spring systems.

In Exercises 25-30, models for concentrations of reactants in simple chemical reactions are developed. These models reappear in Section 5.2 (Exercises 16-20).

2.2 The Geometry of Systems

This section introduces geometric and qualitative ideas that we use throughout the remainder of the text. Vector fields, direction fields, solution curves in the phase plane, equilibrium solutions, and phase portraits are discussed. Vector fields for autonomous systems play the role of slope fields for first-order equations. They are the fundamental geometric representation of a system of differential equations.

DETools

In addition to the general purpose HPGSystemSolver, the special purpose tools DESketchPad and GraphingSolutionsQuiz help students understand the relationship between a solution curve in the phase plane and the corresponding x(t)- and y(t)-graphs.

Comments on selected exercises

Exercises 1-6 involve vector fields and phase portraits for very simple systems.

Exercises 7 and 8 involve vector fields and phase portraits for simple second-order equations.

In Exercises 9 and 10, we provide the direction field, and the student sketches the phase portrait. A description of the long-term behavior of a particular solution is also requested.

Exercises 11 is a "matching" exercise reminiscent of the matching exercise in Section 1.3. This type of exercise is an excellent way to test understanding of geometric concepts such as slope fields and direction fields. If brief explanations are required to justify why a given system corresponds to a particular direction field, students must examine the fields closely.

In Exercise 12, students find equilibrium solutions for a predator-prey system discussed in Section 2.1.

In Exercises 13-18, students are expected to solve for the equilibrium points by hand, but then they produce the phase portraits with the aid of a computer. One important goal is to give them practice using paper-and-pencil techniques and technology in the same exercise. Exercise 19 is a similar exercise for a second-order equation (an example of Duffing's equation).

In Exercise 21, students match solution curves with their corresponding x(t)- and y(t)-graphs.

Exercise 22 is one of our favorite "animated" homework problems.

In Exercises 23-26, students must produce x(t)- and y(t)-graphs from solution curves in the xy-phase plane. Their answers should include a scale on the x,y-axis, but they cannot provide a scale on the t-axis.

Exercise 27 goes the opposite direction as Exercises 23-26. The x(t)- and y(t)-graphs are given, and the student must draw the corresponding solution curve. In theory, doing this exercise is fairly mechanical, but try it before you assign it.

Exercises 28 and 29 are meant to motivate the study of uniqueness. These two exercises are repeated (exactly) as Exercises 12 and 13 in Section 2.4. Logically they make more sense in Section 2.4, but they are also stated here in case you want to get your students thinking about these issues before uniqueness is discussed.

2.3 Analytic Methods for Special Systems

Many students are uncomfortable with the heavy emphasis on geometry in Sections 2.1 and 2.2. Consequently, this section is devoted to certain special systems for which closed-form solutions can be obtained. This is also a good time to repeat the observation that any given function can always be checked to see if it is a solution, and we begin this section with that observation.

Using the method of undetermined coefficients from Section 1.8 or integrating factors from Section 1.9, we also show how one can derive the general solution of a "partially" decoupled linear system. Knowing how to solve this class of systems is useful when we study repeated eigenvalues in Chapter 3, and it provides the students with a systematic method for deriving the solutions to a system.

We also take this opportunity to introduce damped harmonic oscillators in this section. This model is used throughout the text, so the terms and equations presented on pages 193 and 194 are extremely important. Since this section is devoted to techniques that produce closed-form solutions, we also introduce the technique of guessing solutions of the form y(t) = est. The Linearity Principle does not appear until Section 3.1, so we stop short of finding general solutions of the form y(t) = k1 es1 t + k2 es2 t. However, we do make the connection between solutions of the form y(t) = est and the corresponding vector-valued solution Y(t) = (y(t), v(t)) = (est, s est).

DETools

Some of the solutions obtained in this section can be illustrated with MassSpring tool, and HPGSystemSolver can be used to plot phase portraits and x(t)- and y(t)-graphs.

Comments on selected exercises

Exercises 1-4 provide routine practice at checking functions to see if they are solutions.

Exercises 5-12 involve a given partially decoupled system. They illustrate the technique introduced in this section, and they provide additional practice with the geometric ideas introduced in Section 2.2.

Exercises 13-14 involve the analysis of a mass-spring system that involves two springs.

Exercises 15-18 provide the same type of practice for second-order equations. These equations reappear in Exercises 21-24 of Section 3.2.

Exercise 19 is a nonlinear partially decoupled system.

2.4 Euler's Method for Autonomous Systems

Section 2.4 generalizes Euler's method to systems, and we are now able to explain how the graphs and solution curves for the predator-prey equations in Sections 2.1 and 2.2 are obtained. Students usually have little difficulty understanding Euler's method for systems if they have a solid understanding of Section 1.4 and the vector field ideas introduced in Section 2.2.

If Section 1.4 was supplemented with material from Chapter 7, then it is also possible to supplement this section with material from Section 7.3.

As has already been mentioned, it is also convenient to discuss existence and uniqueness theory for systems at this point.

The example of a swaying building is presented in this section because quantitative information determines which model is more appropriate.

DETools

EulersMethodForSystems can be used in class to demonstrate the method, and the students are encouraged to use this tool when they answer Exercises 3-6.

Comments on selected exercises

Exercises 1 and 2 involve computing Euler's method solutions with fairly large step sizes and comparing the results with the direction field and/or actual solutions. The computations are tedious but manageable if done by hand, but the expectation is that the students will use some type of technology such as a spreadsheet to help with the calculations.

Exercises 3-6 were designed to be completed using EulersMethodForSystems. The phase portraits in Exercises 5 and 6 are entertaining.

Exercises 7-9 provide some practice at checking solutions, but their main intent is to illustrate the implications of the Uniqueness Theorem. Exercises 12-15 are more abstract versions of the same question.

Exercises 10-11 illustrate how Euler's method for systems applies to second-order equations.

2.5 The Lorenz Equations

We introduce the Lorenz system here mainly because it is possible to do so. Almost none of our students have seen any modern (i.e., post 1800) mathematics, and they are surprised to learn that there are unanswered questions and that there is active research in mathematics. At this point, we can only describe the Lorenz system and display some numerical solutions. Consequently, this is something of a "golly-gee-whiz" section. Three-dimensional linear systems are discussed in Section 3.8, and the Lorenz system is studied more carefully in Sections 5.5 and 8.5.

If you cover this section, we recommend your mentioning James Gleick's book Chaos. There are also a number of interesting videos that have been produced. They usually do a better job of illustrating the solution curves than we can do with our solvers. One of our favorite videos is Fractals: An Animated Discussion, by Peitgen, Jurgens, Saupe, and Zahlten, published by W. H. Freeman, New York. At this writing, there is also a Java applet written by Patrick Worfolk available at the Geometry Center.

DETools

There are two tools, LorenzEquations and ButterflyEffect, that illustrate the Lorenz attractor.

Comments on selected exercises

Exercises 1-3 and 5 cover aspects of the Lorenz system that can be verified by hand.

Exercise 4 uses ButterflyEffect to illustrate sensitive dependence on initial conditions for the Lorenz attractor.

Review Exercises

Exercises 1-12 are "short answer" exercises. The answers are (usually) one or two sentences. Most (but not all) are relatively straightforward.

Exercises 13-24 are true/false problems. We always expect our students to justify their answers. Exercises 19-24 all refer to the same phase portrait.

Exercise 25 involves solving for equilibrium points and interpreting solutions.

Exercises 26-28 involve decoupled or partially decoupled systems.

In Exercises 29-32, students sketch x(t)- and y(t)-graphs from solution curves in the phase plane.

Exercise 33 corresponds to the simple model of a glider that is given in HMSGlider.

Exercise 34 relates to Euler's method for systems.


Comments on the Labs

All of these labs require technology that is capable of plotting phase portraits. The ability to draw graphs of the coordinate functions is also very useful. HPGSystemSolver is one option that works well.

Lab 2.1 Two Magnets and a Spring

This lab is a computer exploration that involves a bifurcation. It can be started as soon as Section 2.2 is covered, and it pairs up nicely with Lab 4.1.

Lab 2.2 Cooperative and Competitive Species Population Models

This lab can be started as soon as Section 2.2 is covered. It is mostly a computer exploration, but item 2 anticipates certain qualitative ideas that are discussed more specifically in later chapters. Particular attention should be paid to the interpretation of the solutions in physical terms.

Lab 2.3 The Harmonic Oscillator with Modified Damping

Section 2.3 must be covered before this lab can be assigned. This lab is designed to give students practice applying the ideas of this chapter to various second-order equations, both linear and nonlinear ones.

Lab 2.4 A Mass-Spring System with a Rubber Band

This lab is a modification of the standard mass-spring model. It involves both a spring and rubber band. It anticipates the discussion of the Tacoma Narrows Bridge in Chapter 4 and pairs up nicely with Lab 4.3.

Lab 2.5 Active Shock Absorbers

Recently developed magnetorheological fluids allow the damping coefficient of a harmonic oscillator to be adjusted in real time, and this lab examines various models where the damping coefficient is a function of the velocity of the system.