xianjie.m

来自「一些数值计算用的MATLAB代码,适合初学MATLAB和数字计算的人使用,根据M」· M 代码 · 共 26 行

x0=input('请输入区间起点：');
x1=input('请输入区间终点：');
eps=input('精度要求：');
n=input('运算次数：');
test=0;
f0=fun(x0);
f1=fun(x1);
for i=1:n
    x2=x1-(x1-x0)*f1/(f1-f0);
    f2=fun(x2);
    if abs(f2)<eps
        break;
    end
    test=test+1;
    x0=x1;
    x1=x2;
    f0=f1;
    f1=f2;
end
if test<n
    disp(x2);
else
    disp('请重新选择区间进行求解！！！');
end