anchor.diff:          CG1=CBY1
           DO 90 K=2,NM
              CYK=2.0D0*(K-1.0D0)/Z*CG1-CG0
  -           IF (CDABS(CYK).GT.1.0D+290) GO TO 90            
  +           IF (CDABS(CYK).GT.1.0D+290) GO TO 90
              YAK=CDABS(CYK)
              YA1=CDABS(CG0)
              IF (YAK.LT.YA0.AND.YAK.LT.YA1) LB=K
  
positive.code:         END
  
  
          
  C       **********************************
  
          SUBROUTINE KMN(M,N,C,CV,KD,DF,DN,CK1,CK2)
  
negative.code:         END
  
  
  
  C       **********************************
  
          SUBROUTINE AIRYB(X,AI,BI,AD,BD)
  