examp3_29.m

来自「《高等应用数学问题MATLAB求解》的PDF版及习题解答」· M 代码 · 共 10 行

x=[0:0.01:3*pi/2, 3*pi/2];  % 这样赋值能确保 $3\pi/2$ 点被包含在内
y=cos(15*x); plot(x,y)

syms x, A=int(cos(15*x),0,3*pi/2)   % 求取理论值


h0=[0.1,0.01,0.001,0.0001,0.00001,0.000001]; v=[];
for h=h0,
    x=[0:h:3*pi/2, 3*pi/2]; y=cos(15*x); I=trapz(x,y); v=[v; h,I,1/15-I];
end