Динамические матрицы в процедуре, Есть проблема с динамическими данными... Опции

[P-Tigr]

Столкнулся с одной проблемой, 2 дня уже бьюсь... что делать - не знаю...
Если коротко, то дело вот в чем.
Я написал процедуру для подсчета определителя матрицы любого порядка (методом Гаусса). И там над матрицей выполняются различные преобразования (складывания строк и т.д.).
Причем - внимание - по условию задания матрица должна быть обязательно динамической!
Процедура работает нормально, определитель находит, НО в то же время изменяет исходную матрицу, а этого нельзя допустить!! Происходит это, как мне кажется, из-за того, что такие матрицы - ссылки на память, поэтому передать ее в процедуру строго по значению нельзя...

Замучился я с этими динамическими структурами...
Please, умные люди, help me!

[P-Tigr]

Извините пожалуйста, еще 1 вопрос, помогите разобраться!
Почему при вводе некоторых данных, например
1 1 1
1 1 1
1 1 1
в блоке функции нахождения определителя:

for j:=0 to high(AC)-1 do
    begin
      for i:=j+1 to high(AC) do
        begin
          mn1:=AC[i,j]; mn2:=AC[j,j];
// ОШИБКА - в следующей строке   генерируется       
          mnoj:=mn1/mn2;
          for jk:=0 to high(AC) do
            begin
              AC[i,jk]:=AC[i,jk]-AC[j,jk]*mnoj;
            end;
        end;
    end;

в строке
mnoj:=mn1/mn2;

появляется ошибка, связанная с плавающей запятой:
"Invalid floating point operation" или "Floating divizion by zero".
ПОЧЕМУ?

P.S.1. Ошибка появляется и при вводе разных чисел, но правда не всегда - очень странно...
P.S.2. Я присоединил к письму exe-файл в rar-архиве - чтоб удобней было тестировать.

[volvo]

появляется ошибка, связанная с плавающей запятой:
"Invalid floating point operation" или "Floating divizion by zero".
ПОЧЕМУ?

Это не связано с вводимыми данными... Только с тем, что у тебя в отмеченных строках просто напросто происходит выход за границы массива (не забывай, что динамические массивы индексируются с 0 до N-1), и кто знает, какие значения там хранятся...

for j:=0 to high(AC)-1 do
   begin
     for i:=j+1 to high(AC) do  // **** Здесь - i принимает значение на 1 больше чем нужно
       begin
         mn1:=AC[i,j]; mn2:=AC[j,j];

         mnoj:=mn1/mn2; // ОШИБКА

         for jk:=0 to high(AC) do
         // **** Здесь - jk тоже принимает значение на 1 больше, чем допустимо,
         // это тоже может привести к проблемам
           begin
             AC[i,jk]:=AC[i,jk]-AC[j,jk]*mnoj;
           end;
       end;
   end;

Я совсем забыл написать, что я менял границы индексов когда проверял... Привычка...

P.S. Смотри и выше в тексте функции. Там тоже есть выход за пределы массива...