输入：数据点数N, 反应温度TO, 时间间隔H 时间t与总压P的实验数据T),P),反应物起始压力PO X1)0,Y1)=P0,X=T,Y0F2*P)-P0(I=2,N+1) 调用中心差分子程序计算反应速率R(W)(dp/dr) X(I)=ln[T(D)].Y(D)=ln[-R(D)](I=1.N+I) 调用线性回归子程序计算A,B (反应级数S=B,速率常数KS=EXP(A)N) 输出：A,BS,KS 3.源程序(18分) 结束 PROGRAM MAIN DIMENSION X(11).Y(11).V(11),XI(10),YI(10) DOUBLE PRECISIONX,Y.X1.Y1,V.H.A.B,R.S.SS.P0 OPEN(6,FILE-TESTB.TXT.STATUS-UNKNOWN) M=10 N-lI H=1.0E-5 P0=0.8315E05 DATA X1/8.02,12.18,17.30.24.55,33.00,42.50,55.08,68.05,& 90.05.1190/ DATAY1/0.78660.76860.7464,0.71940.6944,0.6701,0.6450,& 0.6244,0.5964,0.569 X(1=0 Y(1=P0 D0101=2N X0=X1-1 Y1-1)=Y1(-1)*1E05 10 YI=2*Y1L-1)-P CALLCF(N,H.X,Y.V) D01001=1N X(D-LOG(Y(D) Y(I)-LOG(ABS(V 100 CONTINUE CALLPK(N.X.Y.A.B.R) S=B3.源程序（18 分） PROGRAM MAIN DIMENSION X(11),Y(11),V(11),X1(10),Y1(10) DOUBLE PRECISION X,Y,X1,Y1,V,H,A,B,R,S,SS,P0 OPEN(6,FILE='TESTB.TXT',STATUS='UNKNOWN') M=10 N=11 H=1.0E-5 P0=0.8315E05 DATA X1/8.02,12.18,17.30,24.55,33.00,42.50,55.08,68.05, & 90.05,119.0/ DATA Y1/0.7866,0.7686,0.7464,0.7194,0.6944,0.6701,0.6450, & 0.6244,0.5964,0.5694/ X(1)=0 Y(1)=P0 DO 10 I=2,N X(I)=X1(I-1) Y1(I-1)=Y1(I-1)*1E05 10 Y(I)=2*Y1(I-1)-P0 CALL CF(N,H,X,Y,V) DO 100 I=1,N X(I)=LOG(Y(I)) Y(I)=LOG(ABS(V(I))) 100 CONTINUE CALL PK(N,X,Y,A,B,R) S=B 输入：数据点数 N, 反应温度 TO, 时间间隔 H 时间 t 与总压 p 总的实验数据 T(I),P(I)，反应物起始压力 P0（I=1， X(1)=0,Y(1)=P0;X(I)=T(I),Y(I)=2*P(I)-P0 (I=2,N+1) 调用中心差分子程序计算反应速率 R(W)（dp/dt） X(I)=ln[T(I)],Y(I)=ln[-R(I)] (I=1,N+1) 调用线性回归子程序计算 A，B （反应级数 S=B，速率常数 KS=EXP(A)）N） 输出： A,B, S,KS 结束