平衡关系式为:G+4)=K (K为Cu(NH)42的不稳定常数) c-x 设 f0)=,+4-K C-x 调用Newton-Raphson法求解方程f(x)=O 取初值x。=0.0 因fx)的导数∫"(x)不便计算,此处用差商代替导数。 2.程序框图(30分) 开始☐ 输入:输入c1,c2,K,x0,EPS 调用Newton-Raphson法子程序求解方程f(x)=0 (其中调用子程序计算f(x)值 输出:x值 结束 3.源程序(30分) PROgram MAIn REALCI.C2.K.X.X0 EXTERNALE DATAC1,C2,K/0.1,0.1,4.6E-14 DATA EPS.E/1E-6.1E-7 COMMON CL.C2.K X0-0.01 CALLNR(X0,F,F.0,EPS,E,X) WRITE(*,*Y[Cu2+]='X.'molL-I' END 平衡关系式为: K c x x c x = − + 1 4 2 ( 4 ) (K 为 Cu(NH3)4 2+的不稳定常数) 设 K c x x c x f x − − + = 1 4 2 ( 4 ) ( ) 调用 Newton-Raphson 法求解方程 f (x) = 0 取初值 x0 = 0.01 因 f (x) 的导数 f ′(x) 不便计算,此处用差商代替导数。 2.程序框图(30 分) 3.源程序(30 分) PROGRAM MAIN REAL C1,C2,K,X,X0 EXTERNAL F DATA C1,C2,K/0.1,0.1,4.6E-14/ DATA EPS,E/1E-6,1E-7/ COMMON C1,C2,K X0=0.01 CALL NR(X0,F,F,0,EPS,E,X) WRITE(*,*)'[Cu2+]= ',X,'molL-1' END 开始 输入:输入 c1,c2,K,x0,EPS 调用 Newton-Raphson 法子程序求解方程 f (x) = 0 (其中调用子程序计算 f (x) 值) 输出:x 值 结束