Please provide results screen capture and python codes, Thanks

(2) Write a Python function that solves a system of linear equations using the Gauss-Seidel method. 1- Your function should take the matrix A, vector b, and e as input. 2-Then, it should check if matrix A is diagonally dominant. If not, it should print "matrix A is not diagonally dominant". 3-If matrix A is diagonally dominant, the function should run the Gauss-Seidel iterations until all the elements of the absolute value of the difference vector are less than e. 4-It should then return the solution of the system of linear equations