Bernoulli Equations - Differential Equations

SECTION 1

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.

SECTION 2

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.

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.

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.

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.

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.

SECTION 3

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 = y¹⁻ⁿ, 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.

SECTION 4

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.

SECTION 5

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.

Value of n	Substitution v = y¹⁻ⁿ	Character	Example Equation
n = 0	v = y (no change needed)	Already linear	dy/dx + 2y = x
n = 1	v = y⁰ = 1 (trivial)	Separable/linear	dy/dx + 2y = 3y → dy/dx = y
n = 2	v = y⁻¹ = 1/y	Most common case	dy/dx + y = xy²
n = 3	v = y⁻² = 1/y²	Common in physics	dy/dx − y = eˣy³
n = 1/2	v = y¹ᐟ² = √y	Fractional exponents	dy/dx + y = y¹ᐟ²

Common values of n in Bernoulli equations and their corresponding substitutions

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.

SECTION 6

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 n

Compare 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 factor

The 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 sides

Integrate: 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 y

Since 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.

SECTION 7

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.

Strengths	Limitations
Converts a nonlinear equation into a linear one, which always has a systematic solution	Only 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 problems	The 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 answer	Division by yⁿ may lose the solution y = 0, which should be checked separately

When to use — and watch out for — the Bernoulli method

✦

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.

SECTION 8

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.

Feature	Bernoulli Equations	Riccati Equations	Exact Equations
General form	dy/dx + Py = Qyⁿ	dy/dx = P + Qy + Ry²	M(x,y)dx + N(x,y)dy = 0
Linearity	Nonlinear (reducible to linear)	Nonlinear (reducible to Bernoulli)	May be linear or nonlinear
Key technique	Substitution v = y¹⁻ⁿ	Requires one known particular solution	Check ∂M/∂y = ∂N/∂x
Difficulty	Moderate	Hard	Moderate

Bernoulli equations compared with other first-order ODE types

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.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

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?

Show answer

Comparing with dy/dx + P(x)·y = Q(x)·yⁿ, we get P(x) = 3, Q(x) = 5, and n = 4. The substitution is v = y¹⁻ⁿ = y¹⁻⁴ = y⁻³ = 1/y³. This means v represents the reciprocal of the cube of y. After substitution, the resulting linear equation in v is dv/dx + (1 − 4)(3)v = (1 − 4)(5), which simplifies to dv/dx − 9v = −15.

Question 2: BASIC CALCULATION

Solve the Bernoulli equation dy/dx − y = −y². Write your answer as y in terms of x and an arbitrary constant C.

Show answer

Here P(x) = −1, Q(x) = −1, and n = 2. Divide by y²: y⁻²(dy/dx) − y⁻¹ = −1. Substitute v = y⁻¹, so dv/dx = −y⁻²(dy/dx), giving −dv/dx − v = −1, or dv/dx + v = 1. The integrating factor is μ = eˣ. Multiplying: d/dx[eˣv] = eˣ. Integrating: eˣv = eˣ + C, so v = 1 + Ce⁻ˣ. Since v = 1/y, the solution is y = 1/(1 + Ce⁻ˣ). Note that y = 0 is also a singular solution.

Question 3: INTERMEDIATE

Solve the initial value problem: dy/dx + (2/x)y = x²y³, with y(1) = 1. Find the particular solution.

Show answer

With P(x) = 2/x, Q(x) = x², n = 3, substitute v = y¹⁻³ = y⁻². Then dv/dx = −2y⁻³(dy/dx). Dividing the original equation by y³ gives y⁻³(dy/dx) + (2/x)y⁻² = x², so (−1/2)(dv/dx) + (2/x)v = x², giving dv/dx − (4/x)v = −2x². The integrating factor is μ = e^(∫−4/x dx) = x⁻⁴. Multiplying: d/dx[x⁻⁴v] = −2x⁻². Integrating: x⁻⁴v = 2x⁻¹ + C, so v = 2x³ + Cx⁴. Since v = y⁻² = 1/y², we have y² = 1/(2x³ + Cx⁴). Applying y(1) = 1: 1 = 1/(2 + C), so C = −1. The particular solution is y² = 1/(2x³ − x⁴), or y = 1/√(2x³ − x⁴).

Question 4: APPLIED

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?

Show answer

Rewrite as dP/dt − 0.5P = −0.001P². This is a Bernoulli equation with P(t) = −0.5 (the coefficient function), Q(t) = −0.001, and n = 2. The substitution is v = P¹⁻² = P⁻¹ = 1/P. After substitution you get a linear equation in v. Solving gives the logistic function P(t) = 500/(1 + Ce⁻⁰·⁵ᵗ). As t → ∞, the exponential term vanishes and P → 500. This means the population approaches a carrying capacity of 500 regardless of the initial population size, which is the hallmark of logistic growth.

Question 5: CRITICAL 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.

Show answer

The key property is that yⁿ is a power function, so when you divide by yⁿ, you get y¹⁻ⁿ, and the chain rule applied to v = y¹⁻ⁿ produces dv/dx = (1−n)y⁻ⁿ(dy/dx). The factor y⁻ⁿ(dy/dx) is exactly what appears after dividing the original equation by yⁿ. With sin(y), dividing by sin(y) produces y/sin(y) and (1/sin(y))·(dy/dx), which cannot be expressed as the derivative of any simple function of y. The substitution works because power functions are closed under differentiation via the chain rule. For n = 2, v = y⁻¹ = 1/y and dv/dx = −y⁻²(dy/dx), giving a clean negative sign. The linear equation becomes dv/dx − P(x)v = −Q(x), which is elegant and straightforward to solve with an integrating factor.

SUMMARY

Bernoulli Equations — Summary

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.