anchor.diff:    470 VMULTC(K)=DMAX1(0.0d0,TEMP*VMULTC(K)+RATIO*VMULTD(K))
         PRINT *, '  3. VMULTC = ', (VMULTC(I),I=1,M+1)
         IF (MCON .EQ. M) RESMAX=RESOLD+RATIO*(RESMAX-RESOLD)
  +      PRINT *, ' RESMAX, MCON, M, ICON = ', 
  +     1    RESMAX, MCON, M, ICON
   C
   C     If the full step is not acceptable then begin another iteration.
   C     Otherwise switch to stage two or end the calculation.
  
positive.code:               TEMP=ALPHA*Z(I,K)+BETA*Z(I,KP)
                Z(I,KP)=ALPHA*Z(I,KP)-BETA*Z(I,K)
    120         Z(I,K)=TEMP
            END IF
            K=K-1
            GOTO 100
  
negative.code:       SUBROUTINE TRSTLP (N,M,A,B,RHO,DX,IFULL,IACT,Z,ZDOTA,VMULTC,
       1  SDIRN,DXNEW,VMULTD,IPRINT)
        IMPLICIT DOUBLE PRECISION (A-H,O-Z)
        DOUBLE PRECISION TEMP
        DIMENSION A(N,*),B(*),DX(*),IACT(*),Z(N,*),ZDOTA(*),
       1  VMULTC(*),SDIRN(*),DXNEW(*),VMULTD(*)
  C
  