\left(J \cdot \left(e^{\ell} - e^{-\ell}\right)\right) \cdot \cos \left(\frac{K}{2}\right) + U\left(\frac{1}{3} \cdot {\ell}^{3} + \left(\frac{1}{60} \cdot {\ell}^{5} + 2 \cdot \ell\right)\right) \cdot \left(J \cdot \cos \left(\frac{K}{2}\right)\right) + Udouble code(double J, double l, double K, double U) {
return ((double) (((double) (((double) (J * ((double) (((double) exp(l)) - ((double) exp(((double) -(l)))))))) * ((double) cos(((double) (K / 2.0)))))) + U));
}
double code(double J, double l, double K, double U) {
return ((double) (((double) (((double) (((double) (0.3333333333333333 * ((double) pow(l, 3.0)))) + ((double) (((double) (0.016666666666666666 * ((double) pow(l, 5.0)))) + ((double) (2.0 * l)))))) * ((double) (J * ((double) cos(((double) (K / 2.0)))))))) + U));
}



Bits error versus J



Bits error versus l



Bits error versus K



Bits error versus U
Results
Initial program 16.9
Taylor expanded around 0 0.4
rmApplied *-commutative0.4
Applied associate-*l*0.4
Final simplification0.4
herbie shell --seed 2020113
(FPCore (J l K U)
:name "Maksimov and Kolovsky, Equation (4)"
:precision binary64
(+ (* (* J (- (exp l) (exp (- l)))) (cos (/ K 2))) U))