w=sum(b')*inf+1;
w=sum(b')+inf;
p={1};
q=1;
while q
r={};
s=[];
for i=find(q)
y=p{i};
z=y(end);
g=q(i);
if g<w(z)
w(z)=g;
f{z}=y;
end
for j=find(a(z,:))
if ~any(j==y)
r=[r {[y j]}];
s=[s g+norm(b(z,:)-b(j,:))];
end
end
end
p=r;
q=s;
end
[i j]=max(w);
c=f{j};