Henrywood and Agarwal, Equation (12)

Average Error: 25.6 → 13.2

Time: 57.1s

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

C

double f(double d, double h, double l, double M, double D) {
        double r8279074 = d;
        double r8279075 = h;
        double r8279076 = r8279074 / r8279075;
        double r8279077 = 1.0;
        double r8279078 = 2.0;
        double r8279079 = r8279077 / r8279078;
        double r8279080 = pow(r8279076, r8279079);
        double r8279081 = l;
        double r8279082 = r8279074 / r8279081;
        double r8279083 = pow(r8279082, r8279079);
        double r8279084 = r8279080 * r8279083;
        double r8279085 = M;
        double r8279086 = D;
        double r8279087 = r8279085 * r8279086;
        double r8279088 = r8279078 * r8279074;
        double r8279089 = r8279087 / r8279088;
        double r8279090 = pow(r8279089, r8279078);
        double r8279091 = r8279079 * r8279090;
        double r8279092 = r8279075 / r8279081;
        double r8279093 = r8279091 * r8279092;
        double r8279094 = r8279077 - r8279093;
        double r8279095 = r8279084 * r8279094;
        return r8279095;
}

↓

double f(double d, double h, double l, double M, double D) {
        double r8279096 = M;
        double r8279097 = D;
        double r8279098 = r8279096 * r8279097;
        double r8279099 = -2.1021038359939756e+79;
        bool r8279100 = r8279098 <= r8279099;
        double r8279101 = 1.0;
        double r8279102 = 0.5;
        double r8279103 = h;
        double r8279104 = cbrt(r8279103);
        double r8279105 = l;
        double r8279106 = cbrt(r8279105);
        double r8279107 = r8279104 / r8279106;
        double r8279108 = 2.0;
        double r8279109 = d;
        double r8279110 = r8279097 / r8279109;
        double r8279111 = r8279108 / r8279110;
        double r8279112 = r8279096 / r8279111;
        double r8279113 = r8279107 * r8279112;
        double r8279114 = r8279113 * r8279113;
        double r8279115 = r8279102 * r8279114;
        double r8279116 = r8279115 * r8279107;
        double r8279117 = r8279101 - r8279116;
        double r8279118 = r8279109 / r8279105;
        double r8279119 = pow(r8279118, r8279102);
        double r8279120 = cbrt(r8279109);
        double r8279121 = r8279120 / r8279104;
        double r8279122 = fabs(r8279121);
        double r8279123 = sqrt(r8279121);
        double r8279124 = r8279122 * r8279123;
        double r8279125 = r8279119 * r8279124;
        double r8279126 = r8279117 * r8279125;
        double r8279127 = 1.3396067446246392e-79;
        bool r8279128 = r8279098 <= r8279127;
        double r8279129 = cbrt(r8279104);
        double r8279130 = r8279129 / r8279105;
        double r8279131 = r8279096 * r8279104;
        double r8279132 = r8279097 / r8279108;
        double r8279133 = r8279109 / r8279132;
        double r8279134 = r8279131 / r8279133;
        double r8279135 = r8279129 * r8279134;
        double r8279136 = r8279135 * r8279135;
        double r8279137 = r8279102 * r8279136;
        double r8279138 = r8279130 * r8279137;
        double r8279139 = r8279101 - r8279138;
        double r8279140 = r8279109 / r8279106;
        double r8279141 = sqrt(r8279140);
        double r8279142 = r8279101 / r8279106;
        double r8279143 = r8279142 / r8279106;
        double r8279144 = sqrt(r8279143);
        double r8279145 = r8279141 * r8279144;
        double r8279146 = r8279145 * r8279124;
        double r8279147 = r8279139 * r8279146;
        double r8279148 = 1.5623144389475482e+166;
        bool r8279149 = r8279098 <= r8279148;
        double r8279150 = cbrt(r8279141);
        double r8279151 = r8279150 * r8279150;
        double r8279152 = r8279150 * r8279151;
        double r8279153 = r8279152 * r8279144;
        double r8279154 = r8279124 * r8279153;
        double r8279155 = r8279103 / r8279105;
        double r8279156 = r8279108 * r8279109;
        double r8279157 = r8279098 / r8279156;
        double r8279158 = pow(r8279157, r8279108);
        double r8279159 = r8279158 * r8279102;
        double r8279160 = r8279155 * r8279159;
        double r8279161 = r8279101 - r8279160;
        double r8279162 = r8279154 * r8279161;
        double r8279163 = r8279096 / r8279156;
        double r8279164 = r8279163 * r8279097;
        double r8279165 = r8279164 * r8279155;
        double r8279166 = r8279164 * r8279165;
        double r8279167 = r8279102 * r8279166;
        double r8279168 = r8279101 - r8279167;
        double r8279169 = r8279109 / r8279103;
        double r8279170 = pow(r8279169, r8279102);
        double r8279171 = r8279119 * r8279170;
        double r8279172 = r8279168 * r8279171;
        double r8279173 = r8279149 ? r8279162 : r8279172;
        double r8279174 = r8279128 ? r8279147 : r8279173;
        double r8279175 = r8279100 ? r8279126 : r8279174;
        return r8279175;
}

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

