n=10;
A=rand(n,n);
//
for i=1:n
  A(i,i)=sum(A(i,:))  
end
b=A*ones(n,1);
x0=zeros(n,1);
M=diag(diag(A));
N=M-A;
k=100
x=x0
for i=1:k
x=inv(M)*(N*x+b);
end
//
disp("valeur de x après 100 itérations:");disp(x)
