ex3-13-1_2.m

matlab程序

% 用梯形法求数值积分

d=pi/1000; t=0:d:3*pi;  nt=length(t);

Y=exp(-0.5*t).*sin(t+pi/6);

Z=trapz(Y)*d    % 用梯形法求数值积分