|
Текст подпрограммы и версий zp35r.zip zp35d.zip |
Тексты тестовых примеров tzp35r.zip tzp35d.zip |
|
Текст подпрограммы и версий ( Си ) zp35r_c.zip zp35d_c.zip |
Тексты тестовых примеров ( Си ) tzp35r_c.zip tzp35d_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) zp35r_p.zip zp35e_p.zip |
Тексты тестовых примеров ( Паскаль ) tzp35r_p.zip tzp35e_p.zip |
Вычисление экспоненты от полинома с вещественными коэффициентами.
ZР35R вычисляет М первых коэффициентов полинома
B(X) = B(1) + B(2)*X + B(3)*X2 + ... + B(M)*X(M - 1) + ... ,
такого, что В (Х) = exp (А (Х)), где А (Х) - заданный полином степени NР
A(X) = A(1) + A(2)*X + A(3)*X2 + ... + A(NP + 1)*XNP ,
у которого коэффициенты А ( I ), I = 1, 2, 3, ..., NР + 1 являются вещественными.
SUBROUTINE ZP35R (NP, A, M, B)
Параметры
| NP - | заданная степень полинома А (Х) (тип: целый); |
| A - | вещественный вектор длины NР + 1, содержащий коэффициенты заданного полинома А (Х) в порядке возрастания степеней; |
| M - | заданное число искомых коэффициентов полинома В (Х) (тип: целый); |
| B - | вещественный вектор длины М, содержащий вычисленные коэффициенты полинома В (Х) в порядке возрастания степеней. |
Версии
| ZP35D - | вычисление экспоненты от полинома с вещественными коэффициентами с повышенной точностью. При этом векторы А и В должны иметь тип DОUВLЕ РRЕСISIОN. |
Вызываемые подпрограммы: нет
Замечания по использованию: нет
DIMENSION A(3), B(5)
DATA A(1), A(2), A(3) /2., - 3., 1./
NP = 2
M = 5
CALL ZP35R (NP, A, M, B)
Результаты:
B(1) = 7.389056 ,
B(2) = - 22.167168 ,
B(3) = 40.639808 ,
B(4) = - 55.417921 ,
B(5) = 61.883345 .