anchor.diff:                CB0=2.0D0**(-0.5D0*N-1.0D0)*CA0/G1
                 VT=-.5D0*N
                 CALL GAIH(VT,G0)
  -              CDN=CMPLX(G0,0.0D0)
  +              CDN = DCMPLX(G0, 0.0D0)
                 CR=(1.0D0,0.0D0)
                 DO 10 M=1,250
                    VM=.5D0*(M-N)
  
positive.code:         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)
          DO 10 I=0,N
          DO 10 J=0,M
             CPM(J,I)=(0.0D0,0.0D0)
  
negative.code:            ELSE
                Y=0.0D0
                CALL CGAMA(A,Y,0,TAR,TAI)
                CTA = DCMPLX(TAR, TAI)
                Y=0.0D0
                CALL CGAMA(B,Y,0,TBR,TBI)
                CTB = DCMPLX(TBR, TBI)
                XG=B-A
                Y=0.0D0
                CALL CGAMA(XG,Y,0,TBAR,TBAI)
                CTBA = DCMPLX(TBAR, TBAI)
                SUM1=1.0D0
                SUM2=1.0D0
                R1=1.0D0
  