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

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

    # севере-западный блок размера ms на ms
    # в полученной матрице J
    # должен совпадать с матрицей JS.
    if not Equal(JS,SubMatrix(J,1..ms,1..ms)) then
      ERROR(`Равенство JS=SubMatrix(G,1..ms,1..ms)
             не выполняется!`);
    fi;

    print(str);print(str);

    # Печать окончательных результатов.
    print(evaln(J)=J);
    print(evaln(T)=T);
    print(evaln(det('T'))=dt);
    print(str);print(str);

fi;

#   "Сценарная часть" работы процедуры завершена.
#   Далее следует завершение работы всей процедуры
#   с возвращением (в виде, допускающем дальнейшее использование)
#   1) матрицы J - (частичной) ж.н.ф. для A;
#   2) матрицы T - перехода от исходного базиса
#      к (частично) жорданову.

RETURN(J,T);

end proc;