DIFFERENTIAL EQUATIONS • FIRST-ORDER DIFFERENTIAL EQUATIONS

Bernoulli Equations

A clever substitution transforms nonlinear differential equations into ones you already know how to solve.

Historical Context & Motivation

In the late 1600s and early 1700s, mathematicians were racing to solve differential equations — equations that describe how quantities change over time. Linear first-order equations were well understood, but many real-world problems produced equations with nonlinear terms, like a variable raised to a power. These equations seemed much harder to crack. The Bernoulli equation emerged from this challenge, offering a brilliant trick: a substitution that converts a nonlinear equation into a familiar linear one.

The Bernoulli family of Switzerland produced several outstanding mathematicians across two generations. Jacob Bernoulli first encountered these types of equations, and his younger brother Johann refined the solution method. Their work laid important groundwork for what would eventually become the systematic study of differential equations.

1687
Newton's Principia Published
Isaac Newton published his laws of motion and calculus-based physics, creating enormous demand for methods to solve differential equations that described physical phenomena.
1695
Jacob Bernoulli Poses the Problem
Jacob Bernoulli proposed a class of nonlinear first-order equations of the form dy/dx + P(x)y = Q(x)yⁿ and challenged the mathematical community to find solutions.
1696
Leibniz and Johann Bernoulli Solve It
Gottfried Wilhelm Leibniz and Johann Bernoulli independently discovered the substitution v = y¹⁻ⁿ that reduces the equation to a linear one, providing a general solution technique.
1700s
Applications Expand
Bernoulli equations found applications in population modeling, fluid dynamics, and mechanical systems, becoming a standard tool in the differential equations toolkit.

The central question that Bernoulli equations address is this: when a differential equation has a nonlinear term involving yⁿ, is there a systematic way to make it linear? The answer, as we'll see, is yes — through a clever change of variable that transforms the equation into something you already have the tools to solve.

Core Principles & Definitions

Before diving into the solution method, you need to understand what makes a Bernoulli equation special and how it differs from the linear equations you may have already studied. A Bernoulli equation is a first-order ordinary differential equation that looks almost linear, except for a power of y on the right side. It has the standard form:

BERNOULLI STANDARD FORM
dy/dx + P(x) · y = Q(x) · yⁿ
where P(x) and Q(x) are known functions of x, and n is a real number. The equation is nonlinear when n ≠ 0 and n ≠ 1.
1

Recognizing the Form

A Bernoulli equation has a yⁿ term on one side. If n = 0, it's just a linear equation. If n = 1, the yⁿ becomes y, and you can combine it with the P(x)y term — still linear. The interesting cases are n ≠ 0, 1.
2

The Key Substitution

The substitution v = y¹⁻ⁿ transforms the nonlinear Bernoulli equation into a first-order linear equation in v. This is the core idea behind the entire method.
3

Solve the Linear Equation

After substitution, you get a standard linear first-order equation dv/dx + (1 − n)P(x)v = (1 − n)Q(x). Solve this using an integrating factor, then back-substitute to find y.
4

Back-Substitute for y

Once you find v(x), remember that v = y¹⁻ⁿ. So y = v^(1/(1−n)). This gives you the solution to the original nonlinear equation.
KEY TAKEAWAY
Think of a Bernoulli equation like a locked door. On its own, it looks too complex to open. But the substitution v = y¹⁻ⁿ acts like a key that converts it into a linear equation — a door you already know how to open. Once you walk through (solve for v), you just translate back to find y.

Visual Explanation

The following diagram shows the complete workflow for solving a Bernoulli equation, from identifying the form to arriving at the final solution. Each step transforms the equation into something more manageable.

The flowchart shows the four main steps: identify the Bernoulli equation, perform the substitution v = y1−n, solve the resulting linear equation for v, and convert back to y. The color-coded arrows guide you through each transformation.

Notice how the diagram moves from a nonlinear equation at the top to a linear equation in the middle. That transformation is the heart of the Bernoulli method. The substitution step (shown in cyan) is where the magic happens — it absorbs the troublesome yⁿ term and produces a clean, standard linear equation that you can solve with methods you already know.

Mathematical Framework

Let's walk through the algebra that makes the Bernoulli substitution work. We start with the standard form and systematically derive the linear equation that results from the change of variable.

Starting Point: The Standard Form

STANDARD FORM
dy/dx + P(x) · y = Q(x) · yⁿ
We require n ≠ 0 and n ≠ 1, since those cases are already linear.

The first manipulation is to divide both sides of the equation by yⁿ. This isolates the nonlinear part and sets up the substitution. After dividing:

DIVIDE BY yⁿ
y⁻ⁿ · (dy/dx) + P(x) · y¹⁻ⁿ = Q(x)
Notice that the right side is now just Q(x) — no y at all! And the left side has y¹⁻ⁿ, which suggests our substitution.

The Substitution

