anchor.diff:  
     close(99)
   
  +  ! General options
  +
  +  verbose = .false.
  +
  +  ftype = 2
  +
     ! Call the solver - First run
   
  -  verbose = .false.
  -
     epsfeas = 1.0D-04
   
     epsopt  = 1.0D-01
  
positive.code: 
    ! Call the solver - Second run
  
    verbose = .false.
  
    epsfeas = 1.0D-04
  
    epsopt  = 1.0D-04
  
    call set_epsfeas(epsfeas)
  
    PEN = 1.0D+10
  
    ftype = 2
  
    call fird(n,x,l,u,me,mi,evalf,evalc,evaljac,verbose,ftype, &
         epsfeas,epsopt,f,feas,fcnt,flag)
  
negative.code:   call set_epsfeas(epsfeas)
  
    PEN = 1.0D+03
  
    call fird(n,x,l,u,me,mi,evalf,evalc,evaljac,verbose,ftype, &
         epsfeas,epsopt,f,feas,fcnt,flag)
  