MathIsimple
ODE-03
First-Order Methods

Exact Equations and Integrating Factors

Learn to identify and solve exact differential equations, find integrating factors to convert non-exact equations, and understand first-order implicit equations including the elegant Clairaut equation.

Learning Objectives

  • Identify and solve exact differential equations
  • Apply the exactness criterion for differential forms
  • Master common total differential formulas
  • Find integrating factors depending on one variable
  • Use the grouped integrating factor method
  • Solve Clairaut equations and find singular solutions
  • Understand the geometric meaning of envelopes
  • Apply integrating factor techniques to applications

1. Exact Differential Equations

A first-order differential equation in the form M(t,x)dt+N(t,x)dx=0M(t,x)dt + N(t,x)dx = 0 is called exactif the left side is the total differential of some function u(t,x)u(t,x). This means we can find the general solution by simply finding this potential function.

Definition 3.1: Exact Equation

A differential equation M(t,x)dt+N(t,x)dx=0M(t,x)dt + N(t,x)dx = 0 is called an exact equation(or total differential equation) if there exists a continuously differentiable function u(t,x)u(t,x) such that:

du(t,x)=M(t,x)dt+N(t,x)dxdu(t,x) = M(t,x)dt + N(t,x)dx

Equivalently, ut=M(t,x)\frac{\partial u}{\partial t} = M(t,x) and ux=N(t,x)\frac{\partial u}{\partial x} = N(t,x).

The key question is: how do we determine if an equation is exact, and how do we find the potential function uu?

Theorem 3.1: Exactness Criterion

Let M(t,x),N(t,x)C(D)M(t,x), N(t,x) \in C(D) where DR2D \subset \mathbb{R}^2 is a simply connected domain, and suppose Mx\frac{\partial M}{\partial x} and Nt\frac{\partial N}{\partial t} are continuous in DD. Then the equation M,dt+N,dx=0M\\,dt + N\\,dx = 0 is exact if and only if:

Mx=Nt\frac{\partial M}{\partial x} = \frac{\partial N}{\partial t}

When this condition holds, the general integral is given by:

u(t,x)=t0tM(s,x)ds+x0xN(t0,s)ds=Cu(t,x) = \int_{t_0}^{t} M(s, x)\,ds + \int_{x_0}^{x} N(t_0, s)\,ds = C
Proof of Theorem 3.1:

Necessity: If the equation is exact, there exists u(t,x)u(t,x) such thatut=M\frac{\partial u}{\partial t} = M and ux=N\frac{\partial u}{\partial x} = N. Since uu has continuous second partial derivatives:

Mx=2uxt=2utx=Nt\frac{\partial M}{\partial x} = \frac{\partial^2 u}{\partial x \partial t} = \frac{\partial^2 u}{\partial t \partial x} = \frac{\partial N}{\partial t}

Sufficiency: Assume Mx=Nt\frac{\partial M}{\partial x} = \frac{\partial N}{\partial t}. We construct u(t,x)u(t,x) as follows. Let:

u(t,x)=t0tM(s,x)ds+φ(x)u(t,x) = \int_{t_0}^{t} M(s, x)\,ds + \varphi(x)

where φ(x)\varphi(x) is a function to be determined. Then ut=M(t,x)\frac{\partial u}{\partial t} = M(t,x). For ux=N\frac{\partial u}{\partial x} = N:

ux=t0tMx(s,x)ds+φ(x)=t0tNt(s,x)ds+φ(x)\frac{\partial u}{\partial x} = \int_{t_0}^{t} \frac{\partial M}{\partial x}(s,x)\,ds + \varphi'(x) = \int_{t_0}^{t} \frac{\partial N}{\partial t}(s,x)\,ds + \varphi'(x)
=N(t,x)N(t0,x)+φ(x)= N(t,x) - N(t_0,x) + \varphi'(x)

Setting φ(x)=N(t0,x)\varphi'(x) = N(t_0,x), we get ux=N(t,x)\frac{\partial u}{\partial x} = N(t,x). ∎

