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

下载文档

网友关注

Grammar Revision Article

百人一首

译林版5B期中测试卷

以学生日语能力需求为前提提高日语教学的效益

旧版初级标准日本语--课文1-49整理(全套)

China is booming, and its hunger for energy is insatiable.

网友关注视频

苏科版八年级数学下册7.2《统计图的选用》

沪教版八年级下次数学练习册21.4（2）无理方程P19

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

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

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

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

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

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

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

冀教版小学数学二年级下册第二单元《租船问题》

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

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

北师大版八年级物理下册第六章常见的光学仪器（二）探究凸透镜成像的规律

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

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

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

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

外研版英语七年级下册module1unit3名词性物主代词讲解

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

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

第19课我喜欢的鸟_第一课时(二等奖)(人美杨永善版二年级下册)_T644386

沪教版八年级下册数学练习册21.4（1）无理方程P18

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

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

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

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

冀教版小学数学二年级下册第二周第2课时《我们的测量》宝丰街小学庞志荣.mp4

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

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

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

精品推荐

营造良好的企业文化

分类导航

生物学

索胜军

此文档贡献者

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

教育资源为主的文档平台

文档分类