数值分析(第五版)计算实习题第四章作业.pdfVIP

第四章： 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.4170 输入： syms x f=inline(sqrt(x).*log(