clear all;
close all;
A = [2,10;2,5;8,4;5,8;7,5;6,4;1,2;4,9];
T = [1,4,7];
B = A(T,:);
for l = 1:2
cltr1 = [];
cltr2 = [];
cltr3 = [];
for i=1:8
    for j = 1:3
d(i,j) = ((B(j,1)-A(i,1))^2 + (B(j,2)-A(i,2))^2) ;
%d(i,j) = ((B(j,1)-A(i,1))^2 + (B(j,2)-A(i,2))^2) ;
%d(i,j) = B(j,:)*A(i,:)';
    end
end

for i = 1:8
    mn = min(d(i,:));
    jj = find(d(i,:) == mn);
    if (jj==1)
        cltr1 = [cltr1,i];
    else if (jj==2)
        cltr2 = [cltr2,i];
    else if (jj==3)
        cltr3 = [cltr3,i]; 
        end
        end
    end
end
m = 0;
n = 0;
for i = 1:size(cltr1,2)
m = m + A(cltr1(i),1);
n = n + A(cltr1(i),2);
end

c(1,1) = m / size(cltr1,2);
c(1,2) = n / size(cltr1,2);

m = 0;
n = 0;
for i = 1:size(cltr2,2)
m = m + A(cltr2(i),1);
n = n + A(cltr2(i),2);
end

c(2,1) = m / size(cltr2,2);
c(2,2) = n / size(cltr2,2);

m = 0;
n = 0;
for i = 1:size(cltr3,2)
m = m + A(cltr3(i),1);
n = n + A(cltr3(i),2);
end

c(3,1) = m / size(cltr3,2);
c(3,2) = n / size(cltr3,2);

B = c;

end