clear all;
citynum=[0.4 0.4439;0.2439 0.1463;0.1707 0.2293;0.2293 0.761;0.5171 0.9414;0.8732 0.6536;0.6878 0.5219;0.8488 0.3609;0.6683 0.2536;0.6195 0.2634];
Vx=[0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000
    0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000
    0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000
    0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000
    0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000
    1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
    0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000
    0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000
    0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
    0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000];
for j=1:10
    for i=1:10
        if Vx(i,j)==1
           G(j,1:2)=citynum(i,:);
           G(j,3)=i;
        end    
    end    
end 
dis2=0;
dis1=0;
for i=1:10
    if i==10
       dis1=dis1+distance(G(i,1),G(1,1),G(i,2),G(1,2));
    else   
    dis2=dis2+distance(G(i,1),G(i+1,1),G(i,2),G(i+1,2));
    end       
end  
dis=dis1+dis2
G=[G;G(1,:)];
plot(G(:,1),G(:,2),'bo-')