自适应simpson积分模板

.
好像还挺好写的
具体思路就是一个区间，判断一下两半积分和整体差值是否过大，如果很大那么分成两边继续递归，否则直接返回
（复杂度未知，常数很大）

#include<math.h>
#include<stdio.h>
#include<string.h>
#define db double
using namespace std;
db A,B,C,D,l,r,t=1/log(10);
inline db F(db x){
	db g=log(fabs(x))*2;
	return (g+sin(A*x*x+B*x)+x)/(2*g*t+cos(C*x*x+D*x));
}
inline db f(db a,db b){  
    return (F(a)+4*F((a+b)/2.)+F(b))*(b-a)/6.;
}
inline db asr(db a,db b,db A){  
    db c=(a+b)/2.,L=f(a,c),R=f(c,b);
    if(fabs(L+R-A)<=1.5e-11) return L+R+(L+R-A)/15.;
    return asr(a,c,L)+asr(c,b,R);
}
int main(){
	scanf("%lf%lf%lf%lf%lf%lf",&A,&B,&C,&D,&l,&r);
	printf("%.5lf\n",asr(l,r,f(l,r)));
}