1.3 Slope Fields and Solution Curves

Consider a differential equation of the form

dydx=f(x, y)
dydx=f(x, y)
(1)

where the right-hand function f(x,y)f(x,y) involves both the independent variable x and the dependent variable y. We might think of integrating both sides in (1) with respect to x, and hence write y(x)=f(x, y(x))dx+Cy(x)=f(x, y(x))dx+C. However, this approach does not lead to a solution of the differential equation, because the indicated integral involves the unknown function y(x) itself, and therefore cannot be evaluated explicitly. Actually, there exists no straightforward procedure by which a general differential equation can be solved explicitly. Indeed, the solutions of such a simple-looking differential equation as y=x2+y2y'=x2+y2 cannot be expressed in terms of the ordinary elementary functions studied in calculus textbooks. Nevertheless, the graphical and numerical methods of this and later sections can be used to construct approximate solutions of differential equations that suffice for many practical purposes.

Slope Fields and Graphical Solutions

There is a simple geometric way to think about solutions of a given differential equation y=f(x,y)y'=f(x,y). At each point (x, y) of the xy-plane, the value of f(x,y)f(x,y) determines a slope m=f(x,y)m=f(x,y). A solution of the differential equation is simply a differentiable function whose graph y=y(x)y=y(x) has this “correct slope” at each point (x, y(x))(x, y(x)) through which it passes—that is, y(x)=f(x,y(x))y'(x)=f(x,y(x)). Thus a solution curve of the differential equation y=f(x,y)y'=f(x,y)—the graph of a solution of the equation—is simply a curve in the xy-plane whose tangent line at each point (x, y) has slope m=f(x,y)m=f(x,y). For instance, Fig. 1.3.1 shows a solution curve of the differential equation y=xyy'=xy together with its tangent lines at three typical points.

FIGURE 1.3.1.

A solution curve for the differential equation y=xyy'=xy together with tangent lines having

  • slope m1=x1y1m1=x1y1 at the point (x1, y1)(x1, y1);

  • slope m2=x2y2m2=x2y2 at the point (x2, y2)(x2, y2); and

  • slope m3=x3y3m3=x3y3 at the point (x3, y3)(x3, y3).

This geometric viewpoint suggests a graphical method for constructing approximate solutions of the differential equation y=f(x,y)y'=f(x,y). Through each of a representative collection of points (x, y) in the plane we draw a short line segment having the proper slope m=f(x,y)m=f(x,y). All these line segments constitute a slope field (or a direction field) for the equation y=f(x,y)y'=f(x,y).

Example 1

Figures 1.3.2 (a)(d) show slope fields and solution curves for the differential equation

dydx=ky
dydx=ky
(2)

with the values k=2, 0.5, 1k=2, 0.5, 1, and 33 of the parameter k in Eq. (2). Note that each slope field yields important qualitative information about the set of all solutions of the differential equation. For instance, Figs. 1.3.2(a) and (b) suggest that each solution y(x) approaches ±± as x+x+ if k>0k>0, whereas Figs. 1.3.2(c) and (d) suggest that y(x)0y(x)0 as x+x+ if k<0k<0. Moreover, although the sign of k determines the direction of increase or decrease of y(x), its absolute value |k||k| appears to determine the rate of change of y(x). All this is apparent from slope fields like those in Fig. 1.3.2, even without knowing that the general solution of Eq. (2) is given explicitly by y(x)=Cekxy(x)=Cekx.

A slope field suggests visually the general shapes of solution curves of the differential equation. Through each point a solution curve should proceed in such a direction that its tangent line is nearly parallel to the nearby line segments of the slope field. Starting at any initial point (a, b), we can attempt to sketch freehand an approximate solution curve that threads its way through the slope field, following the visible line segments as closely as possible.

FIGURE 1.3.2(a)

Slope field and solution curves for y=2yy'=2y.

FIGURE 1.3.2(b)

Slope field and solution curves for y=(0.5)yy'=(0.5)y.

FIGURE 1.3.2(c)

Slope field and solution curves for y=yy'=y.

FIGURE 1.3.2(d)

Slope field and solution curves for y=3yy'=3y.

