Текст подпрограммы и версий ( Фортран )
aej2r.zip , aej2d.zip 		Тексты тестовых примеров ( Фортран )
taej2r.zip , taej2d.zip
Текст подпрограммы и версий ( Си )
aej2r_c.zip , aej2d_c.zip 		Тексты тестовых примеров ( Си )
taej2r_c.zip , taej2d_c.zip
Текст подпрограммы и версий ( Паскаль )
aej2r_p.zip , aej2e_p.zip 		Тексты тестовых примеров ( Паскаль )
taej2r_p.zip , taej2e_p.zip

Подпрограмма:  AEJ2R

Назначение

Вычисление всех собственных значений вещественной матрицы Якоби при помощи неявного Q*L - алгоpитма.

Математическое описание

Исходная вещественная трехдиагональная матрица обладает следующим свойством: произведения пар соответствующих элементов побочных диагоналей все неотрицательны. Вычисленные собственные значения располагаются в возрастающем порядке.

Дж.Х.Уилкинсон, Алгебраическая проблема собственных значений, "Наука", М., 1970.

Использование 

    SUBROUTINE  AEJ2R ( N, A, EV, RAB, IERR) 
   Параметры
N - порядок исходной матрицы (тип: целый);
A - вещественный двумерный массив размерности N на 3, содержащий в последних N - 1 компонентах первого столбца элементы нижней диагонали, во втором столбце - элементы главной диагонали, в первых N - 1 компонентах третьего столбца - элементы верхней диагонали;
EV - вещественный одномерный массив размерности N, содержащий вычисленные собственные значения в возрастающем порядке;
RAB - вещественный одномерный массив размерности 2*N, используемый как рабочий;
IERR - целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы подпрограммы, при этом:
  а) значение IЕRR больше N, если не все попарные произведения соответствующих элементов побочных диагоналей неотрицательны,
  б) значение IЕRR полагается равным номеру собственного значения, для вычисления которого потребовалось более 30 итераций, при этом собственные значения с индексами 1, 2, ..., IЕRR - 1 вычислены правильно и расположены в возрастающем порядке, но они не обязательно являются самыми меньшими из всех N собственных значений. 
   Версии
AEJ2D - вычисление всех собственных значений вещественной матрицы Якоби с повышенной точностью. Массивы А,ЕV, RАВ имеют тип DОUВLЕ РRЕСISIОN. 
   Вызываемые подпрограммы
UTAE10 - подпрограмма выдачи диагностических сообщений при работе подпрограмм АЕJ2R и АЕJ2D. 
   Замечания по использованию
  Подпрограмма АЕJ2R сохраняет исходную информацию.

Пример использования 

      DIMENSION  A(5, 3), EV(5), RAB(10)
      DATA  A /0., 4*1., 5*0., 4*1., 0./
      N = 5
      CALL  AEJ2R (N, A, EV, RAB, IERR)

Результаты:

      EV(1)  =  -1.7320508 , 
      EV(2)  =  -1.0 , 
      EV(3)  =  6.28*10-12 , 
      EV(4)  =  1.0 , 
      EV(5)  =  1.7320508 ;

      IERR  =  0