This function solves a linear system Ax=b using the Gauss-Seidel iteration method. The A-matrix is supplied in Compressed Row Storage (CRS) format. The iteration is terminated when either (1) the user-specified maximum number of iterations has been reached or (2) the norm of successive iterates is less than a user-specified epsilon.
Reference: (1) http://mathworld.wolfram.com/Gauss-SeidelMethod.html
Cheuk (2020). Compressed Row Storage (CRS) Gauss-Seidel Method (https://www.mathworks.com/matlabcentral/fileexchange/45350-compressed-row-storage-crs-gauss-seidel-method), MATLAB Central File Exchange. Retrieved .