Henrywood and Agarwal, Equation (12)

Average Error: 26.7 → 11.1

Time: 1.7m

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.0718122983865245 \cdot 10^{-118}:\\ \;\;\;\;\left(1 - \frac{\left(\frac{\frac{M}{\frac{d}{D}}}{2} \cdot \frac{h}{\ell}\right) \cdot \frac{\frac{M}{\frac{d}{D}}}{2}}{2}\right) \cdot \left({\left(\frac{d}{\ell}\right)}^{\frac{1}{2}} \cdot \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)\right)\\ \mathbf{else}:\\ \;\;\;\;\frac{\left(\left(\left|\sqrt[3]{d}\right| \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}}\right) \cdot \left(\left|\sqrt[3]{d}\right| \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right)\right) \cdot \left(1 - \frac{\frac{\frac{M \cdot D}{d}}{2} \cdot \left(h \cdot \frac{\frac{M \cdot D}{d}}{2}\right)}{\ell \cdot 2}\right)}{\sqrt{\sqrt[3]{h} \cdot \sqrt[3]{h}} \cdot \sqrt{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}\\ \end{array}\]

C

double f(double d, double h, double l, double M, double D) {
        double r5191205 = d;
        double r5191206 = h;
        double r5191207 = r5191205 / r5191206;
        double r5191208 = 1.0;
        double r5191209 = 2.0;
        double r5191210 = r5191208 / r5191209;
        double r5191211 = pow(r5191207, r5191210);
        double r5191212 = l;
        double r5191213 = r5191205 / r5191212;
        double r5191214 = pow(r5191213, r5191210);
        double r5191215 = r5191211 * r5191214;
        double r5191216 = M;
        double r5191217 = D;
        double r5191218 = r5191216 * r5191217;
        double r5191219 = r5191209 * r5191205;
        double r5191220 = r5191218 / r5191219;
        double r5191221 = pow(r5191220, r5191209);
        double r5191222 = r5191210 * r5191221;
        double r5191223 = r5191206 / r5191212;
        double r5191224 = r5191222 * r5191223;
        double r5191225 = r5191208 - r5191224;
        double r5191226 = r5191215 * r5191225;
        return r5191226;
}

↓