BERNOULLI SUBSTITUTION
v = y¹⁻ⁿ → dv/dx = (1 − n) · y⁻ⁿ · (dy/dx)
We define v as y raised to the power (1 − n). Using the chain rule, its derivative has a factor of y⁻ⁿ · (dy/dx), which is exactly the combination that appears in our divided equation.

Now we can express y⁻ⁿ · (dy/dx) as (1/(1 − n)) · (dv/dx). Substituting this into our divided equation and also replacing y¹⁻ⁿ with v gives:

RESULTING LINEAR EQUATION
dv/dx + (1 − n) · P(x) · v = (1 − n) · Q(x)
This is a standard first-order linear equation in v, solvable using the integrating factor μ(x) = e^(∫(1−n)P(x) dx).
💡 Remember the Integrating Factor
To solve dv/dx + R(x)·v = S(x), multiply both sides by the integrating factor μ(x) = e^(∫R(x) dx). The left side collapses to d/dx[μ·v], so integrating both sides isolates v. This step draws directly on the integrating-factor method for linear first-order equations, so it's worth reviewing that technique before working through the example below.

Special Cases & Classifications

Not all Bernoulli equations are equally challenging. The value of n determines the difficulty and the nature of the substitution. Let's examine the most common cases you'll encounter and how they differ.

Common values of n in Bernoulli equations and their corresponding substitutions
Value of nSubstitution v = y¹⁻ⁿCharacterExample Equation
n = 0v = y (no change needed)Already lineardy/dx + 2y = x
n = 1v = y⁰ = 1 (trivial)Separable/lineardy/dx + 2y = 3y → dy/dx = y
n = 2v = y⁻¹ = 1/yMost common casedy/dx + y = xy²
n = 3v = y⁻² = 1/y²Common in physicsdy/dx − y = eˣy³
n = 1/2v = y¹ᐟ² = √yFractional exponentsdy/dx + y = y¹ᐟ²
Solution curves for linear equations (blue) grow steadily, while Bernoulli equation solutions (pink, with n = 2) can grow more dramatically or approach vertical asymptotes. The dashed lines show how different initial conditions produce distinct solution curves.

The graph above illustrates why Bernoulli equations are fundamentally different from linear ones. The linear solution curves (in blue) show smooth, predictable growth. In contrast, the Bernoulli solution curves (in pink) can exhibit much steeper growth — and in some cases may even blow up to infinity in finite time. This kind of nonlinear behavior is precisely why Bernoulli equations require their own solution technique.

Worked Example

Let's solve a complete Bernoulli equation step by step. We'll work through every algebraic detail so you can see exactly how the method unfolds.

Solve: dy/dx + y = xy²
1
Step 1 — Identify P(x), Q(x), and nCompare with the standard form dy/dx + P(x)·y = Q(x)·yⁿ. Here, P(x) = 1, Q(x) = x, and n = 2. Since n ≠ 0 and n ≠ 1, this is indeed a Bernoulli equation that requires substitution.
P(x) = 1, Q(x) = x, n = 2
2
Step 2 — Divide both sides by yⁿ = y²Dividing every term by y² gives: y⁻² · (dy/dx) + y⁻¹ = x. Notice that y⁻¹ = y¹⁻² = y¹⁻ⁿ. This is exactly the form we need for our substitution.
y⁻²(dy/dx) + y⁻¹ = x
3
Step 3 — Make the substitution v = y¹⁻ⁿ = y⁻¹Let v = y⁻¹ = 1/y. Then dv/dx = −y⁻² · (dy/dx), which means y⁻² · (dy/dx) = −dv/dx. Substituting into our divided equation: −dv/dx + v = x, or equivalently dv/dx − v = −x.
dv/dx − v = −x (a linear equation!)
4
Step 4 — Find the integrating factorThe linear equation dv/dx − v = −x has R(x) = −1. The integrating factor is μ(x) = e^(∫(−1)dx) = e⁻ˣ. Multiply both sides by e⁻ˣ: e⁻ˣ · (dv/dx) − e⁻ˣ · v = −x · e⁻ˣ. The left side is d/dx[e⁻ˣ · v].
d/dx[e⁻ˣ · v] = −x · e⁻ˣ
5
Step 5 — Integrate both sidesIntegrate: e⁻ˣ · v = ∫(−x · e⁻ˣ) dx. Using integration by parts (with u = −x and dw = e⁻ˣ dx, so du = −dx and w = −e⁻ˣ), we get ∫(−x · e⁻ˣ) dx = xe⁻ˣ + e⁻ˣ + C = (x + 1)e⁻ˣ + C. So e⁻ˣ · v = (x + 1)e⁻ˣ + C.
v = (x + 1) + Ceˣ
6
Step 6 — Back-substitute for ySince v = 1/y, we have y = 1/v. Therefore the solution is:
y = 1 / (x + 1 + Ceˣ)
Check Your Answer
You can verify the solution by substituting y = 1/(x + 1 + Ceˣ) back into the original equation dy/dx + y = xy². Computing dy/dx using the quotient rule and confirming both sides are equal is a good habit that catches algebra mistakes.