2. if (* M D) < -2.1021038359939756e+79

   1. Initial program 35.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 35.6

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

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

      \[\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 29.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 28.0

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

   8. Simplified 28.0

      \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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 28.0

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

   11. Applied add-cube-cbrt 28.0

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

   12. Applied times-frac 28.0

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

   13. Applied associate-*r* 27.7

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

   14. Simplified 18.3

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

   if -2.1021038359939756e+79 < (* M D) < 1.3396067446246392e-79

   1. Initial program 22.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 22.6

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

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

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

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

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

   8. Simplified 17.4

      \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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 17.5

       \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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 17.5

       \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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 17.4

       \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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 13.0

       \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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 13.0

       \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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 13.0

       \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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 *-un-lft-identity 13.0

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

   18. Applied add-cube-cbrt 13.0

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

   19. Applied times-frac 13.0

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

   20. Applied associate-*r* 9.8

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

   21. Simplified 9.7

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

   23. Applied *-un-lft-identity 9.7

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

   24. Applied add-cube-cbrt 9.8

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

   25. Applied times-frac 9.8

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

   26. Applied associate-*r* 9.1

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

   27. Simplified 9.2

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

   if 1.3396067446246392e-79 < (* M D) < 1.5623144389475482e+166

   1. Initial program 29.5

      \[\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 29.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 29.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 29.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 24.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 24.1

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

   8. Simplified 24.1

      \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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.2

       \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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.2

       \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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.2

       \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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 20.4

       \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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 20.4

       \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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 20.4

       \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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 20.7

       \[\leadsto \left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \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(\left(\sqrt[3]{\sqrt{\frac{d}{\sqrt[3]{\ell}}}} \cdot \sqrt[3]{\sqrt{\frac{d}{\sqrt[3]{\ell}}}}\right) \cdot \sqrt[3]{\sqrt{\frac{d}{\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 1.5623144389475482e+166 < (* M D)

   1. Initial program 40.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 associate-*l* 40.3

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

   4. Simplified 33.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 - \frac{1}{2} \cdot \color{blue}{\left(\left(\frac{h}{\ell} \cdot \left(\frac{M}{d \cdot 2} \cdot D\right)\right) \cdot \left(\frac{M}{d \cdot 2} \cdot D\right)\right)}\right)\]
3. Recombined 4 regimes into one program.

4. Final simplification 13.2

   \[\leadsto \begin{array}{l} \mathbf{if}\;M \cdot D \le -2.1021038359939756 \cdot 10^{+79}:\\ \;\;\;\;\left(1 - \left(\frac{1}{2} \cdot \left(\left(\frac{\sqrt[3]{h}}{\sqrt[3]{\ell}} \cdot \frac{M}{\frac{2}{\frac{D}{d}}}\right) \cdot \left(\frac{\sqrt[3]{h}}{\sqrt[3]{\ell}} \cdot \frac{M}{\frac{2}{\frac{D}{d}}}\right)\right)\right) \cdot \frac{\sqrt[3]{h}}{\sqrt[3]{\ell}}\right) \cdot \left({\left(\frac{d}{\ell}\right)}^{\frac{1}{2}} \cdot \left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right)\right)\\ \mathbf{elif}\;M \cdot D \le 1.3396067446246392 \cdot 10^{-79}:\\ \;\;\;\;\left(1 - \frac{\sqrt[3]{\sqrt[3]{h}}}{\ell} \cdot \left(\frac{1}{2} \cdot \left(\left(\sqrt[3]{\sqrt[3]{h}} \cdot \frac{M \cdot \sqrt[3]{h}}{\frac{d}{\frac{D}{2}}}\right) \cdot \left(\sqrt[3]{\sqrt[3]{h}} \cdot \frac{M \cdot \sqrt[3]{h}}{\frac{d}{\frac{D}{2}}}\right)\right)\right)\right) \cdot \left(\left(\sqrt{\frac{d}{\sqrt[3]{\ell}}} \cdot \sqrt{\frac{\frac{1}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}}\right) \cdot \left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right)\right)\\ \mathbf{elif}\;M \cdot D \le 1.5623144389475482 \cdot 10^{+166}:\\ \;\;\;\;\left(\left(\left|\frac{\sqrt[3]{d}}{\sqrt[3]{h}}\right| \cdot \sqrt{\frac{\sqrt[3]{d}}{\sqrt[3]{h}}}\right) \cdot \left(\left(\sqrt[3]{\sqrt{\frac{d}{\sqrt[3]{\ell}}}} \cdot \left(\sqrt[3]{\sqrt{\frac{d}{\sqrt[3]{\ell}}}} \cdot \sqrt[3]{\sqrt{\frac{d}{\sqrt[3]{\ell}}}}\right)\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}:\\ \;\;\;\;\left(1 - \frac{1}{2} \cdot \left(\left(\frac{M}{2 \cdot d} \cdot D\right) \cdot \left(\left(\frac{M}{2 \cdot d} \cdot D\right) \cdot \frac{h}{\ell}\right)\right)\right) \cdot \left({\left(\frac{d}{\ell}\right)}^{\frac{1}{2}} \cdot {\left(\frac{d}{h}\right)}^{\frac{1}{2}}\right)\\ \end{array}\]

Reproduce

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