ConceptComplete

Higher-Order Constant Coefficient Equations

The methods for second-order linear equations extend naturally to higher-order equations with constant coefficients, providing systematic solution techniques for many physical and engineering problems.

Definitionn-th Order Linear Homogeneous ODE

An n-th order linear homogeneous differential equation with constant coefficients has the form:

any(n)+an1y(n1)++a1y+a0y=0a_ny^{(n)} + a_{n-1}y^{(n-1)} + \cdots + a_1y' + a_0y = 0

where a0,a1,,ana_0, a_1, \ldots, a_n are real constants with an0a_n \neq 0.

Characteristic Equation Method

The solution method parallels the second-order case. Seeking solutions of the form y=erxy = e^{rx} leads to the characteristic equation:

anrn+an1rn1++a1r+a0=0a_nr^n + a_{n-1}r^{n-1} + \cdots + a_1r + a_0 = 0

This is an n-th degree polynomial equation whose roots determine the general solution.

ExampleThird-Order Equation

Solve y6y+11y6y=0y''' - 6y'' + 11y' - 6y = 0.

Characteristic equation:

r36r2+11r6=0r^3 - 6r^2 + 11r - 6 = 0

Factoring: (r1)(r2)(r3)=0(r-1)(r-2)(r-3) = 0

Roots: r1=1r_1 = 1, r2=2r_2 = 2, r3=3r_3 = 3 (all distinct and real)

General solution:

y=c1ex+c2e2x+c3e3xy = c_1e^x + c_2e^{2x} + c_3e^{3x}

Cases for Roots

Case 1: Distinct Real Roots

If the characteristic equation has n distinct real roots r1,r2,,rnr_1, r_2, \ldots, r_n, the general solution is:

y=c1er1x+c2er2x++cnernxy = c_1e^{r_1 x} + c_2e^{r_2 x} + \cdots + c_ne^{r_n x}

Case 2: Repeated Real Roots

If a root rr has multiplicity mm, it contributes m linearly independent solutions:

erx,xerx,x2erx,,xm1erxe^{rx}, \quad xe^{rx}, \quad x^2e^{rx}, \quad \ldots, \quad x^{m-1}e^{rx}

ExampleRepeated Roots

Solve y(4)4y+6y4y+y=0y^{(4)} - 4y''' + 6y'' - 4y' + y = 0.

Notice this is (r1)4=0(r-1)^4 = 0, giving a root r=1r = 1 with multiplicity 4.

General solution:

y=(c1+c2x+c3x2+c4x3)exy = (c_1 + c_2x + c_3x^2 + c_4x^3)e^x

Case 3: Complex Roots

Complex roots occur in conjugate pairs. If r=α±iβr = \alpha \pm i\beta with multiplicity mm, they contribute:

eαxcosβx,xeαxcosβx,,xm1eαxcosβxe^{\alpha x}\cos\beta x, \quad xe^{\alpha x}\cos\beta x, \quad \ldots, \quad x^{m-1}e^{\alpha x}\cos\beta x

eαxsinβx,xeαxsinβx,,xm1eαxsinβxe^{\alpha x}\sin\beta x, \quad xe^{\alpha x}\sin\beta x, \quad \ldots, \quad x^{m-1}e^{\alpha x}\sin\beta x

ExampleMixed Roots

Solve y(4)+2y+y=0y^{(4)} + 2y'' + y = 0.

Characteristic equation: r4+2r2+1=(r2+1)2=0r^4 + 2r^2 + 1 = (r^2 + 1)^2 = 0

Roots: r=±ir = \pm i (each with multiplicity 2)

So α=0\alpha = 0, β=1\beta = 1, and we have:

y=(c1+c2x)cosx+(c3+c4x)sinxy = (c_1 + c_2x)\cos x + (c_3 + c_4x)\sin x

Remark

The number of linearly independent solutions always equals the order of the differential equation. For an n-th order equation, the general solution contains exactly n arbitrary constants.

ExampleMechanical Vibrations

A mass-spring system with three masses leads to a third-order equation. The characteristic roots determine whether the system exhibits:

  • Distinct real roots: Overdamped (no oscillation)
  • Complex roots: Oscillatory motion with exponential envelope
  • Repeated roots: Critical damping transitions

The mathematical structure directly reflects the physical behavior.

Higher-order equations with constant coefficients are completely solvable once the characteristic polynomial is factored, though factoring high-degree polynomials may require numerical methods.