a = normrnd( 0, 1, 3, 3 ) [ L, U ] = lu(a)