ODE: Substitution Methods

Solving differential equations often requires clever transformations to convert a complex, unsolvable form into one we can integrate. For engineers, mastering these substitution methods is non-negotiable; they are the key to modeling everything from heat transfer and fluid dynamics to circuit analysis and population growth. This guide will equip you with a systematic toolkit for identifying and applying the correct substitution to transform an ordinary differential equation (ODE) into a separable or linear equation, which you can then solve with confidence.

The Core Idea: Change of Variable

At its heart, a substitution is a strategic change of variable. You replace the dependent variable $y$ (or sometimes the independent variable $x$) with a new function, say $v$ or $w$, which is defined in terms of $y$ and $x$. The goal is to rewrite the original ODE entirely in terms of this new variable and its derivative, hoping the new equation is of a simpler, solvable type. Success hinges on two skills: correctly executing the derivative chain rule during the substitution and, more importantly, recognizing which substitution to use based on the equation's structure.

Homogeneous Equations: The $v = y/x$ Substitution

A first-order ODE is called homogeneous if it can be written in the form $dy/dx = F(y/x)$. The right-hand side is a function of the combined ratio $y/x$ alone. This structure invites the substitution $v = y/x$, which implies $y = vx$.

To execute this, we differentiate the new relation using the product rule: $dy/dx = v + x (dv/dx)$. We then substitute both expressions into the original ODE:

1. Replace $y/x$ with $v$.
2. Replace $dy/dx$ with $v + x (dv/dx)$.

This yields an equation of the form $v + x \frac{dv}{dx} = F(v)$. Crucially, this new ODE is separable. We can rearrange it as $\frac{dv}{F(v) - v} = \frac{dx}{x}$, which can be solved by direct integration. After integrating and solving for $v$, you must remember to back-substitute $v = y/x$ to express the final solution in terms of the original variables $y$ and $x$.

Example: Solve $x \frac{dy}{dx} = y + \sqrt{x^2 - y^2}$ for $x > 0$. First, rewrite in standard form: $\frac{dy}{dx} = \frac{y}{x} + \sqrt{1 - (\frac{y}{x})^2}$. This is $F(y/x)$. Now, let $v = y/x$, so $y = vx$ and $dy/dx = v + x dv/dx$. Substitute: $v + x \frac{dv}{dx} = v + \sqrt{1 - v^2}$. The $v$ terms cancel, giving $x \frac{dv}{dx} = \sqrt{1 - v^2}$. This is separable: $\int \frac{dv}{\sqrt{1 - v^2}} = \int \frac{dx}{x}$. Integrating yields $\arcsin(v) = \ln|x| + C$. Therefore, $v = \sin(\ln|x| + C)$. Finally, back-substitute: $y = x \sin(\ln|x| + C)$.

The Bernoulli Equation: Linearizing via $w = y^{1-n}$

The Bernoulli equation has the standard form $\frac{dy}{dx} + P(x)y = Q(x)y^n$, where $n$ is any real number except 0 or 1. (If $n=0$, it's linear; if $n=1$, it's separable). This nonlinear equation is transformed into a linear ODE through a specific substitution.

The substitution is $w = y^{1-n}$. We then find $dw/dx$ using the chain rule: $\frac{dw}{dx} = (1-n)y^{-n} \frac{dy}{dx}$. To proceed, multiply the original Bernoulli equation by $(1-n)y^{-n}$. This strategic multiplication gives: $$(1-n)y^{-n}\frac{dy}{dx} + (1-n)y^{-n}P(x)y = (1-n)y^{-n}Q(x)y^n$$ Simplify each term:

• Term 1: $(1-n)y^{-n}\frac{dy}{dx}$ is exactly $\frac{dw}{dx}$.
• Term 2: $(1-n)y^{-n}P(x)y = (1-n)y^{1-n}P(x) = (1-n)w P(x)$.
• Term 3: $(1-n)y^{-n}Q(x)y^n = (1-n)Q(x)$.

The equation is now linear in $w$: $\frac{dw}{dx} + (1-n)P(x)w = (1-n)Q(x)$. You solve this using an integrating factor, then back-substitute $w = y^{1-n}$ to find $y$.

