TheoremComplete

Integrating Factor Theorem

The integrating factor method transforms non-exact equations into exact ones by multiplication with a suitable function. This theorem formalizes when and how such factors can be found.

TheoremIntegrating Factor for Linear Equations

Every linear first-order differential equation:

dydx+p(x)y=q(x)\frac{dy}{dx} + p(x)y = q(x)

where p(x)p(x) is continuous, has an integrating factor:

μ(x)=ep(x)dx\mu(x) = e^{\int p(x)dx}

Multiplying the equation by μ(x)\mu(x) yields:

ddx[μ(x)y]=μ(x)q(x)\frac{d}{dx}[\mu(x)y] = \mu(x)q(x)

which can be integrated directly to obtain the general solution.

This theorem guarantees that every linear first-order ODE can be solved by quadratures (integration), making the class of linear equations completely solvable in principle.

Proof Sketch

Given dydx+p(x)y=q(x)\frac{dy}{dx} + p(x)y = q(x), multiply by μ(x)\mu(x):

μ(x)dydx+μ(x)p(x)y=μ(x)q(x)\mu(x)\frac{dy}{dx} + \mu(x)p(x)y = \mu(x)q(x)

For the left side to equal ddx[μ(x)y]\frac{d}{dx}[\mu(x)y], we need:

μ(x)dydx+μ(x)y=μ(x)dydx+μ(x)p(x)y\mu(x)\frac{dy}{dx} + \mu'(x)y = \mu(x)\frac{dy}{dx} + \mu(x)p(x)y

This requires μ(x)=μ(x)p(x)\mu'(x) = \mu(x)p(x), a separable equation:

dμμ=p(x)dx\frac{d\mu}{\mu} = p(x)dx

Integrating: lnμ=p(x)dx\ln|\mu| = \int p(x)dx, giving μ(x)=ep(x)dx\mu(x) = e^{\int p(x)dx}.

ExampleVariable Coefficient Equation

Solve dydx+2yx=x2\frac{dy}{dx} + \frac{2y}{x} = x^2 for x>0x > 0.

Here p(x)=2/xp(x) = 2/x, so:

μ(x)=e(2/x)dx=e2lnx=x2\mu(x) = e^{\int (2/x)dx} = e^{2\ln x} = x^2

Multiply through:

x2dydx+2xy=x4x^2\frac{dy}{dx} + 2xy = x^4

ddx[x2y]=x4\frac{d}{dx}[x^2 y] = x^4

Integrate:

x2y=x55+Cx^2 y = \frac{x^5}{5} + C

y=x35+Cx2y = \frac{x^3}{5} + \frac{C}{x^2}

TheoremIntegrating Factors for Non-Exact Equations

Consider the equation M(x,y)dx+N(x,y)dy=0M(x,y)dx + N(x,y)dy = 0 which is not exact.

(a) If 1N(MyNx)=h(x)\frac{1}{N}\left(\frac{\partial M}{\partial y} - \frac{\partial N}{\partial x}\right) = h(x) depends only on xx, then:

μ(x)=eh(x)dx\mu(x) = e^{\int h(x)dx}

is an integrating factor.

(b) If 1M(NxMy)=k(y)\frac{1}{M}\left(\frac{\partial N}{\partial x} - \frac{\partial M}{\partial y}\right) = k(y) depends only on yy, then:

μ(y)=ek(y)dy\mu(y) = e^{\int k(y)dy}

is an integrating factor.

ExampleFinding an Integrating Factor

Consider (3xy+y2)dx+(x2+xy)dy=0(3xy + y^2)dx + (x^2 + xy)dy = 0.

Check exactness: M=3xy+y2M = 3xy + y^2, N=x2+xyN = x^2 + xy

My=3x+2y,Nx=2x+y\frac{\partial M}{\partial y} = 3x + 2y, \quad \frac{\partial N}{\partial x} = 2x + y

Not exact. Try case (a):

1N(MyNx)=(3x+2y)(2x+y)x2+xy=x+yx(x+y)=1x\frac{1}{N}\left(\frac{\partial M}{\partial y} - \frac{\partial N}{\partial x}\right) = \frac{(3x+2y)-(2x+y)}{x^2+xy} = \frac{x+y}{x(x+y)} = \frac{1}{x}

This depends only on xx! So μ(x)=e(1/x)dx=x\mu(x) = e^{\int (1/x)dx} = x.

Multiply the original equation by xx:

(3x2y+xy2)dx+(x3+x2y)dy=0(3x^2y + xy^2)dx + (x^3 + x^2y)dy = 0

Now My=3x2+2xy=Nx\frac{\partial M}{\partial y} = 3x^2 + 2xy = \frac{\partial N}{\partial x}. Exact!

Remark

While the theorem provides conditions for finding integrating factors depending on one variable, finding general integrating factors (depending on both xx and yy) is usually very difficult. The cases where μ=μ(x)\mu = \mu(x) or μ=μ(y)\mu = \mu(y) are the most practically useful.

The integrating factor theorem bridges the gap between exact and non-exact equations, showing that many seemingly different equation types can be unified under the exact equation framework.