Ошибка “Invalid floating point operation”

  procedure StsMasID(masID: array of real;szmID:integer;var mo,disp: real);
   var
   i: integer;
   sum: real ;
   begin
     sum:= 0;
     for i := 0 to szmID-1 do
        sum:=sum + masID[i];
      mo := sum/szmID;         // на это ругается
      sum:= 0;   
     for i := 0 to szmID-1 do
        sum:=sum + (masID[i]-mo)*(masID[i]-mo);
      disp:=sum/(szmID-1);
  end;

ошибка в этой строчке  mo := sum/szmID;

function Erlang(L:real; k: integer):real;
  var
   j: integer;
   x,s: real;
 begin
     randomize;
     s:= 0;
  for j:= 1 to k do
      s:=s+ln(random(k));   //Если s = 0 и random(k) = 0, то получите деление на 0
    x:= - 1/(L)*s;
Erlang:=x;
end;

      s:=s+ln(random(k));   //Если s = 0 и random(k) = 0, то получите деление на 0

if random(k) > 0 then begin
s:=s+ln(random(k));
end;

var
j, N : Integer;
...
  for j:= 1 to k do begin
    N:=random(k);
    if N <> 0 then s:=s+ln(N);
  end;

function Erlang(L:real; k: integer):real;
  var
   j, N: integer;
   x,s: real;
 begin
     randomize;
     s:= 0;
  for j:= 1 to k do
  N:=random(k);
    if N <> 0 then s:=s+ln(N);
  end;

mo := sum/szmID;

FrmProizvIS.EditLO.Text:=FloatToStr(Lo);
       FrmProizvIS.EditKO.Text:=FloatToStr(Ko);

function Erlang(L:real; k: integer):real;
  var
   j, N: integer;
   x,s: real;
 begin
     randomize;
     s:= 0;
  for j:= 1 to k do
  N:=random(k);
    if N <> 0 then s:=s+ln(N);
  end;

procedure funcRaspr();   // формирование массива  функции распределения
     var
          i,j,k,n:integer;
        tekZn,intN,intK: real;
         KolInt: integer;              
        mo,disp,maxID,minID,LO,KO: real;
 begin
   if RedTxt then     // если masID редактировался
    begin
     for I := 1 to szmID-1   do
            begin
       masID[i]:= StrToFloat(FrmProizvIS.StrGridID.Cells[i,1]);  //восстановление отредактированного masID
            end;
        RedTxt:=false;
    end;
     KolInt:=StrToInt(FrmProizvIS.EditKolInt.Text);
        ProcUnit.minT(masID,szmID,minID);
        ProcUnit.maxT(masID,szmID,maxID);
         TintFr:= round((maxID-minID)/KolInt)+1;
         FrmProizvIS.EditTint.Text:=IntToStr(TintFr);
      szmFR:=KolInt;
         SetLength(masFR,szmFR);
      for i := 0 to szmFR-1  do
        begin
             intN:=(i)*TintFr;
            intK:=(i+1)*TintFr;
           k:=0;
         for j := 0 to szmID-1  do
            begin
            tekZn:=  masID[j];
          if (intN < tekZn) And (tekZn <= intK) then
               k:=k+1;
            end;
         masFR[i]:= k;
        end;
       n:= FrmProizvIS.StrGridFR.ColCount;
FrmProizvIS.LabelRI.Caption:= 'Распределение по интервалам Tint= '+IntToStr(TintFr)+' сек';
         for I := 0 to n-1  do
           begin
             FrmProizvIS.StrGridFR.Cells[i+1,0]:=IntToStr(i);
             FrmProizvIS.StrGridFR.Cells[i+1,1]:= FloatToStr(masFR[i]);
           end;
       ProcUnit.StsMasID(masID,szmID,mo,disp);
       FrmProizvIS.EditMo.Text:=FloatToStr(mo);
       FrmProizvIS.EditDisp.Text:=FloatToStr(disp);
     Lo:= mo/disp; // вот тут ошибка при нажатии на F8
     Ko:= Lo*mo;
       FrmProizvIS.EditLO.Text:=FloatToStr(Lo);
       FrmProizvIS.EditKO.Text:=FloatToStr(Ko); // ссылается сюда, при наведении каретки на Ko пишет Ko=0
 end;