求面积的数值积分程序exn541 for dx=[2,1,0.5,0.1] %设不同步长 X=0:.1:10;y=-X.*x+115; %取较密的函数样本 plot(x,y),hold on %画出被积曲线并保持 x1=0:dx10;y1=-X1.*x1+115; %求取样点上的y1 %用矩形（欧拉）法求积分，注意末尾去掉一个点 n=length(x1);s=sum(y1(1:n-1))*dx; q=trapz(y1)*dx; %用梯形法求积分 stairs(x1,y1), %画出欧拉法的积分区域 plot(x1,y1) %画出梯形法的积分区域 [dx,s,q],pause(1) hold off,end 求面积的数值积分程序exn541 for dx=[2,1,0.5,0.1] % 设不同步长 x=0:.1:10;y=-x.*x+115; % 取较密的函数样本 plot(x,y),hold on % 画出被积曲线并保持 x1=0:dx:10;y1=-x1.*x1+115; % 求取样点上的y1 % 用矩形（欧拉）法求积分，注意末尾去掉一个点 n=length(x1);s=sum(y1(1:n-1))*dx; q=trapz(y1)*dx; % 用梯形法求积分 stairs(x1,y1), % 画出欧拉法的积分区域 plot(x1,y1) % 画出梯形法的积分区域 [dx,s,q],pause(1), hold off，end