Octave 3.8, jcobi/4

Average Error: 54.4 → 11.0

Time: 37.1s

Precision: 64

\[\alpha \gt -1 \land \beta \gt -1 \land i \gt 1\]

Math

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

↓

\[\begin{array}{l} \mathbf{if}\;i \le 1.963440419475073214093682962109391387285 \cdot 10^{118}:\\ \;\;\;\;\frac{\frac{i}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}}{1} \cdot \frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1}}{\left(\alpha + \beta\right) + i}}}{\frac{\mathsf{fma}\left(i, 2, \alpha + \beta\right)}{\frac{\mathsf{fma}\left(\beta, \alpha, i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}{\mathsf{fma}\left(i, 2, \alpha + \beta\right)}}}\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{i}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}}{1} \cdot \frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1}}{\left(\alpha + \beta\right) + i}}}{\frac{\mathsf{fma}\left(i, 2, \alpha + \beta\right)}{\mathsf{fma}\left(0.25, \alpha, \mathsf{fma}\left(0.5, i, 0.25 \cdot \beta\right)\right)}}\\ \end{array}\]

C

double f(double alpha, double beta, double i) {
        double r225724 = i;
        double r225725 = alpha;
        double r225726 = beta;
        double r225727 = r225725 + r225726;
        double r225728 = r225727 + r225724;
        double r225729 = r225724 * r225728;
        double r225730 = r225726 * r225725;
        double r225731 = r225730 + r225729;
        double r225732 = r225729 * r225731;
        double r225733 = 2.0;
        double r225734 = r225733 * r225724;
        double r225735 = r225727 + r225734;
        double r225736 = r225735 * r225735;
        double r225737 = r225732 / r225736;
        double r225738 = 1.0;
        double r225739 = r225736 - r225738;
        double r225740 = r225737 / r225739;
        return r225740;
}

↓

double f(double alpha, double beta, double i) {
        double r225741 = i;
        double r225742 = 1.9634404194750732e+118;
        bool r225743 = r225741 <= r225742;
        double r225744 = alpha;
        double r225745 = beta;
        double r225746 = r225744 + r225745;
        double r225747 = 2.0;
        double r225748 = r225747 * r225741;
        double r225749 = r225746 + r225748;
        double r225750 = 1.0;
        double r225751 = sqrt(r225750);
        double r225752 = r225749 + r225751;
        double r225753 = 1.0;
        double r225754 = r225752 / r225753;
        double r225755 = r225749 - r225751;
        double r225756 = r225746 + r225741;
        double r225757 = r225755 / r225756;
        double r225758 = r225754 * r225757;
        double r225759 = r225741 / r225758;
        double r225760 = fma(r225741, r225747, r225746);
        double r225761 = r225741 * r225756;
        double r225762 = fma(r225745, r225744, r225761);
        double r225763 = r225762 / r225760;
        double r225764 = r225760 / r225763;
        double r225765 = r225759 / r225764;
        double r225766 = 0.25;
        double r225767 = 0.5;
        double r225768 = r225766 * r225745;
        double r225769 = fma(r225767, r225741, r225768);
        double r225770 = fma(r225766, r225744, r225769);
        double r225771 = r225760 / r225770;
        double r225772 = r225759 / r225771;
        double r225773 = r225743 ? r225765 : r225772;
        return r225773;
}

Error

Bits error versus alpha

Bits error versus beta

Bits error versus i

Derivation

1. Split input into 2 regimes

2. if i < 1.9634404194750732e+118

   1. Initial program 38.6

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

   2. Simplified 34.7

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

   4. Applied associate-/l* 24.8

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

   6. Applied *-un-lft-identity 24.8

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

   7. Applied *-un-lft-identity 24.8

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

   8. Applied times-frac 24.8

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

   9. Applied times-frac 17.6

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

   10. Applied associate-/r* 14.5

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

   11. Simplified 14.4

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

   13. Applied *-un-lft-identity 14.4

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

   14. Applied add-sqr-sqrt 14.4

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

   15. Applied difference-of-squares 14.4

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

   16. Applied times-frac 9.5

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

   if 1.9634404194750732e+118 < i

   1. Initial program 64.0

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

   2. Simplified 63.7

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

   4. Applied associate-/l* 63.7

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

   6. Applied *-un-lft-identity 63.7

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

   7. Applied *-un-lft-identity 63.7

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

   8. Applied times-frac 63.7

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

   9. Applied times-frac 55.5

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

   10. Applied associate-/r* 55.3

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

   11. Simplified 55.3

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

   13. Applied *-un-lft-identity 55.3

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

   14. Applied add-sqr-sqrt 55.3

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

   15. Applied difference-of-squares 55.3

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

   16. Applied times-frac 54.9

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

   17. Taylor expanded around 0 11.9

       \[\leadsto \frac{\frac{i}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}}{1} \cdot \frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1}}{\left(\alpha + \beta\right) + i}}}{\frac{\mathsf{fma}\left(i, 2, \alpha + \beta\right)}{\color{blue}{0.25 \cdot \alpha + \left(0.5 \cdot i + 0.25 \cdot \beta\right)}}}\]

   18. Simplified 11.9

       \[\leadsto \frac{\frac{i}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}}{1} \cdot \frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1}}{\left(\alpha + \beta\right) + i}}}{\frac{\mathsf{fma}\left(i, 2, \alpha + \beta\right)}{\color{blue}{\mathsf{fma}\left(0.25, \alpha, \mathsf{fma}\left(0.5, i, 0.25 \cdot \beta\right)\right)}}}\]
3. Recombined 2 regimes into one program.

4. Final simplification 11.0

   \[\leadsto \begin{array}{l} \mathbf{if}\;i \le 1.963440419475073214093682962109391387285 \cdot 10^{118}:\\ \;\;\;\;\frac{\frac{i}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}}{1} \cdot \frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1}}{\left(\alpha + \beta\right) + i}}}{\frac{\mathsf{fma}\left(i, 2, \alpha + \beta\right)}{\frac{\mathsf{fma}\left(\beta, \alpha, i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}{\mathsf{fma}\left(i, 2, \alpha + \beta\right)}}}\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{i}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}}{1} \cdot \frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1}}{\left(\alpha + \beta\right) + i}}}{\frac{\mathsf{fma}\left(i, 2, \alpha + \beta\right)}{\mathsf{fma}\left(0.25, \alpha, \mathsf{fma}\left(0.5, i, 0.25 \cdot \beta\right)\right)}}\\ \end{array}\]

Reproduce

herbie shell --seed 2019352 +o rules:numerics
(FPCore (alpha beta i)
  :name "Octave 3.8, jcobi/4"
  :precision binary64
  :pre (and (> alpha -1) (> beta -1) (> i 1))
  (/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1)))