Write a Matlab program for implementing Gauss-Seidel iterative method. Check whether the matrix is diagonally dominant at the start of the program. Do not use constant value for the loop (use tolerance to check convergence).