Average Error: 0.2 → 0.2
Time: 6.4s
Precision: 64
\[\left|\frac{1}{\sqrt{\pi}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right|\]
↓
\[\left|\frac{1}{\sqrt{\pi}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left|{\left(\left|x\right|\right)}^{\frac{3}{2}}\right| \cdot \left|{\left(\left|x\right|\right)}^{\frac{3}{2}}\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right|\]
\left|\frac{1}{\sqrt{\pi}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right|
↓
\left|\frac{1}{\sqrt{\pi}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left|{\left(\left|x\right|\right)}^{\frac{3}{2}}\right| \cdot \left|{\left(\left|x\right|\right)}^{\frac{3}{2}}\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right|
double f(double x) {
        double r103121 = 1.0;
        double r103122 = atan2(1.0, 0.0);
        double r103123 = sqrt(r103122);
        double r103124 = r103121 / r103123;
        double r103125 = 2.0;
        double r103126 = x;
        double r103127 = fabs(r103126);
        double r103128 = r103125 * r103127;
        double r103129 = 3.0;
        double r103130 = r103125 / r103129;
        double r103131 = r103127 * r103127;
        double r103132 = r103131 * r103127;
        double r103133 = r103130 * r103132;
        double r103134 = r103128 + r103133;
        double r103135 = 5.0;
        double r103136 = r103121 / r103135;
        double r103137 = r103132 * r103127;
        double r103138 = r103137 * r103127;
        double r103139 = r103136 * r103138;
        double r103140 = r103134 + r103139;
        double r103141 = 21.0;
        double r103142 = r103121 / r103141;
        double r103143 = r103138 * r103127;
        double r103144 = r103143 * r103127;
        double r103145 = r103142 * r103144;
        double r103146 = r103140 + r103145;
        double r103147 = r103124 * r103146;
        double r103148 = fabs(r103147);
        return r103148;
}
↓
double f(double x) {
        double r103149 = 1.0;
        double r103150 = atan2(1.0, 0.0);
        double r103151 = sqrt(r103150);
        double r103152 = r103149 / r103151;
        double r103153 = 2.0;
        double r103154 = x;
        double r103155 = fabs(r103154);
        double r103156 = r103153 * r103155;
        double r103157 = 3.0;
        double r103158 = r103153 / r103157;
        double r103159 = r103155 * r103155;
        double r103160 = r103159 * r103155;
        double r103161 = r103158 * r103160;
        double r103162 = r103156 + r103161;
        double r103163 = 5.0;
        double r103164 = r103149 / r103163;
        double r103165 = r103160 * r103155;
        double r103166 = r103165 * r103155;
        double r103167 = r103164 * r103166;
        double r103168 = r103162 + r103167;
        double r103169 = 21.0;
        double r103170 = r103149 / r103169;
        double r103171 = 1.5;
        double r103172 = pow(r103155, r103171);
        double r103173 = fabs(r103172);
        double r103174 = r103173 * r103173;
        double r103175 = r103174 * r103155;
        double r103176 = r103175 * r103155;
        double r103177 = r103176 * r103155;
        double r103178 = r103177 * r103155;
        double r103179 = r103170 * r103178;
        double r103180 = r103168 + r103179;
        double r103181 = r103152 * r103180;
        double r103182 = fabs(r103181);
        return r103182;
}

Error

Bits error versus x

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Initial program 0.2

    \[\left|\frac{1}{\sqrt{\pi}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right|\]
  2. Using strategy rm

  3. Applied add-sqr-sqrt0.2

    \[\leadsto \left|\frac{1}{\sqrt{\pi}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\color{blue}{\left(\sqrt{\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|} \cdot \sqrt{\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|}\right)} \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right|\]

  4. Simplified0.2

    \[\leadsto \left|\frac{1}{\sqrt{\pi}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\color{blue}{\left|{\left(\left|x\right|\right)}^{\frac{3}{2}}\right|} \cdot \sqrt{\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|}\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right|\]

  5. Simplified0.2

    \[\leadsto \left|\frac{1}{\sqrt{\pi}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left|{\left(\left|x\right|\right)}^{\frac{3}{2}}\right| \cdot \color{blue}{\left|{\left(\left|x\right|\right)}^{\frac{3}{2}}\right|}\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right|\]

  6. Final simplification0.2

    \[\leadsto \left|\frac{1}{\sqrt{\pi}} \cdot \left(\left(\left(2 \cdot \left|x\right| + \frac{2}{3} \cdot \left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{5} \cdot \left(\left(\left(\left(\left|x\right| \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right) + \frac{1}{21} \cdot \left(\left(\left(\left(\left(\left|{\left(\left|x\right|\right)}^{\frac{3}{2}}\right| \cdot \left|{\left(\left|x\right|\right)}^{\frac{3}{2}}\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right) \cdot \left|x\right|\right)\right)\right|\]

Reproduce

herbie shell --seed 2020036 +o rules:numerics
(FPCore (x)
  :name "Jmat.Real.erfi, branch x less than or equal to 0.5"
  :precision binary64
  (fabs (* (/ 1 (sqrt PI)) (+ (+ (+ (* 2 (fabs x)) (* (/ 2 3) (* (* (fabs x) (fabs x)) (fabs x)))) (* (/ 1 5) (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)))) (* (/ 1 21) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))))))