clear;
x= [1.4,1.5
    3.5,5.7
    2.1,8.9
    1.8,2.7
    2.7,1.6
    3.6,6.8
    2.3,4.6
    1.5,5.6
    9.2,8.9
    2.8,6.5];
u=[1,0
    0,1
    1,0
    0,1
    0,1
    0,1
    0,1
    1,0
    1,0
    1,0];
c=[0,0
   0,0];
m=2;
for k=1:1000
    for j=1:2
        up=[0,0];
        down=[0,0];
        for i=1:10
            up=up+(u(i,j)^m*x(i,:));
            down=down+u(i,j)^m;
        end
        c(j,1)=up(1)/down(1);
        c(j,2)=up(2)/down(2);
    end
    oldu=u;
    for i=1:10
        for j=1:2
            u(i,j)=1/((sqrt(sum((x(i,:)-c(j,:)).^2))/sqrt(sum((x(i,:)-c(1,:)).^2)))^(2/(m-1))+(sqrt(sum((x(i,:)-c(j,:)).^2))/sqrt(sum((x(i,:)-c(2,:)).^2)))^(2/(m-1)));      
        end
    end
    if sqrt(sum((u-oldu).^2))<0.001
        c(1,:)
        c(2,:)
        return;
    end
end