2.6 - First-order linear differential equations

2.6 - First-order linear differential equations

Some differential equations cannot be separated, yet become a single product derivative after multiplication by a suitable function. You will identify these linear equations, derive an integrating factor and use an initial condition to select a solution on a valid interval.

Identifying a linear equation

How can you solve xy′+2y=x3xy'+2y=x^3 when the terms involving xx and yy will not separate? The first step is to make the coefficient of the derivative equal to one.

Here and below, y′=dy/dxy'=\mathrm{d}y/\mathrm{d}x. Dividing by xx, on an interval where x≠0x\ne0, gives y′+(2/x)y=x2y'+(2/x)y=x^2. The unknown function and its derivative occur to the first power, with coefficients depending only on xx.

First-order linear equation

A first-order linear differential equation can be written

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

The functions PP and QQ depend only on the independent variable. For example, y′+xy=sin⁡xy'+xy=\sin x is linear; y′+xy2=sin⁡xy'+xy^2=\sin x is not linear in yy.

“Linear” describes the dependence on the unknown function. The coefficient P(x)P(x) itself need not be a linear function of xx.

Worked example 1: making the derivative coefficient one

Question. Write xy′+2y=x3xy'+2y=x^3, for x>0x>0, in standard linear form and identify P(x)P(x) and Q(x)Q(x).

Plan. Divide every term by the nonzero coefficient of y′y'.

y′+2xy=x2.y'+\frac2x y=x^2.

The coefficient of y′y' is now one, so the coefficient of yy is the required PP.

Answer. P(x)=2/xP(x)=2/x and Q(x)=x2Q(x)=x^2, on x>0x>0.

Check. Multiplication by xx restores all three original terms, including the right-hand side x3x^3.

Change only the sign of the term in yy. That sign will change the integrating factor in the next part.

Worked example 2: retaining the signed coefficient

Question. Write xy′−2y=x3xy'-2y=x^3, for x>0x>0, in standard linear form and identify P(x)P(x) and Q(x)Q(x).

Plan. Divide through by xx, retaining the minus sign.

y′−2xy=x2.y'-\frac2x y=x^2.

Comparing with y′+Py=Qy'+Py=Q means that the sign belongs to PP.

Answer. P(x)=−2/xP(x)=-2/x and Q(x)=x2Q(x)=x^2.

Check. At x=2x=2, the normalised equation is y′−y=4y'-y=4. Doubling it gives 2y′−2y=82y'-2y=8, exactly the original equation at that value of xx.

The interval is part of the problem. Dividing by xx cannot establish a solution across x=0x=0; a solution on a negative interval can be treated separately.

Watch out

In a(x)y′+b(x)y=c(x)a(x)y'+b(x)y=c(x), the standard-form coefficient is P=b/aP=b/a, not bb. Divide the right-hand side by aa as well.

Identify both coefficients in this negative-domain example.

The integrating factor

The product rule explains the multiplier we need: (x2y)′=x2y′+2xy(x^2y)'=x^2y'+2xy, so multiplying y′+(2/x)y=x2y'+(2/x)y=x^2 by x2x^2 creates a product derivative.

For a general multiplier μ(x)\mu(x), the product rule gives

(μy)′=μy′+μ′y.(\mu y)'=\mu y'+\mu' y.

After multiplying the standard equation by μ\mu, we need μ′=Pμ\mu'=P\mu. Separating this auxiliary equation gives a convenient choice μ=exp⁡(∫P dx)\mu=\exp(\int P\,\mathrm{d}x).

Integrating factor

For y′+Py=Qy'+Py=Q on an interval where PP and QQ are continuous, choose

μ=e∫P dx,(μy)′=μQ.\mu=e^{\int P\,\mathrm{d}x},\qquad (\mu y)'=\mu Q.

Then

y=∫μQ dx+Cμ.y=\frac{\int\mu Q\,\mathrm{d}x+C}{\mu}.

A nonzero constant multiple of μ\mu works equally well, so no arbitrary constant is needed in the chosen integrating factor. Not in the formula booklet: learn it.

Multiply the whole equation, integrate the right-hand side, then divide the whole result by the integrating factor. This final division also applies to the constant.

Worked example 3: a positive power as integrating factor

Question. Find the general solution of xy′+2y=x3xy'+2y=x^3, for x>0x>0.

Plan. Normalise, multiply by x2x^2, then integrate the product derivative.

y′+2xy=x2μ=e2ln⁡x=x2(x2y)′=x4x2y=x55+C.\begin{aligned} y'+\frac2x y&=x^2 \\ \mu&=e^{2\ln x}=x^2 \\ (x^2y)'&=x^4 \\ x^2y&=\frac{x^5}{5}+C. \end{aligned}

The last line is an equation for x2yx^2y, so divide both terms on the right by x2x^2.

Answer. y=x3/5+C/x2y=x^3/5+C/x^2, for x>0x>0.

Check. Differentiating gives y′=3x2/5−2C/x3y'=3x^2/5-2C/x^3. Thus xy′+2y=(3x3/5−2C/x2)+(2x3/5+2C/x2)=x3xy'+2y=(3x^3/5-2C/x^2)+(2x^3/5+2C/x^2)=x^3.

Continue this lesson.

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

Create free account