\left(9 \cdot {x}^{4} - {y}^{4}\right) + 2 \cdot \left(y \cdot y\right)\sqrt[3]{{\left(9 \cdot {x}^{4} - {y}^{4}\right)}^{3}} + 2 \cdot \left(y \cdot y\right)double f(double x, double y) {
double r58116 = 9.0;
double r58117 = x;
double r58118 = 4.0;
double r58119 = pow(r58117, r58118);
double r58120 = r58116 * r58119;
double r58121 = y;
double r58122 = pow(r58121, r58118);
double r58123 = r58120 - r58122;
double r58124 = 2.0;
double r58125 = r58121 * r58121;
double r58126 = r58124 * r58125;
double r58127 = r58123 + r58126;
return r58127;
}
double f(double x, double y) {
double r58128 = 9.0;
double r58129 = x;
double r58130 = 4.0;
double r58131 = pow(r58129, r58130);
double r58132 = r58128 * r58131;
double r58133 = y;
double r58134 = pow(r58133, r58130);
double r58135 = r58132 - r58134;
double r58136 = 3.0;
double r58137 = pow(r58135, r58136);
double r58138 = cbrt(r58137);
double r58139 = 2.0;
double r58140 = r58133 * r58133;
double r58141 = r58139 * r58140;
double r58142 = r58138 + r58141;
return r58142;
}
Results
Initial program 52.0
rmApplied add-cbrt-cube52.0
Simplified52.0
Final simplification52.0
herbie shell --seed 2019350
(FPCore (x y)
:name "From Rump in a 1983 paper"
:precision binary64
:pre (and (== x 10864) (== y 18817))
(+ (- (* 9 (pow x 4)) (pow y 4)) (* 2 (* y y))))