An nth-order linear ODE has the form y^(n) + a_{n-1}y^(n-1) + ... + a₁y' + a₀y = f(x). The same principles apply: combine n linearly independent homogeneous solutions and add a particular solution. For constant coefficients, the characteristic equation becomes a polynomial of degree n. Higher-order equations arise naturally when modeling complex mechanical and electrical systems.
From your work with undetermined coefficients and variation of parameters, you know the structure of second-order linear ODEs: find two linearly independent solutions to the homogeneous equation, add a particular solution for the forcing term, and the general solution is their combination. Higher-order equations extend this pattern without introducing any fundamentally new ideas — the dimension of the solution space simply grows to match the order.
For an nth-order linear constant-coefficient ODE, the characteristic equation is still found by substituting y = e^(rx) and simplifying. Where a second-order equation gives a quadratic r² + ar + b = 0, an nth-order equation gives a degree-n polynomial: rⁿ + a_{n-1}r^(n-1) + ··· + a₁r + a₀ = 0. The roots of this polynomial — real, complex, or repeated — determine the n linearly independent homogeneous solutions by the same rules you already know. A real distinct root r gives a solution e^(rx). A pair of complex conjugate roots α ± βi gives the pair e^(αx)cos(βx) and e^(αx)sin(βx). A root of multiplicity k gives k solutions: e^(rx), xe^(rx), x²e^(rx), ..., x^(k-1)e^(rx).
The critical requirement is that you collect exactly n linearly independent solutions — one for each root counting multiplicity, and the multiplicity rule ensures this count works out. For example, a fourth-order ODE with characteristic roots r = 1 (simple), r = 1 (so double root at 1 overall... wait, let me redo this) — if the characteristic polynomial is (r−2)²(r² + 1) = 0, the roots are r = 2 (double), r = i, r = −i. This yields four homogeneous solutions: e^(2x), xe^(2x), cos(x), sin(x). The Wronskian test confirms linear independence, though in practice, roots from different factors of the characteristic polynomial are always independent.
Physical applications make the structure concrete. A fourth-order beam equation, or a coupled mass-spring system with two masses, naturally produces fourth-order or coupled second-order ODEs. The same spring-damper analogy applies: characteristic roots with negative real parts give decaying solutions (stable systems), roots on the imaginary axis give pure oscillation, and positive real parts signal instability. Higher-order equations let you model systems with more degrees of freedom — more masses, more coupled components — while the mathematical machinery remains exactly what you already know, scaled up by the degree.