Henrywood and Agarwal, Equation (12)

Average Error: 25.5 → 13.4

Time: 1.5m

Precision: 64

Math

\[\left({\left(\frac{d}{h}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

↓

\[\begin{array}{l} \mathbf{if}\;\left(1 - \frac{h}{\ell} \cdot \left({\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2} \cdot \frac{1}{2}\right)\right) \cdot \left({\left(\frac{d}{\ell}\right)}^{\frac{1}{2}} \cdot {\left(\frac{d}{h}\right)}^{\frac{1}{2}}\right) \le 3.908642621051716 \cdot 10^{+20}:\\ \;\;\;\;\left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\left(\sqrt{\frac{d}{\sqrt[3]{\sqrt[3]{\ell}}}} \cdot \sqrt{\frac{1}{\sqrt[3]{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}}\right) \cdot \sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}}\right)\right) \cdot \left(1 - \frac{h}{\ell} \cdot \left({\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2} \cdot \frac{1}{2}\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\frac{\left(\left(\sqrt{\frac{d}{\sqrt[3]{\ell}}} \cdot \left(1 - \frac{\left(\frac{\frac{M}{\frac{d}{D}}}{2} \cdot \frac{\frac{M}{\frac{d}{D}}}{2}\right) \cdot h}{\ell \cdot 2}\right)\right) \cdot \sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}}\right) \cdot \left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \left|\sqrt[3]{d}\right|\right)}{\sqrt{\sqrt[3]{h} \cdot \sqrt[3]{h}}}\\ \end{array}\]

C

double f(double d, double h, double l, double M, double D) {
        double r5298004 = d;
        double r5298005 = h;
        double r5298006 = r5298004 / r5298005;
        double r5298007 = 1.0;
        double r5298008 = 2.0;
        double r5298009 = r5298007 / r5298008;
        double r5298010 = pow(r5298006, r5298009);
        double r5298011 = l;
        double r5298012 = r5298004 / r5298011;
        double r5298013 = pow(r5298012, r5298009);
        double r5298014 = r5298010 * r5298013;
        double r5298015 = M;
        double r5298016 = D;
        double r5298017 = r5298015 * r5298016;
        double r5298018 = r5298008 * r5298004;
        double r5298019 = r5298017 / r5298018;
        double r5298020 = pow(r5298019, r5298008);
        double r5298021 = r5298009 * r5298020;
        double r5298022 = r5298005 / r5298011;
        double r5298023 = r5298021 * r5298022;
        double r5298024 = r5298007 - r5298023;
        double r5298025 = r5298014 * r5298024;
        return r5298025;
}

↓

double f(double d, double h, double l, double M, double D) {
        double r5298026 = 1.0;
        double r5298027 = h;
        double r5298028 = l;
        double r5298029 = r5298027 / r5298028;
        double r5298030 = M;
        double r5298031 = D;
        double r5298032 = r5298030 * r5298031;
        double r5298033 = 2.0;
        double r5298034 = d;
        double r5298035 = r5298033 * r5298034;
        double r5298036 = r5298032 / r5298035;
        double r5298037 = pow(r5298036, r5298033);
        double r5298038 = 0.5;
        double r5298039 = r5298037 * r5298038;
        double r5298040 = r5298029 * r5298039;
        double r5298041 = r5298026 - r5298040;
        double r5298042 = r5298034 / r5298028;
        double r5298043 = pow(r5298042, r5298038);
        double r5298044 = r5298034 / r5298027;
        double r5298045 = pow(r5298044, r5298038);
        double r5298046 = r5298043 * r5298045;
        double r5298047 = r5298041 * r5298046;
        double r5298048 = 3.908642621051716e+20;
        bool r5298049 = r5298047 <= r5298048;
        double r5298050 = cbrt(r5298034);
        double r5298051 = cbrt(r5298027);
        double r5298052 = r5298050 / r5298051;
        double r5298053 = sqrt(r5298052);
        double r5298054 = r5298052 * r5298052;
        double r5298055 = sqrt(r5298054);
        double r5298056 = r5298053 * r5298055;
        double r5298057 = cbrt(r5298028);
        double r5298058 = cbrt(r5298057);
        double r5298059 = r5298034 / r5298058;
        double r5298060 = sqrt(r5298059);
        double r5298061 = r5298057 * r5298057;
        double r5298062 = cbrt(r5298061);
        double r5298063 = r5298026 / r5298062;
        double r5298064 = sqrt(r5298063);
        double r5298065 = r5298060 * r5298064;
        double r5298066 = r5298026 / r5298057;
        double r5298067 = r5298066 / r5298057;
        double r5298068 = sqrt(r5298067);
        double r5298069 = r5298065 * r5298068;
        double r5298070 = r5298056 * r5298069;
        double r5298071 = r5298070 * r5298041;
        double r5298072 = r5298034 / r5298057;
        double r5298073 = sqrt(r5298072);
        double r5298074 = r5298034 / r5298031;
        double r5298075 = r5298030 / r5298074;
        double r5298076 = r5298075 / r5298033;
        double r5298077 = r5298076 * r5298076;
        double r5298078 = r5298077 * r5298027;
        double r5298079 = r5298028 * r5298033;
        double r5298080 = r5298078 / r5298079;
        double r5298081 = r5298026 - r5298080;
        double r5298082 = r5298073 * r5298081;
        double r5298083 = r5298082 * r5298068;
        double r5298084 = fabs(r5298050);
        double r5298085 = r5298053 * r5298084;
        double r5298086 = r5298083 * r5298085;
        double r5298087 = r5298051 * r5298051;
        double r5298088 = sqrt(r5298087);
        double r5298089 = r5298086 / r5298088;
        double r5298090 = r5298049 ? r5298071 : r5298089;
        return r5298090;
}

