You are an expert in writing code in many different languages. Your goal is to perform code completion for the following code, keeping in mind the rest of the code and the file meta data.
    Metadata details:
    filename: antTSP.m
    path of file: ./repos/opt_matlab/第14章/例14.1
    
    the code of the file until where you have to start completion: function [Shortest_Route,Shortest_Length]=antTSP(city,iter_max,m,Alpha,Beta,Rho,Q)
%city：n个城市的坐标，n×2的矩阵
% iter_max ：最大迭代次数
% m： 蚂蚁个数
%Alpha ：表征信息素重要程度的参数
%Beta ：表征启发式因子重要程度的参数
% Rho ：信息素蒸发系数

    
function [Shortest_Route,Shortest_Length]=antTSP(city,iter_max,m,Alpha,Beta,Rho,Q)
%city：n个城市的坐标，n×2的矩阵
% iter_max ：最大迭代次数
% m： 蚂蚁个数
%Alpha ：表征信息素重要程度的参数
%Beta ：表征启发式因子重要程度的参数
% Rho ：信息素蒸发系数
% Q ：信息素增加强度系数
% R_best 各代最佳路线
[n,D]=city2d(city);
Eta=1./D;          %Eta为能见度因数，这里设为距离的倒数
Tau=ones(n,n);     %Tau为信息素矩阵
Tabu=zeros(m,n);   %存储并记录路径的生成
nC=1;             %计数器
R_best=zeros(iter_max,n);
L_best=inf.*ones(iter_max,1);%%各代最佳路线及长度
while nC<=iter_max    
  Tabu(:,1)=ceil(rand(m,1).*n);
   %m只蚂蚁按概率函数选择下一座城市，完成各自的周游
  for j=2:n
    for i=1:m
        visited=Tabu(i,1:(j-1));  %已访问的城市
        J=zeros(1,(n-j+1)); %待访问的城市
        P=J; Jc=1; %待访问城市的选择概率分布
        for k=1:n
           if length(find(visited==k))==0
               J(Jc)=k;Jc=Jc+1;
           end