\left(9 \cdot {x}^{4} - {y}^{4}\right) + 2 \cdot \left(y \cdot y\right)\sqrt[3]{\left(\left(9 \cdot {x}^{4} - {y}^{4}\right) \cdot \left(9 \cdot {x}^{4} - {y}^{4}\right)\right) \cdot \left(9 \cdot {x}^{4} - {y}^{4}\right)} + 2 \cdot \left(y \cdot y\right)double f(double x, double y) {
double r7088173 = 9.0;
double r7088174 = x;
double r7088175 = 4.0;
double r7088176 = pow(r7088174, r7088175);
double r7088177 = r7088173 * r7088176;
double r7088178 = y;
double r7088179 = pow(r7088178, r7088175);
double r7088180 = r7088177 - r7088179;
double r7088181 = 2.0;
double r7088182 = r7088178 * r7088178;
double r7088183 = r7088181 * r7088182;
double r7088184 = r7088180 + r7088183;
return r7088184;
}
double f(double x, double y) {
double r7088185 = 9.0;
double r7088186 = x;
double r7088187 = 4.0;
double r7088188 = pow(r7088186, r7088187);
double r7088189 = r7088185 * r7088188;
double r7088190 = y;
double r7088191 = pow(r7088190, r7088187);
double r7088192 = r7088189 - r7088191;
double r7088193 = r7088192 * r7088192;
double r7088194 = r7088193 * r7088192;
double r7088195 = cbrt(r7088194);
double r7088196 = 2.0;
double r7088197 = r7088190 * r7088190;
double r7088198 = r7088196 * r7088197;
double r7088199 = r7088195 + r7088198;
return r7088199;
}
Results
Initial program 52.0
rmApplied add-cbrt-cube52.0
Final simplification52.0
herbie shell --seed 2019173
(FPCore (x y)
:name "From Rump in a 1983 paper"
:pre (and (== x 10864.0) (== y 18817.0))
(+ (- (* 9.0 (pow x 4.0)) (pow y 4.0)) (* 2.0 (* y y))))