x/y 44 33 22 11 0 1 2 3 4
44 0 11 22 33 44 55 66 77 88
33 1 0 11 22 33 44 55 66 77
22 2 1 0 11 22 33 44 55 66
11 3 2 1 0 11 22 33 44 55
0 4 3 2 1 0 11 22 33 44
1 5 4 3 2 1 0 11 22 33
2 6 5 4 3 2 1 0 11 22
3 7 6 5 4 3 2 1 0 11
4 8 7 6 5 4 3 2 1 0

FIGURE 1.3.3.

Values of the slope y=xyy'=xy for 4x, y4.4x, y4.

Example 2

Construct a slope field for the differential equation y=xyy'=xy and use it to sketch an approximate solution curve that passes through the point (4,4)(4,4).

Solution

Figure 1.3.3 shows a table of slopes for the given equation. The numerical slope m=xym=xy appears at the intersection of the horizontal x-row and the vertical y-column of the table. If you inspect the pattern of upper-left to lower-right diagonals in this table, you can see that it was easily and quickly constructed. (Of course, a more complicated function f(x, y) on the right-hand side of the differential equation would necessitate more complicated calculations.) Figure 1.3.4 shows the corresponding slope field, and Fig. 1.3.5 shows an approximate solution curve sketched through the point (4,4)(4,4) so as to follow this slope field as closely as possible. At each point it appears to proceed in the direction indicated by the nearby line segments of the slope field.

Although a spreadsheet program (for instance) readily constructs a table of slopes as in Fig. 1.3.3, it can be quite tedious to plot by hand a sufficient number of slope segments as in Fig. 1.3.4. However, most computer algebra systems include commands for quick and ready construction of slope fields with as many line segments as desired; such commands are illustrated in the application material for this section. The more line segments are constructed, the more accurately solution curves can be visualized and sketched. Figure 1.3.6 shows a “finer” slope field for the differential equation y=xyy'=xy of Example 2, together with typical solution curves treading through this slope field.

FIGURE 1.3.4.

Slope field for y=xyy'=xy corresponding to the table of slopes in Fig. 1.3.3.

FIGURE 1.3.5.

The solution curve through (4, 4).(4, 4).

If you look closely at Fig. 1.3.6, you may spot a solution curve that appears to be a straight line! Indeed, you can verify that the linear function y=x1y=x1 is a solution of the equation y=xyy'=xy, and it appears likely that the other solution curves approach this straight line as an asymptote as x+x+. This inference illustrates the fact that a slope field can suggest tangible information about solutions that is not at all evident from the differential equation itself. Can you, by tracing the appropriate solution curve in this figure, infer that y(3)2y(3)2 for the solution y(x) of the initial value problem y=xy, y(4)=4y'=xy, y(4)=4?

FIGURE 1.3.6.

Slope field and typical solution curves for y=xyy'=xy.

Applications of Slope Fields

The next two examples illustrate the use of slope fields to glean useful information in physical situations that are modeled by differential equations. Example 3 is based on the fact that a baseball moving through the air at a moderate speed v (less than about 300 ft/s) encounters air resistance that is approximately proportional to v. If the baseball is thrown straight downward from the top of a tall building or from a hovering helicopter, then it experiences both the downward acceleration of gravity and an upward acceleration of air resistance. If the y-axis is directed downward, then the ball’s velocity v=dy/dtv=dy/dt and its gravitational acceleration g=32 ft/s2g=32 ft/s2 are both positive, while its acceleration due to air resistance is negative. Hence its total acceleration is of the form

dvdt=gkv.
dvdt=gkv.
(3)

A typical value of the air resistance proportionality constant might be k=0.16k=0.16.

Example 3

Falling baseball Suppose you throw a baseball straight downward from a helicopter hovering at an altitude of 3000 feet. You wonder whether someone standing on the ground below could conceivably catch it. In order to estimate the speed with which the ball will land, you can use your laptop’s computer algebra system to construct a slope field for the differential equation

FIGURE 1.3.7.

Slope field and typical solution curves for v=320.16vv'=320.16v.

dvdt=320.16v.
dvdt=320.16v.
(4)

