数值分析（第五版）计算实习题第四章作业(共6页).doc

精选优质文档-倾情为你奉上第四章：1、（1）：复合梯形建立m文件：function t=natrapz(fname,a,b,n) h=(b-a)/n; fa=feval(fname,a);fb=feval(fname,b);f=feval(fname,a+h:h:b-h+0.001*h); t=h*(0.5*(fa+fb)+sum(f); 输入： syms x f=inline(sqrt(x).*log(x);); natrapz(f,eps,1,10)输出：ans = -0.9470输入： syms x f=inline(sqrt(x).*log(x);); natrapz(f,eps,1,100)输出：ans = -0.8157输入： syms x f=inline(sqrt(x).*log(x);); natrapz(f,eps,1,1000)输出：ans = -0.7162复合辛普森