tc_easy_newton.m

来自「％ Nevill等插值算法地实现 不会花费大家的时间的 不希望其他人没有帐号」· M 代码 · 共 20 行

function x=tc_easy_newton(x0)
%简化newton法
%f(x)=x^3-x-1=0,取初值x0=1.5
%%the input format is :
%x0=1.5;x=tc_easy_newton(x0)
f=inline('[x^3-x-1,3*x^2-1]');
ff=f(x0);
A=ff(2);
fprintf('the result is ; ')
while 1
    if ff(1)<1e-4
        return
    end
    x1=x0-ff(1)/A;
    ff=f(x1);
    x0=x1;
    x=x1;
end