输入：数据点数，EPS,T①，C(I)(=1,N)(2分) 曰，NX0T0Y0-C02S2=(4分) 调用Simpson法子程序计算积分S2 (用插值法计算函数值)，S1(I)=S+S2(4分) I=1,NX1①F=S1①/T①Y1①=C①/T④(4分) 调用一元线性回归法计算 k1、k2(4分) 输出：第一步和第二步的反应速率常数k、 (2分) 结束 3.源程序(20分) PROGRAM MAIN (10分) DOUBLEPRECISIONEPS.A.B.X.Y.C.T.S.R.K1.K2.A1.B1.X1.Y1.S2.X2.Y2 DIMENSION X(7.Y(7).C(7)T(7).SI(7.XI(7).YI(7).X2(6).Y2(6) -7 M-6 EPS-1 0E6 OPEN(6.FILE-2009TESTB.Txt',STATUS-'UNKNOWN) DATA TA0.4.7.14.7.24.7.34.7.447.54.71 DATA CA0.0.0170.0.0495.008020.1060.0.1265,0.1360/ Y11=0. X1(1=0 D010=l,N XOETO 10 Y①-CU*2 S24.72)*(0+Y2) X1(2=S24.7 Y1(2=C2y4.7 D015I=3.N A=47 B=T④ 66 3. 源程序（20 分） PROGRAM MAIN （10 分） DOUBLE PRECISION EPS,A,B,X,Y,C,T,S,R,K1,K2,A1,B1,X1,Y1,S2,X2,Y2 DIMENSION X(7),Y(7),C(7),T(7),S1(7),X1(7),Y1(7), X2(6),Y2(6) N=7 M=6 EPS=1.0E-6 OPEN(6,FILE='2009TESTB.Txt',STATUS='UNKNOWN') DATA T/0.,4.7,14.7,24.7,34.7,44.7,54.7/ DATA C/0.,0.0170,0.0495,0.0802,0.1060,0.1265,0.1360/ Y1(1)=0. X1(1)=0 DO 10 I=1,N X(I)=T(I) 10 Y(I)=C(I)**2 S2=(4.7/2)*(0+Y(2)) X1(2)=S2/4.7 Y1(2)=C(2)/4.7 DO 15 I=3,N A=4.7 B=T(I) 输入：数据点数，EPS，T (I), C（I）（I=1，N）（2 分） I=1,N X(I)=T(I) Y(I)=C(I)**2 S2=∫ 4.7 0 2 cHdt （4 分） 调用 Simpson 法子程序计算积分 S2 （用插值法计算函数值），S1（I）＝S+S2 （4 分） ） 输出：第一步和第二步的反应速率常数 k1、k2 （2 分） 结束 I=1,N X1(I)=S1(I)/T(I) Y1(I)=C(I)/T(I) （4 分） 调用一元线性回归法计算 k1、k2 （4 分） ）