|
Текст подпрограммы и версий ( Фортран ) zp34r.zip zp34d.zip |
Тексты тестовых примеров ( Фортран ) tzp34r.zip tzp34d.zip |
|
Текст подпрограммы и версий ( Си ) zp34r_c.zip zp34d_c.zip |
Тексты тестовых примеров ( Си ) tzp34r_c.zip tzp34d_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) zp34r_p.zip zp34e_p.zip |
Тексты тестовых примеров ( Паскаль ) tzp34r_p.zip tzp34e_p.zip |
Вычисление натурального логарифма от полинома с вещественными коэффициентами.
ZР34R вычисляет М первых коэффициентов полинома
B(X) = B(1) + B(2)*X + B(3)*X2 + ... + B(M)*X(M - 1) + ... ,
такого, что В (Х) = ln А (Х), где А (Х) заданный полином степени NР
A(X) = A(1) + A(2)*X + A(3)*X2 + ... + A(NP + 1)*XNP ,
у которого коэффициенты А ( I ), I = 1, 2, 3, ..., NР + 1 являются вещественными и А (1) > 0.
SUBROUTINE ZP34R (NP, A, M, B, IERR)
Параметры
| NP - | заданная степень полинома А (Х) (тип: целый); |
| A - | вещественный вектор длины NР + 1, содержащий коэффициенты заданного полинома А (Х) в порядке возрастания степеней; |
| M - | заданное число искомых коэффициентов полинома В (Х) (тип: целый); |
| B - | вещественный вектор длины М, содержащий вычисленные коэффициенты полинома В (Х) в порядке возрастания степеней; |
| IERR - | целая переменная, значение которой в результате работы подпрограммы полагается равной 65, если А (1) Ј 0. |
Версии
| ZP34D - | вычисление натурального логарифма от полинома с вещественными коэффициентами с повышенной точностью. При этом векторы А и В должны иметь тип DОUВLЕ РRЕСISIОN. |
Вызываемые подпрограммы
| UTZP10 - | подпрограмма выдачи диагностических сообщений при работе подпрограммы ZР34R. |
| UTZP11 - | подпрограмма выдачи диагностических сообщений при работе подпрограммы ZР34D. |
Замечания по использованию: нет
DIMENSION A(5), B(6)
DATA A(1), A(2), A(3), A(4), A(5) /24., - 50., 35., - 10., 1./
NP = 4
M = 6
CALL ZP34R (NP, A, M, B, IERR)
Результаты:
В(1) = 3.178053 , В(2) = - 2.083333 , В(3) = - 0.711556 ,
B(4) = - 0.392554 , B(5) = - 0.269688 , B(6) = - 0.207268 ,
IERR = 0