anchor.diff:  C                                                                     *
   C**********************************************************************
   
  -      SUBROUTINE TT(NMAX,NCHECK)
  +      SUBROUTINE TT(NMAX)
         INCLUDE 'ampld.par.f'
         IMPLICIT REAL*8 (A-H,O-Z)
  -      REAL*8 F(NPN2,NPN2),B(NPN2),WORK(NPN2),
  -     *       QR(NPN2,NPN2),QI(NPN2,NPN2),
  -     *       RGQR(NPN2,NPN2),RGQI(NPN2,NPN2),
  -     *       A(NPN2,NPN2),C(NPN2,NPN2),D(NPN2,NPN2),E(NPN2,NPN2)
  +      REAL*8 QR(NPN2,NPN2),QI(NPN2,NPN2),
  +     *       RGQR(NPN2,NPN2),RGQI(NPN2,NPN2)
         REAL*8 TR1(NPN2,NPN2),TI1(NPN2,NPN2)
         COMPLEX*16 ZQ(NPN2,NPN2),ZW(NPN2)
  -      INTEGER IPIV(NPN2),IPVT(NPN2)
  +      INTEGER IPIV(NPN2)
         COMMON /CT/ TR1,TI1
         COMMON /CTT/ QR,QI,RGQR,RGQI
         NDIM=NPN2
  
positive.code:            SI=-SI
             SIG(N)=SI
      5 CONTINUE
     20 DO 25 I=1,NGAUSS
           I1=NGAUSS+I
           I2=NGAUSS-I+1
           CALL VIG (X(I1),NMAX,M,DV1,DV2)
  
negative.code: 
  C**********************************************************************
  
        SUBROUTINE CONST (NGAUSS,NMAX,X,W,AN,ANN,S,SS,NP,EPS)
        IMPLICIT REAL*8 (A-H,O-Z)
        INCLUDE 'ampld.par.f'
        REAL*8 X(NPNG2),W(NPNG2),X1(NPNG1),W1(NPNG1),
  