计算方法书中的MATLAB实例.pptVIP

用Matlab实现求积分 在Matlab程序编辑器中输入： y0=log(6.0/5.0); fprintf(y[%d]=%f\n,0,y0) n=1; while(1) yl=1.0/n-5*y0; fprintf(y[%d]=%f\n,n,yl) if(n=20) break; end y0=yl; n=n+1; end function x=nabisect(fname,a,b,e) %fname为内嵌函数 表达式；a,b为区间端点；e为输入定义的精度 if nargin4,e=1e-4;end; % nargin为输入变量的个数， 若未输入精度变量，则给定默认精度 fa=feval(fname,a);fb=feval(fname,b); % 计算端点的 函数值 if fa*fb0,error(函数在两端点值必须异号);end x=(a+b)/2 while (b-a)(2*e) fx=feval(fname,x); if fa*fx0,b=x;fb=fx; else a=x;fa=fx;end