load citys_data.mat n = size(citys,1); D = zeros(n,n); for i = 1:n for j = 1:n if i ~= j D(i,j) = sqrt(sum((citys(i,:) - citys(j,:)).^2)); else D(i,j) = 1e-4; end end end m = 50; alpha = 1; beta = 5; rho = 0.1; Q = 1; Eta = 1./D; Tau = ones(n,n); Tab