Предыдущая страница реферата | 16  17  18  19  20  21  22  23  24  25  26 | Следующая страница реферата

function TMainForm.VS_Det(InMass : TVS_MassData): Real;

var

Temp, A: TVS_MassData;

Cols, Rows, Count: Word;

i, j, k: Integer;

begin

Result := 1; //Результат функции по умолчанию

If InMass.N <> InMass.M Then Exit; //если матрица не квадратная - уходим, так как решение методом диагонали

Count := InMass.M ; // Получили размерность исходного массива

SetLength(A.Mass, Count, Count);//Установили размер матрицы

SetLength(Temp.Mass, 1, Count); //Установили размер мартицы

AssignMass(InMass, A); //Во временный массив заносим данные из исходного, чтобы не портить исходный массив

//Поиск и решение

for i := 0 to Count - 2 do {Начало преобразования к верхнему треугольному виду}

begin

for j := i to Count - 1 do {* Поиск }

begin {* нулевых }

Rows := 0; {* строк }

Cols := 0; {* и }

for k := i to Count - 1 do {* столбцов }

begin {* в }

Rows := Rows + Ord(A.Mass[j, k] = 0); {* матрице }

Cols := Cols + Ord(A.Mass[k, j] = 0); {* }

end;{for k := i to Count - 1 do} {* }

if Rows + Cols = 0 then {* }

Break; {* }

if (Cols = Count - i) or (Rows = Count - i) then {* }

begin {* }

Result := 0; {* }

Exit {* }

end {if (Cols = Count - i) or (Rows = Count - i) then }

end; {for j := i to Count - 1 do } {* }

if A.Mass[i, i] = 0 then

for j := i + 1 to Count - 1 do

if A.Mass[j, i] <> 0 then

begin

Result := -Result; {* меняем строку }

Temp.Mass[0] := A.Mass[i]; {* на строку с }

A.Mass[i] := A.Mass[j]; {* первым }

A.Mass[j] := Temp.Mass[0]; {* ненулевым }

Break {* элементом }

end;

for j := i + 1 to Count - 1 do

if A.Mass[j, i] <> 0 then

begin

for k := i + 1 to Count - 1 do

A.Mass[j, k] := A.Mass[j, k] - A.Mass[i, k] * A.Mass[j, i] / A.Mass[i, i];

A.Mass[j, i] := 0

end

end; {Конец преобразования}

Рекомендуем скачать другие рефераты по теме: доклад 6 класс, здоровый образ жизни реферат.





Предыдущая страница реферата | 16  17  18  19  20  21  22  23  24  25  26 | Следующая страница реферата