Программа MULLRD вычисляет вещественные и комплексные корни многочлена с вещественными коэффициентами:

Структура:

Тип:	-	SUBROUTINE
Имена входа для пользователя:	-	MULLRD

Обращение:

CALL MULLRD(A,N,C), где:

А	-	(REAL*8) массив вещественных коэффициентов многочлена длиной N+1, расположенных в порядке убывания степеней;
N	-	(INTEGER) степень многочлена n;
C	-	(COMPLEX*16) массив вычисляемых комплексных корней многочлена длиною N. На выходе массив C(N) содержит N корней.

Метод:

Нахождение комплексных корней многочлена N-й степени c вещественными коэффициентами методом Мюллера.

Литература:

1. D.E.Muller, MTAC, 1965, pp.208-215.

Пример:

       . . .
       IMPLICIT REAL*8 (A-H,O-Z)
       COMPLEX*16 C(4)
       DIMENSION A(5)
       DATA A/1.0D0,-10.0D0,35.0D0,-50.0D0,24.0D0/   ! (x-1)*(x-2)*(x-3)*(x-4)
       . . .
       CALL MULLRD(A,4,C)
       DO I=1,4
         WRITE(*,*) C(I)
       ENDDO
       . . .

Результат:

       (2.999999999999988,0.000000000000000E+000)
       (4.000000000000010,0.000000000000000E+000)
       (1.000000000000001,0.000000000000000E+000)
       (2.000000000000001,0.000000000000000E+000)

---