A=[1 1 0; 1 3 1; 0 1 5];

%M Gauss-Seidel puntos; M=D-E

M=[1 0 0; 1 3 0; 0 1 5]; N=M-A; 
L1=inv(M)*N 
rho=max(abs(eig(L1)))

%M Gauss-Seidel bloques 1; M=D-E

M=[1 0 0; 1 3 1; 0 1 5]; N=M-A; 
L1=inv(M)*N 
rho=max(abs(eig(L1)))

%M Gauss-Seidel bloques 2; M=D-E

M=[1 1 0; 1 3 0; 0 1 5]; N=M-A; 
L1=inv(M)*N
rho=max(abs(eig(L1)))

%{
ejecucion:

ejemplo_m_gauss_seidel

L1 =

         0   -1.0000         0
         0    0.3333   -0.3333
         0   -0.0667    0.0667


rho =

    0.4000


L1 =

         0   -1.0000         0
         0    0.3571         0
         0   -0.0714         0


rho =

    0.3571


L1 =

         0         0    0.5000
         0         0   -0.5000
         0         0    0.1000


rho =

    0.1000
%}