题目描述
题解
辛普森积分裸题
答案要乘2吖,因为椭圆是两半
代码
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
const double eps=1e-10;
int dcmp(double x)
{
if (x<=eps&&x>=-eps) return 0;
return (x>0)?1:-1;
}
int T;
double a,b,l,r,ans;
double F(double x)
{
return b*sqrt(1-((x*x)/(a*a)));
}
double calc(double l,double r)
{
return (F(l)+F(r)+F((l+r)/2)*4)*(r-l)/6;
}
double simpson(double l,double r,double now)
{
double mid=(l+r)/2;
double p=calc(l,mid),q=calc(mid,r);
if (dcmp(now-p-q)==0) return now;
else return simpson(l,mid,p)+simpson(mid,r,q);
}
int main()
{
scanf("%d",&T);
while (T--)
{
scanf("%lf%lf%lf%lf",&a,&b,&l,&r);
ans=simpson(l,r,calc(l,r))*2;
printf("%.3lf\n",ans);
}
}