關於matlab的問題

這個程式能正常作用.他取點的方式為亂數產生,有沒辦法改成為自己設定呢!?
還有蟻群的矩陣一定要用方陣才可以嗎?可改成非方陣嗎?

%蟻周模型,解決TSP問題
%clc
%初始化
format short;
d=[11,2,5;7,8,7;8,6,3];
n=length(d);  %n 城市數目
m=5;  %m 螞蟻數量
Nmax=100;%最大循環次數
%d(i,j) 城市i,j之間的距離,d is a nn matrix

y=zeros(n,n);%y(i,j)=1/d(i,j) 在TSP問題中,啟發信息
for i=1:n
    for j=1:n
        y(i,j)=1/d(i,j);
    end
end
e=1;%信息啟發因子
f=1;%期望啟發因子
Q=20;%
S=ones(n,n);%(i,j)路段初始化起始信息素
for i=1:n
    for j=1:n
         
           S(i,i)=0;
       
    end
end

S1=zeros(n,n);%(i,j)路段信息素增量

notallowed=ones(m,n);%禁忌表,0表示已經訪問過
a=zeros(m,n);%螞蟻循環一次的路徑
for k=1:m
    a(k,1)=1+round(rand
(n-1));%將螞蟻隨機放到n個城市上
end
for k=1:m % 將初始城市放入禁忌表中
    notallowed(k,a(k,1))=0;
end

for N=1:Nmax  %N 循環次數
    t=2;
    L=zeros(1,m);
    while t<=n %重複直至禁忌表滿為止
        for k=1:m
            %計算螞蟻k轉移的概率
            i=a(k,t-1);
            p=zeros(1,n);%p(j)螞蟻k選擇路徑i,j的概率
            for j=1:n
                if notallowed(k,j)~=0
                   u=(S(i,j)^e)(y(i,j)^f);
                   v=0;
                   for w=1:n
                       v=v+(S(i,w)^e)
(y(i,w)^f)*notallowed(k,w);
                   end
                   if v~=0
                      p(j)=u/v;
                   end
                end
            end
            [pk,j]=max§;
            notallowed(k,j)=0;
            L(k)=L(k)+d(i,j);
            a(k,t)=j;
        end
        t=t+1;
    end
    for k=1:m
    L(k)=L(k)+d(a(k,n),a(k,1));
    end
    %一次循環結束,回到起始位置
  
    %更新
    for k=1:m
        for i=1:n-1
           
        end
       
    end
    for i=1:n
        for j=1:n
            if d(i,j)~=inf
               for k=1:m
                 
               end
            end
        end
    end
    for i=1:n
        for j=1:n
            if d(i,j)~=inf
            S(i,j)=(1-rand)*S(i,j)+S1(i,j);
            end
        end
    end
    for k=1:m %將禁忌表中除起始城市,全都置為未訪問
        for t=1:n
            if t~=a(k,1)
               notallowed(k,t)=1;
            end
        end
    end
    S1=zeros(n,n);%(i,j)路段信息素增量清零
    
   
end %循環最大次數結束
[result,k]=min(L)
a(k,:slight_smile: