Distinct Real Roots

If the characteristic equation has distinct real roots [latex]\lambda_1[/latex] and [latex]\lambda_2[/latex], then [latex]e^{\lambda_1 x}[/latex] and [latex]e^{\lambda_2 x}[/latex] are linearly independent solutions to Example “Classifying Second-Order Equations”, and the general solution is given by

[latex]\large{y(x)=c_1e^{\lambda_1 x}+c_2e^{\lambda_2 x}}[/latex],

where [latex]c_1[/latex] and [latex]c_2[/latex] are constants.

For example, the differential equation [latex]y^{\prime\prime}9y^\prime+14y=0[/latex] has the associated characteristic equation [latex]\lambda^2+9\lambda+14=0[/latex]. This factors into [latex](\lambda+2)(\lambda+7)=0[/latex], which has roots [latex]\lambda_1=-2[/latex] and [latex]\lambda_2=-7[/latex]. Therefore, the general solution to this differential equation is

[latex]\large{y(x)=c_1e^{-2}+c_2e^{-7x}}[/latex].

Single Repeated Real Root

Things are a little more complicated if the characteristic equation has a repeated real root, [latex]\lambda[/latex]. In this case, we know [latex]e^{\lambda x}[/latex] is a solution to Example “Classifying Second-Order Equations”, but it is only one solution and we need two linearly independent solutions to determine the general solution. We might be tempted to try a function of the form [latex]ke^{\lambda x}[/latex], where [latex]k[/latex] is some constant, but it would not be linearly independent of [latex]e^{\lambda x}[/latex]. Therefore, let’s try [latex]xe^{\lambda x}[/latex] as the second solution. First, note that by the quadratic formula,

[latex]\large{\lambda=\frac{-b\pm\sqrt{b^2-4ac}}{2a}}[/latex].

But, [latex]\lambda[/latex] is a repeated root, so [latex]b^2-4ac=0[/latex] and [latex]\lambda=\frac{-b}{2a}[/latex]. Thus, if [latex]y=xe^{\lambda x}[/latex], we have

[latex]\large{y^\prime=e^{\lambda x}+\lambda xe^{\lambda x}\text{ and }y^{\prime\prime}=2\lambda e^{\lambda x}+\lambda^2xe^{\lambda x}}[/latex].

Substituting these expressions into a linear second-order differential equation, we see that

[latex]\begin{aligned} ay^{\prime\prime}+by^\prime+cy&=a(2\lambda e^{\lambda x}+\lambda^2e^{\lambda x})+b(e^{\lambda x}+\lambda xe^{\lambda x})+cxe^{\lambda x} \\ &=xe^{\lambda x}(a\lambda^2+b\lambda+c)+e^{\lambda x}(2a\lambda+b) \\ &=xe^{\lambda x}(0)+e^{\lambda x}\left(2\left(\frac{-b}{2a}\right)+b\right) \\ &=0+e^{\lambda x}(0) \\ =0 \end{aligned}[/latex].

This shows that [latex]xe^{\lambda x}[/latex] is a solution to a linear second-order differential equation. Since [latex]e^{\lambda x}[/latex] and [latex]xe^{\lambda x}[/latex] are linearly independent, when the characteristic equation has a repeated root [latex]\lambda[/latex], the general solution to a linear second-order differential equation is given by

[latex]\large{y(x)=c_1e^{\lambda x}+c_2xe^{\lambda x}}[/latex],

where [latex]c_1[/latex] and [latex]c_2[/latex] are constants.

For example, the differential equation [latex]y^{\prime\prime}+12y^\prime+36y=0[/latex] has the associated characteristic equation [latex]\lambda^2+12\lambda+36=0[/latex]. This factors into [latex](\lambda+6)^2[/latex], which has a repeated root [latex]\lambda=-6[/latex]. Therefore, the general solution to this differential equation is

Complex Conjugate Roots

The third case we must consider is when [latex]b^2-4ac<0[/latex]. In this case, when we apply the quadratic formula, we are taking the square root of a negative number. We must use the imaginary number [latex]i=\sqrt{-1}[/latex] to find the roots, which take the form [latex]\lambda_1=\alpha+\beta i[/latex] and [latex]\lambda_2=\alpha-\beta i[/latex]. The complex number [latex]\alpha+\beta i[/latex] is called the conjugate of [latex]\alpha-\beta i[/latex]. Thus, we see that when [latex]b^2-4ac<0[/latex], the roots of our characteristic equation are always complex conjugates.

