Theorem · SuperpositionIf y₁ and y₂ solve a·y″ + b·y′ + c·y = 0, then C₁·y₁ + C₂·y₂ solves it for every pair of constants C₁ and C₂.
The rate of a sum is the sum of the rates, and a constant factor comes out of a rate, so for y = C₁·y₁ + C₂·y₂ the combination a·y″ + b·y′ + c·y splits into C₁·(a·y₁″ + b·y₁′ + c·y₁) + C₂·(a·y₂″ + b·y₂′ + c·y₂).
Each bracket is zero because y₁ and y₂ are solutions, so the whole is C₁·0 + C₂·0 = 0. The proof used two things: that the equation is linear, and that its right side is zero.
□
The second condition matters. The equation y″ + y = 1 is linear, and y = 1 solves it, since y″ = 0 and 0 + 1 = 1. The double y = 2 does not: 0 + 2 = 2, not 1. Superposition belongs to equations whose right side is zero.
Two solutions are not always enough
Take y₁ = e^(−t) and y₂ = 2·e^(−t) for y″ + 3y′ + 2y = 0. Both are solutions, and so is every combination C₁·y₁ + C₂·y₂ = (C₁ + 2·C₂)·e^(−t). The starting value is y(0) = C₁ + 2·C₂ and the starting rate is y′(0) = −C₁ − 2·C₂, which is minus the first. The start y(0) = 1, y′(0) = 0 would need C₁ + 2·C₂ to be 1 and 0 together.
Definition · The WronskianFor two functions y₁ and y₂, the Wronskian is W(t) = y₁(t)·y₂′(t) − y₁′(t)·y₂(t). Its value at t = 0 is the determinant of the coefficients of the pair of equations that fit C₁ and C₂.
For the pair above, y₁(0) = 1, y₂(0) = 2, y₁′(0) = −1 and y₂′(0) = −2. The determinant is 1·(−2) − 2·(−1) = 0, and the Wronskian is zero for every t, because the second function is a multiple of the first.
Theorem · What a nonzero Wronskian guaranteesIf y₁ and y₂ solve a·y″ + b·y′ + c·y = 0 and W(0) ≠ 0, then for every y(0) and y′(0) there is exactly one pair C₁, C₂ with C₁·y₁ + C₂·y₂ starting there.
The equations are C₁·y₁(0) + C₂·y₂(0) = y(0) and C₁·y₁′(0) + C₂·y₂′(0) = y′(0). Multiply the first by y₂′(0) and the second by y₂(0) and subtract, and the C₂ terms cancel: W(0)·C₁ = y₂′(0)·y(0) − y₂(0)·y′(0).
In the same way W(0)·C₂ = y₁(0)·y′(0) − y₁′(0)·y(0). If W(0) is not zero both can be divided by it, which gives one pair of constants and no other. That the combination is then the solution with that start is the fact that two starting values fix one solution, which the course uses and does not prove.
□
The four presets
The lab's four presets give the pair of solutions that the roots determine and print W(0). For y″ + 3y′ + 2y = 0 the pair is e^(−t) and e^(−2t), and the worked example below computes it. For y″ + 4y = 0 the pair is cos(2t) and sin(2t), and W is 2 for every t. For y″ + 4y′ + 4y = 0 the pair is e^(−2t) and t·e^(−2t), and W = e^(−4t), which is 1 at 0.
For two exponentials W(0) is the second root minus the first, with the larger root first, and a repeated root would make it zero: that is why the repeated case needed the second function t·e^(rt). For a complex pair W(0) is β. The fourth preset, y″ − y′ − y = 0, has the roots (1 ± √5)/2, which differ by √5, and the lab prints −√5 as an exact surd, not a rounded number; the sign is the order of the two roots. In all three computed cases W(t) = W(0)·e^(−(b/a)·t), so a Wronskian that starts nonzero never reaches zero; this is stated as a pattern, not proved.
The sign of W(0) carries no information. Swapping the roles of y₁ and y₂ reverses it, and what the guarantee needs is only that it is not zero.