% Gauss-Seidel
A=hilb(5)
A=diag(1./diag(A))*A
U=triu(A,1)
L=tril(A,-1)
I=eye(5,5)
U+I+L
A

