CP 9.1 - The integrating factor method

CP 9.1 - The integrating factor method

A carefully chosen multiplier turns a first-order linear differential equation into one product derivative. You will recognise when the integrating factor method applies, put the equation into the right form, and use constant, algebraic and trigonometric factors to obtain its solution.

Recognising the form and finding the factor

How would you solve y′+2y=2xy'+2y=2x? Integrating y′y' is easy, but the extra 2y2y is an unknown function too. Can one multiplier make both left-hand terms integrate together? The product rule supplies exactly that multiplier.

You already know

The product rule gives (uv)′=uv′+u′v(uv)'=uv'+u'v. Separating variables works when you can rearrange an equation into g(y) dy=f(x) dxg(y)\,\mathrm{d}y=f(x)\,\mathrm{d}x. The equation y′+2y=2xy'+2y=2x does not separate in that way, so it needs a different method.

A first-order linear equation has the form y′+P(x)y=Q(x)y'+P(x)y=Q(x): the unknown function and its first derivative occur to the first power, with coefficients depending only on xx. Thus y′+x2y=sin⁡xy'+x^2y=\sin x is linear; y′+y2=xy'+y^2=x is not. Divide by the coefficient of y′y' before identifying PP, on an interval where that coefficient is non-zero.

Multiply the standard form by a function μ\mu. Its left side becomes μy′+μPy\mu y'+\mu Py. This is (μy)′(\mu y)' precisely when μ′=Pμ\mu'=P\mu. Dividing by μ\mu and integrating suggests ln⁡μ=∫P dx\ln\mu=\int P\,\mathrm{d}x.

Integrating factor

For y′+P(x)y=Q(x)y'+P(x)y=Q(x), on an interval where the coefficients are continuous, an integrating factor is

μ(x)=e∫P(x) dx\mu(x)=\mathrm{e}^{\int P(x)\,\mathrm{d}x}

It gives

ddx(μy)=μQ\frac{\mathrm{d}}{\mathrm{d}x}(\mu y)=\mu Q

You may quote the factor without proof. Choose one antiderivative when calculating it; no arbitrary constant is needed inside the factor. Not in the formula booklet: learn it.

The factor comes from the coefficient of yy after normalising the derivative. Keep the sign of that coefficient.

Worked example 1: a constant coefficient

Question. Put 2y′+6y=4x2y'+6y=4x into standard linear form and find an integrating factor.

Plan. Divide every term by 22, identify PP, and exponentiate its integral.

y′+3y=2xdivide through by 2P(x)=3coefficient of yμ=e∫3 dx=e3xintegrate then exponentiate\begin{aligned} y'+3y&=2x&&\text{divide through by }2\\ P(x)&=3&&\text{coefficient of }y\\ \mu&=\mathrm{e}^{\int3\,\mathrm{d}x}=\mathrm{e}^{3x}&&\text{integrate then exponentiate} \end{aligned}

Why this step? Using 66 before dividing would give e6x\mathrm{e}^{6x}, which would not combine the normalised left side into one product derivative.

Answer. Standard form is y′+3y=2xy'+3y=2x; an integrating factor is e3x\mathrm{e}^{3x}.

Check. (e3xy)′=e3xy′+3e3xy(\mathrm{e}^{3x}y)'=\mathrm{e}^{3x}y'+3\mathrm{e}^{3x}y, exactly the multiplied left side.

The second example changes the derivative coefficient to a variable. Dividing now makes the coefficient of yy a reciprocal function.

Worked example 2: normalising a variable coefficient

Question. For x>0x>0, put xy′+3y=4xxy'+3y=4x into standard form and find an integrating factor.

Plan. Divide by xx, then use the logarithmic antiderivative of 3/x3/x.

y′+3xy=4divide every term by xμ=e∫(3/x) dx=e3ln⁡xx>0=x3inverse exponential and logarithm\begin{aligned} y'+\frac3x y&=4&&\text{divide every term by }x\\ \mu&=\mathrm{e}^{\int(3/x)\,\mathrm{d}x}=\mathrm{e}^{3\ln x}&&x>0\\ &=x^3&&\text{inverse exponential and logarithm} \end{aligned}

Why this step? P(x)=3/xP(x)=3/x, rather than 33. Its integral is 3ln⁡x3\ln x, so the factor is a power of xx, not e3x\mathrm{e}^{3x}.

Answer. y′+3xy=4y'+\frac3x y=4, with integrating factor x3x^3.

Check. Multiplying by x3x^3 gives x3y′+3x2y=(x3y)′x^3y'+3x^2y=(x^3y)'.

