anchor.diff:  
       IF ( FLAG .NE. 0 ) RETURN
   
  -    VQUAD_A=  F + Q(1) ! MODEL IN XNOVO 
  -    VQUAD=  -VQUAD + VQUAD_A  
  +    VQUAD_A = F + Q(1) ! Model in the new point
  +    VQUAD   = VQUAD_A - VQUAD ! Denominator for the sufficient decrease test
  +
       IF (VQUAD .GE. 0.D0 .or. cnorm .gt. xeps)  THEN
          DO I=1, N
             X(I) = XANTIGO(I) 
  
positive.code: 
      ! LOCAL ARRAYS
      REAL(8) :: FF(NPT),D(INN),Y(NPT,N),Q(1+N+N*(N+1)/2), &
           H(NPT+N+1,NPT+N+1),XNOVO(INN),SL(INN),SU(INN), VETOR1(NPT+N+1) 
  
      ! LOCAL SCALARS
      integer :: i,it,j,k,kn,flag
      real(8) :: alfa,beta,c,cnorm,delta,distsq,dsq,fopt,gama, &
           mindelta,rho,rhobeg,rhoend,sigm,sum,tau,tempofinal,tempoinicial
  
      IF ( OUTPUT ) WRITE(*,3000)
  
  
negative.code:     evaljac => evaljac_
      evalhc  => evalhc_
  
      F      = 1.0D300
      IC     = 0
      MAXIC  = maxfcnt ! MAXIMUM NUMBER OF FUNCTION EVALUATIONS
      RHOBEG = RBEG
      RHO    = RHOBEG
      DELTA  = RHO
      RHOEND = REND
      GAMA   = 0.1D0
  
      !     ---------------------------
      !     Feasibility phase - Phase 0
  