newton.m

来自「单变量寻优基本算法的原码」· M 代码 · 共 22 行

%用牛顿法求方程的根
%x(k+1)=x(k)-f(xk)/df(xk)
%不是特别好用，有时不收敛

x1=0;
e=0.0005;

for i=1:200
    y=f(x1);
    x2=x1-y/df(x1);
    dy=df(x2);
    if abs(dy)<e
        break;
    else
        x1=x2;
    end
end

i
x2
x2-x1
f(x2)