At first, we sample

in the

(

is odd) equidistant points around

:
![Rendered by QuickLaTeX.com \[ f_k = f(x_k),\: x_k = x^*+kh,\: k=-\frac{N-1}{2},\dots,\frac{N-1}{2} \]](https://mathgraphs.ru/wp-content/ql-cache/quicklatex.com-ba5d0fab0e07a4362dd6a77cbf5d2e65_l3.png)
where

is some step.
Then we interpolate points

by polynomial
(1) 
Its coefficients

are found as a solution of system of linear equations:
(2) 
Here are references to existing equations: (
1), (
2).
Here is reference to non-existing equation (??).