At first, we sample in the ( is odd) equidistant points around :
where is some step.
Then we interpolate points by polynomial
Its coefficients are found as a solution of system of linear equations:
Here are references to existing equations: (1), (2).
Here is reference to non-existing equation (??).