Текст подпрограммы и версий ( Фортран )
zp40r.zip  zp40d.zip  		Тексты тестовых примеров ( Фортран )
tzp40r.zip  tzp40d.zip 
Текст подпрограммы и версий ( Си )
zp40r_c.zip  zp40d_c.zip  		Тексты тестовых примеров ( Си )
tzp40r_c.zip  tzp40d_c.zip 
Текст подпрограммы и версий ( Паскаль )
zp40r_p.zip  zp40e_p.zip  		Тексты тестовых примеров ( Паскаль )
tzp40r_p.zip  tzp40e_p.zip 

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

Назначение

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

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

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

       B(X)  =  B(1) + B(2)*Y + B(3)*Y2 + ... + B(NP + 1)*YNP ,

получающегося из данного полинома с вещественными коэффициентами 

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

линейным изменением аргумента Y = u * Х + v.

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

    SUBROUTINE  ZP40R (NP, A, B, U, V, RAB) 
   Параметры
NP - заданная степень полинома (тип: целый);
A - вещественный вектор длины NР + 1, содержащий коэффициенты заданного полинома  А (Х) в порядке возрастания степеней;
B - вещественный вектор длины NР + 1, содержащий вычисленные коэффициенты полинома  В (Х) в порядке возрастания степеней;
U, V - заданные значения  u  и  v (тип: вещественный);
RAB - вещественный вектор длины NР + 1, используемый в подпрограмме как рабочий. 
   Версии
ZP40D - вычисление коэффициентов полинома, получающегося из данного полинома с вещественными коэффициентами линейным изменением аргумента с повышенной точностью. При этом векторы А, В, RАВ, а также U и V должны иметь тип DОUВLЕ РRЕСISIОN. 
   Вызываемые подпрограммы
ZP30R - вычисление значения полинома с вещественными коэффициентами от заданного вещественного аргумента. Используется в подпрограмме ZР40R.
ZP30D - вычисление значения полинома с вещественными коэффициентами от заданного вещественного аргумента с повышенной точностью. Используется в подпрограмме ZР40D. 
   Замечания по использованию:  нет

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

       DIMENSION  A(3), B(3), RAB(3)
       DATA  A(1), A(2), A(3) /1., 2., 3./
       NP = 2
       U = 2.
       V = 3.
       CALL  ZP40R (NP, A, B, U, V, RAB)

Результаты:    В(1) = 4.75 ,   В(2) = - 3.5 ,   В(3) = 0.75