Average Error: 24.1 → 11.0
Time: 14.4s
Precision: 64
\[\alpha \gt -1 \land \beta \gt -1 \land i \gt 0.0\]
\[\frac{\frac{\frac{\left(\alpha + \beta\right) \cdot \left(\beta - \alpha\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1}{2}\]
↓
\[\begin{array}{l} \mathbf{if}\;\alpha \le 1.63898582180695576 \cdot 10^{156}:\\ \;\;\;\;\frac{\sqrt[3]{{\left(\left(\alpha + \beta\right) \cdot \left(\frac{1}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}} \cdot \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}}\right) + 1\right)}^{3}}}{2}\\ \mathbf{else}:\\ \;\;\;\;\frac{\left(2 \cdot \frac{1}{\alpha} + 8 \cdot \frac{1}{{\alpha}^{3}}\right) - 4 \cdot \frac{1}{{\alpha}^{2}}}{2}\\ \end{array}\]
\frac{\frac{\frac{\left(\alpha + \beta\right) \cdot \left(\beta - \alpha\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1}{2}
↓
\begin{array}{l}
\mathbf{if}\;\alpha \le 1.63898582180695576 \cdot 10^{156}:\\
\;\;\;\;\frac{\sqrt[3]{{\left(\left(\alpha + \beta\right) \cdot \left(\frac{1}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}} \cdot \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}}\right) + 1\right)}^{3}}}{2}\\

\mathbf{else}:\\
\;\;\;\;\frac{\left(2 \cdot \frac{1}{\alpha} + 8 \cdot \frac{1}{{\alpha}^{3}}\right) - 4 \cdot \frac{1}{{\alpha}^{2}}}{2}\\

\end{array}
double f(double alpha, double beta, double i) {
        double r110281 = alpha;
        double r110282 = beta;
        double r110283 = r110281 + r110282;
        double r110284 = r110282 - r110281;
        double r110285 = r110283 * r110284;
        double r110286 = 2.0;
        double r110287 = i;
        double r110288 = r110286 * r110287;
        double r110289 = r110283 + r110288;
        double r110290 = r110285 / r110289;
        double r110291 = r110289 + r110286;
        double r110292 = r110290 / r110291;
        double r110293 = 1.0;
        double r110294 = r110292 + r110293;
        double r110295 = r110294 / r110286;
        return r110295;
}
↓
double f(double alpha, double beta, double i) {
        double r110296 = alpha;
        double r110297 = 1.6389858218069558e+156;
        bool r110298 = r110296 <= r110297;
        double r110299 = beta;
        double r110300 = r110296 + r110299;
        double r110301 = 1.0;
        double r110302 = 2.0;
        double r110303 = i;
        double r110304 = r110302 * r110303;
        double r110305 = r110300 + r110304;
        double r110306 = r110305 + r110302;
        double r110307 = sqrt(r110306);
        double r110308 = r110301 / r110307;
        double r110309 = r110299 - r110296;
        double r110310 = r110309 / r110305;
        double r110311 = r110310 / r110307;
        double r110312 = r110308 * r110311;
        double r110313 = r110300 * r110312;
        double r110314 = 1.0;
        double r110315 = r110313 + r110314;
        double r110316 = 3.0;
        double r110317 = pow(r110315, r110316);
        double r110318 = cbrt(r110317);
        double r110319 = r110318 / r110302;
        double r110320 = r110301 / r110296;
        double r110321 = r110302 * r110320;
        double r110322 = 8.0;
        double r110323 = pow(r110296, r110316);
        double r110324 = r110301 / r110323;
        double r110325 = r110322 * r110324;
        double r110326 = r110321 + r110325;
        double r110327 = 4.0;
        double r110328 = 2.0;
        double r110329 = pow(r110296, r110328);
        double r110330 = r110301 / r110329;
        double r110331 = r110327 * r110330;
        double r110332 = r110326 - r110331;
        double r110333 = r110332 / r110302;
        double r110334 = r110298 ? r110319 : r110333;
        return r110334;
}

Error

Bits error versus alpha

Bits error versus beta

