11 FORMAT(1X,3F10.4) 4 运行结果(15分) CALL GS(A,N,M,EPS) ORIGINAL LINEAR EOUATIONS (5分) R12=EXPA1.3+24AM23-AL3)}◆X12122◆*2) n74 2568 -2070 R22=EXP(A(2,3+2(AM1,3A2,3)X22)*(X12+2》 2428 -.0066 -2459 WRITE(6.33A(LMLEI.N) A=-1034D+01 (5分) 33F0RMAT5X,A=,D12.4,/5X,B=,D12.40 B=-.7764D+00 WRITE(6,44)R12,R22 WHEN X1=0.200 (5分) 44 FORMAT/5X,WHEN X1=0.200,/5X,'r1=,D12.4/5X,2=,D12.40 rl= 5511D+00 END r2=9536D+00 SUBROUTINE GS(A.N,M,EPS) (10分) D DOUBLE PRECISION A,EPS 【,计算原理（化学原理和计算方法）(15分) DIMENSION A(N,M) 某物质恒压下无相变时，其摩尔恒压热容才Cm与温度T之间的关系可表示为： D050K=1,N P=00 D010=K,N Cn小mol.K=b。+bT+b,T++h,lhT IF(P.GE.ABS(A(I,K)))GOTO 10 P=ABS(A(LK)) bo=12.4,6=2.25,b2=0.34,b3=80.6,b=24.3(5分) 10 CONTINUE 系统的拾变AH-CdT (1)(5分) IF(P.LT.EPS)STOPNO SOLUTION IF(L.EQ.K)GOTO 30 DO20 J=K,M TMP=A(LJ) 由于被积函数Ca小molK=+bT+h,T'+ T +b,lnT有确定的形式，所 A(LJ=A(KJ) A(K,J)=TMP 以可直接利用Sim即son求积分的方法计算可得到(1)式右边的积分值，其中调用自编子程序 20 CONTINUE 30 P=1.0/A(K.K) 计算函数Cp,m的位.然后代入(1)式，可得从1OK升温到50K系统的格变△H.(5分) D040j=K+1,M AKJ=AKJ◆P D0401=K+1,N 2.程序框图(20分)》 A(I,JFA(IJ)-A(LKYA(KJ) 开始 40 CONTINUE 50 CONTINUE 输入：bm=12.4,b1=2.25,b2=0.34, b=80.6,b4=24.3(5分)N D060K=2,N =MK D060J-+1,N A=10K,B=50K,(I=1,N) (5分) A(LM=A(LM)A(LJ)A(J.M) 60 CONTINUE 调用Simpson法子程序计算积分S RETURN END (用自编子程序计算函数C,m值)(5分)6 11 FORMAT(1X,3F10.4) CALL GS(A,N,M,EPS) R12=EXP((A(1,3)+2*(A(2,3)-A(1,3))*X12)*(X22**2)) R22=EXP((A(2,3)+2*(A(1,3)-A(2,3))*X22)*(X12**2)) WRITE(6,33)(A(I,M),I=1,N) 33 FORMAT(/5X,'A=',D12.4,/5X,'B=',D12.4/) WRITE(6,44) R12,R22 44 FORMAT(/5X,'WHEN X1=0.200', /5X,'r1=',D12.4,/5X,'r2=',D12.4/) END SUBROUTINE GS(A,N,M,EPS) （10 分） DOUBLE PRECISION A,EPS DIMENSION A(N,M) DO 50 K=1,N P=0.0 DO 10 I=K,N IF(P.GE.ABS(A(I,K)))GOTO 10 P=ABS(A(I,K)) L=I 10 CONTINUE IF(P.LT.EPS)STOP'NO SOLUTION' IF(L.EQ.K)GOTO 30 DO 20 J=K,M TMP=A(L,J) A(L,J)=A(K,J) A(K,J)=TMP 20 CONTINUE 30 P=1.0/A(K,K) DO 40 J=K+1,M A(K,J)=A(K,J)*P DO 40 I=K+1,N A(I,J)=A(I,J)-A(I,K)*A(K,J) 40 CONTINUE 50 CONTINUE DO 60 K=2,N I=M-K DO 60 J=I+1,N A(I,M)=A(I,M)-A(I,J)*A(J,M) 60 CONTINUE RETURN END 4． 运行结果（15 分） ORIGINAL LINEAR EQUATIONS: （5 分） .0074 .2568 -.2070 .2428 -.0066 -.2459 A= -.1034D+01 （5 分） B= -.7764D+00 WHEN X1=0.200 （5 分） r1= .5511D+00 r2= .9536D+00 D 1. 计算原理（化学原理和计算方法）（15 分） 某物质恒压下无相变时，其摩尔恒压热容才 Cp,m与温度 T 之间的关系可表示为： b T T b Cp /J mol K b b T b T ln 4 2 3 0 1 2 -1 -1 ,m ⋅ ⋅ = + + + + b0＝12.4， b1＝2.25， b2＝0.34， b3＝80.6， b4＝24.3 （5 分） 系统的焓变 H C T T T p d 21 Δ = ∫ ,m （1）（5 分） 由于被积函数 b T T b Cp /J mol K b b T b T ln 4 2 3 0 1 2 -1 -1 ,m ⋅ ⋅ = + + + + 有确定的形式，所 以可直接利用 Simpson 求积分的方法计算可得到（1）式右边的积分值，其中调用自编子程序 计算函数Cp,m 的值，然后代入（1）式，可得从 10K 升温到 50K 系统的焓变 ΔH 。（5 分） 2. 程序框图（20 分） 开始 输入：b0＝12.4， b1＝2.25， b2＝0.34， b3＝80.6， b4＝24.3（5 分）N A=10K, B=50K,（I=1，N） （5 分） 调用 Simpson 法子程序计算积分 S （用自编子程序计算函数Cp,m 值） （5 分）