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

下载文档

网友关注

课堂练习第十章物质结构答案

2017年甘肃省选拔299名普通高校毕业生参加大学生志愿服务西部计划实施方案

十八届六中全会精神学习体会：贯彻落实十六届六中全会精神,做合格党员

心理咨询师三级职业道德满分-攻略

氧化亚铁硫杆菌作用下形成的黄钾铁矾的SEM研究_朱长见

2015版第八章建设中国特色社会主义总布局第二节Word

中阶管理才能训练系列之六《业务工作执行管理》Word

重庆市标准化考试知识点(根据大纲整理)答案

国际环境下构建新疆形象的新闻话语方式研究_郭越

2014年江苏省中学生生物学奥林匹克初赛理论试题解析(动物学)答案

校园环境绿化美化规划

你真伟大Word

网友关注视频

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

《小学数学二年级下册》第二单元测试题讲解

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

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

苏教版二年级下册数学《认识东、南、西、北》

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

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

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

六年级英语下册上海牛津版教材讲解 U1单词

第五单元民族艺术的瑰宝_15. 多姿多彩的民族服饰_第二课时(市一等奖)(岭南版六年级上册)_T129830

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

青岛版教材五年级下册第四单元（走进军营——方向与位置）用数对确定位置（一等奖）

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

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

冀教版小学数学二年级下册第二单元《余数和除数的关系》

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

沪教版八年级下册数学练习册一次函数复习题B组（P11）

外研版英语三起6年级下册（14版）Module3 Unit1

8.练习八_第一课时(特等奖)(苏教版三年级上册)_T142692

第8课对称剪纸_第一课时(二等奖)(沪书画版二年级上册)_T3784187

人教版历史八年级下册第一课《中华人民共和国成立》

小学英语单词

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

3月2日小学二年级数学下册（数一数）

北师大版数学四年级下册第三单元第二节小数点搬家

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

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

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

第12章圆锥曲线_12.7 抛物线的标准方程_第一课时(特等奖)(沪教版高二下册)_T274713

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

精品推荐

营造良好的企业文化

分类导航

生物学

索胜军

此文档贡献者

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

教育资源为主的文档平台

文档分类