Bits error versus i

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 2 regimes

  2. if alpha < 1.6389858218069558e+156

    1. Initial program 16.1

      \[\frac{\frac{\frac{\left(\alpha + \beta\right) \cdot \left(\beta - \alpha\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1}{2}\]
    2. Using strategy rm

    3. Applied *-un-lft-identity16.1

      \[\leadsto \frac{\frac{\frac{\left(\alpha + \beta\right) \cdot \left(\beta - \alpha\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\color{blue}{1 \cdot \left(\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2\right)}} + 1}{2}\]

    4. Applied *-un-lft-identity16.1

      \[\leadsto \frac{\frac{\frac{\left(\alpha + \beta\right) \cdot \left(\beta - \alpha\right)}{\color{blue}{1 \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}}{1 \cdot \left(\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2\right)} + 1}{2}\]

    5. Applied times-frac5.2

      \[\leadsto \frac{\frac{\color{blue}{\frac{\alpha + \beta}{1} \cdot \frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}}{1 \cdot \left(\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2\right)} + 1}{2}\]

    6. Applied times-frac5.2

      \[\leadsto \frac{\color{blue}{\frac{\frac{\alpha + \beta}{1}}{1} \cdot \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}} + 1}{2}\]

    7. Simplified5.2

      \[\leadsto \frac{\color{blue}{\left(\alpha + \beta\right)} \cdot \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1}{2}\]
    8. Using strategy rm

    9. Applied add-cbrt-cube5.2

      \[\leadsto \frac{\color{blue}{\sqrt[3]{\left(\left(\left(\alpha + \beta\right) \cdot \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1\right) \cdot \left(\left(\alpha + \beta\right) \cdot \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1\right)\right) \cdot \left(\left(\alpha + \beta\right) \cdot \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1\right)}}}{2}\]

    10. Simplified5.2

      \[\leadsto \frac{\sqrt[3]{\color{blue}{{\left(\left(\alpha + \beta\right) \cdot \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1\right)}^{3}}}}{2}\]
    11. Using strategy rm

    12. Applied add-sqr-sqrt5.2

      \[\leadsto \frac{\sqrt[3]{{\left(\left(\alpha + \beta\right) \cdot \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\color{blue}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} \cdot \sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}}} + 1\right)}^{3}}}{2}\]

    13. Applied *-un-lft-identity5.2

      \[\leadsto \frac{\sqrt[3]{{\left(\left(\alpha + \beta\right) \cdot \frac{\frac{\beta - \alpha}{\color{blue}{1 \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} \cdot \sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}} + 1\right)}^{3}}}{2}\]

    14. Applied *-un-lft-identity5.2

      \[\leadsto \frac{\sqrt[3]{{\left(\left(\alpha + \beta\right) \cdot \frac{\frac{\color{blue}{1 \cdot \left(\beta - \alpha\right)}}{1 \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} \cdot \sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}} + 1\right)}^{3}}}{2}\]

    15. Applied times-frac5.2

      \[\leadsto \frac{\sqrt[3]{{\left(\left(\alpha + \beta\right) \cdot \frac{\color{blue}{\frac{1}{1} \cdot \frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} \cdot \sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}} + 1\right)}^{3}}}{2}\]

    16. Applied times-frac5.2

      \[\leadsto \frac{\sqrt[3]{{\left(\left(\alpha + \beta\right) \cdot \color{blue}{\left(\frac{\frac{1}{1}}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}} \cdot \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}}\right)} + 1\right)}^{3}}}{2}\]

    17. Simplified5.2

      \[\leadsto \frac{\sqrt[3]{{\left(\left(\alpha + \beta\right) \cdot \left(\color{blue}{\frac{1}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}}} \cdot \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}}\right) + 1\right)}^{3}}}{2}\]

    if 1.6389858218069558e+156 < alpha

    1. Initial program 64.0

      \[\frac{\frac{\frac{\left(\alpha + \beta\right) \cdot \left(\beta - \alpha\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1}{2}\]

    2. Taylor expanded around inf 39.8

      \[\leadsto \frac{\color{blue}{\left(2 \cdot \frac{1}{\alpha} + 8 \cdot \frac{1}{{\alpha}^{3}}\right) - 4 \cdot \frac{1}{{\alpha}^{2}}}}{2}\]
  3. Recombined 2 regimes into one program.

  4. Final simplification11.0

    \[\leadsto \begin{array}{l} \mathbf{if}\;\alpha \le 1.63898582180695576 \cdot 10^{156}:\\ \;\;\;\;\frac{\sqrt[3]{{\left(\left(\alpha + \beta\right) \cdot \left(\frac{1}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}} \cdot \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}}\right) + 1\right)}^{3}}}{2}\\ \mathbf{else}:\\ \;\;\;\;\frac{\left(2 \cdot \frac{1}{\alpha} + 8 \cdot \frac{1}{{\alpha}^{3}}\right) - 4 \cdot \frac{1}{{\alpha}^{2}}}{2}\\ \end{array}\]

Reproduce

herbie shell --seed 2020018 
(FPCore (alpha beta i)
  :name "Octave 3.8, jcobi/2"
  :precision binary64
  :pre (and (> alpha -1) (> beta -1) (> i 0.0))
  (/ (+ (/ (/ (* (+ alpha beta) (- beta alpha)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) 2)) 1) 2))