function [U,value]=MXU(V,D,A,M);
value=eig(A);
null=0;
null_v=zeros(2*M,1);
for I=1:2*M
    for J=(I+1):(2*M)
        if abs(value(J))>abs(value(I))
            null=value(J);
            value(J)=value(I);
            value(I)=null;
            null_v=V(:,J);
            V(:,J)=V(:,I);
            V(:,I)=null_v;
        end
    end
end
U=V;