newton.m

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

x0=input('请输入求根起点：');
eps=input('请输入最大误差要求：');
h=0.01;
f0=fun(x0);
f1=(fun(x0+h)-fun(x0))/h;                           
while(abs(f0)>eps)
    x1=x0-f0/f1;
    x0=x1;
    f0=fun(x0);
    f1=(fun(x0+h)-fun(x0))/h;
end
disp(x1);