\frac{e^{x} - e^{-x}}{2}\frac{\left(\frac{1}{3} \cdot {x}^{3} + \frac{1}{60} \cdot {x}^{5}\right) + 2 \cdot x}{2}double f(double x) {
double r54490 = x;
double r54491 = exp(r54490);
double r54492 = -r54490;
double r54493 = exp(r54492);
double r54494 = r54491 - r54493;
double r54495 = 2.0;
double r54496 = r54494 / r54495;
return r54496;
}
double f(double x) {
double r54497 = 0.3333333333333333;
double r54498 = x;
double r54499 = 3.0;
double r54500 = pow(r54498, r54499);
double r54501 = r54497 * r54500;
double r54502 = 0.016666666666666666;
double r54503 = 5.0;
double r54504 = pow(r54498, r54503);
double r54505 = r54502 * r54504;
double r54506 = r54501 + r54505;
double r54507 = 2.0;
double r54508 = r54507 * r54498;
double r54509 = r54506 + r54508;
double r54510 = 2.0;
double r54511 = r54509 / r54510;
return r54511;
}



Bits error versus x
Results
Initial program 58.1
Taylor expanded around 0 0.6
rmApplied associate-+r+0.6
Final simplification0.6
herbie shell --seed 2019326
(FPCore (x)
:name "Hyperbolic sine"
:precision binary64
(/ (- (exp x) (exp (- x))) 2))