There are infinitely many valid integrating factors: any non-zero constant multiple of μ\mu works too. For example, 7e3x7\mathrm{e}^{3x} gives seven times the same product derivative. That constant disappears into the final arbitrary constant, leaving exactly the same solution family.

Watch out

Wrong: read P=3P=3 directly from xy′+3y=4xxy'+3y=4x. Right: first divide by xx, obtaining P=3/xP=3/x. State the interval when dividing by a variable expression.

Change the coefficient while keeping the normalisation step visible.

Turning the equation into a product derivative

Once you have the factor, multiply the whole equation by it, integrate both sides, and divide to make yy the subject.

The left side is no longer two separate integration problems: it is one derivative. Write that derivative explicitly, so the structure remains visible even if the right-hand integration needs work.

Worked example 3: the opening equation

Question. Find the general solution of y′+2y=2xy'+2y=2x.

Plan. Multiply by e2x\mathrm{e}^{2x}, integrate the product derivative, and use integration by parts on the right.

μ=e∫2 dx=e2x(e2xy)′=2xe2xmultiply every terme2xy=∫2xe2x dxintegrate the derivative=xe2x−12e2x+Cintegration by partsy=x−12+Ce−2xdivide by the factor\begin{aligned} \mu&=\mathrm{e}^{\int2\,\mathrm{d}x}=\mathrm{e}^{2x}\\ (\mathrm{e}^{2x}y)'&=2x\mathrm{e}^{2x}&&\text{multiply every term}\\ \mathrm{e}^{2x}y&=\int2x\mathrm{e}^{2x}\,\mathrm{d}x&&\text{integrate the derivative}\\ &=x\mathrm{e}^{2x}-\frac12\mathrm{e}^{2x}+C&&\text{integration by parts}\\ y&=x-\frac12+C\mathrm{e}^{-2x}&&\text{divide by the factor} \end{aligned}

Why this step? ∫(e2xy)′ dx=e2xy\int(\mathrm{e}^{2x}y)'\,\mathrm{d}x=\mathrm{e}^{2x}y, so you never need to integrate the unknown yy separately.

Answer. y=x−12+Ce−2xy=x-\frac12+C\mathrm{e}^{-2x}, where CC is arbitrary.

Check. y′=1−2Ce−2xy'=1-2C\mathrm{e}^{-2x}, so y′+2y=1−2Ce−2x+2x−1+2Ce−2x=2xy'+2y=1-2C\mathrm{e}^{-2x}+2x-1+2C\mathrm{e}^{-2x}=2x.

The second example changes only the forcing term, from 2x2x to ex\mathrm{e}^x. The coefficient of yy, and therefore the integrating factor, stays the same.

Worked example 4: exponential forcing

Question. Find the general solution of y′+2y=exy'+2y=\mathrm{e}^x.

Plan. Use e2x\mathrm{e}^{2x} again, combine the exponential factors and integrate.

(e2xy)′=e2xex=e3xmultiply throughe2xy=13e3x+Cintegratey=13ex+Ce−2xdivide by e2x\begin{aligned} (\mathrm{e}^{2x}y)'&=\mathrm{e}^{2x}\mathrm{e}^x=\mathrm{e}^{3x}&&\text{multiply through}\\ \mathrm{e}^{2x}y&=\frac13\mathrm{e}^{3x}+C&&\text{integrate}\\ y&=\frac13\mathrm{e}^x+C\mathrm{e}^{-2x}&&\text{divide by }\mathrm{e}^{2x} \end{aligned}

Why this step? The integrating factor depends on P(x)P(x), not on Q(x)Q(x). Changing the right side changes the integration, not the factor.

Answer. y=13ex+Ce−2xy=\frac13\mathrm{e}^x+C\mathrm{e}^{-2x}.

Check. y′+2y=(13+23)ex+(−2C+2C)e−2x=exy'+2y=(\frac13+\frac23)\mathrm{e}^x+(-2C+2C)\mathrm{e}^{-2x}=\mathrm{e}^x.

Keep the arbitrary constant until information about one particular curve is supplied. The next lesson uses conditions to choose it.

Watch out

Wrong: multiply the left side by e2x\mathrm{e}^{2x} but leave the right side unchanged. Right: (e2xy)′=e2xQ(x)(\mathrm{e}^{2x}y)'=\mathrm{e}^{2x}Q(x). The factor multiplies every term.

Use the same factor with a different exponential on the right.

Continue this lesson.

You’re previewing NovaLearn without an account. Create a free account to read the rest.

Create free account