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

下载文档

网友关注

全局变量、局部变量、自动变量、寄存器变量、静态变量、宏定义

信号实验(二)

基于自适应阈值估计的模极大值去噪方法_华春红

双容水箱上的几种液位控制实验及被控对象的数学模型_王志新

光子晶体微腔稀疏波分解复用器的设计与模拟_李旭_乐孜纯

A+Magnetization+Measuring+Instrument—Vibrating+Sample+Magnetometer

浅析矩阵的秩_罗雪梅[1]

2015净水器选哪个品牌好,净水器哪个牌子好

网友关注视频

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

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

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

化学九年级下册全册同步人教版第25集生活中常见的盐（二）

二年级下册数学第二课

二年级下册数学第一课

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

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

北师大版数学四年级下册第三单元第四节街心广场

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

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

小学英语单词

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

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

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

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

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

河南省名校课堂七年级下册英语第一课（2020年2月10日）

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

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

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

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

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

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

8.对剪花样_第一课时(二等奖)(冀美版二年级上册)_T515402

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

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

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

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

七年级下册外研版英语M8U2reading

精品推荐

营造良好的企业文化

分类导航

生物学

贺昌政

此文档贡献者

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

教育资源为主的文档平台

文档分类