w=(b(:,1)+1)*inf;
w=(b(:,1)+1)/0;
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 all(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};