Applications of Differential Equations

A differential equation is not really about $x$ and $y$ — it is a statement about change. It says: "the rate at which this quantity changes right now depends on how much of it there is right now." That single idea turns out to describe cooling coffee, spreading rumors, charging capacitors, decaying isotopes, vibrating bridges, and shrinking populations of fish. When you learn to read a physical situation and translate "the rate of change is proportional to ..." into an equation, you gain one of the most powerful modeling tools in all of science.

This page walks through four of the classic, must-know applications — logistic population growth, Newton's law of cooling, RC electrical circuits, and mechanical vibrations — and shows you not just the formulas but why each equation looks the way it does and how to solve it end to end.

Learning Objectives

• Translate a described physical process ("rate of change is proportional to...") into a differential equation.
• Set up and solve the logistic growth model and interpret the carrying capacity.
• Apply Newton's law of cooling to temperature problems and find unknown times.
• Model RC circuits and compute charging/discharging behavior and time constants.
• Set up the equation for mechanical vibrations and classify the motion (undamped, under-, over-, and critically damped).
• Recognize the shared mathematical structure behind physically different phenomena.

Quick Answer

Differential equations model any system where change depends on the current state. Logistic growth, $\frac{dP}{dt} = rP\left(1 - \frac{P}{K}\right)$, describes populations that grow fast when small but level off at a carrying capacity $K$. Newton's law of cooling, $\frac{dT}{dt} = -k(T - T_s)$, says an object's temperature approaches its surroundings exponentially. RC circuits obey $R\frac{dq}{dt} + \frac{q}{C} = V$, giving exponential charging and discharging governed by the time constant $\tau = RC$. Mechanical vibrations follow $m\frac{d^2x}{dt^2} + c\frac{dx}{dt} + kx = F(t)$, a second-order equation whose damping term $c$ decides whether the system oscillates. The remarkable lesson is that these four unrelated fields share the same handful of equation types.

Where It Came From

The story begins with Isaac Newton and Gottfried Leibniz in the late 1600s. Newton did not invent calculus for its own sake — he invented it because he needed to describe motion. To state his second law, $F = ma$, precisely, he needed a language for instantaneous rate of change, and acceleration is literally the second derivative of position. From the very start, then, calculus and differential equations were the native language of physics.

Through the 1700s a spectacular generation of mathematicians — the Bernoullis, Leonhard Euler, Joseph-Louis Lagrange, and Pierre-Simon Laplace — realized that almost every physical law, once stated as "this rate depends on that quantity," became a differential equation. Euler's work on vibrating strings and beams, and later Jean le Rond d'Alembert's wave equation, showed that the same mathematics governed sound, elasticity, and fluid flow. Newton himself, in Principia, effectively used the law of cooling in his study of how bodies lose heat.

The nineteenth century made the language universal. When physicists needed to describe heat (Fourier), electricity and magnetism (Maxwell), and later quantum systems (Schrödinger), each of them arrived at a differential equation. The logistic equation came from a different need: in 1838 the Belgian mathematician Pierre François Verhulst was dissatisfied with Thomas Malthus's prediction of unlimited exponential population growth, which clearly could not continue forever. Verhulst added a braking term to capture the finite resources of the real world. The deep reason differential equations became the core language of engineering and physics is simple: nature's laws are almost never about how much there is, but about how fast it changes — and that is exactly what a derivative expresses.

Logistic Population Growth

Pure exponential growth, $\frac{dP}{dt} = rP$, assumes each individual reproduces at the same rate no matter how crowded things get. Real populations run into limited food, space, and other resources. Verhulst's fix multiplies the growth rate by a factor that shrinks as the population approaches a maximum sustainable size $K$, called the carrying capacity:

$\frac{dP}{dt} = rP\left(1 - \frac{P}{K}\right)$

When $P$ is small, $1 - P/K \approx 1$ and growth is nearly exponential. As $P \to K$, the bracket $\to 0$ and growth stalls. The solution (by separation of variables and partial fractions) is the S-shaped logistic curve:

$P(t) = \frac{K}{1 + A e^{-rt}}, \qquad A = \frac{K - P_0}{P_0}$

