保捱科技网
您的当前位置:首页复化辛普森公式

复化辛普森公式

来源:保捱科技网
format long

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

end

S=h*(sy1+4*sy2+2*sy3)/3;

> s

s =

17.3222

因篇幅问题不能全部显示,请点此查看更多更全内容