anchor.diff:           PRINT *, '  **CON = ', (CON(I),I=1,M)
            PRINT *, '  **A = ', ((A(I,K),I=1,N),K=1,MP)
            PRINT *, '  **W = ', (W(I),I=1,ITOTAL)
  +         print *, ' '
  +         print *, ' '
         END IF
         IF (IFULL .EQ. 0) THEN
             TEMP=0.0d0
  
positive.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
  
negative.code: C     Further, SIMI holds the inverse of the matrix that is contained in the
  C     first N columns of SIM.
  C
        ITOTAL=N*(3*N+2*M+11)+4*M+6
        IPTEM=MIN0(N,5)
        IPTEMP=IPTEM+1
        NP=N+1
  