计算方法书中的MATLAB实例2程序运行结果:functionx=nabisect(fname,a,b,e)%fname为内嵌函数表达式；a,b为区间端点；e为输入定义的精度ifnargin<4,e=1e-4;end;%nargin为输入变量的个数，若未输入精度变量，则给定默认精度fa=feval(fname,a);fb=feval(fname,b);%计算端点的函数值iffa*fb>0,error('函数在两端点值必须异号');endx=(a+b)/2while(b-a)>(2*e)fx=feval(fname,x);iffa*fx<0,b=x;fb=fx;elsea=x;fa=fx;endx=(a+b)/2end5程序运行结果：用Matlab实现Newton迭代法求解8大家有疑问的，可以询问和交流程序运行结果：用Matlab实现顺序Gauss消去法%回代x=zeros(n,1);x(n)=a(n,n+1)/a(n,n);fork=n-1:-1:1x(k,:)=(a(k,n+1)-a(k,(k+1):n)*x((k+1):n))/a(k,k);end13程序运行结果：用Matlab实现选列主元Gauss消去法解线性方程组ifflag==0,a,endend%回代x=zeros(n,1);x(n)=a(n,n+1)/a(n,n);fork=n-1:-1:1x(k,:)=(a(k,n+1)-a(k,(k+1):n)*x((k+1):n))/a(k,k);end17程序运行结果:用Matlab实现LU分解2023求解例题得一次插值和二次插值如右图所示：26用Matlab实现用Euler格式解常微分方程29Euler格式求解得：用Matlab实现用改进Euler格式解常微分方程32