anchor.diff:          RJAC(2,3) = -Vv*FUGVy(icomp)
   C
           dold=delx
  -c        CALL DLSARG (N, RJAC, LDA, -F, IPATH, delX)
  +c       CALL DLSARG (N, RJAC, LDA, -F, IPATH, delX)
           db = -F
           AJ=RJAC
           call dgesv( N, 1, AJ, lda, ipiv, db, ldb, info )
  
negative.code:     INCX=1
      LDA=3
      ldb=3
  c   IPATH=1
      TCmod=TC
      if(TCGC(1).ne.0.0D0)TCmod=TCGC
      inf=2
  