clc
clear all
nind=80;
fielddr=[0 -20 -5;0 27.001 17.00];
chrom=crtrp(nind,fielddr)
y=objectval(chrom);
maxval(1)=max(y);
avgval(1)=sum(y)/80;
minval(1)=min(y);
maxgen=100;
tic;
for i=2:maxgen
    fit=ranking(y')
    for j=1:80
           m=max(fit);
           if fit(j)>=m
               topval=chrom(j,:);
           end
    end    
    a=select('sus',chrom,fit);
    newch=xovsp(a,0.6);
    %newch1=mutate('mut',newch,2*ones(1,20),0.0001);
    newch1=mutbga(newch,fielddr,[0.0001,1]);
    y=objectval(newch1);
   fit=ranking(y');
   for j=1:80
       m=min(fit);
       if fit(j)<=m
          newch1(j,:)=topval;
       end
    end 
    chrom=newch1;
    y=objectval(newch1);
    maxval(i)=max(y);
    avgval(i)=sum(y)/80;
    minval(i)=min(y);
end
toc;
t=toc;
max(y);
x=1:maxgen;
plot(x,-maxval,'b-',x,-avgval,'r-',x,-minval,'k-')
hold on;
%plot(x,-100*ones(1,maxgen),'w');
hold off;