anchor.diff:  C
           IMPLICIT DOUBLE PRECISION (A-H,O-Z)
           DIMENSION BN(0:N)
  -        TPI=6.283185307179586D0 
  +        TPI=6.283185307179586D0
           BN(0)=1.0D0
           BN(1)=-0.5D0
           BN(2)=1.0D0/6.0D0
  
positive.code: C
  C       ========================================================
  C       Purpose: Compute the expansion coefficients for the
  C                asymptotic expansion of Bessel functions 
  C                with large orders
  C       Input :  Km   --- Maximum k
  C       Output:  A(L) --- Cj(k) where j and k are related to L 
  C                         by L=j+1+[k*(k+1)]/2; j,k=0,1,...,Km
  C       ========================================================
  C
  
negative.code: 35         DL(N)=-0.5D0*X/(N+1.0D0)*UK
             RETURN
          ENDIF
          IF (N.EQ.0) NM=1
          M=MSTA1(X,200)
          IF (M.LT.NM) THEN
             NM=M
  