False Position Linear Interpolator
Calculates chord roots c = [a·f(b) − b·f(a)] / [f(b) − f(a)] connecting interval endpoints.
Determine the root of f(x) = x³ − 5x + 1 = 0 on [0.0, 1.0] using the classical Regula Falsi (false position) method, investigate one-sided endpoint stagnation, evaluate the Illinois acceleration algorithm, and extract the superlinear convergence order α = 1.442 via 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.
Classical Regula Falsi can suffer from severe endpoint stagnation, converging much slower than bisection. The numerical simulation model evaluates the root of f(x) = x³ − 5x + 1 = 0 on [a₀, b₀] = [0.000, 1.000]. Initial values: f(0.0) = +1.000, f(1.0) = −3.000. True root: x* = 0.201637667114... Using the Illinois algorithm: step 1 yields c₁ = 0.250000 (ε₁ = 0.048362); step 2 detects stagnation on a, halves f(a)=0.5, yielding c₂ = 0.203815 (ε₂ = 0.002177); step 3 yields c₃ = 0.2016421 (ε₃ = 4.43 × 10⁻⁶); step 4 yields c₄ = 0.2016376672 (ε₄ = 1.15 × 10⁻¹⁰). Logarithmic coordinates: x = ln(ε_k) = [−3.029, −6.129, −12.327] and y = ln(ε_{k+1}) = [−6.129, −12.327, −22.887]. Ordinary least-squares linear regression on these coordinates yields an empirical slope α = 1.44225 and Y-intercept ln(C) = −1.758 with Pearson correlation coefficient R² = 0.99999998. The empirical order α = 1.442 matches the theoretical Illinois superlinear convergence order α = 3^(1/3) = 1.44225 within 0.003% relative error.
Calculates chord roots c = [a·f(b) − b·f(a)] / [f(b) − f(a)] connecting interval endpoints.
Monitors endpoint retention and halves the function value of stagnant bounds.
Guarantees that f(a_k) · f(b_k) < 0 is strictly maintained at every step.
Linearizes log-errors to measure the transition from linear to superlinear order.
Visualizes secant chords across [a_k, b_k] and demonstrates how Illinois breaks stagnation.
Supplies exact analytical roots for error quantification.
Determine the root of f(x) = x³ − 5x + 1 = 0 on [0.0, 1.0] using the classical Regula Falsi (false position) method, investigate one-sided endpoint stagnation, evaluate the Illinois acceleration algorithm, and extract the superlinear convergence order α = 1.442 via linear regression. Standardized input settings: 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 steps.
Evaluate f(0.0) = +1.000 > 0 and f(1.0) = −3.000 < 0. Bolzano's Theorem guarantees a root in [0.0, 1.0].
Compute c₁ = [0·(−3) − 1·(1)] / [−3 − 1] = −1 / −4 = 0.250000. Evaluate f(c₁) = (0.25)³ − 5(0.25) + 1 = 0.015625 − 1.25 + 1 = −0.234375 < 0. Update [a₁, b₁] = [0.0, 0.25].
Observe that if classical Regula Falsi is applied, the convex curvature of f(x) causes left bound a = 0 to remain fixed while b crawls slowly toward x*, resulting in sluggish linear convergence.
Apply the Illinois algorithm: because the right endpoint replaced b twice consecutively, replace f(a) with f(a)/2 = 0.500. Compute new chord root c₂ = [0·(−0.2344) − 0.25·(0.5)] / [−0.2344 − 0.5] = −0.125 / −0.7344 = 0.2038.
Subsequent Illinois iterations yield c₃ = 0.201642 and c₄ = 0.20163767, rapidly accelerating error reduction.
Perform linear regression of ln(ε_{k+1}) against ln(ε_k). Extract empirical order α = 1.442, matching the theoretical Illinois superlinear convergence order α = 3^(1/3) ≈ 1.442.
Points are your recorded observations. The solid regression line validates Regula Falsi Method mathematical scaling.
| # | Iteration Step k (steps) | False Position Iterate c_k () | 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.