单组分压缩因子.for

运用P-R状态方程

	REAL P,T,Z,TC,PC,TR,AT,K
      REAL W,A,B
	OPEN(3,FILE='INDATE1.TXT')
	READ(*,*) P,T
	open(6,file='计算结果1.txt')
	R=8.314
	AM=0.0
      BM=0.0
      READ(3,*)PC,TC,W
	TR=T/TC
	K=0.3764+1.54226*W-0.26992*W**2	
	AT=(1+K*(1-SQRT(TR)))**2             
      A=AT*0.45724*R*R*TC*TC/PC
      B=0.07880*R*TC/PC
   	AM=AM+A**0.5
      BM=BM+B
   10 CONTINUE
	AM=AM**2
      AA=AM*P/(R*T)**2
	BB=BM*P/(R*T)
   20	X1=X
      F=X1**3-(1-BB)*X1**2+(AA-2*BB-3*BB**2)*X1-(AA*BB-BB**2-BB**3)
	F1=3*X1**2-2*(1-BB)*X1+(AA-2*BB-3*BB**2)
	X=X1-F/F1
	Z=X
	WRITE(6,100)X1,X,Z
	IF(ABS(X-X1).GT.1E-6) GOTO 20
  100 FORMAT(2X,'X1=',F15.7,2X,'X=',F15.7,2X,'Z=',F15.7)  
      END