#include #include #include using namespace std; int main(){ double R, H, D; cin >> R >> H >> D; D *= 0.5; auto area = [&](double d){ return 0.5*H/R*(d*d*(log(d)-log(R+sqrt(R*R-d*d))) + R*sqrt(R*R-d*d)); }; double a = 0.0; double pre = area(D); double preD = D; for(int i=1;i<=1000000;i++){ double r = 0.000001 * i; double d = (1-r)*D + r * R; double a1 = area(0.5*(preD+d)); double a2 = area(d); a += (d-preD)/6*(pre+4*a1+a2); preD = d; pre = a2; } printf("%.8lf\n", 4*a); // double res = 1.0/6*H*R*(D*D*D*log(D) + 2*D*R*sqrt(R*R-D*D) + R*R*R*atan(D/sqrt(R*R-D*D)) - D*D*D*log(sqrt(R*R-D*D)+R)); // printf("%.8lf\n", -res); }