anchor.diff:  C                BY(n-NMIN) --- Yn(x)   ; if indexing starts at 0
   C                NM --- Highest order computed
   C       Routines called:
  -C                MSTA1 and MSTA2 to calculate the starting 
  +C                MSTA1 and MSTA2 to calculate the starting
   C                point for backward recurrence
   C       =====================================================
   C
  
positive.code:         END
  
  
          
  C       **********************************
  
          SUBROUTINE CFS(Z,ZF,ZD)
  
negative.code:         END
  
  
  
  C       **********************************
  
          SUBROUTINE CJYVB(V,Z,VM,CBJ,CDJ,CBY,CDY)
  C
  C       ===========================================================
  C       Purpose: Compute Bessel functions Jv(z), Yv(z) and their
  C                derivatives for a complex argument
  C       Input :  z --- Complex argument
  C                v --- Order of Jv(z) and Yv(z)
  