The result is shown in Fig. 1.3.7, together with a number of solution curves corresponding to different values of the initial velocity v(0) with which you might throw the baseball downward. Note that all these solution curves appear to approach the horizontal line v=200v=200 as an asymptote. This implies that—however you throw it—the baseball should approach the limiting velocity v=200 ft/sv=200 ft/s instead of accelerating indefinitely (as it would in the absence of any air resistance). The handy fact that 60 mi/h=88 ft/s60 mi/h=88 ft/s yields

v=200fts×60 mi/h88 ft/s136.36mih.
v=200fts×60 mi/h88 ft/s136.36mih.

Perhaps a catcher accustomed to 100 mi/h fastballs would have some chance of fielding this speeding ball.

Comment

If the ball’s initial velocity is v(0)=200v(0)=200, then Eq. (4) gives v(0)=32(0.16)(200)=0v'(0)=32(0.16)(200)=0, so the ball experiences no initial acceleration. Its velocity therefore remains unchanged, and hence v(t)200v(t)200 is a constant “equilibrium solution” of the differential equation. If the initial velocity is greater than 200, then the initial acceleration given by Eq. (4) is negative, so the ball slows down as it falls. But if the initial velocity is less than 200, then the initial acceleration given by (4) is positive, so the ball speeds up as it falls. It therefore seems quite reasonable that, because of air resistance, the baseball will approach a limiting velocity of 200 ft/s—whatever initial velocity it starts with. You might like to verify that—in the absence of air resistance—this ball would hit the ground at over 300 mi/h.

In Section 2.1 we will discuss in detail the logistic differential equation

dPdt=kP(MP)
dPdt=kP(MP)
(5)

that often is used to model a population P(t) that inhabits an environment with carrying capacity M. This means that M is the maximum population that this environment can sustain on a long-term basis (in terms of the maximum available food, for instance).

Example 4

Limiting population If we take k=0.0004k=0.0004 and M=150M=150, then the logistic equation in (5) takes the form

dPdt=0.0004P(150P)=0.06P0.0004P2.
dPdt=0.0004P(150P)=0.06P0.0004P2.
(6)

FIGURE 1.3.8.

Slope field and typical solution curves for P=0.06P0.0004P2P'=0.06P0.0004P2.

The positive term 0.06P on the right in (6) corresponds to natural growth at a 6%annual rate (with time t measured in years). The negative term 0.0004P20.0004P2 represents the inhibition of growth due to limited resources in the environment.

Figure 1.3.8 shows a slope field for Eq. (6), together with a number of solution curves corresponding to possible different values of the initial population P(0). Note that all these solution curves appear to approach the horizontal line P=150P=150 as an asymptote. This implies that—whatever the initial population—the population P(t) approaches the limiting population P=150P=150 as tt.

Comment

If the initial population is P(0)=150P(0)=150, then Eq. (6) gives

P(0)=0.0004(150)(150150)=0,
P'(0)=0.0004(150)(150150)=0,

so the population experiences no initial (instantaneous) change. It therefore remains unchanged, and hence P(t)150P(t)150 is a constant “equilibrium solution” of the differential equation. If the initial population is greater than 150, then the initial rate of change given by (6) is negative, so the population immediately begins to decrease. But if the initial population is less than 150, then the initial rate of change given by (6) is positive, so the population immediately begins to increase. It therefore seems quite reasonable to conclude that the population will approach a limiting value of 150—whatever the (positive) initial population.

Existence and Uniqueness of Solutions

Before one spends much time attempting to solve a given differential equation, it is wise to know that solutions actually exist. We may also want to know whether there is only one solution of the equation satisfying a given initial condition—that is, whether its solutions are unique.

