Assignment 20

......!!!!!!!!...................................

......!!!!!!!!...................................

First analyzing the equation qualitatively we draw the following conclusions:

If P is initially 200 then we have at t = 0 the equation

1000 / 200 * dP/dt = 100 - 200 so that

dP/dt = -20,

which tells us that the population initially decreases.

As long as P > 100 the population will continue to decrease, since dP/dt = (100 - P) * (P / 1000) will remain negative as long as P > 100.

As P approaches 100 from 'above'--i.e., from populations greater than 100�the rate of change dP/dt = (100 - P) * P / 1000 approaches 0. So the population cannot decrease to less than 100 if it starts at a value greater than 100.

** The analytical solution will confirm these conclusions: **

This equation is of the form dP/dt = k P ( 1 - P / L). Its solution is obtained by the process in the text (write in the form dP / [ P ( L � P) ] = k / L dt , integrate the left-hand side by separation of variables, etc.).

The solution is P = L / (1 + A e^(-k t) ) with A = (L - P0) / P0.

The present problem has dP/dt = P/1000 * (100 - P) has L = 100 and k = 1/1000, with P0 = 200, so A = (100 - 200) / 200 = -1/2 and its solution is

P = 100 / (1 + (-1/2) e^(-t/1000) ) = 100 / (1 - 1/2 e^(-t/1000) ).

For t = 0 the denominator is 1/2; as t -> infinity e^(-t/1000) -> 0 and the denominator approaches 1.

Thus P is initially 100 / (1/2) = 200, as we already know, and at t increases the population approaches 100 / 1 = 100, with the denominator 1 - � e^(-t/1000) always less than 1 so that the population is always greater than 100.

.................................................

......!!!!!!!!...................................

.................................................

......!!!!!!!!...................................

......!!!!!!!!...................................

.................................................

......!!!!!!!!...................................

......!!!!!!!!...................................

.................................................

......!!!!!!!!...................................

......!!!!!!!!...................................

** dP/dt = P^2 - 6 P is a quadratic function a P^2 + b P + c of P, with a = 1 and b = -6. The function opens upward, with its vertex at P = -b / (2a) = -(-6) / (2 * 1) = 3. At this point P^2 - 6 P = 9 - 18 = -9.

The graph has zeros where P^2 - 6 P = 0, or P(P-6) = 0. The zeros are therefore at P = 0 and P = 6.

Thus between P = 0 and P = 6 the value of dP/dt = P^2 - 6 P is negative. dP/dt = 0 at P =0 and at P = 6 and reaches its minimum at P = 3.

If P(0) = 5 then dP/dt = 5^2 - 6 * 5 = -5 so the population initially decreases. As P decreases it comes closer to P = 3, at which value dP/dt is minimized so that the rate of decrease is greatest, so the P vs. t curve will become steeper and steeper in the downward direction. Up to the P = 3 point the graph of P vs. t will therefore be concave downard.

After P decreases to less than 3 the values of dP/dt begin to increase toward 0, however still remaining negative. The graph of P vs. t will become concave upward, with P approaching zero for large t.

If P(0) = 8 then dP/dt = 8^2 - 6 * 8 = 16 so the population initially increases. As P increases dP/dt also increases rate of increase is increasing, so the P vs. t curve will become steeper and steeper in the upward direction. The graph of P vs. t will therefore remain concave upward. The graph of P vs. t therefore continues to increase at an increasing rate, exceeding all bounds as t -> infinity. **

** This equation is of the form dP/dt = k P (1 - P / L):

P^2 - 6 P = P ( P - 6) = -P(6 - P) = -6 P ( 1 - P/6). So the right-hand side has form k P ( 1 - P / L) with k = -6 and L = 6.

You should know how to derive the solution to dP/dt = k P (1 - P / L), which is P = L / (1 + A e^(-k t) ) with A = (L - P0) / P0.

The solution in this case has A = (6 - P0) / P0.

For P0 = 5 we have A = 1/5 and P = 6 / (1 + 1/5 e^(-(-6) * t) ) = 6 / (1 + 1/5 e^(6 t) ).

The denominator increases without bound as t -> infinity, so that P approaches zero as t increases. **

.................................................

......!!!!!!!!...................................

......!!!!!!!!...................................

** dw/dt = w - wr and dr/dt = -r + wr so that dr/dw = r(w-1) / w(1-r).

For r = 0 (the horizontal axis) this gives us dr/dw = 0, so slope lines along the horizontal axis are horizontal.

As w -> 0 we see that dr/dw -> -infinity so the slope lines are vertical along the vertical axis, with the understanding that the vertical slopes are vertical downward.

As we approach r = 1, w = 1 we see that dr/dw approaches the form 0 / 0 so that the slope near (1, 1) will depend very sensitively on whether we are slightly to the right, slightly to the left, slightly above or slightly below (1, 1). This results in the 'circling' behavior we observe around this point.

Looking at the slope field starting at (3, 1) we see that we move upward and to the left, with the worm population decreasing and the robin population increasing. This continues but with less and less upward movement and more and more leftward movement, indicating a slowing of the growth of the robin population while the worm population continues to decrease quickly.

After awhile the graph is moving directly to the left, perhaps near the point (1.5, 2.5), just after which it begins moving downward and to the left, indicating a decreasing robin population as the worm population continues to be depleted. The movement becomes more and more downward and less and less to the left, indicating that as the robin population declines the worm population is not being decimated as quickly as before.

The graph eventually reaches a vertical downward direction, perhaps around the point (.3, 1), after which it continues moving downward but also to the right. This indicates an increasing worm population as the robins continue to decline.

The graph descends less and less rapidly as the increasing worm population begins to provide food for the robins, who stop dying off as quickly. The graph reaches its 'low point' around maybe (1, .2) before it begins to move upward and to the right once more. This indicates that while the robin population is still small enough to permit an increase in the worm population,there are enough worms to feed a growing robin population.

This increase in both populations continues until the graph returns to the initial point (3, 1), at which point there are enough robins to again begindepleting the worm population and the cycle begins to repeat. **

.......................................................!!!!!!!!...................................

......!!!!!!!!...................................

.................................................

......!!!!!!!!...................................

......!!!!!!!!...................................

.................................................

......!!!!!!!!...................................

......!!!!!!!!...................................

.................................................

&#This looks good. See my notes. Let me know if you have any questions. &#