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)