MATLAB code of the Gauss-Seidel method

Description

eps = 1e-8 ;                x0 = [0;0;0] ;

 A = [2 -1 1; -2 5 -1; 1 -2 4] ;                           b = [-1;1;3] ;

Final result : 

x =

-1.0000
0
1.0000

Iters =

11

A = [-5 1 -2; 1 6 3; 2 -1 -4] ;                                   b = [13; 1; -1] ;

Final result : 

x =
-2.0000
1.0000
-1.0000

Iters =
15