Example 5

  1. Failure of existence The initial value problem

    y=1x,y(0)=0
    y'=1x,y(0)=0
    (7)

    has no solution, because no solution y(x)=(1/x)dx=ln |x|+Cy(x)=(1/x)dx=ln |x|+C of the differential equation is defined at x=0x=0. We see this graphically in Fig. 1.3.9, which shows a direction field and some typical solution curves for the equation y=1/xy'=1/x. It is apparent that the indicated direction field “forces” all solution curves near the y-axis to plunge downward so that none can pass through the point (0, 0).

    FIGURE 1.3.9.

    Direction field and typical solution curves for the equation y=1/x.y'=1/x.

    FIGURE 1.3.10.

    Direction field and two different solution curves for the initial value problem y=2y, y(0)=0.y'=2y, y(0)=0.

  2. Failure of uniqueness On the other hand, you can readily verify that the initial value problem

    y=2y,y(0)=0
    y'=2y,y(0)=0
    (8)

    has the two different solutions y1(x)=x2y1(x)=x2 and y2(x)0y2(x)0 (see Problem 27). Figure 1.3.10 shows a direction field and these two different solution curves for the initial value problem in (8). We see that the curve y1(x)=x2y1(x)=x2 threads its way through the indicated direction field, whereas the differential equation y=2yy'=2y specifies slope y=0y'=0 along the x-axis y2(x)=0y2(x)=0.

Example 5 illustrates the fact that, before we can speak of “the” solution of an initial value problem, we need to know that it has one and only one solution. Questions of existence and uniqueness of solutions also bear on the process of mathematical modeling. Suppose that we are studying a physical system whose behavior is completely determined by certain initial conditions, but that our proposed mathematical model involves a differential equation not having a unique solution satisfying those conditions. This raises an immediate question as to whether the mathematical model adequately represents the physical system.

The theorem stated below implies that the initial value problem y=f(x,y), y(a)=by'=f(x,y), y(a)=b has one and only one solution defined near the point x=ax=a on the x-axis, provided that both the function f and its partial derivative f/yf/y are continuous near the point (a, b) in the xy-plane. Methods of proving existence and uniqueness theorems are discussed in the Appendix.

FIGURE 1.3.11.

The rectangle R and x-interval I of Theorem 1, and the solution curve y=y(x)y=y(x) through the point (a, b).

Remark 1

In the case of the differential equation dy/dx=ydy/dx=y of Example 1 and Fig. 1.3.2(c), both the function f(x,y)=yf(x,y)=y and the partial derivative f/y=1f/y=1 are continuous everywhere, so Theorem 1 implies the existence of a unique solution for any initial data (a, b). Although the theorem ensures existence only on some open interval containing x=ax=a, each solution y(x)=Cexy(x)=Cex actually is defined for all x.

Remark 2

In the case of the differential equation dy/dx=2ydy/dx=2y of Example 5(b) and Eq. (8), the function f(x,y)=2yf(x,y)=2y is continuous wherever y>0y>0, but the partial derivative f/y=1/yf/y=1/y is discontinuous when y=0y=0, and hence at the point (0, 0). This is why it is possible for there to exist two different solutions y1(x)=x2y1(x)=x2 and y2(x)0y2(x)0, each of which satisfies the initial condition y(0)=0y(0)=0.

Remark 3

In Example 7 of Section 1.1 we examined the especially simple differential equation dy/dx=y2dy/dx=y2. Here we have f(x,y)=y2f(x,y)=y2 and f/y=2yf/y=2y. Both of these functions are continuous everywhere in the xy-plane, and in particular on the rectangle 2<x<2, 0<y<22<x<2, 0<y<2. Because the point (0, 1) lies in the interior of this rectangle, Theorem 1 guarantees a unique solution—necessarily a continuous function—of the initial value problem

dydx=y2,y(0)=1
dydx=y2,y(0)=1
(10)

on some open x-interval containing a=0a=0. Indeed this is the solution

y(x)=11x
y(x)=11x

that we discussed in Example 7. But y(x)=1/(1x)y(x)=1/(1x) is discontinuous at x=1x=1, so our unique continuous solution does not exist on the entire interval 2<x<22<x<2. Thus the solution interval I of Theorem 1 may not be as wide as the rectangle R where f and f/yf/y are continuous. Geometrically, the reason is that the solution curve provided by the theorem may leave the rectangle—wherein solutions of the differential equation are guaranteed to exist—before it reaches the one or both ends of the interval (see Fig. 1.3.12).

FIGURE 1.3.12.

