System resolution function with Gauss-Seidel method

System resolution function with Gauss-Seidel method

System resolution function with Gauss-Seidel method
Seidel function is a function created to solve linear systems that
have matrix of scattered coefficients.
This function receives in input the matrix of coefficients A even in
compact notation, vector of terms known, tolerance and maximum
number of iterations so you can write the command
[ sol, res, err, nit ] = seidel(A,b, TOL, NMAX ) or
[ sol, res, err, nit ] = seidel(A,b)
The function returns the vector of solutions sol, the
residue vector res, the relative error err and the number of iterations
made nit and estimation of convergence speed rho
The function stops if it reaches the desired tolerance
or if the number of iterations has reached its maximum,
in this case the function will give a warning
TOL and NMAX input parameters are optional.
CREATED BY GIOVANNI D'AVANZO Università della Campania "Luigi Vanvitelli"
A.A. 2016/2017