dx=0.001;%采样间隔 x=0:dx:1;%数值计算适合于有限区间上,取有限个采样点 Ft=sqrt(2*x.^2-x+5); Sx=dx*cumtrapz(Ft);%计算区间内曲线下图形面积,为小矩形面积累加得 disp(['积分值=' num2str(Sx(end))])%所求定积分值 %图示 plot(x,Ft,'r','MarkerSize',4)%被积分函数图 hold on plot(x,Sx,'.k','MarkerSize',15)%积分图 hold off xlabel('x') legend('被积分函数图','积分图')