anchor.diff:          END
   
   
  -        
  +
   C       **********************************
   
           SUBROUTINE STVH0(X,SH0)
  
positive.code:         END
  
  
          
  C       **********************************
  
          SUBROUTINE CLPN(N,X,Y,CPN,CPD)
  
negative.code: 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                point for backward recurrence
  C       =====================================================
  C
  