MATLAB function: reverse.m
function reverse(ncity,n)
% reverse(ncity,n)
global iorder;
nn=(1+mod((n(2)-n(1)+ncity),ncity))/2;
for j=1:nn
k=1 + mod((n(1)+j-2),ncity);
l=1 + mod((n(2)-j+ncity),ncity);
itmp=iorder(k);
iorder(k)=iorder(l);
iorder(l)=itmp;
end
-----------------------------------...