newtontwo2.m

来自「数值分析方法,主要里面涵盖了计算方法方面matlab程序！」· M 代码 · 共 26 行

%Newton 
clear
clc
eps=5e-6;
delta=1e-6;
N=100;
k=0;
x0=1.0;
while(1)
    x1=x0-func2_2_1(x0)/func2_2_2(x0);
    k=k+1;
    if(k>N|abs(x1)<eps)
        disp('Newton method failed');
        break;
    end
    if abs(x1)<1
        d=x1-x0;
    else
        d=(x1-x0)/x1;
    end
    x0=x1;
    if(abs(d)<eps|abs(func2_2_1(x1))<delta)
        break;
    end
end
fprintf('%f\n',x0);