MATLAB：牛顿算法 - 雅可比矩阵的双边逼近

Question

我对MATLAB中的牛顿算法有疑问。

该算法应该能够执行双侧近似 在没有提供分析雅可比矩阵的情况下雅可比矩阵的技术。

function [x,fx,ef] = newton(f,x,cc)

% convergence criteria
tole = cc(1,1); told = cc(2,1); maxiter = cc(3,1);

% newton algorithm
ef = 0;
for j = 1:maxiter
   [fx,dfx] = f(x);
   xp = x − dfx\fx;
   D = (norm(x−xp)<=tole*(1+norm(xp)) && norm(fx)<=told);
   if D == 1;
      ef = 1; break;
   else
      x = xp;
   end
end