Example: Solve $\frac{dy}{dx} - \frac{4}{x}y = x \sqrt{y}$, where $y > 0$. This is a Bernoulli equation with $n = 1/2$. Our substitution is $w = y^{1 - 1/2} = y^{1/2}$. Then, $\frac{dw}{dx} = \frac{1}{2}y^{-1/2} \frac{dy}{dx}$. Multiply the original ODE by $\frac{1}{2}y^{-1/2}$: $\frac{1}{2}y^{-1/2}\frac{dy}{dx} - \frac{1}{2}y^{-1/2} \cdot \frac{4}{x}y = \frac{1}{2}y^{-1/2} \cdot x \sqrt{y}$. Simplify: $\frac{dw}{dx} - \frac{2}{x} w = \frac{x}{2}$. This is linear in $w$. The integrating factor is $\mu(x) = e^{\int (-2/x) dx} = e^{-2\ln|x|} = x^{-2}$. Multiply through: $x^{-2}\frac{dw}{dx} - 2x^{-3}w = \frac{1}{2x}$. The left side is the derivative of $x^{-2}w$. Integrate: $\frac{d}{dx}(x^{-2}w) = \frac{1}{2x} \implies x^{-2}w = \frac{1}{2}\ln|x| + C$. Thus, $w = \frac{x^2}{2}\ln|x| + Cx^2$. Back-substitute $w = \sqrt{y}$: $\sqrt{y} = \frac{x^2}{2}\ln|x| + Cx^2$, so $y = \left(\frac{x^2}{2}\ln|x| + Cx^2\right)^2$.

Recognizing the Correct Substitution

Choosing the right tool is 90% of the battle. Your analysis should be a quick structural check of the ODE $dy/dx = f(x, y)$.

1. Check for Homogeneous Form: Can you manipulate the equation so all terms are functions of the ratio $y/x$? If yes, use $v = y/x$.
2. Check for Bernoulli Form: Is the equation of the form $y' + P(x)y = Q(x)y^n$? If yes, use $w = y^{1-n}$.
3. Check for "Almost Separable" Forms: Sometimes an equation like $dy/dx = f(ax+by+c)$ can be made separable with the substitution $v = ax+by+c$. Differentiating gives $dv/dx = a + b dy/dx$, allowing you to replace $dy/dx$ with $(1/b)(dv/dx - a)$.
4. General Heuristic: Look for composite arguments in the function $f(x,y)$. The substitution is often the "inner function" of that composite. If $f$ involves $(y/x)$, substitute $v=y/x$. If it involves $y^2/x$, you might try $v=y^2$, but always check if a standard form applies first.

Common Pitfalls

1. Misapplying the Homogeneous Substitution

• Pitfall: Attempting to use $v = y/x$ on an equation that is not homogeneous, leading to a messier, unsimplified equation.
• Correction: Always perform the test: An ODE is homogeneous if replacing $x$ with $tx$ and $y$ with $ty$ yields the original equation multiplied by $t^n$ (often $t^0$). In practice, rewrite the ODE as $dy/dx = F(y/x)$. If you cannot factor $x$ out to get this form, it's not homogeneous.

2. Forgetting the Chain Rule During Substitution

• Pitfall: Substituting $y = vx$ but then writing $dy/dx = dv/dx$, forgetting the product rule term. This is the most common algebraic error.
• Correction: Be meticulous. Write $y = vx$, then differentiate: $\frac{dy}{dx} = \frac{d}{dx}(vx) = x\frac{dv}{dx} + v$. This step is not optional.

3. Neglecting the Back-Substitution

• Pitfall: Solving the transformed ODE for the intermediate variable (like $v$ or $w$) and stopping, leaving the answer in the wrong variable.
• Correction: Your final answer must be an explicit or implicit relation between $y$ and $x$. The last step is always to replace the intermediate variable with its definition in terms of $y$ and $x$.

4. Incorrectly Handling the Bernoulli Exponent $n$

• Pitfall: For the Bernoulli equation $y' + P(x)y = Q(x)y^n$, mistakenly using the substitution $w = y^n$ instead of $w = y^{1-n}$.
• Correction: Memorize the standard transformation: $w = y^{1-n}$. Derive it once to understand why, then apply it consistently. The factor $(1-n)$ that appears is critical for the cancellation that yields a linear equation.

Summary

• Substitution methods are powerful techniques that convert otherwise unsolvable first-order ODEs into separable or linear forms by introducing a new dependent variable.
• For homogeneous equations of the form $dy/dx = F(y/x)$, the substitution $v = y/x$ (with $dy/dx = v + x dv/dx$) always leads to a separable equation in $v$ and $x$.
• The Bernoulli equation, $y' + P(x)y = Q(x)y^n$, is linearized via the substitution $w = y^{1-n}$, which, after multiplication by $(1-n)y^{-n}$, produces a linear ODE in $w(x)$ solvable with an integrating factor.
• Successful application depends on recognizing the equation's structure to choose the correct substitution. Your first step should always be to classify the ODE.
• Avoid common errors by carefully applying the chain rule during differentiation, methodically simplifying the substituted equation, and never forgetting the final back-substitution to return to the original variables.