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

UptoLike

ВУЗ: 

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

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

Рубрика: 

Прил. 1 Коды Maple-процедур 585
# Далее следует ее заполнение.
for i from 1 to n do
if i<=s then
T2[i,i]:=1/sqrt(B1[i,i]);
elif i>s and i<=r then
T2[i,i]:=1/sqrt(-B1[i,i]);
else
T2[i,i]:=1;
fi;
od:
B2:=Transpose(T2).B1.T2;
# Вычисление матрицы,
# отвечающей нормальному виду
# квадратичной формы.
print(`T2=`,T2,`B2=`,B2);
RETURN(B2,T1.T2,[s,t]);
# Возвращаются: матрица нормального вида,
# матрица перехода
# от диагонального вида к нормальному,
# сигнатура (индексы инерции)
# данной квадратичной формы.
end proc;
> save Quadro,"F:/MaplePackages/Quadro.m";
# Сохранение пакета. 
Прил. 1                 Коды Maple-процедур   585

# Далее следует ее заполнение.
for i from 1 to n do
  if i<=s then
    T2[i,i]:=1/sqrt(B1[i,i]);
  elif i>s and i<=r then
    T2[i,i]:=1/sqrt(-B1[i,i]);
  else
    T2[i,i]:=1;
  fi;
od:

B2:=Transpose(T2).B1.T2;
# Вычисление матрицы,
# отвечающей нормальному виду
# квадратичной формы.

print(`T2=`,T2,`B2=`,B2);


RETURN(B2,T1.T2,[s,t]);
# Возвращаются: матрица нормального вида,
# матрица перехода
# от диагонального вида к нормальному,
# сигнатура (индексы инерции)
# данной квадратичной формы.

end proc;


> save Quadro,"F:/MaplePackages/Quadro.m";
# Сохранение пакета.

Страницы