Gauss Seidel with Relaxation

It is a modified version of the Gauss-Seidel method that modifies the solved root with a weighted average for each iteration. The purpose of this is either to make the system convergent or enhances the convergence. It asks the user the coefficient and the RHD (Right-hand side) values matrix, and as well as the relaxation or weighing factor. It also lets the user choose which termination criterion to have, whether a %relative error or number of iterations. It checks if the system is diagonally dominant; if not, it would re-arrange the equations in the most diagonally dominant form possible, ensuring convergence. If the given system cannot be totally diagonally dominant, it alerts the user.