Текст подпрограммы и версий
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 

Подпрограмма:  ZP18R

Назначение

Вычисление числа вещественных корней полинома с вещественными коэффициентами вне и внутри заданного отрезка.

Математическое описание

Подпрограмма 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