7 Differential Equations and Dynamic Models

A progressive guide to modeling change with differential equations, analyzing slope fields and equilibria, solving separable equations, and interpreting exponential and logistic dynamics.

Modeling Change with Differential Equations

A differential equation describes a quantity through its rate of change rather than specifying the quantity directly. If y(t)y(t) is a quantity that varies with time tt, then dydt\frac{dy}{dt} is its instantaneous rate of change. A general first-order model has the form

dydt=f(t,y).\frac{dy}{dt}=f(t,y).

This equation says that the current rate depends on the time and state of the system. Applications include population change, radioactive decay, medication concentration, temperature, and motion.

An supplies the starting value, such as y(0)=y0y(0)=y_0. Together, the differential equation and the form an initial-value problem. The desired solution must satisfy both the rate equation and the specified starting value.

A model is meaningful only under its assumptions. The variables, parameters, units, and time interval must correspond to the real system. For example, a population model should use a meaningful population scale and should not predict negative populations.

Takeaway: A differential equation gives the rule for change, while an identifies the particular trajectory that starts from the stated state.

Slope Fields and Equilibria

For a first-order equation

y′=f(t,y),y'=f(t,y),

a places a short segment at each point (t,y)(t,y). The segment at that point has slope f(t,y)f(t,y), so a solution curve follows the local directions shown by the field.

Read the field qualitatively:

  • Upward-tilting segments indicate y′>0y'>0, so solutions are increasing.

  • Downward-tilting segments indicate y′<0y'<0, so solutions are decreasing.

  • Horizontal segments indicate y′=0y'=0; a solution may have a local extremum there or may be constant.

For an autonomous equation, the slope depends only on the dependent variable:

y′=g(y).y'=g(y).

Therefore, every horizontal line has the same slope pattern. Values satisfying g(y)=0g(y)=0 give equilibrium solutions. If the system starts at one of these values, it remains there.

Consider

y′=y(4−y).y'=y(4-y).

The equilibrium values are y=0y=0 and y=4y=4. When 0<y<40<y<4, the derivative is positive, so solutions increase. When y>4y>4, the derivative is negative, so solutions decrease. When y<0y<0, the derivative is also negative. Thus, y=4y=4 attracts nearby positive solutions, while y=0y=0 repels solutions that start just above it.

Takeaway: A provides a visual and qualitative description of solution behavior before an explicit formula is found.

Solving Separable Equations

A can be written as

dydt=f(t)g(y).\frac{dy}{dt}=f(t)g(y).

For values where g(y)≠0g(y)\neq 0, separate the variables:

1g(y) dy=f(t) dt.\frac{1}{g(y)}\,dy=f(t)\,dt.

Then integrate both sides. A reliable procedure is:

  1. Find constant solutions by solving g(y)=0g(y)=0.

  2. Separate the variables.

  3. Integrate both sides.

  4. Solve explicitly for yy, if possible.

  5. Apply the .

The first step matters because dividing by g(y)g(y) can remove constant solutions.

For example, consider

dydt=2ty,y(0)=3.\frac{dy}{dt}=2ty, \qquad y(0)=3.

For y≠0y\neq 0, separation gives

1y dy=2t dt.\frac{1}{y}\,dy=2t\,dt.

After integration,

ln⁡∣y∣=t2+C.\ln|y|=t^2+C.

Exponentiating and incorporating the sign into the constant gives

y=Cet2.y=Ce^{t^2}.

The gives C=3C=3, so the solution is

y(t)=3et2.y(t)=3e^{t^2}.

For t>0t>0, this solution is increasing, and its rate of increase becomes faster because both tt and yy contribute to y′=2tyy'=2ty.

Takeaway: Separation converts a rate equation into two integrals, but constant solutions and the must be handled explicitly.

Exponential Growth and Decay

The assumes that the rate of change is proportional to the current amount:

dPdt=rP.\frac{dP}{dt}=rP.

Separating and solving gives

P(t)=P0ert,P(t)=P_0e^{rt},

where P0=P(0)P_0=P(0). If r>0r>0, the quantity grows; if r<0r<0, it decays.

For growth, the doubling time is

ln⁡2r.\frac{\ln 2}{r}.

For decay written with a positive decay constant kk, so that r=−kr=-k, the half-life is

ln⁡2k.\frac{\ln 2}{k}.

As an example, if a substance begins with 8080 milligrams and has decay constant k=0.12k=0.12 per hour, then

A(t)=80e−0.12t.A(t)=80e^{-0.12t}.

After 55 hours,

A(5)=80e−0.6≈43.9 mg.A(5)=80e^{-0.6}\approx 43.9\text{ mg}.

The model predicts an amount that approaches zero without reaching it in finite time. Its reliability depends on whether the proportional rate remains approximately constant. Unlimited exponential growth is often unrealistic because resources, space, competition, or other constraints eventually become important.

Takeaway: Exponential behavior is appropriate when the relative rate of change is approximately constant, but its assumptions should be checked over the intended time interval.

Logistic Growth and Long-Term Behavior

The incorporates a limiting level that restricts growth:

dPdt=rP(1−PK),\frac{dP}{dt}=rP\left(1-\frac{P}{K}\right),

where P(t)P(t) is the population, r>0r>0 is the intrinsic growth rate, and K>0K>0 is the . The factor 1−PK1-\frac{P}{K} reduces growth as PP approaches KK.

With initial population P(0)=P0P(0)=P_0, the solution is

P(t)=P0Kert(K−P0)+P0ert.P(t)=\frac{P_0K e^{rt}}{(K-P_0)+P_0e^{rt}}.

The equilibrium solutions are P=0P=0 and P=KP=K. If 0<P<K0<P<K, then P′>0P'>0, so the population increases. If P>KP>K, then P′<0P'<0, so it decreases toward KK. For positive populations, KK is stable, while P=0P=0 is unstable because a small positive population moves away from zero.

When PP is small compared with KK, the correction factor is close to 11, so logistic growth resembles exponential growth. Growth is fastest at P=K/2P=K/2. At that point, the graph changes from accelerating growth to decelerating growth, producing the inflection point of the typical S-shaped curve. As tt becomes large,

P(t)→K.P(t)\to K.

Compared with the , the is more suitable when resources impose a meaningful upper limit. The choice between the models should be based on the system and the time interval, not only on which formula is easier to use.

Takeaway: Logistic growth begins approximately exponentially when the population is small, then slows and approaches the .

Interpreting and Evaluating a Model

A complete differential-equation analysis combines calculation with interpretation. Use the following checklist:

  1. Identify each variable and its units, including the meanings of tt, the dependent variable, and its derivative.

  2. Find equilibria by setting the right-hand side equal to zero.

  3. Determine where the derivative is positive or negative.

  4. Use a , sign analysis, or an explicit solution to describe short-term behavior.

  5. Examine long-term behavior through limits, stable equilibria, or possible unbounded growth.

  6. Check that the solution satisfies the .

  7. Evaluate whether the parameters remain reasonable and whether the model's assumptions apply.

  8. Respect the domain of the original quantity and the time interval for which the model was calibrated.

For the exponential equation P′=rPP'=rP, growth continues without an upper bound when r>0r>0. For the logistic equation P′=rP(1−P/K)P'=rP(1-P/K), positive solutions approach KK. The models can therefore agree when PP is small but differ substantially when PP becomes comparable to KK.

A useful final interpretation should state what the solution means in the original context, not merely report a formula. Include the direction of change, the relevant equilibrium or limit, and the assumptions that make the prediction credible.

Takeaway: The strongest model analysis explains what the mathematics predicts, why it predicts it, and when that prediction should be trusted.