anchor.diff:  C       Purpose: Compute lambda function with arbitrary order v,
   C                and their derivative
   C       Input :  x --- Argument of lambda function
  -C                v --- Order of lambda function 
  +C                v --- Order of lambda function
   C       Output:  VL(n) --- Lambda function of order n+v0
  -C                DL(n) --- Derivative of lambda function 
  +C                DL(n) --- Derivative of lambda function
   C                VM --- Highest order computed
   C       Routines called:
  -C            (1) MSTA1 and MSTA2 for computing the starting 
  +C            (1) MSTA1 and MSTA2 for computing the starting
   C                point for backward recurrence
   C            (2) GAM0 for computing gamma function (|x| ≤ 1)
   C       =========================================================
  
positive.code:         END
  
  
          
  C       **********************************
  
          SUBROUTINE LAMN(N,X,NM,BL,DL)
  
negative.code:         END
  
  
  
  C       **********************************
  
          SUBROUTINE LPMNS(M,N,X,PM,PD)
  