\frac{1 + \left(2 - \frac{\frac{2}{t}}{1 + \frac{1}{t}}\right) \cdot \left(2 - \frac{\frac{2}{t}}{1 + \frac{1}{t}}\right)}{2 + \left(2 - \frac{\frac{2}{t}}{1 + \frac{1}{t}}\right) \cdot \left(2 - \frac{\frac{2}{t}}{1 + \frac{1}{t}}\right)}\frac{1 + \left(2 - \frac{2}{1 + t \cdot 1}\right) \cdot \left(2 - \frac{2}{1 + t \cdot 1}\right)}{2 + \left(2 - \frac{2}{1 + t \cdot 1}\right) \cdot \left(2 - \frac{2}{1 + t \cdot 1}\right)}double f(double t) {
double r31971 = 1.0;
double r31972 = 2.0;
double r31973 = t;
double r31974 = r31972 / r31973;
double r31975 = r31971 / r31973;
double r31976 = r31971 + r31975;
double r31977 = r31974 / r31976;
double r31978 = r31972 - r31977;
double r31979 = r31978 * r31978;
double r31980 = r31971 + r31979;
double r31981 = r31972 + r31979;
double r31982 = r31980 / r31981;
return r31982;
}
double f(double t) {
double r31983 = 1.0;
double r31984 = 2.0;
double r31985 = t;
double r31986 = r31985 * r31983;
double r31987 = r31983 + r31986;
double r31988 = r31984 / r31987;
double r31989 = r31984 - r31988;
double r31990 = r31989 * r31989;
double r31991 = r31983 + r31990;
double r31992 = r31984 + r31990;
double r31993 = r31991 / r31992;
return r31993;
}



Bits error versus t
Results
Initial program 0.1
Simplified0.0
Final simplification0.0
herbie shell --seed 2019323
(FPCore (t)
:name "Kahan p13 Example 2"
:precision binary64
(/ (+ 1 (* (- 2 (/ (/ 2 t) (+ 1 (/ 1 t)))) (- 2 (/ (/ 2 t) (+ 1 (/ 1 t)))))) (+ 2 (* (- 2 (/ (/ 2 t) (+ 1 (/ 1 t)))) (- 2 (/ (/ 2 t) (+ 1 (/ 1 t))))))))