Текст подпрограммы и версий zp18r.zip zp18d.zip |
Тексты тестовых примеров tzp18r.zip tzp18d.zip |
Текст подпрограммы и версий ( Си ) zp18r_c.zip zp18d_c.zip |
Тексты тестовых примеров ( Си ) tzp18r_c.zip tzp18d_c.zip |
Текст подпрограммы и версий ( Паскаль ) zp18r_p.zip zp18e_p.zip |
Тексты тестовых примеров ( Паскаль ) tzp18r_p.zip tzp18e_p.zip |
Вычисление числа вещественных корней полинома с вещественными коэффициентами вне и внутри заданного отрезка.
Подпрограмма ZР18R вычисляет число вещественных корней полинома
A(1) + A(2)*X + A(3)*X2 + ... + A(NP + 1)*XNP
с вещественными коэффициентами вне и внутри заданного отрезка [ХN, ХК] методом Штурма.
Дж.Х.Уилкинсон, Алгебраическая проблема собственных значений, "Наука", М., 1970.
SUBROUTINE ZP18R (NP, A, XN, XK, NB, NR, NA, RAB, NP1, IERR)
Параметры
NP - | заданная степень полинома (тип: целый); |
A - | вещественный вектор длины NР + 1, содержащий коэффициенты исходного полинома в порядке возрастания степеней; |
XN - | заданная нижняя граница отрезка (тип: вещественный); |
XK - | заданная верхняя граница отрезка (тип: вещественный); |
NB - NR NA | целые перемеменные, значения которых в результате работы подпрограммы полагаются равными числу вещественных корней исходного полинома, меньших ХN, лежащих внутри [ХN, ХК] и больших ХК соответственно; |
RAB - | вещественный двумерный массив размера NР1 на NP1, используемый в подпрограмме как рабочий; |
NP1 - | задается равным NР + 1 (тип: целый); |
IERR - | целая переменная, значение которой полагается равным 65, если в пределах ошибок округления заданный полином имеет кратные корни. |
Версии
ZP18D - | вычисление числа вещественных корней полинома с вещественными коэффициентами вне и внутри заданного отрезка с повышенной точностью. При этом вектор А, двумерный массив RАВ, а также ХN и ХК должны иметь тип DОUВLЕ РRЕСISIОN. |
Вызываемые подпрограммы
UTZP10 - | подпрограмма выдачи диагностических сообщений при работе подпрограммы ZР18R. |
UTZP11 - | подпрограмма выдачи диагностических сообщений при работе подпрограммы ZР18D. |
Замечания по использованию: нет
DIMENSION A(5), RAB(5, 5) DATA A(1), A(2), A(3), A(4), A(5) /24., - 50., 35., - 10., 1./ NP = 4 NP1 = 5 XN = 1.5 XK = 3.5 CALL ZP18R (NP, A, XN, XK, NB, NR, NA, RAB, NP1, IERR) Результаты: NВ = 1 , NR = 2 , NА = 1 , IЕRR = 0