# Вторая пачка часть 40

.N*2]OF REAL; MSS = ARRAY [1..N] OF MST; VAR A1 : MSS; I, J, M, K : INTEGER; H : REAL; BEGIN FOR I := 1 TO N DO { ФОРМИРОВАНИЕ ВСПОМОГАТЕЛЬНОЙ МАТРИЦЫ А1 } FOR J := 1 TO N*2 DO IF J<= N THEN A1 [I,J] := A[I,J] ELSE IF I=J-N THEN A1 [I,J] := 1.0 ELSE A1[I,J] :=0.0; TOL := 0; DET := 1; { НАЧАЛЬНЫЕ ЗНАЧЕНИЯ СОNSТ-АНТ} IF N > 50 THEN BEGIN TOL := 2; EXIT; END; IF ABS(A1[N,N]) < 1.0E-10 THEN BEGIN TOL := 1; EXIT; END; FOR I := 1 TO N DO BEGIN IF ABS(A1[I,I])<1.0E-10 THEN BEGIN TOL:= 1; EXIT; END; H := A1[I,I]; DET := DET * H; FOR J := I TO N*2 DO A1[I,J] := A1[I,J] / H; IF I<N THEN FOR J := I+1 TO N DO BEGIN H := A1[J,I]; FOR K := 1 TO N*2 DO A1[J,K] := A1 [J,K] - H*A1[I,K]; END; END; FOR I := N DOWNTO 1 DO FOR K := I-1 DOWNTO 1 DO BEGIN H := A1[K,I]; FOR J:=1 TO 2*N DO A1[K,J] := A1[K,J] - A1[I,J]*H; END; FOR I := 1 TO N DO FOR J := N+1 TO N*2 DO V [I,J-N] := A1[I,J]; END. Формальные параметры про¬цедуры. Входные: А (тип re¬al) - матрица, для которой ищется об¬рат¬¬¬ная; N (тип in¬te¬ger) - размерность матрицы (матрица долж¬на быть квад¬ратной, т.е. размером N?N). Вы¬ход¬ные: tоl (тип re¬al) - це¬лое чис¬ло, ко¬то¬¬¬рое при нор¬мальном завершении про¬¬це¬ду¬ры равно 0. Ес¬ли при об¬ра¬щении мат¬ри-цы один из ди¬а¬го¬наль¬¬ных элементов стал равен 0, то про¬¬¬це¬¬дура воз¬вра¬ща¬¬ет зна¬чение дан¬но¬го па-ра¬мет¬¬¬¬ра 1; V (тип real) - об¬ратная для А мат¬ри¬ца; SUМ- оп¬ре¬¬¬де¬ли¬тель мат¬ри¬цы. В процедуре GAUS_OBR вводится вспо-мо¬га¬тель¬ная ма¬т¬¬ри¬ца А1, размером [N?2*N], в ко¬то¬рой пер¬вая по¬ловина рав¬на введенной мат-ри¬це А, а вто¬рая (от N+1 до N*2) - еди¬¬ни¬чной. Вы¬пол¬няя пре¬об¬ра¬зо¬вания над пер¬вой по¬ло¬ви-ной матрицы А1 (от 1 до N) для при¬ве¬де¬ния ее к виду единичной матрицы, по¬лу¬ча¬ем во второй половине мат¬ри¬цы А1 матрицу, об¬рат¬ную по от-ношению к А, т

Индекс
Элементарные функции Линейные уравнения Нелинейные уравнения Случайные числа