|
Текст подпрограммы и версий ( Фортран ) zp33r.zip zp33d.zip |
Тексты тестовых примеров ( Фортран ) tzp33r.zip tzp33d.zip |
|
Текст подпрограммы и версий ( Си ) zp33r_c.zip zp33d_c.zip |
Тексты тестовых примеров ( Си ) tzp33r_c.zip tzp33d_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) zp33r_p.zip zp33e_p.zip |
Тексты тестовых примеров ( Паскаль ) tzp33r_p.zip tzp33e_p.zip |
Вычисление М первых коэффициентов полинома, который является обратным заданному полиному с вещественными коэффициентами.
ZР33R вычисляет М первых коэффициентов полинома
B(X) = X(- K)*( B(1) + B(2)*X + ... + B(M)*X(M - 1) + ... )
обратного заданному полиному степени NР
A(X) = A(1) + A(2)*X + A(3)*X2 + ... + A(NP + 1)*XNP ,
у которого коэффициенты А ( I ), I = 1, 2, 3, ..., NР + 1 являются вещественными. Здесь К - число первых коэффициентов полинома А (Х), равных 0. Иными словами, ищется полином В (Х) = 1 / А (Х).
SUBROUTINE ZP33R (NP, A, M, B, K)
Параметры
| NP - | заданная степень полинома А (Х) (тип: целый); |
| A - | вещественный вектор длины NР + 1, содержащий коэффициенты заданного полинома А (Х) в порядке возрастания степеней; |
| M - | заданное число искомых коэффициентов полинома В (Х) (тип: целый); |
| B - | вещественный вектор длины М, содержащий вычисленные коэффициенты полинома В (Х) в порядке возрастания степеней; |
| K - | целая переменная, значение которой в результате работы подпрограммы полагается равным числу первых коэффициентов полинома А (Х), равных 0. |
Версии
| ZP33D - | вычисление с повышенной точностью М первых коэффициентов полинома, который является обратным заданному полиному с вещественными коэффициентами. При этом векторы А и В должны иметь тип DОUВLЕ РRЕСISIОN. |
Вызываемые подпрограммы: нет
Замечания по использованию: нет
DIMENSION A(5), B(6)
DATA A /24., - 50., 35., - 10., 1./
NP = 4
M = 6
CALL ZP33R (NP, A, M, B, K)
Результаты:
В(1) = 0.041667 , В(2) = 0.0868056 , В(3) = 0.120081 ,
B(4) = 0.140938 , B(5) = 0.152937 , B(6) = 0.159499