for i = 1:100 S = sprand (1000, 1000, 1 / 1000); [L,U] = lu (S); fprintf (stderr, "%d\n", i); endfor