Integrating factor

Navier–Stokes differential equations used to simulate airflow around an obstruction.

Scope

Natural sciences Engineering
Astronomy Physics Chemistry Biology Geology
Applied mathematics
Continuum mechanics Chaos theory Dynamical systems
Social sciences
Economics Population dynamics

Classification

Types

By variable type
Ordinary Partial Differential-algebraic Integro-differential Fractional Linear Non-linear
Dependent and independent variables Autonomous Complex Coupled / Decoupled Exact Homogeneous / Nonhomogeneous
Features
Order Operator Notation

Relation to processes

Difference (discrete analogue)

Solution

General topics

Lyapunov / Asymptotic / Exponential stability
Rate of convergence
Series / Integral solutions

Solution methods

People

In mathematics, an integrating factor is a function that is chosen to facilitate the solving of a given equation involving differentials. It is commonly used to solve ordinary differential equations, but is also used within multivariable calculus when multiplying through by an integrating factor allows an inexact differential to be made into an exact differential (which can then be integrated to give a scalar field). This is especially useful in thermodynamics where temperature becomes the integrating factor that makes entropy an exact differential.

Use in solving first order linear ordinary differential equations

Integrating factors are useful for solving ordinary differential equations that can be expressed in the form

y'+P(x)y=Q(x)

The basic idea is to find some function $M(x)$ , called the "integrating factor," which we can multiply through our DE in order to bring the left-hand side under a common derivative. For the canonical first-order, linear differential equation shown above, our integrating factor is chosen to be

M(x)=e^{{\int _{{s_{0}}}^{{x}}P(s)ds}}

In order to derive this, let $M(x)$ be the integrating factor of a first order, linear differential equation such that multiplication by $M(x)$ transforms a partial derivative into a total derivative, then:

${\begin{aligned}(1)\qquad &M(x){\underset {{\text{partial derivative}}}{(\underbrace {y'+P(x)y})}}\\(2)\qquad &M(x)y'+M(x)P(x)y\\(3)\qquad &{\underset {{\text{total derivative}}}{\underbrace {M(x)y'+M'(x)y}}}\end{aligned}}$

Going from step 2 to step 3 requires that $M(x)P(x)=M'(x)$ , which is a separable differential equation, whose solution yields $M(x)$ in terms of $P(x)$ :

${\begin{aligned}(4)\qquad &M(x)P(x)=M'(x)\\(5)\qquad &P(x)={\frac {M'(x)}{M(x)}}\\(6)\qquad &\int P(x)dx=\ln M(x)\\(7)\qquad &e^{{\int P(x)dx}}=M(x)\end{aligned}}$

To verify see that multiplying through by $M(x)$ gives

y'e^{{\int _{{s_{0}}}^{{x}}P(s)ds}}+P(x)ye^{{\int _{{s_{0}}}^{{x}}P(s)ds}}=Q(x)e^{{\int _{{s_{0}}}^{{x}}P(s)ds}}

By applying the product rule in reverse, we see that the left-hand side can be expressed as a single derivative in $x$

y'e^{{\int _{{s_{0}}}^{{x}}P(s)ds}}+P(x)ye^{{\int _{{s_{0}}}^{{x}}P(s)ds}}={\frac {d}{dx}}(ye^{{\int _{{s_{0}}}^{{x}}P(s)ds}})

We use this fact to simplify our expression to

\frac{d}{dx}\left(y e^{\int_{s_0}^{x} P(s) ds}\right) = Q(x) e^{\int_{s_0}^{x} P(s) ds}

We then integrate both sides with respect to $x$ , firstly by renaming $x$ to $t$ , obtaining

ye^{{\int _{{s_{0}}}^{{x}}P(s)ds}}=\int _{{t_{0}}}^{{x}}Q(t)e^{{\int _{{s_{0}}}^{{t}}P(s)ds}}dt+C

Finally, we can move the exponential to the right-hand side to find a general solution to our ODE:

y=e^{{-\int _{{s_{0}}}^{{x}}P(s)ds}}\int _{{t_{0}}}^{x}Q(t)e^{{\int _{{s_{0}}}^{{t}}P(s)ds}}dt+Ce^{{-\int _{{s_{0}}}^{{x}}P(s)ds}}

In the case of a homogeneous differential equation, in which $Q(x)=0$ , we find that

y={\frac {C}{e^{{\int _{{s_{0}}}^{{x}}P(s)ds}}}}

where $C$ is a constant.

Example

Solve the differential equation

y'-{\frac {2y}{x}}=0.

We can see that in this case $P(x)={\frac {-2}{x}}$

M(x)=e^{{\int P(x)\,dx}}

M(x)=e^{{\int {\frac {-2}{x}}\,dx}}=e^{{-2\ln x}}={(e^{{\ln x}})}^{{-2}}=x^{{-2}}

(Note we do not need to include the integrating constant - we need only a solution, not the general solution)

M(x)={\frac {1}{x^{2}}}.

Multiplying both sides by $M(x)$ we obtain

{\frac {y'}{x^{2}}}-{\frac {2y}{x^{3}}}=0

{\frac {y'x^{3}-2x^{2}y}{x^{5}}}=0

{\frac {x(y'x^{2}-2xy)}{x^{5}}}=0

{\frac {y'x^{2}-2xy}{x^{4}}}=0.

Reversing the quotient rule gives

\left({\frac {y}{x^{2}}}\right)'=0

{\frac {y}{x^{2}}}=C\,

which gives

y\left(x\right)=Cx^{2}.

General use

An integrating factor is any expression that a differential equation is multiplied by to facilitate integration and is not restricted to first order linear equations. For example, the nonlinear second order equation