?? 單組分壓縮因子.for
字號:
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
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -