1 - \frac{x}{\left(y - z\right) \cdot \left(y - t\right)}1 - \frac{\frac{1}{\sqrt[3]{y - z} \cdot \sqrt[3]{y - z}} \cdot \left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right)}{\sqrt[3]{\sqrt[3]{y - z}} \cdot \sqrt[3]{\sqrt[3]{y - z}}} \cdot \frac{\frac{\sqrt[3]{x}}{y - t}}{\sqrt[3]{\sqrt[3]{y - z}}}double f(double x, double y, double z, double t) {
double r326590 = 1.0;
double r326591 = x;
double r326592 = y;
double r326593 = z;
double r326594 = r326592 - r326593;
double r326595 = t;
double r326596 = r326592 - r326595;
double r326597 = r326594 * r326596;
double r326598 = r326591 / r326597;
double r326599 = r326590 - r326598;
return r326599;
}
double f(double x, double y, double z, double t) {
double r326600 = 1.0;
double r326601 = 1.0;
double r326602 = y;
double r326603 = z;
double r326604 = r326602 - r326603;
double r326605 = cbrt(r326604);
double r326606 = r326605 * r326605;
double r326607 = r326601 / r326606;
double r326608 = x;
double r326609 = cbrt(r326608);
double r326610 = r326609 * r326609;
double r326611 = r326607 * r326610;
double r326612 = cbrt(r326605);
double r326613 = r326612 * r326612;
double r326614 = r326611 / r326613;
double r326615 = t;
double r326616 = r326602 - r326615;
double r326617 = r326609 / r326616;
double r326618 = r326617 / r326612;
double r326619 = r326614 * r326618;
double r326620 = r326600 - r326619;
return r326620;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t
Results
Initial program 0.7
rmApplied *-un-lft-identity0.7
Applied times-frac1.1
rmApplied add-cube-cbrt1.3
Applied *-un-lft-identity1.3
Applied times-frac1.3
Applied associate-*l*1.3
Simplified1.3
rmApplied add-cube-cbrt1.4
Applied *-un-lft-identity1.4
Applied add-cube-cbrt1.4
Applied times-frac1.4
Applied times-frac0.9
Applied associate-*r*0.7
Simplified0.7
Final simplification0.7
herbie shell --seed 2020062
(FPCore (x y z t)
:name "Data.Random.Distribution.Triangular:triangularCDF from random-fu-0.2.6.2, A"
:precision binary64
(- 1 (/ x (* (- y z) (- y t)))))