Higher-Order Homogeneous DEs

Higher-order linear differential equations are those where the highest derivative is greater than 1. A linear DE is homogeneous if the right-hand side (the term without yy or its derivatives) is zero.

General Form

an(x)dnydxn+an1(x)dn1ydxn1++a1(x)dydx+a0(x)y=0a_n(x) \frac{d^ny}{dx^n} + a_{n-1}(x) \frac{d^{n-1}y}{dx^{n-1}} + \dots + a_1(x) \frac{dy}{dx} + a_0(x)y = 0

Existence, Uniqueness, and Superposition

Before solving higher-order linear DEs, we rely on two fundamental theorems that guarantee our methods will work and that our solutions are complete.

Existence and Uniqueness Theorem

For a linear nn-th order IVP defined on an interval II:
an(x)y(n)++a1(x)y+a0(x)y=g(x)a_n(x)y^{(n)} + \dots + a_1(x)y' + a_0(x)y = g(x)
subject to y(x0)=y0,y(x0)=y1,,y(n1)(x0)=yn1y(x_0)=y_0, y'(x_0)=y_1, \dots, y^{(n-1)}(x_0)=y_{n-1}.
If all coefficients ai(x)a_i(x) and g(x)g(x) are continuous on II, and an(x)0a_n(x) \neq 0 for every xx in II, then there exists a unique solution y(x)y(x) to the IVP on that interval. This ensures that when we find a general solution and apply initial conditions, we have found the only possible particular solution.

Superposition Principle (Homogeneous)

If y1,y2,,yky_1, y_2, \dots, y_k are solutions to a homogeneous linear nn-th order DE on an interval II, then any linear combination of these solutions is also a solution:
y=c1y1+c2y2++ckyky = c_1y_1 + c_2y_2 + \dots + c_ky_k
where cic_i are arbitrary constants.

If the nn solutions are linearly independent, this linear combination forms the general solution.

Linear Independence

Before solving, we need to understand that the general solution is a linear combination of linearly independent solutions.

Linear Independence & Wronskian

A set of functions {y1,y2,,yn}\{y_1, y_2, \dots, y_n\} is linearly independent if none can be written as a linear combination of the others (i.e. c1y1+c2y2++cnyn=0c_1 y_1 + c_2 y_2 + \dots + c_n y_n = 0 is only true when all constants cic_i are 00).
Wronskian Test: Compute the determinant W(y1,,yn)W(y_1, \dots, y_n). If W0W \neq 0 for at least one point in the interval, the functions are linearly independent.
W(y1,y2)=det(y1y2y1y2)=y1y2y2y1\begin{aligned} W(y_1, y_2) = \det \begin{pmatrix} y_1 & y_2 \\ y_1' & y_2' \end{pmatrix} = y_1 y_2' - y_2 y_1' \end{aligned}

Constant Coefficients

For equations like ay+by+cy=0ay'' + by' + cy = 0 (where a,b,ca,b,c are constants), we assume a solution of the form y=emxy = e^{mx}. Substituting this yields the Auxiliary Equation (or Characteristic Equation):
am2+bm+c=0am^2 + bm + c = 0

Roots of Auxiliary Equation

Solve for mm using the quadratic formula. The form of the general solution depends on the nature of the roots:
  1. Distinct Real Roots (m1m2m_1 \neq m_2):
    y=c1em1x+c2em2xy = c_1 e^{m_1x} + c_2 e^{m_2x}
  2. Repeated Real Roots (m1=m2=mm_1 = m_2 = m):
    y=c1emx+c2xemxy = c_1 e^{mx} + c_2 x e^{mx}
  3. Complex Conjugate Roots (m=α±βim = \alpha \pm \beta i):
    y=eαx(c1cos(βx)+c2sin(βx))y = e^{\alpha x} (c_1 \cos(\beta x) + c_2 \sin(\beta x))

Spring-Mass System Visualization

Many higher-order homogeneous differential equations relate to physical systems, such as a mass on a spring. This interactive simulation illustrates how the mass moves under different conditions:

Spring-Mass-Damper System

Parameters

1 kg
0.5 Ns/m
10 N/m
1 m
System State
Underdamped
Crit. Damping = 6.32
m
Eq
Loading chart...

Cauchy-Euler Equation

A linear DE with variable coefficients of the form:
anxny(n)+an1xn1y(n1)++a1xy+a0y=0a_n x^n y^{(n)} + a_{n-1} x^{n-1} y^{(n-1)} + \dots + a_1 x y' + a_0 y = 0
is called a Cauchy-Euler equation. Notice the power of xx matches the order of the derivative.
Substitution for Second-Order: ax2y+bxy+cy=0ax^2y'' + bxy' + cy = 0. Assume y=xmy = x^m.

This leads to the auxiliary equation:
am(m1)+bm+c=0am(m-1) + bm + c = 0 am2+(ba)m+c=0am^2 + (b-a)m + c = 0

Roots and Solutions for Cauchy-Euler

  1. Distinct Real Roots: y=c1xm1+c2xm2y = c_1 x^{m_1} + c_2 x^{m_2}
  2. Repeated Real Roots: y=c1xm+c2xmlnxy = c_1 x^m + c_2 x^m \ln|x|
  3. Complex Conjugate Roots (m=α±βim = \alpha \pm \beta i): y=xα[c1cos(βlnx)+c2sin(βlnx)]y = x^{\alpha} [c_1 \cos(\beta \ln|x|) + c_2 \sin(\beta \ln|x|)]
Key Takeaways
  • Existence and Uniqueness guarantees that IVPs have exactly one solution if coefficients are well-behaved.
  • Superposition Principle states that linear combinations of solutions to a homogeneous linear DE form the general solution.
  • Constant Coefficients: Use the auxiliary equation am2+bm+c=0am^2+bm+c=0. Assume y=emxy = e^{mx}.
  • Complex Roots: Lead to sinusoidal solutions with exponential decay/growth (eαxcos(βx)e^{\alpha x}\cos(\beta x)).
  • Repeated Roots: Multiply by an independent variable (xx for constant coeffs, lnx\ln x for Cauchy-Euler) to maintain linear independence.
  • Cauchy-Euler: Use substitution y=xmy=x^m, leading to the auxiliary equation am2+(ba)m+c=0am^2 + (b-a)m + c = 0.