Error

Bits error versus d

Bits error versus h

Bits error versus l

Bits error versus M

Bits error versus D

Derivation

1. Split input into 2 regimes

2. if (* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l)))) < 3.908642621051716e+20

   1. Initial program 17.7

      \[\left({\left(\frac{d}{h}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]
   2. Using strategy rm

   3. Applied add-cube-cbrt 18.0

      \[\leadsto \left({\left(\frac{d}{\color{blue}{\left(\sqrt[3]{h} \cdot \sqrt[3]{h}\right) \cdot \sqrt[3]{h}}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   4. Applied add-cube-cbrt 18.1

      \[\leadsto \left({\left(\frac{\color{blue}{\left(\sqrt[3]{d} \cdot \sqrt[3]{d}\right) \cdot \sqrt[3]{d}}}{\left(\sqrt[3]{h} \cdot \sqrt[3]{h}\right) \cdot \sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   5. Applied times-frac 18.1

      \[\leadsto \left({\color{blue}{\left(\frac{\sqrt[3]{d} \cdot \sqrt[3]{d}}{\sqrt[3]{h} \cdot \sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   6. Applied unpow-prod-down 16.7

      \[\leadsto \left(\color{blue}{\left({\left(\frac{\sqrt[3]{d} \cdot \sqrt[3]{d}}{\sqrt[3]{h} \cdot \sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   7. Simplified 16.7

      \[\leadsto \left(\left(\color{blue}{\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}}} \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   8. Simplified 16.7

      \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \color{blue}{\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}}\right) \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]
   9. Using strategy rm

   10. Applied add-cube-cbrt 16.8

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot {\left(\frac{d}{\color{blue}{\left(\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}\right) \cdot \sqrt[3]{\ell}}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   11. Applied *-un-lft-identity 16.8

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot {\left(\frac{\color{blue}{1 \cdot d}}{\left(\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}\right) \cdot \sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   12. Applied times-frac 16.8

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot {\color{blue}{\left(\frac{1}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{d}{\sqrt[3]{\ell}}\right)}}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   13. Applied unpow-prod-down 12.2

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \color{blue}{\left({\left(\frac{1}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   14. Simplified 12.2

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\color{blue}{\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}}} \cdot {\left(\frac{d}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   15. Simplified 12.2

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \color{blue}{\sqrt{\frac{d}{\sqrt[3]{\ell}}}}\right)\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]
   16. Using strategy rm

   17. Applied add-cube-cbrt 12.2

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{d}{\sqrt[3]{\color{blue}{\left(\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}\right) \cdot \sqrt[3]{\ell}}}}}\right)\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   18. Applied cbrt-prod 12.3

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{d}{\color{blue}{\sqrt[3]{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \sqrt[3]{\sqrt[3]{\ell}}}}}\right)\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   19. Applied *-un-lft-identity 12.3

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{\color{blue}{1 \cdot d}}{\sqrt[3]{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \sqrt[3]{\sqrt[3]{\ell}}}}\right)\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   20. Applied times-frac 12.3

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \sqrt{\color{blue}{\frac{1}{\sqrt[3]{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}} \cdot \frac{d}{\sqrt[3]{\sqrt[3]{\ell}}}}}\right)\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   21. Applied sqrt-prod 11.8

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \color{blue}{\left(\sqrt{\frac{1}{\sqrt[3]{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}} \cdot \sqrt{\frac{d}{\sqrt[3]{\sqrt[3]{\ell}}}}\right)}\right)\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   if 3.908642621051716e+20 < (* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l))))

   1. Initial program 32.7

      \[\left({\left(\frac{d}{h}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]
   2. Using strategy rm

   3. Applied add-cube-cbrt 32.9

      \[\leadsto \left({\left(\frac{d}{\color{blue}{\left(\sqrt[3]{h} \cdot \sqrt[3]{h}\right) \cdot \sqrt[3]{h}}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   4. Applied add-cube-cbrt 33.0

      \[\leadsto \left({\left(\frac{\color{blue}{\left(\sqrt[3]{d} \cdot \sqrt[3]{d}\right) \cdot \sqrt[3]{d}}}{\left(\sqrt[3]{h} \cdot \sqrt[3]{h}\right) \cdot \sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   5. Applied times-frac 33.0

      \[\leadsto \left({\color{blue}{\left(\frac{\sqrt[3]{d} \cdot \sqrt[3]{d}}{\sqrt[3]{h} \cdot \sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   6. Applied unpow-prod-down 24.5

      \[\leadsto \left(\color{blue}{\left({\left(\frac{\sqrt[3]{d} \cdot \sqrt[3]{d}}{\sqrt[3]{h} \cdot \sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   7. Simplified 24.5

      \[\leadsto \left(\left(\color{blue}{\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}}} \cdot {\left(\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   8. Simplified 24.5

      \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \color{blue}{\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}}\right) \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]
   9. Using strategy rm

   10. Applied add-cube-cbrt 24.6

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot {\left(\frac{d}{\color{blue}{\left(\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}\right) \cdot \sqrt[3]{\ell}}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   11. Applied *-un-lft-identity 24.6

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot {\left(\frac{\color{blue}{1 \cdot d}}{\left(\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}\right) \cdot \sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   12. Applied times-frac 24.6

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot {\color{blue}{\left(\frac{1}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{d}{\sqrt[3]{\ell}}\right)}}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   13. Applied unpow-prod-down 22.2

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \color{blue}{\left({\left(\frac{1}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   14. Simplified 22.1

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\color{blue}{\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}}} \cdot {\left(\frac{d}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   15. Simplified 22.1

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \color{blue}{\sqrt{\frac{d}{\sqrt[3]{\ell}}}}\right)\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]
   16. Using strategy rm

   17. Applied associate-*l/ 22.1

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{d}{\sqrt[3]{\ell}}}\right)\right) \cdot \left(1 - \color{blue}{\frac{1 \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}}{2}} \cdot \frac{h}{\ell}\right)\]

   18. Applied frac-times 16.7

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{d}{\sqrt[3]{\ell}}}\right)\right) \cdot \left(1 - \color{blue}{\frac{\left(1 \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot h}{2 \cdot \ell}}\right)\]

   19. Simplified 16.7

       \[\leadsto \left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{d}{\sqrt[3]{\ell}}}\right)\right) \cdot \left(1 - \frac{\color{blue}{h \cdot \left(\frac{M \cdot D}{d \cdot 2} \cdot \frac{M \cdot D}{d \cdot 2}\right)}}{2 \cdot \ell}\right)\]
   20. Using strategy rm

   21. Applied frac-times 16.7

       \[\leadsto \left(\left(\sqrt{\color{blue}{\frac{\sqrt[3]{d} \cdot \sqrt[3]{d}}{\sqrt[3]{h} \cdot \sqrt[3]{h}}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{d}{\sqrt[3]{\ell}}}\right)\right) \cdot \left(1 - \frac{h \cdot \left(\frac{M \cdot D}{d \cdot 2} \cdot \frac{M \cdot D}{d \cdot 2}\right)}{2 \cdot \ell}\right)\]

   22. Applied sqrt-div 15.4

       \[\leadsto \left(\left(\color{blue}{\frac{\sqrt{\sqrt[3]{d} \cdot \sqrt[3]{d}}}{\sqrt{\sqrt[3]{h} \cdot \sqrt[3]{h}}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{d}{\sqrt[3]{\ell}}}\right)\right) \cdot \left(1 - \frac{h \cdot \left(\frac{M \cdot D}{d \cdot 2} \cdot \frac{M \cdot D}{d \cdot 2}\right)}{2 \cdot \ell}\right)\]

   23. Applied associate-*l/ 15.4

       \[\leadsto \left(\color{blue}{\frac{\sqrt{\sqrt[3]{d} \cdot \sqrt[3]{d}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}}{\sqrt{\sqrt[3]{h} \cdot \sqrt[3]{h}}}} \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{d}{\sqrt[3]{\ell}}}\right)\right) \cdot \left(1 - \frac{h \cdot \left(\frac{M \cdot D}{d \cdot 2} \cdot \frac{M \cdot D}{d \cdot 2}\right)}{2 \cdot \ell}\right)\]

   24. Applied associate-*l/ 15.6

       \[\leadsto \color{blue}{\frac{\left(\sqrt{\sqrt[3]{d} \cdot \sqrt[3]{d}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{d}{\sqrt[3]{\ell}}}\right)}{\sqrt{\sqrt[3]{h} \cdot \sqrt[3]{h}}}} \cdot \left(1 - \frac{h \cdot \left(\frac{M \cdot D}{d \cdot 2} \cdot \frac{M \cdot D}{d \cdot 2}\right)}{2 \cdot \ell}\right)\]

   25. Applied associate-*l/ 15.5

       \[\leadsto \color{blue}{\frac{\left(\left(\sqrt{\sqrt[3]{d} \cdot \sqrt[3]{d}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{d}{\sqrt[3]{\ell}}}\right)\right) \cdot \left(1 - \frac{h \cdot \left(\frac{M \cdot D}{d \cdot 2} \cdot \frac{M \cdot D}{d \cdot 2}\right)}{2 \cdot \ell}\right)}{\sqrt{\sqrt[3]{h} \cdot \sqrt[3]{h}}}}\]

   26. Simplified 15.0

       \[\leadsto \frac{\color{blue}{\left(\left|\sqrt[3]{d}\right| \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}} \cdot \left(\sqrt{\frac{d}{\sqrt[3]{\ell}}} \cdot \left(1 - \frac{\left(\frac{\frac{M}{\frac{d}{D}}}{2} \cdot \frac{\frac{M}{\frac{d}{D}}}{2}\right) \cdot h}{\ell \cdot 2}\right)\right)\right)}}{\sqrt{\sqrt[3]{h} \cdot \sqrt[3]{h}}}\]
3. Recombined 2 regimes into one program.

4. Final simplification 13.4

   \[\leadsto \begin{array}{l} \mathbf{if}\;\left(1 - \frac{h}{\ell} \cdot \left({\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2} \cdot \frac{1}{2}\right)\right) \cdot \left({\left(\frac{d}{\ell}\right)}^{\frac{1}{2}} \cdot {\left(\frac{d}{h}\right)}^{\frac{1}{2}}\right) \le 3.908642621051716 \cdot 10^{+20}:\\ \;\;\;\;\left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\left(\sqrt{\frac{d}{\sqrt[3]{\sqrt[3]{\ell}}}} \cdot \sqrt{\frac{1}{\sqrt[3]{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}}\right) \cdot \sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}}\right)\right) \cdot \left(1 - \frac{h}{\ell} \cdot \left({\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2} \cdot \frac{1}{2}\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\frac{\left(\left(\sqrt{\frac{d}{\sqrt[3]{\ell}}} \cdot \left(1 - \frac{\left(\frac{\frac{M}{\frac{d}{D}}}{2} \cdot \frac{\frac{M}{\frac{d}{D}}}{2}\right) \cdot h}{\ell \cdot 2}\right)\right) \cdot \sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}}\right) \cdot \left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \left|\sqrt[3]{d}\right|\right)}{\sqrt{\sqrt[3]{h} \cdot \sqrt[3]{h}}}\\ \end{array}\]

Reproduce

herbie shell --seed 2019151 +o rules:numerics
(FPCore (d h l M D)
  :name "Henrywood and Agarwal, Equation (12)"
  (* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l)))))