程序如下:
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('被积分函数图','积分图')
结果:
积分值=2.2719