计算复杂函数的定积分

设，求。

解：求函数f(x)在[a,b]上的定积分，其几何意义就是求曲线y=f(x)与直线x=a,x=b,y=0所围成的曲边梯形的面积。为了求得曲边梯形的面积，先将积分区间[a,b]分成n等份，每个区间的宽度为h=(b-a)/n，对应的将曲边梯形分成n等份，每个小部分即是一个小曲边梯形。近似求出每个小曲边梯形的面积。

a=0;
b=3*pi;
n=1000;
h=(b-a)/n;
x=a;
s=0;
f=exp(-0.5*x)*sin(x+pi/6);
for i=1:n
    x=x+h;
    f1=exp(-0.5*x)*sin(x+pi/6);
    s=s+(f+f1)*h/2;
    f=f1;
end
s