double f(double a, double b, double c) {
double r4004052 = a;
double r4004053 = b;
double r4004054 = r4004052 + r4004053;
double r4004055 = c;
double r4004056 = r4004054 + r4004055;
double r4004057 = 2.0;
double r4004058 = r4004056 / r4004057;
double r4004059 = r4004058 - r4004052;
double r4004060 = r4004058 * r4004059;
double r4004061 = r4004058 - r4004053;
double r4004062 = r4004060 * r4004061;
double r4004063 = r4004058 - r4004055;
double r4004064 = r4004062 * r4004063;
double r4004065 = sqrt(r4004064);
return r4004065;
}