Tabu search by matlab to solve TSP

% function [best_sofarbestV_sofar]=ftabu（VtableFeatableNtabu_tenurepenal_iter）
%This is a tabu search function.
%The default objective is minimization
%The move is exchange here
%Another function fcost is needed to calculate each move value
%[val]=fcost（Vtable） 
%Vtable is a table includes all available movement combinations.
%Featable is a cell table having the same dimensions as Vtable. For each cell 
%it consistes of a table which including 0 or 1 to indicate the availability 
%of movements between nodes. For example if node1（ij） and node2（uv） 
%cannot be exchanged then Featable{ij）（uv）=0;
XY=[41 94;37 84;54 67;25 62;
    7 64;2 99;68 58;71 44;54 62;
    83 69;64 60;18 54;22 60;83 46;
    91 38;25 38;24 42;58 69;71 71;
    74 78;87 76;18 40;13 40;82 7;
    62 32;58 35;45 21;41 26;44 35;
    4 50];%node coordinates
% XY=[4272;1747;1325;3517;];

SZ=size（XY1）;
% Vtable=[182324911719201021141582526292827161722233051213461;];
Vtable=[231891171920102114158252629282716172223305121346241;];

crnt_tour=Vtable;
short_path=Vtable;
figure （1）;
plot（[XY（short_path（SZ）1）XY（short_path（1）1）][XY（short_path（SZ）2）XY（short_path（1）2）]）;
plot（XY（SZ1）XY（SZ2）‘bd‘‘MarkerFaceColor‘‘g‘‘MarkerSize‘10）;
hold on;

for i=1:SZ-1
    plot（XY（i1）XY（i2）‘bd‘‘MarkerFaceColor‘‘g‘‘MarkerSize‘10）;
    plot（[XY（crnt_tour（i）1）XY（crnt_tour（i+1）1）][XY（crnt_tour（i）2）XY（crnt_tour（i+1）2）]）;
    
end

hold off;




Featable=cell（size（Vtable））;
% for i=1:size（Vtable1）
%     for j=1:size（Vtable2）
%         for u=1:size（Vtable1）
%             for v=1:size（Vtable2）
%                 if i==u && j==v || j==1 || v==1 || Vtable（ij）==0 || Vtable（uv）==0
%                     Featable{ij}（uv）=0;
%                 else
%                     Featable{ij}（uv）=1;
%                 end
%             end
%         end
%     end
% end

for i=1:size（Vtable2）
    for ii=1:size（Vtable2）
        if i~=1 && i~=length（Vtable） && ii~=1 && ii~=length（Vtable） && i~=ii
            Featable{1i}（1ii）=1;
        else
            Featable{1i}（1ii）=0;
        end
    end
end

%Featable generation
N=1000;
penal_iter=60;
tabu_tenure=50;


track=size（N3）;



%************************Function start here*******************************
fprintf（‘\n\n\n**************ftabu start here*****************‘）;
start_time = cputime;

dim1=size（Vtable1）;
dim2=size（Vtable2）;

tabu_list=cell（dim1dim2）;
costtable=cell（dim1dim2）;
frequency = cell（dim1dim2）;

for i=1:dim1
    for j=1:dim2
        tabu_list{ij}=zeros（dim1dim2）;
        costtable{ij}=zeros（dim1dim2）;
        frequency{ij}=zeros（dim1dim2）;
    end
end

for i=1:dim1
        for j=1:dim2
            for u=1:dim1
                for v=1:dim2
                    if Featable{ij}（uv）==0
                        tabu_list{ij}（uv）=10000;