procedure TEvPopulus.Digenesize;
var
  i, j, n, nPopu, nBorn, nDead: Integer;
  fFav, fAch, fRep, fSur, fElm, fAge, fFert, fR: Double;
  P, F: TEvIndividual;
  lBabies: TEvIndividuals;
  L: TStringList;
  FF: TIniFile;
begin
  FF := TIniFile.Create(FileName);

  fFav := Individuals[0].Favorable;
  Individuals.SortByFavorable;
  fAch := Individuals[0].Favorable;
  if fAch > fFav * 1.01 then
    FF.WriteString(‘Specimens’, IntToStrN(LineageNo, 6),
        Individuals[0].Genes[0].AsHex + ‘,’
      + FormatFloat(‘0.000000′, fAch));

  nPopu := Individuals.Count;
  fFav := 0;
  i := 0;
  while (i < nPopu) and (Individuals[i].Experience > 0) do
  begin
    fFav := fFav + Individuals[i].Favorable;
    Inc(i);
  end;
  if i > 0 then fFav := fFav / i;

  fElm := Max(0.0, (0.00001 + nPopu – Population) / (nPopu * 1.0));

  fFert := Fertility * Min(3, Population / nPopu);

  lBabies := TEvIndividuals.Create(False);
  L := TStringList.Create;
  L.Sorted := True;
  L.Duplicates := dupIgnore;
  fAch := 0; fAge := 0;
  nDead := 0;
  i := 0;
  while i < Individuals.Count do
  begin
    P := Individuals[i];
    L.Add(P.Genes[0].AsHex);
    fAch := fAch + P.Favorable;
    fAge := fAge + P.Experience;
    if fFav > 0 then
    begin
      fRep := Power(P.Favorable / fFav, EOR);
      fSur := Power(P.Favorable / fFav, EOS);
    end else
    begin
      fRep := 1;
      fSur := 1;
    end;
    P.Reproductivity := fRep;
    P.Survivability  := fSur;

    fRep := fFert * fRep;
    n := Trunc(fRep);
    fR := Random;
    if fR <= fRep – n then
      Inc(n);
    for j := 0 to n – 1 do
    begin
      fR := Random;
      F := IndividualClass.Create;
      F.Inherit(P, fR <= MutateRate);
      F.Develop;
      lBabies.Add(F);
    end;

    fR := Random;
    if fR * fSur <= fElm then
    begin
      Individuals.Delete(i);
      Inc(nDead);
    end else
      Inc(i);
  end;

  nBorn := lBabies.Count;
  for i := 0 to nBorn – 1 do
    Individuals.Add(lBabies[i]);

  lBabies.Free;

  fAge := fAge / nPopu;
  Achievement := fAch / nPopu;
  Variety := L.Count / nPopu;

  FF.WriteString(‘Progress’, IntToStrN(LineageNo, 6),
      FormatFloat(‘0.000000′, Achievement) + ‘,’
    + FormatFloat(‘0.000000′, Individuals[0].Favorable) + ‘,’
    + FormatFloat(‘0.0000′, Variety) + ‘,’
    + IntToStrN(nBorn, 6) + ‘,’
    + IntToStrN(nDead, 6) + ‘,’
    + FormatFloat(‘#0.0′, fAge));
  FF.Free;

  LineageNo := LineageNo + 1;

  L.Free;
end;