anchor.diff:        SC= COS(DEL)
         SI=0.
         CO =1.0
  -      do 50 I=1,NT
  -      K6=NTOT2-2*I+5
  -      CIRE= A(2*I-1) + A(K6)
  -      CIIM=A(2*I)-A(K6+1)
  -      CNIRE=(-SI*(A(2*I)+A(K6+1))+CO*(A(2*I-1)-A(K6)))
  -      if (SI)62,61,62
  -   62 CNIIM=(A(2*I-1)-A(K6)-CO*CNIRE)/SI
  -      go to 63
  -   61 CNIIM=0.
  -   63 A(2*I-1)=CIRE
  -      A(2*I)=CIIM
  -      A(K6)=CNIRE
  -      A(K6+1)=CNIIM
  -      SIS=SI
  -      SI=SI*SC+CO*SS
  -   50 CO=CO*SC-SIS*SS
  +      do I=1,NT
  +        K6=NTOT2-2*I+5
  +        CIRE= A(2*I-1) + A(K6)
  +        CIIM=A(2*I)-A(K6+1)
  +        CNIRE=(-SI*(A(2*I)+A(K6+1))+CO*(A(2*I-1)-A(K6)))
  +        if (SI)62,61,62
  +  62    CNIIM=(A(2*I-1)-A(K6)-CO*CNIRE)/SI
  +        go to 63
  +  61    CNIIM=0.
  +  63    A(2*I-1)=CIRE
  +        A(2*I)=CIIM
  +        A(K6)=CNIRE
  +        A(K6+1)=CNIIM
  +        SIS=SI
  +        SI=SI*SC+CO*SS
  +        CO=CO*SC-SIS*SS
  +      end do
         KO=NTOT+1
  -      do 40 I=1,KO,2
  -      K1=NTOT2-I+4
  -      AP1RE=A(I)-A(K1+1)
  -      AP2RE=-(A(I+1)+A(K1))
  -      AP1IM=A(I)+A(K1+1)
  -      AP2IM=A(I+1)-A(K1)
  -      A(I)=AP1RE
  -      A(I+1)=AP2RE
  -      A(K1)=AP1IM
  -   40 A(K1+1)=AP2IM
  +      do I=1,KO,2
  +        K1=NTOT2-I+4
  +        AP1RE=A(I)-A(K1+1)
  +        AP2RE=-(A(I+1)+A(K1))
  +        AP1IM=A(I)+A(K1+1)
  +        AP2IM=A(I+1)-A(K1)
  +        A(I)=AP1RE
  +        A(I+1)=AP2RE
  +        A(K1)=AP1IM
  +        A(K1+1)=AP2IM
  +      end do
         NTOP=NTOT2+2
         NT00=NTOT+1
         A(1)=A(NTOT2+3)
         A(2)=A(NTOT2+4)
  -   21 do 52 I=NT00,NTOP,2
  -      A(I)=A(I+2)
  -   52 A(I+1)=A(I+3)
  +   21 do I=NT00,NTOP,2
  +        A(I)=A(I+2)
  +        A(I+1)=A(I+3)
  +      end do
         call  FFT(A,L,INV,S,IFSET,IFERR)
  -      do 20 I=1,NTOT2
  -   20 A(I)=A(I)*FN
  -      do 10 I=2,NTOT2,2
  -   10 A(I)=-A(I)
  +      do I=1,NTOT2
  +        A(I)=A(I)*FN
  +      end do
  +      do I=2,NTOT2,2
  +        A(I)=-A(I)
  +      end do
         return
         end
   C*************************************
  
positive.code: C
        dimension X(1)
        XMAX = 0.
        do 1 I = 1,MX
        XA = ABS(X(I))
        if (XMAX .gt. XA) go to 1
        NXMAX = I
        XMAX = XA
      1 continue
        return
        end
  
negative.code:       A(NN+1)=A(NN-1)
        FN=NTOT
        NTOT3=NTOT2+4
        do I=3,NTOT2,2
          A(I)=0.5* A(I)
          A(I+1)= .5*A(I+1)
        end do
        do I=1,NTOT,2
          K8=NTOT2+2-I
          A(K8)= A(K8-2)
          A(K8+1)=A(K8-1)
        end do
        NTO=NTOT/ 2
        NT=NTO+1
        DEL=3.141592654/FN
  