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:            IF (KF.EQ.2) PSQ=2.0D0*NR**(0.5)
             PX=PSQ-DLOG(PI*PSQ)/(PI*PI*PSQ**3.0)
             PY=DLOG(PI*PSQ)/(PI*PSQ)
             Z=CMPLX(PX,PY)
             IF (KF.EQ.2) THEN
                IF (NR.EQ.2) Z=(2.8334,0.2443)
                IF (NR.EQ.3) Z=(3.4674,0.2185)
  
negative.code:               CB0=2.0D0**(-0.5D0*N-1.0D0)*CA0/G1
                VT=-.5D0*N
                CALL GAIH(VT,G0)
                CDN = DCMPLX(G0, 0.0D0)
                CR=(1.0D0,0.0D0)
                DO 10 M=1,250
                   VM=.5D0*(M-N)
  