Обработка прямоугольных матриц Опции

[ProtasSoft]

Задана большая матрица вещественных чисел размером 1000 на 2000 в текстовом файле.

Необходимо произвести попарную обработку обработку ее строк, т.е.
для каждой пары СТРОКАi и СТРОКАj строк необходимо вычислить функцию от элементов этих строк, например, Ain - элементы СТРОКИi и Bjn - элементы СТРОКИj.
Из результатов вычисления функции F(Ai1,...,Ai1000;Bj1,...,Bj2000) сформировать новую матрицу и записать в новый файл.

Подскажите как мне действовать.
Насколько я понимаю с большими массивами Borland Pascal не работает.

[volvo]

А вот и причина возникающей ошибки (RunTime Error 200 - "Деление на 0"):

function Gamma(x: extended): extended;
{ ... }
begin { gamma function }
  If x=0 then gamma := 1;
  If x>0.0 then begin { ... } end
  else { x<0 } begin
    { ... }
    gam:=gamma(y);{ recursive call }
    for i:=0 to j-1 do
      gam:=gam/(x+i); { <--- Здесь !!! }
    Gamma:=gam
  end { x<0 }
end; { gamma function }

В указанном мной месте при вычислении Funk[10, 10] у тебя одновременно i = 0 и x = 0 (да, да... При х = 0 ты только присваиваешь функции значение 1, но не выходишь из нее !!!)

P.S. Возможный путь решения проблемы:

If x = 0 Then Begin Gamma := 1; Exit End;
{ Далее - по тексту... }