suckysuck

    
N = size(b,1);
s = [];
for i=1:N
  for j=1:N
    s(i,j) = sqrt(sum((b(i,:)-b(j,:)).^2));
  end
end

o = 1;
z = {};

while(size(o,1))
  p = [];
  for n=1:size(o,1)
    w = o(n,:);
    for m=2:N
      if ~ismember(m, w) & a(m, w(end))
        p = [p; w m];
      end
    end
    z = {z{:} w};
  end
  o = p;
end

d = [];
for i=1:numel(z)
  l = 0;
  w = z{i};
  for j=2:nnz(w)
    l = l + s(w(j-1),w(j));
  end
  d = [d l];
end

L = 0;
c = 1;
for n=2:N
  l = -1;
  v = [];
  for m=1:numel(z)
    w = z{m};
    if w(end) == n
      if l<0 | d(m) < l
        v = w;
        l = d(m);
      end
    end
  end
  if l > L
    c = v;
    L = l;
  end
end