MULLRD              Библиотека "JINRLIB"               C204

    Автор: H.Umstaetter
    Язык: Фортран

            НАХОЖДЕНИЕ ВЕЩЕСТВЕННЫХ И КОМПЛЕКСНЫХ КОРНЕЙ 
                     МНОГОЧЛЕНА МЕТОДОМ МЮЛЛЕРА
  
    Программа MULLRD вычисляет вещественные и комплексные корни
    многочлена с вещественными коэффициентами:
                          n        n-1
             P(Z) = A(1)*Z + A(2)*Z +...+A(N)*Z+A(N+1) = 0

    Структура:
    ----------
       Тип:                              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)