S=[-500*ones(1,10);500*ones(1,10)];
pop=genrpop(30,S);

for i= 1:1000
    Naj = testfn3(pop);
    Newpop = selbest(pop,Naj,5);
    Popina = selsus(pop,Naj,25);
    krizenie = crossov(Popina,1,0);
    nove = mutx(krizenie,0.1,S);
    nove1 = muta(nove,0.1,ones(1,10)*5,S);
    pop= [Newpop;nove1];
    Min(i)=min(Naj);
end
pop(1,:)
plot(Min);