Differential Equations: Concepts, Models, and Methods

A progressive guide to modeling change with differential equations, solving separable equations, interpreting slope fields, approximating solutions numerically, and comparing exponential with logistic growth.

Foundations: Functions, Derivatives, and Conditions

A relates an unknown function to one or more of its derivatives. For example, the equation

dydx=3x2\frac{dy}{dx}=3x^2

specifies the rate of change of an unknown function y=f(x)y=f(x). A function is a solution when substitution of the function and its derivatives makes the equation true. Since

ddx(x3+C)=3x2,\frac{d}{dx}(x^3+C)=3x^2,

every function of the form y=x3+Cy=x^3+C is a solution. The constant CC produces a family of solutions rather than one unique function.

An adds a condition such as y(0)=2y(0)=2. That condition identifies one particular member of the solution family.

Differential equations are especially useful for modeling quantities that change. If P(t)P(t) denotes a population, temperature, amount of radioactive material, or concentration, then dP/dtdP/dt represents its instantaneous rate of change.

Takeaway: A describes how a quantity changes, while an initial condition helps determine which solution describes the particular situation.

Building Differential-Equation Models

The main modeling task is to translate a verbal statement about change into a derivative equation. A reliable sequence is:

  1. Define each variable and include its units.

  2. Translate the rate statement into a derivative.

  3. Add an initial condition when one is given.

  4. Solve the equation exactly or approximate its solution.

  5. Interpret the result in context and check its units and behavior.

For example, if a population PP increases at a rate proportional to its current size, the statement becomes

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

Here, kk is a constant of proportionality. More generally, a model can be written as

dPdt=f(t,P),\frac{dP}{dt}=f(t,P),

meaning that the rate can depend on time, the current amount, or both.

A gives the slope or rate of change at each point; it does not directly give the value of the quantity. This distinction explains why graphical tools such as slope fields and numerical procedures such as are useful.

Takeaway: Modeling begins with a precise translation from words about rates into derivatives, followed by solving or approximating and then interpreting the result.

Solving Separable Equations

A can be written in the form

dydx=f(x)g(y).\frac{dy}{dx}=f(x)g(y).

When g(y)g(y) is nonzero, rearrange the equation so that the variables are separated:

1g(y) dy=f(x) dx.\frac{1}{g(y)}\,dy=f(x)\,dx.

Then integrate both sides:

∫1g(y) dy=∫f(x) dx+C.\int\frac{1}{g(y)}\,dy=\int f(x)\,dx+C.

For an

dydx=f(x)g(y),y(x0)=y0,\frac{dy}{dx}=f(x)g(y),\qquad y(x_0)=y_0,

use the following procedure:

  1. Separate the variables.

  2. Integrate both sides and include a constant of integration.

  3. Apply the initial condition.

  4. Solve explicitly when useful, or retain an implicit equation.

  5. Differentiate and substitute back to check the result.

When dividing by an expression involving yy, check whether that expression can equal zero. A constant solution may be lost during the division.

For example, solve

dydx=2xy,y(0)=3.\frac{dy}{dx}=2xy,\qquad y(0)=3.

Separating and integrating gives

1y dy=2x dx,\frac{1}{y}\,dy=2x\,dx,
ln⁡∣y∣=x2+C.\ln|y|=x^2+C.

The solution can be written as y=Cex2y=Ce^{x^2}. Applying y(0)=3y(0)=3 gives C=3C=3, so

y=3ex2.y=3e^{x^2}.

A check confirms that

dydx=3ex2(2x)=2xy.\frac{dy}{dx}=3e^{x^2}(2x)=2xy.

Takeaway: Separation converts one into two integrals, but possible equilibrium or constant solutions must be checked before accepting the final answer.

Visualizing Solution Behavior

A represents the equation

dydx=f(x,y)\frac{dy}{dx}=f(x,y)

by placing a short line segment with slope f(x,y)f(x,y) at each point (x,y)(x,y). A solution curve passing through a point must be tangent to the segment at that point.

To construct a field for

dydx=x−y,\frac{dy}{dx}=x-y,

calculate the slope at selected grid points. At (1,2)(1,2), for instance,

dydx∣(1,2)=1−2=−1.\left.\frac{dy}{dx}\right|_{(1,2)}=1-2=-1.

Draw a short segment with slope −1-1 there, and repeat the process across the grid.

Use the field to identify behavior:

  • Positive slopes indicate that solution curves rise as xx increases.

  • Negative slopes indicate that solution curves fall as xx increases.

  • Zero slopes indicate horizontal tangents.

  • Repeated horizontal levels can reveal constant or equilibrium solutions.

  • The arrangement of segments can show whether curves approach or move away from a level.

For an autonomous equation such as dP/dt=f(P)dP/dt=f(P), the slope depends only on PP, so the pattern repeats horizontally across the field. A solution curve is not one individual segment; it is a smooth path that follows the local directions throughout the plane.

Takeaway: A gives a qualitative picture of solution behavior without requiring an exact formula.

Approximating Solutions Numerically

When an exact solution is inconvenient or unavailable, estimates the solution by taking successive tangent-line steps. For

dydx=f(x,y),y(x0)=y0,\frac{dy}{dx}=f(x,y),\qquad y(x_0)=y_0,

and step size hh, the update rules are

