anchor.diff:        COMMON /CBESS/ J,Y,JR,JI,DJ,DY,DJR,DJI
         NG=NGAUSS*2
         IF (NP.GT.0) CALL RSP2(X,NG,A,EPS,NP,R,DR)
  -      IF (NP.EQ.-1) CALL RSP1(X,NG,NGAUSS,A,EPS,NP,R,DR)
  +      IF (NP.EQ.-1) CALL RSP1(X,NG,NGAUSS,A,EPS,R,DR)
         IF (NP.EQ.-2) CALL RSP3(X,NG,NGAUSS,A,EPS,R,DR)
         IF (NP.EQ.-3) CALL RSP4(X,NG,A,R,DR)
         PI=P*2D0/LAM
  
positive.code: 
  C**********************************************************************
  
        SUBROUTINE TMATR0 (NGAUSS,X,W,AN,ANN,S,SS,PPI,PIR,PII,R,DR,DDR,
       *                  DRR,DRI,NMAX,NCHECK)
        INCLUDE 'ampld.par.f'
        IMPLICIT REAL*8 (A-H,O-Z)
        REAL*8  X(NPNG2),W(NPNG2),AN(NPN1),S(NPNG2),SS(NPNG2),
       *        R(NPNG2),DR(NPNG2),SIG(NPN2),
       *        J(NPNG2,NPN1),Y(NPNG2,NPN1),
       *        JR(NPNG2,NPN1),JI(NPNG2,NPN1),DJ(NPNG2,NPN1),
       *        DY(NPNG2,NPN1),DJR(NPNG2,NPN1),
       *        DJI(NPNG2,NPN1),DDR(NPNG2),DRR(NPNG2),
       *        D1(NPNG2,NPN1),D2(NPNG2,NPN1),
       *        DRI(NPNG2),DS(NPNG2),DSS(NPNG2),RR(NPNG2),
       *        DV1(NPN1),DV2(NPN1)
  
        REAL*8  R11(NPN1,NPN1),R12(NPN1,NPN1),
       *        R21(NPN1,NPN1),R22(NPN1,NPN1),
  
negative.code:             D1(I2,N)=DD1*SI
              D2(I2,N)=-DD2*SI
     25 CONTINUE
        DO 40 I=1,NGSS
             RR(I)=W(I)*R(I)
     40 CONTINUE
  
  