Example 3.1: Solving an Exact Equation

Problem: Solve (2tsinx+3t2x+t4)dt+(t3+t2cosx+x2)dx=0(2t\sin x + 3t^2 x + t^4)dt + (t^3 + t^2\cos x + x^2)dx = 0.

Solution:

First, verify exactness. Let M=2tsinx+3t2x+t4M = 2t\sin x + 3t^2 x + t^4 and N=t3+t2cosx+x2N = t^3 + t^2\cos x + x^2.

Mx=2tcosx+3t2,Nt=3t2+2tcosx\frac{\partial M}{\partial x} = 2t\cos x + 3t^2, \quad \frac{\partial N}{\partial t} = 3t^2 + 2t\cos x

Since Mx=Nt\frac{\partial M}{\partial x} = \frac{\partial N}{\partial t}, the equation is exact.

Taking (t0,x0)=(0,0)(t_0, x_0) = (0, 0), we compute:

u(t,x)=0t(2ssinx+3s2x+s4)ds+0xs2dsu(t,x) = \int_0^t (2s\sin x + 3s^2 x + s^4)\,ds + \int_0^x s^2\,ds
=t2sinx+t3x+t55+x33= t^2\sin x + t^3 x + \frac{t^5}{5} + \frac{x^3}{3}

The general integral is:

t2sinx+t3x+t55+x33=Ct^2\sin x + t^3 x + \frac{t^5}{5} + \frac{x^3}{3} = C

2. Common Total Differential Formulas

Recognizing common total differentials can greatly simplify solving exact equations. Here are the most frequently used formulas:

Essential Total Differential Formulas

1.tdt+xdx=12d(t2+x2)t\,dt + x\,dx = \frac{1}{2}d(t^2 + x^2)
2.xdt+tdxtx=d(ln(tx))\frac{x\,dt + t\,dx}{tx} = d(\ln(tx))
3.xdttdxt2+x2=d(arctantx)\frac{x\,dt - t\,dx}{t^2 + x^2} = d\left(\arctan\frac{t}{x}\right)
4.xdttdxx2=d(tx)\frac{x\,dt - t\,dx}{x^2} = d\left(\frac{t}{x}\right)
5.xdttdxt2=d(xt)\frac{x\,dt - t\,dx}{-t^2} = d\left(\frac{x}{t}\right)
6.tdt+xdxt2+x2=12dln(t2+x2)\frac{t\,dt + x\,dx}{t^2 + x^2} = \frac{1}{2}d\ln(t^2 + x^2)
7.xdttdxt2x2=12dlntxt+x\frac{x\,dt - t\,dx}{t^2 - x^2} = \frac{1}{2}d\ln\left|\frac{t-x}{t+x}\right|
Remark:

These formulas are particularly useful when you can identify parts of a differential equation that match these patterns. Often, the key to solving an equation is rewriting it to reveal these standard forms.

3. Integrating Factors

Most first-order equations are not exact. However, we can often multiply by a suitable function called anintegrating factor to make the equation exact.

Definition 3.2: Integrating Factor

For the equation M(t,x)dt+N(t,x)dx=0M(t,x)dt + N(t,x)dx = 0, a non-zero function μ(t,x)C(D)\mu(t,x) \in C(D) is called an integrating factor if the equation:

μ(t,x)M(t,x)dt+μ(t,x)N(t,x)dx=0\mu(t,x)M(t,x)dt + \mu(t,x)N(t,x)dx = 0

is an exact equation.

Theorem 3.2: Condition for Integrating Factor

A function μ(t,x)\mu(t,x) is an integrating factor for M,dt+N,dx=0M\\,dt + N\\,dx = 0 if and only if:

NμtMμx=μ(MxNt)N\frac{\partial \mu}{\partial t} - M\frac{\partial \mu}{\partial x} = \mu\left(\frac{\partial M}{\partial x} - \frac{\partial N}{\partial t}\right)

This PDE is generally harder to solve than the original equation. However, in special cases, we can find integrating factors that depend on only one variable.

