steffensen.m

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

%------------init------------------
x0=input('近似值：');
delta1=input('精度要求：');
delta2=3*delta1;
%---------------begin-------------
i=1;
while abs(x0-f(x0))>delta1
    printf('这是第%d个解\n',i)
    y=f(x0);
    z=f(y);
    if abs(z-2*y+x0)<delta2
        continue;
    else x1=x0-(y-x0)^2/(z-2*y+x0);
        x0=x1;
    end
    i=i+1;
    x1
end
disp('结果是: ');
x1