function [sol] = AcatspEval(sol,distMatrix,numvars)
% 功能－－－計算適應(yīng)度
% 輸入：sol一個個體
% 輸出：該個體及其適應(yīng)度val

val = sum(diag(distMatrix(sol(1:numvars),[sol(2:numvars) sol(1)])));
sol(numvars+1)=val;