newton.m

来自「matlab牛顿迭代方法,上机试验对于初学者又很大的帮助,简单易懂」· M 代码 · 共 25 行

%牛顿法

x=1.5;
x1=x-fun(x)/gun(x);
if(abs(x1)<1)
    delt=abs(x1-x);
else
    delt=abs((x1-x)/x1);
end

while(delt>1e-6 & abs(gun(x))>1e-6)
    x=x1;
    x1=x-fun(x)/gun(x);
    if(abs(x1)<1)
        delt=abs(x1-x);
    else
        delt=abs((x1-x)/x1);
    end
    if gun(x1)==0
        break
    end         
end

b=x1;
fprintf('%.6f\n',b)