ex3-13-1_2.m

来自「matlab程序」· M 代码 · 共 8 行

% 用梯形法求数值积分

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

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

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