function [V, D] = eigOrder(X)
%EIGORDER Computes eigen values and eigen vectors and order them.
%

% Tomo Eguchi
% 8 February 2002


[nrX, ncX] = size(X);

% find eigenvalues and eigenvectors:
[V1, D1] = eig(X);

% Order eigenvalues and eigenvectors:
orderEigenD = [1:nrX];
eigenDVec = [sum(D1); orderEigenD];
eigenDOrdered = fliplr(rot90(sortrows(eigenDVec')));
D = diag(eigenDOrdered(2,:));
V = reorder(V1, eigenDOrdered(1,:));