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

%M Jacobi puntos; M=D

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

%M Jacobi bloques 1; M=D

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

%M Jacobi bloques 2; M=D

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

%{
ejecucion:

J =

         0   -1.0000         0
   -0.3333         0   -0.3333
         0   -0.2000         0


rho =

    0.6325


J =

         0   -1.0000         0
   -0.3571         0         0
    0.0714         0         0


rho =

    0.5976


J =

         0         0    0.5000
         0         0   -0.5000
         0   -0.2000         0


rho =

    0.3162
%}