Линейная алгебра. Теоремы и алгоритмы. Яцкин Н.И. - 581 стр.

UptoLike

ВУЗ: 

ИвГУ | Иваново

Составители: 

Рубрика: 

Прил. 1 Коды Maple-процедур 581
for i from 1 to n do
Delta[i]:=Determinant(SubMatrix(A,1..i,1..i));
# Вычисление угловых миноров.
print(evaln(Delta[i])=Delta[i]);
if i<n and Delta[i]=0 then
# Если очередной (не последний) угловой минор
# оказывается нулевым, то
# выдается сообщение об ошибке.
ERROR("Условия Якоби не выполнены!");
else
# В противном случае вычисляются
# диагональные элементы диагональной матрицы,
# задающей диагональный вид
# данной квадратичной формы.
mu[i]:=Delta[i]/Delta[i-1];
fi;
od;
for i from 1 to n do
print(evaln(mu[i])=mu[i]);
od;
B:=Matrix(n):
# Заготовка для диагональной матрицы.
# Далее следует ее заполнение.
for i from 1 to n do
B[i,i]:=mu[i]:
od:
T:=Matrix(n,symbol=`t`):
# Заготовка для унитреугольной
# матрицы перехода к диагональному виду.
# Далее следует ее заполнение
# (элементы выше главной диагонали
# остаются неопределенными).
for i from 1 to n do
T[i,i]:=1;
od:
for i from 2 to n do
for j from 1 to i-1 do
T[i,j]:=0;
od;
od; 
Прил. 1                    Коды Maple-процедур     581

for i from 1 to n do
  Delta[i]:=Determinant(SubMatrix(A,1..i,1..i));
  # Вычисление угловых миноров.

  print(evaln(Delta[i])=Delta[i]);

  if i

Страницы