p130-6 实际气体逸度的计算

上传者：索胜军
|
上传时间：2015-05-10
|
密次下载

p130-6 实际气体逸度的计算

PROGRAM MAIN

DOUBLE PRECISION EPS,P,FI,FF,A,B,X0,Y0,X,Y,S DIMENSION X0(14),Y0(14),X(14),Y(14)

OPEN(2,FILE='INPUT.TXT',STATUS='OLD') READ(2,*) N,EPS

READ(2,*) (X0(I),I=1,N)

READ(2,*) (Y0(I),I=1,N)

READ(2,*) M

WRITE(*,22)

22 FORMAT(3X,'P/P0',10X,'FI',10X,'FF/P0'/1X,46(1H-)) DO 111 J=1,M

READ(2,*) P

DO 100 I=1,N

X(I)=X0(I)

100 Y(I)=(Y0(I)-1)/X0(I)

A=1E-4

B=P

CALL SIMP(N,A,B,X,Y,EPS,S)

FI=EXP(S)

FF=FI*B

WRITE(*,11) P,FI,FF

11 FORMAT(1X,3G12.5)

111 CONTINUE

END

SUBROUTINE SIMP(M,A,B,X,Y,EPS,S2)

DOUBLE PRECISION A,B,X,Y,EPS,S2,FA,FB,X0,Y0 DIMENSION X(M),Y(M)

H=B-A

N=1

CALL LGRG2(X,Y,M,A,FA)

CALL LGRG2(X,Y,M,B,FB)

T1=0.5*H*(FA+FB)

10 S=0

DO 20 K=1,N

X0=A+(K-0.5)*H

CALL LGRG2(X,Y,M,X0,Y0)

20 S=S+Y0

T2=0.5*(T1+H*S)

S2=T2+(T2-T1)/3.0

IF(N.EQ.1) GOTO 40

30 D=ABS(S2-S1)

IF(ABS(S2).GT.1.0)THEN

D=ABS((S2-S1)/S2)

END IF

IF(D.LT.EPS) GOTO 50

40 N=N+N

H=0.5*H

T1=T2

S1=S2

GOTO 10

50 RETURN

END

SUBROUTINE LGRG2(X,Y,N,T,Z) DIMENSION X(N),Y(N)

DOUBLE PRECISION X,Y,T,Z,S

Z=0.0

IF(N.LE.0)RETURN

IF(N.EQ.1)THEN

Z=Y(1)

RETURN

END IF

IF(N.EQ.2)THEN

Z=(Y(1)*(T-X(2))-Y(2)*(T-X(1)))/(X(1)-X(2)) RETURN

END IF

IF(T.LE.X(2))THEN

K=1

M=3

ELSE IF(T.GE.X(N-1))THEN

K=N-2

M=N

ELSE

K=1

M=N

10 IF(IABS(K-M).NE.1)THEN

L=(K+M)/2

IF(T.LT.X(L))THEN

M=L

ELSE

K=L

END IF

GOTO 10

END IF

IF(ABS(T-X(K)).LT.ABS(T-X(M))) THEN K=K-1

ELSE

M=M+1

END IF

Z=0.0

DO 30 I=K,M S=1.0

DO 20 J=K,M IF(J.NE.I)THEN

S=S*(T-X(J))/(X(I)-X(J)) END IF

20 CONTINUE Z=Z+S*Y(I) 30 CONTINUE RETURN

END

下载文档

网友关注

宁夏公务员面试结构化面试模拟题答案2.28

宁夏公务员行测每日一练资料分析练习题答案02.11

宁夏公务员行测每日一练资料分析练习题02.08

宁夏公务员申论每周一练答案：乡村振兴人才先行

宁夏公务员行测每日一练言语理解练习题02.07

宁夏公务员行测每日一练资料分析练习题答案02.08

宁夏公务员面试每日一练结构化面试模拟题答案2.26

宁夏公务员面试每日一练结构化面试模拟题答案2.22

宁夏公务员行测言语理解练习题02.27

宁夏公务员行测每日一练判断推理练习题答案02.09

宁夏公务员面试每日一练结构化面试模拟题答案2.5

宁夏公务员行测每日一练言语理解练习题答案02.24

宁夏公务员行测每日一练判断推理练习题02.23

宁夏公务员申论每周一练答案：东北振兴需要软环境

宁夏公务员面试每日一练结构化面试模拟题2.7

宁夏公务员行测每日一练资料分析练习题02.26

宁夏公务员申论每周一练：乡村振兴人才先行

宁夏公务员行测数量关系练习题答案02.28

宁夏公务员面试每日一练结构化面试模拟题答案2.8

宁夏公务员面试结构化面试模拟题2.28

宁夏公务员行测每日一练资料分析练习题02.11

宁夏公务员行测每日一练数量关系练习题02.06

宁夏公务员行测数量关系练习题02.28

宁夏公务员行测每日一练言语理解练习题答案02.07

宁夏公务员面试每日一练结构化面试模拟题答案2.11

宁夏公务员面试每日一练结构化面试模拟题2.11

宁夏公务员行测每日一练数量关系练习题答案02.12

宁夏公务员申论每周一练：东北振兴需要软环境

宁夏公务员行测每日一练判断推理练习题02.22

宁夏公务员行测每日一练判断推理练习题答案02.05

网友关注视频

沪教版牛津小学英语（深圳用）四年级下册 Unit 4

人教版二年级下册数学

七年级英语下册上海牛津版 Unit5

沪教版八年级下册数学练习册20.4（2）一次函数的应用2P8

苏科版数学八年级下册第八章第二节可能性的大小

【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案，安徽省

冀教版英语三年级下册第二课

【部编】人教版语文七年级下册《老山界》优质课教学视频+PPT课件+教案，安徽省

七年级英语下册上海牛津版 Unit9

【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案，安徽省

苏科版数学七年级下册7.2《探索平行线的性质》

外研版八年级英语下学期 Module3

沪教版牛津小学英语（深圳用）四年级下册 Unit 8

30.3 由不共线三点的坐标确定二次函数_第一课时(市一等奖)(冀教版九年级下册)_T144342

19 爱护鸟类_第一课时(二等奖)(桂美版二年级下册)_T3763925

冀教版小学数学二年级下册第二单元《有余数除法的整理与复习》

8 随形想象_第一课时(二等奖)(沪教版二年级上册)_T3786594

化学九年级下册全册同步人教版第18集常见的酸和碱（二）

外研版英语七年级下册module3 unit2第一课时

飞翔英语—冀教版(三起)英语三年级下册Lesson 2 Cats and Dogs

冀教版小学数学二年级下册第二单元《有余数除法的竖式计算》

冀教版小学英语四年级下册Lesson2授课视频

【部编】人教版语文七年级下册《泊秦淮》优质课教学视频+PPT课件+教案，辽宁省

《空中课堂》二年级下册数学第一单元第1课时

第五单元民族艺术的瑰宝_16. 形形色色的民族乐器_第一课时(岭南版六年级上册)_T3751175

冀教版英语五年级下册第二课课程解读

化学九年级下册全册同步人教版第22集酸和碱的中和反应（一）

第五单元民族艺术的瑰宝_16. 形形色色的民族乐器_第一课时(岭南版六年级上册)_T1406126

冀教版小学数学二年级下册第二单元《有余数除法的简单应用》

三年级英语单词记忆下册（沪教版）第一二单元复习

精品推荐

营造良好的企业文化

分类导航

生物学

索胜军

此文档贡献者

2015-05-10
贡献时间
密
下载次数

教育资源为主的文档平台

文档分类