anchor.diff:                ELSE IF (KB.EQ.2) THEN
                    DF(2)=F2
                    FS=-((D(2)-CV)*F2+G(2)*F1)/A(2)
  -              ELSE 
  +              ELSE
                    DF(2)=F2
                    DO 20 J=3,KB+1
                       F=-((D(J-1)-CV)*F2+G(J-1)*F1)/A(J-1)
  
positive.code:         END
  
  
          
  C       **********************************
  C       SciPy: Changed P from a character array to an integer array.
          SUBROUTINE JDZO(NT,N,M,P,ZO)
  
negative.code: C
  C       =======================================================
  C       Purpose: Compute the associated Legendre function
  C                Pmv(x) with an integer order and an arbitrary
  C                nonnegative degree v
  C       Input :  x   --- Argument of Pm(x)  ( -1 ≤ x ≤ 1 )
  C                m   --- Order of Pmv(x)
  