t0 = clock;

[l,u] = lu(W);

D = prod(diag(u))

time = etime(clock,t0)