double f(double d, double h, double l, double M, double D) {
        double r5191227 = 1.0;
        double r5191228 = h;
        double r5191229 = l;
        double r5191230 = r5191228 / r5191229;
        double r5191231 = M;
        double r5191232 = D;
        double r5191233 = r5191231 * r5191232;
        double r5191234 = 2.0;
        double r5191235 = d;
        double r5191236 = r5191234 * r5191235;
        double r5191237 = r5191233 / r5191236;
        double r5191238 = pow(r5191237, r5191234);
        double r5191239 = 0.5;
        double r5191240 = r5191238 * r5191239;
        double r5191241 = r5191230 * r5191240;
        double r5191242 = r5191227 - r5191241;
        double r5191243 = r5191235 / r5191229;
        double r5191244 = pow(r5191243, r5191239);
        double r5191245 = r5191235 / r5191228;
        double r5191246 = pow(r5191245, r5191239);
        double r5191247 = r5191244 * r5191246;
        double r5191248 = r5191242 * r5191247;
        double r5191249 = -3.0718122983865245e-118;
        bool r5191250 = r5191248 <= r5191249;
        double r5191251 = r5191235 / r5191232;
        double r5191252 = r5191231 / r5191251;
        double r5191253 = r5191252 / r5191234;
        double r5191254 = r5191253 * r5191230;
        double r5191255 = r5191254 * r5191253;
        double r5191256 = r5191255 / r5191234;
        double r5191257 = r5191227 - r5191256;
        double r5191258 = cbrt(r5191235);
        double r5191259 = cbrt(r5191228);
        double r5191260 = r5191258 / r5191259;
        double r5191261 = sqrt(r5191260);
        double r5191262 = r5191260 * r5191260;
        double r5191263 = sqrt(r5191262);
        double r5191264 = r5191261 * r5191263;
        double r5191265 = r5191244 * r5191264;
        double r5191266 = r5191257 * r5191265;
        double r5191267 = fabs(r5191258);
        double r5191268 = cbrt(r5191229);
        double r5191269 = r5191258 / r5191268;
        double r5191270 = sqrt(r5191269);
        double r5191271 = r5191267 * r5191270;
        double r5191272 = r5191267 * r5191261;
        double r5191273 = r5191271 * r5191272;
        double r5191274 = r5191233 / r5191235;
        double r5191275 = r5191274 / r5191234;
        double r5191276 = r5191228 * r5191275;
        double r5191277 = r5191275 * r5191276;
        double r5191278 = r5191229 * r5191234;
        double r5191279 = r5191277 / r5191278;
        double r5191280 = r5191227 - r5191279;
        double r5191281 = r5191273 * r5191280;
        double r5191282 = r5191259 * r5191259;
        double r5191283 = sqrt(r5191282);
        double r5191284 = r5191268 * r5191268;
        double r5191285 = sqrt(r5191284);
        double r5191286 = r5191283 * r5191285;
        double r5191287 = r5191281 / r5191286;
        double r5191288 = r5191250 ? r5191266 : r5191287;
        return r5191288;
}

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.0718122983865245e-118

   1. Initial program 27.9

      \[\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 28.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 28.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 28.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 27.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 27.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 27.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 associate-*l/ 27.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(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\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)\]

   11. Applied associate-*l/ 27.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(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \color{blue}{\frac{\left(1 \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}}{2}}\right)\]

   12. Simplified 20.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(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \frac{\color{blue}{\left(\frac{h}{\ell} \cdot \frac{\frac{M}{\frac{d}{D}}}{2}\right) \cdot \frac{\frac{M}{\frac{d}{D}}}{2}}}{2}\right)\]

   if -3.0718122983865245e-118 < (* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l))))

   1. Initial program 26.4

      \[\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 26.7

      \[\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 26.8

      \[\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 26.8

      \[\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 20.4

      \[\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 20.4

      \[\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 20.4

      \[\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 20.5

       \[\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 add-cube-cbrt 20.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(\frac{\color{blue}{\left(\sqrt[3]{d} \cdot \sqrt[3]{d}\right) \cdot \sqrt[3]{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 20.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 {\color{blue}{\left(\frac{\sqrt[3]{d} \cdot \sqrt[3]{d}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{\sqrt[3]{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 15.9

       \[\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{\sqrt[3]{d} \cdot \sqrt[3]{d}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{\sqrt[3]{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 15.9

       \[\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{\sqrt[3]{d}}{\sqrt[3]{\ell}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}}} \cdot {\left(\frac{\sqrt[3]{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 15.9

       \[\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{\sqrt[3]{d}}{\sqrt[3]{\ell}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}} \cdot \color{blue}{\sqrt{\frac{\sqrt[3]{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/ 15.9

       \[\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{\sqrt[3]{d}}{\sqrt[3]{\ell}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{\sqrt[3]{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 13.0

       \[\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{\sqrt[3]{d}}{\sqrt[3]{\ell}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{\sqrt[3]{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 13.0

       \[\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{\sqrt[3]{d}}{\sqrt[3]{\ell}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{\sqrt[3]{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 13.0

       \[\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{\color{blue}{\frac{\sqrt[3]{d} \cdot \sqrt[3]{d}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}} \cdot \sqrt{\frac{\sqrt[3]{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 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}{\frac{\sqrt{\sqrt[3]{d} \cdot \sqrt[3]{d}}}{\sqrt{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}} \cdot \sqrt{\frac{\sqrt[3]{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/ 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}{\frac{\sqrt{\sqrt[3]{d} \cdot \sqrt[3]{d}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}}}{\sqrt{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}}\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 frac-times 12.2

       \[\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 \frac{\sqrt{\sqrt[3]{d} \cdot \sqrt[3]{d}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}}}{\sqrt{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}\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)\]

   25. Applied sqrt-div 11.3

       \[\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 \frac{\sqrt{\sqrt[3]{d} \cdot \sqrt[3]{d}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}}}{\sqrt{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}\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)\]

   26. Applied associate-*l/ 11.3

       \[\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 \frac{\sqrt{\sqrt[3]{d} \cdot \sqrt[3]{d}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}}}{\sqrt{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}\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)\]

   27. Applied frac-times 11.3

       \[\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{\sqrt[3]{d} \cdot \sqrt[3]{d}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}}\right)}{\sqrt{\sqrt[3]{h} \cdot \sqrt[3]{h}} \cdot \sqrt{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}} \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)\]

   28. Applied associate-*l/ 10.9

       \[\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{\sqrt[3]{d} \cdot \sqrt[3]{d}} \cdot \sqrt{\frac{\sqrt[3]{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}} \cdot \sqrt{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}}\]

   29. Simplified 9.1

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

4. Final simplification 11.1

   \[\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.0718122983865245 \cdot 10^{-118}:\\ \;\;\;\;\left(1 - \frac{\left(\frac{\frac{M}{\frac{d}{D}}}{2} \cdot \frac{h}{\ell}\right) \cdot \frac{\frac{M}{\frac{d}{D}}}{2}}{2}\right) \cdot \left({\left(\frac{d}{\ell}\right)}^{\frac{1}{2}} \cdot \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)\right)\\ \mathbf{else}:\\ \;\;\;\;\frac{\left(\left(\left|\sqrt[3]{d}\right| \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}}\right) \cdot \left(\left|\sqrt[3]{d}\right| \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right)\right) \cdot \left(1 - \frac{\frac{\frac{M \cdot D}{d}}{2} \cdot \left(h \cdot \frac{\frac{M \cdot D}{d}}{2}\right)}{\ell \cdot 2}\right)}{\sqrt{\sqrt[3]{h} \cdot \sqrt[3]{h}} \cdot \sqrt{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}\\ \end{array}\]

Reproduce

herbie shell --seed 2019142 
(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)))))