a=1;b=9;
n=40;
h=(b-a)/(2*n);
sy1=fun(a)+fun(b);
sy2=0;
sy3=0;
for j=1:2*n-1
x=a+j*h2;?
if rem(j,2)==0
sy3=sy3+fun(x);
else
sy2=sy2+fun(x);
end
S=h*(sy1+4*sy2+2*sy3)/3;
> s
s =
17.3222
因篇幅问题不能全部显示,请点此查看更多更全内容