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

544                       Коды Maple-процедур          Прил. 1

# Печать массива q.
for i from 1 to s do
  for k from 1 to l[i] do
    print(evaln(q[i][k])=q[i][k]);
  od;

  print(str);
od;

print(str);

#   Задание массива DIAGR столбчатых диаграмм
#   для матриц B[i].
#   Элементами этого массива являются
#   матрицы DIAGR[i] {i=1,...,s).

#   Для их визуализации используются
#   матрицы RIS[i] с противоположным порядком строк.
#   (Это связано с тем, что строки в матрице
#   нумеруются сверху вниз,
#   а в столбчатой диаграмме - наоборот.)

# Сначала диаграммы даются без нумерации векторов
# (заполняются звездочками),
# а затем - с нумерацией по принципу:
# столбцы нумеруются слева направо,
# векторы в столбцах - снизу вверх.
DIAGR:=array(1..s):RIS:=array(1..s);

# Цикл по номеру собственного значения.
for i from 1 to s do

    # Задание (заполненных "пустыми словами" ``)
    # матриц DIAGR[i],
    # каждая из которых содержит
    # l[i] (= показатель стабилизации) строк
    # и d[i][1] (= первый дефект) столбцов.
    DIAGR[i]:=Matrix(l[i],d[i][1],fill=``):

    # Заполнение звездочками тех позиций в DIAGR[i],
    # которые относятся к столбчатой диаграмме.

    # Цикл по номеру строки в матрице DIAGR[i].
    for k from 1 to l[i] do

      # Цикл по номеру столбца в DIAGR[i].
      for j from 1 to p[i][k] do
        # Звездочка в j-м столбце k-й строки матрицы
        # проставляется, если 1<=j<=p[i][k].
        DIAGR[i][k,j]:=`*`;
      od:
    od: