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:            N0=N1
             F0=F1
             N1=NN
  10         F1=F
  20      MSTA2=NN+10
          RETURN
          END FUNCTION
  
negative.code:         IF (X.LT.1.0D-60) THEN
             DO 10 K=0,N
                SY(K)=-1.0D+300
                DY(K)=1.0D+300
  10         END DO
             RETURN
          ENDIF
          SY(0)=-DCOS(X)/X
  