newton.m

来自「Numerical analysis （Gauss、Bisection、jaco」· M 代码 · 共 29 行

function s=Newton(fx,dfx,x0,eps)
if nargin==3
    eps=1.0e-6;
elseif nargin<3
    return
end
x=x0;
f0=eval(fx);
f1=eval(dfx);
if (f0~=0 && f1~=0)
    x1=1;
    i=0;
    while abs(x1-x0)>eps    
%        x1=x0-f0/f1;
%        x=x1;
%        f2=eval(fx);
%        f3=eval(dfx);
        x1=x0;
        x0=x0-f0/f1;
        x=x0;
        f0=eval(fx);
        f1=eval(dfx);
        i=i+1;
        x=[i,x1,x0,f0];
        disp(x);
    end
end