IF(D.LT.EPS)GOTO 50 N-N+N 输出：从1OK升温到50K系统的焓变△H (5分) H=0.5*H T1=12 结束 S1=S2 3.源程序(20分) G0T010 PROGRAM MAIN 50 RETURN DOUBLE PRECISION EPS,TA,TB,S2 (10分) END EXTERNALF 4.运行结果(15分) OPEN(6,FILE-2010TESTD.Txr',STATUS-'UNKNOWN) DH10K->50K=20603.032227Jmol-1 EPS=1.0E-6 TA=10.0 1.计算原理（化学原理和计算方法）(15分) TB=50.0 CALL SIMP(TA,TB,EPS,S2,F) 反应： 2A—→B 10 WRITE(6.1)52 FORMAT('DH(10K->50K-,F12.6,3X,Jmol-1) 已知反应物A和B的初始浓度分别为C。,cm,A的转化率随反应时间()的变化数 END 据。 组分A的速率方程为 FUNCTION F(T) (5分 DOUBLE PRECISION T,F,BO,B1,B2,B3,B4 DATA B0W12.46/,B12.25/,B2/0.34/B3/80.6/.B4/24.3/ v,=de=k,c以 (1)(4分) F=B0+BI*T+B2*T*T+B3/T+B4*LOG(T) dt RETURN END (1)式中，VA为组分A的化学反应速率，ca为反应物A浓度，1为反应时间，n为A的 SUBROUTINE SIMP(A,B,EPS,S2,F) (5分) 反应级数，ka为A的速率常数，dC仙为反应物A浓度随时间的变化率。 DOUBLE PRECISION ABHEPS S2AI H=B-A (1)式两边取对数： T1=O.5H(F(A+F(B》 Inv =Ink +nlnc N=1 s1=0 10 S=0 y=Inv,x Inc,a=Ink,b=n (4分) D020K=1,N A1=A+K-0.5H 20 S=S+F(AI) 得y=a+br T2=0.5*(T1+H*s) S2=T2+(T2-T1)3.0 对应的反应物的浓度©A(,)可由反应关系式求出 cA()= CAd1-0A) (3分 IF(ABS(S2).LE.1.0)D=ABS(S2-S1) 将计算所得同时间1的反应物浓度ca数据通过插值和中心差商法求出反应速率V。·计算 IF(ABS(S2).GT.1.0)THEN 1nV,和1ca后，用一元线性回归子程序计算该反应的级数m,反应的表观速率常数kA·(4分) D=ABS((S2-S1)/S2) END IF 7 3. 源程序（20 分） PROGRAM MAIN DOUBLE PRECISION EPS,TA,TB,S2 （10 分） EXTERNAL F OPEN(6,FILE='2010TESTD.Txt',STATUS='UNKNOWN') EPS=1.0E-6 TA=10.0 TB=50.0 CALL SIMP(TA,TB,EPS,S2,F) 10 WRITE(6,11) S2 11 FORMAT('DH(10K-->50K)=',F12.6,3X,'J mol-1') END FUNCTION F(T) （5 分） DOUBLE PRECISION T,F,B0,B1,B2,B3,B4 DATA B0/12.46/,B1/2.25/,B2/0.34/,B3/80.6/,B4/24.3/ F=B0+B1*T+B2*T*T+B3/T+B4*LOG(T) RETURN END SUBROUTINE SIMP(A,B,EPS,S2,F) （5 分） DOUBLE PRECISION A,B,H,EPS,S2,A1 H=B-A T1=0.5*H*(F(A)+F(B)) N=1 S1=0 10 S=0 DO 20 K=1,N A1=A+(K-0.5)*H 20 S=S+F(A1) T2=0.5*(T1+H*S) S2=T2+(T2-T1)/3.0 IF(ABS(S2).LE.1.0)D=ABS(S2-S1) IF(ABS(S2).GT.1.0)THEN D=ABS((S2-S1)/S2) END IF IF(D.LT.EPS)GOTO 50 N=N+N H=0.5*H T1=T2 S1=S2 GOTO 10 50 RETURN END 4. 运行结果（15 分） DH(10K-->50K)=20603.032227 J mol-1 E 1．计算原理（化学原理和计算方法）（15 分） 反应： 2A ⎯⎯→B 已知反应物 A 和 B 的初始浓度分别为 cA0， cB0， A 的转化率随反应时间（t）的变化数 据。 组分 A 的速率方程为 n k c t c v A A A A dd = − = （1）（4 分） （1）式中， Av 为组分 A 的化学反应速率，cA为反应物 A 浓度，t 为反应时间，n 为 A 的 反应级数，kA为 A 的速率常数，dcA/dt 为反应物 A 浓度随时间的变化率。 （1）式两边取对数： A A A ln v = ln k + nln c 令 y v x c a k b n = ln , 1 = ln A , = ln , = （4 分） 得 y = a + bx 对应的反应物的浓度 cA（ti）可由反应关系式求出 cA（ti）＝ cA0(1-αA) （3 分） 将计算所得同时间 t 的反应物浓度 cA数据通过插值和中心差商法求出反应速率 Av 。计算 ln Av 和 lncA后，用一元线性回归子程序计算该反应的级数 n，反应的表观速率常数 kA。（4 分） 输出：从 10K 升温到 50K 系统的焓变ΔH （5 分） 结束