anchor.diff:  
     end subroutine mevalhc
   
  -
   end module trdf
  
positive.code: 
      IF ( FLAG .NE. 0 ) RETURN
  
      VQUAD_A=  F + Q(1) ! MODEL IN XNOVO 
      VQUAD=  -VQUAD + VQUAD_A  
      IF (VQUAD .GE. 0.D0 .or. cnorm .gt. xeps)  THEN
         DO I=1, N
            X(I) = XANTIGO(I) 
  
negative.code:   ! ******************************************************************
    ! ******************************************************************
  
    SUBROUTINE SIGMA(H,N,NPT,Y,X,VETOR1,SIGM,ALFA,BETA,TAU,IT,DELTA)
  
      implicit none
  