#include<bits/stdc++.h>
#define double long double
#define pi 3.14159265358979323846
using namespace std;
struct point
{
	double x,y;
	point rot(double t)const{return {x*cos(t)-y*sin(t),x*sin(t)+y*cos(t)};}
	double len2()const{return x*x+y*y;}
	double len()const{return sqrt(len2());}
	point operator*(const double t)const{return {t*x,t*y};}
	point operator/(const double t)const{return {x/t,y/t};}
	point operator+(const point o)const{return {x+o.x,y+o.y};}
	point operator-(const point o)const{return {x-o.x,y-o.y};}
};
double det(point i,point j){return i.x*j.y-j.x*i.y;}
double dot(point i,point j){return i.x*j.x+i.y*j.y;}
struct line
{
	point s,t;
};
point line_int(line a,line b)
{
	double s1=det(a.t-a.s,b.s-a.s);
	double s2=det(a.t-a.s,b.t-a.s);
	return (b.s*s2-b.t*s1)/(s2-s1);
}
point project_to_line(point a,line b)
{
	return b.s+(b.t-b.s)*(dot(a-b.s,b.t-b.s)/(b.t-b.s).len2());
}
int main()
{
	int t;
	cin>>t;
	while(t--)
	{
		double la,lb,a,b;
		cin>>la>>lb>>a>>b;
		a*=pi/180,b*=pi/180;
		if(b<=pi/2)
		{
			printf("%.10Lf\n",(la+lb)*(la+lb)*a+lb*lb*b);
			continue;	
		}
		double ans=(la+lb)*(la+lb)*a+lb*lb*b;
		point O={0,0},A={la,0};
		A=A.rot(a);
		point A_={A.x,-A.y};
		point P={lb,0};
		P=P.rot(a+b)+A;
		point Q=project_to_line(O,{P,A});
		if(Q.x<P.x)Q=P;
		point Q_=Q.rot(-2*a);
		if(det(A,Q_)>=0)
		{
			ans+=det(A,Q)-Q.len2()*(acos(Q.x/Q.len())-acos(Q_.x/Q_.len()));
			ans-=det(line_int({O,A},{Q_,A_}),Q_);
			printf("%.10Lf\n",ans);
		}
		else 
		{
			ans+=det(A,Q)-Q.len2()*(acos(Q.x/Q.len())-a);
			printf("%.10Lf\n",ans);
		}
	}
	return 0;
}