anchor.diff:           PRINT *, '  VSIG = ', (VSIG(J),J=1,N), ' -- ', PARSIG
            PRINT *, '  VETA = ', (VETA(J),J=1,N), ' -- ', PARETA
            PRINT *, '  IBRNCH, IFLAG = ', IBRNCH, IFLAG
  +         PRINT *, '  A = ', ((A(I,J),I=1,N),J=1,MP)
         END IF
   C
   C     If a new vertex is needed to improve acceptability, then decide which
  
positive.code:      1        'MAXFUN limit has been reached.')
           GOTO 600
        END IF
        IF (IPRINT .EQ. 3) THEN
           PRINT *, '  SIM = ', (SIM(J,NP),J=1,N)
        END IF
        NFVALS=NFVALS+1
        IF (IPRINT .EQ. 3) THEN
           PRINT *, '  DX = ', (DX(I),I=1,N)
        END IF
        IF (IPRINT .EQ. 3) THEN
           PRINT *, '  BEFORE: ', N, M, (X(I),I=1,N), F, (CON(I),I=1,M)
        END IF
        CALL CALCFC (N,M,X,F,CON)
  
negative.code: C     Partition the working space array W to provide the storage that is needed
  C     for the main calculation.
  C
        ITOTAL=N*(3*N+2*M+11)+4*M+6
        MPP=M+2
        ICON=1
        ISIM=ICON+MPP
  