newton.m

来自「用牛顿法解方程组」· M 代码 · 共 20 行

function X=Newton(X0,tol,nmax,varargin)%输入初值X=X0+1，误差tol，和最大迭代次数，X0为行向量格式。
X=X0'+1;%                               可以编辑F.m和J.m文件来输入方程组和对应的雅克比矩阵。
niter=0;
while ((X(1)-X0(1)) >= tol || (X(2)-X0(2)) >= tol || (X(3)-X0(3)) >= tol ||....
        (X(4)-X0(4)) >= tol) && niter <= (nmax-1)
    niter=niter+1;
    X1=X-inv(J(X))*F(X);
    X0=X;
    X=X1;
end
niter
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Newton([0 0 0 0],1.e-3,100)
%niter = 20
%ans =
%    2.9937
%    0.0317
%   -0.3240
%    0.2607