anchor.diff:  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
  
positive.code:          PRINT *, '  VSIG = ', (VSIG(J),J=1,N), ' -- ', PARSIG
           PRINT *, '  VETA = ', (VETA(J),J=1,N), ' -- ', PARETA
           PRINT *, '  IBRNCH, IFLAG = ', IBRNCH, IFLAG
        END IF
  C
  C     If a new vertex is needed to improve acceptability, then decide which
  
negative.code:      1        'MAXFUN limit has been reached.')
           GOTO 600
        END IF
        NFVALS=NFVALS+1
        IF (IPRINT .EQ. 3) THEN
           PRINT *, '  SIM = ', (SIM(J,NP),J=1,N)
           PRINT *, '  DX = ', (DX(I),I=1,N)
           PRINT *, '  BEFORE: ', N, M, (X(I),I=1,N), F, (CON(I),I=1,M)
        END IF
        CALL CALCFC (N,M,X,F,CON)
  