A First Course in Discrete Dynamical Systems - First Edition - Richard A. Holmgren
A First Course in Discrete Dynamical Systems - First Edition - Richard A. Holmgren
Editorial Board
(North America):
J.H. Ewing
F.W. Gehring
P.R. Halmos
Universitext
Editors (North America): J.H. Ewing, F.W. Gehring, and P.R. Halmos
A First Course in
Discrete Dynamical
Systems
With 55 Figures
Springer-Verlag
New York Berlin Heidelberg London Paris
Tokyo Hong Kong Barcelona Budapest
Richard A. Holmgren
Department of Mathematics
Allegheny College
Meadville, PA 16335-3902
USA
Editorial Board
(North America):
J.H. Ewing F. W. Gehring
Department of Mathematics Department of Mathematics
Indiana University University of Michigan
Bloomington, IN 47405 Ann Arbor, MI 48109
USA USA
P.R. Halmos
Department of Mathematics
Santa Clara University
Santa Clara, CA 95053
USA
987654321
DOl: 10.1007/978-1-4684-0222-3
Preface
How to use this book. This text is suitable for a one-semester course
on discrete dynamical systems. It is based on notes from undergraduate
courses that I have taught over the last few years. The material is intended
for use by undergraduate students with a year or more of college calculus
vi Preface
Richard A. Holmgren
Meadville, Pennsylvania
Contents
Preface v
How to use this book . v
Acknowledgments . vi
1 Introduction 1
1.1 Phase Portraits 4
Exercise Set 1 . . . . 6
5 Sarkovskii's Theorem 41
Exercise Set 5 . . . . . . . 45
AppendixA.
Computer Algorithms 191
A.l Iterating Functions . . . . . . . . . . . . . . . 191
Finding the Value of a Point Under Iteration 191
Tables of Iterates . . . . . . . . . . . . . . . . 192
Controlling the Precision of the Computations. 193
Graphing Iterated Functions. 193
A.2 Graphical Analysis . . 194
A.3 Bifurcation Diagrams. 196
A.4 Julia Sets . . . . . . . 197
Contents xi
References 203
Dynamical Systems. . . . . . . . . . 203
Topics in Mathematics . . . . . . . . 205
General Interest Books on Dynamics 206
Computer Programs and Algorithms 206
Index 209
which our apartment is set. For instance my address might be 2. Each year
the government has decreed that I must move to a new apartment whose
address can be found by cubing my present address and then finding the
additive inverse. That is, I apply the function f(x) = -x3 to my present
address to find my new address. So, I will be living at f(2) = _2 3 = -8
next year. The year after I will be living at f( -8) = j2(2) = _(-8)3 = 512
and after n years I will be living at r(2) = (_1)n23 n • Even ifllive to be
a very old man it is clear that I will never live in the same place twice. In
fact, each year the absolute value of my address will quickly become larger,
and I will move from one side of zero to the other.
Hence, if we start at the point 2, then fn(2) grows without bound in
absolute value and oscillates from one side of 0 to the other. What happens
if we start someplace else? Suppose we begin at the point ~. Then the
next year we will be at fG) = -~. The following year we will be at
f(-~) = f2(~) = 5~2 and after n years we will be at r(~) = (_1)n~.
So each year we will move from one side of 0 to the other, but over time
the move will not be very far.
In general, our address after n years will be r(xo) = (_1)n(xo)3 n ,
where Xo is our first address. If Ixol > 1 we see that
lim
n--+(X)
Ir(xo)1 = lim
n--+ 00
Ixol3n = 00.
r
In absolute value (xo) will increase without bound and the factor (-1) n
will cause r(xo) to oscillate from one side of 0 to the other. On the other
hand, if 0 < Ixol < 1, then lim r(xo) = O. So while r(xo) will still move
n->oo
from one side of 0 to another, its value will approach 0 as n gets larger.
In terms of choosing apartments, an apartment whose address is less than
one in absolute value would be preferable to an apartment whose address
is larger than one since, as time passes, we won't have to move so far each
year.
It remains to see what happens to -1, 0, and 1. Note that f(-1) = 1
and f(1) = -1. Hence, -1 and 1 form a periodic cycle. Thus, if we were
living at 1 and could work out a deal with the person living at -1, then
we could probably leave some of our belongings at each address since we'll
be back every other year, a luxury indeed. Of course, if we lived at 0, then
we'd not have to move at all, since 0 is fixed by f: a stable and perhaps
boring existence.
To recap, a discrete dynamical system consists of a function and its
iterates. If a dynamical system is given, we would like to know where each
point goes as we iterate the function, and what route it takes to get there.
In our example, f(x) = -x3, we have seen a number of possibilities. There
are points which head off to infinity, others which get close to 0, a pair
of points which oscillate back and forth periodically, and a point which
1. Introduction 3
EXAMPLE 1.1.
Suppose we have a population of rabbits which is growing by 10% an-
nually. Let's assume we start with Xo rabbits and denote the population
in the nth year by X n . How many rabbits will we have in n years, or
equivalently, what is the value of xn?
Clearly,
Xl = Xo + .1xo = 1.1xo
X2 = Xl + .1XI = 1.1xI
Xl = p(xo)
X2 = p(xd = (p 0 p)(xo) = p2(xO)
X3 = p(X2) = (p 0 p2)(xO) = p3(xO)
EXAMPLE 1.2.
Let 1(x) = x 2 • The dynamical system we are considering consists of the
domain of 1 and the function itself. The domain is the set of real numbers
which we represent by a line. Note that 0 and 1 are fixed, that is, 1(0) = 0
and 1(1) = 1. We indicate this in our phase portrait by dots at 0 and 1.
We notice that if 0 < x < 1, then In(x) tends towards 0 as n becomes
larger and if x> 1, then In(x) tends towards infinity as n becomes larger.
To represent this, we draw an arrow from 1 towards 0 and another from
one towards positive infinity. Now what happens if x is less than O? The
point -1 becomes fixed at 1 since 1(-1) = 1. This implies r( -1) = 1 for
all n larger than or equal to 1. This is shown by the arrow from -1 to 1.
Points which lie between -1 and 0 are mapped into the interval (0,1) and
1.1. Phase Portraits 5
then they move towards 0 under iteration of the function. Similarly, points
which are to the left of -1 are mapped to the right of 1 and then move
towards infinity under iteration of the map. We represent this by arrows
from the negative portion of the line to the positive portion. All of this
information is encoded in the phase portrait shown in Figure 1.1. 0
EXAMPLE 1.3.
Consider the function f(x) = -x 3 . 0 is fixed and the point 1 goes to
-1 and then returns to 1 on the second iteration. Points which are greater
than 1 in absolute value oscillate from side to side of zero under iteration
of f and grow ever larger in absolute value. Points which are less than 1 in
absolute value oscillate from side to side of zero under iteration of f and
grow ever smaller in absolute value. All of this is encoded in the phase
portrait shown in Figure 1.2. 0
The reader will be given an opportunity to explore and develop more
phase portraits in the exercises.
6 1. Introduction
Exercise Set 1
EXAMPLE 2.2.
Let f : lR -lR be defined by f(x) = x 2. Then
f([-1,2]) = [0,4] and f-1((1,4]) = [-2, -1) U (1,2].
Note that
f(1-1([1,4])) = [1,4], but f-1(1([1,2])) = [-2,-1] U [1,2].
Note also that the use of the notation f- 1 does not necessarily imply that
f is an invertible function. 0
DEFINITION 2.3. A function is one-to-one if there is exactly one element
of its domain assigned to each element of its range.
The "if and only if" in this proposition means that the two statements,
"f : D - C is one-to-one" and "f(x) = f(y) implies x = y" are equivalent.
In other words, if one of the statements is true, then the other is true; if one
of the statements is false, then the other is false. The proof of Proposition
2.4 is not hard and is left as an exercise.
Proposition 2.4 provides us an easy way to test whether or not a function
is one-to-one. For example, consider the function f : lR - lR defined by
f(x) = x 2. Clearly, f is not one-to-one since f(l) = 1 = f(-I). On the
°
other hand, the function g: [0,00) -lR defined by g(x) = x 2 is one-to-one
since if g(x) = g(y), then x 2 = y2, X ~ 0, and y ~ imply x = y. Note that
even though the functions f and 9 have the same rule, they have different
domains and, hence, different properties.
A graphical test for determining whether functions of the real numbers
are one-to-one is the horizontal line test. A function is one-to-one if and
only if every horizontal line crosses the graph of the function at most once.
By looking at the graphs of f(x) = x 2 and p(x) = X 1/ 3 we see that f is not
2. A Quick Look at Functions 11
-2
-I -2
EXAMPLE 2.6.
Construct an example of a function f : [1,3] ---+ [2,4] which is onto but
not one-to-one.
Let us first determine what the graph of f must look like. Since all the
values of f must lie in the codomain, the graph lies between the lines y = 2
and y = 4. As the domain of f is [1,3]' a vertical line will intersect the
graph if and only if it passes through a point on the x-axis between and
including 1 and 3. So the graph of the function lies in the box with corners
(1 , 2), (1 , 4), (3,4), and (3 , 2) as shown in Figure 2.2.
Since f is onto, we know that every horizontal line which passes through
the codomain (thought of as a subset of the y-axis) must intersect the graph
of the function. Also, since f is not one-to-one there must be at least one
12 2. A Quick Look at Functions
(1.4) (3,4)
- -+t--
(3,2)
2 4
horizontal line through the codomain which strikes the graph two or more
times. Obviously, there are many graphs which satisfy these criteria. One
example is shown in Figure 2.2. 0
While it is often sufficient to think of continuous functions as functions
whose graphs are unbroken lines, it is essential to have a more precise
definition. In particular, if a function's domain or codomain is not an
intervai of the real line, then the graph of the function is not always so easy
to draw and not always informative when it can be drawn. For example, if
Z is the set of integers, then the function J : Z --> Z defined by J(z) = 2z
is continuous, but its graph is a set of discrete points.
DEFINITION 2.7. Let J : D --> C and Xo be a point in D. Then J
is continuous at Xo if for every positive number E, there exists a positive
number 8 such that
if x is in D and Ixo - xl < 8, then IJ(xo) - J(x)1 < E.
x I(x) E
• • • j(xo )
D c
EXAMPLE 2.8.
Let Z be the set of integers. We use the definition to show that the
function f : Z -+ Z defined by f (z) = 2z is continuous.
Suppose a is any integer and f is any number greater than o. We must
find a number 8 which is greater than 0 and such that If(z) - f(a)1 < f
whenever z is an integer satisfying Iz - al < 8. Let's try 8 = If z is an !.
integer and Iz - al < !, then z must be equal a. Hence, If(z) - f(a)1 is
equal to 0, which is less than f. 0
Lest we think it is always easy to find a 8 which works , we consider
another example. However, in oder to do that we need to introduce the
Triangle Inequality.
The Triangle Inequality needn't be mysterious. As most of us realize,
it is further to go from New York to Chicago by way of Los Angeles than
it is by traveling directly; this is not surprising. Essentially the Triangle
Inequality tells us that the distance from the point a to the point b is less
than or equal to the distance from a to the point e plus the distance from
e to b. We can express the Triangle Inequality mathematically by
la - bl :::; la - el + Ie - bl (2.1)
where a, band e represent real numbers. Recall that the value la - bl
represents the distance from the point a to the point b. If we let x = a - e
and y = e - b in formula 2.1, then we get
This is th€ form of the Triangle Inequality which is most useful to us; we
use it in virtually every chapter of this text. The reader is well advised to
take the time to prove it and to learn how to use it. Exercise 2.8 at the end
of this chapter provides hints for proving the inequality and practice using
it. Our first use of the the Triangle Inequality is in the following example.
EXAMPLE 2.9.
Show that the function g(x) = x 2 is continuous on the real numbers.
Let a be any real number. We must show g that is continuous at a.
We consider the case where a is not O. Let 10 be any positive real number.
We must find a positive number 8 such that whenever x is a real number
satisfying Ix - al < 8, then If(x) - f(a)1 = Ix2 - a2 1< E. To accomplish
this, we choose 8 so that 8 < lal, 8 < 13:1' and 8 > O. This is possible since
both 10 and lal are positive numbers. If Ix - al < 8, then using the Triangle
Inequality we find that
EXAMPLE 2.11.
I, for x < 1
Let f : [0,2] --+ lR be defined by f(x) = { . Then f(O) = 1
3, for x 2: 1
and f(2) = 3. Even though 2 is between f(O) and f(2) there is no c in
[0,3] such that f(c) = 2. D
DEFINITION 2.12. Let f : D --+ C be a function with range f(D). The
function 9 : f(D) --+ D is an inverse of f if (g 0 f)(x) = x for all x in D.
If f has an inverse, then it is usually denoted f- 1 and we say that f is
invertible.
The following proposition is well-known and its proof is left as an exer-
cise:
PROPOSITION 2.13. A function is invertible if and only if it is one-to-
one.
DEFINITION 2.14. If a function is one-to-one, onto, and continuous, and
its inverse is continuous, then it is a homeomorphism. In this case, we say
the domain and codomain are homeomorphic to one another.
Homeomorphisms play an important role in mathematics in general,
and in dynamical systems, in particular. When two mathematical spaces
or structures are homeomorphic they are, in some sense, the same. Math-
ematically, we describe this by saying that they have the same topolog-
ical properties. The intervals [1,2] and [3,5] are homeomorphic since
9 : [1,2] --+ [3,5] defined by g(x) = 2x + 1 is a homeomorphism. (The
reader should check that 9 is in fact one-to-one, onto, continuous and has a
continuous inverse.) However, the intervals [1,2] and (3,5) are not homeo-
morphic. That is, there is no function f : [1,2] --+ (3,5) which is one-to-one,
16 2. A Quick Look at Functions
Exercise Set 2
2.7 Let f be a function whose domain and codomain are subsets of the
real numbers. Prove f is one-to-one if and only if every horizontal
line intersects the graph of f at most once.
a) Let a and b be real numbers. Prove: la+ bl ~ lal + Ibl. You may
use the fact that -Ial ~ a ~ lal and -Ibl ~ b ~ Ibl or you may
wish to show that la + W
~ (Ial + Ib1)2.
-1 for x < 0, . .
2.12 a) Show that f(x) = { 1
for x > ° IS contmuous.
-1 for x S 0, . .
b) Show that f(x) = { 1 IS not contmuous.
forx>O
2.13 Show that every function whose domain is the integers and whose
codomain is the real numbers is continuous.
2.15 Prove that if a and b are real numbers satisfying a < b, then the
interval [a, b] is homeomorphic to [0,1].
Hint: Find a linear function f : [a, b]-t [0,1].
EXAMPLE 3.5.
a) The interval (a , b) is open.
Let x be a point in (a, b). To show that (a, b) is open, we must find E > 0
so that Iy - xl < E implies y is in (a, b) . Suppose E is the smaller of Ix - al
and Ix - bl. In other words, suppose E is the distance from x to the nearest
endpoint of (a, b). Since a < x < b, E is positive. If Ix - yl < E, then
Ix - yl < E :::; Ix - al·
Since x > a, this implies
-x + a < x - y < x-a.
Solving this for y we find that
a < y < 2x - a.
3. The Topology of the Real Numbers 23
Hence, a < y. To complete the proof that y is in (a, b), we need to show
that y < b. The proof is similar to the one above. The details are left as
an exercise.
Since there are open sets in mathematics, it's not surprising there are
closed sets as well. To define closed sets we use the concepts of convergent
sequences and accumulation points.
DEFINITION 3.6. Let Xl, X2, X3,'" be a sequence of real numbers. The
sequence converges to X if for each E > 0, there exists N such that if k ;::: N,
then Ix-xkl < E. The sequence grows without bound or converges to infinity
if for each M, there exists N such that if k ;::: N, then Xk > M.
Again, we can think of E as a a distance. If the sequence Xl, X2, X3, ...
converges to X, then when k is large, Xk is close to x, that is, IXk - xl < E.
Notice that in the case of convergent sequences we must be able to find an
N which works for each E. The N may change as E changes and it is not
sufficient to find an N which works for one particular E.
EXAMPLE 3.7.
a) The sequence Xl = .9, X2 = .99, X3 = .999, ... converges to 1.
To see this, let E > O. Then there exists an integer N such that l~N < E.
Thus, 11 - Xk I = 16k < E for all k ;::: N as desired.
b) The sequence Xl = 1, X2 = 2, X3 = 1, X4 = 2, X5 = 1, ... does not
converge to any real number.
To demonstrate this, it is sufficient to find a single E for which there is
no N. Choose E = [Link] claim that if Xk is such that IXk - xl < then i,
IXk+! - xl > i·To see this suppose that Xk = 1; it has to be either 1 or 2.
Then Xk+l = 2 and by the Triangle Inequality,
Ixx+! - xl = 12 - xl
= 11 - (x - 1)1
3
;::: 111 -Ix - 11 > 4'
A similar result holds when .Tk = 2. Hence, when E = i, there does not
exist x and N such that IXk - xl < E for all k ;::: N. 0
24 3. The Topology of the Real Numbers
EXAMPLE 3.10.
a) The point 0 is an accumulation point of (0,1]. If E is any positive number,
then either 1 < E, or ~ is an element of (0, 1] and 0 < I~I < E. Hence, there
is a point in (0,1] whose distance from 0 is positive and less than Eo So 0
is an accumulation point of (0,1) by part (b) of Proposition 3.9.
b) The irrational number x = .123456789101112131415 ... is an accumu-
lation point of the set of rational numbers.
To prove this we define a sequence of rational numbers by
Xl =.1
X2 = .12
X3 = .123
X4 = .1234
EXAMPLE 3.12.
a) The interval 1= [2,4] is closed.
We prove this by showing that if x is not in I, then x is not an accumu-
lation point of I. This implies that if x is an accumulation point of I, then
x is in I.
Assume x is not in I. Now either x < 2 or x > 4. If x < 2, then
2 - x is a positive number and the neighborhood of x with radius 2 - x is
(x - (2 - x), x + (2 - x)) or (2x - 2,2). Since this neighborhood doesn't
contain any point of [2,4], x is not an accumulation point of [2,4]. A similar
argument will show that if y > 4, then y is not an accumulation point of
[2,4]. Consequently, [2,4] must contain its accumulation points.
b) The set of rational numbers is not closed in the real numbers since
we demonstrated a sequence of rational numbers which converged to an
irrational number in Example 3.10b. 0
We have seen that the rational numbers are neither open nor closed in
the real numbers so it is not true that all sets must be either open or closed.
In the exercises we ask the reader to show that the interval [1,2) is neither
open nor closed.
If A is a set of real numbers, then the complement of A is the set of
all numbers which are not contained in A. For example, the complement
of [1,2) is (-00,1) U [2,00). Often we can discover properties of a set by
examining its complement. The following proposition employs this tech-
nique.
PROOF. Let A be a subset of the real numbers and let B be the com-
plement of A. We begin by showing that if A is open, then B is closed.
Suppose A is open. If x is an accumulation point of B, then every
neighborhood of x contains an element of B. Consequently, there is no
neighborhood of x which is completely contained in A. As A is open,
Proposition 3.3 implies x is not in A. It follows that x is in B and so
B contains all of its accumulation points. Hence, if A is open, then B is
closed.
It remains to show that B is open when A is closed.
Suppose A is closed. Let x be a real number such that for all E > 0 there
exists y in A satisfying Ix - yl < E. Then x must be in A and it follows
that B is open. Why? 0
26 3. The Topology of the Real Numbers
EXAMPLE 3.17.
The rationals are dense in ~.
It suffices to show that every irrational number is an accumulation point
of the rationals.
Let t = * * .tlt2t3t4 ... be an irrational number expressed in decimal
notation. Define a sequence of rational numbers by
Xl = * * .tl
X2 = * * .tlt2
X3 = * * .tlt2t3
X4 = * * .tlt2t3t4
We can easily prove that this sequence converges to t. Hence, the rationals
are dense in the real numbers. 0
It is interesting to note that the irrational numbers are also dense in the
real numbers. This demonstrates that even though a set is large enough to
be dense, it doesn't necessarily have to be so large that there isn't enough
room left over for other large (that is, dense) sets. In fact, we can find
an infinite collection of subsets of the real numbers, each of which is dense
in the real numbers, such that no two of the sets have points in common.
Next time you are waiting patiently in a doctor's office or for a lecture to
end you might try to do just that.
Exercise Set 3
3.2 Prove that the interval (a, b) is open. (Note that most of the proof
is contained in Example 3.5.)
3.4 Show that the interval [1,2) is neither open nor closed.
3.6 Indicate whether the following sets are open, closed, neither open
nor closed, or both open and closed. Justify your answer.
a) [-1,0) U (O,lj
b) (-l,OjU[O,l)
c) (-l,O)n[O,l)
d) The set of natural numbers, {I, 2, 3, ... }.
e) U(-l,~)
n=l n n
3.11 Prove that lR is both open and closed. Use Proposition 3.13 to
show the empty set 0 is both open and closed. Are there any other
subsets of lR which are both open and closed?
°
g(a) = f(a) - a > and g(b) = f(b) - b < 0. Since g is continuous, the
Intermediate Value Theorem implies there is c in [a, b] such that g(c) = 0.
°
But g(c) = f(c) - c = so f(c) = c and we are done. Figure 4.1 illustrates
this proof. 0
EXAMPLE 4.3.
The function f(x) = 1 - x 2 has a fixed point in the interval [0,1].
To see this we might note that f is a continuous function and the range
of f over [0,1] is contained in [0,1]. Thus, by theorem 4.2, f has a fixed
point in [0,1].
Alternatively, we look at the graph of f shown in Figure 4.2 and note
that it crosses the line y = x in the interval [0,1]. Hence, f has a fixed
point in [0,1].
Finally, we can find the fixed point algebraically by solving the equation
f(x) = x or 1 - x 2 = X . By doing this we find that the fixed points of f
are -~ ± 1 and note that -~ + 1
is in [0,1]. 0
in the interval. In the following theorem we show that there is a fixed point
in the interval when the containment is reversed.
PROOF. Let I = [a, b]. Since f(I) :J I there are c and d in I such that
fCc) = a and fed) = b. If c = a or d = b we are done. If not, then a < c < b
and a < d < b. If we define g(x) = f(x)-x, then g(c) = f(c) -c = a-c < 0
and g(d) = fed) - d = b - d > O. Since g(c) < 0 and g(d) > 0, and g is
continuous, the Intermediate Value Theorem implies there is e between c
and d (and hence in I) satisfying gee) = 0 and fee) = e. The reader is
encouraged to draw a figure which illustrates this proof. 0
EXAMPLE 4.7.
Let g(x) = Ix -11- Then 0 and 1 form a periodic cycle. f(2) = 1 so 2 is
eventually periodic. In fact, every integer is eventually periodic. 0
The function f(x) = -x3 has a fixed point at 0 and a periodic cycle
consisting of 1 and -1. Examination reveals that there are no eventually
periodic points which are not periodic. How should we characterize the rest
of the points? For example, if we start at ! and iterate with f, then we get
the sequence, "2'
1 - 23'
1 2§"'
1 - 227'
1 . . .. S0"21·IS not perlO
. d'lC and never reach es
0, though it does get closer and closer to it. That is, fn (!) converges to 0
as n goes to 00. We say! is forward asymptotic to O.
DEFINITION 4.8. Let f be a function and p be a periodic point of f with
prime period k. Then x is forward asymptotic to p if the sequence x, fk (x),
Pk(x), Pk(X), ... converges to p. In other words, limn_co rk(x) = p.
The stable set of p, denoted by WS(p), consists of all points which are
forward asymptotic to p.
4. Periodic Points and Stable Sets 35
If the sequence lxi, If(x)l, IP(x)l, IP(x)l, ... grows without bound, then
x is forward asymptotic to 00. The stable set of 00, denoted by WS(oo),
consists of all points which are forward asymptotic to 00.
Notice that when we are searching for points in the stable set of a point
with prime period k we must consider the sequence
x, fk(x), f2k(x), f3k(x), ...
not the sequence
EXAMPLE 4.9.
a) Let f(x) = -x3. Then the stable set of 0 consists of all points in the
interval (-1, 1), the stable set of 00 consists of all points in the intervals
(-00, -1) and (1,00), the stable set of -1 contains only the point -1
and the stable set of 1 contains only the point 1. Symbolically we write
WS(O) = (-1,1), WS(oo) = (-00,-1) U (1,00), WS(-l) = {-I}, and
W S(l) = {I}. Notice that it doesn't make sense to discuss the stable set
of any other points since f has no other periodic points.
b) In Example 4.7 we examined the function g(x) = Ix -11- Returning to
that example we see WS(O) is the set of all even integers and W S (l)is the
set of all odd integers. 0
In the preceding example we tacitly assumed that if a point is in the
stable set of one periodic point, then it cannot be in the stable set of a
different periodic point. This idea is important and we prove it in the
following theorem.
THEOREM 4.10. The stable sets of distinct periodic points do not in-
tersect. In other words, if p and q are periodic points and p =I- q, then
WS(p) n WS(q) = 0.
PROOF. Let PI and P2 be periodic points of period kI and k 2 , respec-
tively. We'll show that if WS(PI) n WS(P2) =I- 0, then PI = P2.
Let x be in WS(PI) n W S (P2). Then for each E > 0 there exist NI and
N2 such that n ::::: NI implies IPI - rk1(x)1 < ~ and n ::::: N2 implies
Ip2 - rk2(x)1 < ~. If M is the larger of NI and N 2 , then n ::::: M implies
both IPI - r k1 (x) I < ~ and Ip2 - r k2 (x) I < ~. By the Triangle Inequality
we have
IPI - P21 = IPI - rklk2(x) + rklk2(x) - P21
~ IPI - rklk2(x)1 + Irklk2(x) - P21
E E
< 2" + 2" = E.
36 4. Periodic Points and Stable Sets
Since we have shown the distance between PI and P2 is less than E for each
E > 0, it must be that PI = P2. D
°
the fixed point 0, the periodic points 1 and -1, and points in the stable
sets of or 00. As we continue our investigations we shall find that not all
systems are quite so easy to characterize. We will need new tools and ideas
to classify the dynamics of even a simple function like h(x) = 4x(1 - x).
One of the simplest tools we have at our disposal is graphical analysis.
EXAMPLE 4.11.
We begin our examination of the dynamics of the function f(x) = x 3
by graphing f and the line y = x on the same set of coordinate axes as in
Figure 4.3, and try to determine the itinerary of the point a which lies in
the interval (0,1). Rather than follow the movement of the point on the
x-axis we keep track of its whereabouts on the line y = x. It starts at a
(or, if you prefer, at the point (a, a)). From that point we travel vertically,
until we strike the graph of f. Since we moved vertically, the x value of
the point of intersection must be a and the y value must be f(a). We
now travel horizontally back to the line y = x to arrive at the point f(a).
Repeating the process we travel vertically to the graph of f to arrive at
P
°
the point (f (a), (a)) and then horizontally back to f2 (a) on the y = x
line. Continuing, we see that r(a) approaches as n approaches 00. After
°
trying more points in the interval (0,1), we conclude that every point in
(0,1) will approach under iteration of f(x) = x 3 . Further examination
suggests that W S (0) = (-1,1) and W S (00) = (-00, -1) U (1,00). The
points -1, 0, and 1 are fixed points of f. D
EXAMPLE 4.12.
We now consider the function g( x) = _xl/3. The graph of 9 and the line
y = x are shown in Figure 4.4. We start with a point a which is near and °
°
in the interval (0,1). Note that as n grows the value of gn(a) approaches 1
in absolute value and oscillates from one side of to the other. Note also
4.1. Graphical Analysis 37
(/2(a),f3(a»
~=====:;t:====~:JL!(/t4(a),f5ca»
that -1 and 1 form a periodic cycle with prime period 2. Further, gn (a)
is always positive when n is even. Thus, we conclude that a is in W S (1).
Further analysis suggests that WS(1) = (0,00), WS(-1) = (-00,0), is
fixed, and WS(oo) = 0. D
°
The reader is asked to develop additional examples of graphical analysis
in the following exercises.
Exercise Set 4
4.1 Let g( x) = Ix -11. Find all fixed points and eventually fixed points
of g. How many periodic points does g have? What are they? How
many eventually periodic points does g have? What are they?
4.2 From exercise 4.1 we know that ~ is fixed by g(x) = Ix - 11. Find
WS(~).
4.3 Find two points in [0, 1] that are eventually fixed at °after exactly
3 iterations of h(x) = 4x(1 - x).
4.4 a) Can a function have a point which is eventually fixed after ex-
actly 2 iterations, if it has no points which are eventually fixed
after exactly 1 iteration? If so, give an example; if not, explain
why.
b) Suppose n is a natural number. Can a function have a point
which is eventually fixed after exactly n + 1 iterations, if it has no
points which are eventually fixed after exactly n iterations? If so,
give an example; if not, explain why.
4.6 For each of the following functions, find all periodic points and
describe the stable sets of each. You may wish to use graphical
analysis.
a) q(x) = x2
b) E (x) = eX - 1
c) p( x) = x 2 - .r
d) r(x) = -~ arctan x
Exercise Set 4 39
e)q(x)=x2+i
* f) h(x) = 3.5x(1 - x)
4.7 Show that all periodic points of h(x) = 3.2x(1-x) lie in the interval
[0,1].
4.8 Find a function 1 : [0,1] -+ [0,1] with no fixed points and draw its
graph.
4.9 Find a function 1 : [0,1] -+ [0,1] with exactly three fixed points
and draw its graph.
4.11 Draw a figure that illustrates the proof of Theorem 4.4. Use your
figure to explain the proof.
4.13 Find a continuous function of the real numbers which has a periodic
point of period three and indicate the periodic point.
* What other periodic points does the function have?
Hint: Choose the periodic orbit first. Use the fact that (a, b)
is on the graph of 1 if and only if 1(a) = b and find a continuous
function whose graph goes through the appropriate points.
5
Sarkovskii's Theorem
(4) r-l(An _ l ) = 10
(5) r(An) = h
Since An C h, (5) and Theorem 4.4 imply that r has a fixed point in
An. Of course, this is identical to saying that f has a periodic point p
with period n in An. In the next two paragraphs we use the other four
conditions to show that p has prime period n.
If p is a periodic point of period n in An, then (1) and (3) imply that
p, !(p), P(p), ... , r - 2(p) are in II = [b,cl and (4) implies r-l(p) is in
10 = [a, bl. If p = c, then f(P) = a which is in not in h. Since the only
one of the first n iterates of p which doesn't have to be in h is fn-l(P),
it must be that n = 2. However, this contradicts the fact that the prime
period of c is three. Hence, p must be in the half-open interval [b, c).
If p = b, then n = 3 since p(P) = a which is not in II and the only
one of the first n iterates of p which is not in h is r-l(p). Assume that
n is not three. It follows that p must be in the open interval (b, c). Since
r-l(p) is in 10 = [a, bl, r-l(p) is not equal to p and so p can't have prime
period n - 1. If the prime period of p were less than n - 1, then (3) and
the fact that p is not b or c imply that the orbit of p is contained entirely
in (b, c) and this would contradict (4). Since r- l (p) is not an element of
fk(Ak) = Ao = II
for each k = 1,2, ... ,n - 2 as required by (3). To define A n - 1 consistent
with (4) we note that
r - 1(A n_ 2 ) = f(r- 2 (A n _ 2 )) = f(11) ::) 10 .
Hence, there is A n - 1 C A n- 2 such that f n - 1(A n_ 1) = 10 , Finally,
r(An-d = f(r- 1(A n _ 1)) = f(Io) ::) II
implies there is An C A n- 1 such that r(An) = II as required by (5) and
the proof is complete. 0
As surprising as Theorem 5.1 might be, it is only a special case of a more
general theorem proven by A. N. Sarkovskii in 1964. This interesting and
beautiful result depends only on the continuity of the function in question
and should give pause to anyone who might think that all of the really
good theorems which don't rely on the mathematical equivalent of a B-1
bomber were proven before the turn of the century.
In his theorem, Sarkovskii places an order on the integers which we
define below. The statement of the theorem follows the definition.
DEFINITION 5.2. SARKOYSKII'S ORDERING. Sarkovskii's ordering of
the natural numbers is
3 >- 5 >- 7 >- ... >- 2 . 3 >- 2 . 5 >- 2 . 7 >- ...
. . . >- 22 . 3 >- 22 . 5 >->- 22 . 7 >- . . . >- 2n . 3 >- 2n . 5 >- 2n . 7 >- ...
. . . >- 23 >- 22 >- 2 >- 1.
a >- b indicates a precedes b in the order. When writing the order, all odd
numbers except one are listed in ascending order, then two times every odd,
then four times every odd, and so on. The order is completed by listing the
powers of two in descending order. Every natural number can be found
exactly once in Sarkovskii's ordering.
THEOREM 5.3. [SARKOYSKII, 1964] Suppose that f : lR --+ lR is contin-
uous and that f has a periodic point with prime period n. If n >- m in
Sarkovskii's ordering, then f also has a periodic point with prime period
m.
While the proof of Sarkovskii's Theorem uses no tools which are not
used in the proof of Theorem 5.1, it is somewhat longer. However, in
spirit it is much the same. Essentially we show that if n >- m, then the
existence of a prime period n point implies the existence of a sequence of
closed intervals h, 12 , •.• ,1m , such that f(1i) ::) IH1 for all i < m and
f(1m) ::) h· Then Theorem 4.4 implies there is a period m point in II.
44 5. Sarkovskii's Theorem
If one of the intervals is disjoint from all of the other intervals, then there
is a point with prime period m. The key then is to find these intervals h
Rather than write the proof out in detail here, we refer the reader to the
proof by Huang on pages 91-102 of the April 1992 issue of Mathematics
Magazine. Huang's exposition is clear and complete. A similar proof can be
found in Devaney's text, Introduction to Chaotic Dynamics, which is listed
in the references. Both authors demonstrate that Sarkovskii's Theorem
is sharp by constructing continuous functions which have period 5 points,
but no period 3 points. The technique used allows for the construction of
a continuous function which has a periodic point with prime period n, and
no periodic point with prime period which precedes n in the Sarkovskii
ordering.
In addition to implying the existence of orbits with prime periods lower
than that of a known orbit, Sarkovskii's Theorem can be used to show that
points with certain periods do not exist. For example, in Figure 5.1 the
graphs of h(x) = 3.2x(1 - x), h2 (x) and h4 (x) are shown.
In the first graph we see that h has fixed points at 0 and near .7. In the
second graph, we see the fixed points again, and a period 2 orbit near .5
and .8. In the third graph, we see that the only period 4 points are at 0
and near .7, .5, and .8. That is, h has no prime period 4 orbits. Therefore,
we can conclude that h has no orbits with prime period other than 2 and
1, since those are the only two orbits less than 4 in the Sarkovskii ordering.
We should note that it is safe to assume that all of the periodic points of
h lie in the interval [0, 1], so we do not need to graph any other portion of
the function. The reader should verify that this assumption is valid.
Exercise Set 5 45
Exercise Set 5
5.1 Let hand 10 be as described in the proof of Theorem 5.1. Use the
Intermediate Value Theorem to show that lUo) ~ h, IUd ~ h,
and IUd ~ 10 .
5.3 Show that l(x) = """"-rr4 arctan(x) does not have a point with prime
period 32.
Hint: Use a computer graphics package to look at the graphs
of the iterates of f.
5.5 Let 1 : lR --t lR be continuous, n > 3, and Xl, X2, ... ,Xn be points
such that Xl < X2 < ... < Xn . Show that if l(xi) = Xi+! for
i = 1,2, ... , n - 1 and l(xn) = Xl, then 1 has points with all
prime periods.
x-a X - a
exists. In this case we say f is differentiable at a and call the limit f' (a)
or the derivative of f at a. A function is differentiable if it is differentiable
at each point in its domain.
Recall that a function which is differentiable at a point is also continuous
at that point. An important property of differentiable functions which we
will often exploit is the Mean Value Theorem.
THEOREM 6.2. MEAN VALUE THEOREM. Let f : [a, b] -* IR be differ-
entiable on [a, b]. Then there exists c in [a, b] such that
f(b) - f(a) = f'(c)(b - a).
After examining a few graphs the reader will realize that the Mean Value
Theorem is an intuitively clear result. A proof can be found in any good
calculus text. We use the Mean Value Theorem to get another result about
48 6. Differentiability and Its Implications
fixed points. Theorem 4.2 states that every continuous function which maps
a closed interval back into itself has at least one fixed point. If the function
is differentiable on the closed interval and the value of the derivative is less
than one in absolute value, we get a much stronger result.
THEOREM 6.3. Let I be a closed interval and f : I --+ I be a differen-
tiable function satisfying 1f'(x)1 < 1 for all x in I. Then f has a unique
fixed point in I. Moreover, if x and yare any two points in I and x I- y,
then If(x) - f(y)1 < Ix - yl·
In the hypothesis of Theorem 6.3 we assume f is differentiable on I.
Consequently, f must also be continuous on I and Theorem 4.2 guarantees
that f has a fixed point. The additional differentiability condition guaran-
tees that f has only one fixed point and points in I become closer together
as we iterate the function. In fact, it can be shown that I is contained
in the stable set of the unique fixed point. If we assume f' is continuous,
then this is a corollary of the proof of Theorem 6.5. If f' is not contin-
uous, then the proof requires more sophisticated tools from mathematical
analysis than we have at our disposal. In particular, an understanding of
least upper bounds is useful. Finally, we note that, as with Theorem 4.2,
Theorem 6.3 does not necessarily hold when I is an open interval. (See
exercise 4.10.)
The proof of Theorem 6.3 follows from the Mean Value Theorem.
PROOF. We begin by proving the second part of the theorem. Let x
and y be any two points in I satisfying x I- y. Without loss of generality
we assume x < y. Then our hypotheses assert that f is differentiable on
[x, yj. By the Mean Value Theorem, there exists c in [x, yj such that
EXAMPLE 6.4.
a) Let f(x) = ~x +~. Solving ~x + ~ = x we find that 3 is the only fixed
point of f. Notice also that fl (x) = ~ for all x. Hence, if x is any real
number other than 3, then the Mean Value Theorem implies there exist
real numbers, Cl, C2, and C3 such that
and
Second, in the case where If'(p) I > 1, it is not generally true that once
a point leaves the open interval containing p it is gone for good. When
we examine the logistic function more carefully we will see cases where a
point's iterates return over and over again. As we demonstrated in Chapter
3, we may assume the open set containing p is a neighborhood of p.
Finally, as we noted in Example 6.4, it is true that when the derivative
is negative on the open interval U, then x and f(x) are on opposite sides
of p and when the derivative is positive on U, then x and f(x) are on the
same side of p. The reader is asked to prove this in the exercises.
To sum up, we restate Theorem 6.5.
THEOREM 6.5 (RESTATED). Let f be a 0 1 function and p be a fixed
point of f. Then If'(p) I < 1 implies that there is a neighborhood ofp which
is contained in W'(p) and If'(p)1 > 1 implies that there is a neighborhood
of p all of whose points must leave the neighborhood under iteration of f·
PROOF. Let f be differentiable with a continuous derivative, assume pis
a fixed point of f, and suppose 1f'(p)1 < 1. We wish to find a neighborhood
of p that is contained in WS (p). Our strategy is to find 8 > 0 and E > 0
such that 1f'(x)1 < 1- E for all x in (p- 8,p+ 8). We then apply the Mean
Value Theorem to show that (p - 8,p + 8) C WS(p).
Let E = ~(1 - If'(p)l). Since 1f'(p)1 < 1, E > O. To find 8 we use
the continuity of f'. Since f' is continuous there exists 8 > 0 such that
If'(x) - f'(p) I < E whenever Ip - xl < 8. (The reader should be able to
explain why this is so.) Equivalently, If'(x) - f'(p) I < E whenever x is in
(p - 8,p + 8). As E = ~(l-If'(p)I), this implies that if x is in (p - 8,p+ 8),
then
1!,(x)1 = + !,(p) I
If'(x) - f'(p)
~ If'(x) - f'(p) I + If'(p) I (by the Triangle Inequality)
< E + If'(p) I
= ~(l-If'(p)1) + 1f'(p)1
= ~ + ~1!,(p)1
= 1-(~(1-I!,(p)l))
=1-E.
EXAMPLE 6.6.
a) Let g( x) = x - x 3 . Then 0 is a fixed point of 9 and g' (0) = 1. Graphical
analysis easily demonstrates that (-1,1) C WS(O).
b) Let p(x) = x + x 3 . Then 0 is again a fixed point of p and p'(O) = 1.
However, in this case WS(O) = 0 and WS(oo) contains all real numbers
except O.
c) Let k(x) = eX -1. In this case 0 is a fixed point of k and k'(O) = 1, but
WS(O) = (-00,0] and WS(oo) = (0,00). That is, all points to the left of 0
approach 0 under iteration of k, but all other points grow without bound
under iteration of k.
The graphs of each of the functions in this example are shown in Fig-
ure 6.2. 0
6. Differentiability and Its Implications 53
y=x y=x
y=x
k(x)
p(x)
Example 6.6 demonstrates that nothing definitive can be said about the
behavior of points near a fixed point whose derivative is equal to one in
absolute value.
Recall that p is a periodic point with period n of the function f if p is a
fixed point ofr. In light of this fact, the following definition of a hyperbolic
periodic point is an obvious extension of the definition of hyperbolic fixed
points.
DEFINITION 6.7. Let p be a periodic point of the differentiable function f
with prime period k. Then p is a hyperbolic periodic point if l(fk)'(p)1 =I- l.
If l(fk)'(p)1 = 1, then p is a nonhyperbolic periodic point.
Notice that if p is a periodic point of f with period k, then it is the
derivative of fk at p that is important in determining whether or not p is
hyperbolic rather than the derivative of f at p. In exercise 6.3 we show
that if the prime period of p is larger than 1, then the derivative of f at p
by itself is not sufficient to determine the behavior of points near p under
iteration of the function.
As a corollary to Theorem 6.5 we get
THEOREM 6.S. Let f be a C 1 function and p be a per-iodic point of f
with prime period k. If l(fk)'(p)1 < 1, then there exists an open interval
containing p that is contained in the stable set of p. If I(fk)'(p)1 > 1, then
there exists an open interval containing p such that all points in the inter-val
except p must leave the inter-val under itemtion of fk.
Theorem 6.S is an easy extension of Theorem 6.5 and its proof is left as
an exercise.
We have observed that whenever p is a periodic point of f with prime
period k and the derivative of fk at p is less than one in absolute value,
then points near p are attracted to p under iteration of fk. On the other
hand, if p is a periodic point of f with prime period k and the derivative
54 6. Differentiability and Its Implications
Exercise Set 6
6.1 Find the periodic points of each of the following functions. Indicate
whether the periodic points are repelling, attracting or neither and
find their stable sets.
a) f(x) = -x 3
b) p(x) =x3 -x
c) k(x) = -x3 - x
d) g(x) = e x - 1
e) E(x) = e- X
f) S(x) = sin x
Exercise Set 6 55
g) r(x) = _x 1 / 3
h) A(x) = ~4 arctan x
6.4 Find a fixed point of I(x) = -!x + 3 and use the Mean Value
Theorem to prove that every real number is in the stable set of
this fixed point.
6.5 Prove the second half of Theorem 6.5 using the steps outlined be-
low. The proof depends on finding an € > 0 and an open interval
(p - 8,p + 8) so that 1!,(x)1 > 1 + € whenever x is in (p - 8,p+ 8).
!
a) Let € = (I!' (p) I - 1). Choose an appropriate 8 and use the
Triangle Inequality to show that 1/'(x)1 > 1f'(p)l- € whenever x
is in (p - 8,p + 8).
b) Prove 1f'(p)l- € = 1 + € and conclude 1/'(x)1 > 1 + € whenever
x is in (p - 8,p + 8).
c) Show that if y =1= p and Ik(y) is in (p - 8,p + 8) for all k < n,
then Ir(y) - pi> (1 + €)nly - pI.
d) Use (c) to prove the result.
6.9 Let Xo and a be positive real numbers and consider the sequence
defined by
Xn+l = ~2 (xn+~).
Xn
a) Find a function 1 such that Xn = r(xo).
b) Describe the periodic points of I. Are they hyperbolic? attract-
ing or repelling?
c) What is the limit of the sequence xo, Xl, X2, •.• ?
Note: 1 defines a very useful algorithm for computing square
roots as it has a strongly attracting fixed point. That is, the se-
quence converges quickly to the desired square root. For a discus-
sion of the error estimate see exercise 16 in Chapter 3 of the third
edition of W. Rudin's text Principles 01 Mathematical Analysis.
7
Parametrized Families of
Functions and Bifurcations
EXAMPLE 7.l.
Consider the parametrized family fm(x) = mx. For all values of m
except 1, it is clear that 0 is the only fixed point of fm.
If m < -1, then 0 is a hyperbolic repelling fixed point and all other
points are in the stable set of infinity.
If m = -1, then all points except 0 are periodic points with prime
period 2.
When -1 < m < 1, 0 is an attracting fixed point and WS(O) = R
If m = 1, then all points are fixed points.
Finally, when m > 1, 0 is a repelling fixed point and all other points are
in the stable set of 00. D
58 7. Parametrized Families of Functions and Bifurcations
We note that the dynamics of the family fm(x) = mx are unchanged for
large intervals of parameter values. Then at a particular parameter value
(Le., -1 and 1), the dynamics change suddenly, after which they again
remain constant for a prolonged interval. We call these sudden changes in
dynamics bifurcations. Notice also 0 is the fixed point and that If'(O)1 =
1 when m = 1 or -1. The presence of nonhyperbolic fixed points for
parameter values at which there is a bifurcation is typical.
DEFINITION 7.2. Let fc(x) be a parametrized family of functions. Then
there is a bifurcation at Co if there exists E > 0 such that whenever a and b
satisfy Co - E < a < Co and Co < b < Co + E, then the dynamics of fa (x) are
different from the dynamics of fb(x). In other words, the dynamics of the
function changes when the parameter value crosses through the point Co.
Readers who are used to a more rigorous approach may wonder what
we mean by saying "the dynamics are the same" or the "dynamics are
different". We will rely on an intuitive idea of sameness until we discuss
topological conjugacy in Chapter 10. Then we will say two functions have
the same dynamics if and only if they are topologically conjugate. Inter-
ested readers are encouraged to look ahead. For now, our intuition will be
able to deal well with the functions we examine.
Let's look at a few more simple examples.
EXAMPLE 7.3.
Consider the family of functions Ec(x) = eX +c. We examine the dy-
namics of this family by using graphical analysis. We are all familiar with
the general shape of the exponential function. We should notice that as
the parameter increases in value the shape of the graph remains the same,
but it moves to the left. A careful examination reveals that we have three
cases: the graphs of E(x) and y = x intersect in exactly two points, the
graphs of E(x) and y = x have a single point of tangency, and the graphs
of E(x) and y = x do not intersect. The graphs of typical examples are
shown in Figure 7.1.
If c < -1, then the graphs of E(x) and y = x intersect in two places
as shown in Figure 7.1. We call the two points of intersection a and b,
with a < b. Note that both a and b are fixed points, 0 < E'(a) < 1,
and E'(b) > 1. Hence, a is an attracting hyperbolic fixed point and b is
a repelling hyperbolic fixed point. Graphical analysis demonstrates that
WS(a) = (-00, b) and WS(oo) = (b,oo).
When c = -1, the two fixed points fuse into one fixed point at x = 1.
This fixed point is not hyperbolic since E~ 1 (1) = 1. Graphical analysis
demonstrates that WS(l) = (-00,1] and WS(oo) = (1,00). This case is
illustrated in Figure 7.1.
7. Parametrized Families of Functions and Bifurcations 59
4
The repelling [Link] poin! b
2 3
·1
·2
·2 4
-1
·2
a) b)
0.5
-\ -\
-\ -\
Finally, if c > -1, then the graphs of E(x) and y = x do not in-
tersect and consequently, E(x) has no fixed point. Since E(x) is con-
tinuous, Sarkovskii's Theorem implies that E(x) can have no periodic
points when c> -1. Using graphical analysis we show in Figure 7_1 that
WS(oo) =~. 0
EXAMPLE 7.4.
Consider the family of functions Ak(X) = karctan(x) for parameter
values near 1. When 0 < k < 1, the function has one attracting fixed point
near 0 and all real numbers are in the stable set of 0_ If k = 1, then 0
is still a fixed point but Ai (0) = 1. Graphical analysis indicates that 0
is weakly attracting and WS(O)is still ~ even though 0 is not hyperbolic.
Finally, when k > 1, we find a repelling fixed point at 0 (A~(O) = k > 1)
and two other fixed points have been formed, both of which are attracting.
If the two fixed points are labeled a and b and a < 0 < b, then we see that
WS(a) = (-00,0) and WS(b) = (0,00)_ The graph of each case is shown in
Figure 7.2 along with the line y = x. 0
7. Parametrized Families of Functions and Bifurcations 61
EXAMPLE 7.5.
The function hr (x) = rx( 1- x) exhibits a third type of bifurcation when
r = 1: a trans critical bifurcation. When 0 < r < 1, h has two fixed points,
one of which is less than one and repelling, and another at 0, which is
attracting. When r = 1, these two points merge .to form one fixed point at
o which attracts points which are greater than 0 and repels points which
are less than O. Finally, if 1 < r < 3, we see that 0 becomes a repelling
fixed point and there is an attracting fixed point which is greater than O.
The graphs of these three situations are shown in Figure 7.4. 0
EXAMPLE 7.6.
Continuing the investigation of hr(x) = rx(l - x) that we started in
Example 7.5, we find that hr undergoes another bifurcation when r = 3.
We have seen that hr has a repelling fixed point at 0 and an attracting
fixed point greater than 0 when 1 < r < 3. A quick look at the graph
of h; demonstrates that hr has no other periodic points when r is in this
range. At r = 3 the larger fixed point is at .75 and is weakly attracting:
it is no longer hyperbolic since h' (.75) = -1. When r > 3 the fixed point
is repelling. (The proof of this is left as an exercise.) Graphical analysis
when r = 3.1 indicates that most points in (0,1) are attracted to a period
2 orbit which straddles the larger fixed point. What has happened?
To understand this, we look at the graph of h;
in the vicinity of the
fixed point Pr = r-;l for parameter values slightly less than 3, equal to 3,
62 7. Parametrized Families of FUnctions and Bifurcations
j 1+
a) b)
rr
2
4
t I
I
2 0
ttl
0
-1
-2
1i~ 1t
-3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1.5 2
a) b) c)
-OJ 0.3
~ t t~~{
i i i'f)
0.5
O r-----------~~----------------------------------
·0.5
.1~-----L------~------------~2------------~-------
/
a) b) /
/
0.15 0.15
0.1 0.1
'''V
0.65
0.6 0.6
0.55 / 0.55
0.55 0.6 0.65 0.1 0.15 •
0.55 0.6 0.65 0.1 0.15
c)
0.15 Q2
\
0.1
0.6S
-----
new allracLing period 2
0.6
0.55 Q.
0.55 0.6 0.65 0.1 0.15
and slightly larger than 3. Some examples are shown in Figure 7.6. Notice
that as the parameter value passes 3, the value of (h~)'(Pr) changes from
being less than 1 to being greater than 1. Hence, the fixed point changes
from attracting to repelling. In addition, the continuity of hr requires that
as this happens a period 2 attracting orbit must be added. (The details
of the proof of this fact are also left as an exercise.) This orbit is labeled
by the points al and a2 in Figure 7.6. This type of bifurcation is called a
period doubling bifurcation since a periodic orbit of twice the period of the
original periodic point is added. A bifurcation diagram for the family hr
is shown in Figure 7.5. 0
Exercise Set 7
7.3 In Example 7.6 we stated that a period two orbit must be added
as the derivative of h~ at Pr = r-;:l changed from being less than
one to being greater than one. Prove that this is true.
Hint: First show that (h~(Pr))' > 1 implies there exists Xo > Pr
so that h~(xo) > Xo. Then use the fact that h~(l) = 0 and the
Intermediate Value Theorem to show that there must exist a fixed
point of h~ in the interval (Pr, 1). Using a similar technique show
there exists a fixed point of h~ which is in (!, Pr). How do we know
these points are prime period 2 points of h r ? Can you prove that
they are in the same orbit?
°
Note that the stable set of is still {O, I}, but the stable set of
Pr contains an infinite number of points. However, there is now
an attracting period 2 orbit and "most" of the points in (0,1) are
forward asymptotic to one of the two points in the orbit.
(7) When r ~ 3.45 another period doubling bifurcation occurs and the
period 2 attracting orbit splits into a period 4 attracting orbit and
a period 2 repelling orbit.
(8) In the interval 3.44 < r ~ 4 changes occur rapidly. There is a pe-
riod 3 orbit whenever r > 3.6, so by Sarkovskii's Theorem we know
there are periodic points of all orders. Because of the complexity
of the dynamics when r is in this range we will first discuss what
happens when r > 4.
Suppose r > 4. Then h( ~) > 1 and ~ is forward asymptotic to infinity.
Since h(~) > 1, h(l) = 0, and h(O) = 0, the Intermediate Value Theorem
°
as an exercise the proof that there are infinitely many points which are
eventually fixed at and that there are infinitely many points which are
eventually fixed at Pro We would like to know what other points of [0,1]
are forward asymptotic to something in [0,1].
For each natural number n we define
n=l
An: the set of points that remain in
( ~2 _ vr2 - 4r ~
2r' 2 +
vr2 - 4r)
2r
8. The Logistic Function, Part I 69
and
A= [0 ~2 - vr22r- 4r] U[~2 + vr22r- 4r ' 1].
1 ,
° °
one on the intervals of AI, we note that it is strictly monotone on these
intervals since h'(x) = r(l - 2x) > when x < ~ and h'(x) < when
x > ~. Therefore, if I is one of the two intervals in AI, then h : I ----+ [0,1]
is one-to-one and onto.
To continue the induction argument, we suppose that Ak consist of 2k
disjoint closed intervals. Further, we suppose that if I is one of the closed
intervals in Ak , then hk : I ----+ [0,1] is one-to-one and onto and that either
° °
(hk)'(x) < for all x in I or (hk)'(x) > for all x in [a, b].
Now consider Ak+1' It is easily seen that A k+1 C Ak; the reader should
°
be able to explain why this is so. Assume [a, b] is one of the intervals in Ak
and (hk)'(x) > for all x in [a,b]. Since hk is strictly increasing on [a,b],
hk is continuous, and hk([a, b]) = [0,1], the Intermediate Value Theorem
implies there exist points X2 and X3 such that
(1) a < X2 < X3 < b,
(2) h k ([a,x2]) = [0, ~ _ v'r~;4r],
(3) hk((x 2, x))
3
-- (1.2 _ v'r 2
2 -4r 1.
r' 2
+ v'r2
2-4r) and
r'
(4) hk([X3, b]) = [~ + v'r~;4r, 1].
The first condition implies that the intervals [a, X2] and [X3, b] are disjoint.
The last three conditions imply that hk+l ([a,x2]) = [0,1]' h k+1(x) > 1
70 8. The Logistic Function, Part I
___ 1 __________ _
for all x in (X2, X3), and hk+1([X3, b]) = [0,1]. Thus, the set of points in
[a, b] which are also in A k + 1 consists of the two disjoint closed intervals
[a,x2] and [x3,b]. If x is in [a,x2]' then hk(x) is in [0, ~ - v'r~;4T] and
h'(hk(x)) > O. By assumption (hk)'(x) > 0, so for all x is in [a, X2],
(h k+1)'(x) = h'(hk(x))· (hk)'(x) > O.
A similar argument demonstrates (h k+1)'(x) < 0 for all x in [x3,b].
If we begin the argument in the preceding paragraph by assuming that
(hk)'(x) < 0 for all x in [a, b], then we can use an analogous argument to
show that the points of A k + 1 which are in [a, b] are contained in two disjoint
intervals [a,x2] and [.T3,b], hk+1([a,x2]) = [0,1], hk+1([X3,b]) = [0,1], and
that (h k + 1 )' is either strictly positive or strictly negative on each of the
intervals [a, X2] and [X3, b]. Since [a, b] is an arbitrary interval in Ak it
follows that there are twice as many intervals in A k + 1 than there are in A k .
That is, Ak+1 contains 2(2k) = 2k+1 disjoint closed intervals. Further, we
have shown that if J is one of the intervals in A k+1, then hk+1 : J -'> [0,1]
is onto (by conditions (2) and (4) listed above). Also, hk+1 is one-to-one
on J and we can make the next step in the induction since (h k + 1 ), is either
strictly positive or strictly negative on J. D
n=l
An. In the
proof of Proposition 8.1 we saw that Ak+1 was formed by removing an open
interval from the middle of each interval in Ak . The set of points which
are left after this has been done for each n is A. Students of topology
may recognize that this construction is similar to that of a Cantor Middle
Thirds Set.
8. The Logistic Function, Part I 71
EXAMPLE 8.3.
CONSTRUCTING CANTOR SETS.
Cantor sets are usually constructed by iterative processes. To construct
the Cantor Middle Thirds Set we begin with the interval [0,1] and remove
the open set (!, ~). In the second step we remove the middle third of each
remaining interval. That is, we remove the intervals (!,~) and (~, ~). We
continue the process by removing the middle third of each of the remaining
intervals at each step. The set of points remaining is called the Cantor
Middle Thirds Set. As we will see shortly, this is in fact·a Cantor Set.
H 0 < a < 1, then we can construct similar sets called Cantor Middle-a
sets by removing an interval from the center of each remaining interval
whose length is a times the length of the remaining interval. For example,
in the Middle-i set we remove the middle fifth of each interval at each step.
To be more precise, Let ro = [0,1]. Let r l be the two closed inter-
vals of equal length left when the open interval of length a is removed
from the middle of r o. Define r 2 to be the set of 4 closed intervals
of equal length obtained by removing an open interval whose length is
a x (the length of an interval in r d from the middle of each interval in
r l . We continue to define r n inductively as the set of closed intervals
of equal length formed by removing an open interval whose length is a x
(the length of an interval in r n-I) from the center of the closed intervals
in r n-I. The first few iterations of this construction for a = 3 are shown
in Figure 8.2. The Cantor Middle-a Set is the set r = n~o r n.
Clearly, r is not empty since the endpoints of the open intervals removed
are in r. We shall see in Proposition 8.5 that r is perfect. Students who
have studied count ability will be interested to know that nonempty perfect
subsets in lR must be uncountable. This is proven in the analysis text by
Rudin which is listed in the references. It is also interesting to note that
72 8. The Logistic Function, Part I
113 2/3
only count ably many of the points in r are endpoints of the intervals that
are removed, so there are uncountably many points in r which are not
endpoints. 0
Note that the proof of Proposition 8.5 does not use the fact that we
removed the a-interval from the exact middle of the closed interval. The
fact that the total length of all the intervals approaches zero follows from
the observation that if we remove a fraction a at each stage, then the
length is reduced by a factor of 1 - a so the total length of the remaining
intervals is (1- a)n after n steps. In fact, if at each step we remove an open
interval whose length is a times the length of each closed interval and we
remove this interval in such a way that the endpoints of the open interval
do not coincide with either of the endpoints of the closed interval, then we
obtain a Cantor set. It can also be shown that if a subset of the real line
is closed and bounded, contains no open intervals, and is perfect, then it is
homeomorphic to the Cantor Middle Thirds Set.
We complete this chapter by showing that the set of points which remain
in [0,1] under iteration of hr is a Cantor set whenever r > 2 + J5. This is
true for all r > 4, but it is a bit trickier to prove for r < 4 ::::; 2 + J5. To
do this we use the following lemma:
74 8. The Logistic Function, Part I
LEMMA 8.6. Let r > 2 + v'5 and h(x) = rx(l- x). Then there is € > 0
such that Ih'(x)1 > 1 + € whenever x is in A 1 . Further, the length of each
interval in An is less than (l;€)n.
Note that the value of € found in the lemma depends on the value of r.
As r gets smaller € must get smaller as well.
The proof of the first part of the lemma follows from the fact that if
x is in Al. then Ih~(x)1 is greater than or equal to the absolute value of
the derivative of hr at (~ ± v'r~;:-4r). We prove the second portion of the
lemma by using the first part and the Mean Value Theorem. The proof is
outlined in exercise 8.3.
n=l
An is a Cantor set.
Exercise Set 8
8.1 Show that h(x) = rx(l - x) has infinitely many eventually fixed
points when r ?: 4.
8.4 a) Use Proposition 8.1 and Theorem 4.4 to show that if r > 2+V5,
then hr(x) = rx(l - x) has 2n periodic points with period n.
b) Use Lemma 8.6 and exercise 6.2 to show that the periodic points
of from part (a) are all repelling.
c) Show that if r > 2 + V5, then the periodic points of hr are dense
in A.
8.7 • THE TENT MAP: Let T: [0,1]- [0,1] be the tent map defined
by
T(x) = {2X' for x ~n [~,~]
2 - 2x, for x 10 [2,1].
a) Graph T(x), T2(x), T3(x), and T 4 (x). (You may wish to use
a computer graphing tool. See the appendix for suggestions as to
how one can define T for the computer.)
Exercise Set 8 77
DEFINITION 9.1. Let 2::2 be the set of all infinite sequences of O's and
1 'So This set is called the sequence space of 0 and 1 or the symbol space of
o and 1. More precisely, 2::2 = {(SOS1S2 ... ) I Si = 0 or I}. We will often
refer to elements of 2::2 as points in 2:: 2 .
drs, t] = f
i=O
1
8i ; til.
2i =2
i=O
80 9. Symbolic Dynamics and Chaos
DEFINITION 9.3. Let X be a set and let d be a function from the set of
all ordered pairs of elements of X into the real numbers. If the following
conditions hold for all x, y, and z in X, then d is a metric on X.
a) d[x, y] ~ 0 and d[x, y] = 0 if and only if x = y.
b) d[x,yj = d[y,x].
c) d[x,yj ~ d[x,zj +d[z,yj.
We are already familiar with one metric: the distance between two points
in the real numbers is the absolute value of their difference. For example
the distance between the numbers 2 and 3 is \2 - 3\ = 1 and the distance
between 2 and -5 is \2 - (-5)\ = 7. It should be clear that this distance
function satisfies conditions (a) and (b) of the definition. The reader may
recognize condition (c) in Definition 9.3 as the Thiangle Inequality, which
was proven for the metric on the real numbers in exercise 2.8. In the
exercises at the end of this chapter we show that the distance function in
symbol space is a metric and evaluate several other functions as possible
metrics.
The obvious question now is "Why do we care?". If we refer back to
Chapter 3 we see that all of our topological definitions and the definition
of continuity are based on the idea of distance. So once we have defined a
metric on a set we can define a topology for the set. That is, we can decide
which subsets are open, which are closed, what a convergent sequence is
in the set, and which functions of the set are continuous. All of these
definitions are brought together below.
drs, t) = f
i=O
lSi; til
n 0 00
=L 2i+ L
i=O i=n+l
EXAMPLE 9.12.
a) The distance between 8 = 000000 ... and t = 010101010 ... is
b) If S is the set of all elements of ~2 which begin with the sequence 011,
then S is closed.
To demonstrate this we show that all of the accumulation points of S are
contained in S. If 8 is an accumulation point of S, then every neighborhood
of 8 contains an element of S other than 8. Consequently, there is 8* in
N 1 (8) n S satisfying 8* =I- 8. But then the definition of neighborhood
~
implies d[8*, 8] < d-:> and Lemma 9.11 implies that the first 3 digits of 8 and
8* must agree. Since 8* is in S, the first three digits of 8* are 011. Hence,
the first three digits of 8 are 011 and 8 is in S.
c) The set of all sequences in ~2 which end with an infinite string of O's is
dense in symbol space. That is, if
(J'(808182 .. ·) = 81 8 2 8 3 ....
In other words, the shift map "forgets" the first digit of the sequence. For
example, (J'(01110101 ... ) = 1110101 ....
84 9. Symbolic Dynamics and Chaos
There are several immediate observations we can make about the shift
map.
PROPOSITION 9.14. The shift map is continuous.
PROOF. Let s be an element of E2 and f > O. We must show there is
8 > 0 such that whenever drs, t] < 8, then d[a(s),a(t)] < f. Choose n so
that 2~ < f and let 8 = 2n\1' If drs, t] < 8, then we know from Lemma
9.11 that sand t agree on the first n + 2 digits. So a(s) and a(t) agree on
the first n + 1 digits and d[a(s), a(t)] S 2~ < f as desired. 0
It should be clear by now that two points are close together if and only
if their initial digits agree. The more digits they agree on before they differ
the closer together they are. Of course, this is exactly what Lemma 9.11
claims. Example 9.12 and the proof of Proposition 9.14 illustrate the use
of this idea. We will use this idea once again in the proof of part (a) of
Proposition 9.15 and the reader will be given an opportunity to use it in
proving parts (c) and (d).
PROPOSITION 9.15. The shift map has the following properties:
a) The set of periodic points of the shift map is dense in E 2 •
b) The shift map has 2n periodic points of period n.
c) The set of eventually periodic points of the shift map which are not
periodic is dense in E 2 •
d) There is an element of E2 whose orbit is dense in E 2. That is,
there is s* in E2 such that the set {s*, a(s*), a 2(s*), a 3 (s*), ... } is
dense in E 2 .
e) The set of points which is neither periodic nor eventually periodic
is dense in E 2 •
PROOF. a) Suppose 8 = 808182 ... is a periodic point of a with period
k. Then an(ak(s)) = a n(8). Since an(s) "forgets" the first n digits of 8 we
see that
the initial k digits infinitely often. But, if t = totlt2t3 ... and we choose n
so that in < €, then we can let s = totl ... tntotl ... tntotl . . .. As t and s
agree on the first n + 1 digits, Lemma 9.11 implies drs, tj :::; 2~ < € or s is
in NE (t) and s is a periodic point by construction.
The proofs of parts (b) and (c) are left as exercises.
d) The sequence which begins with 0 1 00 01 10 11 and then includes
all possible blocks of 0 and 1 with three digits (there are eight) followed
by all possible blocks of 0 and 1 with four digits (there are sixteen) and so
forth is called the Morse Sequence. We leave as an exercise the proof that
the orbit of this sequence is dense in ~2.
e) Since the set of non-periodic points includes as a subset the orbit of
the Morse Sequence defined above, the truth of part (e) of this proposition
follows from part (d). D
The shift map has several characteristics worth highlighting. First we
note that the periodic points of a are dense in ~2' This was proven in
the previous proposition. This implies that in every neighborhood of every
point in ~2 there is a periodic point of a. Further, by Proposition 9.9, if
s is in ~2' then there is a sequence of periodic points Sl, S2, S3, •.. which
converges to s. We can think of this property as guaranteeing a certain
amount of regularity in the dynamics of the shift map. Periodic points of
a are everywhere in the domain.
Second, we note that ~2 is well mixed by a. We define this mathemati-
cally as topological transitivity.
DEFINITION 9.16. The function f : D ---> D is topologically transitive if
for all open sets U and V in D, there is x in U and a natural number n
such that r(x) is in V.
We sometimes describe the action of a function which is topologically
transitive as mixing the domain. If U is any open set in the domain of the
function, then some point of U will eventually land in every neighborhood
of every point in the domain under iteration of the function. Proposition
9.15d in conjunction with the following proposition guarantees that the
shift map is topologically transitive.
PROPOSITION 9.17. Let f be a function and suppose there is a point
whose orbit under iteration of f is dense in the domain of f. Then f is
topologically transitive.
PROOF. Let f : D ---> D and let U and V be open sets in D. We
must find x in U and a natural number n such that fn(x) is in V. By
hypothesis, there is a point in D, whose orbit is dense in D. If Xo is such
a point, then we know there exists k such that fk(xO) is in U. If we can
86 9. Symbolic Dynamics and Chaos
E = min{lv - r(xo)l}
°
is greater than and the neighborhood N€/2 (v) doesn't contain any it-
erates of Xo. Thus, the orbit of Xo is not dense in D, a contradiction.
Consequently, V must contain infinitely many iterates of Xo. Since there
are only finitely many natural numbers less than k, there exists m > k such
that fm(xo) is in V. If we let n = m - k, the proof is complete. 0
d[an(t),an(s)] = f
i=O
ISHn; tHnl = f
i=O
~ = 2.
Functions which exhibit the three characteristics just outlined for a are
said to be chaotic.
DEFINITION 9.19. (DEVANEY!) Let D be a metric space. Then the
function f : D -+ D is chaotic if
a) the periodic points of f are dense in D,
b) f is topologically transitive, and
c) f exhibits sensitive dependence on initial conditions.
Restated, this definition implies that a chaotic function observes a cer-
tain amount of regularity and mixes the domain well. Further, even the
smallest changes in initial position may result in dramatically different re-
sults in values under iteration of the function. The regularity is evidenced
by the fact that we can always find a periodic orbit in every neighborhood,
no matter how small, of every point in the domain of the function. On the
other hand, the domain is well mixed by the function since if we choose
any open set, then we can find a point in every other open set that will
eventually end up in the first set under iteration of the function. As we
have seen, a : ~2 -+ ~2 exhibits all of these behaviors and is chaotic.
We complete this section with a new and surprising result which ap-
peared in the April 1992 issue of the American Mathematical Monthly
(pages 332-334) and which was proven by Banks, Brooks, Cairns, Davis,
and Stacey.
THEOREM 9.20. Let f : X -+ X be topologically transitive and suppose
the periodic points of f are dense in X. If X contains an infinite number
of elements, then f exhibits sensitive dependence on initial conditions.
This result implies that when the periodic points of a function are dense
in its domain and the function is topologically transitive, then the function
is chaotic. It is not necessary to check for sensitive dependence on initial
conditions. The assumption that X contain an infinite number of elements
is added to eliminate trivial cases. A function of a set containing only one
element will satisfy the hypotheses ofthe theorem, but the conclusion can't
hold. Counter examples exist for larger finite sets as well.
PROOF. Let f : X -+ X be topologically transitive and suppose the
periodic points of f are dense in X. Let d be the metric on X. We begin
by showing there is 80 > 0 such that for all x in X there is a periodic point
q such that the distance from x to r(q) is greater than or equal to 80 for
all n. In other words, every x in X is at least 80 away from a periodic
IThere are several definitions of chaotic. This definition was introduced by R. De-
vaney in his text, Introduction to Chaotic Dynamical Systems, which was published in
1989 and is listed in the references.
88 9. Symbolic Dynamics and Chaos
Since there are only finitely many points in the orbits of p and q, Do > O.
By definition, all of the points in the orbit of p are at least 2Do units away
from all of the points in the orbit of q. Then by the triangle inequality,
V = nk
i=O
f-i(No(fi(q)))
nf-
k
V = i (N6 (t(q)))
i=O
9. Symbolic Dynamics and Chaos 89
.\ .
d[x , Jkj-m(q)] :::; d[x ,p] + d[p, Jk j (y)] + d[Jkj (y), Jkj-m(q)]. (9.1)
Since we have assumed that d[x,p] < f < 8 and we have shown that
d[Jkj(y),Jkj-m(q)] < 8, inequality (9.1) implies
(9.2)
Exercise Set 9
9.2 Let C be the set of all ordered pairs ofreal numbers. (Le., C is the
Cartesian plane.)
• a) Show that dl [(XI,YI),(X2,Y2)] = V(XI-X2)2 + (YI-Y2)2 is
a metric.
Let y = (1, 2) be a point in C and graph the set
NI(y) = {x in C I dl(x,y) < I}.
b) Show that d2[(XI, YI), (X2' Y2)] = IXI - x21 + IYI - Y21 is a metric.
Let y = (1, 2) be a point in C and graph the set
NI(y) = {x in C I d 2 (x,y) < I}.
9.4 Prove Propositions 9.5, 9.6, 9.7, 9.8, 9.9, and 9.10.
Hint: Use the Triangle Inequality for Proposition 9.6.
Exercise Set 9 91
9.5 Let 8' be the set of all points which don't begin with OIl. Use
Definition 9.4a to prove that 8' is open. Note that we know this
is true since 8' is the complement of the closed set 8 in Example
9.12b.
9.6 Let r = 101010 ... and N~(r) = {s I d[s,r] < n. Describe the set
of points which are in Nl (r).
2
k=l
is an open set.
9.9 Find all the points of period three in I: 2. Which of these points
have prime period three? Which are in the same orbit? Does
Sarkovskii's Theorem apply here? Why or why not?
9.11 Suppose m > 1 and let f : JR. -+ JR. be defined by f(x) = mx.
Show that f exhibits sensitive dependence on initial conditions. Is
f chaotic? Explain.
9.12 • THE DOUBLING MAP ON THE CIRCLE: Let 8 1 be the unit circle
and identify each point on the circle by the radian measure of
the angle between the positive x-axis and the ray beginning at
the origin and passing through the point. We will always measure
angles in a counterclockwise direction. The reader should recognize
this as the usual representation of the unit circle and angles on it.
For example the point (1,0) is labeled 0, and the point (0,1) is
labeled ~. We will also find it useful to identify the point a with
the point oo+2mf where n is any integer. For example, the points 0,
21f, -21f, 41f, ... are all identical. Also, i, 737T , and - 5; are all the
same point on the circle. Again, this is the standard representation
of angles on the unit circle.
We define a metric on 8 1 by letting d[oo,,6] be the length of the
shortest arc on the circle from a to,6. More precisely, if a and ,6
92 9. Symbolic Dynamics and Chaos
9.13 • THE TENT MAP REVISITED: Determine whether or not the tent
map defined by
is chaotic.
Hint: It may help to review the work from exercise 8.7.
We return now to the analysis of the logistic function. Our goal is to prove
that if r > 2 + yI5, then hr(x) = rx(1- x) is chaotic on A. We recall that
A is the set of all numbers in [0,1] which remain in [0,1] under iteration
of h. That is, A = {x I hn(x) is in [0,1] for all n}. By Theorem 9.20, it is
sufficient to show that the periodic points of h are dense in A and that h
is topologically transitive on A. We have already shown that the periodic
points of h are dense on A in exercise 8.4. Unfortunately, proving that h
is topologically transitive on A directly from the definition is a relatively
difficult task. Consequently, we will show instead that the dynamics of h
on A are the same as the dynamics of a on L: 2 • Mathematically speaking,
we say that h on A is topologically conjugate to a on L: 2 .
f
x ~ f(x)
or TilT
E~E r(x) ~ g(r(x)).
g g
96 10. The Logistic Function, Part II: Topological Conjugacy
x y
or
E ----7 E T(X) ----7 g(T(X)).
9 9
f T-l(a) f f(T-l(a))
D ----7 D ----7
r-11 Ir or r-11 Ir
E ----7 E a ----7 g(a).
9 9
UeE VeE
With a little practice this diagram chasing becomes second nature and an
easy way to complete proofs involving commuting functions like topological
conjugacies.
PROOF OF THEOREM 10.3:
a) It suffices to show r- 1 0 9 = for-I. The details are left to the reader.
b) To prove that r 0 fn = gn 0 r we need only consider the following
commutative diagram consisting of n squares. As each of the squares com-
mutes, the whole diagram commutes. Thus, we get the same result by
going around the outside of the diagram in either direction and the truth
of the assertion follows.
D~D~D~ ~D
Alternatively, for those with a more analytical bent, we may use the
Principle of Mathematical Induction. The statement r 0 fn = gn 0 r is true
when n = 1 since r is a topological conjugacy. Now suppose it is true for
n, that is, r 0 fn = gn 0 r. Then it is true for n + 1 since
r 0 r+ 1 = (r 0 r) 0 f = (gn 0 r) 0 f = gn 0 (r 0 f) = gn 0 (g 0 r) = gn+1 0 r.
c) Let p be a periodic point of f with prime period k. We first show
gk(r(p)) = r(p). Since fk(p) = p, we get the commutative diagram
p p
show that for each E > 0, there is M such that whenever n ~ M, then
dE [gkn(r(x)), r(p)] < E.
Suppose E > O. Since r is continuous, there exists 8 > 0 so that if
dD[Y,p] < 8, then dE[r(y), r(p)] < E. Choose M so that if n ~ M, then
dD[fkn(x),p] < 8. Then by continuity,
when n ~ M.
e) The preservation of density of periodic points by a topological conju-
gacy follows easily from part (c) of this proposition and part (d) of Propo-
sition 10.2.
f) The argument for this portion is outlined above.
g) This portion follows from parts (e) and (f) and Theorem 9.20. 0
We complete this chapter by constructing a topological conjugacy from
hr : A -+ A to a : E2 -+ E2 when r > 2 +.J5. As a is chaotic on E 2 ,
Theorem 10.3 then implies that h is chaotic on A.
We begin by defining the function 'I/J : A -+ E 2 • Set 10 = [0, ~ - v'r~;4r]
and h = [~ + jr~;4r, 1]. Recall from Proposition 8.1 that
and
So, if Sn = 0, then Isosl ... Sn = [a, Cl] and if Sn = 1, then Isosl",Sn = [C2, b].
As this is the same way we derived the intervals of An+! in the proof of
Proposition 8.1, we see that Isosl ... Sn is one of the intervals comprising
A n +l .
a) Let S = SOSlS2S3 ... be an element of E 2 . To show that 'lj; is one to
one and onto we must show that 'lj;-l(S) contains exactly one point. But if
x is in 'lj;-l(s), then x is in Isosl ... Sn for all n. Thus,
'lj;-l(s) = n
00
n=O
Isosl",Sn'
We need to show this intersection is not empty and contains exactly one
point. Let Isosl",Sn = [an, bn]. Suppose 'lj;-l(s) contains two points x and
y. Then Ix-yl ~ Ibn -ani for all n since x and yare in each of the intervals
10. The Logistic Function, Part II: Topological Conjugacy 101
[an, bn ]. However, Lemma 8.6 implies that Ian - bnl approaches 0 as n goes
to infinity. So it must be that Ix - yl = 0 and x = y. Consequently, we see
that 'lj;-I(S) contains no more than one point.
It remains to show that 'lj;-I(s) contains at least one point. This follows
immediately from the well known Nested Interval Theorem of mathematical
analysis which states that the intersection of a nested set of closed intervals
is not empty. (See Theorem 2.38 of W. Rudin's Principles of Mathematical
Analysis, 3rd edition.) Recall that the intervals 1SQS1, .. Sn are denoted by
[an, bn ]. The proof ofthe Nested Interval Theorem depends on the principle
that there exists a unique real number which is larger than or equal to
every an and smaller than any other number which is also greater than
all the an. Such a number is called the least upper bound of the set
{al,a2, ... ,an , ... }. Now since an ::; am ::; bm whenever m ~ nand
an ::; bn ::; bm whenever m ::; n, we see that each bm is greater than
every an. So the least upper bound, a, is less than or equal to every bm .
Therefore, an ::; a ::; bn for each n and a is in every [an, bn] = 1SQS1 ... Sn' It
follows that a is in 'lj;-I(s) and we have completed the proof of part (a).
b) Let E > 0 and x be in A. To show 'lj; is continuous at x we must find
8> 0 so that if Ix - yl < 8, then d['lj;(x), 'lj;(y)] < E. If we choose n so that
2~ < E, then Lemma 9.11 implies that it is sufficient to show there exists
a 8 > 0 so that the sequences 'lj;(x) and 'lj;(y) agree on the first n + 1 digits
whenever Ix - yl < 8. In other words, if'lj;(x) = SOSIS2 .. " then x is in
1SQ81S2",Sn and we need to show that y is in1sQs1s2",sn whenever y is in A
and Ix - yl < 8.
Let [a b b1], [a2,b 2], ... , [a2n'H,b2n+1] be the 2n+1 intervals in An+l
!
indexed so that bi - 1 < ai for all i. Set 8 = min{lai - bi - 1 1}. Since the
intervals are not overlapping and finite in number, 8 is positive. Also, if
x and yare in An+l and Ix - yl < 8, then it follows that x and yare in
the same interval of An+1. As 1SQS182 ... 8n is an interval in An+1 and x is in
18Q8182 ... 8n' we see that y must be in 18Q8182 ... 8n and the proof of part (b) is
complete.
c) The proof that 'lj;-1 is continuous is,1eft as an exercise.
d) It remains to show 'lj;oh = ao'lj;. Let x be in A and 'lj;(x) = SOSIS2 ....
As we demonstrated in the proof of part (a), x is the unique point in
18Q8182'" = n00
n=O
18Q8182 ... 8n =
n=O
n
00
n n
00 00
1818283 .• ' = 18182 ... 8n = {x I hk+1(x) is in 18k for all k ::; n}.
n=1 n=O
102 10. The Logistic Function, Part II: Topological Conjugacy
Then h is chaotic on A.
The proof follows from Theorems 10.3 and 10.4 and Proposition 9.15.
Exercise Set 10
10.5 • THE TENT MAP: THE ADVENTURE CONTINUES. Define the tent
map by
T(x) = {2Xl for x in [O,~]
2 - 2x, for x in [~, 1].
a) Show that the map 7(X) = ~ cos(7r(l-x)) + ~ = cos2(~(1-x))
is a topological conjugacy between the map h(x) = 4x(1 - x) on
[0, 1] and the tent map.
Exercise Set 10 103
b) Use part (a) of this exercise, exercise 9.13, and exercise 8.7 to
show that h: [0,1]-+ [0,1] is chaotic.
o·
2
.
3 4
r-values
0.8
fixed point of h~ in the interval (xo, 1). This point is a prime period 2
point. This was proven in exercise 7.2.
Now consider the graph of h~(x) for r = 3.4 shown in Figure 11.2. In
particular, notice the boxed region. If we rotate the graph of h~ on the
interval over which the box sits, then it closely resembles the graph of hs
in the unit square when s = 2.646. By the unit square we mean the box
with corners at (0,0), (0,1), (1,1), and (1,0).
We can draw a square box like the one on the graph of h~ in Figure 11.2
for all parameter values r satisfying 2 < r :s 4. In each case the box has
corners at (!r'r!) ' (!r' r-l)
r '
(r-l r-l)
r' r '
and (r-l
r 'r
!). Note that when r = 2
the upper-right corner of the box coincides with the lower-left corner and
the box is no longer defined. If h~(!) = 2 r (i;;r)is greater than ~, then h2
maps the interval Ir = [~, r;l] to itself. Solving numerically, we see this
occurs when r is less than approximately 3.678.
In Figure 11.3 we consider the graphs of h~(x) for r = 3.4, r = 3.45,
and r = 3.5. In each case the graph is shown on the interval I r . Note that
the interval changes as we change r, but the behavior is identical to that
which we observed when we looked at the graph of h~(x) for r near 3 in
Example 7.6; a period doubling bifurcation occurs. If p is the period two
point, then (h~)'(P) < 1 when r < 3.45. If r = 3.45, then (h~)'(P) ~ 1 and
when r > 3.45, a new attracting period 4 point has appeared.
Looking at the graph of h~ for r = 3.5, we again notice square boxes
in which the graph of h~ is similar to the graph of hs in the unit square.
These boxes are identified in Figure 11.4. This process can be continued
indefinitely. If hr has a periodic point with prime period n, then we can
find regions on the graph of h~ around which we can draw a square box so
that the graph of h~ in the box is similar to the graph of hs on the unit
square. As r changes we will find all of the same dynamics for h~ within
these boxes that we find for hs on the unit square. In particular, since
hs has a period doubling bifurcation when s = 3, h~ will have a period
doubling bifurcation for some value of r.
Recall that h~ : Ir -+ Ir where Ir = [~, r;l] whenever r < 3.678. When
r > 3.679, the graph of h r "sticks out" of the box. That is, there are points
in Ir which are not mapped back into I r . This situation is illustrated in
Figure 11.5. Note that the graph of h~ on Ir for r > 3.679 is similar to the
graph of hs on [0,1] for s > 4 with one significant difference; points which
leave the interval Ir under iteration of h~ when3.679 < r < 4 may return
to the interval, whereas points which leave [0,1] under iteration of hs are
asymptotic to infinity. In Figure 11.5 we show via graphical analysis of h~
that though! leaves Ir when r = 3.8, it returns within a few iterations.
Recall from Chapter 8 that if s > 4, then hs is chaotic on a Can-
tor set which is contained in [0,1]. A similar result holds for h~ on Ir
11. The Logistic Function, Part III 109
0.3
~0~.3-----0-.4----~0.-5-----0.~6----0~.7-
0.3
0.3 0.4 0.5 .6 0.7
O. 0.8
0.6
0.4
0.2
Exercise Set 11
r-values 4
.6
so, what is the period of the new attracting fixed point? At what
parameter value does it appear? Support your statements with
graphs and explanations.
c) Use Feigenbaum's constant, the values Cn, and the fact that
00 1
"'" xn = - - to estimate the value of r at which the period dou-
L I-x
n=O
bling cascade will have been completed.
114 11. The Logistic Function, Part III
11.9 What causes the shadow curves seen in the bifurcation diagrams
of Figure 11.6?
Note: An interesting discussion of the shadow curves can be
found in section 11.4 of the book Fractals and Chaos by Peitgen,
Jurgens, and Saupe, which is listed in the references.
12
Newton's Method
The search for solutions of the equation f (x) = 0 is ancient and methods
that can solve equations of the form ax2 + bx + c = 0 are several thousand
years old. In the 16th century, Italian mathematicians discovered methods
for solving third and fourth degree polynomials. However, it was shown in
the early part of the nineteenth century that there is no general method
for solving polynomials of degree five or higher. Consequently, methods for
estimating solutions of equations as simple as polynomials are necessary.
Isaac Newton developed such a method, which was later refined by Joseph
Raphson, and which we now know as Newton's method or the Newton-
Raphson Method. Newton's method is easy to use and is often taught in
first semester calculus since it only requires knowledge of the derivative.
To illustrate Newton's method, suppose f (x) is differentiable and that Xo
is a reasonable approximation of a solution to the equation f(x) = O. How
can we find a better approximation? We begin by drawing the line tangent
to the graph of f (x) at the point (xo, f (xo) ). The slope of this line is f' (xo).
We label the point of intersection of the tangent line and the x-axis (Xl, 0)
and take Xl as our next approximation. Since f is differentiable, its graph
is reasonably straight on small intervals and if Xo is close to the root, then
Xl is a better approximation than Xo. This is illustrated in Figure 12.l.
J(xo)
Xl = Xo - J/(xo)·
Now if we let Xl be our estimate and apply the process again to find another
estimate X2, then we see that
Xl = Nf(xo)
X2 = Nf(XI) = Nf(Nf(xo)) = NJ(xo)
Xn = N7(xo).
To simplify notation when the function in question is clear, we will suppress
the subscript and write Nf(x) as N(x). We call Nf(x) Newton's function
for J.
EXAMPLE 12.1.
Let p(x) = x3 - x. The roots of pare -1, 0, and 1. Let's see how fast
we can find them using Newton's method. We begin by finding N(x).
N( ) _ _ p(x) _ _ x 3 - X _ 2x 3
X - X p/ (x) - X 3x2 _ 1 - 3x2 _ 1
Xo = .25
Xl = -.0385
X2 = .000114
X3 = -.00000000000299
118 12. Newton's Method
The sequence converges quickly to the nearest root. If our initial value is
.75 we get
Xo = .75
Xl = 1.227
X2 = 1.051
X3 = 1.003
Xo = .45
Xl = -.464
X2 = .567
X3 = -10.156
X4 = -6.793
X5 = -4.562
X6 = -3.091
X7 = -2.135
Xs = -1.536
Xg = -1.192
XlO = -1.038
Xn = -1.002
0.5
-1
The problem we wish to avoid with the first conclusion in Theorem 12.2
is illustrated by the function, f(x) = x 2 - 2x + 1. The only root of f is
1 and f'(l) = O. As Nf(x) =x - J/~~), Nf(l) is not defined since there
is a zero in the denominator. More generally, if r is a root of p( x) and
120 12. Newton's Method
p(x)
Np(x) = x - pl(X)
(x - r)nq(x)
= x-
(x - r)nql(x) + n(x - r)n-lq(x)
-;----:----:7-:---"-..,.:...:.~-:---:-....,......,._
(x - r)nq(x)
= x - -;----:---:;-'-;-~---'-,....::-:.~-____;~
(x - r)n-l((x - r)ql(x) + nq(x))
(x - r)q(x)
= x - -----'--,--....,.....:,-=-:......:..--,---,-
(x - r)ql(x) + nq(x)
Thus, after cancelling the term (x - r)n-l we find that Np(r) is defined
and
N(r)=r-
p
(r-r)q(r)
(r - r)q'(r) + nq(r)
=r-
0+ nq(r)
=r °
°
since q(r) =1= by assumption. Therefore, Np(x) is defined at the root and
has a fixed point there.
On the other hand, if Np(r) = r, then
per)
r - p'er) = r
and INp(r) I =
I Ip(r)p"(r) I
((p' (r))2 = °< 1.
12.1. Newton's Method for Quadratic Functions 121
and
Nor = (2ax + b)2 + A = 4a 2x 2 + 4abx + b2 + A
q 2(2ax+b) 2(2ax+b)
4a x + 4abx + b + b - 4ac
2 2 2 2 2a 2x 2 + 2abx + b2 - 2ac
2(2ax+b) 2ax+b
122 12. Newton's Method
EXAMPLE 12.5.
Let p(x) = x 2 - 1. From Theorem 12.2 we know that -1 and 1 are the
only fixed points of N(x). Since p'(O) = 0, N(O) is undefined. We will show
that WS(1) = (0, (0) and W S( -1) = (-00,0). Graphical analysis of N(x)
justifies this claim though of course it doesn't prove it. The graph of N (x)
with graphical analysis done for two initial points is shown in Figure 12.4.
We begin by proving the claim for WS(1) . We note that if Xo is in
x2 + 1
(0 , 00) , then N(xo) = _0_ _ is in [1,(0) since this is equivalent to saying
2xo
that x6 - 2xo + 1 2:: O. It remains to show that lim Nn(xd
n-+oo
= 1 whenever
12.1. Newton's Method for Quadratic Functions 123
, q(x)q"(x) X2 - 1 . ..
From equation (12.1), N (x) = (q'(X))2 =~. DIfferentIatmg
1
again we find that N"(x) = 3' Since N"(x) > 0 when x is in the interval
x
(1, (0) we know that N' (x) is increasing on this interval. As N' (1) = 0 and
x 2 -1 1
lim N'(x) = x->oo
lim - - 2 -
x->oo 2x 2'
this implies that 0 < N'(x) < ~ for all x in (1, (0). The Mean Value
Theorem implies that there is CI between 1 and Xl so that
EXAMPLE 12.7.
x 2 -1
Let p(x) = x 2 + 1. Then N(x) = - - . To get some idea of what
2x
might happen for this function we look at graphical analysis for N(x) with
initial point .1 in Figure 12.5.
Obviously this isn't as simple as the cases we have examined already.
Looking at Figure 12.5 we are led to wonder if the action of this function
might, in fact, be chaotic. We prove that this is the case by showing that
N is related to the doubling map of the circle.
Recall from exercise 9.12 that the circle is denoted 8 1 and is usually
depicted as the unit circle. Angles or points on the circle are measured
in radians in the counterclockwise direction with the initial point on the
positive x-axis. The doubling map, D : 8 1 ~ 81, is defined by D(()) = 2().
We demonstrated in exercise 9.12 that D was chaotic on 8 1 .
We define a map from the circle with the point 0 deleted to the real
numbers by ¢(x) = cot(~). This map satisfies No ¢ = ¢ 0 D, since by the
half angle identities,
L
/
/
./
./
/ -'"
3 -2 2 3
./ V
./
/'
~V
-~
[7
-3
FIGURE 12.5. Graphical analysis for 100 iterations of
Newton's function for p(x) = x 2 + 1 with .1 used as a
starting point.
Proposition 12.10 coupled with Proposition 12.9 implies that either the
dynamics of Newton's method for a cubic polynomial are the same as the
dynamics of Newton's method for a polynomial in the parametrized family
fc(x) = x 3 + ex + 1, or the dynamics are the same as those of Newton's
method for a polynomial of the form ga(x) = x 3 +ax. The following propo-
sition implies that Newton's method for functions of the latter form must be
topologically conjugate to Newton's method for one of three polynomials:
or
y=x
.\~
. . . ,. :'fi"..
2
....~'. .. ... :.l;... .'
.. ."' ......;.
1.5 .' 0(' - '"!.
~·I·.;~
:;:.,...~ -:. ~~:
'"
0.5
0
- :C
'.
I :..
.
-0.5
..
~
, '.
" .1
-1
-1. 5
-2
-1.75-1.5-1.25 -1 -0.75-0.5-0.25 0
1.5
. ...
1
. .. .. • I •
0.5
o . I t • I •
-------
-0.5
-1
-1.5
-
-2~---1~.~7-5---1~.-5---1~.-2~5---~1----0~.~7~5---0~.5----0~.2~5--~0
0.025
-0.05
EXAMPLE 12.13.
2x 3 - 1
Let f(x) = x 3 -1.265x + 1. Then N(x) = 3x 2 - 26' When iterating
1. 5
N starting at 0, we get the following sequence:
and I is empty, N 2 : I _ I, and I(N 2)'(x)1 < 1 for all x in I. The first
condition ensures that I contains no fixed point of N, and Theorem 6.3,
along with the last two conditions, guarantees that N has an attracting
periodic point with prime period 2 in I. Figure 12.9 suggests that the
interval we are looking for might be [-.03, .03]. Our goal then is to show
that 1= [-.03, .03] satisfies the following three properties:
(1) N(I)nI=0
(2) N2 : I _ I
(3) I(N 2 )'(x) I < 1 for all x in I
. , f(x)f" (x)
In equatIOn (12.1) we calculated that Nf(x) = (f'(x))2 . Hence, for
_ 3 , _ (x 3 - 1.265x + 1)(6x) .
f(x) - x - 1.265x + 1 we have N (x) - (3x 2 _ 1.265)2 . DIfferen-
tiating again and simplifying we find that
0.2
0.1
-0 . 1
-0 . 2
imply IN'(x)1 < .12 when x in I and inequality (12.6) implies IN'(x)1 < 6.4
when x is in [.79, .793], we have
I(N 2 )'(x)1 = IN'(x)IIN'(N(x))1 < (.12)(6.4) = .768 < 1
as required by condition 3.
We have shown that N has an attracting periodic point with prime
period 2 in [-.03, .03]. We leave as an exercise the proof that [-.03, .03] is
contained in the stable set of this periodic point. D
Exercise Set 12
12.1 For each of the following functions use Newton's method to find as
many roots of the equation as you can. Describe the stable set of
each root. Indicate any points you find which are not in the stable
set of any root and explain why they aren't in the stable set of a
root.
a) f(x) = x2 - 3x + 2
b) g(x) = log x
c) s(x) = sec x
d) E(x) = eX
1
e) k(x) = - - 1
x
f) p( x) = X4 + x2
g) r(x) = ..;x
x
h) F(x) = ~1
x +
a
12.2 Draw the graph of function for which Newton's method converges
to infinity whenever the initial value is greater than 2 and for which
it converges to 0 whenever the initial value is less than 2. Find a
formula for such a function. Prove that your formula works; that
is, show that the stable sets of zero and infinity are as desired.
12.3 Investigate Newton's method for f(x) = x 2 • Find the unique fixed
point and prove that its stable set is R
then If(y) - pi < €Ip - yl· Explain why this property suggests the
name super attracting.
h
12.6 A N I NVESTIGATION: ConSl·der t h e sequence di = ai-l - ai were
ai - aHl
the ai are as defined in exercise 12.5. Use a computer to determine
numerically that d i converges. To what number does it converge?
Why is this so? Can you prove that your answer is correct?
12.7 Let p+(x) = x 3 + x and Po(x) = x 3 . Show the only fixed point
of Newton's method for both functions is 0 and that WS(O) is R
(This is Proposition 12.12.)
12.8 Let f(x) = x 3 + ex + 1 where e > O. Show that Nf has exactly one
fixed point and that the stable set of the fixed point contains all
Exercise Set 12 137
real numbers. Show that when e = 0 there is still only one fixed
point and its stable set contains all real numbers except O.
12.9 Let f(x) = x 3 -1.265x + 1. Show that the interval [-.03, .03) is in
the stable set of an attracting periodic point of N f .
Hint: Use the results of Example 12.13.
12.12 a) Figure 12.10 indicates that Newton's method for the function
f(x) = x 3 + ex + 1 goes through a cascade of period doubling
bifurcations as e varies from -1.25 to -1.3. Use a sequence of
graphs of N 2 for parameter values in this range to explain why
this is to be expected. Why do we use the graph of N 2 instead of
N?
b) Using the information from part (a) and our experience with
period doubling cascades in Chapter 11, predict a parameter value
and an interval on which we would expect N 2 to be chaotic. Ex-
plain why this is a good choice. Can you find numerical or graphical
evidence to support your guess? Can you prove that your guess is
correct?
c) If N 2 is chaotic on I, does it follow that N is chaotic on IUN(I)?
Explain.
12.16 Let 8* be the set points of 8 1 which are not eventually fixed at 0
by the doubling map D((}) = 2(}. That is,
x 2 -1
b) Let N(x) = -- be Newton's method for f(x) = x 2 + 1 and
2x
let
lR* = {x in lR I Nn(x) is defined for all n}.
That is, lR* is the set of points which are never mapped to 0 by N.
Show ¢ : 8* -+ lR * defined by ¢ (x) = cot ( ~) is a homeomorphism.
Hint: It may help to use the previous exercise and the discussion
in Example 12.7.
c) Show that D : 8* -+ 8* and N : lR* -+ lR* are topologically
conjugate. Conclude that N is chaotic on lR*.
d) Use part (c) to show that N is chaotic on lR.
12.17 Show that Theorem 12.15 implies that the number of correct digits
to the right of the decimal place approximately doubles with each
iteration of Newton's method.
The reader is probably familiar with the differential equation y' = ky from
calculus. Often it is introduced as a simple model for population growth
or accumulating interest. For example, if p(x) represents the population of
bacteria in a laboratory sample at time x and the population increases by
fixed percentage during each time interval, then it is reasonable to assume
that there is k > 0 so that p' = kp: the rate of change of the size of the
population, as reflected by the derivative, is proportional to the population.
The larger the population, the more individuals that are added to it during
each time interval. If we wish a more sophisticated model that takes into
account the fact that the crowding of a large population may effect growth
negatively, then we might reflect this by adding a quadratic term to get
w' = kw - aw 2 where both k and a are positive. When the population is
small, kw is greater than aw 2 , the derivative is positive, and the population
is growing. As w increases the difference between kw and aw 2 decreases
until we reach a point where aw 2 is larger than kw, the derivative is negative
and the population is shrinking.
Both of these equations are solvable by techniques we learned in calculus.
If ~~ = kp and p(O) = Po, then we can separate the variables to find that
~ dp = k dx. Now we integrate both sides and solve for p.
p
142 13. Numerical Solutions of Differential Equations
log(p) = kx +c
P = Ge kx (13.1)
where G = Po.
. J
Similarly, if we assume that w(O) 2': 0, w(O) =1= ~, separate the variables,
and mtegrate k 2 1 dw, then we can show that
w -aw
kG
w (x) - ----:c---,.- (13.2)
- aG + e- kx
EXAMPLE 13.1.
Let pi = .5p and suppose p(O) = 1. Without solving the differential
equation, estimate the value of p(5).
If there is a solution to the differential equation, then we know the slope
of the tangent line at each point on its graph. For instance, we know
that (0,1) is on the graph since p(O) = 1 and we know that the slope
of the tangent line at (0,1) is .5(p(O)) = .5. Also, since the function is
differentiable its graph is close to the tangent line for a small interval around
O. The y-coordinate of the point on the tangent line with x-coordinate 1
can be found by adding the product of the slope and the change in x to
p(O). Doing so we find that when x = 1, the y-coordinate on the tangent
line is .5(1) + 1 = 1.5. Thus, we estimate that p(l) is 1.5. We illustrate
this in Figure 13.1.
Now the slope of the tangent to the solution at (1, p(l)) is .5(p(1)) or
approximately .5(1.5) = .75. Using the same procedure as before we obtain
a new x value by adding 1 to the current x value and estimate a new y
value by adding the product of the slope and the change in x to the current
y value. That is, we add (.75)(1) = .75 to 1.5 to find that p(2) ~ 2.25. This
is also illustrated in Figure 13.1.
Iterating the procedure, we estimate the slope of the solution at (2, p(2))
as .5(2.25) = 1.125 and add (1.125)(1) = 1.125 to 2.25 to find p(3) ~ 3.375.
13. Numerical Solutions of Differential Equations 143
2.5
(2,2.25)
~~"_bY_'75XlIoP(I)
2
1.5
(1,1.5), slope of this line isp'(1) or .75
10
(5.7.59315)
(4.5.0625)
(0. 1)
o 2 3 5
We continue with this process until we have reached Pn ~ p(x). The first
few steps in the process are illustrated in Figure 13.3.
If we wish to apply the methods of the previous chapters we need to find
a function 9 so that Pn = gn (Po) . Suppose that f is a function of P alone, for
13. Numerical Solutions of Differential Equations 145
f(x,.p,) ·h
f(xflPo) ·!J
" h II
example, f(x,p) = .5p, as in Example 13.1. Then setting g(p) = f(p) ·h+p
we can write equation (13.3) as
Thus, in the case where f depends only on p, Euler's method is the itera-
tion of a function of the real numbers. If f depends on both x and p, then
Euler's method can be modeled as the iteration of a function from ]R2 to
]R2. Unfortunately, we haven't the tools necessary to completely analyze
iterated functions of more than one variable. However, we will see in Ex-
ample 13.5 that even in the latter case we can still get some information
from our understanding of iterated maps of the real line.
EXAMPLE 13.2.
Suppose, as in Example 13.1, that p' = .5p, p(O) = Po, and we wish to
approximate p(x) = Pn by using n steps of size h = ~(x). Then
Pn = (1 + .5h)Pn-l = g(Pn-l)
= g ((1 + .5h)Pn-2) = l(Pn-2)
Therefore, the estimate we get from Euler's method converges to the solu-
tion as the step size approaches 0, and the number of steps goes to infin-
ity. 0
EXAMPLE 13.3.
Consider pi = kp where k is any real number and p(O) = po. As we
have seen, the actual solution is p(x) = poe kx . Suppose we approximate
the solution by using Euler's method with fixed step size h > O. We
denote successive approximations, Po = p(O), PI ~ p(h), P2 ~ p(2h), ....
Using an analysis similar to that in the preceding example, we find that
Pn = kpn-I h + Pn-I· That is, Pn = (1 + kh)npo or Pn = gn(po) where
g(p) = (1 + kh)p.
We note that the only fixed point of 9 is 0 and that Ig' (p) I = 11 + khl
for all t. If k > 0 we see that 19'(p) I > 1 for all p and that gn(po) tends to
(Xl whenever Po i= 0 and remains constant at 0 when Po = o. Qualitatively
this is the same behavior we see in the solutions of the differential equation;
when k > 0, the solution p(x) = poe kx tends to infinity as x increases and
is identically zero when Po = O.
Now consider k < o. If 0 < h < -~, then g'(p) = 11 + khl < 1. Hence,
when h is in this range, then gn(po) tends 0 for all Po since WS(O) = lEt If
h = - ~, then 19' (p) I = 1 and gn (Po) = Po for all n. Finally, if h > - ~, then
gn(po) tends to infinity. Which behavior is correct? Checking the solution
we see that when k < 0, p(x) = poe kx tends to 0 as x goes to infinity so
Euler's approximation has the same qualitative behavior as the solution if
and only if h < - ~. How small the step size h needs to be in order to
guarantee the correct qualitative behavior depends on the absolute value
of k. No single h will work for all k. The larger k is in absolute value the
smaller h must be to exhibit the correct behavior. D
EXAMPLE 13.4.
Let pi = P - p2, p(O) = Po, and h represent a positive step size for an
Euler approximation. Again, we denote the estimated value of p( nh) by
Pn. If g(p) = (p - p2)h + p, then
In exercise 10.4 we saw that all quadratic maps are topologically conjugate
to maps of the form qc(P) = p2 + c via a linear conjugacy. In particular,
there exists T(p) = mp+b so that Toqc = gOT. A few computations reveals
that g(p) = (p - p2)h + p is conjugate to qc(P) = p2 + c where c = 1 ~ h .
2
Since h > 0, there is a positive value of h for which the qualitative behavior
of g(p) is the same as the qualitative behavior of qc(p) = p2 + c for each
c < ~. Recall from exercise 11.5 that when -3/4 < c < 1/4, the only
periodic points of q are a single attracting fixed point and a repelling fixed
point. However, when c < -3/4 we see a cascade of period doubling
bifurcations followed by chaos. Consequently, the qualitative behavior of
the Euler approximation depends on the choice of h and hence c. We now
ask how this corresponds to the qualitative behavior of the solutions of the
differential equation.
From equation (13.2) we know that the solution of p' = p - p2 has the
form p(x) = C C
+ e- X where C = 1-p(O~)'
pO
Therefore, if p(O) is not 1 or 0,
then p(x) tends to 1 as x tends to infinity. If p(O) = 1, then the constant
function p(x) = 1 is a solution which attracts all other solutions and if
p(O) = 0, then the solution is the constant function p = 0 which repels all
other solutions. Returning to the Euler approximation, we are led to ask
if the solutions p = 1 and p = 0 correspond to the attracting and repelling
fixed points we see when -3/4 < c < 1/4.
We note that I-fk is an attracting fixed point of qc and l+fk is
a repelling fixed point. In the exercises we ask the reader to show that
the topological conjugacy maps I-fk to 1 and 1+f4Cto O. So when
-3/4 < c < 1/4, then 1 is an attracting fixed point of Euler's approxi-
mation and 0 is a repelling fixed point. Hence, the qualitative behavior of
the approximation is the same as the qualitative behavior of the solution.
Since c = I-t, -3/4 < c < 1/4 implies 0 < h < 2 and we know that
Euler's approximation displays the correct behavior when h < 2. When
h> 2 the approximation still has fixed points at 0 and 1, but they are both
repelling. As h tends towards infinity, c tends towards negative infinity and
the approximation goes through a series of period doubling bifurcations.
When h > 3, c < -2 and by exercise 10.4, the Euler approximation is
chaotic on a Cantor set with the remaining orbits tending to infinity. A
few typical behaviors are illustrated in Figure 13.4.
As a concluding comment, we note that if Po = 3 and h = 1, then the
Euler approximation does not converge to 1, even though the value of h is
within the range for which 1 is an attracting fixed point. We explore this
apparent contradiction in exercise 13.2. 0
13. Numerical Solutions of Differential Equations 149
1.4
1.2
0.8
0.6
0.4
0.2
0 10 20 30 40
EXAMPLE 13.5.
Now consider the differential equation pi = - xp. By using techniques
similar to those we used earlier, we find that the solution of this equation
is p(x) = Ce- x2 / 2 where C = p(O). However, computing the Euler ap-
proximation in this case is a bit more complicated since the value of the
derivative is dependent on both the value of x and the value of p. We
extend our notation to include this fact by labeling the initial point in the
approximation (O,Po). Then the first iterate of Euler's approximation with
step size h yields (h,pd, the second yields (2h,P2), and so forth, with the
nth iterate being (nh , Pn). Note that as before we increase the first coor-
dinate by the step size h with each iterate. Also, the slope of the tangent
line at the point (nh,Pn) is -(nh)Pn. To determine Pn+l we use equation
(13.3). We have
150 13. Numerical Solutions of Differential Equations
G(x,p) = (x + h, y - xph).
This is more or less the behavior we expected, though the rise in the value
of P for the last few iterations is worrisome. Looking ahead we see things
aren't getting any better: in fact they are substantially worse.
and
Finally,
Clearly, things are getting out of hand. The iterations are swinging wildly
from one side of 0 to the other and the magnitude of the swings is growing.
Perhaps it will settle down again, but it doesn't look good. Taking a
hint from our earlier examples we try to avoid this unpleasant behavior
by using a smaller step size. In the past a small enough step size was
all that we needed to guarantee that the qualitative behavior of the Euler
approximation was the same as that of the actual solution. Letting h = .1
13. Numerical Solutions of Differential Equations 151
PlOO ~ 10- 43
P200 ~ _10- 99
P300 ~ _10- 83
P500 ~ _1011
PlOOO ~ _10 412 .
Exercise Set 13
13.1 a) Complete the steps in Example 13.1 and verify the estimate for
p(5).
b) Let p' = .5p and p(O) = 1. Using Euler's method with 5 steps
each of length .2, estimate p(l) without using a computer.
13.3 For each of the following differential equations, represent the Euler
approximation as the iteration of a function of the real numbers.
Determine for which step sizes the qualitative behavior of the ap-
proximation is the same as the qualitative behavior of the solutions.
a) p' = -!p
b) p' = p - 2p2
c) p' = _p2
Exercise Set 13 153
Note: For part (c), if p(O) f:- 0, then the solution is p(x) = x.!.c
where c = ptO) and if p(O) = 0, then the solution is the constant
function p = O.
13.4 In Example 13.3 we claimed that when p' = kp and k < 0, the
larger k is in absolute value, the smaller the step size must be in
order for the Euler approximation to exhibit the correct qualitative
behavior. Use derivatives to explain why this is so. For which
values of k must we use a step size which is smaller than .01 in
order to guarantee the correct qualitative behavior?
13.5 In Example 13.5 we claim that if x > 2/h, then the growth of the
second coordinate of the function G(x,p) = (x+h, (l-xh)p) under
iteration of G behaves qualitatively like iteration of the function
£(p) = -mp where m > 1. In particular, as n tends to infinity, the
absolute value of the second coordinate of Gn (x, p) tends to infinity
as well. To make this precise, fix h > 0, let (xo,Po) be an initial
point satisfying xo > 2/h, and define £(p) = (l-xoh)p. Prove that
the absolute value of the second coordinate of Gn (xo, Po) is greater
than or equal to l£n(po)1 for all natural numbers n.
Hint: Let g2(X,P) = (1 - xh)p and prove that
Ig2(G n (XO,Po))1 ~ £n+1(po) = 11 - xohln+lpo.
Many interesting and beautiful results in dynamics are seen in the realm
of complex functions. In the last two chapters of this text we examine
Newton's method for complex functions and the dynamics of the quadratic
map qc(z) = z2 + e when z is complex. Both of these subjects yield in-
teresting mathematics, as well as beautiful illustrations. In particular, the
well-known Mandelbrot set is intimately connected with the dynamics of
the quadratic map.
and (14.1 )
This angle is generally called the argument of a+bi and denoted arg(a+bi).
Of course, there are a few problems with the definition of the argument.
First, there are infinitely many angles which satisfy equation (14.1). Two of
these are demonstrated in Figure 14.1. Second, the argument isn't defined
when (a, b) = (0,0). The first of these problems is easy to deal with; we just
need to remember that if two complex numbers have the same modulus and
their arguments differ by an integer multiple of 27f, then they are the same
number. For example, V3 + i has modulus 2 and argument ~ + 2n7f where
n can be any integer. We acknowledge the second problem by not defining
an argument for zero. (There appear to be several good candidates for the
argument of zero, but all of them lead to contradictions in complex function
theory. Hence, we must resign ourselves to not being able to express zero
in polar coordinates. As we will see, this is not an undue hardship.)
Now let z be a nonzero complex number, r be the modulus of z , and a
be the argument of z. Then it should be clear that
EXAMPLE 14.1.
Find the square roots of z = 4e i i .
If we indicate the root by w, then Iwl 2 = 4 and 2 arg(w) = ~. Obviously,
Iwl = 2 and arg(w) = ~ suffice. How do we find the second root?
Recall that we may write arg(z) = ~ + 2mf where n is any integer. So
the equation we really need to solve is 2 arg( w) = ~ + 2mf. Dividing by 2
we arrive at arg( w) = ~ + mr. Writing this as a real number plus integer
powers of 21f we see that arg(w) = ~ + 2n1f or arg(w) =
'IT' 7r' 97l'" .
9;+ 2n1f. So the
roots of z = 4e'42 are WI = 2e s ' and W2 = 2eT'. The complex number z
and its roots are shown in Figure 14.2. 0
Note that this definition is the same as the definition used for func-
tions of the real numbers (Definition 2.7) and functions on metric spaces
(Definition g.4g). In particular, the numbers E and 8 are still positive real
numbers representing distances. When we visualize continuous functions of
the complex plane, we can use the same imagery as we used when consider-
ing continuous functions of the real numbers. The definition implies that if
f is continuous at zo, then for every neighborhood of f(zo), there is a 8 > 0
so that any point in the neighborhood of Zo with radius 8 is mapped inside
the neighborhood of f(zo). That is, f(N6 (zo)) is contained in N€(f(zo)).
We represent this in the complex plane by drawing the neighborhoods as
disks and visualize the neighborhood of Zo being mapped inside the neigh-
borhood of f(zo). This is illustrated in Figure 14.3.
The definition of the derivative for complex functions is also the same
as that of the derivative for real functions.
. f(z) - f(zo)
11m
Z->Zo Z - Zo
exists. In this case we say f is differentiable at Zo and and call the limit
f' (zo) or the derivative of f at Zo. A function is differentiable if it is
differentiable at each point in its domain.
160 14. The Dynamics of Complex FUnctions
Of course, the problem with this definition is that we must define what
. J(z) - J(zo) .
we mean by hm for complex valued functIOns. We recall
Z->Zo z - Zo
from calculus that lim g(x) = L if the value of g(x) is close to L when x
X--+XQ
iJ Jor every E > 0, there exists 0 > 0 so that iJ 0 < Iz - zol < 0, then
IJ(z) - LI < E.
If 0 and E are small, then the expression 0 < Iz - zol < 0 implies that
z is close to but not equal to Zo and IJ(z) - LI < E guarantees that J(z)
is close to L. Drawing a diagram in the complex plane which illustrates
Definition 14.4 can be an instructive exercise and the reader is encouraged
to do so. Notice that we allow z to approach Zo from any direction. This
is a significant difference from the real case in which we can approach
numbers only from the left and the right. As a result of this difference the
properties of differentiable complex functions are significantly different from
those of differentiable real functions. All differentiable complex functions
are analytic which means that their Taylor series expansion is defined and
converges on a neighborhood of any point in the domain. However, the
properties of complex polynomials and rational functions are analogous to
those of real polynomials and rational functions. As these are the only types
14.2. Complex Functions 161
lim I(z)
Z-+Zo
=L
if for every neighborhood N of L, there exists a neighborhood N' of Zo so
that if z is in N' and not equal to Zo, then f(z) is in N.
(£)'
9
(zo) = g(zo)I'(zo) - I(zo)g'(zo).
(g(ZO))2
THEOREM 14.7. THE CHAIN RULE. If 1 is differentiable at Zo and 9
is differentiable at 1(zo), then the composition (g 0 I) is differentiable at Zo
and (g 0 I)'(zo) = g'(f(zo))l'(zo).
_d zn = nzn-l
dz
holds for all natural numbers n.
162 14. The Dynamics of Complex Functions
EXAMPLE 14.9.
If J(z) = z2, then J'(zo) = 2zo for all Zo since
Z2_ Z2
J' (zo) = lim
Z-+Zo z - Zo
0 = lim (z
Z-+Zo
+ zo) (z-zo)
- - = 2zo.
Z - Zo
o
EXAMPLE 14.10.
Let J(z) = az where a is a complex number whose modulus is not 1.
We note that 0 is the only fixed point of J and investigate the orbit of Zo
when Zo -# o.
Clearly, P(zo) = J(azo) = a2zo and in general r(zo) = anzo. If we
write a as lale Oi and recall that an = lalne nOi , then we can write r(zo) as
lalnenOizo or in exponential form
•
/(1)
•
•
/(1)
• •
•
•
PROOF. We prove the theorem for the case 11'(p)1 < 1. The proof of
the case 11'(p)1 > 1 is left as an exercise.
Assume f(p) = p and 11'(p)1 < 1. Since ~(l-IJ'(p)l) > 0, the definitions
of the derivative and limits imply that there exists 8 > 0 such that if
0< Iz - pi < 8, then
Inequalities (14.4) and (14.5) together imply that if 0 < Iz - pi < 8, then
We note that !(1 + If'(p) I) < 1 since 1f'(P) I < 1 and set A = !(1 + I!,(p)l)·
Multiplying both sides of inequality (14.6) by Iz - pi and substituting A for
!(1 + I!,(p)l) we see that when Iz - pi < 8, then
If(z) - pi = If(z) - f(p) I :::; Alz - pI·
By induction we can show that if Iz - pi < 8 and n is any natural number,
then
Ir(z) - pi:::; Anlz - pI·
Since A < 1, this implies that r(z) tends to p as n tends to infinity
whenever z is such that Iz - pi < 8. Thus, the neighborhood of p with
radius 8 is in the stable set of p. 0
Of course, we can generalize Theorem 14.11 to periodic points in the
same way that we generalized the analogous theorem for fixed points of the
real numbers.
THEOREM 14.12. Let f be a differentiable complex function and p be
a periodic point of f with period k. If IUk),(p)1 < 1, then there exists a
neighborhood ofp which is contained in the stable set ofp. If IUk)'(p) I > 1,
then there exists a neighborhood of p such that all points must leave the
neighborhood under iteration of fk .
As with the real case, periodic points which are not hyperbolic do not
have predictable behavior in a neighborhood of the point. In the next
example we examine the behavior of points in C when we iterate a function
of the form f(z) = az where lal = 1. Note that If'(O)1 = 1 and 0 is the
only fixed point of f unless a = 1.
14.3. The Dynamics of Complex FUnctions 165
EXAMPLE 14.14.
Suppose j(z) = e9i z and Zo is a nonzero complex number. Then the
orbit of Zo under iteration of j lies on the circle which is centered at the
origin and whose radius is Izol. If 0 is a rational multiple of 'Jr, then Zo is
a periodic point of j. If the argument of 0 is not a rational multiple of 'Jr,
then Zo is not a periodic point and its orbit is dense on the circle.
To see this, note that jn(zo) = en9izo. Consequently,
EXAMPLE 14.15.
Let q( z) = z2. The only fixed points of q are 0 and 1. Since q' (0) = 0
and q' (1) = 2, 0 is an attracting fixed point and 1 is repelling. If z = refJi,
then q(z) = (refJi)(refJi) = r 2e2fJi . Iterating, we find q2(Z) = r 4e4fJi and in
general qn(z) = r2ne2nfJi. Since Iqn(z)1 = r 2n , we see that z is in WS(O) if
and only if Izl < 1 and z is in WS(oo) if and only if Izl > 1.
If Izl = 1, the dynamics are more interesting. Notice that Izl = 1 implies
that Iqn(z)1 = 1 for all n. Thus, if we denote the circle of points whose
modulus is one by 8 1 , then we can write q : 8 1 -+ 8 1 . (The notation 8 1 is
typically used to denote the unit circle in the plane.) Now if arg(z) = (),
then arg(q(z)) = 2() and we recognize our friend from exercise 9.12: the
doubling map on the circle. Each time we apply the function q(z) = z2 to a
point in 8 1 the modulus of the point remains unchanged and the argument
is doubled. We demonstrated in exercise 9.12 that the doubling map on 8 1
is chaotic so it follows that q(z) = z2 is chaotic on 8 1 as well. Orbits of
typical points in C under iteration of q are shown in Figure 14.5. D
initial pointS
_ - - - - (0.0.1)
have defined this new representation we will treat 00 as a point and include
it in the domain and range of functions.
To define the Riemann Sphere we use stereographic projection, the same
technique which is sometimes used to create maps of the world. A good way
to visualize it is to picture a sphere with a radius of one and center at the
origin in ]R3. We denote this sphere l by S2. The xy-plane represents the
complex plane and intersects the sphere at the equator. This arrangement
is pictured in Figure 14.6 (p. 167).
At the top of the sphere is the point (0,0,1). If we draw the line through
the point at the top and any point in the complex plane, then it intersects
the sphere at exactly one point other than (0,0,1). If we draw a line from
(0,0,1) to a point which is greater than one in modulus, then it intersects
the top half of the sphere. If we draw one through a point which is on the
unit circle, then it intersects at the equator, and if we draw a line from
(0,0,1) through a point which is less than one in modulus, it intersects the
sphere at a point on the bottom half of the sphere. To get a correspondence
between points on the sphere and the points in the complex plane we assign
each complex number to the point on the sphere which is on the line through
(0,0,1) and the number. It should be geometrically obvious that this is
a one-to-one, onto correspondence between all of the points on the sphere
except (0,0,1) and the complex plane. If in addition we assign (0,0,1)
to the point 00, then we have a one-to-one, onto correspondence between
S2 and C U {oo}. We leave as an exercise the verification of the fact
that (0,0,-1) corresponds to 0, (0,1,0) corresponds to i, and (~,~, 0)
corresponds to (~V3 + 2.;2)(1 + i).
Now consider a sequence of points 0, Zl, Z2, Z3, ... which satisfy IZkl = k.
We would like limk--+oo Zk = 00 since it is clear that the points march
towards (0,0,1) on S2. To define this limit we need to define a set of
neighborhoods for each point on S2. Once we have defined neighborhoods
we can also define open sets, convergence of sequences, and continuous
functions without reference to a metric.
DEFINITION 14.16. Let z' be a point on the Riemann Sphere which cor-
responds to Z in the complex plane. The set U' on the Riemann Sphere
°
is a neighborhood of z' if the corresponding set U is a neighborhood of Z
in the complex plane. Let E > and define NE (00) = {z in C I Izi > ~}.
1 We have already used the symbol 8 1 to indicate a sphere (circle) with radius one
and center at the origin in IR2. The symbol 8 2 is used here to indicate a sphere with
radius one and center at the origin in IR3. In general, 8 n indicates a sphere with radius
one and center at the origin in IRn+1. Note that position can be defined in terms of one
variable (the angle) on 8 1 , it takes two variables to describe position on 8 2 , and it takes
n variables to describe a position on 8 n , hence the superscript.
14.4. The Riemann Sphere 169
DEFINITION 14.17. A set is open if every element of the set has a neigh-
borhood that is completely contained in the set.
EXAMPLE 14.21.
Let J(z) = Z2 - 1. A quick calculation reveals that N(z) = z22: 1.
From Theorem 14.20 we know that 1 and -1 are attracting fixed points of
N. We claim that if a+bi is an element ofC and a> 0, then a+bi is in the
stable set of 1. On the other hand, if a + bi is such that a < 0, then it is in
the stable set of -1. Finally, we claim that N is chaotic on the imaginary
axis. We prove this claim by showing that N(z) is topologically conjugate
to q(z) = z2 on the extended plane via a Mobius transformation.
If such a Mobius transformation exists, we should be able to find it by
searching for one that has the required properties. Since any three points
determine a circle and Mobius transformations must map circles to circles,
we look for three significant points of N and the corresponding points of
q to which they must be mapped. The attracting fixed points of N are 1
and -1, so we map them to 0 and 00, the attracting fixed points of q. The
only other fixed point of N is 00 and so it must be mapped to the repelling
fixed point 1 of q. If we are lucky, the coefficient of z in the numerator
of the Mobius transformation we are looking for is not zero, and we can
divide the numerator and denominator of the transformation by it. In this
case, the requisite transformation has the form
T(z) = z + b .
cz+d
Setting T(l) = 0, T( -1) = 00, and T(oo) = 1 and solving for b, c, and d
we find that
T(z) = z -1.
z+l
Verifying that ToN = q 0 T is an easy exercise in algebra. As T is a
homeomorphism, q and N are topologically conjugate.
Since T(O) = -1, T(i) = i, and T(oo) = 1, the image of the imaginary
axis must be the unit circle. (The unit circle is the only circle passing
through -1, i, and 1, and T maps circles and lines to circles and lines.)
Since we know from Example 14.15 that q is chaotic on the unit circle
and that T is a topological conjugacy, it is tempting to conclude that N
is chaotic on the imaginary axis. This is in fact true, but we have to be
a little bit careful in proving it. Theorem 10.3 states that if a function of
a metric space is topologically conjugate to a second function of a metric
space which is chaotic, then the first function is chaotic as well. But T is a
topological conjugacy between functions on the Riemann Sphere, which we
have not defined as a metric space. As the definition of sensitive dependence
on initial conditions uses a metric, the notion of chaotic as we have defined
it doesn't even make sense for functions of the Riemann Sphere, unless we
define a metric there.
172 14. The Dynamics of Complex Functions
We have two alternatives for dealing with this problem. First, we can
define the distance between two points on the Riemann Sphere by
d[(Xl' X2, X3), (Yl, Y2, Y3)] = J (Yl - xd 2 + (Y2 - X2)2 + (Y3 - X3)2.
This is the usual euclidean metric used in m. 3 . It is easy to show that the
subset U of 8 2 is open (as defined in Definitions 14.16 and 14.17) if and
only if for each x in U there exists a positive real number E so that if Y is
in 8 2 and d[x, y] < E, then Y is in U. Thus, we can use this metric to define
chaotic behavior on the Riemann Sphere and Theorem 10.3 holds.
On the other hand, we can avoid these problems by noting that the
periodic points of q are dense on the unit circle and that q is topologically
transitive on the unit circle. These properties are defined only in terms of
open sets and are preserved by topological conjugacy. Thus, the periodic
points of N are dense on the imaginary axis and N is topologically transitive
there as well. Now, since sets that are open in the complex plane are also
open in the extended complex plane, we may conclude that the periodic
points of N are dense and N is topologically transitive on the imaginary
axis, considered as a subset of the complex plane. It follows from Theorem
9.20 that N is chaotic on the imaginary axis.
It remains to show that a + bi is in the stable set of 1 if and only if a > 0
and it is in the stable set of -1 if an only if a < o. But, T(l) = 0 and
T({a+bi I a> O}) = {z Ilzl < I},
which is the stable set of 0 under iteration of q. Since T preserves stable
sets the first assertion holds. A similar observation demonstrates the second
statement. The reader may wish to verify that the proof that topological
conjugacies preserve stable sets does not depend on the existence of a metric
in the spaces in question. 0
When reflecting on the previous example, we note that the imaginary
axis is the perpendicular bisector of the line segment between the roots of
f(z) = Z2 - 1. Also, recall from our work in Chapter 12 that Newton's
method had similar dynamics for large classes of quadratic polynomials.
We are led to ask for which complex quadratic polynomials the perpen-
dicular bisector of the line segment between the roots divides the stable
sets of the roots, and in which of those is Newton's method chaotic on the
bisector. In 1879, the English mathematician Arthur Cayley proved that
the perpendicular bisector of the line segment joining the roots divides the
stable sets of Newton's function for all complex quadratic polynomials with
distinct roots. We can now add that it is chaotic on the bisector as well.
THEOREM 14.22. [CAYLEY, 1879] If the complex quadratic polynomial
q(x) has two distinct roots, then Nq(x) is chaotic on the perpendicular
bisector of the line segment joining the two roots. Further, the stable set
14.5. Newton's Method in the Complex Plane 173
oj each root under iteration oj Nq(x) is the set oj points which lies on the
same side oj the bisector as the root.
Theorem 14.23 is easily proven by showing that it holds for q(z) = Z2 and
then proving that Newton's function for a complex quadratic polynomial
with a single repeated root is topologically conjugate to N q • Theorem
14.23 also generalizes to complex polynomials of higher degree with a single
repeated root.
174 14. The Dynamics of Complex Functions
EXAMPLE 14.24.
Let J(z) = z3 -1. Then N(z) = 2z 3 ~ 1 and has attracting fixed points
h· k·
3z .
at 1, ea ', and ea'. To get some idea of what we mIght expect to find we
conduct a computer experiment. We select a grid of points covering the
square region whose corners are at 2 + 2i, 2 - 2i, -2 - 2i, and -2 + 2i.
We calculate the value of the 100th iterate of each of these points under
iteration of N. If the distance from the 100th iterate to 1 is less than ~,
then we assume the point is in the stable set of 1 and color it light gray. If
the distance from the 100th iterate to e 23" i is less than ~, then we assume
the point is in the stable set of e 2; i and color it dark gray. Finally, if the
distance from the 100th iterate to e 4; i is less than ~, then we assume the
point is in the stable set of e 4; i and color it black. All points which are not
within ~ of one of the roots after 100 iterations of N are left white. The
result of this experiment using a 250 by 250 grid is shown in Figure 14.8.
We can see that the dynamics of Newton's method for this cubic are
much more complicated than those of Newton's method for quadratic func-
tions. In fact, the dynamics of Newton's method for most polynomials of
degree higher than two are not completely understood. The reader is en-
couraged to use experiments similar to the one above to try to understand
the dynamics of Newton's method for cubics. In particular, we would like
to know what happens on the boundary of the stable sets. Specific sugges-
tions for projects can be found in the exercises. 0
14.5. Newton 's Method in the Complex Plane 175
Exercise Set 14
14.2 Find the modulus and argument of each of the following complex
numbers:
a) 1- 4i
b) ! + £i
c) 2
14.4 Find all of the solutions of the following equations and plot them
on the complex plane. What do you notice about the geometry of
the solutions?
a) z2 = 8.
b)z3-8=0.
c) zn - 1 = 0, for n = 2, 3, 4, 5, and 6.
14.6 Verify that the function d[a + bi, c + di] = y'(c - a)2 + (d - b)2 is
a metric on Co
Exercise Set 14 177
14.10 a) Complete the details of the induction argument used in the proof
of the first half of Theorem 14.1l.
b) Let J be a differentiable complex function and p be a fixed point
of J. Prove that if If'(p) I > 1, then there is a neighborhood of p
all of whose points must leave the neighborhood under iteration of
J.
14.11 Let g(z) = az + b where a and b are complex numbers.
a) Find the fixed point of g if it exists. In what case(s) will there
not be a fixed point?
b) Show that if a of. 1, then g(z) = az+b is topologically conjugate
to a function of the form J(z) = cz.
c) In those cases where g has a fixed point, describe its stable set.
d) In those cases where g does not have a fixed point, describe the
dynamics of g.
14.15 Let q(z) = Z2. Show that 1 is a repelling fixed point of q and that
W8(1) is dense on the unit circle in Co Why doesn't this contradict
Theorem 14.11?
14.18 Show that Definition 14.18 and Definition 9.4c are equivalent. To
do this, let Xl, X2, X3, . •• be a sequence of elements in a metric
space, then show that for every E > 0, there exists an integer
N such that d[x, Xk] < E whenever k ~ N if and only if every
neighborhood of X contains an element of the sequence.
14.19 Show that Definition 14.19 and Definition 14.2 are equivalent.
14.24 Prove Cayley's Theorem (Theorem 14.22) which states that if the
complex quadratic polynomial q(x) has two distinct roots, then
N q (x) is chaotic on the perpendicular bisector of the line segment
joining the two roots. Further, the stable set of each root under
iteration of Nq(x) is the set of points which lies on the same side
of the bisector as the root.
Hint: Use one of the suggestions given after the statement of
the theorem in the text.
Iq2(W)1 = Iq(q(w))1
= I(q(W))2 + el
;::: Iq(w)12 - lei
;::: (1c1 2+ lei + 1)2 - lei
= lel 4 + 21el 3 + 31el 2+ lei + 1
;::: 31el 2+ lei + 1.
Continuing by induction we show that
Since the right-hand side of this inequality tends to infinity as n gets larger,
our assertion that w is in the stable set of infinity whenever Iwl > lei + 1
must hold. It follows that either the orbit of a complex number under
iteration of q(z) = z2 + e is bounded by lei + 1 or the number is in the
stable set of infinity. 0
The set of points with bounded orbits is interesting enough to have
earned its own name.
DEFINITION 15.2. The set of points whose orbits are bounded under it-
eration of q(z) = Z2 + e is called the filled Julia set of q and denoted by Kc.
The boundary of a filled Julia set is called the Julia set. The boundary of
a set is the collection of points for which every neighborhood eontains an
element of the set as well an element which is not in the set.
15. The Quadratic Family and the Mandelbrot Set 183
K ..719 •.251
THEOREM 15.5. If 0 is in the filled Julia set of qc(z), then the filled
Julia set is connected. That is, given any two points in the filled Julia set,
there is a path contained in the set which connects those two points. On the
other hand, if 0 is not in the filled Julia set of qc, then the set is a Cantor
set.
Proposition 15.1 and Theorem 15.5 naturally divide the complex num-
bers into two sets: those parameter values c for which the filled Julia set is
connected and the orbit of 0 is bounded under iteration of qc (z) = z2 + c;
and those for which the filled Julia set is totally disconnected and 0 is in
the stable set of infinity. The set we investigate here is the former. More
specifically, we are interested in
M = {c I q~(O) is bounded}.
While we will not prove Theorem 15.5, in exercises 15.4 and 15.5 the
reader is invited to prove that the filled Julia set is a Cantor set and that
q~(O) tends to 00 when Icl > 2. This implies that the set M is contained
inside a disk of radius 2. In exercise 15.3 we show that the set of param-
eter values for which qc has an attracting fixed point is the interior of the
cardioid parametrized by p(O) = ~eei - ~e2ei. By Theorem 15.4, it follows
that this cardioid is contained in M since 0 must be in the stable set of an
attracting fixed point. In fact, if we define
In exercise 15.3 we also show that M2 is the disk with radius ~ and center
at -1. It is interesting to try to find the other sets Mk.
So far we have stated that M is contained inside the circle of radius
2, M contains the interior of the cardioid p(O) = ~el1i - ~e2I1i, and M
contains the interior of the disk with center at -1 and radius ~. To get a
better picture of M we conduct another computer experiment. We know
that M consists of the parameter values for which the orbit of 0 is bounded
and that if Iq~(O)1 > 2 for any n, then c is not in M. Thus, we cover the
disk of radius 2 by a grid, use each point in the grid as a parameter c, and
calculate q~OO(O). If Iq~OO(O)1 ::::: 2, we assume that c is in M and color the
corresponding grid point black. If Iq~OO(O)1 > 2, we leave the grid point
white. The results of such a program are shown in Figure 15.2.
The reader may recognize the sketch of M in Figure 15.2 as the famous
Mandelbrot set, the first images of which were published by Benoit Man-
delbrot in 1980. The Mandelbrot set has been studied intensely since its
discovery and much is known about it, though there is still more to be
learned. It is known that it is connected, and the manner in which the
various bulbs are connected is also understood. The details of this connec-
tion are well beyond the scope of this book. An interesting and accessible
discussion of the properties of the Mandelbrot set can be found in Chapter
14 of the book Chaos and Fractals by Peitgen, Jurgens, and Saupe. It is
the author's understanding that a more complete account of the properties
of the Mandelbrot set will be included in the forthcoming book on complex
dynamics by J. Milnor. Both books are listed in the references.
2r---~----------~----~--~----~--------~
1.5
0.5
-0.5
-1
-1.5
-2~--~----~----~----~--~----~----~--~
-1.5 -1 -0.5 0 0.5 1 1.5 2
Exercise Set 15
15.1 Show that all complex quadratic polynomials are topologically con-
jugate to a polynomial of the form qc(z) = Z2 + c.
Hint: See exercise 10.4.
15.3 a) Show that the set of parameter values c for which qc(z) = Z2 +c
has an attracting fixed point is the interior of the cardioid
1 t 1 2n"
p(B) = _eofl"
- -e ot.
2 4
Hint: We wish to find c so that there exists Zo = re lii such that
z5 + c = Zo and Iq~ (zo) I < 1. Substitute for Zo and solve for c.
b) Show that the set of parameter values for which qc(z) has an
attracting periodic point with prime period two is the interior of
i
the disk with radius and center at -1.
**c) Describe the set of parameter values for which qc(z) = Z2 + c
has an attracting periodic point with prime period k for as many
188 15. The Quadratic Family and the Mandelbrot Set
15.5 LOOKING FOR CANTOR DUST: It has been suggested that Cantor
sets in the complex plane be called Cantor dust since, like a cloud
of dust, no two points are connected yet there is another point
close to every point in the cloud. In this exercise we show that the
filled Julia set of the function qc(z) = Z2 + c forms such a Cantor
cloud when lei> 2.
Let Icl > 2 and define q(z) = Z2 + c. Let D be the circle with
radius Icl and center at the origin.
a) Use exercise 15.4a to show that the filled Julia set is contained
inside D. In fact, Kc lies inside q-n(D) for all n.
b) Show that q-l(D) is the figure eight with crossing at the origin
shown in Figure 15.3.
c) Prove that q-2(D) must be inside q-l(D) and that q-2(D)
consists of two figure eights as demonstrated in Figure 15.3.
d) Prove that q-n(D) consists of 2n figure eights situated as in
Figure 15.3. Further, show that as n goes to infinity the diameter
of the figure eights in q-n(D) converges to O.
Exercise Set 15 189
Below we define f(x) = x 3 and then iterate f three times starting at 1.2.
Clear[f] ;
f [x_] ;= x-3;
192 A. Computer Algorithms
StartingValue = .1;
Firstlteration = 10;
Lastlteration = 20;
i=O;
Y = N[StartingValue];
While[i <= Lastlteration,
If[i >= Firstlteration,Print[i," ",N[y,8]]];
y = h[y] ;
i = i+l]
The output of this program is as follows:
10 0.86414351
11 0.41089828
12 0.84721309
13 0.45305075
14 0.86728519
15 0.40285557
16 0.84197036
17 0.46569696
18 0.87088155
19 0.39356405
20 0.83534986
A.I. Iterating Functions 193
StartingValue = .1;
Firstlteration = 10;
Lastlteration = 20;
SigDigits = 64;
i=O;
y = SetPrecision[StartingValue,SigDigits];
While[i <= Lastlteration,
If [i >= FirstIteration, Print [i, " " ,N [y, 8]]] ;
y = SetPrecision[h[y] ,SigDigits];
i = i +1]
As the function h(x) = 3.5x(1 - x) is not sensitive to small changes in
initial conditions, the output of this program is the same as that of the
preceding program.
for x in [0, ~]
T(x) = {22x,
2x, for x in [~, 1].
xmin O',
xmax 1;
NumberOflterations 2',
194 A. Computer Algorithms
Clear[h];
h[x_] := 3.5x (1-x)
StartingValue = .1;
FirstIt = 0;
LastIt = 20;
xmin = 0;
xmax = 1;
i 0;
y = N[StartingValue] ;
DataTable = {{y,y},{y,h[y]}};
AppendTo[DataTable,{h[y] ,h[y]}] ;
a) b)
First1t = 100
LastIt = 200
Seed = .5
rmin = 2
rmax = 4
Norsteps 250
xmin = 0
xmax = 1
realC = 0
imaginC = 1
modsq = (realC)-2 + (imaginC)-2
xmin = -2
xmax = 2
ymin = -2
ymax = 2
xgrid 250
ygrid = 250
iterations = 25;
There are other methods for computing images of Julia sets. One utilizes
the fact that an equivalent definition of a Julia set for the parameter c is
the closure of the set of repelling periodic points of qc (z) = z2 + c. Another
A.S. The Mandelbrot Set 199
uses the fact that for any positive €, we know the circle with radius 1 + Icl +€
is not in the filled Julia set Kc. Successive inverse images of this circle are
computed. As these inverse images must converge to Kc as the number of
iterations goes to infinity, we can often get a good picture quickly. Both of
these methods are discussed in the book by Peitgen, Jurgens, and Saupe,
which is list ed in the references.
xmin = -2;
xmax 2;
ymin -2;
ymax 2;
xgridpoints = 400;
ygridpoints = 400;
iterations = 75;
mandelbrot = {};
y y + ystep
c = (x + xstep/2) + (y + ystep/2) I;
r = Abs [c]
w = c;
For [i=O, i<=iterations, i++,
If[i == iterations,
AppendTo[mandelbrot,
Rectangle [{x,y},{x+xstep,y+ystep}]]] ;
If[Re[w]-2 + Im[w]-2 > 2,Break[]];
w = q[w,c]
]]]
xmin -2;
xmax 2;
ymin -2;
ymax 2;
xgridpoints = 300;
ygridpoints = 300;
iterations = 50;
stable1 {GrayLevel[.8]};
stable2 {GrayLevel [.4] };
stable3 {GrayLevel[O]};
A.6. Stable Sets of Newton's Method 201
Show[Graphics[{stable1,stable2,stable3}],AspectRatio->1]
The boundaries ofthe stable sets of the fixed points of Newton's function
for a cubic polynomial are the closure of the set of repelling periodic points.
By exploiting this fact, one can write algorithms which will compute these
boundaries efficiently. We again refer the interested programmer to the
book by Peitgen, Jurgens, and Saupe for suggestions as to how this might
be done.
References
Topics in Mathematics.
BAK, J. AND NEWMAN, D. J. Complex Analysis, Springer-Verlag, 1982
BELDING, D. AND MITCHELL, K. Foundations of Analysis, Prentice-Hall,
1991
CHURCHILL, R. V. AND BROWN, J. W. Complex Variables and Applica-
tions, 5th edition, McGraw-Hill, 1990
OSTROWSKII, A. M. Solutions of Equations and Systems of Equations, 2nd
edition, Academic Press, 1966
RUDIN, W. Principles of Mathematical Analysis, 3rd edition, McGraw-Hill,
1976
WILLARD, S. General Topology, Addison-Wesley, 1970
206 References
There are a multitude of other options one can pursue. Probably the
easiest choice for those with some computer experience is to write their
own programs. None of the programs for doing the dynamics experiments
mentioned in this text need be more than 20 or 30 lines long. Many will
have less than 10. To write the programs one needs to be able to input
initial data to a program, construct loops to do the iteration, and output
graphics information to the screen; that's it. Several of the books listed
above in the section on dynamics include computer programs. True Basic
programs which can easily be converted into other languages are included
in appendix B of the text by Devaney. Numerous programs are included
in the book edited by Peitgen and Saupe. Efficient programs in Basic for
computing Mandelbrot and Julia sets are found in the book by Peitgen,
Jurgens, and Saupe.
For those who want a canned package the alternatives are almost endless.
There are probably hundreds of public domain packages which will do the
job. Check local bulletin boards or use Archie and search the internet. By
the time this book appears in print, 1 expect to have located or written
programs for the NeXTstep environment which will do the computations
necessary for the exercises in my text. The software will be compatible with
both the original NeXT hardware and the intel version. Those interested
in getting a copy (free or close to it) may contact me at the Mathematics
Department, Allegheny College, Meadville PA 16335. My internet address
is: rholmgre@[Link]. I accept NeXT mail. Anyone who has already
written such programs for the NeXT or knows of such programs is especially
encouraged to contact me. (I know that there is a program to compute the
Mandelbrot set in the NextDeveloperjDemos folder, thanks. If you have a
NeXT and haven't seen it, have a look; it's a nice program.)
Below 1 have included a reference for a workbook and software which
is intended to serve as a supplement to Devaney's text by the same name
(listed above). While I have never seen it, I believe it contains programs
capable of completing all the computations necessary for completing the
exercises in my text. It requires a MacIntosh computer running system
6.0.5 or higher. To get adequate speed one should have a Mac II and
a math coprocessor. For the exact system requirements, please contact
Addison-Wesley.
GEORGES, J. JOHNSON, D. AND DEVANEY, R. L. A First Course in
Chaotic Dynamical Systems: Laboratory Addison-Wesley, 1992
Index
A subset A of a metric space B is dense in B if every point in B is an accumulation point of A, a point of A, or both . This can also be characterized by the statement that for each point x in B and for every ε > 0, there exists a point y in A such that the distance between x and y is less than ε . Additionally, for every point in B, there exists a sequence of points in A that converges to that point .
Open sets in a metric space are characterized by the property that each point within the set has a neighborhood completely contained in the set . Closed sets contain all their accumulation points . A set is dense in a space if every point in the space is either an accumulation point of the set or a point of the set, implying every open set contains points from the dense subset .
In the complex plane, continuity and differentiability are extended from real functions with several notable differences. Continuity for complex functions involves a function mapping neighborhoods around points effectively into neighborhoods in essence is similar to real functions . For differentiability, one key difference is that in the complex plane, differentiability at a point requires the function to be analytic, meaning it is infinitely differentiable, and its Taylor series converges to the function in a neighborhood . This arises because limits can approach a point from any direction, unlike real functions, confined to left or right direction . Visualizing these concepts in the complex plane often involves mapping through the Riemann Sphere, which allows for an extended view, incorporating the point at infinity and making continuity and differentiability more comprehensible. The Riemann Sphere provides a means to visually interpret mappings, where complex numbers are projected onto a sphere, enabling visualization of transformations or dynamics like Mobius transformations, which maintain properties such as mapping circles to circles . This visualization aids in understanding how complex functions maintain certain geometric properties, unlike functions limited to the real plane .
Qualitative analysis focuses on the behavior and characteristics of solutions, such as stability, periodicity, or tendencies as variables approach infinity. For example, knowing a solution tends to zero gives insights into system behavior regardless of precise values . Quantitative analysis involves exact or approximate solutions at specific points, which can suffer from computational inaccuracies in Euler's method beyond specific step size limits, affecting reliability without always reflecting qualitative behavior accurately . Euler approximations aim for quantitative accuracy, but qualitative traits like convergence to fixed points remain vitally informative .
In a metric space, a point sequence \( x_1, x_2, x_3, ... \) converges to a point \( x \) if, for every \( \epsilon > 0 \), there exists an integer \( N \) such that for all \( k \geq N \), the distance \( d(x, x_k) < \epsilon \). Convergence relates closely to accumulation points; a point \( x \) is an accumulation point of a subset \( S \) if every neighborhood of \( x \) contains another point from \( S \), distinct from \( x \) itself . Accumulation points generalize the idea of limits, providing a foundation for sequences that may not necessarily converge but still "accumulate" around a point . This interplay illustrates how topological concepts in metric spaces, such as open and closed sets, revolve around these distance-based definitions, forming a basis for further properties like continuity and density .
The Euler method is a first-order numerical technique used to approximate solutions to ordinary differential equations (ODEs). It operates by using the derivative to estimate the slope of the solution over small intervals, known as step sizes, and iteratively computes future values based on this estimation. Specifically, if the derivative at a point \( (x_o, p(x_o)) \) is denoted as \( f(x_o, p(x_o)) \), then the next value \( p(x_1) \) is approximated as \( f(x_0, p(x_0)) \cdot h + p(x_0) \), where \( h \) is the step size . The method continues by recalculating at each subsequent point \( x_k = x_{k-1} + h \). Challenges in implementing the Euler method include ensuring stability and accuracy. The method requires a sufficiently small step size to accurately capture the behavior of the solution, particularly for stiff differential equations where larger step sizes can lead to divergence from the true solution . Furthermore, even with appropriately small step sizes, the Euler method may not always exhibit the correct qualitative behavior of the differential equation's solution, especially for complex or rapidly changing equations . Additionally, using very small step sizes can lead to computational inefficiency and numerical errors due to limitations in machine precision .
Bifurcation in dynamical systems refers to a qualitative change in the behavior of a system as a parameter is varied. A saddle-node bifurcation occurs when two fixed points, one stable and one unstable, collide and annihilate each other; this is seen when parameter values make the number of equilibria change . A pitchfork bifurcation involves the splitting of a stable fixed point into one unstable point and two stable points or vice versa, usually occurring symmetrically and resulting in a change from one stable state to multiple states . A trans-critical bifurcation is characterized by the exchange of stability between two fixed points as they intersect; one fixed point becomes stable while the other becomes unstable . These bifurcations represent the critical points where the long-term behavior of a dynamical system changes drastically.
A bifurcation diagram provides a visual representation of how fixed points and periodic points of parameterized functions change as parameters vary. It plots periodic values against parameter values, distinguishing stable (attracting) points with solid lines and unstable (repelling) points with dotted lines . This visual approach helps explain phenomena such as the birth or annihilation of periodic points and illustrates complex changes like pitchfork bifurcations .
Chaotic behavior in complex functions can arise through dynamics such as those seen in Newton's method for complex polynomials where sensitivity to initial conditions and dense winding paths in the complex plane lead to unpredictability. This can be mathematically characterized by the existence of unstable periodic orbits and fractal basin boundaries, often illustrated through topological conjugacies and bifurcations . Mobius transformations that map circles and lines to other circles demonstrate structure amid chaos, while the presence of chaotic sets like Julia sets represents chaotic behavior .
Möbius transformations, defined as T(z) = (az + b) / (cz + d) with ad - bc ≠ 0, are homeomorphisms of the extended complex plane. They preserve angles and map circles to circles, which includes lines as a special case of circle through infinity . They play a crucial role in complex dynamics by assisting in the study of the Riemann sphere and facilitating transformations that simplify the analysis of function behaviors, such as mapping the real axis to a unit circle .