Separation of Variables

Separation of Variables is a technique for solving first-order differential equations of the form

\frac{d y}{d x} = f (x) g (y)

by algebraically gathering all

y

-terms on one side and all

x

-terms on the other, then integrating both sides. Key result: a separable ODE is solved by

\int \frac{d y}{g ( y )} = \int f (x) d x + C

. Applies to: population growth, radioactive decay, Newton's cooling, RC circuits, and any process whose rate depends on the product of a function of the input and a function of the output. Who/when: rooted in the work of Leibniz, the Bernoullis, and Euler in the late 17th–18th centuries.

Separation of variables is the most elementary and widely used method for solving ordinary differential equations. A differential equation is called separable when its derivative can be written as a product of a function depending only on the independent variable and a function depending only on the dependent variable. The method matters because an enormous range of natural phenomena — anything growing or decaying at a rate proportional to its current size — produces separable equations. The surprising fact is that we treat the symbol $\frac{d y}{d x}$ , which is formally a limit of a ratio rather than a true fraction, as if it were a fraction by multiplying through by $d x$ — and remarkably, this informal manipulation always yields the correct answer, justified rigorously by the chain rule.

The Core Insight

If a rate of change factors cleanly into "a part about $x$ " times "a part about $y$ ," then the two variables can be untangled. Once untangled, each side involves only one variable, so ordinary integration — a tool we already possess — finishes the job. The entire difficulty of the differential equation collapses into two separate, independent integrals.

Formal Statement

Separable First-Order ODE

A first-order ordinary differential equation is separable if it can be written in the form

\frac{d y}{d x} = f (x) g (y)

⟳This is the defining form — no derivation, but note that the right side must factor into one pure-

x

function and one pure-

y

function with nothing mixing them.

where:

$x$ = the independent variable,
$y = y (x)$ = the unknown function we seek,
$\frac{d y}{d x}$ = the derivative (rate of change of $y$ with respect to $x$ ),
$f (x)$ = a function depending only on $x$ ,
$g (y)$ = a function depending only on $y$ .

Equivalently, in differential form,

\frac{1}{g ( y )} d y = f (x) d x

⟳Divide both sides of [1] by

g (y)

(assuming

g (y) \neq = 0

) and multiply by

d x

, treating

d y / d x

formally as a quotient of differentials.

Solution by Integration

Given the separable equation [1] with $g (y) \neq = 0$ , the general solution is obtained implicitly by

\int \frac{1}{g ( y )} d y = \int f (x) d x + C

⟳Integrate both sides of [2]; the single arbitrary constant

C

absorbs the two constants of integration from the two indefinite integrals.

where $C$ is an arbitrary constant determined by an initial condition $y (x_{0}) = y_{0}$ .

Why It Holds / Derivation

Justification via the Chain Rule

The "multiply by $d x$ " move looks like illegal fraction abuse. Here is the rigorous version that vindicates it.

Step 1. Start from the separable equation, dividing by

g (y)

\frac{1}{g ( y )} \frac{d y}{d x} = f (x)

⟳Divide both sides of [1] by

g (y)

; this is valid wherever

g (y) \neq = 0

Step 2. Define an antiderivative

G

\frac{1}{g}

, so that

