Henrywood and Agarwal, Equation (12)

Average Error: 26.0 → 12.1

Time: 2.0m

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

C

double f(double d, double h, double l, double M, double D) {
        double r4025414 = d;
        double r4025415 = h;
        double r4025416 = r4025414 / r4025415;
        double r4025417 = 1.0;
        double r4025418 = 2.0;
        double r4025419 = r4025417 / r4025418;
        double r4025420 = pow(r4025416, r4025419);
        double r4025421 = l;
        double r4025422 = r4025414 / r4025421;
        double r4025423 = pow(r4025422, r4025419);
        double r4025424 = r4025420 * r4025423;
        double r4025425 = M;
        double r4025426 = D;
        double r4025427 = r4025425 * r4025426;
        double r4025428 = r4025418 * r4025414;
        double r4025429 = r4025427 / r4025428;
        double r4025430 = pow(r4025429, r4025418);
        double r4025431 = r4025419 * r4025430;
        double r4025432 = r4025415 / r4025421;
        double r4025433 = r4025431 * r4025432;
        double r4025434 = r4025417 - r4025433;
        double r4025435 = r4025424 * r4025434;
        return r4025435;
}

↓

double f(double d, double h, double l, double M, double D) {
        double r4025436 = 1.0;
        double r4025437 = h;
        double r4025438 = l;
        double r4025439 = r4025437 / r4025438;
        double r4025440 = M;
        double r4025441 = D;
        double r4025442 = r4025440 * r4025441;
        double r4025443 = 2.0;
        double r4025444 = d;
        double r4025445 = r4025443 * r4025444;
        double r4025446 = r4025442 / r4025445;
        double r4025447 = pow(r4025446, r4025443);
        double r4025448 = 0.5;
        double r4025449 = r4025447 * r4025448;
        double r4025450 = r4025439 * r4025449;
        double r4025451 = r4025436 - r4025450;
        double r4025452 = r4025444 / r4025438;
        double r4025453 = pow(r4025452, r4025448);
        double r4025454 = r4025444 / r4025437;
        double r4025455 = pow(r4025454, r4025448);
        double r4025456 = r4025453 * r4025455;
        double r4025457 = r4025451 * r4025456;
        double r4025458 = -2.604184409542022e+233;
        bool r4025459 = r4025457 <= r4025458;
        double r4025460 = cbrt(r4025437);
        double r4025461 = cbrt(r4025438);
        double r4025462 = r4025460 / r4025461;
        double r4025463 = r4025441 / r4025444;
        double r4025464 = r4025440 / r4025443;
        double r4025465 = r4025463 * r4025464;
        double r4025466 = r4025462 * r4025465;
        double r4025467 = r4025466 * r4025466;
        double r4025468 = r4025448 * r4025467;
        double r4025469 = r4025462 * r4025468;
        double r4025470 = r4025436 - r4025469;
        double r4025471 = r4025470 * r4025456;
        double r4025472 = 6.854055542649234e+166;
        bool r4025473 = r4025457 <= r4025472;
        double r4025474 = cbrt(r4025444);
        double r4025475 = r4025474 / r4025460;
        double r4025476 = sqrt(r4025475);
        double r4025477 = fabs(r4025474);
        double r4025478 = r4025476 * r4025477;
        double r4025479 = r4025478 * r4025477;
        double r4025480 = r4025474 / r4025461;
        double r4025481 = sqrt(r4025480);
        double r4025482 = r4025479 * r4025481;
        double r4025483 = r4025442 / r4025443;
        double r4025484 = r4025483 / r4025444;
        double r4025485 = r4025484 * r4025484;
        double r4025486 = r4025448 * r4025485;
        double r4025487 = r4025439 * r4025486;
        double r4025488 = r4025436 - r4025487;
        double r4025489 = r4025482 * r4025488;
        double r4025490 = r4025460 * r4025460;
        double r4025491 = sqrt(r4025490);
        double r4025492 = r4025461 * r4025461;
        double r4025493 = sqrt(r4025492);
        double r4025494 = r4025491 * r4025493;
        double r4025495 = r4025489 / r4025494;
        double r4025496 = r4025436 / r4025492;
        double r4025497 = r4025484 * r4025483;
        double r4025498 = r4025496 * r4025497;
        double r4025499 = r4025448 * r4025498;
        double r4025500 = r4025437 * r4025499;
        double r4025501 = r4025461 * r4025444;
        double r4025502 = r4025500 / r4025501;
        double r4025503 = r4025436 - r4025502;
        double r4025504 = r4025475 * r4025475;
        double r4025505 = sqrt(r4025504);
        double r4025506 = r4025476 * r4025505;
        double r4025507 = r4025480 * r4025480;
        double r4025508 = sqrt(r4025507);
        double r4025509 = r4025481 * r4025508;
        double r4025510 = r4025506 * r4025509;
        double r4025511 = r4025503 * r4025510;
        double r4025512 = r4025473 ? r4025495 : r4025511;
        double r4025513 = r4025459 ? r4025471 : r4025512;
        return r4025513;
}

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 3 regimes

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

   1. Initial program 52.3

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

      \[\leadsto \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}{\color{blue}{\left(\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}\right) \cdot \sqrt[3]{\ell}}}\right)\]

   4. Applied add-cube-cbrt 52.4

      \[\leadsto \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{\color{blue}{\left(\sqrt[3]{h} \cdot \sqrt[3]{h}\right) \cdot \sqrt[3]{h}}}{\left(\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}\right) \cdot \sqrt[3]{\ell}}\right)\]

   5. Applied times-frac 52.4

      \[\leadsto \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 \color{blue}{\left(\frac{\sqrt[3]{h} \cdot \sqrt[3]{h}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{\sqrt[3]{h}}{\sqrt[3]{\ell}}\right)}\right)\]

   6. Applied associate-*r* 46.5

      \[\leadsto \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 - \color{blue}{\left(\left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{\sqrt[3]{h} \cdot \sqrt[3]{h}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right) \cdot \frac{\sqrt[3]{h}}{\sqrt[3]{\ell}}}\right)\]

   7. Simplified 34.0

      \[\leadsto \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 - \color{blue}{\left(\frac{1}{2} \cdot \left(\left(\left(\frac{D}{d} \cdot \frac{M}{2}\right) \cdot \frac{\sqrt[3]{h}}{\sqrt[3]{\ell}}\right) \cdot \left(\left(\frac{D}{d} \cdot \frac{M}{2}\right) \cdot \frac{\sqrt[3]{h}}{\sqrt[3]{\ell}}\right)\right)\right)} \cdot \frac{\sqrt[3]{h}}{\sqrt[3]{\ell}}\right)\]

   if -2.604184409542022e+233 < (* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l)))) < 6.854055542649234e+166

   1. Initial program 7.6

      \[\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 8.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 8.2

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

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

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

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

      \[\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 7.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(\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 7.4

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

       \[\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 2.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 2.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 2.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 frac-times 2.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{\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 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

   18. Applied sqrt-div 2.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(\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 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

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

   20. Applied frac-times 2.6

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

   21. Applied sqrt-div 2.6

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

   22. Applied associate-*l/ 2.6

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

   23. Applied frac-times 2.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{\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 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\]

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

   25. Simplified 1.9

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

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

   1. Initial program 50.8

      \[\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 50.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 50.9

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

      \[\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 37.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 37.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 37.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 37.4

       \[\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 37.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{\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 37.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 {\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 33.4

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

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

       \[\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 add-cube-cbrt 33.4

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

   18. Applied *-un-lft-identity 33.4

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

   19. Applied times-frac 33.4

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

   20. Applied associate-*r* 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(\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}{\left(\left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{1}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right) \cdot \frac{h}{\sqrt[3]{\ell}}}\right)\]

   21. Simplified 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(\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}{\left(\frac{1}{2} \cdot \left(\left(\frac{\frac{M \cdot D}{2}}{d} \cdot \frac{\frac{M \cdot D}{2}}{d}\right) \cdot \frac{1}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right)\right)} \cdot \frac{h}{\sqrt[3]{\ell}}\right)\]
   22. Using strategy rm

   23. Applied associate-*r/ 28.4

       \[\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 - \left(\frac{1}{2} \cdot \left(\color{blue}{\frac{\frac{\frac{M \cdot D}{2}}{d} \cdot \frac{M \cdot D}{2}}{d}} \cdot \frac{1}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right)\right) \cdot \frac{h}{\sqrt[3]{\ell}}\right)\]

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

   25. Applied associate-*r/ 26.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{\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{\frac{1}{2} \cdot \left(\left(\frac{\frac{M \cdot D}{2}}{d} \cdot \frac{M \cdot D}{2}\right) \cdot \frac{1}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right)}{d}} \cdot \frac{h}{\sqrt[3]{\ell}}\right)\]

   26. Applied frac-times 24.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(\frac{1}{2} \cdot \left(\left(\frac{\frac{M \cdot D}{2}}{d} \cdot \frac{M \cdot D}{2}\right) \cdot \frac{1}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right)\right) \cdot h}{d \cdot \sqrt[3]{\ell}}}\right)\]
3. Recombined 3 regimes into one program.

4. Final simplification 12.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 -2.604184409542022 \cdot 10^{+233}:\\ \;\;\;\;\left(1 - \frac{\sqrt[3]{h}}{\sqrt[3]{\ell}} \cdot \left(\frac{1}{2} \cdot \left(\left(\frac{\sqrt[3]{h}}{\sqrt[3]{\ell}} \cdot \left(\frac{D}{d} \cdot \frac{M}{2}\right)\right) \cdot \left(\frac{\sqrt[3]{h}}{\sqrt[3]{\ell}} \cdot \left(\frac{D}{d} \cdot \frac{M}{2}\right)\right)\right)\right)\right) \cdot \left({\left(\frac{d}{\ell}\right)}^{\frac{1}{2}} \cdot {\left(\frac{d}{h}\right)}^{\frac{1}{2}}\right)\\ \mathbf{elif}\;\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 6.854055542649234 \cdot 10^{+166}:\\ \;\;\;\;\frac{\left(\left(\left(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}} \cdot \left|\sqrt[3]{d}\right|\right) \cdot \left|\sqrt[3]{d}\right|\right) \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}}\right) \cdot \left(1 - \frac{h}{\ell} \cdot \left(\frac{1}{2} \cdot \left(\frac{\frac{M \cdot D}{2}}{d} \cdot \frac{\frac{M \cdot D}{2}}{d}\right)\right)\right)}{\sqrt{\sqrt[3]{h} \cdot \sqrt[3]{h}} \cdot \sqrt{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}\\ \mathbf{else}:\\ \;\;\;\;\left(1 - \frac{h \cdot \left(\frac{1}{2} \cdot \left(\frac{1}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\frac{\frac{M \cdot D}{2}}{d} \cdot \frac{M \cdot D}{2}\right)\right)\right)}{\sqrt[3]{\ell} \cdot d}\right) \cdot \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(\sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{\ell}} \cdot \frac{\sqrt[3]{d}}{\sqrt[3]{\ell}}}\right)\right)\\ \end{array}\]

Reproduce

herbie shell --seed 2019153 +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)))))