#include <iostream>
#include <cmath>

using namespace std;

const double EPSILON = 1e-9;

struct Vector2 {
  double x, y;

  explicit Vector2(double _x =0, double _y =0): x(_x), y(_y) {}

  Vector2 operator -(Vector2 r) {
    return Vector2(x-r.x, y-r.y);
  }

  double norm() {
    return hypot(x,y);
  }

  double dot(Vector2 r) {
    return x*r.x + y*r.y;
  }
  double cross(Vector2 r) {
    return x*r.y-y*r.x;
  }

};

double CircleArea(double r) {
  return r*r*M_PI;
}

double Theta(double y, double x) {
  return fmod(atan2(y,x)+2*M_PI, 2*M_PI);
}
double Theta2(Vector2 a, Vector2 b) {
  return acos(a.dot(b)/(a.norm()*b.norm()));
}

int main() {
  cout.precision(20);
  ios::sync_with_stdio(false);
  int cases;
  cin >>cases;
  while(cases--) {
    double d, r, v_t, t, v_h;
    cin >> d >> r>> v_t >> t>> v_h;
    double ans =0.0;
    if(fabs(v_h) < EPSILON) {
      // ans = 0;
    } else if(fabs(v_t) < EPSILON) {
      ans = CircleArea(r)/4.0;
    } else {
      double r1 = (d/v_h -t)*v_t;
      if(r+r1 <=d) {
        ans = CircleArea(r)/4.0;
      } else if(r >= d+r1) {
        ans = CircleArea(r)/4.0 - CircleArea(r1)/2.0;
      } else if(r1 >= d+r) {
        //ans = 0.0
      } else {
        double y = (r*r+d*d-r1*r1)/(2*d);
        double x = sqrt(r*r - y*y);
        if(y <=0) {
          //ans = 0.0;
        } else {
          ans = CircleArea(r)/4 - (CircleArea(r)/(2*M_PI)*Theta2(Vector2(x,y), Vector2(0,d)) - d*x/2 + CircleArea(r1)/(2*M_PI)*Theta2(Vector2(0,r1-d), Vector2(x,y-d)));
          if(r1>d) {
            double x2 =sqrt(r1*r1-d*d);
            ans = ans - (CircleArea(r1)/ (2*M_PI)*Theta2(Vector2(0, -d), Vector2(x2 ,-d)) - d*x2/2);
          }
        }
      }
    }
    if(fabs(ans) < EPSILON) cout << "0.0000000000"<<endl;
    else cout << ans << endl;
  }
}