Part 2: Implement the Gauss-Seidel algorithm to solve Ar=b. INPUT: An n x n matrix A, a right-hand side b, an initial guess x0, tolerance TOL and maximum number of iterations Nmax OUTPUT: An approximate solution x or a message of failure STEP 1: set x - x0; STEP 2: For k = 1,2,...,Nmax do STEPS 3-7 STEP 3: Set xold = x; STEP 4: For i =1,2,...,n do STEPS STEP 5: Set SUM = -sum( Ali,j)*x(j), j-1,2,...,1-1) - sum( A(i, j)*xold(j), j=i+1,i+2,...,n) + b(i) STEP 6: Set x(i) = SUM/A(i,i) STEP 7: If max(abs(x-xold))