This computation of does not explain how to compute a derivative
approximation. But if the point **b** is extended to a series of
uniformly spaced points
each a multiple of **h** from **a** and we write the TS expansion as an
expression, we get

Assume the difference operator is centered at the origin. Now formulating this as a matrix equation using terms up to the second derivative, we have