Higher-Order Homogeneous DEs

Learning Objectives

  • Define the general form of a higher-order homogeneous linear differential equation.
  • Apply the Existence and Uniqueness Theorem and the Superposition Principle.
  • Determine linear independence of solutions using the Wronskian.
  • Solve higher-order linear DEs with constant coefficients using the auxiliary equation.
  • Solve Cauchy-Euler equations with variable coefficients.

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 of a Homogeneous Linear DE

The standard representation of an nn-th order linear homogeneous differential equation.

an(x)dnydxn+an−1(x)dn−1ydxn−1+⋯+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

Variables

SymbolDescriptionUnit
y(x)y(x)The unknown function of xxvaries
nnThe order of the highest derivativeunitless
ai(x)a_i(x)Coefficient functions depending only on xxvaries

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 with continuous coefficients and an(x)≠0a_n(x) \neq 0, 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.

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

A determinant used to verify if a set of solutions is linearly independent. If W≠0W \neq 0 for at least one point in the interval, the functions are linearly independent.

W(y1,y2)=det⁡(y1y2y1′y2′)=y1y2′−y2y1′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'

Variables

SymbolDescriptionUnit
WWThe Wronskian determinantvaries
yiy_iSolution functionsvaries
yi′y_i'First derivatives of the solution functionsvaries

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

Auxiliary (Characteristic) Equation

Substituting y=emxy = e^{mx} yields this algebraic equation used to find the roots mm.

am2+bm+c=0am^2 + bm + c = 0

Variables

SymbolDescriptionUnit
a,b,ca, b, cConstant coefficients from the differential equationvaries
mmRoots of the auxiliary equationunitless

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 (m1≠m2m_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. Interact with the spring-mass simulation below to explore how roots of the auxiliary equation dictate the motion of a homogeneous DE.

Spring-Mass System

my′′+cy′+ky=0m y'' + c y' + k y = 0
Underdamped (Complex Conjugate Roots)
1.0
0.5
4.0
Auxiliary Roots
r=−0.25±1.98ir = -0.25 \pm 1.98i
m
y(t) vs t

Cauchy-Euler Equation

A linear DE with variable coefficients where the power of xx matches the order of the derivative is called a Cauchy-Euler equation.

Substitution for Second-Order: For ax2y′′+bxy′+cy=0ax^2y'' + bxy' + cy = 0, we assume y=xmy = x^m.

Cauchy-Euler Auxiliary Equation

The resulting auxiliary equation after substituting y=xmy = x^m into a second-order Cauchy-Euler DE.

am2+(b−a)m+c=0am^2 + (b-a)m + c = 0

Variables

SymbolDescriptionUnit
a,b,ca, b, cConstant coefficients from the original Cauchy-Euler equationvaries
mmRoots of the Cauchy-Euler auxiliary equationunitless

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+c2xmln⁡∣x∣y = c_1 x^m + c_2 x^m \ln|x|
  3. Complex Conjugate Roots (m=α±βim = \alpha \pm \beta i): y=xα[c1cos⁡(βln⁡∣x∣)+c2sin⁡(βln⁡∣x∣)]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 continuous.
  • Superposition Principle states that linear combinations of linearly independent 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, ln⁡x\ln x for Cauchy-Euler) to maintain linear independence.
  • Cauchy-Euler: Use substitution y=xmy=x^m, leading to the auxiliary equation am2+(b−a)m+c=0am^2 + (b-a)m + c = 0.
PrevNext
Quiz Me