double f(double a, double b, double c) {
double r1551562 = a;
double r1551563 = b;
double r1551564 = r1551562 + r1551563;
double r1551565 = c;
double r1551566 = r1551564 + r1551565;
double r1551567 = 2.0;
double r1551568 = r1551566 / r1551567;
double r1551569 = r1551568 - r1551562;
double r1551570 = r1551568 * r1551569;
double r1551571 = r1551568 - r1551563;
double r1551572 = r1551570 * r1551571;
double r1551573 = r1551568 - r1551565;
double r1551574 = r1551572 * r1551573;
double r1551575 = sqrt(r1551574);
return r1551575;
}