anchor.diff:          ENDIF
           IF (KF.EQ.1) X=-0.503
           IF (KF.EQ.2) X=0.577
  -        ZERO=CMPLX(X,Y)
  +        ZERO = DCMPLX(X, Y)
           Z=ZERO
           W=0.0D0
           DO 35 NR=1,NT
  
positive.code:               ELSE
                   CALL GAIH(VA0,GA0)
                   PD=DSQRT(PI)/(2.0D0**(-.5D0*N)*GA0)
                   CDN=CMPLX(PD,0.0D0)
                ENDIF
             ELSE
                XN=-N
  
negative.code:         IMPLICIT DOUBLE PRECISION (D,X,Y)
          IMPLICIT COMPLEX*16 (C,Z)
          DIMENSION CPM(0:MM,0:N),CPD(0:MM,0:N)
          Z = DCMPLX(X, Y)
          DO 10 I=0,N
          DO 10 J=0,M
             CPM(J,I)=(0.0D0,0.0D0)
  