anchor.diff:    integer, parameter :: MMAX = 1000
     integer, parameter :: JCNNZMAX = NMAX * MMAX
     integer, parameter :: HCNNZMAX = NMAX ** 2 * (1 + MMAX)
  -  integer, parameter :: MAXXEL = 3
  +  integer, parameter :: MAXXEL = 30
     integer, parameter :: INN = 1000
   
     ! COMMON SCALARS
  
positive.code:     DELTA   = MAX(DELTA, RHO) ! Correcting delta if necessary
      GAMMA   = 0.1D0
  
      IF (OUTPUT) WRITE(*,1001)
  
      DO I=1,N
         Z(I) = X(I)
         XNOVO(I) = X(I)
  
negative.code:     CALL  PRIMEIROMODELO1 (N,X,Q_,H_,NPT,RHO,Y_,FF_,FLAG) 
  
      IF ( OUTPUT ) WRITE(*,1002) RHO,DELTA,FF_(1),IC,MIN(N,MAXXEL), &
           (X(I), I=1,MIN(N,MAXXEL))
      IF ( FLAG .NE. 0 ) GOTO 31
  
      FZ = FF_(1)         
  
  11  CALL SUBPROBLEMA(N,NPT,Q_,DELTA,D,X,XL,XU,DSQ,M,EQUATN,LINEAR, &
           CCODED,XEPS,FLAG)
  
      ! Actually, we should sum Q(1) to have the correct value
      ! of the model at the points. But, in order to calculate
  