function popparm1=density(parpopsize,conpoppar1,n1)
%
popparm1=conpoppar1;
totalfitness=sum(conpoppar1(:,n1+1));
tataldensity=sum(1./conpoppar1(:,n1+n2));
b=0.75; %
if totalfitness==0
    popparm1(:,n1+3)=b/parpopsize+(1-b)*conpoppar1(:,n1+2)/tataldensity;
else
    popparm1(:,n1+3)=b*conpoppar1(:,n1+1)/totalfitness+((1-b)