Mathematics / Pure Mathematics Differential Equations I 100% Free Open Access
Chapter 2 • Theory & Derivations

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:

$$\frac{dy}{dx} = g(x) h(y) \implies \frac{1}{h(y)} dy = g(x) dx \quad (h(y) \neq 0)$$

Integrating both sides directly yields the one-parameter implicit general solution:

$$\int \frac{dy}{h(y)} = \int g(x) dx + C$$

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$:

$$\frac{dv}{dx} = a + b \frac{dy}{dx} = a + b f(v) \implies \frac{dv}{a + b f(v)} = dx$$

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:

$$v + x \frac{dv}{dx} = F(v) \implies x \frac{dv}{dx} = F(v) - v \implies \frac{dv}{F(v) - v} = \frac{dx}{x}$$

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:

$$d\Phi = \frac{\partial \Phi}{\partial x} dx + \frac{\partial \Phi}{\partial y} dy = M(x, y)dx + N(x, y)dy$$

If $M(x, y)dx + N(x, y)dy = 0$ is exact, its general solution is immediately given by the level curves:

$$\Phi(x, y) = C$$
$$\mathbf{\text{Theorem (Criterion for Exactness): Let } M(x, y) \text{ and } N(x, y) \text{ be continuous with continuous first partial derivatives on } D.}}$$

The differential equation $M(x, y)dx + N(x, y)dy = 0$ is exact if and only if:

$$\frac{\partial M}{\partial y} = \frac{\partial N}{\partial x} \quad \forall (x, y) \in D$$

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:

$$\Phi(x, y) = \int M(x, y) dx + g(y)$$

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)$:

$$\frac{\partial \Phi}{\partial y} = \frac{\partial}{\partial y}\left(\int M(x, y) dx\right) + g'(y) = N(x, y) \implies g'(y) = N(x, y) - \frac{\partial}{\partial y}\left(\int M(x, y) dx\right)$$

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:

$$\frac{\partial}{\partial y}(\mu M) = \frac{\partial}{\partial x}(\mu N) \implies \mu \frac{\partial M}{\partial y} + M \frac{\partial \mu}{\partial y} = \mu \frac{\partial N}{\partial x} + N \frac{\partial \mu}{\partial x}$$

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:

$$\frac{dy}{dx} + P(x) y = Q(x)$$

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:

$$\frac{M_y - N_x}{N} = \frac{P(x) - 0}{1} = P(x)$$

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)$:

$$e^{\int P(x)dx} \frac{dy}{dx} + P(x) e^{\int P(x)dx} y = Q(x) e^{\int P(x)dx}$$

By the product rule of differentiation, the left side is the exact derivative of the product $\mu(x) y$:

$$\frac{d}{dx}\left[ y \cdot e^{\int P(x)dx} \right] = Q(x) e^{\int P(x)dx}$$

Integrating both sides directly with respect to $x$:

$$y \cdot e^{\int P(x)dx} = \int Q(x) e^{\int P(x)dx} dx + C \implies y(x) = e^{-\int P(x)dx} \left[ \int Q(x) e^{\int P(x)dx} dx + C \right]$$

This closed-form formula provides the complete general solution to any first-order linear ODE.

TIERED UNIVERSITY HONORS PROBLEMS

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.

Tier 1 • Foundational Example 2.1: Exact Differential Equation with Trigonometric Potential

Test for exactness and find the general solution of the differential equation:

$$(2x \cos y + 3x^2 y) dx + (x^3 - x^2 \sin y - y) dy = 0$$

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:

$$\frac{\partial M}{\partial y} = -2x \sin y + 3x^2$$


$$\frac{\partial N}{\partial x} = 3x^2 - 2x \sin y$$


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$

$$\Phi(x, y) = \int (2x \cos y + 3x^2 y) dx + g(y) = x^2 \cos y + x^3 y + g(y)$$



Step 3: Differentiate with respect to $y$ and match $N$

$$\frac{\partial \Phi}{\partial y} = -x^2 \sin y + x^3 + g'(y) = N(x, y) = x^3 - x^2 \sin y - y$$

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

Final Answer & Physical Insight

$\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'$).

Tier 2 • Analytical University Exam Example 2.2: Special Integrating Factor $\mu(y)$ for a Non-Exact Equation

Solve the non-exact differential equation:

$$(y + x y^2) dx - x dy = 0$$


by deriving an appropriate integrating factor.

Step 1: Test exactness
$M = y + x y^2$ and $N = -x$.

$$\frac{\partial M}{\partial y} = 1 + 2xy, \quad \frac{\partial N}{\partial x} = -1 \implies M_y - N_x = 2 + 2xy = 2(1 + xy)$$

Since $M_y \neq N_x$, the equation is not exact.

Step 2: Determine integrating factor
Check $\frac{N_x - M_y}{M}$:

$$\frac{N_x - M_y}{M} = \frac{-1 - (1 + 2xy)}{y(1 + xy)} = \frac{-2(1 + xy)}{y(1 + xy)} = -\frac{2}{y}$$

This depends purely on $y$! The integrating factor is:

$$\mu(y) = \exp\left(\int -\frac{2}{y} dy\right) = e^{-2 \ln |y|} = \frac{1}{y^2}$$



Step 3: Multiply ODE by $\mu(y)$

$$\frac{y + x y^2}{y^2} dx - \frac{x}{y^2} dy = 0 \implies \left(\frac{1}{y} + x\right) dx - \frac{x}{y^2} dy = 0$$

Now $M^ = \frac{1}{y} + x$ and $N^ = -\frac{x}{y^2}$. Note $M^_y = -1/y^2 = N^_x$ (exact!).

Step 4: Integrate

$$\Phi = \int \left(\frac{1}{y} + x\right) dx + g(y) = \frac{x}{y} + \frac{x^2}{2} + g(y)$$
$$\frac{\partial \Phi}{\partial y} = -\frac{x}{y^2} + g'(y) = -\frac{x}{y^2} \implies g'(y) = 0 \implies g(y) = C_0$$

Thus $\frac{x}{y} + \frac{x^2}{2} = C$.

Final Answer & Physical Insight

$\mu(y) = \frac{1}{y^2}$, yielding the general solution $\frac{x}{y} + \frac{x^2}{2} = C$ (or $2x + x^2 y = 2Cy$).

Tier 3 • Honors Challenge Example 2.3: Homogeneous Differential Form & Polar Integrating Factor

Solve the homogeneous initial value problem:

$$(x^2 + 3xy + y^2) dx - x^2 dy = 0, \quad y(1) = 0$$


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$:

$$(x^2 + 3x(vx) + v^2 x^2) dx - x^2 (v dx + x dv) = 0$$

Dividing by $x^2 \neq 0$:

$$(1 + 3v + v^2) dx - v dx - x dv = 0 \implies (1 + 2v + v^2) dx - x dv = 0 \implies (v + 1)^2 dx = x dv$$

Separating variables:

$$\frac{dx}{x} = \frac{dv}{(v + 1)^2} \implies \ln |x| = -\frac{1}{v + 1} + C$$

Back-substituting $v = y/x$:

$$\ln |x| = -\frac{1}{\frac{y}{x} + 1} + C = -\frac{x}{x + y} + C \implies \frac{x}{x + y} + \ln |x| = C$$


Method 2: Homogeneous Integrating Factor
$M = x^2 + 3xy + y^2$, $N = -x^2$.

$$Mx + Ny = x(x^2 + 3xy + y^2) - x^2 y = x^3 + 2x^2 y + x y^2 = x(x + y)^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|}$.

Final Answer & Physical Insight

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|}$.