ВУЗ:
Составители:
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
