#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define MP make_pair
mt19937 rnd(time(0));
const long double pi=acos(-1);
long double solve(){
	long double x,y,a,b;cin>>x>>y>>a>>b;
	a=a/180*pi,b=b/180*pi;
	long double ans=y*y*b+(x+y)*(x+y)*a;
	if(b>pi/2){
		long double r=sqrt(x*x+y*y-2*cos(pi-b)*x*y);
		long double c=acos((r*r+x*x-y*y)/(2*r*x)),d=acos((r*r+y*y-x*x)/(2*r*y)); 
//		cerr<<r<<' '<<c<<endl;
		if(r*r+y*y<x*x){
			ans+=2*(0.5*sin(pi-b)*x*y-c/(2*pi)*r*r*pi);
			if(c>2*a){
				auto f=[&](long double l,long double a,long double b){
					long double n=l/sin(a+b)*sin(b);
					return 0.5*n*l*sin(a)-l*l*pi*(b/(2*pi));
				};
				ans-=2*f(r,d,c-2*a);
			}
		}else{
			long double h=sin(pi-b)*x;
			ans+=2*(0.5*sin(pi-b)*cos(pi-b)*x*x-(b-pi/2)/(2*pi)*h*h*pi);
			if(b-pi/2>2*a){
				auto f=[&](long double r,long double a){
//					cout<<r<<' '<<a<<endl;
					return 0.5*r*r*tan(a)-a/(2*pi)*r*r*pi;
				};
				ans-=2*f(h,b-pi/2-2*a);
			}
		}
	} 
	return ans;
}
int main(){
	ios::sync_with_stdio(false);
	int _;cin>>_;
	while(_--) cout<<setprecision(12)<<solve()<<endl;
}