函数文件: function x=newton_Iterative_method(f,n,Initial) x0=Initial; tol=1e-11; x1=x0-Jacobian(f,n,x0)\F(f,x0); while (norm(x1-x0,2)>tol) %数值解的2范数是否在误差范围内 x0=x1; x1=x0-Jacobian(f,n,x0)\F(f,x0); end x=x1;%不动点 function g=Jacobian(f,n,a) %求解任意矩阵的雅可比矩阵 %% s