We define recursively finite differences $ g_k (x) $ of order $ k $ of function $ f $ as follows: $g_0(x)=f(x)$, $g_n(x)=g_{n-1}(x+h_n)-g_{n-1}(x) (n\in\mathbb{N})$.
It is known that all arguments of the function $ f $ (namely, $ x_0, x_0 + h_1, ... $) used in determining $ g_k (x_0) $ lie on the segment $ [a, b] \subset \mathbb {R} $ and the length of this segment is minimally possible, and also that $f \in C ^ {(n-1)} [a, b] $ and there is $ f ^ {(n)} (x) $ at least in the interval $ (a, b) $.
It is necessary to prove that then $\forall k\in \mathbb {N}$ there is a point $ \xi \in [a, b] $ such that $ g_k (x_0) = f ^ {(k)} (\xi) h_1 \cdot h_2 \cdot ... \cdot h_k $. I suppose we should use Cauchy's mean value theorem, but do not have any idea how to do it