function [stpop,D1,N1]=FitnessEvaluate(i_pop,i_n, i_fitInd, o_share1,o_share2);
s=size(i_pop);
tempop=i_pop;
%i_n=s(1,2)-2;
for i=1:s(1,1)
    D(i,1)=0;
    N(i,1)=0;
    for j=1:s(1,1)
        if (i_pop(j,i_n+1)<=i_pop(i,i_n+1)) && (i_pop(j,i_n+2)<=i_pop(i,i_n+2)) && ...
                ((i_pop(j,i_n+1)+i_pop(j,i_n+2))<(i_pop(i,i_n+1)+i_pop(i,i_n+2)))
            D(i,1)=D(i,1)+1;
        end
%         fdij=sqrt((i_pop(j,i_n+1)-i_pop(i,i_n+1))^2+(i_pop(j,i_n+2)-i_pop(i,i_n+2))^2);
%         if fdij<=o_share1
        f1dij=abs(i_pop(j,i_n+1)-i_pop(i,i_n+1));
        f2dij=abs(i_pop(j,i_n+2)-i_pop(i,i_n+2));
%         if f1dij<=o_share1
        if (f1dij<=o_share1) && (f2dij<=o_share2)
            N(i,1)=N(i,1)+1;
        end
    end
    tempop(i,i_n+3)=(1+D(i,1))*N(i,1);
end
%F=tempop(:,i_n+3);
%f=sort(F);
%stpop=sortD(tempop,F,f); 
stpop=tempop;
D1=D;
N1=N;
%D1=sortD(D,F,f);
%N1=sortD(N,F,f);