Z=Z+S*Y④ 30 CONTINUE ETURN END 4.运行结果。(20分) S=3890.804468790690000 S1= 479.166938781738300 The ratio of content of the two substances(S/S1)is: 81199351497 (A/B) C 1.计算原理（化学原理和计算方法）(20分) 反应 其动力学方程式为： d[A]--k,[A]+k:[B] dr d[B =k [A]-k2 [B]-k;[B] dIC]=k[B] dr 利用Runge-Kuta法解此方程组，可得不同时刻的[A,B],[C]值。 2.程序框图(30分) 开始 输入：「A1,B1,C1初值 时间间隔H, 计算精度EPs 计算：积分步长H0 调用变步长Runge--Kuta积分法子程序解动力学方程组计算时刻t 的A,B,[C值 (其中调用计算动力学方程组的子程序计算DY值) 列表输出：10min内每隔0.5min所对应的[A,B,C值 结束 Z=Z+S*Y(I) 30 CONTINUE RETURN END 4.运行结果。（20 分） S= 3890.804468790690000 S1= 479.166938781738300 The ratio of content of the two substances(S/S1) is: 8.1199351497 （A/B） C 1.计算原理（化学原理和计算方法）（20 分） 反应 A ⎯⎯→1 B ⎯⎯→3 C k k 其动力学方程式为： [B] d d[C] [A] [B] [B] d d[B] [A] [B] d d[A] 3 1 2 3 1 2 k t k k k t k k t = = − − = − + 利用 Runge-Kutta 法解此方程组，可得不同时刻的[A]，[B]，[C]值。 2.程序框图（30 分） 开始 输入：[A]，[B]，[C]初值 时间间隔 H， 计算精度 EPS 计算：积分步长 H0 调用变步长 Runge-Kutta 积分法子程序解动力学方程组计算时刻 t 的[A]，[B]，[C]值 （其中调用计算动力学方程组的子程序计算 DY 值） 列表输出：１0min 内每隔 0.5min 所对应的[A]，[B]，[C]值 结束 k2