Текст подпрограммы и версий ( Фортран )
aee4r.zip , aee4d.zip 		Тексты тестовых примеров ( Фортран )
taee4r.zip , taee4d.zip
Текст подпрограммы и версий ( Си )
aee4r_c.zip , aee4d_c.zip 		Тексты тестовых примеров ( Си )
taee4r_c.zip , taee4d_c.zip
Текст подпрограммы и версий ( Паскаль )
aee4r_p.zip , aee4e_p.zip 		Тексты тестовых примеров ( Паскаль )
taee4r_p.zip , taee4e_p.zip

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

Назначение

Вычисление К максимальных собственных значений симметричной трехдиагональной матрицы.

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

Подпрограма АЕЕ4R вычисляет матодом бисекций К максимальных собственных значений симметричной трехдиагональной матрицы.

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

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

    SUBROUTINE  AEE4R ( A, B, N, RU, W, K) 
   Параметры
A - вещественный вектор длины N, содержащий диагональные элементы симметричной трехдиагональной матрицы;
B - вещественный вектор длины N, содержащий в последних N - 1 компонентах поддиагональные элементы симметричной трехдиагональной матрицы;
N - заданный порядок симметричной трехдиагональной матрицы (тип: целый);
RU - вещественный вектор длины К, используемый как рабочий;
W - вещественный вектор длины К, содержащий вычисленные собственные значения в неубывающем порядке;
K - заданое число искомых собственных значений (тип: целый). 
   Версии
AEE4D - вычисление к максимальных собственных значений симметричной трехдиагональной матрицы, заданной с двойной точностью. 
   Вызываемые подпрограммы :   нет 
   Замечания по использованию
  1. 

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

  2. 

В подпрограмме АЕЕ4D параметры А, В, RU, W имеют тип DОUВLЕ РRЕСISIОN.

  3.  Подпрограмма АЕЕ4R может быть использована и для вычисления К максимальных собственных значений комплексной эрмитовой трехдиагональной матрицы. Для этого нужно предварительно привести комплексную эрмитовую трехдиагональную матрицу к вещественной трехдиагональной симметричной матрице преобразованиями подобия (например, обратиться к подпрограмме АFE0С).
  4.  Подпрограмма АЕE4R использует служебную подпрограмму AMA1R.
Подпрограмма АЕE4D использует служебную подпрограмму AMA1D.

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

      DIMENSION  A(50), B(50), RU(5), S(5)
      DATA  A /1., 49., 48*0.0/
      DATA  B /0., 7., 48*0.0/
      N = 50
      K = 5
      CALL  AEE4R (A, B, N, RU, S, K)

Результаты:

      S(1)  =  S(2)  =  S(3)  =  S(4)  =  -2.55*10-12, 
      S(5)  =  50 .