# Elementary differential equations 7th edition - Boyce W.E

ISBN 0-471-31999-6

**Download**(direct link)

**:**

**476**> 477 478 479 480 481 482 .. 486 >> Next

250

Exploration 13.1

3. In the Discrete Tool enter the tent map Tc on the interval 0 < x < 1:

[2cx, 0 < x< 0.5

Tc (x) = c(1 — 2abs(x — 0.5)) = \

|2c(1 — x), 0.5 < x< 1

where the parameter c is allowed to range from 0 to 1. Describe and explain what you see as c is incremented from 0 to 1. [Suggestion: use the Edit option in the Menu box for the bifurcation diagram and set 200 < n < 300 in order to suppress the initial transients.] Any orbits of period 2? Period 3?

Answer questions in the space provided, or on

attached sheets with carefully labeled graphs. A

notepad report using the Architect is OK, too.

Name/Date______

Course/Section

Exploration 13.2. Circle Maps

Another common type of discrete dynamical system is a circle map, which maps the perimeter of the unit circle onto itself. These functions arise when modeling coupled oscillators, such as pendulums or neurons. The simplest types of circle maps are rotations that take the form

Rm(Q) = (Q + m) mod 2n where 0 < Q < 2n and m is a constant.

1. Show that if m = (p/ q)n with p and q positive integers and p/q in lowest terms, then every point has period q.

2. Show that if m = an with a an irrational number, then no point on the circle is periodic.

252

Exploration 13.2

3. What is the long-term behavior of the orbit of a point on the circle if m = an, where a is an irrational number?

Answer questions in the space provided, or on

attached sheets with carefully labeled graphs. A

notepad report using the Architect is OK, too.

Name/Date______

Course/Section

Exploration 13.3. Two-Dimensional Maps and the Discrete Tool

A two-dimensional discrete dynamical system looks like this:

Xn+1 = f(Xn, yn, c) (3)

yn+1 = g(Xn, yn, c)

where f and g are given functions and c is a "place holder” for parameters. For given values of c, xo, and y0, system (3) defines an orbit of points

( Xo, yo ), (X1, y1), (X2, y2 ),...

in the Xy-plane. The two-dimensional tab in the Discrete Tool allows you to explore discrete systems of the form of system (3).

1. Open the Discrete Tool and explore the default system (a version of what is known as the HenonMap):

Xn+1 = 1 + yn - aXn (4)

yn+1 = bXn

where a and b are parameters. For fixed values of the parameters a and b find the fixed points. Are they sinks, sources, or neither? How sensitive is the long-term behavior of an orbit to small changes in the initial point (xo, yo)? What happens if you increment a through a range of values? If you increment b? Any period-doubling sequences? In your judgment, is there any long-term chaotic wandering? [Suggestion: Keep the values of a and b within small ranges of their default values to avoid instabilities.]

254

Exploration 13.3

2. Repeat Problem 1 with the following version of the Henon map:

xn+1 = a — 4 + byn yn+1 = xn

Start with a = 1.28, b = —0.3, x0 = 0, y0 = 0.

Answer questions in the space provided, or on

attached sheets with carefully labeled graphs. A

notepad report using the Architect is OK, too.

Name/Date______

Course/Section

Exploration 13.4. Julia and Mandelbrot Sets and the Discrete Tool

Note that the color schemes for the Julia and Mandelbrot sets in Module 13 differ from those in the discrete tool.

1. Use the Discrete Tool to explore the Mandelbrot set and Julia sets for the complex family fc = z2 + c. What happens to the filled Julia sets as you move c from inside the Mandelbrot set up toward the boundary, then across the boundary and out beyond the Mandelbrot set? Describe how the Julia sets change as you “walk” along the edge of the Mandelbrot set.

256

2. Repeat Problem 1 for the complex family gc = csin

3. Repeat Problem 1 for the family hc = cez.

Exploration 13.4

GLOSSARY

Acceleration The acceleration of a moving body whose position at time t is u( t) is given by

d2u

dfi

Air resistance A body moving through air (or some other medium) is slowed down by a resistive force (also called a drag or damping force) that acts opposite to the body’s velocity. See also “Viscous damping” and “Newtonian damping.”

Amplitude The amplitude of a periodic oscillating function u(t) is half the difference between its maximum and minimum values.

Angular momentum The angular momentum vector of a body rotating about an axis is its moment of inertia about the axis times its angular velocity vector.

This is the analog in rotational mechanics of momentum (mass times velocity) in linear mechanics.

Angular velocity An angular velocity vector, m(t), is the key to the relation between rotating body axes and a fixed coordinate system of the observer. The component mj of the vector m(t) along the jth body axis describes the spin rate of the body about that axis.

Autocatalator This is a chemical reaction of several steps, at least one of which is autocatalytic.

Autocatalytic reaction In an autocatalytic reaction, a chemical species stimulates more of its own production than is destroyed in the process.

Autonomous ODE An autonomous ODE has no explicit mention of the independent variable (usually t) in the rate equations. For example, 2 = x2 is autonomous, but x = x2 + t is not.

**476**> 477 478 479 480 481 482 .. 486 >> Next