ex3_21.m

来自「这里有很多的matlab程序设计的源代码」· M 代码 · 共 12 行

a=0;b=3*pi;
n=1000; h=(b-a)/n;
x=a; s=0; 
f0=exp(-0.5*x)*sin(x+pi/6);
for i=1:n
    x=x+h;
    f1=exp(-0.5*x)*sin(x+pi/6);
    s=s+(f0+f1)*h/2;
    f0=f1;
end
s