Efficient Successive Over-Relaxation Method

This code solves the linear system Ax=b, where A is a symmetric banded matrix, using banded SOR.

Goal: This code solves the linear system Ax=b, where A is a symmetric banded
matrix, using banded SOR. A is first stored in compact storage mode
and next SOR is applied to the compact storage.
Input: This code accepts a symmetric banded matrix A, initial guess x0,
vector b, SOR parameter omega and tolerance epsilon.
Output: This code outputs the SOR-computed solution x, number of iteration
called Iteration and error called Error, which is the Euclidean
distance between kth iterate and k+1 iterate just when Error <
epsilon.