Theorem 3.3: Integrating Factor Depending on One Variable

Case 1: If 1N(MxNt)=G(t)\frac{1}{N}\left(\frac{\partial M}{\partial x} - \frac{\partial N}{\partial t}\right) = G(t) depends only on tt, then:

μ(t)=eG(t)dt\mu(t) = e^{\int G(t)\,dt}

is an integrating factor.

Case 2: If 1M(MxNt)=H(x)-\frac{1}{M}\left(\frac{\partial M}{\partial x} - \frac{\partial N}{\partial t}\right) = H(x) depends only on xx, then:

μ(x)=eH(x)dx\mu(x) = e^{\int H(x)\,dx}

is an integrating factor.

Proof of Theorem 3.3:

We prove Case 2. Suppose μ=μ(x)\mu = \mu(x) depends only on xx. Then μt=0\frac{\partial \mu}{\partial t} = 0, and the condition becomes:

Mdμdx=μ(MxNt)-M\frac{d\mu}{dx} = \mu\left(\frac{\partial M}{\partial x} - \frac{\partial N}{\partial t}\right)

Rearranging:

1μdμdx=1M(MxNt)=H(x)\frac{1}{\mu}\frac{d\mu}{dx} = -\frac{1}{M}\left(\frac{\partial M}{\partial x} - \frac{\partial N}{\partial t}\right) = H(x)

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

Algorithm 3.1: Finding an Integrating Factor

Input: Equation M(t,x)dt+N(t,x)dx=0M(t,x)dt + N(t,x)dx = 0

Output: An integrating factor μ\mu if one exists in simple form

1. Check if equation is already exact: Mx=?Nt\frac{\partial M}{\partial x} \stackrel{?}{=} \frac{\partial N}{\partial t}

2. Compute R=MxNtR = \frac{\partial M}{\partial x} - \frac{\partial N}{\partial t}

3. If R/NR/N depends only on tt: μ(t)=e(R/N)dt\mu(t) = e^{\int (R/N)\,dt}

4. Else if R/M-R/M depends only on xx: μ(x)=e(R/M)dx\mu(x) = e^{\int (-R/M)\,dx}

5. Else try special forms: μ(tx),μ(t2+x2),μ(t/x)\mu(tx), \mu(t^2+x^2), \mu(t/x), etc.

6. return Integrating factor μ\mu

Example 3.2: Finding an Integrating Factor

Problem: Solve x,dt+(t2+t+x2)dx=0x\\,dt + (t^2 + t + x^2)dx = 0.

Solution:

Here M=xM = x and N=t2+t+x2N = t^2 + t + x^2. Check exactness:

Mx=1,Nt=2t+1\frac{\partial M}{\partial x} = 1, \quad \frac{\partial N}{\partial t} = 2t + 1

Not exact. Let's make the substitution y=t+1y = t + 1 (so dy=dtdy = dt):

xdy+(y2y+x2)dx=0x\,dy + (y^2 - y + x^2)dx = 0

Rewrite as:

xdyydx+(x2+y2)dx=0x\,dy - y\,dx + (x^2 + y^2)dx = 0

Using the integrating factor μ=1x2+y2\mu = \frac{1}{x^2 + y^2}:

xdyydxx2+y2+dx=0\frac{x\,dy - y\,dx}{x^2 + y^2} + dx = 0

This gives d(arctanyx)+dx=0d\left(\arctan\frac{y}{x}\right) + dx = 0, so:

arctanyx+x=C    arctant+1x+x=C\arctan\frac{y}{x} + x = C \implies \arctan\frac{t+1}{x} + x = C

4. Grouped Integrating Factor Method

For equations that can be split into several groups, each with its own integrating factor, we can use a powerful technique to find a common integrating factor.

Theorem 3.4: Grouped Integrating Factor Theorem

Suppose the non-exact equation i=1n(Midt+Nidx)=0\sum_{i=1}^{n}(M_i\,dt + N_i\,dx) = 0 can be grouped such that each group Mi,dt+Ni,dx=0M_i\\,dt + N_i\\,dx = 0 has integrating factor μi\mu_i satisfying:

