南开大学计算物理讲义3.ppt

* 4. Secant method 割线, 正切 * Convergence on the root using the secant method * Divergence using the secant method * * =1.618 =0.618 * 例：用割线法求方程 在x = 1.5附近的根，使绝对误差精确到10-4。 按上式计算得： 解：取初值x0 = 1.5，x1 = 1.4，按公式 得迭代格式 x2 = 1.33522；x3 = 1.32541； x4 = 1.32472； x5 = 1.32472 取x* ? 1.32472 * f(x)=exp(x)*ln(x)-x2 f(1)=-1 f(2)=1 x?[1,2] A=1.0 B=2.0 DX=(B-A)/10.0 X0=(A+B)/2.0 WRITE(*,*) ISTEP,X0,DX STOP END Call secant(DL,X0,DX,ISTEP) FUNCTION F(X) F=exp(x)*ALOG(x)-X*X RETURN END DL=1.0E-06 * DO WHILE(ABS(DX).GT.DL) D= F(X1)- F(X0) X2=X1-(X1-X0)* F(X1)/D X0=X1 X1=X2 DX=X1-X0 X1=X0+DX ISTEP=0 I