弦截法的例题及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)