e0319.asv

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

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