Разложение Холецкого, Делфи

[DeeKay]

Собственно, тема в заголовке.

Формула выглядит как:


? I>J

    begin //Собственно формула
      IsLower:=false;
      for i:=1 to InputMatrix.RowCount do
        for j:=1 to InputMatrix.RowCount do
          if i>j then
            for k:=1 to j-1 do
              try
                LMat[i,j]:=(InputMat[i,j]-(LMat[i,k]*LMat[j,k]))/LMat[j,j];
              except
                LMat[i,j]:=(InputMat[i,j]/LMat[i,i]);
              end;
          if i=j then
            begin
              for k:=1 to i-1 do
                LMat[i,i]:=sqrt(InputMat[i,i]-LMat[i,k]*LMat[i,k]);
            end;

          if i<j then
            LMat[i,j]:=0;
    end

InputMat - массив в который записаны значения из стринггрида на форме.

на строчке LMat[i,j]:=(InputMat[i,j]-(LMat[i,k]*LMat[j,k]))/LMat[j,j]; вылазит ЕOverflow, Floating point overflow. нужна помощь

Может заменить сумматор на отдельную переменную?

[Stilet]

А тип элемента матрицы какой?