e0319.m

来自「Matlab7.x程序设计-王建卫-源代码及习题答案」· M 代码 · 共 11 行

a=0;b=2*pi;n=1000;h=(b-a)/n;
x=a:h:b;y=0;
f=cos(x-pi/6).*sin(x+pi/6);
for i=1:n
    s(i)=(f(i)+f(i+1))*h/2; 
    y=y+s(i); 
end
y