μi(Midt+Nidx)=dΦi(t,x)\mu_i(M_i\,dt + N_i\,dx) = d\Phi_i(t,x)

Then for any differentiable function gg, the function μig(Φi)\mu_i g(\Phi_i) is also an integrating factor for the ii-th group equation.

Proof of Theorem 3.4:

For any i=1,,ni = 1, \ldots, n, we compute:

x(μig(Φi)Mi)=x(μiMi)g(Φi)+Φixg(Φi)μiMi\frac{\partial}{\partial x}(\mu_i g(\Phi_i)M_i) = \frac{\partial}{\partial x}(\mu_i M_i)g(\Phi_i) + \frac{\partial \Phi_i}{\partial x}g'(\Phi_i)\mu_i M_i

Since μiMi=Φit\mu_i M_i = \frac{\partial \Phi_i}{\partial t} and μiNi=Φix\mu_i N_i = \frac{\partial \Phi_i}{\partial x}:

=t(μiNi)g(Φi)+μi2MiNig(Φi)= \frac{\partial}{\partial t}(\mu_i N_i)g(\Phi_i) + \mu_i^2 M_i N_i g'(\Phi_i)
=t(μiNi)g(Φi)+Φitg(Φi)μiNi=t(μig(Φi)Ni)= \frac{\partial}{\partial t}(\mu_i N_i)g(\Phi_i) + \frac{\partial \Phi_i}{\partial t}g'(\Phi_i)\mu_i N_i = \frac{\partial}{\partial t}(\mu_i g(\Phi_i)N_i)

Thus μig(Φi)\mu_i g(\Phi_i) is an integrating factor. ∎

Algorithm 3.2: Grouped Integrating Factor Method

Input: Equation that can be split as (Midt+Nidx)=0\sum (M_i\,dt + N_i\,dx) = 0

Output: General integral

1. Split the equation into groups

2. For each group, find μi\mu_i and Φi\Phi_i

3. Choose functions gig_i such that μ1g1(Φ1)=μ2g2(Φ2)==μ\mu_1 g_1(\Phi_1) = \mu_2 g_2(\Phi_2) = \cdots = \mu

4. Multiply original equation by common μ\mu

5. Integrate to find general integral

Example 3.3: Grouped Integrating Factor Method

Problem: Solve (xt+3t2)dt+(1+t3x)dx=0\left(\frac{x}{t} + 3t^2\right)dt + \left(1 + \frac{t^3}{x}\right)dx = 0.

Solution:

Group the equation as:

(xtdt+dx)+(3t2dt+t3xdx)=0\left(\frac{x}{t}\,dt + dx\right) + \left(3t^2\,dt + \frac{t^3}{x}\,dx\right) = 0

Group 1: xtdt+dx=0\frac{x}{t}\,dt + dx = 0

Integrating factor: μ1=t\mu_1 = t, giving xdt+tdx=d(tx)x\,dt + t\,dx = d(tx), so Φ1=tx\Phi_1 = tx.

Group 2: 3t2dt+t3xdx=03t^2\,dt + \frac{t^3}{x}\,dx = 0

Integrating factor: μ2=x\mu_2 = x, giving 3t2xdt+t3dx=d(t3x)3t^2 x\,dt + t^3\,dx = d(t^3 x), so Φ2=t3x\Phi_2 = t^3 x.

Find g1,g2g_1, g_2 such that tg1(tx)=xg2(t3x)tg_1(tx) = xg_2(t^3 x).

Take g1(s)=s2g_1(s) = s^2 and g2(s)=sg_2(s) = s, giving μ=t(tx)2=t3x2\mu = t(tx)^2 = t^3 x^2.

Multiply original equation by μ=t3x2\mu = t^3 x^2:

