function [y]=mutation(x,pmutation,lchrom)
% variable 'x' stands for a simple chromosome
y=x;
for m=1:lchrom
    if flip(pmutation)
      y(m)=xor(x(m),1);
%     else
%         y(m)=x(m);
    end
end