function [y,n]=BDD(x,eps)
%該函數(shù)文件是通過(guò)不動(dòng)點(diǎn)迭代格式求解非線性方程組的近似解
%x為迭代初值，eps為允許的誤差，n記錄迭代次數(shù)
if nargin==1
    eps=1.0e-6;
elseif nargin<1
    error
return
end
x1=gg(x);   %第一次迭代
n=1;
while (norm(x1-x)>=1e-6)&(n<=10000)    %循環(huán)迭代
    x=x1;
    x1=gg(x);
    n=n+1;
end
y=x;