disp(' '), disp('question 2)a)')
M =[ 1 1 1;
0 -1 1;
-2 1 0]
[L,U,P]=lu(M) %avec pivot partiel
%rrefmovie(M)
disp(' '), disp(' '), disp('question 2)b)')
L=tril(sprand(10,10,0.2),-1)+eye(10);
U=triu(sprand(10,10,0.2),1)+eye(10);
disp('Le profil de LU (fig. 3) est réunion de ceux de L et U')
figure(1)
spy(L)
figure(2)
spy(U)
figure(3)
spy(L*U)