At first, we sample in the ( is odd) equidistant points around :

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 (??).