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+xdv/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)wP(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}⟹x^{-2}w=\frac{1}{2}\ln |x|+C$. Thus, $w=\frac{x^2}{2}\ln |x|+C{x}^2$. Back-substitute $w=\sqrt{y}$: $\sqrt{y}=\frac{x^2}{2}\ln |x|+C{x}^2$, so $y={\left(\frac{x^2}{2}\ln |x|+C{x}^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^{\prime }+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+bdy/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^{\prime }+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+xdv/dx$) always leads to a separable equation in $v$ and $x$.
• The Bernoulli equation, $y^{\prime }+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.