2cbrt (problem 3.3.4)

Average Error: 29.9 → 0.5

Time: 5.1s

Precision: binary64

Cost: 32896

Math

\[\sqrt[3]{x + 1} - \sqrt[3]{x} \]

↓

\[\begin{array}{l} t_0 := \sqrt[3]{1 + x}\\ \frac{1}{{t_0}^{2} + \sqrt[3]{x} \cdot \left(t_0 + \sqrt[3]{x}\right)} \end{array} \]

FPCore

(FPCore (x) :precision binary64 (- (cbrt (+ x 1.0)) (cbrt x)))

↓

(FPCore (x)
 :precision binary64
 (let* ((t_0 (cbrt (+ 1.0 x))))
   (/ 1.0 (+ (pow t_0 2.0) (* (cbrt x) (+ t_0 (cbrt x)))))))

C

double code(double x) {
	return cbrt((x + 1.0)) - cbrt(x);
}

↓

double code(double x) {
	double t_0 = cbrt((1.0 + x));
	return 1.0 / (pow(t_0, 2.0) + (cbrt(x) * (t_0 + cbrt(x))));
}

Java

public static double code(double x) {
	return Math.cbrt((x + 1.0)) - Math.cbrt(x);
}

↓

public static double code(double x) {
	double t_0 = Math.cbrt((1.0 + x));
	return 1.0 / (Math.pow(t_0, 2.0) + (Math.cbrt(x) * (t_0 + Math.cbrt(x))));
}

Julia

function code(x)
	return Float64(cbrt(Float64(x + 1.0)) - cbrt(x))
end

↓

function code(x)
	t_0 = cbrt(Float64(1.0 + x))
	return Float64(1.0 / Float64((t_0 ^ 2.0) + Float64(cbrt(x) * Float64(t_0 + cbrt(x)))))
end

Wolfram

code[x_] := N[(N[Power[N[(x + 1.0), $MachinePrecision], 1/3], $MachinePrecision] - N[Power[x, 1/3], $MachinePrecision]), $MachinePrecision]

↓

code[x_] := Block[{t$95$0 = N[Power[N[(1.0 + x), $MachinePrecision], 1/3], $MachinePrecision]}, N[(1.0 / N[(N[Power[t$95$0, 2.0], $MachinePrecision] + N[(N[Power[x, 1/3], $MachinePrecision] * N[(t$95$0 + N[Power[x, 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Initial program 29.9

   \[\sqrt[3]{x + 1} - \sqrt[3]{x} \]

2. Applied egg-rr 29.3

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

3. Taylor expanded in x around 0 0.5

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

4. Final simplification 0.5

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

Alternatives

Alternative 1
Error	15.4
Cost	26624

\[\frac{1}{{\left(\sqrt[3]{1 + x}\right)}^{2} + \left(\sqrt[3]{x \cdot x} + \sqrt[3]{x \cdot \left(1 + x\right)}\right)} \]

Alternative 2
Error	29.9
Cost	25984

\[\sqrt{{\left(\sqrt[3]{1 + x} - \sqrt[3]{x}\right)}^{2}} \]

Alternative 3
Error	29.9
Cost	25920

\[\log \left(e^{\sqrt[3]{1 + x} - \sqrt[3]{x}}\right) \]

Alternative 4
Error	29.9
Cost	13120

\[\sqrt[3]{1 + x} - \sqrt[3]{x} \]

Alternative 5
Error	61.7
Cost	64

\[0 \]

Alternative 6
Error	31.9
Cost	64

\[1 \]

Reproduce

herbie shell --seed 2022313 
(FPCore (x)
  :name "2cbrt (problem 3.3.4)"
  :precision binary64
  (- (cbrt (+ x 1.0)) (cbrt x)))