lit
    JINRLIB

DPOLY2 -

C206

: T.Pomentale :

DPOLY2 c op ooa c N c oc oa A1, A2, ..., An+1:

C206

:

: - SUBROUTINE
: - DPOLY2
: - DCIRCL,DDEFLT,DVP

:

CALL DPOLY2(A,N,ROOT,H,MAXFUN,MODE), :

A - (COMPLEX*16) acc papoc >=N+1, o A(I), I=1,2,...,N+1;
N - (INTEGER) c ooa P(Z);
ROOT - (COMPLEX*16) acc papoc >=N, MODE, 0, . , MODE = 0. ROOT N .
MAXFUN - (INTEGER) acao co pa;
MODE - (INTEGER) c. ROOT;
H - (REAL*8) acc papoc >=N, copa pac H(I) opoc c po ROOT(I), aa c c a op ooa.

:

N<=20. N>20 o oo pa a c oo op. c co pa pa pao a MAXFUN, aac , .

:

  1. T.Pomntal. Homotop itrativ mthods for polynomial quations, J.Inst. Maths. Applics, (1974) 13, 201-213.
:
       IMPLICIT REAL*8 (A-H,O-Z)
       COMPLEX*16 A(5),ROOT(4)
       DIMENSION H(4)
       N=4
       A(1)=( 1.D0, 0.D0)
       A(2)=(-8.D0, 0.D0)
       A(3)=12.D0*(0.117038D-2,-0.1709D0)
       A(4)=( 8.D0, 0.D0)
       A(5)=(-1.D0, 0.D0)
       CALL DPOLY2(A,N,ROOT,H,2500,0)
       WRITE(*,*) (ROOT(I),I=1,N)
       . . .

:

       ROOT(1) = ( 0.12668606459D+00  ,  0.43106652050D-02)
       ROOT(2) = ( 0.98658544797D+00  , -0.17027327325D+00)
       ROOT(3) = (-0.99401348618D+00  , -0.10257573240D+00)
       ROOT(4) = ( 0.78807419736D+01  ,  0.26853834044D+00)


home up e-mail