load HWMar26data;  %Use Matrix A for HW problem 1
clear B C

[idx,C,disterr]=kmeans(A,3);


% Plotting routines
idx1=find(idx==1);
idx2=find(idx==2);
idx3=find(idx==3);

plot(A(idx1,1),A(idx1,2),'r^');
hold on
plot(A(idx2,1),A(idx2,2),'b*');
plot(A(idx3,1),A(idx3,2),'ko');
hold off