function [r,n]=mulNewton_H(x0,eps)
if nargin==1
    eps=1.0e-6;
end
Fx = myfun_H(x0);
dF = detamyfun_H(x0);
r=x0-inv(dF)*Fx;
n=1;
tol=1;
while tol>eps
    x0=r;
    Fx = myfun_H(x0);
    dF = detamyfun_H(x0);
    r=x0-inv(dF)*Fx;                                 %核心迭代公式
    tol=norm(r-x0);
    n=n+1;
    if(n>100000)                                              %迭代步數(shù)控制
        disp('迭代步數(shù)太多，可能不收斂！');
        return;
    end
end