**Example
2.**** **Use successive over relaxation - SOR
iteration to solve the linear system .

Try 10 iterations. Compare the speed of convergence with
Jacobi and Gauss-Seidel iteration.

**Solution 2.**

Enter the 9×9 matrix **A** and the column vector
**B**.

Use 10 iterations of Jacobi iteration.

Use 10 iterations of Gauss-Seidel iteration.

Use 10 iterations of the SOR method with the parameter .

We can compare the approximations with the "true solution."

For this type of linear system the SOR method is known to give fast, accurate results.

**We are done.**

**Aside.**

The solution to the above linear system is
an approximate the solution of the
Laplace's equation

over
the rectangle

with the boundary values

, , for ,

and

and ,
for .

(Remark. There are singularities at the corners of the square and at
these points an average value along adjacent edges is assumed.)

The nine points are the interior points of a 5×5 grid and must be arranged in the following manner.

And the boundary values need to be filled in around these center points as follows.

Then we can plot the solution to the P. D. E.

(c) John H. Mathews 2004