xn+1=xn+h,x_{n+1}=x_n+h,
yn+1=yn+hf(xn,yn).y_{n+1}=y_n+h f(x_n,y_n).

To approximate y(0.3)y(0.3) for

dydx=x+y,qquady(0)=1,\frac{dy}{dx}=x+y,\\qquad y(0)=1,

with h=0.1h=0.1, begin with (x0,y0)=(0,1)(x_0,y_0)=(0,1). The successive calculations are:

  1. At (0,1)(0,1), the slope is 11, so y1=1+0.1(1)=1.1y_1=1+0.1(1)=1.1.

  2. At (0.1,1.1)(0.1,1.1), the slope is 1.21.2, so y2=1.1+0.1(1.2)=1.22y_2=1.1+0.1(1.2)=1.22.

  3. At (0.2,1.22)(0.2,1.22), the slope is 1.421.42, so y3=1.22+0.1(1.42)=1.362y_3=1.22+0.1(1.42)=1.362.

Therefore,

y(0.3)≈1.362.y(0.3)\approx1.362.

A smaller step size generally improves the approximation because each tangent-line step covers a shorter interval, although it requires more calculations. The method can overestimate or underestimate depending on the curvature of the actual solution. Compare the numerical path with the when checking the result.

Takeaway: converts local slope information into a sequence of approximate values.

and Decay

and decay assume that the rate of change is proportional to the current amount:

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

Separating and integrating produces

P(t)=P0ekt,P(t)=P_0e^{kt},

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

For a culture beginning with 500 bacteria and having proportionality constant 0.20.2 per hour,

P(t)=500e0.2t.P(t)=500e^{0.2t}.

After 4 hours,

P(4)=500e0.8≈1112.P(4)=500e^{0.8}\approx1112.

For growth, the doubling time is

Td=ln⁡2k.T_d=\frac{\ln 2}{k}.

For decay, the half-life is

T1/2=ln⁡(1/2)k,T_{1/2}=\frac{\ln(1/2)}{k},

which is positive because k<0k<0.

The exponential model assumes unlimited resources and a constant proportional growth rate. It may describe a short interval well, but it predicts unbounded growth when k>0k>0, which is often unrealistic over a long period.

Takeaway: Exponential models are appropriate when the current amount controls the rate and limiting resources are not yet important.

Logistic Growth and Limited Resources

A modifies by including a , the sustainable level represented by KK:

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

where r>0r>0 is the intrinsic growth rate. The factor 1−P/K1-P/K reduces growth as the population approaches KK.

The sign of the derivative explains the behavior:

  • If 0<P<K0<P<K, then dP/dt>0dP/dt>0, so the population increases.

  • If P>KP>K, then dP/dt<0dP/dt<0, so the population decreases.

  • If P=0P=0 or P=KP=K, then the rate is zero, giving levels.

For a population starting between zero and KK, the curve is typically sigmoidal. Growth begins slowly, becomes fastest near P=K/2P=K/2, and then slows as the population approaches KK. At the midpoint,

dPdt∣P=K/2=rK4.\left.\frac{dP}{dt}\right|_{P=K/2}=\frac{rK}{4}.

The general solution is

P(t)=K1+Ae−rt,P(t)=\frac{K}{1+Ae^{-rt}},

where the initial condition determines AA. If P(0)=P0P(0)=P_0, then

A=K−P0P0,A=\frac{K-P_0}{P_0},

and therefore

P(t)=K1+(K−P0P0)e−rt.P(t)=\frac{K}{1+\left(\frac{K-P_0}{P_0}\right)e^{-rt}}.

For K=1000K=1000, P0=100P_0=100, and r=0.4r=0.4 per year, A=9A=9, so

P(t)=10001+9e−0.4t.P(t)=\frac{1000}{1+9e^{-0.4t}}.

As tt increases, e−0.4te^{-0.4t} approaches zero and P(t)P(t) approaches 1000.

Takeaway: Logistic growth incorporates resource limits: the population grows rapidly near the middle of its range but approaches its over time.

Connecting Models and Representations

The main models differ in the assumptions they make about resources and long-term behavior.

  • The exponential equation is P′=kPP'=kP. It assumes unlimited resources and produces unbounded growth when k>0k>0.

  • The logistic equation is P′=rP(1−P/K)P'=rP(1-P/K). It includes limited resources and approaches the KK.

  • has no finite maximum growth rate within the model.

  • Logistic growth is fastest at P=K/2P=K/2.

  • The has equilibrium levels P=0P=0 and P=KP=K, while the exponential model has the constant solution P=0P=0.

A strong solution process moves among several representations:

  1. Verbal: describe how a quantity changes.

  2. Symbolic: write a such as dy/dx=x+ydy/dx=x+y.

  3. Graphical: inspect a .

  4. Numerical: use to generate approximations.

  5. Algebraic: find an exact solution when possible.

  6. Contextual: interpret the result, including its units and limitations.

For example, the verbal statement that a population grows proportionally to its current size becomes P′=kPP'=kP, which leads to an exponential model. If limited resources matter, the statement must be modified to produce a instead.

Final takeaway: Differential equations connect rates, formulas, graphs, numerical approximations, and real-world interpretations. The most appropriate model depends on what controls the rate of change and on whether the assumptions remain reasonable over the time interval being considered.