Worked example. A fish population in a lake has carrying capacity $K = 1000$, intrinsic growth rate $r = 0.4$ per year, and starts at $P_0 = 100$. Find the population after 5 years.

First, $A = \frac{1000 - 100}{100} = 9$. Then

$P(5) = \frac{1000}{1 + 9 e^{-0.4 \times 5}} = \frac{1000}{1 + 9 e^{-2}}.$

Since $e^{-2} \approx 0.1353$, the denominator is $1 + 9(0.1353) = 1 + 1.218 = 2.218$, so

$P(5) = \frac{1000}{2.218} \approx 451 \text{ fish}.$

The population has more than quadrupled but is nowhere near $K$ yet — it is in the steepest part of the S-curve.

Newton's Law of Cooling

Newton observed that a hot object loses heat at a rate proportional to how much hotter it is than its surroundings. If $T(t)$ is the object's temperature and $T_s$ the (constant) ambient temperature:

$\frac{dT}{dt} = -k(T - T_s), \qquad k > 0.$

The minus sign encodes that a hotter object ($T > T_s$) cools down. Separating variables gives the solution:

$T(t) = T_s + (T_0 - T_s) e^{-kt}.$

The temperature difference decays exponentially; the object asymptotically approaches room temperature.

Worked example. A cup of coffee at $90^\circ\text{C}$ sits in a $20^\circ\text{C}$ room. After 5 minutes it has cooled to $70^\circ\text{C}$. When will it reach $40^\circ\text{C}$?

Here $T_0 = 90$, $T_s = 20$, so $T_0 - T_s = 70$. Use the 5-minute data to find $k$:

$70 = 20 + 70 e^{-5k} \;\Rightarrow\; 50 = 70 e^{-5k} \;\Rightarrow\; e^{-5k} = \tfrac{5}{7}.$

So $-5k = \ln(5/7) = -0.3365$, giving $k = 0.0673$ per minute. Now solve for the time $t$ when $T = 40$:

$40 = 20 + 70 e^{-kt} \;\Rightarrow\; e^{-kt} = \tfrac{20}{70} = 0.2857.$

Then $-kt = \ln(0.2857) = -1.2528$, so $t = \frac{1.2528}{0.0673} \approx 18.6$ minutes. The coffee reaches $40^\circ\text{C}$ about 18 to 19 minutes after it was poured.

RC Electrical Circuits

A resistor ($R$) and capacitor ($C$) in series with a voltage source is one of the most important circuits in electronics. Kirchhoff's voltage law says the source voltage equals the sum of the voltage drops across the resistor ($iR$) and the capacitor ($q/C$). Because current is the rate of change of charge, $i = \frac{dq}{dt}$:

$R\frac{dq}{dt} + \frac{q}{C} = V.$

This is a first-order linear ODE. For a capacitor charging from $q(0) = 0$ toward a battery voltage $V$, the solution is:

$q(t) = CV\left(1 - e^{-t/RC}\right), \qquad \tau = RC.$

The quantity $\tau = RC$ is the time constant: after one $\tau$ the capacitor reaches about $63\%$ of full charge; after $5\tau$ it is over $99\%$ charged, effectively "full." For discharging (no source, initial charge $q_0$), the equation becomes $R\frac{dq}{dt} + \frac{q}{C} = 0$ with solution $q(t) = q_0 e^{-t/RC}$.

Worked example. A circuit has $R = 2000\ \Omega$, $C = 100\ \mu\text{F} = 100 \times 10^{-6}\ \text{F}$, and a $12\ \text{V}$ battery. The time constant is

$\tau = RC = 2000 \times 100 \times 10^{-6} = 0.2 \text{ s}.$

The final charge is $CV = 100 \times 10^{-6} \times 12 = 1.2 \times 10^{-3}\ \text{C} = 1.2\ \text{mC}$. After $t = 0.2\ \text{s}$ (one time constant):

$q = 1.2 \times 10^{-3}\left(1 - e^{-1}\right) = 1.2 \times 10^{-3}(1 - 0.3679) = 7.585 \times 10^{-4}\ \text{C},$

about $0.76\ \text{mC}$, which is indeed $63\%$ of the full $1.2\ \text{mC}$.

Mechanical Vibrations

