anchor.diff:          SUBROUTINE LGAMA(KF,X,GL)
   C
   C       ==================================================
  -C       Purpose: Compute gamma function (x) or ln[(x)]
  -C       Input:   x  --- Argument of (x) ( x > 0 )
  +C       Purpose: Compute gamma function Г(x) or ln[Г(x)]
  +C       Input:   x  --- Argument of Г(x) ( x > 0 )
   C                KF --- Function code
  -C                       KF=1 for (x); KF=0 for ln[(x)]
  -C       Output:  GL --- (x) or ln[(x)]
  +C                       KF=1 for Г(x); KF=0 for ln[Г(x)]
  +C       Output:  GL --- Г(x) or ln[Г(x)]
   C       ==================================================
   C
           IMPLICIT DOUBLE PRECISION (A-H,O-Z)
  
positive.code: C
  C       ========================================================
  C       Purpose: Compute the characteristic value of Mathieu
  C                functions  for q  3m
  C       Input :  m  --- Order of Mathieu functions
  C                q  --- Parameter of Mathieu functions
  C       Output:  A0 --- Initial characteristic value
  
negative.code: C             (2) CV0 for finding initial characteristic
  C                 values using polynomial approximation
  C             (3) CVQM for computing initial characteristic
  C                 values for q ≤ 3*m
  C             (3) CVQL for computing initial characteristic
  C                 values for q ≥ m*m
  C       ======================================================
  C
          IMPLICIT DOUBLE PRECISION (A-H,O-Z)
  