1

%%%牛顿迭代法通用函数function [gen,time]=Newton(f,x0,tol)%x0为迭代初值%tol 为指定误差容限 %nargin:表示这个函数的输入变量个数%nargout:表示这个函数输出变量个数if (nargin==2)  %nargin:表示这个函数的输入变量个数,相当于定义默认误差容限    tol=1.0e-5;end%计算原函数的导数   df=diff(sym(f));   x1=x0;%time：记迭代次数   time=0;   wucha=0.1;%给定一个误差初值 以方便进入循环计算      while (wucha>tol)       time=time+1;       fx=subs(f,x1);        df=subs(df,x1);       gen=x1-fx/df;  %牛顿迭代法的核心公式       gen=vpa(gen);       wucha=abs(gen-x1);       x1=gen;       %x_R(time)=gen;   endend