Data=dlmread('GM4.txt');
Class1=Data(Data(:,end)==1,:);
Class2=Data(Data(:,end)==2,:);

[V,D]=NDET_mex(Data);
[D,ind]=sort(diag(D),'descend');
V=V(:,ind);
New_Data=Data(:,1:end-1)*V;

plot(Class1(:,1),Class1(:,2),'*','MarkerSize',20,'Color',[0 0.4 0]);
hold all;
plot(Class2(:,1),Class2(:,2),'o','MarkerSize',20,'Color',[0 0 0.5]);
hold all;

MIN=-6;
MAX=6;
MEAN=mean(Data);

plot([MIN*V(1,1)+MEAN(1),MAX*V(1,1)+MEAN(1)],[MIN*V(2,1)+MEAN(2),MAX*V(2,1)+MEAN(2)],'-','LineWidth',25,'Color',[0.501 0 0]);
hold all;
plot([MIN*V(1,2)+MEAN(1),MAX*V(1,2)+MEAN(1)],[MIN*V(2,2)+MEAN(2),MAX*V(2,2)+MEAN(2)],'-','LineWidth',10,'Color',[0.501 0 0]);
hold all;