D030=1,N RK(IERK(IH2.DY(I) 2. 程序框图(20分) 30 Y2(I)=Y(I)+HDY(I) CALL F(Y2.DY) 开始 D0401=1,N RK(I=RK(D+DY(D) Y(I=Y(I)+RK(I*H/6 输入：x1,x11,r1,r2 (4分)N) CONTINUE RETuRN X2=1-x1;x22=1-x11 (4分) END 4运行结果(15分，每列7.5分) 调用GAUSS消去法计算A,B(4分) N=2T0= 0T1= 10.0H0=2.00 EPS=.100E-05 常数A,B及x1=0200,x2=0.800代入原方程组计算出相应的常数 t H/a U r11,22(4分) ,00 2.00 234.000 230.000 2.00 1.00 233999 229.997 4.00 229.995 3. 结束 233.99 源程序(20分 6.00 1.00 233.996 229.992 PROGRAM MAIN (10分) 8.00 100 233.995 229.989 DIMENSION A(2 3 10.00 1.00 233.993 229.987 DOUBLE PRECISION A,EPS,X11,X21,R11,R21,X12,X22,R12,R22 INTEGER I,N,M OPEN(6,FILE-2010TESTC.Txt.STATUS-UNKNOWN) 计算原理（化学原理和计算方法）(15分) N=2 按照玛格物斯方程，二组分液态混合物中两个组分的活度系数，上与浓度，2间存在 M=N+1 如下关系： EPS=1.0D-06 X11=0.486 In=[A+2(B-A)xx (1) X21=1-X11 R11=0.813 R21=0.782 n5=[B+2(A-B)x] X12=0.200 22=1-X12 将上述方程组(1)变形为： A(1,M)=LOG(R11) A(2,M)=LOG(R21) lnr=[x(1-2x】A+2xx2B (2) A41,1=X21*2)(1-2*X11) Inr =2x,XA+[x(1-2x2)]B AM2,1=2*X21*X112 AM1,2=2*X11*X21**2 将=0.486， 1=0.813,=0.782,=1-1 代入(2)。买用高临消去法星方程阁 A(2,2=(X11*2)*(1-2X21) 可得到4，B的值.然后将x1=0.200,=1-1代入(1)，可求得相应的丙酮和氯仿的活度系数 WRITE(6,YORIGINAL LINEAR EQUATIONS: n,乃的值. WRITE(6.11(A(JEI.N+11EI.N) 5 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.运行结果（15 分，每列 7.5 分） N= 2 T0= .0 T1= 10.0 H0=2.00 EPS= .100E-05 ------------------------------------------------------------------ t/a H/a U Th ------------------------------------------------------------------ .00 2.00 234.000 230.000 2.00 1.00 233.999 229.997 4.00 1.00 233.997 229.995 6.00 1.00 233.996 229.992 8.00 1.00 233.995 229.989 10.00 1.00 233.993 229.987 ------------------------------------------------------------------ C 1． 计算原理（化学原理和计算方法）（15 分） 按照玛格勒斯方程，二组分液态混合物中两个组分的活度系数 r1，r2与浓度 x1，x2间存在 如下关系： ㏑ 2 1 12 r A B Ax x =+ − [ 2( ) ] （1） ㏑ 2 2 21 r B A Bx x =+ − [ 2( ) ] 将上述方程组（1）变形为： r x x A x x B r x x A x x B ln 2 [ (1 2 )] ln [ (1 2 )] 2 2 2 1 2 2 2 1 2 1 1 2 2 1 2 = + − = − + （2） 将 x1＝0.486， r1＝0.813， r2＝0.782，x2＝1-x1 代入（2），采用高斯消去法解方程组， 可得到 A，B 的值。然后将 x1＝0.200，x2＝1-x1代入（1）,可求得相应的丙酮和氯仿的活度系数 r1，r2的值。 2． 程序框图（20 分） 3． 源程序（20 分） PROGRAM MAIN （10 分） DIMENSION A(2,3) DOUBLE PRECISION A,EPS,X11,X21,R11,R21,X12,X22,R12,R22 INTEGER I,N,M OPEN(6,FILE='2010TESTC.Txt',STATUS='UNKNOWN') N=2 M=N+1 EPS=1.0D-06 X11=0.486 X21=1-X11 R11=0.813 R21=0.782 X12=0.200 X22=1-X12 A(1,M)=LOG(R11) A(2,M)=LOG(R21) A(1,1)=(X21**2)*(1-2*X11) A(2,1)=2*X21*X11**2 A(1,2)=2*X11*X21**2 A(2,2)=(X11**2)*(1-2*X21) WRITE(6,*)'ORIGINAL LINEAR EQUATIONS:' WRITE(6,11)((A(I,J),J=1,N+1),I=1,N) 开始 输入：x1,x11, r1,r2 （4 分）N） X2=1-x1; x22=1-x11 （4 分） 调用 GAUSS 消去法计算 A，B （4 分） ） 输出：r11，r22 （4 分） 结束 常数 A，B 及 x1=0.200， x2=0.800 代入原方程组计算出相应的常数 r11，r22 （4 分）