anchor.diff:    ! Maximum number of iterations
     integer, parameter :: MAXITER = 10000
     ! Maximum number of printing elements
  -  integer, parameter :: MAXNEL  = 3
  +  integer, parameter :: MAXNEL  = 20
     ! Minimum trust region radius
     real(8), parameter :: DELMIN = 1.0D-30
     ! Initial value of RHO
  
positive.code:   implicit none
  
    ! PARAMETERS
    real(8), parameter :: BETA = 1.0D+04
    real(8), parameter :: ALPHA = 1.0D-4
    real(8), parameter :: MU = 1.0D-1
    real(8), parameter :: ETA = 2.5D-1
    ! Restoration reduction factor
  
negative.code:        hxnorm = hynorm
         
         iter = iter + 1
  
         ! Verify convergence conditions
  
  