vEigenVecs=lEigenVecs;

% Display results with classification
while(true)    
    h=figure;scatter(vEigenVecs(:,2),vEigenVecs(:,3),20,vEigenVecs(:,4),'filled');
    k = waitforbuttonpress;
    point1 = get(gca,'CurrentPoint');    % button down detected
    finalRect = rbbox;                   % return figure units
    point2 = get(gca,'CurrentPoint');    % button up detected
    point1 = point1(1,1:2);              % extract x and y
    point2 = point2(1,1:2);
    p1 = min(point1,point2);             % calculate locations
    offset = abs(point1-point2);         % and dimensions
    x = [p1(1) p1(1)+offset(1) p1(1)+offset(1) p1(1) p1(1)];
    y = [p1(2) p1(2) p1(2)+offset(2) p1(2)+offset(2) p1(2)];
    hold on
    axis manual
    plot(x,y,'r','linewidth',5)    
    lDocLabels = cLabels(find( (vEigenVecs(:,2)>p1(1)) & (vEigenVecs(:,3)>p1(2)) & (vEigenVecs(:,2)<p1(1)+offset(1)) & (vEigenVecs(:,3)<p1(2)+offset(2))),1);
    for lk = 1:length(lDocLabels)
        eval(sprintf('dos(''notepad c:\\sndata\\%d.'');',lDocLabels(lk)));
    end;   
    pause
    close(h);
end;