C EQUIV1    SOURCE    CHAT      05/01/12    23:44:44     5004           
C EQUIV1     SOURCE    INSL       24/10/96    
      SUBROUTINE EQUIV1(DEFR0,EPSEQC,EQSTRC,EX,RB,EPO1,ICAL)
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C
C            ********************************************
C            **  CALCUL  APPROCHE DE EPSEQ1 ET EQSTR1 ***
C            ** COURBE DE COMPRESSION  INSA-HYPERBOLE ***
C            ********************************************
C
      IF(ICAL.EQ.0) THEN
        NBIT=200
        PREC=1.D-8
        EPEQ0=DABS(DEFR0)
        EPEQ1=EPO1
        DO 10 I=1,NBIT
          EPEQ=EPEQ0+(EPEQ1-EPEQ0)/2.D0
          SEQ1=EX*(EPEQ-DABS(DEFR0))
          SEQ2=EX*EPEQ/(1.D0+(EX/RB-2.D0/EPO1)*EPEQ+(EPEQ/EPO1)**2)
          CONT=DABS((SEQ1-SEQ2)/SEQ1)
          IF(CONT.LT.PREC) GOTO 20 
          IF(SEQ1.GT.SEQ2) THEN
            EPEQ1=EPEQ
          ELSE
            EPEQ0=EPEQ
          ENDIF
 10     CONTINUE
 20     CONTINUE
      ELSE
        EPEQ=EPO1
        SEQ2=RB
      ENDIF
C      
      IF(EPEQ.GT.EPO1) EPEQ=EPO1
      IF(SEQ2.GT.RB) SEQ2=RB
      EPSEQC=-1.D0*EPEQ
      EQSTRC=-1.D0*SEQ2
      DEFR0 = EPSEQC-EQSTRC/EX
C      
      RETURN
      END