1 - \frac{x}{\left(y - z\right) \cdot \left(y - t\right)}1 - \frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{y - z} \cdot \frac{\sqrt[3]{x}}{y - t}(FPCore (x y z t) :precision binary64 (- 1.0 (/ x (* (- y z) (- y t)))))
(FPCore (x y z t) :precision binary64 (- 1.0 (* (/ (* (cbrt x) (cbrt x)) (- y z)) (/ (cbrt x) (- y t)))))
double code(double x, double y, double z, double t) {
return 1.0 - (x / ((y - z) * (y - t)));
}
double code(double x, double y, double z, double t) {
return 1.0 - (((cbrt(x) * cbrt(x)) / (y - z)) * (cbrt(x) / (y - t)));
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t
Results
Initial program 0.7
rmApplied add-cube-cbrt_binary64_59100.9
Applied times-frac_binary64_58810.7
Final simplification0.7
herbie shell --seed 2021044
(FPCore (x y z t)
:name "Data.Random.Distribution.Triangular:triangularCDF from random-fu-0.2.6.2, A"
:precision binary64
(- 1.0 (/ x (* (- y z) (- y t)))))