newton.m

来自「用Jacobi叠待法解线性方程组 function Jacobi(A,b,n,」· M 代码 · 共 17 行

function Newton(fx,x0,E,N)
%digits(10)
x=0
for n=1:N
    x=x0-subs(fx,x0)/subs(diff(fx),x0);
    fprintf('%f\n',x);
    if abs(det(x-x0))<E
        break;
    end
    if n<N
        n=n+1;
        x0=x;
    else
        break;
    end
end
end