Recursive Difference Engine
Evaluates zero-th through third-order divided differences via two-point quotient recurrences.
Implement the recursive Newton divided difference tableau, compute higher-order difference quotients, and verify the connection between the leading divided difference f[x_0, ..., x_n] and the n-th derivative of the underlying function. Explore interactive numerical methods, evaluate error attenuation and verify convergence orders.
Ready to execute Step = 1.
Adjust Step, then execute computation.
Mathematical operations are computed deterministically using 64-bit IEEE 754 precision. Residuals reflect ordinary least-squares regression across successive observations.
Like all polynomial interpolation on equidistant nodes, high-degree Newton polynomials eventually diverge near boundaries due to Runge's phenomenon. The simulation evaluates intermediate and leading divided differences for the cubic test function across six successive observation trials. The linear regression fit of the response y = c_3 · k verifies the theoretical leading coefficient c_3 = 0.1667 (1/6) with R² = 1.000000.
Evaluates zero-th through third-order divided differences via two-point quotient recurrences.
Maintains the lower triangular tableau matrix of intermediate divided differences without redundant operations.
Evaluates the Newton polynomial form P_n(x) = c_0 + c_1(x-x_0) + c_2(x-x_0)(x-x_1) + ... in O(n) operations.
Compares the highest divided difference f[x_0, ..., x_n] against the theoretical derivative value f^{(n)}(ξ)/n!.
Appends new observation nodes to the existing interpolant without recalculating previously computed table rows.
Monitors the cumulative quotient progression across steps to fit the invariant leading coefficient c_3.
Implement the recursive Newton divided difference tableau, compute higher-order difference quotients, and verify the connection between the leading divided difference f[x_0, ..., x_n] and the n-th derivative of the underlying function. Standardized input settings: 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 step.
Set four distinct nodes x_0, x_1, x_2, x_3 and compute the base function values f[x_i] = f(x_i) in the first column of the tableau.
Compute first-order differences f[x_i, x_{i+1}] = (f[x_{i+1}] - f[x_i]) / (x_{i+1} - x_i) for all adjacent pairs.
Iteratively construct second and third-order columns using the recurrence f[x_i, ..., x_{i+k}] = (f[x_{i+1}, ..., x_{i+k}] - f[x_i, ..., x_{i+k-1}]) / (x_{i+k} - x_i).
Read the top diagonal entries c_k = f[x_0, ..., x_k], which constitute the coefficients of the Newton basis polynomials.
Confirm that for the test cubic polynomial, the third divided difference c_3 = 1/6 ≈ 0.1667 remains constant regardless of node placement.
Log the cumulative quotient response across six observation settings and perform linear regression to verify the theoretical slope c_3 = 0.1667.
Points are your recorded observations. The solid regression line validates Newton Divided Differences mathematical scaling.
| # | Observation Step (k) (step) | Empirical Difference Quotient Value (quotient) | Observation Step Index k | Cumulative Quotient Metric y = c_3 · k | Remove |
|---|
No readings yet. Execute a computation and record it.