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

[volvo]

 const
  NMAX = 20;
  MMAX = 20;

увеличивать не пробовал? У тебя же опять выход за пределы матрицы должен произойти (файл большой, а матрица - маленькая...)

[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;
{ Далее - по тексту... }

[ProtasSoft]

Можно это так обойти

function Gamma(x: extended): extended;
const
pi	= 3.1415926;

var
i,j	: integer;
y,gam	: extended;

begin		{ gamma function }
If x=0 then gamma:=1;
  if x>0.0 then
    begin
      y:=x+2.0;
      gam:=sqrt(2*pi/y)*exp(y*ln(y)+(1-1/(30*y*y))/(12*y)-y);
      Gamma:=gam/(x*(x+1))
    end;
  if  x<0 then
    begin {x<0}
      j:=0;
      y:=x;
      repeat
	j:=j+1;
	y:=y+1.0
      until y>0.0;
      gam:=gamma(y);		{ recursive call }
      for i:=0 to j-1 do
	gam:=gam/(x+i);
      Gamma:=gam
    end; { x<0 }
end;		{ gamma function }

[klem4]

Так к слову : в паскале есть встроенная константа Pi

[volvo]

Так к слову : в паскале есть встроенная константа Pi

:yes: ProtasSoft, ты будешь смеяться, но и функция Power тоже есть (в модуле Math)