Strengths & Limitations

The Bernoulli equation method is powerful, but like any mathematical technique, it has both advantages and constraints. Understanding these helps you know when to reach for this tool and when a different approach might be better.

When to use — and watch out for — the Bernoulli method
StrengthsLimitations
Converts a nonlinear equation into a linear one, which always has a systematic solutionOnly works for equations in the specific form dy/dx + P(x)y = Q(x)yⁿ — not all nonlinear equations qualify
Works for any real value of n (integers, fractions, negatives)Requires knowledge of the integrating factor method as a prerequisite
Produces a general solution with one arbitrary constant, suitable for initial value problemsThe resulting integrals may be difficult or impossible to evaluate in closed form
A well-defined, algorithmic process — follow the steps and you'll arrive at the answerDivision by yⁿ may lose the solution y = 0, which should be checked separately
KEY TAKEAWAY
The Bernoulli method is like a universal adapter for a specific kind of nonlinear equation. Just as an adapter lets you plug a foreign device into a local outlet, the substitution v = y¹⁻ⁿ lets you plug a nonlinear equation into the linear solution method. But it only fits equations shaped the right way — if the 'plug' doesn't match the Bernoulli form, you'll need a different adapter.
Don't Forget y = 0
When you divide by yⁿ in Step 2, you're assuming y ≠ 0. But y = 0 might actually be a valid solution to the original equation. Always check whether y = 0 satisfies the original Bernoulli equation. If it does, it's called a singular solution that you could miss if you only follow the substitution method.

Connection to Advanced Methods

The Bernoulli equation is your first introduction to solving nonlinear differential equations by reducing them to linear ones. This strategy — transforming a hard problem into a simpler one — appears over and over in more advanced mathematics. Understanding where Bernoulli equations sit in the broader landscape of differential equations helps you see the bigger picture.

Bernoulli equations compared with other first-order ODE types
FeatureBernoulli EquationsRiccati EquationsExact Equations
General formdy/dx + Py = Qyⁿdy/dx = P + Qy + Ry²M(x,y)dx + N(x,y)dy = 0
LinearityNonlinear (reducible to linear)Nonlinear (reducible to Bernoulli)May be linear or nonlinear
Key techniqueSubstitution v = y¹⁻ⁿRequires one known particular solutionCheck ∂M/∂y = ∂N/∂x
DifficultyModerateHardModerate

An interesting connection to note: the Riccati equation is a generalization of the Bernoulli equation. If you know one particular solution of a Riccati equation, you can reduce it to a Bernoulli equation, which you then reduce to a linear equation. It's like a chain of transformations, each step making the problem simpler. In advanced courses, you'll also encounter the logistic equation dy/dx = ky(1 − y/M), which is actually a Bernoulli equation with n = 2. This equation models population growth with a carrying capacity, showing that Bernoulli equations have direct real-world applications.

Practice Problems

PROBLEM 1CONCEPTUAL
Consider the equation dy/dx + 3y = 5y⁴. Identify P(x), Q(x), and n. What substitution would you use to convert this into a linear equation? What would the new variable v represent in terms of y?
PROBLEM 2BASIC CALCULATION
Solve the Bernoulli equation dy/dx − y = −y². Write your answer as y in terms of x and an arbitrary constant C.
PROBLEM 3INTERMEDIATE
Solve the initial value problem: dy/dx + (2/x)y = x²y³, with y(1) = 1. Find the particular solution.
PROBLEM 4APPLIED
A population P(t) is modeled by the logistic equation dP/dt = 0.5P − 0.001P², where t is in years. Rewrite this as a Bernoulli equation in standard form, identify n, and describe the substitution you would use. What does the long-term behavior (as t → ∞) of the solution predict about the population?
PROBLEM 5CRITICAL THINKING
Explain why the substitution v = y¹⁻ⁿ works specifically for Bernoulli equations but not for a general nonlinear equation like dy/dx + P(x)y = Q(x)·sin(y). What property of the yⁿ term makes the Bernoulli trick possible? Additionally, show that when n = 2, the Bernoulli substitution and the resulting linear equation have a particularly clean form.

Bernoulli Equations — Summary

A Bernoulli equation has the form dy/dx + P(x)·y = Q(x)·yⁿ, where n ≠ 0 and n ≠ 1. The key substitution v = y¹⁻ⁿ transforms this nonlinear equation into the linear equation dv/dx + (1 − n)P(x)v = (1 − n)Q(x), which can be solved using an integrating factor. After finding v(x), you back-substitute using y = v^(1/(1−n)) to obtain the solution to the original equation.

The method follows four clear steps: identify P(x), Q(x), and n; divide by yⁿ; substitute v = y¹⁻ⁿ to get a linear ODE; and solve and back-substitute. Always check whether y = 0 is a singular solution that may be lost during the division step. Common applications include the logistic equation for population modeling and various problems in physics and engineering.

Varsity Tutors • Differential Equations • Bernoulli Equations