D060J=1N 60 6.22T0,HY20,=1,N 80 FORMAT(/Ix.N-.I2.2X.T0-.F6.1.2X.Tl-F8.12X. H0- G8.3.lX.EPS- 103 234 1X,96(1 3X Dmoll-1 6X'IBVmolL-16X'IC 1X961H-) FORMATUIX 2E8 22X 4G156) 3 FORMAT(/1X,'FAILED TO FIND STABLE SOLUTION IN MAX IT/) 4 END SUBROUTINE FCYDYY (5分) IMPLICIT REAL *8(A-H.O-Z) DEAI RKIK2K2 DIMENSION YUODYGON DATA K10.13/.K2J0.0491/.K3/0.11/ DY(1=K1*Y1)*Y2+K2*Y(3) DYCEKPY0Y(2)+K2*YG) DY(3=K1*YI)Y(2K2Y3K3*Y(3) DY4=K3*Y3) RETURN END SUBROUTINERK4(X0.Y0,Y.H,N.M) (10分) IMPLICIT REAL *8(A-H.O-Z) DIMENSION YO(N).Y(N) DIMENSION Y2(20).DY(20).RK(20) X=X0 D05I=1,N Y(IEYOD) DO 50L=1 M CALL ECYDYY Do101=1N RK(EDY(D 10 Y2(=Y①0.5*H*DY(D CALL F(Y2,DY) DO20I=IN RKI=RK①+2*DYD) 20 Y2=Y(①+0.5*H*DY④ CALL F(Y2,DY) D0301=1N RKI=RK①+2*DY① 30 Y2(D)=Y(D)+H*DY(I) CALL F(Y2,DY) D040=1,N RK(I-RK①HDYO5 DO 60 J=1,N 60 Y0(J)=Y2(J) WRITE(6,22)T0,H,(Y2(I),I=1,N) 80 CONTINUE WRITE(6,44) 11 FORMAT(/1X,'N=',I2,2X,'T0=',F6.1,2X,'T1=',F8.1,2X, & 'H0=',G8.3,1X,'EPS=',G10.3//20X,'Y(I),I=1,2,3,4,...,N'//& 1X,96(1H-)/5X,'t/min',3X,'H/min',6X,'[A]/molL-1', & 6X,'[B]/molL-1',6X,'[C]/molL-1',6X,'[D]/molL-1'& /1X,96(1H-)) 22 FORMAT(1X,2F8.2,2X,4G15.6) 33 FORMAT(/1X,'FAILED TO FIND STABLE SOLUTION IN MAX IT'/) 44 FORMAT(1X,96(1H-)) END SUBROUTINE F(Y,DY) （5 分） IMPLICIT REAL *8(A-H,O-Z) REAL *8 K1,K2,K3 DIMENSION Y(10),DY(10) DATA K1/0.13/,K2/0.0491/,K3/0.11/ DY(1)=-K1*Y(1)*Y(2)+K2*Y(3) DY(2)=-K1*Y(1)*Y(2)+K2*Y(3) DY(3)=K1*Y(1)*Y(2)-K2*Y(3)-K3*Y(3) DY(4)=K3*Y(3) RETURN END SUBROUTINE RK4(X0,Y0,Y,H,N,M) （10 分） IMPLICIT REAL *8(A-H,O-Z) DIMENSION Y0(N),Y(N) DIMENSION Y2(20),DY(20),RK(20) X=X0 DO 5 I=1,N 5 Y(I)=Y0(I) DO 50 L=1,M CALL F(Y,DY) DO 10 I=1,N RK(I)=DY(I) 10 Y2(I)=Y(I)+0.5*H*DY(I) CALL F(Y2,DY) DO 20 I=1,N RK(I)=RK(I)+2*DY(I) 20 Y2(I)=Y(I)+0.5*H*DY(I) CALL F(Y2,DY) DO 30 I=1,N RK(I)=RK(I)+2*DY(I) 30 Y2(I)=Y(I)+H*DY(I) CALL F(Y2,DY) DO 40 I=1,N RK(I)=RK(I)+DY(I)