anchor.diff:          IMPLICIT DOUBLE PRECISION (D,X,Y)
           IMPLICIT COMPLEX*16 (C,Z)
           DIMENSION CPM(0:MM,0:N),CPD(0:MM,0:N)
  -        Z=CMPLX(X,Y)
  +        Z = DCMPLX(X, Y)
           DO 10 I=0,N
           DO 10 J=0,M
              CPM(J,I)=(0.0D0,0.0D0)
  
positive.code:         IMPLICIT DOUBLE PRECISION (X,Y)
          IMPLICIT COMPLEX*16 (C,Z)
          DIMENSION CQM(0:MM,0:N),CQD(0:MM,0:N)
          Z=CMPLX(X,Y)
          IF (DABS(X).EQ.1.0D0.AND.Y.EQ.0.0D0) THEN
             DO 10 I=0,M
             DO 10 J=0,N
  
negative.code:         IMPLICIT DOUBLE PRECISION (X,Y)
          IMPLICIT COMPLEX *16 (C,Z)
          DIMENSION CPN(0:N),CPD(0:N)
          Z = DCMPLX(X, Y)
          CPN(0)=(1.0D0,0.0D0)
          CPN(1)=Z
          CPD(0)=(0.0D0,0.0D0)
  