CALCULUS 2 • DIFFERENTIAL EQUATIONS

Solving Separable DEs — Solving Separable Differential Equations with Initial Conditions

Learn to isolate variables, integrate both sides, and pin down particular solutions using initial conditions.

Historical Context & Motivation

The study of differential equations arose naturally alongside the development of calculus itself, driven by the desire to model physical phenomena such as planetary motion, fluid flow, and the cooling of heated objects. From the moment Newton and Leibniz formalized the derivative in the late seventeenth century, mathematicians recognized that many laws of nature are most naturally expressed as relationships between a function and its rate of change. The challenge, then, was to reverse-engineer the function from such a relationship—what we now call solving a differential equation. Among the earliest and most tractable classes of DEs were those in which the variables could be cleanly separated onto opposite sides of the equation, a technique that remains one of the first tools taught in any differential equations course.

1687
Newton's Principia
Isaac Newton publishes the Principia Mathematica, implicitly solving differential equations to describe gravitational motion and laying the groundwork for the entire field.
1691
Leibniz's Separation Method
Gottfried Wilhelm Leibniz explicitly describes the technique of separating variables in a differential equation, making it one of the earliest systematic solution methods for DEs.
1748
Euler's Integral Formalization
Leonhard Euler rigorously develops integration techniques and initial-value frameworks that allow separated equations to be solved in closed form and matched to specific physical scenarios.
1820s
Cauchy's Existence Theorems
Augustin-Louis Cauchy proves foundational existence and uniqueness theorems for initial value problems, guaranteeing that under suitable conditions a differential equation with an initial condition has exactly one solution.
Modern Era
Computational & Applied Use
Separable DEs remain central in mathematical modeling across biology (population growth), physics (radioactive decay), chemistry (reaction kinetics), and engineering (RC circuits), often serving as the analytical backbone before numerical methods are needed.

The core question that separable equations address is deceptively simple: given a rate of change expressed as a product of two single-variable functions, can we reconstruct the original relationship between the dependent and independent variables? When an initial condition is also specified—for instance, the population at time zero or the temperature at the moment observation begins—we move from a general solution (a family of curves) to a unique particular solution (a single curve) that exactly fits the given scenario. This lesson walks through the full process from identification to final answer.

Core Principles & Definitions

Before diving into the solution procedure, it is essential to understand the structural requirements that make a differential equation separable and to clarify how initial conditions constrain the solution. A first-order ordinary differential equation of the form dy/dx = f(x, y) is called separable if f(x, y) can be factored as a product g(x)·h(y), where g depends only on x and h depends only on y. This factorization is the key structural property that enables us to move all y-dependent terms to one side and all x-dependent terms to the other, reducing the problem to two independent integration tasks.

1

Separability Condition

A DE dy/dx = f(x, y) is separable when f(x, y) = g(x)·h(y). This allows algebraic rearrangement to isolate y-terms with dy and x-terms with dx on opposite sides.
2

General Solution

After separating and integrating both sides, the result includes an arbitrary constant C, producing a family of solution curves. Each value of C yields a different member of this family.
3

Initial Condition

An initial condition y(x₀) = y₀ specifies that the solution curve must pass through the point (x₀, y₀). Substituting this into the general solution determines the unique value of C.
4

Particular Solution

The particular solution is the single member of the family of curves that satisfies the initial condition. It represents the unique trajectory consistent with the given starting data.
5

Domain Awareness

Division by h(y) during separation is only valid where h(y) ≠ 0. Constant solutions where h(y) = 0 (equilibrium solutions) must be checked separately.
KEY TAKEAWAY
Think of the general solution as a stack of transparent map overlays, each showing a slightly different route through a city. The initial condition is like specifying your starting intersection—it selects exactly one overlay (one route) from the entire stack. Without it, you have infinitely many possible paths; with it, only one path remains.

Visual Explanation — Solution Families and Initial Conditions

The violet curves represent the general solution family y = Ce for various values of C. The cyan curve is the particular solution selected by the initial condition y(0) = 1, which forces C = 1. The amber dot marks the initial point (0, 1).

The diagram above illustrates a fundamental idea: integrating a separable DE produces a one-parameter family of curves indexed by the constant C. Each curve is a valid solution to the differential equation, meaning that if you compute dy/dx along any one of them, you recover the original right-hand side 2xy. The initial condition y(0) = 1 picks out a single curve—the one passing through the point (0, 1)—by fixing C = 1. Notice how the curves fan outward from the origin: for C > 0 they grow exponentially upward, for C < 0 they dive downward, and the trivial solution y = 0 (when C = 0) sits along the x-axis as an equilibrium.

Mathematical Framework

