suckysuck2

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 o
  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)
  d = [d 0];
  w = z{i};
  for j=2:nnz(w)
    d(i) = d(i) + s(w(j-1),w(j));
  end
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