This final project discuss a preconditioned Gauss-Seidel method to solve a system of linear equation Ax = b by A which is a strictly diagonally dominant Z-matrix. Preconditioning matrix to be used is P = (I + U)−1, where I is an identity matrix and U is a strictly upper triangular matrix. Using Neumann’s expansion to approximate P, we show that the preconditioning matrix is equivalent …