P130-7 色谱分析的计算

上传者：贺昌政
|
上传时间：2015-05-09
|
密次下载

P130-7 色谱分析的计算

北京化工大学计算化学作业答案

PROGRAM MAIN

DOUBLE PRECISION EPS,A,B,X,Y,A1,B1,X1,Y1,S,S1,SS DIMENSION X(10),Y(10),X1(7),Y1(7)

OPEN(6,FILE='OUTPUT',STATUS='UNKNOWN') N=10

M=7

EPS=1E-5

DATA X/38,40,53,64,75,86,97,108,119,130/ DATA Y/0,3,17,40,69,90,68,43,19,4/

DATA X1/250,258,266,274,282,290,298/ DATA Y1/0.5,6,14,21,13,5,0/

A=38

B=130

A1=250

B1=298

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

CALL SIMP(M,A1,B1,X1,Y1,EPS,S1)

SS=S/S1

WRITE(6,11) S,S1,SS

11 FORMAT(/2X,'S=',F18.10/2X,'S1=',F18.10/2X,'The substance(S/S1)is:',f18.10)

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 ratio of content of two

北京化工大学计算化学作业答案

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

下载文档

网友关注

网友关注视频

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

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

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

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

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

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

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

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

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

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

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

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

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

七年级下册外研版英语M8U2reading

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

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

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

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

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

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

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

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

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

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

二年级下册数学第二课

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

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

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

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

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

精品推荐

营造良好的企业文化

分类导航

生物学

贺昌政

此文档贡献者

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

教育资源为主的文档平台

文档分类