function L=TourLength(tour,D)

    n=numel(tour);
    
    tour=[tour tour(1)];
    
    L=0;
    for k=1:n
        i=tour(k);
        j=tour(k+1);
        L=L+D(i,j);
    end

end