元青花龙纹小盘:复合simpson求积公式设计(C语言)

Nweton—Cotes公式的求积余项表明，求积节点n越大，对应的求积公式精度越高，但由于Nweton—Cotes公式在n>8时数值不稳定，因此不能用增加求积节点数的方法来提高计算精度。实用中常将求积区间[a，b]分成若干个小区间，然后在每个小区间上采用数值稳定的Nweton—Cotes公式求小区间上的定积分，最后把所有小区间上的计算结果相加来作为原定积分的近似值。采用这种方法构造的求积公式就称为复合求积公式。复合求积公式具有计算简单且可以任意逼近所求定积分值的特点，这是Nweton—Cotes公式一般做不到的。常用的复合求积公式有复合梯形公式和复合Simpson公式。

① 复合梯形公式

取等距节点 xk=a+kh ,h=(b-a)/n ,k=0,1,....,n 将积分区间[a,b]n等分，在每个小区间[xk，xk+1] k=0,1,...n-1上用梯形公式做近似计算，就有

得求积公式

（15）是称为复合梯形公式。

通常记（15）的右端为

它称为 T 形值。因为

故复合梯形公式的求积余项为

如果|f"(x)|≤M2，由于a,b是有限数，故若给定计算精度ξ，由（16），令

即只要取h满足（17），及n=(b-a)/h，利用复合求积公式（15）计算，就能得到计算误差小于ξ的定积分近似值。

②复合Simpson公式

取[a,b]上的等距节点 xk=a+kh ,h=(b-a)/n ,k=0,1,...n，将[a,b]n等分，在每个小区间[xk，xk+1]上用Simpson公式做近似计算，就有

式中xk+1/2 = xk+h/2，于是求得下面复合Simpson公式：

利用Simpson公式的求积余项，可以得到复合Simpson公式的求积余项为：

从复合Simpson公式的求积余项可以看出复合Simpson公式比复合梯形公式好，不过前者计算量大些。复合Simpson公式也称为复合抛物线公式。

采用如上构造复合求积公式的方法可以构造出其他类型的复合求积公式，如复合Cotes公式等。实际上复合求积公式本质上是用在求积节点上的分段插值函数代替被积函数f(x)获得求积公式，所以它们归为插值型求积公式，如复合梯形公式是用分段型插值函数取代f(x)参与积分，复合Simpson公式是用分段二次插值函数取代f(x)参与积分等。

数学的实践与认识
MATHEMATICS IN PRACTICE AND THEORY
2005 Vol.35 No.6 P.245-247

--------------------------------------------------------------------------------

经典SIMPSON求积公式的新证明

A New Proof of the Classical Simpson Quadrature Formula

石艳霞 刘证

摘 要：对于定积分近似计算中常使用的经典SIMPSON求积公式介绍一种新的简洁的证明方法并给出误差的最佳估计.
关键词：定积分;分部积分法;SIMPSON求积公式;误差

作者单位：石艳霞（鞍山科技大学理学院,辽宁,鞍山,114044）
刘证（鞍山科技大学理学院,辽宁,鞍山,114044）

参考文献：

[1]菲赫金哥尔茨Γ М著,北京大学高等数学教研组译.微积分学教程:第二卷,第一分册[M].上海:商务印书馆,1955.
[2]胡祖炽.计算方法[M].北京:高等教育出版社,1959.
[3]林成森.数值计算方法[M].北京:科学出版社,2001.
[4]刘玉琏,傅沛仁.数学分析讲义(第三版):上册[M].北京:高等教育出版社,1992.
[5]张筑生.数学分析新讲:第二册[M].北京:北京大学出版社,1999.
[6]Pearce C E M, Pecari c J, Ujevi c N, Varosanec S. Generalizations of some inequalities of Ostrowski-Gruss type[J]. Mathematical Inequalities & Applications, 2000, 3(1): 25-34.

都是扯一堆废话，没有给源代码

拿公式迭代就可以