anchor.diff:          IF (X.LT.1.0D-60) THEN
              DO 10 K=0,N
                 SY(K)=-1.0D+300
  -10            DY(K)=1.0D+300
  +              DY(K)=1.0D+300
  +10         END DO
              RETURN
           ENDIF
           SY(0)=-DCOS(X)/X
  
positive.code:         IF (DABS(X).LT.1.0D-100) THEN
             DO 10 K=0,N
                SJ(K)=0.0D0
  10            DJ(K)=0.0D0
             SJ(0)=1.0D0
             IF (N.GT.0) THEN
                DJ(1)=.3333333333333333D0
  
negative.code:            N0=N1
             F0=F1
             N1=NN
             F1=F
   10     END DO
   20     MSTA1=NN
          RETURN
          END FUNCTION
  