function pop=mutation2(pmutation,crosspop)
[mpop,npop]=size(crosspop);

   for i=1:mpop
       temp=rand;
       if pmutation>temp
           j1 = rand_function(1,npop)
           j2 = rand_function(1,npop) 
           while j1==j2
                j1 = rand_function(1,npop); 
                j2 = rand_function(1,npop); 
           end
           if j1>j2
               j_int=j2;
               j2=j1;
               j1=j_int;
           end
           p=crosspop(i,:);
           p_cho=p(j1:j2);
           for j=j1:j2
           crosspop(i,j)=



for i=1:mpop
       for j=1:npop
           if rand<pmutation
                  if crosspop(i,j)==1
                     crosspop(i,j)=0;
                  else
                     crosspop(i,j)=1;
                  end
           end
      end
end
pop=crosspop;