Skip to content

Mathematics · Ch 10 — Ordinary Differential Equations

First Order Linear Differential Equations

10.7

First Order Linear Differential Equations

Definition. A first-order differential equation of the form

dydx+Py=Q\dfrac{dy}{dx}+Py=Q

— where PP and QQ are functions of xx only — is called linear: no product of yy with dydx\dfrac{dy}{dx} occurs, and yy together with its derivative occur only to the first degree.

Deriving the integrating factor. Consider first the associated homogeneous equation dydx+Py=0\dfrac{dy}{dx}+Py=0. Separating variables, dyy=−P dx\dfrac{dy}{y}=-P\,dx; integrating gives y e∫P dx=Cy\,e^{\int P\,dx}=C. Now examine what happens on differentiating the quantity y e∫P dxy\,e^{\int P\,dx} by the product rule:

ddx ⁣(y e∫P dx)=e∫P dxdydx+y P e∫P dx=e∫P dx(dydx+Py)=Q e∫P dx\dfrac{d}{dx}\!\left(y\,e^{\int P\,dx}\right)=e^{\int P\,dx}\dfrac{dy}{dx}+y\,P\,e^{\int P\,dx}=e^{\int P\,dx}\left(\dfrac{dy}{dx}+Py\right)=Q\,e^{\int P\,dx}

using the original (non-homogeneous) equation in the last step. So the left side is an exact derivative whenever yy solves the linear equation. Integrating both sides with respect to xx gives the general solution

y e∫P dx=∫Q e∫P dx dx+C.y\,e^{\int P\,dx}=\int Q\,e^{\int P\,dx}\,dx+C.

The factor e∫P dxe^{\int P\,dx} is called the integrating factor (I.F.) of the equation.

Remarks.

  1. The solution is often written compactly as y×(I.F.)=∫Q×(I.F.) dx+Cy\times(I.F.)=\displaystyle\int Q\times(I.F.)\,dx+C.
  2. In the integrating factor e∫P dxe^{\int P\,dx}, PP is the coefficient of yy provided the coefficient of dydx\dfrac{dy}{dx} has first been made unity — if the equation is not already in this normalised form, divide through by the coefficient of dydx\dfrac{dy}{dx} first.
  3. The mirror equation dxdy+Px=Q\dfrac{dx}{dy}+Px=Q (with P,QP,Q functions of yy only, no product of xx with dxdy\dfrac{dx}{dy}, both to the first degree) is solved the same way, with the roles of xx and yy exchanged: x e∫P dy=∫Q e∫P dy dy+Cx\,e^{\int P\,dy}=\displaystyle\int Q\,e^{\int P\,dy}\,dy+C.

Worked pattern (Example 10.22). dydx+2y=e−x\dfrac{dy}{dx}+2y=e^{-x}: here P=2, Q=e−xP=2,\ Q=e^{-x}, so I.F.=e2xI.F.=e^{2x}; the solution is ye2x=∫e−xe2xdx+C=ex+Cye^{2x}=\displaystyle\int e^{-x}e^{2x}dx+C=e^{x}+C, i.e. y=e−x+Ce−2xy=e^{-x}+Ce^{-2x}.

Worked pattern (Example 10.23, a trigonometric IF). After dividing through by xx to normalise, P=tan⁡xx−1xP=\dfrac{\tan x}{x}-\dfrac{1}{x}-type expressions arise; ∫P dx\int P\,dx collapses via log⁡(cos⁡x)\log(\cos x) and log⁡x\log x identities to give I.F.=1xcos⁡xI.F.=\dfrac{1}{x\cos x}, and the solution follows.

Worked pattern (Example 10.24, a power-of-trig IF). dydx+2ycot⁡x=3x2cosec⁡2x\dfrac{dy}{dx}+2y\cot x=3x^2\operatorname{cosec}^2x: P=2cot⁡xP=2\cot x, ∫P dx=2log⁡(sin⁡x)=log⁡(sin⁡2x)\int P\,dx=2\log(\sin x)=\log(\sin^2x), so I.F.=sin⁡2xI.F.=\sin^2x; the solution is ysin⁡2x=x3+Cy\sin^2x=x^3+C. …