broy.m

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

function [y,n]=Broy(x0,eps)
%该函数文件是通过Broyden迭代格式求解非线性方程组的近似解
%x为迭代初值，eps为允许的误差，n记录选代次数
if nargin==1
    eps=1.0e-6;
elseif nargin<1
    error
    return
end
A=eye(length(x0));
x1=x0-bf(x0)/A;
n=1;
while (norm(x1-x0)>=1e-6)&(n<=10000)
    x0=x1;
    x1=x0-bf(x0)/A;
    p=x1-x0;
    q=bf(x1)-bf(x0);
    A=A+(q-p*A)'*p/norm(p);
    n=n+1;
end
y=x1;