Текст подпрограммы и версий ( Фортран )
zp45r.zip  zp45d.zip  zp46r.zip  zp46d.zip 
Тексты тестовых примеров ( Фортран )
tzp45r.zip  tzp45d.zip  tzp46r.zip  tzp46d.zip 
Текст подпрограммы и версий ( Си )
zp45r_c.zip  zp46r_c.zip 
Тексты тестовых примеров ( Си )
tzp45r_c.zip  tzp46r_c.zip 
Текст подпрограммы и версий ( Паскаль )
zp45r_p.zip  zp45e_p.zip  zp46r_p.zip  zp46e_p.zip 
Тексты тестовых примеров ( Паскаль )
tzp45r_p.zip  tzp45e_p.zip  tzp46r_p.zip  tzp46e_p.zip 

Подпрограмма:  ZP45R (версия: ZP46R)

Назначение

Вычисление значений ортогональных полиномов Чебышева.

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

Подпрограмма вычисляет значения ортогональных полиномов Чебышева для заданного аргумента Х по рекуррентной формуле:

          T(K)  =  2*X*T(K-1) - T(K-2)  ,   T(1) = 1  ,   T(2) = X  , 

где  Т (К) - значение полинома порядка  К в точке  Х.

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

    SUBROUTINE  ZP45R (X, N, T) 

Параметры

X - заданное значение аргумента (тип: вещественный);
N - длина вектора  Т (тип: целый);
T - вещественный вектор длины  N, в котором содержатся значения полиномов в порядке возрастания степеней, начиная с полинома степени ноль.

Версии

ZP45D - вычисление значений ортогональных полиномов Чебышева для данных двойной точности.
ZP46R -
ZP46D  
вычисление значений смещенных ортогональных полиномов Чебышева. для данных одинарной и двойной точности.

Вызываемые подпрограммы: нет

Замечания по использованию

  Смещенные ортогональные многочлены Чебышева ТS (К, Х) определяются по формуле ТS (К, X) = Т (К, 2 * Х - 1), где  К - порядок полинома,  Х - аргумент и  Т (К, Х) - значение обычного полинома Чебышева порядка  К в точке  Х: Т (К, Х) = cos (К * АRССОS (Х)).

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

       DIMENSION  T(3)
       X = 0.5
       N = 3
       CALL  ZP45R (X, N, T)

Результаты:     T = ( 1., 0.5, - 0.5 )