弦截法的例题及MATLAB程序

用弦截法求方程 在0。4附近的近似根,初始值为 , ,精确至4位有效数字。

建立如下的XianjieMethod。m --75982

function f=XianjieMethod(x0,x1)

x2=x1-chap1_fun2(x1)*(x1-x0)/(chap1_fun2(x1)-chap1_fun2(x0));

eps=1e-4;

n=0;

fprintf('迭代次数  x_n    feval(x_n)n')

fprintf('n=%3。0f   x_%d=%10。5f   %10。6en',n,n,x0,chap1_fun2(x0))

while abs(x1-x0)>eps&(n<600)

    x0=x1;

    x1=x2;

    x2=x1-chap1_fun2(x1)*(x1-x0)/(chap1_fun2(x1)-chap1_fun2(x0));

    n=n+1;

    fprintf('n=%3。0f   x_%d=%10。5f   %10。6en',n,n,x0,chap1_fun2(x0))

end

fprintf('n迭代次数n=%3。0f  x*=%10。5fn',n,x0)

再建立chap1_fun2

function f=chap1_fun2(x)

f=x*(x+1)^2-1;

在命令窗口调用

>> XianjieMethod (0。4,0。6)