Определитель и транспонирование матрицы в Delphi

[Ivan3000]

Привет!Помогите написать программу в delphi,которая будет считать определитель матрицы(nxm),и выполнять транспонирование матрицы(без использования компонента Matrix в Delphi)!

[__STDC__]

Помочь или написать? Транспонирование матрицы вообще три строки кода... определитель лениво писать.. наработки покажите..

[Ivan3000]

транспонирование я сделал,а вот с определителем вообще ничего в голову не идет!Напишите,пожалуйста,код программы для вычисления определителя!
Заранее благодарю!

[__STDC__]

Обратите внимание на LU разложение.. что еще тут можно сказать.. писать я за вас не буду

[Ivan3000]

С определителем ничего не получается!ПОМОГИТЕ,ПОЖАЛУЙСТА!

[__STDC__]

Вот пример нахождения матриц U и L... (Прямой ход метода Гаусса это...)

procedure direct(A:matrix; var U,L:matrix);
var
  i,j,m,k:integer;
  w:extended;
begin
  i := 1;
  j := 1;
  repeat
    for m:=j to n do begin
      w := A[m,j];
      if not(j < 2) then
	for k:=1 to j-1 do
	  w := w-L[m,k]*U[k,j];
      L[m,j] := w;
    end;
    for m:=i+1 to n+1 do begin
      w := A[i,m];
      if not(i < 2) then
	for k:=1 to i-1 do
	  w := w-L[i,k]*U[k,m];
      U[i,m] := w/L[i,i];
    end;
  inc(i);
  inc(j);
  until (i >= m);
end;

дальше понятно что делать?
Единственное что, в таком случае главная диагональ матрицы U находится на главной диагонали матрицы L... соответственно детерминант надо считать из матрицы L

[Ivan3000]

Спасибо большое STDC!!!!

[__STDC__]

Сообщение от Ivan3000

Спасибо большое STDC!!!!

Получилось хоть? Я проверил у себя, на матрице 3*3, у меня сошлись определители =)

[apel]

дайте, пожайлуста, сам проект. чтобы можно было скомпилировать

[BDA]

Сообщение от apel

дайте, пожайлуста, сам проект. чтобы можно было скомпилировать

Вас не смущает, что теме 4 года?