Cтраница 2
Вычислительная часть программы оформлена в виде процедуры UNITIM, входными параметрами которой являются порядок матрицы Р, матрица U, точность расчета EPS. Выходным параметром процедуры является матрица L размерности Р х 2, строки которой содержат действительные и мнимые части найденных собственных значений исходной матрицы. В процедуре UNITIM используются две процедуры SDM и СОМР, первая из которых реализует сложение и вычитание матриц, а вторая - преобразование комплексных чисел из алгебраической в тригонометрическую форму и обратно. [16]
Вычислительная часть программы оформлена в виде процедуры UNITIM, входными параметрами которой являются порядок матрицы Р, матрица [ 7, точность расчета EPS. Выходным параметром процедуры является матрица L размерности Р X 2, строки которой содержат действительные и мнимые части найденных собственных значений исходной матрицы. В процедуре UNITIM используются две процедуры SDM и СОМР, первая из которых реализует сложение и вычитание матриц, а вторая - преобразование комплексных чисел из алгебраической в тригонометрическую форму и обратно. [17]
Интерполяционная таблица формируется подпрограммой TAB, в которую исходная матрица В и интервал поиска собственных значений ALC, AL9 передаются из основной программы с помощью неименованного COMMON-блока. Левая часть характеристического уравнения вычисляется с помощью подпрограммы-функции PN, в которую через COMMON-блок с именем XFAL передаются узлы и коэффициенты полинома Ньютона, найденные собственные значения, номер искомого корня и степень полинома. [18]