计算体积比DN) 调用Newton-Raphson法子程序求解方程f(x)=0 (其中用差高近似f"(x)值)得x 计算平衡时C12的摩尔分数YC12 输出：T,P,D,x,YCI2值 结束 3.源程序(30分) PROGRAM MAD (10分) REALT.P.KPNHCL.NO2,YHCL,YA,D.EPS.X0.YCL2.X.E EXTERNAL KP.NHCL.NO2 OPEN(6.FILE-2005TC.TXT.STATUS-UNKNOWN) 320 YHCL=0.355 DATA EPS,E/IE-6,1E-7/ NO2=11+D) 17 11 FORMAT(/5X.T-F104.K/5X.P-.GI2.4.Pa/5X.VHCL:VO2-.F124& PX, P,F6.4//RESULT:/5X,'y(CL2P,F6.4 END 0分) O.NHCL.NO2 COMMON P.KP.NHCL,NO2 KPP(NHCL-2X 2(NHCL-2*X)**2)*N02-0.5◆X =S-0 10 10 3.源程序（30 分） PROGRAM MAIN （10 分） REAL T,P,KP,NHCL,NO2,YHCL,YA,D,EPS,X0,YCL2,X,E EXTERNAL F COMMON P,KP,NHCL,NO2 OPEN(6,FILE='2005TC.TXT',STATUS='UNKNOWN') T=370 P=1E5 YHCL=0.355 YA=0.645 KP=2.225E-4 DATA EPS,E/1E-6,1E-7/ D=YHCL/(YA*0.21) NHCL=D/(1+D) NO2=1/(1+D) X0=0.01 CALL NR(X0,F,F,0,EPS,E,X) YCL2=X/(1-0.5*X) YCL2=YCL2*(1-0.645*0.79) WRITE(6,11)T,P,D,x,YCL2 11 FORMAT(/5X,'T=',F10.4,'K'/5X,'P=',G12.4,'Pa'/5X,'VHCL:VO2=',F12.4 & /5X,'X=',F6.4//'RESULT:'/5X,'y(CL2)=',F6.4 ) END FUNCTION F(X) （10 分） REAL P,KP,X,S,Q,NHCL,NO2 COMMON P,KP,NHCL,NO2 S=X**2*(X**2) *(1-0.5*X) Q=KP*P*(NHCL-2*X)**2*((NHCL-2*X)**2)*(NO2-0.5*X) F=S-Q 调用 Newton-Raphson 法子程序求解方程 f (x) = 0 （其中用差高近似 f ′(x) 值）得 x 输出：T，P，D，x，YCl2 值 结束 计算体积比 DN） 计算平衡时 Cl2 的摩尔分数 YCl2