The solution curve through the initial point (0, 1) leaves the rectangle R before it reaches the right side of R.

The following example shows that, if the function f(x,y)f(x,y) and/or its partial derivative f/yf/y fail to satisfy the continuity hypothesis of Theorem 1, then the initial value problem in (9) may have either no solution or many—even infinitely many—solutions.

Example 6

Consider the first-order differential equation

xdydx=2y.
xdydx=2y.
(11)

Applying Theorem 1 with f(x,y)=2y/xf(x,y)=2y/x and f/y=2/xf/y=2/x, we conclude that Eq. (11) must have a unique solution near any point in the xy-plane where x0x0. Indeed, we see immediately by substitution in (11) that

y(x)=Cx2
y(x)=Cx2
(12)

satisfies Eq. (11) for any value of the constant C and for all values of the variable x. In particular, the initial value problem

FIGURE 1.3.13.

There are infinitely many solution curves through the point (0, 0), but no solution curves through the point (0, b) if b0b0.

xdydx=2y,y(0)=0
xdydx=2y,y(0)=0
(13)

has infinitely many different solutions, whose solution curves are the parabolas y=Cx2y=Cx2 illustrated in Fig. 1.3.13. (In case C=0C=0 the “parabola” is actually the x-axis y=0y=0.)

Observe that all these parabolas pass through the origin (0, 0), but none of them passes through any other point on the y-axis. It follows that the initial value problem in (13) has infinitely many solutions, but the initial value problem

xdydx=2y,y(0)=b
xdydx=2y,y(0)=b
(14)

has no solution if b0b0.

Finally, note that through any point off the y-axis there passes only one of the parabolas y=Cx2y=Cx2. Hence, if a0a0, then the initial value problem

xdydx=2y,y(a)=b
xdydx=2y,y(a)=b
(15)

has a unique solution on any interval that contains the point x=ax=a but not the origin x=0x=0. In summary, the initial value problem in (15) has

  • a unique solution near (a, b) if a0a0;

  • no solution if a=0a=0 but b0b0;

  • infinitely many solutions if a=b=0a=b=0.

Still more can be said about the initial value problem in (15). Consider a typical initial point off the y-axis—for instance the point (1,1)(1,1) indicated in Fig. 1.3.14. Then for any value of the constant C the function defined by

