ODE: Integrating Factors for Non-Exact Equations

First-order differential equations are the building blocks for modeling dynamic systems in engineering, from heat transfer to circuit analysis. While some equations can be solved directly, many resist simple methods. A powerful technique for solving a specific class of stubborn equations is using an integrating factor to transform a non-exact differential equation into an exact one, thereby unlocking a straightforward path to the solution. Mastering this method is essential for any engineer's mathematical toolkit.

The Foundation: Exact Differential Equations

The journey begins with understanding what makes an equation "exact." A first-order differential equation expressed in differential form $M(x,y)dx+N(x,y)dy=0$ is called an exact equation if there exists a potential function $F(x,y)$ such that its total differential $dF=\frac{\partial F}{\partial x}dx+\frac{\partial F}{\partial y}dy$ matches the equation perfectly. In practical terms, this means: $$\frac{\partial M}{\partial y}=\frac{\partial N}{\partial x}$$ If this condition holds, the solution is given implicitly by $F(x,y)=C$, where $C$ is a constant. You find $F$ by integrating $M$ with respect to $x$ (treating $y$ as constant) and $N$ with respect to $y$ (treating $x$ as constant), then reconciling the results.

Consider the equation $2xydx+(x^2-1)dy=0$. Here, $M=2xy$ and $N=x^2-1$. Checking the condition: $\frac{\partial M}{\partial y}=2x$ and $\frac{\partial N}{\partial x}=2x$. They are equal, so the equation is exact. Integrating gives the solution $F(x,y)=x^{2}y-y=C$.

The Problem: Non-Exact Equations

Most equations you encounter will not satisfy the exactness condition. An equation like $ydx+(2x-y{e}^y)dy=0$ is non-exact because $\frac{\partial M}{\partial y}=1$ while $\frac{\partial N}{\partial x}=2$. This mismatch means there is no immediate potential function $F$. The core idea is to multiply the entire non-exact equation by a cleverly chosen function, $\mu (x,y)$, called an integrating factor. This factor acts as a "mathematical catalyst" that makes the equation exact. $$\mu Mdx+\mu Ndy=0\text{becomes exact.}$$ The challenge, of course, is finding the right $\mu$.

Systematic Search: Integrating Factors Dependent on One Variable

While finding a general integrating factor $\mu (x,y)$ can be difficult, there are systematic formulas for two special—and common—cases: when $\mu$ is a function of $x$ only or $y$ only.

Case 1: $\mu (x)$ - Function of $x$ Only An integrating factor that depends only on $x$ exists if the expression $$\frac{\frac{\partial M}{\partial y}-\frac{\partial N}{\partial x}}{N}$$ is a function of $x$ alone (i.e., all $y$ terms cancel). If this condition holds, the integrating factor is given by: $$\mu (x)=\exp \left(\int \frac{\frac{\partial M}{\partial y}-\frac{\partial N}{\partial x}}{N}dx\right)$$

Example: Solve $(3xy+y^2)dx+(x^2+xy)dy=0$.

1. Check exactness: $M=3xy+y^2$, $N=x^2+xy$. Then $\frac{\partial M}{\partial y}=3x+2y$ and $\frac{\partial N}{\partial x}=2x+y$. Not equal, so non-exact.
2. Test for $\mu (x)$: Compute $\frac{M_y-N_x}{N}=\frac{(3x+2y)-(2x+y)}{x^2+xy}=\frac{x+y}{x(x+y)}=\frac{1}{x}$. This is a function of $x$ only, so $\mu (x)$ exists.
3. Find $\mu (x)$: $\mu (x)=\exp \left(\int \frac{1}{x}dx\right)=e^{\ln |x|}=x$.
4. Multiply the original equation by $x$: $(3x^{2}y+x{y}^2)dx+(x^3+x^{2}y)dy=0$.
5. Verify the new equation is exact and solve.

Case 2: $\mu (y)$ - Function of $y$ Only Similarly, an integrating factor that depends only on $y$ exists if the expression $$\frac{\frac{\partial N}{\partial x}-\frac{\partial M}{\partial y}}{M}$$ is a function of $y$ alone. If so, the integrating factor is: $$\mu (y)=\exp \left(\int \frac{\frac{\partial N}{\partial x}-\frac{\partial M}{\partial y}}{M}dy\right)$$

