newdd.m

来自「有关matlab r2007的实用教程」· M 代码 · 共 25 行

function [y,n]=NewDD(x,eps)
%该函数文件是通过牛顿迭代格式求解非线性方程组的近似解
%x为迭代初值，eps为允许的误差，n记录迭代次数
if nargin==1
    eps=1.0e-6;
elseif nargin<1
    error
    return
end
x1=fx1(x);
x2=-dfx1(x);   %线性方程组导数
x3=inv(x2);
x0=x3*x1';
n=1;
while norm(x0)>=eps   %循环迭代
    x=x0'+x;
    x1=fx1(x);
    x2=-dfx1(x);
    x3=inv(x2);
    x0=x3*x1';
n=n+1;
end
y=x0'+x;