y(x)={x2if x0,Cx2if x>0
y(x)={x2Cx2if x0,if x>0
(16)

is continuous and satisfies the initial value problem

xdydx=2y,y(1)=1.
xdydx=2y,y(1)=1.
(17)

FIGURE 1.3.14.

There are infinitely many solution curves through the point (1, 1).(1, 1).

For a particular value of C, the solution curve defined by (16) consists of the left half of the parabola y=x2y=x2 and the right half of the parabola y=Cx2y=Cx2. Thus the unique solution curve near (1,1)(1,1) branches at the origin into the infinitely many solution curves illustrated in Fig. 1.3.14.

We therefore see that Theorem 1 (if its hypotheses are satisfied) guarantees uniqueness of the solution near the initial point (a, b), but a solution curve through (a, b) may eventually branch elsewhere so that uniqueness is lost. Thus a solution may exist on a larger interval than one on which the solution is unique. For instance, the solution y(x)=x2y(x)=x2 of the initial value problem in (17) exists on the whole x-axis, but this solution is unique only on the negative x-axis <x<0<x<0.

1.3 Problems

In Problems 1 through 10, we have provided the slope field of the indicated differential equation, together with one or more solution curves. Sketch likely solution curves through the additional points marked in each slope field.

  1. dydx=ysin xdydx=ysin x

    FIGURE 1.3.15.

  2. dydx=x+ydydx=x+y

    FIGURE 1.3.16.

  3. dydx=ysin xdydx=ysin x

    FIGURE 1.3.17.

  4. dydx=xydydx=xy

    FIGURE 1.3.18.

  5. dydx=yx+1dydx=yx+1

    FIGURE 1.3.19.

  6. dydx=xy+1dydx=xy+1

    FIGURE 1.3.20.

  7. dydx=sin x+sin ydydx=sin x+sin y

    FIGURE 1.3.21.

  8. dydx=x2ydydx=x2y

    FIGURE 1.3.22.

  9. dydx=x2y2dydx=x2y2

    FIGURE 1.3.23.

  10. dydx=x2+sin ydydx=x2+sin y

    FIGURE 1.3.24.

A more detailed version of Theorem 1 says that, if the function f(x,y)f(x,y) is continuous near the point (a, b), then at least one solution of the differential equation y=f(x,y)y'=f(x,y) exists on some open interval I containing the point x=ax=a and, moreover, that if in addition the partial derivative f/yf/y is continuous near (a, b), then this solution is unique on some (perhaps smaller) interval J. In Problems 11 through 20, determine whether existence of at least one solution of the given initial value problem is thereby guaranteed and, if so, whether uniqueness of that solution is guaranteed.

  1. dydx=2x2y2;y(1)=1dydx=2x2y2;y(1)=1

     

  2. dydx=x ln  y;y(1)=1dydx=x ln  y;y(1)=1

     

  3. dydx=3y;y(0)=1dydx=y3;y(0)=1

     

  4. dydx=3y;y(0)=0dydx=y3;y(0)=0

     

  5. dydx=xy;y(2)=2dydx=xy;y(2)=2

     

  6. dydx=xy;y(2)=1dydx=xy;y(2)=1

     

  7. ydydx=x1;y(0)=1ydydx=x1;y(0)=1

     

  8. ydydx=x1;y(1)=0ydydx=x1;y(1)=0

     

  9. dydx=ln (1+y2);y(0)=0dydx=ln (1+y2);y(0)=0

     

  10. dydx=x2y2;y(0)=1dydx=x2y2;y(0)=1

In Problems 21 and 22, first use the method of Example 2 to construct a slope field for the given differential equation. Then sketch the solution curve corresponding to the given initial condition. Finally, use this solution curve to estimate the desired value of the solution y(x).

  1. y=x+yy'=x+y,   y(0)=0y(0)=0;   y(4)=y(4)=?

  2. y=yxy'=yx,   y(4)=0y(4)=0;   y(4)=y(4)=?

Problems 23 and 24 are like Problems 21 and 22, but now use a computer algebra system to plot and print out a slope field for the given differential equation. If you wish (and know how), you can check your manually sketched solution curve by plotting it with the computer.

  1. y=x2+y21y'=x2+y21,   y(0)=0y(0)=0;   y(2)=y(2)=?

  2. y=x+12y2y'=x+12y2,   y(2)=0y(2)=0;   y(2)=y(2)=?

  3. Falling Parachutist You bail out of the helicopter of Example 3 and pull the ripcord of your parachute. Now k=1.6k=1.6 in Eq. (3), so your downward velocity satisfies the initial value problem

    dvdt=321.6v,v(0)=0.
    dvdt=321.6v,v(0)=0.

    In order to investigate your chances of survival, construct a slope field for this differential equation and sketch the appropriate solution curve. What will your limiting velocity be? Will a strategically located haystack do any good? How long will it take you to reach 95% of your limiting velocity?

  4. Deer Population Suppose the deer population P(t) in a small forest satisfies the logistic equation

    dPdt=0.0225P0.0003P2.
    dPdt=0.0225P0.0003P2.

    Construct a slope field and appropriate solution curve to answer the following questions: If there are 25 deer at time t=0t=0 and t is measured in months, how long will it take the number of deer to double? What will be the limiting deer population?

The next seven problems illustrate the fact that, if the hypotheses of Theorem 1 are not satisfied, then the initial value problem y=f(x,y), y(a)=by'=f(x,y), y(a)=b may have either no solutions, finitely many solutions, or infinitely many solutions.

  1.  

    1. Verify that if c is a constant, then the function defined piecewise by

      y(x)={0for xc,(xc)2for x>c
      y(x)={0(xc)2for xc,for x>c

      satisfies the differential equation y=2yy'=2y for all x (including the point x=cx=c). Construct a figure illustrating the fact that the initial value problem y=2y, y(0)=0y'=2y, y(0)=0 has infinitely many different solutions.

    2. For what values of b does the initial value problem y=2y, y(0)=by'=2y, y(0)=b have (i) no solution, (ii) a unique solution that is defined for all x?

  2. Verify that if k is a constant, then the function y(x)kxy(x)kx satisfies the differential equation xy=yxy'=y for all x. Construct a slope field and several of these straight line solution curves. Then determine (in terms of a and b) how many different solutions the initial value problem xy=y, y(a)=bxy'=y, y(a)=b has—one, none, or infinitely many.

  3. Verify that if c is a constant, then the function defined piecewise by

    y(x)={0for xc,(xc)3for x>c
    y(x)={0(xc)3for xc,for x>c

    satisfies the differential equation y=3y2/3y'=3y2/3 for all x. Can you also use the “left half” of the cubic y=(xc)3y=(xc)3 in piecing together a solution curve of the differential equation? (See Fig. 1.3.25.) Sketch a variety of such solution curves. Is there a point (a, b) of the xy-plane such that the initial value problem y=3y2/3, y(a)=by'=3y2/3, y(a)=b has either no solution or a unique solution that is defined for all x? Reconcile your answer with Theorem 1.

    FIGURE 1.3.25.

    A suggestion for Problem 29.

  4. Verify that if c is a constant, then the function defined piecewise by

    y(x)={+1if xc,cos(xc)if c<x<c+π,1if xc+π

    satisfies the differential equation y=1y2 for all x. (Perhaps a preliminary sketch with c=0 will be helpful.) Sketch a variety of such solution curves. Then determine (in terms of a and b) how many different solutions the initial value problem y=1y2, y(a)=b has.

  5. Carry out an investigation similar to that in Problem 30, except with the differential equation y=+1y2. Does it suffice simply to replace cos(xc) with sin(xc) in piecing together a solution that is defined for all x?

  6. Verify that if c>0, then the function defined piecewise by

    y(x)={0if x2c,(x2c)2if x2>c

    satisfies the differential equation y=4xy for all x. Sketch a variety of such solution curves for different values of c. Then determine (in terms of a and b) how many different solutions the initial value problem y=4xy, y(a)=b has.

  7. If c0, verify that the function defined by y(x)=x/(cx1) (with the graph illustrated in Fig. 1.3.26) satisfies the differential equation x2y+y2=0 if x1/c. Sketch a variety of such solution curves for different values of c. Also, note the constant-valued function y(x)0 that does not result from any choice of the constant c. Finally, determine (in terms of a and b) how many different solutions the initial value problem x2y+y2=0, y(a)=b has.

    FIGURE 1.3.26.

    Slope field for x2y+y2=0 and graph of a solution y(x)=x/(cx1).

    1. Use the direction field of Problem 5 to estimate the values at x=1 of the two solutions of the differential equation y=yx+1 with initial values y(1)=1.2 and y(1)=0.8.

    2. Use a computer algebra system to estimate the values at x=3 of the two solutions of this differential equation with initial values y(3)=3.01 and y(3)=2.99.

    The lesson of this problem is that small changes in initial conditions can make big differences in results.

    1. Use the direction field of Problem 6 to estimate the values at x=2 of the two solutions of the differential equation y=xy+1 with initial values y(3)=0.2 and y(3)=+0.2.

    2. Use a computer algebra system to estimate the values at x=2 of the two solutions of this differential equation with initial values y(3)=0.5 and y(3)=+0.5.

    The lesson of this problem is that big changes in initial conditions may make only small differences in results.

1.3 Application Computer-Generated Slope Fields and Solution Curves

Widely available computer algebra systems and technical computing environments include facilities to automate the construction of slope fields and solution curves, as do some graphing calculators (see Figs. 1.3.2729).

The Expanded Applications site at the URL indicated in the margin includes a discussion of Maple™, Mathematica™, and Matlab™ resources for the investigation of differential equations. For instance, the Maple command

with(DEtools):
DEplot(diff(y(x),x)=sin(x-y(x)), y(x), x=-5..5, y=-5..5);

and the Mathematica command

VectorPlot[{1, Sin[x-y]},{x, -5, 5}, {y, -5, 5}]

produce slope fields similar to the one shown in Fig. 1.3.29. Figure 1.3.29 itself was generated with the Matlab program dfield [John Polking and David Arnold, Ordinary Differential Equations Using Matlab, 3rd edition, Hoboken, NJ: Pearson, 2003] that is freely available for educational use (math.rice.edu/~dfield). This web site also provides a stand-alone Java version of that can be used in a web browser. When a differential equation is entered in the dfield setup menu (Fig. 1.3.30), you can (with mouse button clicks) plot both a slope field and the solution curve (or curves) through any desired point (or points). Another freely available and user-friendly Matlab-based ODE package with impressive graphical capabilities is Iode (www.math.uiuc.edu/iode).

FIGURE 1.3.27.

TI-84 Plus CE™ graphing calculator and TI-Nspire™ CX CAS handheld. Screenshot from Texas Instruments Incorporated. Courtesy of Texas Instruments Incorporated.

Modern technology platforms offer even further interactivity by allowing the user to vary initial conditions and other parameters “in real time.” Mathematica’s Manipulate command was used to generate Fig. 1.3.31, which shows three particular solutions of the differential equation dy/dx=sin(xy). The solid curve corresponds to the initial condition y(1)=0. As the “locator point” initially at (1, 0) is dragged—by mouse or touchpad—to the point (0, 3) or (2,2), the solution curve immediately follows, resulting in the dashed curves shown. The TI-Nspire CX CAS has similar capability; indeed, as Fig. 1.3.28 appears on the Nspire display, each of the initial points (0,b) can be dragged to different locations using the Nspire’s touchpad, with the corresponding solution curves being instantly redrawn.

FIGURE 1.3.28.

Slope field and solution curves for the differential equation

dydx=sin(xy)

with initial points (0, b), b=3, 1, 2, 0, 2, 4 and window 5x, y5 on a TI-89 graphing calculator.

FIGURE 1.3.29.

Computer-generated slope field and solution curves for the differential equation y=sin(xy).

FIGURE 1.3.30.

Matlab dfield setup to construct slope field and solution curves for y=sin(xy).

FIGURE 1.3.31.

Interactive Mathematica solution of the differential equation y=sin(xy). The “locator point” corresponding to the initial condition y(1)=0 can be dragged to any other point in the display, causing the solution curve to be automatically redrawn.

Use a graphing calculator or computer system in the following investigations. You might warm up by generating the slope fields and some solution curves for Problems 1 through 10 in this section.

Investigation A: Plot a slope field and typical solution curves for the differential equation dy/dx=sin(xy), but with a larger window than that of Fig. 1.3.29. With 10x10, 10y10, for instance, a number of apparent straight line solution curves should be visible, especially if your display allows you to drag the initial point interactively from upper left to lower right.

  1. Substitute y=ax+b in the differential equation to determine what the coefficients a and b must be in order to get a solution. Are the results consistent with what you see on the display?

  2. A computer algebra system gives the general solution

    y(x)=x2 tan1(x2CxC).

    Plot this solution with selected values of the constant C to compare the resulting solution curves with those indicated in Fig. 1.3.28. Can you see that no value of C yields the linear solution y=xπ/2 corresponding to the initial condition y(π/2)=0? Are there any values of C for which the corresponding solution curves lie close to this straight line solution curve?

Investigation B: For your own personal investigation, let n be the smallest digit in your student ID number that is greater than 1, and consider the differential equation

dydx=1ncos(xny).
  1. First investigate (as in part (a) of Investigation A) the possibility of straight line solutions.

  2. Then generate a slope field for this differential equation, with the viewing window chosen so that you can picture some of these straight lines, plus a sufficient number of nonlinear solution curves that you can formulate a conjecture about what happens to y(x) as x+. State your inference as plainly as you can. Given the initial value y(0)=y0, try to predict (perhaps in terms of y0) how y(x) behaves as x+.

  3. A computer algebra system gives the general solution

    y(x)=1n[x+2 tan1(1xC)].

    Can you make a connection between this symbolic solution and your graphically generated solution curves (straight lines or otherwise)?

..................Content has been hidden....................

You can't read the all page of ebook, please click here login for view all page.
Reset
18.117.192.247