The choice between checking for $\mu (x)$ or $\mu (y)$ first is often a matter of trial and inspection, looking for simplifications in the test quotients.

Transforming and Solving the Exact Equation

Once you have found a valid integrating factor $\mu$ and multiplied it through the original equation, you have created a new, exact equation. You then solve it using the standard method for exact equations:

1. Integrate: Integrate the new $M_{new}=\mu M$ with respect to $x$, holding $y$ constant. This yields $F(x,y)=\int \mu Mdx+h(y)$, where $h(y)$ is an unknown function of $y$.
2. Differentiate: Differentiate this result with respect to $y$ and set it equal to the new $N_{new}=\mu N$: $\frac{\partial F}{\partial y}=\frac{\partial }{\partial y}\left(\int \mu Mdx\right)+h^{\prime }(y)=\mu N$.
3. Solve for $h(y)$: This equation lets you solve for $h^{\prime }(y)$ and then integrate to find $h(y)$.
4. Write Solution: Assemble the complete implicit solution $F(x,y)=C$.

This process transforms the intractable non-exact equation into a solvable one through systematic integration.

Recognizing Applicability and Alternative Methods

A critical skill is recognizing when the integrating factor method is the right tool versus when another approach is preferable. The method is ideal for equations in the standard form $Mdx+Ndy=0$ that fail the exactness test but where the test quotients for $\mu (x)$ or $\mu (y)$ simplify nicely. Always check for exactness first.

However, other methods may be more efficient. If the equation is separable, that method is always simpler. If it is linear in $y$ (of the form $\frac{dy}{dx}+P(x)y=Q(x)$), the dedicated linear integrating factor $\mu (x)=e^{\int P(x)dx}$ is more direct. The technique covered here is a generalization useful for non-exact, non-linear equations where the special one-variable dependency condition is met. If the test quotients are functions of both variables, finding an integrating factor becomes significantly more complex and may not be worth the effort for introductory problem-solving.

Common Pitfalls

1. Misapplying the Condition Formula: The most common error is mixing up the numerator in the test quotients. Remember the pattern: For $\mu (x)$, it's $(M_y-N_x)/N$. For $\mu (y)$, it's $(N_x-M_y)/M$. A sign error here will lead to an incorrect integrating factor.

Correction: Derive the condition from the requirement that $(\mu M{)}_y=(\mu N{)}_x$. For $\mu (x)$, this leads to ${\mu }^{\prime }N=\mu (M_y-N_x)$, which rearranges to ${\mu }^{\prime }/\mu =(M_y-N_x)/N$.

1. Forgetting to Verify Exactness After Multiplication: After finding and applying your integrating factor, you must verify that the new equation is indeed exact before proceeding to solve it. This quick check can catch algebraic mistakes made in computing $\mu$.

Correction: Always recalculate $\frac{\partial (\mu M)}{\partial y}$ and $\frac{\partial (\mu N)}{\partial x}$ and confirm they are equal.

1. Overlooking Simpler Methods: Immediately jumping to the integrating factor method for a non-exact equation can waste time if the equation is separable or linear.

Correction: Develop a diagnostic checklist: 1) Is it separable? 2) Is it linear? 3) Is it exact? 4) If not, can I find an integrating factor $\mu (x)$ or $\mu (y)$?

1. Incorrect Integration During Solution: When solving the new exact equation, a frequent mistake is mishandling the "constant of integration" $h(y)$, either by omitting it or mis-differentiating it.

Correction: Remember that when you integrate a function of $x$ and $y$ with respect to $x$, the "constant" is actually an arbitrary function of $y$.

Summary

• An integrating factor $\mu (x,y)$ is a multiplier that transforms a non-exact differential equation $Mdx+Ndy=0$ into an exact one, enabling a solution via the potential function $F(x,y)$.
• Systematic formulas exist for finding integrating factors that depend on only one variable: $\mu (x)$ exists if $(M_y-N_x)/N$ is a function of $x$ alone, and $\mu (y)$ exists if $(N_x-M_y)/M$ is a function of $y$ alone.
• The solution process involves calculating the integrating factor, multiplying it through the original equation, verifying the resulting equation is exact, and then finding the potential function $F$ through integration and reconciliation.
• This method is a powerful tool for non-exact equations, but you must first check for simpler forms like separable or linear equations, which have more straightforward solution paths.
• Success hinges on carefully applying the condition formulas, accurately performing the integrations, and methodically solving for the potential function.