Текст подпрограммы и версий ( Фортран )
aee5r.zip , aee5d.zip 		Тексты тестовых примеров ( Фортран )
taee5r.zip , taee5d.zip
Текст подпрограммы и версий ( Си )
aee5r_c.zip , aee5d_c.zip 		Тексты тестовых примеров ( Си )
taee5r_c.zip , taee5d_c.zip
Текст подпрограммы и версий ( Паскаль )
aee5r_p.zip , aee5e_p.zip 		Тексты тестовых примеров ( Паскаль )
taee5r_p.zip , taee5e_p.zip

Функция:  AEE5R

Назначение

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

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

Пусть все собственные значения заданной симметричной трехдиагональной матрицы порядка N пронумерованы в порядке неубывания 

         λ1 ≤ λ2 ≤  ...  ≤ λN-1 ≤ λN

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

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

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

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

Подпрограмма - функция АЕЕ5R не сохраняет исходные векторы А и В.

  2. 

В функции АЕЕ5D праметры А, В имеют тип DОUВLЕ  РRЕСISIОN.

  3.  Подпрограмма - функция АЕЕ5R может быть использована и для вычисления собственного значения комплексной эрмитовой трехдиагональной матрицы по заданному номеpу. Для этого нужно предварительно привести комплексную эрмитову трехдиагональную матрицу к вещественной симметричной трехдиагональной матрице преобразованиями подобия (например, обратиться к подпрограмме АFЕ0С).
  4.  Подпрограмма - функция АЕE5R использует служебную подпрограмму AMA1R.
Подпрограмма - функция АЕE5D использует служебную подпрограмму AMA1D.

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

      DIMENSION  A(50), B(50)
      DATA  A /1., 49., 48*0.0/
      DATA  B /0., 7., 48*0./
      K = 15
      N = 50
      S = AEE5R (A, B, N, K)

Результаты:

      S  =  -2.55*10-12