#include <cstdio>
#include <iostream>
#include <cmath>
using namespace std;
const double PI = 3.14159265358979323846264338327950288419716939937510582097494459230781641;
int main() {
double a, b, c, r;
cin >> a >> b >> c >> r;
double p = (a + b + c) / 2;
double S = sqrt(p * (p-a) * (p-b) * (p-c));
double R = S / p;
double coef = r / R;
double saved = coef * coef * (S - PI*R*R);
printf("%.15f\n", (S - saved) / S);
}
I2luY2x1ZGUgPGNzdGRpbz4KI2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y21hdGg+CnVzaW5nIG5hbWVzcGFjZSBzdGQ7Cgpjb25zdCBkb3VibGUgUEkgPSAzLjE0MTU5MjY1MzU4OTc5MzIzODQ2MjY0MzM4MzI3OTUwMjg4NDE5NzE2OTM5OTM3NTEwNTgyMDk3NDk0NDU5MjMwNzgxNjQxOwoKaW50IG1haW4oKSB7Cglkb3VibGUgYSwgYiwgYywgcjsKCWNpbiA+PiBhID4+IGIgPj4gYyA+PiByOwoJZG91YmxlIHAgPSAoYSArIGIgKyBjKSAvIDI7Cglkb3VibGUgUyA9IHNxcnQocCAqIChwLWEpICogKHAtYikgKiAocC1jKSk7Cglkb3VibGUgUiA9IFMgLyBwOwoJZG91YmxlIGNvZWYgPSByIC8gUjsKCWRvdWJsZSBzYXZlZCA9IGNvZWYgKiBjb2VmICogKFMgLSBQSSpSKlIpOwoJcHJpbnRmKCIlLjE1ZlxuIiwgKFMgLSBzYXZlZCkgLyBTKTsKfQo=