Текст подпрограммы и версий ( Фортран )
aeh6r.zip , aeh6d.zip 		Тексты тестовых примеров ( Фортран )
taeh6r.zip , taeh6d.zip
Текст подпрограммы и версий ( Си )
aeh6r_c.zip , aeh6d_c.zip 		Тексты тестовых примеров ( Си )
taeh6r_c.zip , taeh6d_c.zip
Текст подпрограммы и версий ( Паскаль )
aeh6r_p.zip , aeh6e_p.zip 		Тексты тестовых примеров ( Паскаль )
taeh6r_p.zip , taeh6e_p.zip

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

Назначение

Вычисление всех собственных значений вещественной симметричной матрицы, заданной в компактной форме, с помощью QL - алгоритма.

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

Заданная матрица ортогональными преобразованиями подобия приводится к треугольной симметричной форме, собственные значения которой затем вычисляются с помощью неявного QL - алгоритма.

Уилкинсон, Pайнш. Справочник алгоритмов на языке АЛГОЛ. Линейная алгебра. М.: "Машиностроение", 1976.

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

    SUBROUTINE  AEH6R (N, NV, A, EV, R, IERR) 
   Параметры
N - порядок исходной матрицы (тип: целый);
NV - длина вектора А,  NV ≥ N(N + 1)/2 (тип: целый);
A - вещественный вектор длины NV, содержащий в своих первых N (N + 1)/2 компонентах исходную матрицу, заданную в компактной форме;
EV - вещественный вектор длины  N, содержащий вычисленные собственные значения, расположенные в неубывающем порядке;
R - вещественный вектор длины 2 * N, используемый как рабочий;
IERR - целочисленная переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы подпрограммы; IЕRR имеет тот же смысл, что и в подпрограмме АЕЕ2R. 
   Версии
AEH6D - вычисление с помощью QL - алгоpитма всех собственных значений вещественной симметричной матрицы, заданной в компактной форме, с удвоенной точностью. 
   Вызываемые подпрограммы
AEE2R - вычисление всех собственных значений вещественной симметричной трехдиагональной матрицы при помощи неявного QL - алгоpитма.
AFH5R - приведение вещественной симметричной матрицы, заданной в компактной форме, к симметричной трехдиагональной матрице ортогональными преобразованиями подобия.
UTAE10 - подпрограмма выдачи диагностических сообщений при работе подпрограмм АЕЕ2R и АЕЕ2D. 
   Замечания по использованию
  1. 

Подпрограмма АЕН6R не сохраняет исходный вектор  А.

  2.  В подпрограмме АЕН6D параметры А, ЕV, R имеют тип DОUВLЕ РRЕСISIОN.

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

      DIMENSION  A(6), EV(3), R(6)
      DATA  A /294., - 180., - 25., 108., 240., 231./
      CALL  AEH6R (3, 6, A, EV, R, IERR)

Результаты:

       IERR  =  0

       EV  =  (- 250., 375., 375.)