C     THIS IS THE EXAMPLE PROGRAM DESCRIBING HOW TO USE KDB ROUTINES.
      PROGRAM HC_PVP_E
      IMPLICIT DOUBLE PRECISION (A-H,O-Z) 
      INTEGER NCMP
      CHARACTER*50 CNAME
      CHARACTER*20 FORMULA 
      COMMON/HC_NAME/ICMP(50),CNAME(50),FORMULA(50),ICLASS(50) 
      COMMON/HC_PROP/WT(50),TB(50),TF(50),TC(50),PC(50),VC(50),ZC(50)
	1    ,ACCF(50),WSRK(50),VEST(50),ZRA(50),SOLP(50),VOLP(50)
     2    ,QI(50),RI(50),DM(50) 
      COMMON/HC_KPVP/IEQ_PVP(50),PVP_A(50),PVP_B(50),PVP_C(50),
     1	PVP_D(50),PVP_E(50),PVP_F(50),PVP_G(50),PVP_T1(50),PVP_T2(50)
	DIMENSION TEXP(5),PEXP(5)
	DATA (TEXP(J),J=1,5) /328.03,395.676,457.45,521.1,554.0/
	DATA (PEXP(J),J=1,5) /0.00142D2,0.0781D2,0.7451D2
	1						,3.511D2,6.499D2/
	OPEN(UNIT=5,FILE="RESULT.TXT")
C     NUMBER OF COMPONENT = 1
      NCMP = 1
C     SPECIFY THE ID OF EACH COMPONENT 
      ICMP(1) = 841		! 1-OCTANOL
	CNAME(1) = "1-OCTANOL"
C	VAPOR PRESSURE COEFFICIENT FOR 1-OCTANOL
	IEQ_PVP(1) = 1		
      PVP_A(1)=-2.732603D1
	PVP_B(1)=-1.632739D4
	PVP_C(1)=2.048012D2
	PVP_D(1)=1.224813D-5
	PVP_T1(1)=257.65D0
	PVP_T2(1)=652.50D0
C	PURE PROPERTIES FOR 1-OCTANOL
	TB(1)=468.31D0			
	TC(1)=652.5D0			 
 	PC(1)=2777.0D0
	ACCF(1)=0.587D0
	ICLASS(1)=11

	WRITE(5,*) CNAME(1)
	WRITE(5,*) 'COMPARISON OF EXPERIMENTAL AND ESTIMATED VALUE ',
	1		'(% ERROR)'
	WRITE(5,*) '-----------------------------------------------------'
	WRITE(5,*) ' T(K)   PVP(EXP)   KDBEQN    AMBROSE     RIEDEL'
	WRITE(5,*) '-----------------------------------------------------'
	
	DO I = 1, 5
		CALL HC_PVPEQN(1,TEXP(I),PVP1,IST1)
		CALL HC_PVPAW(1,TEXP(I),PVP2,IST2)
		CALL HC_PVPRD(1,TEXP(I),PVP3,IST3)
		ERR1 = DABS( (PVP1-PEXP(I))/PEXP(I) )*100.D0
		ERR2 = DABS( (PVP2-PEXP(I))/PEXP(I) )*100.D0
		ERR3 = DABS( (PVP3-PEXP(I))/PEXP(I) )*100.D0
		WRITE(5,5) TEXP(I),PEXP(I),ERR1,ERR2,ERR3
	ENDDO
5	FORMAT(1x,F6.2,1X, 4(1X,G10.5)) 
	WRITE(5,*) '-----------------------------------------------------'
	CLOSE(5) 
      STOP 
      END