自适应辛普森积分算法

辛普森积分是数值积分的一种，是中点公式和梯形公式的改进。

——这是我在APIO2017时听某位大神讲课所了解的。

下面就讲讲辛普森积分到底是啥玩意儿。

假定我们要求如下定积分

$$\int_a^bf(x)dx$$
略微懂一点微积分知识的都知道，对于一个黎曼可积的函数，我们要求其在某个闭区间上的定积分，要先求该函数的不定积分，即先求原函数。就是找到一个函数 $F(x)$，使得 $F'(x)=f(x)$，然后根据牛顿—莱布尼茨公式去搞。即 $$\int_a^bf(x)dx=F(b)-F(a)$$
而有时候原函数并不好求，比如要求原的函数长得很复杂，要求的原函数非初等函数。。
这时我们需要“避其锋芒”，采用一种利用“近似”思想的方法。

我们都知道定积分就是曲线$f(x)$下方的面积(假设$f(x)$恒大于0)，我们考虑一种近似求出曲线下方面积的方法——辛普森法。

我们先在区间$[a,b]$中等距地取奇数个点 a=