function kmeans load q1x.dat; a1=round(98*rand+1); a2=round(98*rand+1); miao1=[q1x(a1,1),q1x(a1,2)]; miao2=[q1x(a2,1),q1x(a2,2)]; c=zeros(99,1); sum1=zeros(1,2); sum2=zeros(1,2); for k=1:1 for i=1:99 temp1=norm(q1x(i,:)-miao1)^2; temp2=norm(q1x(i,:)-…