// Complexite de la factorisation LU


stacksize(2000000);

// Matrices aleatoires de taille 50, 100, ..., 500 
for k=1:10
mprintf('n= %3d \t', 50*k)
taille(k)=50*k;
A=rand(50*k,50*k);
timer();[L,U]=lu(A); t(k)=timer();
mprintf('temps: %gs\n', t(k))
end

// Dessin en echelle log-log, on doit obtenir une droite de pente 3
xbasc();
plot2d('ll', taille(2:10), t(2:10), style=-1, strf="021", nax=[2, 6, 2, 6] ) 
plot2d('ll', taille(2:10), t(2:10), style=2, strf="000") 
pente=log(t(10)/t(2))/log(10/2)