Solution of First-Order Equations: Separable, Exact & Special Integrating Factors
Exhaustive treatment of first-order analytical integration: variables separable equations, substitutions reducing to separable forms, homogeneous equations and the Euler substitution v = y/x, exact differential forms, potential function construction, special single-variable and multi-variable integrating factors, and Leibniz first-order linear equations.
§2.1 Separable Differential Equations & Transformations of Dependent Variables
1. Variables Separable Equations
A first-order differential equation is separable if the derivative factorizes into a product of a function of $x$ alone and a function of $y$ alone:
Integrating both sides directly yields the one-parameter implicit general solution:
Caution: Any constant $y = c$ where $h(c) = 0$ is a potential singular or equilibrium solution that must be checked separately, as division by $h(y)$ assumes $h(y) \neq 0$.
2. Equations Reducible to Separable Form via Substitution $v = ax + by + c$
Equations of the form $\frac{dy}{dx} = f(ax + by + c)$, where $b \neq 0$, are rendered separable by defining the substitution $v = ax + by + c$. Differentiating with respect to $x$:
Integrating both sides yields the solution in terms of $v$, and back-substituting $v = ax + by + c$ gives the final answer.
3. Homogeneous Equations of Degree Zero
A function $f(x, y)$ is homogeneous of degree $n$ if $f(tx, ty) = t^n f(x, y)$ for all $t > 0$. An ODE $M(x, y)dx + N(x, y)dy = 0$ is homogeneous if $M$ and $N$ are homogeneous functions of the same degree $n$. In that case, $\frac{dy}{dx} = F\left(\frac{y}{x}\right)$.
The standard transformation is the Euler substitution: $y = v x$, whence $\frac{dy}{dx} = v + x \frac{dv}{dx}$. Substituting into the ODE:
which is completely separated in the variables $v$ and $x$.
§2.2 Exact Differential Equations & Potential Function Construction
1. The Total Differential & The Exactness Criterion
A first-order differential expression $M(x, y)dx + N(x, y)dy$ is an exact differential in an open simply connected domain $D \subset \mathbb{R}^2$ if there exists a continuously differentiable potential function $\Phi(x, y)$ such that:
If $M(x, y)dx + N(x, y)dy = 0$ is exact, its general solution is immediately given by the level curves:
The differential equation $M(x, y)dx + N(x, y)dy = 0$ is exact if and only if:
2. Construction of the Potential Function $\Phi(x, y)$
To find $\Phi(x, y)$, we integrate $\frac{\partial \Phi}{\partial x} = M(x, y)$ with respect to $x$, treating $y$ as a constant:
where $g(y)$ is an arbitrary function of $y$ acting as the integration constant. Differentiating this expression with respect to $y$ and equating to $N(x, y)$:
By the exactness condition $\frac{\partial M}{\partial y} = \frac{\partial N}{\partial x}$, the right-hand side is strictly independent of $x$. Integrating $g'(y)$ with respect to $y$ yields $g(y)$ and completes the potential function $\Phi(x, y) = C$.
§2.3 Special Integrating Factors: Single & Multi-Variable Formulas
1. Concept of the Integrating Factor
If the equation $M(x, y)dx + N(x, y)dy = 0$ is not exact ($\frac{\partial M}{\partial y} \neq \frac{\partial N}{\partial x}$), it can often be made exact by multiplying by an integrating factor $\mu(x, y) \neq 0$ such that:
2. Standard Single-Variable Integrating Factors
- Integrating factor depending only on $x$ ($\mu = \mu(x)$):
Here $\frac{\partial \mu}{\partial y} = 0$ and $\frac{\partial \mu}{\partial x} = \frac{d\mu}{dx}$. The condition reduces to: $$\frac{1}{\mu} \frac{d\mu}{dx} = \frac{M_y - N_x}{N}$$ If $\frac{M_y - N_x}{N} = f(x)$ is a function of $x$ alone, then: $$\mu(x) = \exp\left(\int \frac{M_y - N_x}{N} dx\right)$$ - Integrating factor depending only on $y$ ($\mu = \mu(y)$):
Here $\frac{\partial \mu}{\partial x} = 0$ and $\frac{\partial \mu}{\partial y} = \frac{d\mu}{dy}$. The condition reduces to: $$\frac{1}{\mu} \frac{d\mu}{dy} = \frac{N_x - M_y}{M}$$ If $\frac{N_x - M_y}{M} = g(y)$ is a function of $y$ alone, then: $$\mu(y) = \exp\left(\int \frac{N_x - M_y}{M} dy\right)$$
3. Multi-Variable Integrating Factors & Homogeneous Theorem
- If $\frac{M_y - N_x}{y N - x M} = h(xy)$, then $\mu = \exp\left(\int h(u) du\right)$ where $u = xy$.
- Homogeneous Equation Theorem: If $M dx + N dy = 0$ is a homogeneous equation of degree $n$ and $Mx + Ny \neq 0$, then: $$\mu(x, y) = \frac{1}{Mx + Ny}$$ is always an integrating factor.
§2.4 First-Order Linear Differential Equations (Leibniz Formula)
1. Standard Linear Form
The general first-order linear differential equation is written in canonical form as:
where $P(x)$ and $Q(x)$ are continuous functions on an interval $I$. Rewriting in differential form: $(P(x)y - Q(x))dx + dy = 0$, where $M = P(x)y - Q(x)$ and $N = 1$.
Testing for exactness: $M_y = P(x)$ and $N_x = 0$. Since $M_y \neq N_x$, the equation is not exact, but:
which depends strictly on $x$ alone!
2. Derivation of the Closed-Form General Solution
The integrating factor is $\mu(x) = \exp\left(\int P(x) dx\right)$. Multiplying the standard form by $\mu(x)$:
By the product rule of differentiation, the left side is the exact derivative of the product $\mu(x) y$:
Integrating both sides directly with respect to $x$:
This closed-form formula provides the complete general solution to any first-order linear ODE.
Step-by-Step Solved Examination Problems
Comprehensive analytical derivations, multi-tier solutions (Foundational, Intermediate Exam, and Honors/Proof Challenge) with complete line-by-line verification.
Test for exactness and find the general solution of the differential equation:
Step 1: Test exactness
Here $M(x, y) = 2x \cos y + 3x^2 y$ and $N(x, y) = x^3 - x^2 \sin y - y$.
Compute the partial derivatives:
Since $\frac{\partial M}{\partial y} = \frac{\partial N}{\partial x}$, the equation is exact in all of $\mathbb{R}^2$.
Step 2: Integrate $M$ with respect to $x$
Step 3: Differentiate with respect to $y$ and match $N$
Equating gives $g'(y) = -y \implies g(y) = -\frac{1}{2}y^2$.
Step 4: Formulate the general solution
The potential function is $\Phi(x, y) = x^2 \cos y + x^3 y - \frac{1}{2}y^2 = C$.
$\Phi(x, y) = x^2 \cos y + x^3 y - \frac{1}{2}y^2 = C$ (or $2x^2 \cos y + 2x^3 y - y^2 = C'$).
Solve the non-exact differential equation:
by deriving an appropriate integrating factor.
Step 1: Test exactness
$M = y + x y^2$ and $N = -x$.
Since $M_y \neq N_x$, the equation is not exact.
Step 2: Determine integrating factor
Check $\frac{N_x - M_y}{M}$:
This depends purely on $y$! The integrating factor is:
Step 3: Multiply ODE by $\mu(y)$
Now $M^ = \frac{1}{y} + x$ and $N^ = -\frac{x}{y^2}$. Note $M^_y = -1/y^2 = N^_x$ (exact!).
Step 4: Integrate
Thus $\frac{x}{y} + \frac{x^2}{2} = C$.
$\mu(y) = \frac{1}{y^2}$, yielding the general solution $\frac{x}{y} + \frac{x^2}{2} = C$ (or $2x + x^2 y = 2Cy$).
Solve the homogeneous initial value problem:
using both the substitution $y = vx$ and the homogeneous integrating factor theorem $\mu = 1/(Mx + Ny)$.
Method 1: Euler Substitution $y = vx$
With $y = vx \implies dy = v dx + x dv$:
Dividing by $x^2 \neq 0$:
Separating variables:
Back-substituting $v = y/x$:
Method 2: Homogeneous Integrating Factor
$M = x^2 + 3xy + y^2$, $N = -x^2$.
Therefore, $\mu = \frac{1}{x(x + y)^2}$. Multiplying the ODE by $\mu$ yields the identical potential function differential $d\left[\ln |x| + \frac{x}{x + y}\right] = 0$.
Initial Condition:
At $x = 1, y = 0$: $\frac{1}{1 + 0} + \ln(1) = 1 + 0 = C \implies C = 1$.
Thus $\frac{x}{x + y} + \ln |x| = 1 \implies x + y = \frac{x}{1 - \ln |x|}$.
General solution: $\frac{x}{x + y} + \ln |x| = C$. With $y(1) = 0$, $C = 1$, yielding explicit solution $y(x) = \frac{x \ln |x|}{1 - \ln |x|}$.