ВУЗ:
Составители:
79
Начало
ввод n,
матрицы А
u := 0
i: = 1, n
j
:= i
;
k:=0
2 3
j
≤
n
да
нет
61],[ −> EjiA
нет
j := j+1
k: = j
j: = n+1
да
k< >i
j
:= i, n
C := A[k, j]
A[i, j] := A[k, j]
A[k, j] := C
u: = u+1
1
да
нет
Приведение матрицы
к треугольному виду
k=0
нет
Печать ‘в
столбце’,i,’все
элементы=0’
Exit
Перестановка
строк i и k
да
80
1
k
:= i+1, n
2
C := -A
[
k
,
i
]
/A
[
i
,
i
]
j
:= i
,
n
A[k, j] := A[k, j]+C*A[i, j]
Получение нулей ниже
главной диагонали
3
det := 1
i:=1, n
det := det*A
[
i
,
i
]
u mod 2<>0
det := -det
да
нет
Печать
‘оп
р
е
д
елитель=’
,
det
конец
Рис. 41. Блок-схема вычисления определителя матрицы по методу
Гаусса с выбором главного элемента в столбце.
Начало 1 Получение нулей ниже ввод n, главной диагонали матрицы А k := i+1, n 2 u := 0 Приведение матрицы C := -A[k, i]/A[i, i] к треугольному виду 2 i: = 1, n 3 j := i, n j: = i; k:=0 A[k, j] := A[k, j]+C*A[i, j] нет j≤ n 3 да нет A[i , j ] > 1E − 6 det := 1 да j := j+1 k: = j i := 1, n j: = n+1 да k=0 нет det := det*A[i, i] нет k< >i Печать в нет Перестановка столбце,i,все u mod 2<>0 строк i и k да элементы=0 да j: = i, n Exit det := -det C := A[k, j] A[i, j] := A[k, j] Печать A[k, j] := C определитель=, det u: = u+1 конец Рис. 41. Блок-схема вычисления определителя матрицы по методу 1 Гаусса с выбором главного элемента в столбце. 79 80