examp3_29.m

东北大学 薛定谔的高等应用数学问题MATLAB求解的源代码

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