自适应辛普森积分

问题：求${\int }_a^{b}f(x)dx$
辛普森积分的本质思想是将$f(x)$看作一个二次函数来简化问题
令${\int }_a^{b}f(x)dx\approx {\int }_a^b(A{x}^2+Bx+C)dx$
即：$f(x)\approx A{x}^2+Bx+C$
可以解得：$${\int }_a^{b}f(x)dx\approx \frac{(b-a)}{6}\ast [f(a)+4\ast f(\frac{a+b}{2})+f(b)]$$

因为$|b-a|$越小,误差越小，所以可以递归处理该问题
因为：${\int }_a^{b}f(x)dx={\int }_a^{c}f(x)dx+{\int }_c^{b}f(x)dx$

至于精度问题，由误差分析的知识可得：($\epsilon$为精度误差)
令$c=\frac{a+b}{2}$，$$|{\int }_a^{c}f(x)dx+{\int }_c^{b}f(x)dx-{\int }_a^{b}f(x)dx|<=15\epsilon$$
时即可返回;
注意： 每次递归的时候由于误差的相加性，误差设置必须变为原来的一半

$Code$

inline double f(double x){
	//...
}

inline double simpson(double a, double b){
	return (f(a) + f(b) + 4 * f((a + b) / 2) ) * (b - a) / 6;
}
double asr(double a, double b, double eps, double A){ //A为a,b区间的积分
	double c = (a + b) / 2;
	double l = simpson(a, c), r = simpson(c, b);
	if (fabs(l + r - A) <= 15 * eps)	return l + r + (l + r - A) / 15.0;
	return asr(a, c, eps / 2.0, l) + asr(c, b, eps / 2.0, r);
}