clear all;
clc;

Data=dlmread('GM1_Class3.txt');
Labels=Data(:,end);
Data=Data(:,1:2);
NoClasses=max(unique(Labels));

[V,D]=MDA_mex(Data,Labels);
New_Data=Data*V;

Class1=Data(Labels==1,:);
Class2=Data(Labels==2,:);
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;
if(NoClasses==3)
    Class3=Data(Labels==3,:);
    plot(Class3(:,1),Class3(:,2),'+','MarkerSize',20,'Color',[0.9 0.4 0.2]);
end

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

quiver(MEAN(1),MEAN(2),V(1,1),V(2,1),10,'LineWidth',10,'Color',[0 0.5000 0.5000]);
hold all;
quiver(MEAN(1),MEAN(2),V(1,2),V(2,2),10,'LineWidth',5,'Color',[0 0.5000 0.5000]);