%function popnew=mutation(pop,pmut,n,length)
function popnew=mutation(pop,gen,N,n,length1,length2)
%變異函數
%pmut為變異概率
k=1;
while k<=n
     rk=rand();
     %if rk<pmut
     if rk<0.001+(0.1*gen)./N;    
        pos=fix(rand()*(length1+length2))+1; %隨機產生變異點
        pop(k,pos)=~pop(k,pos); %對變異點進行變異
     end
     k=k+1;
end
popnew=pop;