anchor.diff:          IF (DABS(X).LT.1.0D-100) THEN
              DO 10 K=0,N
                 SJ(K)=0.0D0
  -10            DJ(K)=0.0D0
  +              DJ(K)=0.0D0
  +10         END DO
              SJ(0)=1.0D0
              IF (N.GT.0) THEN
                 DJ(1)=.3333333333333333D0
  
positive.code:               F=(2.0D0*K+3.0D0)*F1/X-F0
                IF (K.LE.NM) SJ(K)=F
                F0=F1
  15            F1=F
             CS=0.0D0
             IF (DABS(SA).GT.DABS(SB)) CS=SA/F
             IF (DABS(SA).LE.DABS(SB)) CS=SB/F0
             DO 20 K=0,NM
  20            SJ(K)=CS*SJ(K)
          ENDIF
          DO 25 K=1,NM
  25         DJ(K)=SJ(K-1)-(K+1.0D0)*SJ(K)/X
          RETURN
          END SUBROUTINE
  
  
negative.code:            N0=N1
             F0=F1
             N1=NN
             F1=F
   10     END DO
   20     MSTA1=NN
          RETURN
          END FUNCTION
  