anchor.diff:          R0=REG
           DO 10 J=1,2*M+IP
   10         R0=R0*J
  -        R=R0    
  +        R=R0
           SUC=R*DF(1)
           SW=0.0D0
           DO 15 K=2,NM
  
positive.code:         END
  
  
          
  C       **********************************
  
          SUBROUTINE SCKB(M,N,C,DF,CK)
  
negative.code:         END
  
  
  
  C       **********************************
  
          SUBROUTINE FCSZO(KF,NT,ZO)
  