function [x1,changex,l1]=drop(x0,l0)

global customer;
global depot;
global demand;    

NumberRoute=length(l0);
  a=ceil(NumberRoute*rand);
  %先計算DD數(shù)組
  %nowdepot為選擇的配送中心
  nowdepot=x0(a,1);
  %計算a路線上的重心
  sum1=0;sum2=0;sum3=0;
    for j=2:l0(a)+1
        sum1=demand(x0(a,j),1)*demand(x0(a,j),3)+sum1;
        sum2=demand(x0(a,j),3)+sum2;
         sum3=demand(x0(a,j),2)*demand(x0(a,j),3)+sum3;
    end
      xcenter=sum1/sum2;
      ycenter=sum3/sum2;
  for j=1:depot
            alldis(j)=distan(demand,xcenter,ycenter,j+customer);
  end
   [dis,n]=min(alldis);
%if  find(x0(:,1)==n)
 % frag=1;
  %else frag=0;
    %end
 %while  frag==0 %看這個最小距離對應(yīng)的配送中心是否開放，如不開放，則將該值刪，從新開始找下一個最小值
    % alldis(n)=100000;
    % [dis,n]=min(alldis);
   %if  find(x0(:,1)==n)
     %  frag=1;
       %else frag=0;
       % end
   %end
while n==x0(a,1)
    alldis(n)=100000;
     [dis,n]=min(alldis);
 end
    
    
x0(a,1)=n;%用選的
l1=l0;
x1=x0;
changex=[1  a n];