欲计算反应器内每隔1.0m转化率及温度分布，需建立微分方程组进行求 解，这是一个常微分方程组的初值问题： =1.1094×10" 1- 1409% -)e d 1000-0.5 m dr=76120x109(100-05x 1-x 1409 -9.9282(T-T)Km 20-0m,x0=0,T。=613K (10分) 调用Runge-Kutta法解上述常微分方程组，可得到反应器内每隔1.0m转化 率及温度分布。 (5分) 2.程序框图(20分) 开始 输入：x;,T初值，TwT0 管长间隔，计算精度EPS (5分 计算：积分步长H(5分) 调用Runge-Kuta法子程序解方程组计算反应器内每隔1.0m转 化率及温度变化。 (5分) (自编子程序计算F(但) 列表输出：L,XT(5分) 结束 3源程序(20分) PROGRAM MAIN (10分) IMPLICIT REAL8(A-H.O-Z) DIMENSION YO(10).Y1(10).Y2(10)10 欲计算反应器内每隔 1.0 m 转化率及温度分布，需建立微分方程组进行求 解，这是一个常微分方程组的初值问题： 1 w 14098 12 1 14098 11 )e 9.9282( ) K m 1000 0.5 1 7.6120 10 ( d d )e m 1000 0.5 1 1.1094 10 ( d d − − − − − − ⋅ − − = × − − = × T T x x z T x x z x T T z0=0 m, x0＝0, K T0 = 613 （10 分） 调用 Runge-Kutta 法解上述常微分方程组，可得到反应器内每隔 1.0 m 转化 率及温度分布。 （5 分） 2.程序框图（20 分） 3.源程序（20 分） PROGRAM MAIN （10 分） IMPLICIT REAL *8(A-H,O-Z) DIMENSION Y0(10),Y1(10),Y2(10) 开始 输入：x；z，T 初值，Tw, T0 管长间隔， 计算精度 EPS （5 分） 计算：积分步长 H（5 分） 调用 Runge-Kutta 法子程序解方程组计算反应器内每隔 1.0 m转 化率及温度变化。 （5 分） （自编子程序计算 F (z)） 列表输出：L，x, T （5 分） 结束