First-Order Linear Differential Equations

[math] \newcommand{\ex}[1]{\item } \newcommand{\sx}{\item} \newcommand{\x}{\sx} \newcommand{\sxlab}[1]{} \newcommand{\xlab}{\sxlab} \newcommand{\prov}[1] {\quad #1} \newcommand{\provx}[1] {\quad \mbox{#1}} \newcommand{\intext}[1]{\quad \mbox{#1} \quad} \newcommand{\R}{\mathrm{\bf R}} \newcommand{\Q}{\mathrm{\bf Q}} \newcommand{\Z}{\mathrm{\bf Z}} \newcommand{\C}{\mathrm{\bf C}} \newcommand{\dt}{\textbf} \newcommand{\goesto}{\rightarrow} \newcommand{\ddxof}[1]{\frac{d #1}{d x}} \newcommand{\ddx}{\frac{d}{dx}} \newcommand{\ddt}{\frac{d}{dt}} \newcommand{\dydx}{\ddxof y} \newcommand{\nxder}[3]{\frac{d^{#1}{#2}}{d{#3}^{#1}}} \newcommand{\deriv}[2]{\frac{d^{#1}{#2}}{dx^{#1}}} \newcommand{\dist}{\mathrm{distance}} \newcommand{\arccot}{\mathrm{arccot\:}} \newcommand{\arccsc}{\mathrm{arccsc\:}} \newcommand{\arcsec}{\mathrm{arcsec\:}} \newcommand{\arctanh}{\mathrm{arctanh\:}} \newcommand{\arcsinh}{\mathrm{arcsinh\:}} \newcommand{\arccosh}{\mathrm{arccosh\:}} \newcommand{\sech}{\mathrm{sech\:}} \newcommand{\csch}{\mathrm{csch\:}} \newcommand{\conj}[1]{\overline{#1}} \newcommand{\mathds}{\mathbb} [/math]

A differential equation which can be written in the form

[[math]] a_1(x) \frac{dy}{dx} + a_0(x)y = b(x), [[/math]]

where [math]a_1, a_0[/math], and [math]b[/math] are given functions of [math]x[/math] and where [math]a_1[/math] is not the zero function, is a first-order linear ditTerential equation. In this section we shall show how to obtain the general solution of equations of this type. Since [math]a_1[/math] is by assumption not the zero function, we can divide both sides of the above equation by [math]a_1(x)[/math]. Setting [math]\frac{a_0(x)}{a_1(x)} = P(x)[/math] and [math]\frac{b(x)}{a_1(x)} = Q(x)[/math], we therefore obtain the differential equation

[[math]] \begin{equation} \frac{dy}{dx} + P(x)y = Q(x), \label{eq11.2.1} \end{equation} [[/math]]

which is the form we shall use in deriving the solution. We shall assume that the functions [math]P[/math] and [math]Q[/math] are continuous, thus assuring ourselves that they have antiderivatives. Let us suppose that the function [math]y = f(x)[/math] is a solution to the differential equation (1). We shall derive a formula which expresses [math]y[/math] in terms of [math]P[/math] and [math]Q[/math] and a constant [math]c[/math] of integration. Conversely, it will be a simple matter to verify that any function [math]y[/math] defined by this formula is a solution to (1). Hence the formula gives the general solution to the differential equation. The derivative of the product of [math]y[/math] and a function [math]\varphi[/math] is given by

[[math]] \begin{equation} \frac{d}{dx} (\varphi(x)y) = \varphi(x) + \varphi'(x)y. \label{eq11.2.2} \end{equation} [[/math]]

Note that the first term on the right has [math]\frac{dy}{dx}[/math] as a factor and the second has [math]y[/math] as a factor, and that the same is true of the two terms on the left side of equation (1). This fact suggests seeking a function [math]\varphi[/math] which has the property that, if both sides of (1) are multiplied by [math]\varphi(x)[/math], then the left side of the resulting equation is the derivative of the product [math]\varphi(x)y[/math]. If both sides of (1) are multiplied by an arbitrary [math]\varphi(x)[/math], the result is

[[math]] \begin{equation} \varphi(x) \frac{dy}{dx} + \varphi(x)P(x)y = \varphi(x)Q(x) . \label{eq11.2.3} \end{equation} [[/math]]

Comparison of this equation with (2) shows that its left side is equal to [math]\frac{d}{dx} (\varphi(x)y)[/math] provided [math]\varphi(x)P(x)y = \varphi'(x)y[/math], which will in turn be true provided

[[math]] \begin{equation} \varphi(x)P(x) = \varphi'(x). \label{eq11.2.4} \end{equation} [[/math]]

However, it is easy to find a function [math]\varphi[/math] which satisfies (4), since, as a differential equation with [math]\varphi[/math] the unknown function, it is separable. Solving it, we obtain

[[math]] \frac{\varphi'(x)}{\varphi(x)} = P(x) [[/math]]

Whence

[[math]] \int \frac{\varphi'(x)}{\varphi(x)} dx = \int P(x) dx , [[/math]]

which implies

[[math]] \ln |\varphi(x)| = \int P(x) dx + c, [[/math]]

and so

[[math]] |\varphi(x)| = e^{\int P(x)dx+c} . [[/math]]

Since we are only seeking a solution to (4), and not the most general form of the solution, we may assume that [math]\varphi(x)[/math] is positive and also ignore the constant of integration. We conclude that if

[[math]] \begin{equation} \varphi(x) = e^{\int P(x) dx}, \label{eq11.2.5} \end{equation} [[/math]]

then the left side of eguation (3) is equal to [math]\frac{d}{dx} (\varphi(x)y)[/math]. With (5), equation (3) therefore becomes

[[math]] \frac{d}{dx} (\varphi(x)y) = \varphi(x)Q(x). [[/math]]

integration yields

[[math]] \varphi(x)y = \int \varphi(x) Q(x) dx + c, [[/math]]

and so

[[math]] y = \frac{1}{\varphi(x)} \Big[ \int \varphi(x) Q(x) dx + c \Big] , [[/math]]

for some real number [math]c[/math]. Replacing [math]\varphi(x)[/math] by [math]e^{\int P(x)dx}[/math], we obtain the promised formula:

Theorem

[[math]] y = e^{- \int P(x)dx} \Big[\int e^{\int P(x)dx} Q(x) dx + c \Big] . [[/math]]

Suppose next that [math]c[/math] is an arbitrary constant and that the function [math]y[/math] is dcfined by (2.1). Then

[[math]] ye^{\int P(x)dx} = \int e^{\int P(x)dx} Q(x) + c. [[/math]]

The derivative of the right side of this equation is [math]e^{\int P(x)dx}Q(x)[/math] and that of the left side is

[[math]] \frac{dy}{dx} e^{\int P(x)dx} + yP(x)e^{\int P(x)dx} . [[/math]]

Hence

[[math]] e^{\int P(x)dx} \Big[\frac{dy}{dx} + P(x)y \Big] = e^{\int P(x)dx} Q(x), [[/math]]

which implies at once that [math]\frac{dy}{dx} + P(x)y = Q(x)[/math]; i.e., [math]y[/math] is a solution to (1). We conclude that formula (2.1) gives the general solution to the differential equation (1). We strongly recommend that no one memorize (2.1). The important fact to remember is that, if the first-order linear differential equation [math]\frac{dy}{dx} + P(x)y = Q(x)[/math] is multiplied through by [math]\varphi(x) = e^{\int P(x)dx}[/math], then the left side of the resulting equation is equal to the derivative of the product [math]\varphi(x)y[/math]. Consequently, the new equation can be integrated to give

[[math]] \varphi(x)y = \int \varphi(x) Q(x) dx + c, [[/math]]

which can then be solved for [math]y[/math]. This function [math]\varphi(x) = e^{\int P(x)dx}[/math], which enables us to change a seemingly nonintegrable sum into the derivative of a product by multiplication, is called an integrating factor.

Example Solve the differential equation

[[math]] x^2 \frac{dy}{dx} - 3xy - 2x^2 = 4x^4 . [[/math]]

To put this in the form of (1), we add [math]2x^2[/math] to both sides and then divide by [math]x^2[/math]. The result is

[[math]] \begin{equation} \frac{dy}{dx} - \frac{3}{x} = 4x^2 + 2, \label{eq11.2.6} \end{equation} [[/math]]

where [math]P(x) = -\frac{3}{x}[/math] and [math]Q(x) = 4x^2 + 2[/math]. An antiderivative of [math]P(x)[/math] is given by

[[math]] \int P(x) dx = \int -\frac{3}{x} dx = -3 \ln |x| = \ln |x^{-3}|, [[/math]]

and it follows that the function

[[math]] \varphi(x) = e^{\int P(x)dx} = e^{\ln |x^{-3}|} = |x^{-3}| [[/math]]

is an integrating factor. Equation (4) shows that if [math]\varphi(x)[/math] is an integrating factor, then so also is [math]-\varphi(x)[/math]. Hence we may drop the absolute values and write simply

[[math]] \varphi(x) = x^{-3}. [[/math]]

Multiplying both sides of (6) by [math]x^{-3}[/math], we obtain

[[math]] x^{-3} \frac{dy}{dx} - 3x^{-4}y = 4x^{-1} + 2x^{-3} , [[/math]]

It is easy to see that the left side of this equation is equal to [math]\frac{d}{dx}(x^{-3}y)[/math]. Hence

[[math]] \frac{d}{dx} (x^{-3} y) = 4x^{-1} + 2x^{-3}, [[/math]]

and so

[[math]] \begin{eqnarray*} x^{-3}y &=& \int (4x^{-1} + 2x^{-3}) dx + c\\ &=& 4 \ln |x| + 2 \frac{x^{-2}}{-2} + c \\ &=& 4 \ln |x| - \frac{1}{x^2} + c, \end{eqnarray*} [[/math]]

where [math]c[/math] is an arbitrary constant. It follows that

[[math]] y = 4x^{3} \ln |x| - x + cx^{3} [[/math]]

is the general solution.

Example Find the general solution of the differential equation

[[math]] \frac{dy}{dx} + 3y = 2 \sin x. [[/math]]

Note that this is a first-order linear differential equation with constant coefficients, but that it is not homogeneous, because the right side is not the zero function. In this example we have [math]P(x) = 3[/math] and [math]Q(x) = 2 \sin x[/math]. Hence

[[math]] \int P(x)dx = \int 3 dx = 3x , [[/math]]

and an integrating factor is

[[math]] \varphi(x) = e^{\int P(x)dx} = e^{3x} . [[/math]]

lt follows that

[[math]] \frac{d}{dx} (e^{3x} y) = 2e^{3x} \sin x, [[/math]]

and so

[[math]] \begin{equation} e^{3x} y = \int 2e^{3x} \sin x dx + c. \label{eq11.2.7} \end{equation} [[/math]]

To evaluate [math]\int 2e^{3x} \sin xdx = 2\int e^{3x} \sin x dx[/math], we use integration by parts twice:

[[math]] \begin{eqnarray*} \int e^{3x} \sin x dx &=& - \int e^{3x} d \cos x \\ &=& -e^{3x} - \cos x + \int \cos x de^{3x} \\ &=& -e^{3x} \cos x + 3 \int e^{3x} \cos x dx. \\ \int e^{3x} \cos x dx &=& \int e^{3x} d \sin x \\ &=& e^{3x} \sin x - \int \sin x de^{3x}\\ &=& e^{3x} \sin x - 3\int e^{3x} \sin x dx. \end{eqnarray*} [[/math]]

Combining these results, we get

[[math]] \int e^{3x} \sin x dx = - e^{3x} \cos x + 3e^{3x} \sin x - 9 \int e^{3x} \sin x dx, [[/math]]

whence

[[math]] 10 \int e^{3x} \sin x dx = e^{3x}(3 \sin x - \cos x), [[/math]]

and so

[[math]] 2\int e^{3x} \sin x dx = \frac{e^{3x}}{5} (3 \sin x - \cos x). [[/math]]

Returning to (7), we have

[[math]] e^{3x}y = \frac{e^{3x}}{5} (3 \sin x - \cos x) + c, [[/math]]

and consequently the general solution of the differential equation is given by

[[math]] \begin{equation} y = \frac{1}{5} (3 \sin x - \cos x) + ce^{-3x} , \label{eq11.2.8} \end{equation} [[/math]]

where [math]c[/math] is an arbitrary constant.

Note that the above solution (8) of the differential equation of Example 2 is the sum of two terms. The second, which is [math]ce^{-3x}[/math], is thc general solution of the homogeneous differential equation [math]\frac{dy}{dx} + 3y = 0[/math]. The first term, [math]\frac{1}{5}(3 \sin x - \cos x)[/math], is one particular solution of the nonhomogeneous differential equation [math]\frac{dy}{dx} + 3y = 2 \sin x[/math]. As we shall see, this situation is typical of the solutions of linear differential equations.

General references

Doyle, Peter G. (2008). "Crowell and Slesnick's Calculus with Analytic Geometry" (PDF). Retrieved Oct 29, 2024.