The solution procedure for a separable differential equation with an initial condition follows a systematic sequence of algebraic manipulation and integration. We begin with a first-order ODE that can be written in the form dy/dx = g(x)·h(y), and our goal is to arrive at an explicit (or implicit) expression for y as a function of x that satisfies both the DE and the constraint y(x₀) = y₀.

STANDARD SEPARABLE FORM
dy/dx = g(x) · h(y)
The right-hand side factors into a function of x alone times a function of y alone. If the DE cannot be written this way, it is not separable.
SEPARATION STEP
(1 / h(y)) dy = g(x) dx
Divide both sides by h(y) and multiply both sides by dx. This is valid only where h(y) ≠ 0. Values of y where h(y) = 0 may yield equilibrium (constant) solutions that must be checked separately.
INTEGRATION
∫ (1 / h(y)) dy = ∫ g(x) dx + C
Integrate both sides independently. The arbitrary constant C encapsulates both constants of integration (we consolidate them into one). The result is the general solution, which may be implicit in y.
APPLYING THE INITIAL CONDITION
y(x₀) = y₀ → substitute to find C
Replace x with x₀ and y with y₀ in the general solution to solve for the specific value of C. Substituting this C back produces the particular solution.
⚠️ When to Check for Lost Solutions
Dividing by h(y) implicitly assumes h(y) ≠ 0. If h(y₀) = 0 for some constant y₀, then y(x) = y₀ for all x is itself a solution. Always verify whether the initial condition places you on one of these singular (equilibrium) solutions before proceeding with the standard separation algorithm.

Step-by-Step Procedure & Flowchart

To consolidate the mathematical framework into a concrete workflow, the following flowchart and accompanying table organize every decision point a solver encounters. Mastering this procedure turns what initially feels like an art—spotting the right algebraic moves—into a repeatable algorithm.

This flowchart summarizes the entire solving procedure. Note the branching at Step 2: if h(y₀) = 0 the initial condition sits on an equilibrium solution, and no further integration is needed.
Summary of the six-step procedure with common pitfalls
StepActionCommon Pitfall
1 — IdentifyFactor RHS as g(x)·h(y)Mistaking sums for products, e.g., x + y is not separable
2 — EquilibriaCheck if h(y₀) = 0Forgetting to check leads to dividing by zero later
3 — SeparateMove all y to one side, all x to the otherForgetting to move dy along with the y-terms
4 — IntegrateEvaluate ∫(1/h(y)) dy and ∫g(x) dxDropping the constant of integration or adding one to each side
5 — Solve for CSubstitute (x₀, y₀) into general solutionArithmetic errors; sign mistakes when exponentiating
6 — SimplifySolve for y explicitly if possibleNot every implicit solution can be solved for y; state domains carefully

Worked Example

Consider the initial value problem dy/dx = (x² + 1)·y, with y(0) = 3. We will solve this step by step, following the procedure outlined in Section 5.

Solve dy/dx = (x² + 1)y, y(0) = 3
1
Step 1 — Identify SeparabilityThe right-hand side is (x² + 1)·y. Here g(x) = x² + 1 and h(y) = y. Since f(x, y) factors as a product of a function of x alone and a function of y alone, the equation is separable.
2
Step 2 — Check Equilibriah(y) = y = 0 gives the equilibrium solution y = 0. Since our initial condition is y(0) = 3 ≠ 0, we are not on this equilibrium and may proceed with separation.
3
Step 3 — Separate VariablesDivide both sides by y (valid since y ≠ 0 near y = 3) and multiply by dx:
(1/y) dy = (x² + 1) dx
4
Step 4 — Integrate Both SidesIntegrate the left side with respect to y and the right side with respect to x: ∫(1/y) dy = ∫(x² + 1) dx → ln|y| = x³/3 + x + C. Since y(0) = 3 > 0 and the exponential function is always positive, we can drop the absolute value near the initial condition and write ln(y) = x³/3 + x + C.
ln(y) = x³/3 + x + C
5
Step 5 — Apply the Initial ConditionSubstitute x = 0 and y = 3: ln(3) = 0 + 0 + C, so C = ln(3).
C = ln 3
6
Step 6 — Write the Particular SolutionSubstituting C = ln 3 back: ln(y) = x³/3 + x + ln 3. Exponentiating both sides: y = ex³/3 + x + ln 3 = eln 3 · ex³/3 + x = 3ex³/3 + x.
y = 3e^(x³/3 + x)
Verification Tip
You can always verify your particular solution by (a) checking that y(0) = 3e⁰ = 3, confirming the initial condition, and (b) computing dy/dx = 3(x² + 1)ex³/3 + x = (x² + 1)·y, confirming it satisfies the DE.

Strengths, Limitations, and Comparisons

