Паскаль вычисление определителя

[Юля23]

Проблема с вычислением определителя матрицы любой размерности...В чем ошибка? как правильно сделать?

uses
  crt;
const
  l=10;
type
  matr=array[0..l,0..l] of integer;
var
 M, A:matr;
  i,j,p,x,k,n,t,r, s: integer;
  D:smallint;

Procedure inputMatrix(var x:matr);
begin
  for i:=0 to n-1 do
    begin
      for j:=0 to n-1 do
        begin
          x[i , j]:=random(10);
          Write(x[i,j]:4);
        end;
      writeln;
    end;
end;

begin
  randomize;
  clrscr;
  writeln('vvedite razmer matrici');
  readln(n);
  writeln;
  writeLn('ishodnaya matrica:');
  inputmatrix(a);

  for i:=0 to n do   {процесс транспонирования}
   for j:=0 to i-1 do
    begin
     k:=A[i,j];
     A[i,j]:=A[j,i];
     A[j,i]:=k;
    end;

  writeln('результат транспонирования: ');
  for i:=0 to n-1 do
   begin
    for j:=0 to n-1 do
     write(A[i,j]:3);
    writeln;
   end;
   WriteLn('ishodnaya matrica');
   inputmatrix(a);
p := n-1 div 2;{поворот на 90}
for i:=0 to p do
for j:=i to n-i do
  begin
    x := A[i,j];
    A[i,j] := A[n-j,i];
    A[n-j,i] := A[n-i,n-j];
    A[n-i,n-j] := A[j,n-i];
    A[j,n-i] := x;
  end;

writeln('Povorot na 90 gradusov po chasovoi strelke:');
for i:=0 to n-1 do
  begin
    for j:=0 to n-1 do
    write(a[i,j]:3);
    writeln;
  end;
  WriteLn('ishodnaya matrica:');{вычисление определителя}
  inputmatrix(a);
  M := a;
for s:=i to n-1 do
 for t:=1 to n do
  M[s,t] := M[s+1,t];
for t:=1 to n-1 do
 for s:=1 to n-1 do
  M[s,t] := M[s,t+1];
 WriteLn('vichislenie opredelitelya');
   D:=Exp((j+1)*Ln(-1))*r[i,j]*M[i,j];
  writeln;
end.

[newerow1989]

Я бы рекомендовал скопировать матрицу в отдельную переменную, а потом преобразовывать

var X,A:matr;
***
begin
   For i:=0 to n-1 do
      For j:=0 to n-1 do
         X[i,j]:=A[i,j];
***
   for i:=0 to n-1 do   {процесс транспонирования}
      for j:=0 to n-1 do
         A[j,i]:=X[i,j];
***
end;

Что за запись:

p := n-1 div 2;{поворот на 90}

Может

p := (n-1) div 2;{в скобках}

О, а это что означает?

WriteLn('ishodnaya matrica:');{вычисление определителя}
inputmatrix(a);//ввод новой матрицы, что ли? (там же стоит x[i,j]:=random(10);)
M := a;//по-моему не получится

Ln(-1) - не существует! И для чего она нужна?

D:smallint;//целые числа
D:=Exp((j+1)*Ln(-1))*r[i,j]*M[i,j];

[Юля23]

О, а это что означает?

WriteLn('ishodnaya matrica:');{вычисление определителя}
inputmatrix(a);//ввод новой матрицы, что ли? (там же стоит x[i,j]:=random(10);)

Это пока незаконченная программа, там потом просто нужно будет меню ввести, которое будет обеспечивать выбор конкретного преобразования, поэтому каждый раз новая матрица.

Ln(-1) - не существует! И для чего она нужна?

D:smallint;//целые числа
D:=Exp((j+1)*Ln(-1))*r[i,j]*M[i,j];

просто есть такая формула вычисления определителя,
D=(-1)^(j+1)*a[i,j]*M[i,j]

[type_Oleg]

Сообщение от Юля23

... просто есть такая формула вычисления определителя,
D=(-1)^(j+1)*a[i,j]*M[i,j]

Просто логарифм от отрицательного числа не берется. По крайней мере в мире действительных чисел.
Да и вообще, всякие экспоненты, логарифмы с индексами - это моветон.

Там же эта степень - просто чтобы знак определить. При четных - минус, нечетных - плюс.
Используйте или if или , если последовательно вычислять последующий из предыдущего - умножайте на (-1), будет +-+-+- ...

[FPaul]

Юля23

Если почитать Wikipedia, то там говорится о способе расчёта приведением матрицы к треугольному виду.

На страничке volvo приведена рекурсивная реализация разложения Лапласа - это то, что пытаешься реализовать ты.