This creates a little bit of a problem for us. If we follow the same process we used for distinct real roots—using the roots of the characteristic equation as the coefficients in the exponents of exponential functions—we get the functions [latex]e^{(\alpha+\beta i)x}[/latex] and [latex]e^{(\alpha-\beta i)x}[/latex] as our solutions. However, there are problems with this approach. First, these functions take on complex (imaginary) values, and a complete discussion of such functions is beyond the scope of this text. Second, even if we were comfortable with complex-value functions, in this course we do not address the idea of a derivative for such functions. So, if possible, we’d like to find two linearly independent real-value solutions to the differential equation. For purposes of this development, we are going to manipulate and differentiate the functions [latex]e^{(\alpha+\beta i)x}[/latex] and [latex]e^{(\alpha-\beta i)x}[/latex] as if they were real-value functions. For these particular functions, this approach is valid mathematically, but be aware that there are other instances when complex-value functions do not follow the same rules as real-value functions. Those of you interested in a more in-depth discussion of complex-value functions should consult a complex analysis text.

Based on the roots [latex]\alpha\pm\beta i[/latex] of the characteristic equation, the functions [latex]e^{(\alpha+\beta i)x}[/latex] and [latex]e^{(\alpha-\beta i)x}[/latex] are linearly independent solutions to the differential equation. and the general solution is given by

[latex]\large{y(x)=c_1e^{(\alpha+\beta i)x}+c_2e^{(\alpha-\beta i)x}}[/latex].

Using some smart choices for [latex]c_1[/latex] and [latex]c_2[/latex], and a little bit of algebraic manipulation, we can find two linearly independent, real-value solutions to a linear second-order differential equation and express our general solution in those terms.

We encountered exponential functions with complex exponents earlier. One of the key tools we used to express these exponential functions in terms of sines and cosines was Euler’s formula, which tells us that

[latex]\large{e^{i\theta}=\cos\theta+i\sin\theta}[/latex]

for all real numbers [latex]\theta[/latex].

Going back to the general solution, we have

[latex]\large{\begin{aligned} y(x)&=c_1e^{(\alpha+\beta i)x}+c_2e^{(\alpha-\beta i)x} \\ &=c_1e^{\alpha x}e^{\beta ix}+c_2e^{\alpha x}e^{-\beta ix} \\ &=e^{\alpha x}\left(c_1e^{\beta ix}+c_2e^{-\beta ix}\right) \end{aligned}}[/latex].

Applying Euler’s formula together with the identities [latex]\cos(-x)=\cos x[/latex] and [latex]\sin(-x)=-\sin x[/latex], we get

[latex]\large{\begin{aligned} y(x)&=e^{\alpha x}[c_1(\cos\beta x+i\sin\beta x)+c_2(\cos(-\beta x)+i\sin(-\beta x))] \\ &=e^{\alpha x}[(c_1+c_2)\cos\beta x+(c_1-c_2)i\sin\beta x] \end{aligned}}[/latex].

Now, if we choose [latex]c_1=c_2=\frac{1}{2}[/latex], the second term is zero and we get

[latex]\large{y(x)=e^{\alpha x}\cos\beta x}[/latex]

as a real-value solution to a linear second-order differential equation. Similarly, if we choose [latex]c_1=-\frac{i}2[/latex] and [latex]c_2=\frac{i}2[/latex], the first term is zero and we get

[latex]\large{y(x)=e^{\alpha x}\sin\beta x}[/latex]

as a second, linearly independent, real-value solution to a linear second-order differential equation.

Based on this, we see that if the characteristic equation has complex conjugate roots [latex]\alpha\pm\beta i[/latex], then the general solution to a linear second-order differential equation is given by

[latex]\large{\begin{aligned} y(x)&=c_1e^{\alpha x}\cos\beta x+c_2e^{\alpha x}\cos\sin\beta x \\ &=e^{\alpha x}(c_1\cos\beta x+c_2\sin\beta x) \end{aligned}}[/latex].

where [latex]c_1[/latex] and [latex]c_2[/latex] are constants.

For example, the differential equation [latex]y^{\prime\prime}-2y^\prime+5y=0[/latex] has the associated characteristic equation [latex]\lambda^2-2\lambda+5=0[/latex]. By the quadratic formula, the roots of the characteristic equation are [latex]1\pm2i[/latex]. Therefore, the general solution to this differential equation is

[latex]\large{y(x)=e^x(c_1\cos2x+c_2\sin2x)}[/latex].