function [V, ID] = fidpev(M,EV,TOL)

   V = []; ID = [];
   if nargin < 3
      TOL = 1e-14;
   end
   
   RKM = rank(M,TOL);
   M = M;
   nc = size(EV,2);
   
   
   i = 1; found = 0; more = 1;
   while (found == 0) & (more == 1)
      mh = [M EV(:,i)];
      if rank(mh,TOL) == (RKM+1)
         found = 1;
	 V = EV(:,i);
	 ID = i;
      else
         i = i+1;
         if i > nc
            more = 0;
         end
      end
   end
end
