Текст подпрограммы и версий ( Фортран )
zp36r.zip  zp36d.zip 
Тексты тестовых примеров ( Фортран )
tzp36r.zip  tzp36d.zip 
Текст подпрограммы и версий ( Си )
zp36r_c.zip  zp36d_c.zip 
Тексты тестовых примеров ( Си )
tzp36r_c.zip  tzp36d_c.zip 
Текст подпрограммы и версий ( Паскаль )
zp36r_p.zip  zp36e_p.zip 
Тексты тестовых примеров ( Паскаль )
tzp36r_p.zip  tzp36e_p.zip 

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

Назначение

Вычисление синуса и косинуса от полинома с вещественными коэффициентами.

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

ZР36R вычисляет М первых коэффициентов полиномов

       B(X)  =  B(1) + B(2)*X + B(3)*X2 + ... + B(M)*X(M - 1) + ... 
  и
       C(X)  =  C(1) + C(2)*X + C(3)*X2 + ... + C(M)*X(M - 1) + ...  , 

таких, что  В (Х) = sin А (Х) и  С (Х) = cos А (Х), где  А (Х) - заданный полином степени NР

       A(X)  =  A(1) + A(2)*X + A(3)*X2 + ... + A(NP + 1)*XNP  , 

у которого коэффициенты  А ( I ),  I = 1, 2, 3, ..., NР являются вещественными.

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

    SUBROUTINE  ZP36R (NP, A, M, B, C) 

Параметры

NP - заданная степень полинома  А (Х) (тип: целый);
A - вещественный вектор длины NР + 1, содержащий коэффициенты заданного полинома  А (Х) в порядке возрастания степеней;
M - заданное число искомых коэффициентов полиномов  В (Х) и  С (Х) (тип: целый);
B, C - вещественные векторы длины М, содержащие вычисленные коэффициенты полиномов  В (Х) и  С (Х) в порядке возрастания степеней.

Версии

ZP36D - вычисление синуса и косинуса от полинома с вещественными коэффициентами с повышенной точностью. При этом векторы А, В и С должны иметь тип DОUВLЕ РRЕСISIОN.

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

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

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

       DIMENSION  A(5), B(4), C(4)
       DATA  A(1), A(2), A(3), A(4), A(5) /24., - 50., 35., - 10., 1./
       NP = 4
       M = 4
       CALL  ZP36R (NP, A, M, B, C)

Результаты:
 
       В(1) = - 0.905578 ,      В(2) = - 21.208950 ,
       B(3) = 1146.819218 ,   B(4) = 7248.058730

       C(1) = 0.424179 ,         C(2) = - 45.278918 , 
       C(3) = - 498.528516 ,   C(4) = 19599.473354