A system of first-order linear ODEs can be written in matrix form: y' = Ay + b(t), where y is a vector, A is a matrix, and b(t) is a forcing vector. This unified framework handles coupled equations and higher-order ODEs (converted to systems) via eigenvalue analysis.
When you studied first-order linear ODEs, each equation involved a single unknown function. Real systems — from predator-prey dynamics to electrical circuits — involve multiple coupled unknowns that influence each other's rates of change. A system of first-order linear ODEs captures exactly this: you have n unknown functions y₁, y₂, ..., yₙ, and each derivative yᵢ' depends linearly on all the others.
The key insight is that such a system can be written as a single matrix equation y' = Ay + b(t), where y is a vector of unknowns and A is a matrix whose entries encode the coupling. This is structurally identical to the scalar equation y' = ay + b(t) that you already know how to solve. The matrix A plays the same role as the scalar constant a — it governs how fast and in what direction the system evolves. This analogy is not just suggestive; the solution methods are direct generalizations.
One of the most powerful applications of this framework is converting higher-order ODEs into first-order systems. A second-order equation y'' = f(t, y, y') can be rewritten by introducing y₁ = y and y₂ = y', turning one second-order equation into two coupled first-order equations. This means every technique you develop for first-order systems — including eigenvalue methods and matrix exponentials — automatically applies to second-order and higher-order problems as well.
The solution to the homogeneous system y' = Ay (with b = 0) is built from the eigenvalues and eigenvectors of A. If A has eigenvalue λ with eigenvector v, then y(t) = e^(λt)v is a solution — you can verify this by differentiating. The general solution is a linear combination of such solutions, one per eigenvalue. The eigenvalue method, which you will study next, makes this systematic: it translates the ODE problem entirely into linear algebra, where you already have the tools to proceed.