clear; clc;

X = xlsread('DATA1B');
X=X(1:13,1:4)
[n,p] = size(X);

k = 3;

idx = [3 9 12];
% 
C = X(idx,:)
% 
temp = zeros(n,1);

it=0
while 1,
    for i =1:n
        for l = 1:k        
            sum=0;
            for j = 1:p
                sum=sum+(X(i,j)-C(l,j))^2;
            end
            d(i,l)=sqrt(sum);
        end
    end
    [z,g]=min(d,[],2);
    if(g==temp)
        break;
    else
        temp = g;
    end
    for i=1:k
        f=find(g==i);
        if f
            C(i,:)=mean(X(find(g==i),:),1);
        end
    end
    it=it+1;
end


y = [X,g];
for i=1:n
    no(i) = norm(X(i,:));
end
y1 = [(1:n)', no',g];

k1 = X(find(g==1),:)
k2 = X(find(g==2),:)
k3 = X(find(g==3),:)
[n1,p1] = size(k1);
[n2,p2] = size(k2);
[n3,p3] = size(k3);
for i=1:n1
    norm1(i) = norm(k1(i,:));
end

x1 = [k1 norm1']

for i=1:n2
    norm2(i) = norm(k2(i,:));
end

x2 = [k2 norm2']

for i=1:n3
    norm3(i) = norm(k3(i,:));
end

x3 = [k3 norm3']

cov1 = cov(k1)
cov2 = cov(k2)
cov3 = cov(k3)
det1 = det(cov1)
det2 = det(cov2)
det3 = det(cov3)
c1 = y1(find(g==1),:);
c2 = y1(find(g==2),:);
c3 = y1(find(g==3),:);

plot(c1(:,1),norm1, 'r*', c2(:,1),norm2, 'b*', c3(:,1),norm3, 'g*')
hold on
grid on
legend({'Cluster 1','Cluster 2', 'Cluster 3'},'Location','northeast')
title('Norm Distance Plot');