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 the differential equation, 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.

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 the differential equation, 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]\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]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.

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, when [latex]b^2-4ac<0[/latex], the roots of our characteristic equation are always complex conjugates.

This creates 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-valued functions, we do not address the idea of a derivative for such functions here.

If possible, we’d like to find two linearly independent real-valued 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-valued functions. For these particular functions, this approach is valid mathematically, but be aware that there are other instances when complex-valued functions do not follow the same rules as real-valued functions.

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]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.

Going back to the general solution, we have:

[latex]\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]\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]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]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]\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.

Summary of Results