Newton's interpolation formula and sums of powers

Abstract In this manuscript, we derive closed formulas for multifold sums of powers of integers by combining the forward Newton interpolation formula with hockey-stick identities for binomial coefficients. We further obtain representations of multifold sums of powers in terms of Stirling numbers of the second kind and Eulerian numbers. Finally, we provide Wolfram Mathematica programs for the efficient verification of the derived identities. Related works Newton's interpolation formula and sums of powers (2025) Sums of powers via central finite differences and Newton's formula (2025) Sums of powers via backward finite differences and Newton's formula (2026) Sums of powers of integers: A complete framework for closed formulas (2026) OEIS https://oeis.org/A131689 — Triangle of numbers T(n,k) = k!*Stirling2(n,k) read by rows, T(n, k) for 0