Matlab- numerical methods help !

Find the solution to the system of equations with coefficient matrices a shown below and vector answer b = [24; - 16; 0; 18; 5] by using Gauss-Seidel iterative methods and compare the answer with the solution of x = a\b. Compare the results by using the build-in solution as the tine solution in the formula for the relative error. a = [9 2 1 4 1 2 8 1 1 3 2 3 9 2 1 0 0 1 7 5 1 0 2 5 9] a = [5 1 1 4 1 2 8 3 1 3 2 3 8 2 1 1 1 1 7 5 1 1 2 5 9] a = [1 2 3 4 5 3 0 4 5 6 2 1 2 3 4 0 0 0 6 5 0 0 0 5 6]