|
Текст подпрограммы и версий ( Фортран ) zp37r.zip zp37d.zip |
Тексты тестовых примеров ( Фортран ) tzp37r.zip tzp37d.zip |
|
Текст подпрограммы и версий ( Си ) zp37r_c.zip zp37d_c.zip |
Тексты тестовых примеров ( Си ) tzp37r_c.zip tzp37d_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) zp37r_p.zip zp37e_p.zip |
Тексты тестовых примеров ( Паскаль ) tzp37r_p.zip tzp37e_p.zip |
Возведение в степень полинома с вещественными коэффициентами.
ZР37R вычисляет М первых коэффициентов полинома
B(X) = XS*(B(1) + B(2)*X + ... + B(M)*X(M - 1) + ...) ,
такого, что В (Х) = (А (Х))R, где А (Х) - заданный полином, степени NР
A(X) = A(1) + A(2)*X + A(3)*X2 + ... + A(NP + 1)*XNP ,
у которого коэффициенты А ( I ), I = 1, 2, 3, ..., NР + 1 являются вещественными. Здесь S = К * R, где К - число первых коэффициентов полинома А (Х), равных 0.
SUBROUTINE ZP37R (NP, A, M, B, R, S)
Параметры
| NP - | заданная степень полинома А (Х) (тип: целый); |
| A - | вещественный вектор длины NР + 1, содержащий коэффициенты заданного полинома А (Х) в порядке возрастания степеней; |
| M - | заданное число искомых коэффициентов полинома В (Х) (тип: целый); |
| B - | вещественный вектор длины М, содержащий вычисленные коэффициенты полинома В (Х) в порядке возрастания степеней; |
| R - | заданный показатель степени (тип: вещественный); |
| S - | вещественная переменная, значение которой в результате работы подпрограммы полагается равным К * R. |
Версии
| ZP37D - | возведение в степень полинома с вещественными коэффициентами с повышенной точностью. При этом векторы А и В должны иметь тип DОUВLЕ РRЕСISIОN. |
Вызываемые подпрограммы: нет
Замечания по использованию: нет
DIMENSION A(5), B(4)
DATA A(1), A(2), A(3), A(4), A(5) /24., - 50., 35., - 10., 1./
NP = 4
M = 4
R = 2
CALL ZP37R (NP, A, M, B, R, S)
Результаты:
В(1) = 576. , В(2) = - 2400. , В(3) = 4180. , B(4) = - 3980. ,
S = 0