(x3t2dt+t3x2dx)+(3t5x2dt+t6xdx)=0(x^3 t^2\,dt + t^3 x^2\,dx) + (3t^5 x^2\,dt + t^6 x\,dx) = 0
=d(13x3t3)+d(12t6x2)=0= d\left(\frac{1}{3}x^3 t^3\right) + d\left(\frac{1}{2}t^6 x^2\right) = 0

The general integral is:

13x3t3+12t6x2=C\frac{1}{3}x^3 t^3 + \frac{1}{2}t^6 x^2 = C

5. First-Order Implicit Equations and Clairaut Equations

When the derivative xx' cannot be explicitly solved, we have an implicit equation F(t,x,x)=0F(t, x, x') = 0. A particularly elegant case is the Clairaut equation.

Definition 3.3: Clairaut Equation

A Clairaut equation has the form:

x=tx+f(x)x = tx' + f(x')

where ff is a given function with f≢0f'' \not\equiv 0. Setting p=xp = x', the equation becomes x=tp+f(p)x = tp + f(p).

Theorem 3.5: Solution of Clairaut Equations

The Clairaut equation x=tx+f(x)x = tx' + f(x') has:

  • General solution: x=Ct+f(C)x = Ct + f(C) — a family of straight lines
  • Singular solution: Given parametrically by t=f(p)t = -f'(p), x=pf(p)+f(p)x = -pf'(p) + f(p)
Proof of Theorem 3.5:

Let p=xp = x'. Differentiate x=tp+f(p)x = tp + f(p) with respect to tt:

p=p+tdpdt+f(p)dpdt=p+(t+f(p))dpdtp = p + t\frac{dp}{dt} + f'(p)\frac{dp}{dt} = p + (t + f'(p))\frac{dp}{dt}

This gives:

(t+f(p))dpdt=0(t + f'(p))\frac{dp}{dt} = 0

Case 1: dpdt=0\frac{dp}{dt} = 0, so p=Cp = C (constant). Substituting back: x=Ct+f(C)x = Ct + f(C).

Case 2: t+f(p)=0t + f'(p) = 0, giving t=f(p)t = -f'(p). Substituting into the original equation:

x=pf(p)+f(p)x = -pf'(p) + f(p)

This parametric form defines the singular solution. ∎

Definition 3.4: Singular Solution

A singular solution of a differential equation is a solution that cannot be obtained from the general solution by any choice of the arbitrary constant, and at every point on its curve, it is tangent to some member of the family of curves representing the general solution.

Geometrically, the singular solution is the envelope of the family of general solution curves.

Example 3.4: Solving a Clairaut Equation

Problem: Solve x=tx14(x)2x = tx' - \frac{1}{4}(x')^2.

Solution:

This is a Clairaut equation with f(p)=14p2f(p) = -\frac{1}{4}p^2.

General solution: x=Ct14C2x = Ct - \frac{1}{4}C^2 (a family of lines).

Singular solution: Since f(p)=12pf'(p) = -\frac{1}{2}p, we have:

t=f(p)=12p    p=2tt = -f'(p) = \frac{1}{2}p \implies p = 2t

Substituting:

x=t(2t)14(2t)2=2t2t2=t2x = t(2t) - \frac{1}{4}(2t)^2 = 2t^2 - t^2 = t^2

The singular solution is x=t2x = t^2 (a parabola).

Geometric interpretation: The parabola x=t2x = t^2 is the envelope of the family of lines x=CtC24x = Ct - \frac{C^2}{4}. Each line is tangent to the parabola.

Remark:

The singular solution often has physical significance. For example, in optics, if each line represents a light ray, the envelope (singular solution) represents a caustic curve where light concentrates.

6. Applications

Example 3.5: Searchlight Mirror Shape

Problem: Find the shape of a mirror such that light rays from a point source at the origin reflect as parallel rays.

Solution:

Let point P(t,x)P(t, x) be on the mirror curve. By the law of reflection, if α\alpha is the angle between the tangent and the reflected ray (parallel to tt-axis), then:

x(t)=tanαx'(t) = \tan\alpha

Since the angle of incidence equals the angle of reflection, OA=OP|OA| = |OP| where AA is on the tt-axis. This gives:

