|
Текст подпрограммы и версий ( Фортран ) zp53r.zip zp53d.zip |
Тексты тестовых примеров ( Фортран ) tzp53r.zip tzp53d.zip |
|
Текст подпрограммы и версий ( Си ) zp53r_c.zip zp53d_c.zip |
Тексты тестовых примеров ( Си ) tzp53r_c.zip tzp53d_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) zp53r_p.zip zp53e_p.zip |
Тексты тестовых примеров ( Паскаль ) tzp53r_p.zip tzp53e_p.zip |
Вычисление значений полинома с вещественными коэффициентами и его производных от заданного вещественного аргумента
Пусть задан полином степени NP с вещественными коэффициентами
A(1) + A(2)*X + A(3)*X2 + ... + A(NP + 1)*XNP
Подпрограмма ZP53R вычисляет значения полинома и его ND - 1 производных от заданного вещественного переменного и помещает эти значения в элементах массива PD (1), PD (2), ..., PD (ND) соответственно
SUBROUTINE ZP53R (NP, A, X, ND, PD)
Параметры
| NP - | заданная степень полинома (тип: целый); |
| A - | вещественный вектор длины NP + 1, содержащий коэффициенты исходного полинома в порядке возрастания степеней; |
| X - | заданное значение аргумента (тип: вещественный); |
| ND - | длина вектора PD, ND > 1 (тип: целый); |
| PD - | вещественный вектор длины ND, содержащий в PD (1) вычисленное значение полинома в точке X, а в PD (2), ..., PD (ND) - вычисленные значения его производных порядков от 1 до ND - 1 в точке X |
Версии
| ZP53D - | вычисление значений полинома с вещественными коэффициентами и его производных от заданного вещественного аргумента в режиме удвоенной точности; при этом параметры A, X и PD должны иметь тип DOUBLE PRECISION |
Вызываемые подпрограммы: нет
Замечания по использованию: нет
DIMENSION A(5), PD(6)
DATA A(1), A(2), A(3), A(4), A(5) /1., 2., 3., 4., 5./
NP = 4
ND = 6
X = 1.0
CALL ZP53R (NP, A, X, ND, PD)
Результаты: PD = ( 15, 40, 90, 144, 120, 0 )