\frac{1}{2} \cdot \log \left(\frac{1 + x}{1 - x}\right)\frac{1}{2} \cdot \left({x}^{5} \cdot \frac{2}{5} + \left(2 \cdot x + \left(x \cdot \left(x \cdot \frac{2}{3}\right)\right) \cdot x\right)\right)double f(double x) {
double r7964313 = 1.0;
double r7964314 = 2.0;
double r7964315 = r7964313 / r7964314;
double r7964316 = x;
double r7964317 = r7964313 + r7964316;
double r7964318 = r7964313 - r7964316;
double r7964319 = r7964317 / r7964318;
double r7964320 = log(r7964319);
double r7964321 = r7964315 * r7964320;
return r7964321;
}
double f(double x) {
double r7964322 = 0.5;
double r7964323 = x;
double r7964324 = 5.0;
double r7964325 = pow(r7964323, r7964324);
double r7964326 = 0.4;
double r7964327 = r7964325 * r7964326;
double r7964328 = 2.0;
double r7964329 = r7964328 * r7964323;
double r7964330 = 0.6666666666666666;
double r7964331 = r7964323 * r7964330;
double r7964332 = r7964323 * r7964331;
double r7964333 = r7964332 * r7964323;
double r7964334 = r7964329 + r7964333;
double r7964335 = r7964327 + r7964334;
double r7964336 = r7964322 * r7964335;
return r7964336;
}



Bits error versus x
Results
Initial program 58.7
Simplified58.7
Taylor expanded around 0 0.2
Simplified0.2
rmApplied distribute-rgt-in0.2
Final simplification0.2
herbie shell --seed 2019128
(FPCore (x)
:name "Hyperbolic arc-(co)tangent"
(* (/ 1 2) (log (/ (+ 1 x) (- 1 x)))))