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

UptoLike

ВУЗ: 

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

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

Рубрика: 

530 Коды Maple-процедур Прил. 1
pr:=[seq(0,k=1..rA)];
# Заготовка для списка номеров главных неизвестных.
for i from 1 to rA do
j:=1;
while JA[i,j]=0 do
j:=j+1;
od;
pr[i]:=j;
od;
# Сформирован список номеров главных неизвестных.
fr:=[];
for k from 1 to n do
if not member(k,pr) then
fr:=[fr[],k];
fi;
od;
# Сформирован список номеров свободных неизвестных.
FA:=Matrix(n,n-rA);
# Заготовка для фундаментальной матрицы,
# соответствующей с.л.у. Ax=0.
# Если rA=n, то W=O, и матрица FA пуста.
# (Выводятся пустой базис и пустой список
# номеров свободных неизвестных.)
if rA<n then
# Неособый случай. Подпространство W решений
# однородной с.л.у. Ax=0 нетривиально.
for h from 1 to n-rA do
FA[ fr[h], h ]:=1;
for i from 1 to rA do
FA[ pr[i], h ]:=-GA[ i, fr[h] ];
od;
od;
# Каждой свободной неизвестной
# соответствует столбец
# в фундаментальной матрице FA -
# базисное частное решение
# однородной с.л.у. Ax=0,
# получающееся, если эта неизвестная равна 1.
fi;
fi;
sys:=GenerateEquations(JA,[seq(x[k],k=1..n)]);
# Генерация однородной с.л.у. по матрице JA. 
530                       Коды Maple-процедур           Прил. 1

  pr:=[seq(0,k=1..rA)];
  # Заготовка для списка номеров главных неизвестных.

  for i from 1 to rA do
    j:=1;
    while JA[i,j]=0 do
      j:=j+1;
    od;
    pr[i]:=j;
  od;
  # Сформирован список номеров главных неизвестных.

  fr:=[];
  for k from 1 to n do
    if not member(k,pr) then
      fr:=[fr[],k];
    fi;
  od;
  # Сформирован список номеров свободных неизвестных.

  FA:=Matrix(n,n-rA);
  # Заготовка для фундаментальной матрицы,
  # соответствующей с.л.у. Ax=0.
  # Если rA=n, то W=O, и матрица FA пуста.
  # (Выводятся пустой базис и пустой список
  # номеров свободных неизвестных.)

  if rA

Страницы