A mass on a spring is the archetype of oscillation. Newton's second law balances three forces: the spring's restoring force ($-kx$, Hooke's law), a damping (friction) force proportional to velocity ($-c\frac{dx}{dt}$), and any external driving force $F(t)$:

$m\frac{d^2x}{dt^2} + c\frac{dx}{dt} + kx = F(t).$

This second-order linear ODE governs an enormous range of systems: car suspensions, building responses to earthquakes, and even electrical LC circuits (which are mathematically identical). The behavior depends entirely on the damping, judged by the discriminant of the characteristic equation $m\lambda^2 + c\lambda + k = 0$:

• Undamped ($c = 0$): pure sinusoidal oscillation forever at natural frequency $\omega_0 = \sqrt{k/m}$.
• Underdamped ($c^2 < 4mk$): oscillation with decaying amplitude — a car door swinging shut with a couple of bounces.
• Critically damped ($c^2 = 4mk$): returns to rest as fast as possible without oscillating — the design goal for a good shock absorber.
• Overdamped ($c^2 > 4mk$): returns to rest slowly, no oscillation — like a screen door with a heavy closer.

Worked example (undamped). A $2\ \text{kg}$ mass on a spring with stiffness $k = 8\ \text{N/m}$, no damping, is pulled $0.1\ \text{m}$ and released from rest. The equation is $2x'' + 8x = 0$, i.e. $x'' + 4x = 0$. The natural frequency is $\omega_0 = \sqrt{4} = 2\ \text{rad/s}$. With $x(0) = 0.1$ and $x'(0) = 0$, the solution is

$x(t) = 0.1 \cos(2t).$

The mass oscillates with amplitude $0.1\ \text{m}$ and period $T = \frac{2\pi}{\omega_0} = \frac{2\pi}{2} = \pi \approx 3.14\ \text{s}$, returning to its starting position every $\pi$ seconds.

Real-World Applications

• Epidemiology. The logistic equation and its cousins (SIR models) predict how diseases spread and level off, guiding vaccination targets.
• Pharmacology. Drug concentration in the bloodstream decays much like Newton's cooling, setting safe dosing intervals.
• Forensic science. Investigators use Newton's law of cooling on a body's temperature to estimate time of death.
• Electronics and signal processing. RC circuits form the timing elements in clocks, filters that remove noise, and the smoothing stages of power supplies.
• Structural and automotive engineering. Vibration analysis prevents resonance disasters (like the Tacoma Narrows Bridge) and tunes suspensions toward critical damping for a smooth, safe ride.
• Ecology and fisheries management. Carrying-capacity models set sustainable harvest quotas.

Common Mistakes

1. Confusing exponential with logistic growth. Misconception: populations grow exponentially forever. Why wrong: resources are finite, so growth must slow. Correction: use the logistic model when there is a carrying capacity; exponential is only a good approximation while $P \ll K$.

2. Getting the sign wrong in cooling. Misconception: writing $\frac{dT}{dt} = k(T - T_s)$ with a plus sign. Why wrong: that predicts a hot object gets hotter, which is nonsense. Correction: the rate must be negative when $T > T_s$, so the equation is $\frac{dT}{dt} = -k(T - T_s)$ with $k > 0$.

3. Mishandling units in RC circuits. Misconception: plugging in $C = 100$ when the capacitance is $100\ \mu\text{F}$. Why wrong: microfarads are $10^{-6}\ \text{F}$; using raw numbers inflates $\tau$ by a million. Correction: always convert to base SI units (farads, ohms, seconds) before computing $\tau = RC$.

4. Assuming vibrations always oscillate. Misconception: a mass-spring system always bounces. Why wrong: enough damping ($c^2 \geq 4mk$) kills all oscillation. Correction: check the discriminant before deciding whether the solution has a $\sin$/$\cos$ part.

Comparison and Connections

Although these phenomena come from biology, thermodynamics, electronics, and mechanics, their equations rhyme. The first three are first-order; vibrations is second-order, which is why only it can oscillate on its own.

Application	Governing equation	Order	Typical behavior
Logistic growth	$\frac{dP}{dt} = rP(1 - P/K)$	1st (nonlinear)	S-curve to $K$
Newton cooling	$\frac{dT}{dt} = -k(T - T_s)$	1st (linear)	Exponential decay to $T_s$
RC circuit	$R\frac{dq}{dt} + q/C = V$	1st (linear)	Exponential charge/discharge
Vibrations	$mx'' + cx' + kx = F$	2nd (linear)	Oscillation, possibly damped

Notice that cooling and the discharging RC circuit are the same equation, $y' = -ky$, wearing different physical clothes. This is the great economy of differential equations: master a handful of equation types and you understand dozens of fields.

Practice Questions

Recall

State the logistic differential equation and identify what $r$ and $K$ represent.

Answer: $\frac{dP}{dt} = rP(1 - P/K)$; $r$ is the intrinsic (small-population) growth rate and $K$ is the carrying capacity, the maximum sustainable population.

Understanding

Explain why Newton's law of cooling has a negative sign but the temperature never goes below the surroundings.

Guidance: The factor $(T - T_s)$ shrinks toward zero as $T \to T_s$. When $T = T_s$ the rate is zero, so cooling stops exactly at ambient — the object asymptotes to $T_s$ rather than overshooting.

Application

A capacitor with $R = 500\ \Omega$ and $C = 2\ \mu\text{F}$ discharges from $q_0 = 4\ \mu\text{C}$. Find the charge after $1\ \text{ms}$.

Answer: $\tau = RC = 500 \times 2\times 10^{-6} = 10^{-3}\ \text{s}$. So $q = 4\times 10^{-6} e^{-1\times 10^{-3}/10^{-3}} = 4\times 10^{-6} e^{-1} \approx 1.47\ \mu\text{C}$.

Analysis

A $1\ \text{kg}$ mass sits on a spring with $k = 16\ \text{N/m}$. What damping coefficient $c$ makes the system critically damped, and why might an engineer want exactly this value?

Answer: Critical damping needs $c^2 = 4mk = 4(1)(16) = 64$, so $c = 8\ \text{N·s/m}$. Engineers target it because the system then returns to equilibrium in the shortest possible time with no overshoot — ideal for shock absorbers and instrument needles.

FAQ

Do I always have to solve the differential equation to use it? No. Often the equation itself answers the question — for example, setting $\frac{dP}{dt} = 0$ finds equilibrium populations ($P = 0$ or $P = K$) without any integration.

Why is $RC$ called a "time constant" — it's a product of ohms and farads, not seconds? Work out the units: ohms $\times$ farads $=$ (volts/amp) $\times$ (coulombs/volt) $=$ coulombs/amp $=$ seconds. The units genuinely come out to time.

How do I know whether to use exponential or logistic growth? If the population is far below any resource limit and you only need a short-term forecast, exponential is fine. If you care about long-term behavior or the population is a sizeable fraction of $K$, use logistic.

What's the difference between the damping cases in real life? Underdamped bounces a few times (a plucked guitar string). Critically damped snaps to rest fastest without bouncing (a good car shock). Overdamped oozes back slowly (a hydraulic door closer set too stiff).

Why does the same equation describe cooling and a discharging capacitor? Both are systems where the rate of decrease is proportional to the current amount — $y' = -ky$. Nature reuses this pattern constantly, which is exactly why learning ODEs pays off across so many subjects.

Quick Revision

• Logistic: $\frac{dP}{dt} = rP(1 - P/K)$, solution $P = \frac{K}{1 + Ae^{-rt}}$, $A = (K - P_0)/P_0$.
• Cooling: $\frac{dT}{dt} = -k(T - T_s)$, solution $T = T_s + (T_0 - T_s)e^{-kt}$.
• RC charging: $q = CV(1 - e^{-t/RC})$; discharging: $q = q_0 e^{-t/RC}$; time constant $\tau = RC$.
• Vibrations: $mx'' + cx' + kx = F(t)$; natural frequency $\omega_0 = \sqrt{k/m}$; critical damping at $c^2 = 4mk$.
• Cooling and RC discharge are both $y' = -ky$ — pure exponential decay.

Related Topics

Prerequisites

• Introduction to Differential Equations
• Separable Equations
• Linear First-Order Equations

Related Topics

• Second-Order Linear Equations
• Differential Equations overview

Next Topics

• Laplace Transforms
• Systems of differential equations and predator-prey models