\left(\frac{\pi}{2} \cdot \frac{1}{b \cdot b - a \cdot a}\right) \cdot \left(\frac{1}{a} - \frac{1}{b}\right)\frac{\frac{0.5 \cdot \left(\frac{\pi}{a} - \frac{\pi}{b}\right)}{b + a}}{b - a}double f(double a, double b) {
double r54462 = atan2(1.0, 0.0);
double r54463 = 2.0;
double r54464 = r54462 / r54463;
double r54465 = 1.0;
double r54466 = b;
double r54467 = r54466 * r54466;
double r54468 = a;
double r54469 = r54468 * r54468;
double r54470 = r54467 - r54469;
double r54471 = r54465 / r54470;
double r54472 = r54464 * r54471;
double r54473 = r54465 / r54468;
double r54474 = r54465 / r54466;
double r54475 = r54473 - r54474;
double r54476 = r54472 * r54475;
return r54476;
}
double f(double a, double b) {
double r54477 = 0.5;
double r54478 = atan2(1.0, 0.0);
double r54479 = a;
double r54480 = r54478 / r54479;
double r54481 = b;
double r54482 = r54478 / r54481;
double r54483 = r54480 - r54482;
double r54484 = r54477 * r54483;
double r54485 = r54481 + r54479;
double r54486 = r54484 / r54485;
double r54487 = r54481 - r54479;
double r54488 = r54486 / r54487;
return r54488;
}



Bits error versus a



Bits error versus b
Results
Initial program 15.0
rmApplied difference-of-squares9.9
Applied associate-/r*9.4
rmApplied associate-*r/9.3
Applied associate-*l/0.3
rmApplied associate-*r/0.3
Applied associate-*l/0.3
Taylor expanded around 0 0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2019353 +o rules:numerics
(FPCore (a b)
:name "NMSE Section 6.1 mentioned, B"
:precision binary64
(* (* (/ PI 2) (/ 1 (- (* b b) (* a a)))) (- (/ 1 a) (/ 1 b))))