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
因篇幅问题不能全部显示,请点此查看更多更全内容
Copyright © 2019- baoaiwan.cn 版权所有 赣ICP备2024042794号-3
违法及侵权请联系:TEL:199 18 7713 E-MAIL:2724546146@qq.com
本站由北京市万商天勤律师事务所王兴未律师提供法律服务