xxt=t2+x2\frac{x}{x'} - t = \sqrt{t^2 + x^2}

Let u=t/xu = t/x. After substitution and integration:

x=2Ct+C2x = 2Ct + \frac{C}{2}

This is a parabola, confirming that parabolic mirrors produce parallel reflected rays from a point source at the focus.

Practice Quiz

Exact Equations and Integrating Factors Quiz
10
Questions
0
Correct
0%
Accuracy
1
For the equation M(t,x)dt+N(t,x)dx=0M(t,x)dt + N(t,x)dx = 0 to be exact, which condition must hold?
Easy
Not attempted
2
Which of the following is a common total differential formula?
Easy
Not attempted
3
What is an integrating factor for the equation Mdt+Ndx=0M\,dt + N\,dx = 0?
Easy
Not attempted
4
Is the equation (2tsinx+3t2x)dt+(t2cosx+t3)dx=0(2t\sin x + 3t^2x)dt + (t^2\cos x + t^3)dx = 0 exact?
Medium
Not attempted
5
If 1M(MxNt)\frac{1}{M}\left(\frac{\partial M}{\partial x} - \frac{\partial N}{\partial t}\right) depends only on xx, then the integrating factor μ\mu depends on:
Medium
Not attempted
6
The general form of a Clairaut equation is:
Medium
Not attempted
7
What type of solutions does a Clairaut equation have?
Medium
Not attempted
8
For the equation xdt+(t2+t+x2)dx=0x\,dt + (t^2 + t + x^2)dx = 0, which substitution helps find an integrating factor?
Medium
Not attempted
9
In the grouped integrating factor method, if equation Midt+Nidx=0M_i\,dt + N_i\,dx = 0 has integrating factor μi\mu_i with first integral Φi\Phi_i, then μig(Φi)\mu_i g(\Phi_i) is also an integrating factor for any:
Hard
Not attempted
10
The singular solution of a Clairaut equation is geometrically the:
Hard
Not attempted

Frequently Asked Questions

How do I know if an equation is exact?

Check if Mx=Nt\frac{\partial M}{\partial x} = \frac{\partial N}{\partial t}. If this equality holds and the domain is simply connected, the equation M,dt+N,dx=0M\\,dt + N\\,dx = 0 is exact. Remember that MMis the coefficient of dtdt and NN is the coefficient of dxdx.

What if I can't find an integrating factor?

Not every equation has a simple integrating factor. Try these approaches in order: (1) Check if R/NR/Ndepends only on tt, (2) Check if R/M-R/M depends only on xx, (3) Try special forms likeμ(tx)\mu(tx) or μ(t2+x2)\mu(t^2+x^2), (4) Try the grouped integrating factor method. If none work, the equation may require a different solution technique.

What is the geometric meaning of the singular solution?

The singular solution is the envelope of the family of curves that form the general solution. Geometrically, at each point of the singular solution curve, exactly one member of the general solution family is tangent to it. This is why the singular solution cannot be obtained by choosing a specific constant in the general solution.

How do I verify my solution to an exact equation?

After finding u(t,x)=Cu(t,x) = C, verify by checking: (1) ut=M\frac{\partial u}{\partial t} = M and (2) ux=N\frac{\partial u}{\partial x} = N. Also, you can differentiate u(t,x)=Cu(t,x) = C implicitly to recover the original differential equation.

When should I use the grouped integrating factor method?

Use this method when: (1) The equation naturally splits into groups, (2) Each group has a recognizable integrating factor, (3) Standard integrating factor methods don't work. The key is finding functionsgig_i such that all μigi(Φi)\mu_i g_i(\Phi_i) are equal.

Can a Clairaut equation have no singular solution?

Yes, if f(p)0f''(p) \equiv 0, meaning f(p)f(p) is linear in pp. In this case, the equation reduces to x=tp+ap+b=(t+a)p+bx = tp + ap + b = (t+a)p + b, which is a family of lines through a single point, and there is no envelope (singular solution).