h=pi/20;
x=0:h:pi/2;
y=sin(x);
z1=sum(y(1:10))*h,  %矩形公式1（內部）
z2=sum(y(2:11))*h,  %矩形公式2（外部）
z=cumsum(y);
z11=z(10)*h,        %矩形公式1（內部）
z12=(z(11)-z(1))*h, %矩形公式2（外部）
z3=trapz(y)*h,      %梯形公式