function pop2=mutation(pop,pm,popsize,n1)
%
%
%
%
%

pop2=pop;
for i=1:popsize
    if(rand<pm)
        mpoint=round(rand*(n1-1))+1;
        pop2(i,mpoint)=abs(pop(i,mpoint)-1);
    end
end