求一段区间函数与x轴的面积问题
const double eps = 1e-7;
double a,b,c,d,L,R;
double f(double x) {
return (c*x+d)/(a*x+b);
}
double simpson(double l,double r) {
return (f(l)+f(r)+4*f((l+r)/2))*(r-l)/6;
}
double asr(double l,double r,double ep,double val) {
double mid = (l+r)/2;
double lval = simpson(l,mid),rval=simpson(mid,r);
if(fabs(lval+rval-val)<=15*ep) {
return lval+rval+(lval+rval-val)/15;
}
return asr(l,mid,ep/2,lval)+asr(mid,r,ep/2,rval);
}
double asme(double l,double r,double ep) {
return asr(l,r,ep,simpson(l,r));
}
void work() {
scanf("%lf%lf%lf%lf%lf%lf",&a,&b,&c,&d,&L,&R);
printf("%lf\n",asme(L,R,eps));
}
/*
const double inf = 30,zero = 1e-9
double f(double x) {
return pow(x,(a/x)-x);
}
double ans = asme(zero,inf,eps);
*/