北京化工大学2005一一2006学年第一学期 《计算化学》期末考试试卷标准答案 1. 计算原理（化学原理和计算方法）(20分) 苯为低沸点组分。按题意将数值代入方程 (1) 式(1)中，F=F。/2,a=2.5,0=0.6 (5分) 整理后得 2.5ln(1-x)-ln(x)+0.7402=0 令f(x)=2.5ln(1-x)-ln(x)+0.7402 (5分) f(x)的一阶导数容易求得 f(x)=-1-xx 251 (5分) 故可调用Newton-Raphson法求解方程f(x)=0 此处，取初值=0.4，可求出x。其中调用函数子程序计算f(x)、f"(x)值。(5分) 程序框图(30分) 开始 输入已知变量：F=F。/2,a=2.5,x0=0.6(5分) 取初值=0.4(5分)N)1 北京化工大学 2005——2006 学年第一学期 《计算化学》期末考试试卷标准答案 A 1． 计算原理（化学原理和计算方法）（20 分） 苯为低沸点组分。按题意将数值代入方程 ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ − − + − = 0 0 0 1 1 ln ln 1 1 ln x x a x x F a F （1） 式（1）中， 0 F F = / 2 ，a＝2.5，x0＝0.6。 （5 分） 整理后得 2.5ln(1 ) ln( ) 0.7402 0 − x x −+ = 令 fx x x ( ) 2.5ln(1 ) ln( ) 0.7402 = −− + （5 分） f ( ) x 的一阶导数容易求得 2.5 1 ( ) 1 f x x x ′ =− − − （5 分） 故可调用 Newton-Raphson 法求解方程 f (x) = 0 此处，取初值＝0.4，可求出 x。其中调用函数子程序计算 f ( ) x 、 f ′(x) 值。（5 分） 2． 程序框图（30 分） 开始 输入已知变量： 0 F F = / 2 ，a＝2.5，x0＝0.6 （5 分） 取初值＝0.4（5 分）N）