*Intermediate Physics for Medicine and Biology*, Russ Hobbie and I discuss the logistic map, a difference equation that can describe phenomena such as population dynamics. We are by no means the first to use the logistic map to illustrate deterministic chaos. Indeed, it has become the canonical example of chaos since Robert May published “Simple Mathematical Models With Very Complicated Dynamics” in 1976 (

*Nature*, Volume 261, Pages 459–467). This paper has been cited nearly 2500 times, implying that it has had a major impact.

Russ and I write the logistic equation as (Eq. 10.36 in our book)

*x*

_{j+1}= a x_{j}(1 – x_{j})where

*x*is the population in the

_{j}*j*generation. Our first task is to determine the equilibrium value for

^{th}*x*.

_{j}The equilibrium valueUsually, when Russ and I say something like “it turns out”, we include a homework problem to verify the result. Homework 34 in Chapter 10 does just this; the reader must prove that the magnitude of the slope is greater than 1 forx*can be obtained by solving Eq. 10.36 withx:_{j+1}= x_{j}= x*

x* = a x* (1 – x*) = 1 – 1/a.

Pointx*can be interpreted graphically as the intersection of Eq. 10.36 with the equationxas shown in Fig. 10.22. You can see from either the graph or from Eq. 10.37 that there is no solution for positive_{j+1}= x_{j}xifais less than 1. Fora= 1 the solution occurs atx*= 0. Fora= 3 the equilibrium solution isx* = 2/3. Figure 10.23 shows how, fora= 2.9 and an initial valuex= 0.2, the values of_{0}xapproach the equilibrium value_{j}x* = 0.655. This equilibrium point is called an attractor.

Figure 10.23 also shows the remarkable behavior that results whenais increased to 3.1. The values ofxdo not come to equilibrium. Rather, they oscillate about the former equilibrium value, taking on first a larger value then a smaller value. This is called a period-2 cycle. The behavior of the map has undergone period doubling. What is different about this value of_{j}a? Nothing looks strange about Fig. 10.22. But it turns out that if we consider the slope of the graph ofxvs_{j+1}xat_{j}x*, we find that foragreater than 3 the slope of the curve at the intersection has a magnitude greater than 1.

*a*greater than 3.

One theme of

*Intermediate Physics for Medicine and Biology*is the use of simple, elementary examples to illustrate fundamental ideas. I like to search for such examples to use in homework problems. One example that has great biological and medical relevance is discussed in Problems 37 and 38 (a model for cardiac electrical dynamics based on the idea of action potential restitution). But when reading May’s review in

*Nature*, I found another example that—while it doesn’t have much direct biological relevance—is as simple or even simpler than the logistic map. Below is a homework problem based on May’s example.

Section 10.8

Problem 33 ½Consider the difference equation

(a) Plot x_{n+1}versus x_{n}for the case of a=3/2, producing a figure analogous to Fig. 10.22.

(b) Find the range of values of a for which the solution for large n does not diverge to infinity or decay to 0. You can do this using either arguments based on plots like in part (a), or using numerical examples.

(c) Find the equilibrium value x* as a function of a, using a method similar to that in Eq. 10.37.

(d) Determine if this equilibrium value is stable or unstable, based on the magnitude of the slope of the x_{n+1}versus x_{n}curve.

(e) For a = 3/2, calculate the first 20 values ofxusing 0.250 and 0.251 as initial conditions. Be sure to carry your calculations out to at least five significant figures. Do the results appear to be chaotic? Are the results sensitive to the initial conditions?_{n}

(f) For one of the data sets generated in part (e), plotxfor 25 values of n, and create a plot analogous to Fig. 10.27. Explain how you could use this plot to distinguish chaotic data from a random list of numbers between zero and one._{n+1}versus x_{n}

## No comments:

## Post a Comment