for i=1:40
    if rand<0.1 %通過變異概率
        M_pos=round(10*rand);
        if M_pos~=0 %若變異位為0則無意義
           pop(i,M_pos)=1-pop(i,M_pos);
       end
     end
     if rand<0.1 %通過變異概率
        M_pos=round(10*rand);
        if M_pos~=0 %若變異位為0則無意義
           pop(i,10+M_pos)=1-pop(i,10+M_pos);
        end
     end
end
 
clear M_pos;
clear i;
clear k;