此模型用到的程序。
M=inf;
a=矩阵（方阵）；
b=a+a';b(find(b==0))=M;
for i=1:53
b(i,i)=0;
end
[D,path] = floyd(b);
fenzu1 = [15 16 17 18 20 21 22 23 24 25 26 27 28 44 46 48 49 50 51]; zuyi = D(fenzu1,fenzu1); fenzu2=[2 5 6 7 8 9 10 11 12 13 14 19 40 41 42 43 45 47 50]; zuer=D(fenzu2,fenzu2); fenzu3=[1 3 4 29 30 31 32 33 34 35 36 37 38 39 50 52 53]
zusan=D(fenzu3,fenzu3);
model: sets: city / 1.. 17/: u; link( city, city): dist, ! 距离矩阵; x; endsets n = @size( city); data: !距离矩阵，它并不需要是对称的; dist =@ole('E:\完全图.xls','dist4');//dist是floyd算法产生的方阵；例如：zuyi ,zuer,zusan; !随机产生，这里可改为你要解决的问题的数据; enddata !目标函数;