2.程序框图(30分) 开始 输入：数据点数N,EPS,S1, T(I),Cp(D)(I=1,N) N) X(I=T(D):Y()-Cp(IT()(I=IN):A=TIB=T2 调用Simpson法子程序计算AB面积S (用插值法计算F),S2-S1+S 输出：S2 结束☐ 3.源程序(30分) PROGRAM MAIN DOUBLE PRECISION EPS.A.B.X.Y.S.KP1.KP DIMENSION X(9).Y(9) N=9 EPS=1 0E.6 OPEN(6.FILE-TESTC.Txl'.STATUS-UNKNOWN) DATA X/15,30,50,100,150,200,30,400,550/ DATA Y/7.2,16.4,21.1,24.3,252,25.8,26.5,27.4,28.9/ A=15 B=500 61=24 D010=1,N X=x 10 Y(I)-Y(I/XI CALLSIMP(N.A.B,X,Y,EPS.S) 92=S+91 WRITE(6 20)A BS WRITE(6.0)B.S2 20 FORMAT(IX.'S(.F6.2.K->.F6.2.K)-.G12.6) 30 FORMAT(1X,'S,F7.2,'K)尸，G12.6,'JK-1mol-1D END SUBROUTINESIMP(M.A.B.X.Y.EPS.S2) DOUBLEPRECISIONA.B.X.Y.EPS.S2.FA.FB.X0.Y0 DIMENSION X(M),Y(M) H=B-A 2.程序框图（30 分） 3.源程序（30 分） PROGRAM MAIN DOUBLE PRECISION EPS,A,B,X,Y,S,KP1,KP2 DIMENSION X(9),Y(9) N=9 EPS=1.0E-6 OPEN(6,FILE='TESTC.Txt',STATUS='UNKNOWN') DATA X/15,30,50,100,150,200,300,400,550/ DATA Y/7.2,16.4,21.1,24.3,25.2,25.8,26.5,27.4,28.9/ A=15 B=500 S1=2.4 DO 10 I=1,N XI=X(I) 10 Y(I)=Y(I)/XI CALL SIMP(N,A,B,X,Y,EPS,S) S2=S+S1 WRITE(6,20) A,B,S WRITE(6,30) B,S2 20 FORMAT(1X,'S(',F6.2,'K-->',F6.2,'K)=',G12.6) 30 FORMAT(1X,'S(',F7.2,'K)=',G12.6,'JK-1mol-1') 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,EPS,S1 ， T(I),Cp(I) （I=1，N） N） X(I)=T(I)；Y(I)=Cp(I)/T(I) (I=1,N)；A=T1,B=T2 调用 Simpson 法子程序计算 A-B 面积 S （用插值法计算 F），S2=S1+S 输出：S2 结束