ada_sum=0;
for i=1:40
    ada_sum=ada_sum+adapt(i);
    e_adapt(i)=0;
end
for i=1:40
    e_adapt(i)=40*adapt(i)/ada_sum;
end
r=1;     i=1;
while (r<=40)&(i<40)
   if e_adapt(r)>=1.0
      while(e_adapt(r)>=1.0)
        
           for k=1:20       
               new_pop(i,k)=pop(r,k);
           end
     
        e_adapt(r)=e_adapt(r)-1;
        i=i+1;
      end 
   end
     r=r+1;
 end
while(i<=40)
     q=1;
     for j=2:40
        if e_adapt(q)<e_adapt(j)
              q=j;
        end
     end   
     for k=1:20       
        new_pop(i,k)=pop(q,k);
      end
      e_adapt(q)=0;
      i=i+1;
  end    
%將選擇產(chǎn)生的新群體復(fù)制給pop種群
for i=1:40
     for j=1:20
        pop(i,j)=new_pop(i,j);
     end
end
 
clear i;
clear j;
clear k;
clear r;
clear q;