From here the unknown is called x, because the system treats position and velocity as a pair, and the equation under study is a·x″ + b·x′ + c·x = 0 with a not zero. Divide by a, and solve for the second rate.
x″ = −(c/a)·x − (b/a)·x′
v = x′
x′ = v
v′ = −(c/a)·x − (b/a)·v
The first line is the original equation, solved for its highest rate. The last two lines are a system of two first-order equations; the second is the first line with v written in place of x′. Nothing has been lost, since v is the rate of x by definition. The system is written as a matrix with one row for each equation; the words trace and determinant are those of Algebra's Systems and Matrices, and the two formulas in the claim below are all this lesson needs of them.
A = 0 1
−c/a −b/a
Theorem · The matrix carries the characteristic equationThe matrix A has trace τ = −b/a and determinant Δ = c/a, and the equation λ² − τ·λ + Δ = 0 is the characteristic equation a·r² + b·r + c = 0 divided by a.
The trace is 0 + (−b/a) and the determinant is 0·(−b/a) − 1·(−c/a), which are −b/a and c/a.
Then λ² − τ·λ + Δ = λ² + (b/a)·λ + c/a. Multiplying by a gives a·λ² + b·λ + c, which is the characteristic quadratic with λ in place of r. The eigenvalues of A, the numbers a later course defines as the roots of this quadratic, are the roots r of the equation.
□
The system and the equation agree
The first preset is x″ + 3x′ + 2x = 0, and the solution of the equation from the earlier lesson with x(0) = 1 and x′(0) = 0 is x = 2·e^(−t) − e^(−2t). Its rate is v = −2·e^(−t) + 2·e^(−2t). Both equations of the system hold exactly.
x′ = −2·e^(−t) + 2·e^(−2t) = v
v′ = 2·e^(−t) − 4·e^(−2t)
−2x − 3v = −4·e^(−t) + 2·e^(−2t) + 6·e^(−t) − 6·e^(−2t)
= 2·e^(−t) − 4·e^(−2t) = v′
The pair (x, v) is a solution of the system, and its first component is the solution of the equation. The pair also starts at (x(0), v(0)) = (1, 0), which is the two initial values of the equation written as a point.
The three presets
The lab takes the matrix as its input and prints the trace, the determinant, the eigenvalues and a name for the picture the system makes in the plane. The matrix for x″ + 3x′ + 2x = 0 has τ = −3 and Δ = 2, and its eigenvalues are −2 and −1. The matrix for x″ + 4x = 0 has τ = 0 and Δ = 4, eigenvalues ±2i. The matrix for x″ + 2x′ + 5x = 0 has τ = −2 and Δ = 5, eigenvalues −1 ± 2i. These are the roots from the earlier lessons, one for one.
The name the lab gives the picture follows the roots: negative real roots, imaginary roots and complex roots with a negative real part are named a stable node, a centre and a stable spiral. What the names mean is the subject of Systems and the Phase Plane; here they are a label on roots already understood.
The lab prints six more tiles that belong to that course: eigenvectors, a general solution written with them, the constants fitted to the start, and two Euler tiles that the next lesson reads. One of them is worth a look now. For the first preset the constants tile says C₁ = −1, C₂ = 2, while “Real Distinct Roots” fitted C₁ = 2, C₂ = −1 to the same start. The two labs number the roots in opposite orders, this one from the smaller eigenvalue up, so here C₁ = −1 goes with e^(−2t), and the first component is 2·e^(−t) − e^(−2t) either way.