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

下载文档

网友关注

福州话学习日常用语200句二栏版本配MP3

宁都县大沽中心小学校园广播6分区U卡机说明书(内页)

2015.4.6时政

杭州市事业单位2014年公开招聘工作人员考试

广西2015年公务员面试题型解读

宁波市事业单位2014年公开招聘工作人员考试

2015年广西公务员面试必考题型

俄罗斯女人名小名或爱称

few, little, a few, a little

玉溪市事业单位2014年公开招聘工作人员考试

申论答题卡_打印版(练习用)

工业设计专业英语第三版第十课翻译

网友关注视频

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

【部编】人教版语文七年级下册《过松源晨炊漆公店（其五）》优质课教学视频+PPT课件+教案，江苏省

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

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

冀教版小学数学二年级下册1

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

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

沪教版八年级下册数学练习册21.3（3）分式方程P17

第4章幂函数、指数函数和对数函数（下）_六指数方程和对数方程_4.7 简单的指数方程_第一课时(沪教版高一下册)_T1566237

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

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

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

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

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

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

每天日常投篮练习第一天森哥打卡上脚 Nike PG 2 如何调整运球跳投手感？

二年级下册数学第三课搭一搭⚖⚖

苏科版数学八年级下册9.2《中心对称和中心对称图形》

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

【部编】人教版语文七年级下册《过松源晨炊漆公店（其五）》优质课教学视频+PPT课件+教案，辽宁省

沪教版八年级下册数学练习册21.3（2）分式方程P15

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

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

3.2 数学二年级下册第二单元表内除法（一）整理和复习李菲菲

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

二年级下册数学第一课

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

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

人教版二年级下册数学

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

精品推荐

营造良好的企业文化

分类导航

生物学

索胜军

此文档贡献者

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

教育资源为主的文档平台

文档分类