Текст подпрограммы и версий ( Фортран )
zp38r.zip  zp38d.zip  		Тексты тестовых примеров ( Фортран )
tzp38r.zip  tzp38d.zip 
Текст подпрограммы и версий ( Си )
zp38r_c.zip  zp38d_c.zip  		Тексты тестовых примеров ( Си )
tzp38r_c.zip  tzp38d_c.zip 
Текст подпрограммы и версий ( Паскаль )
zp38r_p.zip  zp38e_p.zip  		Тексты тестовых примеров ( Паскаль )
tzp38r_p.zip  tzp38e_p.zip 

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

Назначение

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

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

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

        C(X)  =  XK*( C(1) + C(2)*X + ... + C(M)*X(M - 1) + ... )  ,

такого, что  С (Х) = А (Х) * В (Х), где  А (Х) и  В (Х) - заданные полиномы степени NР и L соответственно 

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

        B(X)  =   B(1) + B(2)*X + B(3)*X2 + ... + B(L + 1)*XL  ,

у которых коэффициенты  А ( I ),  I = 1, 2, 3, ..., NР + 1 и   В ( I ),  I = 1, 2, 3, ..., L + 1 являются вещественными. Здесь К = R + S, где R и S представляют собой число первых коэффициентов полиномов  А (Х) и  В (Х), равных 0.

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

    SUBROUTINE  ZP38R (NP, A, L, B, M, C, K) 
   Параметры
NP - заданная степень полинома  А (Х) (тип: целый);
A - вещественный вектор длины NР + 1, содержащий коэффициенты заданного полинома  А (Х) в порядке возрастания степеней;
L - заданная степень полинома  В (Х) (тип: целый);
B - вещественный вектор длины L + 1, содержащий коэффициенты заданного полинома  В (Х) в порядке возрастания степеней;
M - заданное число искомых коэффициентов полинома  С (Х) (тип: целый);
C - вещественный вектор длины М, содержащий вычисленные коэффициенты полинома  С (Х) в порядке возрастания степеней;
K - целая переменная, значение которой в результате работы подпрограммы полагается равным R + S. 
   Версии
ZP38D - вычисление произведения двух полиномов с вещественными коэффициентами с повышенной точностью. При этом векторы А, В и С должны иметь тип DОUВLЕ РRЕСISIОN. 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию:  нет

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

       DIMENSION  A(3), B(3), C(5)
       DATA  A(1), A(2), A(3) /2., - 3., 1./
       DATA  B(1), B(2), B(3) /12., - 7., 1./
       NP = 2
       L = 2
       M = 4
       CALL  ZP38R (NP, A, L, B, M, C, K)

Результаты:

       С(1) = 24. ,    С(2) = - 50. ,    С(3) = 35. ,
       C(4) = - 10. ,  C(5) = 1. ,
       K = 0