\sin x \cdot \frac{\sinh y}{y}\sin x \cdot \left(\frac{1}{6} \cdot {y}^{2} + \left(\frac{1}{120} \cdot {y}^{4} + 1\right)\right)double code(double x, double y) {
return ((double) (((double) sin(x)) * ((double) (((double) sinh(y)) / y))));
}
double code(double x, double y) {
return ((double) (((double) sin(x)) * ((double) (((double) (0.16666666666666666 * ((double) pow(y, 2.0)))) + ((double) (((double) (0.008333333333333333 * ((double) pow(y, 4.0)))) + 1.0))))));
}



Bits error versus x



Bits error versus y
Results
Initial program 0.0
Taylor expanded around 0 0.6
Final simplification0.6
herbie shell --seed 2020121
(FPCore (x y)
:name "Linear.Quaternion:$ccos from linear-1.19.1.3"
:precision binary64
(* (sin x) (/ (sinh y) y)))