main.m

来自「应用matlab应用平台信号采集程序设计与整理。」· M 代码 · 共 18 行

f=inline('(-exp(-x*t)).*sin(2.*(x.^0.5)).*x./(1+x.^2)');
b=1;
a=0;
M=10;
h=(b-a)/(2*M);
s1=0;
s2=0;
t=0:0.1:10;
for K=1:M
    x=a+h*(2*K-1);
    s1=s1+f(x,t);
end
for K=1:(M-1)
    x=a+h*2*K;
    s2=s2+f(x,t);
end
s=h*(f(a,t)+f(b,t)+4*s1+2*s2)/3;
<i