Solution
= Solution
After multiplying the five-point equations by $-1$ if necessary, their matrix $A$ is symmetric positive definite. For any ordering write $A=D+L+L^T$, where $D=4I$. Gauss-Seidel method uses $M=D+L$ and has iteration matrix $-M^{-1}L^T$. In part (a) take $B=L^T$; then $A-B=M$ and
$$
A-B-B^T=D,
$$
which is positive definite. Therefore the iteration matrix has spectral radius below one for every ordering, proving convergence.