Separation of variables is the most elementary analytical technique for first-order ODEs, but it is far from universal. Understanding where the method excels and where it breaks down helps you choose the right approach for a given problem and appreciate why more advanced methods (integrating factors, exact equations, variation of parameters) exist.

Separable vs. non-separable first-order DEs
AspectSeparable DEsNon-separable First-Order DEs
ApplicabilityOnly when dy/dx = g(x)·h(y)Linear, exact, Bernoulli, or numerical methods needed
Conceptual difficultyLow — relies on basic integrationModerate to high — special techniques required
Closed-form solutionsOften obtainable; may be implicitSometimes; many require numerical solution
Risk of lost solutionsYes — dividing by h(y) can lose equilibriaDepends on method; integrating factors avoid this
Physical modelingExponential growth/decay, logistic models, simple mixingForced oscillations, coupled systems, nonlinear dynamics
KEY TAKEAWAY
Separation of variables is like having a universal socket wrench that fits a specific bolt size perfectly: when the bolt (DE) matches, the tool is fast and reliable. But if the bolt is a different size (the DE isn't separable), forcing the wrench will strip the thread. Recognizing separability upfront saves time and prevents algebraic dead ends—and knowing the limitations of the method guides you toward the right alternative tool.

Connection to Advanced Theory

The separable-equation technique is not an isolated trick; it connects deeply to several pillars of differential equations theory that you will encounter in subsequent courses. Below we outline how the ideas you have learned here extend into more advanced territory.

How separable-DE concepts generalize to advanced topics
This LessonAdvanced Extension
Separation of variables in ODESeparation of variables in PDEs (e.g., heat equation u_t = k·u_xx assumes u(x,t) = X(x)·T(t))
Single arbitrary constant Cn-parameter families in nth-order ODEs; boundary-value problems replace initial conditions
Existence & uniqueness via ICPicard–Lindelöf theorem guarantees unique solutions under Lipschitz continuity of f(x, y)
Equilibrium solutions h(y) = 0Stability analysis: linearization, phase portraits, Lyapunov functions classify equilibria as stable, unstable, or saddle
Explicit antiderivatives requiredWhen antiderivatives don't exist in closed form, numerical methods (Euler, Runge–Kutta) approximate solutions

Perhaps the most elegant extension is into partial differential equations, where the same principle—assume the solution factors into functions of individual variables—reduces a PDE in several variables to a system of ODEs, each of which may itself be separable. The Fourier series solutions to the heat and wave equations that you will study in later courses rely precisely on this maneuver. Mastering separation of variables now thus provides the conceptual scaffolding for some of the deepest techniques in applied mathematics.

Practice Problems

PROBLEM 1CONCEPTUAL
Explain why the differential equation dy/dx = x + y is not separable, while dy/dx = xy is separable. What structural property distinguishes the two?
PROBLEM 2BASIC CALCULATION
Solve the initial value problem dy/dx = 3x²y, y(0) = 5.
PROBLEM 3INTERMEDIATE
Solve dy/dx = (1 + y²) sin x, y(0) = 1. Express your answer as y as an explicit function of x.
PROBLEM 4APPLIED
A radioactive substance decays at a rate proportional to the amount present: dN/dt = −λN, where λ > 0 is the decay constant. If 100 grams are present initially and 80 grams remain after 5 years, find the particular solution N(t) and determine the half-life.
PROBLEM 5CRITICAL THINKING
Consider dy/dx = y(1 − y) with y(0) = y₀. (a) Identify all equilibrium solutions. (b) Use partial fractions to solve the IVP for 0 < y₀ < 1 and show that the particular solution approaches 1 as x → ∞. (c) Discuss what happens when y₀ = 0 or y₀ = 1 and explain why uniqueness of solutions prevents trajectories from crossing equilibria.

Lesson Summary

A first-order ODE dy/dx = f(x, y) is separable when f(x, y) factors as g(x)·h(y). The solution method consists of dividing by h(y), multiplying by dx to place all y-terms with dy and all x-terms with dx, and then integrating both sides independently. The result is a general solution containing an arbitrary constant C that represents a family of curves. An initial condition y(x₀) = y₀ selects the unique particular solution from this family by determining C.

Key cautions include checking for equilibrium solutions (values of y where h(y) = 0), which can be lost during the division step, and ensuring that the solution's domain is stated correctly when implicit forms arise. The separable technique generalizes powerfully into partial differential equations and underpins the Fourier-series approach to the heat and wave equations. Mastery of this method—particularly the discipline of applying initial conditions systematically—provides an essential foundation for all subsequent work in differential equations.

Varsity Tutors • Calculus 2 • Solving Separable DEs — Solving Separable Differential Equations with Initial Conditions