t=time();[l,u]=doolittle(Q); disp(time() - t)
0.2687
octave:8> t=time();[l,u]=crout(Q); disp(time() - t)
0.2562
octave:9> t=time();[l,u]=lu(Q); disp(time() - t)
0.1086