Dist = zeros(1,n);
Distx = zeros(n,numel(x)); % <-- ADDED, allocate for saved x vectors
for k=1:n
xp = randperm(numel(x), 2);
x(xp) = x(fliplr(xp)) % vector x permuted in two positions
s = sub2ind(size(d),x(1:end-1),x(2:end ));
Dist(k) = sum(d(s));
Distx(k,:) = x; % <-- ADDED, save the x vector for this iteration
Distance = Dist(1,k) % travaled distance by elements of the vector x
end
[Lowest_Distance,ix] = min(Dist); % <-- MODIFIED, get the index of the min
Lowest_Distancex = Distx(ix,:); % <-- ADDED, pick off x vector corresponding to the min