牛顿法及牛顿下山法求零点

牛顿法

difffun该函数的导函数

a在a附近求零点

ep容忍误差

k迭代次数

function [x_star,k]=MyNewton(fun,difffun,a,ep)


k=1
x_k=a
f_xk=feval(fun,x_k)
fdiff_xk=feval(difffun,x_k)
x_k1=x_k-f_xk/fdiff_xk

    while abs(x_k-x_k1)>ep
        k=k+1
        x_k=x_k1
        f_xk=feval(fun,x_k)
        fdiff_xk=feval(difffun,x_k)
        x_k1=x_k-f_xk/fdiff_xk
    end
    x_star=x_k1

end

牛顿下山法

function [x_star,k]=MyNewton_xiashan(fun,difffun,a,ep)


k=1
x_k=a
f_xk=feval(fun,x_k)
fdiff_xk&