文档介绍：Simpson-Rule-Summary---辛普森法则
LT
考虑积分，如果在区间[a，b]内取等间隔的N份，间隔长度为h，简述矩形（、梯形、Simpson法则）计算积分的i）理论、误差精度分析，和算法计算流程。
解：对于缓变函数我们可以用各个区间中点上函数值作来近似该区间的平均值
其中。
矩形法则：f(x)在区间[a,b]上的积分用矩形求积定义如下
第i个区间对积分的贡献为：
如果围绕该区间中点的邻域内对函数f(x)作泰勒级数展开， 有
其中，和分别表示了f(x)在处的一阶，二阶和三阶导数。相应地，积分在子区间内的值可以表示为
其中第一项是矩形积分的近似值，第二项则由于其中的积分等于零而消除。从而，矩形法则在宽度为h的单个子区间内的最高阶误差由第三项给出
在整个[a,b]区间上的总误差则通过将所有N个子区间的贡献相加得到
其中我们利用了Nh=(b-a)，并且取为f(x)在[a,b]上的二阶导数的均值。
梯形法则：
如果用一个梯形来近似代替每个子区间的面积，梯形的四个顶点分别位于，，与。梯形的面积为
（2-8）
对于整个[a,b]区间，积分值由所有窄带的和给出
Simpson法则
当积分的上下限相对于展开的中心点对称时，积分泰勒展开式中含有f(x)的奇数阶导数的项都将等于零。利用这一性质，我们可以在相邻的两个子区间内对面积作泰勒级数展开，可得
积分的精确度可以达到。利用中心差分，可将f(x)的二阶导数近似表为
将此结果代入前述中心查分表式，则子区间[xi-1,xi+1]内的积分近似为
对于整个[a,b]区间，结果为
假设区间的数目N为偶数，则奇数格点的贡献是偶数格点贡献的2倍。这种权重上的差异来自我们为修正基本方法所得的一级结果所引入的f(x)二阶导数的贡献。
而两端点和的权重则仅为偶数点的一半。Simpson法则的另一优点在于其自然引出了一种算法，即通过迭代使积分达到所需要的精确度。
误差估计---原则上，我们可以通过估计所省略的高阶项大小得出数值计算中的误差。其前提是知道高阶导数（如）的平均值），这样做并不容易。比较相继的两次迭代结果是较为简便的一种方法。例如，如果两次计算值I[a,b]之差小于预设的允许值，则我们认为计算已经达到了所需的精确度。对于Simpson法则，该算法可通过如下方法实施。假设我们在x=a和x=b之间对函数f(x)进行积分，且要求达到精确度。各个节点上的函数值对积分的贡献可分为三个部分，即端点区，奇数点与偶数点的贡献，
其中，。可用，及以如下形式写出积分
由此可见，对于一个新的N，我们可利用上一次循环中已有的
，及的值。为便于进行迭代计算，应从N=开始，其中应为某个较小而合理的区间数目，例如取=6。这是迭代计算的第一步。在下一轮迭代中，将区间的数目翻倍，而令步长h减半，即令N==2。以这种方式改变子区间数目的好处在于，对于一个新的N，来自两个端点的贡献（存储为）是不变的，而偶数点与上一次迭代计算过的偶数点完全相同。也就是说，在新的迭代中的的值就是上一次循环的与之和。所需计算的仅仅是函数在（新的）奇数点上的值。从而我们得到结果
其中。利用上面等式进行变步长的循环计算即。如果计算值与前一循环所得结果之差小于ε，则可认为计算结果已收敛；否则即再次令子区间数目加倍，使N=N3=2N2。然后依此类推，直到计算结果收敛或者达到最大迭代次数。框图给出了算法的主要步骤。
Box 2-2 Program DM_SIMPS
Simpson’s rule integration
Initialization:
Set ε as error tolerance and N as the initial number of points.
Set a maximum for the number of iterations allowed.
Input a as the lower limit of integration and b as the upper limit.
Make an initial calculation of the integral with N points:
Calculate the contributions from end points.
Sum over the contributions from odd points.
Sum over the contributions from even points.
Ca