function [bestpath,minpath]=IATSP(TSPMatrix,CM,Popsize,TC,Pm,Pc)
path=Initpath(TSPMatrix,Popsize);
MemoryLab=InitLab(TSPMatrix,CM);
[ktjqhl,ktyqhl]=qinheliCalculate(path,TSPMatrix);
MemoryLab=GenerateM(path,ktyqhl,CM,TSPMatrix);
for i=1:TC
    path=UGenerateNewpath(path,MemoryLab);
    SelectP=CalculateSelectP(ktjqhl,ktyqhl);
    path=Select(path,SelectP);
    path=GGenerateNewpath(path);
    [ktjqhl,ktyqhl]=qinheliCalculate(path,TSPMatrix);
    MemoryLab=GenerateM(path,ktyqhl,CM);
end
[bestpath,minpath]=Output(path,TSPMatrix);