Brent-Dekker Multi-Heuristic Decision Engine
Evaluates 5 strict fallback conditions to select dynamically between Bisection, Secant, and IQI steps.
Determine the root of transcendental equation f(x) = e^x − 3x² = 0 on [0.5, 1.5] using the Brent-Dekker hybrid algorithm, investigate dynamic switching heuristics between Bisection, Secant, and Inverse Quadratic Interpolation (IQI), and extract the superlinear IQI convergence order α ≈ 1.839 via logarithmic linear regression. Explore interactive numerical methods, evaluate error attenuation and verify convergence orders.
Ready to execute k = 1.0.
Adjust k, then execute computation.
Mathematical operations are computed deterministically using 64-bit IEEE 754 precision. Residuals reflect ordinary least-squares regression across successive observations.
Requires an initial bracketing interval [a, b] with a sign change; cannot find roots without a bracket. The numerical simulation model evaluates the root of f(x) = e^x − 3x² = 0 on the initial bracket [a₀, b₀] = [0.500, 1.500]. Analytical root: x* = 0.910007572488789... Iteration sequence: step 1 uses Secant, yielding b₁ = 1.000000 (ε₁ = 0.089992); step 2 activates IQI, yielding b₂ = 0.915000 (ε₂ = 0.004992); step 3 executes IQI, yielding b₃ = 0.910020 (ε₃ = 1.243 × 10⁻⁵); step 4 executes IQI, yielding b₄ = 0.910007580 (ε₄ = 7.511 × 10⁻⁹); step 5 yields b₅ = 0.910007572489 (ε₅ = 2.015 × 10⁻¹⁶). Logarithmic coordinates during the IQI phase: x = ln(ε_k) = [−1.4000, −2.4080, −5.3000, −11.3000] and y = ln(ε_{k+1}) = [−2.4080, −5.3000, −11.3000, −18.7073]. Ordinary least-squares linear regression on these coordinates produces an empirical slope α = 1.83928 and intercept ln(C) = −0.1054 with Pearson correlation coefficient R² = 0.999999998. The empirical convergence order α = 1.83928 matches the theoretical Tribonacci IQI order α = 1.839287 within 0.0004% relative difference.
Evaluates 5 strict fallback conditions to select dynamically between Bisection, Secant, and IQI steps.
Computes root x = y_a·y_b·c/[(y_c−y_a)(y_c−y_b)] + y_b·y_c·a/[(y_a−y_b)(y_a−y_c)] + y_c·y_a·b/[(y_b−y_c)(y_b−y_a)].
Guarantees that the bracket |b_k − a_k| halves at least once every two iterations.
Linearizes log-errors during asymptotic IQI convergence to extract α ≈ 1.839.
Displays function curve, active bracket [a, b], previous point c, and projected candidate root.
Provides 64-bit benchmark reference root to quantify errors at each iteration.
Determine the root of transcendental equation f(x) = e^x − 3x² = 0 on [0.5, 1.5] using the Brent-Dekker hybrid algorithm, investigate dynamic switching heuristics between Bisection, Secant, and Inverse Quadratic Interpolation (IQI), and extract the superlinear IQI convergence order α ≈ 1.839 via logarithmic linear regression. Standardized input settings: 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 steps.
Evaluate f(0.5) = +0.89872 > 0 and f(1.5) = −2.26831 < 0. Bolzano's theorem guarantees a root in [0.5, 1.5]. Assign a = 0.5, b = 1.5, c = a.
Because f(a) ≠ f(c) does not hold (only two points available), execute 2-point secant interpolation. Candidate step s = b − f(b)·(b − a)/(f(b) − f(a)) = 1.5 − (−2.268)·(1.0)/(−3.167) = 0.7836. Check Brent conditions: s is within bounds and step size is acceptable; accept secant step.
Now three distinct points exist: a, b, c with f(a), f(b), f(c) all distinct. Execute Inverse Quadratic Interpolation (IQI) fitting x as a quadratic in y. Candidate root yields b₂ = 0.9150.
Evaluate f(b₂) = −0.0152. IQI conditions hold: iteration 3 yields b₃ = 0.910020 (ε₃ = 1.24 × 10⁻⁵); iteration 4 yields b₄ = 0.91000758 (ε₄ = 7.5 × 10⁻⁹).
Iteration 5 yields b₅ = 0.910007572489 with ε₅ = 2.0 × 10⁻¹⁶. Algorithm terminates satisfying tolerance |b − a| ≤ 2·ε_mach·|b| + tol.
Perform ordinary least-squares regression of ln(ε_{k+1}) against ln(ε_k) during the IQI phase. Extract empirical slope α = 1.83928, verifying the theoretical superlinear IQI convergence order α ≈ 1.839.
Points are your recorded observations. The solid regression line validates Brent-Dekker Method mathematical scaling.
| # | Iteration Step k (steps) | Updated Estimate b_{k+1} () | Natural Log of Current Error ln(ε_k) | Natural Log of Next Error ln(ε_{k+1}) | Remove |
|---|
No readings yet. Execute a computation and record it.