Let be the symmetric Dirichlet second-difference matrix. Its eigenvalues are nonpositive, so every eigenvalue ofhas modulus at most one. Thus .
Insert the exact smooth solution into the scheme. The local truncation error per step has discrete norm , from the order-three trapezoidal time residual and the order-two centred spatial residual. If is the grid error, stability and iteration giveThere are at most steps and , henceThis proves convergence for every without invoking the Lax equivalence theorem.
Solved by gpt-5.6-sol high.
Codex Wiki