function?Sn=CompoundSimpson（fabn）
h=（b-a）/n;
for?k=0:n???
????x（k+1）=a+k*h;?????
????xk（k+1）=x（k+1）+1/2*h;???
????if?（x（k+1）==0）????
????????x（k+1）=10^（-10）;???????
????end
????if?（xk（k+1）==0）????????????
????????xk（k+1）=10^（-10）;
????end
end
p1=h/6*（f（x（1））+f（x（n+1）））;
for?i=2:n?????
????fx1（i）=h/3*f（x（i））;
end
for?j=1:n?????
????fx2（j）=2/3*h*f（xk（j））;
end
p2=sum（fx1）;
p3=sum（fx2）;
Sn=p1+p2+p3;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2019-03-18?00:43??5\
?????文件?????????271??2018-07-05?20:23??5\CompoundTrapezoid.m
?????文件?????????422??2018-07-05?20:37??5\CompoundSimpson.m