d[A] d =-[A]+[B]-%[A]IC] d(B]-&,[A]-k;[B] C]=-k,IAJIC] D]=&,[AJIC] dt 利用Runge-Kutta法解此方程组，可得不同时刻的[A,B],[C],D]值。 2.程序框图(18分) 开始 输入：[A],B,[C],D初值 时间间隔H 计算精度EPS 计算：积分步长H0 调用变步长Runge-.Kuta积分法子程序解动力学方程组计算时刻t 的ABC,D (其中调用计算动力学方程组的子程序计算DY值)》 列表输出：0-200min内每隔10min所对应的A,B),[C,D]值 3.源程序(18分) PROGRAM MAIN IMPLICITREAL*8(A-H.O-Z) DIMENSION YO(10).Y1(10),Y2(10) OPEN(6,FILE-TESTC.TXT,STATUS-UNKNOWN) DATA Y0/2.0,0.0,1.0.0.0.6*0.0/ N=4 T0-0.0 T1=100.0 H0=10 EPS=1.0E-6 H=HO WRITE(6.11)N.T0.T1.H0.EPS WRITE(6.22)T0.H.(Y0(D).1=1.N) M=IDINT((T1-T0VHO+10) D08010=1,M K=1[A][C] d d[D] [A][C] d d[C] [A] [B] d d[B] [A] [B] [A][C] d d[A] 3 3 2 2 1 2 3 2 1 k t k t k k t k k k t = = − = − = − + − 利用 Runge-Kutta 法解此方程组，可得不同时刻的[A]，[B]，[C]，[D]值。 2.程序框图（18 分） 3.源程序（18 分） PROGRAM MAIN IMPLICIT REAL *8(A-H,O-Z) DIMENSION Y0(10),Y1(10),Y2(10) OPEN(6,FILE='TESTC.TXT',STATUS='UNKNOWN') DATA Y0/2.0,0.0,1.0,0.0,6*0.0/ N=4 T0=0.0 T1=100.0 H0=10 EPS=1.0E-6 H=H0 WRITE(6,11)N,T0,T1,H0,EPS WRITE(6,22)T0,H,(Y0(I),I=1,N) M=IDINT((T1-T0)/H0+10) DO 80 I0=1,M K=1 开始 输入：[A]，[B]，[C]，[D]初值 时间间隔 H， 计算精度 EPS 计算：积分步长 H0 调用变步长 Runge-Kutta 积分法子程序解动力学方程组计算时刻 t 的[A]，[B]，[C]，[D]值 （其中调用计算动力学方程组的子程序计算 DY 值） 列表输出：0-200min 内每隔 10min 所对应的[A]，[B]，[C]，[D]值 结束