DIMENSION Y2(20).DY(20),RK(20) X=XO D051,N 5 Y(IYO(D) D050L=1.M CALL ECYDYY D0101=1N RK(D)-DY(I) 10 Y2=Y+0.5*H*DY CALL F(Y2.DY) D020I=1N RKOD=RKOH2DY(D) 30 Y2(D-Y(H+0.5*H*DY(I) CALLF(Y2.DY D030=1,N RK(D=RK(H2DY(D) Y2I=YID+H*DY① CALLF(Y2,DY) D040I1.N RK①=RK④)+DYO YI=Y④+RK(①*HV6 50 CONTINUE RETURN END 4运行结果。(12分) N=4T0= 0T1=1000H0=100 EpS=100E-05 Y0.=1,23,4N t/min H/min [Al/molL-1 [Bl/molL-1 [Cl/molL-1 [DV/molL-1 10.00 2.00000 000000 1.00000 .000000 10.00 1.5579% 292 98740 125974E.0 20.00 1.2 1.27539 .70206 .977448 .225519E-01 30.00 1.25 1.07941 .889829 969235 307650E-01 40.00 2.50 935702 102655 962253 377469F.01 50.00 825968 11302 956186 .438144E-01 600 73955 12112 95023 4 1770E-01 70.00 669828 1.2761 .94601 .539810E-0 80.00 5.00 612467 132920 941668 ,583320E-01 90.00 5.00 564514 1.37318 937691 623089E-01DIMENSION 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) 40 Y(I)=Y(I)+RK(I)*H/6 50 CONTINUE RETURN END 4.运行结果。（12 分） N= 4 T0= .0 T1= 100.0 H0=10.0 EPS= .100E-05 Y(I),I=1,2,3,4,...,N ---------------------------------------------------------------------------------------------------------------------- t/min H/min [A]/molL-1 [B]/molL-1 [C]/molL-1 [D]/molL-1 --------------------------------------------------------------------------------------------------------------------- .00 10.00 2.00000 .000000 1.00000 .000000 10.00 .63 1.55798 .429421 .987403 .125974E-01 20.00 1.25 1.27539 .702061 .977448 .225519E-01 30.00 1.25 1.07941 .889829 .969235 .307650E-01 40.00 2.50 .935702 1.02655 .962253 .377469E-01 50.00 2.50 .825968 1.13022 .956186 .438144E-01 60.00 2.50 .739551 1.21127 .950823 .491770E-01 70.00 2.50 .669828 1.27619 .946019 .539810E-01 80.00 5.00 .612467 1.32920 .941668 .583320E-01 90.00 5.00 .564514 1.37318 .937691 .623089E-01