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