G^{'} (y) = \frac{1}{g ( y )}

. By the chain rule,

\frac{d}{d x} G (y (x)) = G^{'} (y) \frac{d y}{d x} = \frac{1}{g ( y )} \frac{d y}{d x}

⟳Apply the chain rule to the composite function

G (y (x))

and substitute

G^{'} = 1/ g

Step 3. Combining [4] and [5], the left side of [4] is the

x

-derivative of

G (y (x))

\frac{d}{d x} G (y (x)) = f (x)

⟳Substitute the result of [5] into the left-hand side of [4].

Step 4. Integrate both sides with respect to

x

. Let

F

be an antiderivative of

f

G (y (x)) = F (x) + C

⟳Integrate [6] over

x

; the fundamental theorem of calculus undoes the

\frac{d}{d x}

on the left and produces

F (x) + C

on the right.

Since $G (y) = \int \frac{1}{g ( y )} d y$ and $F (x) = \int f (x) d x$ , equation [7] is exactly the claimed result [3]. The informal "separate and integrate" recipe is therefore a shorthand for this chain-rule argument. $■$

Gottfried Wilhelm Leibniz | 1646–1716

Leibniz, co-inventor of calculus, introduced the very notation $\frac{d y}{d x}$ and the integral sign $\int$ that make separation of variables so visually natural. His differential notation was deliberately designed so that the chain rule and substitution look like fraction arithmetic — the "multiply by $d x$ " step that students perform is a direct legacy of Leibniz's vision that differentials could be manipulated almost algebraically. He treated $d y$ and $d x$ as genuine infinitesimal quantities, a viewpoint later made rigorous by the chain rule and, much later, by non-standard analysis.

Building Intuition

Imagine $\frac{d y}{d x} = f (x) g (y)$ as a recipe for the steepness of a curve at every point of the plane. The factor $g (y)$ says "how strongly does the current value of $y$ influence the slope?" and $f (x)$ says "how strongly does our horizontal position influence the slope?". When these two influences multiply cleanly — never tangling into terms like $x y + sin (x + y)$ — you can ask the two questions separately. Accumulating all the tiny $y$ -changes (weighted by $1/ g (y)$ ) on the left must equal accumulating all the tiny $x$ -pushes (weighted by $f (x)$ ) on the right. Integration is just "totalling up" these tiny contributions, and the constant $C$ records where you started.

Slope field of $dy/dx = -y$ with one separated solution curve $y=e^{-x}$

Worked Examples

Radioactive Decay

Given: A sample of carbon-14 decays at a rate proportional to the amount present:

\frac{d N}{d t} = - k N

, with

k = 1.21 \times 1 0^{- 4} yr^{- 1}

and initial amount

N_{0} = 100 g

. Find

N (t)

and the amount after

5000

years. Step 1 — Separate. Move all

N

to one side, all

t

to the other:

\frac{1}{N} d N = - k d t

⟳Divide both sides of

d N / d t = - k N

N

and multiply by

d t

Step 2 — Integrate both sides.

ln ∣ N ∣ = - k t + C

⟳

\int \frac{1}{N} d N = ln ∣ N ∣

and

\int - k d t = - k t

; combine constants into

C

Step 3 — Solve for $N$ and apply $N (0) = N_{0}$ .

N (t) = N_{0} e^{- k t}

⟳Exponentiate [9] to get

N = e^{C} e^{- k t}

; set

t = 0

to identify

e^{C} = N_{0}

Step 4 — Plug in numbers.

N (5000) = 100 e^{- (1.21 \times 1 0^{- 4}) (5000)} = 100 e^{- 0.605}

N (5000) \approx 54.6 g

Newton's Law of Cooling

Given: A cup of coffee at

9 0^{\circ} C

cools in a room at

2 0^{\circ} C

according to

\frac{d T}{d t} = - k (T - 20)

with

k = 0.05 min^{- 1}

. Find the temperature after

30

minutes. Step 1 — Separate. Let

u = T - 20

d u = d T

\frac{1}{T - 20} d T = - k d t

⟳Divide the cooling equation by

(T - 20)

and multiply by

d t

Step 2 — Integrate.

ln ∣ T - 20∣ = - k t + C

⟳

\int \frac{d T}{T - 20} = ln ∣ T - 20∣

via the substitution

u = T - 20

Step 3 — Solve and apply $T (0) = 90$ .

T (t) = 20 + 70 e^{- k t}

⟳Exponentiate [12] to get

T - 20 = A e^{- k t}

; at

t = 0

T = 90

gives

A = 70

Step 4 — Plug in numbers.

T (30) = 20 + 70 e^{- (0.05) (30)} = 20 + 70 e^{- 1.5} = 20 + 70 (0.2231)

T (30) \approx 35. 6^{\circ} C

Common Misconceptions

Do Not Forget the Constant — Or the Lost Solutions

Two classic errors:

Dropping $C$ . The constant of integration is not optional; it encodes the initial condition. Without it you get one specific curve instead of the whole family of solutions.
Dividing by $g (y)$ without checking $g (y) = 0$ . When you divide [1] by $g (y)$ you silently assume $g (y) \neq = 0$ . Any constant $y$ with $g (y) = 0$ is an equilibrium (singular) solution that the integration procedure may discard. Always check these separately.

A Solution Lost by Careless Division

Consider $\frac{d y}{d x} = y^{2/3}$ . Separating gives $\int y^{- 2/3} d y = \int d x$ , hence $3 y^{1/3} = x + C$ , so $y = (\frac{x + C}{3})^{3}$ . But $g (y) = y^{2/3}$ vanishes at $y = 0$ , and indeed $y \equiv 0$ is also a perfectly valid solution that our separation discarded by dividing by $y^{2/3}$ . Worse, this equation violates uniqueness: through the point $(0, 0)$ both $y \equiv 0$ and $y = (x /3)^{3}$ pass, demonstrating that separation alone does not guarantee a unique solution where $g (y) = 0$ .

Connections

Separation of variables is the gateway technique in the study of Ordinary Differential Equations. It relies entirely on Integration (the inverse of differentiation) and is rigorously justified by the Chain Rule. When an equation is not separable, one turns to neighbouring methods such as Integrating Factors for linear first-order equations or Exact Differential Equations. In physics, the equations it solves are instances of Exponential Growth and Decay. The same name, "separation of variables," also denotes a far more advanced technique for partial differential equations — see Separation of Variables (PDEs) — where a solution is sought as a product $X (x) T (t)$ ; the conceptual spirit of splitting variables is shared, though the machinery differs.

Real-World Applications

Separable differential equations govern a remarkable spread of real phenomena:

Radioactive dating — carbon-14 and uranium-lead decay obey $\frac{d N}{d t} = - k N$ , the basis of archaeological and geological age estimates.
Population biology — the logistic model $\frac{d P}{d t} = r P (1 - \frac{P}{K})$ is separable and predicts S-shaped growth toward a carrying capacity $K$ .
Pharmacokinetics — drug clearance from the bloodstream follows first-order exponential decay, setting safe dosing intervals.
RC electrical circuits — capacitor charging/discharging $\frac{d Q}{d t} = - \frac{Q}{R C}$ yields the time constant $τ = R C$ .
Newton's law of cooling — forensic estimation of time of death from body temperature.
Mixing/dilution problems — concentration of salt in a stirred tank with inflow and outflow.
Terminal velocity — a falling body with drag $\frac{d v}{d t} = g - \frac{c}{m} v$ separates to give the approach to terminal speed.

Historical Development

1671 — Newton, in Method of Fluxions (circulated privately), solves rate equations using series and informal "fluxional" arguments anticipating separation.
1684 — Leibniz publishes the first calculus paper, introducing $d$ and the differential notation that makes the method visually natural.
1691 — Leibniz explicitly uses the phrase and procedure of separating variables in correspondence solving geometric rate problems.
1694 — Johann Bernoulli systematizes the separation technique while solving the catenary and brachistochrone-related curves.
1696 — Jacob Bernoulli applies separation-type reasoning to the equation now bearing his name (the Bernoulli equation), reducible to separable form.
1740s — Leonhard Euler organizes first-order methods, distinguishing separable, linear, and exact types in his foundational analysis texts.
1768–1770 — Euler's Institutiones Calculi Integralis presents separation of variables as a standard, named classification of integrable ODEs.
1820s–1830s — Cauchy formalizes the notion of an initial-value problem and existence of solutions, providing the rigorous backdrop for when separation is valid.
1890s — Picard and Lindelöf prove the existence–uniqueness theorem, clarifying exactly when separation's discarded equilibrium solutions matter.

Subtleties & Limits

Separation of variables is a recognition method, not a universal algorithm: most differential equations are simply not separable, and there is no mechanical procedure to transform an arbitrary equation into separable form. Even when separation succeeds, the resulting integrals $\int \frac{d y}{g ( y )}$ and $\int f (x) d x$ may have no closed-form antiderivative in terms of elementary functions (for instance, $\int e^{- x^{2}} d x$ ), forcing numerical or special-function treatment. Deeper still, the solution emerges implicitly via [7] and may not be solvable for $y$ explicitly; deciding when an implicit relation $G (y) = F (x) + C$ defines a genuine function, and over what domain, connects to the implicit function theorem and to questions of solution existence and uniqueness that remain central in the theory of differential equations.