P4525 【模板】自适应辛普森法 1

代码👇：

#include<bits/stdc++.h>
using namespace std;
double a,b,c,d,l,r;
double f1(double x)
{
	return (c*x)/a+((a*d-b*c)/(a*a))*log(fabs(a*x+b));
}
double f2(double x)
{
	return (c/(2*b))*x*x+(d/b)*x;
}
int main() {
	cin>>a>>b>>c>>d>>l>>r;
	if(a==0.0)
	{
		printf("%.6f",f2(r)-f2(l));
	}
	else
	{
		printf("%.6f",f1(r)-f1(l));
	}
	return 0;
}