输出：总塔板数SW (3分) 结束 3· 源程序(20分) PRogram MAIN (10分) ,A2,B2,X,Y,Y1,Y2,S1,S2,R1,R2,SN 178.0.2 00 B2-0 R1=5.0 R2=40 10 CALL SOD-XO) SN-S14S2 VRE6.20) SI S2 SN 20 FORMAT(IX.S1-.5X.F6.2.5X.'S2-,G12.6//5X,THE NUMBEROFTHEORETICALPLATES-.G12.6) END SUBROUTINE SIMPIMA BX YEPS S2) (5分) DOUBLEPRECISIONA,B.X.Y.EPS.S2.FA.FB.X0,Y0 DIMENSION X(M).Y(M) H=B-A N=I CALL LGRG2CX YMA FA) CALLLGRG2(X.Y.M.B.FB) TI=0.5*H*(FA+FB) 10 g=0 D020K=1.N X0=A4K.051*I CALLLGRG2(X.Y.M.X0.Y0) 20 =g+y0 T2=05T14H 92=T2+T2.T1V30 IF(N.EQ.D)GOTO 40 0 D=ABS(S2-S1) IF(ABS(S2)GT LO)THEN 77 3． 源程序（20 分） PROGRAM MAIN （10 分） DOUBLE PRECISION EPS,A1,B1,A2,B2,X,Y,Y1,Y2,S1,S2,R1,R2,SN DIMENSION X(6),Y(6), Y1(6), Y2(6) N=6 EPS=1.0E-4 OPEN(6,FILE='2008TESTB1.Txt',STATUS='UNKNOWN') DATA X/0.178,0.275,0.372,0.456,0.650,0.844/ DATA Y/0.243,0.382,0.518,0.616,0.795,0.931/ A1=0.4 B1=0.9 A2=0.15 B2=0.4 R1=5.0 R2=4.0 DO 10 I=1,N Y1(I)=1/(Y(I)-X(I)-(0.9-Y(I))/R1) 10 Y2(I)=1/(Y(I)-X(I)-(Y(I)-0.15)/R2) CALL SIMP(N,A1,B1,X,Y1,EPS,S1) CALL SIMP(N,A2,B2,X,Y2,EPS,S2) SN=S1+S2 WRITE(6,20) S1,S2,SN 20 FORMAT(1X,'S1=',5X,F6.2,5X,'S2=',G12.6,//5X,'THE NUMBER OF THEORETICAL PLATES =',G12.6) END SUBROUTINE SIMP(M,A,B,X,Y,EPS,S2) （5 分） DOUBLE PRECISION A,B,X,Y,EPS,S2,FA,FB,X0,Y0 DIMENSION X(M),Y(M) H=B-A N=1 CALL LGRG2(X,Y,M,A,FA) CALL LGRG2(X,Y,M,B,FB) T1=0.5*H*(FA+FB) 10 S=0 DO 20 K=1,N X0=A+(K-0.5)*H CALL LGRG2(X,Y,M,X0,Y0) 20 S=S+Y0 T2=0.5*(T1+H*S) S2=T2+(T2-T1)/3.0 IF(N.EQ.1)GOTO 40 30 D=ABS(S2-S1) IF(ABS(S2).GT.1.0)THEN 输出：总塔板数 SN （3 分） 结束