\frac{e^{x} - e^{-x}}{2}\frac{\left(\sqrt{e^{x}} + \sqrt{e^{-x}}\right) \cdot \left(\frac{1}{24} \cdot {x}^{3} + \left(\frac{1}{1920} \cdot {x}^{5} + x\right)\right)}{2}double code(double x) {
return ((exp(x) - exp(-x)) / 2.0);
}
double code(double x) {
return (((sqrt(exp(x)) + sqrt(exp(-x))) * ((0.041666666666666664 * pow(x, 3.0)) + ((0.0005208333333333333 * pow(x, 5.0)) + x))) / 2.0);
}



Bits error versus x
Results
Initial program 57.9
rmApplied add-sqr-sqrt57.9
Applied add-sqr-sqrt58.0
Applied difference-of-squares58.0
Taylor expanded around 0 0.6
Final simplification0.6
herbie shell --seed 2020057
(FPCore (x)
:name "Hyperbolic sine"
:precision binary64
(/ (- (exp x) (exp (- x))) 2))